0.   引言

地震沉积学是建立在地球物理学、地震地层学、层序地层学及沉积学等多学科基础之上的一门新兴学科，由曾洪流最先提出，之后得益于近年来国内三维地震地球物理解释技术的进步，受到了普遍关注^[1~4]。尤其在井网密度较低地区，地震资料的充分利用与否具有决定性的作用。地震沉积学研究的关键技术手段主要包括90°相位转换技术、地层切片技术和分频解释技术^[1]，其中90°相位转换技术主要使用地震振幅作为岩性的指导。通常使用的标准地震数据大部分都是零相位地震数据，主要因为其子波对称、最大振幅的中心与反射界面一致、分辨率较高。但是零相位地震数据的这种优势主要体现在单层反射界面，如重大的不整合以及厚的板状砂岩体。而对于厚层泥岩夹薄层砂岩的楔形地层，90°相位转换的地震剖面可将反射主瓣提到砂岩层中心，提高分辨率，减少噪音，使反射轴与测井曲线更加吻合，使地震反射具有了地层岩性意义^[5~8]。90°相位转换技术自提出以来便引起很多学者的争论^{[2, 9~11]}，主要由于90°相位转换技术使原本数据体意义模糊，不具有严格的等时性。但是，其操作简单、运算速度快、表现力强，适合海上少井条件的研究。本文旨在利用其表现力强的特点预测砂体展布并通过正演的方法，证明90°相位转换后剖面的地质意义。

砂体展布预测是储层及油气藏勘探的基础和关键步骤之一，而储集砂体精细刻画及分布形态模型的建立一直是研究难点，前人曾用多种手段，如测井、切片技术、反演技术、频谱成像、地震属性、地震相等，识别砂体的展布形态及分布特征^[12~19]。这些技术方法在不同地区获得了显著的成效，但是不同地区限制条件不同导致研究方法的选择也因地而异。本文以黄河口凹陷西洼沙河街二段作为研究区，该区的特点是地震资料品质较好、分辨率高、但钻井少。因此，研究将重点关注以下2个方面：① 90°相位转换地震剖面，分析砂泥岩的纵向展布关系，建立地震地质模型并探寻振幅与岩性的对应关系；② 提取相关属性，讨论沉积体系内部砂组的沉积相带平面分布特征及内部展布特征。

1.   研究区地质概况

研究主要以黄河口凹陷西洼古近系沙河街组二段为例。黄河口凹陷古近纪至现今构造面貌表现为西深东浅、北陡南缓、凹中有隆，整体为北断南超的箕状凹陷，其中郯庐断层西支穿过凹陷中部，形成中央隆起带将黄河口凹陷分为东、西2个次洼^[20~21]。近几年的勘探情况表明，黄河口凹陷是一个富生烃凹陷，主要发育沙三段、沙一段—沙二段和东三段3套主力烃源岩。黄河口凹陷构造上具有典型的断坳叠置特征，其中古近系主要以裂陷作用为主，形成很多半地堑，并可划分为4个裂陷伸展期^[22]，而黄河口凹陷西洼区域的沙河街组二段恰位于裂陷伸展的早期，地形经历了相对短暂的区域隆升。据西洼钻井分层的统计，沙二段整体厚度较薄，最多不超过200 m，以南部垦东—青坨子凸起、莱北低凸起的物源供应为主，为发育大面积砂体提供了充足的背景条件。对于研究区沉积体系内部砂体的展布特征，本文将通过90°相位转换构建地质模型，通过正演模拟技术确定砂体的地质-地震响应特征，并通过特征属性的提取预测砂体的展布范围。

2.   90°相位转换

首先对工区的地震资料进行90°相位转换。转换后的地震资料波形清晰，同相轴特征明显，地震剖面与岩性对应更直观（见图 1）。尽管测井分辨率远高于地震资料的分辨率，但从零相位地震资料与90°相位地震资料的对比可以发现，90°相位地震资料的波谷对应于伽马射线的低值段，而零相位地震资料的波谷对应于伽马射线的高值与低值的拼合段。

图  1  地震剖面与测井对比图

Figure  1.  The comparison and relation between seismic data and well data

下载: 全尺寸图片 幻灯片

2.1   岩性特征统计分析

通常，对于地震振幅的认识和使用主要依赖于不同岩性的波阻抗差异。尽管波阻抗值就是速度和密度的乘积，但是实际应用中波阻抗的影响因素众多，如岩体的深度、温度、压实度，以及岩性、物性甚至于沉积环境及储集相带的分布等。前人通过研究认为，随着地层埋深的增加及地质年代的推移，砂岩和泥岩波阻抗值均增大，在浅层通常泥岩的波阻抗值大于砂岩，但是深层随着压实作用的增强砂岩的波阻抗值大于泥岩^[23~24]。

通过邻近工区及本工区测井数据和资料统计分析得出，研究区目的层深度超过了3000 m，厚度较薄，为150~200 m，层段内岩性波阻抗值增幅受深度影响不明显。研究区目的层内主要发育次生孔隙，砂岩孔隙度-声波阻抗、孔隙度-速度及孔隙度-渗透率的交会分析显示，孔隙度、渗透率等属性特征对本研究层段速度及波阻抗影响不明显，而不同岩性波阻抗和速度差异明显，其中砂岩的波阻抗值大于泥岩（见图 2）。砂泥岩波阻抗可分且差异较大为本区依据振幅和岩性特征建立地质模型提供了理论依据。

图  2  砂泥岩波阻抗与自然伽马交会图

Figure  2.  Cross plot of sandstone and mudstone between impedance and natural Gamma

下载: 全尺寸图片 幻灯片

以研究区及邻近工区多口钻井为基础对岩性发育特征进行详细分析，以便于建立地质剖面模型和正演模拟。统计的主要岩石类型有砂岩、粉砂岩和泥岩。粉砂岩形成的岩层组段相对不多，厚度以薄层为主；砂岩层单层最大厚度超过30 m，有40%的砂岩单层组段厚度超过了4 m；泥岩段厚度以小于2 m的居多，典型的砂泥岩组合分布模式为复合厚互层岩性组合（见图 3）。取其中9口钻井资料进行统计分析，沙河街组二段砂岩的波阻抗值最大，粉砂岩次之，泥岩波阻抗值最小（见表 1）。砂岩速度主要分布在3900~4200 m/s，平均速度为4108 m/s；粉砂岩速度主要分布在3700~4000 m/s，平均速度3900 m/s；泥岩速度主要分布在3500~3900 m/s，平均速度为3767 m/s。

图  3  砂泥岩层单层厚度统计及单井砂泥岩组合分布

Figure  3.  Statistics of sandstone and mudstone thickness distribution frequency and the lithologic distribution characters of a drilling well

下载: 全尺寸图片 幻灯片

表  1  钻井不同岩性速度、密度、波阻抗及单层厚度统计

Table  1.  Statistics of velocity, density, impedance and single thickness of different lithologics in a drilling well

下载: 导出CSV 

| 显示表格

2.2   砂体展布模型的建立及正演模拟

当地震剖面同相轴出现反射变弱、同相轴分叉、产生复波等特征时，可从沉积特征和地质概念模式分析^[25]出发，勾画剖面的地质模型。本文以90°相位转换后的剖面为基准，结合上述对研究区块内砂泥岩的统计分析及岩性组合特征，勾画砂岩层（见图 4）。图 4a为90°相位转换后的地震剖面，红色至蓝色的颜色渐变表述了随振幅能量逐渐变大的过程；图 4b为据地震剖面的同相轴及振幅特征描绘勾画出的砂体剖面展布模型，即剖面地质模型。以砂体剖面展布模型为基准进行正演模拟。

图  4  基于90°相位转换地震剖面建立地质模型

Figure  4.  The geological model based on 90° phase of seismic data

下载: 全尺寸图片 幻灯片

本文采用GMAPlus进行模型正演模拟研究，该技术考虑了介质各向异性，采用速度和密度建立地下波阻抗场，实现纵横向速度的梯度变化，构建地质剖面模型；根据垂直入射原理，自激自收，模拟出地震反射波，得到合成记录道。正演模型砂泥岩速度值的设定主要依据对工区20余口钻井资料的砂泥岩速度及密度的统计和岩性分布规律。正演剖面的分辨率与地质模型及入射子波的频率有关^[26~28]。大部分对于正演模拟技术的应用都是通过提高地震主频来获得高分辨率信息，但是实际的地震资料很难达到理想的高频，因此本文选用与地震资料相近的频率，即20 Hz。

正演模拟是校准地质模型的重要手段，对砂体展布预测具有重要意义。本文对得到的正演合成记录道分别进行零相位、45°相位和90°相位的转换（见图 5）。对比3种相位转换的合成记录道可以看出，同相轴主要沿着砂岩层轮廓展布。其中零相位合成记录剖面中，砂岩主要分布于振幅的波峰和波谷之间（见图 5a）；45°相位合成记录剖面中，砂岩层略向波谷一侧集中，但砂岩轮廓顶面又超越零值界面偏向波峰（见图 5b）；90°相位合成记录剖面，砂岩层整体位于波谷，且随波谷的趋势展布，形状直观明显（见图 5c），再次验证了90°相位转换后砂体对应剖面波谷的特征，即90°相位转换剖面的直观性。图 6对比90°相位合成记录剖面与原始正演模型生成的地震剖面（绿色代表砂岩层的轮廓），发现其相似性较高，从而确立砂岩体剖面模型。说明以90°相位转换后地震剖面为基础直接勾画出的岩性模型剖面可信度较高。

图  5  正演合成记录道（黄色为砂体展布层段）

Figure  5.  The seismic modeling synthesis sections

下载: 全尺寸图片 幻灯片

图  6  原始地震剖面与正演合成地震剖面对比

Figure  6.  Comparison between primitive seismic section and seismic modeling synthesis section of 90° phasing

下载: 全尺寸图片 幻灯片

综上所述，90°相位转换后的地震资料可以更直观地勾画出砂岩层的地质模型，然后通过与原始地震资料的对比，进行正演模型的反复迭代，即可得到更为真实的砂泥岩纵向分布模式，从而预测砂体的展布。

3.   砂体展布分析

地层切片技术是地震沉积学的关键技术手段之一，在层切片的基础上可以对地震数据进行不同种类的属性提取，而以地震属性特征进行沉积相、含油气性、砂体分布等的预测是储层研究的必要基础手段。通常情况下，实际地震剖面的主频较低，单层砂岩很难达到调谐厚度。而同一区域砂体速度比较稳定，因此可以根据地震属性研究小于调谐厚度的砂层厚度变化及边界展布范围^[29]。随着不同地震数据提取方法和算法的产生，通过实践论证的地震属性已经由最初的几种增加到几百种，而这些地震属性对岩性波阻抗的反应不一，很多都不具有物理意义或地质意义，因此必须结合研究内容和地质背景进行筛选和优化，利用最能反应研究区砂体展布状况的地震属性切片进行分析。地震数据的振幅分量是判定岩石物理参数（声波阻抗、反射参数、速度、吸收系数等）的主要因素^[30]。而通过以上对砂泥岩的速度分析以及90°相位转换后的地质模型的建立及正演模拟，发现90°相位转换后的地震剖面的波谷和负振幅值属性与砂岩对应关系较好。因此，对研究区90°相位转换后的地震数据进行负振幅和（见图 7）、平均负谷值振幅（见图 8）等属性的提取，依据振幅值的变化可以清楚地识别三角洲沉积体系中砂体的展布特征。

图  7  属性切片与地震剖面相结合预测砂体展布

Figure  7.  Combine the attribution slice and the seismic section to predict the sandbody's distribution

下载: 全尺寸图片 幻灯片

图  8  沙河街组二段平均负谷值属性及砂体展布平面图

Figure  8.  ANT amplitude and sandbody distribution of Es2

下载: 全尺寸图片 幻灯片

纵向剖面与平面属性切片相结合可建立清晰的砂层结构特征，并能精细地刻画出砂岩层的平面分布样式，如平面上呈多块砂体的拼接，说明砂岩层内部并非都是均质的，剖面上可能是多期砂体叠置的复合体，以及受同沉积断裂的影响，砂体的展布和延伸规律均得以进一步的探讨和分析。此外，依据砂岩层的地球物理相应特征，可以圈定出砂岩体的展布边界。

4.   讨论

通过提取地震属性进行砂体展布范围预测是一种非常普遍的方法。本文不同于以往方法对地震剖面直接进行属性提取，而是首先对地震剖面进行90°相位转换，然后分析转换后的地震剖面与岩性之间的对应关系，勾画砂岩层纵向展布特征即建立地质模型。正演模拟证实，90°相位转换后的地震剖面与砂岩体显现较好的匹配关系，通过正演合成记录道与地震剖面的对比，可确立纵向剖面砂体展布模型。

采用正演模拟技术确立地震剖面的地质模型也有一定的缺点和不足，主要表现在其结果不具有唯一性，厚层泥岩夹单层厚砂岩与夹多层砂泥薄互层对应的振幅值相差不大。但是，其总体的岩性特征是可确定的，结合单层厚度岩性统计分析及地质背景，再通过反复正演模拟调整地质模型，可使其更接近于真实地质情况。砂岩主要对应于90°相位转换后地震剖面的波谷，且负振幅绝对值的增大对应于砂泥岩反射系数的增加。据此特征对90°相位转换的地震剖面进行相应的属性提取，其中属性切片振幅值的波动标志着砂岩层的展布及叠置形态。

关于90°相位转换技术的可行性，笔者认为不同地区的沉积环境及岩性物性条件不同，不可一概而论。而本文研究区地层平缓，目的层位于深层且厚度较薄，岩性区分度较好，是90°相位转换及其相关技术的应用得以实现的优势。

5.   结论

通过钻井和测井资料统计分析，发现研究区砂泥岩速度区分度较好，受岩石物性影响较小，且90°相位地震资料岩性特征显示明显。据此建立的地质模型正演得到砂岩恰好对应于合成记录道的波谷，且波谷的波形能直观反映出砂体的轮廓，证实了90°相位转换地震剖面的岩性意义。

研究90°相位转换后的地震剖面的属性特征，并进行相应的属性提取可展示出砂体的分布范围及轮廓，受断层控制砂岩层更倾向于沿断层方向展布。