大地电磁测深二维地形影响及其校正方法研究

采用有限元法实现ＭＴ二维模型的正演模拟,将地形条件下的二维地电模型的ＭＴ响应分解为地形引起的畸变场和与测点高程位置有关的稳定场两部分进行分析。用有限元法模拟二维地形条件下地电结构的ＭＴ 响应和纯地形的ＭＴ响应,采用比值法消除地形引起的干扰场,然后,将消除地形影响后的视电阻率值通过延拓方法获水平地形条件下二维地电结构的视电阻率分布,完善大地电磁测深视电阻率曲线的地形校正。     

高密度电阻率法地形影响校正

采用有限元法实现高密度电阻率法二维模型的正演模拟,将一定地形条件下的二维地电模型的视电阻率响应,分解为畸变场和正常场二部分进行分析.采用比值法消除由地形引起的畸变场,从而获得水平地形条件下二维地电结构的视电阻率值.在实际资料的校正中,取得了较好的应用效果.     

时间域电法测量地形响应研究

这里以对称四极装置为例,采用有限元法实现各种地形地电条件的时间域二维电阻率正演数值模拟,运用比值法原理提取地形响应参数,通过对比分析不同电性条件和不同地形条件下的地形响应曲线,得到地形响应参数随地形与介质电阻率变化的一系列规律,并进一步完善了比值法地形校正方法。     

基于对偶加权后验误差估计的2.5维直流电阻率自适应有限元正演

采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。     

甚低频电磁法边界元数值模拟及地形影响与改正

本文提出了二维地电模型甚低频电磁法中磁场、电场、极化椭圆倾角、视电阻率异常的边界元数值模拟方法。并根据计算结果讨论了地形影响基本规律与改正方法。所编制的计算程序亦可适用于二维非水平层上满足远区场条件的各类频率电磁测深法。     

AMT数据地形畸变影响与带地形反演效果研究

在山区进行音频大地电磁测深(AMT)时,电磁场易受地形起伏影响发生畸变。为研究起伏地形对AMT数据的影响,采用WinGlink软件基于有限差分算法的正演模块,建立二维起伏地形和平地地电模型,并对比分析了两个理论模型的正演响应。结果表明:TE模式的视电阻率曲线在高频段受干扰程度相对较大,在低频段影响减小甚至消失;TM模式的视电阻曲线变化正好相反;二者相位受地形影响则较小。对含异常体的带地形地电模型进行正演计算,在正演响应中加入2.5%的随机误差模拟实测数据,利用WinGlink软件基于非线性共轭梯度(NLCG)算法的反演模块,分别进行带地形和不带地形的TE、TM和TE+TM模式的二维反演。各反演结果的对比表明,带地形的二维反演结果明显优于不带地形的反演结果,能够较好地反映出异常体的位置、形态及电阻率值。依据得到的认识,对青海祁连隧道工程的AMT实测数据进行带地形二维反演,取得了理想的应用效果。     

点源二维有限元法在三电位电极系正演模拟中的应用

利用点源二维有限元法,针对二维地电模型,进行了三电位电极系视电阻率参数及比值参数的计算,对几种二维地电断面三电位电极系的视电阻率及比值参数的异常规律进行了探讨。这些研究结果,为实际勘探中复杂地电断面的异常解释打下了一定的理论基础,在利用比值参数进一步分辨复杂断面的电性异常方面作了有益的尝试。     

基于二次场二维起伏地形MT有限元数值模拟

通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。     

基于二次插值的线源可控源有限元数值模拟

在准静态近似条件下,采用矩形网格单元和双二次函数插值就频率域二维线源边值问题进行了有限元数值模拟。在二维地电条件下,给出了边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成,最后通过求解复系数方程组得到了地表视电阻率响应。引入伪delta函数模拟线源,消弱了源带来的奇异性。通过与均匀大地以及层状介质模型的解析解对比,平均相对误差分别为0.71%和1.12%。建立了两个异常体模型,数值模拟表明异常响应比较明显,为进一步实现三维可控源电磁法有限元数值模拟提供了基础。     

任意三维地电模型视电阻率异常的边界元法数值解

边界元法(BEM)是一种新的数值方法。由于该法可以降低所研究场问题的维数,因而边界元法较之域型法(FEM和FDM)具有应用简便、数据量少、计算快、精度高等优点。从而,用该法解决了域型法难以实现的三维地电模型视电阻率异常的计算问题。本文论述了用边界无法求解点源场地表水平和起伏下三维地电体位场问题的方法原理和数值处理技术,并给出了若干算例:导电球状矿体上视电阻率数值解与解析解结果;导电球状矿体上不同测线上视电阻率平剖曲线;三维山脊地形下导电椭球状矿体上视电阻率联剖曲线及其地形改正结果等。由本文内视电阻率的边界元法数值解与解析解结果对比的一致性和三维地形校正的显著效果,表明了边界元法是求解任意三维地电模型上位场问题的有效方法。由于用该法实现了对任意三线地质体上空间位场计算,必将推动三维电法勘探工作的开展与深入研究,无疑会对提高电法勘探的地质效果发挥重要作用。     

用边界元法计算点源二维地电体ρs与ηs

边界元法是一种新的微分方程数值解法.本文讨论了用该法求解点电源场中任意起伏地形下二维地电体位场问题的方法原理和数值处理技术,并给出了不同的地电体视电阻率与视极化率的计算结果相应用实例. 结果表明:边界元法是求解电法勘探边值问题和资料处理的一种经济有效的数值计算方法.     

水平地形二维地电断面点源电场的正演

本文给出水平地形下二维地电断面的电阻率及面激发极化异常的数值计算方法。首先用镜象法将半空间问题变为全空间问题,然后用付氏变换将点源二维地电断面的三维电场边值问题变换成二维边值问题。再用格林公式将二维边值问题转变成积分方程。用边界单元法解积分方程,并由付氏反变换,获得三维空间中的电位。用本方法计算了几种模型的电阻率异常及激发极化异常,并与模型实验对比,效果是良好的。与起伏地形二维地电断面的边界单元解法比较,本文所需的计算机内存减少很多,计算时间节省几倍。     

大倾角地形下视电阻率数值模拟

为解决大角度起伏地形即陡峭地形(如河流阶地、冲沟、堤坝、山谷、山脊等)对实测视电阻率的影响问题,这里优化了三维有限元数值模拟中的单元积分形式,使之适合大倾角地形条件下对电阻率的数值模拟.为验证方法的正确性,采用保角变换计算了相近地形在均匀电流场情况下,陡峭地形引起的视电阻率变化特征.计算表明,对于缓倾斜地形,有限元数值...     

CSAMT 2.5维有限元数值模拟

对三维场源二维地电模型的正演计算称为2.5维数值模拟。这里从麦克斯韦方程组出发,1分别求解电磁场的一次场和二次场,将三维场源降为二维;2利用傅氏变换将空间域方程转化为波数域,应用有限元求解波数域电磁场方程,引入无限元解决无穷远边界的收敛问题;3根据电磁场值实虚部的曲线特征,按对数等间隔选取21个波数。编写代码计算均匀半空间与解析解结果对比,电阻率的均方相对误差均小于0.5%,证明其有效性。计算了三种地电模型的电磁场响应,对单一高、低阻体和高低阻组合体的模拟效果真实,异常中心位置基本吻合。结果证明,这里的CSAMT2.5-D正演算法可以模拟较为复杂的地电模型,并取得良好的效果。     

基于非结构化网格的大地电磁2D自适应有限元正演模拟

有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。     

线电流源声频大地电磁测深的二维正演计算及响应特点

论述了如何在大地电磁测深二维正演计算的基础上，结合线电流源声频大地电磁测深（ＣＳＡＭＴ）的边界条件，利用三角形６节点有限元法计算ＴＥ极化方式ＣＳＡＭＴ二维响应的方法；给出了４个模型的计算结果     

三维张量CSAMT响应特征研究

张量CSAMT(可控源音频大地电磁法)能完整地描述三维的地电信息,更适应复杂地质条件下的地质勘探。采用将总场分解为一次场和二次场的策略,用差分近似微分形成大型线性方程组,通过解方程组得到张量CSAMT三维正演结果。通过计算低阻和高阻目标体的响应,发现卡尼亚视电阻率和视阻抗相位对低阻体表现为低值异常,对高阻体表现为高值异常,视阻抗相位达到极值时的频率比视电阻率达到极值时的频率要高。张量CSAMT对低阻体比高阻体更灵敏,在同样条件下,高阻体的异常区位置相比于低阻体会上移。某工区的勘察案例表明,根据卡尼亚视电阻率异常能定性地判断地下异常体的存在,为解决地质类问题提供了依据。本研究为张量CSAMT三维反演提供了基础。     

利用有限单元法模拟二维MT倾子响应

在二维地电模型中,考虑到大地电磁测深(MT)倾子响应依赖于横电(TE)模式的数值模拟问题,笔者采用矩形网格单元和双二次插值对MT的倾子响应进行求解。首先,给出了二维地电条件下的边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成;其次,通过求解复系数方程组得到了每个节点的电场值;最后,采用差分方法求得电场沿横向和垂向的偏导数值,从而计算出二维MT的倾子响应。通过对2个二维地电模型的倾子响应进行数值模拟,获得了倾子响应的实部、虚部和振幅,结果很好地反映了不均匀体的横向分布情况。     

测量磁场水平分量H y的电性源广域电磁测深法

基于层状大地表面水平电偶极源电磁场公式系统,提出了比值法获取H y分量全区视电阻率计算公式。计算了层状地电模型多个收发距的电磁场各场分量E x、H y、H z,并用这三个场分量计算了全区视电阻率,以及E x/H y卡尼亚电阻率,并与模型的MT卡尼亚电阻率频谱曲线进行对比。发现H y全区视电阻率在赤道装置和轴向装置的广大区域对地电模型均有较好的响应,可进行单分量广域电磁测深。轴向装置H y全区视电阻率与赤道装置垂直磁场H z全区视电阻率相似,在低频段与层状模型大地电磁测深卡尼亚电阻率有相似的频率响应特征,在较小收发距条件下对大埋深基底就能有较好的响应。赤道装置和轴向装置测量的水平电场E x分量全区视电阻率在低频区均会进入“饱和”区。对于各电阻率定义方法,合适的收发距是较好地反映出地电特征深度变化的前提,多收发距测量有利于揭示深部电性特征。     

张量CSAMT数据的二维反演

在一单行本中,我们提出了一种ＣＳＡＭＴ正演模拟算法并说明了这种算法在解释野外数据中的应用效果,在本文中,提出了一种用于二维电阻率模型的ＣＳＡＭＴ数据反演算法,这种算法包含一个三维源并可以反演（ＴＥ）模式及（ＴＭ）模式下的数据,用这种算法对假设的数据进行试算,结果表明,不用做平面波校正就能对远区场,过渡区场及近区场进行反演,最终的模型结果显示,原始电阻率模型可被恢复,而不存在任何由上源于效应而产生的