时移电阻率反演模拟研究

电阻率法可用于地表浅层的探测,也可用于对动态地下目标进行监测.对于监测数据的反演,不同数据集的单独反演存在一定缺陷,为此本文在常规电阻率反演算法的基础上,推导了时移电阻率反演公式,实现了时移反演算法程序;为了论证时移反演算法对动态地下目标成像的优越性,建立一组多个正演模型,利用模拟数据进行单独反演和时移反演,并对比两种...     

CSAMT一维全频视电阻率反演

本文采用全频率视电阻率反演,使用2~(-3)~2~(13)Hz频率,并且不做近区和过渡区校正,直接用卡尼亚视电阻率作为反演参数。正演采用虚界面法计算有限长导线在水平层状介质的电磁场分量,反演采用有限内存拟牛顿法。数值模拟采用三层模型和实测数据进行反演,得到了符合地下电性变化的反演模型结果,证实了CSAMT一维全频视电阻率反演具有可实施性。     

高密度电阻率法四电极系横向分辨率探讨

用高密度电阻率法探测地下目标地质体,已经成为近几年常用的方法。随着电法仪器的发展,野外工作方法的改进,以及反演软件的不断更新,高密度电阻率法分辨率已有了长足的提高。这里针对目前市场上常用的RES2DMOD电法模型正演软件和RES2DINV反演软件,对高密度电阻率法四电极系横向分辨率进行了初步探讨,并对野外工作方法选择和反演解释提出部份建议。     

跨孔电阻率法有限元数值模拟及反演

油气资源和矿产资源勘探中,常需要较大的勘探深度和较高的深部分辨能力.充分利用现有钻孔的井地或跨孔电阻率成像技术适应这种要求,但引起的视电阻率响应复杂,从而给资料的正确处理与解释带来一定的困难,为此有必要通过正演数值模拟和反演相结合的手段,对理论模型井地或跨孔装置的响应异常规律进行研究.文章通过设计一个典型地质模型,利用Ansys有限元软件对模型进行2D跨孔电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律.     

起伏地形下复电阻率法25维反演研究

复电阻率法以岩、矿石的频谱参数特性为勘探依据,多种参数的组合解释能为评价异常源的性质提供更为丰富的信息。但其反演问题一直没有得到较好地解决,已经严重地影响和制约了该方法的应用和发展。本文由频率域电磁场的偏微分方程出发,结合Cole-Cole数学模型进行公式推导,并利用有限单元法实现了复电阻率2.5维正演模拟,能够计算起伏地形条件下包含电磁效应的复电阻率响应。在此基础上,提出了利用多个排列的视电阻率和视相位数据的复电阻率反演方法,并在反演方程中加入了Occam法的光滑模型约束,能够对复电阻率参数直接进行大规模的全区反演,有效减弱了SIP反演中欠定性严重的问题。通过理论模型的反演算例表明,该方法能够快速、准确地反演出二维地质断面上所有单元的4种复电阻率参数,从而验证了该反演算法的正确性与稳定性。     

时移可控源音频大地电磁法三维反演研究

时移地球物理监测是观测地下物性参数随时间动态变化的有效方法。地球物理反演结果会受到测量误差、测量环境的噪声污染以及地球物理反演多解性等系统性因素的影响。如果对时移地球物理数据进行不同时刻数据的单独反演,可能会导致不同时刻的反演结果可对比性差,从而影响时移地球物理的监测效果。针对可控源音频大地电磁法在监测领域的实际需求,将不同时刻的观测数据放在一起反演,不同时刻的电阻率模型相互约束,实现了时移可控源音频大地电磁法三维反演。设计三种类型模型进行合成数据的单独反演和时移反演对比试算,验证了时移反演算法的有效性。反演结果表明：通过加入相邻时刻模型约束,时移反演可以更好地聚焦异常体的位置,降低噪声以及测量环境变化对反演结果的影响,防止反演伪影掩盖真实的地下变化。研究成果为时移可控源音频大地电磁法的实际应用奠定了良好基础。     

瞬变电磁全区视电阻率数值计算方法研究

基于均匀半无限介质空间中场的表达式,给出了一维层状介质中心回线装置瞬变电磁测深法的全区视电阻率定义,研究了三种全区视电阻率的数值计算方法。特别是在计算的过程中,利用分段计算,解决了全区视电阻率解的非唯一性问题。同时引入二分查找法,解决了隐函数的反函数求解问题。通过正演建模和烟圈反演,进一步对比分析了三种计算方法,这在瞬变电磁测深资料的定性解释中,正确反映地下一维电性层电性垂向的变化特征具有重要的意义。     

高密度电阻率法的2.5维反演

讨论了电阻率法2.5维正演和反演的算法,并在此基础上编制了高密度电阻率法2.5维反演程序,该程序可用于8种常用电极装置观测结果的反演.对理论和实测数据的反演结果表明,采用的算法正确,程序运行稳定,反演效果很好.     

基于高密度电阻率法对季节性冻土的监测研究

通过建立正演冻土模型,验证了将高密度电阻率法应用于冻土监测的可行性。对高密度电阻率法装置进行特殊处理,将其应用于对季节性冻土的监测,并进行了七个月份的连续监测。结果表明:高密度电阻率法在季节性冻土监测中可行;通过对实测数据进行反演处理,得到电阻率反演剖面,可以看出表层土壤在季节变化中有明显的结冻、解冻行为;并可估算出冬季冻土深度范围。     

直流电阻率测深中二维与三维反演结果的对比与分析

直流电阻率测深中,用二维反演程序对三维地质体进行了反演,并与三维反演结果进行对比和分析。首先,对比二维和三维最小二乘反演在正演模拟算法和先验信息的确定方法;然后,对若干比较典型的模型进行反演实例对比。由于二维程序反演仅考虑单个剖面的电阻率信息,无论是在异常位置、形态及电阻率特性上,其反演精度都比较低;三维电阻率反演综合了多个测深剖面的电阻率信息,其结果与实际模型吻合得非常好。     

直流电阻率三维正演的预条件共轭梯度法研究

三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。     

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

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

水下直流电阻率法数值模拟

利用有限单元法对水下电阻率测深进行数值模拟。与地面测量不同的是，水下电极所处的深度、水体本身和水底地形都会给电阻率数据造成较大影响，给数据的解释带来困难。若先通过计算，分析不同电阻率装置在水下测量时的视电阻率与水深的关系，再用有限单元法对水下的三维地电断面进行数值模拟，讨论水底地形对视电阻率的影响．即可用比较法消除地形影响。     

电阻率图像重建的佐迪(Zohdy)反演及其应用效果

本文讨论了利用佐迪反演方法对地表高密度电阻率法的数据采集结果进行电阻率图像重建。通过对多种复杂地电模型的数值模拟和野外实例的反演分析,表明该方法简单易行,计算成本较低,成像的分辨率高,效果较理想,将其用于地表测量的电阻率层析成像完全可行。作者用Fortran语言编制了实现上述过程的程序,对所有数学模型的正演计算均采用2.5维有限元法。为突出成像效果,本文在结果图示时均采用了伪彩色图像处理技术。     

不同电阻率测深方法对旁侧不均匀体的反映

利用三维有限元数值模拟,重点研究了旁侧三维不均匀体所产生的视电阻率的影响问题,展示了4种常规电阻率测深法在旁侧不均匀体上的异常分布特点;探讨了测深点正下方目标体异常与旁侧不均匀体异常的区别问题,指出旁侧不均匀体的影响是产生电阻率法假异常的主要原因之一.     

水平电偶源r方向电场值应用的理论研究

根据准静态极限条件下,水平电偶极源在均匀半空间表面产生的r方向电场分量E r表达式,分别在一维和三维条件下,对E r进行数值模拟并迭代计算出相应视电阻率,讨论E r的观测范围和应用优势。在一维条件下,通过对四种典型层状模型(两层D、G,三层H、K)的数值计算,并与相同条件下的E x视电阻率进行对比,结果表明:E r比E x受测量方位的影响更小;在三维条件下,选取均匀半空间存在异常体的模型,采用积分方程法对E r进行数值计算,与相同条件下的E x视电阻率进行对比,结果表明:当均匀半空间存在低阻异常体时,E r对低阻体y方向边界的反应比E x效果较好。     

隧道含水构造直流电阻率法超前探测研究

含导水构造和岩溶裂隙水的探查是亟需研究和解决的关键问题。直流电阻率法以其对低阻体反映敏感、分辨率高等优点被引入到隧道掌子面前方含水构造的预报工作中。首先,推导出了无限全空间中球体含水构造直流电法超前探测的解析解公式,并在利用有限元数值方法进行正演分析的基础上,借助“比较法”将含隧道腔体的三维空间中超前探测视电阻率数据转换为无限全空间中的数据,并探讨了由MN电极附近电阻率不均匀所造成的干扰以及识别和去除这种干扰的方法。然后,设计了基于全空间球体含水构造超前探测解析公式的阻尼最小二乘反演方法,对隧道中多组球体含水构造模型进行了较为成功的反演。另外,探讨了其他形状的有限体积的含水构造超前探测异常规律,提出了“等效球体”的概念,亦在全空间球体含水构造超前探测解析公式的基础上设计了快速反演方法。最后,利用物理模型试验验证了有限元数值正演方法和阻尼最小二乘反演技术的可行性和有效性,为实际工程中含水构造的直流电法超前预报奠定了基础。     

成像区外异常体对电阻率层析成像结果的影响

作者研究了成像区外不均匀体对电阻率层析成像结果的影响、其结果表明：成像区外非均匀体的存在对成像结果的影响是明显的;在用地面和钻孔中采集的数据成像时这种影响表现为产生一个相对于钻孔呈大致对称关系的假像;这种影响着区外非均匀体远离成像区而减小,。当区外非均匀体与距其较近的一口井的距离大于井的一半时,这种影响即可忽略。     

三维地电断面电阻率测深有限元数值模拟

用有限单元法进行了电导率分块均匀的三维点源电场电阻率测深的数值模拟.首先给出了三维构造中点源电场的边值问题、变分问题; 然后, 用有限单元法求解变分问题, 将区域剖分成六面体单元, 在单元中进行三线性函数插值, 将变分方程化为线性代数方程组; 最后解方程组, 得各节点的电位值, 进而计算出地表的视电阻率.对几例较典型的地电模型进行试算, 结果表明本方法是行之有效的