地球物m理、地震正演模拟方法.pptVIP

建立二次场方程的积分方程 二次电场可通过将散射电流源Js乘以适当的并矢格林函数G(r，r/)，并对异常体所占的体积做积分而得， 如假设异常体内的电导率为常数?2，则可得到实测电场表达式为 上式是一非齐次的第二类矢量弗雷德霍姆(Fredholm)积分方程式 并矢格林函数 g是对应全空间的标量格林函数 ●P/(x,y,-z) ●P(x,y,z) Q(x,y,z) r r 数值积分，求解 将异常体剖分成N个线性尺寸为?的立方体单元．并假设在每个单元内电场是常数，则积分可用求和式来逼近，再经一系列推导得到分块矩阵方程式 矩阵[M]的每个元素本身就是一个3X3阶矩阵： Γmn是相对有限单元体积电流的并矢格林函数 解方程，求出异常体内每个单元中心处的电场值之后，再加上一次场值可求得异常体外任一点处的电场 。对总电场方程应用法拉第定律，可计算任意点的磁场H(r) 物理模拟方法 基本原理：相似原理 方法：按比例复制地质模型 (通常比例尺为1:100和1:100万之间) ； 模型的物性参数一 般也应按一定的比例改变； 观测装置要微型化。 缩小模型响应能代表野外实际模型响应 模拟准则：使实际模型场与模拟场具有相同的幅值和规律 一、电磁场物理模拟的基本原理 1．频率域电磁场的模拟准则 野外条件下地电体参数方程 模型条件下地电体参数方程 室内模型模拟系统尺寸与野外比例关系为1/l——几何结构比例尺 将比例关系代入野外方程得模型参数描述的野外方程 按照模拟准则，要使模拟结果与实际一致，野外和模型磁场满足的波动方程应完全相同，对比可以写出 模拟准则——参数比例尺 上式两边为响应参数k2r2，亦称综合参数,右端=1，表明：具有相同综合参数的每个系统必定产生相同的电磁响应 ，与?、?、?及l的具体数值无关。因此，模拟准则简化为 2．时间域电磁场的模拟准则 3．模拟模型的种类-4种 导电性溶液 或固体作为均匀介质 导电性更好的覆盖 目标体 目标体 目标体 目标体 4、物理模报实验装置 5．模型材料 超声波模拟 几何相似性 L为几何尺度，?为波长 物理相似性 速度比和密度比相似 纵、横波速度比的相似性 谢谢大家！ 模型应能够反映主要地质构造和岩石、矿物特征，具有代表性或普遍性（共性）、针对性（目的性）、特殊性（特殊问题） 模型不宜太复杂，否则无法建立相应的数学模型；或者计算结果太复杂，难以分析、辨认地质特征与地球物理场特征之间的联系。 地球模型建立的要求： 常用数值计算方法 有限差分法 有限元法 积分方程法 快速离散傅里叶变换法 拟谱法（伪、虚）谱法 射线追踪法 计算速度快 边界刻化好 易于处理有源问题 F域计算 易刻化运动学特性 计算方法选择很重要：不同方法适应性不同，计算效率、机时、周期、费用等各不相同，要综合考虑各种因素合理选用计算方法。有针对性地灵活运用计算方法 微分方程法 正演方法 河南理工大学资源与环境学院 云美厚 地球物理学的问题 正演问题 反演问题 按事物一般原理（或模型）及相关的条件（初始条件、边界条件）来预测事物的结果（可由观测可得 据地球物理场的实际观测值（有时也用理论计算值）定量或定性解释推断地球内部结构（地质体形态和岩层物性）。 基础 目的 应用地球物理学的基本方程式——阻尼标量波动方程 式中，u表示地球物理场的一种，如声场．电磁场的某一分量等；f(x,t)为源函数；x为空间的一个点；t为时间；系数h和g对不同场有不同的物理意义。 位场：在场源外区域满足拉普拉斯方程的物理场称为，如重力场、磁场和稳定电流场 波场：在场源外区域满足波动方程或扩散方程的物理场，如电磁场、弹性波场 选择计算方法，编制计算程序，进行数值计算。 数学模拟方法求解地球物理正演问题的一般步骤： 第一步，地质建模： 据研究对象和问题建立地球模型或地质结构模型； 第二步，数学建模： 据使用的物理手段和地球模型建立相应的数学模型； 第三步，模拟计算： 求解正演问题 地球物理模拟 物理模拟 相似原理 投资大，选材难，结果真实， 数学 模拟法 解析方法 最简捷方便 ，仅适用少数简单模型 数值 模拟法 正演主要工具 效率高，机时少，周期短，费用低。 概念：将描述各种地球物理场的方程或表达式及初、边值条件通过数值方法求出它们的数值解。 模型应能够反映主要地质构造和岩石、矿物特征，具有代表性或普遍性（共性）、针对性（目的性）、特殊性（特殊问题） 模型不宜太复杂，否则无法建立相应的数学模型；或者计算结果太复杂，难以分析、辨认地质特征与地球物理场特征之间的联系。 地球模型建立的要求： 常用数值计算方法 有限差分法 有限元法 积分方程