三维地电断面激发极化法有限元数值模拟

用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.      

频率域激发极化法有限元数值模拟

在频率较低和忽略电磁效应的情况下,利用有限单元方法和Cole-Cole模型对频率域激发极化法进行数值模拟.首先在三维地电条件下,给出电场的边值问题和变分问题,将Cole-Cole模型的频率响应引入到地电模型中,运用有限单元法对模型进行单元剖分、插值、积分和总体合成等,通过解方程最后得到表征频率域极化强度的参数幅频率.通过改变地电模型的参数,得到了不同的幅频率响应曲线.模拟结果与实际情况符合,表明该方法是正确和适用的.     

激电二维数值模拟研究

从点源二维地电问题出发,采用有限单元法进行了地电场进行数值模拟,采用自适应三角剖分来实现起伏地表的模拟,针对双边三极装置,实现了多种模型的正演研究;通过多种模型的正反演,总结异常产生的规律,为激电法的分析提供了有效的信息。     

井—地三维激发极化数据快速反演

为实现井—地激发极化数据的快速反演,分析三维有限元正演四面体网格剖分形成的系数矩阵的元素规律,应用MSR非零元素压缩存储和SSOR-PCG方程求解技术,提高了正演执行效率。在正演基础上,利用Jacobian矩阵与正演方程的简单关系,采用共轭梯度法求解最小二乘目标函数的模型修改量,避开Jacobian矩阵的直接计算和存储,实现了三维井-地激电数据的分钟级快速反演。模型正演合成数据的反演结果表明,能基本再现已知的模型。     

四川那朗地区多金属矿化地质及激电特征

四川那朗地区地处三江成矿带,区内出露三叠系灰岩、板岩及火山岩,有闪长岩岩株侵入其中,地表见铁帽。岩体蚀变分带明显。通过激发极化法工作,由视极化率处理后异常显示,原始数据圈出了激电异常,推测为铁帽下部之矿（化）致异常。     

复杂地形条件下激发极化有限单元法三维数值模拟

野外激电勘查过程中,据激电异常进行钻探验证时,常发生验证目标矿体擦边而过现象.究其原因是在实际勘查中,地形的起伏变化改变了地下电场的分布规律,导致激电异常发生位移及幅值改变.为了改善、克服这种现象,必须从理论上进行研究.运用四面体单元三维有限单元数值正演计算方法对此进行了理论研究,选用典型山谷地形水平层状极化体模型进行数值模拟计算.研究结果表明:①起伏地形不产生假的激电异常,仅对激电异常的形态和幅值产生影响;②采用四面体剖分单元有限单元法研究复杂地形条件对激电异常的影响是行之有效的.     

层状介质背景下地-井方位激电异常分析与解释

为了更好地探究层状介质背景下井旁隐伏矿体的地-井方位激电异常特征,首先研究了地-井三维激电的有限元正演模拟算法,并编制了程序;然后构建三维地电模型并进行正演模拟,讨论层状介质背景下地-井方位激电异常,并分析其异常特征规律;最后通过比值法归一化压制层状介质背景,突出局部异常。理论和算例表明,比值法归一化处理是地-井方位激电观测中消除层状介质背景异常影响的有效方法。     

多极距中梯观测与反演研究

针对传统中梯激电效率高、数据量少等特点,提出多极距中梯观测技术,配合现有的多通道发射、接收机,可在不增加太多工作量的情况下,尽可能采集更多的视电阻率和视极化率数据,以获取更多的地电信息,进行二维反演,实现中梯断面成像.模型正演合成和初步实测数据的反演表明,多极距中梯观测与反演的结果可基本等效于对称四极测深效果,但工作效率成倍提高,可推广应用.     

基于有限单元法的激电参数正演数值模拟

通过利用有限单元法模拟频率域激电法中的重要参数幅频率Fs和相频率Фs,研究Fs、Фs在三维地电条件下的激电响应特征。首先,在有限元方法中,使用六面体和四面体相结合的剖分单元,对大地模型进行离散,通过计算不同频率下的Cole—Cole模型复电阻率响应,替换模型中的体极化模型,最后模拟得到了频率域激电法中的主要参数。通过对比分析,证明了该方法程序简便、准确,也证明了相频率参数可作为一种有效区分时间常数不同的极化体方法。     

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

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