CN114858926B 一种基于频域逆时偏移的全矩阵数据高质量成像方法 （浙江大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN114858926B(45)授权公告日2025.07.04

(21)申请号202210504512.6

(22)申请日2022.05.10

(65)同一申请的已公布的文献号申请公布号CN114858926A

(43)申请公布日2022.08.05

(73)专利权人浙江大学

地址310027浙江省杭州市西湖区浙大路

38号

(72)发明人赵朋纪凯鹏卓超杰叶盛焦晓龙傅建中

(74)专利代理机构杭州知闲专利代理事务所(特殊普通合伙)33315

专利代理师王于敏

(51)Int.CI.

GO1N29/06(2006.01)

G01N29/07(2006.01)

G01N29/44(2006.01)

GO6F17/14(2006.01)

GO6F17/16(2006.01)

(56)对比文件

CN111239246A,2020.06.05

CN113552219A,2021.10.26审查员翁永超

权利要求书2页说明书9页附图5页

(54)发明名称

一种基于频域逆时偏移的全矩阵数据高质量成像方法

(57)摘要

CN114858926B本发明提供一种基于频域逆时偏移的全矩阵数据高质量成像方法，对激发信号、声速分布和全矩阵数据分别进行预处理，得到激发声源矩阵、阻抗矩阵和接收声源矩阵；以前述预处理结果作为输入，对阻抗矩阵进行LU分解处理，求解得到发射波场和接收波场，最后利用成像条件得到成像结果。该成像方法利用同一频率下阻抗矩阵固定的特性，通过LU分解可一次性的求出同一频率下多次激发产生的声压场，大大节省了计算时间，提高了计算效率。采用全波动方程为基础进行波场正演，对波场传播的简化少，适应于任

CN114858926B

□…(C)

16(0)…,P(

CN114858926B权利要求书1/2页

2

1.一种基于频域逆时偏移的全矩阵数据高质量成像方法，其特征在于，包括以下步骤：

(1)输入激发信号、声速分布和全矩阵数据；

(2)对激发信号和全矩阵数据分别进行预处理，分别得到激发声源矩阵和接收声源矩阵；

对激发信号进行如下预处理：

对激发信号进行时域上的傅里叶变换处理，得到激发声源矩阵；

对全矩阵数据进行如下预处理：

先采用窗函数作用于全矩阵数据来滤去直达波，将滤去直达波的全矩阵数据在时间维度上进行傅里叶变换，然后将变换得到的频域中的数据取复数共轭实现逆时的效果，整理得到接收声源矩阵；

(3)对声速分布进行离散化处理，并计算得到阻抗矩阵；

步骤(3)中，阻抗矩阵A(o)的计算公式如下：

A(w)=▽2+w2/c2

其中，▽2表示离散化处理后得到的拉普拉斯算子；の表示频率；c表示声速；

(4)根据得到的激发声源矩阵和接收声源矩阵，采用阻抗矩阵的LU分解，分别求解得到发射波场和接收波场；

步骤(4)的具体步骤如下：

4.1取选择频带范围内任一未遍历频率作为当前频率，以当前频率对应的激发声源矩阵、接收声源矩阵和阻抗矩阵作为输入，对阻抗矩阵进行LU分解，求解得到当前频率对应的激发声压场和接收声压场；

4.2遍历选择频带范围内的所有频率，得到发射波场和接收波场；

(5)采用成像条件处理发射波场和接收波场，得到初始图像，对初始图像进行后处理，得成像结果；

对初始图像进行如下后处理：

先采用拉普拉斯滤波对初始图像进行锐化处理，再对锐化处理后的初始图像进行希尔伯特变换得到图像的包络，最后采用中值滤波滤去椒盐噪声后即得所述成像结果。

2.根据权利要求1所述的基于频域逆时偏移的全矩阵数据高质量成像方法，其特征在于，对声速分布进行离散化处理时，通过在离散化模型的四周添加完美匹配层对阻抗矩阵进行优化。

3.根据权利要求1所述的基于频域逆时偏移的全矩阵数据高质量成像方法，其特征在于，步骤4.1中，根据声源矩阵求解其产生的声压场的计算公式如下：

L(o)U(o)[P?(a),P?(a),…,P(a)]=[S?(a),S?(a),…,S,(a)]

其中，n表示激发次数；L(a)和U(a)表示阻抗矩阵A(a)进行LU分解后得到的两个因子；[S?(a),S?(a),…,S。(a)]表示激发声源矩阵或接收声源矩阵；[P?(a),P?(a),…,P

(o