地下空洞地震瑞雷波的旋转交错网格有限差分数值模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
将真空法和旋转交错网格相结合,实现了地下空洞高阶旋转交错网格波场数值模拟,克服了旋转交错网格在模拟自由边界条件时只在空间二阶差分时稳定的不足。通过和已有普通交错网格下改进的真空法对比,验证了该方法的正确性和稳定性。为了减弱模型人工边界产生的反射,采用了吸收效果较好的完全匹配层技术。模拟了不同模型参数组合情况下的地震波传播以及地震波与地下空洞的相互作用,通过对比不同参数组合的地震记录和波场快照,以及包含空洞和不含空洞地震记录和波场快照,分析了地震波在地下存在空洞时的传播特征以及震源子波频率、地下空洞大小和埋深对地表地震记录的影响。
We combined the improved vacuum formulation and the rotated staggered grid finite difference method to implement the high order rotated staggered finite difference elastic wave-field simulation at the presence of underground cavities.Our method overcame the defect of rotated staggered finite difference when modeling the free surface condition which is only stable at second order spatial FD operators,through the comparison with existing vacuum formulation testifies our methods accuracy and stability.The perfectly Matched Layer(PML)absorbing boundary condition is applied to decrease the reflection from the edges of the numerical model.Different seismograms are calculated with models without voids and models contain voids with different size and depth under different source wavelet frequency.The effects of different cavity size and depth and the source wavelet frequency are evaluated from the snapshots and seismograms with and without the cavity.
引文
[1]CHRISOTS ORFANOS,GEORGE APOSTOLOOULOS.Analysis of different geophysical methods in the detection of an underground opening at a controlled test site[J].Journal of The Balkan geophysical society,2012(15):7-18.
    [2]GILLES GRANGJEAN,DONATIENNE LEPAROUX.The potential of seismic methods for detecting cavities and buried objects:experimentation at a test site[J].Journal of Applied Geophysics,2004:93-106.
    [3]SAENGER E H,BOHLEN T.Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid[J].Geophysics,2004,69:583-591.
    [4]SAENGER E H,GOLD N.,Modeling the propagation of elastic waves using a modified finite difference grid[J].Wave Motion,2000,31(1):77-92.
    [5]GELIS C,D LEPAROUX,J.VIRIEUX,A.Bitri,S.Operto,G.Grandjean.Numerical modeling of surface waves over shallow cavities[J].Journal of Environmental and Engineering Geophysics,2005(10):111-121.
    [6]ZENG C J,XIA R D,MILLER,et al.An improved vacuum formulation for finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities[J].Geophysics,2012,77(1):T1-T9.
    [7]严红勇,刘洋.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟[J].地球物理学报,2012,55(4):1354-1365.YAN H Y,LIU Y.Rotated staggered grid high-order finite-difference numerical modeling for wave propagation in viscoelastic TTI media[J].Chinese J.Geophys,2012,55(4):1354-1365.(In Chinese)
    [8]COLLINO F,TSOGKA C.Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66:1749—1755.
    [9]陈浩,王秀明,赵海波.旋转交错网格有限差分及其完全匹配层吸收边界条件[J].科学通报,2006,51(17):1985-1994.CHEN H,WANG X M,ZHAO H B.Rotated staggered grid finite difference and the perfectly matched layer[J].Chinese Science Bulletin,2006,51(17):1985-1994.(In Chinse)
    [10]CHUNLEI CHU,PAUL L STOFFA.High-order rotated staggered finite difference modelling of 3Delastic wave propagation in general anisotropic media[R].SEG Houston International Expositon and Annual Meeting,2009.
    [11]CUNH C A.Elastic seismic modeling in discontinuous media[J].Geophysics,1993,59:1840-1851.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心