[发明专利]一种确定地质变量分位值的方法和系统

权利要求书

1.一种确定地质变量分位值的方法，其特征在于，所述方法包括以下步骤：

步骤S10：地质变量确定及数据准备，根据具体任务确定待分析的地质变量，针对确定的待分析的地质变量进行数据统计；

步骤S20：建立地质变量变差模型，根据步骤S10得到的地质变量的统计数据，建立变差模型；

步骤S30：变差函数分析，基于步骤S20中的变差模型，交互调整变差关键参数，进行变差分析，得到地质变量的等值线分布图；

步骤S40：工区网格化及编码，将步骤S30中得到的地质变量的等值线分布图进行网格化，并对网格进行编码，得到规则的网格化单元，网格单元包括被等值线覆盖的网格单元和未被等值线覆盖的网格单元；

步骤S41：确定感兴趣的研究范围，包括全区范围、单一局部范围和多个局部范围三种情况；

步骤S42：针对感兴趣区域，选择网格化方式，包括行列数设置、行距及列距设置两种模式，前者是设置固定的行数及列数来划分网格，后者是设置不同的行距和列距来划分网格；

步骤S43：基于前述步骤S41和步骤S42，进行网格划分，记录网格中心空间位置并按行序对网格进行编码；

步骤S50：对步骤S40中得到的网格化单元进行地质变量的均值和方差的计算；

步骤S60：随机模拟抽取网格化单元的变量值及分位值，并进行排序；

其中，步骤S50中，对于所述被等值线覆盖的网格单元，根据覆盖网格单元的等值线的参数值给网格单元进行赋值，基于网格的样本点统计得到均值和方差；对于所述未被等值线覆盖的网格单元，步骤S50包括以下子步骤：

步骤S51：基于网格划分的结果，结合已知网格间空间位置关系，确定网格点模拟插值路径；

步骤S52：采用克里金插值算法计算网格点均值和方差。

2.根据权利要求1所述的确定地质变量分位值的方法，其特征在于，步骤S10中所述的数据统计包括变量参数值和空间坐标。

3.根据权利要求1所述的确定地质变量分位值的方法，其特征在于，步骤S20中所述的变差模型为球状变差模型、高斯变差模型、线性变差模型或指数变差模型。

4.根据权利要求1所述的确定地质变量分位值的方法，其特征在于，步骤S30中所述的关键参数包括块金值、基台值和/或变程值。

5.根据权利要求4所述的确定地质变量分位值的方法，其特征在于，步骤S30中所述地质变量的等值线分布图是通过设置迭代次数、网格参数和分位值来进行校正的。

6.一种确定地质变量分位值的系统，其特征在于，包括：

地质变量确定单元，根据具体任务确定待分析的地质变量，针对确定的待分析的地质变量进行数据统计；

变差模型单元，根据地质变量确定单元得到的地质变量的统计数据，建立变差模型；

变差函数分析单元，基于变差模型单元得到的变差模型，交互调整变差关键参数，进行变差分析，得到地质变量的等值线分布图；

等值线分布图单元，用于工区网格化及编码，将变差函数分析单元中得到的地质变量的等值线分布图进行网格化，并对网格进行编码，得到规则的网格化单元；

均值和方差的计算单元，对等值线分布图单元中得到的网格化单元进行地质变量的均值和方差的计算，网格单元包括被等值线覆盖的网格单元和未被等值线覆盖的网格单元；以及

变量值及分位值确定单元，用于随机模拟抽取网格化单元的变量值及分位值，并进行排序；

其中，等值线分布图单元包括以下模块：

研究范围确定模块，用于确定感兴趣的研究范围，包括全区范围、单一局部范围和多个局部范围三种情况；

网格化模块，用于针对感兴趣区域，选择网格化方式，包括行列数设置、行距及列距设置两种模式，前者是设置固定的行数及列数来划分网格，后者是设置不同的行距和列距来划分网格；以及

网格单元编码模块，基于研究范围确定模块和网格化模块，进行网格划分，记录网格中心空间位置并按行序对网格进行编码；

均值和方差的计算单元中，对于所述被等值线覆盖的网格单元，根据覆盖网格单元的等值线的参数值给网格单元进行赋值，基于网格的样本点统计得到均值和方差；对于所述未被等值线覆盖的网格单元，均值和方差的计算单元包括以下模块：

网格点模拟插值路径模块，用于基于网格划分的结果，结合已知网格间空间位置关系，确定网格点模拟插值路径；以及

均值和方差计算模块，采用克里金插值算法计算网格点均值和方差。