華洋科儀
咨詢熱線:400-901-3696
電化學(xué)噪聲(3)——噪聲阻抗Rn的測(cè)定
  • 發(fā)布時(shí)間 : 2020-07-09 15:13:00
  • 瀏覽量:6414
  1.簡(jiǎn)介
  前兩篇文章(EN:電化學(xué)噪聲(1)和(2))表明與ASTM要求相比,Bio-Logic儀器能夠進(jìn)行可靠的噪聲測(cè)量。本文是對(duì)一個(gè)體系的噪聲阻抗Rn和極化阻抗Rp進(jìn)行比較,并用這個(gè)例子來(lái)介紹EC-Lab的噪聲分析工具。Cottis[1]在一篇文章中討論分析噪聲數(shù)據(jù)的幾種方法。如前所述,其中之一是計(jì)算腐蝕體系的特征值,稱為噪聲阻抗Rn。Rn定義為電位的標(biāo)準(zhǔn)偏差除以電流的標(biāo)準(zhǔn)偏差:
  
  在某些條件下,Rn可被視為等效于極化電阻Rp。許多研究都證明了這種等效性[2,3]。在EC-Lab和EC-Lab Express軟件中,電化學(xué)噪聲分析工具ENA可用于Rn的測(cè)定。如果系統(tǒng)的阻抗不是簡(jiǎn)單的電阻,即阻抗隨頻率變化,則需要使用功率譜密度(PSD)[2]計(jì)算標(biāo)準(zhǔn)電位和電流偏差:
  

  其中,和分別是電位和電流噪聲的PSD,分別用V2/Hz和A2/Hz表示。f是噪聲的頻率。使用ENA工具,可以通過(guò)離散傅立葉變換(DFT)或最大熵法(MEM)計(jì)算PSD。本文將描述如何使用ENA工具執(zhí)行和分析噪聲測(cè)量。


  2.實(shí)驗(yàn)條件
  所有測(cè)試的電解液為ASTM文件中的0.005mol L-1的H2SO4和0.495 mol L-1的Na2SO4溶液。
  2.1 阻抗測(cè)量
  首先,采用阻抗法測(cè)定Rp,工作電極(WE)為1cm2的AISI 430不銹鋼樣品,對(duì)電極(CE)為石墨電極,參比電極(Ref)為飽和甘汞電極(SCE)。按照ASTM G5[5]中所述的程序,用240和600號(hào)砂紙對(duì)AISI樣品進(jìn)行處理。所用技術(shù)為PEIS,條件如圖1所示。在進(jìn)行阻抗測(cè)量之前,將樣品在OCV下浸泡在溶液中1小時(shí),以確保其達(dá)到穩(wěn)定狀態(tài)。




  2.2 噪聲測(cè)量
  使用兩個(gè)1cm2的AISI 430不銹鋼樣品作為WE和CE,一個(gè)SCE(飽和甘汞電極)作為參比電極進(jìn)行噪聲測(cè)量。樣品按照前面描述的程序拋光。使用的技術(shù)是ZRA,條件如圖2所示。電位范圍為[-0.255V;-0.206V],電位測(cè)量分辨率為0.76μV。ZRA技術(shù)的原理是在工作電極(S1)和對(duì)電極(S3)之間保持0V的電壓。在WE和CE之間產(chǎn)生的電流流動(dòng)即為ECN(電化學(xué)電流噪聲)。Ewe/Ece(S1/S3)對(duì)參比電極(S2)的電位即為EPN(電化學(xué)電位噪聲)。建議在表示噪聲幅度時(shí),不要試圖按面積[1]進(jìn)行標(biāo)準(zhǔn)化。
  采樣時(shí)間取決于感興趣問(wèn)題的頻率范圍。根據(jù)Shannon定理,如果感興趣的現(xiàn)象具有頻率f,則采樣頻率必須至少等于2f。所使用的條件如圖2所示。設(shè)置CE的E range為+/-2.5 V。
  測(cè)量是在法拉第籠里進(jìn)行,恒電位儀接地。噪聲的頻率上限和下限分別為fmin=1/(30*60*30*60*10)=30nHz,fmax=1/(2*0.1)=5Hz。我們還使用了一個(gè)5Hz模擬濾波器來(lái)消除高于5Hz頻率的混疊。此篩選器可用于SP-200、SP-240、SP-300和VSP-300。


2ZRA條件


  3.結(jié)果
  3.1 阻抗測(cè)量
  圖3為上述阻抗測(cè)試的Nyquist圖。使用Zfit分析工具可以看出,系統(tǒng)阻抗的等效電路為:R1+Q2/R2。使用最小化算法(Randomize+Simplex)獲得的參數(shù)值為:
  R1=RΩ(電解液電阻)=70±0.3Ω
  R2=Rp(極化電阻)=76 000±500Ω
  C2(用Q2表示的偽雙層電容)=167000±66pF[6]。




3AISI430鋼在0.005mol L-1H2SO40.495 mol L-1Na2SO4溶液中Ecorr下的阻抗Nyquist

3.2 ZRA技術(shù)進(jìn)行噪聲測(cè)量

4中可以看到EPNECN的演化。噪聲實(shí)驗(yàn)的總持續(xù)時(shí)間為30分鐘,我們選擇的數(shù)據(jù)是浸泡20分鐘后獲得的數(shù)據(jù)。ENA工具在軟件的Analysis/Corrosion子菜單中(參見(jiàn)圖5)。



4AISI430鋼在0.005mol L-1H2SO40.495 mol L-1Na2SO4溶液中的EPNECN




  圖5:電化學(xué)噪聲分析EPA工具


  4.噪聲分析
  4.1 Rn測(cè)定
  (a)方法說(shuō)明
  圖6為分析窗口的圖片。分析工具提供了去除EPN或ECN軌跡上的趨勢(shì)或漂移的可能性。趨勢(shì)消除函數(shù)實(shí)際上是從原始時(shí)間軌跡中減去通過(guò)原始軌跡的線性或二階多項(xiàng)式擬合得到的數(shù)值曲線。如圖7所示,趨勢(shì)消除產(chǎn)生的軌跡集中在0左右。



  圖6:ENA窗口


  圖6為ENA窗口和用于計(jì)算Rn的方法。第一種方法使用等式(1),其中使用標(biāo)準(zhǔn)偏差。根據(jù)式(2),還可以使用通過(guò)信號(hào)的傅里葉變換獲得的PSD來(lái)計(jì)算Rn。提出了兩種計(jì)算傅里葉變換的方法:離散傅里葉變換(DFT)和最大熵法(MEM)。
  DFT方法使用一個(gè)函數(shù)來(lái)表示頻域中的時(shí)變離散輸入。分析工具還允許在計(jì)算DFT之前獲得ECN和EPN的時(shí)間軌跡。Windowing(加窗)選項(xiàng)用于去除在計(jì)算DFT時(shí)由時(shí)間記錄開(kāi)始和結(jié)束處的不連續(xù)性產(chǎn)生的信號(hào)的頻率成分。關(guān)于趨勢(shì)消除和加窗的影響的更多細(xì)節(jié)將在附錄中給出。在ENA工具中提出了五種不同的加窗處理方法。DFT和Windowing計(jì)算的表達(dá)式見(jiàn)EC-Lab用戶手冊(cè)的傅里葉變換部分。





  圖7:二階多項(xiàng)式擬合趨勢(shì)去除的ECN
  MEM基于這樣一個(gè)事實(shí),即任何隨機(jī)時(shí)間序列x(t)都可以表示為由白噪聲(所有頻率的振幅恒定的信號(hào))激發(fā)的濾波器(這里稱為自回歸)的輸出。白噪聲被認(rèn)為是解釋該方法名稱的最隨機(jī)信號(hào)。在ENA窗口中(圖6),可以選擇MEM的順序。這個(gè)順序代表濾波器傳遞函數(shù)的極點(diǎn)。如果階數(shù)等于n,則得到的DFT的峰數(shù)為n/2。Bertocci[7]詳細(xì)介紹了MEM在噪聲信號(hào)傅里葉變換計(jì)算中的應(yīng)用。
  (b)結(jié)果及與Rp的比較
  表1為由等式(1)定義的在各種可用條件下的噪聲阻抗。我們可以看到趨勢(shì)消除會(huì)導(dǎo)致巨大的變化。
  表1 使用標(biāo)準(zhǔn)差和所有可用趨勢(shì)消除選項(xiàng)計(jì)算的Rn

  

根據(jù)表2,用DFT方法計(jì)算的Rn值與在相同趨勢(shì)消除條件下(ECN和EPN上的二階多項(xiàng)式)用標(biāo)準(zhǔn)差得到的值大致相同。
  表2 用電位和電流的二階多項(xiàng)式趨勢(shì)和用于計(jì)算PSD的DFT方法計(jì)算Rn


表3為使用Barlett窗口和遞增順序通過(guò)MEM獲得的值。在相同的趨勢(shì)消除和加窗條件下(二階多項(xiàng)式和Barlett),用這種方法得到的值與DFT得到的值接近。有趣的是,在順序號(hào)6之后,Rn似乎達(dá)到了最小值。
  表3 用計(jì)算傅里葉變換的MEM,采用遞增階數(shù)和Barlett加窗獲得的Rn值


4.2 趨勢(shì)消除和加窗的影響
  圖8a中比較原始EPN(圖8a,■)的DFT和通過(guò)使用線性(圖8a,▲)或二階多項(xiàng)式擬合(圖8a,·)獲得的去渲染信號(hào)。當(dāng)頻率低于1Hz時(shí),在原始EPN上可以看到對(duì)數(shù)圖中斜率為-1的直分量。這種非常平滑的成分與漂移有關(guān),文獻(xiàn)[8]中也有相關(guān)介紹。圖9a中可以看出,去趨勢(shì)的影響是從DFT中去除該平滑成分。
  圖8b為加窗EPN的DFT:原始EPN(圖8b,■),使用Hamming窗口(圖8b,▲)和Barlett窗口(圖8b,·)。首先可以看出Hamming窗在去除平滑成分方面不如Barlett窗。最后,圖8c為EPN在以下情況下的DFT:無(wú)去趨勢(shì)和Barlett窗(圖8b,■)、線性擬合去趨勢(shì)和Barlett窗(圖8b,▲)、二階多項(xiàng)式擬合去趨勢(shì)和Barlett窗(圖8a,·)??梢钥闯?,二階多項(xiàng)式擬合去趨勢(shì)化和Barlett窗在最低頻率下提供了振幅最小的DFT。附錄中給出了有關(guān)取消渲染和窗口效果的更詳細(xì)信息。


圖8: 加窗和去趨勢(shì)對(duì)EPN的DFT的影響:
  a) 僅去趨勢(shì);b)僅加窗處理;c)去趨勢(shì)和加窗處理


  1.結(jié)論
  EC-Lab和EC-Lab Express中ZRA技術(shù)可用于噪聲測(cè)量。這項(xiàng)技術(shù)使工作電極和對(duì)電極之間的電壓保持在0V。本文還介紹了相關(guān)的分析工具。計(jì)算值Rn非常依賴于數(shù)據(jù)處理的類型(去趨勢(shì)化和加窗)。結(jié)果表明,這兩種方法具有相似的效果,可以去除低頻成分。必須使用EC Lab和EC Lab Express中可用的不同方法來(lái)交叉檢查Rn值;例如,在上面的示例中,Rn值約為18kΩ,因?yàn)檫@三種方法(標(biāo)準(zhǔn)偏差、DFT和MEM)都得到這個(gè)結(jié)果。這個(gè)Rn值大約是阻抗譜測(cè)定的Rp值的3倍。


  附錄
  本附錄旨在說(shuō)明加窗對(duì)信號(hào)的影響。圖9a為恒電位儀在25s到75s之間記錄的實(shí)際輸入。圖9b為對(duì)該時(shí)間記錄執(zhí)行DFT得到的假定的輸入信號(hào),其中時(shí)間記錄在整個(gè)時(shí)間內(nèi)重復(fù)。
  假定的輸入信號(hào)的不連續(xù)性會(huì)對(duì)信號(hào)產(chǎn)生額外的偽低頻。去趨勢(shì)和加窗的主要作用是去除這些偽頻率。
  通過(guò)使用線性去趨勢(shì),我們?nèi)コ藞D9a中信號(hào)的線性成分,獲得一個(gè)以0為中心的信號(hào),類似于圖7所示。上述時(shí)間記錄的DFT假定輸入如圖10a所示,可以看出不連續(xù)性已經(jīng)消失。窗口函數(shù)是僅在時(shí)間記錄上定義的函數(shù),它的邊界等于0,定義域的中心等于1。圖10b中可以看出,假定的輸入不包含圖9b所示的不連續(xù)性。

9: 時(shí)間記錄的實(shí)際輸入(a)和假定輸入(b


10: a)去趨勢(shì)化和(b)加窗對(duì)實(shí)際輸入信號(hào)的時(shí)間記錄的影響

參考文獻(xiàn)

[1] R. A. Cottis, Corrosion 57, 3 (2001) 265.

[2] U. Bertocci, C. Gabrielli, F. Huet, M. Keddam, P. Rousseau, J. Electrochem. Soc. 144 (1997) 37.

[3] J. F. Chen, W.F. Bogaerts, Corr. Sci. 37 (1995) 1839.

[4] J.R. Kearns, J.R. Scully, eds., “Electrochemical Noise Measurements for Corrosion Applications” ASTM STP 1277, Annual Book of ASTM Standards (West Conshohoscken, PA: ASTM,1996) 446.

[5] ASTM G5-94 Standard Reference Test Method for Making Potentiostatic and Potentiodynamic Anodic Polarization Measurements.

[6] Pseudo capacitance calculation, Applica- tion note #20, http://www.bio-logic.info/potentiostat/notesan.html

[7] U. Bertocci, J. Frydman, C. Gabrielli, F. Huet, M. Keddam, J. Electrochem. Soc. 145 8 (1998) 2780.

[8] U. Bertocci, F. Huet, R. P. Nogueira, P. Rousseau, Corrosion NACE 58 4 (2002) 337.





友情鏈接: 友情鏈接2 友情鏈接3 友情鏈接4 友情鏈接5 友情鏈接6 友情鏈接9 友情鏈接10