时间域偏移成像方法和装置.pdf

(19)中华人民共和国国家知识产权局(12)发明专利(10)授权公告号CN105607117B(45)授权公告日2018.05.15(21)申请号201510577384.8(56)对比文件(22)申请日2015.09.11CN102841379A,2012.12.26,全文.CN102914791A,2013.02.06,全文.(65)同一申请的已公布的文献号GB2379271B,2004.10.27,全文.申请公布号CN105607117AUS2005/0105393A1,2005.05.19,全文.(43)申请公布日2016.05.25马召贵等.基于CPU+GPU联合计算真地表叠(73)专利权人中国石油集团东方地球物理勘探前时间偏移实用化研究.《石油地球物理勘探》有限责任公司.2014,第49卷(第3期),第1栏第3-11行地址072751河北省涿州市范阳西路189号及第2栏第34行-第1栏第10行.专利权人中油油气勘探软件国家工程研究韩复兴等.波前构建法研究现状.《地球物理中心有限公司学进展》.2011,第26卷(第3期),摘要及第1栏第2行-第1栏第3行.(续)(72)发明人岳玉波钱忠平史云燕孙鹏远李建峰徐少波审查员张秀丽(74)专利代理机构北京三友知识产权代理有限公司11127代理人党晓林(51)Int.Cl.G01V1/28(2006.01)G01V1/32(2006.01)权利要求书2页说明书7页附图6页(54)发明名称时间域偏移成像方法和装置(57)摘要本发明提供了一种时间域偏移成像方法和装置，其中，该方法包括：读取一道地震数据；通过双平方根公式，求取该道的震源到成像道粗网格成像点的第一走时，和粗网格成像点到接收点的第二走时，并将所述第一走时和所述第二走时之和作为粗网格成像点的双程走时；通过三次卷积插值公式，根据所述粗网格成像点的双程走时，插值计算出细网格成像点的双程走时；根据所述细网格成像点的双程走时，在输入地震道上拾取能量，作为该输入道的成像结果，最终的偏移结果是所有输入地震道成像结果的叠加。本发明解决了现有技术中偏移叠加剖面中的高频信息被严重破坏的技术问题，达到了有效提高成像结果精度和准确度的技术效果。CN105607117B[转续页]CN105607117B2/2页[接上页](56)对比文件.1981,第ASSP-29卷(第6期),第2栏第ROBERTG.KEYS.CubicConvolution38行-第2栏第35行.InterpolationforDigitalImage许建国.克希霍夫叠前时间偏移走时算法应Processing.《IEEETRANSACTIONSON用效果评价.《中外能源》.2015,第20卷(第5期),ACOUSTICS,SPEECH,ANDSIGNALPROCESSING》第36-40页.2CN105607117B权利要求书1/2页1.一种时间域偏移成像方法，其特征在于，包括：读取一个地震道的地震数据；通过双平方根公式，求取该地震道的震源到成像道粗网格成像点的第一走时，和粗网格成像点到接收点的第二走时，并将所述第一走时和所述第二走时之和作为粗网格成像点的双程走时，所述双平方根公式为：其中，t表示粗网格成像点的双程走时，ts表示震源到成像道粗网格成像点的第一走时，tg表示粗网格成像点到接收点的第二走时，t0表示所述粗网格成像点对应的垂向走时深度，vrms表示所述粗网格成像点处的均方根速度，x表示所述粗网格成像点到共中心点的水平距离，h表示半炮检距；通过三次卷积插值公式，根据所述粗网格成像点的双程走时，插值计算出细网格成像点的双程走时；根据所述细网格成像点的双程走时，在该地震道上拾取能量，并将拾取的能量进行叠加，将叠加结果作为该地震道的成像结果；对所有地震道重复执行上述步骤，得到一个成像道的偏移成像结果。2.如权利要求1所述的方法，其特征在于，读取一个地震道的地震数据，包括：获取待分析区域的CMP道集地震数据；从所述CMP道集地震数据中读取一道地震数据。3.如权利要求1所述的方法，其特征在于，所述三次卷积插值的公式为：3232tn＝Tk-1(-s+2s-s)/2+Tk(3s-5s+2)/23232+Tk+1(-3s+4s+s)/2+Tk+2(s-s)/2其中，tn表示插值计算后得到的细网格成像点的双程走时，Tk-1、Tk、Tk+1、Tk+2表示细网格成像点附近的四个粗网格成像点的走时，s表示细网格步长。4.如权利要求1至3中任一项所述的方法，其特征在于，所述粗网格成像点是按照预定的网格采样间隔，对成像道进行稀疏采样得到的。5.一种时间域偏移成像装置，其特征在于，包括：读取模块，用于读取一个地震道的地震数据；第一双程走时计算模块，用于通过双平方根公式，求取该地震道的震源到成像道粗网