自适应有限元GPR频率域正演方法技术

本发明专利技术公开了一种自适应有限元GPR频率域正演方法，包括设定正演方法的参数；计算当前网格的数值解和局部误差指示值；标记需细化的单元集合；对标记网格进行细化；进行数值求解得到最终的自适应有限元GPR正演结果。本发明专利技术提供的这种自适应有限元GPR正演方法，在后验误差估计策略的基础上，采用自适应有限元方法开展GPR频率域正演，将EPML边界条件引入到频率域正演计算当中，简化了参数优化过程，使得吸收层厚度较小情况下便能达到较好的吸收效果，能够降低节点自由度，将很大程度上节约计算成本。

ADAPTIVE FINITE ELEMENT GPR FREQUENCY DOMAIN FORWARD METHOD

【技术实现步骤摘要】
自适应有限元GPR频率域正演方法
本专利技术具体涉及一种自适应有限元GPR频率域正演方法。
技术介绍
探地雷达(GroundPenetratingRadar，GPR)是利用高频电磁波对地下结构或者物体内部不可见目标体进行定位的一种无损探测技术(Daniels，2004；曾昭发，2010)。通过GPR正演，可以加深对GPR传播规律与反射剖面的认识，提高雷达数据解释水平(Giannopoulos,2005；Irving&Knight，2006)。GPR的数值模拟可分为时间域与频率域两大类。其中时间域模拟方法应用较为广泛，包括时域有限差分(FiniteDifferenceTimeDomain，FDTD)(Chen&Huang,1998；Cassidy&Millington,2009；Diamanti&Giannopoulos，2009；Linetal.，2012；Zhangetal.，2016；Warrenetal.，2016)、时域有限单元(FiniteElementTimeDomain，FETD)(Feng&Wang,2017；Fengetal.本文档来自技高网...

【技术保护点】
1.一种自适应有限元GPR频率域正演方法，包括如下步骤：S1.设定正演方法的参数；S2.根据步骤S1设定的参数，计算当前网格的数值解；S3.计算局部误差指示值；S4.标记需细化的单元集合；S5.对步骤S4标记的网格进行细化；S6.在步骤S5进行细化后的网格基础上，进行数值求解，从而得到最终的自适应有限元GPR正演结果。

【技术特征摘要】
1.一种自适应有限元GPR频率域正演方法，包括如下步骤：S1.设定正演方法的参数；S2.根据步骤S1设定的参数，计算当前网格的数值解；S3.计算局部误差指示值；S4.标记需细化的单元集合；S5.对步骤S4标记的网格进行细化；S6.在步骤S5进行细化后的网格基础上，进行数值求解，从而得到最终的自适应有限元GPR正演结果。2.根据权利要求1所述的自适应有限元GPR频率域正演方法，其特征在于步骤S1所述的参数包括正演频率、模型参数和初始网格。3.根据权利要求2所述的自适应有限元GPR频率域正演方法，其特征在于所述的初始网格，具体为采用Delaunay三角形网格，并对网格进行剖分从而得到非结构化网格剖分结果。4.根据权利要求1～3之一所述的自适应有限元GPR频率域正演方法，其特征在于步骤S2所述的计算当前网格的数值解，具体为采用如下步骤进...

【专利技术属性】
技术研发人员：冯德山，丁思元，王珣，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43