共查询到18条相似文献，搜索用时 296 毫秒
1.
完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果. 相似文献
2.
我国城市地下管线种类繁多、规模巨大.近年来,因城市管网家底不清、位置不明导致的管道挖断事故频发.探地雷达(GPR)是目前最为有效的地下设施探测定位设备,通过对探地雷达回波信号进行反演分析,可获取管道埋深、位置等信息.然而,地下管道等地下圆形目标,传统正演模型中阶梯近似方法会产生一定误差.本文结合辛Euler算法与曲面共形技术,建立探地雷达电磁波在含非金属管道地下结构中传播的精细化数值模型,并通过图形处理器(GPU)加速技术进一步提升模型计算效率.数值模拟结果显示,共形网格技术有效降低了由于阶梯近似造成的虚假绕射波,并通过结合GPU并行计算,在不增加网格密度的前提下,实现了对地下管道电磁响应的高效精细化计算. 相似文献
3.
《地球物理学进展》2017,(4)
探地雷达是一种非常重要的管线探测技术,为了提高地下管线雷达图像特征的认识,确定管线异常体的位置,提高雷达资料的解释精度.论文从Maxwell两个旋度方程出发,推导了二维TM波的差分方程、CFL数值稳定性条件、频散关系.然后,基于Matlab平台编写了探地雷达正演的FDTD程序,应用该FDTD程序开展了管线探测中探地雷达探测效果分析,包括对管线埋藏深度、管线间距、管线内物质、管线材质等影响因素的数值模拟.通过分析雷达正演剖面特征,可以清晰了解并掌握雷达管线探测与各种影响参数之间的关系,对实际地下管线探测可起到指导作用.最后,将GPR应用于武广高速浏阳河隧道管线探测中,GPR准确地定位了PVC通迅电缆的位置在埋深,为工程施工与处置提供了依据. 相似文献
4.
探地雷达正演模拟在真实雷达数据解译及全波形反演中扮演着重要的角色，针对传统探地雷达(Ground Penetrating Radar, GPR)正演模拟计算量巨大、耗时、不利于实时探测等问题，提出一种基于机器学习框架的近实时GPR正演模拟方法.以混凝土中的钢筋探测作为GPR应用场景，混凝土的含水量、钢筋半径及埋地深度作为模型参数，利用时域有限差分数值模拟散射回波信号；运用主成分分析对回波数据进行降维处理得到相应的主成分权值系数，并以此作为机器学习网络的输出；设计了一种基于随机森林的多层循环网络架构和学习策略，不仅充分挖掘学习模型参数和主成分权值系数之间的内在因果关系，也共享主成分间的相互联系，并具有对每个预测主成分完善和修正的功能，以此实现基于机器学习的探地雷达快速正演模拟，与传统机器学习相比，有效提高了正演模拟的精度.在此基础上将两个深度神经网络与随机森林相结合，以回波数据主成分系数为输入，建立了基于机器学习的场景参数预测模型，实现了近实时的埋地目标探测，预测的混凝土含水量最大误差为2%,钢筋埋地深度最大误差为6.7%. 相似文献
5.
《应用地球物理》2020,(1)
探地雷达是一种高效快速无损的浅地表勘探方法,为了提高探地雷达探测的解释精度,需要开展探地雷达高精度的数值模拟方法研究。因此,本文将第二代小波有限元引入探地雷达的正演模拟中,以提升方案构造的第二代小波尺度函数代替多项式函数作为有限元基函数,具有多尺度、多分辨的特性。该算法可根据实际需要,任意改变分析尺度,在大梯度处采用小的分析尺度以提高分析精度,而在小梯度处采用大的分析尺度以提高分析效率,从而有利于捕捉解的局部突变特征,在不改变网格剖分的前提下提高分辨率,实现GPR正问题的高效求解。将该算法应用于具有解析解的线电流辐射源及隧道衬砌不密实模型的数值模拟中,结果表明:第二代小波有限元的求解结果与解析解及常规有限元算法求解结果都能较好地吻合,验证了第二代小波有限元算法的正确性,而且第二代小波有限元可以根据实际问题的需要任意改变分析尺度,不需要网格单元的重新剖分,是一种优于传统单元网格加密和阶次升高的自适应算法,特别适合于求解大梯度、奇异性的复杂GPR正演模拟问题。 相似文献
6.
《地球物理学报》2021,64(6)
常规探地雷达(GPR)逆时偏移大都未考虑高频电磁波在地下高电导率介质中传播的强衰减特性,致使高衰减区域的成像质量低;实际地下结构呈三维空间分布,二维GPR逆时偏移难以将反射波准确归位和绕射波完全收敛于真实位置.为此,本文构建了一种基于电磁波衰减补偿的三维GPR逆时偏移算法.该算法采用三维时域有限差分法计算正传和反传电磁波场,并通过改变反传电磁波方程中包含电导率的衰减项的正负号,人为保持反传时的时间对称性和反转不变性,以补偿正传时衰减的能量;零时刻成像条件用于获取地下三维结构的成像结果.数值试验结果表明:相比于常规GPR逆时偏移算法,基于电磁波衰减补偿的三维GPR逆时偏移可精确补偿电磁波在地下高电导率介质中正传衰减的能量,高电导率区域的成像精度更高,分辨率更好,抗干扰能力更强,其结果更有利于指导后续雷达剖面的解译. 相似文献
7.
基于双二次插值的探地雷达有限元数值模拟 总被引：3，自引：0，他引：3
从探地雷达(GPR)满足的波动方程出发,详细介绍了二维GPR模型单元剖分、二次插值、数值积分和有限元刚度矩阵总体合成的GPR有限元求解过程.为解决数值模拟时截断边界处的超强反射,采用Clay Bout透射边界条件对雷达波进行衰减,进而压制了来自截断边界处的反射波.在满足时间步长与空间网格差分稳定性前提下,采用中心差分法对GPR有限元方程进行离散,并用不完全LU分解预处理的BICGSTAB算法求解系数方程组,然后编制了基于双二次插值的GPR有限元正演模拟matlab程序.运用该程序分别对矩形和"V"字形两个典型地电模型进行正演计算,得到了正演剖面图,将该正演剖面图与基于线性插值的FEM算法的正演剖面图做了对比分析.结果表明基于双二次插值FEM算法相比基于双线性插值FEM算法异常响应更明显,具有更高的模拟精度,更有利于指导雷达剖面的数据解译. 相似文献
8.
探地雷达数值模拟中,时域有限差分法在时间和空间上一般采用二阶精度的中心差分近似(FDTD(2,2)),其形式简单,但数值色散误差较大,在复杂模型模拟时不能很好地反映模型的精细变化.高阶时域有限差分法能很好地改善数值色散带来的误差,提高模拟精度.本文基于三维高阶时域有限差分法的基本原理实现了探地雷达正演模拟,采用单轴各向异性完全匹配层(UPML)作为吸收边界条件,可以有效地吸收外向传播的电磁波,在大大地提高计算效率的同时,也能很好地改善边界的吸收效果.分析对比正演模拟结果,通过三维高阶时域有限差分正演能获得目标体准确电磁响应信息,并能很好的提高模拟精度. 相似文献
9.
应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷达正演方法,并利用该自制程序,对三角形金属体模型进行了正演模拟,得到了其相应的正演合成三维剖视图及切片图,通过对这些模拟结果进行分析,可以加深对三维雷达反射特征的认识,提高探地雷达探测的可靠性、准确度,同时也说明时域多分辨率法在探地雷达三维正演模拟中的有效性. 相似文献
10.
高频电磁波传播速度在水及淤积砂土中影响因素实验研究 总被引：2，自引：1，他引：1
水及淤积土中电磁波传播速度是实现水上探地雷达(GPR)探测数据准确时深转换的关键因素.论文基于探地雷达实测结果分析了水温、盐度、浊度对水中电磁波速的影响和粒径组成、含水量变化对砂土介电常数(电磁波速)的影响,建立了淤积砂土介电常数模型并给出了相应速度计算方法.研究表明水中电磁波速主要受水体盐度影响,随盐度增大而指数减小.淤积砂土介电常数符合Looyenga模型,现场探测时可根据土样三相体积比确定土体介电常数,进而确定土体电磁波速,实现GPR数据时深准确转换. 相似文献
11.
We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research. 相似文献
12.
Realistic numerical modeling of ground penetrating radar (GPR) using the finite-difference time-domain (FDTD) method could greatly benefit from the implementation of subgrids – supporting finer spatial resolution – into the conventional FDTD mesh. This is particularly important, when parts of the computational domain need to be modeled in detail or when there are features or regions in the overall computational mesh with values of high dielectric constant supporting propagation of waves at very short wavelengths. A scheme that simplifies the process of implementing these subgrids into the traditional FDTD method is presented. This scheme is based on the combination of the standard FDTD method and the unconditionally stable alternating-direction implicit (ADI) FDTD technique. Because ADI-FDTD is unconditionally stable its time-step can be set to any value that facilitates the accurate calculation of the fields. By doing so, the two grids can efficiently communicate information across their boundary without requiring to use a costly time-interpolation scheme. This paper discusses the performance of ADI-FDTD subgrids when implemented into the traditional FDTD method, using different communication schemes for the information exchange at the boundary of the two grids. The developed algorithm, can handle cases where the subgrid crosses dielectrically inhomogeneous media. In addition, results from the comparison between the proposed scheme and a commonly employed purely FDTD subgridding technique are presented. 相似文献
13.
复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟. 相似文献
14.
15.
在采用有限差分方法开展探地雷达复杂目标体精细结构模拟时,为了提高计算精度,常采用非均匀网格对目标区域划分小尺寸的网格,以压制离散网格频散现象和保证有限差分方法的稳定性.常规非均匀网格和自适应亚网格技术在网格剖分数量和粗细网格边界处理上难以达到计算效率和计算精度的均衡.本文根据隐形斗篷(invisible cloak)理论,将基于变换光学(Transformation optics)理论应用于有限差分探地雷达数值计算中.该理论的主要思想是基于目标参数变化而保持电磁场的传播不变性,在坐标变换后,Maxwell方程的形式可以维持不变,而使得相对介电常数与磁导率的表达式变得复杂.通过这种方式可以虚拟地扩大目标体所占的网格节点数,减少背景介质区域的网格数,不增加模型空间的网格总数.另外,这种网格划分方式不但提高了计算效率,同时也可以克服亚网格技术边界反射误差的影响.本文推导实现了基于变换光学的二维有限差分方法,通过典型探地雷达模型测试,对比分析了该方法与常规有限差分、变网格有限差分和自适应亚网格有限差分的优缺点.计算结果验证了基于变换光学的有限差分可用于探地雷达目标精细结构模拟,具有较高的计算精度和计算效率. 相似文献
16.
将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性. 相似文献
17.
本文基于声波方程的哈密尔顿系统,构造了一种新的保辛数值格式,简称NSPRK方法.该方法在时间上采用二阶辛可分Runge-Kutta方法,空间上采用近似解析离散算子进行离散逼近.针对本文发展的新方法,我们给出了NSPRK方法在一维和二维情况下的稳定性条件、一维数值频散关系以及二维数值误差,并在计算效率方面与传统辛格式和四阶LWC方法进行了比较.最后,我们将本文方法应用于声波在三层各向同性介质和异常体模型中的波传播数值模拟.数值结果表明,本文发展的NSPRK方法能有效压制粗网格或具有强间断情况下数值方法所存在的数值频散,从而极大地提高了计算效率,节省了计算机内存. 相似文献
18.
提出混合ADI-FDTD亚网格技术开展频散介质GPR正演,即在物性参数变化剧烈局部区域采用细网格剖分ADI-FDTD计算,其他的区域采用粗网格剖分常规FDTD计算,ADI-FDTD突破了CFL条件的限制,可选取与粗网格一致的大时间步长,有效地提高了计算效率.本文首先基于Debye方程,推导了粗网格FDTD及细网格ADI-FDTD频散介质差分格式,着重对粗细两种网格结合的场值交换方式进行了深入探讨,给出了该算法的计算流程.然后以一个薄层模型为例,分别应用粗网格、细网格、混合ADI-FDTD亚网格算法对该模型进行正演,计算资源的占用及模拟精度说明了混合ADI-FDTD亚网格算法的优势.最后,建立频散介质与非频散介质的组合模型,应用3种方法对该模型进行正演,对比3种方法优劣,分析雷达剖面中非频散介质及频散介质中波形特征,有效地指导雷达资料的精确解释. 相似文献