复电阻率法二维有限元数值模拟

伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映.     

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.     

三维地电模型数值模拟中视电阻率真假异常特征分析

本文概要介绍了三维数值模拟中双网格系统的剖分方法,并通过几种典型模型的有限单元法三维数值模拟计算（网格单元电导率分块均匀、电位三线性变化）,充分揭示了对称四极电阻率测深法三维地电模型视电阻率真假异常的分布规律,分析了产生真假异常的原因,并研究了真异常的形态与引起真异常的地质体位置的关系。     

点源两层大地三维地电模型视电阻率边界元解

本文用两层均匀大地人工点源作用下的空间电位分布作为基本解，导出了下伏基岩中陷含三维地质体情况下，关于电位的边界积分方程组，用边界单元法解积分方程求得地表电位，从而求得了一定装置下的视电阻率。本文方法的独特优点是积分只须在异常体边界面上进行了测线可以沿任意方向。用联剖法计算了均匀半空间下含球形异常体的ρｓ曲线，所得结果与其他学者用解析方法所求得的结果一致。用中间梯度法计算下层基岩中含有一长方体，测线     

复杂地形三维直流电阻率有限元数值模拟

系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础.     

电导率连续变化2.5维直流电阻率法有限元数值模拟

在阮百尧、徐世浙(1998,2001)工作基础上,给出四边形内三角网格单元剖分,电导率连续变化2.5维直流电阻率数值模拟方法.在二维起伏地形情况下,有限单元设计为四边形内剖分四个三角网格单元,单元内的电位值和电导率设计为线性变化,将研究区域远边界的第三类边界条件简化为第二类齐次边界条件,重新给出一组适用性较好的Fourier反变换系数,并采用分离式Cholesky分解法快速求解线性方程组.通过对3个地电模型的计算与对比分析,其结果与解析法的均方根误差小于0.2%,计算速度比较常规方法提高了5倍,能有效模拟起伏地形和复杂形态的地电体模型.     

电导率连续变化的三维电阻率测深有限元模拟

针对地下岩矿石的物性参数常常是连续变化的情况,给出了电导率分块连续变化的三维地电断面电阻率测深的有限元数值模拟方法. 首先列出了三维地电断面点源电场的变分问题;然后利用有限单元法解变分问题,采用六面体单元对研究区域进行剖分,对电导率及电位均用三线性函数进行插值;并对2例模型进行了计算,将其结果与二维有限单元法的情况进行了对比,提高了精度.     

时间域瞬变电磁法全区视电阻率的平移算法

深入分析了均匀半空间中瞬变电磁场的解析表达式,发现瞬变响应曲线随地下电导率、发射回线边长与观测时间具有平移伸缩特性,并据此提出了一种新的直接计算全区视电阻率的方法.该法与现有计算全区视电阻率的方法的不同点在于：无需迭代和求解非线性方程,速度快,精度高,算法稳定,且只有惟一解.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

裂隙非均匀分布岩石中的电位及视电阻率计算

计算了在均匀高阻介质中含有低阻裂隙的岩石的电场分布，并由此推导出裂隙非均匀分布的岩石中视电阻率的计算公式．结果表明，裂隙的埋深、规模、走向（平行或垂直视电阻率测线）以及真电阻率等，对视电阻率的观测和计算结果都有不同程度的影响．本文结果对视电阻率法用于地震分析预报提供了一种有效的计算方法．