复杂模型高精度数值模拟研究
|
摘要
通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。
Second-order finite difference and fourth-order finite difference are used to the numerical modeling of wave equation with. the hybrid boundary condition consisted of both the absorbing boundary and the nonreflecting boundary. The computed results indicated the precision of fourth-order finite difference is obviously higher than that of second-order difference and the boundary reflection may be absorbed completely. For the difficulty to simulate the complex model regarding the single seismic record, the single stack section and the blast reflecting surface zero offset section are presented, and the processing of section migration is realized. The seismic section reflects complex structure information more clearly. The seismic modeling shows the satisfactory performance in the case of the complex model.
Second-order finite difference and fourth-order finite difference are used to the numerical modeling of wave equation with. the hybrid boundary condition consisted of both the absorbing boundary and the nonreflecting boundary. The computed results indicated the precision of fourth-order finite difference is obviously higher than that of second-order difference and the boundary reflection may be absorbed completely. For the difficulty to simulate the complex model regarding the single seismic record, the single stack section and the blast reflecting surface zero offset section are presented, and the processing of section migration is realized. The seismic section reflects complex structure information more clearly. The seismic modeling shows the satisfactory performance in the case of the complex model.
引文
[1]贺振华,熊晓军.等时叠加波动方程叠加正演[J].物探化探计算技术,2005,27(3):194.
[2]邵治龙,贺振华,黄德济.正交各向异性介质中地震记录的有限差分数值模拟[J].物探化探计算技术,1998,20(4):300.
[3]张永刚.地震波场数值模拟方法[J].石油物探,2003,42(2):143.
[4]KELLY.Synthetic seismograms-A finited difference approach[J].Geoph,1976,41(2):27.
[5]RENOLDS A C.Boundary Conditions for the Numerical Solution of wave propagation Problems[J].Geoph,1978,43(6):1099.
[6]CLAYTON R,ENGQUISTB.Absorbing boundary Con-ditions for wave-equation migration[J].Geoph,1980,45(5):895.
[7]WANG RUNQIU,JIAXIAOFENG,HUTIANYUE.The Precise Finite Difference Method for seismic Modeling[J].Pplied Geoph,2004,1(2):69.
[8]AHMAD ZAKARIAL,JOHN PENROSE,FRANK THO-MASL,et al.The Two Dimensional Numerical Modeling Of Acoustic Wave Propagation in ShallowWater[J].A-coustical Society Conference,Joodalup,Australia,2006,6:15.
[9]常旭,刘伊克.地震正反演与成像[M].北京:华文出版社,2001.
[10]张关泉,魏素花.声波方程吸收边界条件的稳定性分析[J].计算数学,1998,20(1):103.
[2]邵治龙,贺振华,黄德济.正交各向异性介质中地震记录的有限差分数值模拟[J].物探化探计算技术,1998,20(4):300.
[3]张永刚.地震波场数值模拟方法[J].石油物探,2003,42(2):143.
[4]KELLY.Synthetic seismograms-A finited difference approach[J].Geoph,1976,41(2):27.
[5]RENOLDS A C.Boundary Conditions for the Numerical Solution of wave propagation Problems[J].Geoph,1978,43(6):1099.
[6]CLAYTON R,ENGQUISTB.Absorbing boundary Con-ditions for wave-equation migration[J].Geoph,1980,45(5):895.
[7]WANG RUNQIU,JIAXIAOFENG,HUTIANYUE.The Precise Finite Difference Method for seismic Modeling[J].Pplied Geoph,2004,1(2):69.
[8]AHMAD ZAKARIAL,JOHN PENROSE,FRANK THO-MASL,et al.The Two Dimensional Numerical Modeling Of Acoustic Wave Propagation in ShallowWater[J].A-coustical Society Conference,Joodalup,Australia,2006,6:15.
[9]常旭,刘伊克.地震正反演与成像[M].北京:华文出版社,2001.
[10]张关泉,魏素花.声波方程吸收边界条件的稳定性分析[J].计算数学,1998,20(1):103.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心