综合信息重磁构造解译研究
发布网友
发布时间：2022-04-29 15:53
共1个回答
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
热心网友
时间：2023-10-17 03:28
地球物理信息的综合解译，主要是对重、磁信息进行解译。从统计观点分析航磁图和重力图分别是不同性质、不同时代、不同等级、不同深度磁性体、密度体集合的叠加场，数据处理的目的是对它们的分解，研究不同磁性体、密度体三维空间关系。因为它们是对天然磁场和重力场测量的结果，重、磁信息覆盖全国，有普遍意义。地质图和矿产图均是二维的平面图，只有通过以地质体、矿产资源体为单元的综合信息解译，方能达到三维化认识地质体和矿产资源体的目的。我们对重磁的解译有以下几点地质假设：
1)地质体是无限延深的板状体的集合。从宏观投影分析，地质体的产状近直立。
2)地质体的磁性是均匀的，以感磁场为主，而剩余磁性是次要的。服从地球感磁体的基本特征。
3)地质体集合的重、磁场特征，有很大相关性。
一、重磁环型与线型构造解译
下面以航磁解译工作为例介绍环型与线型构造方法的研究过程(图6-2)：
大兴安岭成矿带铜多金属矿勘查选区
0°，45°，90°，135°4个方向水平一阶导数，共计12个导数，4个有效导数，8个斜交导数。提取4个有效导数，绘制于不同上延高度的平面图上，实现不同上延高度的环型和线型的构造解译
↓不同延拓高度环线构造的上下关联解译
图6-2 环型与线型构造方法的研究流程图
上延高度最大的平面构造最简单，而原平面的构造深度最为复杂。故应先在最大延拓高度的平面图上，进行构造解译，提取出最大深度的环型、线型构造，逐步筛选出不同高度的构造，应用这种方法可以区别不同上延高度的环型和线型构造，形成综合解译的航磁环、线构造格架图(图6-3)。
图6-3 1：20万航磁线型构造解译图
重力和航磁构造解译方法相似，重力场本身是垂直向下，故不用化极，其他工作步骤基本一致。(图6-4)
对航磁异常的研究，首先对比研究化极前后的磁异常特征，是正磁异常、负磁异常、正负跳跃等。其次是否为叠加异常，研究它们的正剩余磁异常、负剩余异常。最后对航磁原平面平剖面图同航磁化极平剖面图相对比，判别是感磁异常还是剩余磁异常特征。
对重力异常的研究，首先注意研究重力异常特征，其次研究是否是叠加重力异常，是正重力剩余异常还是负重力剩余异常。最后研究重、磁异常是否同源，是否吻合。
图6-4 1：20万重力环型构造解译图
二、重磁剩余异常解译
这里的重磁剩余异常是区别于传统意义上的剩余异常，它是利用克里格插值方法对重磁场进行区域场和局部场进行分离，用来突出局部异常，从而更有效的分析重磁场与地质、矿产等成矿因素之间的关系。
克里格方法是矿物学家D.G.Krige(1982)在研究南非金矿分布时提出的一种依据已知点金数据进行空间插值估计金矿分布的一种方法。克里格方法是建立在地质统计学基础之上的一种定量方法，故其数学基础是牢固的。现它已被广泛应用于估计矿产的分布、地下水资源分布等各个方面。该方法依据样本点的空间位置、样本属性及其空间相关性，应用统计学原理来估计未知样本点属性。
克里格法也称空间局部估计或空间局部插值，它建立在变异函数理论及结构分析的基础上，在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克里格法实质上是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏最优估计。从数学角度讲就是对空间分布的数据求线性最优无偏内插估计量的一种方法。更具体的讲，它是根据待估样点或待估地段有限区域内若干已测定的样点数据，在认真考虑了样点的形状、大小和空间相互位置关系，及其与待估样点空间相互位置关系，以及由变异函数提供的结构信息之后，对该待估样点值进行的一种线性无偏最优估计。
设Z(x)是场变量，x是一维、二维或三维空间中的某一个位置，变量Z(x)在x处的值表示为
Z(x)=m(x)+ε(x)+ε′(x) (6.1)
式中：m(x)为描述恒定均值或趋势性有关的结构性成分；ε(x)为与空间变化有关的随机变化项，即区域性变量；ε′(x)为空间上具有零均值及方差σ2与位置无关的高斯噪声项。
Zi(i=1，2，3，…，n)为一组离散的样本数据，现在要对点x0处的变量值Z0进行插值。第一步是确定m(x)，最简单的情况是m(x)等于采样的平均值，距离矢量为h的两点之间的数学期望差为零：
E[z(x)-z(x-h)]=0 (6.2)
假定两点之间的方差只与距离h有关，记为
E[z(x)-z(x+h)]2=E[ε(x)-ε(x+h)]2=2γ(x，h) (6.3)
式中：γ(h)为半方差函数。γ(h)的估算公式如下：
大兴安岭成矿带铜多金属矿勘查选区
变异函数描述随机变量增量间的线性相关程度，是通过平稳增量描述较大尺度的空间范围内随机变量的空间相关结构和随机性的有效工具。Z0的插值表达式为
大兴安岭成矿带铜多金属矿勘查选区
克里格法的原则就是使得权重λi的选择达到Z0为无偏估计，且估计的方差小于其他线性组合产生的方差。无偏估计要求E(Z-Z0)=0，由此得到无偏条件式为
=1，估计方差的计算公式为
大兴安岭成矿带铜多金属矿勘查选区
令
，其中F为n个权系数λi和一个拉格朗日乘子μ的n+1元函数。令F对λi(i=1，2，…，n)和μ的偏导数为零，便得到克里格方程组：
大兴安岭成矿带铜多金属矿勘查选区
式中：j=1，2，…，n-1，n；
=1。由上式计算出权系数λi，便可利用式6.5去估计Z0的值。
三、重磁剩余异常研究
克里格插值方法的优势在于它是一种线性无偏最优插值方法，并且引进了区域化变量的概念，从而使其插值能够结合空间位置关系来进行，尤其是在重磁区域场和局部场的分离，突出局部场方面，比其他插值方法效果更好。
