第2章电测深曲线正演数学模型及程序设计08.doc

第二章 电测深曲线正演数学模型及程序设计2-1 电测深正演数学模型一电测深视电阻率褶积分表达式对称四极时：（视电阻率褶积积分表达式的导出过程）已知： （1）（为了将上式变换成可以用计算机计算的离散形式，对和的自变量取对数，引入新变量和）设 （2）代入 式（意注意推导过程）令： （3）输入信号输出信号滤波器脉冲响应结论：电测深曲线正演计算可采用空间域数字滤波过程来完成。也可写成褶积表达式： 采用对数变量后，其它常用装置系统的视电阻率褶积运算的表达形式分别为：二极装置： 温纳装置： 式中：偶极装置： 式中是为方位装置（时为赤道偶极）；时为径向装置（时为轴向偶极）二电测深视电阻率离散化数学模型1离散式： （4） 采样间隔 第个采样点上核函数值2将（4）式代入（3）式 令，当 时 令： （5）改写成编程使用的正演计算模型为： （6）第j个计算点序列号 第i个滤波点序号 正演滤波系数如果给出滤波系数，再根据递推公式计算出相应点的核函数值，则可由上述公式计算出相应采样点的视电阻率的值。（滤波系数计算过程将在下一节详细说明）3的计算：取 每个对数节取6个点，间隔固定。，由递推公式计算由以上分析可知，只要给出地电断面层参数的值，由递推公式求得一系列核函数采样值。2-2 滤波系数计算方法及正演滤波系数在电测深正演和反演问题中用到三种滤波系数：1）、由核函数计算视电阻率函数值的正演滤波器；2）、由视电阻率函数确定核函数（视电阻率转换函数）的反演滤波器；3）、由一种装置的视电阻率函数确定另一种装置的视电阻率函数的变换滤波器。（正演和反演滤波系数计算方法相同，目前常用的计算方法有三种：）一频率域中的付氏变换法：由上节分析已经得出了各种装置电测深视电阻率函数和核函数在空间域中的褶积积分形式： （1）DPGhosh首次采用付氏变换计算电测深滤波系数，给出满足上式的三对函数，其中之一为： （2）对（2）式中的函数进行付氏变换，分别求出的频谱以及的频谱，则（1）式的褶积可表示为： （3） 为滤波器滤波函数的频谱；式中： 滤波器振幅，相位频率响应。对（3）式进行反付氏变换：（得到滤波器的脉冲响应） （4）式中： 奈奎斯特频率 对 作付氏变换时的采样间隔为了缩短滤波器的长度，减小截断误差的影响，滤波系数取样尽量与脉冲响应的节点（脉冲响应曲线与横轴的交点）重合。为此在用上（4）式计算滤波系数时，取样点向左移动一段距离，使取样点位于的节点处，其位移值为： （5）计算点上的滤波系数时，（4）式应变为： （6）将（5）式中的值代入（6）式得：式中： 截止频率 频率为时的相位响应值。利用付氏变换，可计算各种装置类型的电测深滤波系数。对于反演滤波系数，仍可选择满足褶积公式的成对函数，采用付氏变换法求取。若设计由一种电测深装置转换到另一种装置的视电阻率滤波系数，则只要分别选取这两种装置相应的视电阻率函数，一个为输入，一个为输出，仍可采用付氏变换法计算其变换滤波系数。二最小二乘法基本原理：调整滤波系数，使滤波器实际输出与理想输出之间的偏差平方和最小。视电阻率线性滤波公式： （7） 设待求滤波系数共个，记为。选择输入函数和输出函数变换对作线性组合，作为（7）式的输入函数和输出函数对。记输入函数，由上式可得滤波器的实际输出为： ， （8）设滤波器的理论输出函数度采样值为，则实际输出和理论输出之差的平方和记为： （9）为使目标函数值达到最小，令对所有的滤波系数的偏导数为零，即：， （10）（10）式可变为： （11） 令：， 则（11）式可变为：，也可列出以下方程组形式： （12）该方程组系数矩阵为正定矩阵，称托布里兹矩阵。因此，计算滤波系数问题最终归结为求解方程组。三数值积分法：（课后自己看，讨论）前面讲到由时，有： 条件： 上式存在能写出的滤波器输出函数的精确表达式。确定滤波系数的无穷积分必须存在。滤波系数的精度：1代数和为1：2二层理论曲线计算的位与滤波得到的值比较，3利用Ghosh给出的一组式对函数计算比较。四各种装置的正演滤波系数1对称四极电测深曲线正演滤波系数 表2-2 温纳2二极装置3偶极装置2-3 电测深理论曲线正演计算实用程序数学模型 （1）计算步骤：1输入层参数等。2递推公式计算 3滤波计算：由式计算：4存盘，显示等。程序见P172，附表一。附录一： 程序设计 P173 , 420430 2-4 线性滤波法用于电测深正演计算的其它应用在讨论线性滤波法正演计算视电阻率曲线原理的基础上，这里介绍用同一滤波系数计算不同装置的视电阻率方法；同时为了理论研究的需要，也简要地讨论利用线性滤波法计算电位、电场强度及电流密度的基本方法。一、用二极装置滤波系数计算各种装置视电阻率值二极装置视电阻率函数式 （1）采用对数变量时： （2）因此，二极装置滤波函数称为零阶贝塞尔函数滤波器。经过数学变换，各种装置视电阻率函数都可以表示为某一输入函数和二极装置滤波函数的褶积，因而可用用已知的二极装置滤波系数计算任意装置的电测深曲线。1、对称四极装置 令：（称为偏心率），则有： （3）由此可以看出，对称四极装置视电阻率值是两个二极装置视电阻率和的线性组合。 （4）2、温纳装置 （5）3、偶极装置 由第一章可知偶极装置视电阻率表达式为（式1-86）： （6）（式中为对称四极梯度装置）将二极装置的表达式两边对求导并且同剩以： （7） （二极装置） （对称四极梯度装置） 可变换为 （8）将（8）式代入（6）式： （9）二、用理论对称四极装置滤波系数计算各种装置视电阻率 其滤波函数为称为一阶贝塞尔函数滤波器。 经过数学变换，可把各种装置视电阻率表示为某一输入函数和理论对称四极装置滤波函数的褶积，从而可用理论对称四极滤波系数计算其它装置的视电阻率值（具体过程略，只列出其表达式）。 1、二极装置 （10）式中：将（10）式写成对数形式：2、偶极装置 （书P35）三、线性滤