地震分频技术在扇三角洲前缘砂体预测中的应用——以海拉尔盆地乌尔逊断陷铜钵庙断裂带南屯组为例

时间：2024-07-28

王江，陈沫，王杰，刘霜，周豪

1 中国石油大庆油田有限责任公司勘探开发研究院；2 中国石油大庆油田有限责任公司采油一厂；3 中国石油新疆油田公司采油一厂

0 前言

海拉尔盆地为中、新生代的多旋回、叠合式、断-拗转换型盆地［1］，盆地经过多期改造和建造，断裂发育、构造复杂。乌尔逊断陷位于海拉尔盆地中部断陷带，中石油第四次资源评价的资源量为2.52×108t,现已提交探明储量累计0.38×108t，探明率仅为15.1%，表明乌尔逊断陷具有较大的勘探潜力。随着乌尔逊断陷勘探开发的逐步深入，勘探开发区带由隆起带转向了凹陷区，油藏类型由构造油藏转向了岩性油藏。下白垩统铜钵庙组、南屯组和大磨拐河组等目的层发育扇三角洲前缘砂体，砂岩相带窄、相变快、非均质性强，砂泥岩以薄互层组合，储层厚度薄，纵、横向岩性变化大。砂岩储层横向连续性差，导致区块优选、井位部署及滚动钻井评价难度加大。针对乌尔逊断陷扇三角洲前缘砂泥岩薄互层组合，曾经开展过多轮次地震预测工作，但受地震资料及单一预测方法限制，研究区扇三角洲前缘砂体、湖底扇砂体的地震响应特征不明显，砂体边界不清，预测精度无法满足现阶段岩性油藏勘探评价的地质需求。因此，为了解决岩性油藏薄窄砂体刻画难题，在借鉴前人研究成果的基础上，通过井震结合开展精细储层预测，形成了适合于复杂断陷岩性目标区的地震匹配追踪时频处理解释技术。

地震分频处理技术是地震属性分析中的一种频率域解释方法［2］。国内外学者对分频处理、解释研究可追溯到20 世纪80—90 年代，当时只是应用分频处理资料沿地震反射层的一定时窗内提取地震属性，进行定性的沉积环境、沉积相分析，但受傅里叶变换固定时窗的影响，限制了在薄储层领域的应用。前人［3-5］在分析地震资料时，引入了在时间域和频率域都可以表征信号局域化特征及变时窗特点的小波变换，其分频处理资料得到了明显改善，优势频率地震剖面的储层分辨能力得到了提高，也出现Gabor小波、离散Mallat小波等各种基于小波决定的算法。同傅氏变换相比，小波变换的分频效果更好，但数学运算上不但要求正交的小波函数，同时在一定程度上受地震处理时窗的影响。Stockwell等［6］提出的S 变换时频分析是介于短时傅里叶变换和连续小波变换之间的、应用小波变换的多尺度时频分析方法，S 变换时频分析在克服傅里叶变换不能调节分析窗口频率的同时还可以进行相位校正。应用地震分频资料开展碎屑岩、碳酸盐岩储层预测，以及小断层识别解释和高频沉积旋回研究，取得了不同程度的研究效果［6-9］。

本文应用地震匹配追踪时频处理技术，通过一系列小波构建初始模型，经过迭代优化逼近原始信号的瞬时谱，将时间域的地震信号分解成频率域的离散调谐体，应用频率域的离散数据开展扇三角洲前缘砂体预测，精细刻画扇体的三维轮廓和内部结构，应用研究成果指导岩性油藏井位部署及滚动评价，在乌尔逊断陷铜钵庙断裂带识别的9 个有利岩性体中的7 个实施勘探部署，其中在5 个砂体中获工业油流，在2个砂体中获得低产油流，预测砂体顶界构造深度相对误差均小于0.3%，预测砂岩厚度的绝对误差均小于1 m。

1 地震匹配追踪时频处理技术

地震匹配追踪时频处理是一种基于投影追踪、逐步递推的小波算法，该算法是将地震信号分解为一系列反射时间不同、振幅各异的地震子波，对分解得到的子波进行时频分析，原始地震道频谱就是所有子波频谱叠加的结果。

2006 年Liu［10］基于Ricker 子波匹配追踪算法通过复地震道分析技术得到三瞬地震属性（瞬时振幅、瞬时相位和瞬时频率），以最大瞬时振幅处的地震反射时间作为时间延迟的初始值，得到延时处的瞬时相位和瞬时频率。

主频为fj的零相位Ricker 子波在时间域的表达式为：

式中:fj为地震主频，Hz；t为地震反射时间，s；WR（t,fj）为时间域Ricker子波；j为迭代次数。

在频率域的傅氏频谱为：

式中:f为地震频率，Hz；WR（f,fj）为主频为fj的频率域Ricker子波傅氏频谱；j为迭代次数。

地震信号的稀疏表达式为：

式中:s(t)为时间域带限地震信号；Rs(m)(t)是匹配追踪后的残差；m为通过设置的阈值来判断迭代终止条件而确定的迭代原子个数。原子个数通过时间迟延tj（单位：s）、地震主频fj（单位：Hz）、相关系数aj和相位φj等4个处理参数来控制。

s(t)匹配追踪时频谱u(t,f)的解析式为：

式中：env是WR的瞬时振幅。

2 技术应用及效果

乌尔逊断陷北部地震数据频带范围为8～60 Hz,主频（f）为28Hz,砂岩平均层速度（v）为3 897～4 012 m/s，根据地震可分辨的极限厚度ΔZ=v/(4×f) 计算,可识别的最薄砂体厚度为34～36 m。可见，目前地震资料无法识别研究区20～30 m 厚的砂体,利用叠后地震波阻抗反演可识别的砂岩厚度仅为几十米, 因此, 对于薄层的单砂体, 利用传统的地震波阻抗反演方法难以分辨。在频谱分解之前首先对地震资料进行脉冲反褶积提频处理, 将地震的频宽拓展到8～80 Hz, 提频后的地震带宽比常规地震数据带宽在高频端高出20 Hz,低频端保持不变,主频从原先的28 Hz 拓展到37 Hz, 这样拓频后的地震数据比常规地震数据丰富了很多细节。

2.1 分频处理的优势频率分析

传统的地震分辨率是指时间域内的地震储层分辨能力，分辨的极限厚度是λ/4(λ为波长)，而频谱分析突破了传统地震分辨率的限制，在频率域内对储层进行分辨，每个薄储层产生的地震反射在频率域都有一个特定的优势频率与其对应。图1为纵向砂岩厚度分别为3 m、5 m、10 m、20 m，横向频率f分别为30 Hz、40 Hz、50 Hz、60 Hz、70 Hz 和80 Hz 的时频正演理论模型，从图中可以看到高频调谐能量反映的是薄层的调谐响应，而低频调谐能量是厚层的调谐响应，各单一频率的子波对不同厚度储层的调谐响应是不同的，不同厚度的砂岩对应着不同的优势频率。因此经分频处理得到的单频三维地震资料，在有效地震频带内降低了子波叠覆和鸣振效应，突出了不同地质体对不同频率成分地震资料的敏感响应，而常规叠后地震资料是对整个地震频带范围的地震数据成像，无法满足对所有地质体都达到最佳的成像效果。因此，利用频率域的振幅谱可以识别沉积相带变化、沉积扇体三维轮廓以及薄储层的时间厚度。

图1 不同厚度砂岩的时频正演理论模型Fig.1 Theoretical model of time-frequency forward modeling of sandstone with different thickness

根据对乌尔逊断陷北部主要含油目的层储集特征与测井相的分析，目的层砂体厚度(ΔZ)一般为20～26 m (大磨拐河组20～25 m，南屯组23～30 m)，砂岩平均层速度大磨拐河组为3 897 m/s、南屯组为4 012 m/s。据f=v/(4×ΔZ)计算砂岩储层的最小调谐频率，大磨拐河组为39～48 Hz，南屯组为38～43 Hz，即识别大磨拐河组、南屯组有利含油砂体的地震优势频率理论上为38～48 Hz。图2 为乌尔逊断陷北部S46 井实际钻井的时频分析对比图，可见26 m 砂组对应的优势频率为40 Hz 左右，该频率下表现的调谐响应最强，而8～15 m 砂组对应的优势频率为60～80 Hz。优势频率分析为单频剖面处理、利用优势频率识别砂体提供了理论基础和实践根据。

图2 铜钵庙断裂带S46井不同厚度砂岩的时频分析对比图Fig.2 Time frequency analysis and comparison of sandstone with different thickness of Well S46 in Tongbomiao fault zone

2.2 优势频率的单频剖面特征

针对乌尔逊断陷复杂的地质条件，根据大磨拐河组和南屯组有利含油储层38～48 Hz 的优势频率范围，对乌尔逊断陷北部铜钵庙断裂带三维地震资料应用匹配追踪时频处理技术开展分频处理。首先将地震道与小波库内所有小波互相关，再将从地震道中去除相关性最大的子波投影所得到的残差与小波库互相关，再去除最大相关子波的投影，如此递推迭代，直到剩余能量达到某一给定的门槛值，迭代过程收敛，得到一系列离散频率的调谐数据体。然后通过地震四维时频频谱分析，利用代表某一特定主频的单频能量的单频体、峰值振幅和峰值频率等多种时频地震属性来预测有利储层。

从图3 常规地震偏移剖面（图3a）与42 Hz 单频偏移剖面（图3b）的对比图中可见，42 Hz 单频偏移数据反映的S46 井砂体边界清晰(见图3b 虚线标注部分)，反射结构波组特征清楚，地层接触关系明确，而常规叠后地震资料是对整个地震频带范围的地震数据成像叠加，无法满足所有不同尺度地质体都能达到最佳的成像效果。

图3 铜钵庙断裂带三维常规偏移剖面与单频偏移剖面对比Fig.3 Comparison of 3D seismic conventional migration profile and single frequency migration profile in Tongbomiao fault zone

2.3 分频处理技术的地质效果

地震属性是对地震反射资料的描述和量化，是地震资料中所包含全部地质信息的子集，能揭示出地震数据体中的沉积、岩性和储层等信息［11-15］。分频处理增强了薄储层频率响应，提高了地震资料的分辨率，可精细地反映地震属性异常变化，增强了薄储层的识别能力。

2.3.1 储层分布定性预测效果

在下白垩统铜钵庙组—南屯组沉积时期，乌尔逊断陷东部受北东向控陷断层——铜钵庙断层控制［11］，在半深湖—深湖背景下发育扇三角洲前缘砂体，与早期断层配合形成构造-岩性、岩性油气藏［12］。但是，研究区地震资料分辨率较低，砂体的地震响应特征不明显，叠后、叠前时间偏移资料均无法确定砂体边界及岩性体间的关系。为了验证单频地震资料的预测效果，分别应用常规地震资料和单频地震资料，考虑了储层厚度、位置及层速度，以等时地层格架为约束，在大磨拐河组一段、南屯组一段3砂组（简称南屯组1-3砂组）、南屯组2-1砂组出油层位的上、下30 ms 选取时窗，由线到面提取敏感地震属性。图4a 为乌北铜钵庙断裂带南屯组1-3砂组常规地震均方根振幅属性图，从图中可见，振幅类属性基本可以反映扇三角洲前缘和湖底扇储层的展布规律，但储层边界不清楚。通过与钻井对比发现，常规资料单一属性在预测储层分布上存在局限性。例如，钻井证实W1 井、S11 井、S22 井和S37 井南屯组1-3 砂组岩性以泥岩夹泥质粉砂岩为主，砂岩不发育，应该具有弱振幅、中高频的反射特征，但在常规地震资料均方根振幅属性上，这4口井均落入了强振幅指示的砂岩发育区内。而且，钻井证实的S35 井湖底扇砂体和S45 井扇三角洲前缘砂体间的均方根振幅属性边界不清，也影响了岩性圈闭的识别与刻画。图4b为铜钵庙断裂带南屯组1-3砂组42 Hz 单频地震均方根振幅属性图，通过强振幅属性与钻遇砂岩一致性的对比分析发现：区域内28 口探井和评价井中，有14 口井在南屯组1-3 砂组钻遇砂岩储层，在图4b上显示均方根相对振幅值大于400；另外14 口井南屯组1-3 砂组的砂岩储层不发育，在图4b 上显示均方根相对振幅值小于400。可见单频地震资料预测砂体的符合率为100%，明显优于常规地震资料的预测结果。

图4 铜钵庙断裂带南屯组1-3砂组均方根振幅地震属性图Fig.4 Root mean square amplitude attribute map of 1-3 sand group of Nantun Formation in Tongbomiao fault zone

通过对图4b 单频数据体地震属性精细解释认为：在乌尔逊断陷铜钵庙断裂下降盘发育扇三角洲前缘砂体，中部发育湖底扇砂体，是岩性油藏的发育区。在单频数据振幅属性切片上，扇三角洲前缘砂体和湖底扇砂体均表现为明显的强振幅特征，由此推断乌尔逊断陷北部铜钵庙断裂带受东部物源控制，发育扇三角洲—湖底扇—半深湖（深湖）沉积体系。

2.3.2 储层厚度定量预测效果

根据Partyka等［16］的薄层分析理论，振幅谱伴随频率响应的陷频周期Pt（单位：Hz）与薄层时间厚度t（双程旅行时，单位：s）的关系为t=1/Pt；而t=2△Z/v,其中△Z为薄层厚度（单位：m），v为薄层速度（单位：m/s），这样薄层厚度△Z与Pt的关系为△Z=v/（2×Pt）；将首波峰值频率f0（单位：Hz）=Pt/2 代入时间厚度和薄层厚度等式后，以上两式则变为：

式中首波峰值频率f0可以通过地震资料分频分析得到，这样通过公式（6）就可以计算薄层厚度。

图5a 和图5b 分别为S9 井区南屯组1-3 砂组砂体顶界构造图和根据分频处理资料得到的砂体厚度图，结合地震属性图分析，S9 井区砂体是受东部物源控制在铜钵庙断裂下盘形成的扇三角洲前缘砂体，北西和南东方向受断层控制，沿铜钵庙断裂走向向北东和南西向尖灭，砂体最大厚度可达30 m。针对在南1-3 砂组中预测的4 个砂体部署的13 口探井和评价井均钻遇所预测的储层，钻井实际与钻前预测相符，预测砂体顶界构造深度相对误差均小于0.35%（表1），预测砂岩厚度的绝对误差均小于1 m，相对误差均小于4%，预测精度满足精细储层预测和勘探开发的生产需求。

通过地震分频处理与解释，在铜钵庙断裂带大磨拐河组一段、南屯组1-3 砂组、南屯组2-1 砂组和铜钵庙组共识别有利岩性体9 个，累计面积49.4 km2。目前已钻砂体7个，其中在5个砂体中获工业油流，在2个砂体中获得低产油流，其中S46井在大磨拐河组一段、南屯组和铜钵庙组分别获得4.81 t/d 、11.68 t/d和9.6 t/d的高产工业油流。

图5 铜钵庙断裂带S9井区南屯组1-3砂组砂体顶面构造图与砂岩厚度图Fig.5 Sand body top structural map and sandstone thickness map of 1-3 sand group of Nantun Formation in S9 well block of Tongbomiao fault zone

表1 铜钵庙断裂带南屯组1-3砂组含油砂体预测精度统计表Table 1 Statistical table of prediction accuracy of oil-bearing sandbodies in sand group 1-3 of Nantun Formation in Tongbomiao fault zone