
激光誘導擊穿光譜(LIBS)技術是一種新興的元素檢測和分析技術。當激光作用于樣品表面時,在極短的時間內誘導產生含有樣品物質信息的等離子體。通過對等離子體發(fā)出的光譜進行收集分析,建立光譜與元素一一對應的關系。由此可以檢測出元素周期表中的大部分元素,檢測限可達ppm級。通過標準物質定標后可以進一步得到元素的含量信息。
本研究采用5組國標土壤樣品的光譜數據,通過對多種預處理方法的排列組合來探討光譜預處理的必要性及預處理方法的選擇。面向復雜基質樣品的LIBS光譜預處理方法及定量精度優(yōu)化策略。

我們通過5組國標土壤樣品的實驗得到元素的發(fā)射光譜。但是實驗時可能因為種種原因導致光譜的質量不是很好。一方面可以調整激光參數、采樣方式、脈沖次數、延遲時間等再次得到質量較好的光譜。有了較好質量的光譜是后期計算得到更好結果的前提。通常能夠提高光譜質量的策略是:能量越大光譜信號越強;脈沖次數多,經過累加后可以得到更強的信號;對于基體效應強的樣品,可以采用慢速畫線法進行采樣,可以有效避免單點多次重復采樣帶來的信號衰減問題;適當提前延遲時間可以有效采集到一些信號較弱的元素和一些等離子體壽命較短的元素信息;實驗時通入氦氣可以有效降噪、通入氬氣可以達到增敏的作用;對于質地疏松的樣品可按1:10的比例加入少量粘合劑進行壓片。

另一方面,我們也需要后期對光譜數據進行一些必要預處理。元素含量的計算首先要進行標準曲線擬合。生態(tài)類樣品,例如土壤、植物、中藥材等屬于基質較為復雜的樣品。它們的光譜信號往往會受到雜散光、噪聲、基線漂移等一些因素的干擾,從而影響最終的定量分析結果。因此通常需要在建模前對原始光譜進行預處理。目前已有的光譜預處理方法有很多種,如何選擇合適的預處理方法就顯得尤為重要了。常見的基線校正方法包括自動縮放、均值中心化、帕累托縮放等多種方法。光譜數據預處理算法包括一階導數、二階導數、基線校正、平滑處理、減法運算、歸一化、多元散射校正(MSC)、標準正態(tài)變量變換(SNV)等多種方法。
基線校正方法中,均值中心化是一種基礎的數據平移方法。它通過減去每個變量的均值,使變換后的變量均值為零。能讓光譜數據的差異以零為基準,便于后續(xù)PCA/PLS分析。但由于各波長信號強度天然差異很大,故很少單獨使用。自動縮放的目的是消除波長間量綱差異,讓所有變量被“平等"對待。經過變換后,每個波長對所有樣本的貢獻在數值上變得同等重要。帕累托縮放是一種介于均值中心化和自動縮放之間的折中方案,使用標準差平方根作為縮放因子。可有效放大弱信號,抑制噪聲。能更穩(wěn)健地捕捉數據的真實變化,也能有效的提升微量元素較弱的信號,幫助模型建立更準確的定量關系。是LIBS數據建模的第一選擇,推薦使用。

下面具體介紹一下光譜數據預處理算法對于LIBS光譜的作用及影響。
| 基線校正 (Baseline Correct) | 多種算法組合,旨在扣除由連續(xù)背景輻射產生的非線性基線。 | 此算法對于LIBS光譜的分析至關重要。LIBS信號常疊加在強連續(xù)背景上,基線校正是定量的必要步驟。 |
| 平滑(Smooth) | 通過移動平均、Savitzky-Golay等濾波器,降低白噪聲等高頻隨機噪聲。 | 此算法常用于LIBS光譜的分析,能有效提高信噪比。 |
| 一階導數 (1st Derivative ) | 計算相鄰波長點的強度差,凸顯光譜的斜率變化特征。 | 此算法常用于LIBS光譜的分析,可有效抑制平移的基線漂移,增強峰形分辨能力,常用于分離重疊峰。 |
| 二階導數 (2nd Derivative) | 對一階導數結果再次求導,凸顯光譜的曲率變化。 | 此算法常用于LIBS光譜的分析,能消除線性基線,提供比一階導更高的分辨率,但對噪聲更敏感。 |
| 歸一化(Normalize) | 將光譜按面積(總強度)、背景或內標線進行縮放。 | 此算法常用于LIBS光譜的分析,用于校正激光脈沖能量波動等導致的整體強度變化。 |
| 多元散射校正(MSC) | 通過線性回歸將各光譜校正至“理想"參考光譜。 | 此算法常用于LIBS光譜的分析,用于修正由樣品顆粒度、表面不平整等物理因素引起的乘性和加性散射效應。 |
| 標準正態(tài)變量變換(SNV) | 每條光譜獨立進行標準化(減去均值除以標準差),使數據服從正態(tài)分布。 | 此算法常用于LIBS光譜的分析,用于修正由樣品顆粒度、表面不平整等物理因素引起的乘性和加性散射效應。但計算不依賴參考光譜,更獨立。 |
| 其他算法,如Align, Divide By, Log10, Multiply, Subtract等 | 對于LIBS光譜作用不是很大 | 可根據實際情況按需使用 |
基于上述的說明,處理復雜基質樣品LIBS光譜有效的方法原則是先粗后細,降噪優(yōu)先,再進行適當的校正。本著精準實施預處理策略的原則。處理物理干擾(顆粒效應、能量波動)等,宜采用SNV 或 MSC的算法。消除復雜基線漂移,宜采用平滑和校正基線。增強譜峰分辨率、分離重疊峰,宜采用一階導數和二階導數。應對激光能量波動,宜采用歸一化處理。在軟件里一次性設置好順序,同步運行。可有效提高數據分析效率及分析準度。處理流程可以是:原始光譜 → 平滑 → 基線校正→ 一階導數 → SNV或MSC歸一化 → 完成模型運算 → 檢驗準度。其中的歸一化處理,如果樣品均勻性較好,SNV效果通常不錯;如果樣品差異大且能找到合理的“參考光譜",MSC可能更優(yōu)。當然用戶也可以根據自己的實際情況和需求自行設定模型算法、執(zhí)行順序和更多細節(jié)的設定。
我們的實驗為5組國標土壤樣品的光譜曲線進行擬合,數據預處理算法分別為:Mean-center, Pareto + Smooth, Pareto + Smooth、Baseline Correct, Pareto + Smooth、Baseline Correct、1st Devivative、2nd Devivative、MSC、SNV。由于每個元素對激光的敏感度不一樣,其在光譜上的表征也會不一樣,因此針對不同的元素實際情況也可以采樣不同的算法。比如Na元素的特征光譜是雙峰,Mean-center的單一算法更適合。處理的原則是根據光譜的實際情況采用適當的算法,過多的算法疊加也會導致標準曲線的過度擬合,達不到預期,比如最后一種算法組合就屬于過度擬合。數據結果中的r、Slope和SEP含義分別是擬合相關系數、斜率和預測標準誤差。其中r的值越接近于1表示預測值與真實值的線性相關性越強,結果也會越好。Slope值大于1,說明模型對高濃度樣品預測偏高,對低濃度樣品預測偏低;如果小于1,則情況相反;越接近于1效果越理想。SEP值越小表示預測越精確。它直接反映了模型在預測新樣本時的平均誤差大小,單位與你的元素含量(如wt%)相同。例如加入Baseline Correct算法后就對降低SEP值有著較大的貢獻,比如Pareto + Smooth、Baseline Correct 算法組合里針對Al元素的預處理效果有很大提升。

由數據表可見,采用Mean-center算法后,除Ca和Mg元素,其他元素擬合情況總體較好;其中Na元素,由于其特征光譜是雙峰,Mean-center的單一算法更適合。采用Pareto + Smooth算法組合,所有元素擬合情況較好,r值均大于0.9。在此基礎上增加Baseline Correct算法,r值進一步提升,Al元素的SEP值大幅減小,說明Baseline Correct算法對降低SEP值有著較大的貢獻。最后一種多算法組合由于算法疊加太多,導致了過度的擬合。
綜上所述,前期得到高質量的光譜需要按如下策略進行采集:較強的激光能量,多次脈沖累加,畫線法采樣,適當提前延遲時間,按需通入緩沖氣體;對于質地疏松的樣品可按比例加入粘合劑進行壓片。后期計算推薦的光譜數據預處理算法組合為Pareto(帕累托縮放) + Smooth(5點平滑)、Baseline Correct(基線校正)。根據實際情況可以適當增加1st Devivative、2nd Devivative、MSC、SNV等預處理算法進行進一步的優(yōu)化。或者針對個別特殊的元素進行特殊的算法優(yōu)化方案。但切不可疊加過多的算法,導致過度的擬合。
實驗裝置采用北京澳作生態(tài)儀器有限公司自主研發(fā)的EcoChem 激光光譜元素分析系統(tǒng),該系統(tǒng)采用激光誘導擊穿光譜(LIBS)技術。激光能量200mJ@1064nm,能量輸出0-100%可調;重現率20Hz;光斑大小50-300μm連續(xù)可調;三維全自動工作臺;樣品室可通入氦氣或氬氣。

通過本實驗的分析研究,我們了解了光譜預處理的原則、各種算法的實際作用、建立了LIBS光譜預處理方法及定量精度優(yōu)化策略有效的方法。激光誘導擊穿光譜技術作為一種新興的、快速的、微損的微區(qū)化學分析技術,有著傳統(tǒng)儀器所不具備的更多功能。隨著這一技術的發(fā)展,它必將在元素分析領域發(fā)揮越來越重要的作用。
聯系我們
北京澳作生態(tài)儀器有限公司 公司地址:北京市海淀區(qū)路3號院6號樓1單元一層101A 技術支持:化工儀器網掃一掃 更多精彩
微信二維碼
網站二維碼