CN112946745A - 一种基于引力拟合的转换波剩余静校正方法及系统 - Google Patents
一种基于引力拟合的转换波剩余静校正方法及系统 Download PDFInfo
- Publication number
- CN112946745A CN112946745A CN201911265677.7A CN201911265677A CN112946745A CN 112946745 A CN112946745 A CN 112946745A CN 201911265677 A CN201911265677 A CN 201911265677A CN 112946745 A CN112946745 A CN 112946745A
- Authority
- CN
- China
- Prior art keywords
- wave
- converted
- static correction
- analysis
- residual static
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 125
- 230000003068 static effect Effects 0.000 title claims abstract description 124
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000005484 gravity Effects 0.000 title claims abstract description 49
- 238000004458 analytical method Methods 0.000 claims abstract description 130
- 238000001514 detection method Methods 0.000 claims abstract description 89
- 238000005381 potential energy Methods 0.000 claims abstract description 18
- 230000006870 function Effects 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 239000002245 particle Substances 0.000 description 20
- 238000010586 diagram Methods 0.000 description 13
- 238000003860 storage Methods 0.000 description 12
- 238000012545 processing Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于引力拟合的转换波剩余静校正方法及系统,所述方法包括:根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位;根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位;选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正,本发明基于共接收点道集通过引力拟合函数求取转换波剩余静校正量,改善转换波剩余静校正量的求取准确度。
Description
技术领域
本发明涉及石油地震勘探中地震数据处理技术领域,尤其涉及一种基于引力拟合的转换波剩余静校正方法及系统。
背景技术
近年来,多波勘探得到了越来越多的重视。特别是多波勘探在岩性分析、流体识别、裂缝检测和气云区成像等方面显示出其独特优势,在实际工业生产中得到了越来越多的重视。转换波剩余静校正量的求取是转换波处理的关键步骤,由于近地表的横波速度很低,导致转换波接收点的剩余静校正量比传统纵波的剩余静校正量大得多。横波的静校正量通常要大于相同位置纵波静校正量的2~10倍,与纵波相比,横波实际上不受近地表潜水面的影响,因此纵波和横波静校正量不存在简单的比例关系,使用纵波静校正量乘以某一系数来对横波静校正量进行近似通常是不能用于精细成像的,同时,沿用常规纵波剩余静校正算法进行转换波剩余静校正也具有一定的局限性,例如往往会出现不稳定的“周期跳跃”现象。
发明内容
本发明的一个目的在于提供一种基于引力拟合的转换波剩余静校正方法,基于共接收点道集通过引力拟合函数求取转换波剩余静校正量,改善转换波剩余静校正量的求取准确度。本发明的另一个目的在于提供一种基于引力拟合的转换波剩余静校正系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。
为了达到以上目的,本发明一方面公开了一种基于引力拟合的转换波剩余静校正方法,包括:
根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位;
根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位;
选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。
优选的,所述根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位具体包括:
在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴;
对多段所述同相轴插值得到所述第一分析层位。
优选的,所述根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位具体包括:
通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面;
确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
优选的,所述方法进一步包括:
通过所述第二分析层位对所述转换波叠加剖面进行层拉平。
优选的,所述选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量具体包括:
确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串;
依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
本发明还公开了一种基于引力拟合的转换波剩余静校正系统,包括:
纵波分析层位确定单元,用于根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位;
横波分析层位确定单元,用于根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位;
剩余静校正量确定单元,用于选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。
优选的,所述纵波分析层位确定单元具体用于在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴,对多段所述同相轴插值得到所述第一分析层位。
优选的,所述横波分析层位确定单元具体用于通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面,确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
优选的,所述横波分析层位确定单元进一步用于通过所述第二分析层位对所述转换波叠加剖面进行层拉平。
优选的,所述剩余静校正量确定单元具体用于确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串,依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
本发明还公开了一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,
所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质,其上存储有计算机程序,
该程序被处理器执行时实现如上所述方法。
本发明将纵波叠加剖面转换为转换波叠加剖面,并确定转换波叠加剖面上与纵波叠加剖面的第一分析层位对应的第二分析层位,在第二分析层位上设置预设时间间隔,将第二分析层位上预设时间间隔中的检测数据赋予质量,形成第二分析层位的质子串,然后根据最小势能原理分析质子串得到转换波剩余静校正量,本发明在剩余静校正计算理论的基础上结合最小势能原理,得到转换波剩余静校正量,可改善现有转换波剩余静校正量确定过程中易出现的不稳定的“周期跳跃”现象,且确定方法简单和稳定,计算量低。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体实施例的流程图之一;
图2示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体实施例的流程图之二;
图3示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体实施例的流程图之三;
图4示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体实施例的流程图之四;
图5示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体实施例的流程图之五;
图6示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体例子中纵向叠加剖面的示意图;
图7示出图6中纵向叠加剖面添加扰动后的示意图;
图8示出本发明一种基于引力拟合的转换波剩余静校正方法一个具体例子中静校正后的纵向叠加剖面的示意图;
图9示出本发明一种基于引力拟合的转换波剩余静校正系统一个具体实施例的结构图;
图10示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
根据本发明的一个方面,本实施例公开了一种基于引力拟合的转换波剩余静校正方法。如图1所示,本实施例中,所述方法包括:
S100:根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位。
S200:根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
S300:选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。其中,在地表一致性的假设下,地下构造不复杂且接收点静校正量很大的情况下,剩余静校正量可以近似等于接收点的静校正量,通过剩余静校正量对接收点进行静校正,能够很好的对近地表的检测数据进行修正,以减少误差。其中,地下构造不复杂和静校正量很大的情况的确定为本领域技术人员可以根据经验得到,在实际应用中,也可通过其他计算方式确定下地构造的复杂程度和静校正量大小的程度。
本发明将纵波叠加剖面转换为转换波叠加剖面,并确定转换波叠加剖面上与纵波叠加剖面的第一分析层位对应的第二分析层位,在第二分析层位上设置预设时间间隔,将第二分析层位上预设时间间隔中的检测数据赋予质量,形成第二分析层位的质子串,然后根据最小势能原理分析质子串得到转换波剩余静校正量,本发明在剩余静校正计算理论的基础上结合最小势能原理,采用引力拟合得到转换波剩余静校正量,可改善现有转换波剩余静校正量确定过程中易出现的不稳定的“周期跳跃”现象,且确定方法简单和稳定,计算量低。
在优选的实施方式中,如图2所示,所述S100具体可包括:
S110:在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴。
S120:对多段所述同相轴插值得到所述第一分析层位。
可以理解的是,在实际应用中,可获取共检波点的检测数据,根据检测数据形成纵波叠加剖面,在形成的纵波叠加剖面分段拾取满足预设条件的层位进行分析。优选的,通过采取反应波形特征和能量特征的同相轴,以提高分析结果的准确性。在实际应用中,根据精度和条件等方面的要求,也可采用其他预设条件拾取满足预设条件的层位。对于拾取的离散的同相轴层位,可进一步通过插值的方式形成连续的层位数据得到第一分析层位。
在优选的实施方式中,如图3所示,所述S200具体包括:
S210:通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面。
S220:确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
可以理解的是,需在转换波叠加剖面上分析转换波剩余静校正量,可根据纵波和转换波速度的比值将纵波叠加剖面转换为转换波叠加剖面,进一步的可得到转换波叠加剖面上与第一分析层位对应的第二分析层位,转换波叠加剖面上的第二分析层位分析引力分析的结果即为转换波剩余静校正量,可根据该转换波剩余静校正量对共检波点进行静校正。
在优选的实施方式中,如图4所示,所述S200进一步包括:
S230:通过所述第二分析层位对所述转换波叠加剖面进行层拉平。具体的,通过根据第一分析层位转换得到的第二分析层位对转换波叠加剖面进行层拉平,可对转换波叠加剖面进行地层校正,以提高转换波叠加剖面的准确度。
在优选的实施方式中,如图5所示,所述S300具体包括:
S310:确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串。
S320:依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
可以理解的是,预设时间间隔即为预设的时窗,即确定一个时间范围,在第二分析层位的多个位置上分别确定一个时间范围内的检测数据,然后将多个位置的时窗内的检测数据根据位置形成质子串。
其中,在一个具体例子中,对于质子串中的每个检测数据,当作一个真实存在的质点,且定义质点的质量为Mi,Mi可取该检测数据的振幅能量值,当然也可以取其他对应值进行表征。
每个质点还对应于转换波叠加剖面上的一个位置,可设置与转换波叠加剖面对应的坐标系,将质点的位置通过坐标进行表示。可选的,在一个具体例子中,坐标系可为二维坐标系,质点的位置坐标为(x,y)其中,x可表示分析道的道间距,y可表示时间差。在另一个具体例子中,坐标系也可为三维坐标系,质点的位置坐标为(x,y,z),其中,x可为质点在x方向的坐标,y可为质点在y方向的坐标,z可表示时间差。
在优选的实施方式中,在采用最小势能原理分析所述质子串得到转换波剩余静校正量的过程中,假设所有质点可在时间方向移动,即可在转换波叠加剖面与地层垂直的纵向移动。其中,每个质点在初始位置(x0,y0)都会受到一个弹性约束,假设弹性约束为ki,并且各质点之间均存在万有引力。根据最小势能原理,每个质点在初始位置上会在万有引力的驱动和弹性约束的限制下达到新的平衡位置,令这个位置坐标为(xi,yi)。进一步可形成新坐标位置与各道质点的引力公式。因yi可能移动的位置有N个,所以形成的引力公式也有N个。当引力最大时,质点达到平衡,所以求解形成的N个引力公式形成的引力方程组的最大引力值,可以求得该最大引力值下的质点位置坐标。质点达到新平衡位置的位置坐标与原始位置坐标的差就是该道的剩余静校正量,可将求得的检波点静校正量应用到转换波叠加剖面,对检波点进行静校正。
在另一个具体例子中,检测数据为二维数据时,坐标系为二维坐标系。需要引力拟合的质子串中的n个质点的位置坐标分别为Ai(xi,yi),i=1,2,...,n。令引力拟合得到的最大引力值时各质点的位置坐标为(xi,y′i),并用点集A’表示。即A’的坐标为拟合后的质点的坐标。此时,对于质点Ai具有力平衡公式:
其中,θij为点Ai和点Aj间连线与x轴所形成的夹角,G为引力常量。
令Ki=ki/G,带入上式并整理得到力平衡公式:
上述力平衡公式共有n(正整数)个,未知量y′i亦有n个。求解上述n个非线性方程组即可得到在引力和弹性力作用平衡下的拟合点坐标。优选的,可采用牛顿迭代法求解得到拟合点坐标,则拟合点坐标与初始位置坐标的差值即为转换波剩余静校正量。
在一个具体例子中,检测数据为三维数据时,坐标系为三维坐标系。需要引力拟合的质子串中的n个质点的位置坐标分别为Ai(xi,yi,zi),i=1,2,...,n。令引力拟合得到的最大引力值时各质点的位置坐标为(xi,yi',zi),则力平衡公式为:
同理,上述力平衡公式共有n(正整数)个,未知量y′i亦有n个。求解上述n个非线性方程组即可得到在引力和弹性力作用平衡下的拟合点坐标。优选的,可采用牛顿迭代法求解得到拟合点坐标,则拟合点坐标与初始位置坐标的差值即为转换波剩余静校正量。同理,通过该最小势能原理分析得到转换波剩余静校正量的方法可以进一步适用于任意维度空间下分析得到转换波剩余静校正量的方案。
在一个具体实例中,图6示出一个纵向叠加剖面,图7为图6中加入了扰动的数据形成的叠加剖面,图8示出通过本发明转换波剩余静校正方法进行静校正得到的校正后的叠加剖面,通过图中可以看出,本发明通过全波形剩余静校正处理,即把每一个地震道看作一个质点串,质点分布在该道的所有波峰和波谷上,波峰位置上的质点带正电荷,波谷位置上的质点带负电荷。与之前的引力拟合不同,这里用电荷间的引力和斥力代替之前的万有引力。通过拟合得到的每道时间差即为该道的残差时间。能够很好的改善转换波剩余静校正量的求取准确度,避免现有静校正过程中的“周期跳跃”现象。
基于相同原理,本实施例还公开了一种基于引力拟合的转换波剩余静校正系统。如图9所示,所述系统包括纵波分析层位确定单元11、横波分析层位确定单元12和剩余静校正量确定单元13。
其中,纵波分析层位确定单元11用于根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位。
横波分析层位确定单元12用于根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
剩余静校正量确定单元13用于选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。其中,在地表一致性的假设下,地下构造不复杂且接收点静校正量很大的情况下,剩余静校正量可以近似等于接收点的静校正量,通过剩余静校正量对接收点进行静校正,能够很好的对近地表的检测数据进行修正,以减少误差。其中,地下构造不复杂和静校正量很大的情况的确定为本领域技术人员可以根据经验得到,在实际应用中,也可通过其他计算方式确定下地构造的复杂程度和静校正量大小的程度。
本发明将纵波叠加剖面转换为转换波叠加剖面,并确定转换波叠加剖面上与纵波叠加剖面的第一分析层位对应的第二分析层位,在第二分析层位上设置预设时间间隔,将第二分析层位上预设时间间隔中的检测数据赋予质量,形成第二分析层位的质子串,然后根据最小势能原理分析质子串得到转换波剩余静校正量,本发明在剩余静校正计算理论的基础上结合最小势能原理,采用引力拟合得到转换波剩余静校正量,可改善现有转换波剩余静校正量确定过程中易出现的不稳定的“周期跳跃”现象,且确定方法简单和稳定,计算量低。
在优选的实施方式中,所述纵波分析层位确定单元11具体用于在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴,对多段所述同相轴插值得到所述第一分析层位。
可以理解的是,在实际应用中,可获取共检波点的检测数据,根据检测数据形成纵波叠加剖面,在形成的纵波叠加剖面分段拾取满足预设条件的层位进行分析。优选的,通过采取反应波形特征和能量特征的同相轴,以提高分析结果的准确性。在实际应用中,根据精度和条件等方面的要求,也可采用其他预设条件拾取满足预设条件的层位。对于拾取的离散的同相轴层位,可进一步通过插值的方式形成连续的层位数据得到第一分析层位。
在优选的实施方式中,所述横波分析层位确定单元12具体用于通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面,确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
可以理解的是,需在转换波叠加剖面上分析转换波剩余静校正量,可根据纵波和转换波速度的比值将纵波叠加剖面转换为转换波叠加剖面,进一步的可得到转换波叠加剖面上与第一分析层位对应的第二分析层位,转换波叠加剖面上的第二分析层位分析引力分析的结果即为转换波剩余静校正量,可根据该转换波剩余静校正量对共检波点进行静校正。
在优选的实施方式中,所述横波分析层位确定单元12进一步用于通过所述第二分析层位对所述转换波叠加剖面进行层拉平。具体的,通过根据第一分析层位转换得到的第二分析层位对转换波叠加剖面进行层拉平,可对转换波叠加剖面进行地层校正,以提高转换波叠加剖面的准确度。
在优选的实施方式中,所述剩余静校正量确定单元13具体用于确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串,依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
可以理解的是,预设时间间隔即为预设的时窗,即确定一个时间范围,在第二分析层位的多个位置上分别确定一个时间范围内的检测数据,然后将多个位置的时窗内的检测数据根据位置形成质子串。
其中,在一个具体例子中,对于质子串中的每个检测数据,当作一个真实存在的质点,且定义质点的质量为Mi,Mi可取该检测数据的振幅能量值,当然也可以取其他对应值进行表征。
每个质点还对应于转换波叠加剖面上的一个位置,可设置与转换波叠加剖面对应的坐标系,将质点的位置通过坐标进行表示。可选的,在一个具体例子中,坐标系可为二维坐标系,质点的位置坐标为(x,y)其中,x可表示分析道的道间距,y可表示时间差。在另一个具体例子中,坐标系也可为三维坐标系,质点的位置坐标为(x,y,z),其中,x可为质点在x方向的坐标,y可为质点在y方向的坐标,z可表示时间差。
在优选的实施方式中,在采用最小势能原理分析所述质子串得到转换波剩余静校正量的过程中,假设所有质点可在时间方向移动,即可在转换波叠加剖面与地层垂直的纵向移动。其中,每个质点在初始位置(x0,y0)都会受到一个弹性约束,假设弹性约束为ki,并且各质点之间均存在万有引力。根据最小势能原理,每个质点在初始位置上会在万有引力的驱动和弹性约束的限制下达到新的平衡位置,令这个位置坐标为(xi,yi)。进一步可形成新坐标位置与各道质点的引力公式。因yi可能移动的位置有N个,所以形成的引力公式也有N个。当引力最大时,质点达到平衡,所以求解形成的N个引力公式形成的引力方程组的最大引力值,可以求得该最大引力值下的质点位置坐标。质点达到新平衡位置的位置坐标与原始位置坐标的差就是该道的剩余静校正量,可将求得的检波点静校正量应用到转换波叠加剖面,对检波点进行静校正。
在一个具体例子中,检测数据为二维数据时,坐标系为二维坐标系。需要引力拟合的质子串中的n个质点的位置坐标分别为Ai(xi,yi),i=1,2,...,n。令引力拟合得到的最大引力值时各质点的位置坐标为(xi,y′i),并用点集A’表示。即A’的坐标为拟合后的质点的坐标。此时,对于质点Ai具有力平衡公式:
其中,θij为点Ai和点Aj间连线与x轴所形成的夹角,G为引力常量。
令Ki=ki/G,带入上式并整理得到力平衡公式:
上述力平衡公式共有n(正整数)个,未知量y′i亦有n个。求解上述n个非线性方程组即可得到在引力和弹性力作用平衡下的拟合点坐标。优选的,可采用牛顿迭代法求解得到拟合点坐标,则拟合点坐标与初始位置坐标的差值即为转换波剩余静校正量。
在另一个具体例子中,检测数据为三维数据时,坐标系为三维坐标系。需要引力拟合的质子串中的n个质点的位置坐标分别为Ai(xi,yi,zi),i=1,2,...,n。令引力拟合得到的最大引力值时各质点的位置坐标为(xi,yi',zi),则力平衡公式为:
同理,上述力平衡公式共有n(正整数)个,未知量y′i亦有n个。求解上述n个非线性方程组即可得到在引力和弹性力作用平衡下的拟合点坐标。优选的,可采用牛顿迭代法求解得到拟合点坐标,则拟合点坐标与初始位置坐标的差值即为转换波剩余静校正量。同理,通过该最小势能原理分析得到转换波剩余静校正量的方法可以进一步适用于任意维度空间下分析得到转换波剩余静校正量的方案。
由于该系统解决问题的原理与以上方法类似,因此本系统的实施可以参见方法的实施,在此不再赘述。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备,具体的,计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述的由客户端执行的方法,或者,所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图10,其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图10所示,计算机设备600包括中央处理单元(CPU)601,其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中,还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605:包括键盘、鼠标等的输入部分606;包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607;包括硬盘等的存储部分608;以及包括诸如LAN卡,调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器610上,以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地,根据本发明的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括有形地包含在机器可读介质上的计算机程序,所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分609从网络上被下载和安装,和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (12)
1.一种基于引力拟合的转换波剩余静校正方法,其特征在于,包括:
根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位;
根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位;
选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。
2.根据权利要求1所述的转换波剩余静校正方法,其特征在于,所述根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位具体包括:
在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴;
对多段所述同相轴插值得到所述第一分析层位。
3.根据权利要求1所述的转换波剩余静校正方法,其特征在于,所述根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位具体包括:
通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面;
确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
4.根据权利要求3所述的转换波剩余静校正方法,其特征在于,所述方法进一步包括:
通过所述第二分析层位对所述转换波叠加剖面进行层拉平。
5.根据权利要求1所述的转换波剩余静校正方法,其特征在于,所述选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量具体包括:
确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串;
依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
6.一种基于引力拟合的转换波剩余静校正系统,其特征在于,包括:
纵波分析层位确定单元,用于根据纵波共检波点检测数据形成的纵波叠加剖面得到第一分析层位;
横波分析层位确定单元,用于根据纵波和转换波速度将所述纵波叠加剖面转换为转换波叠加剖面,并得到转换波叠加剖面上与所述第一分析层位对应的第二分析层位;
剩余静校正量确定单元,用于选取所述第二分析层位预设时间间隔内各分析道的检测数据形成质子串,根据最小势能原理分析所述质子串得到转换波剩余静校正量,根据所述转换波剩余静校正量对检波点的检测数据进行静校正。
7.根据权利要求6所述的转换波剩余静校正系统,其特征在于,所述纵波分析层位确定单元具体用于在根据纵波共检波点检测数据形成的纵波叠加剖面上选取多段反映波形特征和能量特征的同相轴,对多段所述同相轴插值得到所述第一分析层位。
8.根据权利要求6所述的转换波剩余静校正系统,其特征在于,所述横波分析层位确定单元具体用于通过纵波和转换波速度比值将所述纵波叠加剖面转换为转换波叠加剖面,确定所述转换波叠加剖面上与所述第一分析层位对应的第二分析层位。
9.根据权利要求8所述的转换波剩余静校正系统,其特征在于,所述横波分析层位确定单元进一步用于通过所述第二分析层位对所述转换波叠加剖面进行层拉平。
10.根据权利要求6所述的转换波剩余静校正系统,其特征在于,所述剩余静校正量确定单元具体用于确定所述第二分析层位多个位置预设时间间隔内检测数据的振幅能量值得到多个位置的质子串,依次改变每个检测数据在时间方向上的位置形成与其他检测数据的引力形成目标函数,确定所述目标函数取得最大值时各检测数据的位置与预设时间间隔内的初始位置的差为转换波剩余静校正量。
11.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,
所述处理器执行所述程序时实现如权利要求1-5任一项所述方法。
12.一种计算机可读介质,其上存储有计算机程序,其特征在于,
该程序被处理器执行时实现如权利要求1-5任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911265677.7A CN112946745B (zh) | 2019-12-11 | 2019-12-11 | 一种基于引力拟合的转换波剩余静校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911265677.7A CN112946745B (zh) | 2019-12-11 | 2019-12-11 | 一种基于引力拟合的转换波剩余静校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112946745A true CN112946745A (zh) | 2021-06-11 |
CN112946745B CN112946745B (zh) | 2024-06-25 |
Family
ID=76226357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911265677.7A Active CN112946745B (zh) | 2019-12-11 | 2019-12-11 | 一种基于引力拟合的转换波剩余静校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112946745B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8410296D0 (en) * | 1983-04-20 | 1984-05-31 | Chevron Res | Seismic processing |
US4611311A (en) * | 1983-04-20 | 1986-09-09 | Chevron Research Company | Method of seismic processing involving converted P- or S-wave data |
US5508914A (en) * | 1994-09-28 | 1996-04-16 | Western Atlas International, Inc. | Method for calculating static corrections for seismic data |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
RU22830U1 (ru) * | 2001-07-26 | 2002-04-27 | Миколаевский Эрнест Юлианович | Автоматизированная система поиска, разведки и оценки эксплуатационных свойств залежей и месторождений полезных ископаемых и прогноза тектонических и литофизических свойств геологических сред |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN104155694A (zh) * | 2014-08-28 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种反射转换横波共检波点叠加剖面的剩余静校正方法 |
CN104656139A (zh) * | 2013-11-21 | 2015-05-27 | 中国石油天然气集团公司 | 一种转换波的静校正方法及设备 |
CA2879033A1 (en) * | 2014-01-28 | 2015-07-28 | Cgg Services Sa | Methods and apparatus for determining residual static corrections using individual ranges |
CN104914469A (zh) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | 一种转换横波中的静校正方法及装置 |
CN106443791A (zh) * | 2016-11-11 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN108957548A (zh) * | 2018-05-23 | 2018-12-07 | 中国地质调查局油气资源调查中心 | 一种多波多分量联合观测地震页岩气富集区预测技术 |
CN109254322A (zh) * | 2017-07-14 | 2019-01-22 | 中国石油天然气股份有限公司 | 一种基于遗传的剩余静校正模型道建立方法 |
-
2019
- 2019-12-11 CN CN201911265677.7A patent/CN112946745B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611311A (en) * | 1983-04-20 | 1986-09-09 | Chevron Research Company | Method of seismic processing involving converted P- or S-wave data |
GB8410296D0 (en) * | 1983-04-20 | 1984-05-31 | Chevron Res | Seismic processing |
US5508914A (en) * | 1994-09-28 | 1996-04-16 | Western Atlas International, Inc. | Method for calculating static corrections for seismic data |
US6188964B1 (en) * | 1999-09-14 | 2001-02-13 | Ut-Battelle, Llc | Method for using global optimization to the estimation of surface-consistent residual statics |
RU22830U1 (ru) * | 2001-07-26 | 2002-04-27 | Миколаевский Эрнест Юлианович | Автоматизированная система поиска, разведки и оценки эксплуатационных свойств залежей и месторождений полезных ископаемых и прогноза тектонических и литофизических свойств геологических сред |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN104656139A (zh) * | 2013-11-21 | 2015-05-27 | 中国石油天然气集团公司 | 一种转换波的静校正方法及设备 |
CA2879033A1 (en) * | 2014-01-28 | 2015-07-28 | Cgg Services Sa | Methods and apparatus for determining residual static corrections using individual ranges |
CN104914469A (zh) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | 一种转换横波中的静校正方法及装置 |
CN104155694A (zh) * | 2014-08-28 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种反射转换横波共检波点叠加剖面的剩余静校正方法 |
CN106443791A (zh) * | 2016-11-11 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN109254322A (zh) * | 2017-07-14 | 2019-01-22 | 中国石油天然气股份有限公司 | 一种基于遗传的剩余静校正模型道建立方法 |
CN108957548A (zh) * | 2018-05-23 | 2018-12-07 | 中国地质调查局油气资源调查中心 | 一种多波多分量联合观测地震页岩气富集区预测技术 |
Non-Patent Citations (6)
Title |
---|
BYCHKOV,S等: "Near-Surface Correction on Seismic and Gravity Data", JOURNAL OF EARTH SCIENCE, vol. 26, no. 06, pages 851 - 857, XP035563460, DOI: 10.1007/s12583-015-0546-7 * |
J・T・贝拉蒂;王维佳;: "常规纵波、转换波以及水平偏振横波反射法的野外对比研究", 勘探地球物理进展, no. 01, pages 49 - 59 * |
井西利: "大剩余静校正量求解的两步法", 石油地球物理勘探, vol. 38, no. 01, pages 22 - 28 * |
佟恺林;李瑞;唐建明;: "纵波构造控制的共检波点叠加道相关算法研究与应用", 石油地球物理勘探, vol. 51, no. 04, pages 691 - 698 * |
宋玉梅;: "三维转换波静校正技术", 新疆石油地质, vol. 32, no. 06, pages 660 - 663 * |
熊晶璇等: "百分比速度扫描剖面在速度提取中应用", 物探化探计算技术, vol. 37, no. 03, pages 355 - 360 * |
Also Published As
Publication number | Publication date |
---|---|
CN112946745B (zh) | 2024-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107817525B (zh) | 一种近地表常速填充方法及装置 | |
CN107783186B (zh) | 一种近地表速度建模方法及装置 | |
CN109711051B (zh) | 一种考虑滑床岩体结构特征的桩顶位移非线性预测方法 | |
CN117268332B (zh) | 一种建筑物不均匀沉降的监测方法及系统 | |
CN111323814B (zh) | 基于岩石物理模板定量确定砂岩储层含水饱和度的方法 | |
CN112946745B (zh) | 一种基于引力拟合的转换波剩余静校正方法及系统 | |
CN112379436A (zh) | 混合震源激发地震资料静校正量计算方法及装置 | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
JP7438005B2 (ja) | 模擬地震動作成方法 | |
CN112241022B (zh) | 基于射线密度生成层析反演模型速度界面的方法和装置 | |
CN111337973B (zh) | 地震数据重建方法、系统 | |
CN113392372A (zh) | 气井全生命周期无阻流量确定方法及系统 | |
CN112415598A (zh) | 近地表地层品质因子模型的建立方法和装置 | |
CN112946736A (zh) | 一种三维观测系统重建方法及系统 | |
CN113900141A (zh) | 油气分布预测方法及装置 | |
CN112578458A (zh) | 叠前弹性阻抗随机反演方法、装置、存储介质及处理器 | |
CN109100794B (zh) | 一种时窗加权的相干速度反演方法及系统 | |
CN111596346B (zh) | 弹性波速度反演方法和装置 | |
CN113204855A (zh) | 地层缩短速率确定方法、系统、设备及可读介质 | |
CN113031071B (zh) | 地震波长波长静校正方法及装置 | |
CN112180441B (zh) | 转换波初始速度建模方法及装置 | |
CN108254787B (zh) | 波的走时获得方法及装置、成像方法及装置 | |
CN111721292B (zh) | 一种利用地磁的室内定位方法 | |
Guo et al. | A wavenumber-domain iterative approach for apparent density mapping of an undulant layer and its application in central South China |
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 |