2. 西安交通大学电子与信息工程学院, 陕西西安 710049
2. School of Electronic and Information Engineering, Xi'an Jiaotong University, Xi'an 710049, China
如何获得高分辨率的地震资料始终是地震勘探领域的重点研究内容之一。除了改进采集方式获得高精度宽频资料外, 对所获得的资料进行高分辨率处理也是一条重要途径。近年来, 国内外学者提出了多种提高分辨率的处理方法, 如反褶积[1-2]、谱白化[3-4]、宽带约束反演[5]、反Q滤波[6-8]、井地联合处理[9-13]、叠前分角度处理[14]等, 均取得了明显的应用效果。反褶积[1-2]类方法通过对地震记录中的地震子波进行压缩达到提高分辨率的目的, 是目前提高叠后地震资料分辨率重要且有效的手段。目前该类方法已经发展了很多实用算法, 但为了实现稳定的反褶积, 需作子波最小相位假设。然而, 由于实际资料中的地震子波多为混合相位, 该方法的应用效果受到限制, 同时, 为了减少噪声对反褶积算子的影响, 发展了时频域反褶积方法, 并取得一定的进展[15]。
工业界常用的谱白化方法[3-4]具有操作简单、鲁棒性强等优点, 缺点是保真性较差。宽带约束反演方法[5]是将地震资料、测井信息和先验地质知识有机结合起来的叠后迭代反演方法, 可以大幅度提高地震剖面的分辨率, 尤其是在井资料比较密集的情况下, 效果明显优于反褶积方法。但是此类方法强烈依赖于初始模型, 反演结果往往具有多解性。反Q滤波方法[6-8]通过补偿地震波在地下介质中传播时产生的衰减效应来提高地震信号分辨率, 其应用效果直接取决于地层品质因子Q值模型和Q值估计精度, 如何较为准确地估计Q值是该方法需要解决的难题。CHOPRA等[9]利用井中不同深度VSP资料获得地震波的幅度和频率衰减特性, 成功地恢复了地震资料中被衰减的高频分量。程亮等[13]给出了一种井地联合的提高地震资料分辨率方法, 并利用实际资料验证了方法的有效性。张明振[14]根据叠前道集小角度数据及大角度数据同源的机理, 提出一种压制大角度数据调谐效应, 进而提高叠加数据分辨率的方法。
利用小波变换提高地震资料分辨率也一直是业界的研究热点。吴爱娣等[16]应用二维连续小波变换提高地震资料分辨率, 方法是利用信号及噪声在时间、频率、空间及波数域的差异来提高地震资料信噪比, 同时通过调整各尺度信号能量大小提高地震资料分辨率。陈文超等[17]在连续小波域大幅度提高了过井地震剖面的分辨率, 指出井震资料匹配程度限制了该方法的应用效果。王西文等[18]将井旁地震道及合成地震记录在小波域分频展开, 提出了一种参考合成地震记录优化实际地震数据小波重构系数的方法, 该方法能明显地提高地震资料分辨率。此外, 王西文等[19]构造了能提高地震资料分辨率、提取瞬时属性的导数小波函数方法, 取得了明显的应用效果。邱庚香等[20]对地震记录进行二进小波多分辨率分解, 通过补偿高频分量来提高地震记录分辨率, 然而, 该方法的应用效果受制于原始资料信噪比。SMITH等[21-22]利用连续小波变换拓展地震资料的有效带宽, 明显提高了地震资料的分辨率。该方法不仅能对信号的高频段进行补偿, 而且能对信号的低频段进行补偿, 在提高地震资料分辨率的同时, 使地震信号保真度明显高于常用的谱白化等方法。陈文超等[23]进一步发展了连续小波域提高地震资料分辨率的方法, 利用时变信号在时间-尺度域的多分辨特性, 基于宽频带高分辨率地震信号振幅谱能量近似水平的假设, 在连续小波变换域拓展地震资料的有效带宽。
本文研究了连续小波域自适应地震资料频谱拓展方法, 给出了方法的具体实现步骤及主要参数选择策略, 利用准噶尔盆地高密度采集三维叠前时间偏移数据验证了方法的有效性。
1 方法原理地震波在地下介质中传播时, 高频信号相对低频信号吸收衰减更快, 形成了地震信号主频低、带宽窄的特点。若将地震信号分解到连续小波域, 则表现为高频信号的小尺度小波系数能量相对较弱。基于宽频带高分辨率地震信号振幅谱能量近似水平的假设, 在连续小波域补偿传播损失的高频能量, 可以恢复地震资料的高分辨率特征。
首先, 对单道地震信号做连续小波变换。将单道地震数据分解到时间-尺度域, 即获得拓频前地震信号的小波变换系数:
$ W(\tau, s)=\int\limits_{-\infty}^{\infty} f(t) \frac{1}{\sqrt{s}} \varPsi^{*}\left(\frac{t-\tau}{s}\right) \mathrm{d} t $ | (1) |
式中, f(t)为单道地震信号, Ψ(t)为母小波, τ和s为平移因子及尺度因子。连续小波变换在时频平面(即时间-尺度平面)内具有良好的局部化特性和冗余特性, 为频率域拓展频带提供了可能性。连续小波变换的时频分辨率受不确定性原理的约束, 同时与尺度相关:尺度增大时, 时间分辨率较低, 频率分辨率较高; 反之, 尺度减小时, 时间分辨率较高, 频率分辨率较低。
其次, 在原始地震信号的振幅谱中选择一个基准频率, 计算所需要拓展频率的范围。由于信号的有效频宽按照倍频变化时, 时间域观察到的信号分辨率变化最明显, 因此, 频谱拓展范围的计算采用倍频程的方式。图 1为拓展频率范围选择示意图, A表示所选择的高频端基准频率, 即原始地震数据有效频宽的高截频对应的频率值, B为A/2, B到A之间的频率段为与拓展高频段对应的基础频率。下文例子中, 基准频率A(即高频基准频率)约为20 Hz, 则B约为10 Hz。拓展的高频分量的频率称为谐波(一次谐波、二次谐波等), 基础频率段为B(即A/2)到A之间, 一次谐波频率段为A到2A之间, 二次谐波为2A到4A之间, 以此类推。利用基础频率段的振幅谱预测谐波的振幅谱, 即通过计算基础频率段的能量密度校正谐波的能量密度, 使得谐波能量密度与基础频率段能量密度基本一致, 从而达到拓展带宽的目的。低频端拓展频率范围选择与高频端类似, 图 1虚线即为拓展低频端的基准频率(假设为C)。下文例子中, 低频基准频率约为5 Hz(即低截频对应的频率值)。拓展的低频分量的频率称为次谐波(一次次谐波等), 其中一次次谐波频率段在C/2到C之间。低频端(次谐波)频谱拓展与高频端(谐波)频谱拓展类似。以上校正过程均需在连续小波域完成。
最后, 在时频域(即时间-尺度域)得到已拓展带宽的信息, 方法是:将谐波、次谐波频率段对应的时间域尺度下的小波系数W(τ, s)乘以一个相应的权值因子, 得到加权后的小波系数
$ \tilde{f}(t)=\frac{1}{C_{\varPsi}} \int\limits_{0-\infty}^{\infty} \int\limits_{-\infty}^{\infty} \frac{1}{\sqrt{s}} \widetilde{W}(\tau, s) \varPsi\left(\frac{t-\tau}{s}\right) \frac{\mathrm{d} s \mathrm{d} \tau}{s^{2}} $ | (2) |
式中:CΨ为满足可允许性条件的常数, 其数值由母小波决定。另外, 本文采用了Morlet小波作为小波母函数。
2 实际资料应用 2.1 沙漠区M6三维工区M6三维工区位于准噶尔盆地腹部沙漠地区。勘探目的层包括侏罗系(J)、三叠系(T)、二叠系(P)、石炭系(C), 以二叠系、石炭系为主。地表沙丘对地震波有很强的吸收衰减作用, 严重降低了地震分辨率; 同时, 目的层埋深大, 进一步增加了目的层识别的难度。该工区的地震资料为可控震源高密度高效采集资料, 采用宽频带震源激发(扫描频率为1.5~84.0 Hz)、宽方位接收(主要目的层段观测系统横纵比约为0.78)、高覆盖次数(大约1 000次)以及基于表层精细调查的吸收补偿等。在数据处理阶段, 采用分类、分阶段、多域等多种手段压制各种噪声[24-26], 提升了地震资料的信噪比。但该工区叠前时间偏移数据存在低频信息丰富、主频低、高频弱的问题, 影响了地质目标的识别。为此, 采用本文方法进行了提高分辨率处理。
图 2对比了过M6井线地震剖面应用本文方法处理前后的结果, 图 2a为处理前的叠前时间偏移剖面, 图 2b为本文方法提高分辨率处理后的结果。可以看出, 应用连续小波变换自适应拓频方法后, 剖面时间方向分辨率得到了明显提升, 且保持了较高的信噪比。时间分辨率的提高使白垩系(K)、侏罗系(J)、三叠系(T)、二叠系(P)、石炭系(C)地震反射接触关系更加清楚, 断裂特征改善, 断点清晰。尤其是在大幅度提高侏罗系强反射同相轴分辨率的同时, 改善了石炭系顶界与二叠系不整合接触关系的成像, 为后续高精度属性分析提供了良好的数据, 也为石炭系构造研究提供了丰富的信息。
图 3是图 2对应数据的多道归一化振幅谱, 其中红线为应用本方法前的归一化振幅谱, 绿线为本文方法拓频后数据的归一化振幅谱。可以看出, 本文方法明显拓展了地震记录的有效频宽(包括高频段及低频段), 尤其是对低频段有效能量具有良好的保持能力, 为后续与反演相关的处理环节提供了有力保障。
以下参照原始数据的多道归一化振幅谱具体说明本文方法的参数选择问题。如图 3所示, 原始资料高截频约20 Hz, 即基准频率A=20 Hz。经过分析发现, 原始数据有较高的信噪比, 且在60 Hz频段还存在有效信息, 因此可将原始数据高频截止频率由20 Hz扩展至80 Hz。应用本文方法, 即可扩展一次谐波段(A至2A)及二次谐波段(2A至4A)的高频能量。原始数据低频能量的拓展与高频能量拓展类似。在利用本文方法进行实际资料处理时, 可根据具体情况对频谱拓展后的数据进行低通滤波, 以保持数据的信噪比。例如, 虽然本文方法可将高频端拓展到80 Hz, 但考虑到60 Hz以上基本没有有效信息, 因此可利用高截频为60 Hz的滤波器进行低通滤波。由图 3可以看出, 如果低频按2 Hz计算, 则提高分辨率处理后数据频宽约为5个倍频程。也就是说, 该区域利用低频震源勘探得到了宽频成果数据(5个倍频程), 这与目前工业界关于宽频带高分辨率地震数据采集、处理的认识是一致的。
2.2 玛湖地区MZ三维工区MZ三维工区位于准噶尔盆地玛湖西斜坡中段, 为玛湖地区高密度三维高分辨率勘探攻关的重点区域, 主要地质目标是三叠系和二叠系。对该目标的勘探主要存在以下难题:①地震资料主频低、频带窄、整体分辨率低; ②储层薄导致岩性目标边界难以落实; ③断裂不清晰影响圈闭识别。为获得高品质宽频原始地震资料, 该工区三维地震采用井炮采集, 并开展了精细的近地表激发岩性与深度调查。然而, 经攻关处理得到的叠前时间偏移数据分辨率还是难以满足地质目标研究的需求。
为此, 我们抽取过MZ4井主测线的叠前时间偏移数据进行了小波域提高分辨率处理试验。图 4对比了过MZ4井主测线三维地震剖面应用本文方法处理前后的结果, 可以看出, 本文方法在明显提高三叠系和二叠系等主要目的层反射同相轴分辨率的同时, 保持了较高的信噪比。与本文方法处理前的叠前时间偏移结果相比, 处理后的结果具有三个显著特点:①三叠系及以上地层分辨率明显提高; ②改善了断裂成像质量, 断点更加清晰干脆; ③三叠系白碱滩组(T3b)、百口泉组(T1b)油层所在砂层组的空间展布变得清晰, 有利于识别圈闭边界。
图 5为本文方法应用前后过MZ4井线归一化平均振幅谱。可以看出, 经本文方法处理后的结果频宽拓展明显, 35~60 Hz范围的信号能量得到提高, 主频从20 Hz提高到35 Hz左右。同时, 本文方法处理后较好地保持了低频信号的能量, 说明本文方法具有显著的保持振幅能力。
本文研究了连续小波域提高分辨率处理方法, 对准噶尔盆地高密度宽频采集的三维叠前时间偏移数据进行了频带拓展, 同时拓展了原始数据的高、低频有效信号, 取得以下成果与认识:
1) 沙漠区M6工区的实际资料拓频处理使白垩系、侏罗系、三叠系、二叠系、石炭系地震反射接触关系更加清楚, 断裂特征改善, 断点清晰, 尤其是在大幅度提高侏罗系强反射同相轴分辨率的同时, 改善了石炭系顶界与二叠系不整合接触关系成像。
2) 玛湖地区MZ工区实际资料拓频处理较好地提高了侏罗系至三叠系地层分辨率, 改善了断裂成像, 使断点更加清晰干脆, 三叠系白碱滩组、百口泉组油层所在砂层组的空间展布变得清晰, 有利于识别圈闭边界。
3) 本文方法对于宽频采集的地震资料有良好的适应能力, 采用基于低频震源的陆上宽频采集技术并结合高分辨率处理方法, 获得宽频地震资料切实可行。
4) 目前常规勘探技术用于沙漠、黄土源探区或者深层目标很难进一步拓展高频信息, 建议采用宽频采集技术, 即尽量提高低频有效信息能量, 获得宽频成果资料。另外, 将叠后拓频方法推广至叠前, 配合保真去噪、高精度速度估计及各向异性校正等, 也可进一步拓展地震资料的高频有效信息。
[1] |
YILMAZ O. Seismic data processing[M]. USA: Society of Exploration Geophysicists, 1987: 82-155.
|
[2] |
黄绪德. 反褶积与地震道反演[M]. 北京: 石油工业出版社, 1992: 11-45. HUANG X D. Deconvolution and seismic trace inversion[M]. Beijing: Petroleum Chemical Press, 1992: 11-45. |
[3] |
边国柱, 张立群. 地震数据的谱白化处理[J]. 石油物探, 1986, 25(2): 32-39. BIAN G Z, ZHANG L Q. Spectral whitening of seismic data[J]. Geophysical Prospecting for Petroleum, 1986, 25(2): 32-39. |
[4] |
颜中辉, 方刚, 徐华宁, 等. 希尔伯特谱白化方法在海洋地震资料高分辨率处理中的应用[J]. 海洋地质与第四纪地质, 2018, 38(4): 212-220. YAN Z H, FANG G, XU H N, et al. The application of Hilbert spectral whitening method to high resolution processing of marine seismic data[J]. Marine Geology & Quaternary Geology, 2018, 38(4): 212-220. |
[5] |
彭传平.宽带约束反演方法实现及应用研究[D].西安: 长安大学, 2008 PENG C P.Implementation and application of broadband constrained inversion method[D].Xi'an: Chang'an University, 2008 http://cdmd.cnki.com.cn/Article/CDMD-11941-2009066663.htm |
[6] |
WANG Y H. A stable and efficient approach of inverse Q filtering[J]. Geophysics, 2002, 67(2): 657-663. DOI:10.1190/1.1468627 |
[7] |
WANG Y H. Inverse-filter for seismic resolution enhancement[J]. Geophysics, 2006, 71(3): V51-V60. DOI:10.1190/1.2192912 |
[8] |
CARDIMONA S. Waveform inversion and digital filter theory[J]. Geophysics, 1991, 56(4): 534-536. DOI:10.1190/1.1443069 |
[9] |
CHOPRA S, ALEXEEV V, SUDHAKAR V, et al. High-frequency restoration of surface seismic data[J]. The Leading Edge, 2003, 22(8): 730-738. DOI:10.1190/1.1605071 |
[10] |
孔庆丰. 利用井间地震提高地面地震资料分辨率[J]. 物探与化探, 2010, 34(3): 353-357. KONG Q F. The study and application of crosswell seismic technology to improving surface seismic resolution[J]. Gophysical and Geochemican Exploration, 2010, 34(3): 353-357. |
[11] |
严红勇, 刘洋, 于茜. 利用VSP提高地面地震资料的分辨率[J]. 石油天然气学报, 2010, 32(3): 65-68. YAN H Y, LIU Y, YU Q. Improving the resolution of poststack surface seismic data by VSP[J]. Journal of Oil and Gas Technology, 2010, 32(3): 65-68. DOI:10.3969/j.issn.1000-9752.2010.03.014 |
[12] |
张涛, 林承焰, 张宪国. 利用测井曲线反褶积方法提高薄层识别能力[J]. 石油勘探与开发, 2010, 37(5): 579-582. ZHANG T, LIN C X, ZHANG X G. Improve thin beds discrimination using the log curves deconvolution method[J]. Petroleum Exploration and Development, 2010, 37(5): 579-582. |
[13] |
程亮, 孔庆丰, 宋建国. 步进迭代法井地联合地震资料拓频处理[J]. 石油地球物理勘探, 2017, 52(4): 696-703. CHENG L, KONG Q F, SONG J G. The frequency spectrum extension for borehole joint seismic data by marching iteration[J]. Oil Geophysical Prospecting, 2017, 52(4): 696-703. |
[14] |
张明振. 基于地震叠前大角度道集的调谐压制提高分辨率方法[J]. 科学技术与工程, 2017, 17(11): 169-174. ZHANG M Z. The method of resolution improvement through tuning suppression based on large angle pre-stack seismic data[J]. Science Technology and Engineering, 2017, 17(11): 169-174. DOI:10.3969/j.issn.1671-1815.2017.11.027 |
[15] |
高静怀, 汪玲玲, 赵伟. 基于反射地震记录变子波模型提高地震记录分辨率[J]. 地球物理学报, 2009, 52(5): 1289-1300. GAO J H, WANG L L, ZHAO W. Enhancing resolution of seismic traces based on the changing wavelet model of the seismogram[J]. Chinese Journal of Geophysics, 2009, 52(5): 1289-1300. DOI:10.3969/j.issn.0001-5733.2009.05.018 |
[16] |
吴爱娣, 牟永光. 二维小波变换在提高地震资料分辨率中的应用[J]. 石油物探, 1997, 36(1): 36-42. WU A D, MOU Y G. The use of 2-D wavelet transform in improving the resolution of seismic data[J]. Geophysical Prospecting for Petroleum, 1997, 36(1): 36-42. |
[17] |
陈文超, 高静怀. 小波变换用于提高过井地震资料分辨率的方法研究[J]. 西安交通大学学报, 1998, 32(12): 44-47. CHEN W C, GAO J H. Enhancement of the resolution of the seismogram across a well using wavelet transform technique[J]. Journal of Xi'an Jiaotong University, 1998, 32(12): 44-47. |
[18] |
王西文, 高静怀, 李幼铭. 高分辨地震资料处理中导数小波函数的构造[J]. 石油物探, 2000, 39(2): 64-71. WANG X W, GAO J H, LI Y M. Construction of derivative wavelet function in high resolution data processing[J]. Geophysical Prospecting for Petroleum, 2000, 39(2): 64-71. DOI:10.3969/j.issn.1000-1441.2000.02.007 |
[19] |
王西文, 刘全新, 高静怀, 等. 地震资料在小波域的分频处理与重构[J]. 石油地球物理勘探, 2001, 36(1): 78-85. WANG X W, LIU Q X, GAO J H, et al. Frequency-shared processing and reconstruction of seismic data in wavelet domain[J]. Oil Geophysical Prospecting, 2001, 36(1): 78-85. DOI:10.3321/j.issn:1000-7210.2001.01.012 |
[20] |
邱庚香, 熊彬, 梁锦文. 应用小波变换提高浅层地震记录分辨率[J]. 桂林理工大学学报, 2003, 23(2): 159-163. QIU G X, XIONG B, LIANG J W. Application of wavelet transformation in resolution of shallow seismic record[J]. Journal of Guilin Institute of Technology, 2003, 23(2): 159-163. DOI:10.3969/j.issn.1674-9057.2003.02.004 |
[21] |
SMITH M, GARY P, JAIME S, et al. Extending seismic bandwidth using the continuous wavelet transform[J]. First Break, 2008, 26: 97-102. |
[22] |
SMITH M.Bandwidth extension using harmonics[R].Salvador: Eleventh International Congress of Geophysical Society, 2009
|
[23] |
陈文超, 王伟, 李瑞萍.连续小波域地震信号自适应频谱拓展方法[C]//中国石油学会2015年物探技术研讨会论文集.宜昌: 中国石油学会物探专业委员会, 2015: 649-652 CHEN W C, WANG W, LI R P.The method to expand the frequency spectrum bandwidth adaptively for seismic data in continuous wavelet domain[C]//Proceedings of the 2015 Symposium on Geophysical Exploration Technology of the Chinese Petroleum Association.Yichang: Society of Petroelum Geophysicists, 2015: 649-652 http://cpfd.cnki.com.cn/Article/CPFDTOTAL-SYDQ201505001161.htm |
[24] |
王西文, 赵邦六, 吕焕通, 等. 地震资料相对保真处理方法研究[J]. 石油物探, 2009, 48(4): 319-331. WANG X W, ZHAO B L, LV H T, et al. Relative amplitude preservation method for seismic data processing[J]. Geophysical Prospecting for Petroleum, 2009, 48(4): 319-331. DOI:10.3969/j.issn.1000-1441.2009.04.001 |
[25] |
王伟, 姜呈馥, 张丽霞, 等. 一种边界保持随机噪声衰减的各向异性拉普拉斯滤波技术[J]. 中国煤炭地质, 2015, 27(9): 61-66. WANG W, JIANG C F, ZHANG L X, et al. An anisotropic laplace filtering technology with edge random noise attenuation retaining[J]. Coal Geology of China, 2015, 27(9): 61-66. DOI:10.3969/j.issn.1674-1803.2015.09.15 |
[26] |
CHEN J Y, CHEN W C, WANG X K, et al.A modified Cadzow filtering algorithm based on anti-picket fence effect[C]//2017 CGS/SEG Annual Meeting.Qindao: China University of Petroleum (East China), 2017: 17-20
|