陳安娜¹，曹發和^{1, *}，張鑒清¹，曹楚南^{1, 2}

       1浙江大學化學系，杭州，中國，310027

　　2金屬腐蝕與防護國家重點實驗室，沈陽，中國，110016

　　* 通訊作者：Email: nelson_cao@zju.edu.cn

　　作者簡介

　　陳安娜：2009年7月畢業于浙江大學化學系，獲學士學位；同年9月保送至浙江大學化學系電化學與功能材料研究所攻讀碩士研究生，師從曹發和副教授。目前研究興趣為典型金屬的模擬大氣腐蝕研究和電化學監檢測技術在腐蝕失效過程中的應用，特別是電化學噪聲技術。針對電化學噪聲信號隨機波動且具有一定確定性的特性，應用包括散粒噪聲理論和混沌理論等多種數據提取和分析方法，獲取典型金屬在模擬大氣條件下腐蝕過程和金屬電沉積晶體生長機理的有效信息，主要致力于在時域、頻域、混沌理論（包括分形理論）及其他非線性時間序列分析方法的建立，對動力學狀態空間中的波動數據的熵值、關聯維數、分形維數、Lyapunov指數等性質進行考察，并結合電化學基本原理，明晰腐蝕機制和相應腐蝕反應動力學參數。

　　相關論文如下：

　　1.Study of pitting corrosion on mild steel during wet-dry cycles by electrochemical noise analysis based on chaos theory, Anna Chen, Fahe Cao, Xiaoning Liao, Wenjuan Liu, Liyun Zheng, Jianqing Zhang, Chunan Cao, Corrosion Science, submitted.

　　2.Application of shot noise analysis for corrosion process of zinc alloy （ZnAl4Cu1） under dry-wet cycles, Anna Chen, Fahe Cao, Wenjuan Liu, Liyun Zheng, Zhao Zhang, Jianqing Zhang, Chunan Cao, Trans. Nonferrous. Soc. China, in press.

　　3.Characterization of Corrosion during Dry-Wet Test by Current Noise Fractal Analysis,Anna Chen, Fahe Cao, Jianqing Zhang, “2010年全國腐蝕與電化學及測量方法大會”收錄。

陳安娜

       摘要：本文用電化學噪聲法檢測低碳鋼（Q235）在0.1 M NaCl和0.1 M Na₂SO₄溶液干濕循環腐蝕試驗中亞穩態點蝕的發生特性。電流噪聲（ECN）的混沌分析中最大Lyapunov指數（LLE）正值的百分率越高，亞穩態點蝕的數目越多；LLE數值越正，亞穩態蝕點分布越均勻。文章對監控的5天15個干濕循環內，每小時15分鐘后11段電流噪聲進行先降噪后計算LLE的處理。LLE計算結果可預測在Na₂SO₄溶液中，低碳鋼表面將發生大量亞穩態蝕點的產生和修復過程，并且由較大的LLE正值顯示這些蝕點的分布較在NaCl溶液中均勻。同時，電流噪聲LLE計算結果表明Q235在濕循環中比干循環中產生更多可修復的亞穩態點蝕且濕循環中蝕點的分布更加均勻。

　　關鍵詞：電化學噪聲，亞穩態點蝕，低碳鋼，混沌分析

　　1 引言

　　眾所周知，亞穩態蝕點的成核、生長、滅亡將引起電位、電流的波動。腐蝕介質中，鈍性金屬（如鐵、鋁）表面亞穩態蝕點的生長是迅速的、不可避免的。具有自鈍化特性的金屬，對孔蝕的敏感性較高。Y.F. Cheng認為電位暫態峰緩慢恢復是由鈍化膜電容的充放電引起的^[1]。而電流暫態峰的來源類型可以由功率密度譜衰減斜率的不同分為6種^[2]。其中，純鐵和低碳鋼上亞穩態點蝕通常引發電流先迅速上升后指數衰減回到基線。通過理論計算，此類電流峰功率密度譜衰減斜率為-2 mA²/Hz²。Speckert等人^[3]的研究表明，Ti-Al合金表面亞穩態點蝕將產生正向電流峰，Fe-Al合金表面亞穩態點蝕能產生雙向電流峰。Dong等人采用電化學噪聲研究了模擬混凝土孔隙液中浸泡的碳鋼表面亞穩態孔蝕產生和發展與銹層的關系，該研究發現銹層可以引導亞穩態點蝕的生成，生成速率在出現宏觀蝕孔以后達到最大值^[4]。#p#分頁標題#e#

　　為了原位監檢測鋼鐵等金屬在潮差區、泥土內、微生物影響下的腐蝕行為，電化學噪聲法已被廣泛采用。傳統的電化學噪聲分析方法主要包括時域分析、頻域分析（常用快速傅里葉變換法FFT,最大熵值法MEM和小波變換法WT）等。另一方面，由于腐蝕的非平穩性和不可預知性，所表現出的電化學噪聲具有混沌特性。混沌是指在確定性系統中出現類似隨機波動的過程，普遍認為混沌系統應該具有以下幾個特征：非線性、確定性、對初始條件的依賴性等。李雅普諾夫指數（LE）是刻畫混沌系統的主要特征之一，它表示耗散體系相空間中相體積收縮過程中幾何特征變化的物理量。一維映射只有一個LE,它可能大于、等于、小于零。對于正的LE來說，就意味著信息量的損失。存在LE大于0的情況下，體系才被認為具有混沌的性質。

　　本文采用電化學噪聲法，基于混沌分析參數LLE，研究了在0.1 M NaCl和0.1 M Na₂SO₄溶液干濕循環腐蝕試驗中，低碳鋼Q235表面亞穩態點蝕發生及分布的情況。

　　2 實驗方法

　　2.1 干濕循環實驗

　　用于海洋工程等自然環境中的鋼材經常暴露在短暫濕環境和持久性風干環境相交替的物理環境中。本試驗設計碳鋼Q235材料在0.1 M NaCl和0.1 M Na₂SO₄溶液中進行模擬自然環境干濕循環實驗，實驗裝置可見文獻^[5]。Q235的組成為：C 0.19%， Si 0.22%， Mn 0.56%， P 0.0086， S 0.22， 其余為鐵。試驗中一個循環為8小時，包括1個小時溶液噴淋階段和7個小時自然風干階段，分別稱為濕循環和干循環。

　　2.2 電化學噪聲測量

　　電化學噪聲測量采用三電極體系。工作電極和對電極均為Q235鋼，試樣為底部工作面積為0.20 cm²的圓柱，用銅導線鉚接好后，非工作面用聚甲基丙烯酸甲酯封成電極。實驗前用400#至1000#砂紙逐級打磨，拋光至鏡面，用95%的乙醇擦拭除油，水洗，冷風吹干后放入干燥箱保存。參比電極為飽和甘汞電極（SCE）。實驗采用美國Gamry公司PCI4-G300恒電位儀ESA410軟件實時并且同時測量雙電極偶合時電位、電流的波動。采樣頻率為5 Hz。

　　2.3 銹層表征

　　干濕循環結束后用XRL3400光學顯微鏡和SRION-100電子掃描顯微鏡對工作電極表面形貌進行觀察分析。電極表面產物晶相和組成分別由PANalytical XPert PRO MPD X射線衍射（XRD）儀和JDbin-yvon LabRam HRUV Raman光譜儀表征。XRD分析采用Cu-Kα靶，功率為40 kV×40 mA, 掃描速度為 2 °/min。

　　3 電化學噪聲數據處理

　　3.1 電流數據和頻域譜獲取

　　本文取干濕循環每小時15分鐘后11段電流時間序列做數據分析，每段時間序列長204.8s，即包括1024個數據點。一般亞穩態點蝕引起的電流峰持續時間為幾秒至幾十秒，因此所選時間長度完全可以收集亞穩態點蝕從發生到滅亡的信息。碳鋼在兩種電解液中交替進行15個濕循環和15個干循環，因此在單種溶液中，進行混沌分析的時間序列總數是1320個。由于電流噪聲來源廣泛，包括熱噪聲、散粒噪聲、閃爍噪聲、儀器噪聲等，為了得到與點蝕有關的確定性的電流組分，我們引入了奇異值分解降噪的辦法。

　　矩陣奇異值分解屬于矩陣對角化分解的一種類型，它是通過正交變換，將矩陣A分解為A=USV^T的形式。其中U和V為酉矩陣，S為和A具有相同維數的對角矩陣，其中的元素λ稱為奇異值。用表示無噪聲污染的信號矩陣，構造，S_r代表無噪聲的信號空間，S₀是代表噪聲空間。將λ_i(i>r)置為零，僅利用U、S、V的前r列計算得到無噪聲污染的矩陣。本文選取r=4，選取依據可參閱文獻^[6]。

　　文中采用的頻譜轉換方法為傅里葉變換和小波變換法^[7]。根據特征的小波系數實現對電化學噪聲信號的多尺度分辨空間傳播性的檢測。

　　3.2 相空間重構及LLE計算

　　相空間重構是根據有限的數據來重構吸引子以研究系統動力行為的方法，其基本意思是：系統中任一分量的演化都是由與之相互作用的其它分量所決定的。構造一個等價的狀態空間，只需考察一個分量，并將在某個固定的時間延遲點上的測量作為新維處理。重復這一過程并測量相對于不同時間延遲的量，就可以產生出許多這樣的點，它可以將動力系統的許多行為保留下來。重構相空間需要選擇兩個重要參數，時間延遲t和嵌入維數m， 本文中應用自相關法計算t， Cao方法計算m^[8]。

　　對電流時間序列X₁，X₂ …，X_n相空間重構后有：

　　Wolf^[9]給出了在一維數據中提取Lyapunov指數的方法。以初始點為基點V₀，選取最近的點V_η0，構成初始向量，V₀與V_η0間的距離可記為L₀，取時間步長為t，初始向量沿軌線向前演化得到一新向量，其相應基點與端點間的歐氏距離可記為L₁’，如此直至所有相點，然后取各指數增長率的平均值為最大Lyapunov指數估計值。#p#分頁標題#e#

　　4 結果與討論

　　4.1 Q235干濕循環下的ECN

　　兩個相同材料的低碳鋼樣品在0.1 M NaCl和0.1 MNa₂SO₄溶液中干濕循環，由于兩電極在干濕循環下擴散層厚度的不同，引起電極間極化差異不同。根據Klapper等人對電化學噪聲的定義，由電流的絕對值，可以判斷得出濕循環中兩電極間不對稱性大，干循環中不對稱性小（如圖1）。Na₂SO₄溶液中，兩個樣品的不對稱性在前六個小時逐漸增大并達到相對穩定。NaCl溶液里第一個濕循環中兩電極不對稱性最大。這可解釋為Cl^-離子對鐵的局部吸附作用。Fe的晶格邊緣由于Cl^-的侵蝕首先出現缺陷，然后鐵的氧化反應均勻發生。

　　圖1. Q235在（a） 0.1 M Na₂SO₄和（b） 0.1 M NaCl中干濕循環電位、電流圖。

　　如3.1節所述，選取局部的電流序列去除漂移、SVD降噪后做FFT分析得到的頻譜均具有特征頻率， 可認為經過降噪的電流值具有一定周期性。小波分解結果表明波動能量集中在高階晶胞，引起電流波動的反應為慢反應。鐵的溶解在研究范圍內受到傳質控制。中性溶液中，與鐵溶解相偶合的陰極反應為O₂的還原和銹層中β-FeOOH和γ-FeOOH還原為Fe3O4^[10]的反應。同時，Zou等人指出β-FeOOH具有更強的還原性。

　　自相關所得的時間延遲t在干循環中大于濕循環，表明干循環中，低碳鐵表面亞穩態點蝕的產生、滅亡周期大于濕循環中。這與干循環中產物累積有關。電解液在濕循環中的流動，降低了擴散層厚度，O₂容易擴散到碳鋼表面，而干循環中，產物在鋼鐵表面累積，O₂的擴散受到抑制，FeOOH的還原反應成為消耗鐵氧化產生電子的主要反應。后者反應的速率慢，延時長。

　　圖2. 點蝕成核生長模型與相關的ECN曲線示意。

　　4.2 亞穩態點蝕分布模型

　　LLE表現了電流在多維空間的演化情況。波形試驗表明，如果波動是純周期性的（單頻正弦波），無不確定性，則LLE接近于0。如果電流是趨勢上升或下降的，LLE將小于0。如果在單頻正弦波上施加一定比例的白噪聲，LLE將為正值，在閾值以下，LLE數值隨白噪聲強度增加而增加。假設一個亞穩態點蝕將引起電流某單頻正弦波，考慮兩個亞穩態蝕點在Q235表面發生的情況，如圖2。當第二個蝕點在第一個蝕點邊緣或在第一個蝕孔內“深挖”，腐蝕微電極具有大陰極小陽極的面積結構，陽極電流密度很大，第二個單頻正弦波強度應與第一個相似。疊加的波形應如圖2C。如第二個蝕點在距離第一個蝕點較遠的鈍化膜的另一缺陷處發生，其再鈍化能力比第一種情況強，其延續時間比較短，且強度小（圖2B）。

　　在氯化鈉溶液中ECN的LLE正數的百分率小于硫酸鈉溶液中。這表明當介質中含有氯離子時，氯離子優先吸附在碳鋼鈍化膜表面，蝕孔形成后，平衡便受到破壞，溶解占優勢。同時，氯化鈉溶液中LLE正數數值小于硫酸鈉，表明NaCl溶液中蝕點分布的局部性。SEM圖片可驗證這一結果。碳鋼在NaCl溶液中干濕循環48小時后可觀察到明顯的蝕孔。Na₂SO₄溶液中，樣品表面積累了較厚的腐蝕產物。由EDS和XRD譜判斷Q235在Na₂SO₄溶液中干濕循環24小時的表層腐蝕產物為γ-FeOOH， Fe₂O₃和FeSO₄，在NaCl溶液中的腐蝕產物為無定型的鐵的氧化物。循環48小時以后產物的原位Raman光譜圖表明NaCl溶液中的腐蝕產物主要為Fe₃O₄， 由此推斷NaCl溶液中的無定型產物可能為β-FeOOH。XRD和Raman光譜數據因為篇幅原因沒有列出。產物類型的差別影響了點蝕的發生和分布。模型如圖2。

　　同時，LLE計算顯示濕循環中ECN的LLE大于0的百分率大于干循環中，表示濕循環中將發生更多的亞穩態點蝕。且LLE值大于干循環中，如表1。現象可解釋為介質處于流動狀態，有利于溶解氧的輸送，亞穩態點蝕雖然發生多，但是再鈍化也較容易，蝕點均勻分布。干循環中不溶產物在金屬表面沉積，局部裸露的金屬表面容易發生穩定的孔蝕。

　　表1. 低碳鋼干濕循環中電流噪聲 （ECN） 的LLE值計算結果#p#分頁標題#e#

　　5 結論

　　LLE計算結果表明Q235在Na₂SO₄溶液干濕循環條件下，表面將發生大量亞穩態蝕點的產生和修復過程，并且由較大的LLE正值顯示這些蝕點的分布較在NaCl溶液中均勻。同時，電流噪聲LLE計算結果驗證了Q235在濕循環中比干循環中產生更多可修復的亞穩態點蝕且濕循環中蝕點的分布更加均勻的結論。

　　致謝

　　本研究得到國家自然科學基金（50801056）資助。

　　References（參考文獻）

　　[1] Y.F. Cheng, M. Wilmott, J.L. Luo. Analysis of the role of electrode capacitance on the initiation of pits for A516 carbon steel by electrochemical noise measurements. Corrosion Science, 1999, 41（7）： 1245-1256.

　　[2] Y.F. Cheng, J.L. Luo, M. Wilmott. Spectral analysis of electrochemical noise with different transient shapes. Electrochimica Acta, 2000,  45（11）： 1763-1771.

　　[3] L. Speckert, G.T. Burstein. Combined anodic/cathodic transient currents within nucleating pits on Al-Fe alloy surfaces. Corrosion Science, 2011, 53（2）： 534-539.

　　[4] Z.H. Dong, W. Shi, X.P. Guo. Initiation and repassivation of pitting corrosion of carbon steel in carbonated concrete pore solution, Corrosion Science, 2011, 53（4）： 1322-1330.

　　[5] Y.Y. Shi, Z. Zhang, J.X. Su, F.H. Cao, J.Q. Zhang. Electrochemical noise study on 2024-T3 Aluminum alloy corrosion in simulated acid rain under cyclic wet-dry condition. Electrochimica Acta, 2006, 51（23）： 4977-4986.

　　[6] B. Malki, A. Legris, J.L. Pastol, D. Gorse. Comprehensive dissolution current noise analysis during stress corrosion cracking of Cu3Au alloys, Journal of the Electrochemical Society, 1999, 146（10）： 3702-3710.

　　[7] F.H. Cao, Z. Zhang, J.X. Su, Y.Y. Shi, J.Q. Zhang. Electrochemical noise analysis of LY12-T3 in EXCO solution by discrete wavelet transform technique. Electrochimica Acta, 2006, 51（7）： 1359-1364.

　　[8] L.Y. Cao, A.S. Soofi. Nonlinear deterministic forecasting of daily dollar exchange rates. International Journal of Forecasting, 1999, 15（4）： 421-430.

　　[9] A. Wolf, J.B. Swift, H.L. Swinney, J.A. Vastano. Determining Lyapunov Exponents from a Time-Series. Physica D, 1985, 16（3）： 285-317.

　　[10] Y. Zou, J. Wang, Y.Y. Zheng. Electrochemical techniques for determining corrosion rate of rusted steel in seawater. Corrosion Scienc, 2011, 53（1）： 208-216.