起伏地表二维声波方程地震波场模拟与分析
|
摘要
为了模拟和分析起伏地表上的地震波场,本文采用一种坐标映射的方法将起伏地表映射成为水平表面,推导出映射后坐标系下的二维声波方程,然后在映射后坐标系下采用有限差分法对起伏地表地震波场进行数值模拟。计算结果表明该方法能够正确模拟起伏地表上的地震波场。通过对波形、能量、频谱、相位等方面的分析表明:起伏地表对地震波场的动力学特征影响很大。
In order to do simulation and analysis of seismic wavefield on relief surface, the paper used a coordinate mapping method to map relief surface onto horizontal surface and deduced 2-D acoustic wave equation on post-mapping coordinate, then, used finite-difference algorithm for numeric simulation of seismic wavefield on relief surface after coordinate mapping. The calculated results showed the method could correctly simulate the seismic wavefield on relief surface. It is shown by analyses in aspects of waveform, energy, spectrum and phase that the relief surface significantly affects the dynamic feature of seismic wavefield.
In order to do simulation and analysis of seismic wavefield on relief surface, the paper used a coordinate mapping method to map relief surface onto horizontal surface and deduced 2-D acoustic wave equation on post-mapping coordinate, then, used finite-difference algorithm for numeric simulation of seismic wavefield on relief surface after coordinate mapping. The calculated results showed the method could correctly simulate the seismic wavefield on relief surface. It is shown by analyses in aspects of waveform, energy, spectrum and phase that the relief surface significantly affects the dynamic feature of seismic wavefield.
引文
[1]Tessmer E,Kosloff D and Behle A.3-D elastic mod-eling with surface topography by a chebychev spectral method.Geophysics,1994,59(3):464~473
[2]Hesthol mSand Ruud B.2Dfinite-difference elastic wave modeling including surface topograph.Geophys-ical Prospecting,1994,42:371~390
[3]Nielsen F,Berg P and Skovgaard O.Using the pseudo spectral technique on curved grids for2Dacoustic for-ward modeling.Geophysical Prospecting,1994,42:321~341
[4]赵景霞,张叔伦,孙沛勇.曲网格伪谱法二维声波模拟.石油物探,2003,42(1):1~5
[5]Tessmer E,Kosloff Dand Behle A.E-lastic wave propagation si mulation in the presence of surface to-pography.Geophysics J Int,1992,108,621~632
[6]Hesthol m S and Ruud B.3-Dfinite-difference elastic wave modeling including surface topography.Geo-physics,1998,63(2):613~622
[7]Hesthol m S.Elastic wave modeling with free sur-faces:stability of long si mulations.Geophysics,2003,68(1):314~321
[8]张永刚.地震波场数值模拟方法.石油物探,2003,42(2):143~148
[9]崔兴福,张关泉.地震波方程人工边界的两种处理方法.石油物探,2003,42(4):452~455
[10]马永旺,詹正彬,顾汉明.伪谱法地震波传播数值模拟.石油物探,2001,40(2):42~48
[11]李文杰,魏修成,刘洋.声波正演中一种新的吸收边界条件——双重吸收边界条件.石油物探,2004,43(6):528~531
[12]熊小兵,贺振华.相位移加有限差分法波动方程正演模拟.石油物探,1998,37(3):22~28
[13]刘洋,李承楚,牟永光.任意偶数阶精度有限差分法正演模拟.石油地球物理勘探,1998,31(1):1~10
[14]王润秋等.地震波正演模拟的任意差分精细积分方法.石油地球物理勘探,2003,38(3):252~257
[15]朱生旺等.波动方程非规则网格任意阶精度差分法正演.石油地球物理勘探,2005,40(2):149~153
[2]Hesthol mSand Ruud B.2Dfinite-difference elastic wave modeling including surface topograph.Geophys-ical Prospecting,1994,42:371~390
[3]Nielsen F,Berg P and Skovgaard O.Using the pseudo spectral technique on curved grids for2Dacoustic for-ward modeling.Geophysical Prospecting,1994,42:321~341
[4]赵景霞,张叔伦,孙沛勇.曲网格伪谱法二维声波模拟.石油物探,2003,42(1):1~5
[5]Tessmer E,Kosloff Dand Behle A.E-lastic wave propagation si mulation in the presence of surface to-pography.Geophysics J Int,1992,108,621~632
[6]Hesthol m S and Ruud B.3-Dfinite-difference elastic wave modeling including surface topography.Geo-physics,1998,63(2):613~622
[7]Hesthol m S.Elastic wave modeling with free sur-faces:stability of long si mulations.Geophysics,2003,68(1):314~321
[8]张永刚.地震波场数值模拟方法.石油物探,2003,42(2):143~148
[9]崔兴福,张关泉.地震波方程人工边界的两种处理方法.石油物探,2003,42(4):452~455
[10]马永旺,詹正彬,顾汉明.伪谱法地震波传播数值模拟.石油物探,2001,40(2):42~48
[11]李文杰,魏修成,刘洋.声波正演中一种新的吸收边界条件——双重吸收边界条件.石油物探,2004,43(6):528~531
[12]熊小兵,贺振华.相位移加有限差分法波动方程正演模拟.石油物探,1998,37(3):22~28
[13]刘洋,李承楚,牟永光.任意偶数阶精度有限差分法正演模拟.石油地球物理勘探,1998,31(1):1~10
[14]王润秋等.地震波正演模拟的任意差分精细积分方法.石油地球物理勘探,2003,38(3):252~257
[15]朱生旺等.波动方程非规则网格任意阶精度差分法正演.石油地球物理勘探,2005,40(2):149~153
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心