(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 20221093720 5.7 (22)申请日 2022.08.05 (71)申请人 中国地质科 学院地质力学研究所 地址 100000 北京市海淀区民族大 学南路 11号 (72)发明人 王文磊　袁长江　乔东旭　王文君　 (74)专利代理 机构 北京超凡宏宇专利代理事务 所(特殊普通 合伙) 11463 专利代理师 董艳芳 (51)Int.Cl. G06F 30/20(2020.01) G06F 17/18(2006.01) G06Q 50/02(2012.01) G06F 119/02(2020.01) (54)发明名称 分析成矿元素相关性的地理加权回归模型 创建方法及装置 (57)摘要 本申请提供了一种分析成矿元素相关性的 地理加权回归模 型创建方法及装置， 所述创建方 法包括： 将目标研究区域的多个采样点中的任意 一个采样点作为回归分析点， 构建以该回归分析 点为椭圆中心的多个椭圆搜索区域； 根据空间加 权U统计量法确定该回归分析点的目标椭圆搜索 区域； 获取目标椭圆搜索区域中包括的所有采样 点； 根据预先确定的每个采样点到该回归分析点 的目标权重， 确定出多个回归系数； 根据回归系 数创建对应的地理加权回归模型。 采用本申请提 供的技术方案不仅可以提升矿区中对应的预定 成矿元素的信息提取效果， 还可用于定量化分析 和描述矿区中预定成矿元素和预定地质活动敏 感元素的共生关系， 也提高了地理加权回归模型 的适用性。 权利要求书3页 说明书23页 附图16页 CN 115292933 A 2022.11.04 CN 115292933 A 1.一种分析成矿元素相关性的地理加权回归模型创建方法， 其特征在于， 所述创建方 法包括： 获取目标研究区域中的多个采样点对应地区的样点信息； 其中， 所述样点信息包括该 采样点对应地区的空间位置坐标、 预定成矿元素的浓度以及至少一种预定地质活动敏感元 素的浓度； 将所述多个采样点中的任意一个采样点作为 回归分析点， 根据 预先确定的椭圆搜索参 数构建以该回归分析点 为椭圆中心的多个椭圆搜索区域； 根据空间加权U统计量法在该回归分析点的多个椭圆搜索区域中确定该回归分析点的 目标椭圆搜索区域； 获取所述目标椭圆搜索区域中包括的所有采样点； 根据所述目标椭圆搜索区域中的每个采样点对应地区的样点信息以及预先确定的每 个采样点到该回归分析点的目标权重， 确定出所述目标椭圆搜索区域的该回归分析点对应 的多个回归系数； 根据所述回归系数创建以该回归分析点为椭圆中心的所述目标椭圆搜索区域对应的 地理加权回归模型， 以用于 分析以该回归分析点为椭圆中心的所述目标椭圆搜索区域中各 个采样点对应地区的每种预定地质活动敏感元素的浓度分别对对应的目标成矿元素的影 响程度。 2.根据权利要求1所述的创建方法， 其特 征在于， 通过以下步骤确定的椭圆搜索参数： 根据每个采样点对应地 区的空间位置坐标， 确定椭圆搜索区域的长半轴范围和短半轴 范围； 根据所述椭圆搜索区域的长半轴范围和短半轴范围确定所述椭圆搜索区域的等效半 径的值域； 根据所述等效半径的值域以及预先设置的步长， 将所述等效半径的值域按照所述步长 进行拆分， 确定出多个等效半径； 根据预先设置的压缩系数的值域以及预设压缩系数间隔， 将所述压缩系数的值域按照 所述预设压缩系数间隔进行拆分， 确定出多个压缩系数； 根据预先设置的方位角的值域以及预设方位间隔， 将所述方位角的值域按照所述预设 方位间隔进行拆分， 确定出多个方位角； 将所述多个等效半径、 多个压缩系数以及多个方位角确定为椭圆搜索参数。 3.根据权利要求1所述的创建方法， 其特征在于， 所述根据 预先确定的椭圆搜索参数构 建以该回归分析点 为椭圆中心的多个椭圆搜索区域的步骤， 包括： 根据预先确定的椭圆搜索参数， 获取所述椭圆参数中的多个等效半径、 多个压缩系数 以及多个方位角； 在所述多个等效半径中， 将初始等效半径作为当前等效半径， 针对每个压缩系数， 构建 以该回归分析点为椭圆中心在当前等效半径下的该压缩系数的每个方位角对应的椭圆搜 索区域； 确定是否具有下一个等效半径； 若具有， 则将下一个等效半径更新为当前等效半径， 继续构建下一个等效半径作为当 前等效半径的多个椭圆搜索区域；权　利　要　求　书 1/3 页 2 CN 115292933 A 2若不具有， 则得到以所述 回归分析点为椭圆中心的所有等效半径中每个等效半径对应 的多个椭圆搜索区域。 4.根据权利要求1所述的创建方法， 其特征在于， 所述在该回归分析点的多个椭圆搜索 区域中确定该回归分析点的目标椭圆搜索区域的步骤， 包括： 在该回归分析点的多个椭圆搜索区域中， 针对椭圆搜索参数中的每个等效半径， 确定 该等效半径对应的多个椭圆搜索区域中每 个椭圆搜索区域对应的空间加权U统计量； 在该等效半径下的每个椭圆搜索区域对应的空间加权U统计量中， 将空间加权U统计量 的绝对值 最大的椭圆搜索区域确定为该等效半径的最优椭圆搜索区域； 根据所述椭圆搜索参数中的每个等效半径， 将每个等效半径的最优椭圆搜索区域对应 的空间加权U统计量的绝对值进 行排序， 将绝对值最大的空间加权U统计量对应的最优椭圆 搜索区域确定为该回归分析点的目标椭圆搜索区域。 5.根据权利要求4所述的创建方法， 其特征在于， 通过以下步骤确定每个椭圆搜索区域 对应的空间加权U统计量： 针对每个椭圆搜索区域， 获取该椭圆搜索区域中的所有采样点中每个采样点对应地 区 的空间位置坐标和预定成矿元 素的浓度； 根据该椭圆搜索区域对应的椭圆参数， 确定出 该椭圆搜索区域的长半轴和短 半轴； 根据所述该椭圆搜索区域的长半轴和短半轴以及该椭圆搜索区域中的所有采样点中 每个采样点对应地区的空间位置坐标， 确定出每 个采样点到所述回归分析点的目标权 重； 根据每个采样点对应地 区的预定成矿元素的浓度和每个采样点对应的目标权重， 确定 该椭圆搜索区域对应的空间加权U统计量。 6.根据权利要求5所述的创建方法， 其特征在于， 所述根据 所述该椭圆搜索区域的长半 轴和短半轴以及该椭圆搜索区域中的所有采样点中每个采样点对应地区的空间位置坐标， 确定出每 个采样点到所述回归分析点的目标权 重的步骤， 包括： 根据所述该椭圆搜索区域的长半轴和短半轴以及该椭圆搜索区域中的所有采样点中 每个采样点对应地区的空间位置坐标， 确定每 个采样点对应的权 重参数； 针对每个采样点， 确定该采样点对应的权重参数是否在大于零并且小于等于一的范围 内； 若是， 则将该采样点对应的权重参数进行平方得到第一权重参数， 将一与所述第一权 重参数的差再进行平方， 得到该采样点到所述回归分析点的目标权 重； 若否， 则将零确定为该采样点到所述回归分析点的目标权 重。 7.根据权利要求1所述的创建方法， 其特征在于， 所述根据 所述目标椭圆搜索区域中的 每个采样点对应地区的样点信息以及预先确定的每个采样点到该回归分析点的目标权重， 确定出所述目标椭圆搜索区域的该回归分析点对应的多个回归系数的步骤， 包括： 根据所述目标椭圆搜索区域中的每个采样点对应地区的样点信息中的预定成矿元素 的浓度， 将所述每个采样点对应地区的预定成矿元素的浓度组成的列向量确定为地理加权 回归模型的因变量列向量； 根据所述目标椭圆搜索区域中的每个采样点对应地区的样点信息中的至少一种预定 地质活动敏感元素的浓度， 将所述每个采样点对应地区的至少一种预定地质活动敏感元素 的浓度组成的矩阵确定为 地理加权回归 模型的自变量矩阵；权　利　要　求　书 2/3 页 3 CN 115292933 A 3