[發明專利]一種放射性綜合信息異常遠景區圈定方法

有效

• 申請號: 201611217118.5
• 申請日: 2016-12-26
• 公開（公告）號: CN108241176B
• 公開（公告）日: 2019-08-13
• 發明（設計）人: 李必紅;程紀星;楊龍泉;趙丹;張玉燕
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G01V5/00
• 分類號: G01V5/00
• 代理公司: 核工業專利中心 11007
• 代理人: 包海燕

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繪制的等值線圖，圈定異常遠景區。