VSP数据逆时偏移方法研究
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
相比于地面地震,VSP数据地震勘探具有高分辨率和高信噪比的特点,目前常规的VSP数据偏移方法难于获得理想的成像效果,逆时偏移方法(RTM)基于双程波动方程,迄今为止被认为是精度最高的地震资料偏移算法.文中依据VSP数据采集特点,从二维声波波动方程入手,研究基于高阶有限差分数值算法的VSP数据逆时偏移,数值计算采用节约部分存储的完全匹配层吸收边界条件(PML),在有效抑制地震波场的边界反射以及保持波场精度的情况下,很大程度上节约了波场存储.文中也采用了互相关成像条件以及拉普拉斯低频噪声压制技术,实现了VSP数据逆时偏移的高精度成像.测试表明,高精度的逆时偏移方法可有效成像VSP地震数据,可实现对小断层构造、大角度构造,以及井周小构造的精细刻画.
Compared with ground seismic exploration,VSP seismic exploration has the characteristics of high resolution and high signal to noise ratio.At present,conventional VSP data migration methods are difficult to obtain the ideal imaging effect.Based on the two-way wave equation,the reverse-time migration(RTM)has become the most accurate migration method so far.In this paper,on the basis of VSP data acquisition characteristics,we deal with VSP data by high order finite difference algorithm during the RTM based on 2D acoustic wave equation.We adopt the perfect matched layer absorbing boundary condition(PML)by the way of setting checkpoints,so that we saved the wave fields storage,in the meantime,effectively suppressed the boundary reflection of seismic wave fields and kept the precision.We also used the cross-correlation imaging condition and laplacian low frequency noise suppression technique to make the high resolution image with reverse-time migration of VSP data come true.The test indicated that high precision reverse-time migration could validly image VSP data and nicely depict small faults,large angle structure and the local area around the well.
引文
Baysal E,Kosloff D D,Sherwood J W C.1983.Reverse timemigration[J].Geophysics,48(11):1514-1524.
    Chang W F,McMechan G A.1986.Reverse time migration of offsetvertical seismic profiling data using the excitation time imagingcondition[J].Geophysics,51(1):67-84.
    Chang W F,McMechan G A.1987.Elastic reverse-time migration[J].Geophysics,52(10):1365-1375.
    Chen K Y.2010.Study on perfectly matched layer absorbingboundary condition[J].Geophysical Prospecting for Petroleum,49(5):472-477.
    Dillon P B.1988.Vertical seismic profile migration using theKirchhoff integral[J].Geophysics,53(6):786-799.
    Dong C H,Wang S X.2010.Finite-difference method with variabletime-steps for reverse-time migration[J].Progress in Geophys,25(4):1298-1305,doi:10.3969/j.issn.1004-2903.2010.04.018.
    Gal’perin E I.1974.Vertical seismic profiling[M].Society ofExploration Geophysicists special Publication.
    Gao G,He Z H,Huang D J,et al.2011.The attenuation factoranalysis perfect matched layer absorbing boundary condition[J].Geophysical Prospecting for Petroleum,50(5):430-433.
    Gao J H,Yang S L.2007.On the method of quality factorsestimation from zero offset VSP data[J].Chinese J.Geophys.(in Chinese),50(4):1198-1209.
    Gao J H,Yang S L,Wang D X.2008.Quality factor extractionusing instantaneous frequency at envelope peak of direct wavesof VSP data[J].Chinese J.Geophys.(in Chinese),51(3):853-861.
    Guo N M,Wu G C.2012.High-order finite difference acousticequation based on PML boundary with variable grids forreverse-time migration[J].Oil geophysics prospecting,47(2):256-265.
    Hu H,Liu Y K,Chang X,et al.2013.Analysis and application ofboundary treatment for the computation of reverse-timemigration.Chinese J,Geophys.(in Chinese),56(6):2033-2042,doi:10.6038/cjg20130624.
    Li B,Liu H W,Liu G F,et al.2010.Computational strategy ofseismic 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 orderfinite different pre-stack reverse time migration and GPUimplementation[J].Chinese J.Geophys.(in Chinese),53(7):1725-1733,doi:10.3969/j.issn.0001-5733.2010.07.024.
    Liu H W,Liu H,Zou Z,et al.2010b.The problems of denoise andstorage 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,Liu H,Li B,et al.2011.Pre-stack reverse timemigration for rugged topography and GPU accelerationtechnology[J].Chinese J.Geophys.(in Chinese),54(7):1883-1892,doi:10.3969/j.issn.0001-5733.2011.07.022.
    Liu S W,Wang H Z,Chen S C,et al.2012.Joint imaging methodof VSP upgoing and downgoing reflection wave[J].Chinese J.Geophys.(in Chinese),3126-3133,doi:10.6038/j.issn.0001-5733.2012.09.030.
    Liu Y,Li C C,Mou Y G.1998.Finite-difference numericalmodeling of any even-order accuracy[J].Oil geophysicalprospecting,33(1):1-10.
    McMechan G A.1983.Migration by extrapolation of timedependent boundary values[J].Geophysical Prospecting,31:413-420.
    Seeman B,Horowicz L.1983.Vertical seismic profiling:separationof upgoing and downgoing acoustic waves in a stratified medium[J].Geophsics,48(5):555-568.
    Wang J,Li Z C,Tao L.2012.The research on imaging condition ofreverse time migration[J].Progress in Geophysics,27(3):1173-1182.
    Wang J,Yang C C,Liu H H,et al.2006.VSP multiwave imagewith Kirchhoff Method[J].Progress in Geophysics,21(3):845-855,doi:10.6038/j.issn.1004-2903.2012.03.042.
    Wang S D.2003.Absorbing boundary condition for acoustic waveequation by perfectly matched layer[J].Oil GeophysicalProspecting,38(1):31-34.
    Whitmore N D.1983.Iterative depth migration by backward timepropagation[A].53rd SEG Meeting,Expanded Abstracts[C],382-385.
    Xu Y.2008.Prestack reverse-time migration by the grid method[J].Progress in Geophysics,23(3):839-845.
    Yu Z,Wang Y C,He J.2010.A stable inverse Q filtering methodfor VSP[J].Progress in Geophysics,25(5):1676-1684,doi:10.3969/j.issn.1004-2903.2010.05.021.
    Zhang S,Zhu H L,Wang J Y,et al.2011.Walk-away VSP and3D-VSP self-adapted static correction method[J].Progress inGeophysics,26(5):1755-1762,doi:10.3969/j.issn.1004-2903.2011.05.030.
    Zhu H L,Zhang S.2012.A high fidelity VSP method of wavefieldseparation[J].Progress in Geophysics,(4):1557-1564,doi:10.6038/j.issn.1004-2903.2012.04.032.
    Zhu G M.1988.Vertical seismic profiling method[M].Beijing:Petroleum industry press.
    Zou Y Y,Xu Y X,Sha C.2009.Review on the forward modelingand inversion of vertical seismic profile[J].Progress inGeophysics,24(1):145-153.
    陈可洋.2010.完全匹配层吸收边界条件研究[J].石油物探,49(5):472-477.
    董春晖,王尚旭.2010.变时间步长有限差分法逆时偏移[J].地球物理学进展,25(4):1298-1305,doi:10.3969/j.issn.1004-2903.2010.04.018.
    高刚,贺振华,黄德济,等.2011.完全匹配层人工边界条件中的衰减因子分析[J].石油物探,50(5):430-433.
    高静怀,杨森林.2007.利用零偏移VSP资料估计介质品质因子方法研究[J].地球物理学报,50(4):1198-1209.
    高静怀,杨森林,王大兴.2008.利用VSP资料直达波的包络峰值处瞬时频率提取介质品质因子[J].地球物理学报,51(3):853-861.
    郭念民,吴国忱.2012.基于PML边界的变网格高阶有限差分声波方程逆时偏移[J].石油地球物理勘探,47(2):256-265.
    胡昊,刘伊克,常旭,等.2013.逆时偏移计算中的边界处理分析及应用[J].地球物理学报,56(6):2033-2042,doi:10.6038/cjg20130624.
    李博,刘红伟,刘国峰,等.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.
    刘红伟,刘洪,李博,等.2011.起伏地表叠前逆时偏移理论及GPU加速技术[J].地球物理学报,54(7):1883-1892,doi:10.3969/j.issn.0001-5733.2011.07.022.
    刘守伟,王华忠,陈生昌,等.2012.VSP上下行反射波联合成像方法研究[J].地球物理学报,3126-3133,doi:10.6038/j.issn.0001-5733.2012.09.030.
    刘洋,李承楚,牟永光.1998.任意偶数阶精度有限差分数值模拟[J].石油地球物理勘探,33(1):1-10.
    王娟,李振春,陶丽.2012.逆时偏移成像条件研究[J].地球物理学进展,27(3):1173-1182,doi:10.6038/j.issn.1004-2903.2012.03.042.
    王珺,杨长春,刘海河,等.2006.克希霍夫法VSP多波联合成像[J].地球物理学进展,21(3):845-855.
    王守东.2003.声波方程完全匹配层吸收边界[J].石油地球物理勘探,38(1):31-34.
    徐义.2008.格子法在起伏地表叠前逆时偏移深度偏移中的应用[J].地球物理学进展,23(3):839-845.
    余振,王彦春,何静.2010.一种稳定的VSP反Q滤波方法[J].地球物理学进展,25(5):1676-1684,doi:10.3969/j.issn.1004-2903.2010.05.021.
    朱海龙,张山.2012.高保真VSP波场分离方法研究[J].地球物理学进展,(4):1557-1564,doi:10.6038/j.issn.1004-2903.2012.04.032.
    张山,朱海龙,王家映,等.2011.Walk-away VSP和3D-VSP自适应静校正方法研究[J].地球物理学进展,26(5):1755-1762,doi:10.3969/j.issn.1004-2903.2011.05.030.
    邹延延,徐义贤,沙椿.2009.VSP正反演综述[J].地球物理学进展,24(1):145-153.
    朱光明.1988.垂直地震剖面方法[M].北京:石油工业出版社.

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