激光光譜法同步測量尿素水溶液液膜厚度與濃度

楊薈楠，郭曉龍，楊斌，蘇明旭，蔡小舒

?

激光光譜法同步測量尿素水溶液液膜厚度與濃度

楊薈楠，郭曉龍，楊斌，蘇明旭，蔡小舒

(上海理工大學能源與動力工程學院，上海市動力工程多相流動與傳熱重點實驗室，上海200093)

對尿素水溶液液膜厚度和濃度進行定量分析在工業應用中至關重要，而傳統的測量方法只能對液膜厚度和濃度分別進行測量。提出了一種用于尿素水溶液液膜厚度和濃度高精度同步在線測量的新方法。通過靈敏度分析和內部交叉驗證，選取了測量精度及靈敏度最高的一對波長組合(7040.75 cm-1與 6761.00 cm-1)。由內部交叉驗證可知，該方法測量液膜厚度誤差小于1.40%，濃度誤差小于5.76%。

溶液；測量；膜；濃度；厚度

引 言

液膜形成的現象廣泛地存在于各種工業過程中，例如在電子元器件噴霧冷卻中器件表面液膜的形成[1]，高效清潔滅火技術中細水霧碰壁液膜的形成[2]，水平管降膜蒸發器中水滴撞擊換熱管表面液膜的形成[3]，以及選擇性催化還原系統(SCR)中汽車尾氣排放管上的尿素水溶液液膜的形成等[4]，對液膜多個參數（厚度、濃度等）進行定量分析不僅能更好地了解液膜形成和蒸發這個極其復雜的物理過程的本質，也對優化所涉及的各種工業過程具有重要的意義。

傳統的測量方法只能對尿素水溶液液膜的厚度與濃度分別進行測量。液膜厚度的測量主要分為電學和光學兩大類。傳統的電學膜厚測量方法包括電容法和電導法，電導法基本原理是基于探針區域液膜電阻的變化從而測得液膜厚度[5-6]。電容法測量膜厚的基本原理為：兩塊平行金屬極板之間形成電容，當其間的介質或平板距離發生變化時，引起輸出電壓的變化[6]。而電學法測量精度和測量范圍易受到如液膜波動形態、平板表面液膜覆蓋情況和探針結構參數等因素的影響[7-9]。光學測量方法，特別是基于激光光譜的測量方法，因其具有無干擾、靈敏度高、響應快等優點在液膜測量方面更具優勢[10-12]。如Yang等[13]利用可調諧半導體激光吸收光譜技術對純水水膜蒸發過程中水膜的厚度(5～1000mm)、溫度(20～75℃)以及水膜上方水蒸氣的溫度進行了同步測量，Greszik等[14-15]利用拉曼散射法和激光誘導熒光技術獲得了透明石英片上純水水膜的二維分布信息，但測量結果需進行標定且設備昂貴。溶液濃度常見的測量方法主要包括折光法、密度法和超聲波聲速法等。折光法利用介質的濃度與光的折射率的相關性對溶液的濃度進行測量，雙波長分光光度計方法[16]也可實現尿素溶液濃度無接觸測量，但不能對液膜厚度進行同步測量。雖然基于密度法和超聲波聲速法而研發的濃度計已成熟地運用在如水泥脫硝工藝中的尿素水溶液濃度在線測量，但這兩種方法均需與溶液接觸，因此它們不能滿足對液膜進行高精度檢測的需求。

本文提出了一種基于激光光譜同步測量尿素水溶液液膜厚度和濃度的新方法。利用傅里葉變換紅外光譜儀對不同濃度(質量分數為5%～50%)的尿素水溶液進行測量，建立基于比爾-朗伯定律的液膜厚度及濃度同步反演模型，通過對尿素水溶液吸收率與濃度的關系進行優化分析，選取兩個波長，從而實現對尿素水溶液液膜厚度及濃度高精度和高敏感度的同步測量。并對該波長組合進行交叉驗證，獲得該方法的測量誤差。

1 理論與建模

液態水在近紅外區域存在一個寬帶的吸收光譜[17]，Yang等[18]精確地測定了液態水在該范圍內(5800～7800 cm-1)不同溫度下(25～75℃)的紅外吸收截面。Halbout等[19]測定了尿素單晶體在近紅外區域的透射曲線，然而現有研究中并沒有精確的尿素水溶液的吸收率。因此，本文將質量分數為5%～50%的尿素水溶液放置在光程為1 mm的石英樣品池中，用高分辨率的傅里葉紅外變換光譜儀(0.25 cm-1)對其進行掃描，測得了不同濃度下尿素水溶液的吸收率(,)。如圖1所示，隨著濃度的增加，尿素水溶液的吸收譜線向著波數小的方向偏移。因此，可基于比爾-朗伯定律，建立雙波長尿素水溶液液膜厚度與濃度計算的數學模型。

1.1 尿素水溶液液膜厚度與濃度計算

由比爾-朗伯定律可知，特定波長處的透射率為

(V)=(t/0)=exp(-(,)d) (1)

式中，t和0分別是透射光和入射光的光強，d是吸收介質的長度。由于尿素水溶液中溶質和溶劑在近紅外區域內均有吸收，本文只考慮溶液（溶質與溶劑）的透射率隨濃度、波數及厚度d的變化，從而建立雙波長透射率與濃度的數學模型以確定液膜的厚度和濃度。吸收率(,)是和的函數，由圖1所獲得的不同濃度下尿素水溶液的吸收率，可在一對選定波長下對(,)進行關于的線性擬合

(v,) =A+Bc,=1,2 (2)

式中，A和B分別是擬合系數。因此，透射率是兩個未知參數和d的函數。

計算雙波長(波長1和波長2)透射率的對數比值可得

由此可見，1/2僅是尿素水溶液吸收率的函數，可轉換為尿素水溶液濃度的函數

因此，溶液濃度可由比值1/2和擬合的系數A和B推出

將式(5)計算得到的值代入式(1)可得

1.2 靈敏度分析

為了實現對尿素水溶液液膜厚度和濃度高精度及高靈敏度的同步測量，需要對雙波長進行優選分析。對尿素水溶液吸收率求關于濃度的導數(d/d)，如圖2所示。選取d/d的極值點與零點所對應的波長位置進行組合，即最小值所在波長位置(7040.75 cm-1)，最大值所在波長位置(6719.75 cm-1)，以及兩個零點所對應的波長位置(6613.25 cm-1和6761.00 cm-1)，這4個波長位置在不同濃度下的吸收率及其關于濃度的線性擬合曲線如圖3所示。本文對4個波長位置進行組合，形成3組波長對(7040.75 cm-1與6613.25 cm-1，7040.75 cm-1與6719.75 cm-1及7040.75 cm-1與6761.00 cm-1)，并對每組波長對進行靈敏度分析。

本文所建立的數學模型對濃度測量的靈敏度取決于1/2對濃度的導數(d1/2/d)。由式(4)可得

圖4是所選取的3對波長組合在不同濃度下的d1/2/d值。如圖所示，3對波長d1/2/d的值均隨著濃度的增加而變大，即隨著濃度的增加，本文所建立的雙波長測量模型的靈敏度增加。在不同的尿素水溶液濃度下，7040.75 cm-1與6761.00 cm-1的波長組合所對應的d1/2/d的值均為最大，可見該對波長組合對濃度測量的靈敏度最高。

2 最優波長對選擇

為獲得雙波長尿素水溶液液膜厚度和濃度反演數學模型的測量精度，本文采用內部交叉驗證對1.2節中所提到的3對波長組合進行研究，從而實現最優波長對的選擇。模型內部交叉驗證是指每個參與建模的樣品都要被其他樣品建立的模型預測一次[20]。本文所有樣本點是不同濃度下(質量分數為5%～50%)每對波長所對應的透射率，依次剔除兩波長位置在選定濃度下的透射率，用其余濃度所對應的透射率按照式(2)獲得吸收率對應的擬合系數A和B，并根據本文建立的數學模型來預測所選定濃度處尿素水溶液液膜的厚度值和濃度值，并與厚度和濃度的真實值進行比較。

通過比較3對波長組合交叉驗證所獲得的預測值與真實值，可獲得本文所提出新方法的測量誤差。如圖5所示，3對波長組合獲得的預測值和真實值均吻合良好。其中，波長組合(7040.75 cm-1與 6761.00 cm-1)計算所得的尿素水溶液液膜參數與真實值相對誤差最小，液膜厚度的最大相對誤差為1.40%，濃度最大相對誤差為5.76%。可見，該對波長組合為最優波長對，這與本文1.2節中的靈敏度分析結果一致。

3 結 論

本文提出了一種基于激光光譜法同步測量尿素水溶液液膜參數的新方法，通過靈敏度分析以及內部交叉驗證選取了一對最優波長組合(7040.75 cm-1與 6761.00 cm-1)，可實現對尿素水溶液液膜厚度和濃度高精度、高靈敏度的測量。由靈敏度分析可知，該對波長組合對濃度測量的靈敏度最高；由內部交叉驗證可得，該對波長組合的液膜厚度和濃度預測值與真實值相對誤差最小，因此，本文將該波長組合選定為最優波長對。由于尿素水溶液的吸收系數比純水小，結合前期工作，本文提出的新方法對尿素水溶液液膜厚度測量的有效范圍可覆蓋5～1000mm，而厚度測量的上限則需要通過實驗進行進一步驗證。進一步工作還將把本文提出的新方法應用到尿素水溶液液膜的動態變化測量中，并對其他水溶液液膜進行研究。

符 號 說 明

Ai,Bi——吸收率線性擬合系數 c——溶液濃度（質量分數），% dl——液膜厚度，cm It,I0——分別為透射光、入射光的光強 k——吸收率，cm-1 v——波數，cm-1 t——透射率 下角標 t——透射光 0——入射光

References

[1] Pautsch A G, Shedd T A, Nellis G F. Thickness measurements of the thin film in spray evaporative cooling [J]., 2004, 1: 70-76

[2] Mawhinney J R, Richardson J K. A review of water mist fire suppression research and development [J]., 1996, 33: 54-90

[3] Schuseberger P. Heat transfer in horizontal tube falling film evaporators//IDA World Congress on Desalination and Water Reuse[C]. Atlantis, The Palm-Dubai, 2009

[4] Gieshoff J, Pfeifer M, Schafer-Sindlinger A, Spurk P C, Garr G, Leprince T, Crocker M. Advanced urea SCR catalysts for automotive applications [J]., 2001: 2001010514

[5] Guo Bin (郭斌), Li Huixiong(李會雄), Guo Dupeng (郭篤鵬). Measurement and analysis of the thickness distributions of falling liquid films around horizontal tubes [J].(工程熱物理學報), 2011, 32:71-74

[6] Cui Jie (崔潔), Chen Xueli (陳雪莉), Wang Qingli(王清立), Wang Fuchen(王輔臣), Gong Xin(龔欣). Double-parallel conductance probe for measuring thickness of liquid film in new-type cyclone separator [J].(化工學報), 2009, 60(6): 1487-1493

[7] Zhang Feng(張鋒), Geng Jiao(耿皎), Ma Shaoling(馬少玲), Zhang Zhibing(張志炳). Measurement of the falling film thickness by capacitance method [J].（測量與設備), 2005, 4: 33-34

[8] Wang R, Lee B A, Lee J S, Kim K Y, Kim S. Analytical estimation of liquid film thickness in two-phase annular flow using electrical resistance measurement [J]., 2012, 36: 2833-2840

[9] Li Weidong (李衛東), Yang Jian (楊健), Li Rongxian (李榮先), Zhou Fangde (周方德), Experimental study and statistical analysis of interfacial waves in stratified gas-liquid flow [J].()(化工學報), 2001, 52(6): 495-498

[10] Zhang Huishu (張會書), Yuan Xigang(袁希鋼), Kalbassi Mohammad Ali. Laser induced fluorescence technique for measuring liquid distribution in structured packing [J].(化工學報), 2014, 65(9): 3331-3339

[11] Mouza A A, Vlachos N A, Paras S V, Karabelas A J. Measurement of liquid film thickness using a light absorption method [J]., 2000, 28:355-359

[12] Hurlburt E T, Newell T A. Optical measurement of liquid film thickness and wave velocity in liquid film flows [J]., 1996, 21: 357-362

[13] Yang H, Greszik D, Wlokas I, Dreier T, Schulz C. Tunable diode laser absorption sensor for the simultaneous measurement of water film thickness, liquid- and vapor-phase temperature [J]., 2011, 104:21-27

[14] Greszik D, Yang H, Dreier T, Schulz C. Laser-based diagnostics for the measurement of liquid water film thickness [J].., 2011, 50: A60-A67

[15] Greszik D, Yang H, Dreier T, Schulz C. Measurements of water film thickness by laser-induced fluorescence and Raman imaging [J]., 2010, 102:123-132

[16] Guo Nairong(過乃蓉), Luo Qingyao(羅慶堯), Zeng Yun’e(曾云鶚). Dual-wavelength spectrophotometer [J].(冶金分析), 1982, 17: 50-57

[17] Hale G M, Querry M R. Optical constants of water in the 200-nm to 200-mm wavelength region [J]., 1973, 12: 555-563

[18] Yang H, Greszik D, Dreier T, Schulz C. Simultaneous measurement of liquid water film thickness and vapor temperature using near-infrared tunable diode laser spectroscopy [J]., 2010, 99: 385-390

[19] Halbout J M, Blit S, Donadson W, Tang C L. Efficient phase-matched second-harmonic generation and sum-frequency mixing in urea [J]., 1979, 15: 1176-1180

[20] Yan Yanlu(嚴衍祿), Chen Bin(陳斌), Zhu Dazhou(朱大洲). Near Infrared Spectroscopy—Principles, Technologies and Applications(近紅外光譜分析的原理、技術與應用)[M]. Beijing：China Light Industry Press，2013: 103-105

Simultaneous measurement of film thickness and concentration of urea-water solution by laser spectroscopy

YANG Huinan，GUO Xiaolong，YANG Bin，SU Mingxu，CAI Xiaoshu

(School of Energy and Power Engineering, Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

Measurement of film thickness and concentration of urea-water solution is very important in many practical applications. However, these two important parameters can only be determined separately by conventional measurement methods. Here, a novel online measurement method by laser spectroscopy was developed to obtain the film thickness and concentration of urea-water solution simultaneously. By sensitivity analysis and cross-validation, an optimal wavelength pair (7040.75 cm-1and 6761.00 cm-1) was chosen to determine film thickness and concentration. As shown from cross-validation, uncertainties of the developed method were smaller than 1.40% for thickness measurement and 5.76% for concentration measurement.

solution; measurement; film; concentration; thickness

2014-07-03.

YANG Bin, yangbin@usst.edu.cn

10.11949/j.issn.0438-1157.20141006

TK 31

A

0438—1157（2015）02—0759—05

國家自然科學基金項目（51306123）；博士點基金聯合資助課題新教師類 (20133120120008)；上海高校青年教師培養資助計劃項目。

2014-07-03收到初稿，2014-11-25收到修改稿。

聯系人：楊斌。第一作者：楊薈楠（1983—），女，博士，講師。

supported by the National Natural Science Foundation of China (51306123), the Joint Specialized Research Found for the Doctoral Program of Higher Education (20133120120008) and the Foundation for Training Young Teacher in University of Shanghai.