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

权利要求书

1.一种构建用于矿产勘查的地理加权回归模型的方法，其特征在于，所述方法包括：

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

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

通过以下公式确定距离权重矩阵的估算模型：

其中，σ₁、σ₂代表椭圆筛选区域对应的协方差矩阵的对角线上的方差，σ₁₂为椭圆筛选区域对应的协方差矩阵的协方差；椭圆筛选区域对应的协方差矩阵为：其中，R为椭圆筛选区域的旋转矩阵，S为椭圆筛选区域的缩放矩阵，其中，a为椭圆筛选区域的长半轴，b为椭圆筛选区域的短半轴，为θ转换后的弧度数，θ为旋转角度，转换公式：

确定所述椭圆筛选区域的长半轴参数和短半轴参数，以及所述目标样本点与所述椭圆筛选区域内除所述目标样本点之外的其他样本点之间的距离；

将所述长半轴参数、所述短半轴参数和所述距离输入至所述距离权重矩阵的估算模型中，得到各个所述椭圆筛选区域对应的目标样本点的距离权重矩阵；

基于加权最小二乘法，建立回归系数估算模型；

将各个所述距离权重矩阵，以及所述目标样本点包括的属性特征和对应的成矿作用输入至所述回归系数估算模型中，拟合确定各个所述椭圆筛选区域对应的估计参数，其中，所述属性特征包括控矿要素特征和位置特征，所述控矿要素特征包括构造强度特征、中酸性侵入岩特征和地层特征；

基于各个所述椭圆筛选区域对应的估计参数，确定所述目标样本点对应的多个待定地理加权回归模型；

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

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

2.根据权利要求1所述的方法，其特征在于，通过以下步骤构建多个椭圆筛选区域：

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

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

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

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

3.根据权利要求1所述的方法，其特征在于，所述回归系数估算模型为：

其中，代表回归系数估算矩阵，(u_i，v_i)为第i个样本点的空间几何坐标，X表示包含控矿要素特征的矩阵，Y表示包含成矿作用的矩阵，W(u_i，v_i)表示n×n的距离权重矩阵，其主对角线上的元素为给定位置i周围的样本点j的权重，矩阵中其他位置的元素为0。

4.根据权利要求1所述的方法，其特征在于，所述方法还包括：

获取矿产勘查的实际工区所包括的采样点的属性特征，其中，所述属性特征包括控矿要素特征和位置特征，所述控矿要素特征包括构造强度特征、中酸性侵入岩特征和赋矿地层特征；

基于每个所述采样点的属性特征，确定每个所述采样点对应的地理加权回归模型；

从每个所述采样点的地理加权回归模型中获取估计参数，将获取到的估计参数通过插值运算得到各个所述采样点对应的成矿作用及各向异性特征；

基于所述目标工区所包括的各个采样点对应的成矿作用及各向异性特征，确定所述目标工区的矿产勘查结果。