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

权利要求书

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绘制的等值线图，圈定异常远景区。