CN109212598A - 基于直达波反演的三维空间二次定位方法 - Google Patents
基于直达波反演的三维空间二次定位方法 Download PDFInfo
- Publication number
- CN109212598A CN109212598A CN201811292082.6A CN201811292082A CN109212598A CN 109212598 A CN109212598 A CN 109212598A CN 201811292082 A CN201811292082 A CN 201811292082A CN 109212598 A CN109212598 A CN 109212598A
- Authority
- CN
- China
- Prior art keywords
- direct wave
- solution
- localization method
- dimensional space
- coordinate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000004807 localization Effects 0.000 title claims abstract description 24
- 239000000523 sample Substances 0.000 claims abstract description 8
- 239000013535 sea water Substances 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 abstract description 6
- 238000007689 inspection Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000010415 tropism Effects 0.000 description 1
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
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
本发明提供一种基于直达波反演的三维空间二次定位方法,该基于直达波反演的三维空间二次定位方法包括:步骤1,加载观测***,并进行直达波的拾取;步骤2,求取实际炮检距;步骤3,列出关于检波点的坐标方程,并进行坐标方程的求解;步骤4,去除坐标方程的无效解;步骤5,求取检波点二次定位坐标;步骤6,进行坐标替换。该基于直达波反演的三维空间二次定位方法相对常规的二维方法结果更合理,精确更高,效果可靠,流程及参数设置简单,运算速度快,操作简单易实现。
Description
技术领域
本发明涉及油气勘探地震资料处理技术领域,特别是涉及到一种基于直达波反演的三维空间二次定位方法。
背景技术
现有的检波点二次定位方法主要有声波定位方法和初至波定位方法,声波二次定位是一种只能由记录的炮点位置推算检波点位置的定位***,在实际应用中,声波二次定位存在一定的局限性,声波二次定位方法所需设备昂贵,数据采集慢且回收率低,影响生产效率,无法在实际生产中大量推广使用,并且无法对推算出的检波点位置的准确性进行验证。目前使用较多的是初至波二次定位方法,该方法利用圆圆相交的原理,但是该原理是基于平面二维空间的算法,实际生产中,炮点、检波点并不能正好在一个平面,这样就会产生一定的误差,使得计算出来的检波点坐标的准确性降低,并且,该方法无法求得检波点的水深。
在浅海和海陆过渡带,检波器是通过抛掷或机械放缆方式沉入海底,在受到海流、潮汐、涌浪的强大冲击时,检波器很难完全固定在海底,会出现漂移现象,而海底电缆又无法实现实时定位,造成检波点实际位置与设计的坐标不一致。如果对存在偏差的检波点坐标进行后续地震资料处理,会影响资料的同向性,进一步影响地震资料的分辨率和偏移归位的准确性,最终必然会影响到整个资料处理的质量。为此我们发明了一种新的基于直达波反演的三维空间二次定位方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种方法效果可靠、操作简单易实现的基于直达波反演的三维空间二次定位方法。
本发明的目的可通过如下技术措施来实现:基于直达波反演的三维空间二次定位方法,该基于直达波反演的三维空间二次定位方法包括:步骤1,加载观测***,并进行直达波的拾取;步骤2,求取实际炮检距;步骤3,列出关于检波点的坐标方程,并进行坐标方程的求解;步骤4,去除坐标方程的无效解;步骤5,求取检波点二次定位坐标;步骤6,进行坐标替换。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,加载观测***时,对野外采集的地震数据和测量成果进行合并。
在步骤1中,进行直达波拾取时,对加载观测***的单炮用海水速度进行线性动校正,再进行直达波拾取,假设要求取的某个检波点D总共接收N炮信号,N≥4,第i炮在这个检波点的直达波时间为ti,第i炮的水深为zi,第i炮的坐标为(xi,yi,-zi)。
在步骤2中,记野外测定的海水的速度为v,则检波点P到第i炮的距离为v*ti。
在步骤3中,根据检波点D与炮点的距离v*ti,设检波点的坐标为(x,y,-z),z为检波点水深,列出检波点P关于第i炮的表达式:
(x-xi)2+(y-yi)2+(-z-(-zi))2=(v*ti)2
(式1)
其中,炮点坐标(xi,yi,-zi)和海水的速度为v通过野外测量得到,为已知数,ti为室内拾取直达波获得,这样就成了N个关于(x,y,z)的三元二次方程。
在步骤3中,求三个联立方程的解,此时会有三种情况产生,一种是无解,第二种是超过两个解,第三种是两个解;根据三球相交的原理,求得的解应该为两个,所以去掉产生前两种解的方程组,只保留两个解的方程组,这样组合得到个方程组因此最多可以求得2个解。
在步骤4中,采用均方差法进行计算,对所有的解求取平均值,然后再求取均方差,去除均方差较大的解,剩余的解即为有效解。
在步骤5中,采用均方差法计算获得实际检波点的坐标,对所有的有效解求取平均值,然后再求取均方差;去除均方差较大的解,对于剩余的解再求平均值,该值即是检波点二次定位的坐标。
在步骤6中,将原来不准确的检波点坐标用新求取的坐标进行替换,实现检波点二次定位,并运用线性动校正检查。
本发明中的基于直达波反演的三维空间二次定位方法,从三维空间综合考虑炮点、检波点的坐标,相对常规的二维方法结果更合理,精确更高,效果可靠。该方法流程及参数设置简单,运算速度快,操作简单易实现。
附图说明
图1为本发明的一具体实施例中检波点漂移野外原始单炮的示意图;
图2为本发明的一具体实施例中直达波拾取图;
图3为本发明的一具体实施例中二次定位前线性动校正图;
图4为本发明的一具体实施例中二次定位后线性动校正图;
图5为本发明的基于直达波反演的三维空间二次定位方法的一具体实施例的流程图;
图6为本发明的一具体实施例中二次定位前叠加剖面的示意图;
图7为本发明的一具体实施例中采用二次定位后的叠加剖面的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
直达波是地震波由炮点经过海水直接传到检波点的波场,直达波可以在单炮上通过人工拾取出来;海水的速度测出来比较容易,并且在某一工区基本是恒定的。利用直达波和海水的速度,可以直接求出炮点到检波点的距离。炮点的坐标可以通过卫星定位***直接测出来,这样就得到一个以炮点为球心,以求出的炮点到检波点的距离为半径的球面。对于一个检波点,用同样的方法,求取以第二个炮点、第三个炮点、第四个炮点为球心的球面,根据球面相交的原理,四个球面相交可以确定一个点,这个点就是检波点的坐标。
如图5所示,图5为本发明的基于直达波反演的三维空间二次定位方法的流程图。
步骤1,加载观测***。对野外采集的地震数据和测量成果进行合并。
步骤2,直达波拾取。对加载观测***的单炮用海水速度进行线性动校正(水速约为1500米/秒,此时的速度不用很精确),再进行直达波拾取,这里假设要求取的某个检波点D总共接收N炮信号(这里N≥4,因为一是实际野外施工不会小于4,二是小于4就不能满足本发明的基本原理),第i炮在这个检波点的直达波时间为ti,第i炮的水深为zi,第i炮的坐标为(xi,yi,-zi)。
步骤3,求取实际炮检距。记野外测定的海水的速度为v(海水的速度近似为1500米/秒左右,但每个地区会因含盐度、温度和水深的影响而略有不同),则检波点P到第i炮的距离为v*ti。
步骤4,列关于检波点坐标方程。根据检波点D与炮点的距离v*ti,设检波点的坐标为(x,y,-z),z为检波点水深,可以列出检波点P关于第i炮的表达式:
(x-xi)2+(y-yi)2+(-z-(-zi))2=(v*ti)2
(式1)
其中,炮点坐标(xi,yi,-zi)和海水的速度为v通过野外测量得到,为已知数,ti为室内拾取直达波获得,这样就成了N个关于(x,y,z)的三元二次方程。
步骤5,组合方程求解。根据本发明的基本原理,在无任何误差的情况下,4个炮点就可以求得一个检波点坐标。但是实际生产中,由于测量误差、直达波拾取误差的存在,方程很可能无解。解决方法是,只求三个联立方程的解,此时会有三种情况产生,一种是无解,第二种是超过两个解,第三种是两个解。根据三球相交的原理,求得的解应该为两个,所以去掉产生前两种解的方程组,只保留两个解的方程组,当然这两个解中,只有一个可能是正确的解。这样可以组合得到个方程组因此最多可以求得2个解。
步骤6,去除无效解。根据解的分布,去掉离散点。采用均方差法进行计算,对所有的解求取平均值,然后再求取均方差。去除均方差较大的解(大约是解总个数的50%),剩余的解即为有效解。
步骤7,求取检波点二次定位坐标。同样采用均方差法计算获得实际检波点的坐标,对所有的有效解求取平均值,然后再求取均方差。去除均方差较大的解(约去除总数的20%),对于剩余的解再求平均值,该值即是检波点二次定位的坐标。
步骤8,坐标替换。将原来不准确的检波点坐标用新求取的坐标进行替换,实现检波点二次定位。运用线性动校正检查,如图3、图4所示。
在应用本发明的一具体实施例中,以XX油田ZH地区三维地震资料为目标靶区,应用本方法对该资料进行处理,以验证本方法的效果,具体流程图见图5。该实际资料采用16线4炮观测***采集,地震资料时间长度7000ms,时间采样间隔为1ms,采样点数7000,每线道数为240道。采用上述方法对该资料进行处理。
1)首先进入步骤1,加载观测***,从单炮的初至能看出是否存在检波点漂移的单炮。在道号130附近,存在初至扭曲现象,说明道号附近的130检波点有漂移,如图1所示。
2)然后依据步骤2,用水速做的线性动校正单炮上进行直达波拾取,获得每一道的直达波时间。如图2所示,500ms附近的白色线为拾取的直达波时间。
3)依据步骤3-4,可以利用无漂移的单炮获取海水的速度,也可以从野外测量获得海水的速度,计算实际炮检距。再根据公式1列出所有的方程。
4)依据步骤5-6,对所有组合成方程组再求解,对所有的解去除无效解。
5)依据步骤7,利用均方差法去除偏离解,对解求取平均值得到二次定位坐标。
6)依据步骤8,应用二次定位坐标对原坐标进行替换,实现二次定位。图3为加原始的坐标的线性动校正图,可以看出,利用海水的速度做动校正,直达波没有拉平,说明坐标不准确。图4是二次定位后的线性动校正图,直达波拉平,说明得到了正确的检波点坐标。
7)图6是二次定位前对地震数据做的叠加,可以看出,由于检波点坐标不准,造成资料不同向,叠加效果较差,图7是采用本文提出二次定位后的叠加剖面,从剖面上看,同向轴连续性得到了明显增强,地震剖面效果更理想,说明本方法确实可行,效果突出。
Claims (9)
1.基于直达波反演的三维空间二次定位方法,其特征在于,该基于直达波反演的三维空间二次定位方法包括:
步骤1,加载观测***,并进行直达波的拾取;
步骤2,求取实际炮检距;
步骤3,列出关于检波点的坐标方程,并进行坐标方程的求解;
步骤4,去除坐标方程的无效解;
步骤5,求取检波点二次定位坐标;
步骤6,进行坐标替换。
2.根据权利要求1所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤1中,加载观测***时,对野外采集的地震数据和测量成果进行合并。
3.根据权利要求1所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤1中,进行直达波拾取时,对加载观测***的单炮用海水速度进行线性动校正,再进行直达波拾取,假设要求取的某个检波点D总共接收N炮信号,N≥4,第i炮在这个检波点的直达波时间为ti,第i炮的水深为zi,第i炮的坐标为(xi,yi,-zi)。
4.根据权利要求3所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤2中,记野外测定的海水的速度为v,则检波点P到第i炮的距离为v*ti。
5.根据权利要求4所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤3中,根据检波点D与炮点的距离v*ti,设检波点的坐标为(x,y,-z),z为检波点水深,列出检波点P关于第i炮的表达式:
(x-xi)2+(y-yi)2+(-z-(-zi))2=(v*ti)2
(式1)
其中,炮点坐标(xi,yi,-zi)和海水的速度为v通过野外测量得到,为已知数,ti为室内拾取直达波获得,这样就成了N个关于(x,y,z)的三元二次方程。
6.根据权利要求5所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤3中,求三个联立方程的解,此时会有三种情况产生,一种是无解,第二种是超过两个解,第三种是两个解;根据三球相交的原理,求得的解应该为两个,所以去掉产生前两种解的方程组,只保留两个解的方程组,这样组合得到个方程组因此最多可以求得个解。
7.根据权利要求1所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤4中,采用均方差法进行计算,对所有的解求取平均值,然后再求取均方差,去除均方差较大的解,剩余的解即为有效解。
8.根据权利要求1所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤5中,采用均方差法计算获得实际检波点的坐标,对所有的有效解求取平均值,然后再求取均方差;去除均方差较大的解,对于剩余的解再求平均值,该值即是检波点二次定位的坐标。
9.根据权利要求1所述的基于直达波反演的三维空间二次定位方法,其特征在于,在步骤6中,将原来不准确的检波点坐标用新求取的坐标进行替换,实现检波点二次定位,并运用线性动校正检查。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811292082.6A CN109212598B (zh) | 2018-10-31 | 2018-10-31 | 基于直达波反演的三维空间二次定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811292082.6A CN109212598B (zh) | 2018-10-31 | 2018-10-31 | 基于直达波反演的三维空间二次定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109212598A true CN109212598A (zh) | 2019-01-15 |
CN109212598B CN109212598B (zh) | 2021-05-04 |
Family
ID=64998360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811292082.6A Active CN109212598B (zh) | 2018-10-31 | 2018-10-31 | 基于直达波反演的三维空间二次定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109212598B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110780350A (zh) * | 2019-10-25 | 2020-02-11 | 中国海洋大学 | 一种海底地震检波器动态定位方法、***及检波器 |
CN112162322A (zh) * | 2020-09-28 | 2021-01-01 | 中海油田服务股份有限公司 | 一种地震勘探采集***的同步性质控方法和装置 |
CN113703050A (zh) * | 2021-08-05 | 2021-11-26 | 青岛海洋地质研究所 | 一种深海地震垂直缆二次定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101086534A (zh) * | 2006-06-07 | 2007-12-12 | 中国石油集团东方地球物理勘探有限责任公司 | 检波点二次定位方法 |
CN101852867A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种矢量合成检波点二次定位方法 |
CN101930080A (zh) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
KR101693310B1 (ko) * | 2015-11-20 | 2017-01-06 | 한국지질자원연구원 | 평면파를 이용한 이방성 매질에서의 전파형 역산방법 |
-
2018
- 2018-10-31 CN CN201811292082.6A patent/CN109212598B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101086534A (zh) * | 2006-06-07 | 2007-12-12 | 中国石油集团东方地球物理勘探有限责任公司 | 检波点二次定位方法 |
CN101852867A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种矢量合成检波点二次定位方法 |
CN101930080A (zh) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
KR101693310B1 (ko) * | 2015-11-20 | 2017-01-06 | 한국지질자원연구원 | 평면파를 이용한 이방성 매질에서의 전파형 역산방법 |
Non-Patent Citations (1)
Title |
---|
吕双等: "基于直达波的二次定位软件开发及应用", 《内蒙古石油化工》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110780350A (zh) * | 2019-10-25 | 2020-02-11 | 中国海洋大学 | 一种海底地震检波器动态定位方法、***及检波器 |
CN110780350B (zh) * | 2019-10-25 | 2021-10-01 | 中国海洋大学 | 一种海底地震检波器动态定位方法、***及检波器 |
CN112162322A (zh) * | 2020-09-28 | 2021-01-01 | 中海油田服务股份有限公司 | 一种地震勘探采集***的同步性质控方法和装置 |
CN113703050A (zh) * | 2021-08-05 | 2021-11-26 | 青岛海洋地质研究所 | 一种深海地震垂直缆二次定位方法 |
CN113703050B (zh) * | 2021-08-05 | 2022-06-17 | 青岛海洋地质研究所 | 一种深海地震垂直缆二次定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109212598B (zh) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212598A (zh) | 基于直达波反演的三维空间二次定位方法 | |
Guan et al. | An effective method for submarine buried pipeline detection via multi-sensor data fusion | |
CN110260044B (zh) | 一种海底管道定位方法 | |
CN108896040B (zh) | 天空海一体化水下潜器惯性/重力组合导航方法和*** | |
US10520631B2 (en) | Magnetic field measurement via streamer cables | |
CN109085655B (zh) | 一种水下平台重力测量方案与验证方法 | |
CN106443776B (zh) | 一种基于时间切片法的海底地震仪重定位方法 | |
CN106546954A (zh) | 一种深海超短基线定位精度检验方法 | |
CN103926627A (zh) | 水下载体地磁三分量测量方法 | |
CN107664758A (zh) | 基于长基线或超短基线组网的深海导航定位***及方法 | |
CN109738902A (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN111854699A (zh) | 一种基于无人机航测河道崩岸过程的监测方法 | |
CN103353612B (zh) | 一种地下目标物体的测量定位设备及测量定位方法 | |
CN110146924A (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
CN110309581B (zh) | 一种水下潜标位置综合校准测量点快速优化布局方法 | |
CN109633540B (zh) | 一种磁源的实时定位***及实时定位方法 | |
CN101266153B (zh) | 测绘工程类陀螺全站仪精度评定方法 | |
CN112378376B (zh) | 一种基于传感阵列和测斜仪的海底变形联合监测方法 | |
CN109782355B (zh) | Obs检波点漂移的检测方法及装置 | |
CN108761470B (zh) | 一种基于拖缆形态方程解析的目标定位方法 | |
CN109116416A (zh) | 基于三维方式的倾斜海底检波器二次定位方法 | |
CN109579845A (zh) | 一种基于电场探测阵列的船舶跟踪定位方法 | |
Zhao et al. | Study on underwater navigation system based on geomagnetic match technique | |
Li et al. | Compensation method for the carrier magnetic interference of underwater magnetic vector measurement system | |
CN109490868A (zh) | 一种基于分布式垂直线列阵的海上目标运动分析方法 |
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 |