5 P/ f* L0 q( z3 v 引言4 Y. E) d/ E1 }) K
地球上的油气资源约有50%储藏在海洋中，全球已探明的海洋石油储量约为540亿吨，其探明程度远低于陆地油气资源。因此，海洋油气资源具有巨大的开发潜力，是石油勘探未来的希望。
. _9 a h) q* [ N 本文的主要研究内容是海洋地震勘探双检拖缆系统中的几个关键技术，这些技术的研究成果主要用于指导双检拖缆系统的实际研制，除此以外也可将这些研究成果用于其他包含大规模传感器网络的系统中。
$ i; ]3 b* {6 K# h 一、海上拖缆地震勘探技术概述
J1 |" l" A$ p- P% z 地震勘探技术是海洋物探最实用的技术手段，按检波器所处位置的差异，主要分为海上拖缆地震勘探技术和海底电缆地震勘探技术。与海底电缆地震勘探技术相比，海上拖缆具有施工灵活、成本低、效率高的优点，是目前海洋地震勘探的主要方法。海上拖缆地震勘探技术用于油气勘探领域，在海洋工程、科学研究等方面也应用广泛。
3 l5 z& B d7 i$ B 海上拖缆地震勘探系统的结构示意图如图1.1所示。海洋地震勘探通常是以物探船的形式来进行实际勘探作业，物探船在海上利用人工震源激发地震波，地震波传播到海底，经海底地质反射，反射波被水下拖缆中的检波器阵列接收。拖缆采集系统将检波器得到的模拟信号进行量化得到数字信号；传输系统将采集系统得到的地震波数字信号通过电缆、光纤等技术上传到位于船上的上位机中，同时还需要对信号进行实时处理和存储以便于后续的数据分析。由于海底不同的地质结构会导致声波在其中的传播阻抗不同，根据接收到的地震波信号进行声波阻抗反演，即可获得海底的地质信息。 ( I4 }+ |- h' X8 d
' Q0 n5 x* O! K9 [, F+ {5 b
图1.1海上拖缆地震勘探系统结构图 3 u" K. V+ N7 ?0 {8 I
二、海上拖缆地震勘探技术发展5 O6 ]( M+ H& w3 i
常规海上拖缆技术的主要工作方式是通过一条水平缆中的水听器阵列采集地震波信号，这种方式下拖缆布置方便、成本低、效率高。但是，在实际作业中拖缆需要置于一定的沉放深度，由于海水的声阻抗比空气大，当海底反射回来的有效反射波信号向上经海水传播到空气中时，会在海水界面发生全反射。
1 E1 m5 R* B% q" z3 ]: P) v 这些经全反射反向传播的地震波信号同样会被拖缆中的水听器捕获接收，从而形成对有效反射波信号的干扰，这种海面反射波也称为鬼波。鬼波最大的危害是对地震波数据造成陷频效应:当拖缆沉放深度较浅时，地震波数据频带较宽，高频性能优异，但是低频部分则缺失严重；拖缆沉放深度较深时，低频成分相对丰富，但是频带变窄，这极大地制约了采集的地震资料质量。为解决上述问题，一系列用于海洋拖缆地震勘探的新技术得以发展并被运用。 . f; S% T* C- l. D
斜缆采集技术(又称为变深度缆采集)最早由C.Ray于1982年提出，指的是通过控制拖缆中检波器的沉放深度随检波器偏移距离的增大而增大，使拖缆倾斜布置在水中，如图1.2所示。
* l r% P% v/ L r$ ]" H3 z 由于检波器深度的不同，接受信号的频带也是不同的，对应的探测岩层深度也不同。因此，斜缆采集数据的频带宽，不同深度岩层均能获得较好的成像。同时，针对斜缆采集的地震资料，利用镜像偏移和联合反褶积技术，可以有效实现鬼波压制。
' z& o* r" b8 C: ~5 p; f/ l # @, ?: Q8 n/ V0 j
图1.2斜缆采集技术示意图
+ r7 X/ Y7 e a& x, T7 m. l1 j 上下双缆采集技术，是将两条平行拖缆布置在同一水平位置、不同深度上实现地震资料采集，如图1.3所示。 p* K, I n- ]7 s% f
一方面，上下双缆在不同的深度，相较于单缆采集，接受信号的频带更宽，在低频部分也可以获得丰富的信息；另一方面，针对上下双缆，可以通过稀疏反演技术、F-K域波场沿拓技术等方法实现鬼波压制。 3 d6 b* J4 b( B) W& Z% B0 X. O" N
% p+ ^! ]7 r, e8 Z6 J0 c* I9 ?- r2 P' a
图1.3上下双缆采集技术示意图 % s* q, C) y* q$ w
双检采集技术，指的是将压力检波器和垂直速度检波器同时布置在同一个检测点中，实现同一检测点上压力和垂直速度测量。对上行波，压力检波器和垂直速度检波器的极性是相同的；对下行波，也就是海面反射产生的鬼波，压力检波器和垂直速度检波器的极性是相反的。 ( [; U- a- f4 s, q" `8 \
据此，可以对压力检波器和垂直速度检波器得到的地震资料进行合成，实现有效的鬼波压制，提高分辨率。在双检拖缆采集技术中，除压力检波器外，还采用了垂直速度检波器来探测垂直方向的速度并用于双检数据合成。四分量拖缆采集技术可看成双检采集技术的延伸和发展。
/ w v+ P- U) m 四分量拖缆中使用MEMS三分量检波器代替水检拖缆的垂直速度检波器，从而获得横向空间上的数据。因此在双检采集去鬼波的基础上，还可以根据横向空间数据，利用重构算法来提高横向空间分辨率，大幅提高了地震资料的信噪比。
/ n4 B w7 J9 b5 i; Z4 v3 l9 L, A 3 ~7 W- K$ t6 r. i: ~ ?6 u% c
图1.4四分量拖缆示意图(水听器+三分量MEMS检波器) 0 {2 P# B3 a% M# \% U
在上面介绍的海上拖缆采集新技术中，斜缆采集技术需要利用外部的控制系统实现拖缆的倾斜布置，布置难度大，布置成本高；上下缆采集技术，需要的拖缆数更多，相应地采集成本也更高，且上下缆的相对位置对海底地震勘探分辨率影响较大，勘探效果易受拖缆漂移情况的影响。
& I4 D9 H: c% e5 X' @ 相较于上述两种技术，双检拖缆采集技术只是在常规单水检采集系统的基础上，将采集系统由单一的水检改为双检检波器，从拖缆的布置难度和采集成本的角度来考虑都相对较容易实现。除此以外，对双检数据合成去鬼波算法的研究也己经相对成熟。 1 v& O3 C8 \+ G- [9 O
三、海洋物探双检拖缆技术
* C% M! F, t' _! | 3.1双检拖缆技术概述 " B2 @3 A8 E) j/ y# u, Z" A, g$ @
双检采集技术最先是在20世纪80年代末被提出的，不过在当时的条件下只能用运动传感器来测量速度分量，实际使用效果并不理想。在2007年，Carlson等学者提出用万向架技术来固定垂直速度传感器，从而可以准确地测量出垂直方向上速度分量。至此，双检拖缆采集技术才正式进入了实用阶段。
, c" s& m% ]4 ]" `3 h 目前，市面上能看到应用双检拖缆技术最成熟的是挪威PGS公司推出的Geo Streamer双检拖缆采集系统，其中陆检部分使用的是基于MEMS技术的垂直速度检波器。双检拖缆系统示意图如图1.5所示，检波器组中H表示水检压力检波器(Hydrophone)，G表示陆检速度检波器(Geo phone)。 , L. C4 B6 q6 B) g3 y- w' @
区别于传统海上拖缆系统，双检拖缆中用集成水检和陆检的检波器组代替单独的水检检波器，两种检波器被固定在拖缆上同一检波点处。水检压力检波器利用压电效应记录检波点处由地震波引起的水声压力波场变化，灵敏度高；陆检速度检波器则是通过机电转换记录检波点处速度波场的变化，适应范围广。对两种检波器的输出信号利用一
& M8 j" t' x) ] 定的算法进行数据合成，就可得到消除鬼波后的有效反射波信号。
x! Z% k* p6 v; X1 P. |& b2 t
7 r. ]; A: N8 q. T: P0 V 图1.5双检拖缆地震勘探技术示意图
9 s4 ~" a& I8 |1 K: @" Z" d0 C* {9 A 3.2双检拖缆技术的特点 ( Q- b, p6 X9 T* |
1)消除鬼波
" e& O; f( y, U6 o$ t# i% E3 z 双检拖缆技术中利用集成水检和陆检的检波器组在同一时刻记录该检波点处的地震波场，由陆检速度检波器记录垂直方向上由地震波引起的质点运动速度，由水检检波器记录同一点由质点运动引起的水声压力变化。 1 M8 M9 X' ~' P7 c9 H% |
水检记录的水压分量是无方向的标量信息，只有极性的概念；陆检记录的垂直速度分量是矢量，速度方向与检波器极性对应。对于上行波来说水检和陆检记录的极性相同，而对于下行波水检和陆检记录的极性相反。上行波在两种检波器记录的极性全为正或者全为负；下行波水检记录极性为正时，陆检记录到的极性则为负。 $ J5 G( q2 ^1 \5 y: A
两种检波器在地震波场中所处位置相同，因而反射周期也完全一致。将这二者记录的地震波场相加则下行波得以压制，上行波得以加强；相减则会只得到下行波。综上所述，利用水检记录的水压分量与陆检记录的垂直速度分量在对于不同类型地震波信号表现出的特征差异，可以有效消除鬼波对于有效地震波的干扰，提高采集系统的信噪比，改善记录的地震波数据质量。
+ t7 q* q% v3 W 2)拓宽频带 6 g, W. {/ J) M$ k" s
由于双检拖缆中集成水检和陆检的检波器组，可以同时对地震波场中的水压分量和垂直速度分量进行记录，利用特定的数据处理算法可以实现对鬼波的压制，从而消除由鬼波造成的陷频效应。双检拖缆中水检水听器、陆检速度检波器和二者合并道数据的频谱曲线如图1.6所示。 $ F; H6 R- \. L; M! v1 C. a. [
通过分析此图可得到以下结论:陆检速度检波器记录的数据低频成分相对丰富，但是在中高频段振幅迅速下降；水检水听器记录的数据中高频端振幅较大，然而低频成分相对丢失严重；二者合并后刚好可以产生互补，水陆检合并数据低频成分丰富、频谱曲线相对光滑。综上所述，双检技术实现了提升低频成分、拓展频带宽度的目的。
- a3 i1 R8 N; z7 X5 U) r: W; N' X ; F% N1 ]/ p: o' L
图1.6双检拖缆检波器幅频响应曲线
& U0 ]- W. j9 ?9 A 3.3双检消除鬼波的原理 ) n1 N8 z1 K, P
双检拖缆采集技术中，水检记录的是地震波场引起的海水疏密程度的变化，陆检记录的是地震波场引起的检波点处质点在垂直方向上的运动。以弹簧的伸缩来模拟地震波在海水中的传播，考虑如图1.7所示的情况，假设纵波的传播方向由下向上。压缩图中位于初始状态的弹簧，弹簧变为压缩状态。 ( W6 n! T& @6 d. t* A
如图中压缩部分所示，在弹簧上会产生一个由下向上的压缩波场，同时弹簧末尾的质点在垂直方向上的运动同样是由下向上，此时纵波传播方向与质点运动方向相同。考虑相反的情况，拉伸图中位于初始状态的弹簧，弹簧变为拉伸状态。如图中膨胀部分所示，在弹簧上会产生一个由下向上的膨胀波场，同时弹簧末尾的质点在垂直方向上的运动变为由上向下，此时纵波传播方向与质点运动方向相反。 ; n2 F. o9 m; v5 f! R
综上所述，可得到以下结论:纵波传播方向与质点运动方向相同会产生压缩波场；纵波传播方向与质点运动方向相反会产生膨胀波场。
9 R* Z# h7 D3 P- ]6 S3 a4 J2 \! T B. \ ) `3 m4 u& x+ Q' c8 A6 g% J" n+ ~
图1.7纵波传播方向与质点运动方向的关系 9 \+ Q/ b* v/ z1 J
水检记录的是水压随着海水疏密程度的变化，当海水由疏变密开始压缩水听器，水听器记录的极性为负，海水由密变疏时水听器记录的极性为正。陆检记录质点在垂直方向上的运动，质点运动方向发生变化时陆检检波器记录的极性也会随之改变。质点运动方向为上时，陆检记录的极性为负；质点运动方向为下时，陆检记录的极性为正。由于海平面反射系数接近-1，海底有效反射波经海平面反射后，极性反转，同时由上行波转变为下行波。 ' `5 U, ^. ^$ `0 J
因此，如下图1.8所示，对于水检检波器来说，其记录的有效反射波和鬼波极性相反；对于陆检检波器来说，其记录的有效反射波和鬼波极性相同。即水、陆检波器接收到的海底有效反射波极性相同，而鬼波信号极性相反。 " \- O) _' f# T. i
3 _1 }7 _1 y" E% U# s, v: _9 @
图1.8双检数据合成去鬼波示意图
$ ]* \' R, ^% F4 Z) K# r 双检合成消除鬼波信号的基本原理，正是利用水检检波器和陆检检波器对于有效反射波和鬼波记录极性规则的差异，采用适当的双检数据合成算法合成数据就可得到加强后的有效反射波信号，实现消除鬼波的目的。 1 t4 K3 W u7 ?+ t( N
3.4双检数据合成的算法 ( I# u$ E' R+ {. y$ z
海洋物探双检拖缆技术中，水检部分采集到的压力波信号P可分为两部分:海底有效反射波pup及海面反射鬼波pdown，所以有公式(1.1)。再对公式(1.1)左右两边同时做频率-波数域的二维傅里叶变换可得到公式(1.2)。其中，w表示角频率；kx表示沿拖缆方向的波数；z表示拖缆的沉放深度（ 即检波器深度 ）。
4 n' k! r/ Y7 D* G) h 1 R/ L2 h: j6 w
（1.1）
& q* Y2 [( A' q3 w
0 H$ X4 v( i( _: f （1.2)
- G4 i1 q5 F* S 由于海面反射系数接近于-1，根据频率-波数域的二维声波方程，即可得二维波场沿深度轴方向的波场延拓公式 + h0 D+ w' p$ G; U4 O
, N5 p! s# A9 g% F7 z' Z; S
有效反射波pup及鬼波pdown可表示为式(1.3)所示。其中，kz表示不垂直方向的波数，P0(W，kx，0)为海平面处波场。在不考虑y 方向的情况下，kx、kz和波数k的关系为式(1.4)。由式(1.2)和(1.3)可以得到pup和p之间在频率-波数域的关系式，公式(1.5)。
& B" \0 H4 W8 A5 F
5 ~/ ]8 p2 X8 @: e6 e5 D! [ (1.3)
: M- f) o! M! u- Z p , J2 z" D" j9 Q5 m1 m! h+ d
(1.4)
2 x& g: H) i! u( ~% f# d 1 Z3 N; S" v' |$ g
(1.5) ) e$ n R s7 l* Z+ T9 O4 e/ p
根据牛顿第二运动定律，水检测得的压力波场p和陆检测得的速度波场vz有公式(1.6)的关系，即声波扰动时的动力学方程。式中，ρ表示海水的密度。对式(1.6)两边做频率-波数域的二维傅里叶变换，依据傅里叶变换的性质，得到式(1.7)。 * _$ x C4 Z# z" I) v
4 ~( t3 T5 ]* G1 [' u0 x （1.6)
4 _3 W: z+ Z* T+ n0 W; n! P
: s1 Q9 N. g1 N" R+ L (1.7)
E! ?- e+ F, @7 p' e J8 N9 v 将等式(1.3)代入等式(1.7)可得到频率-波数域速度波场vz与海底有效反射波pup的关系式，公式(1.8)。
9 V0 X, ^3 M7 y/ t/ e
$ ~) w0 i7 A; K( k" @: [1 O (1.8) , S- F- T7 C4 ~, c: A
联立等式(1.2)和(1.8)，可得到方程组(1.9)。 % g% g1 d* \$ [3 A2 y
9 t/ ]: A) B1 V4 b6 h (1.9) ! E# q, {8 M6 d: w& l
通过解析上面方程组就可以将有效反射波pup和海面反射鬼波pdown分离开来，达到消除鬼波的目的。因为这种算法是将双检数据变换到频率-波数域进行处理，所以通常称之为f-k域波场分离法。 & N4 [) m1 Z) _% w6 e
四、双检拖缆中的关键技术3 |! F f: k4 `+ n4 m
海洋物探双检拖缆不仅是海上拖缆系统的简单延续，也是海上地震勘探发展的新阶段。与传统的单检波器拖缆相比，双检拖缆因工作原理和系统结构的不同，在设计与实现过程仍需解决传统拖缆不曾遇到的各种挑战。双检拖缆的主体研究内容包括:
+ N5 G. m, |& V- [" H. ^2 c 多种不同类型检波器输出信号的调理及读出，多检波器对应的多倍数据长距离传输，多通道多检波器釆集过程的同步，水下拖缆数据的汇聚，上位机与拖缆系统的数据交互方式，整个系统的远距离供电，数据的存储与分析，双检
; ~8 ^8 H0 J) E0 N u+ O, H9 F 数据合成算法的模拟与实现，拖缆姿态的控制，上位机的实时控制与状态监测软件，震源激发子系统等等。 3 v2 i* y. Q$ Z/ R$ o( W! _7 }3 [
根据上文对于双检拖缆技术基本原理以及所需研究内容的分析，结合实验室现有的条件，明确本文拟解决的关键技术可分为以下部分: $ |( v% ?8 R6 s0 i9 u; R
4.1数据采集
P4 o7 T. s ?0 F0 o 双检拖缆中采集系统的性能决定了整个勘探过程的质量和精度。采集系统部分主要用来实现地震波信号采集的功能，包括地震检波器、信号调理电路和信号读出电路三部分。由于双检采集过程中需要采集地震波同一时刻同一位置的不同表征分量，例如压力波场、速度波场等，所以采集系统中需要采用不同类型的检波器作为基本的传感单元。 ( C: z2 B v; p7 X% a2 F3 D9 i
传统拖缆的采集系统中往往只需要采用水检压电式水听器来将水声压力波场转化为电信号，而双检拖缆中除此之外还应该包含用于检测速度波场变化的陆检检波器。在陆上石油勘探中，主要是采用动圈式速度检波器(如SM-4、20DX等)来进行转化地震波振动为电信号的工作。
+ H5 D+ I+ I( ~2 t: A/ b% P 近些年来，在国外的地震勘探市场上，基于MEMS技术制成的的加速度传感器也同样广泛运用于陆地石油勘探，MEMS加速度传感器具有大动态范围、高灵敏度、低噪声、低零漂等优点。
. E4 Z1 v3 ]" f; t) z' g3 F 因此本文研究的采集系统中检波器部分将包含压电式水听器、动圈式速度传感器和MEMS加速度传感器。
: v" K6 J8 {) c/ a 对于不同的传感器来说，其测量的地震波分量不同，工作条件和工作原理也都不尽相同，这造成了不同传感器输出信号特征是不同的，所以每种传感器都应设计相应的信号调理电路。地震勘探中通常要求采集系统的动态范围高达120dB，如何实现大动态范围的高精度波形数字化技术，也是数据采集部分研究的重点。
1 z' @! P' n/ w$ x 4.2数据传输
$ R5 n) b4 }9 @# y- U$ R 双检拖缆中传输系统的性能好坏决定着采集系统得到的地震波数据能否被正常传输至上位机，传输系统将分布在整条拖缆各处的采集电路连接为一个整体。
" p$ x% y& C3 z 双检拖缆相较于传统拖缆系统带来最直接的变化就是检波器数量成倍的增力口，传统拖缆系统中每个采集节点只对应于一个水检检波器，而本文研究的双检拖缆系统中每个采集节点至少包含一个水检检波器、一个陆检速度检波器和一个陆检加速度检波器。
- q7 p, T! ^" @# d, O/ l+ l7 N 若考虑到将双检拖缆延伸至四分量拖缆，每个采集节点需要一个水检检波器和一个三分量陆检检波器即总共四个地震波检波器。检波器数量增加的同时带来的是数据传输速率需求增大。参考国内外海上拖缆地震勘探技术支持厂商的相关产品，拖缆的有效采集长度基本维持在12km以上，道间距普遍为12.5m，最大采样率一般能达到0.5ms，同时采样精度为24bit。 + W4 F. y v$ r( C8 v3 M6 a$ l
以此为依据来讨论双检拖缆技术中传输系统性能的需求。双检拖缆系统中各采集电路的数据是通过上传至传输电路，再由传输电路通过级联的方式逐级向上传输，这样越靠近前级，传输电路需满足的数据传输率越会呈线性增大的趋势。若考虑实际工作过程中数据的实时传输，最前级传输电路(与船上系统直接相连)需满足的最小数据传输速率为(12000/12.5)X(1000/0.5)x24x3=138.24M6ps。
b G$ O3 T1 Q V( `. ]* X, G" i 以此为基础考虑数据传输协议层的损耗以及冗余设计，要求的数据传输率将达到250Mbps左右。综上所述，双检拖缆的数据传输部分主要研究内容为传输电路间高速数据的长距离可靠传输。
( l: V$ j4 ^& s: d) P" }1 A( I) _ 随着检波器数量和拖缆数量的增加，上传至船上系统的数据量也会随之倍增。如果仍使用传统拖缆技术中基于读出机箱中多块采集板卡实现的数据读出方法，所需的机箱数和连接线缆数量将是个很庞大的数字。如何简单有效得将汇聚到船上系统的数据上传至上位机也将是数据传输部分研究的重难点。 3 b* R" n' c; `3 @/ j
4.3系统同步 ; F4 G" D- T" r7 ^( H S
水下拖缆系统的同步是传统拖缆设计过程中需解决的关键性问题，对于双检拖缆来说，这个问题同样存在。相对于传统拖缆，双检拖缆采集系统中检波器数量更多，地震波采样得到的数据量更大，系统的同步难度也更高。双检拖缆水下有效采集段的长度普遍能达到几千米，在该采集段内布置了大量的地震波检波器，对于这样大规模的传感器网络来说，如何保证多个采集通道间的高精度同步一直是海洋地震勘探技术中的关键性问题。 ) B6 J' S. E! k6 ]: R" ]' @3 P
在双检拖缆系统中，系统的同步可以分为命令同步和时钟同步两部分。命令同步，即命令在传输过程中出现的传输延时问题，传统的拖缆设计常采用延时补偿的方式来解决，逐级增加命令的延迟时间可以保证命令信号在每一级的同步性。延迟时间测量比较精确的方法是利用命令回环技术。
; U7 N, e; x/ Z$ f5 `, Y, ] 对于命令同步技术的研究已经相对成熟，本文对此部分仅进行简单介绍不做深入讨论。对于时钟同步问题的解决，传统拖缆设计中常采用的方法是基于PLL锁相环的时钟分发技术和基于SerDes串并转换的时钟恢复技术。这两种技术都存在着技术线路复杂、功耗高等问题，对于新型双检拖缆的设计已不再适用。所以，时钟同步技术将会是双检拖缆系统同步的研究重点。
8 Y' m0 X: g6 H* i& ^; s) P$ Z; F& G 总结 b5 Q% z' d3 p, Y
本文的研究内容主要是双检拖缆系统中的几个关键技术，分别是数据采集、数据传输和系统同步三个方面。这些技术的研究成果主要用于指导双检拖缆系统的实际研制工作，以外也可将这些技术成果用于其他包含大规模传感器网络的系统中。
: S* X5 t* Y5 q
5 u6 P0 ^3 O/ \1 _* o8 k( _, i6 N( o! m# T! m& r3 P' H
| 
