共查询到20条相似文献，搜索用时 15 毫秒
1.
随着重力梯度全张量测量技术的日趋成熟,重力梯度全张量数据的三维反演技术日益受到重视与关注.全张量数据反演与重力数据反演一样仍然面临着严重的多解性问题.本文将基于地质统计学的协同克里金方法应用于重力梯度全张量数据三维反演,建立了密度约束下的多变量协同克里金联合反演方程,以降低反演的多解性.模型试验表明密度信息的加入能够有效降低反演的多解性,提高反演结果的分辨率,尤其是纵向分辨率能够得到显著提高.最后对美国德克萨斯州一个岩丘区所获得实际资料的应用表明了本文方法的实用性. 相似文献
2.
基于GPU并行的重力、重力梯度三维正演快速计算及反演策略 总被引：4，自引：2，他引：2
利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度Uxx、Uxy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴. 相似文献
3.
随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性. 相似文献
4.
《应用地球物理》2017,(2)
随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU(Graphics Processing Unit)和预处理分解技术(Preconditioning methods)在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATesla C2050 GPU的并行SSOR-ICCG算法和在2.0GHz CPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。 相似文献
5.
重力梯度数据相对于传统重力数据,能够更细致、准确地描述地球浅部构造和研究矿产资源分布等信息.本文采用共轭梯度算法,在加权密度域求解重力梯度数据三维聚焦反演最优化问题,以恢复地下三维密度分布,目标函数包括数据不拟合函数和最小支撑稳定函数.首先,在推导目标函数对加权密度的一阶导数时,为了得到更合理的计算公式,我们考虑变加权函数中含有密度变量;此外,本文通过密度上下限约束,改善了传统聚焦反演中聚焦因子选取困难的问题.新算法获得的反演结果,对聚焦因子的选择约束较少,相比传统聚焦算法,能够更容易的获得理想结果.将方法应用于理论模型验证其有效性和正确性,并应用本文方法处理文顿盐丘地区的航空全张量重力梯度数据,得到了与已知地质信息匹配的密度分布,表明本文方法具有处理实际数据的能力. 相似文献
6.
7.
重力梯度全张量数据具有高灵敏度反映探测目标局部微弱变化的特点,可以对单一地震数据在建立速度模型过程中的不确定性,尤其是各向异性速度模型建立过程中的精细结构进行补充修正.本文主要研究重震两类数据与各自所反映的物理模型关系以及密度模型与速度模型之间的相互联系和制约关系,从而建立重力全张量数据与地震速度模型间的信息互补关系和模型参数间的解析表达式.在此基础上,提出多参数加权正反演算法流程,逐层异常分离技术和约束条件下迭代反演模型修正方法,并实现模块编程和验证.实验数据和结果表明,重力全张量数据和地震数据多参数融合能够明显提高地震速度建模精度,减小模型解释的不确定性. 相似文献
8.
9.
重力数据的物性反演面临着严重的多解性问题,降低多解性的有效手段是加入约束条件.而边界识别、深度估计及成像方法可获取地质体的水平位置、深度范围等几何参数信息,本文将基于数据本身挖掘的地质体几何参数信息约束到物性反演中,以降低反演的多解性.通过引入基于深度信息的深度加权函数及基于水平位置的水平梯度加权函数建立优化约束条件,有效地提高了反演结果的横向及纵向分辨率.重力梯度数据包含更多的地质体空间特征信息,将优化约束反演方法应用到全张量数据的反演中,模型试验表明本文方法反演结果与理论模型更加吻合.最后对美国路易斯安那州文顿盐丘实测航空重力梯度数据的应用表明,本文方法在其他地球物理、地质资料不足的情况下获得更可靠的反演结果. 相似文献
10.
航空电磁数据的三维解释由于数据量大需要有高效的反演算法作为支撑.本文利用两种目前主流的数值优化技术（非线性共轭梯度和有限内存的BFGS法）实现了三维频率域航空电磁反演,并进一步比较了两种方法的有效性和运算效率.在反演过程中,为了更好地反演异常体的空间位置,模型方差矩阵中的光滑系数在反演起始阶段取值较大;当数据拟合差下降趋于平缓时,再利用较小的光滑因子约束反演过程来实现聚焦和获得精确的反演结果.理论数据反演表明这两种优化策略具有相似的内存需求,但是有限内存的BFGS技术比非线性共轭梯度法在计算时间和模型反演分辨率上具有一定的优越性,因此有限内存BFGS法更适合于求解大规模三维反演问题. 模型试验进一步表明目前主流的迭代法求解技术不适合大规模航空电磁数据反演,未来移动平台多源电磁数据快速正反演可通过引入矩阵分解技术来实现. 相似文献
11.
起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性. 相似文献
12.
13.
重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合. 相似文献
14.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 相似文献
15.
大地电磁三维快速松弛反演并行算法研究(英文) 总被引：8，自引：1，他引：7
为了进一步提高大地电磁三维快速松弛反演的计算效率,在深入分析大地电磁三维快速松弛反演算法的基础上,结合MPI自身的优越性,确定了并行计算的思路,实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维快速松弛反演并行程序运行结果正确,效率提高明显。此成果在普通微机上实现,推动了大地电磁三维反演技术的实用化,可为其它地球物理三维正反演研究所借鉴。 相似文献
16.
本文实现了地面矩形大定源三维频率域反演.矩形大定源三维模型响应计算采用交错网格有限差分技术.正演的微分方程为异常电场满足的非齐次Helmholtz方程,方程右手边源项中的大定源产生的背景格林函数由虚界面法结合虚框法计算.频率域三维反演采用非线性共轭梯度反演技术.反演的数据类型为垂直磁场的频率域响应Hz的实部和虚部分量.数值结果表明,(1)三维模型正演模拟数值结果与前人一致,为三维反演奠定基础;(2)针对两个三维导电模型,分别进行了三维反演数值试算.反演结果可以清晰恢复出异常体的电阻率和位置信息,表明地面矩形大定源三维频率域非线性共轭梯度反演具有可行性.本文研究的意义在于,在电磁响应时频转换技术的基础上,如果将野外实测的瞬变电磁数据变换为对应的频率响应,则结合本文提出的三维反演技术,可以为矩形大定源瞬变电磁数据的三维解释提供一个新的思路. 相似文献
17.
本文利用Extrapolation Tikhonov正则化算法处理重力梯度数据三维密度反演的线性不适定问题。与Tikhonov正则化方法相比,Extrapolation Tikhonov正则化方法减小了因正则化参数的引入而带来的反演结果误差,提高了预测数据与观测数据之间的拟合精度。同时为了消除位场数据反演时位置函数快速衰减对反演结果的影响,本文提出了基于重力梯度全张量特征向量法的深度加权函数,模型试验证明了该深度加权函数能有效识别异常体密度分布特征。对澳大利亚Kauring地区实测重力梯度数据进行反演,并和已有研究成果对比分析。结果表明该反演方法能够较好的获取地下异常体的密度分布信息。 相似文献
18.
基于深度学习的重力异常与重力梯度异常联合反演 总被引：1，自引：0，他引：1
高效高精度的反演算法在重力大数据时代背景下显得尤为重要,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于深度学习的重力异常及重力梯度异常的联合反演方法.文中首先提出了一种基于网格点几何格架的重力异常及重力梯度异常的空间域快速正演算法,这为本文深度学习反演算法的实现奠定了基础;其次对大量的不同密度模型进行正演计算获得样本数据集;然后设计了一种端到端的深度学习网络结构(GraInvNet),再利用样本数据对该网络结构进行训练;最后进行反演预测.组合模型试验表明,多维度数据联合反演相比单一分量反演其结果更“聚焦”,且与模型边界高度吻合,并且对于复杂模型的姿态与物性预测具有极为显著的优势,以及对于含噪声数据的反演,其质量也不会降低;Vinton岩丘实测重力数据也验证了文中方法的有效性;从而证明了深度学习在重力数据的高效高精度反演方面具有的巨大潜力. 相似文献
19.
实现大地电磁法快速三维反演的关键在于找到快速计算灵敏度矩阵的方案. 本文在对大地电磁三维张量阻抗表达式进行深入分析的基础上，获得了三维快速松弛反演算法的灵敏度函数表达式，解决了三维快速松弛反演的核心问题；为了减小反演解的多解性，针对三维问题定义了最小构造函数，实现了求最小构造的三维快速松弛反演算法. 对二维棱柱体、三维棱柱体理论模型的大地电磁合成数据进行了反演试算，反演结果与理论模型相吻合. 对日本Kayabe地区实测资料进行了XY模式反演，反演得到的地电模型较好地反映了地热资源的分布. 反演试算结果表明大地电磁三维快速松弛反演具有计算速度快、结果稳定可靠等特点；使三维快速松弛反演在普通微机上得以实现. 相似文献
20.
常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好. 相似文献