位场处理与解释技术位场转换及处理.pptx

位场转换及处理位场转换及处理 主要是指空间换算，例如向上向下解析延拓；地形改正；不同分量主要是指空间换算，例如向上向下解析延拓；地形改正；不同分量之间的换算，如磁异常的垂直分量化为水平分晕；不同磁化方向之间的之间的换算，如磁异常的垂直分量化为水平分晕；不同磁化方向之间的换算，如化到磁极；以及各阶导数换算、局部异常和区域异常的划分、换算，如化到磁极；以及各阶导数换算、局部异常和区域异常的划分、滤波等等；滤波等等；位场转换和处理的目的可以归纳为以下几方面：位场转换和处理的目的可以归纳为以下几方面：将复杂异常化为简单异常，以满足某些解释方面的需要。将复杂异常化为简单异常，以满足某些解释方面的需要。将实测异常分解及变换，从而可更方便地利用信息，为解释提供更多将实测异常分解及变换，从而可更方便地利用信息，为解释提供更多的手段，提高解释的效果。的手段，提高解释的效果。突出异常的有用信息，压制干扰，区分异常的性质，及提供产状等。突出异常的有用信息，压制干扰，区分异常的性质，及提供产状等。位场转换和处理的方法很多，总的可分为空间域和频率域位场转换和处理的方法很多，总的可分为空间域和频率域两大类。两大类。3.1 等效源法等效源法 在场论中，场和场源具有唯一的对应关系，但在实际中，观测到的场在场论中，场和场源具有唯一的对应关系，但在实际中，观测到的场只是整个场的一部分，再加上观测误差和随机干扰，使其对观测值的解释只是整个场的一部分，再加上观测误差和随机干扰，使其对观测值的解释不可避免地出现多解性，即是说有多种可能的场源分布与观测场在一定的不可避免地出现多解性，即是说有多种可能的场源分布与观测场在一定的误差范围内对应。误差范围内对应。与观测场对应的场源，如果不是真正的场源，我们称它为等效源。与观测场对应的场源，如果不是真正的场源，我们称它为等效源。3.1.1 基本原理基本原理 在观测场对应的多个场源中，选择一组最简单的场源，例如按一定位置分布或在观测场对应的多个场源中，选择一组最简单的场源，例如按一定位置分布或不固定位置的点荷、线荷、磁偶极子、偶极线等，用最优化方法确定它们的质量或不固定位置的点荷、线荷、磁偶极子、偶极线等，用最优化方法确定它们的质量或磁量，使它们产生在观测面上的场值与实测场值相吻合，利用这组等效源就能很方磁量，使它们产生在观测面上的场值与实测场值相吻合，利用这组等效源就能很方便地作各种位场变换，即是说观测值的各种互换值可以等价地认为是等效源产生的便地作各种位场变换，即是说观测值的各种互换值可以等价地认为是等效源产生的各种变换值。各种变换值。曲化平、向上、向下延拓，异常不同分量之间的换，求变换磁化方向的磁异常（包曲化平、向上、向下延拓，异常不同分量之间的换，求变换磁化方向的磁异常（包括化到磁极），以及垂向一次、二次导数等等，均可以用等效源正演计算来求得。括化到磁极），以及垂向一次、二次导数等等，均可以用等效源正演计算来求得。等效源法的特点：等效源法的特点：1.把各种繁杂的位场转换，变成一个简单的正演计算处理，计算过程简单，便于统一处理。2.不丢掉边部测点，条件好的情况下，可适当外推；3.对于地形起伏较大的观测面，作位场转换的效果仍较好；4.由于等效源产生的场仅在观测面内与真实场源产生的场在一定误差范围内吻合，因而用等效源进行位场的转换，特别是向下延拓时就必然只限于一定范围之内，而不是整个空间，这就是该法的局限；等效源方法可以下两步进行：等效源方法可以下两步进行：1.选择等效源模型，并用最优化方法求取等效源的质量或磁荷量；选择等效源模型，并用最优化方法求取等效源的质量或磁荷量；2.用等效源来计算各种位场转换值。用等效源来计算各种位场转换值。3.1.2 等效源的求取等效源的求取 以二维剖面垂直磁异常转换为例：以二维剖面垂直磁异常转换为例：如图如图1所示，选择等效源模型为磁偶极子，磁化方所示，选择等效源模型为磁偶极子，磁化方向的倾角为向的倾角为A，那么第，那么第i个等效源产生在第个等效源产生在第k个观测个观测点上的垂直磁场为：点上的垂直磁场为：记：记：则：则：对于点磁荷：对于点磁荷：叠加求和：叠加求和：按照等效源的定义有：按照等效源的定义有：目标函数：目标函数：从而归结为求解多元函数的极值问题，可用第二章所述的最优化方法求解。需从而归结为求解多元函数的极值问题，可用第二章所述的最优化方法求解。需要注意的是，当等效源位置固定时，要注意的是，当等效源位置固定时，Zmk为自变量为自变量m1，m2，,mn。的线性函。的线性函数，可用直接解线性方程组的方法求解，但是我们发现，此线性方程组的系数数，可用直接解线性方程组的方法求解，但是我们发现，此线性方程组的系数矩阵矩阵CTC的条件不好易于奇异，所以我们采用以下最优化方法解此方程：的条件不好易于奇异，所以我们采用以下最优化方法解此方程：1、最速下降法、最速下降法迭代格式：迭代格式：最速下降方向：最速下降方向：最佳步长：最佳步长：2、阻尼最小二乘法、阻尼最小二乘法做法类同于第二章所述，设初值为做法类同于第二章所述，设初值为 为增量，求解式为增量，求解式(3.1.9)得得 ，用，用 作为新的初值，继续求解新的作为新的初值，继续求解新的 ，迭代下去，直到满足精度要求为止。从式，迭代下去，直到满足精度要求为止。从式(3.1.2)可知，可知，Zm为为m的线性函数，因而的线性函数，因而C矩阵矩阵(即即Zm对对m的一阶偏导数矩阵的一阶偏导数矩阵)为常数矩阵，所为常数矩阵，所以以CTC在整个迭代过程小是不变的，它的各元素仅与等效源的位置及磁化方向和测在整个迭代过程小是不变的，它的各元素仅与等效源的位置及磁化方向和测点的位置有关，而与变量点的位置有关，而与变量m无关这是一个突出优点，可以减少计算丁作量。无关这是一个突出优点，可以减少计算丁作量。用最速下降法求取的等效源各分量之间相差较小，且有一定的规律，作位用最速下降法求取的等效源各分量之间相差较小，且有一定的规律，作位场转换的效果通常较好，但计算比较费时间；场转换的效果通常较好，但计算比较费时间；用阻尼最小二乘法术取的等效源各分量之间相差较大，作位场转换的效果，用阻尼最小二乘法术取的等效源各分量之间相差较大，作位场转换的效果，特别是向下解析延拓的效果较差，但计算较快。特别是向下解析延拓的效果较差，但计算较快。3.1.3 利用等效源进行各种位场转换利用等效源进行各种位场转换 如果已求得了等效源，将它当作等效的场源，那么对于观测值的各种位场转如果已求得了等效源，将它当作等效的场源，那么对于观测值的各种位场转换值，在一定空间范围和一定误差范围内可以等价地认为与等效源所产生的相换值，在一定空间范围和一定误差范围内可以等价地认为与等效源所产生的相应场值相同，因而利用等效源作一系列的正演计算即可得到各种转换值，下面应场值相同，因而利用等效源作一系列的正演计算即可得到各种转换值，下面举几个例子来说明。举几个例子来说明。不管是解析延拓还是曲化平实质上都是求观测面以上或以下某一位置上的不管是解析延拓还是曲化平实质上都是求观测面以上或以下某一位置上的场值，当利用观测值求出等效源以后，要求空间任意一点的场值，可以利场值，当利用观测值求出等效源以后，要求空间任意一点的场值，可以利用式用式（1）向上、向下解析延拓及曲化平）向上、向下解析延拓及曲化平（2）磁异常不同分量之间的换算磁异常不同分量之间的换算水平分量的系数表达式：水平分量的系数表达式：（3）求转换磁化方向后的场值求转换磁化方向后的场值（4）求观测值的垂向一次、二次导数求观测值的垂向一次、二次导数垂向一次导数垂向一次导数垂向二次导数垂向二次导数3.1.4 等效源法的其他应用等效源法的其他应用（1）复杂情况下的化极效果复杂情况下的化极效果 在小比例尺大范围的化磁极资料处理中，会遇到如下的情况，即在一个大范在小比例尺大范围的化磁极资料处理中，会遇到如下的情况，即在一个大范围内，不同位置的地质体的磁化方向不同。这种情况对于常用的化磁极的方法围内，不同位置的地质体的磁化方向不同。这种情况对于常用的化磁极的方法包括频率域的化极方法包括频率域的化极方法)来说是难以实现的。考虑到等效源法各自等效源之间来说是难以实现的。考虑到等效源法各自等效源之间具有一定的独立性，模拟变磁化方向比较容易，于是用等效源理论模型进行了具有一定的独立性，模拟变磁化方向比较容易，于是用等效源理论模型进行了变磁化方向的化极计算，并取得了满意的结果变磁化方向的化极计算，并取得了满意的结果 图图3.5示出一示出一f变磁化方向化极的效果。图变磁化方向化极的效果。图3.5(a)中为三个不同磁化方向的板中为三个不同磁化方向的板状体状体(磁倾角分别为磁倾角分别为30,45,0)所产生的磁异常。这里排列了两层等效源，如图所产生的磁异常。这里排列了两层等效源，如图中中“”所示，让它们的磁倾角在三个板状体附近接近于真实倾角，通过拟合所示，让它们的磁倾角在三个板状体附近接近于真实倾角，通过拟合山曲线所得的等效源再进行化极计算得到如图山曲线所得的等效源再进行化极计算得到如图3.5 (b所示的虚线。图中实线为所示的虚线。图中实线为这三个板状体垂直磁化理论曲线。两者对比可见，除极值部分有较小差异外，这三个板状体垂直磁化理论曲线。两者对比可见，除极值部分有较小差异外，曲线基本吻合，最大相对误差为曲线基本吻合，最大相对误差为4.2%，均方差为，均方差为3.32nT（2）等效源反演效果等效源反演效果32 直接解拉氏方程位场转换直接解拉氏方程位场转换采用边界条件：采用边界条件：得到这组新边界条件下的解析解：得到这组新边界条件下的解析解：3.2.1 以曲面为边界的拉氏方程解以曲面为边界的拉氏方程解 重磁场的位和各分量在没有场源的空间，是满足拉普拉斯方程的，一切问题总是重磁场的位和各分量在没有场源的空间，是满足拉普拉斯方程的，一切问题总是可以归纳为按具体的边界条件求拉氏方程的解。以往考虑该问题的范围是整个上半可以归纳为按具体的边界条件求拉氏方程的解。以往考虑该问题的范围是整个上半空间，但实际测量大多是在地面上进行的，而地面常是起伏的曲面，而且异常值仅空间，但实际测量大多是在地面上进行的，而地面常是起伏的曲面，而且异常值仅存在于局部空间，因此考虑整个上半空间没有必要，而且往往会引起具体异常边部存在于局部空间，因此考虑整个上半空间没有必要，而且往往会引起具体异常边部的各种计算减低精度，所以改用以下的边界条件更为合适。的各种计算减低精度，所以改用以下的边界条件更为合适。方程（方程（3.2.）在边界条件（）在边界条件（3.2.1）的解：）的解：其中：其中：方程（方程（3.2.8）在边界条件（）在边界条件（3.2.2）的解：）的解：方程（方程（3.2.9）在边界条件（）在边界条件（3.2.3）的解：）的解：（1）当测量在一个平面内进行时当测量在一个平面内进行时式中：式中：于是：于是：采样梯形积分公式采样梯形积分公式（2）测点分布于曲面上测点分布于曲面上极小条件：极小条件：3.2.1 二维情况二维情况分离变量法求解：分离变量法求解：于是：于是：采样梯形积分公式采样梯形积分公式极小条件：极小条件：1、影响计算精度的主要因素、影响计算精度的主要因素3.2.3 有关技术方法有关技术方法（2）测区场值零边界的确定，）测区场值零边界的确定，关键是区域场的消除关键是区域场的消除（3）所得拉氏方程基本级数解的系数的近似程度）所得拉氏方程基本级数解的系数的近似程度。随着离实测线和面的距离增加，。随着离实测线和面的距离增加，其近似程度将降低，离得越远，降低愈多。其近似程度将降低，离得越远，降低愈多。(1)在边界条件在边界条件(3.2.1)和和(3.2.2)中，中，x0，x=Lx，y=0，y=Ly定义了一组直立边界面，而对于某一高度定义了一组直立边界面，而对于某一高度z=z0而言，边界面而言，边界面上的位场值为零的假设就不成立了，因为异常为零的边界，上的位场值为零的假设就不成立了，因为异常为零的边界，随着高度的增加，要向外扩展，应为随着高度的增加，要向外扩展，应为x=z+a1,x=z+a2,y=z+b1,y=z+b2的斜面为位或场的零边界，如图的斜面为位或场的零边界，如图3.10所示。所示。这样一来，使拉氏方程的解十分复杂和困难，即使求得了一这样一来，使拉氏方程的解十分复杂和困难，即使求得了一些数字解，也将会有较大的误差，而且不同场源的异常，其些数字解，也将会有较大的误差，而且不同场源的异常，其零边界随高度向外扩展的规律也不相同，零边界随高度向外扩展的规律也不相同，不能淮确地不能淮确地确定，只能近似地取得。而取直立边界使问题变得十分简单，确定，只能近似地取得。而取直立边界使问题变得十分简单，当取当取Lx、Ly比实际边界稍大时，近似的程度也是可行的。比实际边界稍大时，近似的程度也是可行的。2实际达到的精度实际达到的精度 上述方法系统的误差来源，其中一部分，如计算数学上的近似性，是上述方法系统的误差来源，其中一部分，如计算数学上的近似性，是任何数据处理方法都有的，另一部分则是所述方法系统所特有的。但这任何数据处理方法都有的，另一部分则是所述方法系统所特有的。但这绝不是说该法的计算精度差。它的精度是较高的，其任何一种功能均能绝不是说该法的计算精度差。它的精度是较高的，其任何一种功能均能完全满足数据处理的要求。只要给出足够的计算量，一定可以得到所需完全满足数据处理的要求。只要给出足够的计算量，一定可以得到所需的计算精度。所得的基本解是一个级数，在不超过实测点数目的范围内，的计算精度。所得的基本解是一个级数，在不超过实测点数目的范围内，显然项数愈高，计算精度也愈高。取级数项数超过测点数是不行的，因显然项数愈高，计算精度也愈高。取级数项数超过测点数是不行的，因为系数要实测数据确定。由确定系数的方程组可知，当测点数少于项数为系数要实测数据确定。由确定系数的方程组可知，当测点数少于项数即系数的数目时，会发生方程式少于未知数的情况，其解是不能唯一确即系数的数目时，会发生方程式少于未知数的情况，其解是不能唯一确定的。实际情况是定的。实际情况是:为了满足数据处理的要求，所取级数项数可以比测点为了满足数据处理的要求，所取级数项数可以比测点数小得很多。对一般变化比较平缓的异常数据，级数项数和测点数的比数小得很多。对一般变化比较平缓的异常数据，级数项数和测点数的比有有1=10左右就可保证边部亦有相当高的计算精度。对于梯度大，如强地左右就可保证边部亦有相当高的计算精度。对于梯度大，如强地形的磁测资料等情况，用形的磁测资料等情况，用1:3左右也够了。再提高计算精度的潜力还相当左右也够了。再提高计算精度的潜力还相当大，这是问题的一个方面大，这是问题的一个方面:级数项数取得愈高，精度也愈高。级数项数取得愈高，精度也愈高。但这里所说的精度愈高，是仅只对用以确定级数系数的实测数据所在的邻近但这里所说的精度愈高，是仅只对用以确定级数系数的实测数据所在的邻近区域而言的。对比较远的区域，其精度则反而要降低。这是由于级数项数取区域而言的。对比较远的区域，其精度则反而要降低。这是由于级数项数取得愈高，对实测数据的模拟就愈细，因而函数也愈适应于实测数据区域的近得愈高，对实测数据的模拟就愈细，因而函数也愈适应于实测数据区域的近旁。其远处的适应能力，就反而要有所降低了。这在所述方法系统中是经过旁。其远处的适应能力，就反而要有所降低了。这在所述方法系统中是经过多次试算证明了的。在自然界，这样一种规律是有普遍意义的，对局部的过多次试算证明了的。在自然界，这样一种规律是有普遍意义的，对局部的过分适应，往往要以对于全局适应性的降低为代价的。分适应，往往要以对于全局适应性的降低为代价的。3.3 位场转换数值积分法位场转换数值积分法 从重磁勘探中所介绍的积分插值法可知，我们可以将无限平面的积分，变为有从重磁勘探中所介绍的积分插值法可知，我们可以将无限平面的积分，变为有限个环带的积分，最后变成在各圆环上的场值平均值，乘以某一系数，叠加求和来限个环带的积分，最后变成在各圆环上的场值平均值，乘以某一系数，叠加求和来得到各种延拓值。要求圆环上的场值平均值，只能在圆环上取一定数量的点的场值得到各种延拓值。要求圆环上的场值平均值，只能在圆环上取一定数量的点的场值来求平均，而这些点必须能代表场的特点，所以它不可能只取落在圆环上的测点，来求平均，而这些点必须能代表场的特点，所以它不可能只取落在圆环上的测点，而要大量的取一些非测点的点，这些点的场值只能靠插值的方法获得。从而可以看而要大量的取一些非测点的点，这些点的场值只能靠插值的方法获得。从而可以看出积分插值法的计算方法是反复地使用插值法，由于插值的误差可能给结果带来较出积分插值法的计算方法是反复地使用插值法，由于插值的误差可能给结果带来较大的影响，对此必然会想到是否可以直接采用测点的场值，用方形面积来代替无限大的影响，对此必然会想到是否可以直接采用测点的场值，用方形面积来代替无限平面的积分，因而相应的计算格式也就会发生变化，这就是本节要介绍的内容。平面的积分，因而相应的计算格式也就会发生变化，这就是本节要介绍的内容。除此以外，积分插值法的基本原理是，对磁场沿垂直方向的变化规律用一个四除此以外，积分插值法的基本原理是，对磁场沿垂直方向的变化规律用一个四次多项式来描述，正因为这样，使场源处的奇点性不能被反映出来，那么如果用一次多项式来描述，正因为这样，使场源处的奇点性不能被反映出来，那么如果用一个有理分式来描述，不是比多项式更为合理吗个有理分式来描述，不是比多项式更为合理吗?在场源处，有理分式可以反映出奇在场源处，有理分式可以反映出奇点的特性，并能利用这些奇点为我们提供场源的某些信息。点的特性，并能利用这些奇点为我们提供场源的某些信息。331 平面向上延拓平面向上延拓应用梯形积分公式得：应用梯形积分公式得：332 用有理分式计算向下延拓用有理分式计算向下延拓