首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 453 毫秒
1.
点源二维电法正演的有限单元法   总被引:3,自引:3,他引:3  
点源二维任意地电断面电阻率法的正演计算是电法勘探中的一个重要问题。对于这个问题,我们已经实现了有限差分方法的计算,获得了较好的效果。有限单元法是解决电法正演计算又一重要途径,在1971年Coggon首先将有限单元法引进到电法正演的模拟计算中来,以后Rijo等人又发展了这一方法,有限单元法进行电阻率法正演计算的优点,是比较容易拟合实际复杂的地电断电和起伏不平的地形。在实现点源二维地电断面视电阻率异常正演计算的有限单元法中,我们将混合边界条件引入求解,并采朋了若干优化措施,使有限单元法计算视电阻率异常的精度和速度均有进一步提高,获得较好的计算效果。在文献中我们简单地介绍了这个方法,下面我们较详细  相似文献   

2.
变步长激发极化法有限元数值模拟   总被引:1,自引:0,他引:1  
用改变剖分步长的方法对三维地电模型进行激发极化法正演模拟计算.首先,在给出三维地电模型异常电位的边值问题与变分问题的基础上,通过对三维激电模型进行变步长的六面体剖分,将异常体部分用较短步长的单元进行加密剖分,并在各个单元上进行插值,于是得到其总刚度矩阵;然后用变带宽方式对其进行编号存储,并选用成熟的LDLT方法求解大型稀疏矩阵,运用等效视电阻率方法求取地表视极化率值.通过用中间梯度装置对几例典型的激电模型进行了正演计算,结果表明,该方法能够达到满意的精度和效果,从而为激发极化正演解释提供有益的帮助.  相似文献   

3.
将一种快速数值模拟方法用于激电法正演模拟中,利用在低电阻率差情况下,积分方程法模拟时阻抗矩阵的非对角线项可忽略,而只需计算矩阵主对角值这一关键点。这里详细阐明了三维地电断面激电法快速模拟方法,推导了求解过程,并以此为基础编制了计算程序。实例试算结果说明,该模拟方法在计算速度、计算精度上都收到了较满意的结果。  相似文献   

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

5.
积分方程法三维模拟井地电法并行算法研究   总被引:1,自引:1,他引:1  
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。  相似文献   

6.
直流电法正演计算中的有限单元法(第三讲)   总被引:1,自引:1,他引:1  
为了解决在复杂地电条件下电法勘探中场的计算,近十年来,国内外发展了电法勘探正演计算的数值模拟方法,或者称为计算机模拟,其中一个比较活跃的领域就是有限单元法。电法勘探中场的正演计算就是给定地电条件和供电电源的位置,计算地下和地面的电位分布。例如对电阻率法就是给定地下的电阻率(或导电率)的分布和地形形状,计算电位的分布,进而对于具体的供电和测量电极装置计算视电阻率值,这种正演计算是反演计算的基础。但遗憾的是直流电法正演的解析计算往往很困难,对许多地电剖面情况甚至不可能求解  相似文献   

7.
从点源场的边值问题出发,利用有限元推导了井地电法和高密度电法2.5维的详细解法,并编制了相应的程序。利用均匀半空间的地电模型进行了正演计算,并与地下真实电阻率进行了对比分析,验证了算法的正确性。在此基础上对均匀介质中存在对称高阻体和边缘高阻体的模型,进行正演模拟,分别得到其井地电法的视电阻率断面图和高密度电法的视电阻率断面图。结果表明,结合两种方法可以对地下存在的高阻异常体进行有效精确地探测。  相似文献   

8.
二维高密度电阻率法数值模拟与应用   总被引:2,自引:0,他引:2  
金聪 《地质与勘探》2014,50(5):984-990
高密度电阻率法作为工程勘察领域广泛使用的勘察方法之一,与传统电阻率法相比,具有效率高、信息丰富等优点。本文首先从高密度电法原理出发,介绍了两种常见的装置(温施装置和温纳装置)的特点。然后依据实测资料建立了地电理论模型,并通过理论正反演数值模拟,分析了两种装置下模型的电性响应特征,并与实际资料与进行了比较,进而进行地层分层和地质解释。结果表明,温施装置进行地质分层效果更佳,通过高密度电法正演模拟指导反演可提高地层分层和地质解释的精度。  相似文献   

9.
为探究坑道聚焦直流激电法电流场分布特性及其探测影响因素,利用有限单元法对空间场正常电位和异常电位进行求解,基于Comsol软件构建均质三维地电探测模型,采用四面体网格自适应算法进行剖分,对比分析数值求解精度。研究聚焦电流场分布变化规律,确定聚焦效应电流比系数取值范围。沿掘进断面的正前方构造异常地质的三维地电模型,利用等效电阻率法正演模拟掘进断面面积、干扰体、聚焦效应电流比系数等因素对激电效应参量的影响。结果表明:聚焦直流激电法对掘进断面前方不良地质体的探测具有较好的敏感性,坑道腔体内干扰体对目标异常体探测的影响可以忽略,旁侧干扰异常体距坑道底板距离越远,对目标异常体探测影响越小;增大掘进断面面积和聚焦效应电流比系数,可有效增加勘探距离。此研究可为聚焦直流激电法反演提供依据和参考,对推动聚焦电法勘探理论的发展具有重要意义。  相似文献   

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

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

12.
电阻率/激电测深二维人机交互正演模拟   总被引:1,自引:0,他引:1  
在有限元点源二维电阻率正演基础上,针对传统对称四极装置其电极距跨度大且不等间隔分布等特点,提出将整个正演计算区域划分为目标区、观测区和扩展区三个区域,在观测区域先采用等间隔剖分,然后在区域中插入实际电极作为正演计算节点的剖分方式,同时将无穷远边界条件简化技术引入到正演计算中,并结合人机交互与可视化技术实现了大极距电阻率/激电测深的二维正演模拟。最后通过对相关模型的试算结果表明,作者在本文中提出的针对传统对称四极装置电阻率/激电测深的二维正演模拟技术是有效可行的。  相似文献   

13.
在三维电阻率的正反演计算中,快速、准确的正演计算是反演的关键。而正演计算往往涉及到求解大型线性方程组Ax=b的问题,通过Lanczos迭代构造出对称三对角阵方程组,并采用正交分解法进行求解,与传统算法相比,此算法占用内存少、收敛速度快、且稳定;针对大型稀疏矩阵的特点,采用简单地记录矩阵的非零元素值及其所在行、列值的方法,来存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。通过理论分析和点电源三维地电场计算实例,阐述该法是地电三维正演计算的有效方法。  相似文献   

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

15.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。  相似文献   

16.
高密度电法是一种以地球物理为基础,集电剖面和电测深为一体,采用高密度布点,进行二维地电断面测量,研究地下介质体电阻率差异的一种勘查技术.在此,以厦门至成都公路四川境川黔界至纳溪段高速公路C9合同段大田头大桥为例,介绍了高密度电法的基本原理及其在勘查中的应用效果.  相似文献   

17.
时间谱电阻率法的二维正演算法   总被引:5,自引:0,他引:5       下载免费PDF全文
孟永良  罗延钟 《地球科学》2000,25(6):656-662
时间谱电阻率(TSR)法的二维正演, 是计算可极化二维地电构造上三维电流偶极源的电场瞬变响应, 属所谓2.5维时间域电磁场数值模拟问题, 是目前国际上未妥善解决的计算地球物理疑难问题.针对现有算法的局限性, 建立了新的算法, 其特点是: (1)采用二维有限单元算法, 在矩形网格中增加两对角线形成三角网格剖分, 同时用高斯消元法消除矩形网格中心结点的待求未知量.这样, 既可较准确地模拟任意二维复杂地电断面, 又可节省计算量.(2)采用直接计算二次场的新算法, 只需计算电场和磁场沿地电构造走向两个分量的一次场, 因而, 不但计算精度较高, 而且不显著增加计算量.(3)采用G -S变换法作逆拉氏变换, 并利用拉氏变换延迟定理在倍增的时间间隔中插值, 从而实现对密集采样时间瞬变过程的快速计算.(4)能对可极化和导电大地(即同时包括IP和EM效应).   相似文献   

18.
随机岩溶介质的电阻率正演模拟   总被引:1,自引:1,他引:0  
受岩溶发育的影响,岩溶地区电阻率分布具有很大的不确定性。采用传统的正演方法在一定程度上不能满足岩溶探测的需要。为了更真实地反映岩溶介质的电性特点,本文提出了一种根据随机过程的谱分解理论来构造随机介质模型的方法,并以此建立了岩溶地区随机地电模型。对理论和实测地电模型的正演计算表明,随机介质能更有效模拟实际地电结构,正演结果更接近实测的异常形态。因此,岩溶地区的电阻率正演模拟计算采用随机介质模型是十分有必要的,一方面有助于对异常的认识,另一方面,随机介质模型可以更好地模拟岩溶地区的地电模型。   相似文献   

19.
可控源音频大地电磁法(CSAMT)正演具有计算耗时、精度低的特点;常规的反演方法边界处理简单,占用内存大,搜索方向也难以精准,严重影响反演效果。为了实现CSAMT快速高精度正演,改善边界处理方法以减少内存占用,优化搜索方向以减少迭代反演次数,提出基于边界约束有限内存的拟牛顿CSAMT一维反演方法。首先,在CSAMT正演中基于高斯积分与并行算法实现了一维介质模型条件下有限长导线源激发的水平电场快速高精度计算。其次,在反演方法中,分别构建了基于相对误差和绝对误差的两种目标函数,并引入光滑模型约束,采取自适应正则化策略更新正则化因子,并采用基于边界约束的有限内存拟牛顿算法(LBFGS_B算法),实现了CSAMT一维准确快速反演。以含低、高阻薄层的多层地电模型的水平电场Ex振幅为反演数据进行直接反演和模型检验。结果表明,该方法能够反映出地层电阻率随深度变化的趋势,并对高、低阻薄层有较好的分辨能力,低阻薄层分辨效果更佳。实际CSAMT资料的反演剖面的分辨率优于常规的连续介质反演方法,地质效果明显,表明该方法具有良好的应用前景。   相似文献   

20.
在徐世浙等(1995)工作的基础上,给出了二维地电断面电阻率测深的变分问题以及电导率分块线性变化时的有限元方法数值解,并计算了电导率随深度变化时的地电断面的测深曲线.结果与解析解相比,精度很高.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号