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

2.
在使用有限单元法解决了井-井2.5维正演问题后,采用矩形剖分方式将井间剖面离散成若干个点,再将正演所得的视电阻率数据采用一种新颖的方式分配重组到井间每个离散点上,绘制出了井间剖面电阻率近似成像图。该近似成像方法避开了反演成像中庞大的内存需求和计算量,计算效率高,并且成像结果表明该方法具有较好的效果,是一种行之有效的井间视电阻率资料处理手段。  相似文献   

3.
复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。  相似文献   

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

5.
不平地形条件下高密度电阻率法的2.5维反演   总被引:6,自引:0,他引:6  
讨论了不平地形条件下电阻率法2.5维正演和反演的算法,并在此基础上编制了高密度电阻率法2.5维反演程序,该程序可用于九种常用电极装置观测结果的反演.对理论和实测数据的反演结果表明,本文的算法正确,程序运行稳定,反演效果很好.  相似文献   

6.
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。  相似文献   

7.
首先以频点为并行粒度,完成了MT三维正演并行计算,然后将该算法引入共轭梯度法MT三维反演过程中的正演求数据的残差、反演方程右端项的求取、共轭梯度求解模型增量这三个计算环节中,使得MT三维反演中的主要计算量实现并行化。对编制的MT三维正、反演并行计算程序,在多核工作站上通过理论地电模型进行了正演试算和反演验证,并在串行和并行两种计算方式下进行了计算效率对比。  相似文献   

8.
基于异常复电位2.5维CR有限元数值模拟   总被引:3,自引:0,他引:3  
论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交流双点电源复电阻法的变分问题;然后利用有限单元法求解变分问题,在矩形网格中增加2条对角线形成三角网格对研究区域进行剖分,在单元中对复电位和岩、矿石复电导率进行线性插值,将变分方程化为线性代数方程组;最后,求解线性方程组,得到各节点的异常复电位值,并根据选定的测量装置计算出地表的视复电阻率,进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法的有效性后构建模型算例,给出了其不同测量装置下的计算结果。  相似文献   

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

10.
复电阻率法在油气资源、矿产勘查中发挥了重要的作用,为了深刻认识复电阻率法异常特征变化规律,本文对复电阻率2.5D正演问题展开研究。首先直接给出复电阻率法2.5D有限元正演所满足的变分问题,并详细地推导相应的刚度矩阵的计算过程。引入Fractal模型作为等效模型研究频谱激电异常特征。对单元内的复电导率及复电位均进行线性插值,而后,采用不完全LU分解的稳定双共轭梯度算法求解有限元线性方程组,获得异常复电位值。设计3个典型的地电模型验证了本文算法的正确性及精确性,并分析了不同装置下,不同频率的2.5D复电阻率异常响应特征。数值模拟结果表明,采用Fractal模型研究激发极化异常特征是可行、有效的;不同装置、不同频率下的复电阻率法异常特征有着显著的差异。  相似文献   

11.
大地的感应激发极化效应有时会在瞬变电磁晚期响应上产生负值变号现象。常规实电阻率瞬变电磁反演由于没有考虑激发极化效应,对于观测数据负值部分的解释一直存在问题,这影响了反演解释的真实性和可信度。对此,本文首先将Cole-Cole复电阻率模型加入正演计算,并验证了计算的准确性。然后在阻尼最小二乘法的基础上加入Occam平滑约束来构建反演方程,能够同时反演出零频电阻率、充电率、时间常数以及频率相关系数,实现了一维瞬变电磁复电阻率反演算法。最后,建立具有不同程度激发极化效应的典型层状理论模型进行反演试算与结果分析,结果表明,在激发极化效果明显的低阻高极化地层中,复电阻率反演效果更好。与实电阻率反演结果的对比说明,瞬变电磁复电阻率反演既可以达到实电阻率的常规反演效果,也能解决实电阻率无法实现的负值拟合问题。  相似文献   

12.
电阻率图像重建的佐迪(Zohdy)反演及其应用效果   总被引:19,自引:0,他引:19  
王兴泰  李晓芹 《物探与化探》1996,20(3):228-233,208
本文讨论了利用佐迪反演方法对地表高密度电阻率法的数据采集结果进行电阻率图像重建。通过对多种复杂地电模型的数值模拟和野外实例的反演分析,表明该方法简单易行,计算成本较低,成像的分辨率高,效果较理想,将其用于地表测量的电阻率层析成像完全可行。作者用Fortran语言编制了实现上述过程的程序,对所有数学模型的正演计算均采用2.5维有限元法。为突出成像效果,本文在结果图示时均采用了伪彩色图像处理技术。  相似文献   

13.
大回线瞬变电磁法全区视电阻率的逆样条插值计算   总被引:7,自引:2,他引:5  
给出了一种方便快捷的瞬变电磁全区视电阻率数值计算方法--逆样条插值法。对几例由2.5维正演程序获得的理论数据的计算结果表明,逆样条插值法用于全区视电阻率的定义是可行的。与此同时,该方法的实现也为瞬变电磁自动反演中初始模型的构制奠定了基础。  相似文献   

14.
为探究三维地质体对二维反演结果地影响以及二维反演中极化模式的选择问题,设计了一系列二维、三维地电模型进行正演计算,并对正演结果进行二维反演计算。研究结果表明:相位受到三维畸变的影响较小,视电阻率较大;当三维异常体模型走向延伸较小时,TE模式的反演结果误差较大,TM模式的反演结果相对较好,但是反演得到了虚假的下覆构造,TE+TM联合反演结合了TE和TM模式反演的优点,可以较好地反映异常体的赋存位置及电阻率值;当三维异常体模型走向延伸长度增大时,三维数据的二维反演结果是可靠的,用TM和TE+TM模式反演都比较合理。  相似文献   

15.
开展了不同观测方式的井地2.5D直流电阻率反演研究。①从2.5维直流电阻率满足的边值问题出发,采用变分原理结合节点线性基函数推导了2.5维井地DC满足的积分弱解形式;②构建了二阶最大平滑稳定泛函的2.5维井地直流电阻率正则化目标函数,采用共轭梯度算法对正则化目标函数进行最优化求解,并采用逐步衰减正则化因子的求解策略来提高反演的稳定性;③设计了均匀半空间模型得到的数值解与解析解对比,电位的相对误差在2%以内,阐述正演算法的正确性和高精度。另外,分析了不同观测装置的2.5维井地直流电阻率异常体特征,并对不同观测方式对2.5维井地DC理论数据进行反演研究。研究结果表明,井中数据的引入提高了2.5DC对纵向探测的分辨率能力,同时提高了2.5DC反演有效性以及准确性。  相似文献   

16.
Zohdy方法在三维电阻率测深数据反演中的应用   总被引:2,自引:0,他引:2  
三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果.  相似文献   

17.
在寻找金属和非金属矿产,解决水文,工程地质问题和地热资源的勘探中,电阻率法是一种重要的物探方法,但由于该方法计算复杂,正演、反演问题往往停留在定性阶段,从而大大限制了该方法的应用效果和解决问题的能力,我们有必要逐步实现电阻率法资料解释的定量化。如果我们已知一个地区的地电资料,对给定的装置和供电电流值计算电异常值或视电阻率值,这个问题即是电法勘探的正演问题,它是电阻率法定量解释的基础。由于地电条件的复  相似文献   

18.
一种新的复电阻率频谱参数的递推反演方法   总被引:11,自引:1,他引:11  
从Cole-Cole模型出发,导出了一种新的复电阻率谱参数的递推反演方法。该方法不需要给定参数初始值,反演过程中也无需任何中间量的介入,简单方便。尤其是只需要四个不同频率的复电阻率,从而提高了复电阻率法野外的工作效率。反演结果表明,这种方法具有较高的反演精度。  相似文献   

19.
二维直流电阻率测深曲线的快速反演   总被引:5,自引:1,他引:5  
本文采用作者提出的一维直流电阻率测深的直接反演公式,结合二维正演和渐近迭代反演,推出了一种二维直流电阻率测深的快速反演方法,经实际应用,效果很好。  相似文献   

20.
在山区进行音频大地电磁测深(AMT)时,电磁场易受地形起伏影响发生畸变。为研究起伏地形对AMT数据的影响,采用WinGlink软件基于有限差分算法的正演模块,建立二维起伏地形和平地地电模型,并对比分析了两个理论模型的正演响应。结果表明:TE模式的视电阻率曲线在高频段受干扰程度相对较大,在低频段影响减小甚至消失;TM模式的视电阻曲线变化正好相反;二者相位受地形影响则较小。对含异常体的带地形地电模型进行正演计算,在正演响应中加入2.5%的随机误差模拟实测数据,利用WinGlink软件基于非线性共轭梯度(NLCG)算法的反演模块,分别进行带地形和不带地形的TE、TM和TE+TM模式的二维反演。各反演结果的对比表明,带地形的二维反演结果明显优于不带地形的反演结果,能够较好地反映出异常体的位置、形态及电阻率值。依据得到的认识,对青海祁连隧道工程的AMT实测数据进行带地形二维反演,取得了理想的应用效果。  相似文献   

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

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