一种基于反射信号重构的薄层介质参数反演方法

技术特征：
1.一种基于反射信号重构的薄层介质参数反演方法，其特征在于，该方法用于在未知先验信息的情况下，反演地下薄层介质参数，所述先验信息包括探测场景材料和实际厚度，所述方法包括如下步骤：首先通过广义反射系数的定义以及地下分层的褶积模型重构出反射信号模型，该模型包含了层状介质的物理参数；然后构建优化问题，所述优化问题的代价函数为模型和实际反射信号之间的误差平方，通过遗传算法求解所述优化问题，使所述误差平方最小化；所述通过遗传算法求解所述优化问题在初值范围内进行优化。2.如权利要求1所述的一种基于反射信号重构的薄层介质参数反演方法，其特征在于，该方法适用于公路路面、桥梁内部分层结构以及其他包含薄层的分层介质探测场景，同时也适用于非薄层的介质探测。3.如权利要求1所述的一种基于反射信号重构的薄层介质参数反演方法，其特征在于，所述通过广义反射系数的定义以及地下分层的褶积模型重构出反射信号模型，具体为：地下各介质层都是线性、均匀、各向同性的，且为非磁性无损材料；则第i层和第i 1层之间界面的反射系数r
i,i 1
和透射系数分别为τ
i,i 1
：电磁波的衰减常数为：其中：ω表示角频率；ε
ri
表示第i层的相对介电常数，ε
ri 1
表示相对介电常数；σ
i
表示第i层的电导率；μ
i
表示第i层的磁导率；定义广义反射系数为接收到的每层分界面的反射波振幅与发射波的比值，记为r
i,i 1
；由于第一层顶部上方为空气介质，因此有r
0,1
＝r
0,1
，r
0,1
表示第0层和第1层之间界面的反射系数；得第i个广义反射系数r
i,i 1
为：式中d
k
为第k层的层厚；k取值为1～i；α
k
为第k层的电磁波的衰减常数；r
k,k 1
表示第k层和第k 1层之间界面的反射系数；层状介质的反射系统简化为线性时不变lti系统，则接收信号y
r
(t)为：式中，x(t)为t时刻的雷达发射信号，x(t-t
i
)为t-t
i
时刻的雷达发射信号；n为总分层个数；n(t)表示高斯白噪声。
4.如权利要求1所述的一种基于反射信号重构的薄层介质参数反演方法，其特征在于，所述构建优化问题，所述优化问题的代价函数为模型和实际反射信号之间的误差平方，通过遗传算法求解所述优化问题，使所述误差平方最小化，具体为：最优化问题为：其中，o(l)为包含参数的代价函数，l为待优化的参数矩阵；o(l)表达式为：其中：y
r
(t)为模型信号；y(t)为实际反射信号；|| ||2为l2范数；所述优化问题的物理意义为：求取最优参数，使得重构的反射信号与实际接收的反射信号拟合度最高。5.如权利要求4所述的一种基于反射信号重构的薄层介质参数反演方法，其特征在于，所述通过遗传算法求解所述优化问题在初值范围内进行优化，具体为：将介质层双程发射时间t
i
与电性参数ε
i
和σ
i
作为优化对象，在给定初值范围的条件下得到参数最优解；利用广义反射系数频谱估计参数，得到反演初值：通过对广义反射系数序列的频谱进行分析，初步估计出待反演参数的值，并以此作为遗传算法的输入参数初值，在采用遗传算法求解所述优化问题在初值范围内进行优化；其中分层介质第一层底面和第二层底面的反射系数序列为r(t)：r(t)＝r
12
δ(t-t1) r
23
δ(t-t2)其中r
12
和r
23
为广义反射系数，δ(t-t1)和δ(t-t2)分别为t1，t2时刻的冲激响应；对反射系数系列作傅里叶变换即得到反射系数序列谱r(f)：其中r(f)为时间零点在中间层的中点位置处的广义反射系数频谱；j为复数符号；将反射系数序列谱r(f)分为幅度谱|r(f)|和相位谱θ(f)：θ(f)＝arctan[tan(πfδt)r
o
/r
e
]-2πft
1-πfδt其中：t1为第一层介质的电磁波双程反射时延，δt＝t
1-t0为第二层介质的电磁波双程反射时延，r
e
和r
o
为广义反射系数对进行奇偶分解的结果：对幅度谱|r(f)|进行求导并令其等于0：
解得频率为n为正整数，n＝1,2,...；f为幅度谱极大值点和极小值点所对应的频率，且为周期性分布，即相邻的两个极大/小值的间隔为幅度谱的凹陷周期δf：利用时间厚度与幅度谱凹陷周期之间的反比关系，利用反射系数频谱的周期来估计出时间厚度；θ(f)分为两个部分，第一部分包含参数δt和反射系数奇偶分量，影响相位谱的振幅大小；第二部分包含参数δt和t1，影响相位谱的变化速度；反射系数奇偶分量包含介电常数和与电导率相关的衰减常数，则通过相位谱估计出三层介质的介电常数和中间层的电导率；得到的反演初值包括上述估计的时间厚度、三层介质的介电常数和中间层的电导率。

技术总结
本发明公开了一种基于反射信号重构的薄层介质参数的方法，适用于含有薄层的分层介质模型，能够解决薄层的时延估计问题，且该方法精度和收敛速度均较高，能够获得更好的薄层参数反演性能，是一种高效准确的反演方法。该方法用于在未知先验信息的情况下，反演地下薄层介质参数，先验信息包括探测场景材料和实际厚度，方法包括如下步骤：首先通过广义反射系数的定义以及地下分层的褶积模型重构出反射信号模型，该模型包含了层状介质的物理参数。然后构建优化问题，优化问题的代价函数为模型和实际反射信号之间的误差平方，通过遗传算法求解优化问题，使误差平方最小化；通过遗传算法求解优化问题在初值范围内进行优化。求解优化问题在初值范围内进行优化。求解优化问题在初值范围内进行优化。


技术研发人员：杨小鹏 李奕璇 刘仁杰 兰天 渠晓东
受保护的技术使用者：北京理工大学
技术研发日：2022.03.03
技术公布日：2022/7/9