邹桂红
- 作品数：3 被引量：2H指数：1
- 供职机构：中国石油大学（北京）地球物理与信息工程学院更多>>
- 发文基金：中国石油大学(北京)科研基金更多>>
- 相关领域：理学天文地球更多>>
- 迭代法中压缩对角存储的应用框架被引量：1
- 2013年
- 在许多有限元、有限差分的应用中,待解方程组的系数矩阵是大型稀疏带状阵,方程组的求解一般使用迭代法。与其它存储方式相比,压缩对角存储由于不存矩阵元素的行列索引,对内存的使用最为节省。使用压缩对角存储,沿对角线操作可以完成高效的矩阵—向量乘。由于在以往的文献中,没有提及按行、列操作的算法,压缩对角存储的应用范围受到一定的限制。用行、列、对角版的矩阵—向量乘代表普通意义的行、列、对角方向操作模式,通过等价矩阵推导和伪代码,给出了一个包括按行、列操作算法在内的应用框架。这里运用C语言实现不完全Chol-esky分解共轭梯度法解方程,阐述了如何在实际编程中使用这个算法框架。经理论与实验分析表明,对角压缩存储应用于框架中是高效的,因为与使用常用的一维行索引存储格式所编程序相比,同样迭代次数的耗时减少了约25%。
- 邹桂红梁华庆尹洪东
- 关键词：迭代法
- 基于四叉树网格的2.5D直流电阻率法有限元数值模拟被引量：1
- 2014年
- 在VC环境下,使用C语言实现了基于四叉树网格的2.5D直流电阻率法有限元数值模拟程序。3个典型模型的算例被用来证明四叉树网格方案适合于复杂地电模型的数值模拟:分别使用四叉树网格和等距全网格对水平两层模型进行了计算。后一方案与前一方案相比,在数值模拟精度上没有明显的差别,而使用的网格数却有大幅度的减少;计算了覆盖层下倾斜高阻板状体模型,高阻板状体的形态、埋深清晰地显示在偶极-偶极视电阻率断面图上;仿真了中间梯度法测量2D地形下埋藏直立高阻板状体的情形,利用比较法对计算结果进行地形改正后,高阻体的埋深被计算出来。
- 邹桂红梁华庆耿敏
- 关键词：有限元法电阻率法
- 基于八叉树网格的直流电阻率法有限元数值模拟
- 2013年
- 在VC环境下,使用C语言,实现了基于八叉树网格的三维直流电阻率法有限元数值模拟程序。三个典型模型的算例用来证明八叉树网格方案适合于复杂地电模型的数值模拟。数值模拟了球体模型。随着网格对球面边界的不断逼近,数值模拟的精度也逐步提高。计算了覆盖层下倾斜板状高阻体模型,在偶极-偶极视电阻率断面上,能清晰地看见高阻体模型的形态、埋深。仿真了联合剖面法测量三维地形下埋藏低阻椭球体的情况。由于地形的影响,在视电阻率曲线上不能发现低阻异常的存在,利用比较法进行地形改正后,低阻异常得到凸显。
- 邹桂红梁华庆耿敏
- 关键词：有限元法直流电阻率法
