[发明专利]用于海山表面空间数据的变差函数模拟的地质统计学方法

说明书

本发明公开了一种用于海山表面空间数据的变差函数模拟的地质统计学方法，定义3D‑Surface区域化变量；利用基于“距离‑坡度”或“距离‑相对水深”的变差函数计算：将变程拟合为坡度的连续函数；进行KRIGING估值。

技术领域

本发明涉及地质统计学，变差函数，海洋地质、矿产资源定量评价领域，是一项关于空间数据的“变差函数”及地质统计学的新计算方法，主要用于海山表面空间数据的变差函数模拟及地质统计学插值。

背景技术

目前，将海山表面产出的矿产资源数据为模型定义三维表面(3D-Surface)区域化变量：二维地理坐标上连续取值，三维空间上只有唯一值.某地理坐标点只对应第三维空间的唯一值，即水深数据不像3D空间数据那样纵向上连续取值，某地理点上只有唯一值；另外三维表面空间上的参数分布不受方向控制，只受第三维空间的相对变化控制；传统变差函数计算方法不再适合这类空间数据；需要新的变差函数计算方法与之相适应。

传统“距离-方向”变差函数，不适合以海山表面矿产资源为数据模型的3D-Surface空间数据，因为海山矿产资源数据分布不受方向控制，而受水深控制。

传统地质统计学(KRIGING)估值方法，也不再适合3D-Surface空间数据的空间差值。

定义3D-Surface区域化变量；三维表面区域化变量有其特殊性，导致传统方法不能有效，对诸如海山表面空间数据等3D-Surface区域化变量，进行变差函数计算，也不能进行KRIGING估值。

本案定义基于“距离-坡度”或“距离-相对水深”的变差函数概念和计算方法；这个模型方法比传统基于“距离-方向”的变差函数，更加适合3D-Surface区域化变量的变差函数计算。

KRIGING估值，传统方法只是将水深当成纵向坐标，本案将点对相对水深或坡度看成点对变化性的自变量。

发明内容

本发明的目的是针对现有技术的不足，提出一种用于海山表面空间数据的变差函数模拟的地质统计学方法。

为实现本发明之目的，本发明采用以下技术方案予以实现：

一种用于海山表面空间数据的变差函数模拟的地质统计学方法，利用如下公式进行计算：

假设有n个数据信息点Z(x_i，y_i，dp_i)，i＝1，...，n

当头点Z(x_h，y_h，dp_h)确定，按照一定的坡度和距离搜寻尾点，按如下步骤：

空间距离：ds′_ht＝sqart[(x_h-x_t)²+(y_h-y_t)²+(dp_h-dp_t)²]

水平距离：ds_ht＝sqart[(x_h-x_t)²+(y_h-y_t)²] (3)

由于坡度计算要用到水平距离，因此用水平距离代替空间距离作为距离的度量参数，下文中所涉及的距离均指水平距离。

坡度：

n个实验点，共有个坡度数据，最大坡度值g_max，最小坡度值g_min