1. 引言
地震资料中同向轴的变化与地层岩性的变化有着密切关系，岩性的变化与沉积相的分布相关联，因此根据地震资料提取地震属性划分地震相对沉积相的分布具有一定的指导意义 [1] [2] [3] 。过去地震相的划分主要通过“相面法”对地震反射特征进行归类总结，但单一的地震属性只能反应某种特定的沉积相的对应特征，以此为依据划分沉积相往往具有多解性和不确定性，将多种属性融合分析可以更为立体地综合反应沉积相的地震反射特征，从而更加精确地建立地震相与沉积相之间的对应关系 [4] 。
准噶尔盆地玛湖一区有多口井钻至二叠系风城组，区域内主要有扇三角洲平原、扇三角洲内前缘、扇三角洲外前缘三个主要亚相，钻遇不同沉积相的井油气显示差别较大，分布于前缘相带的井往往具有良好的开发价值或油气显示。但目前对于研究区及周边地区二叠系风城组的相带边界划分及地震相与沉积相的对应关系研究较少，本文借助钻井较多的玛湖一区，利用Geoeast软件的属性分析模块，分析不同方法划分地震相的效果，建立研究区周边区域地震相和沉积相对应模板，为该区域的沉积相划分工作提供思路。
2. 区域地质概况
研究区位于准噶尔盆地西部中拐凸起的东北坡，整块工区基本在玛湖凹陷里(图1)。研究区经历多期构造旋回影响，发生了多次抬升剥蚀作用，最终在白垩纪基本定型，形成了如今北西高，向南东逐渐变为缓坡的斜坡构造模式。
该区域地层保存较为完整，自石炭系至新近系均有发育，其中石炭系不整合及白垩系不整合为区域内较为明显的标志层，研究区二叠系自下而上发育风城组(P1f)、夏子街组(P2x)，其中风城组分为风一段(P1f1)、风二段(P1f2)、风三段(P1f3)三个亚段，风一段地层厚度向斜坡上倾方向逐渐减薄，风城组地层呈现出明显的“下削上超”的特征。研究区内坡度整体呈现西高东低的形态，发育大型缓坡浅水三角洲相沉积体系，接受中拐扇的物源供给，同时砂体的分布受到古地貌的控制古地貌低点为砂体卸载和前缘相发育的有利区，同时海西期形成的北西——南东向逆断层延伸距离远、断距大，为风城组油气聚集提供有利通道，印支期形成的大部分为断距较小的正断层，延伸距离短，一定程度上调整了油藏的分布，使得风城组具有较好的油气富集条件，具有较大的勘探前景 [5] 。
Figure 1. Structural position diagram of the study area
图1. 研究区构造位置图
3. 玛湖一区地震相模式识别
常用的地震相属性分析包括：振幅、频率、几何、相位属性等 [6] ，从单属性分析的角度出发，振幅及频率属性大致可以在范围内对地震相特征进行划分，对研究区风城组二段底界向上40 ms为时窗提取瞬时振幅属性及瞬时频率属性(图2、图3)，通过剖面对比分析可以发现：研究区地震反射特征可以归为以下几类。
Figure 2. 40 ms root mean square amplitude slice on the Fengcheng formation
图2. 风城组上40 ms均方根振幅切片
Figure 3. 40 ms instantaneous frequency section on the Fengcheng formation
图3. 风城组上40 ms瞬时频率切片
I-I类地震相是玛湖1区块东边附近的振幅较强、频率低、连续性较好地震相，对应滨浅湖的地震反射特征，通过MAH48井剖面可以发现，在分界处的地震反射特征差别较大，玛湖48井右侧三角洲外前缘表现为振幅多为中低值、频率较低、连续反射的地震反射特征，MAH48井三角洲平原–内前缘相表现为振幅值较高、频率高、断续反射的地震反射特征(图4)。
II-II类地震相表现为振幅中值、频率低、较连续反射地震相，三角洲外前缘反射表现为地层向南东逐渐增厚并出现上超现象，通过MAH48井上可以发现在分界处的地震反射特征差别较大。KE207井附近可见部分断续波谷反射，向左地震反射频率逐渐减小、振幅增强，且振幅强弱横向变化，体现出三角洲前缘的沉积特征(图5)。
Figure 4. Deentary entary section of well MH48 (Section positions are shown in Figure 3)
图4. 过MH48井沉积相解释剖面(剖面位置见图3)
II-III类地震相表现为频率较低、振幅中值的断续的地震反射特征，部分地区表现为中频杂乱反射，对应三角洲平原–内前缘反射。
III-II类地震相表现为频率较低、振幅较强的连续地震反射特征，对应三角洲外前缘反射。
Figure 5. Interpretation section of well MAH 26-KE 207 (Section positions are shown in Figure 3)
图5. 过MAH26-KE207井沉积相解释剖面(剖面位置见图3)
III-II类地震相表现为杂乱反射，与三角洲平原–内前缘强水动力及砾岩沉积相对应，MAH26井中上部呈现出明显的砾岩沉积的杂乱空白反射，向KE207井延伸过程中随着搬运距离增加，岩石粒度降低，分选变好，在三角洲外前缘沉积砂砾岩，地震反射特征与连井测井沉积相解释结果较为符合(图5)。
Figure 6. Interpretation section of well KE81-MAH28 (Section positions are shown in Figure 3)
图6. 过KE81-MAH28井沉积相解释剖面(剖面位置见图3)
IV-III类地震相表现为频率高值、振幅较强的杂乱的地震反射特征，对应三角洲平原–内前缘反射(图6)。
通过井震对比分析，并结合振幅及频率属性切片可建立研究区六类地震相与沉积相的对应关系，认为I-I类地震相代表了滨浅湖反射特征，II-II、II-III类地震相为三角洲外前缘反射特征，III-II、III-II、IV-III类地震相为三角洲平原–内前缘反射特征，表明该区块风城组二段振幅及频率属性与沉积相分布具有一定的对应关系，利用振幅及频率属性特征可以在一定程度上预测沉积相的分布 [7] 。
4. 纹理属性分析
Table 1. Principle and meaning table of each texture attribute
表1. 各纹理属性的原理及含义表
纹理是表现图像或物体表面的起伏程度和粗糙程度的一种属性，提取地震纹理属性的方法有统计分析法、模型分析法、结构分析法、信号处理法等，本次主要研究统计方法中的基于灰度共生矩阵的纹理分析方法 [8] 。
利用灰度共生矩阵进行纹理分析的方法最早由Haralick提出，之后Love等利用纹理属性分析构造引起的较大尺度的纹理特征，Gao和West等将灰度共生矩阵提取纹理属性的方法引入三维地震数据，从多角度分析地质体并识别沉积环境 [9] [10] 。现如今提取出的纹理属性已可用于裂缝检测、构造解释、河道预测、地震相及沉积相分析、储层预测等多个方面，因其受地震资料品质的影响较小而获得了较好的应用效果。
单纯的灰度共生矩阵不能直接用于纹理描述，而是需要将其进行属性计算，转化为具有地质含义的属性体进行解释。在Geoeast中可以利用灰度共生矩阵提取14种纹理属性，其中熵，均匀度，能量，对比度，相关性和方差是其中常用的六种纹理属性 [11] [12] 。以下介绍各种纹理属性的含义(表1)。
从结果看，其中纹理方差属性对于三角洲平原–内前缘与三角洲外前缘的相带边界刻画较为模糊，平原–内前缘相内属性值相对复杂，干扰信息较多，MH28周围内前缘相与外前缘相边界较为清晰(图7(a))。纹理均匀度属性平原–内前缘相内干扰信息较多，平原–内前缘与外前缘相带边界混杂，MH26扇体左侧平原–内前缘与外前缘难以区分，内前缘与外前缘相边界较为模糊(图7(b))。纹理能量属性内前缘与外前缘边界较为清晰，内部干扰信息较少(图7(c))。纹理熵和相关度属性图像较为模糊，各项带边界干扰信息较多，MH7附近扇体与JW26附近扇体边界模糊(图7(d)~(e))。纹理对比度属性冲积扇与平原相边界清晰，且平原相内部干扰信息少，可以清晰识别内前缘相和外前缘相边界(图7(f))。
综合以上纹理属性，认为纹理方差属性、纹理对比度属性、纹理能量属性比较适合该三维区块风城组地震相边界识别工作。
5. 结论
1) 研究区扇三角洲平原–内前缘主要为中弱振幅中低频较断续地震相和中振幅低频弱连续地震相，扇三角洲外前缘主要为中低频中强振幅较连续地震相和中强振幅连续地震相，滨浅湖主要为中强振幅低频较连续地震相，各沉积相带对应的地震相具有不同的特征，可建立较为准确的对应关系模板。
2) 在多种纹理属性对比分析的方法下能够清晰识别各地震相的边界，纹理差异性属性和纹理熵属性对于玛湖一区沉积相划分具有较好的指导意义。