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

专利信息
申请号: 201010619550.3 申请日: 2010-12-22
公开(公告)号: CN102176008A 公开(公告)日: 2011-09-07
发明(设计)人: 朱必波;郑乐一;丛卫华;蒋飚;张峰山;傅翔毅;张卫华 申请(专利权)人: 中国船舶重工集团公司第七一五研究所
主分类号: G01S7/539 分类号: G01S7/539;G01S15/89
代理公司: 杭州九洲专利事务所有限公司 33101 代理人: 陈继亮
地址: 310012 浙*** 国省代码: 浙江;33
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 三维 地层 成像 方位 滤波 方法
【权利要求书】:

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

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

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

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

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

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

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

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

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

(2)获得的接收回波矢量、阵列流形矩阵、目标空间场矩阵三者的方程等式将VHV作特征分解,VHV=PΛPH,其中λi是第i个特征值,ei是第i个归一化特征向量,该过程通过对V进行奇异值分解来实现,假设特征值以降序排列:把特征值分为两组:非零特征值(λ1,λ2,…,λfix)和零特征值两组特征值的个数分别为Lfix和Lfree=Lθ-Lfix,特征向量P分为两个子阵,Pfix(Lθ×Lfix)和Pfree(Lθ×Lfree),分别对应于非零特征值和零特征值;

(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′都进行子空间展开为:

ZfixZfree=ΛfixΛfreeSfixSfree=ΛfixSfixΛfreeSfree]]>

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

Zfix=ΛfixSfixSfix=Λfix-1Zfix]]>

zfix,l=λlsfix,lsfix,l=1λlzfix,l,]]>l=1,...,Lfix

计算s′fix=Λ-1fix·z′fix,其中Λ-1fix=diag(λ1-1,λ2-1,…,λfix-1)

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

S=PHS^S^=PS=PfixPfreeSfixSfreePfixSfix+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)要对选取的各个距离分层分别运算。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国船舶重工集团公司第七一五研究所,未经中国船舶重工集团公司第七一五研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201010619550.3/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top