一种碳酸岩型稀土矿富集区判定方法

本发明涉及矿产地质勘查，特别是涉及一种碳酸岩型稀土矿富集区判定方法。

背景技术：

1、现有碳酸岩型稀土矿富集区判定手段多依靠单一稀土总量或轻重稀土比值开展异常圈定，仅能单一表征稀土元素富集程度，无法同步耦合氟、磷等关键成矿挥发分元素，同时忽略区域镧镱比值的整体分异波动干扰，容易将无矿弱蚀变碳酸岩、低分异围岩误判为矿化异常区；传统评价模型缺少标准化阈值标定流程，人工划定异常边界主观性强，插值生成的空间栅格存在大量孤立噪点，且未结合稀土总量背景值构建定量耦合评价系数，无法实现富集程度梯度分级；此外常规均匀网格采样模式造成岩体核心区采样密度不足、外围无用区域采样冗余，后续圈定的富集区域未结合断裂、碱性岩接触带控矿要素分级，最终输出的靶区缺少量化优先级排序依据，导致勘查靶区精准度低、钻探工程布设成本高、资源潜力评估偏差大。

技术实现思路

1、鉴于此，本发明为了解决技术背景中存在的问题，提出了一种碳酸岩型稀土矿富集区判定方法。具体包括如下内容：

2、一种碳酸岩型稀土矿富集区判定方法，包在目标勘查区域布设网格采样点位，采集碳酸岩岩芯或露头岩石样品，对各样品开展la-icp-ms原位元素测试，获取各样品稀土总量r、轻稀土l、重稀土h、氟元素f、磷元素p实测含量，采集各采样点位三维空间坐标数据，将全部参数导入空间分析模块构建多参数数据集，基于数据集计算稀土分异判别指数，计算公式：

3、;

4、式中、分别为第i件样品镧、镱元素含量，n为参与统计的样品总数，为区域镧镱比值均值，选取同一区域无矿碳酸岩样品数据集确定基准阈值，提取目标区域内满足>的采样点位形成初始异常点位集，采用距离权重插值法对初始异常点位开展空间插值得到连续二维面数据，再构建稀土富集耦合系数，计算公式：

5、;

6、式中k为梯度调节系数，为区域碳酸岩稀土背景含量，依据数值区间划分梯度单元，将连续空间上同一梯度单元相连的闭合区域划定为碳酸岩型稀土矿潜在富集区。

7、在本发明的一些实施例中，所述网格采样点位布设时采用变间距网格模式，在碳酸岩岩体边界300m范围内布设50m×50m加密网格，岩体外围300m至1000m区间布设150m×150m常规网格，岩体外围超出1000m区域布设300m×300m稀疏网格，采样时每件样品采集新鲜无风化碳酸岩岩块，剔除蚀变程度高于45%的样品，每件样品选取3组不同矿物颗粒开展原位元素测试，取算术平均值作为该样品最终元素含量数据。

8、在本发明的一些实施例中，计算基准阈值时，选取区域内远离已知矿化地段、蚀变强度低于15%、无热液脉体穿插的碳酸岩样品不少于35件，统计该组样品全部数值，采用箱型图剔除离散异常值后，取上四分位数作为基准阈值。

9、在本发明的一些实施例中，所述权重插值法插值幂指数设置为2.2，插值邻域搜索采用八邻域模式，插值网格分辨率与野外采样最小网格间距保持一致，插值完成后采用3×3均值滤波消除插值产生的孤立噪点栅格。

10、在本发明的一些实施例中，所述梯度调节系数k取值区间设定0.0002～0.0008，采用未矿化碳酸岩样品稀土总量中位数，按照数值由低至高依次划分四级梯度单元，梯度单元边界采用50m缓冲区平滑处理。

11、在本发明的一些实施例中，划定所述潜在富集区之后，叠加区域断裂构造、碱性岩接触带空间矢量图层，对位于断裂交汇部位、碳酸岩-碱性岩接触带两侧200m范围内的闭合区域标记为一级富集靶区，其余闭合区域标记为二级富集靶区。

12、在本发明的一些实施例中，对所述一级富集靶区提取栅格单元内均值、面积、平均岩体埋深参数，按照参数排序完成靶区优先级排序，输出分级富集区矢量成果文件。

13、采用上述技术方案，具有如下有益效果：

14、本方法采用变间距网格采样，匹配碳酸岩岩体不同区域矿化发育程度，在矿化核心区加密采样不漏失小型矿化条带，外围稀疏采样直接减少野外采样样品数量与测试成本；

15、通过多组矿物颗粒原位测试取平均值，直接消除矿物微区元素分布不均造成的la-icp-ms测试数据偏差；构建稀土分异判别指数同步耦合轻重稀土分异、稀土总量、氟磷成矿挥发分、区域镧镱比值背景波动四类关键参数，直接规避单一元素指标评价产生的假异常；稀土富集耦合系数结合逻辑斯蒂函数校正区域稀土背景值，直接实现稀土分异程度与稀土总量的定量耦合，四级梯度单元搭配50m缓冲区平滑边界，直接输出边界规整、潜力分级清晰的潜在富集区；叠加断裂、碱性岩接触带矢量图层分级标记一、二级靶区，直接依托控矿地质要素区分高、低潜力富集区块；提取均值、面积、岩体埋深参数完成一级靶区量化排序，直接输出可直接用于钻探工程部署的分级矢量成果，无需人工二次整理评价数据。

1.一种碳酸岩型稀土矿富集区判定方法，其特征在于，包在目标勘查区域布设网格采样点位，采集碳酸岩岩芯或露头岩石样品，对各样品开展la-icp-ms原位元素测试，获取各样品稀土总量r、轻稀土l、重稀土h、氟元素f、磷元素p实测含量，采集各采样点位三维空间坐标数据，将全部参数导入空间分析模块构建多参数数据集，基于数据集计算稀土分异判别指数，计算公式：

2.根据权利要求1所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，所述网格采样点位布设时采用变间距网格模式，在碳酸岩岩体边界300m范围内布设50m×50m加密网格，岩体外围300m至1000m区间布设150m×150m常规网格，岩体外围超出1000m区域布设300m×300m稀疏网格，采样时每件样品采集新鲜无风化碳酸岩岩块，剔除蚀变程度高于45%的样品，每件样品选取3组不同矿物颗粒开展原位元素测试，取算术平均值作为该样品最终元素含量数据。

3.根据权利要求1所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，计算基准阈值时，选取区域内远离已知矿化地段、蚀变强度低于15%、无热液脉体穿插的碳酸岩样品不少于35件，统计该组样品全部数值，采用箱型图剔除离散异常值后，取上四分位数作为基准阈值。

4.根据权利要求1所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，所述权重插值法插值幂指数设置为2.2，插值邻域搜索采用八邻域模式，插值网格分辨率与野外采样最小网格间距保持一致，插值完成后采用3×3均值滤波消除插值产生的孤立噪点栅格。

5.根据权利要求1所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，所述梯度调节系数k取值区间设定0.0002～0.0008，采用未矿化碳酸岩样品稀土总量中位数，按照数值由低至高依次划分四级梯度单元，梯度单元边界采用50m缓冲区平滑处理。

6.根据权利要求1所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，划定所述潜在富集区之后，叠加区域断裂构造、碱性岩接触带空间矢量图层，对位于断裂交汇部位、碳酸岩-碱性岩接触带两侧200m范围内的闭合区域标记为一级富集靶区，其余闭合区域标记为二级富集靶区。

7.根据权利要求6所述的一种碳酸岩型稀土矿富集区判定方法，其特征在于，对所述一级富集靶区提取栅格单元内均值、面积、平均岩体埋深参数，按照参数排序完成靶区优先级排序，输出分级富集区矢量成果文件。