[发明专利]一种海洋环境与声场不确实性的表征和传递的快速计算方法

权利要求书

1.一种海洋环境与声场不确实性的表征和传递的快速计算方法，其特征在于：

一、海洋环境不确实性和声场不确实性的表征：

(1)海洋环境不确实性表征

用向量b表征环境参量，用b₀表示我们利用测量或反演方法获得的已知环境参量，作为不确实环 境参量的均值，不确实环境参量采用参量不确实性域来表征，它表示了环境参量落在b₀附近的 一个领域之中，用概率密度函数p(b)表征不确实环境参量的分布，表征了环境参量不确实性的 范围，下面对不确实环境参量作多项式混沌展开随机谱表征；

步骤一：将不确实性环境参量空间表示成随机变量b；

步骤二：将不确实环境参量b分解为确实性部分和随机扰动部分：

b(r，z)＝b₀+δb(r，z)    (1)

其中b₀为确定性部分，δb为随机扰动部分，r，z分别表示距离和深度；

步骤三：根据环境参量b的概率概率密度分布p(b)，对随机扰动部分作多项式混沌展开：

δb(r,z;θ)=Σs=0∞αs(r,z)Λs(ξ(θ)),---(2)]]>

式中Λ_s(ξ(θ))是对应环境参量b概率密度分布的随机多项式混沌基函数，θ为随机参量空间 α_s(r，z)是随机扰动环境参量投影到随机多项式混沌基函数得到的确定性系数；

步骤四：对无穷多阶多项式混沌展开作有限截取

对不同分布的环境参量利用对应的多项式混沌基函数展开，展开的阶数是有限的；

步骤五：多项式混沌项数确定

随机环境参量b由多个元素构成，b＝[b₁b₂...b_n]，根据最小均方误差准则，n维随机变量b所 必需的多项式项数为(n+S)！/n！S！；其中S为多项式混沌的阶数；

环境参量b选取为水层的声速梯度，则

b(r，z，θ)＝c(r，z，θ)＝c₀(z)+δc(r，z，θ)    (3)

(2)声场不确实性表征

①确实的环境参量b₀驱动声传播模型产生确实的声场p₀；

②不确实环境参量空间对应不确实声场空间

③用概率密度函数p(p)表征不确实声场的分布；

④表征了声场不确实性的范围；

同样，对不确实声场也作多项式混沌展开随机谱表征；

步骤一：不确实声场作多项式混沌展开：

p(r,z;θ)=Σq=0∞γq(r,z)Λq(ξ(θ))---(4)]]>

步骤二：对无穷多阶多项式混沌展开作有限截取

根据最小均方误差准则，对无穷多阶多项式混沌展开作有限阶Q的截取：

p(r,z;θ)≈Σq=0Qγq(r,z)Λq(ξ(θ))---(5)]]>

式中p(r，z；θ)为预测的随机声压场；

二、环境参量不确实性到声场不确实性的传递方法：

将环境参量与预测声场的多项式随机谱表征与确实性波方程相结合，得到嵌入不确实性的随机波 方程，实现环境参量不确实性到声场不确实性的传递，具体步骤如下：

(1)确定性波方程到随机波方程

步骤一：基于简正波传播模型的确定性波方程

▿2p+κ2p=0---(6)]]>

式中表示拉普拉氏算子，κ是波数，κ＝ω/c(z)，ω为角频率；

步骤二：随机扰动声速作二项式展开

对随机扰动声速作二项式展开，并保留一阶项，则不确实声速的波数可近似为，

κ~2=ω2c02(1+δc)2≈κ2(1-2δc(r,z;θ)c0(z))---(7)]]>

步骤三：确定性波方程推广到随机波方程

将近似的波数代入到确定性波方程，将确定性波方程推广到随机波方程：

∂2p(r,z;θ)∂r2+∂2p(r,z;θ)∂z2+κ2p-2κ2δc(r,z;θ)c0(r,z)p=0---(8)]]>

将式(3)和式(4)代入(8)式，并对等式两边同时乘以Λ₁(ξ(θ))，然后通过系总平均将各项投影到基空间，

可得

Σq=0Q∂2γq(r,z)∂r2<ΛqΛ1>+Σq=0Q∂2γq(r,z)∂z2<ΛqΛ1>+κ2Σq=0Qγq(r,z)<ΛqΛ1>]]>

(9)

-2κ2c0(z)<Λ12>Σq=0QΣs=0Qαs(r,z)γq(r,z)<ΛqΛsΛ1>=0]]>

将混沌基函数的正交性应用于上述方程，得到一组声场不确实性系数的耦合差分方程，

Σq=0Q∂γ12(r,z)∂r2+Σq=0Q∂2γ1(r,z)∂z2+κ2Σq=0Qγ1(r,z)-2κ2c0(z)<Λ12>Σq=0QΣs=0Qαs(r,z)γq(r,z)<ΛqΛsΛ1>=0---(10)]]>

以高斯分布随机声速进行下面公式的推导，δc(r，z，θ)＝α₀He₀+α₁He₁+…＝σξ，ξ是均值为 0，方差为1的高斯随机变量；将δc代入(10)式给出第1个不确实系数的耦合方程：

∂2γ1(r,z)∂r2+∂2γ1(r,z)∂z2+κ2γ1(r,z)-2κ2σc0<He12>Σq=0∞γq(r,z)<HeqξHe1>=0---(11)]]>

利用Hermite基函数的回归特性和正交性，计算耦合项的系总平均，得到系数γ₁的偏微分回归 方程：

∂2γ1(r,z)∂r2+∂2γ1(r,z)∂z2+Dγ1(r,z)-Elγl-1(r,z)-El+1γl+1(r,z)=0---(12)]]>

式中D≡κ²，E≡2κ²σ/c₀；用向量-矩阵形式，上述偏微分方程可化简为

∂2γ∂r2+∂2γ∂z2+Aγ=0---(13)]]>

式中，A为环境参量不确实性到声场不确实性的耦合方程，A是对称的三角阵；

(2)求解随机波方程得到不确实声场

步骤一：耦合矩阵A特征值分解

A=Σq=1QλqgqgqT---(14)]]>

式中，上标T表示转置，G＝[{g₁}，{g₂}，…]为耦合矩阵A的特征向量矩阵，GG^-1＝I，(13) 式可以表示为

∂2G-1γ∂r2+∂2G-1γ∂z2+G-1AGG-1γ=0---(15)]]>

令上式可以简化为

∂2γ^∂r2+∂2γ^∂z2+Ωγ^=0---(16)]]>

式中，Ω＝G^-1AG为对角阵，对角阵元素为特征值λ₁，λ₂，…，利用特征值和特征向量分解技术，得 到第1个确定系数的波方程

∂2γ^1∂r2+∂2γ^1∂z2+λ1γ^1=0,---(17)]]>

步骤二：确定系数求解

等式(17)与标准波方程(6)有相同的形式，因而求解过程与标准的波方程求解过程是一样的；的 解表示为

γ^u(r,z)=i8πre-iπ4au0Σm=1Mψm(zs)ψm(z)krmeikrmr---(18)]]>

式中k_rm、ψ_m(z)分别是水平波数和模深度函数，由简正波传播模型求解得到，a_u0表示矩阵A的特征 向量的逆的第一列中的第u个元素，在声速梯度不确实情况下，模深度函数没有受到不确实声速的影响， 但声速梯度不确实导致了水平波数的不确实，其中k_zm是深度波数，根据γ₁与的 关系式，系数γ₁是的线性组合

γ1=Σu=0Qg1uγ^u(r,z)---(19)]]>

步骤三：Q个多项式混沌基函数的确定性系数加权累加得到不确实性声场

最后，在远场近似下的不确实性声场可表示为：

p(r,z;θ)=Σq=0Q[Σu=0Qgquγ^u(r,z)]Λq(ξ(θ))]]>

(20)

=Σq=0Q[Σu=0Qgqu(i8πre-iπ4au0Σm=1Mψm(zs)ψm(z)λu-kzm24eiλu-kzm2r)]Λq(ξ(θ))]]>

若SSP是确实的，则上式中只有γ₀是不等于0的，此时式(20)收敛到声速确知时的声场简正模 表达式；

步骤四：不确实性声场的均值和相关函数

由多项式混沌展开的不确实声场是随机声场，其均值可表示为

<p(r,z;θ)>=γ^0(r,z)---(21)]]>

相关函数为

Rpp(r1,z1;r2,z2)=Σq=1Qγ^q*(r1,z1)γ^q(r2,z2)---(22).]]>