 点击上方“溪流之海洋人生”即可订阅哦浅地层剖面仪（Sub－BottomProfiler）是利用声波探测海底地形地貌和浅部地层结构、构造的海洋仪器，在海底地质调查、海底资源开发、海洋工程选址等应用中发挥着极为重要的作用。根据声波的产生方式，常见的浅地层剖面仪可分为压电陶瓷式、电磁式、电火花式和参量阵式这4种类型。压电陶瓷式是利用不同矿物晶体的压电效应产生声波，有固定频率和线性调频这2种，具有较高的分辨率但穿透能力较弱；电磁式是利用电磁效应使金属板产生位移从而产生声波，其分辨率与穿透能力适中；电火花式是通过电极间放电气化海水产生大能量声波，其穿透深度大但分辨率略有下降；参量阵式是利用小孔径基阵发射高频信号，经声学非线性效应得到低频、宽频带、无旁瓣高指向性的次频声波，可有效解决线性浅地层剖面仪因受低频工作机制限制而存在的横向分辨率低、混响与噪声抑制能力差等技术缺陷，适用于全海深探测，工作效率高，数据质量好，在浅地层剖面仪探测领域具有独特优势。
4 X0 H7 q# e1 N" s: `. P2 q 本文以参量阵浅地层剖面仪为背景，提出一种高精度、高采样率、宽动态范围和高运算能力的采集处理系统，该系统充分利用现场可编程逻辑门阵列（FPGA）在高速采集领域和数字信号处理器（DSP）在复杂算法处理领域中的优势，采用FPGA＋DSP架构形式的硬件平台和匹配滤波算法处理方式，经试验验证，系统运行稳定功耗低，在其他水声设备中也极具借鉴价值。 3 j6 A r9 w9 c" Y+ H
一、原理简介 : n6 L* f& X8 w4 M/ M
参量阵浅地层剖面仪是以声学非线性理论为基础，利用高频高增益的小孔径基阵发射含有ω1和ω2这2个相近频率成分的宽带原频信号，在介质的非线性作用下，声场中将产生低频无旁瓣高指向性的宽带差频信号（|ω1－ω2|）及高频超声和频信号（ω1＋ω2）等。根据声学经验公式可知，声吸收系数与频率的平方成正比，因此在声波传播过程中，频率较高的和频信号将快速衰减，频率较低的差频信号将持续作用。海底底质结构以层状分布，当声波向下传播时，一部分在介质分层处发生反射，另一部分经过透射后继续向下传播，在下一介质分层处同样进行反射和透射，依次循环递进，其示意图如图1所示。  . M3 d v- u3 f* W/ D# S
图1 参量浅工作原理示意图
, L9 f0 m* O6 s% d2 _9 k 假设水体是第一层介质，其密度和声波传播速度分别为ρi和vi；下一层界面为第二种介质，其密度和声波传播速度分别为ρ2和v2。声波的反射强度与反射系数R密切相关，数学表达式为： 2 a0 T/ G& ]# A% b% N: v% D
R＝（ρ2/V2－ρ1V1）/（ρ2V2＋ρ1V1），⑴
2 ]$ N) b9 E# I, n2 l) S0 P 因此当差频信号穿过海底沉积层时，由于不同介质层间的物理差异，会在界面处反射产生强弱不同的回波信号。设备接收机接收回波信号后经放大滤波，再进行特定的算法处理，结果输出到显示控制端，即可绘制出一串浓淡不一的反映地层结构和构造的像素点。 & b' A' X- Y" J
二、系统设计 . Q2 M2 z" a$ o
⒈系统架构设计
1 _4 V' F+ b# g' r# V4 i 参量阵浅地层面仪主要由基阵、功放单元、采集处理单元及显控单元组成。采集处理单元生成发射信号输出至功放单元，功放单元进行功率放大处理后作用于基阵，基阵将电信号转换为声信号发射出去；然后基阵再将反射回来的回波信号转换为电信号输送至采集处理单元，通过一系列的信号调理后进行采集和数据处理；数据处理过程中，频率较高的原频信号用于精确测量水深，频率较低的差频信号，因其具有很强的沉积层穿透能力而用于探测海底浅部地层结构信息；所有的处理结果都输出至显控单元，进行实时显示。根据上述采集处理单元的功能需求，综合考虑目前各种主流处理芯片在信号采集处理及控制等领域的优劣势，系统硬件平台采用FPGA＋DSP架构形式来实现，其框图如图2所示。 
+ v, A1 W& ^6 ]3 A% R3 p% Q( q. k 图2 采集处理系统硬件框图
( i9 E5 h* g+ `/ _3 t FPGA单元主要实现信号增益控制，高速AD信号采集、滤波处理，并将处理后的数据经高速数据总线（uPP，通用并行端口）传输至DSP单元；DSP单元是系统的核心处理单元，实现发射信号产生，地层部面和海底深度信息解析等功能。 6 W0 i7 V& L. H1 n- V$ t! S2 z
⒉回波信号预处理
, d# \4 E8 q& |4 U8 m5 ^* r) e0 ] 在进行信号采集之前，需要对基阵接收到的原频回波信号和差频回波信号进行低噪声放大和滤波处理，以得到高信噪比的模拟信号。信号预处理的流程如图3所示。  2 m& L' h1 `& @9 l: j9 x- J
图3 回波信号预处理流程图 # M; y: N) Y8 s2 |
根据电路噪声理论，放大电路自身的噪声主要由芯片供电电源噪声、电流噪声以及接收回路的电阻热噪声等噪声源引起。因此在低噪声放大电路设计中，前级接收回路中尽可能不串接电阻，第一级放大电路尽量选择电压噪声和电流噪声较小的器件。本系统采用低噪声JFET分立器件设计差分放大电路作为第一级放大电路。低噪声JFET的最大优点是其电流噪声很低，其电流噪声密度一般为fA/VHz级别，因此使用在传感器内阻于kQ级别的放大电路中，由电流噪声和传感器等效内阻形成的噪声几乎可以忽略不计。 " |7 M4 H" H0 h& [, A( p: u \
本系统原频信号接收通道的中心频率为100kHz，带宽为30kHz，即85～115kHz，设计需实现带外倍频程40dB的衰减。由于上限截止频率与下限截止频率比不大于2，该滤波器属于窄带滤波器，根据中心频率和带宽计算出滤波器Q值为3.3，且考虑到带内纹波及带外衰减要求，使用两级四阶MFB拓扑结构的巴特沃斯带通滤波电路，其仿真结果如图4所示。  - } X0 q) u9 }4 [1 @5 N
图4 MFB带通滤波电路仿真 ! \3 K" U7 Q" U6 `
差频信号接收通道的带宽为2～16kHz，其上限截止频率和下限截止频率分别为2kHz和16kHz，设计需实现带外倍频程40dB的衰减。由于上下限截止频率比为8，因此该滤波器属于宽带滤波器。宽带滤波器的低通和高通滤波器需分别设计，然后通过将低通和高通滤波器级相联实现带通滤波器。考虑到带内纹波及带外衰减要求，使用Sally－ken电路拓扑结构的巴特沃斯带通滤波电路，其仿真结果如图5所示。 
, j6 O- C/ ?' t8 Q$ T 图5 Sally－ken带通滤波电路仿真结果
! P9 k& O: ~; z. ]. @# S, ? ⒊高速数据采集传输
( P* e1 O* c% h: | 本系统设计中，选用AD7760实现AD采集功能。AD7760是一款高性能，24位Σ－△型模数转换器，融合了宽输人带宽（>110dB）、高速采样与Z－△转换技术的优势，非常适合高速数据采集应用。AD7760内部的Σ－△调制器以极高的频率对输入的模拟信号进行采样，并对2个抽样之间的差值进行量化编码。调制器内部积分器对于输入信号表现为低通滤波器，对于量化噪声则表现为高通滤波器，这样可以将大部分噪声能量移出信号频带之外，调制后的数字滤波器可以消除较大的带外量化噪声。
5 d5 u: f) `' W" L0 Z M AD7760内部使用3个串联FIR滤波器，第一个滤波器以20MHz最大频率接收调制器数据，然后进行1/4抽取，以5MHz频率输出数据；第二个滤波器提供2×至32×的可选抽取率，也可以将其完全旁路；第三个滤波器提供2×的固定抽取率，它是用户可编程的，此滤波器也可以旁路。根据AD7760芯片数据手册及硬件配置方式，数据输出速率与晶振频率有如下关系： ; m. k* f6 c5 U7 ]* ~
ODR＝MCLK/(2·F1·F2·F3)，⑵
+ t+ H( u ], d: F 式中：ODR为数据输出速率；MCLK为输入晶振频率；F1为滤波器1抽取率；F2为滤波器2抽取率；F2为滤波器3抽取率。滤波器旁路设置时，抽取率取值为1。
9 d4 i( i2 R8 _' h6 o. a/ j+ u: @$ s 根据设计需求，原频信号采样率不低于600kHz，因此选用38.4MHz时钟进行1/32抽取，输出数据频率为600kHz；差频信号采样率不低于96kHz，因此选用24.576MHz时钟进行1/128抽取，输出数据频率为96kHz，AD采样频率设计如表1所示。 % ]/ {4 i' g, W. z {" V
表1 AD采样频率设计 
* W! \$ ~6 V( O- K FPGA单元获取AD采样数据后通过高速数据总线（uPP接口）将数据传输至DSP单元。设计中，综合考虑数据更新速度及FPGA芯片存储容量，将采集到的数据以2K个点划分为一个单元，使用乒乓结构进行缓存（见图6），采用流水线形式通过uPP总线传输至DSP单元中。  $ e1 y) a8 d0 z" e' O3 n* k; V! A
图6 乒乓缓存设计 ! U8 v' i4 }& x3 G( A9 A
原频和差频数据缓存写入及读写时序关系如图7所示。其中，T1为原频数据缓存写满时间，T2为2K缓存空间读空时间，T3为差频数据缓存写满时间。  : }/ Z$ N# ~4 }/ I. l
图7 原频及差频数据缓存读写示意图 0 ^; J) f O. T& r: G
由于原频信号（采样率600kHz）和差频信号（采样率96kHz）同时采样，因此需要使用分时复用的方式来完成uPP通道数据传输。原频RAM写人速度是差频RAM写入速度的6.25倍，分时复用uPP通道时，缓存RAM的读出速度（2T2）需远大于写人速度（T1）才能保证数据不会发生覆盖。
; a% H: }+ x3 ` uPP时钟设计为48MHz，则原频通道缓存RAM写满耗时（T1）约为3.4ms；读出耗时（T2）约为42us；差频通道缓存RAM写满耗时（T3）约为21ms；读出耗时（T2）约为42us；T1远大约2T2。满足设计需求，实现了不同频不同源信号的同步采集传输。
3 k1 ^6 t5 l$ K: v ⒋系统软件设计 . ^# }! n/ N/ G. S: [( U
系统处理软件主要实现发射信号产生、地层剖面和海底深度信息解析、系统控制及通信等功能。
' M! T+ R# t: f4 G& m- A ⑴发射信号产生 ; b2 M) s7 j7 m/ x
本系统设计中，可以由显控单元设置发射信号形式和频率等信息，采集处理单元接收到显控单元的配置信息后，生成相应的发射信号并输出至功放单元。 9 U8 }- ^+ _+ n O @; u2 {
发射信号形式使用调幅信号，即将单频（CW）或线性调频（LFM）等信号调制到原频载波信号。穿透深度和探测分辨率是浅地层剖面仪的2个重要指标，根据声学理论，单频信号的脉冲宽度越宽，穿透能力越强，但探测分辨率越低；线性调频信号的探测分辨率主要是由信号带宽决定，与脉冲宽度无关，因此可通过增加脉冲宽度的方式来增加探测深度，而不影响探测分辨率。
1 T' ]( B3 L8 _, n- d6 H 本系统中使用LFM信号实现发射信号生成。f1为起始频率；△f为频率变化量；fs为采样频率；N为信号长度。信号表达式为：  - T3 P* Y* j, @, ^ H$ O2 Z
原频信号的频率为fo，信号表达式为：
) M w& E+ z3 O5 ]$ r xo(n)＝cos(2πfon/fs)，⑷
- [# i8 c3 E5 K2 J8 E 则发射信号的形式为： 5 K6 v; |% P. q( F" n" ~
X2(n)＝Xo(n)×x1(n)
& }& m& f$ V$ r* C0 R, Z7 t 本系统设计中，原频信号频率为100kHz，差频信号频率范围为2～16kHz可选。 ; U8 g% t& G" E! n( x8 D8 o
图8给出了一组发射信号。其起始频率为1kHz；采样频率为1000kHz；频率变化量为2.5kHz。该信号经声场非线性效应，即可产生2～7kHz的差频信号。  / T7 O; L7 y' ]2 f [
图8 发射信号示意图 + u5 k2 ?2 ^+ I
⑵地层剖面处理
( v1 p: B/ v1 K% a9 ~& ^ 差频信号的低频特性决定了其具有很强的沉积层穿透能力，差频信号穿透地层剖面后经反射和散射作用产生回波信号，采集回波信号并进行匹配滤波处理，即可得到海底浅部地层结构信息，匹配滤波处理方法如图9所示。  5 w9 k( U* S$ i6 r4 |; y" M
图9 地层剖面处理框图 ) J; n. W1 X% B: Y/ e" a
其中，参考信号为标准水听器在发射换能器下接收到的直达波，y(n)为含有地层剖面信息的结果数据。
( K* e/ T- {0 c0 h9 r0 a" Y3 W) j- F ⑶海底深度测量
- o; J! l+ A5 }- V 原频信号频率较高，其分辨率也较高，因此设计中使用原频信号来精确测量海底深度信息。对原频回波信号滑动积分，进行能量检测，找到回波能量最大值位置，即为海底深度的位置。 7 Y* c1 Y- t* r$ }$ D
⑷系统软件设计流程 & i& ?/ U$ f; M7 o W6 A; g
采集处理单元接收到启动命令后，根据显控单元设置的周期、配置定时器参数、启动定时器、输出发射同步信号、并启动信号采集、执行测量算法流程、将处理结果数据、姿态数据打包发送至显控单元，工作流程如图10所示。  & Q" x6 s* O1 t6 A& z
图10 系统工作流程
0 f' u; ?. {- j8 q- i 三、试验结果及分析 : T, t0 v! r7 [2 U
本采集处理系统已成功应用于某型参量阵浅地层剖面仪中，并于2022年11月在宁波北仑海域开展了海上试验。试验中，若配置原频信号为100kHz的单频信号，差频信号为5～13kHz的LFM信号，其探测效果如图11所示。  7 [7 s" t) G0 `+ w" t/ H3 [
图11 穿透深度和垂直分辨率测量结果
) `: |/ s: J0 L0 w& q 可知，该参量阵浅地层剖面仪运行稳定可靠，探测效果理想，海底界面及地层分层清晰，成像连续一致，设备在该测试海域最大穿透深度可达39.29m，设备垂向分辨率可达10cm。 2 }# [: \5 ]2 x) i' W& ^9 a; B
四、结语
$ Z9 q1 F" C. J; C/ x 本文以参量阵浅地层剖面仪为应用背景，充分利用FPGA、DSP及高速AD在采集处理领域的优势，提出了一种高精度、高采样率、宽动态范围和高运算能力的采集处理系统，并已成功应用于XX型设备中。经试验验证，系统测量精度高，性能稳定可靠；此外系统的硬件架构形式、信号预处理方式、不同频同源信号的同步采集等技术均对其他水声设备研制有极好的借鉴意义。
" }$ S: G1 L, ~, n$ Q 1 % i! k2 e w. C( z7 T# x2 r, u2 T
END 7 ~) q" }6 e4 ]/ Y
1
+ f: @& @/ W$ [8 ?% G7 ] 【作者简介】文/冯宏 郭冉 程何小，来自中国船舶集团有限公司第七一五研究所。第一作者冯宏，女，1987年出生，硕士，高级工程师，研究方向为信号处理机设计。本文受基金项目赞助，国家重点研发计划项目（2022YFC2803600）。文章来自《舰船科学技术》（2024年第21期），用于学习与交流，参考文献略，版权归作者及出版社共同拥有，转载也请备注“由”溪流之海洋人生“微信公众平台编辑与整理。  fill=%23FFFFFF%3E%3Crect x=249 y=126 width=1 height=1%3E%3C/rect%3E%3C/g%3E%3C/g%3E%3C/svg%3E) 9 q0 K2 R. a# _" z1 V, y8 K B
相关阅读推荐 海洋论坛▏参量阵浅地层剖面探测技术在海底管道检测中的应用海洋论坛▏参量阵浅地层剖面的处理方法及在识别浅层天然气水合物中的应用海洋论坛▏参量浅剖现状及其在海洋工程中的应用海洋论坛▏SES-2000标准浅地层剖面仪在埋深障碍物探测中的应用行业资讯▏南方海洋推出一体式和分体式两款国产参量阵型浅地层剖面仪知识窗▏参量阵技术在水声工程中的应用及其相关问题研究
; _0 }) h% g: w I, d$ h 公众号 ' L: H* \; J0 L- r1 h
溪流之海洋人生 0 Y' y$ `9 `6 g0 u
微信号▏xiliu92899 4 Q3 A% U( ^5 h& E- g1 `
用专业精神创造价值 & G: h& h& B0 i {- ]: Y5 O
用人文关怀引发共鸣 您的关注就是我们前行的动力 1 U6 A4 p$ S& m/ O% ?
投稿邮箱▏191291624@qq.com
6 i0 X! X: Y! u$ d: H7 i( \! F4 g; R% D" X5 C U
" I% ?9 \. ]# h; ~. e% a, w/ s' m1 {, I' l1 h) ?0 I1 V
# A, V* B A1 p& G; C8 d: Y/ U | 
