共查询到20条相似文献，搜索用时 156 毫秒
1.
从点源场的边值问题出发,利用有限元推导了井地电法和高密度电法2.5维的详细解法,并编制了相应的程序。利用均匀半空间的地电模型进行了正演计算,并与地下真实电阻率进行了对比分析,验证了算法的正确性。在此基础上对均匀介质中存在对称高阻体和边缘高阻体的模型,进行正演模拟,分别得到其井地电法的视电阻率断面图和高密度电法的视电阻率断面图。结果表明,结合两种方法可以对地下存在的高阻异常体进行有效精确地探测。 相似文献
2.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献
3.
三维地电断面激发极化法有限元数值模拟 总被引：5，自引：0，他引：5
用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的. 相似文献
4.
井地电阻率法是在井中供电地表测量的一种方法。在井地电阻率法电法实际工作中电流源通常采用点状电流源和线状电流源,其中点状电流源适用于裸眼井中供电的情况,而线状电流源则适用于套管井供电的情况。结合野外实践工作,主要研究以下几个方面并得到了相应的结论:(1)电流源为点状时,针对电流源分别在异常体的上面、里面和下面3种模型进行正演数值模拟计算,结果表明地表的视电阻率响应变化很大;(2)电流源为垂直线源时,针对电流源分别在异常体的上部、延伸至异常体中间以及贯穿异常体3种模型进行正演数值模拟计算,结果表明地表视电阻率响应的幅值变化逐渐变大;(3)电流源为倾斜线源时,线源倾斜角度不断加大对应地表视电阻率的影响较大,随着线源倾斜角的加大,地表响应的幅值变小。这些结论为野外工作中设置井地电法仪器参数以提高其观测精度提供了理论依据,为初步判断异常体的埋深做出定性的分析,同时为井地电法三维反演奠定了基础。 相似文献
5.
激发极化法具有观测参数多、观测装置灵活以及对电子导电矿物反应灵敏等特点，在矿产资源勘探中优势明显。为提高激电法对中深部矿产资源的探测能力，这里借助井地观测空间与钻孔揭露的地质信息，系统研究了井中激电阵列观测方法、正演模拟与反演成像技术。首先，设计了两种井中激电阵列观测方法；然后，推导了连续介质模型井中激电2.5维有限元正演模拟的理论公式，给出了井中激电多装置联合约束反演的线性方程；最后，对复杂地电模型的模拟数据和水槽实验数据进行反演试算，验证了本文研究的井中阵列激电法在中深部矿产资源勘探的可行性。 相似文献
6.
充水溶洞特征的高密度电阻率法反演分析研究 总被引：3，自引：3，他引：0
在高密度电法勘探中,不同性质的溶洞其反演电阻率特征有差异,溶洞为空洞时反演电阻率表现为高阻,充满水时表现则为低阻,而充水溶洞的电阻率反演特征很难判断。文章通过构建不同充水程度的溶洞模型进行正演模拟,利用最小二乘法对模型的地电断面进行反演,对比不同反演成像模型的地电响应特征。结果表明,充水溶洞的电阻率随溶洞充水量的增加而减小,低阻的范围随着水量的增加而不断扩大,其可分辨性强。野外探查资料进一步验证了模拟结果,所获得的认识对地下溶洞的电法探测分析具有一定的指导意义。 相似文献
7.
8.
直流电法正演计算中的有限单元法(第三讲) 总被引：1，自引：1，他引：1
周熙襄 《物探化探计算技术》1982,(4)
为了解决在复杂地电条件下电法勘探中场的计算,近十年来,国内外发展了电法勘探正演计算的数值模拟方法,或者称为计算机模拟,其中一个比较活跃的领域就是有限单元法。电法勘探中场的正演计算就是给定地电条件和供电电源的位置,计算地下和地面的电位分布。例如对电阻率法就是给定地下的电阻率(或导电率)的分布和地形形状,计算电位的分布,进而对于具体的供电和测量电极装置计算视电阻率值,这种正演计算是反演计算的基础。但遗憾的是直流电法正演的解析计算往往很困难,对许多地电剖面情况甚至不可能求解 相似文献
9.
10.
本刊编辑部 《物探化探计算技术》1984,(3)
本刊编辑部备有点源二维有限单元法、有限差分法和点源三维积分方程法的Fortran语言程序及BD—200语言程序各一套。这些程序适用于复杂地电断面的正演计算,也可用于反演计算。有限单元法程序还适用于计算复杂地形条件下、任意剖面(包括主剖面和旁侧剖面)的各种装置的视电阻率和视极化率值。该程序具有功能强,操作简便等特点。其计算精度和速度都超过了国外同类型程序的水平,并于1982年由地矿部物化探局评审通过。凡需要该程序以及咨询的用户,请来函向成都地质学院内《物化探计算技术》编辑部联系。 相似文献
11.
基于电偶源的体积分方程法三维电磁反演 总被引：5，自引：0，他引：5
利用体积分方程对电偶源激发时均匀导电半空间频率域三维电磁场进行正反演计算。计算中参考Eaton等提出的方法将雅可比矩阵分为线性和非线性偏微分项求解,然后利用阻尼最小二乘技术进行最优化迭代拟合。这种反演方法具有更广的适用范围和更高的计算精度,能够准确、稳定地得到较高分辨率的地下三维异常体电阻率分布图像。 相似文献
12.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献
13.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格（aggregation-based algebraic multigrid, AGMG）算法引入大地电磁三维正演模拟中。首先从准静态条件下的麦克斯韦方程出发,利用交错网格有限体积法进行离散,并采用第一类Dirichlet边界条件形成大型稀疏复线性方程组;然后阐述AGMG算法的粗化策略和套迭代技术,并实施3种不同的AGMG求解算法：1）传统的V循环AGMG算法;2）AGMG预处理共轭梯度（AGMG-CG）;3）AGMG预处理广义共轭残差法（AGMG-GCR）。最终实现大地电磁法三维正演模拟。对典型地电模型进行正演模拟,并与已有的大地电磁三维正反演程序（ModEM）进行结果对比,以验证本文算法的准确性。另外,不同剖分网格和极化方式正演模拟结果与准残量最小化（QMR）迭代算法的对比表明,AGMG预处理求解算法（AGMG-CG、AGMG-GCR）不仅能够改善算法的稳定性,而且能够快速有效地求解正演问题;其中AGMG-GCR迭代次数更少,求解速度更快,误差衰减曲线更光滑,在144×152×104网格剖分情况下,相对于现有ModEM程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
14.
在航空电磁探测方法应用中,地下探测环境复杂,地下介质物性参数存在各向异性特征。如果采用常规各向同性介质模型,在资料解释过程中会产生严重的偏差。本文基于矢量有限元法开展三维任意各向异性频率域航空电磁响应模拟计算研究。首先,将总场分解成一次场和二次场,对空气介质的均匀全空间进行一次场解析计算,同时利用矢量有限元法对二次电场的双旋度方程进行求解。为提高求解速度,采用共享内存直接求解器PARDISO对大规模稀疏矩阵并行计算,大大提高了三维模型的计算速度。之后开展了4种典型目标模型的三维各向异性介质模拟：围岩各向同性-目标体各向异性（绕z轴旋转）模型;围岩各向同性-目标体各向异性（绕x轴旋转）模型;围岩各向异性-目标体各向异性（绕z轴旋转）模型;围岩各向异性-目标体各向异性（绕x轴旋转）模型。对比分析了不同模型在不同旋转角度情况下,磁场实虚分量的变化特征,进而总结了各向异性参数对航空电磁响应的影响规律和识别方法。 相似文献
15.
基于钻孔数据的三维可视化快速建模技术的研究 总被引：1，自引：0，他引：1
为了实现三维建模过程的智能化以及自动化,以及使复杂地质体三维建模过程变得简单,针对CAD图形资料编写了能够直接将地表、钻孔等地质信息读入到数据库的DXF接口程序,提出了一种利用GA-Kriging插值手段对地层属性进行插值并最终实现三维可视化快速建模的方法,并且在此基础上基于VTK图形库实现了三维建模、等值线绘制、钻孔显示、切片以及特定地层的提取以及体积的计算等功能。在编写Kriging插值算法时,为了实现块金常数、变程等常数取值的自动化、最优化,应用GA智能算法求解。最后,通过地铁隧道的一个工程实例,对建模方法以及相应的算法模块进行了验证,表明方法具有一定的优越性。 相似文献
16.
电偶源频率电磁测深三维地电模型有限元正演 总被引：9，自引：1，他引：8
在有源变频测深(包括可控源音频大地电磁测深CSAMT)的三维有限元正演中,应用了吸收边界条件和边元有限元算法,计算精度为3%左右,基本上可满足电偶极源三维地电模型情况下频率电磁测深响应特征的分析。同时由于广义变分概念的引入,可使泛函的变分原理应用于有耗媒质电磁场问题。 相似文献
17.
18.
19.
傅氏变换有限差分(即FFD)混合波场延拓法波动方程正演模拟精度高,但三维FFD因其计算复杂、数据存储巨大而难以推广。这里提出三维FFD一步法正演模拟并行算法,把一个串行的大规模频率循环分割成并行的多个小规模频率循环,并发执行多个波场延拓来提高运算速度,并通过由给定的资源设定并行数,以及适时重排三维数据,大大增强了算法的适应性和时效性。通过对比并行和非并行正演实验,验证了算法的有效性。 相似文献
20.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献