[发明专利]一种放射性综合信息异常远景区圈定方法

摘要

本发明属于铀矿勘查中放射性地球物理勘查领域，具体涉及一种放射性综合信息异常远景区圈定方法。采用布尔逻辑提取γ能谱实测数据、氡及其子体实测数据、基于上述数据得到的其他数据的各个数据图层的异常面，以异常面和已知矿体在地表投影面的重叠区与已知矿体在地表投影面的比值为分值，对各个数据图层进行综合集成，利用集成后的综合信息异常来圈定异常远景区来预测深部铀矿资源。本发明提供了一种放射性综合信息异常远景区圈定方法，能够准确圈定深部铀矿化异常，为深部铀矿找矿和工程勘探提供参考依据。

主权项

1.一种放射性综合信息异常远景区圈定方法，其特征在于：包括以下步骤：步骤1收集工作区地质资料、铀矿体已知矿体面资料、γ能谱实测数据、氡及其子体实测数据；步骤2对γ能谱实测数据、氡及其子体实测数据进行整理，形成含公里网的空间散点数据，包括：含公里网的γ能谱铀空间散点数据、含公里网的γ能谱钍空间散点数据、含公里网的γ能谱钾空间散点数据，含公里网的土壤氡空间散点数据；步骤3将步骤2所述各含公里网的空间散点数据进行投影，并进行插值计算，形成空间栅格数据，包括：γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据；步骤4根据步骤3所述空间栅格数据，计算铀钍比、铀钾比、氡衬度；步骤5计算γ能谱铀空间栅格数据平均值MU和γ能谱铀空间栅格数据均方差SU、γ能谱钍空间栅格数据平均值MT和γ能谱钍空间栅格数据均方差ST、γ能谱钾空间栅格数据平均值MK和γ能谱钾空间栅格数据均方差SK、土壤氡空间栅格数据平均值MR和土壤氡空间栅格数据均方差SR、铀钍比平均值MUT和铀钍比均方差SUT、铀钾比平均值MUK和铀钍比均方差SUK、以及氡衬度平均值MRC和氡衬度均方差SRC；步骤6将步骤3所述γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据，以及步骤4所述铀钍比、铀钾比、氡衬度数据绘制成等值图并识别铀异常矢量面、钍异常矢量面、钾异常矢量面、氡异常矢量面、铀钍比异常矢量面、铀钾比异常矢量面、氡衬度异常矢量面；步骤6中，采用布尔逻辑法，识别异常矢量面，具体识别方法如下：将γ能谱铀空间栅格数据≥MU+1.5×SU的区域识别为铀异常矢量面；将γ能谱钍空间栅格数据≥MT+1.5×ST的区域识别为钍异常矢量面；将γ能谱钾空间栅格数据≥MK+1.5×SK的区域识别为钾异常矢量面；将γ能谱氡空间栅格数据≥MR+1.5×SR的区域识别为氡异常矢量面；将铀钍比≥MUT+1.5×SUT的区域识别为铀钍比异常矢量面；将铀钾比≥MUK+1.5×SUK的区域识别为铀钾比异常矢量面；将氡衬度≥MRC+1.5×SRC的区域识别为氡衬度异常矢量面；步骤7获取步骤1所述已知矿体面与铀异常矢量面的重叠区面积AU、已知矿体面与钍异常矢量面的重叠区面积AT、已知矿体面与钾异常矢量面的重叠区面积AK、已知矿体面与氡异常矢量面的重叠区面积AR、已知矿体面与铀钍比异常矢量面的重叠区面积AUT、已知矿体面与铀钾比异常矢量面的重叠区面积AUK、已知矿体面与氡衬度异常矢量面的重叠区面积ARC；将已知矿体面面积记为Am，计算各重叠区面积与已知矿体面面积的比值：CU＝AU/Am；CT＝AT/Am；CK＝AK/Am；CR＝AR/Am；CUT＝AUT/Am；CUK＝AUK/Am；CRC＝ARC/Am；步骤8将步骤3所述γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据，以及步骤4所述铀钍比、铀钾比、氡衬度数据分别拉伸或压缩至[0,1]区间内，γ能谱铀空间栅格数据拉伸或压缩至[0,1]区间内之后，乘以CU，得到GU；γ能谱钍空间栅格数据拉伸或压缩至[0,1]区间内之后，乘以CT，得到GT；γ能谱钾空间栅格数据拉伸或压缩至[0,1]区间内之后，乘以CK，得到GK；土壤氡空间栅格数据拉伸或压缩至[0,1]区间内之后，乘以CR，得到GR；铀钍比拉伸或压缩至[0,1]区间内之后，乘以CUT，得到GUT；铀钾比拉伸或压缩至[0,1]区间内之后，乘以CUK，得到GTK；氡衬度拉伸或压缩至[0,1]区间内之后，乘以CRC，得到GRC；计算Q＝(GU+GT+GK+GR+GUT+GTK+GRC)/(CU+CT+CK+CR+CUT+CUK+CRC)，对Q进行滤波和滑动平均处理得到栅格数据S；将栅格数据S绘制成等值线图，统计栅格数据S的平均值MS和均方差SS；对于S值≥MS+0.5×SS的区域为圈定为深部成矿有利区；结合已知矿床、矿点资料，利用步骤8绘制的等值线图，圈定异常远景区。