CN115793058A - 局部路径频变复旅行时的计算方法、装置、设备及介质 - Google Patents
局部路径频变复旅行时的计算方法、装置、设备及介质 Download PDFInfo
- Publication number
- CN115793058A CN115793058A CN202310079400.5A CN202310079400A CN115793058A CN 115793058 A CN115793058 A CN 115793058A CN 202310079400 A CN202310079400 A CN 202310079400A CN 115793058 A CN115793058 A CN 115793058A
- Authority
- CN
- China
- Prior art keywords
- ray
- grid
- travel time
- frequency
- equation
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种局部路径频变复旅行时的计算方法、装置、设备及介质,涉及油气勘探领域,方法包括:获取目标区域的地震速度模型、品质因子模型、炮点位置及对应的地震波速度,并对目标区域进行规则网格化以得到网格;基于地震速度模型构建程函方程;求解程函方程,得到横向方向和深度方向射线参数及非吸收衰减旅行时;基于横向方向和深度方向射线参数以及炮点位置,确定网格的射线传播角度和长度;基于品质因子模型及网格的射线传播角度和长度,确定耗散旅行时;基于非吸收衰减旅行时和耗散旅行时,得到频变复旅行时。基于此,本发明实现了含流体吸收衰减介质的频变复旅行时的计算,进而能实现含流体吸收衰减介质偏移成像和层析成像。
Description
技术领域
本发明涉及油气勘探领域,尤其涉及一种局部路径频变复旅行时的计算方法、装置、设备及介质。
背景技术
油气资源储存在地下岩石的孔隙中,而地下岩石对应的介质类型可视作为含流体吸收衰减介质类型,而在此种类型的介质下的旅行时存在频变和衰减,故此种旅行称之为频变复旅行时。
而频变复旅行时的计算因不满足最小旅行时原则,故无法采用常规的旅行时计算方法,也因此,当需采用射线方法观测含流体吸收衰减介质中的油气资源时,将因射线方法所需的频变复旅行时无法计算导致无法观测。也因此,如何计算介质内的频变复旅行时以完成观测是现阶段需克服的难点之一。
发明内容
有鉴于此,本发明提供一种局部路径频变复旅行时的计算方法、装置、设备及介质,以解决如何计算介质内的频变复旅行时以完成油气勘探的问题。
第一方面,本发明实施例提供一种局部路径频变复旅行时的计算方法,包括:
获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格;
基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况;
基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;
基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度;
基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减;
基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述程函方程包括:
所述基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时,包括:
对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
获取所述扰动量的离散差分格式,并基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度,包括:
基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
可选的,在本发明实施例提供的一种可行方式中,所述基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,包括:
基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,通过局部球面非线性插值计算方式和快速扫描计算方式确定深度方向和横向方向的网格的耗散旅行时,其中,所述炮点位置所在的网格的耗散旅行时为0。
可选的,在本发明实施例提供的一种可行方式中,所述网格对应的频变复旅行时的计算公式包括:
第二方面,本发明实施例提供一种局部路径频变复旅行时的计算装置,包括:
获取模块,用于获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格;
构建模块,用于基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况;
第一得到模块,用于基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;
确定模块,用于基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度;
第二得到模块,用于基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减;
第三得到模块,用于基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述程函方程包括:
所述第一得到模块,包括:
分解子模块,用于对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
扫描子模块,用于基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
参数获取子模块,用于基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述确定模块,包括:
角度得到子模块,用于基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
长度得到子模块,用于基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
第三方面,本发明实施例提供一种计算机设备,包括存储器以及处理器,存储器存储有计算机程序,计算机程序在处理器上运行时执行如第一方面中任一种公开的局部路径频变复旅行时的计算方法。
第四方面,本发明实施例提供一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序在处理器上运行时执行如第一方面中任一种公开的局部路径频变复旅行时的计算方法。
本发明提供的局部路径频变复旅行时的计算方法中,计算机设备首先获取目标区域的地震速度模型、品质因子模型、炮点位置及炮点位置对应的地震波速度,并对目标区域进行规则网格化操作以得到对应的多个网格;然后,基于地震速度模型构建程函方程;之后,基于炮点位置及地震波速度求解程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;接着,基于横向方向射线参数、深度方向射线参数及炮点位置,确定每个网格的射线传播角度和射线路径长度;随后,基于品质因子模型及每个网格的射线传播角度和射线路径长度,确定每个网格对应的耗散旅行时;最后,基于非吸收衰减旅行时和每个网格对应的耗散旅行时,得到每个网格对应的频变复旅行时。基于此,本发明实现了含流体吸收衰减介质的频变复旅行时(也即局部路径频变复旅行时)的计算,进而能实现含流体吸收衰减介质偏移成像和层析成像,从而满足了现阶段油气勘探的需要。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对本发明保护范围的限定。在各个附图中,类似的构成部分采用类似的编号。
图1示出了本发明实施例提供的局部路径频变复旅行时的计算方法的流程示意图;
图2示出了本发明实施例提供的地震速度模型图;
图3示出了本发明实施例体提供的品质因子模型图;
图4示出了本发明实施例提供的非吸收衰减旅行时示意图;
图5示出了本发明实施例提供的耗散旅行时示意图;
图6示出了x=2500m且z=1200m位置的频变复旅行时的实部示意图;
图7示出了x=2500m且z=1200m位置的频变复旅行时的虚部示意图;
图8示出了本发明实施例提供的局部路径频变复旅行时的计算装置的结构示意图。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
在下文中,可在本发明的各种实施例中使用的术语“包括”、“具有”及其同源词仅意在表示特定特征、数字、步骤、操作、元件、组件或前述项的组合,并且不应被理解为首先排除一个或更多个其它特征、数字、步骤、操作、元件、组件或前述项的组合的存在或增加一个或更多个特征、数字、步骤、操作、元件、组件或前述项的组合的可能性。
此外,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
除非另有限定,否则在这里使用的所有术语(包括技术术语和科学术语)具有与本发明的各种实施例所属领域普通技术人员通常理解的含义相同的含义。所述术语(诸如在一般使用的词典中限定的术语)将被解释为具有与在相关技术领域中的语境含义相同的含义并且将不被解释为具有理想化的含义或过于正式的含义,除非在本发明的各种实施例中被清楚地限定。
为更清楚地说明本发明提供的各个实施例,故在此对现有技术进行简要说明,具体如下:
在油气勘探领域中,基于地震波正演技术的地震勘探方法是较为常见的油气勘探方式之一。而地震波正演技术是通过波形、振幅和旅行时等地震波信息来模拟地震波在真实地下介质中的传播情况。
进一步的,地震波正演技术主要分为波动方程方法、有限频方法及射线方法三大类,其中,射线方法因计算效率较高,且能够适应大部分地下沉积介质,故为较常用的地震波正演技术。而在射线方法主要是模拟地震波的振幅和旅行时信息,其中的旅行时信息常用于层析成像和叠前深度成像等成像系统。
更进一步的,当地下岩石的孔隙储存有油气资源时,则此种地下沉积介质的介质类型可视作含流体吸收衰减介质类型,而此介质类型中传播的地震波存在频变和衰减,故其对应的旅行时信息(也即频变复旅行时)将表征地震波的频变和衰减。
但是,常规旅行时的计算方法多是基于最小旅行时原则来完成,而因频变复旅行时不满足最小旅行时原则,导致频变复旅行时无法直接采用常规旅行时的计算方法来得到。
并且,因含流体吸收衰减介质类型的层析成像和叠前深度成像中需采用频变复旅行时,故导致现有技术下的频变复旅行时难以通过射线方法获取,使得油气勘探受阻。
实施例1
参照图1,示出了本发明实施例提供的第一种局部路径频变复旅行时的计算方法的流程示意图,本发明实施例提供的局部路径频变复旅行时的计算方法包括:
步骤S110,获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格。
可以理解的是,地震速度模型用于描述在不同深度和不同横向位置下的地震波传播速度,品质因子模型用于描述在不同深度和不同横向位置下的地震波传播传播介质的品质因子,炮点位置代表地震波被激发的点位,炮点位置对应的地震波速度表示地震波被激发后其在炮点位置的速度。为更好地说明本发明实施例提供的地震速度模型和品质因子模型,请参考图2和图3,图2示出了本发明实施例提供的地震速度模型图,图3示出了本发明实施例体提供的品质因子模型图。
还可以理解的是,本发明实施例中的目标区域即指代介质类型为含流体吸收衰减介质类型的区域,也指代需进行油气勘测的区域。
需说明的是,本发明实施例中的计算机设备将基于射线方法完成的频变复旅行时的计算。并且,为使获取到的频变复旅行时能较为容易地在成像和层析中使用,故将本发明实施例的计算机设备将基于规则网格进行射线追踪。也因此,本发明实施例的计算机设备将对待检测区域(也即目标区域)进行规则网格化,以确定目标区域对应的多个规则网格,从而在后续步骤中确定每个网格中的射线路径以得到频变复旅行时。
步骤S120,基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况。
也即,本发明实施例的计算机设备将基于程函方程确定从炮点位置产生的地震波在不同网格下的非吸收衰减旅行时,也即未因介质类型而产生频变和衰减的旅行时信息。
可选的,在本发明实施例提供的一种可行方式中,前述程函方程包括:
步骤S130,基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
也即,本发明实施例的计算机设备将基于预设代码/程序对程函方程进行求解,从而得到横向方向射线参数和深度方向射线参数,进而能基于横向方向射线参数和深度方向射线参数确定射线在每个网格的射线传播角度。
同时,计算机设备还将得到非吸收衰减旅行时,从而确定任意一个网格下的不因储存介质而发生频变和衰减的旅行时信息。
可选的,在一种可行方式中,本发明实施例中的计算机设备内存储有预设程序/代码,其用于对程函方程执行有限差分求解方法,进而得到前述横向方向射线参数、前述深度方向射线参数及前述非吸收衰减旅行时。
可选的,计算机设备基于有限差分求解方法进行求解时,因炮点位置附近的等时面曲率较大,易使得差分求解结果引起较大误差,故在本发明实施例提供的一种可行方式中,所述程函方程包括:
进而,上述步骤S130,包括:
对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
获取所述扰动量的离散差分格式,并基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
也即,本发明实施例为降低因等时面曲率较大而导致的误差,故于程函方程中引入了因式分解,以得到非吸收衰减旅行时的解析部分和扰动量,具体可参考下式:
需说明的是,观测点/计算点可视作网格中心的位置。
而为求解扰动量,本发明实施例中的计算机设备还采用了扰动量的离散差分格式,具体可参考下式:
其中:
进一步的,为求解程函方程,本发明实施例中的计算机设备将以采用快速扫描计算方式对观测点的横向方向和深度方向分别进行预设次数次来回扫描运算以完成一次高斯-塞德尔迭代,并当迭代到扰动量稳定不变时停止迭代,进而基于稳定的扰动量完成求解。在一种可行方式中,上述预设次数为4。
由此,本发明实施例的计算机设备将基于前述解析部分和前述扰动量的离散差分格式求解所述程函方程,以降低由等时面曲率导致的误差。
步骤S140,基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度。
也即,本发明实施例将基于横向方向射线参数和深度方向射线参数确定每个网格的射线传播角度,并在确定射线传播角度后,基于射线传播角度确定网格内的局部射线长度,也即射线路径长度。
可选的,在本发明实施例提供的一种可行方式中,上述步骤S140具体包括:
基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
不难想到的是,射线传播角度的计算可参考下式:
可以理解的是,在程函方程的求解过程中,射线将从起始点指向当前计算点(也即网格位置),故射线的终点坐标为当前计算点,而网格内射线的起点坐标位于网格的边界。进而,基于已知晓的网格位置和射线传播角度,即可确定出网格内的射线传播路径(也即射线路径长度)。可选的,在本发明实施例提供的一种可行方式中,射线路径长度计算过程包括:(1)确定网格内射线的起点坐标。(2)基于起点坐标和终点坐标构建网格内的射线路径,也即局部射线路径。
还可以理解的是,因不同网格内的射线传播角度不同,故各个网格内射线的起点坐标的确定方式亦存在差异,如当射线传播角度的绝对值小于等于45度且大于等于0度时(也即0°≤||≤45°),则起点坐标的横向坐标和深度坐标的计算过程为:
还可以理解的是,在确定起点坐标和终点坐标后,网格内的射线路径长度可表达为:
步骤S150,基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减。
可以理解的是,相比于不描述频散和衰减的非吸收衰减旅行时,本发明实施例中的耗散旅行时将由地下流体引起的地震波频散效应和振幅衰减。进而在后续过程中将基于非吸收衰减旅行时和耗散旅行时完成各个网格内的频变复旅行时(也即局部路径频变复旅行时)的计算。
可选的,在本发明实施例提供的一种可行方式中,上述步骤S160具体包括:
基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,通过局部球面非线性插值计算方式和快速扫描计算方式确定深度方向和横向方向的网格的耗散旅行时,其中,所述炮点位置所在的网格的耗散旅行时为0。
也即,本发明实施例将基于扫描计算方式,对横向和深度方向分别来回扫描计算,从而各个网格的耗散旅行时。
更为具体的,本发明的申请人还考虑到等时面在局部区域具有球面非线性特征,故采用基于传播角度取值的两点非线性插值法以及扫描计算方式来确定各个网格的耗散旅行时。不难理解的是,确定各个网格耗散旅行时的过程为一种递归过程,进而,将炮点位置所在的网格的耗散旅行时为0以作为初始值后,本发明实施例中的计算机设备将基于以初始值完成递归计算,最终得到所有网格点的耗散旅行时。
进一步的,在递归计算的每一步中,当前网格位置的耗散旅行时依赖于此网格位置周围的其它网格的耗散旅行时,其它网格的耗散旅行时将由基于对应的递归过程得到。
示范性的,设在一个当前网格内,射线传播角度的绝对值小于等于45°并且大于等于0°,射线的起点左右相邻的两个网格的横向坐标分别为和,两个网格对应的耗散旅行时分别为和,和已通过对应的递归过程计算得到,则当前网格对应的起始耗散旅行时(也即局部射线起始耗散旅行时)为:
而当射线传播角度的绝对值小于等于90°且大于45°时,射线的起点上下相邻的两个网格的深度坐标分别为和,两个网格对应的耗散旅行时分别为和,和已通过对应的递归过程计算得到,则当前网格对应的起始耗散旅行时为:
步骤S160,基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
也即,本发明实施例将利用非吸收衰减旅行时和各个网格对应的耗散旅行时,确定非吸收衰减旅行时在发生衰减和频变后的相应结果,即频变复旅行时。可选的,在本发明实施例提供的一种可行方式中,所述网格对应的频变复旅行时的计算公式包括:
基于此,本发明实施例实现了含流体吸收衰减介质的频变复旅行时(也即局部路径频变复旅行时)的计算,进而能实现含流体吸收衰减介质偏移成像和层析成像,从而满足了现阶段油气勘探的需要。
为更清楚地说明本发明实施例提供的局部路径频变复旅行时的计算方法,故在此给出某地区的局部路径频变复旅行时的计算实例,具体如下:
某地区的地震速度模型为v(x,z),品质因子模型为q(x,z),v(x,z)和q(x,z)可参考图2和图3,此地区的炮点位于x=2500m且z=0m的位置,炮点位置对应的地震波速度为1500m/s。
基于v(x,z)和炮点位置构建程函方程,得到式(1),即:
因炮点位置及炮点位置附近的等时面曲率较大,直接对炮点位置进行有限差分将引入较大误差,故为降低等时面曲率带来的误差,引入非吸收衰减旅行时T对应的因式分解部分,则得到式(2),即:
将式(2)带入式(1),得到式(3),即:
其中:
因此地区被划分为规则网格,故每个网格为大小相同的矩形,进而计算点分布于各个网格点上,也因此,计算点即为规则网格点。同时,针对任一个网格,网格内的(局部)射线的终点为当前计算点,射线穿过网格到达当前计算点,也即,射线的起点位于网格的边界。
由此,将式(4)和(5)带入式(3),得到式(6),即:
需说明的是,在计算的过程中,首先将炮点位置的非吸收衰减旅行时设置为0;接着,采用快速扫描方式朝着计算点的横向和深度方向分别进行4次扫描以完成一次高斯-塞德尔迭代;当迭代到不发生改变后,也即,得到的稳定值后,停止迭代,并利用的稳定值计算横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
也因此,网格内的局部射线路径,即射线路径长度d1为:
耗散旅行时采用快速扫描计算方式完成,具体的,将炮点位置处的耗散旅行时设置为0以作为初始值,在基于初始值,分别对计算点的横方向和深度方向进行递归的来回扫描计算,最终得到所有网格点的耗散旅行时。在递归计算的每一步,当前计算点的耗散旅行时未知,而网格周围的其它网格点的耗散旅行时已经由先前完成的递归运算得到,故当前网格内射线的起点的相邻网格点的耗散旅行时已知。同时,等时面在局部区域具有球面非线性特征,由此,当0°≤||≤45°时,当前网格点的左右相邻的两个网格点的横向坐标分别为和两个网格对应的耗散旅行时分别为和,则当前网格的起始耗散旅行时(局部射线起始耗散旅行时)为:
为更好地说明此实例中的耗散旅行时,请参考图5,示出了本发明实施例提供的耗散旅行时示意图。
由此,针对任意一个网格,其对应的局部路径频变复旅行时为:
为更好地说明此实例中的频变复旅行时,请参考图6和图7,图6示出了x=2500m且z=1200m位置的频变复旅行时的实部示意图,图7示出了x=2500m且z=1200m位置的频变复旅行时的虚部示意图。
实施例2
与本发明实施例1提供的局部路径频变复旅行时的计算方法相对应的,本发明实施例2还提供一种局部路径频变复旅行时的计算装置,参照图8,示出了本发明实施例提供的局部路径频变复旅行时的计算装置的结构示意图,本发明实施例提供的局部路径频变复旅行时的计算装置200,包括:
获取模块210,用于获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格;
构建模块220,用于基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况;
第一得到模块230,用于基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;
确定模块240,用于基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度;
第二得到模块250,用于基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减;
第三得到模块260,用于基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述程函方程包括:
所述第一得到模块,包括:
分解子模块,用于对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
扫描子模块,用于基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
参数获取子模块,用于基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
可选的,在本发明实施例提供的一种可行方式中,所述确定模块,包括:
角度得到子模块,用于基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
长度得到子模块,用于基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
可选的,在本发明实施例提供的一种可行方式中,所第二得到模块还用于基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,通过局部球面非线性插值计算方式和快速扫描计算方式确定深度方向和横向方向的网格的耗散旅行时,其中,所述炮点位置所在的网格的耗散旅行时为0。
可选的,在本发明实施例提供的一种可行方式中,所述网格对应的频变复旅行时的计算公式包括:
本申请实施例提供的局部路径频变复旅行时的计算装置200能够实现实施例1对应的局部路径频变复旅行时的计算方法的各个过程,且能达到相同的技术效果,为避免重复,这里不再赘述。
本发明实施例还提供一种计算机设备,包括存储器以及处理器,存储器存储有计算机程序,计算机程序在处理器上运行时执行如实施例1中所述的局部路径频变复旅行时的计算方法。
本发明实施例还提供一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序在处理器上运行时执行如实施例1中所述的局部路径频变复旅行时的计算方法。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和结构图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,结构图和/或流程图中的每个方框、以及结构图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明各个实施例中的各功能模块或单元可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或更多个模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是智能手机、个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种局部路径频变复旅行时的计算方法,其特征在于,包括:
获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格;
基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况;
基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;
基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度;
基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减;
基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
2.根据权利要求1所述的局部路径频变复旅行时的计算方法,其特征在于,所述程函方程包括:
所述基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时,包括:
对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
获取所述扰动量的离散差分格式,并基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
3.根据权利要求1所述的局部路径频变复旅行时的计算方法,其特征在于,所述基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度,包括:
基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
4.根据权利要求1所述的局部路径频变复旅行时的计算方法,其特征在于,所述基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,包括:
基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,通过局部球面非线性插值计算方式和快速扫描计算方式确定深度方向和横向方向的网格的耗散旅行时,其中,所述炮点位置所在的网格的耗散旅行时为0。
6.一种局部路径频变复旅行时的计算装置,其特征在于,包括:
获取模块,用于获取目标区域的地震速度模型、品质因子模型、炮点位置及所述炮点位置对应的地震波速度,并对所述目标区域进行规则网格化操作,得到所述目标区域对应的多个网格;
构建模块,用于基于所述地震速度模型构建程函方程,其中,所述程函方程表示从所述炮点位置产生的地震波沿不同射线方向在各个所述网格的传播情况;
第一得到模块,用于基于所述炮点位置及所述地震波速度求解所述程函方程,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时;
确定模块,用于基于所述横向方向射线参数、所述深度方向射线参数及所述炮点位置,确定每个所述网格的射线传播角度和射线路径长度;
第二得到模块,用于基于所述品质因子模型及每个所述网格的所述射线传播角度和所述射线路径长度,确定每个所述网格对应的耗散旅行时,其中,所述耗散旅行时表示由地下流体引起的地震波频散效应和振幅衰减;
第三得到模块,用于基于所述非吸收衰减旅行时和每个所述网格对应的耗散旅行时,得到每个所述网格对应的频变复旅行时。
7.根据权利要求6所述的局部路径频变复旅行时的计算装置,其特征在于,所述程函方程包括:
所述第一得到模块,包括:
分解子模块,用于对所述程函方程中的非吸收衰减旅行时进行因式分解,得到所述非吸收衰减旅行时的解析部分和扰动量,其中,所述解析部分表示网格的位置到所述炮点位置的距离与所述地震波速度的比值;
扫描子模块,用于基于高斯-塞德尔迭代法对观测点的深度方向和横向方向进行预设次数的快速扫描计算,得到所述扰动量对应的稳定值;
参数获取子模块,用于基于所述解析部分和所述扰动量的稳定值,得到横向方向射线参数、深度方向射线参数及非吸收衰减旅行时。
8.根据权利要求6所述的局部路径频变复旅行时的计算装置,其特征在于,所述确定模块,包括:
角度得到子模块,用于基于所述横向方向射线参数和所述深度方向射线参数的比值的反正切函数值,得到每个所述网格的射线传播角度;
长度得到子模块,用于基于每个所述网格对应的位置及每个所述网格内的射线传播角度,得到每个所述网格内的射线路径长度。
9.一种计算机设备,其特征在于,包括存储器以及处理器,所述存储器存储有计算机程序,所述计算机程序在所述处理器上运行时执行如权利要求1-5任一项所述的局部路径频变复旅行时的计算方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序在处理器上运行时执行如权利要求1-5任一项所述的局部路径频变复旅行时的计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310079400.5A CN115793058B (zh) | 2023-02-08 | 2023-02-08 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310079400.5A CN115793058B (zh) | 2023-02-08 | 2023-02-08 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115793058A true CN115793058A (zh) | 2023-03-14 |
CN115793058B CN115793058B (zh) | 2023-04-18 |
Family
ID=85430457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310079400.5A Active CN115793058B (zh) | 2023-02-08 | 2023-02-08 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115793058B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5062086A (en) * | 1990-08-27 | 1991-10-29 | Conoco Inc. | Calculation of raypaths and wavepaths from traveltime tables for the tomographic estimation of transmission velocities |
US20080075346A1 (en) * | 2006-09-27 | 2008-03-27 | Ziosoft, Inc. | Image processing method |
US20110075516A1 (en) * | 2009-09-25 | 2011-03-31 | Halliburton Energy Services, Inc. | Seismic Imaging Systems and Methods Employing Tomographic Migration-Velocity Analysis Using Common Angle Image Gathers |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN107966729A (zh) * | 2016-10-19 | 2018-04-27 | 中国石油化工股份有限公司 | 一种三维tti介质射线追踪方法及系统 |
CN109100786A (zh) * | 2018-06-21 | 2018-12-28 | 中国石油天然气集团有限公司 | 深度域品质因子的确定方法和装置 |
CN109725345A (zh) * | 2018-11-15 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种初至波正演模拟方法及装置 |
CN109917454A (zh) * | 2019-02-19 | 2019-06-21 | 中国石油天然气集团有限公司 | 基于双基准面的真地表叠前深度偏移成像方法及装置 |
US20190302296A1 (en) * | 2018-03-27 | 2019-10-03 | Mehdi Aharchaou | Directional Q Compensation with Sparsity Constraints and Preconditioning |
CN110658558A (zh) * | 2019-09-25 | 2020-01-07 | 中国石油化工股份有限公司 | 吸收衰减介质叠前深度逆时偏移成像方法及系统 |
CN111123354A (zh) * | 2019-12-30 | 2020-05-08 | 中国石油大学(北京) | 基于频变反射振幅衰减预测致密气层的方法及设备 |
CN112904415A (zh) * | 2021-01-20 | 2021-06-04 | 中国石油大学(北京) | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
CN112987088A (zh) * | 2021-02-22 | 2021-06-18 | 成都理工大学 | 一种渗流介质地震横波数值模拟和成像方法 |
CN115201896A (zh) * | 2022-02-17 | 2022-10-18 | 成都理工大学 | 吸收衰减介质逆时偏移方法、装置、成像方法、介质 |
-
2023
- 2023-02-08 CN CN202310079400.5A patent/CN115793058B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5062086A (en) * | 1990-08-27 | 1991-10-29 | Conoco Inc. | Calculation of raypaths and wavepaths from traveltime tables for the tomographic estimation of transmission velocities |
US20080075346A1 (en) * | 2006-09-27 | 2008-03-27 | Ziosoft, Inc. | Image processing method |
US20110075516A1 (en) * | 2009-09-25 | 2011-03-31 | Halliburton Energy Services, Inc. | Seismic Imaging Systems and Methods Employing Tomographic Migration-Velocity Analysis Using Common Angle Image Gathers |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN107966729A (zh) * | 2016-10-19 | 2018-04-27 | 中国石油化工股份有限公司 | 一种三维tti介质射线追踪方法及系统 |
US20190302296A1 (en) * | 2018-03-27 | 2019-10-03 | Mehdi Aharchaou | Directional Q Compensation with Sparsity Constraints and Preconditioning |
CN109100786A (zh) * | 2018-06-21 | 2018-12-28 | 中国石油天然气集团有限公司 | 深度域品质因子的确定方法和装置 |
CN109725345A (zh) * | 2018-11-15 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种初至波正演模拟方法及装置 |
CN109917454A (zh) * | 2019-02-19 | 2019-06-21 | 中国石油天然气集团有限公司 | 基于双基准面的真地表叠前深度偏移成像方法及装置 |
CN110658558A (zh) * | 2019-09-25 | 2020-01-07 | 中国石油化工股份有限公司 | 吸收衰减介质叠前深度逆时偏移成像方法及系统 |
CN111123354A (zh) * | 2019-12-30 | 2020-05-08 | 中国石油大学(北京) | 基于频变反射振幅衰减预测致密气层的方法及设备 |
CN112904415A (zh) * | 2021-01-20 | 2021-06-04 | 中国石油大学(北京) | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
CN112987088A (zh) * | 2021-02-22 | 2021-06-18 | 成都理工大学 | 一种渗流介质地震横波数值模拟和成像方法 |
CN115201896A (zh) * | 2022-02-17 | 2022-10-18 | 成都理工大学 | 吸收衰减介质逆时偏移方法、装置、成像方法、介质 |
Non-Patent Citations (4)
Title |
---|
DIRK GAJEWSKI;杨勤勇;: "用射线法研究弱衰减各向异性介质中的矢量波场" * |
杜增利;施泽进;尹成;刘文革;姜庆国;: "折射初至波射线追踪方法研究" * |
王玉风;孟祥宁;安岩;: "Q值补偿方法研究" * |
黄兴国;孙建国;孙章庆;王乾龙;: "基于复程函方程和改进的快速推进法的复旅行时计算方法" * |
Also Published As
Publication number | Publication date |
---|---|
CN115793058B (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2690373C (en) | Method for velocity analysis using waveform inversion in laplace domain for geophysical imaging | |
Lee et al. | Parabolic equation development in recent decade | |
US6721694B1 (en) | Method and system for representing the depths of the floors of the oceans | |
US8731838B2 (en) | Fresnel zone fat ray tomography | |
US7460437B2 (en) | Seismic data processing method and system for migration of seismic signals incorporating azimuthal variations in the velocity | |
EA025919B1 (ru) | Системы и способы моделирования трехмерных геологических структур | |
EA008733B1 (ru) | Способ сейсморазведочных построений в сложных геологических формациях | |
US20140257700A1 (en) | System and method for estimating uncertainty for geophysical gridding routines lacking inherent uncertainty estimation | |
CN114428324B (zh) | 叠前高角度快速傅里叶变换地震成像方法、系统、设备 | |
CN104755961A (zh) | 基于蒙特卡洛反向投影的波束反演 | |
CN106896403A (zh) | 弹性高斯束偏移成像方法和系统 | |
US20170045634A1 (en) | Method of using a parabolic equation model for range-dependent seismo-acoustic problems | |
WO2014057440A1 (en) | Wave equation processing | |
US20140129188A1 (en) | System and method for analysis of seismic images | |
CN111624658A (zh) | 深度域成像模拟方法及系统 | |
US20160018541A1 (en) | System and method for rock property estimation of subsurface geologic volumes | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
CN102830431A (zh) | 真地表射线追踪自适应插值方法 | |
CN115793058B (zh) | 局部路径频变复旅行时的计算方法、装置、设备及介质 | |
AU739128B2 (en) | A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration | |
CN104570091B (zh) | 一种获得初至波射线的方法和装置 | |
CN112596103A (zh) | 射线追踪方法、装置和电子设备 | |
CN108254787B (zh) | 波的走时获得方法及装置、成像方法及装置 | |
US8259532B2 (en) | Hybrid modeling in the tau-p domain | |
RU2794968C1 (ru) | Моделирование сейсмических скоростей |
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 |