基于异常复电位2.5维CR有限元数值模拟

论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交流双点电源复电阻法的变分问题;然后利用有限单元法求解变分问题,在矩形网格中增加2条对角线形成三角网格对研究区域进行剖分,在单元中对复电位和岩、矿石复电导率进行线性插值,将变分方程化为线性代数方程组;最后,求解线性方程组,得到各节点的异常复电位值,并根据选定的测量装置计算出地表的视复电阻率,进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法的有效性后构建模型算例,给出了其不同测量装置下的计算结果。

基于Fractal模型的复电阻率法2.5D有限元数值模拟

复电阻率法在油气资源、矿产勘查中发挥了重要的作用,为了深刻认识复电阻率法异常特征变化规律,本文对复电阻率2.5D正演问题展开研究。首先直接给出复电阻率法2.5D有限元正演所满足的变分问题,并详细地推导相应的刚度矩阵的计算过程。引入Fractal模型作为等效模型研究频谱激电异常特征。对单元内的复电导率及复电位均进行线性插值,而后,采用不完全LU分解的稳定双共轭梯度算法求解有限元线性方程组,获得异常复电位值。设计3个典型的地电模型验证了本文算法的正确性及精确性,并分析了不同装置下,不同频率的2.5D复电阻率异常响应特征。数值模拟结果表明,采用Fractal模型研究激发极化异常特征是可行、有效的;不同装置、不同频率下的复电阻率法异常特征有着显著的差异。

基于GEMTIP模型的CR三维有限元数值模拟

本文提出了基于广义等效介质激电(GEMTIP)模型的复电阻率三维有限元数值模拟方法,首先介绍了GEMTIP模型及在不同参数影响下的复电阻率变化特征,而后给出双点电源的复电阻率法三维正演所满足的变分问题,采用六面体单元对研究区域进行网格剖分,在单元内对复电位、复电导率均进行线性插值,将变分方程转为线性方程组;采用不完全LU分解的稳定双共轭梯度算法(BICGSTAB算法)求解线性方程组,得到异常复电位值.设计一个两层极化大地模型,将数值解与解析解进行对比,验证了文中算法的正确性.设计两个典型地电模型的数值模拟结果表明,采用GEMTIP模型研究宏观尺度下的激电效应是有效的,可为进一步研究复电阻率异常响应特征提供一种新的途径.