第4章电测深曲线数字解释原理08.docVIP

第四章 稳定电流场数值模拟方法简介 ① 一维地电断面：水平层状均介质 正演法：拟合核函数 直接拟合视电阻率 DZ法 其他快速反演算法 ② 二、三维地电断面：正演：（70年代 Coggon引入） 有限单元法 有限差分（差商代替微商） 积分方程（三维，面电荷密度） 边界单元 ③ 二维反演 高密度电阻率法（工程，环境） 电阻率CT成像： 引起施工方法的变革： 井←→井，井←→地面 概述：电法勘探以测定人工和天然地下电场或电磁场为基础。 首先应解决正演问题 地质条件电性资料—→电磁场分布 方法：1、经典方法解析求解：—→参数简单规则的地电条件 2、数值模拟（计算机）：有限单元法、有限差分法、积分方程及边界元法。 几种方法简介： 有限差分法： 从电场所满足的偏微分方程和边界条件出发，将微分方程转变为差分方程。 实际上，用差商代替微商： 步骤：首先将求解区域剖分成许多小正方形或小长方形的网格（网格剖分），以网格节点上的电场参数值来表征电场的空间分布；然后，用网格节点上电位函数的差商来近似代替该点的偏导数（微商），由此得到一个关于网格节点电位值的高阶线性方程组，即差分方程；解此方程，就可算出网格节点上的场参数值。 有限差分法基本原理 为使问题简单，以一元函数为例说明差分的概念。设自变量有一增量，则相应的函数的增量为： （1） 称为函数的差分（一阶差分）。其中称为函数的微分。表示的高阶无穷小。 通常，自变量的增量又称为自变量的微分，记作，于是 （2） 则 （3） 说明函数的微分与自变量的微分之商等于该函数的导数，因此，导数又称为微商。 当很小时，差分与微分之间的差异将很小，因此，一阶差分除以增量的商，即一阶差商： （4） 对于一般情况，将自变量离散为、、，则称为的一阶向前差商公式为： 的一阶向前差商公式为： （5） 的一阶向后差商公式为： （6） 的一阶中心差商公式为： （7） 由上述三个公式还可得到二阶中心差商公式： （8） 式（5）~式（8） 有限单元法： 有限单元法是二十世纪五十年代中期发展起来的一种数值计算方法。最初主要用于结构和应力分析。七十年代初，J.H.Coggon首先将其用于电法勘探。他从电磁场总能量最小原理出发，建立了用有限单元法进行电（磁）场模拟的算法，并实际计算了二维地电条件下的直流点电源场和线电源场中的电磁法异常。后来，L.R.Rijo进一步完善了这种数值计算方法，使之成为计算二维地电条件下（点源）电阻率法和激发极化法异常的有效方法。我国从七十年代中期开始在电法勘探中做有限元数值模拟，主要是研究轴对称空间二维问题。后来，又引进和发展了Coggon—Rijo的二维地电条件下点源场的有限元算法。 基本原理： 有限元法是根据电场分布服从最小能量原理，将给定边界条件下求解电位的微分方程等价地转换成求解泛函的极值问题，然后将连续的求解域离散化，即将求解区域剖分为一些在节点处相互连接的网格单元，进而在各单元上近似地将变分问题离散化，得到由空间各点未知电位值组成的高阶线性方程组，最后用计算机求解该方程组，可确定各点的电位分布。 有限元法实质上是一种求解场的变分问题的数值方法。它不是直接去求解场的方程，而是去求解一个能量取最极值的问题（变分问题）。 稳定电流场的Laplace方程或Poisson方程的数值解，即三维椭圆型微分方程的边值问题 （2-8） 上述微分方程边值问题等价于下面泛函的极小值问题： （2-9） 泛函问题 1、泛函： 例：函数： 自变量，因变量。即自变量的每一个值都对应的某个值，称是的函数。 若J又是的函数，每一个对应一个J值，则称J是的泛函： 例：A B点之间曲线方程为 则曲线弧长为： 显然J的值取决于曲线具体形式。 可见：泛函就是函数的函数。 2、泛函的变分 （类似于函数的微分） ① 一般定义： 自变量的增量变化成为自变量的变分。 由此引起泛函的增量为： 若泛函的增量可表示为：（泰勒级数展开） ……………① 其中：关于是奇次线性的，当为无穷小时，为高阶无穷小。 式中泛函增量的线性部分即为泛函增量的主部，把泛函增量的主部定义为泛函的变分，记作： 即： 注：与函数增量的区别： Lagrange泛函变分的定义（拉格郎日） 定义：泛函的变分是对的导数在时的值，即： 变分运算基本方法： ⅰ、变分运算与微分求导运算顺序可交换： 即：， ⅱ、变分运算与积分运算顺序可交换：