双相介质地震波场数值模拟的迭积微分算子及其PML边界条件
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
将基于计算数学中F0rsyte广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度——应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了合流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.
The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation.To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media.To attenuate the reflection artifacts caused by the artificial boundary, the Perfectly Matched Layer (PML) absorbing boundary is also considered in the method to deal with the boundary problem due to its advantage of automatically handling large-angle emission.This paper constructs the constitutive relationship for two-phase media, and further derives the first-order velocity-stress equation for 2D two-phase media.Numerical modeling using the CFPD method is carried out in the above-mentioned media.The results modeled in Biot two-phase media can better explain the liquid pore characteristics and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.
引文
[1] Biot M A Theory of propagation of elastic waves in a fluidsaturated porous solid.Ⅰ.Lower-frequency range, Ⅱ.Higher-frequency range[J].J.Acoust.Soc.Am., 1956, 28:179~191.
    [2] Zhu X, McMechan G A.Numerical simulation of seismic response of poroelastic reservoirs using Biot theory.Geophysics, 1991 , 56 (3) :328~339.
    [3] 刘洋,李承楚.双相各向异性介质中弹性波传播伪谱法数值模拟研究[J].地震学报,2000,22:132~138. Liu Y, Li C C.Study of elastic wave propagation in two phase anisotropic media by numerical modeling of pseudospectral method[J].Acta Seismologica Sinica ,2000 , 22 : 132 ~138.
    [4] 杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报,2002,45(2) :575~583. Yang D H.Finite element method of the elastic wave equation and wavefield simulation in two-phase anisotropic media[J].Chinese J.Geophys.(in Chinese), 2002,45(2) :575~583.
    [5] 王秀明,张海澜,王东.利用高阶交错网格有限差分算法模拟地震波在非均匀孔隙介质中的传播[J].地球物理学报,2003,46(6) :842~849. Wang X M, Zhang H L , Wang D.Modeling of seismic wave propagation in heterogeneous poroelastic media using a high-staggered finite-difference method[J].Chinese J.Geophys.(in Chinese) ,2003,46(6) :842~849.
    [6] 裴正林.三维双相各向异性介质中弹性波方程交错网格高阶有限差分法模拟[J].石油大学学报:自然科学版,2004,28(5) :23~29.Pei Z L.A staggered-grid high-order difference method for modeling seismic wave propagation in 3D dual-phase anisotropic media[J].Journal of the University of Petroleum,China(Edition of Natural Science), 2004 ,28(5) :23~29.
    [7] 王东,张海澜,王秀明.部分饱和孔隙岩石中声波传播数值研究[J].地球物理学报,2006,49(2) :524~532. Wang D, Zhang H L, Wang X M.A numerical study of acoustic wave propagation in partially saturated poroelastic rock[J].Chinese J.Geophys.(in Chinese), 2006, 49 (2) :524~532.
    [8] 卢明辉,巴晶,杨慧珠.双相介质分界面上弹性波的反射与透射[J].地球物理学进展,2007,22(5) :1439~1445. Lu M H, Ba J, Yang H Z.Reflection and transmission of elastic waves from a boundary of two-phase media[J].Progress in Geophysics (in Chinese), 2007, 22 (5) : 1438~1445.
    [9] 巴晶,卢明辉,胡彬,等.黏弹双相介质中的松弛骨架模型[J].地球物理学报,2008,51(5) :1527~1537. Ba J, Lu M H, Hu B, et al.The skeleton-relaxed model for poroviscoelastic media[J].Chinese J.Geophys.(in Chinese),2008, 51(5) :1527~1537.
    [10] 孟庆生,何樵登,樊玉清,等.基于BISQ机制三维双相EDA介质中的反射与透射[J].地球物理学进展,2005,20(2) :323~327. Meng Q S, He Q D ,Fan Y Q , et al.The reflection and transmission in 3D double-phase EDA media based on BISQ mechanisms[J].Progress in Geophysics (in Chinese), 2005,20(2) :323~327.
    [11] Jose M, Carcione, Gerard C, Herman, etal., Review Article:Seismic modeling[J].Geophysics, 2002,67 ( 4 ) : 1304~1325.
    [12] 郑洪伟,李廷栋,高锐,等.数值模拟在地球动力学中的研究进展[J].地球物理学进展,2006,21(2) :360~369. Zheng H W, Li T D, Gao R, et al.The advance of numerical simulation in geodynamics[J].Progress in Geophysics(in Chinese), 2006,21(2) :360~369.
    [13] Zhou B, Greenhalgh S A.Seismic scalar wave equation modeling by a convolutional diffentiator[J].Bulletin of Seismological Society of America,1982(1) :289~303
    [14] 张中杰,滕吉文,杨顶辉.声波与弹性波场数值模拟中的迭积微分算子法[J].地震学报,1996,18(1) :63~69.Zhang Z J, Teng J W, Yang D H.Acoustic and elastic wave modeling by a convolutional differentiator[J].Acta Seismoligica Sinica, 1996,18(1) : 63~69.
    [15] 戴志阳,孙建国,查显杰.地震波场模拟中的褶积微分算子法[J].吉林大学学报(地球科学版),2005,35(4) :520~524. Dai Z Y, Sun J G, Li X J.Seismic wave field modeling with convolution differentiator algorithm[J].Journal of JiLin University (Earth Science Edition).2005,35(4) : 520~524.
    [16] 戴志阳,孙建国,查显杰.地震波混合阶迭积算法模拟[J].物探化探计算技术,2005,27(2) :111~114. Dai Z Y, Sun J G, Li X J.Seismic wave field modeling with hybrid order convolution differentiator[J].Computing Techniques for Geophysics and Geochemical Exploration,2005, 27(2) :111~114.
    [17] 谢靖.物探数据处理的数学方法[M].北京:地质出版社,1981:97~104. Xie J.Numerical method for geophysical data-process[M].Beijing: Geological Publishing House, 1981 : 97 ~ 104.
    [18] 程冰洁,李小凡,龙桂华.基于广义正交多项式褶积微分算子 的地震波场数值模拟方法[J].地球物理学报,2008,51(2) :531~537.Cheng B J,Li X F, Long G H.Seismic waves modeling by convolutional Forsyte polynomial differentiator method[J].Chinese J.Geophys.(in Chinese), 2008,51(2) :531~537.
    [19] 程冰洁,李小凡.2. 5维地震波场褶积微分算子法数值模拟[J].地球物理学进展,2008,23(4) :1099~1105. Cheng B J ,Li X F.2. 5D modelling of elastic waves using the convolutional differentiator method[J].Progress in Geophysics (in Chinese), 2008,23(4) :1099~1105.
    [20] Clayton R, Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations[J].Bulletin of Seismological Society of America, 1997, 67: 1529~1540.
    [21] Cerjan C, Kosloff D, et al.A non-reflection boundary condition for discrete acoustic and elastic wave equation[J].Geophysics, 1985,50:705~708.
    [22] Berenger J P.Three-dimensional perfectly matched layer for the absorption electromagnetic waves[J].Comp.Phys, 1996,127:363~379.
    [23] Francis C, Chrysoula T.Application of the perfectly matched absorbing lamodel to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics, 2001, 66(1) :294~307.
    [24] 赵海波,王秀明,王东,陈浩.完全匹配层吸收边界在孔隙介质弹性波模拟中的应用[J].地球物理学报,2007,50(2) :581~591. Zhao H B, Wang X M, Wang D, Chen H.Applications of the boundary absorption using a perfectly matched layer for elastic wave simulation in poroelastic media[J].Chinese J.Geophys.(in Chinese), 2007, 50(2) :581~591.
    [25] 徐义,张剑锋.地震波数值模拟的非规则网格PML吸收边界[J].地球物理学报,2008,51(5) :1520~1526. Xu Y, Zhang J F.An irregular-grid perfectly matched layer absorbing boundary for seismic wave modeling[J].Chinese J.Geophys.(in Chinese), 2008, 51(5) :1520~1526.
    [26] Zeng Y Q,He J Q,Liu Q H.The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics,2001,66(4) :1258~1266.

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