基于小波边缘分析与井—震联合建模的波阻抗反演技术在陆梁隆起带储层预测中的应用

doi:10.11720/wtyht.2023.1580

基于小波边缘分析与井—震联合建模的波阻抗反演技术在陆梁隆起带储层预测中的应用

史全党¹, 孔令业¹, 吴超^,¹, 丁艳雪¹, 刘泽民¹, 于雪¹, 王江^,²

1.中国石油新疆油田公司采气一厂,新疆克拉玛依 834000

2.中国石油大庆油田有限责任公司勘探开发研究院,黑龙江大庆 163712

Application of wave impedance inversion technology based on wavelet edge analysis and combined well-seismic modeling in reservoir prediction of Luliang uplift zone

SHI Quan-Dang¹, KONG Ling-Ye¹, WU Chao^,¹, DING Yan-Xue¹, LIU Ze-Min¹, YU Xue¹, WANG Jiang^,²

1. PetroChina Xinjiang Oilfield Company Gas Production Plant No.1,Karamay 834000,China

2. Exploration and Development Research Institute of Daqing Oilfield Co.,Ltd.,Daqing 163712,China

通讯作者: 王江(1966-),男,博士,高级工程师,主要从事地震勘探数据处理、解释方法技术研究工作。Email:wangjiaang@petrochina.com.cn

责任编辑: 叶佩

收稿日期: 2022-12-16 修回日期: 2023-09-8

基金资助:

国家重大科技专项“大型油气田及煤层气开发”(2011ZX05001-001-004)

Received: 2022-12-16 Revised: 2023-09-8

作者简介 About authors

吴超(1991-),男,高级工程师,主要从事油田勘探及油藏评价工作。Email:wuchao12@petrochina.com.cn

摘要

针对准噶尔盆地陆梁隆起断裂构造复杂、储层横向非均质性强、含气砂体小而薄预测难的问题,应用小波边缘分析与井—震联合建模的波阻抗反演技术,从地震记录中直接提取反映构造和岩性变化的地震属性特征参数,与测井声波阻抗数据一起建立初始模型并参与波阻抗模型的扰动修改。此反演技术弥补了井间插值建模过程中井间高频成分的缺失和井间岩性的局部变化细节,避免了常规波阻抗反演过程中由于初始模型不准确而产生的错误信息,提高了地震资料识别“小而薄”砂体的分辨能力。结果表明:受东部克拉美丽山物源控制,在陆梁隆起带DX14井区梧桐沟组形成扇三角洲—半深湖(深湖)沉积体系;发育一系列扇三角洲前缘砂体。实际钻井结果与钻前预测对比表明,井点处预测砂岩厚度的绝对误差小于0.60 m,相对均方误差小于2.84%,预测精度满足精细储层预测要求。研究成果可为精细油藏描述和井位部署提供依据。

关键词： 小波边缘分析; 井—震联合建模; 波阻抗反演; 储层预测; 地震属性特征参数; DX14井区; 陆梁隆起带

Abstract

The Luliang uplift zone in the Junggar Basin exhibits intricate fault structures,laterally heterogeneous reservoirs,and small-scale and thin gas-bearing sand bodies,presenting challenges in reservoir prediction.Hence,this study applied the wave impedance inversion technology based on wavelet edge analysis and well-seismic joint modeling to directly extract seismic attributes' characteristic parameters that reflect structural and lithological changes from seismic records.These seismic attribute characteristic parameters were used to build the initial model together with acoustic impedance logs and participated in disturbance modification of the wave impedance model.This inversion technology counteracted the lack of inter-well high-frequency components and the inter-well local lithologic changes during the inter-well interpolation modeling and avoided error information caused by the inaccurate initial model in the conventional wave impedance inversion,thus improving the resolution of seismic data in identifying small-scale and thin sand bodies.The results show that under the control of the provenance of the Kelameili Mountain in the east,a fan delta-semideep (deep) lacustrine sedimentary system formed in the Wutonggou Formation in the DX14 well area of the Luliang uplift zone,hosting many fan delta-front sand bodies.The comparison between the actual drilling results and the pre-drilling prediction results indicates that the absolute error and relative mean error of sandstone thickness prediction at the well site were less than 0.60 m and less than 2.84%,respectively,suggesting that the prediction accuracy meets the requirements of fine-scale reservoir prediction.The research results provide a guide for fine-scale reservoir description and well deployment.

Keywords： wavelet edge analysis; well-seismic joint modeling; impedance inversion; reservoir prediction; seismic attribute characteristic parameter; DX14 well area; Luliang uplift zone

PDF (3342KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

史全党, 孔令业, 吴超, 丁艳雪, 刘泽民, 于雪, 王江. 基于小波边缘分析与井—震联合建模的波阻抗反演技术在陆梁隆起带储层预测中的应用[J]. 物探与化探, 2023, 47(6): 1425-1432 doi:10.11720/wtyht.2023.1580

SHI Quan-Dang, KONG Ling-Ye, WU Chao, DING Yan-Xue, LIU Ze-Min, YU Xue, WANG Jiang. Application of wave impedance inversion technology based on wavelet edge analysis and combined well-seismic modeling in reservoir prediction of Luliang uplift zone[J]. Geophysical and Geochemical Exploration, 2023, 47(6): 1425-1432 doi:10.11720/wtyht.2023.1580

0 引言

随着含油气盆地精细构造解释、精细储层预测和精细油藏描述工作的深入,对地震资料分辨率的需求越来越高^[1⇓⇓-4],精细储层预测的手段也一直以地震波阻抗反演为主^[5⇓⇓-8]。根据目标区的构造背景、沉积环境制定合理的技术对策和利用波阻抗反演方法来减小波阻抗反演结果的多解性、提高储层预测能力一直是地震波阻抗反演和储层精细预测的攻关方向^[9-10]。姚逢昌等^[11]通过对基于模型的波阻抗反演方法进行分析,认为多解性是基于测井约束性框架模型波阻抗反演方法所固有的特性,取决于初始模型和实际地下构造特点、沉积特征的符合程度。现有的测井约束波阻抗反演方法由于各自不同的方法原理,均存在一定局限性^[12-13],地震高频分量主要依赖于初始波阻抗模型从测井曲线中获取,初始建模主要是在解释层位控制下的井间插值,没有考虑断层在初始建模中的作用,横向上也没有充分利用地震信息。当目标区构造条件复杂、断层发育、沉积环境和岩性横向变化较大时,利用钻井数据井间插值或外推得到的初始模型低频背景较好,而反映岩性横向井间局部变化的高频成分获取较难,距离插值井点越远,产生的误差越大,减小了波阻抗反演方法对初始模型高频成分的应用程度,降低了波阻抗反演储层预测精度。

基于小波边缘分析与建模的波阻抗反演(acoustic impedance inversion based on wavelet edge analysis and modelling)技术(简称AIW),利用小波边缘分析方法从地震记录中直接提取地震属性特征参数,在断层参与的高密度层控构造框架约束下,同测井声波阻抗数据一起建立初始模型,相互约束、相互补充,充分利用地震数据横向分辨率高、纵向控制层位密度大的优势,避免了常规波阻抗反演过程中初始模型建立不准而产生的地质影响。许崇宝等^[14]、崔永福等^[15]、谢裕江等^[16]应用AIW波阻抗反演技术对准东煤田和新疆红沙泉区块侏罗系西山窑组含煤地层的煤层厚度、塔里木盆地英买力地区志留系沥青砂岩以及中国MOU气田盒8段河流相心滩砂岩进行了预测,弥补了常规波阻抗反演建模过程中井间高频成分缺失,提高了储层反演精度,但应用的区块大多构造简单、地层横向变化不大,很少针对构造复杂、断裂发育、沉积环境多变的断陷盆地以及低分辨率、低信噪比目标区进行研究。

本文应用小波边缘分析与井—震联合建模的波阻抗反演技术,考虑准噶尔盆地陆梁隆起带3D地震的波组反射特性以及地下地震波场的分布特点,通过小波边缘分析提取反映岩性局部变化的地震特征参数,在断层参与的高密度层控构造框架约束下,参与初始模型的建立和波阻抗的反演迭代,确定波阻抗反演的关键环节、优化克拉美丽气田DX14井区储层预测流程,不但精细刻画了克拉美丽气田DX14井区扇三角洲前缘砂体的分布,也提高了有利储层的预测精度。

1 区域地质概况

克拉美丽气田位于准噶尔盆地陆梁隆起东南部的滴南凸起中段,滴南凸起整体为一大型西倾的鼻状构造,受准噶尔盆地多期的构造建造和改造影响,断裂发育、构造复杂。滴南凸起中段是准噶尔盆地腹部最早展开勘探的地区之一,至今已发现了石炭系、二叠系、侏罗系和白垩系多个油气藏,其中二叠系梧桐沟组油气藏埋深3 500~4 000 m,而且构造破碎、断裂发育,多发育为断层控制的小断块和小断鼻构造,气藏类型主要为构造—岩性气藏。主要含油气层为梧桐沟组一段,岩性为细砂岩和粉砂岩,储集空间类型主要为剩余粒间孔及粒内溶孔。根据电性及储层特征,自下而上可划分为两段(P₃wt₂、P₃wt₁),油层主要发育在P₃wt₁,砂岩储层不但低孔、特低渗、非均质性强,而且储层横向连续性差,纵向变化大,使得滴南凸起中段DX14井区的地震资料不但构造成像效果差、分辨率低,而且砂体的地震响应特征不明显。

通过滴南凸起中段DX14井区原始地震资料分析,在目的层段梧桐沟组存在弱—空白以及低分辨率的地震反射,地震资料主频为30 Hz,如果地层平均速度按3 900 m/s计算,滴南凸起中段DX14井区3D地震资料可识别砂体的最小厚度为32.5 m,而滴西DX14井区含气砂体厚度一般为2.0~15 m,所以常规地震资料的分辨能力是无法满足精细构造解释和精细储层预测的地质需求,因此,应在精细井—震一体化初始模型建立的基础上开展高精度的波阻抗反演和精细储层预测。

2 技术方法

2.1 多尺度小波边缘分析与地震特征参数建模

小波边缘分析是一个理想的地震特征参数提取方法,利用小波边缘分析从地震数据本身提取地震特征参数得到地震数据局部的变化信息,通过地震数据的小波变换函数的一阶微分极大值和二阶微分零值进行处理,对处理后的梯度点进行判断,确定其为边缘点。

对地震数据 $d$ 进行小波变换得到 $W_{s} (d)$ ,计算公式为

(1)${W_s(d)=\frac1s\int_{-\infty}^\infty b\psi(\frac{t-x}s)\mathrm{d}x},$

对式(1)两侧求偏导,且令

(2)

\frac{\partial W_{s} (d)}{\partial t} = 0

(3)

\frac{\partial^{2} W_{s} (d)}{\partial t^{2}} = 0

式中:d为地震数据;W_s(d)为地震数据的小波变换;s为小波尺度;x为平移量,s;ψ为母小波函数; t 为时间,s;b为时间延迟x秒的地震数据,小波尺度s控制小波函数的伸缩,平移量x控制小波函数的平移。

联合式(1)、式(2)、式(3),求解此方程组,可得到地震记录中不同尺度、不同级别的地震特征参数边缘点(断层和岩性体的边缘)。

这样,通过小波变换把地震信号分解为多个不同尺度的分量,利用其良好的局部时频分辨力提取不同精度、不同奇异度的边缘点,最后连接起来就形成综合边缘。提取的地震特征参数的边缘点不但反映了原始剖面中不易识别的地质特征和地下岩性纵、横向的局部变化,将综合边缘点控制的地震特征参数同测井曲线一起参与波阻抗反演建模,克服了井间插值反演建模过程中无法充分利用地震信息的不足,获取了反映岩性横向井间局部变化的高频成分,从而也补充了波阻抗反演过程中高、低频地震分量。

2.2 AIW波阻抗反演算法

AIW波阻抗反演采用非线性全局优化的非常快速模拟退火法(VFSA)作为基本反演方法,该算法通过模型反复扰动修改不但可以避免线性反演算法强烈依赖初始模型而落入局部极值的弊端,且可以高精度求得反演问题的全局最优解^[17],使得AIW波阻抗反演结果与地下地质情况更加符合,其算法如下:

1)将归一化相似系数定义为AIW波阻抗反演的目标函数E^[18],来计算初始框架模型g₀的目标函数值E(g₀),计算公式为

(4)

E = 1 - \frac{\sum_{i = 1} (a_{i} - ā_{i}) (a'_{i} - ā'_{i})}{{[\sum_{i = 1} (a_{i} - ā_{i})^{2}]}^{1 / 2} {[\sum_{i = 1} (a'_{i} - ā'_{i})^{2}]}^{1 / 2}}

式中:E为目标函数值;a_i、ā_i为实际地震道的相对振幅和平均相对振幅,a'_i、ā'_i为合成地震道的相对振幅和平均相对振幅。

2)联合地震特征参数和测井数据一起来约束模型扰动中的模型变量的变化区间[A_i,B_i],修改初始模型g₀获取新模型g,并根据关系式(4)计算新模型g的目标函数值E(g)以及新模型目标函数E(g)与当前模型目标函数E(g₀)差ΔE,即ΔE = E(g)-E(g₀)。

3)当ΔE<0时,确定新模型g为最终的初始模型,如果确定了新模型g,置g₀=g,E(g₀)=E(g)。

4)当ΔE>0,则按依赖于温度的似Cauchy分布产生新模型,计算公式为:

(5)

g'_{i} = g_{i} + y_{i} (B_{i} - A_{i})

(6)

y_{i} = T s g n (u - 0.5) [{(1 + 1 / T)}^{|2 u - 1|} - 1]

式中: $g'_{i}$ 为修改后的第i模型, $g'_{i} \in [A_{i}, B_{i}]$ ;g_i为当前的第i模型,g_i∈[A_i,B_i];y_i为扰动因子;T为温度,K;sgn为符号函数;u为均匀分布的随机量,且u∈[0,1]。

5)同一温度下,多次重复步骤3)和4),缓慢降低温度,重复步骤3)~5),直至收敛条件ΔE<0,反复迭代反演,获取高精度的AIW波阻抗反演数据体。

2.3 AIW波阻抗反演的关键环节

AIW波阻抗反演的关键在于利用小波边缘分析方法从地震记录中直接提取地震属性特征参数,由于提取的地震属性特征参数不但反映地下岩性横向局部变化,而且直接参与了波阻抗反演过程中初始模型的建立和反复迭代反演,弥补了井间插值建模过程中井间高频成分的缺失和井间岩性的局部变化,避免了常规波阻抗反演过程中初始模型建立不准确而产生的影响。主要关键环节为:

1)利用测井和地震资料开展井震精细对比,确定不同波组的地质表征,通过纵向细划地层对比单元、横向精细对比,对三维地震资料开展精细的构造解释,得到精细的构造和沉积格架,联合测井曲线建立波阻抗低频模型;

2)通过小波边缘分析提取反映地下岩性横向局部变化的地震属性特征参数,结合精细标定的测井曲线以及应用解释层位、解释断层和测井曲线建立的波阻抗低频模型重新建立新的AIW波阻抗模型,进一步获取波阻抗高、低频成分;

3)提取的地震属性特征参数同钻井声波测井曲线中速度或声波阻抗数据一起进行模型扰动和修改,利用非线性全局优化的非常快速模拟退火法,反复迭代反演,获取高精度的AIW波阻抗反演数据体。

3 应用实例

3.1 井—震精细标定与储层敏感性分析

为了建立精细的波阻抗反演初始框架模型,以DX14井区梧桐沟组一段上、下分布的两套含气砂层为对比标志,对DX14井区梧桐沟组地层进行纵向细分、横向精细对比,划小纵向地层对比单元。根据地层旋回、地震层序特征,基于井—震精细标定技术进一步将梧桐沟组一段纵向细划了3个十米级地层层序旋回的砂组(P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ 、P₃w $t_{1}^{3}$ ),其上部的P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ 砂组是DX14井区梧桐沟组的主力含气层。从地层纵向细化、横向精细对比出发,以石炭系顶界不整合界面为地震资料时间—深度转换的标准层, 在地震剖面上精细标定各砂组位置, 以砂组为单元开展精细构造追踪、解释。梧桐沟组沉积时期,受东部克拉美丽山物源控制,在DX14井区发育一系列扇三角洲前缘砂体,扇三角洲前缘砂体具有高阻、低伽马、较高波阻抗的特征,与上覆和下伏地层具有较大的波阻抗差,表现为中强地震反射特征。

通过对DX14井区内所有钻井的测井资料进行曲线重采样、去野值、归一化和标准化预处理,对DX14井区主要含气层段(P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ )的测井曲线响应进行砂、泥岩识别的敏感性分析,发现砂岩表现为高电阻、低伽马、低声波和低密度的特征,利用波阻抗、密度和电阻率等电性曲线无法区分含气层段(P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ )的砂、泥岩(图1a)。为了精细刻画DX14井区砂岩分布,以声波曲线作为低频背景、自然伽马曲线作为高频成分重构拟声波曲线,来反演和预测梧桐沟组一段(P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ )砂体的展布情况。经伽马拟声波测井曲线重构,DX14井区梧桐沟组一段(P₃w $t_{1}^{1}$ 、P₃w $t_{1}^{2}$ )的砂岩均表现为高波阻抗特征,而泥岩则表现为低波阻抗特征,砂岩和泥岩能够很好地区分出来,砂、泥岩间波阻抗界限为6 700 [(kg·m^-3)·(m·s^-1)](图1b)。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 波阻抗重构前(a)、后(b)的电阻率与波阻抗交汇

Fig.1 Intersection of resistivity and wave impedance before(a) and after(b) wave impedance reconstruction

3.2 构造框架模型

以往波阻抗反演的构造框架模型建立主要是依赖于地震解释的区域性组、段级标准反射层位,而很少应用一些层间的砂组和油层反射层位,同时也完全不考虑断层对反演结果的影响,在构造简单、地层横向变化不大的凹陷性盆地这种影响较小,但在构造复杂、断层发育、地层横向变化较大的断陷盆地就会产生相当大的影响:

1)由于断陷盆地地层横向变化,在大的区域性标准反射层控下内插的井间测井信息就会发生构造误差,为了减小内插测井信息的构造误差,开展层序地层研究,纵向细化地层单元,开展砂组级层序对比,以砂组级的地震反射层面来实现层控构造框架建模,减小区域性标准反射层控下井间内插测井信息发生的构造误差;

2)由于断陷盆地复杂的断裂结构,井间内插时也没有考虑断层两侧内插测井数据的突变,不但影响断层两侧的岩性接触关系,降低断层两侧波阻抗反演精度,而且降低了波阻抗反演剖面对断层的识别能力,为了解决这一问题,将断层转化为地震层位,明确断层与纵向细化地层单元间的接触关系,同地震反射层位一起控制测井信息的井间内插,克服了断层两侧由于井间测井数据内插而引起的测井数据、地震属性特征参数误差,实现了断层两侧内插测井数据以及地震属性特征参数突变。

3.3 小波边缘分析

小波边缘分析具有多尺度分辨的特点,可以有效地区分地震剖面中不同尺度、不同级别的地震特征参数边缘点,利用小波边缘分析获取的地震特征参数中的边缘点可以是断层、岩性体的边缘。

图2和图3分别为DX14井区三维L1036地震剖面以及经小波边缘分析得到的地震特征参数边缘检测特征点剖面。从图3可见,经过小波多尺度边缘分析处理后,纵向地层有效分开,相邻道的不连续性明显,微断裂的断点被突出、地层微小错动特征清晰,反映的地层接触关系更准确,提取的地震特征参数的边缘点可以反映原始剖面中不易识别的地质特征,将边缘点控制的地震特征参数同测井曲线一起参与波阻抗反演建模,不但补充了波阻抗高、低频成分,而且把地层横向间的接触关系完全考虑在波阻抗反演的整个过程中,并且在迭代过程中应用此参数对其进行扰动限制,使AIW波阻抗反演模型反映的地下地质情况更加真实,减少了波阻抗反演的多解性。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 DX14井区三维L1036地震剖面

Fig.2 L1036 seismic profile in DX14 well area

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 DX14井区三维L1036边缘检测特征点剖面

Fig.3 L1036 edge detection feature point profile in DX14 well area

3.4 井—震一体化初始模型建立

在波阻抗反演框架模型的建立过程中,低频分量主要来源于两种途径:一是将测井数据中的高频成分滤去;二是直接应用地震处理的速度谱资料,目前大多数波阻抗反演的低频成分主要是滤去声波测井曲线中的高频成分获得。

同常规波阻抗反演不同,DX14井区AIW波阻抗反演初始模型建立过程中突出井—震一体化、构造—岩性统一:

1)应用钻井拟声波测井曲线获得速度或者阻抗的低频成分,然后在断层参与的高密度层控构造框架控制下横向沿着纵向细化的层序界面外推。为了纵向细化低频模型,保证低频模型井间构造精度,纵向选择梧桐沟组一段3个砂组和梧桐沟组二段1个砂组的底界以及百口泉组底界5个地震反射层做为主要的控制层,当遇到断层时,按该断层转换后的层位与纵向细化的层序界面间的接触关系处理,断层上盘层序界面与断层关系按超覆处理,下盘层序界面与断层关系按剥蚀处理,建立准确的DX14井区AIW波阻抗反演的低频阻抗模型。

2)通过井—震结合精细标定钻井有利储层和纵向细化的层序界面,提取拟波阻抗曲线的低频阻抗分量,在建立的断层参与的高密度层控构造框架控制下内插和外推低频阻抗建立低频阻抗模型后,根据AIW 波阻抗反演基本原理,依据小波边缘分析技术从经过拓频处理的提频三维地震资料中提取反映岩性局部变化的地震属性特征参数,在断层参与的高密度层控构造框架控制下,用小波边缘检测方法得到的边缘信息结合测井沿层位外推的波阻抗来估算反演的初始波阻抗模型,同拟声波阻抗的低频阻抗模型一起建立AIW波阻抗反演的最终初始模型,并且在迭代过程中应用此参数对其进行扰动限制,会使反演结果收获意想不到的效果。

由于断层参与的高密度层控构造框架减小了层间测井信息内插引起的层间构造误差、实现了断层两侧内插测井数据的突变,通过井—震精细标定,依据小波边缘分析技术从双向拓频地震数据中提取地震属性特征参数不但补充了井间内插造成的高频成分缺失,也加入了反映岩性变化信息的地震属性参数,实现井—震一体化、构造—岩性统一的波阻抗反演初始模型。

图4和图5分别为DX14井区3D地震L1036线测井声波阻抗内插外推得到的低频模型以及小波边缘分析得到的地震属性特征参数联合测井声波阻抗、构造框架建立的井—震统一的AIW波阻抗反演模型。对比图4和图5可见,井—震统一的AIW波阻抗反演模型不但包含测井声波阻抗内插外推的低频模型井间横向内插的低频成分,同时也引入了反映井间岩性变化的地震属性数据的小波边缘分析结果,在图4中,由于井间声波阻抗横向连续内插,无法体现井间岩性的变化和地层间的叠置关系,而在井—震统一的AIW波阻抗反演模型图5中,地层纵向上岩性变化可以通过图4中的低频成分来区分,而横向上岩性的变化、地层间接触关系以及断层、断点则通过地震属性特征参数来表征,同测井声波阻抗内插外推的低频模型相比,井—震统一的AIW波阻抗反演模型更能反映地下的实际地质情况。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 L1036线测井阻抗低频模型

Fig.4 Low frequency logging impedance model of L1036 line

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 井震统一的AIW波阻抗反演模型

Fig.5 AIW wave impedance inversion model for well seismology unification

4 应用效果

在滴南凸起DX14井区,梧桐沟组一段是梧桐沟组的主力含气层系,上部砂岩发育,储层较厚,下部泥岩增多。梧桐沟组一段沉积时期,构造活动相对增强,地貌高差构造变化较大,受东部克拉美丽山物源控制,在DX14井区发育一系列扇三角洲前缘砂体,砂体顶部和底部主要发育泥岩、粉砂质泥岩沉积,储层结构具有典型的“泥包砂”特点,单层砂岩厚度一般2.0~7.0 m,前缘砂体累计厚度一般大于为5.0~18.0 m。应用AIW波阻抗反演方法,充分考虑滴南凸起DX14井区复杂的构造特点、多变的沉积环境形成的“小而薄”扇三角洲前缘砂体以及三维地震资料波组反射特征、地下全方位地震波场的分布特征与变化因素,依据小波边缘分析技术从经过拓频处理的三维地震资料中提取反映岩性局部变化的地震属性特征参数,在断层参与的高密度层控构造框架控制的基础上同测井数据一起参与初始模型建立,并反复迭代修改反演道的波阻抗模型,最终得到滴南凸起DX14井区高精度的地震反演波阻抗资料。

图6为DX14井区在同一个构造地质框架模型基础上,分别采用AIW波阻抗反演和常规基于模型波阻抗反演方法得到的过DX14井波阻抗反演剖面,从图6可见,AIW波阻抗反演剖面的分辨率高于常规基于模型波阻抗反演结果和常规叠后地震剖面,AIW波阻抗反演结果反映的砂层间相互接触关系不但清晰,而且明显反映出了地层岩相、岩性的横向变化,地层岩性信息也更加丰富。通过过DX14井的波阻抗反演剖面也可以看出梧桐沟组一段的砂体横向分布范围较小、不连续,呈现出扇三角洲前缘砂体“小而薄”的特点,符合在DX14井区发育扇三角洲前缘砂体的地质认识。纵向上准确预测了DX14井梧桐沟组一段发育的5套砂体,其中3 532.00~3 536.00 m井段厚度为4.0 m的含气砂体同3 542.00~3 550.00 m井段厚度为8.0 m的含油砂体呈明显的叠置关系,而3 522.00~3 524.00 m井段厚度为2.0 m的含气砂体、3 562.00~3 572.00 m井段厚度为10.0 m的含气砂体则呈现明显的透镜状砂体,向上倾和下倾方向尖灭,反演结果并与实钻结果高度吻合,而常规基于模型波阻抗反演结果不但储层分辨率不高,而且无法反映砂体间的叠置关系。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 AIW波阻抗反演(a)与常规基于模型波阻抗反演(b)剖面对比

Fig.6 Comparison between AIW wave impedance inversion section and conventional model-based wave impedance inversion section

利用AIW波阻抗反演和常规基于模型波阻抗反演结果,分别对梧桐沟组一段下砂组砂体顶、底界面进行精细井震对比、追踪解释,构造成图,以梧桐沟组一段下砂组顶、底地震反射界面为层段时窗界限,按照砂岩和泥岩间波阻抗界限6 700 [(kg·m^-3)·(m·s^-1)]预测砂体分布结果如图7,从中可见常规基于模型波阻抗反演资料分辨率较低,砂体的地震响应特征不明显,砂体边界及岩性体间的关系不清(图7b)。而利用AIW波阻抗反演结果可以明显看出:受东部克拉美丽山物源控制,在DX14井区发育一系列扇三角洲前缘砂体,砂体与断裂相配合在DX14井区形成断层—岩性和岩性油气藏(图7a)。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 DX14井区梧桐沟组一段P₃w $t_{1}^{2}$ AIW波阻抗反演(a)与常规基于模型波阻抗反演(b)波阻抗平面

Fig.7 P₃wt²₁ AIW wave impedance inversion(a) and conventional model-based wave impedance inversion(b) wave impedance plan of the first member of Wutonggou Formation in DX14 well area

目前DX14井区完成各类钻井133口,其中直井76口,水平井57口,以砂岩和泥岩间波阻抗界限6 700 [(kg·m^-3)·(m·s^-1)]为标准对梧桐沟组一段P₃w $t_{1}^{2}$ 砂岩进行预测,通过与实钻砂岩厚度对比,井点处预测砂岩厚度的绝对误差均小于0.60 m,相对均误差小于2.84%。近期完钻的D408井钻前预测梧桐沟组一段P₃w $t_{1}^{2}$ 砂岩厚度为12.6 m,通过D408井钻后对比,该井梧桐沟组一段P₃w $t_{1}^{2}$ 实钻砂岩厚度为12.1 m,绝对误差为0.5 m,证明AIW

波阻抗反演结果的准确度较高,预测精度完全满足精细储层预测、精细油藏描述和精细勘探开发的地质需求。

5 结论

1)应用小波边缘分析技术建立的初始波阻抗模型更加接近真实地下地质情况,提高了反演的真实性,利用所得的特征点信息结合测井波阻抗对反演的约束起到了很好的效果,为储层的研究提供了指导作用。

2)AIW波阻抗反演技术将测井波阻抗、地震特征边缘信息作为约束条件建立的初始模型,弥补了常规反演过程中缺少井间的局部构造、岩性变化的地震信息,增加了波阻抗反演过程中高频成分,同常规反演相比,AIW波阻抗反演结果不但能够准确反映井间岩性变化,而且更接近实际地质情况。

3)在断层参与的高密度层控构造框架控制下,测井波阻抗和地震属性特征参数联合建模实现了井—震一体化、构造—岩性统一,但应用的地震属性特征参数比较单一,多属性联合建模仍是波阻抗反演的研究方向。

4)利用基于小波边缘分析的井—震联合建模AIW波阻抗高反演方法,对准噶尔盆地陆梁隆起DX14井区梧桐沟组一段含气砂岩储层进行精细反演和预测,刻画了扇三角洲前缘砂体,预测砂体厚度的相对误差均小于2.84%,预测精度高,完全适用于准噶尔等复杂断陷盆地以及复杂目标区构造—岩性等各种隐蔽性油气藏的精细储层预测和精细油藏描述。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

G F

, Zhou

, Zhao

Potential risks of spectrum whitening deconvolution—Compared with well-driven deconvolution

[J] .Petroleum Science, 2009, 6(2):146-152.

DOI:10.1007/s12182-009-0023-y URL [本文引用: 2]

[2]

王江, 涂国田, 王杰.

基于载波调制的高分辨率地震双向拓频技术及其应用

[J]. 石油物探, 2021, 60(6):954-963.

DOI:10.3969/j.issn.1000-1441.2021.06.009 [本文引用: 1]

针对断裂复杂、沉积相变化快,含油砂体分布范围小、厚度薄,预测难的问题,将无线电领域的频率域“载波调制”理念引入到地震资料处理中,用基于最大熵(Max Entropy)谱估计的最小二乘上、下三角阵分解(LS-LUD)算法,研发了ButHRS地震高分辨率双向拓频处理技术,将改进型Butterworth带通子波谱的包络作为目标函数,通过L1模计算方法得到谱调制加权因子,将短时窗数据的LS-LUD谱分解得到的振幅包络向改进型Butterworth带通子波谱的包络形态逼近,得到短时窗地震数据高分辨率谱。该技术不但在低频和高频两个方向上拓展了地震频带,有效压制了地震子波旁瓣,而且在提高地震资料分辨率的同时,保持了地震资料的信噪比、相对振幅关系和时频特性,提高了地震资料识别小薄砂体的能力。将该方法应用于乌尔逊断陷铜钵庙断裂带的储层预测,经过三维地震拓频处理,目的层频带拓宽了40Hz(由10～55Hz拓宽到8～95Hz),精细刻画了铜钵庙断裂带南屯组扇三角洲前缘砂体的三维轮廓和有利储层分布,并通过综合分析认为铜钵庙断裂带受东部物源控制,发育扇三角洲—湖底扇—半深湖(深湖)沉积体系。预测结果符合地质认识,与钻井结果吻合度较高。

Wang

, Tu

G T

, Wang

High-resolution seismic bidirectional frequency extension based on carrier modulation

[J]. Geophysical Prospecting for Petroleum, 2021, 60(6):954-963.

DOI:10.3969/j.issn.1000-1441.2021.06.009 [本文引用: 1]

<p>This work was motivated by the need to overcome the issues related to complex faulting,rapid sedimentary transformation,narrow distribution and low thickness of oil-bearing sand bodies,and the difficulties in their prediction.To overcome these issues,for the first time,the concept of “carrier modulation” in the radio frequency domain was introduced to seismic data processing,and a ButHRS high-resolution seismic bidirectional frequency extension method was developed using the LS-LUD algorithm based on the estimation of the maximum entropy spectrum.In this method,the envelope of the improved Butterworth band signal spectrum is taken as the objective function,and the spectral modulation weighting factor is obtained through the L1-mode calculation method.The amplitude envelope of the LS-LUD spectrum decomposed from short-time window data is approximated to the envelope of the improved Butterworth band signal spectrum,thereby obtaining a high-resolution spectrum for the short duration seismic data.The HRSBFE method not only expands the seismic frequency band in both the low-frequency and high-frequency domains but also effectively suppresses the sub-wave lobe of the seismogram.This allows the signal-to-noise ratio,relative amplitude relationship,and time-frequency characteristics of the seismic data to be maintained,while improving their resolution.This is especially beneficial for identifying small and thin sand bodies.The proposed method was used for reservoir prediction in the Tongbomiao fault zone in the Wuerxun fault depression.The frequency band of the target layer could be widened by 40Hz (from '10~55Hz’ to ‘8~95Hz’) via frequency extension processing of three-dimensional (3D) seismic data.Based on these data,the 3D profile of fan delta front sand bodies and the distribution of favorable reservoirs in the Nantun Formation in the Tongbomiao structural belt could be described.It was found that fan delta-sub-lacustrine fan semi-deep- and deep-lake sedimentary systems developed in the Tongbomiao tectonic belt controlled by the eastern source area.The predicted results were in accordance with the geological trends,and exhibited good consistency with data from drilled cores.</p>

[3]

王华忠, 盛燊.

走向精确地震勘探的道路

[J]. 石油物探, 2021, 60(5):693-708,720.

DOI:10.3969/j.issn.1000-1441.2021.05.001 [本文引用: 1]

勘探地震领域已扩展到复杂地表、复杂构造、复杂油藏和深层目标等。地震数据采集技术发展到了宽方位、高密度、宽频带(“两宽一高”)地震数据采集阶段;地震波成像技术发展到了贝叶斯(Bayes)参数估计理论下的全波形反演(FWI)和最小二乘逆时深度偏移成像(LS_RTM);油藏描述发展到了综合信息利用和最佳判定阶段;地震勘探技术已经发展至全新的阶段。横向缓变的层状介质假设、地表一致性假设、射线理论波传播和Zoeppritz方程界定了上一代精确地震勘探的方法技术及其适用性,上一代精确地震勘探以高分辨率地震子波作为成像处理的核心目标,并据此开展薄层油气藏的识别、描述与评价。而描述任意介质中地震波传播的波动理论和贝叶斯参数估计理论构成了新一代高精度地震勘探的理论基础。“两宽一高”的地震数据采集技术和更高精度的子波处理;基于高维、字典基和稀疏特征表达的信号处理技术(解决去噪、数据规则化、数据压缩、去混叠等问题)、建立更精确速度和Q值模型以及估计宽带反射系数的特征波反演成像技术、宽带波阻抗成像技术和基于信息综合的人工智能油藏描述技术代表了走向精确地震勘探的未来方向。

Wang

H Z

, Sheng

Pathway toward accurate seismic exploration

[J]. Geophysical Prospecting for Petroleum, 2021, 60(5):693-708,720.

DOI:10.3969/j.issn.1000-1441.2021.05.001 [本文引用: 1]

The domain of seismic exploration has expanded to cover complex topographies,complex subsurface structures,complex reservoirs,and deep geological targets.The new generation of accurate seismic exploration methods entail wide-azimuth,wide-frequency band,and high-density seismic data acquisition.In these methods,full-wave inversion,and least-square reverse time migration methods based on Bayes estimation theory have become the standard for imaging,and reservoir characterization is implemented with integrated information utilization and optimal inference.Evidently,seismic exploration has entered a new era.For the last-generation accurate seismic exploration,the assumption of consistence for layers and surfaces,the ray-theory wave propagation,and Zoeppritz equations define the scope of seismic data processing and its applicability.A high-resolution seismic wavelet is the key objective of seismic data acquisition and imaging methods,based on which the identification,characterization,and evaluation of thin reservoirs is carried out.In new-generation methods,the theoretical basis is enriched by a number of tools,namely:the Bayesian estimation theory and wave equations that describe wave propagation in complex media.The seismic data acquisition with wide-azimuth,wide-frequency band,and high-density sampling;accurate seismic wavelet processing;amplitude and phase preserving signal processing based on high dimension,dictionary bases and sparse representation;characteristic wave inversion for accurate velocity model building;Q model building and wideband reflectivity estimation;wideband impedance model building;and reservoir characterization with integrated information utilization and artificial intelligence,all of these methods provide novel pathways for accurate seismic exploration.

[4]

纪永祯, 张渝悦, 朱立华,

等.

基于SBL-WVD 的地震高分辨率时频分析

[J]. 石油物探, 2020, 59(1):80-86,107.

DOI:10.3969/j.issn.1000-1441.2020.01.009 [本文引用: 1]

时频分析是地震数据处理和解释过程中重要的数学工具,其精度和分辨率决定了后续处理和解释成果的质量。提出了一种结合贝叶斯学习方法(sparse bayesian learning,SBL)和魏格纳威利分布(wigner-ville distribution,WVD)的两步高分辨率时频分析方法。第一步基于构建的雷克子波库和贝叶斯学习方法将地震数据分解为子波的线性组合;第二步通过求取子波的魏格纳威利分布获得地震数据的时频分布。其中,贝叶斯最大后验概率和第二型最大似然概率通过迭代求解。贝叶斯学习方法可以用最少数量的、具有不同主频和相位的雷克子波重构地震数据,并同时有效压制随机噪声。求取、分解子波的魏格纳威利分布可有效避免交叉项干扰,分辨率高。模拟数据和实际数据实验结果均验证了方法的正确性和有效性。与常规基于Gabor变换和匹配追踪算法的时频分析方法相比,该方法具有更高的精度和分辨率,有利于后续处理和解释研究。

Y Z

, Zhang

Y Y

, Zhu

L H

et al.

High-resolution seismic time-frequency analysis based on sparse Bayesian learning combined with Wigner-Ville distribution

[J]. Geophysical Prospecting for Petroleum, 2020, 59(1):80-86,107.

DOI:10.3969/j.issn.1000-1441.2020.01.009 [本文引用: 1]

The precision and resolution of time-frequency analysis determines the quality of processing and interpretation.Here,a two-step high-resolution time-frequency analysis method is presented,which combines sparse Bayesian learning (SBL) and Wigner-Ville distribution (WVD).First,on the basis of a built Ricker wavelet library and SBL,the seismic record is decomposed as a linear combination of wavelets;secondly,the WVD of the wavelets are calculated to obtain the time-frequency distribution of the seismic record.By iteratively solving a Bayesian maximum posterior and a type-II maximum likelihood problem,the Bayesian learning method can find the minimum number of Ricker wavelets with different frequency and phase that are needed to reconstruct the seismic record.By doing so,the method can effectively suppress random noise.Furthermore,the WVD of the decomposed wavelets can effectively eliminate cross-term interference,thereby achieving high-resolution imaging.An application to both synthetic and actual data samples showed higher precision and resolution of the proposed method in comparison with that of conventional time-frequency analysis methods based on Gabor transformation and the MP algorithm,which facilitates subsequent processing and interpretation.

[5]

肖张波, 雷永昌, 于骏清,

等.

基于宽频资料的扩展弹性阻抗反演方法在陆丰22洼陷低勘探区古近系岩性预测中的应用

[J]. 物探与化探, 2022, 46(2):392-401.