浅水声速剖面的反演方法与实验验证.docxVIP

浅水声速剖面的反演方法与实验验证 浅水声速剖面的反演方法与实验验证 2000年5月第18卷第2期西北工业大学学报 Journal of N o rthw estern Po lytechnical U niversity M ay 2000V o l . 18N o. 2 浅水声速剖面的反演方法与实验验证 沈远海, 马远良, 屠庆平, 姜小权, 朱治富 (西北工业大学航海工程学院, 陕西西安　710072) 摘　要:提出一种在浅水域进行声速剖面反演的方法。通过将某地浅水声速剖面用经验正交函数 (EO F ) 进行表示, 利用一个很短的垂直线列阵接收声源信号, 估计信号的多途到达角结构(DOA ) , 以到达角结构的误差为代价函数, 对实测的声速剖面进行了反演。实验结果表明:所提出的方法是可行的。 关　键　词:反演, 声速剖面, 经验正交函数(EO F ) , 中图分类号:O 422. 2　　　文献标识码:A :2(0220212204　　声场环境反演(inversi on ) (acou stic tom ograp , 。通常, 对声置外, 还包括海洋中的声速剖面(它反映了海水中温度, 盐度等特性的变化) , 海底底质参数和海深等参量。声速剖面的反演问题是声场反演中一个具有代表性的问题[1～4]。利用声场来反演声速剖面首先是M unk 和W un sch 提出来的。在深海环境下, 他们考察了在射线理论下利用相对到达时延TDOA 等参量来定位, 而后发展了匹配场的方法来搜索声速剖面参量, 利用线性处理器, 采用网格法对扫描空间进行搜索。为了达到相当的精度, 需要尽量减少扫描间隔, 这无疑会增加计算量。本文引入了一种简便的高分辨的代价函数, 利用具有很好全局收敛性的模拟退火方法, 在国内某湖泊的实验所测得的数据, 对声速剖面来进行反演。对于一般的低频问题, 简正波模型具有计算量小、简便、精确的特点, 这是低频[4]采用简正波的原因。而对于本文所处理的高频问题, 确定声场采用射线模式则更加适合。另外, 本文对匹配准则加以修改, 利用接收点处垂直平面内的多途到达角(DOA ) 结构进行匹配处理。 声速剖面的EO F 表示[3, 6] 设对声速剖面做M 次采样C m =C (x m , y m ) , m =1, …, M , 每个声速剖面具有N 个采样点, 定义声 场协方差矩阵为 R c 作特征分解, 有 声场中任意一点的声速剖面可以表示为 C (x , y , z ) =C 0(x , y , z ) + ∑Α(x , y ) F (z ) (3) 式中, Αn 是坐标(x , y ) 的函数, 特征向量F n 是深度z 的函数, C 0(x , y , z ) 为样本的平均声速剖面。从国内某湖泊和南海的大量实测数据进行了验证, 可以认为在误差允许的范围内, 浅水的声速剖面可用EO F 函数近似地表示出来。这样, 不同的声速剖面和不同的EO F 函数的系数存在一一对应的关系, 反演不同的声速剖面最终变为搜索最优的系数。1. 2　信号特征参量估计[4] 基于射线理论模型的特征匹配方法, 充分利用了水声场中接受点的时间、空间信息, 使利用小孔径基阵完成反演成为可能, 这一切建立在能够有效精确分辨各入射本征声线的基础上。因此, 多途入射信 收稿日期:1998-06-01基金项目:国家自然科学基金资助 作者简介:沈远海(1974-) , 男, 湖北省人, 西北工业大学博士, 现在船舶系统工程部工作, 从事匹配场处理技术的研究。 第2期　　　　　　　　　　　沈远海等:浅水声速剖面的反演方法与实验验证?213? 号特征参量的高分辨估计是反演方法的关键。 对于同一信号源的不同本征声线而言, 这些入射信号都是相关的。因此, 所采用的DOA 估计算法是能够处理相关密集分布入射信源的高分辨DOA 估计算法。利用P illai 提出的不减少阵列孔径损失的向前向后平滑技术[4], 可以有效分辨入射信源。1. 3　模拟退火全局寻优[2, 5] 常规模拟退火方法包括几个步骤:(1) 构造系统能量函数 该能量函数应能够正确反映参量空间中不同参量对系统的贡献, 并对最希望反演的参量有较高的敏感度, 当所有参数都估计正确时, 能量函数应为最小值。采用函数为 式中, Υm 为[0, 1]内均匀分布随机数, T 为退火温度。 (4) 给定系统的退火进程, 即系统的降温策略 采用[5]中的逆线性的冷却进程 (j +