[发明专利]一种构建用于矿产勘查的地理加权回归模型的方法及装置

说明书

本发明提供了一种构建用于矿产勘查的地理加权回归模型的方法及装置，通过在目标工区选取目标样本点构建多个椭圆筛选区域，并利用加权最小二乘法，对目标样本点包括的属性特征和对应的成矿作用进行拟合，得到目标样本点对应的多个待定地理加权回归模型，然后从中选取出最小的残差平方和对应的待定地理加权回归模型为最优地理加权回归模型，这样一来，以基于水系沉积物地球化学数据的矿产勘查为例，可以表征与成矿相关的各种地质要素的各向异性特征的影响，可以客观地反映复杂的地质过程及其产物的空间分布，有利于提高对矿产勘查的分析精度。

技术领域

本发明涉及矿产勘查技术领域中具有各向异性特征的事件、过程，具体而言，涉及一种构建用于矿产勘查的地理加权回归模型的方法及装置。

背景技术

在成矿过程和矿产资源勘察过程中，常采用地理加权回归分析模型(Geographically Weighted Regression，GWR)对给定区域内的少量样本进行局部加权最小二乘回归，估计地理变量之间的局部关系，揭示局部区域详细的变化。

以成矿地质事件为例，由于成矿过程具有复杂性、各向异性、不可逆性、自相似性等非线性特征，传统的地理加权回归模型认为距回归点相等距离的采样点对于回归点有相同影响力，无法客观地反映复杂的地质过程及其产物的空间分布，并且无法考虑与成矿相关的各种地质要素的各向异性(anisotropy)特征，使得对控矿作用强度及空间分布的分析精度较低，影响矿产勘查精度。

发明内容

有鉴于此，本发明的目的在于提供一种构建用于矿产勘查的地理加权回归模型的方法及装置，可用于变量间各向异性特征属性的描述；并将该地理加权回归模型在基于水系沉积物地球化学矿产勘查中进行应用示范，可以解决与成矿相关的各种地质要素的各向异性特征所带来的影响，有利于提高对矿产勘查的分析精度。

第一方面，本申请提供了一种构建用于矿产勘查的地理加权回归模型的方法，所述方法包括：

获取目标工区所包括的样本点，以及每个样本点包括的属性特征和对应的成矿作用；

以任意一个样本点为目标样本点，以所述目标样本点的所在位置为几何中心，构建多个椭圆筛选区域；

基于所述椭圆筛选区域的长半轴参数和短半轴参数，以及所述目标样本点与所述椭圆筛选区域内除所述目标样本点之外的其他样本点之间的距离，确定各个所述椭圆筛选区域对应的目标样本点的距离权重矩阵；

基于各个所述距离权重矩阵，以及所述目标样本点包括的属性特征和对应的成矿作用，拟合确定各个所述椭圆筛选区域对应的估计参数，得到所述目标样本点对应的多个待定地理加权回归模型；

确定所述目标样本点对应的各个所述待定地理加权回归模型的残差平方和；

从确定出的所有残差平方和中得到最小的残差平方和对应的待定地理加权回归模型为最优地理加权回归模型。

优选地，通过以下步骤构建多个椭圆筛选区域：

获取预设增加步长和预设增加角度值，其中，所述预设增加角度值为所述椭圆筛选区域绕着所述目标样本点旋转的旋转角度变化量；

在所述椭圆筛选区域的长半轴参数和短半轴参数上分别循环增加所述预设增加步长以及在所述椭圆筛选区域的旋转角度参数上循环增加所述预设增加角度值，得到多个待定长半轴参数、待定短半轴参数和待定旋转角度参数；

直至所述待定长半轴参数达到预设长轴参数阈值，所述待定短半轴参数达到预设短轴参数阈值，且所述待定旋转角度参数达到预设旋转角度参数阈值，停止循环；

基于获得的多个待定长半轴参数和待定短半轴参数，以及多个待定旋转角度参数，确定所述目标样本点对应的多个椭圆筛选区域。