共查询到17条相似文献，搜索用时 947 毫秒
1.
2.
《物探化探计算技术》2019,(6)
利用重力异常反演地球内部的介质密度分布,是研究地球内部精细结构的有效方法之一,这里利用EIGEN-6C2超高阶卫星重力场模型,计算龙门山地区的自由空间重力异常、布格重力异常以及剩余重力异常,在此基础上,基于Parker-Oldenburg模型以及粒子群算法反演该地区的三维介质密度分布,发现反演结果与Crust1.0密度模型所提供的结果基本一致,结果说明,使用Parker-Oldenburg模型结合粒子群算法反演地壳的三维密度分布,具有可行性和实用性。 相似文献
3.
利用共轭梯度方法的激发极化三维快速反演 总被引：2，自引：0，他引：2
利用共轭梯度方法实现了激发极化(IP)三维快速反演。首先,利用共轭梯度方法反演电位数据,得到地下的三维电阻率模型,(由于避免了直接求偏导数矩阵,因此反演计算速度很快。)然后,以此电阻率模型为地下电导率分布,再反演IP数据得到三维极化率分布理论模型。试算结果表明其效果较好。 相似文献
4.
《物探化探计算技术》2015,(6)
利用阻尼型高斯牛顿法求解反演目标函数,并将粗糙度矩阵约束和最小支撑泛函约束以及参考模型和模型参数界限约束引入到目标函数中,同时结合预条件共轭梯度法进行反演迭代,实现了激发极化数据三维聚焦反演快速计算。通过对几个典型模型的试算,并与传统光滑模型约束的反演结果进行比较,表明该反演算法的可靠性和稳定性较好,并且算法速度快,占用内存低且易于并行化。 相似文献
5.
6.
基于不同地球物理数据对地下结构不同的敏感性和数据分布,联合地球物理反演可以减少反演的非唯一性和提高反演模型的可靠性。在研究中,利用噪声成像得到的地震面波相速度,并结合当地的布格重力异常数据对四川地区的岩石圈速度结构进行了联合成像研究。地震面波频散数据主要对地下岩石的横波速度敏感,而重力数据对地下介质的岩石密度有很好的约束性,为了将两种数据归并统一到同一个联合反演系统中,利用了地震波速度和岩石密度之间的经验关系。基于面波和重力联合反演成像算法,得到了四川地区岩石圈的三维横波速度模型。该模型不仅与地表已观测到的地质特征有较好的吻合性,而且能够比较好地拟合面波数据和重力数据,新的模型对于研究四川地区的地震灾害有着重要的参考意义。 相似文献
7.
遗传算法是近些年来产生和发展的一种模拟生物进化过程的自适应启发式全局优化的搜索算法。它不完全依赖于初始猜测,且具有全局收敛的特点,可以被用来解决各种复杂的实际问题,如工程优化设计,人工智能和决策系统,以及地球物理反演等。尽管遗传算法是一种效率很高的全局优化算法,但许多仿真结果表明,它具有计算时间长,局部搜索能力弱的缺点。而共轭梯度法属于非启发式全局优化搜索方法,收敛速度快,但容易陷入局部极值,且严重依赖初始猜测。根据遗传算法和共轭梯度法的特点,这里提出了一种混合遗传算法,用来进行地球物理反演。该算法既具有遗传算法的全局收敛性,又有共轭梯度法的快速收敛性,经实际应用,取得了良好的效果。 相似文献
8.
基于粗细网格的有限元并行分析方法 总被引：2，自引：0，他引：2
并行计算己成为求解大规模岩土工程问题的一种强大趋势。探讨了粗细网格与预处理共轭梯度法结合的并行有限元算法。从多重网格刚度矩阵推得有效的预处理子。该算法在工作站机群上实现。用地基处理时土体强夯的数值模拟分析进行了数值测试,对其并行性能进行了详细分析。计算结果表明：该算法具有良好的并行加速比和效率,是一种有效的并行算法。 相似文献
9.
由于地球物理反演存在固有的多解性问题,如何减小多解性,提高地球物理方法反演结果的可靠性和准确性是当今地球物理领域研究的热点。为了减小大地电磁与地震反射波走时资料反演的多解性,在Gallardo和Meju研究的交叉梯度理论的基础上,实现大地电磁与地震反射波走时数据的二维联合反演算法;设计理论地电模型和速度模型进行合成数据的单独反演和联合反演试算,并验证大地电磁和地震反射波走时数据二维联合反演算法的正确性和有效性。反演结果表明,利用交叉梯度方法耦合大地电磁电阻率参数和地震波速度参数能使联合反演过程中的这两个模型在更新的时候相互约束,联合反演结果比单一反演结果更接近真实模型,从而使反演结果更可靠,减少了反演解的多解性。 相似文献
10.
电阻率三维复杂结构的快速反演 总被引：3，自引：1，他引：2
运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。 相似文献
11.
为提高重力及重力梯度数据的正演效率,笔者引入基于MPI (Message Passing Interface)和OpenMP (Open Multi-Processing)的并行计算,通过对比分析不同数据规模的网格数和模型体个数对并行效率和加速比的影响,得出随着正演数据规模的增加,并行效率和加速比均得到提高。同时对比了基于MPI和基于OpenMP的两种并行方式的性能,结果表明,重力及重力梯度数据正演的并行计算中MPI的提速能力优于OpenMP,且在较大规模数据的正演计算中基于MPI的并行效率优于基于OpenMP的并行效率。 相似文献
12.
病态线性方程组的简单迭代解法 总被引：9，自引：0，他引：9
针对地球物理反问题中经常碰到的病态线性方程组的求解问题。本文提出了一种简单迭代（ＳＩ）算法,从理论上证明了解序列收敛且收敛到方程组的真解,然后给出了几个算例,将计算结果与对付病态问题能力很强的ＣＧ类算法的结果进行了对比,结果表明：ＳＩ算法具有极强的抗病态能力,计算精度明显高于ＣＧ类算法,但计算速度稍低于后者 相似文献
13.
基于预优共轭梯度算法的反演方法在求解磁化强度大规模欠定线性方程组时可以改善方程组的条件数,并提高算法的纵向分辨率。该方法对磁性体形状无要求,且收敛速度快、稳定性好,能较直观地反映磁性体特征。本文通过理论模型试验结果分析证明了基于预优共轭梯度的磁化强度反演方法的稳定可靠性,北衙铁金矿区万硐山矿段L32线磁异常反演结果表明了该方法应用于实际资料处理与解释的有效性和可靠性。 相似文献
14.
We developed a three-dimensional(3D) conjugate gradient inversion algorithm for in-verting magnetotelluric impedance tensor measurements.In order to show the importance of including diagonal components of magnetotelluric impedance tensor in 3D inversion,synthetic data were inverted using the 3D conjugate gradient inversion,and the inversion results were compared and analyzed.The results from the 3D inversion of synthetic data indicate that both the off-diagonal and the diagonal components are required in in... 相似文献
15.
随着传感器技术的发展,重力场与重力张量场测量技术发展迅速,为实现地下密度分布精细反演提供了数据保障。正演是反演的基础,解决任意密度分布复杂地质体重力场与重力张量正演高效、高精度计算问题,是实现重力高效、精细反演、人机交互反演解释的关键。针对起伏地形和任意密度分布这种复杂条件下二维重力场及重力张量场高效高精度正演问题,这里提出了一种空间波数混合域正演算法,其关键环节包括:①结合新的矩形二度体组合模型波数域表达式和一维Gauss-FFT算法,提出了一种任意密度分布和起伏地形下重力场及重力张量高效、高精度正演算法;②采用新的二维正演算法,计算观测最高点和最低点之间多个不同高度水平网格重力场及重力张量,结合三次样条插值方法,实现了起伏地形上重力场及重力张量场高效、高精度正演。模型算例结果表明,新方法具有高效、高精度的显著特点。 相似文献
16.
随着岩土工程规模的不断扩大、复杂性的增加以及计算参数的多样化和计算精度的提高, 人们对于计算机计算能力的要求越来越高, 然而单处理器无法满足这类大规模计算.从数据输入、区域分解、线性方程组的迭代求解、后处理等方面详细阐述高性能计算平台上并行有限元求解大规模岩土工程的关键问题.提出了利用MPI2的新特性进行海量数据的分段并行读入, 采用ParMetis软件并行地进行区域分解, 实现了前处理过程的完全并行化; 采用基于Jacobi预处理技术的预处理共轭梯度法(PCG)进行线性方程组的并行迭代求解; 采用Paraview软件实现了后处理的并行可视化.在深腾7000系统上对某隧道工程的三维开挖过程进行了数值模拟, 对其并行性能进行了分析和评价, 验证了采用的区域分解算法和系统方程组的求解方法的可行性, 并且具有较高的加速比和并行效率. 相似文献
17.
抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。 相似文献