地球物理学报  2018, Vol. 61 Issue (1): 284-292   PDF    
地震反Q滤波应用于碳酸盐岩储层特征描述
周衍, 饶莹     
中国石油大学(北京)油气资源与探测国家重点实验室, 北京 102249
摘要:塔里木盆地塔河地区的奥陶系碳酸盐岩储层埋藏深,其深度超过5000 m,储集体类型以裂缝、缝洞为主,储层的纵向和横向上都具有极强的非均质性.这些特点导致反射地震数据的信噪比低、分辨率低,从而使得碳酸盐岩储层特征描述的难度增大.针对这些问题,本文运用稳定化反Q滤波方法对塔河地区的三维地震数据进行提高分辨率处理,同步实现地震振幅补偿和子波相位校正.处理结果与原始地震数据进行地震剖面面貌、储层内幕信息、断裂缝洞刻画等方面的直接对比,一方面展示反Q滤波方法在提高分辨率、提高弱反射能量、以及增强反射同相轴的空间连续性等方面的优势,同时展示稳定化反Q滤波方法能够提高对奥陶系碳酸盐岩缝洞储集体的特征描述精度.
关键词: Q滤波      储层特征描述      碳酸盐岩储层      缝洞储集体刻画      奥陶系     
Seismic inverse Q filtering for carbonate reservoir characterization
ZHOU Yan, RAO Ying     
China University of Petroleum(Beijing), State Key Laboratory of Petroleum Resources & Prospecting, Beijing 102249, China
Abstract: The Ordovician carbonate reservoir in Tahe oilfield, Tarim basin, is buried very deep with a depth being larger than 5000 m. The reservoir is featured by fracture-vuggies, and the formation has extremely strong heterogeneity in both the vertical and horizontal direction.These features lead to low signal-to-noise ratio and low resolution of the seismic reflection data, and in turn to the difficulty in characterization of carbonate reservoirs.Aiming at these problems, here we use a stabilized inverse-Q filter to process the three-dimensional (3D) seismic data in Tahe oilfield, conducting seismic amplitude compensation and waveform phase correction, simultaneously.Comparing the inverse-Q filtered result directly with the raw seismic data in terms of seismic profiling, interior information within the reservoir, characterization of faults and fracture-vuggies, we demonstrate the advantages of stabilized inverse-Q filtering in enhancing seismic resolution, strengthening weak signals, and improving lateral continunity of reflections, and reveal that this method can improve the accuracy of fracture-vuggy characterization in the Ordovician carbonate reservoirs.
Key words: Inverse Q filter    Reservoir characterization    Carbonate reservoir    Fracture-vuggy illustration    Ordovician    
0 引言

地震波在地层中传播时,由于地层的滤波作用,造成地震波的振幅衰减,不同频率的的平面波其衰减系数不同,而且不同频率的平面波传播的速度也不一样,会造成地震波的相位畸变,即速度频散效应.我们通常将振幅衰减和频散效应统称为地层的Q滤波作用.反Q滤波处理就是对实际的地震数据进行振幅补偿和相位校正,以期消弥地层的Q滤波作用(Wang,2008a).

关于地层的Q滤波作用,最初的数学模型由Kolsky等(Kolsky, 1953, 1956Futterman, 1962)提出:

(1)

式中ω是平面波的角频率,α(ω)是衰减系数,v(ω)是相速度,而Q(ω)则代表了Q对频率的依赖性.在地震分析中,因为地震数据一般记录为实数值,所以ω只有正频率.Wang(2008a)指出(1)式为低衰减的近似表达式,在Q≫1时有效,可应用于像勘探地震这种地层衰减作用比较小的情形.地球物理文献中所谓不同的Q数学模型只在于衰减系数α(ω)和相速度v(ω)的定义上的差别.修正版Kolsky模型由Wang等(2004a)提出.该修正模型可以与其它关于Q的数学模型的衰减系数α(ω)和相速度v(ω)进行直接对比.在地震勘探领域由于其频率较低、频带狭窄,人们可以忽略Q对频率的依赖性而运用所谓的常数Q模型(Kjartansson,1979).但是,从(1)式可以看出,衰减系数α(ω)和相速度v(ω)对频率的依赖性仍然存在.这种频率的依赖性一方面可以被充分利用从而从地震数据中估算Q值,另一方面对反Q滤波处理造成了巨大的困难.

Q滤波可以看作是地震波传播的反过程.Robinson(1979)指出反Q滤波的算法类似于地震反褶积的运算过程.在地震波传播的正过程中,振幅衰减和相位畸变都是频率和传播时间的指数函数.由于表示相位畸变的函数指数是虚数,相位校正的函数指数仍然是虚数,其数学实现时是无条件稳定的.Hargreaves等(1991)提出了通过傅里叶变换进行相位校正的反Q滤波方法,其Q模型是一个Q常数值.Bano(1996)将常数值Q模型的相位反Q滤波方法应用到层状的Q模型上,每层的Q为常数值.

但是,由于表示振幅衰减的函数指数为实数,所对应的振幅补偿的函数指数也是实数,而且是正实数,在数学实现时极其不稳定,将会大大抬升地震数据的噪音,处理之后的地震数据信噪比将急剧下降. Hale(1982)曾建议利用级数展开作近似的高频振幅补偿,达到一定频率时即终止补偿处理.Wang(2002)首次提出稳定化的反Q滤波方法,该方法将数学反问题的求解原理应用到地震处理当中,实现了稳定化的振幅补偿.从此,反Q滤波在地球物理领域发展到了一个崭新时代,可以实现振幅补偿和相位校正的同步处理(Chen et al., 2014; Zhang et al., 2014, 2015; 张固澜等, 2014, 2015),甚至可以实现含反Q滤波的地震偏移成像(Wang and Guo, 2004bWang,2008b; Bai et al., 2016).

稳定化的反Q滤波方法属于“三高”处理技术,即在高保真、高信噪比前提之下的高分辨率处理.Wang(2003)指出反Q滤波处理是否提高地震数据的分辨率取决于频带宽度的变化和信噪比的变化两个因素的综合,并且给出了量化公式.稳定化反Q滤波方法已经成功应用于碎屑岩地层,开展砂体描述及油气预测(郭建等, 2007Gan et al., 2009; 甘其刚等, 2012).但是,本文将稳定化反Q滤波方法应用到碳酸盐岩储层的地震数据上,对碳酸盐岩储层的缝洞储集体进行特征描述.研究区域为塔里木盆地的塔河油田.该地区的奥陶系碳酸盐岩油气藏目标储层埋深大(5000 m以上),储集体类型以裂缝、缝洞为主,储层的纵向、横向非均质强,地震数据信噪比低、分辨率低,使得储层预测难度大.针对塔河地区的这些问题,我们运用稳定化反Q滤波方法对塔河地区的三维地震数据进行处理,提高资料的分辨率,提高弱反射的能量,增强反射同相轴的空间连续性,同时突出地震异常,从而提高对奥陶系碳酸盐岩储层缝洞储集体的预测精度.

1 稳定化反Q滤波方法

稳定化反Q滤波的重要特点即在于该方法的稳定性.正因为如此,地震反Q滤波方法能够同步进行振幅补偿和相位校正(Wang,2008a).

正如前文提到的,平面波的振幅衰减因子为指数函数:

(2)

式中ω是平面波的角频率,因为地震数据一般记录为实数值,所以ω只有正频率,ωh是频率调谐参数,与地震频带的最高有效频率相关(Wang and Guo, 2004a),τ是以秒为单位地震反射波的双程走时,τ=0时则是地震观测的地表,γ(τ)为一无量纲算子:

(3)

振幅补偿因子是振幅衰减因子的逆.Wang(2006)提出的稳定化振幅补偿因子表达式如下:

(4)

其中σ2为稳定化参数,它是抑制数据背景噪声的关键参数,用以保证处理之后地震数据的信噪比得以提高.Wang(2006)明确给出了稳定化参数σ2与反Q滤波增益限G之间的关系公式

(5)

式中增益限G的度量单位为分贝.

前文曾经提到过将数学反问题的求解原理应用到稳定化振幅补偿因子.但是,首先应当注意到(4)式所表述的稳定化过程跟常规的最小二乘求逆方法有不同之处,即稳定化参数σ2不仅像最小二乘方法那样应用到分式的分母上,而且应用到了分子部分.其物理意义可以这样理解:当Q-1趋于0时,根据(2)式的表述指数函数β(τ, ω)应当趋于1,在(4)式中,由于分子和分母同时应用了稳定化参数σ2,据此式得到的广义逆Λ(τ, ω)趋于1,正是“β(τ, ω)趋于1”的逆.Wang(2006)指出稳定化振幅补偿因子的一个重要特点是其可逆性,因为Λ(τ, ω)趋于1,Λ(τ, ω)的逆也趋于1.如果采用常规的最小二乘方法,此处给出的渐近解Λ(τ, ω)将趋于0,这个渐近解在数学上是不可逆的.

结合无条件稳定的相位校正,平面波的反Q滤波基本表达式为

(6)

表达式中的指数函数为相位校正因子,i是虚数单位,U(0, ω)为地表所记录的地震数据(τ=0),是反Q滤波处理的输入数据.对所有平面波结果求和即可得到时域的地震样点:

(7)

该表达式称之为成像条件.

三维地震反Q滤波是一个数据驱动的处理过程,即地层衰减系数模型从三维地震数据体中估算得到.衰减系数Q-1是指地震波在地下介质传播中总能量的耗损,是地层介质的内在物理属性.引起地震波能量衰减的因素包括介质中固体与固体、固体与流体、流体与流体接触界面之间波动接触过程的能量损耗.从地震数据中估算衰减系数的方法请参见Wang(2004).

图 1用显示实际数据的方式展示了三维地震反Q滤波的处理流程,包括从三维数据体中估算Q值,建立一个三维的衰减系数(Q-1)模型,然后依据此模型进行反Q滤波.图 1中的a图是从三维地震数据体中抽出的任意线剖面,b图则为从三维衰减系数模型中抽出的沿该任意线的衰减系数剖面,c图是相对应的反Q滤波处理之后的任意线地震剖面.图中的横向距离是沿任意线的累计距离(km).地震反Q滤波的处理明显提高了地震分辨率.下节我们将仔细分析反Q滤波前后的地震数据及其储层特征细节.

图 1 三维地震数据、衰减系数和反Q滤波数据 (a)从三维地震数据体中抽取的任意剖面;(b)相对应的衰减系数任意剖面;(c)相对应的反Q滤波之后的剖面. Fig. 1 3D seismic data, attenuation coefficients, and inverse Q filtered seismic data (a) Seismic profile arbitrarily extracted from 3D seismic data cube; (b) The attenuation profile corresponding to the arbitrary seismic profile; (c) The corresponding seismic profile after inverse Q filtering.
2 碳酸盐岩储层的实际应用

塔里木盆地是在前震旦纪变质岩基底之上发育的大型复合沉积盆地,由古生代克拉通盆地和中新生代前陆盆地叠合而形成.受区域构造活动和海平面升降的影响,形成三个沉积期:震旦—泥盆纪海相沉积,石炭—二叠纪海陆交互相沉积和三叠—第四纪陆相沉积.塔河油田位于塔里木盆地中部.塔河油田构造处于塔里木盆地沙雅隆起阿克库勒凸起西南斜坡.该区地层层系三叠系以上为陆相碎屑岩沉积,石炭系为海陆交互相碎屑岩夹薄层灰岩,奥陶系为海相碳酸盐岩.其中奥陶系碳酸盐岩经风化剥蚀形成了众多低幅度的风蚀残丘,裂缝及溶蚀孔洞发育,储层的非均层性极强.本文将围绕塔河油田奥陶系目标储层中缝洞储集体的特征描述这一主题,展示稳定化反Q滤波方法在碳酸盐岩储层段的实际三维地震数据上的应用效果.

2.1 地震剖面整体面貌

将稳定化反Q滤波方法对塔河油田的三维地震资料进行补偿处理,从地震剖面的整体面貌来看,在数据的保幅性、信噪比提高和横向纵向分辨率增强等方面优势明显.

该区的目标储层段在图 2中的3480 ms往下的奥陶系碳酸盐岩层.稳定化反Q滤波方法的特点如图 2所示,主要表现在以下几个方面:断层、断点更加清楚;对次级残丘、岩溶沟谷的刻画更加清晰;奥陶系顶风化层反射成像连续性好、信噪比高、能量强,奥陶系内幕反射波能量增强,“串珠状”反射能量聚焦.

图 2 目标储层段的三维地震数据 (a)反Q滤波前的地震剖面; (b)反Q滤波之后的地震剖面.该区奥陶系碳酸盐岩储层的目标层段在时窗3480~3900 ms之间. Fig. 2 3D seismic data with in the target window of receivers (a) Seismic profile before inverse Q filtering; (b) Seismic profile after inverse Q filtering. The target of the Ordovician carbonate reservoir is within a time window between 3480 and 3900 ms.
2.2 “串珠状”反射及不整合面内部细节

塔河油田奥陶系碳酸盐岩储层缝洞储集体的特征指示之一是由缝洞体造成的“串珠状”地震反射.即使与逆时偏移得到的三维地震数据体相对比,我们也可以发现地震数据体经反Q滤波处理之后的空间分辨率明显提高,从图 3中可以识别出更多的“串珠状”反射.

图 3 奥陶系碳酸盐岩层缝洞储集体形成“串珠状”地震反射 (a)从反Q滤波前数据体中抽取的地震剖面; (b)从反Q滤波后数据体中抽取的地震剖面. Fig. 3 Fracture-vuggies in the Ordovician carbonate reservoir generate a series of "string-bead" reflections (a) Seismic profile extracted from the seismic data cube before inverse Q filtering; (b) Seismic profile extracted from the inverse-Q filtered seismic data.

同时注意图 4中椭圆标定的位置,可以看到奥陶系顶层不整合面反射的内部细节得以清晰呈现.图中(左侧时间)3500ms的反射为奥陶系碳酸盐岩储层的顶部位置.

图 4 奥陶系顶层不整合面的内部细节 (a)从反Q滤波前数据体中抽取的地震剖面; (b)从反Q滤波后数据体中抽取的地震剖面. Fig. 4 The unconformity details at the top of the Ordovician carbonate reservoir (a) Seismic profile extracted from the seismic data cube before inverse Q filtering; (b) Seismic profile extracted from the inverse-Q filtered seismic data.
2.3 断裂特征和断裂显示

稳定化反Q滤波处理在对地震振幅补偿的同时也对子波的相位进行了校正.由于这是高保真、高信噪比处理,反Q滤波之后的地震成像效果得以改善,不仅奥陶系顶部的不整合面反射还包括奥陶系中下统内幕的反射界面更加清楚.因为地震绕射波的归位表现合理,进一步提高了三维地震数据的空间横向分辨率,使得断裂的纵向延伸及其断点的位置更加清晰(图 5).

图 5 地震剖面展示的断裂及其断点位置 (a)从反Q滤波前数据体中抽取的地震剖面; (b)从反Q滤波后数据体中抽取的地震剖面. Fig. 5 Fault-fractures and fracturing locations revealed from seismic profiles (a) Seismic profile extracted from the seismic data cube before inverse Q filtering; (b) Seismic profile extracted from the inverse-Q filtered seismic data.

倾角属性可以作为碳酸盐岩地层无明显断距的断裂解释依据.通过倾角属性可以从沿层面识别扰曲、褶皱、凸起、差异压实等特征.这些特征不仅在常规解释时无法追踪,在地震相干属性体上也呈现为连续高相干的特征.在图 6的倾角属性切片上,颜色较深处为大断层,颜色较弱部位为次级别的断层或小裂缝.根据反Q滤波之后的三维地震体计算得到的倾角属性在突出主断裂的同时,也能更好地显示出分支断裂和次级断裂.

图 6 根据三维地震数据计算的倾角属性 (a)根据反Q滤波前地震数据计算的倾角属性; (b)根据反Q滤波后地震数据计算的倾角属性. Fig. 6 Dipping angle attributes derived from 3D seismic data (a) Dipping angle attribute extracted from seismic data before inverse Q filtering; (b) Dippong angle attributes extracted from inverse-Q filtered seismic data.
2.4 裂缝及缝洞特征描述

依据地震振幅的变化我们可以对碳酸盐岩储层的缝洞体进行预测.当碳酸盐岩储层的缝洞发育时并充填油、水或其它物质,缝洞与周围介质存在较大的波阻抗差异,从而形成相对强的地震反射特征.利用振幅变化可以刻画这些局部地震波场的变化,振幅变化大的部位很有可能是裂缝、溶洞的发育带(图 7).

图 7 依据地震振幅变化预测碳酸盐岩缝洞体 (a)利用反Q滤波前的数据进行预测; (b)利用反Q滤波之后的数据进行的预测. Fig. 7 Detection of carbonate fracture-vuggies based on seismic amplitude variation (a) Detection using seismic data cube before inverse-Q filtering; (b) Detection using inverse-Q filtered seismic data cube.

致密的碳酸盐岩与上部低速介质接触时,两层的波阻抗差异较大,会形成较强的地震反射.当碳酸盐岩顶界面发育缝洞系统时,波阻抗相对降低,与上覆介质的波阻抗差异减小形成强背景上的相对弱反射特征.但是,反Q滤波之后的地震数据体上相对弱的反射也得以加强,同时反Q滤波处理突出了地震振幅异常,“串珠状”反射显示更加明显,振幅异常的边界刻画明确.

图 8展示的是预测的缝洞体平面分布与地震剖面上“串珠状”反射的对应关系.缝洞体平面图上的红色曲线是任意选取的一条剖面.结合地震剖面对缝洞体进行分析,使得碳酸盐岩储层油气有利缝洞带的预测更加精确.

图 8 预测的缝洞体平面分布与地震剖面上“串珠状”反射的对应关系 (a)反Q滤波前的缝洞体预测及其地震剖面; (b)反Q滤波之后的缝洞体预测及其地震剖面. Fig. 8 Spatial distribution of detected fracture-vuggies and the corresponding "string-bead" reflections in seismic profiles (a) Fracture-vuggy detection and the seismic profile before inverse-Q filtering; (b) Fracture-vuggy detection and the seismic profile after inverse-Q filtering.

由于奥陶系碳酸盐岩储层极强的非均质性,除了存在有大断裂之外,发育更多的是微裂缝.常规的地震相干属性可以预测大断裂但是不能精细刻画微裂缝.我们采用基于数据本征结构分析的方法计算三维地震数据的相干属性,经过线性增强等一系列处理,可以直观地反映裂缝的发育程度和发育特征,如图 9中的蓝色显示.从不同角度来看三维缝洞体雕刻,都可以清晰看到利用反Q滤波之后的地震数据体雕刻出来的裂缝更加精细.因为反Q滤波处理之后的地震数据体中“串珠状”反射能量更强,利用该数据体进行雕刻的缝洞体(图 9中的黄色显示)直观数量更多,同时裂缝与“串珠状”反射的对应关系也更好.

图 9 从奥陶系顶部俯视的裂缝(蓝色)和溶洞(黄色)的三维雕刻体 (a)反Q滤波之前地震数据的缝洞雕刻; (b)反Q滤波处理之后地震数据的缝洞雕刻. Fig. 9 Three-dimrensional illustration of faults (in blue) and fracture-vuggies (in yellow), looking down from the top of the Ordovician formation (a) Illustration of faults and fracture-vuggies using seismic data before inverse-Q filtering; (b) Illustration of faults and fracture-vuggies using inverse-Q filtered seismic data.
3 结论

针对奥陶系碳酸盐岩缝洞型储层埋藏深、非均质性强、分辨率低、信噪比低等问题,采用稳定化反Q滤波方法对塔河油田三维地震数据进行处理.处理结果与处理之前的数据相比表明,稳定化反Q滤波方法能够提高地震分辨率,清晰显示储层构造特征,突出地震“串珠状”反射.利用反Q滤波处理之后的三维地震数据进行储层特征描述,奥陶系碳酸盐岩储层的层间小断裂以及缝洞储集体成像清晰准确.

参考文献
Bai M, Chen X H, Wu J, et al. 2016. Q-compensated migration by Gaussian beam summation method. Journal of Geophysics and Engineering, 13(1): 35-48. DOI:10.1088/1742-2132/13/1/35
Bano M. 1996. Q-phase compensation of seismic records in the frequency domain. Bulletin of the Seismological Society of America, 86(4): 1179-1186.
Chen Z B, Chen X H, Wang Y H, et al. 2014. Estimation of Q factors from reflection seismic data for a band-limited and stabilized inverse Q filter driven by an average-Q model. Journal of Applied Geophysics, 101: 86-94. DOI:10.1016/j.jappgeo.2013.12.003
Futterman W I. 1962. Dispersive body waves. Journal of Geophysical Research, 67(13): 5279-5291. DOI:10.1029/JZ067i013p05279
Gan Q Q, Xu D, Tang J M, et al. 2009. Seismic resolution enhancement for tight-sand gas reservoir characterization. Journal of Geophysics and Engineering, 6(1): 21-28. DOI:10.1088/1742-2132/6/1/002
Hale D. 1981. An inverse-Q filter. Standford Exploration Project Report, 26: 231-243.
Hargreaves D, Calvert A J. 1991. Inverse Q filtering by Fourier transform. Geophysics, 56(4): 519-527. DOI:10.1190/1.1443067
Kjartansson E. 1979. Constant Q wave propagation and attenuation. Journal of Geophysical Research:Solid Earth, 84(B9): 4737-4748. DOI:10.1029/JB084iB09p04737
Kolsky H. 1953. Stress Waves in Solids. Oxford: Clarendon Press.
Kolsky H. 1956. The propagation of stress pulses in viscoelastic solids. Philosophical Magazine, 1(8): 693-710. DOI:10.1080/14786435608238144
Robinson J C. 1979. A technique for the continuous representation of dispersion in seismic data. Geophysics, 44(8): 1345-1351. DOI:10.1190/1.1441011
Wang Y H. 2002. A stable and efficient approach of inverse Q filtering. Geophysics, 67(2): 657-663. DOI:10.1190/1.1468627
Wang Y H. 2003. Quantifying the effectiveness of stabilized inverse Q filtering. Geophysics, 68(1): 337-345. DOI:10.1190/1.1543219
Wang Y H. 2004. Q analysis on reflection seismic data. Geophysical Research Letters, 31(17): L17606.
Wang Y H, Guo J. 2004a. Modified Kolsky model for seismic attenuation and dispersion. Journal of Geophysics and Engineering, 1(1): 187-196.
Wang Y H, Guo J. 2004b. Seismic migration with inverse Q filtering. Geophysical Research Letters, 31(21): L21608.
Wang Y H. 2006. Inverse Q-filter for seismic resolution enhancement. Geophysics, 71(3): V51-V60. DOI:10.1190/1.2192912
Wang Y. 2008a. Seismic Inverse Q filtering. Oxford: Blaclwell.
Wang Y H. 2008b. Inverse-Q filtered migration. Geophysics, 73(1): S1-S6. DOI:10.1190/1.2806924
Zhang G L, Wang X M, He Z H, et al. 2014. Impact of Q value and gain-limit to the resolution of inverse Q filtering. Journal of Geophysics and Engineering, 11(4): 045011. DOI:10.1088/1742-2132/11/4/045011
Zhang G L, Wang X M, He Z H. 2015. A stable and self-adaptive approach for inverse Q-filter. Journal of Applied Geophysics, 116: 236-246. DOI:10.1016/j.jappgeo.2015.03.012
Zhang G L, He Z H, Wang X M, et al. 2014. Seismic wave dispersion effects and inverse Q-filter phase compensation. Chinese Journal of Geophysics, 57(5): 1655-1663. DOI:10.6038/cjg20140528
Zhang G L, Lin J, Wang X M, et al. 2015. A self-adaptive approach for inverse Q-filtering. Chinese Journal of Geophysics, 58(7): 2525-2535. DOI:10.6038/cjg20150727
甘其刚, 许多, 王仰华. 2012. 利用稳定反Q滤波技术提高储层预测精度-以川西新场气田上三叠统须家河组气藏为例. 天然气工业, 32(3): 42–44.
郭建, 王成彬, 胡中平, 等. 2007. Q补偿技术在提高地震分辨率中的应用-以准噶尔盆地Y1井区为例. 石油物探, 46(5): 509–513, 530.
张固澜, 贺振华, 王熙明, 等. 2014. 地震波频散效应与反Q滤波相位补偿. 地球物理学报, 57(5): 1655–1663, 125. DOI:10.6038/cjg20140528
张固澜, 林进, 王熙明, 等. 2015. 一种自适应增益限的反Q滤波. 地球物理学报, 58(7): 2525–2535. DOI:10.6038/cjg20150727