准P波方程紧致交错网格井间地震波场模拟及边界条件
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义。基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程。重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟。数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果。一阶准P波方程能以足够的精度描述VTI介质中准P波特征。完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法。
Study on the build-up process and the spreading mechanism of cross-hole seismic wave field is very important for guiding actual cross well seismic exploration.Based on the first-order quasi-P wave equation in transversely isotropic media with a vertical symmetry axis(VTI),the finite difference scheme of the equation which was solved in the compact staggered grid was constructed.The compact staggered grid is implicit difference scheme which has the characteristic of unconditional stability.The perfectly matched layer(PML) absorbing boundary condition which fit the equation were thoroughly studied.Based on the studies above the forward modeling was realized for the first-order quasi-P wave equation in VTI medium.In summary,compact staggered-grid finite difference scheme could have a high degree of accuracy to simulate the spreading process of quasi-P wave in VTI medium.The first-order quasi-P wave equation could describe the characteristic of quasi-P wave in VTI medium with a high accuracy.PML is one of the best absorption boundary conditions which can effectively deal with the problem of article boundary.
引文
[1]周建宇,何惺华,李安夏,等.罗家地区井间地震方法与效果[J].石油地球物理勘探,2001,36(6):745.
    [2]孔庆丰,王延光,左建军,等.井间地震有限角度叠加方法研究与应用[J].石油地球物理勘探,2007,42(3):256.
    [3]李万万,裴正林.井间地震弹性波传播特征数值模拟[J].物探与化探,2008,32(2):207.
    [4]ZHANG J F.VERSCHUUR D J.Elastic wave prop-agation in heterogeneous anisotropic media using thelumped finite-element method[J].Geophysics,2002,67(2):625.
    [5]祝贺君,张伟,陈晓非.二维各向异性介质中地震波场的高阶同位网格有限差分模拟[J].地球物理学报,2009,52(6):1536.
    [6]汪利民,徐义贤,江贵.井中和井间地震波场正演模拟[J].石油物探,2009,48(2):146.
    [7]韩令贺,何兵寿.VTI介质一阶准P波方程正演模拟及边界条件[J].石油地球物理勘探,2010,45(6):819.
    [8]ALKHALIFAH T.An acoustic wave equation foranisotropic media[J].Geophysics,2000,65(4):1239.
    [9]HESTHOLM S.Acoustic VTI modeling using high-order finite-differences[J].Geophysics,2009,74(5):T67.
    [10]MADARIAGA R.Dynamics of an expanding circularfault[J].BSSA,1976,66(3):639.
    [11]LELE S K.Compact finite difference schemes withspectral-like resolution[J].J.Comput Phys,1992,103(1):16.
    [12]王书强,杨顶辉,杨宽德.弹性波方程的紧致差分方法[J].清华大学学报:自然科学版,2002,42(8):1128.
    [13]马廷福,金涛,曹福军,等.三维双调和方程的高阶紧致差分格式及其多重网格方法[J].兰州理工大学学报,2010,36(3):142.
    [14]葛永斌,刘国涛.三维Helmholtz方程的高阶隐式紧致差分方法[J].工程数学学报,2010,27(5):853.
    [15]NAGARAJAN S.LELE S K,FERZIGER J H.Arobust high-order compact method for large eddysimulation[J].Comput.Phys,2003,191(2):392.
    [16]BENDIKS.A staggered compact finite difference for-mulation for the compressible Navier-Stokes equa-tions[J].Comput.Phys,2005,208(2):675.
    [17]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411.
    [18]DU Q Z,LIN B,HOU B.Numerical modeling ofseismic wavefields in transversely isotropic mediawith a compact staggered-grid finite differencescheme[J].Applied Geophysical,2009,6(1):42.
    [19]COLLINO F,TSOGKA C.Application of the per-fectly matched absorbing layer model to the linearelastodynamic problem in an isotropic heterogeneousmedia[J].Geophysics,2001,66(1):294.
    [20]何兵寿,张会星.VTI介质中准P波方程的数值解法[J].煤炭学报,2006,31(4):446.
    [21]TARIQ A.An acoustic wave equation for anisotropicmedia[J].Geophysics,2000,65(4):1239.
    [22]吴律.层析基础及其在井间地震中的应用[M].北京:石油工业出版社,1997.

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