李 偉 洪厚勝,2 朱曼利 孫宏韜
(1. 南京工業大學生物與制藥工程學院,江蘇 南京 211816;2. 南京匯科生物工程設備有限公司,江蘇 南京 210009;3. 南京工業大學化學與分子工程學院,江蘇 南京 211816)
機械蒸汽再壓縮系統的流程模擬及性能研究
李 偉1洪厚勝1,2朱曼利3孫宏韜1
(1. 南京工業大學生物與制藥工程學院,江蘇 南京 211816;2. 南京匯科生物工程設備有限公司,江蘇 南京 210009;3. 南京工業大學化學與分子工程學院,江蘇 南京 211816)
為使得整個機械蒸汽再壓縮(MVR)蒸發系統的穩定性、節能性更好,借助Aspen Plus軟件學習版,根據MVR運行原理,構建了MVR性能分析模型,并通過改變MVR節點的參數,模擬研究了蒸發量、補充新鮮蒸汽量與進料溫度、蒸發壓強的關系;蒸發量、沸點與進料濃度(NaCl質量分數)、蒸發壓強的關系;蒸發器換熱量、COP與蒸發壓強、壓縮機壓力升的關系。通過分析模擬結果得出:適當減小壓縮比,可以提高蒸發系統的COP;原料液應該加熱到沸點后,再進入蒸發器中進行換熱;從節能效果看,MVR蒸發系統更適合在低溫低壓下運行;在蒸發前,應對濃度較大的原料液進行適當的稀釋。
機械蒸汽再壓縮系統;Aspen Plus;性能研究;節能
Abstract:Mechanical vapor recompression (MVR) system has been more and more popular because of its good energy-saving effect and small footprint, but how to make the stability and energy-saving of whole MVR evaporation system are better still need to further research on its performance. With the help of Aspen Plus software version, the MVR performance analysis model was set according to the operation principle of MVR to study the relationship between the amount of evaporation, the amount of fresh steam and feed temperature, evaporation pressure, and the relationship between evaporation and boiling temperature and feed concentration (NaCl mass fraction) and evaporation pressure, the relationship between heat transfer rate of evaporator, COP and evaporation pressure, pressure increase of compressor by changing the parameters of MVR node. Make conclusions by analyzing the simulation results: the COP can be improved by reducing the compression ratio properly; the feed should be heated to the boiling point temperature and then into the evaporator for heat transfer; MVR evaporation system is more suitable for operation at low temperature and low pressure from the view of energy saving; the feed with high concentration should be dilute before evaporation.
Keywords: mechanical vapor recompression system; Aspen Plus; performance analysis; energy saving
蒸發濃縮作為基本的化工單元操作,早已經廣泛應用于化工、食品、制藥、原子能等領域中[1-2]。但蒸發作為耗能大戶也一直面臨著新技術的改革,機械蒸汽再壓縮系統是繼單效蒸發、多效蒸發、TVR之后第四代蒸發技術,相比于前三代,MVR系統有著節能效果明顯、蒸發溫度低、占地面積小、自動化程度高等諸多優點[3-5]。已經廣泛應用于海水淡化[6-8]、廢水處理[9-10]、中藥濃縮[11-12]、固體干燥[13]等領域。
MVR蒸發系統工藝流程見圖1,其工作原理[14-16]:原料液可以利用蒸汽冷凝水的余熱以及新鮮蒸汽的熱量進行預熱,使其達到飽和溫度,然后進入蒸發器中與殼程中的蒸汽進行換熱蒸發,生成的二次蒸汽經過氣液分離器進行分離提純后,進入蒸汽壓縮機,經由蒸汽壓縮機壓縮后的高品位的二次蒸汽經過冷卻水消除過熱度后,作為加熱物料可以替代新鮮的生蒸汽進入到蒸發器的殼程中,冷凝換熱。同時,由于蒸汽壓縮機性能的問題,通常需要通入一定量的新鮮蒸汽來保持整個系統的穩定換熱,該系統完全利用了二次蒸汽所蘊含的潛能,并且壓縮機的抽吸作用會給系統產生一定的真空度,省去了前三代蒸發器中的用于冷凝二次蒸汽的冷凝器。
機械蒸汽再壓縮系統主要有蒸汽壓縮機、蒸發器、氣液分離器,預熱器、真空泵以及其他附屬設備[17-18]。在該系統中,蒸汽壓縮機是最重要的設備,其工作效率的高低直接關系到MVR系統節能的效果。目前中國用于MVR系統上性能良好的壓縮機主要依賴于進口[19]。根據蒸發量的不同需要,應用最多的是離心式壓縮機和羅茨式壓縮機。蒸發器的種類有很多,目前應用較廣泛、技術成熟的蒸發器是:中央循環管式蒸發器、外熱式蒸發器。由于食品行業物料具有熱敏性的特點,一般選用豎直管降膜式蒸發器較多[20]。水平管降膜蒸發器憑借著高傳熱效率、設備體積小等優點已經越來越受到青睞,但適用于潔凈的物料,易于結垢的物料則有難清理結垢的缺點[21-22],蒸發出來的二次蒸汽含有部分液態水,直接進入蒸汽壓縮機會造成壓縮機的損壞。因此,需要性能良好的氣液分離器,適合MVR系統的氣液分離器主要是懸液分離器和重力沉降分離器;預熱器則主要選用結構簡單、換熱系數高的板式換熱器。

1. 原料罐 2. 一級預熱器 3. 二級預熱器 4. 原料液泵 5. 冷凝水泵 6. 冷凝水罐 7. 濃縮液泵 8. 蒸發器 9. 蒸汽壓縮機
圖1 MVR蒸發系統工藝流程
Figure 1 The process flow of MVR evaporation system
國內外也已經對MVR系統作了大量的模擬研究。顧承真等[23]以水作為工質,利用Aspen Plus軟件建立了分析模型,并研究了循環蒸汽量、補充水的量與進料溫度、冷凝液溫度、蒸汽壓縮比以及蒸發壓強等之間的變化關系,并提出了最佳工藝操作條件。劉立等[24]構建了MVR循環工藝模型,并將實際運行中的蒸汽冷凝放熱量以及電動機功率進行了驗證,得出了所構建模型的可靠性,并以此進行了MVR系統放大分析。Han[25]利用Aspen Plus基于MVR技術的工作原理構建了單級和雙級的零排放的海水淡化系統,并與傳統的蒸發技術進行了能耗對比。這些研究足可以看出Aspen Plus在分析MVR系統上有很大的應用。
MVR蒸發系統的性能受到許多節點的控制,而且每個節點之間的參數都有一定的關聯度,本試驗擬通過流程模擬的方法構建了MVR工藝流程,通過控制進料溫度、進料濃度、蒸發壓強、蒸汽的壓力升等節點的參數,討論其對MVR蒸發系統性能的影響,旨在為離子型溶液的MVR蒸發過程優化提供借鑒,得到能夠模擬MVR蒸發過程的方法,指導MVR系統的工程放大。
1.1 MVR系統理論分析
MVR系統存在物料衡算與能量衡算,物料衡算由式(1)、(2)計算:
F=D+W,
(1)
Fxf=Wxw,
(2)
式中:
F——原料液的質量流量,kg/h;
D——二次蒸汽的質量流量,kg/h;
W——濃縮液的質量流量,kg/h;
xf——原料液中NaCl的質量分數,%;
xw——濃縮液中NaCl的質量分數,%。
能量衡算由式(3)計算:
Dr0=Fc0(t-t0)+Dr1+QL,
(3)
式中:
r0——壓縮后二次蒸汽的汽化熱,kJ/kg;
c0——原料液的比熱容,kJ/(kg·℃);
t——原料液的沸點,℃;
t0——原料液的初始溫度,℃;
r2——二次蒸汽的汽化熱,kJ/kg;
QL——系統的熱量損失,kJ/h。
蒸發器中產生的二次蒸汽經過蒸汽壓縮機壓縮后變為過熱蒸汽,其中,壓縮后的過熱溫度[26]可由式(4)計算:

(4)
式中:
Td——壓縮后過熱蒸汽的溫度,K;
Ts——料液在蒸發壓強下的飽和溫度,K;
k——絕熱系數;
ε——壓縮比;
ηv——容積效率,%;
τ——壓縮機溫度偏差,K。
一定工況下,蒸汽壓縮機在壓縮過程中所消耗的理論功率[27]由式(5)計算:
(5)
式中:
Wth——壓縮機的理論功率,kW;
qvs——二次蒸汽的體積流量,m3/s;
p1——蒸發壓強,kPa;
m——多變系數。
由于壓縮機壓縮二次蒸汽屬于多變過程,同時,壓縮機本身存在機械效率,所以實際運行過程中,壓縮機所消耗的總功率按式(6)計算:

(6)
式中:
WT——壓縮機的實際功率,kW;
ηe——機械效率,%;
ηm——多變效率,%。
過熱蒸汽在進入蒸發器以前,需要對其進行飽和化處理,所需要冷卻水的質量流量由式(7)計算:
Dhw+qshs=(D+qs)hsa,
(7)
式中:
hw——過熱蒸汽的單位比焓,kJ/kg;
qs——冷凝水的質量流量,kg/h;
hs——冷凝水的單位比焓,kJ/kg;
hsa——飽和蒸汽的單位比焓,kJ/kg。
MVR屬于熱泵技術,評價熱泵工作過程性能好壞的參數便是轉為高溫高品位的熱源獲得的熱量Q1與為此所消耗的機械功WT之間的比值大小,工程上稱這個比值為供熱系數[28],用符號COP來表示:

(8)
式中:
COP——供熱系數;
Q1——蒸發器換熱量,kW。
1.2 MVR系統模型的構建
Aspen Plus 是基于穩態化模擬、優化、靈敏度分析和經濟評價的大型化工流程模擬軟件,主要由三部分組成:物性數據庫、單元操作模塊、系統實現策略,可用于各種操作過程的模擬及從單個操作單元到整個工藝流程的模擬[29]。本文根據MVR工藝原理以及Aspen Plus模擬的特點,選用通過熱物流連接的2個Heater與兩相閃蒸器flash2共同組成蒸發器及氣液分離器,選用compr作為蒸汽壓縮機,選用一個Heater將壓縮后的過熱蒸汽冷凝到飽和蒸汽,整個MVR系統分析模型見圖2[30-31],物料選用一定質量分數的NaCl溶液,考察進料溫度、蒸發壓強、進料濃度、二次蒸汽經過壓縮機壓縮后的溫升以及補充新鮮蒸汽的量對蒸發量、壓縮機功率、COP的影響,為保證模擬結果的合理性。現作以下假設:
(1) 整個MVR系統的熱量損失主要集中在蒸發器的換熱過程,選擇總換熱量的5%作為整個系統的熱量損失。
(2) 以產生一定量的二次蒸汽作為初始條件,忽略啟動時所需要的新鮮蒸汽。
(3) 蒸汽壓縮機的多變效率、機械效率均取為90%。
(4) 為達熱平衡,補充0.5 MPa的新鮮飽和蒸汽。

圖2 MVR系統模擬圖
2.1 蒸發量、補充新鮮蒸汽量與進料溫度、蒸發壓強的關系
選取原料的進料量為5 000 kg/h,初始階段產生的二次蒸汽量為1 000 kg/h,原料液中NaCl的質量分數為10%,選取不同的進料溫度,在不同的蒸發壓力下對蒸發量、補充新鮮蒸汽量進行模擬,結果見圖3、4。

圖3 進料溫度、蒸發壓強與蒸發量的關系
Figure 3 The relation between feed temperature, evapora-tion pressure and evaporation rate

圖4 進料溫度、蒸發壓強與補充新鮮蒸汽量的關系
Figure 4 The relation between feed temperature, evapora-tion pressure and Supplementary fresh vapor
由圖3可知,在同一蒸發壓強下,蒸發量隨著進料溫度的增加成線性增加,在同一進料溫度下,蒸發量隨蒸發壓強的增大而減小。由圖4可知,為達到熱平衡,在同一蒸發壓強下,隨著進料溫度的增加,所需要補充的新鮮蒸汽量呈線性減少,在同一進料溫度下,需補充的新鮮蒸汽量隨著蒸發壓強的增大而逐漸增大。主要原因是:① 原料液在進入蒸發器后需要先加熱到相應壓強下的沸點溫度才可以進行蒸發,在蒸汽換熱量一定的前提下,進入蒸發器的原料液溫度越低,本身所攜帶的熱量越少,則需要加熱到沸點所需要的熱量越多,則用于蒸發的剩余熱量會越少,產生的二次蒸汽量越少,為維持系統一定的蒸發量,所需要補充的新鮮蒸汽的量會相應地減少;② 在一定的進料溫度下,系統的蒸發壓強越低,其所對應的飽和沸點溫度越低,則會有更多的熱量用于蒸發過程,產生的二次蒸汽越多,所補充的新鮮蒸汽量越少。
2.2 蒸發量、沸點與進料濃度(NaCl質量分數)、蒸發壓強的關系
料液的進料量為5 000 kg/h,初始階段產生的二次蒸汽量為1 000 kg/h,進料溫度為飽和溫度,選取不同的進料濃度,在不同的蒸發壓強下對蒸發量、溶液的沸點進行模擬,結果見圖5、6。

圖5 進料濃度、蒸發壓強與蒸發量的關系
Figure 5 The relation between feed concentration, evaporation pressure and evaporation rate

圖6 進料濃度、蒸發壓強與沸點的關系
Figure 6 The relation between feed concentration, evaporation pressure and boiling temperature
由圖5可知,在同一蒸發壓強下,隨著進料濃度的增加,料液的蒸發量逐漸減小,且蒸發壓強為40 kPa時蒸發量的下降速度要小于蒸發壓強為60,80 kPa時的。由圖6可知,隨著進料濃度的增加,沸點呈上升趨勢,在同一進料濃度下,沸點隨著蒸發壓強的增大而升高。主要原因在于在進料溫度與系統的換熱量一定的情況下,進料濃度越高,溶液的沸點越高,溶液達到沸騰需要的熱量會越多,因此,用于蒸發的熱量會相應地減小,產生的二次蒸發量減小。另外,隨著進料濃度的增加蒸發壓強為40 kPa的沸點升速要小于蒸發壓強為60,80 kPa下的。
2.3 壓縮機功率、冷凝器熱負荷與蒸發壓強,蒸汽壓力升的關系
料液的進料量為5 000 kg/h,初始階段產生的二次蒸汽量為1 000 kg/h,進料溫度為飽和溫度,進料濃度為10%,選取二次蒸汽經過壓縮機壓縮后不同的壓力升,在不同的蒸發壓強下對壓縮機功率、冷凝器熱負荷進行模擬,結果見圖7、8。
由圖7、8可知,兩者的變化趨勢差別不大,在一定的蒸汽壓力升下,壓縮機功率、冷凝器熱負荷隨著蒸發壓強的增大而減小,在一定的蒸發壓強下,壓縮機功率、冷凝器熱負荷隨著蒸汽壓力升的增大而增大。主要原因是:① 當蒸汽的壓力升一定、壓縮機的處理量不變時,隨著蒸發壓強的增大,壓縮機所對應的壓縮比會相應減小,壓縮后二次蒸汽的能量會變小,過熱溫度會減小,被冷卻水冷卻成飽和蒸汽所放出的熱量同樣會減小;② 在同一蒸發壓強下,蒸汽的壓力升越大,相對應壓縮機的壓縮比會越大,壓縮后蒸汽的熱量會變大,需要壓縮機做更多的功,壓縮機所消耗的功率越大,則經過壓縮后的過熱蒸汽的溫度越高,為冷凝到飽和狀態,所消耗冷卻水的量也相應地增加。所以,兩者具有相同的趨勢。

圖7 蒸發壓強、壓力升與壓縮機功率的關系Figure 7 The relation between evaporation pressure, pressure increase of vapor and compressor power

圖8 蒸發壓強、壓力升與熱負荷的關系
Figure 8 The relation between evaporation pressure, pressure increase of vapor and thermal load
2.4 蒸發器換熱量、COP與蒸發壓強、蒸汽壓力升的關系
進料量為5 000 kg/h,初始二次蒸汽量為1 000 kg/h,原料液的濃度為10%,沸點進料,熱損失保持換熱總量的5%,在不同的蒸發壓強,不同的蒸汽壓力升下,對蒸發器的換熱量、COP進行模擬,結果見圖9、10。

圖9 蒸發壓強、壓力升與換熱量的關系
Figure 9 The relation between evaporation pressure, pressure increase of vapor and heat exchange

圖10 蒸發壓強、壓力升與COP的關系
Figure 10 The relation between evaporation pressure, pressure increase of vapor andCOP
由圖9可知,在同一蒸汽壓力升下,蒸發器的換熱量隨蒸發壓強的增大而減小;在同一蒸發壓強下,隨著壓力升的減小而增大。主要原因是:① 在一定的壓力升下,蒸發壓強越大,壓縮機相對應的壓縮比越小,其被壓縮后的飽和蒸汽壓越大,相對應的汽化熱越小,如保持產生的二次蒸汽量不變,則蒸發器的換熱量會相應地減小;② 在同一蒸發壓強下,由于二次蒸汽的進料量一定,蒸汽的壓力升越大,二次蒸汽被壓縮后的飽和蒸汽壓越大,相對應的汽化熱值減小,蒸發器的換熱量減小。由圖10可知,在同一壓力升下,隨蒸發壓強的增大,相對應的壓縮比減小,COP值變大,因此適當地增大蒸發壓強,可以降低能耗,提高蒸發系統的COP。在不同的蒸發壓強下,壓力升越大,COP的增幅越大,在不同的壓強下,壓力升為30 kPa時的COP增幅最大,也就是說,壓縮比越小,COP值越大,但由于需要維持一定的蒸發量,因此,壓縮比應該有一個合理的取值范圍。
通過Aspen Plus軟件對影響MVR蒸發系統的幾個節點進行了模擬,主要通過進料溫度、蒸發壓強、蒸發量、補充新鮮蒸汽量、壓縮機壓力升、進料濃度、蒸發器換熱量、冷凝器熱負荷、壓縮機功率、COP性能評價指標來對MVR蒸發系統的性能進行研究。并對模擬結果進行了理論分析,結論如下:
(1) 適當地減小壓縮比,可以提高蒸發系統的COP,減少壓縮機所消耗功率,減少冷凝水的用量。理論上,壓縮比越小,系統越節能、效率越高。但由于系統需要維持一定的蒸發量,因此,壓縮比應該有一個合理的取值范圍。
(2) 為維持系統一定的蒸發量,更好地將系統熱能用于蒸發,減少新鮮蒸汽的消耗量,原料液應該加熱到沸點后,再進入蒸發器中進行換熱。
(3) MVR蒸發系統更適合在低溫低壓下運行,同時需要結合物料的熱敏性來考慮。
(4) 物料的濃度越大,沸點升越大,越不易蒸發,在蒸發前,可以對濃度較大的原料液進行適當的稀釋后再進行蒸發。
[1] 姚玉英, 黃鳳廉, 陳常貴, 等. 化工原理: 上冊[M]. 天津: 天津科學技術出版社, 2012: 290-298.
[2] ZHOU Ya-su, SHI Cheng-jun, DONG Guo-qiang. Analysis of a mechanical vapor recompression wastewater distillation system[J]. Desalination, 2014, 353: 91-97.
[3] ETTOUNEY H. Design of single-effect mechanical vapor compression[J]. Desalination, 2006, 190(1/2/3): 1-15.
[4] SHEN Jiu-bing, XING Zi-wen, ZHANG Kai, et al. Develop-ment of a water-injected twin-screw compressor for mechanical vapor compression desalination systems[J]. Applied Thermal Engineering, 2016, 95: 125-135.
[5] 李偉, 朱曼利, 洪厚勝. 機械蒸汽再壓縮系統的設計及其在食品工業中的應用[J]. 食品與機械, 2016, 32(7): 223-226.
[6] HAN Dong, HE Wei-feng, YUE Chen, et al. Study on desalination of zero-emission system based on mechanical vapor compression [J]. Applied Energy, 2017, 185: 1 490-1 496.
[7] 劉鵬, 王永青. 機械壓汽蒸餾海水淡化技術的研究和發展狀況[J]. 機電技術, 2011, 34(4): 161-165.
[8] 焦冬生. 機械壓汽蒸餾海水淡化系統的可用能分析[J]. 太陽能學報, 2008, 29(10): 1 197-1 203.
[9] LI Yu-long, WU Hong, LIANG Xin-gong, et al. Experimental study of waste concentration by mechanical vapor compression technology[J]. Desalination, 2015, 361: 46-52.
[10] 梁林, 韓東, 彭濤. 機械蒸汽再壓縮硫酸銨廢水處理系統的分析[J]. 化學工程, 2012, 40(8): 74-78.
[11] 王谷洪, 郭亮, 周齊. 淺析MVR濃縮技術在中藥方面的應用[J]. 現代制造, 2015(35): 22-24.
[12] 張功臣. MVR蒸發器的節能特點及其在中藥濃縮中的應用[J]. 現代制造, 2015(8): 25-28.
[13] 楊德明, 葉夢飛, 譚建凱, 等. 機械蒸汽再壓縮(MVR)熱泵技術的應用進展[J]. 常州大學學報: 自然科學版, 2015, 27(1): 76-80.
[14] 梁林, 韓東. 蒸汽機械再壓縮蒸發器的實驗[J]. 化工進展, 2009, 28(S1): 358-360.
[15] YANG Jun-ling, ZHANG Chong, ZHANG Zhen-tao, et al. Study on mechanical vapor recompression system with wet compression single screw compressor[J]. Applied Thermal Engineering, 2016, 103: 205-211.
[16] 石成君, 周亞素, 孫韶, 等. 機械蒸汽再壓縮蒸發技術高鹽度廢水處理系統的性能分析[J]. 水處理技術, 2013, 39(12): 63-68.
[17] LIANG Lin, HAN Dong, MA Ran, et al. Treatment of high-concentration wastewater using double-effect mechanical vapor recompression[J]. Desalination, 2013, 314(4): 139-146.
[18] 李偉, 朱曼利, 洪厚勝. 機械蒸汽再壓縮技術(MVR)研究現狀[J]. 現代化工, 2016(11): 28-31.
[19] 王力威, 莊景發, 楊魯偉, 等. 單螺桿水蒸汽壓縮機驅動的MVR系統性能實驗研究[J]. 中國科學院大學學報, 2015, 32(1): 38-45.
[20] 顧承真, 洪厚勝, 顏旭, 等. 羅茨壓縮機驅動MVR熱泵系統的節點分析[J]. 現代化工, 2014, 34(4): 125-129.
[21] 路慧霞, 馬曉建. 水平管外降膜蒸發的傳熱實驗[J]. 化工進展, 2009, 28(2): 203-205.
[22] 何茂剛, 王小飛, 張穎. 制冷用水平管降膜蒸發器的研究進展及新技術[J]. 化工學報, 2008, 59(S2): 23-28.
[23] 顧承真, 閔兆升, 洪厚勝. 機械蒸汽再壓縮蒸發系統的性能分析[J]. 化工進展, 2014, 33(1): 30-35.
[24] 劉立. 機械蒸汽再壓縮式降膜蒸發系統的設計和性能研究[D]. 天津: 河北工業大學, 2014: 37-44.
[25] HAN Dong. Study on zero-emission desalination system based on mechanical vapor recompression technology[J]. Energy Procedia, 2015, 75: 1 436-1 444.
[26] 蘇春模. 羅茨鼓風機及其使用[M]. 長沙: 中南工業大學出版社, 1999: 66-75.
[27] 劉立, 張繼軍, 劉燕, 等. 機械蒸汽再壓縮式熱泵用于降膜蒸發系統的研究[J]. 現代化工, 2014, 34(9): 128-132.
[28] 鄭聰. 新型熱泵操作方式的研究[D]. 天津: 天津大學, 2009: 19-20.
[29] 張治山, 楊超龍. Aspen Plus在化工中的應用[J]. 廣東化工, 2012, 39(3): 77-78.
[30] 郝冬青, 沙作良, 王彥飛, 等. 低溫多效海水淡化系統的Aspen Plus模擬[J]. 天津科技大學學報, 2011, 26(1): 47-50.
[31] 陳軍. 低溫多效蒸發海水淡化系統工藝流程模擬及優化[D]. 北京: 北京化工大學, 2013: 19-33.
Process simulation and performance study of mechanical vapor recompression system
LI Wei1HONGHong-sheng1,2ZHUMan-li3SUNHong-tao1
(1.CollegeofLifeScienceandPharmaceuticalEngineering,NanjingTechUniversity,Nanjing,Jiangsu211816,China; 2.NanjingHuikeBioengineeringEquipmentCorporation,Nanjing,Jiangsu210009,China; 3.CollegeofChemistryandMolecularEngineering,NanjingTechUniversity,Nanjing,Jiangsu211816,China)
10.13652/j.issn.1003-5788.2017.08.016
李偉,男,南京工業大學在讀碩士研究生。
洪厚勝(1965—),男,南京工業大學教授,博士。 E-mail: hhs@njtech.edu.cn
2017—05—04