[发明专利]基于散射声场分离算法在浅海信道中测量方法

说明书

本发明公开了基于散射声场分离算法在浅海信道中测量方法，包括以下步骤：建立基于弹性结构的浅海信道散射声场模型，采用双层柱面阵列进行声压数据采样，提取浅海信道中两同轴柱形全息面上的复声压数据。提取两柱形全息面数据，计算全息面声压角谱。采用柱面声场分离技术对两柱形测量全息面声压数据进行分解，分离得到散射声场。重构散射声场，计算声场重构误差。对重建面的声压角谱进行波数域加窗，对加窗后的重建声压角谱进行傅里叶逆变换获得重建面声压。通过有限元软件对声场进行仿真，获取更为精准的两柱形全息面复声压数据，通过声场分离算法，有效的分离得到入射声和散射声，重构散射声场。

技术领域

本发明涉及声全息领域中浅海信道中目标识别的测量方法，特别涉及基于散射声场分离算法在浅海信道中测量方法。

背景技术

声全息是将全息照相理论引用到声学领域而形成的声成像技术，根据成像距离的不同分为：近场声全息、远场声全息、常规声全息。与传统的全息技术相比，近场声全息由于测量面位置的选取，使其既包含了远场传播波成分，又包含了近场倏逝波成分，从而突破瑞利距离的限制,可以精确的识别噪声源和可视化空间声场。由于这样的优点和特点，使得近场声全息在声源的分析方面具有广阔的应用前景，成为噪声源识别和定位的重要方法之一。近场声全息技术的核心是其全息变换算法。其基本原理是利用基阵测量得到的目标辐射声场的声压或振速等声学量，利用不同的重构算法，进行声场的可视化重建或预测。

常规的近场声全息对环境要求比较严格，要求声源都在测量面的同一侧，另一侧为自由场环境。而在实际工程的测量环境下，这一点是很难保障的。此时进行重建的结果会因为存在其他目标的反射、散射等干扰，重建精度难以保障，甚至会出现完全影响判断。因此，对于非自由场的NAH技术研究日益成为关键问题。

对于目标声散射问题，最原始的是积分方程法，其基本原理是Helmholtz积分公式，理论上可计算任意复杂形状目标的散射声场，但积分方程本身的求解比较困难，但是其有明显的缺点，解的不稳定性和巨大的计算量。Rayleigh提出分离变量，结合边界条件求解简正级数，该方法可对无限长圆柱、球等给出严格解析解，仅适用于规则形状目标，此外，该方法的解为无穷级数形式，导致解的收敛性比较差。已经提出的严格理论解都只适用于一些简单目标。于是，对于复杂目标的散射声场计算，出现了多种数值解法和近似解法。比如：T矩阵法、物理声学计算方法、亮点模型和板块元理论等方法。有限元法是一种十分常见的数值计算方法，通过对区域离散化近似，线性叠加来完成计算，但在单独处理散射声场问题时，尤其是复杂、大尺寸目标和高频计算时，计算速度慢，对计算机硬件的要求高。相较于有限元方法，边界元法计算精度高，计算量减小，但面临大阶次矩阵的求逆计算，因而对具有复杂形状的目标建模十分困难。

发明内容

本发明的目的在于提供基于散射声场分离算法在浅海信道中测量方法，通过声场分离算法，有效的分离得到入射声和散射声，重构散射声场。为解决目标散射声场特性分析和测量问题提供新的研究途径，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：

基于散射声场分离算法在浅海信道中测量方法，包括以下步骤：

步骤一：建立基于弹性结构浅海信道散射声场模型，采用双层柱面阵列进行声压数据采样，提取浅海信道中两同轴柱形全息面上的复声压数据；

步骤二：提取两柱形全息面数据，计算全息面声压角谱；

步骤三：采用柱面声场分离技术对两柱形测量全息面声压数据进行分解，分离得到散射声场；

步骤四：重构散射声场，计算声场重构误差。对重建面的声压角谱进行波数域加窗。对加窗后的重建声压角谱进行傅里叶逆变换获得重建面声压。

进一步地，步骤一结合待分析的模型及有限空间三维尺寸，建立弹性结构浅海信道声场模型，并进行网格划分，水域网格按照一个波长内不少于六个点的规则，采用自由四面体网格进行网格划分。