留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于HHT方法的地层Q值补偿方法研究

颜中辉 徐华宁 李攀峰 刘鸿 杨佳佳 陈珊珊 杨传胜

颜中辉, 徐华宁, 李攀峰, 刘鸿, 杨佳佳, 陈珊珊, 杨传胜. 基于HHT方法的地层Q值补偿方法研究[J]. 海洋地质与第四纪地质, 2020, 40(1): 175-183. doi: 10.16562/j.cnki.0256-1492.2018100901
引用本文: 颜中辉, 徐华宁, 李攀峰, 刘鸿, 杨佳佳, 陈珊珊, 杨传胜. 基于HHT方法的地层Q值补偿方法研究[J]. 海洋地质与第四纪地质, 2020, 40(1): 175-183. doi: 10.16562/j.cnki.0256-1492.2018100901
YAN Zhonghui, XU Huaning, LI Panfeng, LIU Hong, YANG Jiajia, CHEN Shanshan, YANG Chuangsheng. Research of compensation method for formation Q value based on HHT method[J]. Marine Geology & Quaternary Geology, 2020, 40(1): 175-183. doi: 10.16562/j.cnki.0256-1492.2018100901
Citation: YAN Zhonghui, XU Huaning, LI Panfeng, LIU Hong, YANG Jiajia, CHEN Shanshan, YANG Chuangsheng. Research of compensation method for formation Q value based on HHT method[J]. Marine Geology & Quaternary Geology, 2020, 40(1): 175-183. doi: 10.16562/j.cnki.0256-1492.2018100901

基于HHT方法的地层Q值补偿方法研究


doi: 10.16562/j.cnki.0256-1492.2018100901
详细信息
    作者简介:

    颜中辉(1986—),男,硕士,助理研究员,主要从事地震资料处理、多波多分量研究,E-mail:zhonghui4564573@163.com

    通讯作者: 徐华宁,男,博士,教授级高工,主要从事海洋地球物理调查研究,E-mail:xuhuaning1973@163.com
  • 基金项目:  山东省自然科学基金培养基金项目“基于HHT变换的多波地震资料去噪应用研究”(ZR2015PD006);山东省自然基金博士基金 “地质约束下的重磁三维反演在南黄海盆地崂山隆起的应用”(ZR2016DB33);青岛海洋科学与技术国家实验室开放基金“多分量地震波逆时偏移拉伸校正及其并行”(QNLM2016ORP0206);国家自然科学基金青年基金“东海北部外陆架晚更新世以来两期古三角洲的时空展布特征及物源属性研究”(41706069)
  • 中图分类号: P631.4

Research of compensation method for formation Q value based on HHT method

More Information
  • 摘要: 常规的基于叠后地震资料的地层Q值补偿方法受数据平均效应的影响,Q值估算精度较低,导致保真性差,而HHT时频分析可以反映出信号的局部信息,基于以上原因,提出一种基于HHT方法的地层Q值补偿方法。首先应用HHT分频技术将浅层、中深层地震数据分成不同的频率范围,得到不同目标层频率尺度的地震数据;然后在HHT时频域内通过谱比法提取Q值;最后分别进行地层Q值振幅补偿和相位补偿。该方法针对不同的目标地层将地震数据分解成不同的频率尺度,避免了不同频率成分之间的干扰,而且可针对性处理目标层的频率信息,提高了地层Q值提取的稳定性。模型试算和实际资料应用表明,该方法能够有效地对振幅和相位进行补偿,使中深层反射波形得到恢复,主频提高,频带展宽,地震资料的品质得到提升。
  • 图  1  模型数据Q补偿前后的对比

    a. 补偿前的地震数据,b. 补偿后的地震数据。

    Figure  1.  Comparison of model seismic data before and after Q compensation

    a. Model data before compensation,b. Model data after compensation.

    图  2  模型数据Q补偿前后的波形显示对比

    a. 补偿前的波形显示,b. 补偿后的波形显示。

    Figure  2.  Waveform comparison of model seismic data before and after Q compensation

    a. Waveform before compensation,b. Waveform after compensation.

    图  3  基于HHT方法提取的地层Q值曲线

    Figure  3.  Formation Q-value curve extracted based on HHT method

    图  4  模型数据Q补偿前后的频谱对比

    Figure  4.  Spectrum comparison before and after model data Q compensation

    图  5  加噪的模型数据Q补偿前后的对比

    a. 补偿前的地震数据,b. 补偿后的地震数据。

    Figure  5.  Comparison of model seismic data with noise before and after Q compensation

    a. Model data before compensation,b. Model data after compensation.

    图  6  加噪的模型数据Q补偿前后的频谱对比

    Figure  6.  Spectrum comparison of model seismic data with noise before and after Q compensation

    图  7  基于HHT方法提取的地层Q值剖面

    Figure  7.  Formation Q-value profile extracted with HHT method

    图  8  实际地震数据Q补偿结果对比

    a. 补偿前的地震数据,b. 补偿后的地震数据。

    Figure  8.  Comparison of actual seismic data before and after Q compensation

    a. Seismic data before compensation, b. Seismic data after compensation.

    图  9  浅层数据Q补偿结果对比

    a. 补偿前的地震数据,b. 补偿后的地震数据。

    Figure  9.  Comparison of shallow data before and after Q compensation

    a. Seismic data before compensation,b. Seismic data after compensation.

    10  中深层数据Q补偿结果对比

    a. 补偿前的地震数据,b. 补偿后的地震数据。

    10.  Comparison of deep data before and after Q compensation results

    a.Seismic data before compensation,b. Seismic data after compensation.

    11  实际数据Q补偿前后频谱对比

    11.  Spectrum comparison before and after Q compensation

    12  浅层数据Q补偿前后频谱对比

    12.  Spectrum comparison of shallow data before and after Q compensation

    13  中深层数据 Q补偿前后频谱对比

    13.  Spectrum compensation before and after Medium-deep data Q-compensation

    14  噪音成分Q补偿前后频谱对比

    14.  Spectrum comparison before and after noise component Q-compensation

  • [1] 李振春, 王清振. 地震波衰减机理及能量补偿研究综述[J]. 地球物理学进展, 2007, 22(4):1147-1152. [LI Zhenchun, WANG Qingzhen. A review of research on mechanism of seismic attenuation and energy compensation [J]. Progress in Geophysics, 2007, 22(4): 1147-1152. doi: 10.3969/j.issn.1004-2903.2007.04.021
    [2] 王珺, 杨长春, 乔玉雷. 用稳定高效的反Q滤波技术提高地震资料分辨率[J]. 地球物理学进展, 2008, 23(2):456-463. [WANG Jun, YANG Changchun, QIAO Yulei. Improve the seismic resolution with a stable and efficient inverse Q filter [J]. Progress in Geophysics, 2008, 23(2): 456-463.
    [3] Hale D. Q-adaptive deconvolution[C]//Proceedings of 1982 SEG Annual Meeting. Dallas, Texas: Society of Exploration Geophysicists, 1982: 133-158.
    [4] Bickel S H, Natarajan R R. Plane-wave Q deconvolution [J]. Geophysics, 1985, 50(9): 1426-1439. doi: 10.1190/1.1442011
    [5] Hargreaves N D, Calvert A J. Inverse Q filtering by Fourier transform [J]. Geophysics, 1991, 56(4): 519-527. doi: 10.1190/1.1443067
    [6] Bano M. Q-phase compensation of seismic records in the frequency domain [J]. Bulletin of the Seismological Society of America, 1996, 86(4): 1179-1186.
    [7] Wang Y H. A stable and efficient approach of inverse Q filtering [J]. Geophysics, 2002, 67(2): 657-663. doi: 10.1190/1.1468627
    [8] Wang Y H. Quantifying the effectiveness of stabilized inverse Q filtering [J]. Geophysics, 2003, 68(1): 337-345. doi: 10.1190/1.1543219
    [9] Wang Y H. Inverse Q-filter for seismic resolution enhancement [J]. Geophysics, 2006, 71(3): V51-V60. doi: 10.1190/1.2192912
    [10] Yan H Y, Liu Y. Estimation of Q and inverse Q filtering for prestack reflected PP-and converted PS-waves [J]. Applied Geophysics, 2009, 6(1): 59-69. doi: 10.1007/s11770-009-0009-y
    [11] 刘财, 刘洋, 王典, 等. 一种频域吸收衰减补偿方法[J]. 石油物探, 2005, 44(2):116-118. [LIU Cai, LIU Yang, WANG Dian, et al. A method to compensate strata absorption and attenuation in frequency domain [J]. Geophysical Prospecting for Petroleum, 2005, 44(2): 116-118. doi: 10.3969/j.issn.1000-1441.2005.02.005
    [12] Zhang C J, Ulrych T J. Seismic absorption compensation: a least squares inverse scheme [J]. Geophysics, 2007, 72(6): R109-R114. doi: 10.1190/1.2766467
    [13] 王本锋, 陈小宏, 李景叶, 等. 基于反演的稳定高效衰减补偿方法[J]. 地球物理学报, 2014, 57(4):1265-1274. [WANG Benfeng, CHEN Xiaohong, LI Jingye, et al. A stable and efficient attenuation compensation method based on inversion [J]. Chinese Journal of Geophysics, 2014, 57(4): 1265-1274. doi: 10.6038/cjg20140423
    [14] 叶秋焱, 崔宏良, 叶玮, 等. 一种基于叠前数据的综合Q值求取及应用[J]. 石油地球物理勘探, 2017, 52(2):304-308. [YE Qiuyan, CUI Hongliang, YE Wei, et al. Integrated Q-value estimation based on prestack seismic data [J]. Oil Geophysical Prospecting, 2017, 52(2): 304-308.
    [15] Gao J H, Yang S L, Wang D X, et al. Quality factors estimation using wavelet’s envelope peak instantaneous frequency[C]//2009 SEG Annual Meeting. Houston, Texas: SEG, 2009: 2457-2461.
    [16] Wang X J, Yin X Y, Wu G C. Estimation of layer quality factors from pre-stack seismic data[C]//Proceedings of the 80th Annual International Meeting. Denver, Colorado: SEG, 2010: 192-196.
    [17] Zhang F C, Zhang Y Y, Liu H Q. Seismic wave absorption attenuation compensation based on the matching pursuit algorithm[C]// (2015)SEG Technical Program Expanded Abstracts: 5179-5183.
    [18] 王文闯, 李合群. 混合法Q吸收补偿[J]. 石油地球物理勘探, 2012, 47(2):244-248. [WANG Wenchuang, LI Hequn. A hybrid method for Q-compensation [J]. Oil Geophysical Prospecting, 2012, 47(2): 244-248.
    [19] 李金丽, 李振春, 管路平, 等. 地震波衰减及补偿方法[J]. 物探与化探, 2015, 39(3):456-465. [LI Jinli, LI Zhenchun, GUAN Luping, et al. Seismic wave attenuation and compensation method [J]. Geophysical and Geochemical Exploration, 2015, 39(3): 456-465. doi: 10.11720/wtyht.2015.3.04
    [20] Battista B M, Knapp C, McGee T, et al. Application of the empirical mode decomposition and Hilbert-Huang transform to seismic reflection data [J]. Geophysics, 2007, 72(2): H29-H37. doi: 10.1190/1.2437700
    [21] Huang N E. Computer implicated empirical mode decomposition method, apparatus, and article of manufacture: US, U.S. Patent, 5983162[P]. 1996.
    [22] Mijovic B, De Vos M, Gligorijevic I, et al. Source separation from single-channel recordings by combining empirical-mode decomposition and independent component analysis [J]. IEEE Transactions on Biomedical Engineering, 2010, 57(9): 2188-2196. doi: 10.1109/TBME.2010.2051440
    [23] 颜中辉, 栾锡武, 王赟, 等. 基于经验模态分解的分数维地震随机噪声衰减方法[J]. 地球物理学报, 2017, 60(7):2845-2857. [YAN Zhonghui, LUAN Xiwu, WANG Yun, et al. Seismic random noise attenuation based on empirical mode decomposition of fractal dimension [J]. Chinese Journal of Geophysics, 2017, 60(7): 2845-2857. doi: 10.6038/cjg20170729
    [24] 梁岳, 顾汉明, 姚知铭. 改进的希尔伯特-黄变换在储层预测中的应用[J]. 石油物探, 2016, 55(4):606-615. [LIANG Yue, GU Hanming, YAO Zhiming. Application of improved Hilbert-Huang transform in reservoir prediction [J]. Geophysical Prospecting for Petroleum, 2016, 55(4): 606-615. doi: 10.3969/j.issn.1000-1441.2016.04.016
    [25] 颜中辉, 方刚, 徐华宁, 等. 希尔伯特谱白化方法在海洋地震资料高分辨率处理中的应用[J]. 海洋地质与第四纪, 2018, 38(4):212-220. [YAN Zhonghui, FANG Gang, XU Huaning, 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.
    [26] Braga L L S, Moraes F S. High-resolution gathers by inverse Q filtering in the wavelet domain [J]. Geophysics, 2013, 78(2): V53-V61. doi: 10.1190/geo2011-0508.1
    [27] 杨凯, 刘伟. 基于希尔伯特—黄变换的地震资料高分辨率处理[J]. 工程地球物理学报, 2015, 12(1):22-26. [YANG Kai, LIU Wei. The marine seismic data high resolution processing based on Hilbert-Huang transform [J]. Chinese Journal of Engineering Geophysics, 2015, 12(1): 22-26. doi: 10.3969/j.issn.1672-7940.2015.01.005
    [28] 王小杰, 印兴耀, 吴国忱. 基于叠前地震数据的地层Q值估计[J]. 石油地球物理勘探, 2011, 46(3):423-428. [WANG Xiaojie, YIN Xingyao, WU Guochen. Estimation of stratigraphic quality factors on pre-stack seismic data [J]. Petroleum Geophysical Prospecting, 2011, 46(3): 423-428.
    [29] 李合群, 孟小红, 赵波. 地震数据Q吸收补偿应用研究[J]. 石油地球物理勘探, 2010, 45(2):190-195, 229. [LI Hequn, MENG Xiaohong, ZHAO Bo. Application studies on Q absorption and compensation for seismic data [J]. Petroleum Geophysical Prospecting, 2010, 45(2): 190-195, 229.
    [30] 杨学亭, 刘财, 刘洋, 等. 基于连续小波变换的时频域地震波能量衰减补偿[J]. 石油物探, 2014, 53(5):1001-1441. [YANG Xueting, LIU Cai, LIU Yang, et al. The attenuation compensation of seismic wave energy in time-frequency domain based on the continuous wavelet transform [J]. Geophysical Prospecting for Petroleum, 2014, 53(5): 1001-1441.
    [31] 安静娟, 于亮. 一种改进的Q值求取方法[J]. 长江大学学报: 自科版, 2014, 11(2):90-91, 101. [AN Jingjuan, YU Liang. An improved method for obtaining Q value [J]. Journal of Yangtze University: Natural Science Edition, 2014, 11(2): 90-91, 101.
  • [1] 王小杰, 颜中辉, 刘欣欣, 杨佳佳.  基于小波分频的Q值补偿技术在东海中深层油气勘探中的应用 . 海洋地质与第四纪地质, 2019, 39(6): 200-206. doi: 10.16562/j.cnki.0256-1492.2019062101
    [2] 韦振权, 张莉, 帅庆伟, 易海, 钱星, 雷振宇, 王衍棠.  平衡剖面技术在台湾海峡盆地西部构造演化研究中的应用 . 海洋地质与第四纪地质, 2018, 38(5): 193-201. doi: 10.16562/j.cnki.0256-1492.2018.05.019
    [3] 徐云霞, 文鹏飞, 张宝金, 李福元.  琼东南海域OBS纵波资料成像处理关键技术 . 海洋地质与第四纪地质, 2018, 38(5): 185-192. doi: 10.16562/j.cnki.0256-1492.2018.05.018
    [4] 路彩晨, 贾佳, 高福元, 夏敦胜.  全新世安塞剖面的磁学特征变化历史及其受控因子分析 . 海洋地质与第四纪地质, 2018, 38(5): 178-184. doi: 10.16562/j.cnki.0256-1492.2018.05.017
    [5] 颜中辉, 方刚, 徐华宁, 刘俊, 施剑, 潘军, 王建强.  希尔伯特谱白化方法在海洋地震资料高分辨率处理中的应用 . 海洋地质与第四纪地质, 2018, 38(4): 212-220. doi: 10.16562/j.cnki.0256-1492.2018.04.019
    [6] 李三忠, 索艳慧, 郭玲莉, 戴黎明, 周立宏, 楼达.  渤海湾盆地大歧口凹陷变换构造与板内变形差异 . 海洋地质与第四纪地质, 2017, 37(4): 98-109. doi: 10.16562/j.cnki.0256-1492.2017.04.006
    [7] 冷传旭, 袁鸿洁, 徐翠玲, 郝娅楠, 赵广涛.  对数比变换在因子分析法提取东亚冬季风敏感粒级中的应用——以南黄海中部泥质区H07孔为例 . 海洋地质与第四纪地质, 2017, 37(1): 151-161. doi: 10.16562/j.cnki.0256-1492.2017.01.018
    [8] 刘俊, 吴淑玉, 施剑.  南黄海地震资料叠前去噪技术应用 . 海洋地质与第四纪地质, 2015, 35(6): 165-173. doi: 10.16562/j.cnki.0256-1492.2015.06.017
    [9] 蔡廷禄, 贾建军, 汪亚平.  河口海岸和近海沉积物的粒度资料同化技术 . 海洋地质与第四纪地质, 2014, 34(1): 185-193. doi: 10.3724/SP.J.1140.2014.01185
    [10] 李武, 朱怀平.  钻井中油气地球化学勘探技术及其应用 . 海洋地质与第四纪地质, 2012, 32(1): 159-166. doi: 10.3724/SP.J.1140.2012.01159
    [11] 张光学, 张明, 杨胜雄, 雷新华, 徐华宁, 刘学伟, 梁金强, 沙志彬.  海洋天然气水合物地震检测技术及其应用 . 海洋地质与第四纪地质, 2011, 31(4): 51-58. doi: 10.3724/SP.J.1140.2011.04051
    [12] 赵铁虎, 张训华, 冯京.  海底油气渗漏浅表层声学探测技术 . 海洋地质与第四纪地质, 2010, 30(6): 149-156. doi: 10.3724/SP.J.1140.2010.06149
    [13] 杨丹, 扈传昱, 于培松, 武光海, 倪建宇, 张海生, 姚龙奎, 卢冰.  中、东太平洋沉积岩样低成熟烃组成与演化及其受控因子对比 . 海洋地质与第四纪地质, 2010, 30(6): 99-106. doi: 10.3724/SP.J.1140.2010.06099
    [14] 肖军, 肖应凯, 刘丛强, 贺茂勇.  海洋生物碳酸盐δ11B-pH技术影响因素 . 海洋地质与第四纪地质, 2010, 30(4): 215-222. doi: 10.3724/SP.J.1140.2010.04215
    [15] 石正国, 刘晓东, 程肖侠.  东亚季风和南亚季风的反相位关系及其与ENSO的联系:年际尺度和轨道尺度 . 海洋地质与第四纪地质, 2009, 29(5): 83-88. doi: 10.3724/SP.J.1140.2009.05083
    [16] 徐继尚, 李广雪, 曹立华, 杨荣民.  海底管道综合探测技术及东方1-1管道不稳定因素 . 海洋地质与第四纪地质, 2009, 29(5): 43-50. doi: 10.3724/SP.J.1140.2009.05043
    [17] 张荣强, 吴时国, 周雁, 余朝华, 冯德勇, 于正军.  平衡剖面技术及其在济阳坳陷桩海地区的应用 . 海洋地质与第四纪地质, 2008, 28(6): 135-142. doi: 10.3724/SP.J.1140.2008.06135
    [18] 张明, 彭朝旭, 沙志彬.  天然气水合物准三维地震调查导航定位技术 . 海洋地质与第四纪地质, 2008, 28(6): 101-106. doi: 10.3724/SP.J.1140.2008.06101
    [19] 业渝光, 张剑, 胡高伟, 刁少波, 刘昌岭.  天然气水合物超声和时域反射联合探测技术 . 海洋地质与第四纪地质, 2008, 28(5): 101-107.
    [20] 钟广见, 文鹏飞, 罗文造, 赵庆献, 陆敬安.  电缆沉放深度对地震资料品质的影响 . 海洋地质与第四纪地质, 2007, 27(6): 131-134.
  • 加载中
图(14)
计量
  • 文章访问数:  578
  • HTML全文浏览量:  270
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-09
  • 修回日期:  2019-04-16
  • 网络出版日期:  2020-01-15
  • 刊出日期:  2020-02-01

基于HHT方法的地层Q值补偿方法研究

doi: 10.16562/j.cnki.0256-1492.2018100901
    作者简介:

    颜中辉(1986—),男,硕士,助理研究员,主要从事地震资料处理、多波多分量研究,E-mail:zhonghui4564573@163.com

    通讯作者: 徐华宁,男,博士,教授级高工,主要从事海洋地球物理调查研究,E-mail:xuhuaning1973@163.com
基金项目:  山东省自然科学基金培养基金项目“基于HHT变换的多波地震资料去噪应用研究”(ZR2015PD006);山东省自然基金博士基金 “地质约束下的重磁三维反演在南黄海盆地崂山隆起的应用”(ZR2016DB33);青岛海洋科学与技术国家实验室开放基金“多分量地震波逆时偏移拉伸校正及其并行”(QNLM2016ORP0206);国家自然科学基金青年基金“东海北部外陆架晚更新世以来两期古三角洲的时空展布特征及物源属性研究”(41706069)
  • 中图分类号: P631.4

摘要: 常规的基于叠后地震资料的地层Q值补偿方法受数据平均效应的影响,Q值估算精度较低,导致保真性差,而HHT时频分析可以反映出信号的局部信息,基于以上原因,提出一种基于HHT方法的地层Q值补偿方法。首先应用HHT分频技术将浅层、中深层地震数据分成不同的频率范围,得到不同目标层频率尺度的地震数据;然后在HHT时频域内通过谱比法提取Q值;最后分别进行地层Q值振幅补偿和相位补偿。该方法针对不同的目标地层将地震数据分解成不同的频率尺度,避免了不同频率成分之间的干扰,而且可针对性处理目标层的频率信息,提高了地层Q值提取的稳定性。模型试算和实际资料应用表明,该方法能够有效地对振幅和相位进行补偿,使中深层反射波形得到恢复,主频提高,频带展宽,地震资料的品质得到提升。

English Abstract

  • 随着复杂油气藏勘探精度的提高,精细的地震属性反演、油气储层的定量描述都需要高分辨率的地震资料。而地层的吸收衰减作用是影响地震资料分辨率的一个主要因素[1],从而导致了有效反射波频带变窄,相位产生畸变,分辨率和信噪比降低,从而影响了地震资料的成像精度。Q补偿作为一种确定性的处理技术[2],不仅可以补偿地震反射信号的频率和振幅衰减,同时能改善地震子波的相位特性,使中深层的高频部分得到加强,频谱被拓宽,从而有效提高地震资料的质量。

    反Q滤波补偿方面也有许多学者进行了研究,基本可以分为3类:用级数展开作近似高频补偿的反Q滤波、基于波场延拓的反Q 滤波和其他反Q滤波方法。最早Hale提出Futterman模型进行反Q滤波的方法[3],该方法基于预测误差滤波分析迭代求出Q值,然后进行Q补偿,其缺点是计算量大。Bickel和Natarajan根据复函数分析技术提出一种反Q滤波方法,此方法基于平面波假设,通过复数对平面波传播进行描述,从而得到时变的Q补偿[4]。Hargreaves和Calvert根据常Q模型得到相位Q补偿方法[5],其做法是通过波场外推来计算常Q模型,从而对相位进行补偿,缺点是未考虑振幅的影响。Bano利用常Q模型发展到层常Q值模型的相位Q补偿方法,优势在于补偿结果能稳定地改善相位的畸变,而且准确性比较高,但同样忽略了振幅的影响[6]。WANG基于等间隔形式的Q模型提出一种稳定有效的反Q滤波方法,此方法以波场延拓理论为基础,并对Q值求取的抗噪性和稳定性进行深入讨论,并不断完善,形成了最常用的反Q滤波算法[7-9]。严红勇和刘洋在WANG方法的基础上,发展成沿射线路径进行波场延拓从而进行Q补偿的方法,主要是针对多分量资料进行处理[10]。刘财等提出基于频率域分时窗吸收补偿的反Q滤波方法[11];Zhang和Ulrych在反Q滤波过程中引入了最小二乘方法和Bayasian原理,从而弥补了Q求取的不稳定性,并在实际数据中得到了很好的效果[12];王本锋等结合反演的思想正则化方法提出了一种新的补偿方法,该方法基于波场延拓的正Q滤波方程,对有效频带内的频率分量进行分析,提高了计算效率,补偿结果稳定高效[13]。叶秋焱利用广义S变换求取时频域表层Q值的方法,通过表层Q值与井控深层Q值建立综合Q体,从而对地震数据进行Q补偿[14]

    虽然地层Q值补偿方法的理论和算法比较成熟,但是在实际地震资料处理应用中一直存在其稳定性和抗噪性不足的问题[15-19],尤其是地层Q值的准确性。本文在HHT方法的基础上,基于EMD自适应分解和HHT时频谱局部刻画能力[20-23],通过将地震资料分成不同的频率段,在HHT时频域计算不同目标层的地层Q值,从而进行Q值补偿。确保了Q值提取精度和Q值补偿的准确性。

    • HHT方法的核心内容是固有模态分解(EMD)[24-25],其主要思路是将具有复杂尺度成分的非平稳、非线性信号分解成一组较优Hilbert变换性能的固有模态函数(IMF)之和。我们给定复杂信号X(t),经过EMD自适应分解可表示为:

      $ X\left( t \right) = \sum\nolimits_{i = 1}^n {{c_i}\left( t \right) + {r_n}\left( t \right)} $

      (1)

      式中,ci(t)是固有模态函数,rn(t)称为余量。EMD分解其优势在于信号分解过程是自适应的。得到的各个不同尺度IMF都是平稳窄带信号,从而使得波形轮廓对称更好。实际上,EMD方法主要是从频率特征尺度进行分解,第一步得到的是信号中频率特征尺度最小的IMF分量;第二步得到频率特征时间尺度较大的IMF分量;最后分离出来的IMF具有频率特征时间尺度最大的特征。

      原信号X(t)分解成多个不同尺度IMF后,基于复数道构建技术,可以构成解析信号Z(t),

      $ Z\left( t \right) = X\left( t \right) + iY\left( t \right) $

      (2)

      $ Z\left( t \right)\text{可写为:}\;\;\;\;\;\;Z\left( t \right) = a\left( t \right){{\rm{e}}^{i\theta \left( t \right)}} $

      (3)

      Y(t)为X(t)经过Hilbert变换的信号,每一个IMF分量经过Hilbert变换后,构造解析函数可以表示为:

      $ H\left( t \right) = \sum\nolimits_{j + 1}^n {{a_j}\left( t \right)\exp \left[ {i\int {{w_j}\left( t \right){\rm{d}}t} } \right]} $

      (4)

      其中,aj(t)为分解的第j阶IMF分量,wj(t)为特定的瞬时频率,对照上式看出,HHT变换实际是傅里叶变换的扩展,优势在于容易反映局部振幅和频率的信息。去除残余函数,展开后得到Hilbert幅值谱,记作:

      $ H\left( {w,t} \right) = Re\sum\nolimits_{j + 1}^n {{a_j}\left( t \right){{\rm{e}}^{i\int {{w_j}\left( t \right){\rm{d}}t} }}} $

      (5)

      基于HHT分频技术的地层Q值提取方法如下:首先基于HHT时频分析技术对目的层进行分析,选择适合目的层的优势频率范围;然后利用EMD分解技术对地震数据根据频率尺度进行分解,选择只含有目的层优势频率范围的地震数据重构;最后在HHT时频域内通过谱比法计算地层Q值[26]。本文方法在实际地震数据对地层Q补偿应用中的具体步骤为:将实际地震数据分成浅层、中深层,分别对这两套数据进行HHT时频分析,提取局部的频率属性特征,选择适合不同层的频率范围,基于EMD分解技术对这两套数据进行分解、选择、重构处理[27],得到具有不同频率尺度的地震资料;基于复数道技术,在HHT时频域内通过谱比法计算Q值。此做法的优势在于不同频率成分地震信号估算的Q值避免了造成相互影响,同时局部的Q值异常求得更加准确,补偿效果更佳。HHT时频域内谱比法是假定地震子波为零相位子波,其估算地层Q 值的公式为:

      $ \ln \frac{{W\left( {a,{t_1}} \right)}}{{W\left( {a,t} \right)}} = - \alpha \left( {\frac{{{t_1} - t}}{Q}} \right) $

      (6)

      其中,$\alpha = {\rm{\pi }}\left( {\dfrac{{{f_1} + \tau f}}{{\tau + 1}}} \right)$

      式中,τ表示能量衰减率,为一常数,f1表示零相位子波的视频率;f表示平面波频率;tt1表示平面波从地面开始传播到目的层上下界面所用的时间。W(at)表示HHT时频谱。通过公式(6)进行线性回归可以求得地层Q值[28]

    • 地震波在地下传播时,由于信号的中心频率要向较低的频率方向移动,因此找出Q和中心频率的关系至关重要。以此通过Q补偿可以消除地震波在地下介质中传播的频散吸收,提高地震数据的垂向分辨率[29-31]

      在均匀介质中传播的平面波频域表达式为:

      $ U\left( f \right) = {U_0}\left( f \right)\exp \left[ {{{i}}\left( {2{\rm{\pi }}ft - kx} \right)} \right] $

      (7)

      式中,U0(f)表示震源脉冲响应;k为波数;x为波前面沿着射线路径从震源传播到检波器的距离。假如介质是粘弹性,则k为复数,表达式如下:

      $ k\left( f \right) = \frac{{2{\rm{\pi }}f}}{{c\left( f \right)}} = \frac{{2{\rm{\pi }}f}}{{v\left( f \right)}} - i\alpha \left( f \right) $

      (8)

      其中,c(f)为复速度,v(f)为相速度,α(f)为衰减系数。

      在Kolsky-Futterman模型中,相速度和衰减系数对应的关系满足:

      $ \alpha \left( f \right) = \frac{{\rm{\pi }}}{{v\left( f \right)Q}}f $

      (9)

      $ \frac{1}{{v\left( f \right)}} = \frac{1}{{{v_r}}}\left( {1 - \frac{1}{{{\rm{\pi }}Q}}\ln \left| {\frac{f}{{{f_r}}}} \right|} \right) $

      (10)

      式中,fr代表参考频率;vr代表参考相速度。

      于是将公式(8)至公式(10)代入公式(7)中,展开后对应表达式:

      $ U\left( f \right) = {U_0}\left( f \right)\exp \left[ { - \left( {\frac{{{\rm{\pi }}ft}}{Q} + {{i}}\frac{{2ft}}{Q}\ln \left| {\frac{f}{{{f_r}}}} \right|} \right)} \right] $

      (11)

      式中,$t = \dfrac{x}{{{v_r}}}$

      定义$G\left( {f,t} \right) = \exp \left[ { - \left( {\dfrac{{{\rm{\pi }}ft}}{Q} + {{i}}\dfrac{{2ft}}{Q}\ln \left| {\dfrac{f}{{{f_r}}}} \right|} \right)} \right]$大地滤波算子,则反Q滤波算子为:

      $ {G^{ - 1}}\left( {f,t} \right) = A\left( {f,t} \right)P\left( {f,t} \right) = {\rm{exp}}\left( {\frac{{{\rm{\pi }}ft}}{Q} + {{i}}\frac{{2ft}}{Q}{\rm{ln}}\left| {\frac{f}{{{f_r}}}} \right|} \right) $

      (12)

      式中,$P\left( {f,t} \right) = {\rm{exp}}\left( {i\dfrac{{2ft}}{Q}{\rm{ln}}\left| {\dfrac{f}{{{f_r}}}} \right|} \right)$是相位补偿算子,$A\left( {f,t} \right) = {\rm{exp}}\dfrac{{{\rm{\pi }}ft}}{Q}$是振幅补偿算子。

      由此推导出地震波衰减补偿在时频域的表达式:

      $ {U_0}\left( {f,t} \right) = U\left( {f,t} \right)A\left( {f,t} \right)P\left( {f,t} \right) $

      (13)

      最后进行HHT反变换,得到补偿后的时间域地震信号。

    • 基于上述理论基础,首先进行模型参数的试验,通过层状介质的衰减模拟分析,以及Q值求取与补偿研究对理论进行验证。模型参数如下:对应层的速度依次为1 500、2 000、3 000 m/s,Q值依次为120、60、120。正演模拟采用基于单程波动方程的分步傅里叶波场延拓方法,采用的地震子波为雷克子波,频率为30 Hz,时间采样间隔1 ms。图1a为模拟得到的地震记录,图2a为放大的波形显示。由于地下介质吸收衰减作用,体现出振幅减小和波形畸变的现象,从放大的波形分析,子波波形不再是零相位,即分辨率逐渐降低。

      图  1  模型数据Q补偿前后的对比

      Figure 1.  Comparison of model seismic data before and after Q compensation

      图  2  模型数据Q补偿前后的波形显示对比

      Figure 2.  Waveform comparison of model seismic data before and after Q compensation

      图3是基于小波分频技术提取的地层Q值曲线,与模型参数基本一致。图1b图2b是经过本文方法进行补偿后的结果。对比图1a图2a可以看出,通过Q补偿处理的不同深度地层反射波形均恢复到了原始的雷克子波形态,对应每一层地震信号的振幅和相位都得到了恢复。图4是经过Q补偿前后的振幅谱对比,蓝色为补偿前的频谱,红色为补偿后的频谱。由图可以看出,补偿后的频谱得到展宽。由此可以说明,通过本文方法对地层进行Q值补偿,能够达到使地震子波波形恢复原状,主频提高,频带变宽,提高分辨率的作用。

      图  3  基于HHT方法提取的地层Q值曲线

      Figure 3.  Formation Q-value curve extracted based on HHT method

      图  4  模型数据Q补偿前后的频谱对比

      Figure 4.  Spectrum comparison before and after model data Q compensation

      为了验证本文Q补偿方法的抗噪性,通过对模型数据加入8%的随机噪声,同样的步骤对加噪的数据进行Q补偿测试。图5为补偿前后的数据对比,图6为补偿前后的频谱对比。从图中的对比分析,补偿后的地震记录有效信息得到加强恢复,噪声的部分基本没变。也就是说本文方法在信号补偿的同时,也能确保信噪比。

      图  5  加噪的模型数据Q补偿前后的对比

      Figure 5.  Comparison of model seismic data with noise before and after Q compensation

      图  6  加噪的模型数据Q补偿前后的频谱对比

      Figure 6.  Spectrum comparison of model seismic data with noise before and after Q compensation

    • 本文选择某海域地震资料实际数据对本文方法进行应用效果分析。选取某一段信噪比较好的资料,对其浅层、中深层进行Hilbert分频处理,选择适合不同层的频率区间在HHT时频域计算地层的Q值,进而进行反Q滤波,实现对地层振幅和相位的补偿。图7对应的是根据本文方法求得的Q值剖面,图8a是原始的叠后剖面,从图中可看出Q值剖面跟地层吻合较好。图8b是整个剖面进行补偿后的效果,可以看出分辨率有了显著提高,同时剖面的信噪比也能得到很好的保持。图9为浅层经过补偿后的对比图,补偿后同相轴的连续性有了明显的改善。图10是中深层补偿后的对比图,同样的,补偿后的剖面分辨率更高,地层信息更加丰富。

      图  7  基于HHT方法提取的地层Q值剖面

      Figure 7.  Formation Q-value profile extracted with HHT method

      图  8  实际地震数据Q补偿结果对比

      Figure 8.  Comparison of actual seismic data before and after Q compensation

      图  9  浅层数据Q补偿结果对比

      Figure 9.  Comparison of shallow data before and after Q compensation

      图  10  中深层数据Q补偿结果对比

      Figure 10.  Comparison of deep data before and after Q compensation results

      图11是整个补偿前后振幅谱对比图,图12 为浅层补偿前后频谱对比图。对比可知,Q补偿后的数据频谱有了较大的改善,从浅层到深层表现出高频能量提高,数据主频范围向高频移动,且频谱的形态基本没发生变换,说明补偿是有效的,同时频带变宽,对应的剖面高频成分也有很明显的补偿效果,分辨率有所提高。尤其对中深层的补偿,从图13分析可知,原始资料的高频部信息缺失严重,通过Q补偿后,高频的有效信息被恢复。由于实际数据为海上资料,补偿后高频部分振幅能量得到加强,相应的会导致海上虚反射引起的陷波点效应加剧,海上常规资料处理的方法是通过鬼波压制解决。

      图  11  实际数据Q补偿前后频谱对比

      Figure 11.  Spectrum comparison before and after Q compensation

      图  12  浅层数据Q补偿前后频谱对比

      Figure 12.  Spectrum comparison of shallow data before and after Q compensation

      图  13  中深层数据 Q补偿前后频谱对比

      Figure 13.  Spectrum compensation before and after Medium-deep data Q-compensation

      为了验证本文方法在提高分辨率的同时,不会损失信噪比,提取噪音部分进行频谱分析,如图14所示,Q补偿后的频谱基本没变换,也就是说噪音的成分没有加强。从补偿的剖面来看,存在信噪比降低的表象,其实只是整个资料主频向高频移动的表现。

      图  14  噪音成分Q补偿前后频谱对比

      Figure 14.  Spectrum comparison before and after noise component Q-compensation

    • (1)HHT变换具有很高的局部分析能力,因此,基于HHT时频谱的地层Q值提取方法能得到更加准确的地层Q值,同时EMD分解算法是基于频率尺度进行分解,可针对不同目标区频率区间进行地震数据的Q补偿,从而达到提高资料分辨率的目的。

      (2)通过本文方法进行Q值补偿,在有效提高地震资料主频,展宽频带,提高分辨率的同时,确保了资料的信噪比。从模拟数据到实际地震数据的测试表明本文方法是可行且有效的。

      (3)由于EMD分解的过程存在模态混叠效应,同时Q值提取的正确性对资料信噪比有一定的要求。因此,在处理之前应对地震数据进行去噪预处理,以达到更佳的补偿效果。

参考文献 (31)

目录

    /

    返回文章
    返回