一种基于数据正则化的地球化学异常圈定方法及系统技术方案

本发明专利技术涉及一种基于数据正则化的地球化学异常圈定方法及系统。该方法包括：获取研究区的勘查地球化学数据；在地理信息系统GIS平台上，根据所述勘查地球化学数据生成地球化学数据点文件，并获取所述研究区的GIS格式地层文件；根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件；基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；计算所述正则化后的地球化学数据的异常下限；根据所述异常下限圈定地球化学异常。本发明专利技术能够有效去除或过滤由地层岩性造成的假异常，识别弱异常，以及提高选取的异常下限的合理性。以及提高选取的异常下限的合理性。以及提高选取的异常下限的合理性。

【技术实现步骤摘要】
一种基于数据正则化的地球化学异常圈定方法及系统


[0001]本专利技术涉及勘查地球化学领域，特别是涉及一种基于数据正则化的地球化学异常圈定方法及系统。

技术介绍

[0002]勘查地球化学是重要的矿产勘查方法之一，在矿产勘查中占有重要的地位。对勘查地球化学数据的处理结果将直接影响对工作区成矿地质作用的认识以及对勘查靶区的圈定，进而直接影响矿产勘查的结果。
[0003]前人对勘查地球化学数据的处理方法主要是将全部数据采用距离平方反比或者克里格法进行插值后，生成等值线图，进而采用一些通用的方法如均值+2倍的标准差或者分形方法确定异常下限。其不足之处是对工作区内的地质背景差异考虑不够，不能充分反映各岩性或地质体造成的地球化学差异，进而可能导致生成一些假异常或不能识别弱异常；另外，对其异常下限确定的方法，多以假设数据符合正态分布为前提，而实际上，地球化学数据常不呈正态分布，导致选取异常下限的合理性差，进而导致靶区圈定可靠性低。

技术实现思路

[0004]本专利技术的目的是提供一种基于数据正则化的地球化学异常圈定方法及系统，以解决现有的勘查地球化学数据的处理方法不能充分反映各岩性或地质体造成的地球化学差异，导致生成一些假异常或不能识别弱异常，选取的异常下限的合理性低的问题。
[0005]为实现上述目的，本专利技术提供了如下方案：
[0006]一种基于数据正则化的地球化学异常圈定方法，包括：
[0007]获取研究区的勘查地球化学数据；所述勘查地球化学数据包括采样点的坐标以及所述研究区内各个地球化学元素的含量；
[0008]在地理信息系统GIS平台上，根据所述勘查地球化学数据生成地球化学数据点文件，并获取所述研究区的GIS格式地层文件；
[0009]根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件；
[0010]基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；
[0011]计算所述正则化后的地球化学数据的异常下限；
[0012]根据所述异常下限圈定地球化学异常。
[0013]可选的，所述获取研究区的勘查地球化学数据，之后还包括：
[0014]利用每个地球化学元素对应的元素检测限的1/2替换元素含量为空值或者低于元素检测限的值。
[0015]可选的，所述根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件，具体包括：
[0016]将所述地球化学数据点文件以及所述GIS格式地层文件进行空间相交分析，生成带有地层属性的地球化学点文件。
[0017]可选的，所述基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据，具体包括：
[0018]基于所述带有地层属性的地球化学点文件，以不同地层属性为类别，生成各个类别的箱线图，并计算每个所述箱线图内的中位数、中位数绝对偏差、第三四分位数以及四分差；
[0019]利用所述中位数以及所述四分差对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据。
[0020]可选的，所述利用所述中位数以及所述四分差对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据，具体包括：
[0021]利用公式对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；其中，Z
ij
为正则化后的地球化学数据的值，x
ij
为j元素含量值，Median
j
为j元素的中位数，IQR
j
为j元素的四分差；i为地层类别，j为元素类别。
[0022]可选的，所述基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据，之后还包括：
[0023]对所述正则化后的地球化学数据的值进行插值处理，生成元素含量等值线图。
[0024]可选的，所述计算所述正则化后的地球化学数据的异常下限，具体包括：
[0025]利用公式X＝Q3+1.5
×
IQR计算所述正则化后的地球化学数据的异常下限；其中，X为异常下限，Q3为第三四分位数，IQR为四分差。
[0026]一种基于数据正则化的地球化学异常圈定系统，包括：
[0027]勘查地球化学数据获取模块，用于获取研究区的勘查地球化学数据；所述勘查地球化学数据包括采样点的坐标以及所述研究区内各个地球化学元素的含量；
[0028]地球化学数据点文件以及GIS格式地层文件生成模块，用于在地理信息系统GIS平台上，根据所述勘查地球化学数据生成地球化学数据点文件，并获取所述研究区的GIS格式地层文件；
[0029]带有地层属性的地球化学点文件生成模块，用于根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件；
[0030]正则化处理模块，用于基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；
[0031]异常下限确定模块，用于计算所述正则化后的地球化学数据的异常下限；
[0032]地球化学异常圈定模块，用于根据所述异常下限圈定地球化学异常。
[0033]可选的，还包括：
[0034]替换模块，用于利用每个地球化学元素对应的元素检测限的1/2替换元素含量为空值或者低于元素检测限的值。
[0035]可选的，所述正则化处理模块，具体包括：
[0036]箱线图生成单元，用于基于所述带有地层属性的地球化学点文件，以不同地层属性为类别，生成各个类别的箱线图，并计算每个所述箱线图内的中位数、中位数绝对偏差、
第三四分位数以及四分差；
[0037]正则化处理单元，用于利用所述中位数以及所述四分差对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据。
[0038]根据本专利技术提供的具体实施例，本专利技术公开了以下技术效果：本专利技术提供一种基于数据正则化的数据处理的异常圈定方法及系统，通过对带有地层属性的地球化学点文件中不同地层属性内的地球化学数据进行正则化处理，充分考虑了研究区内的地质背景差异，反映了各岩性或地质体的造成的地球化学差异，有效去除或过滤由地层岩性造成的假异常，识别弱异常，进而指导矿产勘查，确定异常下限的方法更为科学合理，因此，本专利技术所提供的圈定方法及系统具有原理简单、方法科学、操作方便、运行效率高等特点。
附图说明
[0039]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于数据正则化的地球化学异常圈定方法，其特征在于，包括：获取研究区的勘查地球化学数据；所述勘查地球化学数据包括采样点的坐标以及所述研究区内各个地球化学元素的含量；在地理信息系统GIS平台上，根据所述勘查地球化学数据生成地球化学数据点文件，并获取所述研究区的GIS格式地层文件；根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件；基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；计算所述正则化后的地球化学数据的异常下限；根据所述异常下限圈定地球化学异常。2.根据权利要求1所述的基于数据正则化的地球化学异常圈定方法，其特征在于，所述获取研究区的勘查地球化学数据，之后还包括：利用每个地球化学元素对应的元素检测限的1/2替换元素含量为空值或者低于元素检测限的值。3.根据权利要求1所述的基于数据正则化的地球化学异常圈定方法，其特征在于，所述根据所述地球化学数据点文件以及所述GIS格式地层文件生成带有地层属性的地球化学点文件，具体包括：将所述地球化学数据点文件以及所述GIS格式地层文件进行空间相交分析，生成带有地层属性的地球化学点文件。4.根据权利要求1所述的基于数据正则化的地球化学异常圈定方法，其特征在于，所述基于所述带有地层属性的地球化学点文件，对不同地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据，具体包括：基于所述带有地层属性的地球化学点文件，以不同地层属性为类别，生成各个类别的箱线图，并计算每个所述箱线图内的中位数、中位数绝对偏差、第三四分位数以及四分差；利用所述中位数以及所述四分差对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据。5.根据权利要求4所述的基于数据正则化的地球化学异常圈定方法，其特征在于，所述利用所述中位数以及所述四分差对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据，具体包括：利用公式对不同的地层属性内的地球化学数据进行正则化处理，确定正则化后的地球化学数据；其中，Z
ij
为正则化后的地球化学数据的值，x
ij
为j元素含量值，Median
j
为j元素的中位数，IQR
j

【专利技术属性】
技术研发人员：孙莉，肖克炎，高阳，
申请(专利权)人：中国地质科学院矿产资源研究所，
类型：发明
国别省市：