国家高技术研究发展计划(2012AA09A20105)
- 作品数:4 被引量:33H指数:3
- 相关作者:戴世坤张钱江陈龙伟李昆赵东东更多>>
- 相关机构:中南大学中国石油大学(北京)更多>>
- 发文基金:国家高技术研究发展计划国家自然科学基金更多>>
- 相关领域:天文地球更多>>
- 频率域可控源电磁法2.5D正反演被引量:19
- 2013年
- CSEM有MT不可比拟的优势,得到国内外学者多年的大力发展,但是,其数据处理水平尚停留在1D和2D近似解释阶段。2.5D相对来说更加符合实际情况(3D),要想进一步发展该方法,必须开展CSEM 2.5D或3D正演和反演研究。采用有限元法实现了CSEM 2.5D正演模拟,在验证其代码有效的前提下,计算了水平地形和地堑地形的电磁场响应。采用共轭梯度法实现了对场值直接进行反演的算法。从试算低阻和高阻模型的结果可知,反演算法精度高,计算速度快,初步达到了实用化的程度。
- 戴世坤王顺国张钱江薛东川
- 关键词:正演反演共轭梯度
- 起伏地形条件下二维声波频率域全波形反演(英文)被引量:1
- 2015年
- 频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和HeSSian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。
- 张钱江戴世坤陈龙伟李昆赵东东黄兴兴
- 关键词:起伏地形有限单元法
- 基于高斯牛顿法的二维直流电阻率法的快速反演被引量:6
- 2015年
- 从二维线源问题出发,对二维直流电阻率法高效、高精度正反演方法进行研究。在正演数值模拟中,引入直接解法求解器求解线性方程组,既保证了起伏地形条件下有限元法正演数值模拟的计算精度和计算效率,又为反演算法中"拟正演"快速回代求解提供了条件。结合高效、高精度的正演算法,采用高斯牛顿法对电阻率进行反演成像。在弱非均匀介质前提下,基于近似海森矩阵主对角线元素严格占优的特点,采用舍弃海森矩阵非对角线元素的策略,提高整个反演计算的效率。最后,利用合成数据对反演算法的有效性进行检验。结果表明:给出的反演算法稳定、快速,结合偶极-偶极装置和三极装置,能有效地反演出异常体的形状、大小和位置。
- 赵东东张钱江戴世坤陈龙伟李昆
- 关键词:直流电阻率法起伏地形反演
- 层状单轴各向异性介质中长导线源电磁响应的快速计算被引量:7
- 2016年
- 在场源附近,长导线源不能视作电偶极子源。长导线源电磁响应的常用计算方法是对场源进行离散,再对离散电偶极子电磁场进行叠加,这种方法计算简单,但计算效率较低。针对这一问题,提出了有限单元形函数积分算法。该算法从电偶极子场的解析表达式出发,将长导线源进行剖分,通过在每个剖分单元上构造形函数来求解积分表达式。通过模型算例对传统电偶极子叠加算法和形函数积分算法的计算精度及计算效率进行了测试。从测试结果可以看出,无论在场源中远区还是近区,形函数积分算法只需要将长导线源剖分成两个单元就能达到很高精度,而电偶极子叠加算法需要将长导线源离散成大量的电偶极子才能达到相同数值精度,相同数值精度条件下形函数积分算法具有更高的计算效率,且离场源越近计算效率优势越明显。
- 欧阳芳戴世坤张钱江陈龙伟陈轻蕊