Numerical simulation and dispersion curve verification of three-dimensional background noises
HUO Keyu1, SHAO Guangzhou1, WANG Guoshun1, BAI Shuai1, WU Hua2
1. School of Geological Engineering and Geomatics, Chang'an University, Xi'an, Shaanxi 710054, China; 2. School of Science, Chang'an University, Xi'an, Shaanxi 710064, China
Abstract：In recent years, background noise imaging technology using vibrations generated by human production activities and natural earthquakes as the source has gradually become a research hotspot. Whether it is to verify the effectiveness of the background noise imaging method and the dispersion curve extraction method based on the known model, or to optimize the parameters of the background noise data acquisition method, it is necessary to obtain the required theoretical synthetic data of background noises through wave field numerical simulation. Three-dimensional (3D) numerical simulation is an important research content in background noise imaging. In view of the randomness of background noises in time and space, this paper introduces the random source into the finite difference numerical simulation of the 3D Rayleigh wave field to effectively calculate background noise data. The dispersion curve of the noise recorded by simulation is extracted by the two-trace method based on Aki formula and compared with the theoretical dispersion curve, so as to verify the feasibility and effectiveness of noise simulation. The test results of the theoretical model show that when the source parameters are within a reasonable range, the dispersion curve of the synthetic background noise data and the theoretical dispersion curve are well fitted;at the same time, the dispersion curve extracted from the 3D simulation data is more accurate than that by the two-dimensional simulation data.
OKADA H. Estimates of an S wave velocity distribution using long-period microtremors[J]. Geophysical Bulletin of Hokkaido University, 1983, 42:119-143.
[2]
TSAI V C, MOSCHETTI M P. An explicit relationship between time-domain noise correlation and spatial autocorrelation (SPAC) results[J]. Geophysical Journal International, 2010, 182(1):454-460.
[3]
EKSTRÖM G, ABERS G A, WEBB S C. Determination of surface-wave phase velocities across US array from noise and Aki's spectral formulation[J]. Geophysical Research Letters, 2009, 36(18):L18301.
[4]
刘军. 基于射线追踪的微地震模型多波场正演模拟[D]. 山东青岛:中国石油大学(华东), 2009.LIU Jun. Research on Multi-Wave Field Forward Simulation of Micro-Seismic Model Based on Ray Tracing[D]. China University of Petroleum (East China), Qingdao, Shandong, 2009.
[5]
容娇君, 张固澜, 郭晓玲, 等. 压裂/微地震事件地面响应信号模拟[J]. 石油天然气学报, 2010, 32(4):247-250, 432.RONG Jiaojun, ZHANG Gulan, GUO Xiaoling, et al. Surface response signal simulation of fracturing and microseismic events[J]. Journal of Oil and Gas Technology, 2010, 32(4):247-250, 432.
[6]
ZHAO X, YOUNG R P. Numerical modeling of seismicity induced by fluid injection in naturally fractured reservoirs[J]. Geophysics, 2011, 76(6):WC167WC180.
[7]
朱恒, 王德利, 时志安, 等. 地震干涉技术被动源地震成像[J]. 地球物理学进展, 2012, 27(2):496-502.ZHU Heng, WANG Deli, SHI Zhian, et al. Passive seismic imaging of seismic interferometry[J]. Progress in Geophysics, 2012, 27(2):496-502.
[8]
张唤兰. 微地震数值模拟及震源定位方法研究[D]. 陕西西安:西安科技大学, 2014.ZHANG Huanlan. Numerical Simulation of Microseismic and Study on Source Location Method[D]. Xi'an University of Science and Technology, Xi'an, Shaanxi, 2014.
[9]
唐杰, 方兵, 蓝阳, 等. 压裂诱发的微地震震源机制及信号传播特性[J]. 石油地球物理勘探, 2015, 50(4):643-649.TANG Jie, FANG Bing, LAN Yang, et al. Focal mechanism of micro-seismic induced by hydrofracture and its signal propagation characteristics[J]. Oil Geophysical Prospecting, 2015, 50(4):643-649.
[10]
徐宗博. 高频背景噪声波场模拟与面波成像[D]. 湖北武汉:中国地质大学(武汉), 2016.XU Zongbo. High-Frequency Ambient Seismic Noise Simulation and Surface-Wave Imaging[D]. China University of Geosciences(Wuhan), Wuhan, Hubei, 2016.
[11]
WANG F, XU G, QU X, et al. Numerical simulation of microseismic wave propagation for underground gas storage[C]. International Geophysical Conference, 2017, 827-830.
[12]
雷朝阳, 刘怀山. 被动源成像中透射与反射响应关系研究[J]. 中国海洋大学学报(自然科学版), 2017, 47(6):112-118.LEI Chaoyang, LIU Huaishan. Study on relation between transmission and reflection responses of passive seismic imaging[J]. Periodical of Ocean University of China(Nature Science), 2017, 47(6):112-118.
[13]
张晟瑞. 水力压裂微地震正演模拟分析及应用研究[D]. 黑龙江大庆:东北石油大学, 2018.ZHANG Shengrui. A Study on Analysis and Application of Hydraulic Fracturing Microseismic Forward Modeling[D]. Northeast Petroleum University, Daqing, Heilongjiang, 2018.
[14]
KNOPOFF L. A matrix method for elastic wave problems[J]. Bulletin of the Seismological Society of America, 1964, 54(1):431-438.
[15]
李远林. 被动源面波法在渭河盆地地层结构分层中的应用研究[D]. 陕西西安:长安大学, 2020.LI Yuanlin. Application of Passive Surface Wave Method in Stratigraphic Classification of Weihe Basin[D]. Chang'an University, Xi'an, Shaanxi, 2020.
[16]
MADARIAGA R. Dynamics of an expanding circular fault[J]. Bulletin of the Seismological Society of America, 1976, 66(3):639-666.
[17]
彭更新, 刘威, 郭念民, 等. 基于时空域交错网格有限差分法的应力速度声波方程数值模拟[J]. 石油物探, 2022, 61(1):156-165, 173.PENG Gengxin, LIU Wei, GUO Nianmin, et al. A time-space domain dispersion-relationship-based staggered-grid finite-difference scheme for modeling the stress-velocity acoustic wave equation[J]. Geophysical Prospecting for Petroleum, 2022, 61(1):156-165, 173.
[18]
唐超, 文晓涛, 王文化. 基于最小范数优化交错网格有限差分系数的波动方程数值模拟[J]. 石油地球物理勘探, 2021, 56(5):1039-1047.TANG Chao, WEN Xiaotao, WANG Wenhua. Numerical simulation of wave equations based on minimum-norm optimization of staggered-grid finite-difference coefficients[J]. Oil Geophysical Prospecting, 2021, 56(5):1039-1047.
[19]
全红娟, 朱光明, 王晋国, 等. 地震波震源加载方式对波场特征的影响[J]. 西北大学学报(自然科学版), 2012, 42(6):902-906.QUAN Hongjuan, ZHU Guangming, WANG Jinguo, et al. The influence of loading method of seismic source on wave field characteristic[J]. Journal of Northwest University (Natural Science Edition), 2012, 42(6):902-906.
[20]
COURANT R, FRIEDRICHS K, LEWY H. On the partial difference equations of mathematical physics[J]. IBM Journal of Research and Development, 1967, 11(2):215-234.
[21]
GRAVES R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences[J]. Bulletin of the Seismological Society of America, 1996, 86(4):1091-1106.
[22]
邵广周, 岳亮, 李远林, 等. 被动源瑞利波两道法提取频散曲线的质量控制方法[J]. 物探与化探, 2019, 43(6):1297-1308.SHAO Guangzhou, YUE Liang, LI Yuanlin, et al. A study of quality control of extracting dispersion curves by two-channel method of passive Rayleigh waves[J]. Geophysical and Geochemical Exploration, 2019, 43(6):1297-1308.
[23]
李欣欣. 主动源与被动源瑞利波联合成像方法研究[D]. 陕西西安:长安大学, 2017.LI Xinxin. Study on the Joint Tomographic Method for Active and Passive Source Rayleigh Wave Data[D]. Chang'an University, Xi'an, Shaanxi, 2017.
[24]
BENSEN G D, RITZWOLLER M H, BARMIN M P, et al. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements[J]. Geophysical Journal International, 2007, 169(3):1239-1260.
[25]
BOSCHI L, WEEMSTRA C, VERBEKE J, et al. On measuring surface wave phase velocity from station-station cross-correlation of ambient signal[J]. Geophysical Journal International, 2013, 192(1):346-358.
