共查询到16条相似文献，搜索用时 171 毫秒
1.
《物探化探计算技术》2015,(3)
针对目前地面直流电阻率法所面临的勘探深度有限的问题,开展了基于TV井地2.5D直流电阻率正则化反演研究。首先从井地直流电满足的边值问题出发,采用三角单元二次插值有限元法推导了2.5D井地电阻率正演公式。为了提高计算效率,对边界采用近似处理、以及利用图论理论的矩阵重排与填入元分析方法,实现了稀疏矩阵的直接解法,大大提高了计算效率。为了能够提高反演结果的稳定以及对异常边界分辨能力,采用L曲线法求取正则化因子、TV的稳定因子以及共轭梯度方法进行反演目标函数的求解,研究表明,反演结果高效、准确。最后对不同井地观测方式的反演结果进行对比分析,得出井地直流电阻率异常特征和分布规律。 相似文献
2.
直流电阻率测深二维反演中,正则化参数的选取影响反演结果分辨率及反演过程稳定性。利用主动约束平衡正则化因子,进行直流电阻率光滑约束最小二乘二维自适应反演,改善直流电阻率测深二维反演的分辨率与稳定性。在反演迭代过程中,正则化因子根据模型参数的空间展布函数进行自适应计算、正则化参数的自适应计算。模拟数据反演结果验证了该方法的有效性与可行性,反演结果能准确地反映地下模型的真实电性结构。 相似文献
3.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明：基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。 相似文献
4.
二维大地电磁正则化共轭梯度法反演算法 总被引：2，自引：0，他引：2
针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。 相似文献
5.
6.
直流电阻率法因效率高、成本低等特点广泛应用于有色金属、煤田等矿产资源勘探和地质调查行业。但在井–地、地–井电阻率法勘探过程中,钻孔因素(井液电阻率和钻孔孔径)对不同探测模式下视电阻率响应的影响尚未明确,是否影响视电阻率资料的解释是值得探讨的一个问题。因此,提出采用混合网格有限元法实现直流电阻率三维正演。给出异常电位法满足的边值问题及有限元变分问题,应用三棱柱和四面体混合网格实现对计算区域的快速离散,并建立2种网格的线性插值基函数和单元系数矩阵;采用SSOR–PCG迭代算法求解异常电位满足的大型线性方程组,得到各观测点的响应。在保证计算精度的前提下,应用混合网格有效地离散钻孔地电模型,探讨钻孔因素对井–地、地–井观测方式视电阻率数据的影响特征。对于井–地观测方式:钻孔附近视电阻率受钻孔因素影响最大,严重影响视电阻率资料的合理解释,但随着测点收发距的增加,视电阻率响应逐渐趋于围岩电阻率;而且随着发射源深度的增加,钻孔对地表视电阻率响应的影响也逐渐减小。对于地–井观测方式:钻孔因素对浅部数据影响较大,对深部数据影响小,相比于井液电阻率,视电阻率响应更易受钻孔孔径的影响。直流电阻率法混合网格有限元三维正演算法的提出,对实际井–地、地–井电阻率法勘探具有理论指导意义,可结合钻孔信息进行正演模拟分析,选择合适的收发距,有效压制钻孔因素对实测资料的影响。 相似文献
7.
地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段之一,正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容。稳定因子的主要功能是对模型的解空间进行限制,以减少多解性,求得稳定解;为更好的定位地下目标体,研究工作采用了最小支持稳定因子,并给出了统一表达形式。正则化因子在模型误差与数据误差之间起到平衡作用,研究工作在＂L曲线法＂自动选取正则化因子的基础上,提出一种修正的计算正则化因子算法。以2.5维直流电阻率为例,运用非线性共轭梯度法,进行了试算研究,结果表明了修正的＂L曲线法＂算法快速、稳定,最小支持稳定因子具有较好的聚焦特性,有利于实现陡变边界反演。 相似文献
8.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
9.
直流电阻率法是一种较常用的地球物理勘探方法,常用于浅层勘探,因此受地形影响较大.其具有较高的勘探分辨率,但不同的布极方式会对灵敏度产生不同的影响.因此,研究不同布极方式下的带地形直流电阻率法反演规律和特点具有一定的意义.这里对带地形的点源直流电阻率法采用有限体积法进行二维正演模拟,然后将该模拟结果引入正则化Tikhon... 相似文献
10.
11.
为了解地下介质物理特性随时间的变化,采用地球物理监测数据的成像方法已经逐渐得到应用。一种常用的途径就是通过独立反演不同时刻采集的地球物理数据,来重构间歇影像,这使得实测的数据信息得不到充分利用。因此,这里提出了一种同时反演所有不同时刻观测数据的新方法,即通过一次反演过程来获得所有不同时刻的地下模型。在反演过程中,通过引入时间和空间正则化因子来改善反演的稳定性,减少反演的多解性。最后,以直流电阻率监测数据反演为例,验证所提算法的有效性。 相似文献
12.
电阻率法可用于地表浅层的探测,也可用于对动态地下目标进行监测。对于监测数据的反演,不同数据集的单独反演存在一定缺陷,为此本文在常规电阻率反演算法的基础上,推导了时移电阻率反演公式,实现了时移反演算法程序;为了论证时移反演算法对动态地下目标成像的优越性,建立一组多个正演模型,利用模拟数据进行单独反演和时移反演,并对比两种方法的结果。研究表明,尽管两种算法都能圈定动态地下目标体,但时移反演算法可以消除不同观测数据集中所包含的随机误差,减少假异常的出现。 相似文献
13.
为消除TEM定源回线装置增大接收区域后明显的边框效应,且提高TEM的分层能力,达到提高对不同电性地层精细分辨的目的,采用带约束的正则化反演方法。在反演算法中,拉格朗日因子采用CMD正则化因子的自适应调节方案,自主调整。为保证反演过程的高效、稳定且结果收敛性,使反演模型每一层的厚度按照等对数间隔进行离散,减少矩阵方程的病态性;同时对反演电阻率值采用可行的上下限约束,使迭代的结果不至于出现极值或错值,降低反演的多解性。通过对多层理论模型的计算,证明了约束反演算法的正确性和可靠性,将其用于实测瞬变电磁法勘探数据的反演计算中,反演电阻率纵向分布规律与已知地层电性资料吻合,解释的异常区得到了钻探验证,结果可靠。 相似文献
14.
大定源瞬变电磁一维自适应正则化反演 总被引：1，自引：0，他引：1
自适应正则化反演对初始模型要求较低,直接对瞬变电磁响应数据进行反演。正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应的得到,从而快速获得地下的地电结构。本文采用自适应正则化反演算法对大定源回线瞬变电磁一维层状模型进行反演,使用均匀半空间作为初始模型,首次采用模型对深度的二阶导数极小的模型约束,通过典型理论模型的反演计算,证明了TEM自适应正则化一维反演算法拟合效果好、精度高,由于不用反复搜索正则化因子,并且收敛速度快,体现了良好的稳定性和可靠性。 相似文献
15.
Seawater intrusion (SWI) is a complex process, where 3D modeling is often necessary in order to monitor and manage the affected aquifers. Here, we present a synthetic study to test a joint hydrogeophysical inversion approach aimed at solving the inverse problem of estimating initial and current saltwater distribution. First, we use a 3D groundwater model for variable density flow based on discretized flow and solute mass balance equations. In addition to the groundwater model, a 3D geophysical model was developed for direct current resistivity imaging and inversion. The objective function of the coupled problem consists of data misfit and regularization terms as well as a coupling term that relates groundwater and geophysical states. We present a novel approach to solve the inverse problem using an alternating direction method of multipliers (ADMM) to minimize this coupled objective function. ADMM enables to treat the groundwater and geophysical part separately and thus use the existing software with minor changes. To further reduce the computational cost, the sensitivities are derived analytically for the discretized system of equations, which allows us to efficiently compute the gradients in the minimization procedure. The method was tested on different synthetic scenarios with groundwater and geophysical data represented by solute mass fraction data and direct current resistivity data. With the ADMM approach, we were able to obtain better estimates for the solute distribution compared to just considering each data separately, solving the problem with a simple coupled approach or by a direct substitution of the coupling constraint. 相似文献
16.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。 相似文献