CN112505750B - 一种深拖多道地震拖缆姿态确定方法及处理终端 - Google Patents
一种深拖多道地震拖缆姿态确定方法及处理终端 Download PDFInfo
- Publication number
- CN112505750B CN112505750B CN202011297248.0A CN202011297248A CN112505750B CN 112505750 B CN112505750 B CN 112505750B CN 202011297248 A CN202011297248 A CN 202011297248A CN 112505750 B CN112505750 B CN 112505750B
- Authority
- CN
- China
- Prior art keywords
- detector
- seismic
- travel time
- streamer
- deep
- 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
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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
- G01V1/201—Constructional details of seismic cables, e.g. streamers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3817—Positioning of seismic devices
- G01V1/3835—Positioning of seismic devices measuring position, e.g. by GPS or acoustically
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1423—Sea
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种深拖多道地震拖缆姿态确定方法及处理终端,所述方法包括如下步骤:步骤1:对每一个炮道集,拾取出每一个检波器上的直达波旅行时、海面反射波旅行时和海底反射波旅行时;步骤2:获得地震波在海水中的传播平均速度;步骤3:计算出每一个检波器相对于震源的水平距离和高度差;步骤4:按对应的方程组求解出每个检波器的三维坐标,得到每一个检波器的三维坐标,从而相当于得到拖缆的姿态。本发明不需要二维假设,可以获得拖缆三维的空间形态。同时利用直达波、海底反射波、海面反射波的旅行时来确定检波器与震源之间的相对位置,提高相对位置的准确性。
Description
技术领域
本发明涉及海洋勘探中拖缆姿态确定技术领域,具体涉及一种深拖多道地震拖缆姿态确定方法及处理终端。
背景技术
深拖多道地震是一种将震源和拖缆布设在距离海底上方100米左右进行海底调查的技术,如图1所示,图1中带有震源的拖***于距离海底100m(米)左右的上方,拖体与拖缆连接,拖缆上设置有若干检波器,拖缆呈理想的水平直线状态,拖缆经过拖体后通过电缆与位于海表面的地震船连接。在常规的海面拖缆地震勘探中,由于勘探***与深海海底之间存在一个厚厚的水层,导致地震波的菲涅尔带较大,影响地震数据的分辨率。而深拖多道地震正好能够减小菲涅尔带,进而提高分辨率。对于海面拖缆地震,震源的位置由GPS精确确定,拖缆的形态和位置由水鸟、罗盘、尾标等控制。
深拖地震勘探存在两个难题:(1)连接深拖***和拖船(即图1中的地震船)的电缆长达1km以上,且需要根据距离海底的高度进行调整,由此深拖***的深度和相对拖船的距离及方向在整个测线上并不是固定的。(2)由于受拖曳速度、洋流、海水浮力变化等因素的影响,深拖多道地震的拖缆的几何形态在拖曳过程中无法保持为理想的水平直线,而是处于一种不断变化的非线性状态。图2为实际深拖拖缆采集时的形态,其拖缆(即图中的接收缆)在三维空间中呈弯曲状态,拖缆上的各个检波器并不在一条理想的水平直线上。
现有的技术中,通过声学超长基线技术能够较好地确定深拖震源的位置,但深拖拖缆的几何形态却难以确定,而深拖拖缆的集合形态不能准确确定也即意味着后续利用拖缆采集的数据也导致出现偏差,因此需要确定拖缆的三维空间上的姿态,也即需要确定出拖缆上各个检波器的三维坐标。
发明内容
针对现有技术的不足,本发明的目的之一提供一种深拖多道地震拖缆姿态确定方法,其能够解决拖缆姿态确定的问题;
本发明的目的之二提供一种处理终端,其能够解决拖缆姿态确定的问题。
实现本发明的目的之一的技术方案为:一种深拖多道地震拖缆姿态确定方法,包括如下步骤:
步骤1:对每一个炮道集,拾取出每一个检波器上的直达波旅行时td、海面反射波旅行时tss和海底反射波旅行时tsf;
步骤2:获得地震波在海水中的传播平均速度v;
步骤3:根据直达波旅行时td、海面反射波旅行时tss、海底反射波旅行时tsf和地震波在海水中的传播平均速度v计算出每一个检波器相对于震源的水平距离d和高度差h;
步骤4:按方程组①求解出每个检波器的三维坐标(x′,y′,h′),得到每一个检波器的三维坐标,从而相当于得到拖缆的姿态,
式中,(xs,ys,H)表示震源的三维坐标,θ表示检波器所在拖缆的方位角。
进一步地,所述计算出每一个检波器相对于震源的水平距离d和高度差h,其具体实现过程按公式②和公式③计算得到:
式中,(x1,0)表示检波器接收到震源的海面反射波在海面的反射点的坐标,(x2,h2)表示检波器的坐标,(x3,h3)表示检波器接收到震源的海底反射波在海底的反射点的坐标,H表示震源到海面的垂直距离。
实现本发明的目的之二的技术方案为:一种处理终端,其包括:
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行所述的深拖多道地震拖缆姿态确定方法的步骤。
本发明的有益效果为:本发明不需要二维假设,可以获得拖缆三维的空间形态。同时利用直达波、海底反射波、海面反射波的旅行时来确定检波器与震源之间的相对位置,提高相对位置的准确性。将海底反射波的反射点深度当做未知数来求解,通过对比求解的结果以及多波束深度数据,可以评估和质控计算结果的准确性。
附图说明
图1为深拖多道地震观测示意图;
图2为实际深拖多道拖缆采集时的形态;
图3为较佳实施例的流程示意图;
图4为求解检波器与震源相对位置的示意图;
图5为求解检波器的三维坐标时的示意图;
图6为处理终端的示意图。
具体实施方式
下面,结合附图以及具体实施方案,对本发明做进一步描述。
如图2-图5所示,一种深拖多道地震拖缆姿态确定方法,包括如下步骤:
步骤2:获得地震波在海水中的传播平均速度v。可以通过分析实际勘探采集时所在区域的海面多道地震数据来获得地震波在海水中传播的平均速度,其为现有技术,在此不进行赘述。
步骤3:计算出每一个检波器相对于震源的水平距离d和高度差h。水平距离也即是检波器和震源两个点之间在水平方向上的距离,而高度差也即是以垂直于海面的竖直方向上检波器和震源两个点之间的深度差。
参考图4,图中的i表示第i个检波器,求解检波器相对于震源的水平距离d和高度差h具体包括,
过震源点和检波点的连线作垂直于海面的平面(如图2所示),沿着该平面与海面的交线建立x轴,在该平面内沿着垂直于x轴的方向建立y轴,从而得到该平面内的二维直角坐标系。设震源的坐标为(0,H),H也即是震源到海面的垂直距离,为已知量,检波器的坐标为(x2,h2),检波器接收到震源的海面反射波在海面的反射点的坐标为(x1,0),检波器接收到震源的海底反射波在海底的反射点的坐标为(x3,h3)。
根据直达波旅行时的计算公式,可得到:
x2 2+(h2-H)2=(tdv)2
根据海面反射波旅行时的计算公式,可得到:
x1 2+H2+(x2-x1)2+h2 2=(tssv)2
根据海底反射波旅行时的计算公式,可得到:
x3 2+(h3-H)2+(x2-x3)2+(h3-h2)2=(tsfv)2
根据海面反射波入射角等于反射角,可得到:
根据海底反射波入射角等于反射角,可得到:
将以上公式联合起来,即得到方程组①:
方程组①中,x1,x2,h2,x3,h3均为未知数,其余均为已知数,所以可以求解方程组,也即属于一个正定方程。考虑到方程组中各个方程均是非线性的,本实施例采用最小二乘最优化的方法求解五个未知数。最终得到检检波器相对于震源的水平距离d和高度差h:
d=x2
h=H-h2
步骤4:参考图5,图中的i表示第i个检波器,根据检波器相对于震源的水平距离d、高度差h、震源的三维坐标(xs,ys,H)和拖缆上的方位传感器所记录的方位角度θ,记录的方位角度θ表示检波器所在拖缆的方位角,也即表示检波器与震源的连线投影到海平面的投射线相对于以海面为轴线的夹角,按方程组②求解出每个检波器的三维坐标(x′,y′,h′),得到每一个检波器的三维坐标,从而相当于得到拖缆的姿态。
按照上述步骤即可得到每一个检波器的三维坐标,从而得到了整个拖缆的三维空间形态,得到拖缆的三维空间形态之后,就可以对深拖多道地震数据进行后续处理了。
本发明不需要二维假设,可以获得拖缆三维的空间形态。同时利用直达波、海底反射波、海面反射波的旅行时来确定检波器与震源之间的相对位置,提高相对位置的准确性。将海底反射波的反射点深度当做未知数来求解,通过对比求解的结果以及多波束深度数据,可以评估和质控计算结果的准确性。
如图6所示,本发明还提供一种处理终端100,其包括:
存储器101,用于存储程序指令;
处理器102,用于运行所述程序指令,以执行所述深拖多道地震拖缆姿态确定方法的步骤。
本说明书所公开的实施例只是对本发明单方面特征的一个例证,本发明的保护范围不限于此实施例,其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说,可根据以上描述的技术方案以及构思,做出其它各种相应的改变以及变形,而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。
Claims (2)
1.一种深拖多道地震拖缆姿态确定方法,其特征在于,包括如下步骤:
步骤1:对每一个炮道集,拾取出每一个检波器上的直达波旅行时td、海面反射波旅行时tss和海底反射波旅行时tsf;
步骤2:获得地震波在海水中的传播平均速度v;
步骤3:根据直达波旅行时td、海面反射波旅行时tss、海底反射波旅行时tsf和地震波在海水中的传播平均速度v计算出每一个检波器相对于震源的水平距离d和高度差h;
步骤4:按方程组①求解出每个检波器的三维坐标(x′,y′,h′),得到每一个检波器的三维坐标,从而相当于得到拖缆的姿态,
式中,(xs,ys,H)表示震源的三维坐标,θ表示检波器所在拖缆的方位角,
所述计算出每一个检波器相对于震源的水平距离d和高度差h,其具体实现过程按公式②和公式③计算得到:
式中,(x1,0)表示检波器接收到震源的海面反射波在海面的反射点的坐标,(x2,h2)表示检波器的坐标,(x3,h3)表示检波器接收到震源的海底反射波在海底的反射点的坐标,H表示震源到海面的垂直距离。
2.一种处理终端,其特征在于,其包括:
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行如权利要求1所述的深 拖多道地震拖缆姿态确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011297248.0A CN112505750B (zh) | 2020-11-18 | 2020-11-18 | 一种深拖多道地震拖缆姿态确定方法及处理终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011297248.0A CN112505750B (zh) | 2020-11-18 | 2020-11-18 | 一种深拖多道地震拖缆姿态确定方法及处理终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112505750A CN112505750A (zh) | 2021-03-16 |
CN112505750B true CN112505750B (zh) | 2021-06-08 |
Family
ID=74958080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011297248.0A Active CN112505750B (zh) | 2020-11-18 | 2020-11-18 | 一种深拖多道地震拖缆姿态确定方法及处理终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112505750B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112987080B (zh) * | 2021-04-22 | 2022-10-25 | 自然资源部第一海洋研究所 | 海底多波震源及海底探测*** |
CN113703050B (zh) * | 2021-08-05 | 2022-06-17 | 青岛海洋地质研究所 | 一种深海地震垂直缆二次定位方法 |
CN114518597B (zh) * | 2022-03-08 | 2023-05-02 | 中国海洋大学 | 基于海底反射与其电缆虚反射校正地震不平电缆的方法 |
CN114460649B (zh) * | 2022-04-14 | 2022-06-24 | 自然资源部第一海洋研究所 | 一种深海近底拖曳式多道地震接收阵列形态重建方法 |
CN115308800B (zh) * | 2022-08-29 | 2024-07-12 | 南方海洋科学与工程广东省实验室(广州) | 一种利用海底反射波走时和地形数据定位海底地震仪的方法及处理终端 |
CN115826056B (zh) * | 2023-02-20 | 2023-06-16 | 山东科技大学 | 深拖式高分辨率多道地震拖缆水听器阵列高精度定位方法 |
CN115980856B (zh) * | 2023-03-21 | 2023-06-30 | 自然资源部第一海洋研究所 | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3831136A (en) * | 1972-05-05 | 1974-08-20 | Chevron Res | Method of initiating and collecting seismic data related to strata underlying bodies of water using a continuously moving seismic exploration system located on a single boat |
CN1664617A (zh) * | 2004-03-03 | 2005-09-07 | Pgs美洲公司 | 组合海洋地震拖缆中的压力传感器和质点运动传感器的信号的*** |
CN106443792A (zh) * | 2016-10-14 | 2017-02-22 | 中国船舶工业***工程研究院 | 一种水下多缆声学网络定位方法 |
CN209433021U (zh) * | 2019-01-30 | 2019-09-24 | 自然资源部第二海洋研究所 | 一种基于拖缆搭载的海洋地震探测*** |
CN111273344A (zh) * | 2020-03-02 | 2020-06-12 | 广州海洋地质调查局 | 一种基于续至折射波的层析反演方法及处理终端 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8010291B2 (en) * | 2008-11-08 | 2011-08-30 | Westerngeco L.L.C. | Processing measurement data that involves computing a derivative of the measurement data that represents a subterranean structure |
-
2020
- 2020-11-18 CN CN202011297248.0A patent/CN112505750B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3831136A (en) * | 1972-05-05 | 1974-08-20 | Chevron Res | Method of initiating and collecting seismic data related to strata underlying bodies of water using a continuously moving seismic exploration system located on a single boat |
CN1664617A (zh) * | 2004-03-03 | 2005-09-07 | Pgs美洲公司 | 组合海洋地震拖缆中的压力传感器和质点运动传感器的信号的*** |
CN106443792A (zh) * | 2016-10-14 | 2017-02-22 | 中国船舶工业***工程研究院 | 一种水下多缆声学网络定位方法 |
CN209433021U (zh) * | 2019-01-30 | 2019-09-24 | 自然资源部第二海洋研究所 | 一种基于拖缆搭载的海洋地震探测*** |
CN111273344A (zh) * | 2020-03-02 | 2020-06-12 | 广州海洋地质调查局 | 一种基于续至折射波的层析反演方法及处理终端 |
Non-Patent Citations (6)
Title |
---|
Dynamic analysis of deep-towed seismic array based on relative-velocity-element-frame;Xiangqian Zhu et al.;《Ocean Engineering》;20201027;第1-14页 * |
Mistie Problem and Correction Method of Single Channel Seismic Data in Shallow Sea;Xi Chen et al.;《2nd International Symposium on Resource Exploration and Environmental Science》;20181231;第1-8页 * |
Source and receiver geometry corrections for deep towed multichannel seismic data;R.Wali et al.;《GEOPHYSICAL RESEARCH LETTERS》;19990701;第26卷(第13期);第1993-1996页 * |
东沙海域复杂海况下单源单缆长排列地震资料处理;邢涛等;《热带海洋学报》;20200731;第39卷(第4期);第91-99页 * |
海洋短排列多道反射地震数据观测***重定义与沉放深度校正;张圣等;《海洋学报》;20200331;第42卷(第3期);第140-148页 * |
深拖式多道高分辨率地震探测***在南海首次应用;魏峥嵘等;《石油地球物理勘探》;20201031;第55卷(第5期);第965-972页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112505750A (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112505750B (zh) | 一种深拖多道地震拖缆姿态确定方法及处理终端 | |
US8760967B2 (en) | Generating an angle domain common image gather | |
Singh et al. | Microbathymetric mapping from underwater vehicles in the deep ocean | |
WO2009048683A2 (en) | Controlling seismic source elements based on determining a three-dimensional geometry of the seismic source elements | |
CN113640808B (zh) | 浅水区海底电缆埋深探测方法及装置 | |
US11802986B2 (en) | Hybrid ocean bottom seismic receiver and streamer seismic data acquisition using wide towed sources | |
US9116257B2 (en) | Method and apparatus for sea current aided, enhanced azimuth seismic data acquisition | |
US10520631B2 (en) | Magnetic field measurement via streamer cables | |
CN107664758A (zh) | 基于长基线或超短基线组网的深海导航定位***及方法 | |
CN111551942B (zh) | 一种基于解卷积算法的水下自主航行器回坞方法 | |
RU2608301C2 (ru) | Система и способ 3d исследования морского дна для инженерных изысканий | |
CN112904428B (zh) | 海洋浅地层剖面探测***及方法 | |
CN108761470B (zh) | 一种基于拖缆形态方程解析的目标定位方法 | |
CN112904316B (zh) | 一种针对海面波动的机载激光测深数据折射误差改正方法 | |
CN110907937A (zh) | 一种基于“t”型阵的掩埋物合成孔径三维成像方法 | |
CN115308800B (zh) | 一种利用海底反射波走时和地形数据定位海底地震仪的方法及处理终端 | |
US8477560B2 (en) | Determining a position of an object of a positioning network associated with a marine survey arrangement | |
CN113703050B (zh) | 一种深海地震垂直缆二次定位方法 | |
US8634270B2 (en) | Determining sea conditions in marine seismic spreads | |
KR102425484B1 (ko) | 연근해 3차원 고주파 지층탐사 자료처리방법 | |
CN116609834B (zh) | 一种基于海洋垂直缆地震勘探的数据处理方法 | |
EP2778717A2 (en) | Wavefield Regularization by 3-D Wavefield Decomposition for Geophysical Data | |
Guandong et al. | Method and application of offset vector summation in node repositioning based on refractions | |
Hou et al. | Multi-beam line observation optimization strategy based on random variable unitary iterative model | |
Wang | Ocean Exploration Model and Analysis Based on Multi-Beam Sounding Technology |
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 |