Surface structure survey by microtremor based on F-J method
ZHAO Rongrong1, YANG Zhentao2,3, REN Chenghao4, ZHOU Xin2, WANG Yi1, CHEN Wei1
1. PetroChina Southwest Oil & Gasfield Company Exploration Division, Chengdu, Sichuan 610041, China; 2. Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China; 3. Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology (Southern University of Science and Technology), Shenzhen, Guangdong 518055, China; 4. GME and Geochemical Surveys, BGP Inc., Zhuozhou, Hebei 072751, China
Abstract:3D seismic exploration is typically employed in areas where micro logging and other static correction methods cannot be implemented, such as towns and roads. Recently researchers have tried to replace existing static correction methods with microtremor surveys and yielded some achievements. At present, in petroleum exploration, microtremor survey mainly employs surface wave dispersion curves to invert formation velocity structures. This method is still in its infancy in petroleum exploration. As the core technology, dispersion imaging methods can be classfied into many types, such as the frequency-Bessel(F-J) transformation method proposed in 2019, which has attracted the attention of geophysicists. Compared with the spatial autocorrelation method(SPAC), the F-J method has great advantages in adopting ambient noise for higher-mode surface wave imaging. To better study the applicability of F-J method, this paper first analyzes the influence of different observation systems and acquisition param-eters on F-J imaging using synthetic data. Then, the effect of different observation systems and acquisition time on FJ dispersion imaging is analyzed through the field ambient noise data of an area in Sichuan. In practical application, the dispersion spectrum obtained by subsection optimization superposition can effectively improve the imaging quality of the dispersion spectrum and provide a data basis for the accurate extraction of dispersion curves. Finally, logging data comparison shows that multi-mode joint inversion using fundamental mode and higher-mode dispersion curves can accurately invert formation structures, thus making it possible to obtain surface structures economically and efficiently.
赵东. 被动源面波勘探方法与应用[J]. 物探与化探, 2010, 34(6):759-764.ZHAO Dong. Passive surface waves:methods and applications[J]. Geophysical and Geochemical Exploration, 2010, 34(6):759-764.
[2]
AKI K. Space and time spectra of stationary stochastic waves, with special reference to microtremors[J]. Bulletin of the Earthquake Research Institute, 1957, 35(3):415-456.
[3]
CHO I, TADA T, SHINOZAKI Y. A new method of microtremor exploration using miniature seismic arrays:quick estimation of average shear velocities of the shallow soil[J]. BUTSURI-TANSA, 2008, 61(6):457-468..
[4]
LING S, OKADA H. An extended use of the spatial autocorrelation method for the estimation of geological structure using microtremors[C]. Proceedings of the 89th SEGJ Conference, Japanese, 1993, 44-48.
[5]
OKADA H, SUTO K, ASTEN M W. The microtremor survey method[J]. Geophysical Monograph Series No 12, SEG, Tulsa, 2003, https://doi. org/10.1190/1.9781560801740.
[6]
ASTEN M, ROBERTS J. Analysis of esg2006 blindtrial microtremor data using the MMSPAC method[J]. Journal of Telemedicine & Telecare, 2006, 13(1):7-14.
[7]
徐佩芬, 李传金, 凌甦群, 等. 利用微动勘察方法探测煤矿陷落柱[J]. 地球物理学报, 2009, 52(7):1923-1930.XU Peifen, LI Chuanjin, LING Suqun, et al. Mapping collapsed columns in coal mines utilizing microtremor survey methods[J]. Chinese Journal of Geophysics, 2009, 52(7):1923-1930.
[8]
徐佩芬, 李世豪, 凌甦群, 等. 利用SPAC法估算地壳S波速度结构[J]. 地球物理学报, 2013, 56(11):3846-3854.XU Peifen, LI Shihao, LING Suqun, et al. Application of SPAC method to estimate the crustal S-wave velocity structure[J]. Chinese Journal of Geophysics, 2013, 52(7):1923-1930.
[9]
苏悦, 杨振涛, 杨博, 等. 基于矢量波数变换法的主动源瑞雷波多模式提取方法在近地表地层结构探测中的应用研究[J]. 北京大学学报(自然科学版), 2020, 56(3):427-435.SU Yue, YANG Zhentao, YANG Bo, et al. Application research of active source Rayleigh wave multi-mode extraction method based on vector wavenumber transformation method in near surface stratigraphic structure detection[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2020, 56(3):427-435.
[10]
杨振涛, 陈晓非, 潘磊, 等. 基于矢量波数变换法(VWTM)的多道Rayleigh波分析方法[J]. 地球物理学报, 2019, 62(1):298-305.YANG Zhentao, CHEN Xiaofei, PAN Lei, et al. Multichannel analysis of Rayleigh waves based on the vector wavenumber transformation method (VWTM)[J]. Chinese Journal of Geophysics, 2019, 62(1):298-305.
[11]
罗银河, 夏江海, 刘江平, 等. 基阶与高阶瑞利波联合反演研究[J]. 地球物理学报, 2008, 51(1):242-249.LUO Yinhe, XIA Jianghai, LIU Jiangping, et al. Joint inversion of fundamental and higher mode Rayleigh waves[J]. Chinese Journal of Geophysics, 2008, 51(1):242-249.
[12]
夏江海, 高玲利, 潘雨迪, 等. 高频面波方法的若干新进展[J]. 地球物理学报, 2015, 58(8):2591-2605.XIA Jianghai, GAO Lingli, PAN Yudi, et al. New findings in high-frequency surface wave method[J]. Chinese Journal of Geophysics, 2015, 58(8):2591-2605.
[13]
WANG J, WU G, CHEN X. Frequency-Bessel transform method for effective imaging of higher-mode Rayleigh dispersion curves from ambient seismic noise data[J]. Journal of Geophysical Research:Solid Earth, 2019, 124(4):3708-3723.
[14]
李雪燕, 陈晓非, 杨振涛, 等. 城市微动高阶面波在浅层勘探中的应用:以苏州河地区为例[J]. 地球物理学报, 2020, 63(1):247-255.LI Xueyan, CHEN Xiaofei, YANG Zhentao, et al. Application of high-order surface waves in shallow exploration:An example of the Suzhou river, Shanghai[J]. Chinese Journal of Geophysics, 2020, 63(1):247-255.
[15]
林伯香, 孙晶梅, 刘清林. 层析成像低速带速度反演和静校正方法[J]. 石油物探, 2002, 41(2):136-140.LIN Boxiang, SUN Jingmei, LIU Qinglin. Tomography for LVZ velocity inversion and statics[J]. Geophysical Prospecting for Petroleum, 2002, 41(2):136-140.
[16]
林伯香, 孙晶梅, 徐颖, 等. 几种常用静校正方法的讨论[J]. 石油物探, 2006, 45(4):367-372.LIN Boxiang, SUN Jingmei, XU Ying, et al. Static correction approaches being frequently applied[J]. Geophysical Prospecting for Petroleum, 2006, 45(4):367-372.
[17]
李庆春, 邵广周, 刘金兰, 等. 瑞雷面波勘探的过去、现在和未来[J]. 地球科学与环境学报, 2006, 28(3):74-77.LI Qingchun, SHAO Guangzhou, LIU Jinglan, et al. Past, present and future of Rayleigh surface wave exploration[J]. Journal of Earth Sciences and Environment, 2006, 28(3):74-77.
[18]
CHEN X. Seismogram synthesis in multi-layered halfspace part Ⅰ. Theoretical formulations[J]. Earthquake Research in China, 1999, 13(2):149-174.
[19]
CHEN X. A systematic and efficient method of com-puting normal modes for multilayered half-space[J]. Geophysical Journal International, 1993, 115(2):391-409.
[20]
LUCO J E, APSEL R J. On the green's functions for a layered half-space. Part Ⅰ[J]. Bulletin of the Seismological Society of America, 1983, 73(4):909-929.
[21]
KENNETT B. Seismic wave propagation in stratified media[J]. Geophysical Journal of the Royal Astronomical Society, 1986, 86(1):219-220.
[22]
HISADA Y. An efficient method for computing green's functions for a layered half-space with sources and receivers at close depths[J]. Bulletin of the Seismological Society of America, 1994, 84(5):1456-1472.
[23]
SÁNCHEZ-SESMA F J, CAMPILLO M. Retrieval of the Green's function from cross correlation:the canonical elastic problem[J]. Bulletin of the Seismological Society of America, 2006, 96(3):1182-1191.
[24]
SATO H, FEHLER M C, MAEDA T. Green's Function retrieval from the cross-correlation function of Random waves//SATO H, FEHLER M C, MAEDA T. Seismic Wave Propagation and Scattering in the Heterogeneous Earth(2nd ed)[M]. Springer, Berlin, Heidelberg, 2012, 401-450.
[25]
张碧星, 肖柏勋, 杨文杰, 等. 瑞利波勘探中"之"形频散曲线的形成机理及反演研究[J]. 地球物理学报, 2000, 43(4):557-567.ZHANG Bixing, XIAO Baixun, YANG Wenjie, et al. Mechanism of zigzag dispersion curves in Rayleigh exploration and its inversion study[J]. Chinese Journal of Geophysics, 2000, 43(4):557-567.