基于EMD噪声互相关提取时域格林函数研究.docVIP

基于EMD噪声互相关提取时域格林函数研究 　　摘要本文在前人研究的基础上，基于简正波理论，给出了海洋环境噪声互相关函数(NCF)及时域格林函数(TDGF)的公式，分析了NCF与TDGF的相似性；基于青岛胶州湾航道附近采集的海洋环境噪声数据，分析相关的时空特性；提出了经验模态分解（EMD）与波束形成相结合的环境噪声互相关技术，提高了相关峰信噪比，成功提取了海表噪声场的直达垂直波至的时域到达结构，与理论上的时域到达结构相符；并进一步反演得到声速结果，与同步观测的声速测量结果十分吻合。 　　【关键词】时域格林函数 噪声互相关函数经验模态分解 　　1 引言 　　近年来，国外许多研究证明，对两个接收点记录到的长时间的海洋环境噪声进行互相关处理，得到的长时间互相关函数可近似表示这两点之间的格林函数。2004年，Roux, Kuperman首次通过实验证明，使用海洋环境噪声可近似地估计两点之间的时域格林函数（TDGF）的波前到达结构。随后，国外许多研究机构对使用NCF提取TDGF的技术及其应用进行了研究。而国内对于该方面的研究基本都处于起步阶段，对这方面的研究报道也不多见。 　　实际上，在海洋中每个独立噪声源产生的声场中，声波经过长距离传播后在不同的两个点接收到的噪声信号是相关的，只是这种微小的相关成分，被整个海洋中广泛分布的非相关声场所淹没。然而，对两个接收点记录的海洋环境噪声进行长时间的互相关处理可以将这种相关成分提取出来，得到的相干波前与格林函数的结构一致。本文的主要工作是利用海洋环境噪声提取格林函数，将经验模态分解(EMD)与波束形成相结合的方法应用到格林函数提取中以提高信噪比,优化提取方法，并进一步进行声速反演。 　　2 原理与方法 　　2.1 NCF与TDGF的相似性 　　本文采用了如下海洋模型，仅有一层水体覆盖于半无限大的海底之上，且水体、海底的密度以及其中的声速分别为ρ1，c1(z)，ρ2，c2(z)，海面为压力释放界面，海底为刚性海底。 　　图2.2.1为 噪声场模型示意图，噪声源平面位于海面下深度z’处，海洋中的两个点分别位于(r1,z1),(r2,z2)处 　　图中海洋的声学环境可看作一个波导。海洋当中位于深度z1和z2，水平间距为R的两个场点之间的格林函数可以用简正波形式表示为： 　　Gω(R,z1,z2)=(knR) 　　（2-1-1） 　　其中，ρ为声源所在处的海水密度，S(ω)为声源的频谱。对式（2-1-1）进行逆傅立叶变换，就可以得到时域格林函数（TDGF）的表达式： 　　Gt(R,z1,z2)= dωGt(R,z1,z2) exp(iωt)（2-1-2） 　　上面已经得到两点之间时域格林函数（TDGF）的表达式。若表面噪声源被等效为位于深度 处的一个布满点噪声源的平面，海洋环境噪声互谱密度的表达式可以表示为： 　　Gω(R,z1,z2)= 　　 　　 　　（2-1-3） 　　其中，Q2(ω)为噪声源的能量谱，z为噪声源深度，αn为模态衰减系数。 　　将式（2-1-1）与式（2-1-3）进行比较，可以很明显的看出，两个式子求和号之前的项为标量不影响格林函数和互谱密度函数的结构，二者求和号之后的项在结构上具有相似性，故而理论上可以采用海洋环境噪声场中两点之间的长时间互相关函数来提取格林函数。 　　2.2 经验模态分解（EMD）方法简介 　　本文使用经验模态分解方法，在对原始噪声数据进行经验模态分解，将其作为一种特殊滤波方法，提取噪声数据中的不同频段成分，再进行互相关处理，来分析其对互相关结果的影响及有效性。 　　经验模态分解可以在时域范围内将复杂的原始信号分解为一系列固有模态函数（Intrinsic Mode Functions，IMF），每一个固有模态函数需要满足以下两个条件： 　　(1)在固有模态函数的时间序列中，函数的过零点的数量与极值点的数量必须相等，或者最多只能相差一个， 　　(2)在任意时间点上，由该函数的局部极大值和极小值所构成的包络的均值必须等于零。 　　经验模态分解方法的本质是根据原始信号自身的时间尺度特征对信号进行分解从而获得信号的固有振动模态，然后根据这些固有振动模态来分解信号的时间序列。其步骤在此不再赘述。 　　3 实验数据分析 　　2012年11月15日，中国海洋大学海洋声学研究室在青岛胶州湾海域航道附近进行了海洋环境噪声采集实验。本文对海洋环境噪声时空特性进行分析，利用噪声互相关技术提取格林函数，并进一步进行声速反演。 　　3.1 数据预处理结合波束形成的噪声互相关 　　为了提取信噪比尽可能高的噪声互相关波至，需要进行单阵元数据预处理。单阵元数据预处理的内容主要包含以下几个方面： 　　（1）重新降频采样。满足采样定理的前提下进行降频采样，降低数据处理难度； 　　（2）去均值与带通滤