本发明属于铀矿勘查中放射性地球物理勘查领域，具体涉及一种放射性综合信息异常远景区圈定方法。
背景技术：
根据目前的钻探揭露，深部铀矿具有较大的找矿潜力，开展深部铀矿勘查技术研究与创新有着重要的意义。
深部隐伏铀矿中所含放射性核素在向地表传递的过程中不断衰变产生新的核素，这些放射性信息在传递过程中经历了围岩及盖层的层层阻隔，到达地表时异常信息十分微弱，部分信息被地表信息掩盖，且单一的信息异常容易产生多解性，给深部铀矿化异常解释带来较大困难。
采用多元信息综合异常提取是放射性地球物理方法找矿的发展趋势。通过多放射性异常信息集成提取综合信息异常来圈定异常远景靶区，对寻找深部铀矿具有重要的意义。
技术实现要素：
本发明需要解决的技术问题为：提供一种放射性综合信息异常远景区圈定方法，能够有效指示深部铀矿化信息，指导深部铀矿找矿。
本发明的技术方案如下所述：
一种放射性综合信息异常远景区圈定方法，包括以下步骤：
步骤1
收集工作区地质资料、铀矿体已知矿体面资料、γ能谱实测数据、氡及其子体实测数据；
步骤2
对γ能谱实测数据、氡及其子体实测数据进行整理，形成含公里网的空间散点数据，包括：含公里网的γ能谱铀空间散点数据、含公里网的γ能谱钍空间散点数据、含公里网的γ能谱钾空间散点数据，含公里网的土壤氡空间散点数据；
步骤3
将步骤2所述各含公里网的空间散点数据进行投影，并进行插值计算，形成空间栅格数据，包括：γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据；
步骤4
根据步骤3所述空间栅格数据，计算铀钍比、铀钾比、氡衬度；
步骤5
计算γ能谱铀空间栅格数据平均值mu和γ能谱铀空间栅格数据均方差su、γ能谱钍空间栅格数据平均值mt和γ能谱钍空间栅格数据均方差st、γ能谱钾空间栅格数据平均值mk和γ能谱钾空间栅格数据均方差sk、土壤氡空间栅格数据平均值mr和土壤氡空间栅格数据均方差sr、铀钍比平均值mut和铀钍比均方差sut、铀钾比平均值muk和铀钍比均方差suk、以及氡衬度平均值mrc和氡衬度均方差src；
步骤6
将步骤3所述γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据，以及步骤4所述铀钍比、铀钾比、氡衬度数据绘制成等值图并识别铀异常矢量面、钍异常矢量面、钾异常矢量面、氡异常矢量面、铀钍比异常矢量面、铀钾比异常矢量面、氡衬度异常矢量面；
步骤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]区间内，之后依次乘以cu、ct、ck、cr、cut、cuk、crc，计算结果依次记为gu、gt、gk、gr、gut、gtk、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绘制的等值线图，圈定异常远景区。
作为优选方案：步骤1中，所述γ能谱实测数据为地面γ能谱实测数据或航空γ能谱实测数据。
作为优选方案：步骤3中，基于gis处理平台将步骤2所述各含公里网的空间散点数据进行投影，并进行克里金插值，形成空间栅格数据。
作为优选方案：步骤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的区域识别为氡衬度异常矢量面。
本发明的有益效果为：
本发明的一种放射性综合信息异常远景区圈定方法，能够准确圈定深部铀矿化异常，为深部铀矿找矿和工程勘探提供参考依据。
附图说明
图1为采用本发明的放射性综合信息异常远景区圈定方法获取的放射性综合信息异常等值线及异常远景区图。
具体实施方式
下面结合说明书附图和实施例对本发明的一种放射性综合信息异常远景区圈定方法进行详细说明。
本发明的一种放射性综合信息异常远景区圈定方法，采用布尔逻辑提取γ能谱实测数据、氡及其子体实测数据、基于上述数据得到的其他数据的各个数据图层的异常面，以异常面和已知矿体在地表投影面的重叠区与已知矿体在地表投影面的比值为分值，对各个数据图层进行综合集成，利用集成后的综合信息异常来圈定异常远景区来预测深部铀矿资源。
具体而言，本发明的方法包括以下步骤：
步骤1
收集工作区地质资料、铀矿体已知矿体面资料、γ能谱实测数据、氡及其子体实测数据。所述γ能谱实测数据可以为地面γ能谱实测数据或航空γ能谱实测数据，如相山火山盆地放射性地面γ能谱实测数据。
步骤2
对γ能谱实测数据、氡及其子体实测数据进行整理，形成含公里网的空间散点数据，包括：含公里网的γ能谱铀空间散点数据、含公里网的γ能谱钍空间散点数据、含公里网的γ能谱钾空间散点数据，含公里网的土壤氡空间散点数据。
步骤3
基于gis处理平台将步骤2所述各含公里网的空间散点数据进行投影，并进行克里金插值，形成空间栅格数据，包括：γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据。
步骤4
根据步骤3所述空间栅格数据，计算铀钍比、铀钾比、氡衬度，计算方法为本领域技术人员公知常识。
步骤5
计算γ能谱铀空间栅格数据平均值mu和γ能谱铀空间栅格数据均方差su、γ能谱钍空间栅格数据平均值mt和γ能谱钍空间栅格数据均方差st、γ能谱钾空间栅格数据平均值mk和γ能谱钾空间栅格数据均方差sk、土壤氡空间栅格数据平均值mr和土壤氡空间栅格数据均方差sr、铀钍比平均值mut和铀钍比均方差sut、铀钾比平均值muk和铀钍比均方差suk、以及氡衬度平均值mrc和氡衬度均方差src。
步骤6
将步骤3所述γ能谱铀空间栅格数据、γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据，以及步骤4所述铀钍比、铀钾比、氡衬度数据绘制成等值图。
采用布尔逻辑法，识别异常矢量面，具体识别方法如下：
将γ能谱铀空间栅格数据≥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所述γ能谱铀空间栅格数据拉伸或压缩至[0,1]区间内，之后乘以cu，计算结果记为gu。同理，将步骤3所述γ能谱钍空间栅格数据、γ能谱钾空间栅格数据、土壤氡空间栅格数据，以及步骤4所述铀钍比、铀钾比、氡衬度数据分别拉伸或压缩至[0,1]区间内，之后依次乘以ct、ck、cr、cut、cuk、crc，计算结果依次记为gt、gk、gr、gut、gtk、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绘制的等值线图，圈定图1所示异常远景区。
上面结合附图和实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做各种变化。本发明中未作详细描述的内容均可以采用现有技术。