共查询到18条相似文献，搜索用时 125 毫秒
1.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
2.
地面核磁共振正则化反演方法研究 总被引：2，自引：1，他引：1
不同的模型约束对基于正则化反演技术的地面核磁共振测深数据反演影响很大。利用理论数据讨论了平坦模型和光滑模型约束及与之相关的正则化参数对反演结果的影响规律。模拟表明,无论是高信噪比数据,还是低信噪比数据,反演的含水量在小正则化参数时,对模型约束的依赖性不大,主要受正则化参数影响。大正则化参数时反演分辨率低,而小正则化参数在保证数据拟合精度提高的同时,增加了模型的分辨率,能获得更可靠的含水量分布。 相似文献
3.
建立在常规的反演后分析数据分辨矩阵和模型分辨矩阵的基础上,自动选取最优数据进行自适应反演,以提高反演解释的准确性.另一方面,在电阻率变化相对剧烈的二维反演中,采用圆滑约束正则化反演方法经常导致反演不收敛或是结果不理想,为此提出加入自适应粗糙度参数α、β、γ.在圆滑约束反演的迭代改正过程中,同时确立突变边界位置,提高模型参数在此边界上电位梯度,以达到最佳的反演成像效果. 相似文献
4.
付德俊 《物探化探计算技术》2019,(3):365-370
开展了不同观测方式的井地2.5D直流电阻率反演研究。①从2.5维直流电阻率满足的边值问题出发,采用变分原理结合节点线性基函数推导了2.5维井地DC满足的积分弱解形式;②构建了二阶最大平滑稳定泛函的2.5维井地直流电阻率正则化目标函数,采用共轭梯度算法对正则化目标函数进行最优化求解,并采用逐步衰减正则化因子的求解策略来提高反演的稳定性;③设计了均匀半空间模型得到的数值解与解析解对比,电位的相对误差在2%以内,阐述正演算法的正确性和高精度。另外,分析了不同观测装置的2.5维井地直流电阻率异常体特征,并对不同观测方式对2.5维井地DC理论数据进行反演研究。研究结果表明,井中数据的引入提高了2.5DC对纵向探测的分辨率能力,同时提高了2.5DC反演有效性以及准确性。 相似文献
5.
《物探化探计算技术》2021,(1)
为提高电测深曲线反演过程的稳定性和自动性,基于广义交叉验证自适应正则化方法,给出了电测深曲线反演的全自动迭代方法。首先介绍了引入光滑约束的线性反演方程的构建方法,综合奇异值分解法与广义交叉验证自适应正则化方法求解线性反演方程,接着给出了反演过程中初始模型的构建方法和反演偏导数矩阵的差分计算方法,并采用C++语言实现了电测深一维反演算法,最后通过对模拟数据和实测数据进行反演试算,验证了反演算法的有效性 相似文献
6.
大定源瞬变电磁一维自适应正则化反演 总被引：1，自引：0，他引：1
自适应正则化反演对初始模型要求较低,直接对瞬变电磁响应数据进行反演。正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应的得到,从而快速获得地下的地电结构。本文采用自适应正则化反演算法对大定源回线瞬变电磁一维层状模型进行反演,使用均匀半空间作为初始模型,首次采用模型对深度的二阶导数极小的模型约束,通过典型理论模型的反演计算,证明了TEM自适应正则化一维反演算法拟合效果好、精度高,由于不用反复搜索正则化因子,并且收敛速度快,体现了良好的稳定性和可靠性。 相似文献
7.
8.
直流电阻率测深中二维与三维反演结果的对比与分析 总被引：1，自引：0，他引：1
直流电阻率测深中,用二维反演程序对三维地质体进行了反演,并与三维反演结果进行对比和分析。首先,对比二维和三维最小二乘反演在正演模拟算法和先验信息的确定方法;然后,对若干比较典型的模型进行反演实例对比。由于二维程序反演仅考虑单个剖面的电阻率信息,无论是在异常位置、形态及电阻率特性上,其反演精度都比较低;三维电阻率反演综合了多个测深剖面的电阻率信息,其结果与实际模型吻合得非常好。 相似文献
9.
《物探与化探》2021,(5)
长导线源半航空瞬变电磁正则化反演正则项通常采用L2范数,其拟合结果较光滑,不能有效刻画层界面信息。针对层状介质陡变模型实现正则项为L1范数的反演算法,采用迭代重加权最小二乘法将原问题转化为L2正则化子问题求解,解决L1范数存在不可导问题;采用OpenMP技术对雅可比矩阵并行计算,提高了反演速度;对自适应正则化因子分段迭代法的调整策略进行分析并改进,改进后的自适应正则化因子调整策略更适合半航空瞬变电磁L1正则反演算法。最后对电阻率进行反演并与Occam反演结果作比较,结果表明L1正则反演充分迭代后能够突出符合真实模型的电性界面,反演电阻率与模型真实值更接近。 相似文献
10.
在中心回线瞬变电磁反演过程中,为保证反演结果横向连续性,可施加横向约束,将多个测点组合成一个整体进行反演,在相邻测点之间建立相关联系,从而保证反演结果的连续性。本文横向约束使用横向加权、垂向粗糙度约束来增强相邻测点之间联系,提高反演精度。在反演过程中引入自适应正则化因子,正则化因子通过计算每次迭代的数据目标函数和模型目标函数的范数自适应得到,可改善反演过程中对初始模型的依赖。最后使用该算法对基于非结构时间域有限元三维正演数据及实测数据进行反演,与实际结果相符,表明该算法可增强反演结果的连续性且具有很好的抗噪能力。 相似文献
11.
A direct inversion scheme for deep resistivity sounding data using artificial neural networks 总被引：1，自引：0，他引：1
Initialization of model parameters is crucial in the conventional 1D inversion of DC electrical data, since a poor guess may
result in undesired parameter estimations. In the present work, we investigate the performance of neural networks in the direct
inversion of DC sounding data, without the need ofa priori information. We introduce a two-step network approach where the first network identifies the curve type, followed by the
model parameter estimation using the second network. This approach provides the flexibility to accommodate all the characteristic
sounding curve types with a wide range of resistivity and thickness. Here we realize a three layer feed-forward neural network
with fast back propagation learning algorithms performing well. The basic data sets for training and testing were simulated
on the basis of available deep resistivity sounding (DRS) data from the crystalline terrains of south India. The optimum network
parameters and performance were decided as a function of the testing error convergence with respect to the network training
error. On adequate training, the final weights simulate faithfully to recover resistivity and thickness on new data. The small
discrepancies noticed, however, are well within the resolvability of resistivity sounding curve interpretations. 相似文献
12.
Investigating 2-D MT inversion codes using real field data 总被引：1，自引：0，他引：1
Reza Ghaedrahmati Ali Moradzadeh Nader Fathianpour Seong Kon Lee 《Arabian Journal of Geosciences》2014,7(6):2315-2328
There are currently a significant number of two-dimensional (2-D) and three-dimensional (3-D) inversion codes available for magnetotelluric (MT) data. Through various 2-D inversion algorithms suggested so far, the classical Occam's inversion, the data space Occam's inversion, the nonlinear conjugate gradient (NLCG) method, and the Gauss–Newton (GN) method are fundamental driving methods to find optimum earth models, and OCCAM, DASOCC, NLCG, and MT2DInvMatlab are possible candidates one can find in the public domain that implement these algorithms for 2-D MT inversions, respectively. In this study, we investigate the pros and cons (strength and weakness) of these codes to help one use them efficiently in practical works and, as an introductory guide, further develop (sophisticate or extend) them, especially for the 3-D case. To achieve this goal, we applied each one of the four aforementioned codes on a profile of real MT field dataset. Then, further investigations have been done by performing several inversion tests to see how each code can find the appropriate model to reconstruct the subsurface resistivity structure. Numerical experiments show that the two parameters, regularization and target misfit, in addition to the main criteria of inversion (such as the forward and the sensitivities calculation method, and the type of inversion algorithm), are very important to produce the expected model in inversion. The regularization parameter that acts to trade off between model norm and data misfit can affect the inversion process in terms of both the computational efficiency and the accuracy of the obtained model. Also, lack of insufficient precision to choose the target misfit can lead the inversion to produce and reach an incorrect model. 相似文献
13.
为消除TEM定源回线装置增大接收区域后明显的边框效应,且提高TEM的分层能力,达到提高对不同电性地层精细分辨的目的,采用带约束的正则化反演方法。在反演算法中,拉格朗日因子采用CMD正则化因子的自适应调节方案,自主调整。为保证反演过程的高效、稳定且结果收敛性,使反演模型每一层的厚度按照等对数间隔进行离散,减少矩阵方程的病态性;同时对反演电阻率值采用可行的上下限约束,使迭代的结果不至于出现极值或错值,降低反演的多解性。通过对多层理论模型的计算,证明了约束反演算法的正确性和可靠性,将其用于实测瞬变电磁法勘探数据的反演计算中,反演电阻率纵向分布规律与已知地层电性资料吻合,解释的异常区得到了钻探验证,结果可靠。 相似文献
14.
The effectiveness of inversion apparent resistivity data to determine accurately the true resistivity distribution over 2D structures has been investigated using a common inversion scheme based on smoothness-constrained nonlinear least-squares optimization with enhancing horizontal resolution (EHR) technique by numerical simulation. The theoretical model generates in RES2DMOD software at specific distance and depth using Wenner, Wenner–Schlumberger, and pole–dipole arrays were inverted. The inversion model was compared with the original 2D model in RES2DINV software. The study model includes horizontal layering, vertical resolution, and horizontal two layers with different resistivity. Also, the response to variations in data density of these arrays was investigated. The study shows the best array suitable to be used in the survey was chosen for real data acquisition at the actual site. Subsequently, the results from borehole were used to verify the results of 2D resistivity imaging method with and without EHR technique. Saturated zone (0–40 Ω-m) was found scattered at the depth of 10–20 m. The borehole is located at 63 m at 2D resistivity imaging survey which shows at depth 10–20 m is sandy silt. Highly weathered sandstone was found at 6 m depth with resistivity value of 800 Ω-m and SPT N value of 20. The bedrock was found at 27 m depth with resistivity value of 3,000 Ω-m and SPT N value of 50. The application of 2D resistivity imaging with EHR technique indicate the ability of the proposed approach in terms of density, depth, and resistivity value of anomalous and layer in a computationally and numerically efficient manner and to exhibit good performance in the data inversion. 相似文献
15.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明：基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。 相似文献
16.
为了解地下介质物理特性随时间的变化,采用地球物理监测数据的成像方法已经逐渐得到应用。一种常用的途径就是通过独立反演不同时刻采集的地球物理数据,来重构间歇影像,这使得实测的数据信息得不到充分利用。因此,这里提出了一种同时反演所有不同时刻观测数据的新方法,即通过一次反演过程来获得所有不同时刻的地下模型。在反演过程中,通过引入时间和空间正则化因子来改善反演的稳定性,减少反演的多解性。最后,以直流电阻率监测数据反演为例,验证所提算法的有效性。 相似文献
17.
可控源音频大地电磁法（CSAMT）正演具有计算耗时、精度低的特点;常规的反演方法边界处理简单,占用内存大,搜索方向也难以精准,严重影响反演效果。为了实现CSAMT快速高精度正演,改善边界处理方法以减少内存占用,优化搜索方向以减少迭代反演次数,提出基于边界约束有限内存的拟牛顿CSAMT一维反演方法。首先,在CSAMT正演中基于高斯积分与并行算法实现了一维介质模型条件下有限长导线源激发的水平电场快速高精度计算。其次,在反演方法中,分别构建了基于相对误差和绝对误差的两种目标函数,并引入光滑模型约束,采取自适应正则化策略更新正则化因子,并采用基于边界约束的有限内存拟牛顿算法（LBFGS_B算法）,实现了CSAMT一维准确快速反演。以含低、高阻薄层的多层地电模型的水平电场Ex振幅为反演数据进行直接反演和模型检验。结果表明,该方法能够反映出地层电阻率随深度变化的趋势,并对高、低阻薄层有较好的分辨能力,低阻薄层分辨效果更佳。实际CSAMT资料的反演剖面的分辨率优于常规的连续介质反演方法,地质效果明显,表明该方法具有良好的应用前景。 相似文献
18.
瞬变电磁法（TEM）以其对低阻体反映灵敏而被广泛用于煤矿水文地质勘探中,当煤层埋藏较深,上覆地层纵向上有多种电性层存在时,瞬变电磁法晚期经验公式计算结果已不能反映实际情况。通过深入研究OCCAM反演算法,认为选择使误差目标函数达到最小的正则化参数,在反演迭代步长过大时进行压缩迭代步长的控制方法,以及最光滑的模型约束,可以提高反演分辨率。在瞬变电磁数据处理中使用加入了上述参数、控制和约束的反演技术,经对理论模型数据进行反演,反演结果与模型的电性分布规律吻合,表明算法正确;在华北煤田深部采场的瞬变电磁法探测应用中,OCCAM反演分辨出了煤层上覆多种电性地层,效果较好。 相似文献