?鉛鋅礦石稀散元素鎘含量的高光譜反演研究

引言

稀散元素即稀有分散元素，又稱分散元素，是指在地殼中豐度低（一般為１０－９級(jí)）、多分散賦存的元素，常形成伴生、共生礦產(chǎn)資源，很少形成獨(dú)立的礦床類型，主要包括鎵（Ｇａ）、鍺（Ｇｅ）、鎘（Ｃｄ）、銦（Ｉｎ）、錸（Ｒｅ）、鉈（Ｔｌ）等８種元素。稀散元素因其獨(dú)特、優(yōu)異的物理、化學(xué)性能而被廣泛應(yīng)用于國民經(jīng)濟(jì)各行業(yè)，是全球高新技術(shù)產(chǎn)業(yè)中不可或缺的關(guān)鍵性礦物原料，是世界各主要經(jīng)濟(jì)體重點(diǎn)勘探的戰(zhàn)略性礦產(chǎn)資源。稀散元素具有超常富集的特征和成礦專屬性（包括礦床類型的專屬性和賦存礦物的專屬性。

本文以西藏某鋅礦集區(qū)的鉛鋅礦石為測(cè)試樣品， 利用地物光譜儀及相關(guān)軟件進(jìn)行鎘（Ｃｄ）元素的光譜采集與預(yù)處理，分析鎘（Ｃｄ）元素的光譜特征與規(guī)律，篩選特征波段，開展隨機(jī)森林（ＲＦ）、人工神經(jīng)網(wǎng)絡(luò)（ＡＮＮ）、支持向量機(jī) （ＳＶＭ）模型的機(jī)器學(xué)習(xí)與反演，評(píng)價(jià)基于高光譜的鎘（Ｃｄ）元素含量反演模型，為鎘（Ｃｄ）等稀散金屬元素含量反演、快速無損檢測(cè)提供參考，為研究區(qū)等高海拔區(qū)域的鎘（Ｃｄ）等稀散金屬礦產(chǎn)資源勘探、評(píng)價(jià)提供思路，為進(jìn)一步的稀散金屬礦產(chǎn)資源高光譜快速掃面勘探研究提供學(xué)術(shù)支持。

實(shí)驗(yàn)部分

2.1 樣品處理

鎘（Ｃｄ）的原子序數(shù)是４８，位于元素周期表第五周期ⅡＢ族，單質(zhì)為銀白色金屬。鎘（Ｃｄ）是典型的稀散元素，具有親硫化和親石性，多伴生于中低溫的鉛鋅礦床中，鉛鋅礦石（方鉛礦－閃鋅礦礦石）是鎘（Ｃｄ）的主要載體，是研究鎘（Ｃｄ）元素地球化學(xué)性質(zhì)及光譜特征的主要對(duì)象。為了保證樣品的有效性和代表性，野外均勻選取礦集區(qū)不同類型的鉛鋅礦石作為實(shí)驗(yàn)樣品，帶回實(shí)驗(yàn)室進(jìn)行挑選、研磨、過篩等預(yù)處理，控制樣品粒徑小于７５μｍ。每件試樣縮分成２份，１份用于化學(xué)分析，１份留作副樣并用于光譜測(cè)試。

2.2 光譜測(cè)試

為避免外界光線等因素的干擾，在夜間封閉的實(shí)驗(yàn)環(huán)境中測(cè)試。測(cè)試前，對(duì)地物光譜儀進(jìn)行白板校正，測(cè)試過程中每隔１０ｍｉｎ左右進(jìn)行一次白板校正。為消除隨機(jī)噪聲的干擾，將光譜試樣置于專用光譜測(cè)試樣品杯中混勻整平，利用專用便攜式地物光譜儀探頭固定裝置保持地物光譜儀探頭垂直試樣平面；每個(gè)試樣分別測(cè)試５次，每次采集后，微調(diào)探頭固定裝置，以便繼續(xù)測(cè)試同一樣品的不同平面微區(qū)。

圖1 樣品光譜采集

結(jié)果與結(jié)論

3.1 實(shí)測(cè)鎘元素含量統(tǒng)計(jì)分析

化學(xué)分析樣品３２件，樣品實(shí)測(cè)的鎘元素含量最大值為３３.９００×１０－６，最小值為１０．１００×１０－６，平均值為１７．９７５×１０－６，樣品變異系數(shù)為４０．３１０％（表１）。對(duì)比鎘元素在地幔、 地殼及各類巖石中的豐度值，實(shí)測(cè)鉛鋅礦石樣品的鎘元素含量異常特征明顯，具有較強(qiáng)的空間異質(zhì)性，適合開展元素含量的光譜反演研究。

表１礦石實(shí)測(cè)鎘元素含量統(tǒng)計(jì)分析

3.2 光譜數(shù)據(jù)預(yù)處理

光譜數(shù)據(jù)預(yù)處理是光譜分析、特征波段選擇、模型構(gòu)建、反演預(yù)測(cè)等的基礎(chǔ)，通過光譜數(shù)據(jù)的預(yù)處理，可有效降低光譜數(shù)據(jù)的噪聲及冗余信息，增強(qiáng)光譜數(shù)據(jù)的有用信息及特征波段。光譜測(cè)試結(jié)束后，利用地物光譜儀配套的光譜數(shù)據(jù)處理軟件進(jìn)行光譜數(shù)據(jù)的均值預(yù)處理和可視化分析，將去除異常值后的光譜數(shù)據(jù)作為鎘元素的實(shí)測(cè)光譜數(shù) 據(jù)，利用該數(shù)據(jù)生成鎘元素的原始光譜曲線（圖2），進(jìn)行光譜數(shù)據(jù)質(zhì)量評(píng)價(jià)和光譜特征分析。

圖2 鉛鋅礦石的鎘元素原始光譜曲線

總體上，實(shí)測(cè)礦石樣品中稀散金屬鎘元素的光譜反射率較高，反射率分布相對(duì)集中，波譜異常明顯，波譜特征穩(wěn)定。在反射率特征方面，所測(cè)樣品中有２９件樣品的反射率集中分布于４０％～６０％區(qū)間，２件樣品的反射率主要分布于６０％～７０％區(qū)間，１件樣品的反射率主要分布于３０％～４０％區(qū)間。在波段特征方面，３５０～７００ｎｍ波段區(qū)間為樣品反射率上升區(qū)間，但波譜上升斜率有所差異。其中，１件樣品的反射率迅速增高后于４００ｎｍ附近開始出現(xiàn)拋物線式的下降特征；４件樣品的反射率上升后于 ４５０ｎｍ附近出現(xiàn)平緩的變化趨勢(shì)；８件樣品的反射率上升趨勢(shì)較緩，于７５０ｎｍ左右出現(xiàn)寬緩波峰后微弱下降；１９件樣品的反射率上升斜率較大，于５００ｎｍ左右形成不對(duì)稱波峰。大部分樣品在８００～８５０ｎｍ區(qū)間出現(xiàn)了一個(gè)寬緩的反射率波峰曲線，并在１０００ｎｍ附近出現(xiàn)小幅度的反射率陡降。在１４２０，１９２０，２２００，２３５０和２４５０ｎｍ附近均形成明顯的反射率波谷曲線（吸 收 峰），尤其以１４２０，１９２０和２２００ｎｍ處的尖棱狀波谷最為明顯。在２２９０～２５００ｎｍ區(qū)間所有樣品的反射率均為右傾下降趨勢(shì)。

3.3 鎘元素特征波段篩選

特征波段選擇是在波段預(yù)處理的基礎(chǔ)上，通過提取原始波譜數(shù)據(jù)中信息量大、相關(guān)性小的波段子集來降低光譜數(shù)據(jù)維度，提高模型反演的精度和穩(wěn)定性。利用軟件中的功能對(duì)鎘元素原始光譜曲線進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析，并選取一階微分（ＦＤ）、二階微分（ＳＤ）、倒數(shù)的對(duì)數(shù)（ＡＴ）、倒數(shù)對(duì)數(shù)的一階微分 （ＡＦＤ）四種數(shù)據(jù)變換方法進(jìn)行光譜數(shù)據(jù)的處理，突出光譜曲線的吸收特征，以便更加直觀和定量地開展鎘元素光譜曲線的變化規(guī)律分析，總結(jié)光譜數(shù)據(jù)特征，進(jìn)行鎘元素特征 波段篩選。基于相關(guān)性、信息量和特征波段的準(zhǔn)則，利用軟件，結(jié)合皮爾森相關(guān)性系數(shù)分析進(jìn)行特征波段的選擇，篩選出對(duì)建模貢獻(xiàn)大的波段作為建立反演模型的最佳變量。皮爾森相關(guān)性系數(shù)ｒ是反映兩個(gè)變量的相關(guān)性程度，ｒ的取值范圍為［－１，１］，ｒ的絕對(duì)值越接近于１，則變量間的相關(guān)性越強(qiáng)，反之則弱。對(duì)鎘元素原始光譜曲線分別進(jìn)行ＦＤ，ＳＤ，ＡＴ，ＡＦＤ變換和皮爾森相關(guān)性系數(shù)ｒ分析，不同變換方式篩選出的特征波段及相關(guān)性系數(shù)見表２。

表２ 不同光譜變換方式的特征波段及相關(guān)性系數(shù)

圖3 SD變換曲線

分析結(jié)果表明，在各類光譜變換方式中，獲得較大相關(guān)性系數(shù)的波段區(qū)間為３５３～２２８３ｎｍ，在２２９０～２５００ｎｍ波段區(qū)間內(nèi)相關(guān)性系數(shù)很小。在光譜數(shù)據(jù)降維與特征波段篩選變換中，ＳＤ光譜變換處理的效果最好，ＳＤ光譜變換處理共篩選出１５個(gè)特征波段，覆蓋可見光及近紅外波段（圖3）。在０．０１級(jí)別相關(guān)性顯著的特征波段有１０個(gè)（６６５、７０９、７８１、９０７、１１４６、１１８２、１２０２、１５６３、 １９２０和１９８４ｎｍ），主要分布在近紅外波段區(qū)間；在０．０５級(jí)別相關(guān)性顯著的特征波段有５個(gè)（３７９、６０３、７７１、８４８ 和 １６２３ｎｍ）。

結(jié)論

鉛鋅礦石樣品中的稀散元素鎘含量較高，鎘元素在可見光－近紅外波段的光譜特征明顯，波形穩(wěn)定，波譜反射率較高且集中分布于４０％～６０％區(qū)間，在１４２０、１９２０和２２００ ｎｍ 處形成有吸收峰，適合開展稀散元素含量的高光譜反演研究。

推薦：

便攜式地物光譜儀 iSpecField-NIR/WNIR

專門用于野外遙感測(cè)量、土壤環(huán)境、礦物地質(zhì)勘探等領(lǐng)域的最新明星產(chǎn)品，由于其操作靈活、便攜方便、光譜測(cè)試速度快、光譜數(shù)據(jù)準(zhǔn)確是一款真正意義上便攜式地物光譜儀。

審核編輯 黃宇