2018/12/29 作者/EWG1990儀器學(xué)習(xí)網(wǎng)

本節(jié)涉及氣相色譜儀的幾種常用檢測器(TCD、FID、ECD、FPD、NPD等)的數(shù)據(jù)采集和數(shù)據(jù)分析，這類信號(hào)不包含質(zhì)譜、光譜類的輔助信息，其關(guān)鍵在于如何提高對(duì)峰的檢測精度和準(zhǔn)確度。

氣相色譜儀檢測器輸出的信號(hào)非常快速(可以視為連續(xù)信號(hào))，信號(hào)強(qiáng)度非常低(小至10+A)，同時(shí)信號(hào)是模擬電信號(hào)，因此色譜儀輸出電信號(hào)無法用簡單的方法進(jìn)行定性、定量處理，首先要把模擬電信號(hào)用記錄儀記錄下來或把模擬信號(hào)轉(zhuǎn)換成數(shù)字信號(hào)儲(chǔ)存下來，然后根據(jù)不同分析要求再做處理，以獲得有關(guān)被分析組分的定性、定量結(jié)果或其他信息。氣相色譜的數(shù)據(jù)記錄處理方式一般可分為三種:臺(tái)式記錄儀、數(shù)據(jù)處理機(jī)和色譜工作站也是數(shù)據(jù)處理的不同發(fā)展階段。

早期色譜數(shù)據(jù)的記錄使用臺(tái)式記錄儀，采用差式電位計(jì)、同步電機(jī)、記錄紙對(duì)模擬信號(hào)進(jìn)行顯示和記錄，然后通過手動(dòng)進(jìn)行保留時(shí)間和峰面積的計(jì)算，由于操作煩瑣，目前已完全被淘汰了。

色譜數(shù)據(jù)處理機(jī)是一種專門用于色譜數(shù)據(jù)分析的專用計(jì)算機(jī)，色譜處理機(jī)采用A/D轉(zhuǎn)換器(模數(shù)轉(zhuǎn)換器)將接收到的氣相色譜檢測器傳過來的信號(hào)轉(zhuǎn)換成數(shù)字信號(hào)存儲(chǔ)到專用計(jì)算機(jī)中，分析人員可以通過鍵盤輸入專用計(jì)算機(jī)指令來確定各種參數(shù)進(jìn)行各種定量計(jì)算，進(jìn)而打印出分析報(bào)告，但色譜譜圖存儲(chǔ)數(shù)量有限，結(jié)果也不是非常直觀，制約了色譜數(shù)據(jù)處理機(jī)的推廣和應(yīng)用。

隨著計(jì)算機(jī)技術(shù)的發(fā)展，色譜工作站已成為色譜數(shù)據(jù)處理的主流，其在組成和工作原理與數(shù)據(jù)處理機(jī)基本相同，但不同的是它是計(jì)算機(jī)程序，可在計(jì)算機(jī)操作系統(tǒng)下工作，存儲(chǔ)和處理能力比數(shù)據(jù)處理機(jī)大得多，同時(shí)作為廠家專門配套的工作站，還可以實(shí)現(xiàn)氣相色譜運(yùn)行情況同步顯示，并且具備氣相色譜參數(shù)的實(shí)時(shí)控制。

本節(jié)將著重介紹濾噪和色譜峰檢測等方法。

一、噪聲平滑處理

色譜數(shù)據(jù)處理系統(tǒng)接收到的數(shù)字信號(hào)(經(jīng)過模數(shù)轉(zhuǎn)換的檢測器信號(hào))，包含兩部分信息部分是檢測器對(duì)色譜柱流出物響應(yīng)產(chǎn)生的電信號(hào)，另一部分是色譜系統(tǒng)(包括進(jìn)樣系統(tǒng)、分離系統(tǒng)、檢測系統(tǒng)、放大系統(tǒng)和模數(shù)轉(zhuǎn)換系統(tǒng)等)產(chǎn)生的一切與色譜流出物無關(guān)的信號(hào)，通常這類信號(hào)被稱為“系統(tǒng)噪聲”。系統(tǒng)噪聲疊加在樣品有關(guān)的真實(shí)信號(hào)之上，使色譜數(shù)據(jù)處理系統(tǒng)對(duì)色譜峰的檢測判別以及進(jìn)一步的數(shù)據(jù)處理帶來了不利因素，所以有必要盡可能地將信號(hào)中含的噪聲部分排除掉，以便更好地進(jìn)行譜圖數(shù)據(jù)處理。

更多的儀器學(xué)習(xí)資料，在儀器學(xué)習(xí)網(wǎng)

應(yīng)該注意的是譜圖中的系統(tǒng)噪聲與數(shù)據(jù)采集頻率密切相關(guān)，同一條件下，噪聲信號(hào)的大小正比于采集速率的1/2次方，隨著采集頻率的増加，噪聲的大小也隨之增加。因此適當(dāng)降低數(shù)據(jù)采集頻率可以降低原始譜圖的噪聲信號(hào)，為保證峰面積積分準(zhǔn)確，每個(gè)完整色譜峰信號(hào)采集點(diǎn)數(shù)不應(yīng)小于20，否則可能導(dǎo)致積分結(jié)果的偏差

系統(tǒng)噪聲通常始終有高頻的、雜亂的信號(hào)，去除系統(tǒng)噪聲的方法一般為濾波或平滑。信號(hào)濾波技術(shù)可分為兩大類，一類是硬件濾波技術(shù)，如RC電路濾波，另一類為軟件濾波技術(shù)也稱為數(shù)字濾波。硬件濾波器存在動(dòng)態(tài)范圍小、濾波效果不好、適應(yīng)性差等問題。軟件濾波是用計(jì)算機(jī)程序?qū)?shù)字信號(hào)進(jìn)行平滑、濾波處理，通過預(yù)定的規(guī)則進(jìn)行運(yùn)算或變換，從而達(dá)到去除噪聲的效果。軟件濾波可以根據(jù)不同的實(shí)際情況，改變算法中的各種參數(shù)，來達(dá)到最佳的濾波效果。

軟件數(shù)字濾波的應(yīng)用廣泛，使用方便，算法容易實(shí)現(xiàn)，越來越得到廣泛的應(yīng)用和重視。軟件數(shù)字濾波方法常用的有均值濾波器、 Savitzky-Golay濾波器。

1.均值濾波器( mean filter)

均值濾波器也叫移動(dòng)平均法，首先確定一個(gè)固定的時(shí)間窗口(n)，窗口內(nèi)數(shù)據(jù)(Y1)的平均值替換窗口中心的原始數(shù)據(jù)點(diǎn)，達(dá)到數(shù)據(jù)平滑的作用：

均值濾波器可以有效地壓制信號(hào)抖動(dòng)，尤其是脈沖式信號(hào)。通常采用的是三點(diǎn)中心或五點(diǎn)中心法，窗口設(shè)定越大，濾波效果越好，但譜圖的畸變也越大。通常而言窗口設(shè)定應(yīng)小于典型峰寬1/10，在此條件下，峰的畸變不是非常大，當(dāng)窗口大于典型峰寬15，峰高會(huì)被明顯變低，畸變嚴(yán)重，可能直接影響峰檢測和積分。采用四點(diǎn)中心對(duì)稱的移動(dòng)平均法將會(huì)產(chǎn)生個(gè)“窗口”中心原先沒有的數(shù)據(jù)點(diǎn)，可以提高峰谷點(diǎn)的判斷精度。

2. Savitzky- Golay濾波器

Savitzky- Golay濾波器是1964年由 Savitzky和 Golay］提出的，根據(jù)其原理可稱為移動(dòng)窗口最小二乘多項(xiàng)式平滑，也稱為最小二乘平滑，是分析化學(xué)中應(yīng)用最廣泛的平滑算法。Savitzky- Golay濾波或稱 Savitzky- Golay平滑，利用中心點(diǎn)以及其前后m-1個(gè)數(shù)據(jù)點(diǎn)進(jìn)行多項(xiàng)式的最小二乘擬合，定出多項(xiàng)式中的各項(xiàng)系數(shù)。然后，用多項(xiàng)式的值來代替實(shí)驗(yàn)值以達(dá)到數(shù)據(jù)平滑的目的。原則上說，選取的多項(xiàng)式次數(shù)越高，數(shù)據(jù)組內(nèi)的點(diǎn)數(shù)越多，所計(jì)算出的數(shù)據(jù)平滑性也就越好，但在應(yīng)用中還需將此法對(duì)實(shí)際峰形畸變的影響因素加以考慮。應(yīng)用五點(diǎn)二次法平滑和七點(diǎn)三次法平滑處理較為理想。

五點(diǎn)二次法平滑:yi=(-3xi-2+12xi+17xi+12xi+1-3xi+2)/35

七點(diǎn)三次法平滑:yi=[7xi+6(xi+1+xi-1)+3(xi+2+xi-2)-2(xi+3+xi-3)]/21

其中，x是原始數(shù)據(jù)，y是經(jīng)過平滑處理后的數(shù)據(jù)。從本質(zhì)上說，移動(dòng)窗口多項(xiàng)式平滑就是利用窗口內(nèi)部各個(gè)點(diǎn)之間的加權(quán)來計(jì)算平滑后的新值。

在平滑濾波過程中，應(yīng)該注意的是:對(duì)信號(hào)濾波，或多或少地會(huì)對(duì)有效信號(hào)產(chǎn)生影響，一般來說，隨著信噪比的增加，峰形的畸變也會(huì)相應(yīng)地增大，因此在選用各種濾波方法進(jìn)行濾波時(shí)，應(yīng)遵守這樣一個(gè)原則:在峰形畸變保持在一定允許范圍內(nèi)，盡可能地提高信噪比。

色譜譜圖的濾波一直是科學(xué)研究的熱點(diǎn)問題，一些新的算法也被用于色請(qǐng)語圖的平滑，其中快速傳里葉變換和小波變換被認(rèn)為是理想的濾波算法，無論是傅里葉變換還是小波變換，其實(shí)質(zhì)都是一樣的，就是將信號(hào)在時(shí)間域和頻率域之間相互轉(zhuǎn)換，從看似復(fù)余的時(shí)間域的數(shù)據(jù)進(jìn)行變換，將其轉(zhuǎn)化分解成不同頻率域的信號(hào)疊加，由于信號(hào)往往在頻率域有比在時(shí)間域更加簡單和直觀的特性，根據(jù)色譜數(shù)據(jù)處理的經(jīng)驗(yàn)，噪聲的頻率最高，色譜峰次之，基線的頻率最低，將噪聲所在的頻率域的信號(hào)進(jìn)行平滑，然后通過進(jìn)一步反變換到時(shí)間域上，完成對(duì)請(qǐng)圖的平滑，由于小波變換、傳里葉變換運(yùn)算所需的數(shù)據(jù)量和運(yùn)算量較大，在使用時(shí)一般不能直接完成，需要根據(jù)實(shí)際的噪聲情況進(jìn)行必要的調(diào)整，因此沒有真正用于色譜工作站，但其良好的無歧視的濾波能力，尤其是針對(duì)微弱信號(hào)的獲取可具有較理想的效果，使它成為數(shù)據(jù)處理研究的重點(diǎn)。

二、色譜峰的檢測判別方法

峰檢測的目標(biāo)是判斷色譜流出峰的起始點(diǎn)、最高點(diǎn)和結(jié)束點(diǎn)，它的準(zhǔn)確性直接影響后續(xù)的定性、定量計(jì)算。隨著計(jì)算機(jī)技術(shù)的發(fā)展，色譜數(shù)據(jù)處理機(jī)和工作站廣泛應(yīng)用于色譜的數(shù)據(jù)處理，基于軟件的峰檢測自動(dòng)判定方法獲得快速的發(fā)展，已產(chǎn)生了各種不同的色譜峰處理方法，如幅值法、統(tǒng)計(jì)門限法、面積靈敏度法與基線靈敏度法、一階導(dǎo)數(shù)判別法、修正的一階導(dǎo)數(shù)判別法、二階導(dǎo)數(shù)判別法、增加斜率法等。

階導(dǎo)數(shù)法是目前工作站中最為常用的色譜峰檢測方法，主要是依據(jù)色譜峰信號(hào)的斜率變化(一階導(dǎo)數(shù))進(jìn)行峰起點(diǎn)、落點(diǎn)和峰頂點(diǎn)的判別。色譜峰信號(hào)的一階導(dǎo)數(shù)如圖1所示，由圖中可以看出在峰的起點(diǎn)處，其相應(yīng)的一階導(dǎo)數(shù)由零變成正，在峰頂點(diǎn)處，相應(yīng)的一階導(dǎo)數(shù)由正變成負(fù)，在峰結(jié)東處一階導(dǎo)數(shù)由負(fù)變成零。

圖1 色譜峰及其相應(yīng)的一階導(dǎo)數(shù)示意圖

通常情況下采樣頻率固定。為了計(jì)算方便，通常直接采用相鄰信號(hào)的差值作為斜率，由于噪聲的存在，基線處的一階導(dǎo)數(shù)并不為零，而是在零點(diǎn)附近上下波動(dòng)，所以，在一階導(dǎo)數(shù)檢測時(shí)須確定一個(gè)閾值。當(dāng)信號(hào)增加，一階導(dǎo)數(shù)值大于預(yù)定的閾值，確定峰的起點(diǎn)。當(dāng)導(dǎo)數(shù)值由正變負(fù)時(shí)為峰頂點(diǎn)。當(dāng)峰頂點(diǎn)確定后，若信號(hào)的一階導(dǎo)數(shù)值的絕對(duì)值小于閾值時(shí)為峰的結(jié)東點(diǎn)。閾值可以根據(jù)基線中噪聲偏差設(shè)定，通常設(shè)定值為基線噪聲的3倍。

傳統(tǒng)的一階導(dǎo)數(shù)的缺陷在于對(duì)于對(duì)稱峰形檢測較為準(zhǔn)確，但對(duì)于拖尾峰而言，由于值的存在可能導(dǎo)致峰結(jié)束點(diǎn)的判定提前。另一方面，對(duì)于那些峰寬較大的小峰可能會(huì)出現(xiàn)漏檢的現(xiàn)象。為了改善信噪比較低的扁平峰的檢測，1982年J.L. Excoffier等提出了修正一階導(dǎo)數(shù)法用于色譜峰檢測。對(duì)于相同峰高的色譜峰而言，峰越寬，一階導(dǎo)數(shù)的值越小，由于值的限定，可能導(dǎo)致某些峰肉眼可見，但未獲得檢出。

由于噪聲的影響，相鄰信號(hào)差值的波動(dòng)往往會(huì)大于色譜峰前沿的斜率，可能對(duì)判斷造成影響。修正一階導(dǎo)數(shù)法采用相隔幾個(gè)點(diǎn)的差值來代替相鄰兩點(diǎn)的差作為斜率。在此情況下斜率的最大值隨兩點(diǎn)間距的增大而增大，當(dāng)兩點(diǎn)的距離為峰寬的一半時(shí)，數(shù)值達(dá)到最大振幅圖2為同一譜圖的一階導(dǎo)數(shù)和修正一階導(dǎo)數(shù)的結(jié)果，表明修正一階導(dǎo)數(shù)明顯優(yōu)于一階導(dǎo)數(shù)。

圖2 色譜譜圖的一階導(dǎo)數(shù)和修正一階導(dǎo)數(shù)比較

A-原始譜圖；B-一階導(dǎo)數(shù)；C-修正一階導(dǎo)數(shù)

修正一階導(dǎo)數(shù)的峰檢測過程及判別依據(jù)與一階導(dǎo)數(shù)法相同，同樣需要根據(jù)噪聲情況來確定峰判斷的閾值，但從圖9-2可以看出相同的基線噪聲條件下，噪聲對(duì)修正一階導(dǎo)數(shù)干擾明顯小于一階導(dǎo)數(shù)，窗口越大干擾越小，其原理與采集頻率-噪聲相關(guān)性相同。

二階導(dǎo)數(shù)也可用于色譜峰的判別，從實(shí)際應(yīng)用效果而言，用二階導(dǎo)數(shù)法的檢測靈敏度要高于一階導(dǎo)數(shù)法，甚至在基線出現(xiàn)線性漂移的情況下也不受影響。但基線噪聲對(duì)信號(hào)導(dǎo)數(shù)的影響隨導(dǎo)數(shù)階數(shù)的增加而增大，會(huì)對(duì)程序上的峰判斷產(chǎn)生不利影響，在實(shí)際工作站中較少使用。

E. F.G. Woerlee等提出增加斜率檢素法(TS法)，如圖3所示。

圖3 ITS法的起峰判斷

將起始點(diǎn)(B1)固定，與其后連續(xù)四個(gè)信號(hào)點(diǎn)(B2、B3、B4、B3)分別計(jì)算其斜率值。當(dāng)四個(gè)斜率值連續(xù)上升，則確定起始點(diǎn)(B1)為峰起點(diǎn)。這些關(guān)系可用下列公式表示:

將上述關(guān)系式加以整理可轉(zhuǎn)換為:

B3-2B2+B1>0＆＆2B4-3B3+B1>0＆＆3B5-4B4+B1>0

我們可以看出，實(shí)際上這種檢測方法是二階導(dǎo)數(shù)的一種變形。當(dāng)數(shù)據(jù)點(diǎn)滿足二階導(dǎo)數(shù)檢測法的條件時(shí)，也必能滿足ITS法，相反，能滿足ITS法的不一定能滿足二階導(dǎo)數(shù)法。ITS法對(duì)噪聲的敏感性要低于二階導(dǎo)數(shù)法，同時(shí)具有更高的檢測靈敏度。

更多的儀器學(xué)習(xí)資料，在儀器學(xué)習(xí)網(wǎng)

傳統(tǒng)的一維氣相色譜通常的峰容量比較低(通常小于1000現(xiàn)實(shí)樣品中組分?jǐn)?shù)則成千上萬，其分離能力遠(yuǎn)不能滿足使樣本的峰獲得完全的分離，通常的色譜處理系統(tǒng)一般對(duì)肩峰的情況不予考慮，將肩峰并入主峰進(jìn)行計(jì)算。這種處理方法將給定量結(jié)果帶來較大偏差。重疊峰分析也是色譜數(shù)據(jù)處理的研究重點(diǎn)，盧佩章教授等通過EMG模型進(jìn)行全譜圖曲線擬合的方法解決了重疊峰的定量問題?。但擬合法的色譜峰的模型參數(shù)以及色譜參數(shù)隨保留時(shí)間的變化對(duì)擬合結(jié)果至關(guān)重要，而且運(yùn)算復(fù)雜，同時(shí)缺乏質(zhì)譜等其他信息進(jìn)行軸助驗(yàn)證，難以實(shí)現(xiàn)廣泛的使用，在通常的數(shù)據(jù)處理軟件(工作站)中較少使用。然而，從實(shí)際的數(shù)據(jù)出發(fā)采用二階導(dǎo)數(shù)可能實(shí)現(xiàn)重疊峰尤其是肩峰的拆分。圖9-4為一后肩峰譜圖及其一階導(dǎo)數(shù)和二階導(dǎo)數(shù)，理論上，后肩峰可根據(jù)以下條件進(jìn)行判斷:有后肩峰時(shí)，從主峰頂點(diǎn)到峰結(jié)東點(diǎn)，信號(hào)的一階導(dǎo)數(shù)保持正號(hào)；而二階導(dǎo)數(shù)出現(xiàn)了正負(fù)變換，由此可判斷肩峰的存在。

圖4 后肩峰的色譜圖及其一階導(dǎo)數(shù)和二階導(dǎo)數(shù)

但如前文所述，信號(hào)的二階導(dǎo)數(shù)對(duì)系統(tǒng)噪聲十分敏感，直接采用上述方式進(jìn)行判斷容易出現(xiàn)錯(cuò)誤。為了減小二階導(dǎo)數(shù)的影響，可以采用查找一階導(dǎo)數(shù)極值的方法進(jìn)行肩峰判斷，從圖9.4中可以看出存在三個(gè)拐點(diǎn)，即一階導(dǎo)數(shù)出現(xiàn)三個(gè)符號(hào)為負(fù)的極值點(diǎn)，兩個(gè)極小、一個(gè)極大。為避免由于噪聲的存在對(duì)肩峰判斷的影響，當(dāng)三個(gè)極值點(diǎn)找到后，極大值和極小值之間的差值需大于一個(gè)預(yù)先給定的門限值。此門限值根據(jù)噪聲大小來確定，同時(shí)相鄰極值點(diǎn)間的時(shí)間間隔也需滿足一定條件，只有滿足以上的系列條件才能確定出肩峰。

三、峰基線的確定

峰高是峰尖至峰底(或基線)的距離，峰面積是色譜峰與峰底(或基線)所圍成的面積。峰高和峰面積是色譜定量的基本參數(shù)。前面論述了色譜峰的起點(diǎn)、終點(diǎn)的判定方法，只有色諧峰曲線與基線間的整個(gè)面積才表示該組分的總量，故峰面積的測量同時(shí)還取決于基線的判斷。因此，要準(zhǔn)確地測量峰高和峰面積，另一關(guān)鍵是基線的確定。

圖9-5表示了不同情況下確定色譜基線的四種方法。

(1)水平基線校正法見圖5(a)，將始點(diǎn)水平畫至最后一個(gè)色譜峰的始點(diǎn)，折上與該峰的終點(diǎn)相連的線，作為色譜的基線。此法用于基線基本平穩(wěn)的情況。

(2)連帶法或膠帶法見圖5(b)，由第一個(gè)色譜峰的始點(diǎn)與最后一個(gè)色譜峰的終點(diǎn)連一直線，色譜基線則由低于此直線的各峰谷相連而成。此法適用于基線波動(dòng)不大的情況(3)峰谷法見圖5(c)，把峰谷點(diǎn)之間的連線作為基線。此法適用于基線波動(dòng)較大的情況。

(4)空白基線法見圖5(d)，是在不進(jìn)樣的情況下先作一次色譜分析，其譜圖稱為空白基線。然后在完全相同的色譜條件下進(jìn)樣分析，此時(shí)譜圖與空白基線譜圖相減，稱為空白基線法。此方法主要用于程序升溫等基線漂移較為嚴(yán)重的情況，對(duì)儀器的重復(fù)性要求很嚴(yán)否則會(huì)產(chǎn)生很大誤差。

圖5 確定色譜基線的幾種方法

譜圖基線必須根據(jù)色譜圖的分離情況進(jìn)行綜合考慮。在使用微處理機(jī)和計(jì)算機(jī)時(shí)可以將色譜圖調(diào)出，看譜圖上基線的判斷、峰的起落點(diǎn)和未分離峰的切割是否合理，若不合理，就必須重新設(shè)定參數(shù)，或通過手動(dòng)的方式，再次處理數(shù)據(jù)直至滿意為止。

兩峰重合時(shí)，在數(shù)據(jù)處理中通常采用谷點(diǎn)作垂直線的方法進(jìn)行切分。當(dāng)谷點(diǎn)小于小峰的半峰高或者兩峰等高時(shí)，誤差不大。應(yīng)當(dāng)注意的是當(dāng)谷點(diǎn)在小峰的半峰高以上時(shí)，小峰面積的測量誤差會(huì)隨大峰和小峰的峰高比值的增加而顯著增加。Proksch等建立了氣相色譜重合峰峰面積的校正系數(shù)表。通過重合峰的峰高比、分率等信息可在表中査到校正系數(shù)，用于面積折算。但其只有在峰形符合高斯分布、兩個(gè)峰之間存在峰谷、小峰在大峰之前流出等條件下才適用。

更多的儀器學(xué)習(xí)資料，在儀器學(xué)習(xí)網(wǎng)

免責(zé)聲明：文章僅供學(xué)習(xí)和交流，如涉及作品版權(quán)問題需要我方刪除，請(qǐng)聯(lián)系我們，我們會(huì)在第一時(shí)間進(jìn)行處理。

---

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的怎么进行数据平滑滤波_气相色谱数据处理方法 EWG1990仪器学习网的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。