基于ADM谱反演的高分辨率裂缝预测技术研究及应用

张　华,贺振华,李亚林,张晓斌,何光明,李德珍

(1.“油气藏地质与开发工程”国家重点实验室,四川成都610059;2.成都理工大学地球物理学院,四川成都610059;3.中国石油天然气集团公司川庆钻探工程有限公司地球物理勘探公司,四川成都610213;4.中国石油天然气集团公司山地地震技术试验基地,四川成都610213)



基于ADM谱反演的高分辨率裂缝预测技术研究及应用

张华1,2,3,4,贺振华1,2,李亚林3,4,张晓斌3,4,何光明3,4,李德珍3,4

(1.“油气藏地质与开发工程”国家重点实验室,四川成都610059;2.成都理工大学地球物理学院,四川成都610059;3.中国石油天然气集团公司川庆钻探工程有限公司地球物理勘探公司,四川成都610213;4.中国石油天然气集团公司山地地震技术试验基地,四川成都610213)

从谱反演基本原理出发,推导出多层稀疏反射系数目标函数,利用基追踪交替方向法(Alternating Direction Method,ADM)谱反演算法求取高精度反射系数,与谱模拟解析法提取的宽频子波褶积重构,得到保持低频信息的高分辨率剖面,利用改进型特征值相干体算法对此剖面进行裂缝预测。理论测试和实际资料应用结果表明,ADM谱反演高分辨率技术能得到极性、位置准确性较高的反射系数,而重构后的宽频剖面分辨率明显提高、同相轴连续性更好、断点更清晰、断层更明了,相干体切片精度更高、细节刻画更清晰,可作为裂缝预测的重要技术手段。

分辨率;ADM;裂缝预测;相干体;反射系数

相干体裂缝特征描述是裂缝性油气藏勘探的重要手段之一,而获取用于裂缝预测所需高分辨率数据体以及裂缝参数的求取一直是裂缝油气藏勘探中重点研究的问题。

薄层反射系数谱反演[1-3]是近些年发展起来的一种能有效分辨薄层、提高分辨率的方法。该方法的主要特点是不需要任何先验信息和数学假设条件,利用反射系数奇偶分解理论和部分频率谱信息反演地层稀疏反射系数,进一步提高分辨低于调谐厚度的地震薄层能力。PURYEAR等[4]、CHOPRA等[5]对谱反演算法理论进行了详细推导,模型测试和实际资料应用结果证明该方法能稳健地分辨薄层厚度。YUAN等[6]详细分析了谱反演的不适定性,并提出了快速的粒子群Levenberg-Marquardt混合反演算法,提高了反演的精度和效率。柴新涛等[7]利用最小二乘QR分解(Least Squares QR,LSQR)算法进行谱反演,得到较好的模型测试效果。陈祖庆等[8]利用地震资料的谱信息和基追踪算法进行稀疏脉冲反射系数谱反演,模型试验结果表明该方法能分辨原始数据中无法识别的薄层。

三维相干体技术是20世纪90年代后期兴起的一项十分有效的地震解释技术,该技术主要从相邻地震道相互之间的相干性出发,给出一种定量描述,能突出断层处地震波的变化情况,形成相干值异常区域,从而预测断裂展布。BAHORICH等[9]提出在振幅上进行互相关相干分析C1算法,但该算法计算出的断点非相干性点模糊,噪声大、能量弱。MARFURT等[10]提出基于相似系数相干C2算法,具有稳定性好、断层刻画精度高的特点,但横向分辨率低的劣势无法改变。GERSZTENKORN等[11]提出基于本征结构分析的C3算法,利用Hilbert变换计算相干体,提高其稳定性和抗噪性,但该算法用于高陡复杂地区时,如不考虑倾角则效果不如第二代算法的效果好,如考虑倾角则因计算效率低而无法应用于生产。张军华等[12]、王玉学等[13]、刘振峰等[14]、齐晴等[15]对C3算法及其优化算法在断层解释、裂缝预测方面进行了讨论和改进,取得不错的应用效果。

本文在前人研究的基础上,利用基追踪算法是全局优化算法,能使地震信号的稀疏分解更加稳定的优势[16],进行基于基追踪算法理论的交替方向法(ADM)谱反演算法研究,求取高精度反射系数,再利用该反射系数与含低频的宽频子波重构得到分辨率更高的三维地震数据体,在该数据体上利用改进的C3相干算法进行高精度断层识别和裂缝预测,避免了解释人员在断层解释和组合的随意性,使断层解释、裂缝预测的精度大大提高。

1　ADM谱反演的高频拓展原理

1.1稀疏谱反演基本原理

根据不含噪声的褶积模型,地震合成记录在频率域可表示为:

(1)

式中:S(f),W(f)和R(f)分别表示频率域地震记录、地震子波、反射系数。

稀疏谱反演利用频率域地震记录与子波的丰富信息,采用合适的反演算法消除子波影响,从记录中得到所需的反射系数。利用反射系数奇偶分解原理可得到频率域下谱反演目标函数表达式:

(2)

式中:ae和ao分别代表反射系数R(f)的偶、奇分量比例。

为建立更符合实际地层情况的反射系数模型,必须推导多层反射系数模型下的目标函数,由反射系数奇偶分解原理得到多层反射系数偶、奇分量,表示为:

(3)

式中:T表示反射系数时间厚度；re和ro分别代表反射系数r(t)的偶、奇分量。

对(3)式进行傅里叶变换得到其频率域实部与虚部表达式:

(4)

利用公式(2)和公式(4)得到多层反射模型的目标函数表达式为:

(5)

式中:Δt表示时移量。

为了便于求解,将公式(5)写成矩阵形式:

(6)

其中,

Ae=

Ao=

be=[re(t1)re(t2)…re(tN)]T

bo=[ro(t1)ro(t2)…ro(tN)]T

(7)

1.2ADM谱反演算法

求解目标函数((6)式)的算法有很多,如共轭梯度法、模拟退火法、匹配追踪法、基追踪法等,这些算法及混合算法各有优劣,其中,基追踪去噪算法由于其数值计算过程中表现出的稳健性、高效性和精确性越来越受到地球物理学者的青睐,本文采用以基追踪去噪理论为基础的ADM谱反演算法来求解目标函数最优解。

由于该目标函数((6)式)有多个全局最优解,为了降低求解问题时的多解性和不确定性,需要结合实际地震反射系数为稀疏的这一假设条件。因此,加上一个稀疏化约束,将目标函数改为范数求解形式为:

(8)

根据基追踪去噪算法理论,可将目标函数(8)的求解问题转化为:

(9)

式中:υ为极小权重值;G为复数集合；n为向量x的维数。

如果引入变量r∈Gm,那么公式(9)等价于公式(10):

(10)

式中:m为向量r的维数。

公式(10)对应的增广拉格朗日子问题可表示为:

(11)

式中:yT表示乘子的共轭转置运算;β为罚参数。

利用公式(11)进行基追踪ADM谱反演算法,实现步骤如下。

1) 令k=0对rk,xk,yk赋初始值,并给定υ,β常数值,然后进行算法运算,如果满足终止准则,则完成运算,否则,进行步骤2);

2) 令x=xk,y=yk,求解r的子问题得到rk+1:

(12)

3) 令r=rk+1,y=yk则关于x的极小化问题公式(11)等价于公式(13):

(13)

那么可以通过公式(14)近似求解来完成对公式(13)的精确求解:

(14)

式中:τ为大于零的邻近参数,gk为x=xk时二次项的梯度,则公式(14)通过公式(15)求解得到xk+1:

(15)

式中:“°”表示逐元素相乘;Shrink(·)表示一维收缩算子;sign表示符号函数。

4) 令x=xk+1,r=rk+1,求解y的子问题得到yk+1:

(16)

式中:γ为大于零的常数。

5) 令k=k+1,重复步骤1)至步骤4)进行迭代运算。

1.3保持低频信息的宽频子波提取

为得到保持低频信息的高分辨率资料,需对求取的反射系数与保持低频信息的宽频子波重构。

零相位子波是具有较高分辨率的子波,但采用不同主频的零相位子波,主频的变化会使其低频相应发生变化,重构后无法保持提频前的低频信息,

所以必须获取能保持低频信息的宽频子波振幅谱,从而得到零相位宽频子波。

我们利用谱模拟解析法[17]从提频前地震数据的振幅谱进行高频拓宽,把该宽频振幅谱作为零相位子波振幅谱,采用傅里叶反变换得到真正具有提频前低频信息的宽频子波。

假设地震子波的振幅谱可以表示为下面光滑的解析函数:

(17)

式中:f表示频率;L为多项式阶数;an为依赖于实际地震记录的待定常数。

对于确定的L,将计算出的an(n=0,1,2,…,L)代入(17)式,可以得到一条光滑的拟合曲线|w(f)|,此曲线即是我们所认为的子波振幅谱。

利用公式(17)求出提频前地震数据振幅谱作为子波振幅谱(图1a),然后通过交互扩展该振幅谱高频部分,低频振幅保持不变,得到宽频子波振幅谱(图1b),以该振幅谱与零相位谱重构经过傅里叶反变换后得到保持低频信息的宽频子波(图1c)。

图1　保持低频信息的宽频子波提取a 输入数据振幅谱及拟合曲线; b 拓宽高频拟合曲线得到的振幅谱; c 拓展变频拟合曲线得到的振幅谱与零相位谱重构得到的宽频子波

2　改进型特征值法相干体裂缝预测

断裂、裂缝对碳酸盐岩储层有明显的改善作用,便于解释人员利用预测裂缝手段进行裂缝型储层预测。

相干数据体技术主要利用相邻地震道的相似性原理,描述地层及岩性的横向不均匀性。图2中,地震信号上半部连续性好的同相轴对应于连续性较好的水平地层、倾斜地层等地质体,下半部连续性差的同相轴对应于连续性较差的断层、裂缝等地质体。采用相干数据体算法,对地震数据体的不连续性进行系统分析,使解释人员能有效识别构造和断层的分布情况,避免解释的随意性。

图2　相干体技术描述裂缝构造原理

为了尽可能压制噪声干扰,提高相干体计算数据的分辨率,GERSZTENKORN等[11]提出了第三代特征值相干体计算法。

首先利用MARFURT等[10]在第二代多道相似性相干体计算中定义的椭圆和矩形分析窗口(图3)中,以坐标为(x,y)的分析点为中心,时窗中心t=nΔt,则时间窗口内含2k+1个元素的分析道与窗口内的J道数据排列成矩阵:

(18)

图3　分析窗口示意a 以分析点(红点)为中心的椭圆分析窗口,主轴为a,径向轴为b,主轴方位角为φa; b 以分析点(红点)为中心的矩形窗口,主轴为a,径轴为b,主轴方位角为φa

(19)

设λj为分析窗口内地震数据协方差矩阵公式(19)的特征值,且有λ1≥λ2≥…≥λJ,则在以分析点为中心的分析窗口内,沿视倾角(p,q)的J道C3特征值相干算法公式可表示为:

(20)

则在以分析点为中心的分析窗口内,沿视倾角(p,q=0)的J道C3特征值相干算法公式表示为:

(21)

则沿着所有倾角的C3特征值相干计算算法公式可表示为:

(22)

第三代特征值相干算法C3在子空间中进行计算,比前两代相干算法能更好地压制噪声,但计算量大、效率慢，影响了实际应用时的经济效益。

(23)

3　数值实验及实际资料应用

3.1断层模型谱反演试验

利用一个断层模型(图4a)进行ADM谱反演实验,图4b是对图4a断层模型正演得到30Hz主频的合成记录,图4c为本文方法反演的反射系数,图4d为主频为50Hz重构记录。分析图4可知,ADM谱反演算法能求出极性和位置准确性较高的反射系数,而且通过重构能得到断层分辨更清晰的合成记录。

图4　断层模型反演实验a 断层模型; b 30Hz合成记录; c 谱反演反射系数;d 50Hz重构合成记录

3.2Marmousi模型谱反演试验

利用SEG提供的Marmousi模型进行谱反演反射系数实验,以验证本文方法对纵、横向均存在变化的地震数据反演的适应性。

图5a为Marmousi模型反射系数,图5b是35Hz主频的合成记录,图5c为本文方法反演的反射系数。分析图5可知,ADM谱反演算法能精确求出含断裂带、深部盐丘两侧高速体等复杂数据体的反射系数,证明了该方法的正确性和对复杂模型较好的适应性。

3.3实际资料应用

将基于ADM谱反演重构后的高分辨率实际地震数据用于裂缝预测中,以验证该方法在储层预测中的有效性。对图6a叠前时间域偏移剖面目的层进行ADM谱反演实验,得到反射系数剖面(图6b),该反射系数与测井系数(蓝色曲线)匹配较好,各层位置和极性求取较佳。用提取的宽频子波进行高分辨率重构,得到高分辨率剖面(图6c)，其分辨率明显比图6a的高,且浅、中、深层的同相轴(如白框处)连续性更好、断点更清晰、断层更明了,更有利于后续解释人员进行裂缝预测。分析高分辨率处理前、后的振幅谱(图7)可知,提频后频带得到明显拓宽,而且由于使用保持低频的宽频子波重构,提频后振幅谱的低频信息得到完整保留,为后续反演等解释提供了丰富的频带信息。

为了分析高分辨率处理前、后的效果,分别提取了20ms的C3相干切片(图8),分析图8可知,利用ADM谱反演重构的高分辨率剖面比处理前剖面能获取分辨率更高的相干切片效果,且切片上裂缝走势更清晰、裂缝刻画细节更精细。

利用过井成像测井资料(图9)进一步验证高分辨率处理前、后裂缝发育所在断层的真实性。从图9中可以看出,在2278～2283m段有1个裂缝发育段,在2279m处裂缝发育,通过合成记录精细标定裂缝发育带相当于须三底向下20ms,对最小负曲率属性数据体进行须三底向下20ms切片(图10a),与过井相干剖面(图10b)对比后可知,过井点处发育黑色异常,也即存在裂缝发育带,这也印证了提频后剖面(图6c)箭头所指位置存在明显断层,证明本文方法得到的高分辨率剖面为后续裂缝预测提供了精度更高的依据。

图5　Marmousi模型反演实验a Marmousi模型反射系数; b 35Hz合成记录; c 谱反演反射系数

图6　提频前实际地震数据(a)、谱反演反射系数剖面(b)和宽频子波重构后高分辨率剖面(c)

图7　提频前(a)、后(b)数据振幅谱

图8　提频前(a)、后(b)数据须三底向下20ms相干切片

图9　过井测井剖面

图10　须三底向下20ms最小负曲率属性(a)和过井数据相干剖面(b)

4　结论

本文研究了ADM谱反演高分辨率技术在裂缝预测中应用情况。模型试验和实际资料应用结果表明,利用ADM谱反演算法求取精度高、极性位置正确、能量强弱分明的地层反射系数；与谱模拟解析法提取的宽频子波褶积重构,得到分辨率高、同相轴连续性好、断层清晰、低频信息丰富的剖面；在此剖面上能得到精度更高、细节刻画更清晰的C3相干切片，便于断层和裂缝的有效识别,避免了解释人员在断层解释和组合时的随意性,大大提高了断层解释、裂缝预测的精度。

[1]PORTNIAGUINE O,CASTAGNA J P.Inverse spectral decomposition[J].Expanded Abstracts of 74thAnnual Internat SEG Mtg,2004:1786-1789

[2]PORTNIAGUINE O,CASTAGNA J P.Spectral inversion:lessons from modeling and Boonesville case study[J].Expanded Abstracts of 75thAnnual Internat SEG Mtg,2005:1638-1641

[3]CHOPRA S,CASTAGNA J P,PORTNIAGUINE O.Thin-bed reflectivity inversion[J].Expanded Abstracts of 76thAnnual Internat SEG Mtg,2006:2057-2061

[4]PURYEAR C I,CASTAGNA J P.Layer-thickness determination and stratigraphic interpretation using spectral inversion:theory and application[J].Geophysics,2008,73(2):R37-R48

[5]CHOPRA S,CASTAGNA J,XU Y.Relative acoustic impedance application for thin-bed reflectivity inversion[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009:3554-3558

[6]YUAN S Y,WANG S X,TIAN N.Swarm intelligence optimization and its applicationin in geophysical data inversion[J].Applied Geophysics,2009,6(2):166-174

[7]柴新涛,李振春,韩文功,等.基于LSQR算法的谱反演方法研究[J].石油物探,2012,51(1):11-18

CHAI X T,LI Z C,HAN W G,et al.Spectral inversion method analysis based on the LSQR algorithm[J].Geophysical Prospecting for Petroleum,2012,51(1):11-18

[8]陈祖庆,王静波.基于压缩感知的稀疏脉冲反射系数谱反演方法研究[J].石油物探,2015,54(4):459-466

CHEN Z Q,WANG J B.A spectral inversion method of sparse-spike reflection coefficients based on compressed sensing[J].Geophysical Prospecting for Petroleum,2015,54(4):459-466

[9]BAHORICH M S,FARMER S L.3-D seismic coherency for faults and stratigraphic features:the coherence cube[J].The Leading Edge,1995,14(10):1053-1058

[10]MARFURT K J,KIRLIN R L,FARMER S L,et al.3-D seismic attributes using a semblance-based coherency algorithm[J].Geophysics,1998,63(4):1150-1165

[11]GERSZTENKORN A,MARFURT K J.Eigenstructure-based coherence computations as an aid to 3-D structural and stratigraphic mapping[J].Geophysics,1999,64(5):1468-1479

[12]张军华,王月英,赵勇.C3相干体在断层和裂缝识别中的应用[J].地震学报,2004,26(5):560-564

ZHANG J H,WANG Y Y,ZHAO Y.Application of the third generation of coherent cube in recognizing faults and fractures[J].Acta Seismologica Sinica,2004,26(5):560-564

[13]王玉学,韩大匡,刘文岭,等.相干体技术在火山岩预测中的应用[J].石油物探,2006,45(2):192-196

WANG Y X,HAN D K,LIU W L,et al.The appli-cation of coherence technology in the volcanic reservoir prediction[J].Geophysical Prospecting for Petroleum,2006,45(2):192-196

[14]刘振峰,曲寿利,孙建国,等.地震裂缝预测技术研究进展[J].石油物探,2012,51(2):191-198

LIU Z F,QU S L,SUN J G,et al.Progress of seismic fracture characterization technology[J].Geophysical Prospecting for Petroleum,2012,51(2):191-198

[15]齐晴,曲寿利,李振春,等.也门S2区块叠前、叠后裂缝综合预测方法应用研究[J].石油物探,2014,53 (5):603-608

QI Q,QU S L,LI Z C,et al.Fracture prediction by using pre-stack and post-stack seismic techniques in

S2 block,Yemen[J].Geophysical Prospecting for Petroleum,2014,53(5):603-608

[16]李海山.基于稀疏表示理论的地震信号处理方法研究[D].东营:中国石油大学,2013

LI H S.Seismic data processing methods based on sparse representation theory[J].Dongying:China University of Petroleum,2013

[17]郭廷超,曹文俊,陶长江,等.时变谱模拟反褶积方法研究[J].石油物探,2015,54(1):36-42

GUO T C,CAO W J,TAO C J,et al.Research on time-varying spectral modeling deconvolution method[J].Geophysical Prospecting for Petroleum,2015,54(1):36-42

(编辑：顾石庆)

Research and application of high resolution fracture prediction technology based on ADM spectral inversion

ZHANG Hua1,2,3,4,HE Zhenhua1,2,LI Yalin3,4,ZHANG Xiaobin3,4,HE Guangming3,4,LI Dezhen3,4

(1.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,Chengdu610059,China;2.CollegeofGeophysics,ChengduUniversityofTechnology,Chengdu610059,China;3.GeophysicalExplorationCompany,ChuanqingDrillingEngineeringCo.Ltd.,CNPC,Chengdu610213,China;4.MountainGeophysicalTechnologyTestCenter,CNPC,Chengdu610213,China)

How to obtain the high resolution data volume and high precision coherence slices have been the key research issues in the exploration of fractured reservoirs.Therefore,the application of high resolution fracture prediction technology based on Alternating Direction Method(ADM) spectrum inversion is studied.From the basic principle of spectral inversion,the target function of the multi-layer sparse reflection coefficient is derived.The high precision reflection coefficient is obtained with the ADM spectrum inversion of basis pursuit algorithm.The high resolution section with the low frequency information is reconstructed by the convolution of the high precision reflection coefficient with the broadband zero phase wavelet extracted from spectrum simulation analytical method.The fracture prediction for this high resolution section using the improved eigenvalue algorithm is conducted.The theoretical experiments and actual data application tests show that fracture prediction technology based on ADM spectrum inversion can get the reflection coefficient with the accurate polarity and position as an important technique for the fracture prediction.The reconstructed broadband sections are characterized by higher resolution,better continuous events and more clear breakpoints,and the coherent body slice has higher accuracy with fine reservoir characterization.

resolution,ADM,fracture prediction,coherent body,reflection coefficient

2015-09-25;改回日期：2016-04-12。

张华(1981—),男,博士在读,工程师,主要从事地震资料处理、反演方法研究工作。

中国石油天然气集团公司科学研究与技术开发项目(2013E-38-08)资助。

P631

A

1000-1441(2016)05-0737-09

10.3969/j.issn.1000-1441.2016.05.013

This research is financially supported by the Science Research and Technology Development Project of China National Petroleum Corporation (Grant No.2013E-38-08).