一种火山岩型铀矿三维空间定位预测方法

摘要

本发明属于铀矿勘查技术领域，具体涉及一种火山岩型铀矿三维空间定位预测方法,步骤1、获取工作区物化探数据；步骤2、圈定放射性异常、地球化学元素异常和高磁异常；步骤3、对圈定的物化探异常区与已知矿体面进行相交分析；步骤4：数据归一化处理，获得综合信息分值；步骤5：圈定综合信息分值异常区；步骤6：获得电阻率低阻区和陡变区，根据音频大地电磁测深数据圈定低阻异常区和电阻率陡变区；步骤7：进行铀成矿空间三维定位预测。本发明方法以物化探异常为依据，综合深部铀矿化信息和成矿环境实现了火山岩型铀矿成矿三维空间定位预测，为深部火山型铀矿勘查和资源量预测提供了技术支撑。

说明书

技术领域

本发明属于铀矿勘查技术领域，具体涉及一种火山岩型铀矿三维空间定位预测方法。

背景技术

火山岩型铀矿是我国一种重要的热液型铀矿，该类铀矿浅层资源已基本查明；由于其深部成矿环境较为复杂，深部铀资源勘查还存在较多的技术“瓶颈”，致使深部铀矿定位及资源量预测精度较低。在火山岩型铀矿成矿定位预测方面，国内外目前现有技术主要依据放射性异常和地球化学元素异常方法圈定成矿有利区，基本都是平面预测，但是随着科学技术的进步和多学科联合攻关，设计一种能实现多维综合勘查与空间定位的铀矿勘查方法势在必行，目前三维数值处理和可视化技术的飞速发展，使得建立三维可视化空间定位预测方法成为可能。

因此，根据上述现有技术的不足，亟需设计一种火山岩型铀矿三维空间定位预测方法以实现通过多维综合勘查与空间定位进行铀矿勘查的技术。

发明内容

本发明目的是针对现有技术的不足，提供一种火山岩型铀矿三维空间定位预测方法，用于解决现有铀矿勘查技术由于多采用平面预测的方法对成矿环境复杂的地区铀矿勘测准确度较低、深部火山岩型铀矿定位难的技术问题。

本发明的技术方案：

一种火山岩型铀矿三维空间定位预测方法，包括如下步骤：

步骤1、获取工作区物化探数据；

步骤2、圈定放射性异常、地球化学元素异常和高磁异常，对步骤1中获取的数据进行处理，圈定出土壤氡浓度异常区、γ能谱含量异常区、土壤地球化学元素异常区及高磁异常区；

步骤3、对圈定的物化探异常区与已知矿体面进行相交分析；

步骤4：数据归一化处理，获得综合信息分值；

步骤5：圈定综合信息分值异常区；

步骤6：获得电阻率低阻区和陡变区，根据音频大地电磁测深数据圈定低阻异常区和电阻率陡变区；

步骤7：进行铀成矿空间三维定位预测。

所述步骤1还包括：选定工作区并通过收集或实测获取覆盖工作区的物化探数据和已知矿床矿点数据资料，覆盖工作区的物化探数据和已知矿床矿点数据资料包括:高精度磁测数据、音频大地电磁数据、土壤氡浓度数据、地面伽马能谱数据、地球化学元素含量数据、已知矿体面等。

所述步骤二和步骤三中的γ能谱包括：U、Th、K；所述步骤二和步骤三中的土壤地球化学元素包括：U、Mo、Be元素。

所述步骤3还包括：通过与已知矿体面的相交分析，获得土壤氡浓度异常区、γ能谱含量异常区、土壤地球化学元素异常区及高磁异常区与已知矿体面的重叠区面积，并将这些重叠区的面积分别除以铀矿体面的面积，获得相应比值。

所述步骤四还包括：将步骤3中的土壤氡浓度、γ能谱含量、土壤地球化学元素含量及磁测数据按照如下公式(1)归一化至[0,1]区间内，再按如下公式(2)进行综合，获得综合打分值；

上述公式(1)中，X

上述公式(2)中，k

所述步骤5还包括：根据步骤4处理的综合信息分值数据，采用累积频率法，按照80％累频确定综合信息分值的异常下限，根据异常下限圈定综合信息分值异常区。

所述步骤7还包括：采用三维制图软件，以公里网坐标为三维体坐标系统，分别导入步骤6中的电阻率陡变区和步骤5中圈定的综合信息分值异常区，将综合信息分值异常区沿深度H方向延伸形成综合信息异常柱体，将综合信息异常柱体和电阻率陡变区重叠的三维空间体划定为火山岩型铀矿成矿三维空间靶区，并根据低阻异常的走向预测矿体规模。

本发明的有益效果：

本发明设计的一种火山岩型铀矿三维空间定位预测方法，解决了深部火山岩型铀矿定位难的问题，本发明方法以物化探异常为依据，综合深部铀矿化信息和成矿环境实现了火山岩型铀矿成矿三维空间定位预测，为深部火山型铀矿勘查和资源量预测提供了技术支撑。

附图说明

图1为本发明实施例中河元背工作区预测铀矿体空间位置示意图

具体实施方式

下面结合附图及实施例对本发明的一种铀钼矿床靶区圈定方法进行详细说明。

一种火山岩型铀矿三维空间定位预测方法，包括如下步骤：

步骤1、获取工作区物化探数据；选定工作区并通过收集或实测获取覆盖工作区的物化探数据和已知矿床矿点数据资料，覆盖工作区的物化探数据和已知矿床矿点数据资料包括:高精度磁测数据、音频大地电磁数据、土壤氡浓度数据、地面伽马能谱数据、地球化学元素含量数据、已知矿体面等。

步骤2、圈定放射性异常、地球化学元素异常和高磁异常，对步骤1中获取的数据进行处理，圈定出土壤氡浓度异常区、γ能谱(U、Th、K)含量异常区、土壤地球化学元素(U、Mo、Be)异常区及高磁异常区；

步骤3、对圈定的物化探异常区与已知矿体面进行相交分析；通过与已知矿体面的相交分析，获得土壤氡浓度异常区、γ能谱含量异常区、土壤地球化学元素异常区及高磁异常区与已知矿体面的重叠区面积，并将这些重叠区的面积分别除以铀矿体面的面积，获得相应比值。

步骤4：数据归一化处理，获得综合信息分值；将步骤3中的土壤氡浓度、γ能谱含量、土壤地球化学元素含量及磁测数据按照如下公式(1)归一化至[0,1]区间内，再按如下公式(2)进行综合，获得综合打分值；

上述公式(1)中，X

上述公式(2)中，k

步骤5：圈定综合信息分值异常区；根据步骤4处理的综合信息分值数据，采用累积频率法，按照80％累频确定综合信息分值的异常下限，根据异常下限圈定综合信息分值异常区。

步骤6：获得电阻率低阻区和陡变区，根据音频大地电磁测深数据圈定低阻异常区和电阻率陡变区；

步骤7：进行铀成矿空间三维定位预测：采用三维制图软件，以公里网坐标为三维体坐标系统，分别导入步骤6中的电阻率陡变区和步骤5中圈定的综合信息分值异常区，将综合信息分值异常区沿深度H方向延伸形成综合信息异常柱体，将综合信息异常柱体和电阻率陡变区重叠的三维空间体划定为火山岩型铀矿成矿三维空间靶区，并根据低阻异常的走向预测矿体规模。

实施例：

步骤1、获取工作区物化探数据：

选择相山盆地，收集已知矿体面数据资料，通过实测获取物化探数据资料。在相山盆地按照1:5万比例尺开展高精度物化探勘查工作，主要包括音频大地电磁(AMT)、高精度磁法、氡气测量、伽马能谱测量、土壤地球化学测量。物化探测线方向NW313°，线距500m，点距100m，共设计测量剖面51条。在河元背、居隆庵、邹家山和云际等4个重点工作区加密至线距100m，点距50m。通过数据处理，形成了相山火山盆地和重点工作区物化探基础图件，这些图件包括：高精度磁测ΔT异常图、土壤氡浓度异常图、地面伽马能谱异常图(U、Th、K)、土壤地球化学元素异常图(U、Mo、Be)以及音频大地电磁二维反演电阻率断面和地质解译图。

步骤2、圈定放射性异常、地球化学元素异常和高磁异常：

对步骤1中获取的数据进行处理，圈定出土壤氡浓度异常区、γ能谱(U、Th、K)含量异常区、土壤地球化学元素(U、Mo、Be)异常区及高磁异常区。

步骤3、对圈定的物化探异常区与已知矿体面进行相交分析：

通过与已知矿体面的相交分析，获得土壤氡浓度异常区、γ能谱(U、Th、K)含量异常区、土壤地球化学元素(U、Mo、Be)异常区及高磁异常区与已知矿体面的重叠区面积，并将这些重叠区的面积分别除以铀矿体面的面积，获得相应的比值。

步骤4：数据归一化处理，获得综合信息分值：

将步骤3中的土壤氡浓度、γ能谱(U、Th、K)含量、土壤地球化学元素(U、Mo、Be)含量及磁测数据按照公式(1)归一化至[0,1]区间内，再按公式(2)进行综合，获得综合打分值；

上述公式(1)中，X

上述公式(2)中，k

步骤5：圈定综合信息分值异常区：

根据步骤4处理的综合信息分值数据，采用累积频率法，按照80％累频确定综合信息分值的异常下限，根据异常下限圈定综合信息分值异常区。

步骤6：获得电阻率低阻区和陡变区：

根据音频大地电磁测深数据圈定低阻异常区和电阻率陡变区。

步骤7：进行铀成矿空间三维定位预测：

采用三维制图软件(GOCAD)，以公里网坐标(X，Y，H)为三维体坐标系统，分别导入步骤6中的电阻率陡变区和步骤5中圈定的综合信息分值异常区，将综合信息分值异常区沿深度H方向延伸形成综合信息异常柱体，将综合信息异常柱体和电阻率陡变区重叠的三维空间体划定为火山岩型铀矿成矿三维空间靶区，并根据低阻异常的走向预测矿体规模。

上面结合实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可采用现有方法技术。