基于水平阵浅海海水声速剖面匹配场反演方法.pdfVIP

基于水平阵的浅海海水声速剖面匹配场反演方法 李风华张仁和李整林 (中国科学院声学研究所声场声息国家重点实验室，北京100080) Sound ProfileInversionWithAHorizontalLine in Speed Array Shallow、№ter LI ZHANGRen—he，LI Feng-hua， Zheng—lin 100080 of J LaboratoryofAcoustics，InsdtmeofAcoustics，ChineseAcademySciences，Beijingchiha) (National 1引言 是海洋声学研究的一个重要内容。该方法与传统的海水声速剖面直接或间接测量技术(如声速仪、温度计、 CTD或XBT等)相比的一个重要特点是通过声学这个积分探头可以实现较大尺度内海水声速变化的监测。 Tolstoy[2]等人提出了用匹配场方法反演海水声速剖面，并指出在介质水平缓变海域，对于计算声场 而言应用多个不同的声速剖面与用其“平均声速剖面”是一样的，可以把介质当作水平不变的，直接反演 介质的“平均声速剖面”，这样就避开了介质水平变化的复杂性。但是以往的研究往往采用垂直阵实现海 水声速剖面的反演，本文采用水平阵进行海水声速剖面的反演研究。 2水平阵匹配场反演海水声速分布方法 在利用匹配场进行声学反演时，需要解决三个问题：代价函数，待反演参数的选取及搜索算法。在本 文中，代价函数为 M盯N l cosT(c(z))=l∑∑∑Q『，。∽碡，。‘以，cG蚓∞)一 lm=ll=m“i=1 l 其中：Q『．，∽)=ql∽k。’∽)为水平阵上第，和第m号阵元接收到的声信号的互功率谱 阵元对应的声信号的互功率谱，且f为水昕器阵的阵元个数，Ⅳ为频率点数，上标“}”表示复数共轭．K 是归一化系数，满足足=、／l∑∑∑9，。∽)9，‘∽』∑∑∑正．。(，：，c(z鹏，。+∽，c0叫。 1 l=l Y lm=l，…m+l lira=1，-m+l 海水声速剖面是深度的函数，为了减少待反演参数，可以采用经验正交函数表示海水声速剖面 c-)=羽+∑a；g；0) (2) f 其中ciz)为平均海水声速分布，而gi(z)为第f阶经验正交函数，aj为待反演的经验正交函数的系数。 利用匹配场技术反演海水声速剖面往往是一个多维反演问题，所以搜索算法的选取是一个重要的问 题。目前已发展了包括模拟退火、遗传算法等在内的大量全局搜索方法。但是为了解决目前模拟退火和遗 传算法等搜索方法的稳定性与可靠性问题，本文采用基于微扰理论的穷举法。其优点在于：1)计算了所 有可能的待反演参量对应的代价函数，反演结果不会陷入局部极值：2)反演结果是确定性的，不会随不 同次的反演得到不同的结果；3)可以获得全局极值与局部极值的分布信息，该信息可以为其它的分析(比 如信噪比、海洋环境失配等)提供基础。4)基于声场微扰理论的穷举反演方法可以大大提高反演速度。 151 3实验结果 2003年8月中国科学院声学研究所声场声信息国家重点实验室进行了一次浅海声学实验，其实验示意 图如图1。长度为45米的32元水平接收阵布放于海底，气枪声源到接收器的距