油气勘探

基于最大熵准则的Wigner-Ville分布与微型古河道刻画方法及应用

  • 徐天吉 ,
  • 程冰洁 ,
  • 牛双晨 ,
  • 秦正晔 ,
  • 王贞贞
展开
  • 1.电子科技大学资源与环境学院,成都 611731;
    2.电子科技大学长三角研究院(湖州),浙江湖州 313000;
    3.成都理工大学“油气藏地质及开发工程”国家重点实验室,成都 610059;
    4.成都理工大学“地球探测与信息技术”教育部重点实验室,成都 610059;
    5.四川中成煤田物探工程院有限公司,成都 610072
徐天吉(1975-),男,四川射洪人,博士,电子科技大学资源与环境学院研究员,主要从事多波地震勘探、储集层地质力学、人工智能等方法研究与教学。地址:成都市高新区(西区)西源大道2006号创新中心,邮政编码:611731。E-mail: xutianji@uestc.edu.cn

收稿日期: 2021-01-29

  修回日期: 2021-11-09

  网络出版日期: 2021-11-25

基金资助

国家自然科学基金面上项目“基于应力诱导各向异性与岩石力学性质的深层页岩储层可压裂性评价方法研究”(42074160); 国家自然科学基金面上项目“基于各向异性介质弹性参数的页岩TOC及脆性预测方法”(41574099); 四川省科技计划“页岩储层智能化含气性预测方法”(2020JDRC0013)

A microscopic ancient river channel identification method based on maximum entropy principle and Wigner-Ville Distribution and its application

  • XU Tianji ,
  • CHENG Bingjie ,
  • NIU Shuangchen ,
  • QIN Zhengye ,
  • WANG Zhenzhen
Expand
  • 1. School of Resources and Environment, University of Electronic Science and Technology of China, Chengdu 611731, China;
    2. Yangtze Delta Region Institute of University of Electronic Science and Technology of China, Huzhou 313000;
    3. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu 610059, China;
    4. College of Geophysics, Chengdu University of Technology, Chengdu 610059, China;
    5. Sichuan Zhongcheng Institute of Coalfield Geophysical Engineering, Chengdu 610072

Received date: 2021-01-29

  Revised date: 2021-11-09

  Online published: 2021-11-25

摘要

针对窄细古河道精细刻画难题,利用最大熵准则增强Wigner-Ville分布的聚焦特性,在有效提升地震信号时频分辨能力的基础上,建立了一种微型古河道识别新方法。基于最大熵功率谱与自回归模型(AR)功率谱等效的原理,首先利用Burg算法和Levinson-Durbin递推规则,求取AR模型的预测误差、自回归系数等参数;然后,在自相关函数一阶导数为0的条件下,计算地震信号的Wigner-Ville分布,获取微型古河道最大熵准则约束下的Wigner-Ville时频功率谱(MEWVD)。通过仿真地震信号和窄薄模型数值模拟信号实验分析,发现MEWVD既能有效避免Wigner-Ville分布的交叉项干扰,还能获得比短时傅里叶变换(STFT)、连续小波变换(CWT)等信号分析方法更加精准的频谱特征;同时,还证实了利用不同频率的MEWVD,可以有效识别不同尺度的窄细古河道。将该方法应用于四川盆地中江气田侏罗系沙溪庙组(J2s33-2小层)气藏,准确地识别出宽度小于500 m、砂岩厚度小于35 m的窄细古河道的宽度、走向等空间信息,可为井位部署、水平井压裂选段等提供依据。图7表2参31

本文引用格式

徐天吉 , 程冰洁 , 牛双晨 , 秦正晔 , 王贞贞 . 基于最大熵准则的Wigner-Ville分布与微型古河道刻画方法及应用[J]. 石油勘探与开发, 2021 , 48(6) : 1175 -1186 . DOI: 10.11698/PED.2021.06.09

Abstract

In view of the problem of fine characterization of narrow and thin channels, the maximum entropy criterion is used to enhance the focusing characteristics of Wigner-Ville Distribution. On the basis of effectively improving the time-frequency resolution of seismic signal, a new method of microscopic ancient river channel identification is established. Based on the principle of the equivalence between the maximum entropy power spectrum and the AR model power spectrum, the prediction error and the autoregression coefficient of AR model are obtained by using the Burg algorithm and Levinson-Durbin recurrence rule. Under the condition of the first derivative of autocorrelation function being 0, the Wigner-Ville Distribution of seismic signal is calculated, and the Wigner-Ville Distribution time-frequency power spectrum (MEWVD) is obtained under the maximum entropy criterion of the microscopic ancient river channel. Through analysis of emulational seismic signal and numerical simulation signal of narrow thin model, it is found that MEWVD can effectively avoid the interference of cross term of Wigner-Ville Distribution, and obtain more accurate spectral characteristics than STFT and CWT signal analysis methods. It is also proved that the narrow and thin river channels of different scales can be identified effectively by using MEWVD of different frequencies. The method is applied to the third member of Jurassic Shaximiao Formation (J2s33-2) gas reservoir of the Zhongjiang gas field in Sichuan Basin. The spatial information of width and direction of narrow and thin river channel with width less than 500 m and sandstone thickness less than 35 m is accurately identified, providing basis for well deployment and horizontal well fracturing section selection.

参考文献

[1] 刘成川, 曹廷宽, 卜淘. 薄窄河道致密砂岩气藏高效开发技术对策[J]. 大庆石油地质与开发, 2020, 39(2): 157-165.
LIU Chengchuan, CAO Tingkuan, BU Tao. High-efficiency developing techniques for the tight sandstone gas reservoir with thin and narrow channels[J]. Petroleum Geology & Oilfield Development in Daqing, 2020, 39(2): 157-165.
[2] 梁宏刚, 先伟, 李文平. 塔河油田超深层、薄储层、窄河道砂体识别与描述技术[J]. 地质科技情报, 2015, 34(5): 180-183.
LIANG Honggang, XIAN Wei, LI Wenping. Identification and description technique for ultra-deep, thin reservoir, narrow, channel sand, Tahe oil field[J]. Geological Science and Technology Information, 2015, 34(5): 180-183.
[3] WANG Y. Inverse Q-filter for seismic resolution enhancement[J]. Geophysics, 2006, 71(3): 51-60.
[4] 李胜强, 刘振东, 严加永, 等. 高分辨深反射地震探测采集处理关键技术综述[J]. 地球物理学进展, 2020, 35(4): 1400-1409.
LI Shengqiang, LIU Zhendong, YAN Jiayong, et al. Review on the key techniques for high resolution deep reflection seismic acquisition and processing[J]. Progress in Geophysics, 2020, 35(4): 1400-1409.
[5] MICHAL S, GARY P, JAIME S, et al. Extending seismic bandwidth using the continuous wavelet transform[J]. First Brieak, 2008, 26(6): 97-102.
[6] CHOPRA S, CASTAGNA J P, XU Y. Thin-bed reflectivity inversion and some applications[J]. First Break, 2009, 27(5): 27-34.
[7] SUN H Y, DEMANET L. Low-frequency extrapolation with deep learning[R]. Anaheim: 88th Society of Exploration Geophysicists International Exposition and Annual Meeting, 2018.
[8] 李祺鑫, 罗亚能, 张生, 等. 高分辨率波阻抗贝叶斯序贯随机反演[J]. 石油地球物理勘探, 2020, 55(2): 389-398.
LI Qixin, LUO Yaneng, ZHANG Sheng, et al. High-resolution Bayesian sequential stochastic inversion[J]. Oil Geophysical Prospection, 2020, 55(2): 389-398.
[9] 陈彦虎, 毕建军, 邱小斌, 等. 地震波形指示反演方法及其应用[J]. 石油勘探与开发, 2020, 47(6): 1149-1158.
CHEN Yanhu, BI Jianjun, QIU Xiaobin, et al. A method of seismic meme inversion and its application[J]. Petroleum Exploration and Development, 2020, 47(6): 1149-1158.
[10] GABOR D. Theory of communication[J]. Journal of the Institute of Electrical Engineers of Japan, 1946, 93(26): 429-457.
[11] 黄昱丞, 郑晓东, 栾奕, 等. 地震信号线性与非线性时频分析方法对比[J]. 石油地球物理勘探, 2018, 53(5): 975-989.
HUANG Yucheng, ZHENG Xiaodong, LUAN Yi, et al. Comparison of linear and nonlinear time-frequecy analysis on seismic signals[J]. Oil Geophysical Prospection, 2018, 53(5): 975-989.
[12] WAERDEN B L van der. Sources of quantum mechanics[M]. Amsterdam: North-Holland Publishing Co, 1967.
[13] PEI S C, HUANG S G. STFT with adaptive window width based on the chirp rate[J]. IEEE Transactions on Signal Processing, 2012, 60(8): 4065-4080.
[14] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings Mathematical Physical & Engineering Sciences, 1998, 454(1971): 903-995.
[15] TORRES M E, COLOMINAS M A, SCHLOTTHAUER G, et al. A complete ensemble empirical mode decomposition with adaptive noise[R]. Prague: IEEE International Conference on Acoustics, 2011.
[16] DAUBECHIES I, LU J, WU H T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool[J]. Applied and Computational Harmonic Analysis, 2011, 30(2): 243-261.
[17] MALLAT S G, ZHANG Z F. Matching pursuits with time-frequency dictionaries[J]. IEEE Transaction on Signal Processing, 1993, 41(12): 3397-3415.
[18] WIGNER E P. On the quantum correction for thermodynamic equilibrium[J]. Physical Review, 1932, 40(5): 749-759.
[19] CHOI H I, WILLIAMS W J. Improved time-frequency representation of multicomponent signals using exponential kernels[J]. IEEE Transactions on Acoustics Speech & Signal Processing, 2002, 37(6): 862-871.
[20] SHANNON C E. A mathematical theory of communication[J]. Mobile Computing and Communications Review, 1948, 5(1): 3-55.
[21] JAYNES E T. Information theory and statistical mechanicsⅡ[J]. Physical Review, 1957, 108(2): 171-190.
[22] BURG J P. A new analysis technique for time series data[R]. Enschede, the Netherlands: NATO Advanced Study Institute on Signal Processing with Emphasis on Underwater Acoustics, 1968.
[23] BOASHASH B. Estimating and interpreting the instantaneous frequency of a signal: Part I: Fundamentals[J]. Proceedings of the IEEE, 1992, 80(4): 520-538.
[24] BOS A V D. Alternative interpretation of maximum entropy spectral analysis[J]. IEEE Trans on Information Theory, 1971, 17(4): 493-494.
[25] LEVINSON N. The Wiener RMS (root mean square) criterion in filter design and prediction[J]. Journal of Mathematical Physics, 1946, 25(1/2/3/4): 261-278.
[26] KAY S M, MARPLE S L. Spectrum analysis a modern perspective[J]. Proceedings of the IEEE, 1981, 68(11): 1380-1419.
[27] MARPLE S L. A new autoregressive spectrum analysis algorithm[J]. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1980, 28(4): 441-454.
[28] BURG J P. Maximum entropy spectrum analysis[D]. Stanford, California: Stanford University, 1975.
[29] 谢林江, 尹东. 基于最大相关熵准则的鲁棒度量学习算法[J]. 计算机系统应用, 2018, 27(10): 146-153.
XIE Linjiang, YIN Dong. Robust metric learning based on maximum correntropy criterion[J]. Computer Systems & Applications, 2018, 27(10): 146-153.
[30] 马永生, 蔡勋育, 赵培荣, 等. 四川盆地大中型天然气田分布特征与勘探方向[J]. 石油学报, 2010, 31(1): 347-354.
MA Yongsheng, CAI Xunyu, ZHAO Peirong, et al. Distribution and further exploration of the large-medium sized gas fields in Sichuan Basin[J]. Acta Petrolei Sinica, 2010, 31(1): 347-354.
[31] 魏国齐, 杨威, 刘满仓, 等. 四川盆地大气田分布、主控因素与勘探方向[J]. 天然气工业, 2019, 39(6): 1-12.
WEI Guoqi, YANG Wei, LIU Mancang, et al. Distribution rules, main controlling factors and exploration directions of giant gas fields in the Sichuan Basin[J]. Natural Gas Industry, 2019, 39(6): 1-12.
文章导航

/