Parameterized random boundary condition and its application on seismic RTM
ZHANG Li-Mei1, CHENG Gu1, 2
1.School of Earth Sciences and Engineering,Sun Yat-Sen University,Guangzhou 510275,China; 2.Guangdong Key Laboratory of Geological Process and Mineral Resources Exploration,Guangzhou 510275,China
Abstract:Pre-stack seismic reverse time migration has important application in high accuracy subsurface imaging.However enormous memory is needed in RTM.Using random boundary in RTM can reduce memory need of RTM effectively.In this paper,random media expressed by correlation function,correlation length,root mean square velocity perturbation is used,the meanings and values of parameters are analyzed and the optimal parameter group(exponential correlation function,the correlation length being equal to wave length and root mean square velocity perturbation being 30%) are obtained.Numerical simulation shows that using random boundary with optimal parameter group can destroy the coherence of artificial boundary reflection wave fields effectively and obtain high accuracy imaging result with less memory at the same time.
[1] Whitmore N D.Iterative depth migration by backward time propagation[C]//Expanded Abstracts of the 53 rd Annual SEG Meeting.Las Vegas,Nevada:Society of Exploration Geophysicists,1983:382-385. [2] McMechan G A.Migration by extrapolation of time-dependent boundary values[J].Geophysical Prospecting,1983,31(3):413-420. [3] Baysal E,Kosloff D D,Sherwood J W C.Reverse time migration[J].Geophysics,1983,48(11):1514-1524. [4] Griewank A.Achieving logarithmic growth of temporal and spatial complexity in reverse automatic differentiation[J].Optimization Methods and software,1992,1(1):35-54. [5] Symes W W.Reverse time migration with optimal checkpointing[J].Geophysics,2007,72(5):SM213-SM221. [6] Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms: Numerical results[J].Geophysics,1986,51(7):1387-1403. [7] Tan S R,Huang L J.Reducing the computer memory requirement for 3D reverse-time migration with a boundary-wavefield extrapolation method[J].Geophysics,2014,79(5):S185-S194. [8] Clapp R G.Reverse time migration with random boundaries[C]//Expanded Abstracts of the 79 th Annual SEG Meeting.Houston,Texas:Society of Exploration Geophysicists,2009. [9] Shen X,Clapp R G.Random boundary condition for low-frequency wave propagation[C]//Expanded Abstracts of the 81 st Annual SEG Meeting.San Antonio,Texas:Society of Exploration Geophysicists,2011,30:2962-2965. [10] Shen X K,Clapp R G.Random boundary condition for memory-efficient waveform inversion gradient computation[J].Geophysics,2015,80(6):R351-R359. [11] Clayton R,Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations[J].Bulletin of the seismological society of America,1977,67(6):1529-1540. [12] Cerjan C,Kosloff D,Kosloff R,et al.A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50(4):705-708. [13] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of computational physics,1994,114(2):185-200. [14] 李博,刘红伟,刘国峰,等.地震叠前逆时偏移算法的CPU/GPU实施对策[J].地球物理学报,2010,12:2938-2943. [15] 胡昊,刘伊克,常旭,等.逆时偏移计算中的边界处理分析及应用[J].地球物理学报,2013,56(6):2033-2042. [16] 石颖,柯璇,张莹莹.逆时偏移边界条件与存储策略分析[J].地球物理学进展,2015,30(2):581-585. [17] Shi Y,Wang Y H.Reverse time migration of 3D vertical seismic profile data[J].Geophysics,2016,81(1):S31-S38. [18] 张慧宇,王立歆,孔祥宁,等.适于大规模数据的逆时偏移实现方案[J].物探与化探,2015,39(5):978-984. [19] Ikelle L T,Yung S K,Daube F.2-D random media with ellipsoidal autocorrelation functions[J].Geophysics,1993,58(9):1359-1372. [20] 奚先,姚姚.二维随机介质及波动方程正演模拟[J].石油地球物理勘探,2001,36(5):546-552. [21] Kerner C.Anisotropy in sedimentary rocks modeled as random media[J].Geophysics,1992,57(4):564-576. [22] Korn M.Seismic waves in random media[J].Journal of applied geophysics,1993,29(3-4):247-269. [23] Ergintav S,Canitez N.Modeling of multi-scale media in discrete form[J].Journal of seismic exploration,1997,6(1):77-96. [24] Flatté S M,Wu R S.Small-scale structure in the lithosphere and asthenosphere deduced from arrival time and amplitude fluctuations at NORSAR[J].Journal of Geophysical Research: Solid Earth,1988,93(B6):6601-6614. [25] 姚姚,奚先.区域多尺度随机介质模型及其波场分析[J].石油物探,2004,43(1):1-7. [26] 奚先,姚姚.二维横各向同性弹性随机介质中的波场特征[J].地球物理学进展,2004,19(4):924-932. [27] 奚先,姚姚.非平稳随机介质模型[J].石油地球物理勘探,2005,40(1):71-75. [28] 陈可洋,刘洪林,杨微,等.随机介质模型的改进方法及应用[J].大庆石油地质与开发,2008,27(5):124-126,131. [29] 王金山,陈可洋,吴清岭,等.随机介质模型的一种构造方法[J].物探与化探,2010,34(2):191-194. [30] 李灿苹,王南萍,李志宏.Von Karman 型自相关函数模拟随机介质[J].物探与化探,2010,34(1):98-102. [31] 殷学鑫,刘洋.二维随机介质模型正演模拟及其波场分析[J].石油地球物理勘探,2011,46(6):862-872. [32] Xie X B.Seismic wave scattering in 3D random media:A finite-difference simulation and slowness domain analysis[C]//Expanded Abstracts of the 83 rd Annual SEG Meeting.Houston,Texas:Society of Exploration Geophysicists,2013:3428-3432. [33] 顾元,朱培民,李辉,等.二维叠后地震数据的平稳随机介质参数估计[J].地球物理学报,2014,57(7):2291-2301. [34] Frankel A,Clayton R W.Finite difference simulations of seismic scattering:Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity[J].Journal of Geophysical Research:Solid Earth,1986,91(B6):6465-6489. [35] 吴如山.地震波散射:理论与应用[J].地球物理学进展,1989,4(4):1-23. [36] 王保利,高静怀,陈文超,等.地震叠前逆时偏移的有效边界存储策略[J].地球物理学报,2012,55(7):2412-2421.