地震反演是储层预测的核心技术之一, 该技术基于地震原始数据, 综合运用地质、测井等多学科的相关信息来获得地下各层位波阻抗的分布, 是一门综合性的地球物理勘探技术^[1]。测井声波时差曲线是反映储层特征最重要也是最常规的资料, 在实际地震反演中利用声波时差曲线开展储层波阻抗反演已取得了非常显著的地质效果^[2], 应用实钻井测井数据, 在沉积背景约束下建立地质模型, 实现井间地震资料的约束反演, 以此补充地震资料带宽的不足, 从而实现优势互补^[3]。波阻抗反演是地震储层反演预测研究中的核心技术, 因此地震反演通常也被称为波阻抗反演^[4], 波阻抗反演是高分辨率地震资料处理的最终表达形式, 在圈闭评价、油气藏描述及井位论证中均起到了至关重要的作用^[5]。

苏里格气田苏X区块是典型的岩性气藏, 属三角洲河流相沉积, 其砂体空间变化规律复杂。该区块经过多年开发, 天然气富集区井网基本完善, 其外围区盒8段辫状河道沉积储层受成藏条件及非均质性影响, 气水分布规律复杂, 含水层较多; 山西组曲流河道砂体宽度窄, 厚度薄, 横向变化快, 连片性差, 砂体地震识别及平面展布规律预测难。研究区储层以砂泥岩薄互储层为主, 相同沉积单元内的薄层单砂体甚至在百米范围内会呈现横向突变的空间几何特征, 储层非均质性强, 且测井声波阻抗曲线对该类砂体响应不明显, 即使在气田深度开发阶段的密集井网条件下, 仅依靠井点信息或地震信息预测河道砂体难度亦较大^[6], 在外围区砂体精细识别与刻画的难度更大。因此, 利用气田开发区综合井网密集和地震资料横向分辨率高的优势, 探索密集井网条件下井震联合, 应用地震资料拓频技术、贝叶斯反演技术及地震数据体分频特征重构技术3种技术手段相结合的高分辨率地震多子波反演方法, 开展储层预测工作, 精细刻画砂体轮廓及内幕、砂体叠置关系, 以此夯实老气田挖潜的重要基础工作。

苏里格气田位于鄂尔多斯盆地中北部, 处于石嘴山和苏里格两个三角洲沉积体系交界处, 区域构造处于鄂尔多斯盆地陕北斜坡西北部, 属内蒙古自治区鄂尔多斯市的鄂托克旗和鄂托克前旗所辖。该气田目的层以石盒子组盒8段和山西组山1段砂体为主, 主要目的层盒8段底界为一宽缓的西倾单斜, 坡降为3~10 m/km, 单斜背景上发育有多排北东-南西走向的低缓鼻隆构造, 幅度在10~20 m之间。通过实际钻探研究证实, 该区发育的低缓鼻隆构造对天然气聚集控制作用较弱, 天然气藏展布主要受砂体和物性控制, 属于典型致密砂岩岩性类型气藏^[7]。苏X区块位于苏里格气田的中西部, 与苏里格气田总体构造特征基本一致, 表现为由北东向南西倾斜的单斜构造, 局部发育一些低缓的鼻隆构造, 幅度在8~15 m之间。研究区内完成二维地震888.3 km/48条, 测网密度相对均匀, 达到1 km× 2 km, 但地震资料受分辨率、信噪比影响, 品质较差, 无法满足精细预测储层展布的需求。

地震反演效果的好坏主要取决于两个条件：一是地震资料品质的好坏, 包括信噪比、分辨率及保真度, 这决定了最终反演结果的分辨率和可信度; 二是测井资料的岩电组合特征明晰与否, 包括测井声波时差、密度及纵波阻抗曲线能否很好地反映储层和围岩的差异, 这决定了波阻抗反演结果能否有效地区分储层与非储层^[8]。在进行地震反演工作前, 首先要对这两方面的工作进行详细分析和处理, 做好高分辨率地震反演的基础工作, 以确保最终反演结果精准有效。

2.1 地震数据拓频处理

苏X区块叠后地震资料信噪比低, 目的层同相轴连续性差, 目的层段的地震资料主频在15 Hz左右, 有效频带宽度10~30 Hz。理论上, 地震垂向分辨率为调谐厚度, 即H=λ /4, 其中λ =V/F^[9]（H为地层厚度, λ 为地震子波长度; V为目的层平均速度; F为目的层地震资料主频）。研究区目的层段的地震层速度约为3 000 m/s, 因此, 地震资料垂向分辨率大约为50 m, 只能识别大于50 m的砂层组, 远不能满足实际生产需要。针对研究区地震资料主频低、频带较窄的客观情况, 首先采用反Q滤波补偿技术对地震数据进行振幅补偿和子波相位一致性校正; 其次, 对反Q滤波剖面开展宽频成像分辨率提高技术处理, 适当提高地震资料主频, 拓展地震数据高频和低频, 增加地震有效频带宽度, 从而提高地震资料品质。

2.1.1 反Q滤波补偿技术

因为地震波在传播过程中存在被吸收后逐渐衰减的情况, 所以地震子波具有时变和空变的特性, 不同程度地制约了提频效果。而叠后提频主要针对子波频率进行处理, 经过大量的试验验证, 在经Q滤波校正了子波一致性后的数据体上提频效果更佳。

反Q滤波补偿是在求得的较为精确的品质因子Q场的基础上, 补偿由于地层滤波作用所造成的振幅和频率的衰减。这里采用基于Gabor变换的反Q滤波补偿方法对地震记录进行Gabor变换, 再结合振幅补偿算子及相位补偿算子在时频域进行反Q滤波补偿, 对所得的波场利用逆Gabor变换可得到经由反Q滤波处理的地震记录^[10], 并且针对不同信噪比的实际数据, 可以通过输入补偿因子来控制噪音, 在提高分辨率的同时, 兼顾地震资料的信噪比。反Q滤波补偿对地震分辨率的提升极为有限, 其主要功能是补偿振幅和子波相位一致性, 从而为后续提高分辨率处理奠定基础。

2.1.2 宽频成像分辨率提高技术

宽频成像分辨率提高技术是通过对原始地震数据应用连续小波变换（CWT）时间序列分析, 利用连续小波域的多分辨率特征, 按照谐频分析的思路, 可以计算出地震子波的谐波信息（基波频率的整数倍）。将谐波信息加到原始地震数据中, 可以有效拓展地震数据的频宽, 进而更好地刻画地下的地层反射, 其优点在于：（1）从原理上突破了传统褶积模型的限制, 不再需要估计地震子波; （2）利用优势频段提取谐波分量, 保证了拓频结果具有较高的信噪比; （3）谐波成分不改变原始信号的相位, 有效保持了地震数据的相对振幅关系; （4）不仅高频端拓频能力突出, 还能够拓展低频端。

经过宽频成像分辨率提高处理后, 地震数据主频提高了10 Hz, 有效频带宽度拓宽了15 Hz, 分辨率获得提高的同时, 前后地震相位特征、波组关系、信噪比、振幅横向变化都得到了保留（图1）。盒8段和山西组分辨率提高后, 出现与砂岩对应的稳定地震轴, 地震资料对薄砂层的响应特征更为明显, 为地震反演储层预测奠定了资料基础。

图1

	Figure Option ViewDownloadNew Window
	图1 宽频成像分辨率提高前后剖面对比

2.2 测井岩石物理分析

测井岩石物理分析是地震反演储层预测的基础, 通过分析岩石弹性参数与储层物性参数之间的关系, 找到能够较好反映岩性或储层物性的敏感参数^[11]。储层敏感参数的选取决定了反演方案的制定及反演结果的可靠性。

2.2.1 盒8段测井岩石物理分析

盒8段地层总体表现为砂岩高波阻抗、泥岩低波阻抗特征, 但砂岩含气或含水后, 其波阻抗降低, 与泥岩无法区分。苏X区块只有叠后二维数据, 因此只能开展叠后地震反演储层预测, 这需要寻找敏感的储层参数。测井自然伽马参数在本区可以识别岩性, 选择自然伽马曲线的高频部分和波阻抗曲线的低频部分进行重构, 获得重构波阻抗, 并选择对子波不敏感的贝叶斯地震反演技术进行本区储层预测。

从盒8段原始波阻抗-重构波阻抗-自然伽马交会图板（图2）可见, 原始波阻抗无法区分盒8段的砂、泥岩, 但是重构波阻抗砂、泥岩界限清晰, 因此可以用自然伽马拟声波重构波阻抗进行苏X区块二维工区多子波反演预测本区储层。

图2

	Figure Option ViewDownloadNew Window
	图2 盒8段原始波阻抗-重构波阻抗-自然伽马交会图板

2.2.2 山西组测井岩石物理分析

山西组地层岩性主要由浅灰色含砾中砂岩、细砂岩、灰色泥岩和煤层组成, 其中砂岩表现为相对低自然伽马、高波阻抗特征; 泥岩表现为相对高自然伽马、高波阻抗特征; 煤层则表现为低自然伽马、低波阻抗特征。由此可知, 单一参数无法同时区分3种岩性, 通过自然伽马拟声波重构波阻抗发现, 重构波阻抗后砂岩表现为高波阻抗, 而泥岩和煤层表现为低波阻抗, 因此自然伽马拟声波重构波阻抗可以用来进行山西组的地震反演预测。

为了验证重构波阻抗识别岩性的可靠性, 对工区内的34口井进行了应用分析对比, 煤层和泥岩在重构波阻抗上均表现为低波阻抗, 砂岩表现为高波阻抗。从山西组原始波阻抗-重构波阻抗-自然伽马交会图板（图3）可以看出, 重构后的波阻抗砂岩和煤层、泥岩的门槛为10 600 g/cm³· m/s, 证实了重构波阻抗识别岩性的可靠性。

图3

	Figure Option ViewDownloadNew Window
	图3 山西组原始波阻抗-重构波阻抗-自然伽马交会图板

通过分析岩石弹性参数与储层物性参数间的关系, 找到了能够较好反映岩性的敏感参数, 采用声波曲线重构的方法, 进行自然伽马拟声波重构波阻抗反演, 较好地提高了岩性区分识别能力, 从而解决了复杂储层描述存在的问题。

研究区目的层石盒子组与山西组属于典型三角洲河流相沉积, 储层为典型的砂泥岩互层, 砂体横向变化快, 单砂体厚度差异大, 厚度范围2~20 m。二维地震常规反演方法的分辨率与地震分辨率较为接近, 难以识别5 m左右的储层。多子波反演是在稀疏脉冲反演方法的基础上演化而来, 以贝叶斯反演算法为依据, 区别于常规稀疏脉冲反演单一界面假设, 多子波反演是以薄层顶底一对反射系数为假设, 考虑了薄层调谐因素, 反演结果分辨率更高, 并且适用于二维地震工区。

3.2 分频特征重构

将多个分频地震反演结果进行分析, 优选出能反映特定地质目标的分频地震反演数据体, 然后在频率域进行重构融合, 最终形成一个具有一定带宽和主频的高分辨率地震反演数据体, 使得最终反演结果呈现更高的纵向分辨率, 对薄储层有更突出的表征。该方法不需要任何先验模型、反射系数的分布假设, 能有效提高地震资料的分辨率, 更有利于薄层岩性的识别和厚度估算。

本文利用分频操作所得15、25、35 Hz分频地震数据体（分别表征低频特征、中频特征及高频特征数据）与钻井现场模型数据结合进行贝叶斯反演, 得到对应的地震反演数据体。在此基础上, 对3个地震反演数据体进行频率域的融合, 其权重系数分别是0.3、0.3、0.4, 最终得到高分辨率多子波反演数据体（图4）。

图4

	Figure Option ViewDownloadNew Window
	图4 苏X区块地震反演剖面

多子波地震反演作为一种高分辨率的二维叠后地震反演储层预测方法, 其优点在于受地震资料品质的影响相对较小, 可突破地震频带宽度的限制, 能获得相对较高分辨率的地层波阻抗等属性资料。

由高分辨率多子波地震反演与常规稀疏脉冲地震反演剖面对比（图5）可知, 对比井段纵向上的砂体分布特征与井资料吻合性较好, 表明自然伽马拟声波多子波地震反演更能预测出砂体的具体发育位置。从图5中可以看出, 相比常规稀疏脉冲地震反演, 多子波地震反演纵向分辨率大幅提高, 同时横向分辨率也具有较大的提升, 井震吻合度更高, 可以区分出砂体的分布, 总体砂体横向变化更自然, 尖灭点更清晰, 无论厚层还是薄层, 井间变化都很合理, 高分辨率多子波地震反演效果明显优于常规稀疏脉冲反演。

图5

	Figure Option ViewDownloadNew Window
	图5 高分辨率多子波地震反演与常规稀疏脉冲地震反演剖面对比

应用高分辨率多子波地震反演技术对苏X区块主要目的层盒8下亚段和山1段进行了储层分析。图6为盒8下亚段多子波地震反演储层平面图和砂岩厚度统计, 可见砂体主要呈近北西-南东向条带状展布, 研究区中部发育一个南北向的砂岩条带; 储层整体为中部厚, 东西两侧薄, 储层厚度一般为10~40 m, 最厚处在研究区南部, 厚度约40 m。图7为山1段多子波地震反演储层平面图和砂岩厚度统计, 显示储层主要呈近北西-南东向的条带状展布, 研究区中部发育一个北西-南东向砂岩条带, 储层整体为中部厚, 东西两侧薄, 厚度一般为10~30 m, 最厚处在研究区东南部, 厚度约30 m。根据反演结果提取目的层砂岩厚度, 结合验证井和新钻井资料来验证多子波地震反演结果的预测精度, 可以达到最佳程度的吻合。

图6

	Figure Option ViewDownloadNew Window
	图6 盒8下亚段多子波地震反演储层平面图和砂岩厚度统计

图7

	Figure Option ViewDownloadNew Window
	图7 山1段多子波地震反演储层平面图和砂岩厚度统计

针对二维地震资料采用反Q滤波补偿和宽频成像分辨率提高技术等地震资料处理方法, 可有效拓宽原始地震数据的频带宽度, 提高主频, 提升地震资料品质, 这是精细刻画储层特征的基础。精细进行测井岩石物理分析, 明确储层地震波阻抗响应特征, 选取对储层敏感的测井曲线, 进行拟声波重构波阻抗, 这是做好储层地震反演的关键。对提频后的地震数据体进行分频操作, 将生成的各分频地震数据体与测井模型数据结合, 分别进行贝叶斯地震反演, 反演后的数据体通过分频特征重构技术进行频率域重构融合, 最终得到高分辨率多子波地震反演数据体, 即高分辨率多子波地震反演技术。采用高分辨率多子波地震反演对苏X区块进行储层预测, 取得良好效果, 相对二维常规稀疏脉冲地震反演, 高分辨率多子波地震反演在纵横向分辨率上提升明显, 砂体的叠置关系及井间砂体的尖灭展现更清晰, 反演结果更可靠, 更符合地质规律, 能够为有利层位的优选和评价提供依据。

（编辑 唐艳军）