塔里木地区勘探地震正演模拟研究
|
摘要
以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过串行程序并行化,大大减少了三维正演模拟的耗时,完成了大数据量三维正演模拟.
On the basis of the geological background of Tarim foreland basin,this paper realizes parallel three-dimensional forward modeling in Kuqa area by using three-dimensional arbitrary difference precise integration algorithm.Compared with the conventional differential method,the three-dimensional arbitrary difference precise integration algorithm has greatly improved the precision by using the time domain local integral semi-analytical method to obtain the recursion operator of wave equations.In this paper,according to the actual examples analysis,the stability of calculation gets better owing to the stable factor constraints;by using an improved adaptive absorbing boundary and the parallelization of serial program,the time consumption of three-dimensional forward modeling is greatly reduced and a large amount of three-dimensional forward modeling is finished.
On the basis of the geological background of Tarim foreland basin,this paper realizes parallel three-dimensional forward modeling in Kuqa area by using three-dimensional arbitrary difference precise integration algorithm.Compared with the conventional differential method,the three-dimensional arbitrary difference precise integration algorithm has greatly improved the precision by using the time domain local integral semi-analytical method to obtain the recursion operator of wave equations.In this paper,according to the actual examples analysis,the stability of calculation gets better owing to the stable factor constraints;by using an improved adaptive absorbing boundary and the parallelization of serial program,the time consumption of three-dimensional forward modeling is greatly reduced and a large amount of three-dimensional forward modeling is finished.
引文
[1] Grary S.Martin,Robert Wiley,Kurt J.Marfurt.Marmousi2:An elastic upgrade for Marmousi.The Leading Edge,2006,25(2) :156～166
[2] Paul Sava.Subsalt Exploration and Development Imaging:Interpretation,and Drilling-What have we learned?The Leading Edge,2006,25(11) :1370～1376
[3] 刘礼农,刘洪,李幼铭.SEG/EAGE盐丘和推覆体模型的波动方程三维叠前深度偏移成像.地球物理学报,2004,47(2) :312～320 Liu L N,Liu H,Li Y M.Wave-equation 3-D prestack depth migration for the SEG/EAGE salt and overthrust model.Chinese J.Geophys.(in Chinese),2004,47(2) :312～320
[4] 孙建国.复杂地表条件下地球物理场数值模拟方法评述.世界地质,2007,26(3) :345～357 Sun J G.Methods for numerical modeling of geophysical fields under complex topographical conditions:a critical review.Global Geology(in Chinese),2007,26(3) :345～357
[5] 钟万勰.子域精细积分及偏微分方程数值解.计算结构力学及其应用,1995,12(3) :253～260 Zhong W X.Subdomain precise integration method and numerical solution of partial differential equations.Computational Structural Mechanics and Applications(in Chinese),1995,12(3) :253～260
[6] 强士中,王孝国,唐茂林等.任意差分精细积分法及数值稳定性分析.应用数学与力学,1999,20(3) :256～262 Qiang S Z,Wang X G,Tang M L.On the arbitrary difference precise integration method and its numerical stability.Applied Mathematics and Mechanics(in Chinese),1999,20(3) :256～262
[7] 贾晓峰,王润秋,胡天跃.求解波动方程的任意差分精细积分法.中国地震,2003,(3) :236～242 Jia X F,Wang R Q,Hu T Y.The Arbitrary Difference Precise Integration Method for Solving the Seismic Wave Equation.Earthquake Research in China(in Chinese),2003,(3) :236～242
[8] 王润秋,张明,张崴.三维精细积分法地震正演及边界条件.石油地球物理勘探,2007,42(2) :202～207 Wang R Q,Zhang M,Zhang W.3 D fine integration seismic forward modeling and boundary condition.Oil Geophysical Prospecting(in Chinese),2007,42(2) :202～207
[9] 牟永光,裴振林.三维复杂介质地震数值模拟.北京.石油工业出版社,2005． 12～14 Mou Y G,Pei Z L.Seismic Numerical Modeling for 3-D Complex Media.Beijing:Petroleum Industry Press(in Chinese),2005． 12～14
[10] Toldi J L,Hale D.Data-dependent absorbing side boundaries.Stanford Exploration Project Report,1982． 111～121
[11] Clayton R W,Engquist B.Absorbing boundary conditions for wave equation migration.Geophysics,1980,45:895 ～904
[2] Paul Sava.Subsalt Exploration and Development Imaging:Interpretation,and Drilling-What have we learned?The Leading Edge,2006,25(11) :1370～1376
[3] 刘礼农,刘洪,李幼铭.SEG/EAGE盐丘和推覆体模型的波动方程三维叠前深度偏移成像.地球物理学报,2004,47(2) :312～320 Liu L N,Liu H,Li Y M.Wave-equation 3-D prestack depth migration for the SEG/EAGE salt and overthrust model.Chinese J.Geophys.(in Chinese),2004,47(2) :312～320
[4] 孙建国.复杂地表条件下地球物理场数值模拟方法评述.世界地质,2007,26(3) :345～357 Sun J G.Methods for numerical modeling of geophysical fields under complex topographical conditions:a critical review.Global Geology(in Chinese),2007,26(3) :345～357
[5] 钟万勰.子域精细积分及偏微分方程数值解.计算结构力学及其应用,1995,12(3) :253～260 Zhong W X.Subdomain precise integration method and numerical solution of partial differential equations.Computational Structural Mechanics and Applications(in Chinese),1995,12(3) :253～260
[6] 强士中,王孝国,唐茂林等.任意差分精细积分法及数值稳定性分析.应用数学与力学,1999,20(3) :256～262 Qiang S Z,Wang X G,Tang M L.On the arbitrary difference precise integration method and its numerical stability.Applied Mathematics and Mechanics(in Chinese),1999,20(3) :256～262
[7] 贾晓峰,王润秋,胡天跃.求解波动方程的任意差分精细积分法.中国地震,2003,(3) :236～242 Jia X F,Wang R Q,Hu T Y.The Arbitrary Difference Precise Integration Method for Solving the Seismic Wave Equation.Earthquake Research in China(in Chinese),2003,(3) :236～242
[8] 王润秋,张明,张崴.三维精细积分法地震正演及边界条件.石油地球物理勘探,2007,42(2) :202～207 Wang R Q,Zhang M,Zhang W.3 D fine integration seismic forward modeling and boundary condition.Oil Geophysical Prospecting(in Chinese),2007,42(2) :202～207
[9] 牟永光,裴振林.三维复杂介质地震数值模拟.北京.石油工业出版社,2005． 12～14 Mou Y G,Pei Z L.Seismic Numerical Modeling for 3-D Complex Media.Beijing:Petroleum Industry Press(in Chinese),2005． 12～14
[10] Toldi J L,Hale D.Data-dependent absorbing side boundaries.Stanford Exploration Project Report,1982． 111～121
[11] Clayton R W,Engquist B.Absorbing boundary conditions for wave equation migration.Geophysics,1980,45:895 ～904
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心