[工学]4位场处理与解释技术位场转换及处理.ppt

采样梯形积分公式 （2） 测点分布于曲面上 极小条件： 3.2.1 二维情况 分离变量法求解： 于是： 采样梯形积分公式 极小条件： 1、影响计算精度的主要因素 3.2.3 有关技术方法 （2）测区场值零边界的确定，关键是区域场的消除 （3）所得拉氏方程基本级数解的系数的近似程度。随着离实测线和面的距离增加，其近似程度将降低，离得越远，降低愈多。 (1)在边界条件(3.2.1)和(3.2.2)中，x＝0，x=Lx，y=0，y=Ly定义了一组直立边界面，而对于某一高度z=z0而言，边界面上的位场值为零的假设就不成立了，因为异常为零的边界，随着高度的增加，要向外扩展，应为x=αz+a1 ,x=βz+a2, y= αz+b1 ,y=βz+b2的斜面为位或场的零边界，如图3.10所示。这样一来，使拉氏方程的解十分复杂和困难，即使求得了一些数字解，也将会有较大的误差，而且不同场源的异常，其零边界随高度向外扩展的规律也不相同，α，β不能淮确地确定，只能近似地取得。而取直立边界使问题变得十分简单，当取Lx、Ly比实际边界稍大时，近似的程度也是可行的。 实际达到的精度 上述方法系统的误差来源，其中一部分，如计算数学上的近似性，是任何数据处理方法都有的，另一部分则是所述方法系统所特有的。但这绝不是说该法的计算精度差。它的精度是较高的，其任何一种功能均能完全满足数据处理的要求。只要给出足够的计算量，一定可以得到所需的计算精度。所得的基本解是一个级数，在不超过实测点数目的范围内，显然项数愈高，计算精度也愈高。取级数项数超过测点数是不行的，因为系数要实测数据确定。由确定系数的方程组可知，当测点数少于项数即系数的数目时，会发生方程式少于未知数的情况，其解是不能唯一确定的。实际情况是:为了满足数据处理的要求，所取级数项数可以比测点数小得很多。对一般变化比较平缓的异常数据，级数项数和测点数的比有1=10左右就可保证边部亦有相当高的计算精度。对于梯度大，如强地形的磁测资料等情况，用1:3左右也够了。再提高计算精度的潜力还相当大，这是问题的一个方面:级数项数取得愈高，精度也愈高。 但这里所说的精度愈高，是仅只对用以确定级数系数的实测数据所在的邻近区域而言的。对比较远的区域，其精度则反而要降低。这是由于级数项数取得愈高，对实测数据的模拟就愈细，因而函数也愈适应于实测数据区域的近旁。其远处的适应能力，就反而要有所降低了。这在所述方法系统中是经过多次试算证明了的。在自然界，这样一种规律是有普遍意义的，对局部的过分适应，往往要以对于全局适应性的降低为代价的。 3.3 位场转换数值积分法 从重磁勘探中所介绍的积分插值法可知，我们可以将无限平面的积分，变为有限个环带的积分，最后变成在各圆环上的场值平均值，乘以某一系数，叠加求和来得到各种延拓值。要求圆环上的场值平均值，只能在圆环上取一定数量的点的场值来求平均，而这些点必须能代表场的特点，所以它不可能只取落在圆环上的测点，而要大量的取一些非测点的点，这些点的场值只能靠插值的方法获得。从而可以看出积分插值法的计算方法是反复地使用插值法，由于插值的误差可能给结果带来较大的影响，对此必然会想到是否可以直接采用测点的场值，用方形面积来代替无限平面的积分，因而相应的计算格式也就会发生变化，这就是本节要介绍的内容。 除此以外，积分插值法的基本原理是，对磁场沿垂直方向的变化规律用一个四次多项式来描述，正因为这样，使场源处的奇点性不能被反映出来，那么如果用一个有理分式来描述，不是比多项式更为合理吗?在场源处，有理分式可以反映出奇点的特性，并能利用这些奇点为我们提供场源的某些信息。 3．3．1 平面向上延拓 应用梯形积分公式得： 3．3．2 用有理分式计算向下延拓 位场数据处理与解释 位场转换与处理 中国石油大学（华东） 唐杰 位场转换及处理 主要是指空间换算，例如向上．向下解析延拓；地形改正；不同分量之间的换算，如磁异常的垂直分量化为水平分晕；不同磁化方向之间的换算，如化到磁极；以及各阶导数换算、局部异常和区域异常的划分、滤波等等； 位场转换和处理的目的可以归纳为以下几方面： 将复杂异常化为简单异常，以满足某些解释方面的需要。 将实测异常分解及变换，从而可更方便地利用信息，为解释提供更多的手段，提高解释的效果。 突出异常的有用信息，压制干扰，区分异常的性质，及提供产状等。 位场转换和处理的方法很多，总的可分为空间域和频率域两大类。 3.1 等效源法 在场论中，场和场源具有唯一的对应关系，但在实际中，观测到的场只是整个场的一部分，再加上观测误差和随机干扰，使其对观测值的解释不可避免地出现多解性，即是说有多种可能的场源