岑 唯,王 剛,廖 平
(中南大學(xué),高性能復(fù)雜制造國家重點實驗室,湖南長沙 410083)
目前,便攜式的重金屬濃度檢測儀大多使用單片機、DSP等微處理器采集和處理檢測信號[1]。但是這類處理器計算能力有限,在預(yù)處理數(shù)據(jù)、擬合實驗曲線和求峰電流時稍顯吃力。若用電腦處理檢測數(shù)據(jù),又會降低便攜性。因此提出一種以安卓手機為上位機,基于三次樣條插值算法和導(dǎo)數(shù)尋峰算法的便攜式重金屬濃度檢測系統(tǒng)。
陽極溶出伏安法是便攜式重金屬濃度檢測系統(tǒng)中最常用的方法之一[2],本文以陽極溶出法為原理基礎(chǔ)設(shè)計檢測系統(tǒng)。本系統(tǒng)采用以STM32單片機為核心的下位機采樣電化學(xué)反應(yīng)產(chǎn)生的電流信號[3],再通過手機對采樣的離散信號進行插值處理后求出溶出峰峰電流。系統(tǒng)實驗表明本便攜式檢測系統(tǒng)能準確還原溶出伏安曲線、準確識別溶出峰出現(xiàn)的位置并計算出峰電流,可廣泛應(yīng)用于室外重金屬濃度快速檢測領(lǐng)域。
陽極溶出法通過在三電極體系上施加掃描電壓促進電極和溶液組成的電解池發(fā)生反應(yīng)完成檢測,可以分為電解富集和電解溶出兩個反應(yīng)過程。在富集階段,根據(jù)待測重金屬的種類給參比電極施加特征負電位E,經(jīng)過時間t0后,待測重金屬通過陰極還原富集到工作電極上。在溶出階段,給參比電極施加由正到負的掃描電位,使富集在工作電極上的重金屬重新溶解,記錄并繪制溶出階段工作電極上產(chǎn)生的電流和所施加掃描電位的伏安特性曲線。最后通過算法識別伏安特性曲線上溶出峰的位置和高度判斷重金屬的種類和濃度[4-5]。 常見重金屬離子Cu2+、Pb2+、Cd2+、Zn2+的陽極溶出伏安特性曲線如圖1所示。
圖1 陽極溶出伏安特性曲線示意圖
Cu2+、Pb2+、Cd2+、Zn2+溶出峰出現(xiàn)的特征電位約為-0.17 V、-0.6 V、-0.8 V、-1.2 V。以汞電極進行陽極溶出法測試為例,根據(jù)Ramdies-Sevcik公式[6],溶出峰峰電流為
(1)
式中:K為常數(shù);n為參與上述電極反應(yīng)的電子數(shù);D為被測物質(zhì)在溶液和金屬汞中的擴散系數(shù);η為黏度;ω為富集時攪拌溶液的速度;v為電位掃描速度;r為汞膜的半徑;t為富集時間;C為待測溶液中重金屬濃度。
在控制公式中其他影響因素不變的情況下,溶出峰電流的大小只與離子濃度有關(guān),且成一定比例關(guān)系。只需通過對一系列標準濃度的重金屬溶液進行標定實驗即可得到峰電流和濃度之間的線性關(guān)系式。
因此只要確定伏安特性曲線溶出峰出現(xiàn)的位置和溶出峰峰電流的大小即可確定待測溶液中重金屬的種類和濃度。
本系統(tǒng)采用絲網(wǎng)印刷電極CL10組成陽極溶出法所需的三電極體系,其結(jié)構(gòu)如圖2所示。絲網(wǎng)印刷電極是一種微小型化學(xué)電極,其采用絲網(wǎng)印刷技術(shù)制作,將工作電極(WE)、對電極(CE)和參比電極(RE)集成于一體。只需將電極插到電路板的接口,再將工作區(qū)域浸入待測溶液中即可開始檢測。
圖2 絲網(wǎng)印刷電極
本系統(tǒng)硬件系統(tǒng)框圖如圖3所示。上位機用于實驗數(shù)據(jù)處理和結(jié)果展示,下位機用于采集并回傳實驗數(shù)據(jù)。上位機為下位機供電的同時通過USB通信。
圖3 硬件結(jié)構(gòu)框圖
掃描電壓模塊用于給恒電位模塊提供各種電壓波形,主要由DAC8831、電壓跟隨器和有源濾波器組成,其電路結(jié)構(gòu)如圖4所示。
圖4 掃描電壓電路簡圖
恒電位模塊是本檢測系統(tǒng)的核心電路之一,其主要由運算放大器OP07搭建組成,其電路簡圖如圖5所示。
圖5 恒電位電路簡圖
由瞬時極性法易知,U2和U3組成了串聯(lián)負反饋,U3滿足“虛斷”和“虛短”。因此有
V1=0
(2)
(3)
式中:V1為U3反相輸入端的電位,V;Vscan為掃描電壓模塊的輸出電位,V;VR為參比電極的電位,V。
聯(lián)立解得式(2)、式(3):
(4)
工作電極在微電流檢測電路中通過運放接“虛地”,所以VW=0。工作電極與參比電極的電位差為
UWR=VW-VR=R3/R6·Vscan
(5)
令R3=R6即可實現(xiàn)UWR=Vscan。
微電流檢測模塊[7]用于將工作電極產(chǎn)生的微弱電流信號轉(zhuǎn)換成0~2.5 V的電壓信號。為提高信噪比,本儀器采用多級放大、邊放大邊濾波的方案,電路結(jié)構(gòu)如圖6所示。第一級電路為I/V轉(zhuǎn)換電路[8],運放采用低偏置電流的ICL7650。利用多路復(fù)用開關(guān)用軟件的方式實現(xiàn)可變增益放大電流信號,并采用兩個多路復(fù)用開關(guān)并聯(lián)去除多路復(fù)用開關(guān)導(dǎo)通壓降帶來的影響。多路開關(guān)1的輸出電壓為
圖6 微電流檢測電路簡圖
V1=Iw(Rg+Ron)
(6)
式中:Rg為跨阻放大電路的增益;Ron為多路開關(guān)的導(dǎo)通電阻;Ω;Iw為流經(jīng)工作電極的電流,μA。
在Rg較小時,Ron的影響將不可忽略。多路開關(guān)2的輸出端接在有源濾波器的輸入端,無電流流過多路開關(guān)2,其輸出電壓為式(7)。從多路開關(guān)2輸出端引出電壓至后級電路即可消除Ron的影響。
V2=IwRg
(7)
第二級電路為壓控電壓源二階低通濾波電路,用于濾除高頻干擾,電路特征頻率為50 Hz,電壓放大倍數(shù)為2。第三級為反相放大電路,放大倍數(shù)為20。微電流檢測電路的總增益為40Rg。第四級為電壓抬升電路,利用同相求和電路將前級輸出電壓抬升1 V。
ADC采用16位芯片LTC2488,其輸出接口為SPI,可直接與單片機直連。該芯片具有4個輸入通道(2個差分輸入通道),在采樣微電流檢測電路輸出電壓的同時,實時監(jiān)測參比電極的電壓。
根據(jù)采樣數(shù)據(jù)還原溶出伏安曲線是準確測量峰電流的前提。系統(tǒng)用三次樣條插值法將離散數(shù)據(jù)還原出陽極溶出伏安曲線,并用導(dǎo)數(shù)尋峰算法確定峰值和基線所在位置。峰值到基線的豎直距離即是峰電流值,由峰電流值最終確定重金屬的濃度。
樣條插值法是數(shù)值逼近的重要分支,其在保留低次插值多項式優(yōu)點的同時,還能提高插值函數(shù)的光滑性。本系統(tǒng)選用應(yīng)用最廣泛且具有連續(xù)二階導(dǎo)數(shù)的三次樣條插值函數(shù)。對于表1給定的一組離散數(shù)據(jù),其中a=x0 表1 離散數(shù)據(jù)表 S(x)的求解步驟如下: 列出每個子區(qū)間的表達式: Si(x)=ai0+ai1x+ai2x2+ai3x3 (8) 其中i=1,2,…,n。因此確定S(x)需要確定4n個待定系數(shù)。 根據(jù)連續(xù)性條件得到如下方程: S(xi-0)=S(xi+0) (9) S′(xi-0)=S′(xi+0) (10) S″(xi-0)=S″(xi+0) (11) S(xi)=yi (12) 其中式(9)~式(11)中,i=1,2,…,n-1;式(12)中,i=0,1,2,…,n。該方程組共有4n-2個方程。最后加上自然邊界條件: S″(x0)=S″(xn)=0 (13) 即可得到4n個方程。解出這些方程即可確定S(x)。 本系統(tǒng)用異常數(shù)據(jù)剔除算法將某濃度硫酸鋅溶液的檢測數(shù)據(jù)預(yù)處理后,用三次樣條插值法進行插值,插值結(jié)果如圖7所示。曲線1為去除異常點后的插值曲線,曲線2為直接插值的曲線。曲線1光滑連續(xù)且溶出峰上無多余干擾峰,與理論曲線更接近,說明三次樣條插值法可以高度還原伏安曲線。 圖7 剔除異常點前后插值曲線對比 三次樣條插值函數(shù)圖像具有光滑的特性,因此本文使用導(dǎo)數(shù)法尋找峰值。同時檢測多種重金屬離子時,溶出伏安曲線上將出現(xiàn)多個溶出峰。為了更準確識別到各溶出峰的位置和高度,首先根據(jù)重金屬的特征溶出電位將曲線分為[-1.3,-1.1]、[-0.9,-0.7]、[-0.7,-0.5]、[-0.3,-0.1]4個區(qū)間,如圖1所示。然后分別求出4個區(qū)間內(nèi)一階導(dǎo)數(shù)為0、二階導(dǎo)數(shù)小于0且數(shù)值最大的點,取該點為此區(qū)間內(nèi)的溶出峰頂點,該點的橫坐標值為溶出峰出現(xiàn)的電位值,通過該值確定重金屬的種類。 以溶出峰頂點為分界線,在曲線的定義域內(nèi)分別向左向右搜尋一個電流值為0左右且離該頂點最近的一個極小值點。取過這兩個點的直線為該溶出峰的基線。 在室溫25 ℃的條件下,各取等量質(zhì)量濃度為200、500、800、1 000 μg/L的鉛離子溶液加入電解池中,并加入等量500 μg/L的鉍離子溶液。 通過手機設(shè)置下位機掃描電壓參數(shù)如表2所示。 表2 掃描電壓參數(shù)表 設(shè)置富集時間設(shè)置為3 min,為模擬野外無攪拌器的情況,富集過程中未做攪拌。靜息時間為30 s。將絲網(wǎng)印刷電極CL10的工作區(qū)域端浸入待測溶液中,另一端接到下位機接口,然后開始檢測。檢測結(jié)果如表3所示。 表3 檢測數(shù)據(jù)表 由表3可知,溶出峰電位都在-600 mV附近,與理論相符合。用最小二乘法對濃度和峰電流進行擬合,得到本系統(tǒng)測量鉛離子的線性擬合方程為 y=0.002 4x+1.569 式中:y為峰電流,μA;x為鉛離子濃度,μg/L。 線性相關(guān)度為R2=0.998 6。實驗結(jié)果表明在一定濃度范圍內(nèi),使用本系統(tǒng)進行重金屬濃度檢測時得到的溶出峰峰電流值和濃度值具有良好的線性關(guān)系。 本文提出了基于安卓的便攜式重金屬濃度檢測系統(tǒng),系統(tǒng)中下位機負責采集絲網(wǎng)印刷電極的電流信號并上傳到手機。利用手機處理器的強大性能采用三次樣條插值法還原出伏安特性曲線,然后用導(dǎo)數(shù)法尋找溶出峰和基線的位置,最終準確計算出溶出峰峰值電流的大小確定溶液中重金屬的濃度。實驗結(jié)果表明:本檢測系統(tǒng)能夠高度還原溶出伏安特性曲線、精準識別峰值點和基線的位置,具有較高的峰電流檢測精度、便攜性和開發(fā)價值。3.2 導(dǎo)數(shù)尋峰算法
4 實驗驗證
5 結(jié)束語