逆时偏移边界条件与存储策略分析
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震波场模拟是叠前逆时偏移方法的核心计算单元,因此,边界条件和存储方案的选择是高精度以及高效率逆时偏移算法不可回避的问题.本文就随机边界条件和吸收边界条件在逆时偏移算法中的应用效果以及相应的四种存储策略展开具体分析,给出了不同存储策略下的计算成本和存储量需求.文中优选随机边界条件和吸收边界条件的有效边界存储方案,引入GPU并行加速技术对SEG起伏地表模型进行测试,详细讨论分析了两种方法的计算效率和计算精度.测试表明,随机边界存储的计算成本要低于有效边界存储的计算成本,而受边界散射的影响,随机边界逆时偏移在浅层伴有噪音干扰.针对具体数据情况,合理地运用GPU加速技术、选择有效的边界条件以及存储方法,是改善逆时偏移成像精度和效率的有效途径.
Seismic wave field modeling is the computing core of pre-stack reverse time migration.Therefore,for improving the accuracy and computational efficiency of reverse time migration,effective boundary condition and storage strategies are the key factors.In this paper,we analyzed the application effect by reverse time migration and four kinds of storage strategies with random and absorbing boundary conditions,summarized the calculation cost and storage requirement for different boundary and storage strategies.We have tested the reverse time migration algorithm on SEG rugged topography model with random boundary condition and absorbing boundary condition with its efficient boundary storage scheme chosen,combined with GPU(Graphic Processing Unit) acceleration technology.The result shows that the computational cost of random boundary storage is lower than which of efficient boundary strategy,but affected by the boundary scattering,the reverse time migration profile with random boundary would get noise interference in shallow layers.The key scheme to improve the accuracy and computational efficiency of reverse time migration is choosing compatible boundary condition,storage strategy,and GPU acceleration technology according to the specific data.
引文
Baysal E,Kosloff D D,Sherwood J W C.1983.Reverse time migration[J].Geophysics,48(11):1514-1524,doi:10.1190/1.1441434.
    Clapp R G.2009.Reverse time migration with random boundaries[C].//79rd Annual International Meeting,SEG Expanded Abstracts,2809-2813.
    Deng F,McMechan G A.2008.Viscoelastic true-amplitude prestack reverse-time depth migration[J].Geophysics,73(4):S143-S155.
    Five ways to avoid storing source wavefield snapshots in 2D elastic prestack reverse time migration.Geophysics,80(1):S1-S18.
    Hu H,Liu Y K,Chang X,et al.2013.Analysis and application of boundary treatment for the computation of reverse-time migration[J].Chinese J.Geophys.,(in Chinese),56(6):2033-2042,doi:10.6038/cjg20130624.
    Li B,Liu G F,Liu H.2009.A method of using GPU to accelerate seismic pre-stack time migration[J].Chinese J.Geophys.,(in Chinese),52(1):245-252.
    Li B,Liu H W,Liu G F,et al.2010.Computational strategy of seismic pre-stack reverse time migration on CPU/GPU[J].Chinese J.Geophys.,(in Chinese),53(12):2938-2943,doi:10.3969/j.issn.0001-5733.2010.12.017.
    Liu H W,Li B,Liu H,et al.2010a.The algorithm of high order finite difference pre-stack reverse time migration and GPU implementation[J].Chinese J.Geophys.,(in Chinese),53(7):1725-1733,doi:10.3969/j.issa 0001-5733.2010.07.024.
    Liu H W,Liu H,Zou Z.et al.2010b.The problems of denoise and storage in seismic reverse time migration[J].Chinese J.Geophys.,(in Chinese),53(9):2171-2180,doi:10.3969/j.issn.0001-5733.2010.09.017.
    Liu H W,Ding R W,Liu L,et al.2013.Wavefield reconstruction methods for reverse time migration[J].Journal of Geophysics and Engineering,10:1-6.
    Liu S Z,Shi Y,Guo X B.et al.2014.Investigation on reverse-time migration of VSP data[J].Progress in Geophysics.(in Chinese),29(5):2211-2218,doi:10.6038/pg20140533.
    Liu S W,Wang H Z,Chen S C,et al.2013.Implementation strategy of 3D reverse time migration on GPU/CPU clusters.Chinese].Geophys.,(in Chinese),56(10):3487-3496,doi:10.6038/cjg20131023.
    McMechan G A.1983.Migration by extrapolation of timedependent boundary values[J].Geophysical Prospecting,31:413-420.
    Nguyen B D,McMechan G A.2015.Symes W W.2007.Reverse time migration with optimal checkpointing[J].Geophysics,72(5):SM213-SM221,doi:10.1190/1.2742686.
    Wang B L,Gao J H,Chen W C,et al.2012.Efficient boundary storage strategies for seismic reverse time migration[J].Chinese J.Geophys.,(in Chinese),55(7):2412-2421,doi:10.6038/j.issn.0001-5733.2012.07.025.
    Wang W H,Ke X,Pei J Y.2013.Application investigation of perfectly matched layer absorbing boundary condition[J].Progress in Geophysics,(in Chinese),28(5):2508-2514,doi:10.6038/pg20130529.
    Whitmore N D.1983.Iterative depth migration by backward time propagation[C].53rd Annual International Meeting,SEG Expanded Abstracts,382-385.
    Yan R,Xie X B.2012.An angle-domain imaging condition for elastic reverse time migration and its application to angle gather extraction[J].Geophysics,75(5):S105-S115.
    胡昊,刘伊克,常旭,等.2013.逆时偏移计算中的边界处理分析及应用[J].地球物理学报,56(6):2033-2042,doi:10.6038/cjg20130624.
    李博,刘国峰,刘洪.2009.地震叠前时间偏移的一种图形处理器提速实现方法[J].地球物理学报,52(1):245-252.
    李博,刘红伟,刘国峰,等.2010.地震叠前逆时偏移算法的CPU/GPU实施对策[J].地球物理学报,53(12):2938-2943,doi:10.3969/j.issn.0001-5733.2010.12.017.
    刘红伟,李博,刘洪,等.2010a.地震叠前逆时偏移高阶有限差分算法及GPU实现[J].地球物理学报.53(7):1725-1733,doi:10.3969/j.issn.0001-5733.2010.07.024.
    刘红伟,刘洪,邹振,等.2010b.地震叠前逆时偏移中的去噪与存储[J].地球物理学报,53(9):2171-2180,doi:10.3969/j.issn.0001-5733.2010.09.017.
    刘诗竹,石颖,郭雪豹,等.2014.VSP数据逆时偏移方法研究[J].地球物理学进展,29(5):2211-2218,doi:10.6038/pg20140533.
    刘守伟,王华忠,陈生昌,等.2013.三维逆时偏移GPU/CPU机群实现方案研究[J].地球物理学报,56(10):3487-3496.doi:10.6038/cjg20131023.
    王保利,高静怀,陈文超,等.2012.地震叠前逆时偏移的有效边界存储策略[J].地球物理学报,55(7):2412-2421,doi:10.6038/j.issn.0001-5733.2012.07.025.
    王维红,柯璇,裴江云.2013.完全匹配层吸收边界条件应用研究[J].地球物理学进展,28(5):2508-2514.doi:10.6038/pg20130529.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心