二维电阻率正演.pptVIP

地球物理正演程序设计 报 告 人: 陈辉 单 位: 东华理工大学 专 业: 地球探测与信息技术 2009年9月 地球物理正演 给定地质体的形状大小，空间位置及物性参数的条件下，在不同场源（天然源和人工源）的作用下，求取空间和时间上任意一点上的物理场值 主要方法 解析计算法 有限差分法 有限元法 边界元法 积分方程法 基本思想 有限代替无限 离散代替连续 求解线性方程组得到近似解 基本步骤 求解区域离散化 地球物理方程离散化 边界条件离散化 建立线性方程组 求解线性方程组 输出地球物理场近似解 关键问题 区域剖分 区域坐标的剖分，物性参数的剖分 线性方程组的建立 系数的求取，右端项计，算矩阵的存储 求解线性方程组 选取合适的方法 需要的知识 地球物理方法的基本理论 地球物理方程，边界条件 计算方法的基本理论 有限差分，有限元，积分方程 线性方程组的求解方法 松弛迭代，共轭梯度 实施步骤 公式推到 程序编写 结果验证 二维直流电阻率有限差分正演 公式推导 二维直流电法的基本方程 欧姆定律： 电位定义： 二维直流电法的基本方程（续） 边界条件 第一类边界条件（Dirichlet） 第二类边界条件(Neuman) 混