首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明:基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。  相似文献   

2.
长导线源半航空瞬变电磁正则化反演正则项通常采用L2范数,其拟合结果较光滑,不能有效刻画层界面信息。针对层状介质陡变模型实现正则项为L1范数的反演算法,采用迭代重加权最小二乘法将原问题转化为L2正则化子问题求解,解决L1范数存在不可导问题;采用OpenMP技术对雅可比矩阵并行计算,提高了反演速度;对自适应正则化因子分段迭代法的调整策略进行分析并改进,改进后的自适应正则化因子调整策略更适合半航空瞬变电磁L1正则反演算法。最后对电阻率进行反演并与Occam反演结果作比较,结果表明L1正则反演充分迭代后能够突出符合真实模型的电性界面,反演电阻率与模型真实值更接近。  相似文献   

3.
为了适应地形起伏的实际地质情况,开展带地形的最小二乘二维反演研究。鉴于大地电磁(MT)反演的不适定问题,引入Tikhonov的正则化方法,从而获得关于总目标函数的方程,利用光滑约束最小二乘法求解总目标函数方程。由于正则化因子值与反演精度以及稳定性相关,采用主动约束平衡方法获取最优化的正则化因子,以确保反演精度和稳定性都达到最佳。与此同时,利用电磁场互易定理以节省反演迭代过程求解雅可比矩阵的计算时间。构建了若干地质构造模型进行试算,分别讨论TE、TM模式以及二者联合模式的反演结果,并与前人研究工作对比以说明本文方法的反演效果。  相似文献   

4.
直流电阻率法是一种较常用的地球物理勘探方法,常用于浅层勘探,因此受地形影响较大.其具有较高的勘探分辨率,但不同的布极方式会对灵敏度产生不同的影响.因此,研究不同布极方式下的带地形直流电阻率法反演规律和特点具有一定的意义.这里对带地形的点源直流电阻率法采用有限体积法进行二维正演模拟,然后将该模拟结果引入正则化Tikhon...  相似文献   

5.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。  相似文献   

6.
COMSOL是一个多物理场仿真软件,笔者尝试将其应用于大地电磁法的数值模拟,研究起伏地形对大地电磁响应的影响及解决方案。本文由麦克斯韦方程组出发,施加第一类边界条件,将地形起伏数字化为参数化曲面,导入COMSOL进行计算。在对COMSOL的精度进行了验证的基础上,分析不同频点所需内存以及自由度,然后选择合适的频点进行研究。将三维起伏地形的模拟结果与无地形起伏时的结果进行对比,总结出山谷、山峰地形对三维大地电磁响应的影响,发现在XY模式下,地形对视电阻率的影响大于相位,在视电阻率等值线图中能够清晰地看出起伏地形的轮廓,而在YX模式下,地形对相位的影响大于对视电阻率的影响。最后,利用“比值法”进行地形校正,校正前由于起伏地形影响,在异常体附近会出现部分视电阻率数值跳跃现象,通过校正有效抑制了起伏地形对异常体附近视电阻率的影响,校正后的视电阻率曲线整体上移,并且异常体周围的视电阻率值也趋于背景场值,获得了较好的效果。  相似文献   

7.
复杂地形下高密度激电法2.5维有限单元法数值模拟   总被引:1,自引:0,他引:1  
目前三维有限元模拟计算量大,计算效率低,对计算机要求较高,2.5维模拟是三维问题的简化,较好地克服了上述问题。在高密度激电法模拟中,模拟的点数一般较大,要求程序计算效率高,因此对高密度激电法2.5维模拟研究是必要的。首先给出了2.5维稳定电流场的边值问题及对应的变分问题;为了能更好模拟复杂地形对极化率异常的影响,采用三角形网格,有利于模拟复杂地形又有利于将其应用于反演计算。推导了基于连续电性介质的2.5维稳定电流场的有限元法,利用等效电阻率法,编制了2.5维高密度激电法有限元模拟程序,计算了水平地形及起伏地形下的极化率异常模型,并验证了方法的可行性,分析得出了山谷地形比山脊地形对极化率观测的影响大的结果。  相似文献   

8.
针对目前地面直流电阻率法所面临的勘探深度有限的问题,开展了基于TV井地2.5D直流电阻率正则化反演研究。首先从井地直流电满足的边值问题出发,采用三角单元二次插值有限元法推导了2.5D井地电阻率正演公式。为了提高计算效率,对边界采用近似处理、以及利用图论理论的矩阵重排与填入元分析方法,实现了稀疏矩阵的直接解法,大大提高了计算效率。为了能够提高反演结果的稳定以及对异常边界分辨能力,采用L曲线法求取正则化因子、TV的稳定因子以及共轭梯度方法进行反演目标函数的求解,研究表明,反演结果高效、准确。最后对不同井地观测方式的反演结果进行对比分析,得出井地直流电阻率异常特征和分布规律。  相似文献   

9.
为解决大角度起伏地形即陡峭地形(如河流阶地、冲沟、堤坝、山谷、山脊等)对实测视电阻率的影响问题,这里优化了三维有限元数值模拟中的单元积分形式,使之适合大倾角地形条件下对电阻率的数值模拟.为验证方法的正确性,采用保角变换计算了相近地形在均匀电流场情况下,陡峭地形引起的视电阻率变化特征.计算表明,对于缓倾斜地形,有限元数值...  相似文献   

10.
直流电阻率测深二维反演中,正则化参数的选取影响反演结果分辨率及反演过程稳定性。利用主动约束平衡正则化因子,进行直流电阻率光滑约束最小二乘二维自适应反演,改善直流电阻率测深二维反演的分辨率与稳定性。在反演迭代过程中,正则化因子根据模型参数的空间展布函数进行自适应计算、正则化参数的自适应计算。模拟数据反演结果验证了该方法的有效性与可行性,反演结果能准确地反映地下模型的真实电性结构。  相似文献   

11.
重磁异常的物性线性反演多是欠定问题求解,求解不仅具有明显的不稳定性,还会导致其垂向分辨率大大降低,反演的"趋肤效应"明显.针对重力数据反演中的这一问题,引入重加权正则化共轭梯度法的约束反演,即在目标函数中加入正则化方程加以约束,以深度加权函数对核矩阵进行补偿,在反演过程中实现了正则化因子的自适应选择.引入了PRP公式,...  相似文献   

12.
二维大地电磁正则化共轭梯度法反演算法   总被引:2,自引:0,他引:2  
针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。  相似文献   

13.
曲面位场的正则化线性规划法直接反演技术   总被引:4,自引:0,他引:4  
讨论了可用于起伏地形上重磁位场数据直接反演的、正则化的线性规划反演方法,给出了求解模型、求解方法及相关的技术措施。通过理论模型试算,证明了方法的可行性,在某铜镍矿区的应用验证了方法的实用性。  相似文献   

14.
高密度电阻率法在岩溶区找水具有较好的效果,但岩溶区地形复杂,高密度电法受地形起伏影响较大。为了提高岩溶区高密度电法找水的勘探效果,掌握地形起伏对高密度电法影响的规律,本文采用高密度电阻率法微测系统,利用铜柱体模拟地下岩溶管道,研究地形起伏对高密度电法探测岩溶管道的影响。研究结果表明:高密度联合剖面法及高密度α2装置受谷地地形影响较大,容易在谷地中间产生虚假的低阻异常,矩形AMN和滚动MNB装置受地形影响较小;受谷地地形影响,高密度电法能勘探的有效异常体深度减小;同一电极距下,平地模型及谷地模型反演异常体顶部埋深均小于真实异常,横向宽度均大于真实异常,平地模型反演异常体的顶部埋深大于谷地模型反演异常体,横向宽度小于谷地模型反演异常体。  相似文献   

15.
航空大地电磁法以其机动灵活、效率高等优点,适用于在地势崎岖的偏远环境中开展快速普查作业.在航空大地电磁探测中,起伏地形会对倾子响应产生影响,另外由于倾子是无量纲传输函数,只包含电性结构的空间变化信息,而与电阻率的数值水平无关,单独利用倾子数据进行反演时,反演多解性更大.为解决上述问题,首先采用四边形网格模拟地形起伏;其...  相似文献   

16.
为了解起伏地表对复电阻率法数据的影响,发展了起伏地表条件下复电阻率法2.5维有限元数值模拟方法。经数值模拟结果与解析解对比表明,该模拟方法在计算精度上得到了较满意的结果;几种常见纯地形异常的数值结果表明:起伏地形对复电阻率法数据的畸变特征与地形坡度大小以及地形分布范围有很强的相关性,坡度越大,地形影响越为剧烈;坡度角为45°的山脊地形、山谷地形和斜坡地形所引起的复电阻率法实分量异常分别达到背景值的21.2倍、3.9倍和3.4倍,虚分量异常分别达到背景值的2.3倍、2.2倍和2.1倍,振幅异常分别达到背景值的21.2倍、3.9倍和3.4倍;而相位受地形影响非常微弱。这些工作为复电阻率法的地形改正和反演奠定了基础。  相似文献   

17.
正演是反演技术的基础,正演速度和求解反演问题的系数矩阵存放一直是起伏地形下重、磁三维反演的关键技术问题。这里提出了一种起伏地形下重磁快速正演计算方法,其计算原理是根据反演在垂向的剖分层数,利用水平地形正演计算形成二个不同大小刻度标尺矩阵,然后在模型空间,使用分段线性插值的方式,直接计算出起伏地形观测点的正演值。该方法的主要特点是在保持很高的计算精度下计算速度可提高二倍,且节省计算内存,适合起伏地形下重磁三维反演技术研究。  相似文献   

18.
针对太原市蒙山煤矿采空区地形起伏较大、需弯线测量等因素制约的特殊场地条件,提出了阵列形式观测的高密度电阻率法和单点测量方式的瞬变电磁法相结合的综合物探方法.前者二维异常反演分辨率较高,可弥补瞬变电磁法浅部探测盲区,后者具有探测深度大的优势,采用一维异常反演方式可弥补高密度电阻率法受剖面折线展布和地形起伏影响的缺点.由此查明了该矿区采空区的位置及分布情况,后期钻探验证表明所采用的综合物探方法是行之有效的.该综合物探方法组合模式实现了两种物探方法的优势互补,可为同类场地条件下探测煤层采空区提供有意义的借鉴.  相似文献   

19.
地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段之一,正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容。稳定因子的主要功能是对模型的解空间进行限制,以减少多解性,求得稳定解;为更好的定位地下目标体,研究工作采用了最小支持稳定因子,并给出了统一表达形式。正则化因子在模型误差与数据误差之间起到平衡作用,研究工作在"L曲线法"自动选取正则化因子的基础上,提出一种修正的计算正则化因子算法。以2.5维直流电阻率为例,运用非线性共轭梯度法,进行了试算研究,结果表明了修正的"L曲线法"算法快速、稳定,最小支持稳定因子具有较好的聚焦特性,有利于实现陡变边界反演。  相似文献   

20.
为消除TEM定源回线装置增大接收区域后明显的边框效应,且提高TEM的分层能力,达到提高对不同电性地层精细分辨的目的,采用带约束的正则化反演方法。在反演算法中,拉格朗日因子采用CMD正则化因子的自适应调节方案,自主调整。为保证反演过程的高效、稳定且结果收敛性,使反演模型每一层的厚度按照等对数间隔进行离散,减少矩阵方程的病态性;同时对反演电阻率值采用可行的上下限约束,使迭代的结果不至于出现极值或错值,降低反演的多解性。通过对多层理论模型的计算,证明了约束反演算法的正确性和可靠性,将其用于实测瞬变电磁法勘探数据的反演计算中,反演电阻率纵向分布规律与已知地层电性资料吻合,解释的异常区得到了钻探验证,结果可靠。   相似文献   

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

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