电阻率法和激发极化法不同装置组合数据三维非线性共轭梯度反演并行算法研究

来源 :中国地质大学(北京) | 被引量 : 0次 | 上传用户:ffgghhaz
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
直流电阻率法和激发极化法在浅层水文、找矿、工程、环境等领域应用广泛,而浅层的地质目标体一般表现为三维电性结构,一维和二维的数据反演不能完整的恢复三维电性结构,为了避免在反演中三维体不均匀性的影响,需要采用三维数据观测和反演。在找矿方面,基于矿山中钻井的相对位置,我们可以将井中装置与地表装置相结合,即形成地表、地井、井地和井间四种观测方式相结合的复合采集方式。由于三维地球物理问题不仅是欠定的、非线性的,在进行三维反演时还存在耗时长的特点,而且线性反演算法往往使得反演陷入局部极小,所以本文采用非线性共轭梯度反演算法来规避使用线性反演时所带来陷入局部极小的问题,并且在该反演算法中加入MPI和GPU并行算法来提高计算效率。针对直流电阻率法和激发极化法的正演问题,本文采用有限差分法求解点电源三维地电场。针对正演算法中的解方程,采用不完全Cholesky分解的共轭梯度法和不完全LU分解的稳定双共轭梯度法;针对大型系数矩阵,本文采用行压缩存储的格式,不仅减小了内存消耗,还提高了计算的效率。反演算法采用非线性共轭梯度法,非线性共轭梯度法的优势就在于不用求雅克比矩阵,而是直接求取雅克比矩阵或者其转置和一个向量的乘积。该方法不仅避免了存储雅克比矩阵,还大大减少了计算量,每次反演迭代计算一次正演,两次“拟正演”。在并行计算方面,由于MPI含有将多个独立任务分配给多进程同时计算的特性,即可利用MPI进行多个点电源地电场的计算。而在计算正演和“拟正演”过程中,都存在解方程的过程,反演过程的大部分时间都消耗在解方程中,因此,利用GPU快速运算的特性,在解方程时,调用基于GPU开发的CULA Sparse函数库参与解方程运算。CULA Sparse函数库包含多种解方程方法和预处理方法,在本文中,我们采用不完全LU分解的共轭梯度解方程。最后,分别设计单一采集装置和多个将地表、地井、井地和井间相结合的装置,通过模型算例验证实现的非线性共轭梯度并行算法的有效性,并说明了多装置相结合采集的数据参与反演比使用单一装置采集的数据参与反演的结果更好。
其他文献
本文将NGDC-720全球卫星岩石圈磁场(LMF)模型应用于中国大陆及邻区,探讨了不同分量的磁场及其梯度在不同高度上的分布特征;将该区域地形、活动构造、地壳厚度、地表热流等与LM
针对独立学院在发展过程中内部控制制度建设存在的问题及构建进行探讨,结合内部控制相关理论提出关于独立学院控制环境、控制活动、风险评估、监督与评价等方面的建议。
对“什么是比较文学”这一关乎比较文学存在依据的重要命题,韦斯坦因超越了重视考证的法国学派和眼界宽泛的美国学派之间的争辩,采取了一种中间路线,他从辨析定义出发,主张回溯比
以甲基丙烯酸羟乙酯和肉桂酰氯为原料,合成了甲基丙烯酰氧乙基肉桂酸酯(CEMA),进一步选用甲基丙烯酸二甲氨基乙酯(DMAEMA)和CEMA为主要单体制备了系列共聚物水凝胶。通过红外
由于岩体中包含了数量众多、分布形式复杂的结构面,从而使岩体力学性质具有尺度效应和高度的不均匀性、不连续性和各向异性。岩体的力学参数取值一直是岩体工作者关注的热点
针对肿瘤放射治疗中的加速器机械碰撞问题,提出一种增强现实环境下基于Unity3D的碰撞检测方法。基于三维建模软件3ds Max对医用直线加速器、治疗床等设备进行准确建模,基于病
目的对肠造口患者进行为期3个月的延续性护理,对3个月内患者自我管理能力的影响因素进行分析;对护理人员在肠造口相关知识、健康宣教、临床实践能力等方面进行分析,描述目前
(一)十五大报告明确指出公有制和公有制的实现形式是两个不同的概念。“公有制经济不仅包括国有经济和集体经济,还包括混合所有制经济中的国有成分和集体成分。”这是从企业
ue*M#’#dkB4##8#”专利申请号:00109“7公开号:1278062申请日:00.06.23公开日:00.12.27申请人地址:(100084川C京市海淀区清华园申请人:清华大学发明人:隋森芳文摘:本发明属于生物技
新媒体传播技术的发展大大改变了媒体和受众之间的关系,使得信息传播呈现出虚拟与现实融合、时空互补的优势,在传播过程中人的价值不断提升,人本媒体的趋势更为明显。本文就