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

第二章 电测深曲线正演数学模型及程序设计 §2-1 电测深正演数学模型 一﹑电测深视电阻率褶积分表达式 对称四极时：（视电阻率褶积积分表达式的导出过程） 已知： （1） （为了将上式变换成可以用计算机计算的离散形式，对和的自变量取对数，引入新变量和） 设 （2） 代入① 式（意注意推导过程） 令： （3） ——输入信号 ——输出信号 ——滤波器脉冲响应 结论：电测深曲线正演计算可采用空间域数字滤波过程来完成。 也可写成褶积表达式： 采用对数变量后，其它常用装置系统的视电阻率褶积运算的表达形式分别为： 二极装置： 温纳装置： 式中： 偶极装置： 式中是为方位装置（时为赤道偶极）；时为径向装置（时为轴向偶极） 二﹑电测深视电阻率离散化数学模型 1﹑离散式： （4） ——采样间隔 ——第个采样点上核函数值 2﹑将（4）式代入（3）式 令，，当 时 令： （5） 改写成编程使用的正演计算模型为： （6） 第j个计算点序列号 第i个滤波点序号 正演滤波系数 如果给出滤波系数，再根据递推公式计算出相应点的核函数值，则可由上述公式计算出相应采样点的视电阻率的值。 （滤波系数计算过程将在下一节详细说明） 3﹑的计算： 取 每个对数节取6个点，间隔固定。 ， 由递推公式计算 由以上分析可知，只要给出地电断面层参数的值，由递推公式求得一系列核函数采样值。 §2-2 滤波系数计算方法及正演滤波系数 在电测深正演和反演问题中用到三种滤波系数： 1）、由核函数计算视电阻率函数值的正演滤波器； 2）、由视电阻率函数确定核函数（视电阻率转换函数）的反演滤波器； 3）、由一种装置的视电阻率函数确定另一种装置的视电阻率函数的变换滤波器。 （正演和反演滤波系数计算方法相同，目前常用的计算方法有三种：） 一﹑频率域中的付氏变换法： 由上节分析已经得出了各种装置电测深视电阻率函数和核函数在空间域中的褶积积分形式： （1） D．P．Ghosh首次采用付氏变换计算电测深滤波系数，给出满足上式的三对函数，其中之一为： （2） 对（2）式中的函数进行付氏变换，分别求出的频谱以及的频谱，则（1）式的褶积可表示为： ∴ （3） 为滤波器滤波函数的频谱； 式中： 滤波器振幅，相位频率响应。 对（3）式进行反付氏变换：（得到滤波器的脉冲响应） （4） 式中： 奈奎斯特频率 对 作付氏变换时的采样间隔 为了缩短滤波器的长度，减小截断误差的影响，滤波系数取样尽量与脉冲响应的节点（脉冲响应曲线与横轴的交点）重合。为此在用上（4）式计算滤波系数时，取样点向左移动一段距离，使取样点位于的节点处，其位移值为： （5） 计算点上的滤波系数时，（4）式应变为： （6） 将（5）式中的值代入（6）式得： 式中：—— 截止频率 ——频率为时的相位响应值。 利用付氏变换，可计算各种装置类型的电测深滤波系数。对于反演滤波系数，仍可选择满足褶积公式的成对函数，采用付氏变换法求取。 若设计由一种电测深装置转换到另一种装置的视电阻率滤波系数，则只要分别选取这两种装置相应的视电阻率函数，一个为输入，一个为输出，仍可采用付氏变换法计算其变换滤波系数。 二﹑最小二乘法 基本原理：调整滤波系数，使滤波器实际输出与理想输出之间的偏差平方和最小。 视电阻率线性滤波公式： （7） 设待求滤波系数共个，记为。 选择输入函数和输出函数变换对作线性组合，作为（7）式的输入函数和输出函数对。 记输入函数，由上式可得滤波器的实际输出为： ， （8） 设滤波器的理论输出函数度采样值为，则实际输出和理论输出之差的平方和记为： （9） 为使目标函数值达到最小，令对所有的滤波系数的偏导数为零，即： ， （10） （10）式可变为： （11） 令：， 则（11）式可变为： ， 也可列出以下方程组形式： （12） 该方程组