共查询到20条相似文献,搜索用时 218 毫秒
1.
《物探化探计算技术》2021,(1)
为提高电测深曲线反演过程的稳定性和自动性,基于广义交叉验证自适应正则化方法,给出了电测深曲线反演的全自动迭代方法。首先介绍了引入光滑约束的线性反演方程的构建方法,综合奇异值分解法与广义交叉验证自适应正则化方法求解线性反演方程,接着给出了反演过程中初始模型的构建方法和反演偏导数矩阵的差分计算方法,并采用C++语言实现了电测深一维反演算法,最后通过对模拟数据和实测数据进行反演试算,验证了反演算法的有效性 相似文献
2.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。 相似文献
3.
4.
5.
针对三维地电场正演数值计算过程中形成的超大规模稀疏线性方程组,在分析此类线性方程组的一般解法基础上,着重阐述一种适宜求解此类方程组的Lanczos迭代过程与算法原理。同时,当地下介质的电性差异较大时,形成系数矩阵A的条件数就很大,可对算法进行适当改进。讨论采用不完全Cholesky分解方法进行预条件处理,经过条件数改善后,形成新的线性方程组系数矩阵,就会变为一个近似的单位矩阵。经改进后的Lanczos算法,将提高数值计算稳定性,从而加快迭代收敛速度,为提高反演质量提供基础。 相似文献
6.
7.
8.
由于离散运算,Radon变换方法中存在截断效应,传统的最小二乘约束反演并不能有效抑制这种效应。介绍了一种频域空间稀疏约束算法,在反演迭代过程中,根据前一次迭代的结果,通过贝叶斯原理将加权矩阵与前一次迭代的结果联系起来,得到新的加权矩阵;然后求解这个加权矩阵方程,得到频率域的稀疏解。比较了用最小二乘反演和压制截断效应的Radon域变换的效果,后者效果明显优于前者。模拟了一个5层的水平地层的地质模型,用这种压制截断效应的Radon变换进行了VSP资料的上下行波场的分离,得到了较好的效果。 相似文献
9.
电阻率三维复杂结构的快速反演 总被引:3,自引:1,他引:2
运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。 相似文献
10.
这里介绍一种加速遗传算法,该算法通过提高计算效率改善了标准遗传算法,加速遗传算法在速度反演应用中具有更高的实用价值。加速遗传算法在每次加速循环中进行一次迭代,用部份最佳个体的边界值优化搜索范围,大大地提高了计算效率。在迭代完成后,进行混沌优化以提高反演精度。二种算法的反演计算结果说明,加速遗传算法性能优于标准遗传算法,无论针对简单速度模型还是复杂速度模型,加速遗传算法都是适用的。 相似文献
11.
12.
13.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引:1,自引:0,他引:1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
14.
V.V. Plotkin 《Russian Geology and Geophysics》2018,59(3):324-333
The paper presents the results of a joint inversion of magnetotelluric sounding (MTS) curves taken at several stations without preliminary selection and normalization and distorted by the presence of lateral electrical-conductivity inhomogeneities in the medium. In the calculations, we used synthetic MTS data for a three-dimensional model. Preparation and interpretation of data are carried out by the Trefftz method using a numerical model of the field and MTS curve distortions. To solve the inverse problem and optimize the subsurface model, we used a nonlinear least-squares method and an iterative process with calculation of the sensitivity matrix and its singular decomposition. The target functional is determined by the discrepancies between the model and synthetic experimental apparent-resistivity curves corresponding to the elements of the impedance tensors on the lateral diagonal. The reliability of the reconstructed subsurface model is characterized by the dispersion of the deviations of its parameters from the parameters of the known model used for the preparation of synthetic experimental data. The joint consideration of distorted apparent-resistivity curves at several stations increases the reliability of interpretation results. The obtained solution to the inverse problem is approximate and can be used as a starting model for more complex algorithms and programs. 相似文献
15.
16.
为消除TEM定源回线装置增大接收区域后明显的边框效应,且提高TEM的分层能力,达到提高对不同电性地层精细分辨的目的,采用带约束的正则化反演方法。在反演算法中,拉格朗日因子采用CMD正则化因子的自适应调节方案,自主调整。为保证反演过程的高效、稳定且结果收敛性,使反演模型每一层的厚度按照等对数间隔进行离散,减少矩阵方程的病态性;同时对反演电阻率值采用可行的上下限约束,使迭代的结果不至于出现极值或错值,降低反演的多解性。通过对多层理论模型的计算,证明了约束反演算法的正确性和可靠性,将其用于实测瞬变电磁法勘探数据的反演计算中,反演电阻率纵向分布规律与已知地层电性资料吻合,解释的异常区得到了钻探验证,结果可靠。 相似文献
17.
《Russian Geology and Geophysics》2015,56(9):1357-1365
Computer inversion of measurement data based on the equation of particle transport is a topical component of the modern nuclear-geophysical technologies. In this work we have elaborated an approach and iterative methods for this inversion according to the gamma-gamma method data. In the framework of the “successive approximations according to characteristic interactions” approach, characteristic interactions have been established and iterative methods for the solution of two problems have been developed, namely, the problem of evaluation of the formation density and the problem of evaluation of two parameters—the formation density and the mud cake thickness. The numerical experiments confirm the theoretical concepts on the properties of iterative processes and convergence of methods. 相似文献
18.
岩石物理建模是叠前反演前重要的基础工作。以致密砂岩储层为目标,剖析了岩石物理建模涉及的主要过程。将测井评价与岩石物理建模视为迭代过程,测井评价提供岩石物理建模的输入数据,岩石物理建模控制测井评价参数的质量。在测井参数具有较高质量的情况下,优选岩石物理建模方法,优化关键参数,较为准确地预测了横波速度等弹性参数,为叠前弹性参数反演提供了可靠的基础数据。利用岩石物理解释模板对叠前反演纵横波速度体进行了解释,提高了致密砂岩气层的识别精度,在鄂尔多斯盆地取得了较好的应用效果。 相似文献
19.
提出了一种适用于大回线源瞬变电磁法框内任意测点垂直Z分量的一维反演方法。该算法采用CMD正则化因子优化方案、光滑度函数和自然边界条件对反演过程进行了约束,减小了初值的依赖度,保证了反演过程的自主性和反演结果的稳定性;在反演过程中,在保证模型残差沿梯度下降的同时,对模型修正量的最大值进行了限制,使反演模型不易陷入局部极小,保证了反演结果全局的收敛性;同时反演的灵敏度矩阵采用牛顿法进行更新,减少了整个反演计算的工作量。最后通过五层的双低阻层模型、双高阻层模型和实测资料证明了该算法可靠。 相似文献