[发明专利]一种三维地层成像的相控方位滤波方法

权利要求书

1.一种三维地层成像的相控方位滤波方法，其特征是：针对相同距离，不同角度、单快拍、多相干源信号，采用最大后验概率准则获得在空间各扫描方位上的源信号最优波形估计，根据各方位上信号的能量大小来判断目标的方位；包括步骤如下：

1)数据预处理模块：根据实际需要，设定系统参数，包括信号类型及脉宽，目标深度范围，分层距离带，方位角度范围，最小方位可分辨角度，选取到达回波数据的时间样本点，将每层到达的地层回波信号均匀分成L_θ个方位栅格，方位栅格集合不大于发射换能器发射指向性宽度；

2)阵列矩阵变换模块：包括等式方程步骤、参数关系式步骤两部分，通过转化获得接收回波矩阵空间、接收回波的转化矩阵空间、目标空间场矩阵、加权后的目标空间场矩阵、特征空间向量、特征值、阵列流形矩阵，建立上述参数之间方程表达式及参数关系式；

3)MAP准则估计模块：以最大后验概率准则来估计方位上的回声强度，根据最大后验概率准则，建立代价函数和约束方程；

4)迭代收敛模块：根据代价函数和约束方程，转化为无约束最小化求解问题，反复迭代使代价函数最小，获得分层距离带上的方位栅格回波信号强度的最优估计值。

2.根据权利要求1所述的三维地层成像的相控方位滤波方法，其特征是：

(1)地层回波信号均匀分成L_θ个方位栅格L_θ＞＞M阵元数目，s₁(t)，s₂(t)，...，为不同方位栅格的源信号，满足远场条件，m∈(1，2，...，M)，l∈(1，...，L_θ)为一系列的源信号，1，…，L_θ源信号，在各阵元，1，…M阵元，产生的相位差的频域表达式，n_1...M(t)噪声序列为零均值高斯过程，阵元间噪声相互独立，阵列各传感器为各向同性，无互耦以及通道不一致干扰；方程等式进一步表达如下：

Z(t)=Σl=1Lθv(φl)sl(t)=V(φ1,...,φLθ)·S(t)+N(t)⇒Z=V(S+ΔS)=VS^]]>

其中Z(t)＝[z₁(t)，…，z_M(t)]^T是M×1的阵元接收信号向量，是L_θ×1的源信号向量，是第l个源信号在M阵元上的M×1的阵列流形向量，是M×L_θ的阵列流形矩阵，N(t)＝[n₁(t)，…，n_M(t)]^T是M×1的噪声向量，为混入不相干噪声后的目标空间场矩阵，φ_l＝2πf_cdcosθ_l/c，其中f_c为中心频率，d为阵元间距，c为水中声速，以最大后验概率准则(MAP)来估计θ_l方位上的回声强度以方位输出谱有峰值的表示该方位有强散射目标，而峰值的高低表示对该方位目标散射强度的估计值，选取代价函数和约束方程

(2)获得的接收回波矢量、阵列流形矩阵、目标空间场矩阵三者的方程等式将V^HV作特征分解，V^HV＝PΛP^H，其中λ_i是第i个特征值，e_i是第i个归一化特征向量，该过程通过对V进行奇异值分解来实现，假设特征值以降序排列：把特征值分为两组：非零特征值(λ₁，λ₂，…，λ_fix)和零特征值两组特征值的个数分别为L_fix和L_free＝L_θ-L_fix，特征向量P分为两个子阵，P_fix(L_θ×L_fix)和P_free(L_θ×L_free)，分别对应于非零特征值和零特征值；

(3)采用矩阵等效替换，距离分层上各阵元接收回波矢量引入阵列流形矩阵共轭值和特征空间向量共轭信息，方程等式左边获得接收回波的转化矩阵空间Z′，将等式右边特征空间向量共轭信息和目标空间场矩阵的乘积，方程等式右边获得转化后的加权后的目标空间场矩阵S′，进一步矩阵等效推导如下：

⇒PHVHZ=ΛPHS^⇒QZ=ΛPHS^(Q=PHVH)⇒Z′=ΛS′(Z′=QZ,S′=PHS^)]]>

(4)获得对每个时刻t的接收数据矢量Z(t)，计算Z′，通过上述转化，利用非零特征值子空间和零特征值子空间相互正交的特性，Z′_fix和S′_fix为非零特征值子空间，Z′_free和S′_free为零特征值子空间，将接收回波的转化矩阵空间Z′和加权后的目标空间场矩阵S′都进行子空间展开为：

Zfix′Zfree′=ΛfixΛfreeSfix′Sfree′=ΛfixSfix′ΛfreeSfree′]]>

(5)获得Z′和S′的子空间中主要能量分布的Z′_fix和S′_fix关系式，每个元素可写为：

Zfix′=ΛfixSfix′⇒Sfix′=Λfix-1Zfix′]]>

zfix,l′=λlsfix,l′⇒sfix,l′=1λlzfix,l′,]]>l＝1，...，L_fix

计算s′_fix＝Λ^-1_fix·z′_fix，其中Λ^-1_fix＝diag(λ₁^-1，λ₂^-1，…，λ_fix^-1)

(6)加权后的目标空间场矩阵中能量分布S′_fix由接收回波的转化矩阵空间Z′得到，依次获得目标空间场矩阵；

S′=PHS^⇒S^=PS′=PfixPfreeSfix′Sfree′PfixSfix′+PfreeSfree′]]>

初始化向量S′_free，所有元素置零，计算计算对角矩阵

(7)对每个时刻t的接收数据矢量Z(t)，为得到的最优化估计结果，需要保证代价函数Γ最小，将代价函数和约束方程的的解可看作一个关于S′_free的无约束最小化求解问题，即求取S′_free，满足

有：∂Γ∂s^l=2s^l|s^l|2,]]>写成矩阵形式为：∂Γ∂S^=2T-1S^]]>

⇒∂Γ∂Sfree′=PfreeH∂Γ∂S^2PfreeHT-1S^=0,]]>其中T=diag(|s1|2,|s2|2,···,|sLθ|2)]]>

Sfree′=-(PfreeHT-1Pfree)-1(PfreeHT-1Pfix)Sfix′]]>

⇒=(PfreeHTPfix)(PfixHTPfix)-1Sfix′]]>

计算Sfree′=(PfreeHTPfix)(PfixHTPfix)-1Sfix′;]]>

(8)重复(6)至(7)步，直到S′_free收敛，最后得的估计量；其中(1)至(3)步只做一遍，而(4)至(8)要对选取的各个距离分层分别运算。