摘要
瑞雷波勘探技术具有快速、无损、衰减小、抗干扰能力强等优点,广泛应用于实际的工程勘探。在实际的工程勘探中提取的瑞雷波频散曲线只有一条,是正演中多个模式频散曲线的拟合,但是在正演中无法直接拟合成一条曲线。因此通过计算研究瑞雷波的能量分布,能够直接计算出各个模式频散曲线的能量相对值。通过与有限差分模拟法对比,实验表明:该算法能够实现不同模式频散曲线的能量相对值的计算,为频散曲线拟合提供了直接的数据基础。
Rayleigh wave exploration technology has the advantages of fast, lossless, low attenuation, strong anti-interference ability, and is widely used in practical engineering exploration. There is only one dispersion curve of Rayleigh wave in the actual engineering exploration, which is the fitting of the multiple modes dispersion curve in the forward process. But it can not be directly fitted into a curve in forward modeling. Therefore, by studying the calculation of Rayleigh wave's energy distribution, we can directly calculate the energy's relative value of each mode dispersion curve. Compared with the finite difference method, the experimental results show that the algorithm can calculate the different modes' relative value of dispersion curve, which provides a direct data foundation for the dispersion curve fitting.
引文
[1] Socco L V, Foti S, Boiero D. Surface-wave analysis for build-ing near-surface velocity models—Established approaches andnew perspectives[J]. Geophysics, 2010, 75(5):75A83-75A102.
[2] Zhang S, Li X, Jeong H, et al. Modeling nonlinear Rayleighwave fields generated by angle beam wedge transducers—Atheoretical study[J]. Wave Motion, 2016, 67:141-159.
[3] Bao X, Dalton C A, Ritsema J. Effects of elastic focusing onglobal models of Rayleigh wave attenuation[J]. GeophysicalJournal International, 2016, 207(2):1062-1079.
[4] Ma S, Audet P. Seismic velocity model of the crust in thenorthern Canadian Cordillera from Rayleigh wave dispersiondata[J]. Canadian Journal of Earth Sciences, 2017, 54(2):163-172.
[5] MaranòS, Hobiger M, F?h D. Retrieval of Rayleigh wave ellip-ticity from ambient vibration recordings[J]. Geophysical Journal International, 2017, 209(1):334-352.
[6] Palomeras I, Villase?or A, Thurner S, et al. Lithospheric struc-ture of Iberia and Morocco using finite-frequency Rayleighwave tomography from earthquakes and seismic ambient noise[J]. Geochemistry, Geophysics, Geosystems, 2017, 18(5):1824-1840.
[7] Beilina L. Domain decomposition finite element/finite differ-ence method for the conductivity reconstruction in a hyperbol-ic equation[J]. Communications in Nonlinear Science and Nu-merical Simulation, 2016, 37:222-237.
[8] Xu Y, Xia J, Miller R D. Numerical investigation of implemen-tation of air-earth boundary by acoustic-elastic boundary ap-proach[J].Geophysics,2007,72(5):SM147-SM153.
[9]辛维,闫子超,梁文全,等.用于弹性波方程数值模拟的有限差分系数确定方法[J].地球物理学报, 2015, 58(7):2486-2495.
[10] Pal J, Ghorai A P. Comparison and Analysis of SurfaceWaves Propagation in Initially Stressed Dry Sandy Layer Us-ing Conventional and Time Dependent Finite DifferenceScheme[J]. International Journal of Applied Engineering Re-search, 2016, 11(10):7119-7124.
[11]刘建宙.有限差分法在瑞雷波波场正演模拟中的应用[D].中国地质大学(北京), 2014.