CN106959152A - 一种测量海水中声波速度的方法 - Google Patents

Abstract

本发明公开了一种测量海水中声波速度的方法，其利用基于海洋水平缆或斜缆采集得到的海洋地震勘探数据直接计算得出。通过先构建海水中一次波和相应的海平面虚反射之间关系的平面波模型；然后，提取若干道地震数据，对每一道海洋地震数据，分别按公式计算出不同方向来的一次波与虚反射的叠加能量E，依据绘制出的E随速度V_W的变化曲线的极大值点所对应的速度确定出海水中声波速度V_w。本发明直接利用海洋地震勘探的数据，不需要额外的观测，即可获得真实准确的海水中声波速度，测量效率高，省时、省力、省费用；并且，由于所获得的海水中声波速度就是地震数据采集时的速度，因而，本发明的方法更适合于地震数据的处理中使用。

Description

一种测量海水中声波速度的方法

技术领域

本发明涉及一种测量海水中声波传播速度的方法，具体涉及一种利用海洋地震勘探数据获取海水中声波传播速度的方法，属于海洋勘探技术领域。

背景技术

海洋地震勘探中，一般采用拖缆进行地震数据采集，拖缆在海水面下一定的深度上，沿着拖缆以一定的间隔布设检波器接收地震波场，若所有检波器的深度是相等的，则为水平缆采集；若各检波器的深度不同，则为斜缆采集，为了避开海浪等对地震波场的影响，检波器的深度一般不小于5m。

由于流体中只能传播纵波，不能传播波横波，因此地震波在海水中是以声波形式传播的。海水中地震波的速度就是海水中声波的速度。

地震波可以看成是不同方向传播的平面波。海洋地震勘探中，可以把地震波分为一次波和虚反射，其中一次波是指从电缆下方传播到接收点的地震波，虚反射指的是地震波传播到海水水面被反射向下到达接收点的地震波。

海水中声波速度是一个非常重要的物理量，在海洋勘探中有着重要的用途，其大小受到温度、盐度和压强(水深)的影响，随着空间位置和时间不断变化。在海洋油气勘探中海水中声波速度是海洋地震资料处理的关键参数之一。

如何精确求取声波在海水中的传播速度是一个重要的研究课题。

早在1827年，法国数学家Sturm和瑞士物理学家Colladon在日内瓦湖进行过水中声速测量的试验，并得到了极为正确的淡水声速。在上世纪初，水中声速值的精确测量为水声的实际应用奠定了扎实的基础：人们在灯塔上安装好潜水钟和雾号，根据在船上测量钟声在水中传播和雾号声在空气中传播的时间间隔测量航船与灯塔间的距离。从量纲上考虑，将声在海水中传播的时间转换为距离时都必须是精确的声速值。不同时间在同一海域，回声测深仪性能会由于“午后效应”的影响而发生变化，所谓“午后效应”，即回声测距设备，每到下午性能就下降的现象。由于太阳辐射，深表层水温于下午升高，形成负梯度型声速分布，使声波向下折射，因此上午能够探测到的目标，午却因声波折射而处于声波阴影区之中，信号显著减弱，甚至完全消失。回声测深仪因此成为精确测量海水声速的重要推动原因，而这种需要对于声纳应用中获得精确的目标距离数据显得尤为突出。

现有技术中，海水中的声速测量方法一般采取两种方式：一种是直接测量法，另一种是间接测量法。其中，直接测量法一般使用的设备为“声速仪”，即通常利用收发换能器在固定距离内测量声速，同时以压力传感器及温度补偿装置测量水深。根据获取声速的方法不同，通常情况下又分为环鸣法、驻波干涉法、脉冲叠加法以及相位法等。

间接测量法则是通过水文仪器，如温盐深仪等来测量海水的温度、盐度和深度，再利用这些环境测量与声速的经验公式，进而通过计算达到声速剖面。

以上所描述的目前国内的海水中声速的测量方法在很大程度上都依赖于对海洋环境要素的变化规律的掌握。但是，大自然的变化万千，随着时间的不同以及各种细微的变化，各项参数都各有不同。

对于海洋地震勘探来说，如何寻求一种能够测量地震数据采集时海水中声波速度方法非常重要，可以更加精确地提供地震资料处理所需要的还是声速信息，提高地震资料处理的精度和效果。

发明内容

本发明的目的是一种测量海水中声波速度的方法，其基于海洋水平缆或斜缆采集的地震数据直接计算获取，具有测量效率高、结果真实、准确，省时省力、成本低等特点。

本发明为实现上述目的所采用的技术方案是，一种测量海水中声波速度的方法，其特征在于，利用海洋地震勘探数据直接计算得出，具体步骤如下：

第一步，构建海水中一次波和相应的海平面虚反射之间关系的平面波模型；

第二步，提取若干道地震数据，对每一道海洋地震数据，分别按下式(3)计算出不同方向来的一次波与虚反射的叠加能量E：

上式(3)中：

s(t)为地震数据在t时刻的样点值；

h为地震波接收点距离海平面的垂直高度；

α为平面波与竖直方向的夹角；

τ为虚反射比一次波晚到达接收点时间；

V_w为变量，表示声波速度设定值；

第三步，将计算出的不同方向来的一次波与虚反射的叠加能量E数值作为纵坐标、以与其对应的声波速度数值作为横坐标，绘制成E随速度V_W的变化曲线，该变化曲线的极大值点所对应的声波速度设定值即为测量出的海水中声波速度V_w。

优选为，上述的测量海水中声波速度的方法，其特征在于，

上式(2)中：

l为虚反射等相位面比一次波多传播的距离，满足下式(1)：

l＝2hcosα (1)。

进一步优选，上述的测量海水中声波速度的方法，其特征在于，第二步由计算机执行，计算流程如下：

(1)、提取一道地震数据；

(2)、从地震数据道头中提取地震波接收点的深度h；

(3)、设定表示平面波方向的角度α范围、初始速度Vs、速度增量Δv和终止速度Ve；

(4)、用公式(3)计算出不同速度下的叠加能量E，并保存E和Vw数值，用于绘制E随速度V_W的变化曲线；

(5)、设定速度增加量V_w＝V_w+Δv；

(6)、判断V_w是否大于Ve，若是，则根据步骤(4)保存的E与V_w数值，以叠加能量E数值作为纵坐标、以与其对应的声波速度数值V_w作为横坐标，绘制E随速度V_W的变化曲线；

若否，则返回步骤(4)；

(7)、当E随速度V_W的变化曲线绘制完成后，找出变化曲线的极大值点，该极大值点所对应的速度V_w，即为海水中声波速度V_w。进一步优选，上述地震数据是基于海洋水平缆或斜缆采集得到的。

上述技术方案直接带来的技术效果是，通过构建海水中一次波和相应的海面虚反射之间的平面波模型；建立一次波与虚反射的时差关系(该时差关系与电缆深度、海水中声波速度相关)；构建一次波与虚反射的叠加能量与海水中声波速度的计算关系式；并基于“使叠加能量最强的速度就是真正的海水声波速度”的规则设计了一种通过速度扫描来确定海水中声波速度的方法。

上述技术方案通过直接利用海洋地震勘探的数据，不需要任何额外的观测，即可获得真实准确的海水中声波速度，并且，所获得的海水中声波速度就是地震数据采集时的速度。该速度更适合于地震数据的处理中使用，同时也能用于一般的海水速度测量目的。

即，本发明相对于现有技术，具有测量效率高，省时、省力、省费用等有益效果。

附图说明

图1为一次波和相应的海面虚反射之间关系的平面波模型；

图2为本发明的能量曲线计算流程图；

图3为海洋中声波速度模型；

图4为图3所速度模型的一道合成地震记录，接收点水深20m；

图5为叠加能量随速度变化的曲线，速度范围为1200m/s到2100m/s，每隔20m/s计算一个能量，极值点在1500m/s；

图6为一道实际的海洋地震数据，接收点深度为20m；

图7为图6所示数据的叠加能量随速度变化的曲线，其中，速度范围为1200m/s到资料2100m/s，每隔20m/s计算一个能量，极值点在1480m/s。

附图标记说明：

1、接收点，2、镜像点，3、一道地震记录，4、接收电缆，5、海平面，6、虚反射等相位面多传播的距离，7、虚反射，8、一次波，9、时差，10、平面波的等相位面，11、虚反射从海平面到镜像点的射线。

具体实施方式

下面结合附图，对本发明进行详细说明。

图1为海水中一次波和相应的海面虚反射之间关系的平面波模型。

如图1所示，图中5为海平面，4为接收电缆，1为接收电缆上的任意一个接收点，该接收点距海平面的距离为h(h不小于5m)；

接收点1相对于海平面的镜像点2，位于海平面之上，距海平面的垂直距离也为h；

8为两组来自不同方向的一次平面波，简称一次波，其中一组以α角入射到海平面的一次波等相位面10；

一次波在海平面被反射，形成虚反射7。

虚反射7从海平面向下传播到达接收点1的射线的长度与虚反射从海平面向上到镜像点的射线11的长度是相等的一次波先到达接收点，然后对应的虚反射到达接收点。

虚反射到达接收点的时间等同于一次波以海水声波速度继续传播到达镜像点2的时间。

图1中，数字6表示的是虚反射等相位面多传播的距离；一道地震3上面有四个地震子波波形，代表两组一次波及其对应的虚反射，数字9表示一次波与其对应的虚反射到达接收点的旅行时之差τ。

如图1所示，虚反射等相位面多传播的距离6(虚反射等相位面比一次波多传播的距离)为：

l＝2hcosα (1)

式(1)中，h为接收点在海面一下的深度，α为平面波与竖直方向的夹角。

根据如图1所示的海水中一次波和相应的海面虚反射之间关系的平面波模型，建立一次波与虚反射的时差关系为：

式中V_w为海水中声波的速度，τ为虚反射比一次晚到达接收点时间。

由于虚反射是一次波在海平面上的反射，海平面的反射系数为“-1”。因此，虚反射与对应的一次波相比，形状相同、极性相反。一次波与对应的虚反射相减，数值变为一次波的两倍(负负得正)，从能量的角度来看，相减后信号的能量(样点数值的平方和)则是一次波能量的四倍。

地震数据是由一次波与相应的虚反射组成的，根据一次波与虚反射的关系提出相减信号的能量计算公式。

考虑到虚反射与一次波存在时差τ，根据上式(1)和上式(2)把不同方向α来的一次波减去虚反射后平方求和，提出下式所示的相减信号叠加能量公式：

上式(3)中，s(t)为地震数据在t时刻的样点值。

在速度正确时，公式(3)中的平方项是一次波减去虚反射(负负得正)样点后数值、成为一次波样点数值两倍后的平方，虚反射加强了一次波能量，理想条件下是四倍。

在速度与海水速度不一致时计算的时间延迟τ与实际情况不符合，公式(3)中的平方项不再是一次波减去对应的虚反射，而是一次波减去一个随机的样点数值后的平方，能量得不到显著加强。

因此，根据公式(3)计算叠加能量E，在速度取值为海水中声波速度时E达到极大值，在速度偏大或偏小时E为较小的数值；

对一道海洋地震数据，用公式(3)计算不同速度下的叠加能量E，得到E随速度V_W的变化曲线，极值点对应的速度就是海水中声波速度。

图2为本发明的能量曲线计算流程图。如图2所示，能量曲线计算流程如下：

(1)、提取一道如图4所示的20m深接收点合成地震记录，如图3所示，其速度模型的合成地震道包含一次波和虚反射；

(2)、从地震数据道头中提取地震波接收点的深度h；

(3)、设定表示平面波方向的角度α范围、初始速度Vs、速度增量Δv和终止速度Ve；

(4)、用公式(3)计算出不同速度下的叠加能量E，并保存E和V_w数值，用于绘制E随速度V_w的变化曲线；

(5)、设定速度增加量V_w＝V_w+Δv；

(6)、判断V_w是否大于Ve，若是，则根据步骤(4)保存的E与V_w数值保留步骤(4)计算出的叠加能量E数值，并，以该叠加能量E数值作为纵坐标、以与其对应的声波速度数值V_w作为横坐标，绘制E随速度V_w的变化曲线；

若否，则返回步骤(4)；

(7)、当E随速度V_w的变化曲线绘制完成后，找出变化曲线的极大值点，该极大值点所对应的速度V_w，即为海水中声波速度V_w。图3为海洋中声波速度模型。

如图3所示，在1000米深的海水中，声波传播速度约为1500m/s，海水以下的地层中，声波的传播速度约2500m/s。

如图5所示，确定极大值点对应的速度为1500m/s，与图3所示模型中海水中声波速度(1500m/s)一致。

该数值与图3中海水中声波速度完全吻合。

证明，本发明的测量海水中声波速度的方法，所获得的声波速度的结果是正确的、可靠的。

下面结合实施例，对本发明作进一步的说明。

实施例1：

(1)取某海洋地震勘探数据一道，如图6所示；

(2)设定角度α范围为0°-90°、初始速度V_w＝1200m/s、速度增量Δv＝20m/s和终止速度V_end＝2100m/s；

(3)用公式(3)计算E；

(4)绘制出能量E随速度V_w变化曲线，如图6所示，确定能量E随速度Vw变化曲线中极大值点对应的速度为1480m/s。

该数值与中海油服天津分公司用声纳在该海域测量后提供的该海区的声波速度1490m/s，数值非常接近。

表明：本发明的测量海水中声波速度的方法，完全可以替代仪器的现场测量方法。从而，更加方便快捷，既省时省力，有可大幅提高测量效率、节省现场测量费用。

Claims (4)

1.一种测量海水中声波速度的方法，其特征在于，利用海洋地震勘探数据直接计算得出，具体步骤如下：

第一步，构建海水中一次波和相应的海平面虚反射之间关系的平面波模型；

第二步，提取若干道地震数据，对每一道海洋地震数据，分别按下式(3)计算出不同方向来的一次波与虚反射的叠加能量E：

$E={\mathrm{\Σ}}_{\mathrm{\α}}{\mathrm{\Σ}}_t{\left(s\right(t)-s(t-\mathrm{\τ}\left)\right)}^2={\mathrm{\Σ}}_{\mathrm{\α}}{\mathrm{\Σ}}_t{\left(s\right(t)-s(t-\frac{2hcos\mathrm{\α}}{V_W}\left)\right)}^2---\left(3\right)$

上式(3)中：

s(t)为地震数据在t时刻的样点值；

h为地震波接收点距离海平面的垂直高度；

α为平面波与竖直方向的夹角；

τ为虚反射比一次波晚到达接收点时间；

V_w为变量，表示声波速度设定值；

第三步，将计算出的不同方向来的一次波与虚反射的叠加能量E数值作为纵坐标、以与其对应的声波速度数值作为横坐标，绘制成E随速度V_W的变化曲线，该变化曲线的极大值点所对应的声波速度设定值即为测量出的海水中声波速度V_w。

2.根据权利要求1所述的测量海水中声波速度的方法，其特征在于，

$\mathrm{\τ}=\frac{l}{V_W}=\frac{2hcos\mathrm{\α}}{V_W}---\left(2\right)$

上式(2)中：

l为虚反射等相位面比一次波多传播的距离，满足下式(1)：

l＝2hcosα (1)。

3.根据权利要求1所述的测量海水中声波速度的方法，其特征在于，第二步由计算机执行，计算流程如下：

(1)、提取一道地震数据；

(2)、从地震数据道头中提取地震波接收点的深度h；

(3)、设定表示平面波方向的角度α范围、初始速度Vs、速度增量Δv和终止速度Ve；

(4)、用公式(3)计算出不同速度下的叠加能量E，并保存E和Vw数值，用于绘制E随速度V_W的变化曲线；

(5)、设定速度增加量V_w＝V_w+Δv；

(6)、判断V_w是否大于Ve，若是，则根据步骤(4)保存的E与V_w数值，以叠加能量E数值作为纵坐标、以与其对应的声波速度数值V_w作为横坐标，绘制E随速度V_W的变化曲线；

若否，则返回步骤(4)；

(7)、当E随速度V_W的变化曲线绘制完成后，找出变化曲线的极大值点，该极大值点所对应的速度V_w，即为海水中声波速度V_w。

4.根据权利要求1至3任一所述的测量海水中声波速度的方法，其特征在于，所述地震数据是基于海洋水平缆或斜缆采集得到的。