基于OpenMP的重力张量并行正演     

余涛  朱自强  鲁光银  曹书锦 《物探化探计算技术》2013,(4):446-451,372

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路.  相似文献   

花岗岩地区浅部地下水井位的快速确定          下载免费PDF全文

杨天春  张正发  许德根  曹书锦 《水文地质工程地质》2017,(5):20-20

文章探讨在花岗岩地区地下水井位的快速定位问题,结合工程实践,采用三极电阻率测深法和天然电场选频法开展综合探测研究。首先实施三极电阻率测深剖面工作,对采集结果进行二维有限元反演研究,获得二维测深断面上低阻导电体的分布规律;在测深剖面上进一步开展天然电场选频法剖面勘探工作,得到低阻含水构造的水平精确位置,并进一步开展选频测深法的实践。后期的成井结果表明:三极电测深勘探成果的剖面异常比较直观,选频法则在异常体的地表水平定位方面比较精确,施工非常方便,两种方法相结合对地下水进行综合勘探是一种行之有效的办法。  相似文献   

基于单元细分H-自适应有限元全张量重力梯度正演   总被引：1，自引：0，他引：1  

曹书锦  朱自强  鲁光银  曾思红  郭文斌 《地球物理学进展》2010,25(3)

在重力梯度各阶张量正演中,难于求取复杂地质体的解析式,转而利用已知解析解的简单规则几何形体近似剖分地质体,其会引起较大的几何拟合误差.本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元分析给出了满足地下复杂地质体正演的型函数,引入体积分的误差指示器用以估计全张量重力梯度误差,提出了适用于全张量重力梯度正演的自适应计算的策略及迭代算法,数值算例证明了本文方法的正确性和有效性.  相似文献   

地质雷达目标体体积计算   总被引：1，自引：0，他引：1  

朱自强  曹书锦 《物探化探计算技术》2009,31(2):117-120

采用地质雷达探测地质体,很难直接从探测图像资料上获取目标体体积。这里基于频域雷达方程,提出了一种根据地质雷达探测影像宽度,电磁波长,以及目标体埋深等参数来计算散射截面,进而求取被测目标体等效体体积的计算模型,并将其应用于工程实例。计算结果表明,该模型准确可靠,是一种快速求取目标体体积的新思路。  相似文献   

基于粒子群算法的重力张量反演研究     

孙中科  朱自强  鲁光银  曹书锦 《物探化探计算技术》2013,35(2)

重力梯度张量的定义是对重力位求二阶导数,相比于传统的布格重力异常,它在反映地下密度异常分布上有着更高的灵敏度,且能够进一步准确直接地反映目标体的边界.但对单个张量分量作反演时,可能会丢失一部份有用的信息,从而造成反演结果的误差.而全张量反演是将重力的五个梯度张量联合起来进行反演,这样做可以综合更丰富的场源信息.相比于传统做法中的布格重力异常反演和单重力张量分量反演,得到的反演结果不仅有了更高的分辨率,在识别目标体特征上也有更好的效果.粒子群算法是一种基于群体智能的优化迭代算法,这里利用粒子群算法对重力张量单分量、布格重力异常和全张量分别进行反演,并对结果进行简要分析.  相似文献   

基于OpenMP的重力张量并行正演CSCD     

余涛朱自强鲁光银曹书锦 《物探化探计算技术》2013,(4):446-451

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。  相似文献