黃福祥，栗大超，伍 鵬，張晶鑫，宋 冰，徐可欣

(天津大學(xué)精密測試技術(shù)及儀器國家重點實驗室，天津 300072)

表面等離子體共振(surface plasmon resonance，SPR)技術(shù)[1-6]是近年發(fā)展起來的先進檢測技術(shù)，具有靈敏度高、無需標記、穩(wěn)定快速、便捷實時等特點[7-8]．筆者根據(jù)表面等離子共振檢測葡萄糖濃度時信號和噪聲的特點，對 SPR檢測技術(shù)的數(shù)據(jù)處理算法進行了研究，采用固定點數(shù)的動態(tài)基線法和加權(quán)質(zhì)心法處理SPR數(shù)據(jù)，提高了SPR系統(tǒng)的檢測分辨率，將改進后的數(shù)據(jù)處理方法應(yīng)用于 SPR葡萄糖濃度檢測系統(tǒng)的信號處理，獲得了很好的實驗效果．

1 SPR檢測葡萄糖時的信號特征

SPR系統(tǒng)葡萄糖濃度檢測原理如下：當入射光在傳感器棱鏡端面發(fā)生全反射時，產(chǎn)生的倏逝波激發(fā)棱鏡底面所鍍金屬膜附近的自由電子，在棱鏡和葡萄糖溶液的界面產(chǎn)生表面等離子共振波．當倏逝波的P分量與表面等離子波發(fā)生諧振時，倏逝波的光能大部分耦合到表面等離子波中，使得棱鏡中全反射光強急劇下降，形成吸收峰．此時全反射的角度稱為共振角，檢測共振角出現(xiàn)的位置就可以推測金屬膜表面葡萄糖溶液的濃度情況．

SPR葡萄糖檢測系統(tǒng)的輸出是由光電接收器接收到的光強度轉(zhuǎn)換而成的電信號，SPR輸出曲線如圖1所示，橫軸為 SPR光電接收器陣列的像元數(shù)，縱軸代表接收器接收到的光強度．當金屬膜表面葡萄糖溶液折射率發(fā)生變化時，輸出 SPR曲線會發(fā)生整體的移動；通過計算SPR曲線的移動，可以推算出葡萄糖溶液折射率的變化．由圖 1可以看出，信號曲線噪聲大，吸收峰不尖銳，信號不能滿足測量的要求．引起噪聲的原因包括輸出光源能量的波動、光源的散彈噪聲、光電接收器的噪聲和電路噪聲等，實驗中應(yīng)予以濾除．

圖1 典型的SPR輸出曲線Fig.1 Typical output curve of SPR sensor

2 SPR檢測葡萄糖時的數(shù)據(jù)處理方法

SPR檢測葡萄糖溶液的數(shù)據(jù)處理方法如下：對原始信號進行濾波；選擇計算基線；對基線以下的數(shù)據(jù)進行分析，尋找最低點的位置，用曲線最低點的變化來表征SPR共振角的變化．

2.1 SPR數(shù)據(jù)預(yù)處理

為了提高分析精度，首先需要對信號進行預(yù)處理，預(yù)處理包括2個步驟．

(1)采用多次平均的方法來減少偶然誤差，將時間間隔很近的 5組數(shù)據(jù)平均，并通過比較，剔除掉突變的奇異點．

(2)采用移動平均濾波，通過每個像元左右各點的累加平均，采用 5點移動平均平滑濾波，降低某些個別像元噪聲的影響．

2.2 動態(tài)基線法

基線是在 SPR測量輸出曲線上人為設(shè)定的直線，基線以下的點用于最終的計算，基線以上的點不參與計算．基線選擇過大，會包含冗余的甚至錯誤的測量點；基線選擇過小，則包含的測量點不足，信息量不夠；兩者都會造成最終的曲線擬合結(jié)果與實際情況產(chǎn)生偏差．SPR檢測葡萄糖系統(tǒng)的輸出會受外界光強變化、電路系統(tǒng)噪聲、機械結(jié)構(gòu)振動等因素影響，所以采用固定的基線進行分析就會產(chǎn)生誤差，需要采用動態(tài)基線法．本文分析和討論了 3種動態(tài)基線調(diào)整的方法．固定半峰寬度方法基于表面等離子共振峰的半峰寬不變的假設(shè)，基線為峰值和最低值的平均值；固定面積比動態(tài)基線分析中，移動基線位置，使基線以上包含面積和基線以下面積之比為固定常數(shù)；固定點數(shù)基線法動態(tài)調(diào)整基線，使基線下包含的點數(shù)為固定值．

在實驗中，建立了 SPR響應(yīng)的等效曲線，如圖 2所示．通過模擬和仿真，分析和討論了不同算法對于測量精度的影響．SPR響應(yīng)模擬曲線可表示為

式中：xi為第i個像元的位置；x0為共振角出現(xiàn)位置；Ii為像元 i的接收光強；Imax為最大光強；Idip為最大光強和最小光強差值，即信號深度；a(0 ＜ a ＜ 1)為信號的不對稱性；ω為共振峰的寬度．Imax、Idip、a和ω取值分別為 1.1、0.7、0.586和 3.9，使響應(yīng)曲線更接近實際測量的信號．如圖 2所示，隨 x0取值不同，曲線發(fā)生均勻移動．

圖2 模擬的SPR響應(yīng)曲線Fig.2 Simulation of SPR response curves

在由式(1)產(chǎn)生的模擬信號中引入不同幅值的白噪聲，比較固定點數(shù)法、固定半峰寬度法和固定面積比法在噪聲環(huán)境下的分辨率．如圖 3所示，采用固定點數(shù)的動態(tài)基線法，標準差最小，抗干擾能力最好．通過模擬計算，當固定測量點數(shù)在 30～40之間時，測量重復(fù)性最好．

圖3 噪聲下用各動態(tài)基線法計算的質(zhì)心標準差Fig.3 Standard deviation of centroid calculated using the dynamic baseline method with noise

2.3 尋找曲線最低點算法

SPR檢測葡萄糖濃度系統(tǒng)中，尋找曲線最低點的算法包括多項式擬合法、質(zhì)心法和加權(quán)質(zhì)心法．多項式擬合法尋找輸出光強最小的像元位置，然后用最小值附近若干點作為分析對象，進行多項式最小二乘擬合，計算最低點位置；質(zhì)心法通過計算基線以下部分數(shù)據(jù)幾何中心位置的移動，來判定曲線的移動，由于 SPR響應(yīng)曲線的非對稱性，質(zhì)心法計算出來的結(jié)果往往不是曲線的最低點；加權(quán)質(zhì)心法通過加權(quán)把每一個像點上的值都累加起來，其中以每個點到基線距離為權(quán)重．

2.3.1 跟蹤能力

采用不同的 x0取值，分析和比較隨 x0取值不同，采用二階擬合法、濾波二階擬合法、質(zhì)心法和加權(quán)質(zhì)心法分別計算得到最低點或質(zhì)心的變化情況．圖 4中，橫軸表示輸入 x0發(fā)生的變化，縱軸表示計算得到的最低點或質(zhì)心的變化．由圖可見，二階擬合法、濾波擬合法和加權(quán)質(zhì)心法能夠較好地跟蹤曲線的變化．質(zhì)心法計算結(jié)果波動較大．

圖4 算法跟蹤能力比較Fig.4 Comparison of algorithm tracking

2.3.2 分辨率

為驗證以上算法的分辨率，令 x0取值產(chǎn)生微小變化，比較各算法的結(jié)果．由圖 5可以看出，使用加權(quán)質(zhì)心法，當 x0發(fā)生 1 0-6變化時，計算得到的質(zhì)心能夠跟蹤輸入量的變化，線性較好，跟蹤精度較高．

圖5 算法分辨率比較Fig.5 Comparison of algorithm resolution

2.3.3 抗干擾能力

對模擬的信號引入幅值 0.05的白噪聲，在白噪聲下產(chǎn)生 10組數(shù)據(jù)，比較各種算法在噪聲環(huán)境下的分辨率．結(jié)果如圖 6所示，加權(quán)質(zhì)心法標準差最小，抗干擾能力最好，二階擬合法標準差最大．

圖6 算法抗干擾能力比較Fig.6 Comparison of algorithm anti-jamming capability

綜合以上分析可見，加權(quán)質(zhì)心法在跟蹤能力、分辨率、抗干擾能力方面都優(yōu)于其他的分析方法，適合于SPR葡萄糖濃度檢測系統(tǒng)，是有效簡便的計算方法．

3 實驗結(jié)果分析

綜合以上分析，采用固定點數(shù)的動態(tài)基線法和加權(quán)質(zhì)心法處理SPR數(shù)據(jù)，提高了SPR系統(tǒng)的檢測分辨率．為驗證算法的效果，使用SPR系統(tǒng)測量去離子水，記錄折射率隨時間變化的情況．保持室溫為 26,℃，用恒溫水浴系統(tǒng)控制去離子水進樣溫度，通過微量注射泵保證測量流速為20 μL/min，監(jiān)測1,h．

當僅采用質(zhì)心法處理數(shù)據(jù)時，測量得到去離子水的折射率在 1.332,985～1.333 000之間，對其做方差統(tǒng)計得到測量標準差為 7 × 1 0-6IU(IU 為折射率)，如圖 7(a)所示．用改進后的算法計算，測量得到去離子水的折射率在 1.332,994～1.333,002之間，對其做方差統(tǒng)計得到測量標準差為 4 × 1 0-6IU，如圖 7(b)所示．利用改進后的算法，能得到更好的測量重復(fù)性．

使用改進后的算法處理 SPR葡萄糖測量的數(shù)據(jù)．配置濃度分別為0.005～200,g/L的葡萄糖水溶液，在室溫 28,℃下進行實驗，實驗用水為去離子水，葡萄糖為分析純葡萄糖．每次更換不同濃度的葡萄糖溶液時都通入去離子水洗滌管道，實驗結(jié)果如圖 8所示．由圖 8(a)可以看出，隨測量葡萄糖溶液濃度的增大，SPR響應(yīng)曲線向左移動，SPR共振角位置向左移動，通過 SPR共振角位置可以計算溶液折射率，測量溶液折射率增大．圖 8(b)表示不同葡萄糖溶液的折射率變化．

圖7 改進算法重復(fù)性分析Fig.7 Repeatability analysis of the improved algorithm

圖8 測量不同濃度葡萄糖溶液的折射率變化Fig.8 Refractive index change when measuring different Fig.8 concentrations of glucose solution

4 結(jié) 語

使用改進后固定計算點數(shù)的動態(tài)基線法和加權(quán)質(zhì)心法結(jié)合算法有效地提高了 SPR系統(tǒng)的分辨率，使用改進后的算法處理 SPR葡萄糖檢測中得到的數(shù)據(jù)，獲得了很好的測量精度，驗證了理論分析與模擬仿真的正確性．

［1］Kurihara Kazuyoshi，Suzuki Koji. Theoretical understanding of an absorption-based surface plasmon sensor based on Kretchmann’s theory[J].Anal Chem，2002，74：696-701.

［2］Lam W W，Chu L H，Wong C L，et al．A surface plasmon resonance system for the measurement of glucose in aqueous solution[J].Sensors and Actuators B：Chemical，2005，105：138-143.

［3］王海明，錢凱先．表面等離子共振技術(shù)在生物分互作用研究中的應(yīng)用[J]. 浙江大學(xué)學(xué)報：工學(xué)版，2003，37(3)：354-357.

Wang Haiming，Qian Kaixian. Principle of surface plasmon resonance biosensor technology and its application in biomolecular interacting analysis[J].Journal of Zhejiang University：Engineering Science，2003，37(3)：354-357(in Chinese).

［4］Stenberg E，Persson B，Roos H. Quantitative determination of surface concentration of protein with surface plasmon resonance using radiolabeled proteins[J].J Colloid InterfSci，1991，143：513-526.

［5］Spadavecchia J， Manera M G， Quaranta F．Surface plasmon resonance imaging of DNA based biosensors for potential applications in food analysis[J].BiosensBioelectron，2005，21(6)：894-900.

［6］Marquette C A，Blum L J．State of the art and recent advances in immunoanalytical systems[J].Biosens Bioelectron，2006，21：1424-1433.

［7］Johnston K S，Yee S S，Booksh K S. Calibration of surface plasmon resonance refractometers using locally weighted parametric regression[J].Anal Chem，1997，69(10)：1851-1884.

［8］Kurihara Kazuyoshi，Nakamura Kaori，Suzuki Koji.Asymmetric SPR sensor response curve-fitting equation for the accurate determination of SPR resonance angle[J].Sensors and Actuators B：Chemical，2002，86(1)：49-57.