logo

您所在位置网站首页 > 海量文档  > 学术论文 > 自然科学论文

直流电阻率三维正演的预条件共轭梯度法研究.pdf 1页

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
·710· 中国地球物理2011 (19)地质调查和矿产勘查地球物理 直流电阻率三维正演的预条件共轭梯度法研究 3一Ddirectcurrent forward the resistivitymodelingby method conjugategradient precondition 杨金凤+邓居智陈辉 JuzhiChenHui Jinfeng+Deng Vang 东华理工大学放射性地质与勘查技术国防重点科学实验室 江西抚州 344000 1.引言 直流电阻率法因其具有灵活性高、勘探费用低等特点被广泛的应用于矿产资源勘查、地质调查等领域。 如何快速、准确的求解离散后的线性方程组是直流电阻率三维正演的核心问题,对三维反演等后续工作具 有重要意义。本文基于有限差分的数值计算方法,引入对称超松弛预条件共轭梯度(SSOR.CG)法,研究直 流电阻率三维有限差分正演问题,并通过异常场法来消除总场中由点电源导致的奇异性。 2.方法原理 为了提高计算精度和计算速度以及避免由源奇异性引起的差分计算复杂性,采用解析计算和数值计算 相结合的方法,即异常场法,推导出不含电流源项的异常场方程。对求解区域进行六面体剖分,然后在网 格点上利用七点差分格式对变形后的积分方程进行离散得到节点电位满足的线性方程组。边界条件采用第 一类边界条件和混合边界条件。 对该方程组的系数矩阵采用对称对角行存储格式进行存储,即用两个一维数组存储。其中一个实型数 组存储矩阵上三角中的非零元,另一个整型数组存储上一个数组中的元素在矩阵中的位置索引。然后采用 对称超松弛预条件共轭梯度(SSOR.CG)法求解。预处理技术的关键是通过合适的预调件矩阵的选择,来降 低系数矩阵的条件数。本算法首先将系数矩阵分解成一个上三角矩阵和下三角矩阵及单位矩阵,然后将预 条件矩阵定义为含松弛因子的一个上三角阵与其转置的乘积。对于松弛因子的选择,根据一定的迭代精度 试算了不同模型最终确定了一个最佳的松弛因子。对于该算法的另一关键步是通过高斯消去法先顺带后回 带最终可求出预处理矩阵与任一列向量的乘积,该计算过程仅涉及少量非零元素的操作,均可以非常简单、 快速地完成。 3.模型验证及结论 首先设计了水平层状介质模型和垂直接触模型进行试算,结果表明采用对称超松弛预条件共轭梯度法 求解有限差分正演中的大型线性方程组,与高斯直接法相比速度显著得到提高,而且内存需求大大减少, 验证了该算法的可行性与有效性。然后设计了典型的高阻(低阻)正方体异常模型进行计算,结果表明将 解析计算和数值计算相结合并利用对称超松弛预条件共轭梯度法求解该正演计算中的异常电位,可以消除 点电源奇异性影响,在保持快速计算的基础上,显著提高了计算精度。 本研究由“东华理工大学研究生创新基金(DYCAl0015)”资助。 参考文献 M P removal:arefinementof 【1】LowryT’AllenB,ShiveN.Singularity resistivitymodeling 1989,54(6):766.--77.

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

“原创力文档”前称为“文档投稿赚钱网”,本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有【成交的100%(原创)】。原创力文档是网络服务平台方,若您的权利被侵害,侵权客服QQ:3005833200 电话:19940600175 欢迎举报,上传者QQ群:784321556