共查询到20条相似文献，搜索用时 562 毫秒
1.
2.
起伏地形下复电阻率法25维反演研究 总被引：1，自引：0，他引：1
复电阻率法以岩、矿石的频谱参数特性为勘探依据,多种参数的组合解释能为评价异常源的性质提供更为丰富的信息。但其反演问题一直没有得到较好地解决,已经严重地影响和制约了该方法的应用和发展。本文由频率域电磁场的偏微分方程出发,结合Cole-Cole数学模型进行公式推导,并利用有限单元法实现了复电阻率2.5维正演模拟,能够计算起伏地形条件下包含电磁效应的复电阻率响应。在此基础上,提出了利用多个排列的视电阻率和视相位数据的复电阻率反演方法,并在反演方程中加入了Occam法的光滑模型约束,能够对复电阻率参数直接进行大规模的全区反演,有效减弱了SIP反演中欠定性严重的问题。通过理论模型的反演算例表明,该方法能够快速、准确地反演出二维地质断面上所有单元的4种复电阻率参数,从而验证了该反演算法的正确性与稳定性。 相似文献
3.
复电阻率法在油气资源、矿产勘查中发挥了重要的作用,为了深刻认识复电阻率法异常特征变化规律,本文对复电阻率2.5D正演问题展开研究。首先直接给出复电阻率法2.5D有限元正演所满足的变分问题,并详细地推导相应的刚度矩阵的计算过程。引入Fractal模型作为等效模型研究频谱激电异常特征。对单元内的复电导率及复电位均进行线性插值,而后,采用不完全LU分解的稳定双共轭梯度算法求解有限元线性方程组,获得异常复电位值。设计3个典型的地电模型验证了本文算法的正确性及精确性,并分析了不同装置下,不同频率的2.5D复电阻率异常响应特征。数值模拟结果表明,采用Fractal模型研究激发极化异常特征是可行、有效的;不同装置、不同频率下的复电阻率法异常特征有着显著的差异。 相似文献
4.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献
5.
受激发极化效应的影响,大地中岩(矿)石电阻率是一个与频率、电化学性质等因素有关的复电阻率。接收到的电磁场信号中,包含了电磁感应效应以及激发极化效应双重因素的影响。将cole-cole复电阻率模型引入二维线源可控源正演理论,加载伪delta函数,实现了含激发极化效应的线源可控源有限元数值模拟。计算了不同极化参数的二维地质体线源可控源的影响特征。结果表明:线源可控源响应对IP效应反映明显,视电阻率出现增高的特征;低阻高极化相对于高阻高极化模型受IP效应影响更大,从低阻中提取IP效应更加乐观。该研究为实现从复杂地质条件下CSAMT信号提取IP效应奠定了基础。 相似文献
6.
7.
多年来很多学者开展了复电阻率电磁场正演模拟研究,然而大多数都是针对电性源进行的。鉴于此,这里采用体积分方程法进行了磁源激发条件下的复电阻率三维正演模拟。通过分析正演结果,得出了磁源激发极化法可行的结论,且磁场响应是由电磁感应和激电效应共同产生的,并引入磁荷概念总结了这两种效应的电磁响应规律。 相似文献
8.
在地震勘探中,地震正演模拟是非常重要的技术。与时间域正演相比,频率域正演速度快,计算效率高。如何高效准确地完成频率域正演计算是目前该领域的一个重要问题。数值频散问题和如何提高计算效率降低求解分解阻抗内存占用量一直是频率域正演所需要解决的问题。与传统的直接法求解阻抗矩阵的频率域正演方法不同,本文采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演,并提出了一种波数补偿的表达式来压制数值频散现象。经过简单模型和复杂模型的数值测试,采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演能够有效降低计算时间,且采用波数补偿的频率域正演方法能够有效压制数值频散现象,提高波场模拟精度。 相似文献
9.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
10.
电阻率/激电测深二维人机交互正演模拟 总被引：1，自引：0，他引：1
在有限元点源二维电阻率正演基础上,针对传统对称四极装置其电极距跨度大且不等间隔分布等特点,提出将整个正演计算区域划分为目标区、观测区和扩展区三个区域,在观测区域先采用等间隔剖分,然后在区域中插入实际电极作为正演计算节点的剖分方式,同时将无穷远边界条件简化技术引入到正演计算中,并结合人机交互与可视化技术实现了大极距电阻率/激电测深的二维正演模拟。最后通过对相关模型的试算结果表明,作者在本文中提出的针对传统对称四极装置电阻率/激电测深的二维正演模拟技术是有效可行的。 相似文献
11.
12.
从点源场的边值问题出发,利用有限元推导了井地电法和高密度电法2.5维的详细解法,并编制了相应的程序。利用均匀半空间的地电模型进行了正演计算,并与地下真实电阻率进行了对比分析,验证了算法的正确性。在此基础上对均匀介质中存在对称高阻体和边缘高阻体的模型,进行正演模拟,分别得到其井地电法的视电阻率断面图和高密度电法的视电阻率断面图。结果表明,结合两种方法可以对地下存在的高阻异常体进行有效精确地探测。 相似文献
13.
14.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献
15.
积分方程法三维模拟井地电法并行算法研究 总被引：1，自引：1，他引：1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
16.
基于二次插值的线源可控源有限元数值模拟 总被引：3，自引：0，他引：3
在准静态近似条件下,采用矩形网格单元和双二次函数插值就频率域二维线源边值问题进行了有限元数值模拟。在二维地电条件下,给出了边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成,最后通过求解复系数方程组得到了地表视电阻率响应。引入伪delta函数模拟线源,消弱了源带来的奇异性。通过与均匀大地以及层状介质模型的解析解对比,平均相对误差分别为0.71%和1.12%。建立了两个异常体模型,数值模拟表明异常响应比较明显,为进一步实现三维可控源电磁法有限元数值模拟提供了基础。 相似文献
17.
射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10～300 kHz,其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度,并进一步影响反演成像分辨率。针对这一问题,提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法,利用Julia的分布式计算将各频点的计算发送到不同进程进行求解,从而达到加速计算的目的,同时在计算中考虑位移电流的影响,提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应,分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明：当浅部存在高阻覆盖层时,基于准静态假设条件计算的RMT视电阻率和相位响应偏高,且频率越高、覆盖层电阻率越大,响应偏差越大;对于煤炭采空区模型,RMT法能有效反映异常体位置,但忽略位移电流会在采空区及其附近引起较大计算误差;起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应,尤其是地形拐角处;2种不同规模的并行算例对比证明了并行算法的高效性,且随着求解问题规模增大,并行算法效率也随之... 相似文献
18.
基于异常复电位2.5维CR有限元数值模拟 总被引：3，自引：0，他引：3
论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交流双点电源复电阻法的变分问题;然后利用有限单元法求解变分问题,在矩形网格中增加2条对角线形成三角网格对研究区域进行剖分,在单元中对复电位和岩、矿石复电导率进行线性插值,将变分方程化为线性代数方程组;最后,求解线性方程组,得到各节点的异常复电位值,并根据选定的测量装置计算出地表的视复电阻率,进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法的有效性后构建模型算例,给出了其不同测量装置下的计算结果。 相似文献
19.
从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序。然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性。在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图。结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置。 相似文献
20.
从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序.然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性.在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图.结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置. 相似文献