利用层析术进行海底地质断层勘探

2020-12-15杨士莪

哈尔滨工程大学学报 2020年10期

杨士莪

(1.哈尔滨工程大学水声技术重点实验室,黑龙江哈尔滨 150001; 2.海洋信息获取与安全工业和信息化部重点实验室(哈尔滨工程大学),黑龙江哈尔滨 150001; 3.哈尔滨工程大学水声工程学院,黑龙江哈尔滨 150001)

实际海底介质中的断层多为不平整的连续曲面[1]，通常海洋学方法多使用沿测线测量方法，效率不高；若需要了解一定海区范围内海底断面情况，需要设计往返多次的测量航线，既费时又耗资不菲，有必要考虑有效的改进方法。若采用水声学中的层析术进行勘探，则可以一次性地完成一定海区面积内各点处海底断面的深度和形状探测，不仅大大提高了探测效率，还可节约所需海试费用。

1 原理方法介绍

利用坐底声呐阵,按照层析术方法进行海底断层测量的原则方法如图1所示，在预定海域周边不同点处布设矢量坐底声呐阵,并在海域范围内不同位置点投放爆炸信号。设海深为H、各坐底声呐位置Si=(xi,yi,H)、各爆炸点位置Bj=(xj,yj,zj)等均为已知。海水中不同点处声速可利用其平均值及带有待定系数的各阶经验正交函数之和表示。

注：坐底声呐,爆炸点。图1 声源与坐底声呐位置示意Fig.1 Sketch of the positions of the sound sources and bottom sonars

若海底介质中断层曲面的函数表示为z=Γ(x,y)，海底介质密度与纵波波速和横波波速均可设定为待求的坐标函数。根据地震观测所得海底介质中断层的常见形式，为适当减少利用层析术测定断层曲面形状进行最优化计算时的待定参数数目，可建议采用以下2种可能形式之一:

1)将探测海域中断层设想为若干块不同断面，每块均具有空间平面形式：akx+bky+zk=0；为分割各断层区间，可再利用若干设定的空间点(xl,yl,zl)，以每个空间平面的四顶角坐标点划定该断层所在区域。这种断层一般位于海底介质中相对较深位置。

利用爆炸声信号进行实验时，此时各接收点一般至少可接收到:1)自爆炸点传来的直达脉冲；2)经一次海面反射后到达海底的脉冲；3)经海底反射后又经一次海面反射最终到达接收点的脉冲、以及经海底界面折射，并再经海底介质中断层界面反射后，最终到达海底界面的纵波和横波脉冲,共计5个信号。为计算简便起见,具体分析时可采用射线声学方法。设自爆炸点发出的声线为L0(α,β),其中(α,β)为表征声线射出方向的空间角的量,则声程方程的完全解可写为[2]：

V(L,x,y,z,α,β)=0

(1)

且有声线的空间方程式为：

(2)

式中：xj为爆炸点坐标；xp为自爆炸点出发的声线到达海底或海面点处的坐标。对于自爆炸点出发；直达接收点xi的声线，可取xp=xi，其中xi为接收点坐标值，由式(1)与式(2) 3个方程式联立求解,即可得出相应的声程数值、自爆炸点出发到达接收点的声线的出射角(L1,α0,β0)。

而对于经一次海面反射后到达接收点的声线，则首先可由爆炸点坐标xj与设想的海面处声反射点位置xp=(xp,yp,0)，列出相应的声程方程式与声线方程式：

(3)

(4)

此时共计6个方程式，6个未知量Lp、L1、xp、yp、αp、βp，即可联立求解得到各未知量数值，而声波传播时间，则可由声程Lp+L1求得。

同理，对自爆炸点经海底反射后又经一次海面反射最终到达接收点的脉冲，可假设声线在海底的反射点为xp=(xp,yp,H)，而声线射出方向为(α1,β1)，当该声线在海底反射后，到达海面坐标点xq=(xq,yq,0)处又经海面反射，最终到达接收点xi。此时可写出各段的声程方程和声线方程式：

(5)

(6)

(7)

以上共计9个方程式，未知数L1、L2、L3、xp、yp、xq、yq、α1、β1也一共9个，不难联立求解各未知数，而总传播时间可由3段的声程和求得。式(5)～(7)中α′、β′、α″、β″系为书写简便，仿前节用以分别代表对应声程的坐标偏导数。

对经过海底界面折射或进入海底介质的声波,由于纵波和横波的折射系数不同,因而需要分别讨论。设海底介质密度和拉米参数分别为ρ2、λ、γ，各参数均可为深度的函数，为简单计，不妨设为随深度线形变化。首先考虑纵波的传播，虽然当介质随深度变化时纵波和横波将会产生耦合，但只要耦合系数不大,并不会引起解算上的困难。另外,由于对纵波的分析方法与对横波的分析方法完全相同,下文中为书写简便计，将采用肩标“单撇”来表示对纵波分析的各量,而采用肩标“双撇”表示对横波分析的各量,不再分别重复书写。

已知海底介质中断层曲面方程为z=Γ(x,y),并设其在断层上下的介质密度与Lame参数分别表示为：

(8)

(9)

(10)

仿照前述方法，此时按声程方程与声线方程式并取二次声源位于海底界面折射点，可同样获得在上层海底介质中纵波的声程方程完全解V′(L′,x,y,z,α′,β′)，其中α′,β′则依折射定理，有：

(11)

由于考虑到希望在海底介质中断面上反射的纵波，能直接到达接收点，因而此时爆炸脉冲入射海底界面上的点，将不同于经一次海底、海面反射到达接收点的爆炸脉中的反射点，而将是xr=(xr,yr,H)。故有：

(12)

若此纵波射到海底介质断层上(x′t,y′t,z′t)点处，则该点处法线方向为：

(13)

(14)

所以

(15)

将各声程梯度项转化为空间方向角α′1、β′1,自断层上反射点发出的声线出射角即为α′1、β′1,根据此反射声应直接到达接收点xi条件,由海底介质中的声程方程和声线方程,以及xi点坐标值,可直接解出自断层反射点至接收点的声程,及此声线在断层面上反射点处的出射角数值:

(16)

至此问题即已完全解决。自爆炸点经海底介质中断层反射到达接收点的传播时间，可由从爆炸点至海面，由海面至断层反射点及自断层反射点至接收点的声程之和计算得到。

对于在海底界面折射并经海底介质中断层反射后到达接收点的横波信号，因为计算方法完全类似，在此即不重复。若使用4个坐底声呐，施行5次爆炸，各坐底声呐可接收到沿5条不同路径由爆炸点到接收点的爆炸脉冲，则共有100个信号，可列出100个方程式，此时若取描述介质空间参数变化的30个参数，描述海底介质断面位置的20个参数，则将有共50个待定系数，不难借助最优化方法，求得各待定系数的最可几值，完成借助层析术对断层位置的勘探结果。