CN115980856A - 适用复杂海底界面的深拖式多道地震炮检精确定位方法 - Google Patents
适用复杂海底界面的深拖式多道地震炮检精确定位方法 Download PDFInfo
- Publication number
- CN115980856A CN115980856A CN202310272911.9A CN202310272911A CN115980856A CN 115980856 A CN115980856 A CN 115980856A CN 202310272911 A CN202310272911 A CN 202310272911A CN 115980856 A CN115980856 A CN 115980856A
- Authority
- CN
- China
- Prior art keywords
- deep
- towed
- seismic
- point
- shot
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 112
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005284 excitation Effects 0.000 claims abstract description 80
- 238000004364 calculation method Methods 0.000 claims abstract description 52
- 239000013535 sea water Substances 0.000 claims abstract description 36
- 239000000523 sample Substances 0.000 claims abstract description 30
- 238000005457 optimization Methods 0.000 claims abstract description 29
- 238000012512 characterization method Methods 0.000 claims abstract description 18
- 230000001360 synchronised effect Effects 0.000 claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 41
- 230000035515 penetration Effects 0.000 claims description 9
- 238000000926 separation method Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了适用复杂海底界面的深拖式多道地震炮检精确定位方法,包括:计算深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量;计算各炮集检波点接收到的直达波旅行时理论计算值、海面反射波旅行时理论计算值和海底反射波旅行时理论计算值;构建复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数;求解目标函数,完成深拖式多道地震探测系统的炮检精确定位。本发明不依赖水听器阵列配置的姿态传感器和深度传感器数据,能够同步获取各炮集阵列几何形态表征参数和海水声速的优化反演结果,有效应对地震波旅行时拾取的不确定性和任意复杂海底界面条件。
Description
技术领域
本发明涉及海洋地球物理勘探技术领域,具体涉及适用复杂海底界面的深拖式多道地震炮检精确定位方法。
背景技术
深拖式多道地震探测技术是一种将地震高频震源和地震接收电缆均置于近海底位置拖曳的、适用于深水工作环境的高分辨率多道地震探测方式,自20世纪80年代以来,世界各国研发出来的深拖式多道地震探测系统已经在天然气水合物系统评估、游离气和气烟囱表征、海底麻坑识别、滑坡灾害评价等方面取得了成功的应用效果。
准确的炮检位置关系是深拖式多道地震数据成像质量的前提条件,对于深拖式高分辨率地震探测系统而言,微小的炮检定位误差引起的地震记录双曲线失真会造成速度分析速度谱聚焦度变差,进而导致叠加处理时出现相消干扰,影响地震成像剖面的分辨率和信噪比。近些年,为了提高深拖式多道地震探测系统数据成像质量,充分挖掘系统的探测潜力,许多学者相继提出和发展了多种适用于深拖式多道地震探测系统的炮检位置反演方法。现有技术多依赖于深拖多道地震接收电缆配置的姿态传感器或深度传感器数据,并且没有考虑海水声速的横向不均匀性及海底界面的复杂性,中国专利(CN114460649A)公开了一种深海近底拖曳式多道地震接收阵列形态重建方法,但是该方法也没有考虑海水声速的横向不均匀性及海底界面的复杂性,以上问题都会导致炮检定位结果不够准确,影响最终地震成像剖面的信噪比和分辨率,甚至出现构造假象。
发明内容
针对现有技术存在的问题,本发明提供适用复杂海底界面的深拖式多道地震炮检精确定位方法,其仅基于直达波、海底反射波和海面反射波走时信息,不依赖于水听器阵列配置的姿态传感器和深度传感器数据,且对任意复杂海底界面均适用,能够有效提高线列阵拖体的几何形态反演精度和可靠性。
适用复杂海底界面的深拖式多道地震炮检精确定位方法,所述方法包括:
分别计算深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量;
分别计算各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值、海面反射波旅行时理论计算值和海底反射波旅行时理论计算值;
构建复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数;
求解目标函数,实现深拖式多道地震炮检精确定位。
在一种可能的实现方式中,深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量的计算公式如下:
其中,表示深拖震源距离深拖地震接收电缆拖曳点的水平方向偏移量,表示深拖震源距离深拖地震接收电缆拖曳点的垂直方向偏移量,表示深拖地震接收电缆拖曳点和首个检波点沿深拖地震接收电缆的间隔距离,表示相邻检波点沿深拖地震接收电缆的间隔距离,表示深拖地震接收电缆分割成的第节线段的俯仰角,表示深拖地震接收电缆分割成的第1节线段的俯仰角;以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第1个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示第1个检波点相对深拖震源激发点位置的垂直方向偏移量,,为深拖地震接收电缆的检波点总个数。
在一种可能的实现方式中,各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值和海面反射波旅行时理论计算值的计算公式如下:
其中,表示当前炮集深拖震源激发点位置的深拖震源入水深度观测值,表示当前炮集深拖震源激发点位置的海水声速;以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,,为深拖地震接收电缆的检波点总个数。
在一种可能的实现方式中,各炮集中深拖地震接收电缆各检波点接收到的海底反射波旅行时理论计算值的计算公式如下:
其中,表示当前炮集深拖震源激发点位置的深拖震源入水深度观测值,表示当前炮集深拖震源激发点位置的海水声速;以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示当前炮集第个海底界面水深变化曲线离散点的横向坐标,表示当前炮集第个海底界面水深变化曲线离散点的水深值,表示当设定第个海底界面水深变化曲线离散点为反射点时,第个检波点接收到的海底反射波传播路径的距离,,为海底界面水深变化曲线离散点的总个数;表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值;,为深拖地震接收电缆的检波点总个数。
在一种可能的实现方式中,复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数如下:
其中,表示使其后的表达式达到最小值时变量的取值,待优化参数,分别表示深拖地震接收电缆分割成的各节线段俯仰角,构成接收阵列几何形态表征参数,表示当前炮集深拖震源激发点位置的海水声速优化参数,表示当前炮集深拖震源激发点位置的海水声速,表示深拖地震接收电缆的检波点序号,,为深拖地震接收电缆的检波点总个数;表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的直达波旅行时实际观测值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时实际观测值,表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海底反射波旅行时实际观测值,表示目标函数最优解,,表示接收阵列几何形态表征参数最优解,表示当前炮集深拖震源激发点位置的海水声速优化参数最优解。
在一种可能的实现方式中,基于公式(7)获得目标函数最优解,将代入公式(1)中的,将代入公式(1)中的,求解出深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量,实现深拖式多道地震炮检精确定位。
基于以上发明内容,相对于现有技术,本发明实现了以下技术效果:
(1)本发明采用最短路径算法对深拖式多道地震探测系统接收到的海底反射波旅行时进行准确求取,不需要对海底界面做任何近似表示,能够适用于任意复杂海底界面情况。
(2)本发明基于直达波、海面反射波和海底反射波旅行时共同构建深拖式多道地震探测系统炮检定位目标函数,能够应对地震波旅行时拾取不确定性带来的影响,能够有效提高反演算法的稳定性。
(3)本发明考虑海水声速横向非均质性对深拖多道地震接收电缆阵列几何形态反演精度的影响,将海水声速优化参数引入目标函数,采用非线性优化算法实现了深拖多道地震接收阵列几何形态表征参数和海水声速优化参数的同步反演,能够有效提高深拖式多道地震探测系统炮点-检波点相对位置关系的反演精度和可靠性。
本发明所述方法不依赖水听器阵列配置的姿态传感器和深度传感器数据,能够同步获取各炮集的阵列几何形态表征参数和海水声速的优化反演结果,能够有效应对地震波旅行时拾取的不确定性和任意复杂海底界面条件。
附图说明
为了更清楚说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本申请实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法的流程示意图;
图2为本发明实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法中深拖式多道地震探测系统接收到的直达波和海面反射波传播路径示意图;
图3为本发明实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法中深拖式多道地震探测系统接收到的海底反射波传播路径示意图;
图4为本发明实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法中作业测线区段对应的海底界面水深变化曲线;
图5为本发明实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法中深拖式多道地震探测系统炮点-检波点精确定位结果。
具体实施方式
为了更好的理解本发明的技术方案,下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
参见图1,为本申请实施例提供的适用复杂海底界面的深拖式多道地震炮检精确定位方法的流程示意图。如图1所示,所述方法具体包括:
步骤S1:基于深拖式多道地震探测系统的工程设计参数,对深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量、垂直方向偏移量进行表征;
步骤S2:根据地震射线理论,结合深拖震源入水深度观测值和深拖震源离底高度观测值,对各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值、海面反射波旅行时理论计算值和海底反射波旅行时理论计算值进行表征;
步骤S3:构建复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数;
步骤S4:求解目标函数,完成深拖式多道地震探测系统的炮检精确定位。
需要特别说明的是,在上述深拖式多道地震炮检精确定位方法前,首先需要进行数据采集。具体包括:对实际炮集记录进行人工拾取获得深拖地震接收电缆检波点接收到的直达波旅行时实际观测值、海面反射波旅行时实际观测值和海底反射波旅行时实际观测值;读取深度计的值获取各炮集深拖震源激发点位置的深拖震源入水深度观测值,读取高度计的值获取各炮集深拖震源激发点位置的深拖震源离底高度观测值。
进一步地,所述步骤S1具体为:
在只考虑二维情况下,将深拖震源和深拖地震接收电缆组成的线列阵拖体以各检波点位置为节点分割为多节线段,基于各节线段长度及其俯仰角,对深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量、垂直方向偏移量进行表征,获得深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量、垂直方向偏移量的表达式:
其中,表示深拖震源距离深拖地震接收电缆拖曳点的水平方向偏移量,表示深拖震源距离深拖地震接收电缆拖曳点的垂直方向偏移量,表示深拖地震接收电缆拖曳点和首个检波点沿深拖地震接收电缆的间隔距离,表示相邻检波点沿深拖地震接收电缆的间隔距离,表示深拖地震接收电缆分割成的第节线段的俯仰角,表示深拖地震接收电缆分割成的第1节线段的俯仰角;以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第1个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示第1个检波点相对深拖震源激发点位置的垂直方向偏移量,,为深拖地震接收电缆的检波点总个数。
本实施例中,深拖式多道地震探测系统的工程设计参数:深拖地震接收电缆的检波点总个数为48,即;深拖震源距离深拖地震接收电缆拖曳点的水平方向偏移量为2.0m,深拖震源距离深拖地震接收电缆拖曳点的垂直方向偏移量为0.6m;相邻检波点沿深拖地震接收电缆的间隔距离为3.125m;深拖地震接收电缆拖曳点和首个检波点沿深拖地震接收电缆的间隔距离为12.5m;分别表示深拖地震接收电缆分割成的各节线段的俯仰角。
进一步地,所述步骤S2具体为:
参见图2,为本发明实施例提供的深拖式多道地震探测系统直达波和海面反射波传播路径示意图。如图2所示,基于地震射线理论可以得出炮检位置和直达波、海面反射波旅行时的关系,由此获得各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值和海面反射波旅行时理论计算值,采用以下公式计算:
其中,表示当前炮集深拖震源激发点位置的深拖震源入水深度观测值,表示当前炮集深拖震源激发点位置的海水声速,以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,,为深拖地震接收电缆的检波点总个数。
区别于水平海底界面和倾斜海底界面,针对复杂海底界面情况,由于受海底反射波传播路径确定问题的影响,不能用简单的几何关系对其海底反射波旅行时进行表征。由各炮集深拖震源激发点位置的深拖震源入水深度观测值和深拖震源离底高度观测值经平滑处理后相加,获得深拖式多道地震探测系统作业测线区段对应的海底界面水深变化曲线。选取当前炮集深拖震源和深拖地震接收电缆下方的海底界面水深信息,横向范围从深拖震源正下方位置向深拖地震接收电缆方向截取一定距离,将该段海底界面水深变化曲线按照横向间隔0.5m离散化,获得离散化后的第个海底界面水深变化曲线离散点的水深值,为海底界面水深变化曲线离散点的序号,则当前炮集海底界面第个海底界面水深变化曲线离散点的横向坐标为
其中,,为海底界面水深变化曲线离散点的总个数。对于深拖地震接收电缆的每个检波点,连接深拖震源、海底界面水深变化曲线离散点和深拖地震接收电缆检波点,可以获得个可能的海底反射波传播路径,也就是说,表示当设定第个海底界面水深变化曲线离散点为反射点时,第个检波点接收到的海底反射波传播路径的距离。如图3所示,根据最短路径算法,确定其中最短传播距离对应的海底反射波传播路径为真实的海底反射波传播路径,由此获得各炮集中深拖地震接收电缆各检波点接收到的海底反射波旅行时理论计算值的表达式:
其中,表示当前炮集深拖震源激发点位置的深拖震源入水深度观测值,表示当前炮集深拖震源激发点位置的海水声速。以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示当前炮集第个海底界面水深变化曲线离散点的横向坐标,表示当前炮集第个海底界面水深变化曲线离散点的水深值,则表示当前炮集第个海底界面水深变化曲线离散点的纵坐标,表示当设定第个海底界面水深变化曲线离散点为反射点时,第个检波点接收到的海底反射波传播路径的距离,,为海底界面水深变化曲线离散点的总个数。表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值。,为深拖地震接收电缆的检波点总个数。
本实施例中,以第75炮为例,第75炮炮集深拖震源激发点位置的深拖震源入水深度观测值为1415.52m,深拖震源离底高度观测值为98.24m,以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值,,48为深拖地震接收电缆的检波点总个数。
图4为本实施例中一共221个炮集深拖震源激发点位置的深拖震源入水深度观测值和深拖震源离底高度观测值经平滑处理后相加,获得的深拖式多道地震探测系统作业测线区段对应的海底界面水深变化曲线,如图4所示,圆点标记表示第75炮对应的海底界面水深值。
以第75炮为例,选取当前炮集深拖震源和深拖地震接收电缆下方的海底界面水深信息,横向范围从深拖震源正下方位置向深拖接收电缆方向截取一定距离(本实施例中深拖地震接收电缆有效工作长度150m,本实施例中将一定距离设为100m),将该段海底界面水深变化曲线按照横向间隔0.5m离散化,获得离散化后的海底界面水深变化曲线离散点的水深值,为海底界面水深变化曲线离散点的序号,则当前炮集海底界面第个海底界面水深变化曲线离散点的横向坐标为
其中,。201为海底界面水深变化曲线离散点的总个数。对于深拖地震接收电缆的每个检波点,连接深拖震源、海底界面水深变化曲线离散点和深拖地震接收电缆检波点可以获得201个可能的海底反射波传播路径,根据最短路径算法,确定其中最短传播距离对应的海底反射波传播路径为真实的海底反射波传播路径。
进一步地,所述步骤S3具体为:
当前炮集深拖震源激发点位置的海水声速对于深拖地震接收阵列几何形态反演精度影响较大,考虑到海水声速的横向不均匀性,使在1480-1500m/s之间优化取值,将海水声速也作为优化参数在目标函数求解中同步反演。因此,在上述各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值、海面反射波旅行时理论计算值和海底反射波旅行时理论计算值的表达式的基础上,构建复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数如下:
其中,表示使其后的表达式达到最小值时变量的取值,待优化参数,分别表示深拖地震接收电缆分割成的各节线段俯仰角,构成接收阵列几何形态表征参数,表示当前炮集深拖震源激发点位置的海水声速优化参数,表示当前炮集深拖震源激发点位置的海水声速,表示深拖地震接收电缆的检波点序号,,为深拖地震接收电缆的检波点总个数,表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的直达波旅行时实际观测值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时实际观测值,表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海底反射波旅行时实际观测值,表示目标函数最优解,,表示接收阵列几何形态表征参数最优解,表示当前炮集深拖震源激发点位置的海水声速优化参数最优解。
本实施例中深拖地震接收电缆的检波点总个数为48,待优化参数,其中,分别表示深拖地震接收电缆分割成的各节线段俯仰角,是当前炮集深拖震源激发点位置的海水声速优化参数。当前炮集第个检波点接收到的直达波旅行时实际观测值、海面反射波旅行时实际观测值和海底反射波旅行时实际观测值通过对实际炮集记录进行人工拾取获得。
进一步地,所述步骤S4具体为:
基于公式(7),利用粒子群优化算法求解目标函数,获得目标函数最优解,将代入公式(1)中的,将代入公式(1)中的,求解出深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量,实现深拖式多道地震炮检精确定位。
在本实施例中,基于公式(7),利用粒子群优化算法求解目标函数,获得目标函数最优解,将代入公式(1)中的,将代入公式(1)中的,求解出深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量,实现深拖式多道地震炮检精确定位。如图5所述,为本实施例第75炮的深拖式多道地震探测系统炮点-检波点位置信息反演结果,即深拖式多道地震炮检精确定位结果,其中星形标记为深拖震源位置,圆形标记为深拖地震接收电缆48个检波点的位置;横坐标表示以实施例第1炮深拖震源位置为参照原点计算的横向距离,定义拖曳作业行进方向为负方向;纵坐标表示深拖震源及深拖地震接收电缆检波点的入水深度。
综上所述,本发明采用最短路径算法对深拖式多道地震探测系统接收到的海底反射波旅行时进行准确求取,方法不需要对海底界面做任何近似表示,能够适用于任意复杂海底界面情况;基于直达波、海面反射波和海底反射波旅行时共同构建深拖式多道地震探测系统炮检定位目标函数,能够应对地震波旅行时拾取不确定性带来的影响,能够有效提高反演算法的稳定性;考虑海水声速横向非均质性对深拖多道地震接收电缆阵列几何形态反演精度的影响,将海水声速优化参数引入目标函数,采用非线性优化算法实现了深拖多道地震接收阵列几何形态表征参数和海水声速优化参数的同步反演,有效提高了深拖式多道地震探测系统炮点-检波点相对位置关系的反演精度和可靠性。
以上所述,仅为本申请的具体实施方式,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。本申请的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.适用复杂海底界面的深拖式多道地震炮检精确定位方法,其特征在于,所述方法包括:
分别计算深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量;
分别计算各炮集中深拖地震接收电缆各检波点接收到的直达波旅行时理论计算值、海面反射波旅行时理论计算值和海底反射波旅行时理论计算值;
构建复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数;
求解目标函数,实现深拖式多道地震炮检精确定位。
2.根据权利要求1所述的适用复杂海底界面的深拖式多道地震炮检精确定位方法,其特征在于,深拖地震接收电缆各检波点相对深拖震源激发点位置的水平方向偏移量和垂直方向偏移量的计算公式如下:
其中,表示深拖震源距离深拖地震接收电缆拖曳点的水平方向偏移量,表示深拖震源距离深拖地震接收电缆拖曳点的垂直方向偏移量,表示深拖地震接收电缆拖曳点和首个检波点沿深拖地震接收电缆的间隔距离,表示相邻检波点沿深拖地震接收电缆的间隔距离,表示深拖地震接收电缆分割成的第节线段的俯仰角,表示深拖地震接收电缆分割成的第1节线段的俯仰角;以深拖震源激发点位置为原点,表示第个检波点相对深拖震源激发点位置的水平方向偏移量,表示第1个检波点相对深拖震源激发点位置的水平方向偏移量,表示第个检波点相对深拖震源激发点位置的垂直方向偏移量,表示第1个检波点相对深拖震源激发点位置的垂直方向偏移量,,为深拖地震接收电缆的检波点总个数。
4.根据权利要求1所述的适用复杂海底界面的深拖式多道地震炮检精确定位方法,其特征在于,各炮集中深拖地震接收电缆各检波点接收到的海底反射波旅行时理论计算值的计算公式如下:
6.根据权利要求2所述的适用复杂海底界面的深拖式多道地震炮检精确定位方法,其特征在于,复杂海底界面条件下深拖式多道地震接收阵列几何形态表征参数和海水声速优化参数同步优化反演的目标函数如下:
其中,表示使其后的表达式达到最小值时变量的取值,待优化参数,分别表示深拖地震接收电缆分割成的各节线段俯仰角,构成接收阵列几何形态表征参数,表示当前炮集深拖震源激发点位置的海水声速优化参数,表示当前炮集深拖震源激发点位置的海水声速,表示深拖地震接收电缆的检波点序号,,为深拖地震接收电缆的检波点总个数;表示当前炮集第个检波点接收到的直达波旅行时理论计算值,表示当前炮集第个检波点接收到的直达波旅行时实际观测值,表示当前炮集第个检波点接收到的海面反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海面反射波旅行时实际观测值,表示当前炮集第个检波点接收到的海底反射波旅行时理论计算值,表示当前炮集第个检波点接收到的海底反射波旅行时实际观测值,表示目标函数最优解,,表示接收阵列几何形态表征参数最优解,表示当前炮集深拖震源激发点位置的海水声速优化参数最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310272911.9A CN115980856B (zh) | 2023-03-21 | 2023-03-21 | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310272911.9A CN115980856B (zh) | 2023-03-21 | 2023-03-21 | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115980856A true CN115980856A (zh) | 2023-04-18 |
CN115980856B CN115980856B (zh) | 2023-06-30 |
Family
ID=85976491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310272911.9A Active CN115980856B (zh) | 2023-03-21 | 2023-03-21 | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115980856B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4992991A (en) * | 1988-10-31 | 1991-02-12 | Exxon Production Research Company | Marine seismic receiving system employing multiple depth hydrostreamer cable |
US20040076077A1 (en) * | 2001-01-25 | 2004-04-22 | Johan Robertsson | Method of processing marine seismic data and a method of seismic surveying |
US20130041616A1 (en) * | 2010-04-06 | 2013-02-14 | Cggveritas Services | Method of determining the position of a detector disposed at the bottom of the sea |
CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
US20170248721A1 (en) * | 2014-09-10 | 2017-08-31 | Cgg Services | Wave-field reconstruction using a reflection from a variable sea surface |
CN107576986A (zh) * | 2017-08-30 | 2018-01-12 | 中国石油集团东方地球物理勘探有限责任公司 | 一种确定海底地震反射波旅行时的方法和装置 |
CN110780350A (zh) * | 2019-10-25 | 2020-02-11 | 中国海洋大学 | 一种海底地震检波器动态定位方法、系统及检波器 |
CN112505750A (zh) * | 2020-11-18 | 2021-03-16 | 广州海洋地质调查局 | 一种深拖多道地震拖缆姿态确定方法及处理终端 |
CN114460649A (zh) * | 2022-04-14 | 2022-05-10 | 自然资源部第一海洋研究所 | 一种深海近底拖曳式多道地震接收阵列形态重建方法 |
CN115308800A (zh) * | 2022-08-29 | 2022-11-08 | 南方海洋科学与工程广东省实验室(广州) | 一种利用海底反射波走时和地形数据定位海底地震仪的方法及处理终端 |
-
2023
- 2023-03-21 CN CN202310272911.9A patent/CN115980856B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4992991A (en) * | 1988-10-31 | 1991-02-12 | Exxon Production Research Company | Marine seismic receiving system employing multiple depth hydrostreamer cable |
US20040076077A1 (en) * | 2001-01-25 | 2004-04-22 | Johan Robertsson | Method of processing marine seismic data and a method of seismic surveying |
US20130041616A1 (en) * | 2010-04-06 | 2013-02-14 | Cggveritas Services | Method of determining the position of a detector disposed at the bottom of the sea |
US20170248721A1 (en) * | 2014-09-10 | 2017-08-31 | Cgg Services | Wave-field reconstruction using a reflection from a variable sea surface |
CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
CN107576986A (zh) * | 2017-08-30 | 2018-01-12 | 中国石油集团东方地球物理勘探有限责任公司 | 一种确定海底地震反射波旅行时的方法和装置 |
CN110780350A (zh) * | 2019-10-25 | 2020-02-11 | 中国海洋大学 | 一种海底地震检波器动态定位方法、系统及检波器 |
CN112505750A (zh) * | 2020-11-18 | 2021-03-16 | 广州海洋地质调查局 | 一种深拖多道地震拖缆姿态确定方法及处理终端 |
CN114460649A (zh) * | 2022-04-14 | 2022-05-10 | 自然资源部第一海洋研究所 | 一种深海近底拖曳式多道地震接收阵列形态重建方法 |
CN115308800A (zh) * | 2022-08-29 | 2022-11-08 | 南方海洋科学与工程广东省实验室(广州) | 一种利用海底反射波走时和地形数据定位海底地震仪的方法及处理终端 |
Non-Patent Citations (2)
Title |
---|
XIANGQIAN ZHU等: "Dynamic analysis of deep-towed seismic array based on relative-velocity-element-frame", OCEAN ENGINEERING * |
张圣等: "海洋短排列多道反射地震数据观测系统重定义与沉放深度校正", 海洋学报 * |
Also Published As
Publication number | Publication date |
---|---|
CN115980856B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3324215B1 (en) | Method of first arrival picking of multi-channel seismic survey data | |
CN114460649B (zh) | 一种深海近底拖曳式多道地震接收阵列形态重建方法 | |
EP3078991B1 (en) | Method for swell effect and mis-tie correction in high-resolution seismic data using multi-beam echo sounder data | |
CN112505750B (zh) | 一种深拖多道地震拖缆姿态确定方法及处理终端 | |
US10520631B2 (en) | Magnetic field measurement via streamer cables | |
CN111854704A (zh) | 一种海洋地球物理综合调查系统 | |
AU2012208972B2 (en) | Quality-based steering methods and systems for 4D geophysical surveys | |
CN113640808B (zh) | 浅水区海底电缆埋深探测方法及装置 | |
CN105259570A (zh) | 一种海上时移地震接收点位移校正方法 | |
CN111487678B (zh) | 一种确定高分辨率小多道地震最小偏移距和系统延迟的分析方法 | |
KR20120076952A (ko) | 굴절법 해양탄성파 탐사를 위한 obc 타입 스트리머 장치 | |
CN109212598B (zh) | 基于直达波反演的三维空间二次定位方法 | |
CN109116416B (zh) | 基于三维方式的倾斜海底检波器二次定位方法 | |
US20130046472A1 (en) | Method of determining the relative position of two detectors at the bottom of the sea | |
CN108693560B (zh) | 一种基于互相关道的散射波成像方法及系统 | |
CN112180452A (zh) | 基于探地雷达和三维速度谱的地下管线埋深估计方法 | |
US9423521B2 (en) | Method of determining the position of a detector disposed at the bottom of the sea | |
CN115980856B (zh) | 适用复杂海底界面的深拖式多道地震炮检精确定位方法 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
CN113703050B (zh) | 一种深海地震垂直缆二次定位方法 | |
EP0297852A2 (en) | Method for real time display of marine seismic survey data coverage | |
CN112904348B (zh) | 一种三维探测方法、装置、设备和存储介质 | |
CN109490962A (zh) | 一种浅地层剖面非零偏移距的消除方法 | |
CN113189646B (zh) | 一种去除拖拽式浅剖地层异常起伏的方法 | |
WO2008037965A1 (en) | Receiver orientation in an electromagnetic survey |
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 |