靜態(tài)糧堆冷卻與干燥通風(fēng)溫濕度模擬研究

張曉靜,王遠成,潘鈺, 尉堯方, 高帥

(山東建筑大學(xué) 教育部可再生能源建筑利用技術(shù)實驗室, 山東 濟南 250101)

?

靜態(tài)糧堆冷卻與干燥通風(fēng)溫濕度模擬研究

張曉靜,王遠成*,潘鈺, 尉堯方, 高帥

(山東建筑大學(xué) 教育部可再生能源建筑利用技術(shù)實驗室, 山東 濟南 250101)

小麥儲藏過程中，溫度和水分是安全儲糧的關(guān)鍵因素，通風(fēng)過程中小麥堆的溫度和水分的模擬研究對于安全儲糧具有指導(dǎo)意義。文章基于局部熱質(zhì)平衡原理采用Fortran程序建立了一套新的數(shù)學(xué)模型，圍繞靜態(tài)糧堆冷卻與干燥通風(fēng)過程，通過對不同通風(fēng)情況下糧堆溫度和水分變化的模擬研究，闡明了糧堆通風(fēng)過程中溫度和水分變化的一般規(guī)律。結(jié)果表明:對于小麥堆的干燥過程，小麥堆的溫度先升高達到峰值26 ℃后又逐漸趨近于進風(fēng)溫度，隨著通風(fēng)時間的增加，小麥堆水分普遍降低；對于小麥的冷卻過程，小麥堆的溫度隨著進風(fēng)狀態(tài)的改變而改變；通過模擬結(jié)果以及實驗結(jié)果的對比分析，文中模型對溫度的模擬結(jié)果的誤差最大只有2 ℃，測量和預(yù)測水分的最大差值為1.19%。

Fortran程序；儲糧通風(fēng)；傳熱傳質(zhì)

0　引言

糧食儲藏過程中，水分和溫度是特別重要的兩個參數(shù)。環(huán)境溫度的季節(jié)性變化，會導(dǎo)致糧堆內(nèi)熱量傳遞和水分遷移，繼而引起糧堆局部的溫度和水分的升高，從而導(dǎo)致微生物和害蟲的生長，使得糧倉中的糧食變得不安全[1]。糧堆的呼吸作用會消耗糧食的干物質(zhì)，同時呼吸過程會產(chǎn)生熱量和水分，如果糧堆通風(fēng)情況不好，有可能引起糧堆的“出汗”甚至“結(jié)頂”[2]。

經(jīng)過長期的實踐和研究，人們發(fā)現(xiàn)當(dāng)小麥堆溫度低于 15 ℃(小麥安全儲藏的溫度)，小麥堆水分在 12%～12.5%(安全水分)時，可以有效地避免蟲害的發(fā)生，抑制糧堆中生物體的生命活動，延緩儲糧品質(zhì)的劣變[3]。通風(fēng)過程中糧堆內(nèi)部的流動是一個非常復(fù)雜的過程，它與糧堆熱物性參數(shù)和糧食生物特性等多種因素有關(guān)，而且涉及到流體力學(xué)、傳熱學(xué)[4],工業(yè)通風(fēng)[5]、生物學(xué)以及多孔介質(zhì)內(nèi)部的流動傳遞理論[6]。在機械通風(fēng)方面，我國近幾年一些技術(shù)力量較好的國儲庫在研究人員指導(dǎo)下進行了研究和試驗，提出原始水分在16.5%～18.0% 左右的糧食，利用機械通風(fēng)完全可均勻地將其干燥于安全水分之內(nèi)[7]。就倉機械通風(fēng)在小麥儲藏中具有降溫效果顯著、費用較低等特點，在確保儲糧安全方面，發(fā)揮著至關(guān)重要的作用。

為了掌握通風(fēng)過程中儲糧生態(tài)系統(tǒng)變化規(guī)律，儲藏工作者進行了大量的實地測試工作，以研究特定的天氣條件下谷物通風(fēng)系統(tǒng)的性能和空氣流速以及糧堆內(nèi)部溫濕度變化規(guī)律。但是通風(fēng)實驗需要投入較大的人力物力，成本較高，實驗結(jié)果還不具有可重復(fù)性。為了減少工作量，Sharma開發(fā)了仿真模型來描述在可變參數(shù)通風(fēng)下熱與質(zhì)量傳遞現(xiàn)象[8]。由于模型的限制，這些模型不能模擬谷物在整個通風(fēng)段的溫度?；跓崃亢唾|(zhì)量守衡定律，Jia等模擬了就倉通風(fēng)過程中糧倉內(nèi)儲藏小麥的溫度變化，但是沒有模擬小麥堆水分的變化[9]。王平等對平房倉橫向通風(fēng)降溫進行了模擬研究[10]，呂宗旺等以CFD技術(shù)對糧倉進行了降溫模擬[11]，王遠成等對大型房式倉地上籠糧堆溫度和水分變化規(guī)律[12]和圓筒倉內(nèi)自然對流對糧堆熱濕傳遞的影響[13]以及倉儲糧堆內(nèi)熱濕耦合傳遞[14]進行了數(shù)值模擬研究。

以上發(fā)展的模型雖然能比較準(zhǔn)確地模擬通風(fēng)過程中溫度和水分的變化，但是一般都存在計算周期長，占用計算機內(nèi)存大的缺陷。而且在設(shè)定初始條件的時候，大都需要知道所模擬糧種的導(dǎo)熱系數(shù)等參數(shù)，而糧種的這些熱物性參數(shù)受糧食水分以及糧溫的影響比較大，目前國內(nèi)對這方面的統(tǒng)計數(shù)據(jù)尚不完善。文章通過集總參數(shù)的方法結(jié)合Fortran程序設(shè)計了一套新的數(shù)學(xué)模型。該模型不針對某一特種糧食，而是普遍適用于各種儲糧品種，只需要在原程序的基礎(chǔ)上修改糧食的比熱容即可。文章對儲糧過程中小麥堆的溫度和水分變化進行模擬，獲得在不同通風(fēng)條件下小麥堆的溫度和水分，通過與試驗結(jié)果的比較也證明這個模型的準(zhǔn)確性與可行性。

1　模型的建立

1.1物理模型的建立

文章將糧堆劃分為相等厚度的谷物薄層，這些糧層間的空氣流動的方向如圖1所示，其中，ΔH為空間增量，Ti為進風(fēng)溫度，Wi為進風(fēng)空氣濕度比率，To為出風(fēng)溫度，Tw為小麥堆溫度，Wo為出風(fēng)空氣濕度比率。建立仿真模型的目的是預(yù)測在規(guī)定的時間內(nèi)每一層糧層最終的溫度和水分含量。從上一層排出的空氣狀態(tài)被用作輸入下一層糧層的條件，此過程被重復(fù)，直到一段時間的增量上整個谷物柱上的模擬完成。然后該模型以這段時間增量末期的谷物狀態(tài)為初始條件，開始模擬在進風(fēng)狀態(tài)下下一個時間增量的第一層的小麥狀態(tài)，然后使用和第一次時間增量相同的程序，模擬了在第二個時間增量的整個小麥堆。不斷循環(huán)直到完成小麥的整個通風(fēng)過程。

圖1　糧層間空氣流動方向示意圖

1.2數(shù)學(xué)模型的建立

1.2.1焓濕關(guān)系

Brooke提出焓濕關(guān)系式來計算空氣的屬性[15]。

當(dāng)干球溫度Tdb大于或等于0 ℃時，飽和蒸氣壓Ps1(g/cm2)由式(1)計算為

Ps1=exp[58.8858-12301.69/(491.69+1.8Tdb)-5.16923ln(491.69+1.8Tdb)]

(1)

當(dāng)干球溫度Tdb低于0 ℃時，飽和蒸氣壓Ps2(g/cm2)由式(2)計算為

Ps2=exp[27.6452-11286.1489/(491.69+1.8Tdb)-0.46057ln(491.69+1.8Tdb)]

(2)

在給定濕度下的情況下，氣體壓力為Pv=(RH)PS，其中，RH為進風(fēng)的相對濕度。

空氣的濕度比率W=0.6219Pv/(Patm-Pv)，Patm是標(biāo)準(zhǔn)大氣壓，g/cm2。

空氣的比體積vsa=0.007573(491.69+1.8Tdb)(1+1.6055W)。

空氣的質(zhì)量流速ma=Qa/vsa；Qa是空氣的體積流量，m3/h?？諝獾谋葻醕pa=0.2405+0.44W，自由水分的蒸發(fā)潛熱ls=597.768-0.56983T， T是糧堆的溫度，℃

空氣焓h由ASHRAE基本手冊中的式(3)計算為[16]h=4.267+0.240Tdb+W(597.386+0.444Tdb)

(3)

1.2.2小麥基本屬性

Othmer等[17]推導(dǎo)出了小麥達到平衡時空氣的相對濕度RHe(小數(shù)) 由式(4)計算為RHe=exp[(l/ls)ln(2.04816Ps)+c]/2.04816Ps

(4)

式中：l是小麥中水分的蒸發(fā)潛熱，kcal/kg水;l/ls=1+23exp(-0.470 M);c=-3.34×104M-4.0;M是小麥的水分，干基百分比l/ls為1+23exp(-0.40M)。

Viravanichai[18]研究了熱量的以下關(guān)系，硬紅春小麥在溫度范圍為-33.5 至21.8 ℃的小麥的比熱cpw(kcal/kg·℃)值由式(5)～(9)表示為

cpw1=0.250+0.00743M(-33.5 ～-24.1 ℃)

(5)

cpw2=0.279+0.00580M(-21.4 ～-10.8 ℃)

(6)

cpw3=0.297+0.00708M(-10.8 ～3.6 ℃)

(7)

cpw4=0.245+0.01058M(0.6 ～8.9 ℃)

(8)

cpw5=0.273+0.00933M(8.9 ～21.8 ℃)

(9)

1.2.3熱濕傳遞公式

作為中國改革開放的前沿，云南具有優(yōu)越的區(qū)位優(yōu)勢，是“一帶一路”、孟中印緬經(jīng)濟走廊、中國—中南半島經(jīng)濟走廊的重要節(jié)點，是中國唯一可以同時從陸上溝通東南亞、南亞的省，并可以通過中東聯(lián)接歐洲、非洲，開拓西向貿(mào)易通道的省份。獨特的區(qū)位優(yōu)勢，決定了云南將在“一帶一路”建設(shè)中發(fā)揮不可替代的作用。

2.3.1熱平衡公式

(1) 根據(jù)局部熱平衡，可以得到式(10)計算為

maΔtcpa(Ti-To)+mgcpw(Tw-To)+(hi-ho)maΔt l/ls=0

(10)

式中：ma為空氣的質(zhì)量流速，kg干空氣/h；mg為糧層中的小麥質(zhì)量，kg；hi、ho分別為進風(fēng)的焓，kcal/kg干空氣。

對于干燥過程，可變形為

空氣損失的熱=糧堆獲得的熱+蒸發(fā)潛熱

(2) 水分平衡公式

根據(jù)質(zhì)量守恒原理，可以得到式(11)計算為

(Wi-Wo)100maΔt+(Mo-Mi)md=0

(11)

式中：md是糧層中的固體質(zhì)量，kg；Mo是在時間增量Δt內(nèi)模擬的小麥堆水分，干基百分比；Mi是小麥堆的初始水分，干基百分比；Δt是時間增量，h。

對于干燥過程，可變形為

空氣獲得的水分=糧堆失去的水分

文中模型以空氣逐時的溫度和相對濕度、糧堆的深度、糧層的厚度、空氣流速、谷物的初始水分和初始溫度為初始輸入數(shù)據(jù)，各種基礎(chǔ)屬性由上述公式計算得到。

比較Pvi和Pg、Tw和Ti的關(guān)系，存在下面四種可能：

加熱干燥小麥 Ti>Tw,Pvi

冷卻加濕小麥 Ti>Tw,Pvi>Pg

冷卻通風(fēng)小麥 Ti

加熱加濕小麥 Ti>Tw,Pvi>Pg

在以上四個過程中，都達到了熱質(zhì)平衡。為了發(fā)展數(shù)學(xué)模型，文中模型假設(shè)在空氣溫度下首先發(fā)生水分傳遞，然后發(fā)生熱量傳遞。假設(shè)達到平衡時，小麥堆的出風(fēng)溫度和小麥堆的溫度是相同的。同樣的，小麥堆的出風(fēng)空氣濕度比率和小麥堆的比濕達到平衡。而實際在小麥堆通風(fēng)過程中，溫度的變化對小麥堆的影響比水分變化快。因此，當(dāng)已經(jīng)達到溫度平衡時，水分平衡還沒有達到。鄰近通風(fēng)結(jié)束階段時，小麥堆的水分和溫度與空氣都達到平衡。

當(dāng)水分傳遞未達到平衡時，出風(fēng)空氣濕度比率可以由式(12)表示為

(12)

R(R以百分比表示)定義為：糧層內(nèi)實際傳遞的水分與在平衡狀態(tài)下傳遞水分之比。

干燥和調(diào)質(zhì)實驗是在溫度為6.8～25.0 ℃的范圍內(nèi)進行的，通過反復(fù)的實驗，R值應(yīng)為80。在-19.0～0 ℃的溫度范圍調(diào)質(zhì)測試使用的R為經(jīng)驗數(shù)值。

R值在小麥堆糧層厚度和空氣流速一定的情況下取決于空氣和小麥的狀態(tài)。

1.3數(shù)值方法和初始條件

文章模擬了Sharma[8]在實驗室做的實驗工況，得出了120cm深的糧堆干燥通風(fēng)和冷卻通風(fēng)時的溫度和濕度,實驗條件見表1。

表1　用于通風(fēng)模型的變量

2　通風(fēng)實驗

通風(fēng)實驗的原理圖如圖2所示。實驗設(shè)備可以被分成三部分：儲存糧食的垂直通風(fēng)柱，調(diào)節(jié)空氣的空氣室，調(diào)節(jié)穿過糧食的空氣的空氣分配室。下面是這個設(shè)備的簡介，更多的細節(jié)Sharma已經(jīng)給出[8]。

圖2　通風(fēng)實驗原理圖

通風(fēng)測試使用的小麥放在由丙烯腈-丁二烯-苯乙烯(ABS)制造的直徑為10.16cm的觀察管內(nèi)。觀察管的高度為24.13cm，糧層的厚度為15.2cm。沿著觀察管設(shè)置了6個直徑為1.27cm的孔洞，通過孔洞里抽樣檢測小麥，沿著觀察管軸線分布有熱電偶。Brooks玻璃管流量計(型號1110)安裝在通風(fēng)柱和排風(fēng)支管之間的流動管線上，用于測量空氣流速。系統(tǒng)采用機械通風(fēng)，處理過的空氣在風(fēng)機的作用下經(jīng)過空氣室穿過全孔板對糧食進行通風(fēng)。

3　模擬結(jié)果分析

將文章利用程序計算出來的結(jié)果(程序值)與實驗值和Sharma的模擬結(jié)果(模擬值)[8]作比較，得到以下圖形。

3.1工況1下模擬結(jié)果分析

工況1是在進風(fēng)參數(shù)恒定的狀態(tài)下進行模擬的，模擬結(jié)果如圖3、4所示。

由模擬結(jié)果可以發(fā)現(xiàn)，隨著通風(fēng)時間的增加，小麥堆溫度逐漸升高，達到峰值后逐漸趨于進風(fēng)溫度(如圖3(a)所示)，在解吸濕過程中小麥釋放潛熱，使得底層小麥堆的溫度在通風(fēng)的前幾個小時里高于測量溫度而在通風(fēng)后期預(yù)測溫度又低于測量溫度，這個差值最大可達到2 ℃。同時由于局部熱平衡的作用，隨著糧層高度的增加，溫度變化越來越緩和(如圖3(b)所示)。通風(fēng)15h時，糧層內(nèi)水分出現(xiàn)一個峰值(如圖4(a)所示)，隨著通風(fēng)時間的增加，峰值小幅降低且峰值后水分含量普遍降低(如圖4(b)所示)，說明當(dāng)溫度達到平衡時，水分尚未達到平衡。

3.2工況2下模擬結(jié)果分析

工況2是在進風(fēng)參數(shù)隨時間不斷變化的狀態(tài)下進行模擬的，模擬結(jié)果如圖5、6、7所示。

圖3　小麥堆各層溫度隨時間變化圖(a)小麥堆第1層；(b)小麥堆第5層

圖4　不同時刻下小麥堆各層水分變化圖(a)15 h ；(b)25 h

圖5　不同R值對小麥堆水分含量的影響圖(80 h)

低R值表征水分傳遞較少，即從空氣和小麥中帶走少量顯熱。而高R值表征水分傳遞較多，即從空氣和小麥中帶走大量顯熱。由模擬結(jié)果可知最終模擬R值等于80最切合實際(如圖5所示)。通風(fēng)前期，第一糧層的預(yù)測溫度比測量溫度要高。這個差值最大到2 ℃(如圖6(a)所示)。這是因為當(dāng)空氣流過糧層時，它的濕度趨向于平衡濕度。因此，在較上面的糧層里，空氣的濕度與平衡時的濕度差值很小。這就會使得達到平衡時水分傳遞的實際速率比最大傳遞速率的80%大，即上層糧層的R值應(yīng)該比80大，這也是通風(fēng)前期預(yù)測水分低于測量水分(如圖7所示)的原因。模型預(yù)測的小麥堆溫度取決于在吸附過程中空氣潛熱的增加。在通風(fēng)的前幾個小時內(nèi)，最高可傳送水分的百分比是最低的。因此，R值應(yīng)該比80小。這可能是在通風(fēng)初期的幾個小時里預(yù)測溫度比測量溫度低(如圖6(b)所示)的原因。

圖6　小麥堆各層溫度隨時間變化圖(a)小麥堆第1層；(b)小麥堆第7層

圖7　80 h后小麥堆各層水分變化圖

在這些實驗中測量和預(yù)測水分的最大差值為1.19%。產(chǎn)生差值的原因有可能和以下因素有關(guān)：(1)由于在每個實驗中小麥和空氣的狀態(tài)不同，取相同的R值是不準(zhǔn)確的；(2)R值的變化預(yù)示著水分含量的變化；(3)糧層的邊緣效應(yīng)，進風(fēng)相對濕度和小麥堆水分含量的測量誤差。

4　結(jié)論

通過上述研究可知：

(1) 儲糧通風(fēng)過程中溫度與濕度是相互耦合的。

(2) 對于小麥堆的干燥過程，小麥堆的溫度先升高達到峰值后(26 ℃)又逐漸降低并趨近于進風(fēng)溫度(22 ℃)，隨著通風(fēng)時間的增加，小麥堆水分峰值降低且峰值后水分含量普遍降低；對于小麥的冷卻過程，小麥堆的溫度隨著進風(fēng)狀態(tài)的改變而改變，水分含量與實驗結(jié)果比較吻和。

(3) 通過與前人模擬結(jié)果以及實驗結(jié)果的比較發(fā)現(xiàn)，文中模型具有更高的準(zhǔn)確性，溫度模擬結(jié)果的誤差最大只有2 ℃，測量和預(yù)測水分的最大差值為1.19%。文中模型在局部平衡的基礎(chǔ)上首次提出了R值，使水分傳遞更接近于實際情況。

(4) 通過與測量結(jié)果的比較，驗證了該模型的準(zhǔn)確性與可行性。該模型的突出優(yōu)點就是占用內(nèi)存小，運行速度快，操作簡單，普適性強，與之前的數(shù)學(xué)模型相比，有更好的應(yīng)用價值。

[1]楊廣靖,任云虹,賈金，等. 我國糧食儲藏的現(xiàn)狀及發(fā)展趨勢[J].糧食加工，2012， 37(1):60-63.

[2]王遠成, 武傳欣, 段海峰. 糧食儲藏過程中糧堆內(nèi)部熱濕耦合傳遞過程對糧食品質(zhì)變化的影響研究[C].中國糧油儲藏學(xué)會第六屆年會報告，德宏：中國糧油儲藏學(xué)會,2009.

[3]王遠成, 段海峰, 張來林.就倉通風(fēng)時糧堆內(nèi)部熱濕耦合傳遞過程的數(shù)值預(yù)測[J].河南工業(yè)大學(xué)學(xué)報,2009,30(6):76-77.

[4]陶文銓.數(shù)值傳熱學(xué)[M].西安:交通大學(xué)出版社,1998.

[5]孫一堅，沈恒根. 工業(yè)通風(fēng)[M].北京：中國建筑工業(yè)出版社,2010.

[6]林瑞泰.多孔介質(zhì)傳熱傳質(zhì)引論[M].北京:科學(xué)出版社,1995.

[7]湯文良，狄育慧.關(guān)于我國糧倉的通風(fēng)技術(shù)及低溫儲藏技術(shù)[J].潔凈與空調(diào)技術(shù),2013,3(1):44-46.

[8]SharmaS.C. .Simulationofheatandmasstransferduringaerationofwheatandrapeseedbu1ks[J].UnpublishedM.Sc.Thesis,UniversityofManitoba,Winnipeg, 1973(31):104-105.

[9]JiaC.C.，SunD.W.，CaoC.W..Computersimulationoftemperaturechangesinawheatstoragebin[J]．JournalofStoredProductsResearch. 2001, 37( 1) : 165 -177.

[10]王平，周焰，曹陽，等.平房倉橫向通風(fēng)降溫技術(shù)研究[J]. 糧油倉儲科技通訊2011( 2)：20.

[11]呂宗旺，馮黎明，孫福艷.糧倉通風(fēng)CFD技術(shù)研究[J].河南工業(yè)大學(xué)學(xué)報，2013,34(6):105-109.

[12]王遠成,張忠杰,吳子丹，等.計算流體力學(xué)技術(shù)在小麥儲藏中的應(yīng)用[J].中國糧油學(xué)報，2008,27(5):86-90.

[13]王遠成, 亓偉, 張中濤,等. 圓筒倉內(nèi)自然對流對糧堆熱濕傳遞的影響研究[J].水動力學(xué)研究與進展A,2015,25(4):487-496.

[14]王遠成,白忠權(quán),張中濤，等.倉儲糧堆內(nèi)熱濕耦合傳遞的數(shù)值模擬[J].中國糧油學(xué)報,2015,3(11):64-66.

[15]BrookerD.B..Mathematicalmodelofthepychometrics[J].Chart.Trans.ASAE,1967,11( 10): 558-560.

[16]Ashrae.HandbookofFundamentals[M].NewYork:AmericanSocietyofHeating,RefrigeratingandAir-conditioningEngineer,1972.

[17]OthmerD.F,HuangH. .Correlatingvapourpressureandlatentheatdata[J],Ind.Eng.Chem. 1940,32: 841-846.

[18]ViravanichaiS..EffectofMoisturecontentandtemperatureonspecificheatofwheat[D].Winnipeg:Vniverstityofmanitoba, 1971.

(學(xué)科責(zé)編：吳芹)

Simulation of temperature and moisture in static grainstorageduringcoolanddryventilation

Zhang Xiaojing,Wang Yuancheng*, Pan Yu, et al.

(Key Laboratory of Renewable Energy Utilization Technologies in Building of the National Education Ministry, Shandong Jianzhu University, Jinan 250101, China)

Thetemperatureandmoisturecontentiscriticaltoachievingsafewheatstorage.Theresearchershavedesignedavarietyofmathematicalmodelstopredictthetemperatureandmoisturecontentofwheatduringventilation,whichgenerallytakesalotofmoneyandtime.TheproblemcanbesolvedeffectivelybyanewmathematicalmodelbasedonlocalheatandmassbalanceprincipleandFortranprogram.Itmakesthesimulationresearchofgraintemperatureandmoistureindifferentventilationconditionsandexploresthegenerallawofgraintemperatureandmoisturevariation.Theresultsshow:thetemperatureofwheatrisespeakfirstandthenapproachesthetemperatureofincomingairgraduallyduringthedryingventilation,andthemoistureofwheatdecreasesgenerallywiththeincreaseofventilationtime;thetemperatureofthewheatchangeswiththeconditionofincomingairduringcoolingventilation.Themodelhashigheraccuracycomparedwiththeprevioussimulationresultsandexperimentalresults,anditexists2 ℃errorintemperatureand1.19%errorinmoistercontentonly.Themodelisuniversalandpracticalfordifferentkindsofgrainanditisverysignificantfortheguidanceofwheatstorage.

Fortranprogramcalculation;ventilationofwheatstorage;heatandmasstransfer

2015-11-06

國家自然基金項目(51276102)；國家糧食公益專項項目(201313001)；國家糧食公益專項項目(2015449-001-03)

張曉靜(1990-)，女，在讀研究生，主要從事多孔介質(zhì)熱質(zhì)傳遞等方面的研究.E-mail:930065394@qq.com

*:王遠成(1963-)，男,教授,博士，主要從事復(fù)雜系統(tǒng)的傳熱傳質(zhì)等方面的研究.E-mail:wycjn1@163.com

1673-7644(2016)01-0047-06

TS205

A