CN109917454A - 基于双基准面的真地表叠前深度偏移成像方法及装置 - Google Patents
基于双基准面的真地表叠前深度偏移成像方法及装置 Download PDFInfo
- Publication number
- CN109917454A CN109917454A CN201910121989.4A CN201910121989A CN109917454A CN 109917454 A CN109917454 A CN 109917454A CN 201910121989 A CN201910121989 A CN 201910121989A CN 109917454 A CN109917454 A CN 109917454A
- Authority
- CN
- China
- Prior art keywords
- double
- true earth
- datum levels
- depth migration
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于双基准面的真地表叠前深度偏移成像方法及装置,包括:基于微测井先验信息对拾取的初至波进行网格层析反演得到近地表速度模型;在近地表速度模型中融入中深层速度场获得初始速度场;基于初始速度场对井炮地震数据进行克希霍夫双基准面真地表偏移得到共成像点道集;根据共成像点道集利用双基准面网格层析法对初始速度场进行优化迭代获得最终迭代速度场;基于最终迭代速度场对井炮地震数据进行克希霍夫双基准面真地表偏移得到最终的叠前深度偏移成像结果;双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。由于该方案设置了双基准面,可以提高地震资料的成像精度。
Description
技术领域
本发明涉及地震资料处理技术领域,特别涉及一种基于双基准面的真地表叠前深度偏移成像方法及装置。
背景技术
中国西部探区速度结构复杂,叠前时间偏移技术的前提条件已经不成立。随着勘探程度的深入,叠前深度偏移已经开始广泛应用。西部探区地表条件一般比较恶劣,多为大沙漠、山地地形。为了获得更多更强的地震反射信息,很多地震采集工程采用炸药震源,在潜水面以下激发,在地表埋置检波器接收,即激发点和接收点不在一个面上。以往的叠前深度偏移技术,由于无法适应浅层速度和地表地形的剧烈变化,是通过静校正技术对地震数据做时移,同时将真实的深度域速度模型替换为平滑地表、浅层填充替换速度的速度模型,这样做破坏了实际的波场传播路径,在旅行时计算时会引起误差。为了克服这种缺陷,已经发展了一些真地表偏移技术,但是目前的真地表偏移技术都是基于单个偏移基准面的,即需要首先通过静校正将炮点从激发点校正到地表,再从地表开始偏移。当激发点深度较深时,这种校正方法仍然会破坏实际的波场传播路径,偏离了真地表偏移成像技术的初衷,如图1所示,从激发点S到检波点G的真实的传播路径为实线型箭头线,而目前的偏移方法将静校正处理放在偏移前引起了波场畸变,导致从激发点S到检波点G的路径为虚线型箭头线,这样会造成旅行时的计算不准确,进而造成地震资料的成像精度不高。
发明内容
本发明实施例提供了一种基于双基准面的真地表叠前深度偏移成像方法及装置,提高了地震资料的成像精度。
该基于双基准面的真地表叠前深度偏移成像方法包括:
基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集;
根据共成像点道集,利用双基准面网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
该基于双基准面的真地表叠前深度偏移成像装置包括:
近地表速度模型确定模块,用于基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
初始速度场确定模块,用于在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
真地表偏移模块,用于基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集;
速度优化迭代模块,用于根据共成像点道集,利用双基准面网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
所述真地表偏移模块还用于:基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述基于双基准面的真地表叠前深度偏移成像方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述基于双基准面的真地表叠前深度偏移成像方法的计算机程序。
在本发明实施例中,基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型,将微测井先验信息作为约束条件可以提高模型反演的精度;通过将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上,设置了双基准面,然后进行克希霍夫双基准面真地表偏移、初始速度场进行优化迭代、再次克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果,这样可以提高地震资料的成像精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种真地表偏移射线路径图;
图2是本发明实施例提供的一种基于双基准面的真地表叠前深度偏移成像方法流程图;
图3是本发明实施例提供的一种微测井约束前后的层析反演近地表模型(下)、与VSP速度对比(上)示意图;
图4是本发明实施例提供的一种中国西部大沙漠地区井炮激发资料的激发接收位置示意图,其中,横坐标表示道序号,纵坐标表示深度;
图5是本发明实施例提供的一种双基准面偏移中的炮点偏移面(下)和检波点偏移面(上)示意图,其中,纵坐标表示深度;
图6是本发明实施例提供的一种常规的CMP面偏移输入道集(下)和双基准面偏移输入道集(上)示意图,其中,横坐标表示道序号,纵坐标表示深度;
图7是本发明实施例提供的一种常规的CMP面偏移速度(下)和双基准面偏移速度(上)示意图,其中,横坐标表示速度,纵坐标表示深度;
图8是本发明实施例提供的一种常规的CMP面偏移CRP道集(下)和双基准面偏移CRP道集(上)示意图,其中,横坐标表示道序号,纵坐标表示深度;
图9是本发明实施例提供的一种常规的CMP面偏移叠加剖面(下)和双基准面偏移叠加剖面(上)示意图,其中,横坐标表示道序号,纵坐标表示深度;
图10是本发明实施例提供的一种基于双基准面克希霍夫偏移CIP道集的双基准面网格层析剩余深度拾取(下,其中,横坐标表示道序号,纵坐标表示深度)和速度更新量(上)示意图;
图11是本发明实施例提供的一种基于双基准面的真地表叠前深度偏移成像装置结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明实施例中,提供了一种基于双基准面的真地表叠前深度偏移成像方法,如图2所示,该方法包括:
步骤201:基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
步骤202:在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
步骤203:基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集(一般是指CMP偏移后得到的CRP道集);
步骤204:根据共成像点道集,利用网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
步骤205:基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
在本发明实施例中,步骤201具体实现如下:
真地表叠前深度偏移成像技术的基础在于近地表速度建模,由于层析反演方法的目标函数是正演走时与实际走时之差,就是说,当根据现有模型计算的正演走时与实际走时之差足够小时,就认为该模型是收敛的。走时是慢度(和速度相对,是速度的倒数)和深度的乘积,就是说,同样的走时,可能对应多种速度和深度的组合,即反演的结果是多解的。但是,传统的层析反演方法由于缺少深度的约束,得到的速度和深度模型往往都偏大,仅可以用于计算静校正量,并不能用于叠前深度偏移速度建模。因此,为了从多解中得到更接近真实情况的解,需要对速度和深度做限制。所谓“先验信息”,就是“预先知道的信息”,在本发明中,微测井就是预先知道的信息,通过微测井调查,可以知道工区里的某些点的速度和深度值,将这些点的信息作为约束条件,在反演过程中对结果加以限制,可以大大提高模型反演的精度,能用于深度偏移融合建模。
在本发明实施例中,步骤202里的“融入”可以理解为“拼接”,对近地表反演得到的模型,根据反演使用的初至波包含的偏移距范围,有一个认为相对可靠的反演深度,将地表到这个深度的速度场,拼接一个中深层的速度场,形成一个真地表偏移的初始速度场。这个拼接可以有过渡的斜坡。
在本发明实施例中,虽然微测井约束的层析反演方法可以得到在趋势上与实际地质情况吻合较好的近地表速度,但是由于近地表速度很低,所以相对较小的误差就能导致较大的时差问题。由于层析本质上还是基于射线的方法,是对波动方程的近似,另外由于反演网格大小的限制等原因,最终反演的结果难免存在一定的误差,这种误差导致的时差就是一种高频静校正问题。基于此,本发明提出的基于双基准面的真地表叠前深度偏移成像方法还可以包括:
步骤203-1:在对井炮地震数据进行克希霍夫双基准面真地表偏移之前,对所述井炮地震数据进行高频剩余静校正,获得静校正数据。
也就是说,层析反演得到的速度模型不能完成解决静校正中的高频问题,需要通过初至波剩余静校正计算这部分高频校正量,应用到井炮地震数据上,作为反演速度模型误差的一个补偿。
下面从具体的实现公式上来说明本发明方法。
首先,基于起伏地表的波前重建射线追踪和旅行时计算。层析反演技术和偏移成像技术,最重要的组成部分都是射线追踪和旅行时计算。目前主流的旅行时计算方法包括基于费马原理、基于程函方程有限差分、基于波前重建的三类方法,前面两种方法速度快,但是不适应速度和界面的剧烈变化,对于近地表层析和真地表偏移来说,必须采用基于波前重建的射线追踪和旅行时计算方法。波前构建法以运动学射线追踪和动力学射线追踪组成的方程组为基本出发点,实现波前射线路径、走时以及振幅的计算,其方程组为:
其中,xi是第i个网格位置的坐标分量,v是波的传播速度,pi是第i个网格表示的慢度分量,τ是时间;Q、P、V均是2×2阶的矩阵,Q是从射线参数坐标到射线中心坐标的转换矩阵,P是从射线参数坐标到射线中心坐标的慢度矢量分量的转换矩阵,V是速度的二阶导数;表示射线路径矩阵。
这里指的是射线路径求取和旅行时的计算方法原理。在本发明中包含近地表速度反演和偏移成像的内容,这两个部分都需要做射线路径求取和旅行时计算。以近地表反演为例,对给定的炮点(震源)和检波点(R)位置,在初至速度模型的条件下,需要在某种物理原则条件下(最小走时、最大能量等),通过某种数学方法(波前重建等)找出一条路径,根据该路径和对应网格的速度可以求出旅行时,如果该旅行时和实际走时有差异,则修改初始速度模型。偏移成像则是上述的前半部分,即根据速度和炮检点位置计算传播路径和走时,再根据走时把地震数据求和的过程。求路径和走时的方法有很多,但是精度各有不同。本发明要强调的是,由于风化侵蚀等作用,近地表的速度往往变化剧烈,必须用精度较高的波场构建的方法来求取路径和走时。
然后,基于微测井等先验信息约束的初至波层析反演。初至波反演的基本公式为:
其中,表示路径矩阵,与深度相关;表示慢度矩阵,是速度的倒数;是走时矩阵,即初至时间;j代表炮检对的序号,一共M个炮检对,即M个走时,j=1,2…,M;N表示网格的个数,i=1,2…,N。对一个炮检对j来说,矩阵A中的第j行是它的射线路径,每个网格的路径乘以对应网格的慢度,得到第j个炮检对的走时。
最后,对井炮地震数据进行克希霍夫双基准面真地表偏移,所使用的是三维克希霍夫积分公式,可以写成如下形式:
其中,U是地震数据的振幅,它可以看做是炮检点坐标和时间的函数;W是振幅加权因子;是炮点坐标;是检波点坐标;ts是炮点到成像点的旅行时;tr是检波点到成像点的旅行时;t是炮点到检波点的旅行时;Ω是偏移孔径。
对给定的成像点坐标x,对偏移孔径Ω内的某一地震道,可以根据该道的炮检点坐标,计算炮点到成像点再到检波点的走时,将该走时处的振幅值按时间的偏导数加权取出,将孔径Ω内所有到按此处理再求和,作为该成像点x处的成像结果,既是克希霍夫偏移的过程。
上述公式的逻辑关系是:在公式(1)的基础上,可以求出公式(2)中的矩阵A,再根据公式(2)做速度反演得到近地表速度模型或者做中深层的速度更新;同样在公式(2)求得的模型下,根据公式(1)可以求出公式(3)中的走时T,在根据公式(3)对地震数据求和即可完成偏移成像。
下面通过具体的效果图来说明本发明提出的基于双基准面的真地表叠前深度偏移成像方法的优点。其中,图3至图10是效果图,横纵坐标的清楚与否不影响附图效果的展示。
图3的下图为微测井约束前后的层析反演近地表模型,下左图为约束前模型,下右图为约束后模型,可以看到约束后的模型,其沙丘的速度特征更符合地质认识:速度更低,成层性更好。图3的上图为约束后模型速度与VSP速度对比,可以看到约束反演得到的近地表速度与VSP速度吻合度非常高。
图4为中国西部大沙漠地区井炮激发资料的激发接收位置示意图,炮点位于潜水面以下,检波点位于地表,两者高差最大可能达到上百米。
图5为双基准面偏移中的两个偏移基准面,下为炮点偏移面,由炮点激发高程得到,上为检波点偏移面,由地表高程得到。
图6左为常规的CMP面偏移输入道集,上为双基准面偏移输入道集。可以看到由于双基准面偏移为真正意义的真地表偏移,所以该方法的输入道集只需应用近地表层析反演无法解决的少量高频问题,基本无需应用静校正,因此初至上能看到明显的静校正问题,这些问题将在深度偏移中以非地表一致性的方式得到解决。
图7下为常规的CMP面偏移速度,上为双基准面偏移速度。可以看到双基准面偏移是基于真地表的,对应的偏移速度包含了真实的近地表速度,速度值最低可低至500m/s。
图8下为常规的CMP面偏移CRP道集,上为双基准面偏移CRP道集。图9下为常规的CMP面偏移叠加剖面,上为双基准面偏移叠加剖面。可以看到虽然双基准面偏移的输入道集没有做任何静校正,但是静校正问题都在深度偏移中得到了解决,而且从道集上看,浅层的信噪比得到了明显的提高,深层成像也得到了一定的改善。
图10是基于双基准面克希霍夫偏移CIP道集的双基准面网格层析剩余深度拾取和速度更新量,可以看到速度更新量是收敛的,可以用于偏移速度优化迭代。
基于同一发明构思,本发明实施例中还提供了一种基于双基准面的真地表叠前深度偏移成像装置,如下面的实施例所述。由于基于双基准面的真地表叠前深度偏移成像装置解决问题的原理与基于双基准面的真地表叠前深度偏移成像方法相似,因此基于双基准面的真地表叠前深度偏移成像装置的实施可以参见基于双基准面的真地表叠前深度偏移成像方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图11是本发明实施例的基于双基准面的真地表叠前深度偏移成像装置的一种结构框图,如图11所示,包括:
近地表速度模型确定模块1101,用于基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
初始速度场确定模块1102,用于在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
真地表偏移模块1103,用于基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集;
速度优化迭代模块1104,用于根据共成像点道集,利用双基准面网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
所述真地表偏移模块1103还用于:基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
在本发明实施例中,所述微测井先验信息包括速度和深度。
在本发明实施例中,还包括:
静校正模块,用于在对井炮地震数据进行克希霍夫双基准面真地表偏移之前,对所述井炮地震数据进行高频剩余静校正,获得静校正数据。
在本发明实施例中,所述近地表速度模型确定模块1101具体用于:
按照如下公式基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型:
其中,xi是第i个网格位置的坐标分量,v是波的传播速度,pi是第i个网格表示的慢度分量,τ是时间;Q、P、V均是2×2阶的矩阵,Q是从射线参数坐标到射线中心坐标的转换矩阵,P是从射线参数坐标到射线中心坐标的慢度矢量分量的转换矩阵,V是速度的二阶导数;表示路径矩阵,表示慢度矩阵;是走时矩阵;M表示走时的个数,j=1,2…,M;N表示网格的个数,i=1,2…,N。
在本发明实施例中,所述真地表偏移模块1103具体用于:
按照如下公式对井炮地震数据进行克希霍夫双基准面真地表偏移:
其中,U是地震数据的振幅;W是振幅加权因子;是炮点坐标;是检波点坐标;ts是炮点到成像点的旅行时;tr是检波点到成像点的旅行时;t是炮点到检波点的旅行时;Ω是偏移孔径。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述基于双基准面的真地表叠前深度偏移成像方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述基于双基准面的真地表叠前深度偏移成像方法的计算机程序。
综上所述,为了在偏移时,得到更准确的旅行时,本发明提出了一种基于激发、接收双基准面的真地表叠前深度偏移成像方法及装置,即在偏移计算旅行时时,将炮点和检波点分别置于其真实激发和接收的高程位置。基于该技术,形成了以微测井等先验信息约束的层析反演近地表速度建模、双基准面克希霍夫叠前深度偏移成像、双基准面网格层析速度优化的真地表叠前深度偏移配套处理技术。该配套处理技术可以完全将静校正问题融入到叠前深度偏移解决,解决了以往将静校正处理放在偏移前进行引起的波场畸变问题,可以使得旅行时的计算准确,从而提高地震资料的成像精度。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种基于双基准面的真地表叠前深度偏移成像方法,其特征在于,包括:
基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集;
根据共成像点道集,利用双基准面网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
2.如权利要求1所述的基于双基准面的真地表叠前深度偏移成像方法,其特征在于,所述微测井先验信息包括速度和深度。
3.如权利要求1所述的基于双基准面的真地表叠前深度偏移成像方法,其特征在于,在对井炮地震数据进行克希霍夫双基准面真地表偏移之前,还包括:
对所述井炮地震数据进行高频剩余静校正,获得静校正数据。
4.如权利要求1所述的基于双基准面的真地表叠前深度偏移成像方法,其特征在于,按照如下公式基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型:
其中,xi是第i个网格位置的坐标分量,v是波的传播速度,pi是第i个网格表示的慢度分量,τ是时间;Q、P、V均是2×2阶的矩阵,Q是从射线参数坐标到射线中心坐标的转换矩阵,P是从射线参数坐标到射线中心坐标的慢度矢量分量的转换矩阵,V是速度的二阶导数;表示路径矩阵,表示慢度矩阵;是走时矩阵;M表示走时的个数,j=1,2…,M;N表示网格的个数,i=1,2…,N。
5.如权利要求1所述的基于双基准面的真地表叠前深度偏移成像方法,其特征在于,按照如下公式对井炮地震数据进行克希霍夫双基准面真地表偏移:
其中,U是地震数据的振幅;W是振幅加权因子;是炮点坐标;是检波点坐标;ts是炮点到成像点的旅行时;tr是检波点到成像点的旅行时;t是炮点到检波点的旅行时;Ω是偏移孔径。
6.一种基于双基准面的真地表叠前深度偏移成像装置,其特征在于,包括:
近地表速度模型确定模块,用于基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型;
初始速度场确定模块,用于在所述近地表速度模型中融入中深层速度场,形成具有深度偏移的初始速度场;
真地表偏移模块,用于基于初始速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到共成像点道集;
速度优化迭代模块,用于根据共成像点道集,利用双基准面网格层析方法对所述初始速度场进行优化迭代,获得最终迭代速度场;
所述真地表偏移模块还用于:基于所述最终迭代速度场,对井炮地震数据进行克希霍夫双基准面真地表偏移,得到最终的叠前深度偏移成像结果;
其中,所述双基准面为将射线追踪的起点放置于激发点高程面上,射线追踪的终点放置于接收点高程上。
7.如权利要求6所述的基于双基准面的真地表叠前深度偏移成像装置,其特征在于,所述微测井先验信息包括速度和深度。
8.如权利要求6所述的基于双基准面的真地表叠前深度偏移成像装置,其特征在于,还包括:
静校正模块,用于在对井炮地震数据进行克希霍夫双基准面真地表偏移之前,对所述井炮地震数据进行高频剩余静校正,获得静校正数据。
9.如权利要求6所述的基于双基准面的真地表叠前深度偏移成像装置,其特征在于,所述近地表速度模型确定模块具体用于:
按照如下公式基于微测井先验信息,对拾取的初至波进行网格层析反演,得到近地表速度模型:
其中,xi是第i个网格位置的坐标分量,v是波的传播速度,pi是第i个网格表示的慢度分量,τ是时间;Q、P、V均是2×2阶的矩阵,Q是从射线参数坐标到射线中心坐标的转换矩阵,P是从射线参数坐标到射线中心坐标的慢度矢量分量的转换矩阵,V是速度的二阶导数;表示路径矩阵,表示慢度矩阵;是走时矩阵;M表示走时的个数,j=1,2…,M;N表示网格的个数,i=1,2…,N。
10.如权利要求6所述的基于双基准面的真地表叠前深度偏移成像装置,其特征在于,所述真地表偏移模块具体用于:
按照如下公式对井炮地震数据进行克希霍夫双基准面真地表偏移:
其中,U是地震数据的振幅;W是振幅加权因子;是炮点坐标;是检波点坐标;ts是炮点到成像点的旅行时;tr是检波点到成像点的旅行时;t是炮点到检波点的旅行时;Ω是偏移孔径。
11.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5任一所述基于双基准面的真地表叠前深度偏移成像方法。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至5任一所述基于双基准面的真地表叠前深度偏移成像方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910121989.4A CN109917454B (zh) | 2019-02-19 | 2019-02-19 | 基于双基准面的真地表叠前深度偏移成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910121989.4A CN109917454B (zh) | 2019-02-19 | 2019-02-19 | 基于双基准面的真地表叠前深度偏移成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109917454A true CN109917454A (zh) | 2019-06-21 |
CN109917454B CN109917454B (zh) | 2020-10-09 |
Family
ID=66961773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910121989.4A Active CN109917454B (zh) | 2019-02-19 | 2019-02-19 | 基于双基准面的真地表叠前深度偏移成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109917454B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624647A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111751875A (zh) * | 2020-07-07 | 2020-10-09 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp带限角度积分叠前时间偏移方法和装置 |
CN112270031A (zh) * | 2020-11-10 | 2021-01-26 | 黄河勘测规划设计研究院有限公司 | 一种多沙河流水库的双泥沙侵蚀基准面排沙系统 |
CN113031069A (zh) * | 2021-03-02 | 2021-06-25 | 吉林大学 | 一种岩溶地区的多元信息约束智能层析静校正方法 |
CN113568040A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5570321A (en) * | 1994-03-03 | 1996-10-29 | Atlantic Richfield Company | Seismic velocity model optimization method using simulated annearling to determine prestack travel-times |
US20050038604A1 (en) * | 2003-08-14 | 2005-02-17 | Marin Irina A. | Method and apparatus for kinematically linking multiple seismic domains |
US20090257308A1 (en) * | 2008-04-11 | 2009-10-15 | Dimitri Bevc | Migration velocity analysis methods |
CN102109615A (zh) * | 2010-11-19 | 2011-06-29 | 恒泰艾普石油天然气技术服务股份有限公司 | 地震资料处理真地表条件下的振幅补偿方法和装置 |
CN102590857A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 真地表起伏叠前深度域双程波成像方法 |
CN102830431A (zh) * | 2012-08-14 | 2012-12-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 真地表射线追踪自适应插值方法 |
CN103472483A (zh) * | 2013-09-27 | 2013-12-25 | 中国石油化工股份有限公司 | 基于真地表或浮动基准面的速度建模方法 |
CN103869368A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种无表层调查资料约束的大炮初至综合建模静校正方法 |
CN105137481A (zh) * | 2015-08-24 | 2015-12-09 | 中国石油天然气集团公司 | 一种调整用于叠前偏移成像的地震数据能量一致性的方法 |
CN106443774A (zh) * | 2016-11-16 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种提高起伏地表叠前深度偏移成像精度的方法 |
CN106932823A (zh) * | 2015-12-31 | 2017-07-07 | 上海青凤致远地球物理地质勘探科技有限公司 | 一种光滑基准面选择技术 |
CN107783186A (zh) * | 2017-09-29 | 2018-03-09 | 中国石油天然气集团公司 | 一种近地表速度建模方法及装置 |
CN107870355A (zh) * | 2017-11-06 | 2018-04-03 | 西南交通大学 | 一种复杂地形条件下的克希霍夫型波束偏移方法 |
CN109085644A (zh) * | 2018-07-30 | 2018-12-25 | 中国石油化工股份有限公司 | 基于双射线走时的真地表成像方法 |
-
2019
- 2019-02-19 CN CN201910121989.4A patent/CN109917454B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5570321A (en) * | 1994-03-03 | 1996-10-29 | Atlantic Richfield Company | Seismic velocity model optimization method using simulated annearling to determine prestack travel-times |
US20050038604A1 (en) * | 2003-08-14 | 2005-02-17 | Marin Irina A. | Method and apparatus for kinematically linking multiple seismic domains |
US20090257308A1 (en) * | 2008-04-11 | 2009-10-15 | Dimitri Bevc | Migration velocity analysis methods |
CN102109615A (zh) * | 2010-11-19 | 2011-06-29 | 恒泰艾普石油天然气技术服务股份有限公司 | 地震资料处理真地表条件下的振幅补偿方法和装置 |
CN102590857A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 真地表起伏叠前深度域双程波成像方法 |
CN102830431A (zh) * | 2012-08-14 | 2012-12-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 真地表射线追踪自适应插值方法 |
CN103869368A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种无表层调查资料约束的大炮初至综合建模静校正方法 |
CN103472483A (zh) * | 2013-09-27 | 2013-12-25 | 中国石油化工股份有限公司 | 基于真地表或浮动基准面的速度建模方法 |
CN105137481A (zh) * | 2015-08-24 | 2015-12-09 | 中国石油天然气集团公司 | 一种调整用于叠前偏移成像的地震数据能量一致性的方法 |
CN106932823A (zh) * | 2015-12-31 | 2017-07-07 | 上海青凤致远地球物理地质勘探科技有限公司 | 一种光滑基准面选择技术 |
CN106443774A (zh) * | 2016-11-16 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种提高起伏地表叠前深度偏移成像精度的方法 |
CN107783186A (zh) * | 2017-09-29 | 2018-03-09 | 中国石油天然气集团公司 | 一种近地表速度建模方法及装置 |
CN107870355A (zh) * | 2017-11-06 | 2018-04-03 | 西南交通大学 | 一种复杂地形条件下的克希霍夫型波束偏移方法 |
CN109085644A (zh) * | 2018-07-30 | 2018-12-25 | 中国石油化工股份有限公司 | 基于双射线走时的真地表成像方法 |
Non-Patent Citations (2)
Title |
---|
张建磊 等: "GeoEast积分法叠前深度偏移处理技术", 《石油工业计算机应用》 * |
马召贵 等: "基于CPU+GPU联合计算真地表叠前时间偏移实用化研究", 《石油地球物理勘探》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624647A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111624647B (zh) * | 2020-06-05 | 2022-06-24 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111751875A (zh) * | 2020-07-07 | 2020-10-09 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp带限角度积分叠前时间偏移方法和装置 |
CN111751875B (zh) * | 2020-07-07 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp带限角度积分叠前时间偏移方法和装置 |
CN112270031A (zh) * | 2020-11-10 | 2021-01-26 | 黄河勘测规划设计研究院有限公司 | 一种多沙河流水库的双泥沙侵蚀基准面排沙系统 |
CN112270031B (zh) * | 2020-11-10 | 2021-06-15 | 黄河勘测规划设计研究院有限公司 | 一种多沙河流水库的双泥沙侵蚀基准面排沙系统 |
CN113031069A (zh) * | 2021-03-02 | 2021-06-25 | 吉林大学 | 一种岩溶地区的多元信息约束智能层析静校正方法 |
CN113568040A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
CN113568040B (zh) * | 2021-07-20 | 2024-01-26 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109917454B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109917454A (zh) | 基于双基准面的真地表叠前深度偏移成像方法及装置 | |
US6687618B2 (en) | Typing picks to horizons in migration velocity analysis | |
CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN107843922A (zh) | 一种基于地震初至波和反射波走时联合的层析成像方法 | |
CN106526674A (zh) | 一种三维全波形反演能量加权梯度预处理方法 | |
CN105301636B (zh) | 速度模型的建立方法和装置 | |
CN109581495A (zh) | 山地表层速度模型构建方法及系统 | |
EA008733B1 (ru) | Способ сейсморазведочных построений в сложных геологических формациях | |
CN105277978A (zh) | 一种确定近地表速度模型的方法及装置 | |
CN107817525B (zh) | 一种近地表常速填充方法及装置 | |
CN106199704B (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
CN107817516A (zh) | 基于初至波信息的近地表建模方法及系统 | |
CN104199103A (zh) | 转换波静校正方法及装置 | |
CN108445532B (zh) | 一种深度域反偏移方法及装置 | |
CN108845350A (zh) | 反演二维速度模型的方法及装置 | |
CN104614762B (zh) | 疏松砂岩气藏边界确定方法及装置 | |
CN108957539B (zh) | 层析偏移速度分析中的射线追踪方法及装置 | |
CN104597496B (zh) | 一种二维地震资料中速度的三维空间归位方法 | |
US9804281B2 (en) | Migration velocity analysis method for VSP data | |
CN109581494A (zh) | 叠前偏移方法及装置 | |
CN109188522A (zh) | 速度场构建方法及装置 | |
CN107193047A (zh) | 非地表一致性时变动态剩余静校正方法和装置 | |
CN107229074A (zh) | 低勘探程度海相探区的变速成图方法和装置 | |
Li et al. | Novel strategies for complex foothills seismic imaging—Part 1: Mega-near-surface velocity estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |