田宇，常樹帥，劉丹，陳江源，魏永強

（1.核工業(yè)航測遙感中心，河北 石家莊，050002；2.國家核應(yīng)急航空監(jiān)測中心，河北 石家莊，050002；3.中核集團鈾資源地球物理勘查技術(shù)中心重點實驗室，河北 石家莊，050002）

2019 年，中國核工業(yè)地質(zhì)局下達了“二連盆地地質(zhì)資料綜合開發(fā)利用及編圖”項目，該項目的工作內(nèi)容包括航磁、航放數(shù)據(jù)集成處理，并進行綜合成礦預測。在對已知砂巖型鈾礦與航放數(shù)據(jù)進行疊加分析時，發(fā)現(xiàn)并非所有鈾礦都位于航放高場內(nèi)，還有部分鈾礦位于航放高場邊緣；而在數(shù)字高程模型（DEM）數(shù)據(jù)上，絕大部分鈾礦處于地形由高向低變化部位或處于相對坡底部位（圖1）。受微地形的影響，已知鈾礦區(qū)航放鈾含量分布特征更為復雜。為分析產(chǎn)生上述情況的原因，以及充分挖掘航放數(shù)據(jù)與DEM 數(shù)據(jù)所蘊含的與成礦有關(guān)信息，筆者試圖通過對航放鈾含量與微地形進行相關(guān)分析，在航放資料處理解釋方面進行探討，以期深化對航放異常的認識，更好地提高航放測量的地質(zhì)效果。

圖1 二連盆地烏蘭察布坳陷航放鈾含量與DEM 數(shù)據(jù)疊合影像圖Fig.1 The superimposed image of airborne uranium content and DEM data in Ulanchabu depression of Erlian basin

1 研究基礎(chǔ)

1.1 微地形定義及其研究意義

微地形簡而言之就是地表的微弱起伏，即相對高程變化小的一種起伏地表，在此特指1 km 范圍內(nèi)的地形變化不超過50 m。對于砂巖型鈾礦而言，其主要產(chǎn)于長時期接受沉降且成礦后構(gòu)造運動并不強烈的盆地內(nèi)，因其地表起伏并不強烈，所以應(yīng)用了微地形概念用以研究其與砂巖型鈾礦成礦之間的關(guān)系。

砂巖型鈾礦是外生后成鈾礦床，根據(jù)砂巖型鈾礦成礦模式，含鈾含氧水沿砂巖層向前運移過程中，不斷消耗氧而使鈾沉淀下來形成主要鈾礦體，或深部的油氣通過斷裂到達砂巖層，使原有的氧化環(huán)境變成還原環(huán)境沉淀出部分鈾礦體［1-3］。對于外生成礦來說，無論是什么類型的礦床，盡管其具體的成礦作用各不相同，但總的來說都是以水作為載體對鈾礦質(zhì)從源、遷移、匯（成礦）到保存4 個環(huán)節(jié)相互嵌套的動力學過程［4］，而地表水及潛水的發(fā)育及流動與原始地形地貌密切相關(guān)。

我國北方主要產(chǎn)鈾盆地砂巖型鈾礦的成礦年齡都比含礦層位的沉積時代年輕很多，成礦時代大部分集中在古近紀和新近紀，部分盆地的成礦時代延伸至了第四紀，具有“幕式”及多階段特征［5-6］。因此，可以推斷現(xiàn)代地形地貌對于砂巖型鈾礦的富集成礦有潛在的影響。

近年來，相關(guān)學者也曾利用數(shù)字高程模型（DEM）數(shù)據(jù)進行過鈾成礦規(guī)律的研究，如核工業(yè)二七〇研究所張萬良等人基于DEM 的面積-高程曲線計算，研究了硬巖型鈾礦相山鈾礦田地貌形體與鈾礦空間分布關(guān)系，并得出了鈾礦空間分布與地貌形體有一定關(guān)系的結(jié)論，即鈾礦床礦點的分布與地表侵蝕程度有關(guān)［7］。東華理工大學孫岳等人研究了烏茲別克斯坦中卡茲庫姆地區(qū)砂巖型鈾礦與DEM 數(shù)據(jù)之間的空間關(guān)聯(lián)，利用DEM數(shù)據(jù)提取了高程、坡度、坡向、地形起伏度等各類地貌因子，并進行了定量分析，研究結(jié)果表明中卡茲庫姆地區(qū)砂巖型鈾礦發(fā)育于Ⅳ級夷平面附近、集中于低海拔和低地形起伏度的有利區(qū)域，并根據(jù)地貌因子圈定了砂巖型鈾礦賦存的有利區(qū)域［8］。

1.2 相關(guān)分析理論基礎(chǔ)

相關(guān)分析在此指相關(guān)系數(shù)分析。相關(guān)系數(shù)是一個界于-1 與＋1 之間的無量綱指數(shù)，反映兩個數(shù)據(jù)集合之間的線性相關(guān)程度。相關(guān)系數(shù)大于零，表示兩變量正相關(guān)，因變量有隨自變量同步增加或同步減少的趨勢。相關(guān)系數(shù)小于零，表示兩變量負相關(guān)，因變量有隨自變量增加而減少（或因變量隨自變量減少而增加）的趨勢。公式為：

式中：r為相關(guān)系數(shù)值；x1，x2為相關(guān)分析的兩變量；n為對比取樣點個數(shù)，即匹配數(shù)。

相關(guān)分析是放射性物探工作中常用的一種數(shù)理統(tǒng)計方法，用以研究某一變量與其他變量之間的關(guān)系。例如巖石或礦石中鈾與釷之間的關(guān)系，放射性元素與其他伴生元素之間的關(guān)系，兩個鉆孔的測井資料的對比、地層剖面對比、鉆孔柱狀圖的對比等。

2 砂巖型鈾礦航放鈾含量與微地形相關(guān)分析

既然已知鈾礦區(qū)受微地形的影響，航放鈾含量可能顯示為高場，也可能顯示為高場的邊緣，使鈾元素分布特征復雜化。那么究竟已知鈾礦區(qū)航放鈾含量與微地形存在怎樣的關(guān)系就值得深入探究。以二連盆地已知的砂巖型鈾礦為研究對象，采用相關(guān)性分析法對航放鈾含量及微地形的相關(guān)性進行了探索實驗，以期總結(jié)已知鈾礦航放鈾含量與微地形的特征關(guān)系，并據(jù)此對成礦有利區(qū)的劃分提供指導。

2.1 相關(guān)系數(shù)提取方法

2.1.1 航放數(shù)據(jù)來源

本次采用的數(shù)據(jù)來源于核工業(yè)航測遙感中心于1979—2018 年間在二連盆地飛行獲取的航放數(shù)據(jù)，因不同時期使用的測量儀器、比例尺、測線方向等均不盡相同，所以對各測區(qū)的數(shù)據(jù)進行了連片及調(diào)平處理，并形成了全區(qū)的網(wǎng)格文件。

2.1.2 地形數(shù)據(jù)來源

地形數(shù)據(jù)來源于中國科學院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺（http：//www.gscloud.cn），以平臺中的ASTER GDEM V2版數(shù)據(jù)為基礎(chǔ)，其空間分辨率為30 m，并將獲取的地形數(shù)據(jù)轉(zhuǎn)換為網(wǎng)格文件。

2.1.3 相關(guān)系數(shù)提取

首先，生成南北向500 m 線距、50 m 點距的測線文件，并據(jù)此提取了測線上的航放鈾含量及地形高程數(shù)據(jù)；其次，利用自編軟件，逐測線進行航放鈾含量與地形高程相關(guān)系數(shù)計算，將計算結(jié)果進行網(wǎng)格化，從而形成全區(qū)的相關(guān)系數(shù)網(wǎng)格文件，進而提取已知鈾礦位置的相關(guān)系數(shù)。因使用滑動窗口法求取每條測線的相關(guān)系數(shù)，所以滑動窗口的寬度選擇對計算結(jié)果有較大的影響，本次分別選取了3 km、5 km、7 km、10 km 窗口寬度的相關(guān)系數(shù)進行實驗，用以分析最佳的滑動窗口寬度。

2.2 二連盆地典型礦床分析

2.2.1 腦木根地區(qū)

腦木根礦床及137、701 礦化點均處于高海拔地區(qū)向低海拔地區(qū)過渡的斜坡帶上，其中腦木根礦床略接近于坡底（圖2A）。從航放鈾含量剖面上可以看出（圖2B），腦木根礦床處于平緩背景場向高場變化的梯度帶上，背景場的鈾含量在1.8 ×10-6左右，坡底位置對應(yīng)的鈾含量最高，可達5.7×10-6；137 及701 礦化點處于高場向背景場的過渡地段，并有微弱的局部異常顯示，鈾含量分別為2.24×10-6、2.39×10-6。在相關(guān)系數(shù)曲線上，當滑動窗口寬度為3 km時，曲線形態(tài)較為雜亂，規(guī)律性差；滑動窗口在5～10 km曲線上均表現(xiàn)為負相關(guān)關(guān)系，尤以7 km及10 km窗口曲線對區(qū)域性相關(guān)性反映的最為明顯，相關(guān)系數(shù)介于-0.8～-0.5之間（圖2C、D、E、F）。

圖2 腦木根地區(qū)航放鈾含量與微地形相關(guān)系數(shù)綜合剖面圖Fig.2 Comprehensive sections of correlation coefficient between airborne uranium content and micro-topography in Naomugen area

2.2.2 努和廷地區(qū)

努和廷礦床及217、240 礦化點同樣處于高海拔地區(qū)向低海拔地區(qū)過渡的斜坡帶上，其中努和廷礦床略接近于坡底（圖3A）。在航放鈾含量剖面上（圖3B），努和廷礦床及217、240 礦化點均表現(xiàn)為低背景場中的局部弱異常，其鈾含量分別為1.49×10-6、1.86×10-6、1.75×10-6，異常特征不顯著。在滑動窗口3 km 及5 km 相關(guān)曲線上努和廷礦床從正相關(guān)關(guān)系向不相關(guān)關(guān)系變化，217 及240礦化點則表現(xiàn)為明顯的負相關(guān)關(guān)系；在5～10 km滑動窗口相關(guān)曲線上，努和廷礦床則表現(xiàn)為負相關(guān)關(guān)系，其中滑動窗口10 km曲線將217礦化點的相關(guān)系數(shù)由負相關(guān)轉(zhuǎn)變?yōu)榱苏嚓P(guān)，這是由于滑動窗口寬度的變大對局部曲線形態(tài)造成的影響；尤以7 km 窗口曲線對區(qū)域性相關(guān)關(guān)系反映最為明顯，其相關(guān)系數(shù)介于-0.6～-0.3（圖3C、D、E、F）。

圖3 努和廷地區(qū)航放鈾含量與微地形相關(guān)系數(shù)綜合剖面圖Fig.3 Comprehensive sections of correlation coefficient between airborne uranium content and micro-topography in Nuheting area

綜上分析，二連盆地典型鈾礦床航放鈾含量與微地形的相關(guān)系數(shù)主要介于負相關(guān)-不相關(guān)之間，即已知鈾礦往往分布于局部地形低、航放鈾含量高的地區(qū)或局部地形高、航放鈾含量低的地區(qū)，說明這些地段分布較多的活性鈾或說明鈾表生遷移能力強。3 km 及5 km 窗口曲線較為凌亂，而7 km 及10 km 對區(qū)域性相關(guān)關(guān)系反映最為顯著，由此認為7 km 或10 km 滑動窗口圈定的負相關(guān)區(qū)域可以作為成礦有利區(qū)。

利用10 km 滑動窗口編制了二連盆地航放鈾含量與微地形相關(guān)系數(shù)影像圖（圖4），從圖中可以看出，絕大多數(shù)鈾礦處于航放鈾含量與微地形的負相關(guān)區(qū)內(nèi)，且負相關(guān)區(qū)的分布區(qū)域與主要的含礦砂體分布區(qū)域具有較強的關(guān)聯(lián)性。

圖4 二連盆地航放鈾含量與微地形相關(guān)系數(shù)影像圖Fig.4 Image map of the correlation coefficient between airborne uranium content and micro-topography in Erlian basin

3 航放鈾含量與微地形負相關(guān)區(qū)形成機理分析

實踐證明，雖然砂巖型鈾礦床是隱伏盲礦，其上覆蓋較厚（特別是有致密的隔水泥巖），但當其埋深較小時（一般小于500 m），仍然能在地表或空中觀測到弱放射性異常顯示，即高精度航放測量能探測到這種微弱變化［9］。航放測量之所以能夠探測到微弱的異常，是由于航放主要探測鈾與氡的共同衰變子體214Bi的γ 射線，其母體的運移機制適用于超微量金屬的地氣運移機理和氡的接力遷移理論［10-15］（圖5）。通過地氣攜帶的“納米級”鈾等微粒在地表局部形成富集，從而被航放測量探測到，而富集在地表的微粒容易受到機械搬運及化學搬運的影響，由高程相對較高地區(qū)向較低地區(qū)進行橫向運移，這也正是部分航放異常區(qū)域與鈾礦體在地表投影位置有一定偏移的原因。

圖5 砂巖型鈾礦與航放鈾含量、微地形關(guān)系模式圖（據(jù)參考文獻［15］修改）Fig.5 The relationship between sandstone-type uranium deposits and airborne uranium content and micro-topography（modified after reference［15］）

通過對砂巖型鈾礦與航放鈾含量、微地形的綜合分析可以看出，他們之間存在三種關(guān)系（圖4）。一是鈾礦處于斜坡帶上，由超微量金屬的地氣運移機理可知，在地表富集的鈾等放射性微粒可形成航放異常，但同時在受到表生作用的影響下，放射性微粒會隨地形變化而遷移，在航放異常曲線上則表現(xiàn)為非對稱形態(tài)，即地形高的一側(cè)較陡，地形低的一側(cè)較緩，而在相關(guān)系數(shù)曲線上會展現(xiàn)出較弱的負相關(guān)關(guān)系；二是礦床處于坡頂位置，地表雖有放射性微粒富集，但由于表生作用強烈，放射性微粒大多由地形高的區(qū)域向地形低的區(qū)域進行了遷移，此時在航放鈾含量曲線上會出現(xiàn)略低于其背景場的特征，在等值線平面圖上則處于高場或偏高場的邊部，在相關(guān)系數(shù)曲線上表現(xiàn)為負相關(guān)關(guān)系；三是礦體處于坡底，此種情況下上升至地表的放射性微粒直接在坡底富集，并且還會接受周邊鈾元素的遷移富集，此時在航放鈾含量曲線上會有明顯的異常顯示，在相關(guān)系數(shù)曲線上則表現(xiàn)為明顯的負相關(guān)關(guān)系。

4 負相關(guān)在其他盆地砂巖型鈾礦的驗證

為了驗證其他盆地的砂巖型鈾礦是否也處于航放鈾含量與微地形的負相關(guān)帶內(nèi)，選取了我國北方的鄂爾多斯盆地及松遼盆地內(nèi)的典型砂巖型鈾礦進行了驗證分析。

4.1 鄂爾多斯盆地大營-柴登壕地區(qū)

大營-柴登壕地區(qū)內(nèi)的礦床及礦化點所在位置的地形較其周邊變化不大，約有50 m 左右的高差（圖6A）。在航放鈾含量剖面上（圖6B），72 礦化點有微弱的異常顯示，鈾含量為2.79×10-6，大營、柴登壕礦床及5029 礦化點均處于背景場中，無明顯異常顯示，背景場的鈾含量在2.0×10-6左右。從相關(guān)系數(shù)曲線上可以看出，滑動窗口為3 km 時，72 礦化點及大營礦床表現(xiàn)為負相關(guān)關(guān)系，柴登壕礦床及5029 礦化點表現(xiàn)為不相關(guān)（相關(guān)系數(shù)接近于0），在5～10 km 窗口曲線上，除柴登壕礦床外，其他礦床、礦化點均表現(xiàn)為負相關(guān)關(guān)系，柴登壕礦床則表現(xiàn)為不相關(guān)（圖6C、D、E、F）。

圖6 大營-柴登壕地區(qū)航放鈾含量與微地形相關(guān)系數(shù)綜合剖面圖Fig.6 Comprehensive sections of correlation coefficient between airborne uranium content and micro-topography in Daying-Chaidenghao area

4.2 松遼盆地錢家店-大林地區(qū)

從錢家店-大林地區(qū)航放鈾含量與微地形相關(guān)關(guān)系綜合剖面圖可知，錢2 及錢4 礦床所在位置的地形較為平緩，大林礦床所在位置的地形較其周邊約有30 m 左右的高差（圖7A）。錢2 及錢4 礦床處于背景場中、無明顯異常顯示，鈾含量在1.5×10-6左右，大林礦床則有明顯的異常顯示，鈾含量達3.05×10-6（圖7B）。在相關(guān)系數(shù)曲線上，錢2 礦床在3～10 km 窗口曲線上均表現(xiàn)為負相關(guān)關(guān)系，錢4 及大林礦床在3 km 及5 km 窗口曲線表現(xiàn)為不相關(guān)，在7 km 及10 km 窗口曲線上則表現(xiàn)為負相關(guān)關(guān)系（圖7C、D、E、F）。

圖7 錢家店-大林地區(qū)航放鈾含量與微地形相關(guān)系數(shù)綜合剖面圖Fig.7 Comprehensive sections of correlation coefficient between airborne uranium content and micro-topography in Qianjiadian-Dalin area

從鄂爾多斯盆地及松遼盆地典型砂巖型鈾礦航放鈾含量與微地形相關(guān)分析結(jié)果可以看出，已知砂巖型鈾礦區(qū)域同樣出現(xiàn)了航放鈾含量與微地形的負相關(guān)區(qū)，說明了負相關(guān)區(qū)對砂巖型鈾礦的賦存具有一定的指示意義，這為砂巖型鈾礦成礦預測提供了十分有參考價值的信息。

5 已知鈾礦航放鈾含量與微地形相關(guān)性統(tǒng)計分析

通過提取二連、鄂爾多斯及松遼盆地砂巖型鈾礦床、礦（化）點、異常點的相關(guān)系數(shù)進行統(tǒng)計分析，進一步說明利用相關(guān)系數(shù)劃分成礦有利區(qū)的有效性。

本次使用10 km 窗口進行微地形與鈾含量的相關(guān)關(guān)系計算，據(jù)此提取了二連、鄂爾多斯及松遼盆地170 個砂巖型鈾礦床、礦（化）點、異常點的相關(guān)系數(shù)。經(jīng)統(tǒng)計，鈾礦床、礦點、礦化點及異常點的相關(guān)系數(shù)平均值分別為-0.22、-0.20、-0.20及-0.17（表1）。

表1 鈾礦航放鈾含量與微地形相關(guān)系數(shù)統(tǒng)計表Table 1 Correlation coefficient between airborne uranium content of uranium mineralization and micro-topography

從相關(guān)系數(shù)分布直方圖中可以看出（圖8），相關(guān)系數(shù)-1～0 之間的鈾礦床有12 個，占總數(shù)的80 %；相關(guān)系數(shù)-1～0 之間的鈾礦點有11個，占總數(shù)的73.30 %；相關(guān)系數(shù)-1～0 之間的鈾礦化點有28 個，占總數(shù)的59.6 %；相關(guān)系數(shù)-1～0 之間的鈾異常點有62 個，占總數(shù)的66.7 %；同時，處于負相關(guān)區(qū)內(nèi)的礦床、礦（化）點等共有113 個，占總數(shù)的66.5%。據(jù)此說明絕大多數(shù)鈾礦處于航放鈾含量與微地形的負相關(guān)區(qū)域內(nèi)，表明利用負相關(guān)區(qū)域劃分成礦有利區(qū)具有一定有效性。負相關(guān)的鈾礦床及礦點所占比例顯著高于負相關(guān)的鈾礦化點及異常點所占比例，推斷為礦化點及異常點的規(guī)模較小、鈾遷移量少，導致其在航放鈾含量曲線上反映不明顯所致。

圖8 相關(guān)系數(shù)分布直方圖Fig.8 Histogram of correlation coefficient

6 結(jié)論

本論文僅是針對砂巖型鈾礦空間分布與航放鈾含量、微地形之間關(guān)系的一次探索性研究工作。研究結(jié)果表明，砂巖型鈾礦與航放鈾含量、微地形的負相關(guān)區(qū)具有空間關(guān)聯(lián)，對尋找砂巖型鈾礦具有一定的指示意義；相關(guān)分析最佳的滑動窗口寬度為7 km 或10 km；當前利用航放數(shù)據(jù)進行砂巖型鈾礦成礦預測時主要依靠航放弱信息提取技術(shù)，相關(guān)分析可作為一種補充手段與弱信息提取技術(shù)結(jié)合預測。