CN105242311B - 预估初至拾取的时窗的方法和装置 - Google Patents
预估初至拾取的时窗的方法和装置 Download PDFInfo
- Publication number
- CN105242311B CN105242311B CN201510689844.6A CN201510689844A CN105242311B CN 105242311 B CN105242311 B CN 105242311B CN 201510689844 A CN201510689844 A CN 201510689844A CN 105242311 B CN105242311 B CN 105242311B
- Authority
- CN
- China
- Prior art keywords
- shot
- point
- time
- offset
- sequence
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001514 detection method Methods 0.000 claims description 63
- 238000012163 sequencing technique Methods 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 8
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供的一种预估初至拾取的时窗的方法和装置。一个实施方式中,所述方法包括:获取预定地区内的至少两个炮点的位置信息及至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;将至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;根据至少两个炮点中每个炮点的第一对应关系以及炮检距序列对应获取至少两个炮点中每个炮点的时间序列;根据待测炮点的位置信息以及对至少两个炮点的位置信息和时间序列进行第一内插运算,得到待测炮点的时间序列;按照第二预定规则对待测炮点的时间序列进行运算,从而得到待测炮点的时窗。本发明实现了提供能提高初至拾取的准确率的一种预估初至拾取的时窗的方法的目的。
Description
技术领域
本发明涉及石油天然气开采领域,尤其涉及一种预估初至拾取的时窗的方法和装置。
背景技术
地震勘探的原理是:在地面上某点放置激发源后,例如:打井放炮,使激发源爆炸产生的地震波向下传播。地震波遇到地层(速度与密度的乘积有差异)的分界面时,通常会发生反射;同时另一部分地震波还会继续向下传播,碰到相似的地层界面后还会产生反射和透射,即一部分地震波的能量反射回地面,另一部分继续向下传播。与此同时,地面上的检波器把来自各个地层分界面的反射波引起地面震动的情况记录下来,然后根据地震波从地面开始向下传播的时刻和地层分界面反射波到达地面的时刻,得出地震波从地面向下传播到达地层分界面,又反射回地面的总时间,再用别的方法测定出地震波在岩层中传播的速度,最后就可以得到地层分界面的埋藏深度了。
初至波是激发源信号经过地层之后最先返回地面的信号,初至波简称初至。初至时间是指初至波到达检波器的时间。在有效信号返回地面之前检波器所接受的信号主要是一些环境噪音。初至拾取就是从检波器所接受到的信号中拾取初至波。该初至拾取的准确性是地震勘探的一项基础而又重要的工作,同时也是地震勘探的重要难题。现有的初至拾取方法主要有手动拾取和自动拾取。手动拾取较自动拾取工作量大,效率低。目前初至拾取的自动拾取方法主要有:相关法、能量比值法、分形维数法、神经网络算法、基于边缘检测算法等等。为了降低拾取难度,一般拾取方法是需要定义一个一定宽度的时窗,使初至时间包含在时窗内,拾取时仅需要考虑时窗范围内的时间,从而使得初至时间得到限制,在一定程度上提高了拾取的准确度。
现有技术中定义时窗的方法是:在地震记录资料中选取一个炮点,称为第一炮点,获取第一炮点的第一大地坐标;选取该第一炮点的至少两个地震道作为第一检测点,获取该第一检测点的炮检距和初至时间,对该至少两个第一检测点的炮检距和初至时间进行线性内插运算,从而得到第一炮点的每一个炮检距所对应的初至时间。将每一个炮检距所对应的初至时间增加一个预定量的值得到每一个炮检距所对应的最大初至时间,将每一个炮检距所对应的初至时间减小一个预定量的值得到每一个炮检距所对应的最小初至时间,如图1所示,该图1的横坐标表示炮检距,纵坐标表示时间。该第一炮点的每一个地震道的时窗50是最小初至时间与最大初至时间之间的范围。所以该第一炮的所有地震道可以根据各自的炮检距从该时窗上得到每一地震道的时窗。对于其他炮点来说,如图2所示,根据各自的炮检距从该第一炮的时窗上也可以得到每一地震道的时窗51。
由于不同位置的炮点,其地表和地下地质条件不同,所以地震波在地下的传播速度不同,所以地震波到达具有相同炮检距的检波器的初至时间不同,所以现有技术的时窗的定义中,由于将第一炮中定义的时窗直接应用于另一炮中,容易造成其他炮点的初至时间不包括在第一炮的时窗内,从而不能在第一炮的时窗内拾取到初至波,从而降低了初至拾取的准确率,从而影响了地震勘探工作的顺利进行。
发明内容
本发明的目的是提供了一种预估初至拾取的时窗的方法及装置,能够提高初至拾取的准确率。
为实现上述目的,本申请提供的技术方案如下。
本申请一方面提供了一种预估初至拾取的时窗的方法,其用于预估位于预定地区内的待测炮点的时窗,其包括:获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;将所述至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;根据所述至少两个炮点中每个炮点的第一对应关系以及所述炮检距序列对应获取所述至少两个炮点中每个炮点的时间序列;根据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的时间序列;按照第二预定规则对所述待测炮点的时间序列进行运算,从而得到所述待测炮点的时窗。
优选地,步骤对所述每个炮点的位置时间序列进行第一内插运算,从而得出该预定地区内的待测炮点的位置信息所对应的初至时间序列包括:获取所述每个炮点的位置时间序列的第K项,构成第K个散点序列;对所述第K个散点序列进行第一内插运算,得到第K个曲面;根据所述待测炮点的位置信息和第K个曲面对应获取第K个初至时间。
优选地,所述第一内插包括二次曲线内插法和/或多次曲线内插法。
优选地,在步骤按照第二预定规则对所述待测炮点的初至时间序列进行运算,从而得到所述待测炮点的时窗之后,还包括:根据所述待测炮点的时窗以及所述炮检距序列对应获取所述待测炮点的炮检距时窗序列。
优选地,所述第一预定的规则为将所述至少两个炮点中所有的炮检距按照从小到大的顺序进行排序,如果所述至少两个炮点所有的炮检距中存在两个相同的炮检距,那么则删除该两个相同的炮检距中的一个,从而得到项数为K的炮检距序列。
优选地,步骤获取所述预定地区内的至少两个炮点的位置信息及每个炮点的每一个炮检距和与其对应的初至时间的第一对应关系,包括:获取所述预定地区内的至少两个炮点的位置信息;选取所述每个炮点的至少两个地震道作为所述每个炮点的检测点,获取所述每个炮点的每个检测点的炮检距和其对应的初至时间;根据所述每个炮点的每一个炮检距、所述每个炮点的每个检测点的炮检距和其对应的初至时间,通过第二内插运算,获取所述每个炮点的每一个炮检距所对应的初至时间;根据所述每个炮点的每一个炮检距和与其对应的初至时间获取所述第一对应关系。
优选地,所述第二内插为线性内插法。
本申请另一方面提供了一种预估初至拾取的时窗的装置,其包括:第一获取模块,用于获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;排序模块,用于将所述至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;第二获取模块,用于根据所述至少两个炮点中每个炮点的第一对应关系以及所述炮检距序列对应获取所述至少两个炮点中每个炮点的时间序列;第一运算模块,用于根据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的时间序列;第二运算模块,用于按照第二预定规则对所述待测炮点的时间序列进行运算,从而得到所述待测炮点的时窗。
本发明提供的一种预估初至拾取的时窗的方法的有益效果是:本发明通过对预定地区内的多个炮点的位置时间序列进行第一内插运算,从而得出该预定地区内的待测炮点的位置信息所对应的初至时间序列,然后按照第二预定规则对待测炮点的初至时间序列进行运算,从而得到待测炮点的时窗。也就是说,通过对预定地区内的多个炮点的时窗进行第一内插运算,从而得到待测炮点的时窗,从而使得预定地区的所有炮点的时窗都能内插到合适的位置,从而区别于现有技术中的将第一炮中定义的时窗直接应用于另一炮中,从而容易造成其他炮点的初至时间不包括在第一炮的时窗内,从而不能在第一炮的时窗内拾取到初至波的情况,从而实现了本发明能提供提高初至拾取的准确率的一种预估初至拾取的时窗的方法的目的。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术的示意图;
图2为现有技术的示意图;
图3为本发明实施方式的预估初至拾取的时窗的方法的一种流程图;
图4为本发明实施方式的预估初至拾取的时窗的方法的一种流程图;
图5为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图6为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图7为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图8为本发明实施方式的预估初至拾取的时窗的方法的一种流程图;
图9为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图10为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图11为本发明实施方式的预估初至拾取的时窗的方法的一种示意图;
图12为本发明实施方式的预估初至拾取的时窗的装置的一种流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,本发明提供一种预估初至拾取的时窗的方法,其用于预估位于预定地区内的待测炮点的时窗,其包括:S14:获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;S16:将所述至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;S18:根据所述至少两个炮点中每个炮点的第一对应关系以及所述炮检距序列对应获取所述至少两个炮点中每个炮点的时间序列;S20:根据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的时间序列;S22:按照第二预定规则对所述待测炮点的时间序列进行运算,从而得到所述待测炮点的时窗。
本发明通过对预定地区内的多个炮点的位置时间序列进行第一内插运算,从而得出该预定地区内的待测炮点的位置信息所对应的初至时间序列,然后按照第二预定规则对待测炮点的初至时间序列进行运算,从而得到待测炮点的时窗。也就是说,通过对预定地区内的多个炮点的时窗进行第一内插运算,从而得到待测炮点的时窗,从而使得预定地区的所有炮点的时窗都能内插到合适的位置,从而区别于现有技术中的将第一炮中定义的时窗直接应用于另一炮中,从而容易造成其他炮点的初至时间不包括在第一炮的时窗内,从而不能在第一炮的时窗内拾取到初至波的情况,从而实现了本发明能提供提高初至拾取的准确率的一种预估初至拾取的时窗的方法的目的。
如图3所示,在本实施方式中,以西部某地区作为预定地区,预估位于该预定地区内的待测炮点的时窗。
如图3所示,在本实施方式中,步骤S14:获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系。
如图4所示,在本实施方式中,该步骤S14包括:S141:获取预定地区内的至少两个炮点的位置信息,从而可以通过运算该预定地区内的至少两个炮点的时窗得到该预定地区内的待测炮点的时窗。
在本实施方式中,在该预定地区的地震记录资料上选取两个炮点。该两个炮点为第一炮点和第二炮点。该第一炮点和第二炮点是按照选取的时间先后顺序所定义的。获取该第一炮点和第二炮点的位置信息。该第一炮点和第二炮点为该预定地区的两个炮点,所以通过运算该第一炮点和第二炮点的时窗从而得到待测炮点的时窗,即为根据该预定地区内两个具有不同位置信息的炮点的时窗运算另一个位于该预定地区内的炮点的时窗,从而区别于现有技术中的将第一炮中定义的时窗直接应用于另一炮中的情况。在本实施方式中,该位置信息为大地坐标。大地坐标是大地测量中以参考椭球面为基准面的坐标。该位置信息记为(Ei,Ni),其中i为炮点的序号,i=1,2,3,…m。在本实施方式中,i=1,2。
在本实施方式中,在预定地区的地震记录资料中,选择1031号炮点为第一炮点,该1031号炮点的位置信息为(5906.0,241532.0)。选择154号炮点为第二炮点,该154号炮点的位置信息为(5955.9,241431.8)。
如图4所示,在本实施方式中,该步骤S14还包括:S143:选取至少两个炮点的至少两个地震道作为至少两个炮点的检测点,获取每个检测点的炮检距和其对应的初至时间。
在本实施方式中,在地震记录资料上分别选取第一炮点和第二炮点的检测点,该第一炮点的检测点称为第一检测点,该第一检测点包括至少两个地震道。该第二炮点的检测点称为第二检测点,该第二检测点包括至少两个地震道,该第一检测点和第二检测点分别用于运算该第一炮点和第二炮点的时窗。该第一检测点和第二检测点的选取需要根据该第一检测点和第二检测点是否能够从地震记录资料中获取得到其初至时间所确定,从而能够获取第一检测点和第二检测点的炮检距和其对应的初至时间。该第一检测点和第二检测点的个数至少两个,从而通过该已知初至时间的第一检测点和第二检测点运算该第一炮点和第二炮点的其他地震道的初至时间将会因为选取的检测点的个数越多而运算的初至时间越准确。在本实施方式中,选取该第一炮点的3个地震道作为第一检测点。选取该第二炮点的3个地震道作为第二检测点。
根据至少两个炮点的每一个炮检距、至少两个炮点的每个检测点的炮检距和每个检测点的初至时间,通过第二内插运算,获取至少两个炮点的每一个炮检距所对应的初至时间;
如图4所示,在本实施方式中,该步骤S14还包括:S145:根据至少两个炮点的每一个炮检距、至少两个炮点的每个检测点的炮检距和每个检测点的初至时间,通过第二内插运算,获取至少两个炮点的每一个炮检距所对应的初至时间。
如图4所示,在本实施方式中,该步骤S14还包括:S147:根据至少两个炮点的每一个炮检距和至少两个炮点的每一个炮检距所对应的初至时间获取第一对应关系。
在本实施方式中,该每个炮点的每个检测点的炮检距和其对应的初至时间可以从地震记录资料中获取。该每个炮点的每一个炮检距可以通过现场测量的方法获取。所以根据每个炮点的每一个炮检距、每个炮点的每个检测点的炮检距和其对应的初至时间,通过第二内插运算,能够获取每个炮点的每一个炮检距所对应的初至时间,从而能够根据每个炮点的每一个炮检距和与其对应的初至时间获取第一对应关系。
如图5所示,在本实施方式中,该第一对应关系可以在以炮检距为横坐标,时间为纵坐标的坐标轴上表示出来。该第一对应关系还可以表示成如下的序列:
炮检距:xi1,xi2,…xin;
时间:ti1,ti2,…tin;
其中,i为炮点,i=1,2,3,…m。在本实施方式中,i=1,2。n为检测点的个数,n>=2。在本实施方式中,n=3。xin为第i个炮点的第n个检测点;tin为第i个炮点的第n个检测点的初至时间。上述序列表示:当第i个炮点的第1个检测点的炮检距为xi1时,该第i个炮点的第1个检测点所对应的初至时间为ti1;当第i个炮点的第2个检测点的炮检距为xi2时,该第i个炮点的第2个检测点所对应的初至时间为ti2;…当第i个炮点的第n个检测点的炮检距为xin时,该第i个炮点的第n个检测点所对应的初至时间为tin。
在本实施方式中,该第一炮点的第一对应关系可以表示成如下的序列:
炮检距:404,3828,…7377;
时间:270,1642,…2702;
由上述序列可以知道:当该第一炮点的第1个检测点的炮检距为404时,该炮检距为404的检测点所对应的初至时间为270;当该第一炮点的第2个检测点的炮检距为3828时,该炮检距为3828的检测点所对应的初至时间为1642;…当该第一炮点的第n个检测点的炮检距为7377时,该炮检距为7377的检测点所对应的初至时间为2702,从而可以根据该序列,通过第二内插运算,获取该第一炮点的每一个炮检距所对应的初至时间。
该第一炮点的其他地震道的炮检距为0,20,40,60,80,100,…10000。所以可以通过第一炮点的每个检测点的炮检距和其对应的初至时间,进行第二内插运算,得到该第一炮点的其他地震道的炮检距所对应的初至时间。该第二内插为线性内插法、非线性内插法,优选为线性内插法。
如图6所示,根据第一炮点的每一个炮检距和与其对应的初至时间获取的第一对应关系为第一折线53。
在本实施方式中,该第二炮点的第一对应关系可以表示成如下的序列:
炮检距:428,2526,…6825;
时间:222,944,…2070;
由上述序列可以知道:当该第二炮点的第1个检测点的炮检距为428时,该炮检距为428的检测点所对应的初至时间为222;当该第二炮点的第2个检测点的炮检距为2526时,该炮检距为2526的检测点所对应的初至时间为944;…当该第二炮点的第n个检测点的炮检距为6825时,该炮检距为6825的检测点所对应的初至时间为2070,从而可以根据该序列,通过第二内插运算,获取该第二炮点的每一个炮检距所对应的初至时间。
该第二炮点的其他地震道的炮检距为0,20,40,60,80,100,…10000。所以可以通过第二炮点的每个检测点的炮检距和其对应的初至时间,进行第二内插运算,得到该第二炮点的其他地震道的炮检距所对应的初至时间。该第二内插为线性内插法、非线性内插法,优选为线性内插法。
如图7所示,根据第二炮点的每一个炮检距和与其对应的初至时间获取的第一对应关系为第二折线55。
如图3所示,在本实施方式中,步骤S16:将至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;第一预定的规则为将至少两个炮点中所有的炮检距按照从小到大的顺序进行排序,如果至少两个炮点所有的炮检距中存在两个相同的炮检距,那么则删除该两个相同的炮检距中的一个,从而得到项数为K的炮检距序列。该项数为K的炮检距序列可以表示成如下的形式:
炮检距序列为:x1,x2,x3,…xk;
其中,xk为炮检距,k为炮检距个数,k=n1+n2+…+nm。x1<x2<x3<x4<…<xk。
如图3所示,在本实施方式中,步骤S18:根据至少两个炮点中每个炮点的第一对应关系以及炮检距序列对应获取至少两个炮点中每个炮点的时间序列。
在本实施方式中,按照炮检距序列,对每个炮点的第一对应关系中进行第三内插运算,从而得到每个炮点的时间序列:ti1,ti2,…tik;其中,i为炮点,i=1,2,3,…m;k为炮检距个数,k=n1+n2+…+nm。tik为第i炮点的第k个炮检距所对应的时间。该第三内插为线性内插法、非线性内插法,优选为线性内插法。根据图5所示,在炮检距:xi1,xi2,…xin中找到x2,然后根据第一内插运算得到x2所对应的t2,从而得到炮检距为x2所对应的时间t2。
在本实施方式中,根据第一炮点的第一对应关系以及炮检距序列对应获取第一炮点的时间序列。该第一炮点的第一对应关系可以表示成如下的序列:
炮检距:404,3828,…7377;
时间:270,1642,…2702;
从而根据炮检距序列:x1,x2,x3,…xk;进行第三内插运算获取第一炮点的时间序列:t11,t12,…t1k;k为炮检距个数,k=n1+n2+…+nm。t1k为第一炮点的第k个炮检距所对应的时间。
在本实施方式中,根据第二炮点的第一对应关系以及炮检距序列对应获取第一炮点的时间序列。该第二炮点的第一对应关系可以表示成如下的序列:
炮检距:428,2526,…6825;
时间:222,944,…2070;
从而根据炮检距序列:x1,x2,x3,…xk;进行第三内插运算获取第二炮点的时间序列:t21,t22,…t2k;k为炮检距个数,k=n1+n2+…+nm。T2k为第二炮点的第k个炮检距所对应的时间。
如图3所示,在本实施方式中,步骤S20:根据待测炮点的位置信息以及对至少两个炮点的位置信息和时间序列进行第一内插运算,得到待测炮点的时间序列。该第一内插包括二次曲线内插法和/或多次曲线内插法。
如图8所示,在本实施方式中,步骤S20包括:S201:根据至少两个炮点中每个炮点的位置信息以及至少两个炮点中每个炮点的时间序列对应获取每个炮点的位置时间序列。
在本实施方式中,每个炮点的位置信息为(Ei,Ni),其中i为炮点的序号,i=1,2,3,…m。在本实施方式中,i=1,2。每个炮点的初至时间序列为ti1,ti2,…tik;其中,i为炮点,i=1,2,3,…m;k为炮检距个数,k=n1+n2+…+nm。tik为第i炮点的第k个炮检距所对应的初至时间。所以每个炮点的位置时间序列为:
(Ei,Ni,ti1),(Ei,Ni,ti2),(Ei,Ni,ti3),…(Ei,Ni,tik)
其中i为炮点的序号,i=1,2,3,…m。在本实施方式中,i=1,2。k为炮检距个数,k=n1+n2+…+nm。
在本实施方式中,第一炮点的位置时间序列为:
(5906.0,241532.0,270),(5906.0,241532.0,1642),…(5906.0,241532.0,2702)
在本实施方式中,第二炮点的位置时间序列为:
(5955.9,241431.8,222),(5906.0,241532.0,944),…(5906.0,241532.0,2070)。
如图8所示,在本实施方式中,步骤S20还包括:S203:获取每个炮点的位置时间序列的第K项,构成第K个散点序列,从而形成空间上的若干个散点,从而k个时间形成k个散点序列:
散点1:(E1,N1,t11)、(E2,N2,t21)、(E3,N3,t31)、…(Em,Nm,tm1)、
散点2:(E1,N1,t12)、(E2,N2,t22)、(E3,N3,t32)、…(Em,Nm,tm2)、
散点3:(E1,N1,t13)、(E2,N2,t23)、(E3,N3,t33)、…(Em,Nm,tm3)、
…
散点k:(E1,N1,t1k)、(E2,N2,t2k)、(E3,N3,t3k)、…(Em,Nm,tmk)、
如图8所示,在本实施方式中,步骤S20包括:S205:对第K个散点序列进行第一内插运算,得到第K个曲面:对散点1序列进行第一内插运算,该散点1序列:(E1,N1,t11)、(E2,N2,t21)、(E3,N3,t31)、…(Em,Nm,tm1)能形成第1个曲面;对散点2序列进行第一内插运算,该散点2序列:(E1,N1,t12)、(E2,N2,t22)、(E3,N3,t32)、…(Em,Nm,tm2)能形成第2个曲面;对散点3序列进行第一内插运算,该散点3序列:(E1,N1,t13)、(E2,N2,t23)、(E3,N3,t33)、…(Em,Nm,tm3)能形成第3个曲面;…对散点K序列进行第一内插运算,该散点K序列:(E1,N1,t1k)、(E2,N2,t2k)、(E3,N3,t3k)、…(Em,Nm,tmk)能形成第K个曲面。
如图8所示,在本实施方式中,步骤S20包括:S207:根据待测炮点的位置信息和第K个曲面对应获取待测炮点的第K个时间。对于该预定地区内的任一待测炮点,其位置信息为(EK,NK),其中K为炮点的序号,根据该位置信息(EK,NK)可以对第1曲面进行第一内插运算,从而可以从第1个曲面获取第1个时间t1;根据该位置信息(EK,NK)可以对第2曲面进行第一内插运算,从而可以从第2个曲面获取第2个时间t2;根据该位置信息(EK,NK)可以对第3曲面进行第一内插运算,从而可以从第3个曲面获取第3个时间t3;…根据该位置信息(EK,NK)可以对第k个曲面进行第一内插运算,从而可以从第k个曲面获取第k个时间tk,从而得出该预定地区内的待测炮点的位置信息所对应的时间序列:t1,t2,…tk;k为曲面的个数。
在本实施方式中,根据第一炮点的位置时间序列:(5906.0,241532.0,270),(5906.0,241532.0,1642),…(5906.0,241532.0,2702)和第二炮点的位置时间序列:(5955.9,241431.8,222),(5906.0,241532.0,944),…(5906.0,241532.0,2070)。进行第一内插运算,从而得出该预定地区内的待测炮点的位置信息所对应的时间序列:t1,t2。
步骤S22:按照第二预定规则对待测炮点的时间序列进行运算,从而得到待测炮点的时窗。该第二预定规则为在时间序列上增加预定值,得到最大初至时间序列,在时间序列上减少预定值,得到最小初至时间序列。该待测炮点的时窗为最小初至时间与最大初至时间之间的范围。该预定值例如为200ms。该预定值的两倍为时窗的宽度。从而通过定义待测炮点的时窗,使得待测炮点的初至时间包含在时窗内,拾取时仅需要考虑时窗范围内的时间,从而使得初至时间得到限制,在一定程度上提高了拾取的准确度。
如图9所示,按照第二预定规则对第一炮点的时间序列进行运算,从而得到第一炮点的时窗57。
如图10所示,按照第二预定规则对第二炮点的时间序列进行运算,从而得到第二炮点的时窗59。
如图11所示,待测炮点的时窗61为:(t1-Δt~t1+Δt),(t2-Δt~t2+Δt),…(tk-Δt~tk+Δt);
本发明通过对预定地区内的多个炮点的位置时间序列进行第一内插运算,从而得出该预定地区内的待测炮点的位置信息所对应的初至时间序列,也就是说,可以通过对预定地区内的多个炮点的时窗进行第一内插运算,从而得到待测炮点的时窗,从而使得预定地区的所有炮点的时窗都能内插到合适的位置,从而区别于现有技术中的将第一炮中定义的时窗直接应用于另一炮中,从而容易造成其他炮点的初至时间不包括在第一炮的时窗内,从而不能在第一炮的时窗内拾取到初至波的情况,从而实现了本发明能提供提高初至拾取的准确率的一种预估初至拾取的时窗的方法的目的。
步骤S22:根据待测炮点的时窗以及炮检距序列对应获取待测炮点的炮检距时窗序列。根据炮检距序列:x1,x2,x3,…xk和待测炮点的时窗(t1-Δt~t1+Δt),(t2-Δt~t2+Δt),…(tk-Δt~tk+Δt);该炮检距时窗序列可以表示为:
〔x1,(t1-Δt~t1+Δt)〕,〔x2,(t2-Δt~t2+Δt)〕,〔x3,(t3-Δt~t3+Δt)〕,…〔xk,(tk-Δt~tk+Δt)〕。从而根据该炮检距时窗序列,可以得到该待测炮点的不同的炮检距所对应的时窗。
基于同一构思,本申请实施方式还提供了一种预估初至拾取的时窗的装置,如下面的实施方式所述。由于预估初至拾取的时窗的装置解决问题的原理与预估初至拾取的时窗的方法相似,因此预估初至拾取的时窗的装置的实施可以参见上述预估初至拾取的时窗的方法的实施,重复之处不再赘述。以下所使用的,术语“模块”或者“单元”可以实现预定功能的软件和/或硬件的组合。尽管以下实施方式所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
如图12所示,预估初至拾取的时窗的装置可以包括:第一获取模块101,排序模块103,第二获取模块105,第一运算模块107,第二运算模块109。
预估初至拾取的时窗的装置,其包括:
第一获取模块101,可以用于获取预定地区内的至少两个炮点的位置信息及至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;
排序模块103,可以用于将至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;
第二获取模块105,可以用于根据至少两个炮点中每个炮点的第一对应关系以及炮检距序列对应获取至少两个炮点中每个炮点的时间序列;
第一运算模块107,可以用于根据待测炮点的位置信息以及对至少两个炮点的位置信息和时间序列进行第一内插运算,得到待测炮点的时间序列;
第二运算模块109,可以用于按照第二预定规则对待测炮点的时间序列进行运算,从而得到待测炮点的时窗。
以上仅为本发明的几个实施例,本领域的技术人员依据申请文件公开的内容可以对本发明实施例进行各种改动或变型而不脱离本发明的精神和范围。
Claims (7)
1.一种预估初至拾取的时窗的方法,其用于预估位于预定地区内的待测炮点的时窗,其特征在于,其包括:
获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;
将所述至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;
根据所述至少两个炮点中每个炮点的第一对应关系以及所述炮检距序列对应获取所述至少两个炮点中每个炮点的时间序列;
根据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的时间序列;
按照第二预定规则对所述待测炮点的时间序列进行运算,从而得到所述待测炮点的时窗其中,
所述据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的初至时间序列包括:
根据所述至少两个炮点中每个炮点的位置信息以及所述至少两个炮点中每个炮点的时间序列对应获取每个炮点的位置时间序列;
获取所述每个炮点的位置时间序列的第K项,构成第K个散点序列;
对所述第K个散点序列进行第一内插运算,得到第K个曲面;
根据所述待测炮点的位置信息和第K个曲面对应获取所述待测炮点的第K个时间。
2.根据权利要求1所述的预估初至拾取的时窗的方法,其特征在于:所述第一内插包括二次曲线内插法和/或多次曲线内插法。
3.根据权利要求2所述的预估初至拾取的时窗的方法,其特征在于:在步骤按照第二预定规则对所述待测炮点的初至时间序列进行运算,从而得到所述待测炮点的时窗之后,还包括:根据所述待测炮点的时窗以及所述炮检距序列对应获取所述待测炮点的炮检距时窗序列。
4.根据权利要求3所述的预估初至拾取的时窗的方法,其特征在于:所述第一预定的规则为将所述至少两个炮点中所有的炮检距按照从小到大的顺序进行排序,如果所述至少两个炮点所有的炮检距中存在两个相同的炮检距,那么则删除该两个相同的炮检距中的一个,从而得到项数为K的炮检距序列。
5.根据权利要求4所述的预估初至拾取的时窗的方法,其特征在于:步骤获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系,包括:
获取所述预定地区内的至少两个炮点的位置信息;
选取所述至少两个炮点的至少两个地震道作为所述至少两个炮点的检测点,获取所述每个检测点的炮检距和其对应的初至时间;
根据所述至少两个炮点的每一个炮检距、所述至少两个炮点的每个检测点的炮检距和所述每个检测点的初至时间,通过第二内插运算,获取所述至少两个炮点的每一个炮检距所对应的初至时间;
根据所述至少两个炮点的每一个炮检距和所述至少两个炮点的每一个炮检距所对应的初至时间获取所述第一对应关系。
6.根据权利要求5所述的预估初至拾取的时窗的方法,其特征在于:所述第二内插为线性内插法。
7.一种预估初至拾取的时窗的装置,其用于预估位于预定地区内的待测炮点的时窗,其特征在于,其包括:
第一获取模块,用于获取所述预定地区内的至少两个炮点的位置信息及所述至少两个炮点中每个炮点的炮检距与初至时间之间的第一对应关系;
排序模块,用于将所述至少两个炮点中所有的炮检距按照第一预定规则进行排序,得到炮检距序列;
第二获取模块,用于根据所述至少两个炮点中每个炮点的第一对应关系以及所述炮检距序列对应获取所述至少两个炮点中每个炮点的时间序列;
第一运算模块,用于根据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的时间序列;
第二运算模块,用于按照第二预定规则对所述待测炮点的时间序列进行运算,从而得到所述待测炮点的时窗;其中,
所述据所述待测炮点的位置信息以及对所述至少两个炮点的位置信息和所述时间序列进行第一内插运算,得到所述待测炮点的初至时间序列包括:
根据所述至少两个炮点中每个炮点的位置信息以及所述至少两个炮点中每个炮点的时间序列对应获取每个炮点的位置时间序列;
获取所述每个炮点的位置时间序列的第K项,构成第K个散点序列;
对所述第K个散点序列进行第一内插运算,得到第K个曲面;
根据所述待测炮点的位置信息和第K个曲面对应获取所述待测炮点的第K个时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510689844.6A CN105242311B (zh) | 2015-10-22 | 2015-10-22 | 预估初至拾取的时窗的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510689844.6A CN105242311B (zh) | 2015-10-22 | 2015-10-22 | 预估初至拾取的时窗的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242311A CN105242311A (zh) | 2016-01-13 |
CN105242311B true CN105242311B (zh) | 2018-01-02 |
Family
ID=55040011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510689844.6A Active CN105242311B (zh) | 2015-10-22 | 2015-10-22 | 预估初至拾取的时窗的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242311B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106610889B (zh) * | 2016-12-28 | 2019-08-09 | 北京奇艺世纪科技有限公司 | 一种数据业务的预估方法及装置 |
CN107576983B (zh) * | 2017-11-02 | 2019-06-04 | 中国石油集团东方地球物理勘探有限责任公司 | 提高单炮地震记录初至拾取准确度的方法和资料评价方法 |
CN112241021B (zh) * | 2019-07-16 | 2023-05-26 | 中国石油天然气集团有限公司 | 初至拾取方法及装置 |
CN112083497A (zh) * | 2020-08-27 | 2020-12-15 | 中国石油天然气集团有限公司 | 一种模型约束下的初至拾取方法及装置 |
CN113109869A (zh) * | 2021-03-30 | 2021-07-13 | 成都理工大学 | 一种页岩超声测试波形初至的自动拾取方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609160B (zh) * | 2008-06-19 | 2011-08-03 | 中国石油天然气集团公司 | 一种抗背景噪音强的地震数据初至波自动拾取方法 |
CN101630015B (zh) * | 2008-07-16 | 2012-01-11 | 中国石油天然气集团公司 | 一种提高初至波拾取精度和效率的方法 |
CN102081167B (zh) * | 2009-12-01 | 2012-08-15 | 中国石油天然气集团公司 | 一种三维vsp数据初至波拾取方法 |
US9103935B2 (en) * | 2013-02-04 | 2015-08-11 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
CN103592698B (zh) * | 2013-10-31 | 2016-03-09 | 中国石油天然气集团公司 | 一种基于地震属性的观测系统评价方法及装置 |
CN104133239A (zh) * | 2014-07-29 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种地震单炮记录目的层时窗的定义方法 |
-
2015
- 2015-10-22 CN CN201510689844.6A patent/CN105242311B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105242311A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105242311B (zh) | 预估初至拾取的时窗的方法和装置 | |
US11209561B2 (en) | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function | |
RU2621482C2 (ru) | Системы и способы для опережающего измерения удельного сопротивления с использованием информации опорной скважины | |
CA2823710C (en) | Methods and systems regarding models of underground formations | |
US20190227185A1 (en) | System and method for predicting fault seal from seismic data | |
EP3394647B1 (en) | A methodology for building realistic numerical forward stratigraphic models in data sparse environment | |
NO334218B1 (no) | Behandling av målinger på lydbølgeformer for å bestemme langsomheten | |
CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
US10024991B2 (en) | Method and device for processing seismic signals | |
WO2017165341A2 (en) | Method and device for estimating sonic slowness in a subterranean formation | |
CN110703354B (zh) | 一种基于井震结合的岩性边界识别方法 | |
EP3631524B1 (en) | System and method for predicting fault seal from seismic data | |
NO20121473A1 (no) | System for modellering av geologiske strukturer | |
CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
Dickmann | 3D Tunnel Seismic Prediction: A next generation tool to characterize rock mass condition ahead of the tunnel face | |
AU2019312924B2 (en) | System and method for seismic amplitude analysis | |
NO20180869A1 (en) | Seismic well ties using blocking schemes | |
CN112649848B (zh) | 利用波动方程求解地震波阻抗的方法和装置 | |
NO20190217A1 (en) | Correcting biases in microseismic-event data | |
CN111312272A (zh) | 一种用于减少近井声数据集中噪声信号的产品、方法及系统 | |
CN111077573A (zh) | 一种确定地层弹性参数的方法、装置及系统 | |
CN111965706B (zh) | 地震反演方法及装置 | |
US20240168189A1 (en) | Systems, devices, and methods for generating average velocity maps of subsurface formations | |
US20240311444A1 (en) | Sequence stratigraphic interpretation of seismic data | |
CN115951408A (zh) | 地下地层的声波探测的慢度提取方法、装置、介质及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |