楊興林,楊歡歡,張禮松,劉明遠(yuǎn)
(1.江蘇科技大學(xué)能源與動力工程學(xué)院,江蘇鎮(zhèn)江212003)
(2.江西中船航海儀器有限公司,江西九江332000)
(3.中國人民解放軍92337部隊(duì),遼寧大連116023)
輔助冷卻的復(fù)合式地源熱泵系統(tǒng)性能分析
楊興林1,楊歡歡1,張禮松2,劉明遠(yuǎn)3
(1.江蘇科技大學(xué)能源與動力工程學(xué)院,江蘇鎮(zhèn)江212003)
(2.江西中船航海儀器有限公司,江西九江332000)
(3.中國人民解放軍92337部隊(duì),遼寧大連116023)
文中提出兩種輔助冷卻的復(fù)合式地源熱泵系統(tǒng),以江南鎮(zhèn)江地區(qū)一棟典型辦公大樓為例,利用TRNSYS軟件建立復(fù)合式地源熱泵系統(tǒng)仿真模型,并與傳統(tǒng)地源熱泵系統(tǒng)的運(yùn)行性能及能耗進(jìn)行對比分析,結(jié)果表明:設(shè)有冷水機(jī)組的輔助冷卻復(fù)合式地源熱泵系統(tǒng)的土壤溫升更小、運(yùn)行更穩(wěn)定高效、系統(tǒng)能耗更低、具有更好的經(jīng)濟(jì)性.
復(fù)合式地源熱泵系統(tǒng);性能系數(shù);能耗
地源熱泵系統(tǒng)以其環(huán)保、高效節(jié)能、環(huán)境與經(jīng)濟(jì)效益顯著、系統(tǒng)運(yùn)行穩(wěn)定、能源可再生利用、維護(hù)量小等優(yōu)點(diǎn)而得以推廣[1].隨著地源熱泵技術(shù)的應(yīng)用范圍越來越廣泛,系統(tǒng)問題逐漸顯現(xiàn).在夏季通過地埋管換熱器向地下巖土的放熱量與在冬季通過地埋管換熱器從地下巖土的取熱量相差過大,導(dǎo)致地埋管換熱器在運(yùn)行數(shù)年后,地下巖土的平均溫度將逐年升高或降低,從而影響地埋管周圍土壤溫度場的恢復(fù),導(dǎo)致地源熱泵空調(diào)系統(tǒng)性能不斷下降,甚至不能運(yùn)行.因此需要在傳統(tǒng)地源熱泵系統(tǒng)中增加輔助散熱裝置或輔助加熱裝置,平衡地埋管換熱器的換熱負(fù)荷,使地下巖土平均溫度維持在有利于地源熱泵系統(tǒng)高效節(jié)能運(yùn)行的范圍內(nèi),即采用復(fù)合式地源熱泵空調(diào)系統(tǒng)[2-3].
選取江南地區(qū)鎮(zhèn)江的一棟典型辦公大樓為研究對象,共12層,占地面積1 715 m2,建筑總面積10 296 m2.利用TRNSYS[4]軟件計(jì)算全年逐時(shí)動態(tài)負(fù)荷,獲得全年累計(jì)冷負(fù)荷為2 381 845 kW,全年累計(jì)熱負(fù)荷814 894 kW,全年累計(jì)冷熱負(fù)荷比為2.92∶1.針對鎮(zhèn)江地區(qū)的氣候特征與該辦公建筑的負(fù)荷特點(diǎn),設(shè)計(jì)地源熱泵系統(tǒng)時(shí)應(yīng)高度重視熱平衡問題,因此文中提出了兩種輔助冷卻的復(fù)合式地源熱泵系統(tǒng).
第一種是聯(lián)合冷卻塔的復(fù)合式地源熱泵系統(tǒng)[5],該系統(tǒng)是在地源熱泵系統(tǒng)的基礎(chǔ)上添加了輔助散熱設(shè)備——冷卻塔,主要由冷卻塔、地埋管、熱泵機(jī)組3個(gè)部分組成.冷卻塔只用于制冷工況,機(jī)組流出的高溫冷卻水先由冷卻塔冷卻,隨后經(jīng)初步冷卻后的冷卻水通過地埋管換熱器進(jìn)行循環(huán)排熱.冬季工況時(shí),地埋管換熱器直接與機(jī)組中循環(huán)的流體介質(zhì)進(jìn)行換熱.圖1給出了該復(fù)合地源熱泵系統(tǒng)原理.

圖1 地源熱泵與冷卻塔復(fù)合式系統(tǒng)原理Fig.1 Principle of the hybrid system of ground source heat pump and cooling tower
該系統(tǒng)的優(yōu)點(diǎn)是:① 室外濕球溫度較低,冷卻塔可以獨(dú)立運(yùn)行,充分發(fā)揮冷卻塔的優(yōu)勢;②冷卻塔的獨(dú)立運(yùn)行使地埋管周圍能夠及時(shí)散熱,土壤溫度不會升高;③地埋管周圍土壤溫度在供冷季后期時(shí)會升高,從而降低系統(tǒng)的效率,此時(shí)冷卻塔獨(dú)立運(yùn)行可以讓系統(tǒng)整體達(dá)到最優(yōu)[6].
第二種是聯(lián)合冷水機(jī)組的復(fù)合式地源熱泵系統(tǒng)[7],該系統(tǒng)是在地源熱泵的基礎(chǔ)上增加了輔助冷源——冷水機(jī)組和輔助散熱設(shè)備——冷卻塔,其主要由地埋管、熱泵機(jī)組、冷水機(jī)組、冷卻塔等組成.冷水機(jī)組聯(lián)合地源熱泵的復(fù)合式系統(tǒng)中,地源熱泵供冷供熱系統(tǒng)與常規(guī)冷水機(jī)組供冷系統(tǒng)在室內(nèi)側(cè)連接在一起,但兩者獨(dú)立運(yùn)行,在夏季,可同時(shí)運(yùn)行兩系統(tǒng)為室內(nèi)供冷;在冬季,僅采用地源熱泵系統(tǒng)為室內(nèi)供熱.由于管網(wǎng)阻抗穩(wěn)定,因此室外側(cè)系統(tǒng)運(yùn)行也較為穩(wěn)定.系統(tǒng)中地埋管換熱器和冷卻塔均能單獨(dú)運(yùn)行,所以即使地埋管換熱器間歇運(yùn)行仍可實(shí)現(xiàn)空調(diào)系統(tǒng)的連續(xù)運(yùn)行,從而使土壤溫度得到恢復(fù),更有利于系統(tǒng)的節(jié)能.其系統(tǒng)原理如圖2.

圖2 地源熱泵與冷水機(jī)組復(fù)合式系統(tǒng)原理Fig.2 Principle of the hybrid ground source heat pump system and water chiller
該系統(tǒng)的優(yōu)點(diǎn)是:① 輔助冷卻塔與地埋管系統(tǒng)的并聯(lián)連接方式較優(yōu);② 將地埋管和常規(guī)冷卻塔空調(diào)系統(tǒng)獨(dú)立設(shè)置,可以減小熱泵機(jī)組冷凝器設(shè)計(jì)溫度與冷卻塔進(jìn)出水溫度之間的差異;③ 系統(tǒng)中不再需要使用板式換熱器,從而減少了熱量損失和初投資[8].
根據(jù)計(jì)算的建筑負(fù)荷,對系統(tǒng)設(shè)備和組件進(jìn)行設(shè)計(jì)計(jì)算選型[9].地埋管總數(shù)為535個(gè),總長度為51 517.2 m,每米管長換熱量為61.9 W/m.熱泵機(jī)組采用環(huán)保制冷劑為HFC134a的半封閉式雙螺桿熱泵機(jī)組,制冷量為2 756 kW,輸入功率為557 kW,蒸發(fā)器流量520 m3/h,冷凝器流量580 m3/h,制熱量為2 732 kW,輸入功率為550 kW,蒸發(fā)器流量580 m3/h,冷凝器流量520 m3/h,通過與壓縮機(jī)的減載系數(shù)相匹配,可實(shí)現(xiàn)10%~100%內(nèi)能量的無級調(diào)節(jié).選用型號為ZXZ-N300T的冷卻塔,冷卻能力為 1.5×106kcal/h.選用型號為 LSBLX1600T的冷水機(jī)組,制冷量為1 600 kW,功率為199 kW.
本工程是冷負(fù)荷遠(yuǎn)大于熱負(fù)荷的辦公大樓,主要特點(diǎn)是白天使用率比較高,針對這一特點(diǎn)采取不同的控制策略.對于熱泵機(jī)組,按照機(jī)組進(jìn)水溫度控制,在夏季制冷季節(jié),機(jī)組進(jìn)水溫度低于7℃時(shí),關(guān)閉熱泵機(jī)組,機(jī)組進(jìn)水溫度高于12℃,開啟熱泵機(jī)組;在冬季制熱季節(jié),機(jī)組進(jìn)水溫度低于40℃時(shí),關(guān)閉熱泵機(jī)組,機(jī)組進(jìn)水溫度高于45℃,開啟熱泵機(jī)組.對于冷卻塔,文中則采取時(shí)間上的控制方法,夏季下班時(shí)間18:00至第二天上班時(shí)間8:00開啟冷卻塔,開始向土壤進(jìn)行蓄冷.對于冷水機(jī)組,根據(jù)建筑冷熱負(fù)荷比2.92∶1,按照負(fù)荷率50%分別分配給冷水機(jī)組和地源熱泵機(jī)組.
利用TRNSYS軟件對復(fù)合式地源熱泵系統(tǒng)的設(shè)計(jì)和模擬過程中,主要有以下幾個(gè)關(guān)鍵模塊:數(shù)據(jù)導(dǎo)入部件、水-水熱泵機(jī)組、地埋管換熱器、冷卻塔、冷水機(jī)組、水泵、控制系統(tǒng)、圖形顯示部件和結(jié)果輸出部件等模塊.圖3~5分別給出了常規(guī)地源熱泵系統(tǒng)仿真模型、地源熱泵與冷卻塔復(fù)合式系統(tǒng)仿真模型、地源熱泵與冷水機(jī)組復(fù)合式系統(tǒng)仿真模型.由于該辦公大樓為大型建筑,分區(qū)多,末端連接復(fù)雜,控制較多,系統(tǒng)較難收斂,因此文中主要研究該建筑對應(yīng)熱泵系統(tǒng)的長期運(yùn)行參數(shù),故將動態(tài)負(fù)荷輸入至負(fù)荷處理模塊Type682中處理.

圖3 常規(guī)地源熱泵系統(tǒng)仿真模型Fig.3 Model of the conventional ground source heat pump system

圖4 地源熱泵與冷卻塔復(fù)合式系統(tǒng)仿真模型Fig.4 Model of the hybrid system of ground source heat pump and cooling tower

圖5 地源熱泵與冷水機(jī)組復(fù)合式系統(tǒng)仿真模型Fig.5 Model of the hybrid ground source heat pump system and water chiller
分別對常規(guī)地源熱泵系統(tǒng)(Only-GSHP系統(tǒng))、冷卻塔與地源熱泵的復(fù)合式系統(tǒng)(CT-GSHP系統(tǒng))以及冷水機(jī)組和地源熱泵的復(fù)合式系統(tǒng)(WCHP-GSHP系統(tǒng))進(jìn)行為期10 a的模擬,對其運(yùn)行性能及能耗進(jìn)行分析.
4.1 土壤熱平衡分析
表1給出了3種系統(tǒng)的土壤溫升情況,土壤平均溫升由大到小依次為Only-GSHP系統(tǒng)、CT-GSHP系統(tǒng)、WCHP-GSHP系統(tǒng),這是由于一方面冷卻塔輔助散熱;另一方面冷水機(jī)組承擔(dān)了部分負(fù)荷,地埋管換熱器承擔(dān)的熱量越少,導(dǎo)致土壤的溫升越小.地源熱泵系統(tǒng)運(yùn)行10 a后,3個(gè)系統(tǒng)的土壤溫升分別為27.21℃,4.17℃及2.24℃,結(jié)果說明3種系統(tǒng)都存在一定程度的溫升,但地源熱泵加冷卻塔加冷水機(jī)組系統(tǒng)的溫度上升值是可以接受的.地下土壤有吸熱和散熱的能力,若在熱泵的一個(gè)運(yùn)行年內(nèi),從土壤中吸收的熱量與排至土壤的熱量的差值不超過土壤的自動恢復(fù)能力,土壤溫度將不會大范圍的偏離設(shè)置地源熱泵系統(tǒng)前的土壤初始溫度,造成土壤熱堆積或冷坑現(xiàn)象.

表1 3種系統(tǒng)的土壤溫升表Table 1 Soil temperaturesrises with three kinds of systems
4.2 系統(tǒng)運(yùn)行性能分析
從圖6中可以看出,3種系統(tǒng)的制熱工況下平均性能系數(shù)COP增加,但是制冷工況下平均性能系數(shù)COP降低.對于CT-GSHP系統(tǒng),單純考慮系統(tǒng)運(yùn)行COP時(shí),其制熱效率COP最低,這是由于隨著冷卻塔等設(shè)備的增加以及土壤熱量平衡問題導(dǎo)致其能耗增加.WCHP-GSHP系統(tǒng)平均性能COP能保持較高的數(shù)值,這是由于利用冷水分擔(dān)一部分負(fù)荷,且使得土壤處于一定的熱平衡環(huán)境下,促使熱泵機(jī)組處于高效的運(yùn)行環(huán)境;隨著系統(tǒng)的持續(xù)運(yùn)行,系統(tǒng)的平均COP降低.Only-GSHP系統(tǒng)的COP從第1年的2.93運(yùn)行10 a后降至2.53,這是土壤熱不平衡越來越嚴(yán)峻,土壤熱堆積問題嚴(yán)重從而大幅度降低制冷工況下系統(tǒng)的運(yùn)行性能,而CT-GSHP系統(tǒng)維持在2.94左右,WCHP-GSHP系統(tǒng)維持在3.34左右,因此可以看出WCHP-GSHP系統(tǒng)具有較高的COP以及穩(wěn)定性.

圖63種系統(tǒng)的地源熱泵系統(tǒng)COPFig.6 COP of three kinds of ground heat pump systems
4.3 系統(tǒng)能耗分析
文中的系統(tǒng)能耗包括熱泵機(jī)組的能耗、冷卻塔的能耗、冷水機(jī)組的能耗,冷凍水泵的能耗、冷卻水泵的能耗以及冷卻塔泵的能耗,因?yàn)橥粋€(gè)建筑物內(nèi)的風(fēng)機(jī)能耗基本上是相等的,所以在進(jìn)行能耗對比分析時(shí)刻不必考慮風(fēng)機(jī)能耗.
從圖7a)和b)中可以看出地源熱泵系統(tǒng)冬夏季的逐年運(yùn)行能耗分別與熱泵機(jī)組冬夏季逐年運(yùn)行能耗變化規(guī)律相似,但是變化幅度不同,全年的系統(tǒng)能耗與全年的機(jī)組能耗變化規(guī)律相似,從圖7c)中可看出,3種系統(tǒng)平均全年累計(jì)運(yùn)行能耗(E)分別為5.5×106kW·h,4.6×106kW·h和2.8×106kW·h.

圖7 3種地源熱泵系統(tǒng)能耗Fig.7 Energy consumption of three kinds of ground source heat pump systems
4.4 經(jīng)濟(jì)性分析
系統(tǒng)的初投資主要包括地埋管和機(jī)組的費(fèi)用.地埋管的費(fèi)用主要包括鉆孔費(fèi)用和管材費(fèi),平均為100元/m.機(jī)組主要是熱泵機(jī)組、冷卻塔、冷水機(jī)組,其他零部件忽略不計(jì).由于文中3個(gè)系統(tǒng)的地埋管數(shù)量一樣,所以在初投資方面,聯(lián)合冷水機(jī)組運(yùn)行的復(fù)合式地源熱泵系統(tǒng)比只聯(lián)合冷卻塔運(yùn)行的復(fù)合式熱泵系統(tǒng)多出冷水機(jī)組的價(jià)格.
系統(tǒng)的運(yùn)行費(fèi)用主要是機(jī)組和系統(tǒng)的能耗(表2),只聯(lián)合冷卻塔的復(fù)合式地源熱泵系統(tǒng)運(yùn)行費(fèi)用比傳統(tǒng)的節(jié)省16.27%,聯(lián)合冷水機(jī)組的復(fù)合式地源熱泵系統(tǒng)運(yùn)行費(fèi)用比傳統(tǒng)的節(jié)省49.05%的費(fèi)用.

表2 3種系統(tǒng)運(yùn)行費(fèi)用Table 2 Operational costs of three kinds of systems
初投資和運(yùn)行費(fèi)用是兩個(gè)不同性質(zhì)的經(jīng)濟(jì)指標(biāo),不能全面反映各方案經(jīng)濟(jì)性,要全面比較、評價(jià)各方案的經(jīng)濟(jì)性就需要一個(gè)綜合的經(jīng)濟(jì)指標(biāo).文中通過計(jì)算各方案初投資及壽命周期內(nèi)各年的運(yùn)行費(fèi)用,運(yùn)用動態(tài)費(fèi)用年值A(chǔ)C法對各方案的經(jīng)濟(jì)性進(jìn)行比較和分析:

式中:i為銀行利率,取10%;n為使用壽命,各個(gè)方案中地源熱泵系統(tǒng)皆取50 a,鍋爐系統(tǒng)取10 a;C0為初投資,元;C為年運(yùn)行費(fèi)用,元.
表3為3種系統(tǒng)的經(jīng)濟(jì)評價(jià).從表中可看出,聯(lián)合冷水機(jī)組的復(fù)合式地源熱泵系統(tǒng)的動態(tài)費(fèi)用年值最少,比傳統(tǒng)的減少41.52%,最具經(jīng)濟(jì)性.

表3 3種系統(tǒng)經(jīng)濟(jì)評價(jià)Table 3 Economic analysis for three kinds of systems
文中提出了聯(lián)合冷卻塔運(yùn)行的復(fù)合式地源熱泵系統(tǒng)和聯(lián)合冷水機(jī)組運(yùn)行的復(fù)合式地源熱泵系統(tǒng),以江南鎮(zhèn)江地區(qū)一棟典型辦公大樓為例,利用TRNSYS軟件建立復(fù)合式地源熱泵系統(tǒng)仿真模型,并與傳統(tǒng)地源熱泵系統(tǒng)的運(yùn)行性能及能耗進(jìn)行對比分析,結(jié)論如下:
1)系統(tǒng)運(yùn)行10 a后,3種系統(tǒng)的土壤溫升依次為27.21℃,4.17℃及2.24℃,都存在一定的土壤熱失衡現(xiàn)象,尤其是傳統(tǒng)地源熱泵系統(tǒng)土壤熱失衡現(xiàn)象尤為嚴(yán)重,不能保證機(jī)組的正常運(yùn)行條件,說明在冷負(fù)荷占主導(dǎo)的建筑需要輔助冷卻散熱裝置.
2)長期運(yùn)行后,傳統(tǒng)地源熱泵系統(tǒng)COP由2.78降至2.56,聯(lián)合冷卻塔運(yùn)行的復(fù)合式地源熱泵系統(tǒng)維持在2.94左右,聯(lián)合冷水機(jī)組運(yùn)行的復(fù)合式地源熱泵系統(tǒng)維持在3.34左右;可以看出聯(lián)合冷水機(jī)組運(yùn)行的復(fù)合式地源熱泵系統(tǒng)具有較高的COP以及穩(wěn)定性.
3)3種系統(tǒng)平均全年累計(jì)運(yùn)行能耗分別為5.5×106kW·h,4.6×106kW·h和2.8×106kW·h.由于聯(lián)合冷水機(jī)組運(yùn)行的復(fù)合式地源熱泵系統(tǒng)增加了冷水機(jī)組,按照50%的負(fù)荷率分配給熱泵機(jī)組和冷水機(jī)組,大大降低了能耗,可見該系統(tǒng)運(yùn)行能耗最低且比較穩(wěn)定.
4)從初投資和年運(yùn)行費(fèi)用上來看,輔助冷水機(jī)組的復(fù)合式地源熱泵系統(tǒng)的動態(tài)費(fèi)用年值最少,有更好的經(jīng)濟(jì)性,該系統(tǒng)有更好的推廣價(jià)值和應(yīng)用前景.
References)
[1] 劉洪,戴文浩.地源熱泵技術(shù)與建筑節(jié)能[J].節(jié)能技術(shù),2011,25(2):176.Liu Hong,Dai Wenhao.Ground source heat pump technology and building energy conservation[J].Energy Conservation Technology,2011,25(2):176.(in Chinese)
[2] 閡科峰.夏熱冬冷地區(qū)復(fù)合式地源熱泵適宜性研究[J].建筑熱能通風(fēng)空調(diào),2014,33(5):57-59.Min Kefeng.The suitability study of hybrid ground source heat pump in hot summer and cold winter region.[J].Building Energy&Environment,2014,33(5):57-59.(in Chinese)[3] 趙忠超,豐威仙.太陽能輔助地源熱泵供暖實(shí)驗(yàn)研究[J].江蘇科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,29 (3):261-265.Zhao Zhongchao,F(xiàn)eng Weixian.Experimental study of a solar assisted ground source heat pump for heating[J].Journal of Jiangsu University of Science and Technology(Natural Science Edition),2014,29(3): 261-265.(in Chinese)
[4]Frederic K,Joseph V,Kevyn J.Development and validation of a new TRNSYS type for the simulation of external building walls containing PCM[J].Energy and Buildings,2010,42(7):1004-1009.
[5] 李曉,王偉紅.冷卻塔輔助排熱的復(fù)合地源熱泵系統(tǒng)及其運(yùn)行控制[J].制冷與空調(diào),2013,27(5): 466-468.Li Xiao,Wang Weihong.The composite GSHP system equipped with cooling tower for removing heat and its operation control[J].Refrigeration and Air Conditioning,2013,27(5):466-468.(in Chinese)
[6] 羅敏.輔助冷卻塔復(fù)合式地源熱泵系統(tǒng)關(guān)鍵技術(shù)問題研究[D].重慶:重慶大學(xué),2013:24-26.
[7] 李玉街,王琪瑋.中央空調(diào)冷水機(jī)組能源效率與負(fù)荷率動態(tài)匹配節(jié)能技術(shù)[J].智能建筑電氣技術(shù),2012(4):1729-1275.Li Yujie,Wang Qiwei.Energy-saving technologies of dynamic adaptation between energy efficiency and load rate for the chillers of central air-conditioning[J].Electrical Technology of Intelligent Buildings,2012 (4):1729-1275.
[8]Scott Hackel,Amanda Pertzborn.Effective design and operation of hybrid ground-source heat pumps:three case studies[J].Energy and Buildings,2013,43 (12):3497-3504.
[9]中華人民共和國建設(shè)部,中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 50366-2009.地源熱泵系統(tǒng)工程技術(shù)規(guī)范[S].北京:中國建筑工業(yè)出版社,2009.
[10] 徐永軍.重慶地區(qū)復(fù)合式地源熱泵系統(tǒng)的應(yīng)用及控制策略研究[D].重慶:重慶大學(xué),2011:8-9.
(責(zé)任編輯:曹 莉)
Performance analysis on hybrid ground source heat pump systems with auxiliary cooling
Yang Xinglin1,Yang Huanhuan1,Zhang Lisong2,Liu Mingyuan3
(1.School of Energy and Power Engineering,Jiangsu University of Science and Technology,Zhenjiang Jiangsu 212003,China)
(2.Jiangxi CSSC Nautical Instrument Ltd.Co.,Jiujiang Jiangxi 332000,China)
(3.No.92337 Unit of the Chinese People's Liberation Army,Dalian Liaoning 116023,China)
Two kinds of hybrid ground source heat pump systems with auxiliary cooling are put forward.A typical office building in Zhenjiang is taken as an example to build the model for hybrid ground source heat pump system using the TRNSYS software.The performance and energy consumption of the new hybrid ground source heat pump system are compared with those of the traditional ground source heat pump system.The results show that hybrid ground source heat pump systems with auxiliary cooling equipped with water chillers have lower soil temperature rises,more stable and efficient performance,lower system energy consumption and better economy.
HGSHP;coefficient of performance; energy consumption
U664.86
A
1673-4807(2015)06-0540-06
10.3969/j.issn.1673-4807.2015.06.007
2015-06-29
楊興林(1964—),男,教授,博士,研究方向?yàn)樾滦涂照{(diào)系統(tǒng)設(shè)計(jì)及性能.E-maill:hcyangxl2010@163.com
楊興林,楊歡歡,張禮松,等.輔助冷卻的復(fù)合式地源熱泵系統(tǒng)性能分析[J].江蘇科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,29(16):540-545.
江蘇科技大學(xué)學(xué)報(bào)(自然科學(xué)版)2015年6期