本發(fā)明屬于放射性物探數(shù)據(jù)校正
技術(shù)領(lǐng)域：
，具體涉及一種地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)調(diào)平校正方法。
背景技術(shù)：
：目前，在鈾礦勘查領(lǐng)域中，普遍采用地面伽瑪能譜測(cè)量獲取地表土壤或巖石中的鈾含量，通過圈定鈾異常區(qū)來預(yù)測(cè)鈾成礦靶區(qū)。當(dāng)勘查區(qū)面積較大時(shí)，需要將勘查區(qū)劃分為幾個(gè)分測(cè)區(qū)，并在不同年度來完成測(cè)量任務(wù)。由于地面伽瑪能譜測(cè)量受氣候條件等影響，導(dǎo)致不同年度相鄰分測(cè)區(qū)的地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)存在一定差別，因此需要對(duì)各分測(cè)區(qū)的地面伽瑪能譜測(cè)量的鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正，從而有效地去除氣候條件等對(duì)鈾含量的影響，獲取整個(gè)勘查區(qū)地表土壤或巖石中鈾含量分布特征。技術(shù)實(shí)現(xiàn)要素：本發(fā)明需要解決的技術(shù)問題為：現(xiàn)有技術(shù)難以對(duì)不同年度相鄰測(cè)區(qū)的地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)進(jìn)行有效調(diào)平校正。本發(fā)明的技術(shù)方案如下所述：一種地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)調(diào)平校正方法，包括以下步驟：步驟1.鈾含量重復(fù)測(cè)量：將不同年度在分測(cè)區(qū)進(jìn)行鈾含量原始測(cè)量時(shí)設(shè)置的測(cè)點(diǎn)作為原始測(cè)點(diǎn)，之后在同一年度對(duì)所有分測(cè)區(qū)進(jìn)行重復(fù)鈾含量測(cè)量，重復(fù)鈾含量測(cè)量過程中設(shè)置的測(cè)點(diǎn)作為重復(fù)測(cè)點(diǎn)，每個(gè)分測(cè)區(qū)中的重復(fù)測(cè)點(diǎn)均為該分測(cè)區(qū)原始測(cè)點(diǎn)中的一部分，即在該分測(cè)區(qū)部分原始測(cè)點(diǎn)同位置處設(shè)置重復(fù)測(cè)點(diǎn)，重復(fù)測(cè)點(diǎn)在該分測(cè)區(qū)均勻分布；步驟2.計(jì)算鈾含量背景值：分別對(duì)每個(gè)分測(cè)區(qū)計(jì)算重復(fù)測(cè)量時(shí)鈾含量背景值和原始測(cè)量時(shí)鈾含量背景值；步驟3.計(jì)算鈾含量調(diào)平校正系數(shù)：根據(jù)同一分測(cè)區(qū)重復(fù)測(cè)量時(shí)鈾含量背景值和原始測(cè)量時(shí)鈾含量背景值相等的關(guān)系、重復(fù)測(cè)點(diǎn)測(cè)得數(shù)據(jù)與同位置經(jīng)調(diào)平校正后的原始測(cè)點(diǎn)測(cè)得數(shù)據(jù)之間的誤差為最小的關(guān)系，計(jì)算調(diào)平校正系數(shù)；步驟4.對(duì)分測(cè)區(qū)鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正：采用調(diào)平校正系數(shù)對(duì)分測(cè)區(qū)鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正。作為優(yōu)選方案：步驟2中，采用“逐步替代剔除法”分別對(duì)每個(gè)分測(cè)區(qū)計(jì)算重復(fù)測(cè)量時(shí)鈾含量背景值和原始測(cè)量時(shí)鈾含量背景值，具體包括以下步驟：步驟2.1.計(jì)算重復(fù)測(cè)量時(shí)鈾含量背景值采用下式計(jì)算重復(fù)測(cè)點(diǎn)鈾含量平均值、鈾含量標(biāo)準(zhǔn)偏差和鈾含量異常下限值：式中，MY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量平均值；SY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量標(biāo)準(zhǔn)偏差；AY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量異常下限值；Yi(k)為第k次迭代采用的第i個(gè)重復(fù)測(cè)點(diǎn)的鈾含量值，其中，Yi(1)為第i個(gè)重復(fù)測(cè)點(diǎn)測(cè)得的鈾含量值；n為該分測(cè)區(qū)重復(fù)測(cè)點(diǎn)個(gè)數(shù)；式(1)計(jì)算中，迭代計(jì)算終止時(shí)的重復(fù)測(cè)點(diǎn)鈾含量平均值即為該分測(cè)區(qū)重復(fù)測(cè)量時(shí)鈾含量背景值BY；步驟2.2.計(jì)算原始測(cè)量時(shí)鈾含量背景值采用下式計(jì)算原始測(cè)點(diǎn)鈾含量平均值、鈾含量標(biāo)準(zhǔn)偏差和鈾含量異常下限值：式中，MX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量平均值；SX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量標(biāo)準(zhǔn)偏差；AX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量異常下限值；Xi(k)為第k次迭代采用的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)的鈾含量值，其中，Xi(1)為第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值；式(2)計(jì)算中，迭代計(jì)算終止時(shí)的原始測(cè)點(diǎn)鈾含量平均值即為該分測(cè)區(qū)原始測(cè)量時(shí)鈾含量背景值BX。作為優(yōu)選方案：步驟3具體包括以下步驟：根據(jù)原始測(cè)量時(shí)鈾含量背景值與重復(fù)測(cè)量時(shí)鈾含量背景值相等的關(guān)系建立函數(shù)關(guān)系式：令BY＝aBX+b(3)式中，a、b均為調(diào)平校正系數(shù)；采用式(4)對(duì)與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值Xi(1)進(jìn)行校正：Xci＝aXi(1)+bi＝1，2,...,n(4)式中，Xci為校正后的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值；采用式(5)計(jì)算Xci與第i個(gè)重復(fù)測(cè)點(diǎn)測(cè)得的鈾含量值Yi(1)的誤差E：E=Σi=1n(Xci-Yi(1))2i=1,2,...,n---(5)]]>將式(3)、式(4)代入式(5)得到式(6)：E=Σi=1n[a(Xi(1)-BX)+BY-Yi(1)]2---(6)]]>根據(jù)重復(fù)測(cè)點(diǎn)測(cè)得數(shù)據(jù)與經(jīng)調(diào)平校正后同位置的原始測(cè)點(diǎn)測(cè)得數(shù)據(jù)之間的誤差為最小的關(guān)系，計(jì)算校正系數(shù)：令∂E∂a=0---(7)]]>根據(jù)式(7)計(jì)算調(diào)平校正系數(shù)a，將a值代入式(3)計(jì)算調(diào)平校正系數(shù)b。作為優(yōu)選方案：步驟4具體包括以下步驟：根據(jù)式(4)計(jì)算出校正后的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值即完成了對(duì)該分測(cè)區(qū)鈾含量數(shù)據(jù)的調(diào)平校正。作為優(yōu)選方案：步驟1中，對(duì)于每個(gè)分測(cè)區(qū)，該分測(cè)區(qū)重復(fù)測(cè)點(diǎn)數(shù)量均大于等于該分測(cè)區(qū)原始測(cè)點(diǎn)數(shù)量的20％。作為優(yōu)選方案：步驟1中，對(duì)所有分測(cè)區(qū)的重復(fù)鈾含量測(cè)量均采用地面 伽瑪能譜測(cè)量方法。本發(fā)明的有益效果為：本發(fā)明的一種地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)調(diào)平校正方法，能夠?qū)Σ煌甓认噜彍y(cè)區(qū)地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)進(jìn)行有效調(diào)平校正，獲取整個(gè)勘查區(qū)的鈾含量分布特征，指導(dǎo)鈾礦勘探工作。具體實(shí)施方式下面結(jié)合實(shí)施例對(duì)本發(fā)明的一種地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)調(diào)平校正方法進(jìn)行詳細(xì)說明。一種地面伽瑪能譜測(cè)量鈾含量數(shù)據(jù)調(diào)平校正方法，包括以下步驟：步驟1.鈾含量重復(fù)測(cè)量：將現(xiàn)有技術(shù)中不同年度在分測(cè)區(qū)進(jìn)行鈾含量原始測(cè)量時(shí)設(shè)置的測(cè)點(diǎn)作為原始測(cè)點(diǎn)，之后在同一年度對(duì)所有分測(cè)區(qū)采用現(xiàn)有技術(shù)中的測(cè)量方法進(jìn)行重復(fù)鈾含量測(cè)量，重復(fù)鈾含量測(cè)量過程中設(shè)置的測(cè)點(diǎn)作為重復(fù)測(cè)點(diǎn)，每個(gè)分測(cè)區(qū)中的重復(fù)測(cè)點(diǎn)均為該分測(cè)區(qū)原始測(cè)點(diǎn)中的一部分，即在該分測(cè)區(qū)部分原始測(cè)點(diǎn)同位置處設(shè)置重復(fù)測(cè)點(diǎn)，重復(fù)測(cè)點(diǎn)在該分測(cè)區(qū)均勻分布；步驟2.計(jì)算鈾含量背景值：采用“逐步替代剔除法”分別對(duì)每個(gè)分測(cè)區(qū)計(jì)算重復(fù)測(cè)量時(shí)鈾含量背景值和原始測(cè)量時(shí)鈾含量背景值；步驟3.計(jì)算鈾含量調(diào)平校正系數(shù)：根據(jù)相同分測(cè)區(qū)重復(fù)測(cè)量時(shí)鈾含量背景值和原始測(cè)量時(shí)鈾含量背景值相等的關(guān)系、重復(fù)測(cè)點(diǎn)測(cè)得數(shù)據(jù)與經(jīng)調(diào)平校正后的原始測(cè)點(diǎn)測(cè)得數(shù)據(jù)之間的誤差為最小的關(guān)系，計(jì)算調(diào)平校正系數(shù)；步驟4.對(duì)分測(cè)區(qū)鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正：采用校正系數(shù)對(duì)分測(cè)區(qū)鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正。實(shí)施例1步驟1.鈾含量重復(fù)測(cè)量在同一年度對(duì)所有分測(cè)區(qū)采用現(xiàn)有技術(shù)中的地面伽瑪能譜測(cè)量方法進(jìn)行 重復(fù)鈾含量測(cè)量，對(duì)于每個(gè)分測(cè)區(qū)，重復(fù)測(cè)點(diǎn)數(shù)量均需大于等于該分測(cè)區(qū)原始測(cè)點(diǎn)數(shù)量的20％。步驟2.計(jì)算鈾含量背景值針對(duì)每一個(gè)分測(cè)區(qū)，分別計(jì)算該分測(cè)區(qū)內(nèi)所有重復(fù)測(cè)點(diǎn)和與其同位置的原始測(cè)點(diǎn)的背景值。步驟2.1.計(jì)算重復(fù)測(cè)量時(shí)鈾含量背景值采用下式計(jì)算重復(fù)測(cè)點(diǎn)鈾含量平均值、標(biāo)準(zhǔn)偏差和異常下限值：式中，MY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量平均值；SY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量標(biāo)準(zhǔn)偏差；AY(k)為第k次迭代得到的重復(fù)測(cè)點(diǎn)鈾含量異常下限值；Yi(k)為第k次迭代采用的第i個(gè)重復(fù)測(cè)點(diǎn)的鈾含量值，其中，Yi(1)為第i個(gè)重復(fù)測(cè)點(diǎn)測(cè)得的鈾含量值；n為該分測(cè)區(qū)重復(fù)測(cè)點(diǎn)個(gè)數(shù)。式(1)計(jì)算中，當(dāng)k＝p時(shí)，不再存在大于AY(p)的Yi(p)數(shù)據(jù)，迭代計(jì)算終止。迭代計(jì)算終止時(shí)的鈾含量平均值MY(p)即為該分測(cè)區(qū)重復(fù)測(cè)量時(shí)鈾含量背景值BY。步驟2.2.計(jì)算原始測(cè)量時(shí)鈾含量背景值采用下式計(jì)算原始測(cè)點(diǎn)鈾含量平均值、標(biāo)準(zhǔn)偏差和異常下限值：式中，MX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量平均值；SX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量標(biāo)準(zhǔn)偏差；AX(k)為第k次迭代得到的原始測(cè)點(diǎn)鈾含量異常下限值；Xi(k)為第k次迭代采用的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)的鈾含量值，其中，Xi(1)為第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值.式(2)計(jì)算中，當(dāng)k＝q時(shí)，不再存在大于AX(q)的Xi(q)數(shù)據(jù)，迭代計(jì)算終止。迭代計(jì)算終止時(shí)的鈾含量平均值MX(q)即為該分測(cè)區(qū)原始測(cè)量時(shí)鈾含量背景值BX。步驟3.計(jì)算鈾含量調(diào)平校正系數(shù)根據(jù)原始測(cè)量時(shí)鈾含量背景值與重復(fù)測(cè)量時(shí)鈾含量背景值相等的關(guān)系建立函數(shù)關(guān)系式：令BY＝aBX+b(3)式中，a、b均為調(diào)平校正系數(shù)。采用式(4)對(duì)與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值Xi(1)進(jìn)行校正：Xci＝aXi(1)+bi＝1，2,...,n(4)式中，Xci為校正后的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值。采用式(5)計(jì)算Xci與第i個(gè)重復(fù)測(cè)點(diǎn)測(cè)得的鈾含量值Yi(1)的誤差E：E=Σi=1n(Xci-Yi(1))2i=1,2,...,n---(5)]]>將式(3)、式(4)代入式(5)得到式(6)：E=Σi=1n[a(Xi(1)-BX)+BY-Yi(1)]2---(6)]]>根據(jù)重復(fù)測(cè)點(diǎn)測(cè)得數(shù)據(jù)與經(jīng)調(diào)平校正后同位置的原始測(cè)點(diǎn)測(cè)得數(shù)據(jù)之間的誤差為最小的關(guān)系，計(jì)算校正系數(shù)：令∂E∂a=0---(7)]]>根據(jù)式(7)計(jì)算調(diào)平校正系數(shù)a，將a值代入式(3)計(jì)算調(diào)平校正系數(shù)b。步驟4.對(duì)分測(cè)區(qū)鈾含量數(shù)據(jù)進(jìn)行調(diào)平校正根據(jù)式(4)計(jì)算出校正后的與第i個(gè)重復(fù)測(cè)點(diǎn)同位置的原始測(cè)點(diǎn)測(cè)得的鈾含量值即完成了對(duì)該分測(cè)區(qū)鈾含量數(shù)據(jù)的調(diào)平校正。其他分測(cè)區(qū)采用與上述方法逐一進(jìn)行調(diào)平校正。當(dāng)前第1頁1 2 3