基于双二次插值的探地雷达有限元数值模拟

从探地雷达(GPR)满足的波动方程出发,详细介绍了二维GPR模型单元剖分、二次插值、数值积分和有限元刚度矩阵总体合成的GPR有限元求解过程.为解决数值模拟时截断边界处的超强反射,采用Clay Bout透射边界条件对雷达波进行衰减,进而压制了来自截断边界处的反射波.在满足时间步长与空间网格差分稳定性前提下,采用中心差分法对GPR有限元方程进行离散,并用不完全LU分解预处理的BICGSTAB算法求解系数方程组,然后编制了基于双二次插值的GPR有限元正演模拟matlab程序.运用该程序分别对矩形和"V"字形两个典型地电模型进行正演计算,得到了正演剖面图,将该正演剖面图与基于线性插值的FEM算法的正演剖面图做了对比分析.结果表明基于双二次插值FEM算法相比基于双线性插值FEM算法异常响应更明显,具有更高的模拟精度,更有利于指导雷达剖面的数据解译.     

井地电位法中变物性参数异常响应的模拟与分析

为了研究地下异常体几何及物性参数变化对井地电位法地表响应的影响规律,对垂直线源激励下的变物性参数地下异常体,采用三维有限单元法进行了正演模拟。经模拟结果表明,在浅部时异常形态清晰,随深度的增加,异常响应的形态变得模糊;地表异常响应对厚度反应不灵敏;对水平规模大小变化反应非常灵敏。这些规律将对油气田剩余油分布和注水推进前沿的确定,以及井地电位仪器分辨率的设计提供指导作用,并为井地电位法反演成像奠定相关基础。     

基于速度估计的改进型线性变换有限差分偏移在探地雷达中的应用

为了解决传统线性变换有限差分偏移算法对变速介质的高敏感性问题,对加入了速度变量的标量波动方程进行线性变换,在横向非均匀介质条件下,详细推导了包含混合因子项的线性变换差分偏移方程,实现了横向非均匀介质中陡倾角反射界面的正确归位,同时改进型算法可以灵活地选择混合因子项,这样就增强了原算法对复杂地质问题的适应能力,在此基础上...