复杂介质地震波传播模拟中边界元法与有限差分法的比较研究
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震波场数值模拟是理论地球物理学和勘探地球物理学的重要研究手段.在众多数值模拟方法中边界元法和有限差分法是两种典型的地震波传播模拟计算方法.边界元法是一种半解析-半数值的边界型方法,它显式地利用边界连续条件,沿着地层边界进行离散,具有降维、高精度和自动满足远场辐射条件的优点;有限差分法是一种典型的基于微分的区域型方法,它隐式地使用边界连续条件,以空间网格形式进行离散和数值逼近,具有高效、实用和容易数值实现的优点.本文以一个半圆形均匀Valley模型和两个非均匀断裂/断层模型为例,从计算精度、计算效率、频散特性以及适用性等方面对这两种方法进行了比较研究.数值计算结果表明:边界元法可以精确地几何描述有内部断点、断面的复杂构造,能够精确地模拟内部不规则界面之间波的反射/传播;有限差分法不能以足够的精度描述几何断点和内部不规则边界.边界元法在高频时计算量大于有限差分法,有限差分法则需要更小的网格间距以压制数值频散.因此,在处理内部非均质和高频计算时,有限差分法更有效;在处理内部不规则边界、断点、大尺度等问题时,边界元法比有限差分法更有优势.
Numerical modeling of seismic wave-field plays an important role in both theoretical geophysics and exploration geophysics.With the development of geophysics and numerical simulation technology,boundary element method(BEM) and finite-difference method(FDM) have been two basic methods in the study of seismic wave propagation.BEM is a semi-analytical and semi-numerical boundary type method,which explicitly uses the boundary continuity conditions for displacement and traction across interfaces,discretizes model along stratigraphic boundaries,and has advantages of spatial dimensionality reduction,high accuracy and easy fulfillment of radiation conditions at infinity.FDM is a typical domain type method,which is characteristic of the implicit use of boundary continuity conditions,discretizes model in the form of space grids,and has advantages of high efficiency and easy numerical implementation.To take full advantage of BEM and FDM,the numerical accuracy,computational efficiency,dispersion effect and applicability of the two methods are comprehensively compared in this paper.A homogeneous Valley model and two heterogeneous complex fault models are applied.Comparisons are implemented in frequency domain as well as in time domain.The results show that each method has specific merits and drawbacks.BEM is more accurate in the simulation of reflection/transmission across internal irregular interfaces,FDM cannot geometrically describe internal boundaries with enough accuracy.At high frequencies,BEM is computationally intensive since a considerable amount of matrix must be inverted,FDM is more efficient but it needs a much larger number of discretization grid points to decrease numerical dispersion.Therefore,FDM is a better choice when dealing with the problems of heterogeneity and high frequency computations,while BEM has potential advantages for solving fault point,internal irregular boundaries and large-scale problems.
引文
[1]Alterman Z,Karal F C.Propagation of seismic wave in layeredmedia by finite difference methods[J].Bull.Seism.Soc.Am.,1968,58(1):367-398.
    [2]Alford R M,Kelly K R,Boore D M.Accuracy of finite-difference modeling of the acoustic wave equation[J].Geophysics,1974,39:834-842.
    [3]Robertsson J O A,Blanch J O,Symes W W.Viscoelasticfinite difference modeling[J].Geophysics,1994,59:1444-1456.
    [4]Juhlin C.Finite-difference elastic wave propagation in 2Dheterogeneous transversely isotropic media[J].Geophys.Prosp.,1995,43:843-858.
    [5]Igel H,Mora P,et al.Anisotropic wave propagation throughfinite-difference grids[J].Geophysics,1995,60(4):1203-1216.
    [6]侯安宁,何樵登.各向异性介质中弹性波动高阶差分法及其稳定性的研究[J].地球物理学报,1995,38(2):243-251.Hou A N,He Q D.Study of an elastic wave high-orderdifference method and its stability in anisotropic media[J].Chinese J.Geophys.(in Chinese),1995,38(2):243-251.
    [7]Ozdenvar T,McMechan G A.Algorithms for staggered-gridcomputations for poroelastic,elastic,and scalar equations[J].Geophysical prospecting,1997,45:403-420.
    [8]VirieuxJ.SH-wave propagation in heterogeneous media:velocity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1957.
    [9]Graves R W.Simulating seismic wave propagation in 3Delasticmedia using staggered-grid finite-difference[J].Bull.Seism.Soc.Am.,1996,86:1091-1106.
    [10]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.Dong L G,Ma Z T,Cao J Z,et al.A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese J.Geophys.(in Chinese),2000,43(3):411-419.
    [11]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.Dong L G,Ma Z T,Cao J Z.A study on stability of thestaggered-grid high-order difference method of first-orderelastic wave equation[J].Chinese J.Geophys.(inChinese),2000,43(6):856-864.
    [12]裴正林,牟永光.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2003,27(6):17-21.Pei Z L,Mu Y G.A staggered-grid high-order differencemethod for modeling seismic wave propagation ininhomogeneous media[J].Journal of the University ofPetroleum,China(Edition of Natural Science)(in Chinese),2003,27(6):17-21.
    [13]程冰洁,李小凡,徐天吉.非均匀介质中交错网格高阶有限差分法数值模拟[J].物探化探计算技术,2006,28(4):294-298.Chen B J,Li X F,Xu T J.High-order difference modeling ofheterogeneous with staggered grid[J].ComputingTechniques for Geophysical and Geochemical Exploration(inChinese),2006,28(4):294-298.
    [14]陈可洋.标量声波波动方程高阶交错网格有限差分法[J].中国海上油气,2009,21(4):232-236.Chen K Y.High-order staggered-grid finite difference schemefor scalar acoustic wave equation[J].China Offshore Oil andGas(in Chinese),2009,21(4):232-236.
    [15]刘洪林,陈可洋,杨微,等.高阶交错网格有限差分法纵横波波场分离数值模拟[J].地球物理学进展,2010,25(3);877-884.Liu H L,Chen K Y,Yang W,et al.Numerical modeling ofP-and S-wave field separation with high-order staggered-gridfinite-difference scheme[J].Progress in Geophys.(inChinese),2010,25(3):877-884.
    [16]Saenger E H,Gold N,Shapiro S A.Modeling the propagation ofelastic waves using a modified finite-difference grid[J].WaveMotion,2000,31:77-92.
    [17]Saenger E H,Bohlen T.Finite-difference modeling of viscoelasticand anisotropic wave propagation using the rotated staggeredgrid[J].Geophysics,2004,69:583-591.
    [18]Jastram C,Tessemer E.Elastic modeling on a grid with verticallyvarying spacing[J].Geophysical prospecting,1994,42:357-370.
    [19]Oprsal I A,Zahradnik J.Elastic finite-difference method forirregular grids[J].Geophysics 1999,64:240-250.
    [20]Robertsson J O.A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence oftopography[J].Geophysics,1996,61:1921-1934.
    [21]Berenger J.A perfectly matched layer for the absorption ofsorption of electromagnetic wave[J].J ComputPhys,1994,114:185-200.
    [22]徐义,张剑锋.地震波数值模拟的非规则网格PML吸收边界[J].地球物理学报,2008,51(5):1520-1526.Xu Y,Zhang J F.An irregular-grid perfectly matched layerabsorbing boundary for seismic wave modeling[J].ChineseJ.Geophys.(in Chinese),2008,51(5):1520-1526.
    [23]Lysmer J,Drake L A.A finite element method for seismology[M],Methods in computational physics II,Seismolgy.Academic Press,1972,181-216.
    [24]Marfurt K J.Accuracy of finite-difference and finite-elementmodeling of the scalar and elastic wave equations[J].Geophysics,1984,49:533-549.
    [25]Seron F J,Sanz F J,Kindelan M,et al.Finite-element methodfor elastic wave propagation[J].Comm.Appl.NumericalMethods,1990,6(2):359-368.
    [26]Padovani E,Priolo E,Seriani.Low and high-order finiteelement method:experience in seismic modeling[J].J.Comp.Acoust.,1994,2(1):371-422.
    [27]Sarma G S,Mallick,Gadhinglajkar V R.Nonreflecting boundarycondition in finite element formulation for an elastic waveequation[J].Geophysics,1998,63(3):1006-1016.
    [28]周辉,徐世浙,刘斌,等.各向异性介质中波动方程有限元法模拟及其稳定性[J].地球物理学报,1997,40(6):833-841.Zhou H,Xu S Z,Liu B,et al.Modeling of wave equation inanisotropic media using finite element method and its stabilitycondition[J].Chinese J.Geophys.(Acta GeophysicaSinica)(in Chinese),1997,40(6):833-841.
    [29]张美根,王妙月,李小凡,等.各向异性弹性波场的有限元数值模拟[J].地球物理学进展,2002,17(3):384-389.Zhang M G,Wang M Y,Li X F,et al.Finite elementforward modeling of anisotropic elastic waves[J].Progressin Geophys.(in Chinese),2002,17(3):384-389.
    [30]杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报,2002,45(4):575-583.Yang D H.Finite element method of the elastic waveequation and wavefield simulation in two-phase anisotropicmedia[J].Chinese J.Geophys.(in Chinese),2002,45(4):575-583.
    [31]Brebbia C A.The boundary element method for engineer[M].New York:Wiley,1978:189.
    [32]Jaswon M A.Integral equation methods in potential theory I[J].Proc.Roy.Soc.Lond.,1963,275:23-32.
    [33]Kawase H.Time-domain response of a semi-circular canyonfor incident SV,P,and Rayleigh waves calculated by thediscrete wave number boundary element method[J].Bull.Seism.Soc.Am.,1988,78:1415-1437.
    [34]Bouchon M,Campillo M,Gaffet S.A boundary integral equation-discrete wave number representation method to study wavepropagation in multilayered media having irregular interfaces[J].Geophysics,1989,54:1134-1140.
    [35]Sanchez-Sesma F J et al.Diffraction of elastic waves by three-dimensional surface irregularities.Part II[J].Bull.Seism.Soc.Am.,1989,79:101-112.
    [36]符力耘,牟永光.弹性波边界元法正演模拟[J].地球物理学报,1994,37(4):521-529.Fu L Y,Mou Y G.Boundary element method for elasticwave forward modeling[J].Chinese J.Geophys.(inChinese),1994,37(4):521-529.
    [37]Fu L Y,Wu R S.A hybrid BE-GS method for modelingregional wave propagation[J].Pure and AppliedGeophysics,2001,158:1251-1277.
    [38]Fu L Y.Seismogram synthesis for piecewise heterogeneousmedia[J].Geophys.J.Int.,2002,150:800-808.
    [39]Fu L Y.Numerical study of generalized Lipmann-Schwingerintegral equation including surface topography[J].Geophysics,2003,68(2):665-671.
    [40]Fu L Y,Bouchon M.Discrete wavenumber solutions to numericalwave propagation in piecewise media heterogeneous media-ITheory of two-dimensional SH case[J].Geophys.J.Int.,2004,157:481-498.
    [41]Fu L Y,Wu R S.Infinite boundary element absorbing boundaryfor wave propagation simulations[J].Geophysics,2004,65:625-637.
    [42]Fu L Y.Rough surface scattering:Comparison of variousapproximation theories for 2DSH waves[J].Bull.Seism.Soc.Am.,2005,95:646-663.
    [43]Yu G X,Fu L Y,Yao Z X.Comparison of different BEM+Born series modeling schemes for wave propagation incomplex geological structures[J].Geophysics,2010,75:71-82.
    [44]Fujiwara H.The fast multipole method for solving integralequations of three-dimensional topography and basin problems[J].Geophysical Journal International,2000,140:198-210.
    [45]Grasso E,Chaillat S,Bonnet M,Semblat JF.Application ofthe multi-level time-harmonic fast multipole BEM to 3-Dvisco-elastodynamics[J].Engineering Analysis withBoundary Elements,2012,36:744-758.
    [46]Yazdchi M,Khalili N,Valliappan S.Non-linear seismic behaviorof concrete gravity dams using coupled finite-boundaryelement technique[J].International Journal for NumericalMethods in Engineering,1999,44(6):101-130.
    [47]Gatmiri B,Arson C,Nguyen K V.Seismic site effects by anoptimized 2D BE/FE method I.Theory,numericaloptimization and application to topographical irregularities[J].Soil Dyn.And Earthq.Eng.,2008,28(8):632-645.
    [48]Ohtsuki A,Harumi K.Effects of topography and subsurfaceinhomogeneities on seismic SV waves[J].Earthquake EngStruct Dyn,1983,11:441-462.
    [49]Bielak J,Xu J F.Earthquake ground motion and structuralresponse in alluvial valleys[J].J Geotechn Geoenviron Eng,1999,125(5):413-423.
    [50]Pedensen H A,Campillo M,Sanchez-Sesma F J.Azimuthdependent wave amplification in alluvial valleys[J].SoilDyn.And Earthq.Eng.,1995,14(4):289-300.
    [51]Barros F C P,Luco J E.Amplification of obliquely incidentwaves by a cylindrical valley embedded in a layered half-space[J].Earthquake Eng Struct Dyn,1995,24:163-175.
    [52]Pao Y H,Varatharajulu V.Huygens’principle,radiationconditions,and integral formulas for the scattering of elasticwaves[J].Journal of the Acoustical Society of America,1976,59(6):1361-1371.
    [53]Trifunac M D.Surface motion of a semicylindrical alluvialvalley for incident plane SH waves.Bull.Seism.Soc.Arrrer.,1971,61:1755-1770.
    [54]Ge Z,Chen X F.Wave propagation in irregularly layeredelastic models:a boundary element approach with a globalreflection/transmission matrix propagator[J].Bull.Seism.Soc.Am.,2007,97:1025-1031.
    [55]曹军,盖增喜,张坚,陈晓非.用于不规则界面多层介质的地震波方法比较研究:不规则地形问题[J].地球物理学报,2004,47(3):495-503.Cao J,Ge Z X,Zhang J,et al.A comparative study onseismic wave methods for multi-layered media with irregularinterfaces Part 1:irregular topography problem[J].ChineseJ.Geophys.(in Chinese),2004,47(3):495-503.
    [56]Marcel Frehner,Stefan M.Schmalholz,Erik H.Saenger,Holger Steeb.Comparison of finite difference and finiteelement methods for simulating two-dimensional scattering ofelastic waves[J].Physics of the Earth and PlanetaryInteriors,2008,171:112-121.
    [57]皮红梅,蒋先芝,刘财,等.波动方程数值模拟的三种方法及对比[J].地球物理学进展,2009,24(2):391-397.Pi H M,Jiang X Z,Liu C,et al.Three methods of waveequation forward modeling and comparisions[J].Progress inGeophys.(in Chinese),2009,24(2):391-397.
    [58]Liu Yang,Mrinal K.Sen.A new time-space domain high-order finite-difference method for the acoustic wave equation[J].Journal of Computational Physics,2009,228:8779-8806.
    [59]冯英杰,杨长春,吴萍.地震波有限差分模拟综述[J].地球物理学进展,2007,22(2):487-491.Feng Y J,Yang C C,Wu P.The review of the finite-difference elastic wave motion modeling[J].Progress inGeophys.(in Chinese),2007,22(2):487-491.
    [60]殷文,印兴耀,吴国忱,等.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报,2006,49(2):561-568.Yin W,Yin X Y,Wu G C,et al.The method of finitedifference of high precision elastic wave equations in thefrequency domain and wavefield simulation[J].Chinese J.Geophys.(in Chinese),2006,49(2):561-568.

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