三维高密度电阻率E SCAN法有限元模拟异常特征研究

工程电法勘探中，常需要探测埋深不大、规模较小、分布复杂的目标体，对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大，可对目标体进行多方位观测，是目前工程探测中常用的方法之一.本文对三维高密度电阻率E SCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型，利用有限单元法进行正演计算，并分析其异常特征和分布规律.算例结果表明，E SCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测，即可分辨地质体电性特征及水平位置，但观测结果难以分辨地质体纵向延伸.     

三维高密度电阻率ESCAN法有限元模拟异常特征研究

工程电法勘探中，常需要探测埋深不大、规模较小、分布复杂的目标体，对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大，可对目标体进行多方位观测，是目前工程探测中常用的方法之一.本文对三维高密度电阻率ESCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型，利用有限单元法进行正演计算，并分析其异常特征和分布规律.算例结果表明，ESCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测，即可分辨地质体电性特征及水平位置，但观测结果难以分辨地质体纵向延伸.     

三维电阻率测深数据Zohdy近似反演方法

Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率．为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算．在电阻率调整公式中加入一个迭代系数,能够加快收敛速度．并对加5％随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致．     

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

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

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.     

延庆台地电阻率铁质干扰的有限元模拟

近年来延庆台地电阻率观测多次出现因环境干扰而引起的较为显著的异常变化,尤其是观测环境中铁质干扰源对地电阻率观测的幅度和形态的影响.本文依据台站测区地下高密度电性结构剖面和电测深结果建立了三维数值模型,分析了几种铁质对延庆台地电阻率观测的影响幅度和形态,结果表明,在扣除地电阻率基准值后,干扰幅度和形态与实际观测值的异常变化较吻合,从而从数值分析的角度确定了地下铁质对延庆地电阻率观测的干扰.      

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

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

三维直流电阻率自然单元-无限元耦合数值模拟

直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难,本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型,并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程;然后介绍了二、三维自然单元法及其形函数的构建过程;再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合;最后通过4个数值模型,结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果,验证了算法的正确性及有效性.数值结果表明,本文所提算法满足精度要求,能有效减小计算区域,无需二次更新刚度矩阵,对复杂模型适应性强.     

基于非结构网格有限元的三维井地电阻率法任意各向异性正演模拟

地层介质的电各向异性增加了井地电阻率法响应的复杂性,开展基于电各向异性介质模型的井地电阻率法响应规律研究对于正确解释各向异性显著地区的观测数据特征至关重要.针对垂直线源井地电阻率法的任意各向异性响应模拟问题,本文提出了一种基于非结构网格有限元三维正演算法,通过引入3×3的对称正定张量来表征任意各向异性的电导率,采用非结构四面体网格有限元方法来离散电位的边值问题,通过将垂直线源等效为一系列点源问题,进而实现了任意各向异性介质中井地电阻率法的高效数值计算.通过与三个地电模型的解的对比,验证了本文数值解算法的精度和有效性.针对线源远离和垂直穿过异常体的两类模型,分别考察了当围岩或异常体为各向异性介质时的井地视电阻率响应特征.结果表明,对于各向异性地层,围岩和异常体的主轴电阻率值和旋转角均会对井地视电阻率的幅值及分布产生显著的影响.研究结果对于提高井地电阻率法的认识和资料解释水平具有重要的理论和实际意义.     

巷道超前探测的并行Monte Carlo方法及电阻率各向异性影响

直流电阻率法通过观测单极-双极装置视电阻率曲线最小值位置估计掌子面前方低阻异常体距离,在巷道超前探测及矿井水害预测与防治工作中有广泛应用,并形成多种经验预测模型.本文基于直流电法巷道超前探测的各向异性电阻率三维有限元数值模拟,给出超前探线性预测模型.然而,实际井下地质结构复杂,异常大小不定且任意分布,而且可能存在有电阻率各向异性,使得目前由简单模型实验或数值模拟获得的预测公式不确定性较大,而且可靠性难于评价.Monte Carlo方法用随机化思想解决非确定性问题,我们将该方法与并行算法结合,对大量不确定的随机模型进行三维数值模拟,以验证各预测模型的准确度及可靠性.结果表明,本文的超前探线性预测模型较其他预测模型更为准确及可靠.但介质各向异性的影响很大,尤其当围岩电阻率具有各向异性时,所有预测模型的准确度及可靠性较差,给安全生产带来很大隐患.     

矿井瞬变电磁法全空间视电阻率解释方法研究

从全空间中心回线源响应公式出发,推导了全空间晚期和全区视电阻率公式;通过与半空间公式比较和对实测资料的分析,研究了全空间视电阻率公式的应用效果.结果显示,全空间视电阻率为半空间视电阻率的(5/2)^2/3倍,其晚期视电阻率与真值的相对误差小于0.56%,而相应的半空间视电阻率与真值误差均大于42%.与常规解释方法相比,采用全空间全区视电阻率换算方法能更好地圈定低阻异常范围,且电阻率差异大.因此,该公式更适用于矿井瞬变电磁法资料处理与解释,并具有更好的效果.     

基于交叉梯度约束的电阻率法和背景噪声法三维联合反演研究

电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率.     

坑道直流电阻率超前聚焦探测新方法研究

研究和提出了一种坑道直流电阻率超前聚焦探测新方法.该方法通过在掌子面上设立不同功能的环状电极组,使一次场电流具有像探照灯一样的聚焦功能,可有效探测坑道掘进前方不良地质体的存在,达到超前预报的目的.本文首先介绍了直流电阻率超前聚焦探测的基本原理,并给出了两种聚焦观测方式.为了验证所提方法的可行性,通过轴对称电性介质二维异常电位有限元数值模拟方法,对聚焦观测条件下几例坑道模型进行了模拟计算.结果表明,该方法对异常反映明显,可进行即时解释,建议在坑道超前预报中逐步推广使用.     

宝昌地震台地电阻率电性结构特征

使用高密度电法和电测深法,对内蒙古宝昌地震台地电阻率场地进行实地测量,获得该场地电性结构特征参数,判定该台地电测区结构属于KH型,并结合钻孔柱状图,获得该台地下介质水平层状电性结构模型,发现地层厚度呈不均匀分布,总体较为连续、完整,可为今后该台地电阻率异常提取、异常识别及观测精度提高等提供基础数据。     

高密度电阻率法勘探长测线多排列数据连接处理

本文通过对高密度电阻率法勘探中的长测线多排列数据进行连接处理,实现了同一测线采用同一套参数的规范化处理,克服了因不同排列采用不同参数处理而导致假异常和信息丢失等弊端,同时解决了在成像时,因各个电性剖面的色标值不一致,剖面坐标与实际桩号对应不直观,导致对整个区域的电性趋势认识产生误导的问题.从实例处理来看,能为整个勘探区提取更丰富的电性信息,不但突出了整条测线的地下电性趋势,更客观、更全面地把握和认识勘探地区的电性分布规律,而且从宏观上把握了一些比较隐蔽的地质现象.     

Resistivity tomography study on samples with water-bearing structure

IntroductionThestudyonthevariationoftheeIectricaIpropertiesofsimulatingsamplesatcertainphysi-calconditionsinthelaboratoryisanimportantmethodtounderstandthecomposition,movementandevolutionofthematerialsintheseismicsourceregionduringseismogenicprocessandisaIsoanimPortantsupplementtothefieldexperimentofeathquakemonitoringbymeansofelectricaImethod.ThepreviouslymadesimuIationexperimentsindicatethatthevalueofapparentresistivityofwatCrsaturatedintactrocksamplesmaydecreaseoneorderofmagnitudebeforefr…     

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

金属大棚对周口台地电阻率干扰定量分析

周口地电阻率观测2015年出现异常变化,因异常期间该台场地存在较多干扰源。为正确评价该次异常,在电测深结果的基础上建立三维数值模型,采用有限元数值计算,定量评估该台测区周围金属蔬菜大棚等干扰源对地电阻率观测的影响量。结果表明:金属大棚可能是周口台本次地电阻率异常的主要影响因素。