具有不规则自由表面的粘弹性介质中声波模拟的有限差分算法
|
摘要
声波在具有不规则自由表面的粘弹性介质传播的数值模拟工作在地震勘探,地震预测中非常重要,特别是当模型有起伏的自由界面和较强的衰减特性时更是如此。基于镜像方法和直接方法,本文发展了一种二维有限差分算法,可以模拟不规则自由表面引起的声波散射问题。该方法将自由表面条件与速度应力方程结合求解粘弹性波动方程,在垂直和水平自由段及拐点处设置相应的边界条件。该方法假设自由表面穿过剪切应力和相应参数的网格点。为了提高计算精度,分别计算了水平和垂直方向上应力镜像值,而对质点速度,采用先水平方向后垂直方向分析进行镜像计算和更新。将粘弹性水平自由表面镜像方法和不同倾斜度平滑自由表面的改进方法的数值结果进行对比,并验证了算法的精度。
Modelling acoustic wave propagation in heterogeneous viscoelastic media with free-surfaces plays an important role in seismic exploration,earthquake prediction,especially for the model with a rugged free surface and strong attenuation.Based on imaging and direct methods,we have developed a new improved two-dimensional finite-difference algorithm for simulating acoustic wave scattering from an irregular free surface.The free surface conditions are combined with a velocity-stress formulation for full viscoelastic wave equations.The boundary conditions are implemented in the vertical and horizontal free surface segments and at their corners.The method assumes that the free surface passes through the grid points of shear stress and the corresponding shear memory variables.In order to increase accuracy,imaging is performed for stress components in both horizontal and vertical directions independently.To update particle velocities,imaging and updating are first carried out in horizontal direction,and then,in vertical direction.Numerical results are compared for a viscoelastic flat horizontal free surface with the imaging method and those for flat free surface of various slope angles with the improved method,and are shown to be in good agreement.
Modelling acoustic wave propagation in heterogeneous viscoelastic media with free-surfaces plays an important role in seismic exploration,earthquake prediction,especially for the model with a rugged free surface and strong attenuation.Based on imaging and direct methods,we have developed a new improved two-dimensional finite-difference algorithm for simulating acoustic wave scattering from an irregular free surface.The free surface conditions are combined with a velocity-stress formulation for full viscoelastic wave equations.The boundary conditions are implemented in the vertical and horizontal free surface segments and at their corners.The method assumes that the free surface passes through the grid points of shear stress and the corresponding shear memory variables.In order to increase accuracy,imaging is performed for stress components in both horizontal and vertical directions independently.To update particle velocities,imaging and updating are first carried out in horizontal direction,and then,in vertical direction.Numerical results are compared for a viscoelastic flat horizontal free surface with the imaging method and those for flat free surface of various slope angles with the improved method,and are shown to be in good agreement.
引文
[1]TESSMER E,KOSLOFF D,BEHLE A.Elastic wave propagation simulation in the presence of surface topography[J].Geophys.J.Internat.,1992,108(2):621-632.
[2]HESTHOLM S O,RUUD B O.2D finite-difference viscoelastic wave modeling including surface topography[J].Geophysical Prospecting,2000,48(2):341-373.
[3]ROBERTSSON J O A.A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J].Geophysics,1996,61(6):1921-1934.
[4]KOMATITSCH D,VILOTTE J P.The spectral element method:an efficient tool to simulate the seismic response of2D and3D geological structures[J].Bull.Seism,Soc.Am.,1998,88(2):368-392.
[5]KISHORE N N,SRIDHAR I,IYENGAR N G R.Finite element modelling of the scattering of ultrasonic waves by isolated flaws[J].NDT&E International,2000,33(5),297-305.
[6]DAUKSHER W,EMERY A F.Accuracy in modeling the acoustic wave equation with Chebyshev spectral finite elements[J].Finite Elements in Analysis and Design,1997,26(2):115-128.
[7]SERIANI G.3-D large-scale wave propagation modelling by spectral element method on Cray T3E[J].Comp.Meth.Appl.Mech.and Eng.,1998,164(1-2):235-247.
[8]TAL-EZER H.Spectral methods in time for hyperbolic problems[J].SIAM J.Numer.Anal.,1986,23(1):11-26.
[9]LEVANDER A R.Fourth-order finite-difference P-SV seismograms[J].Geophysics,1988,53(11):1425-1436.
[10]GRAVES R.Simulation seismic wave propagation in3D elastic media using staggered-grid finite differences[J].Bull.Seism.Soc.Am.,1996,86(4):1091-1106.
[11]WANG X M,ZHANG H L.Modeling of elastic wave propagation on a curved free surface using an improved finite-difference algorithm[J].Science in China Series G.Physics Mechanics and Astronomy,2004,47(5):633-648.
[12]ROBERTSSON J O A,BLANCH J O,SYMES W W.Viscoelastic finite-difference modeling[J].Geophysics,1994,59(9):1444-1456.
[2]HESTHOLM S O,RUUD B O.2D finite-difference viscoelastic wave modeling including surface topography[J].Geophysical Prospecting,2000,48(2):341-373.
[3]ROBERTSSON J O A.A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J].Geophysics,1996,61(6):1921-1934.
[4]KOMATITSCH D,VILOTTE J P.The spectral element method:an efficient tool to simulate the seismic response of2D and3D geological structures[J].Bull.Seism,Soc.Am.,1998,88(2):368-392.
[5]KISHORE N N,SRIDHAR I,IYENGAR N G R.Finite element modelling of the scattering of ultrasonic waves by isolated flaws[J].NDT&E International,2000,33(5),297-305.
[6]DAUKSHER W,EMERY A F.Accuracy in modeling the acoustic wave equation with Chebyshev spectral finite elements[J].Finite Elements in Analysis and Design,1997,26(2):115-128.
[7]SERIANI G.3-D large-scale wave propagation modelling by spectral element method on Cray T3E[J].Comp.Meth.Appl.Mech.and Eng.,1998,164(1-2):235-247.
[8]TAL-EZER H.Spectral methods in time for hyperbolic problems[J].SIAM J.Numer.Anal.,1986,23(1):11-26.
[9]LEVANDER A R.Fourth-order finite-difference P-SV seismograms[J].Geophysics,1988,53(11):1425-1436.
[10]GRAVES R.Simulation seismic wave propagation in3D elastic media using staggered-grid finite differences[J].Bull.Seism.Soc.Am.,1996,86(4):1091-1106.
[11]WANG X M,ZHANG H L.Modeling of elastic wave propagation on a curved free surface using an improved finite-difference algorithm[J].Science in China Series G.Physics Mechanics and Astronomy,2004,47(5):633-648.
[12]ROBERTSSON J O A,BLANCH J O,SYMES W W.Viscoelastic finite-difference modeling[J].Geophysics,1994,59(9):1444-1456.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心