梁坤峰,任峴樂,賈雪迎,王林,阮春蕾,張林泉
(河南科技大學(xué)制冷與空調(diào)技術(shù)研究所,河南洛陽 471003)
冰蓄冷輻射空調(diào)系統(tǒng)運(yùn)行優(yōu)化與能耗分析
梁坤峰*,任峴樂,賈雪迎,王林,阮春蕾,張林泉
(河南科技大學(xué)制冷與空調(diào)技術(shù)研究所,河南洛陽 471003)
基于溫濕度獨(dú)立控制思路,提出了一種將冰蓄冷、毛細(xì)管輻射及地源耦合的冰蓄冷輻射空調(diào)系統(tǒng),以實(shí)現(xiàn)電力負(fù)荷“削峰填谷”、降低能耗和更好的熱舒適性。建立了冰蓄冷輻射空調(diào)系統(tǒng)的運(yùn)行優(yōu)化模型,以杭州地區(qū)某辦公樓為實(shí)例,進(jìn)行了負(fù)荷計(jì)算、運(yùn)行方案優(yōu)化以及能耗分析,結(jié)果表明冰蓄冷輻射空調(diào)系統(tǒng)與常規(guī)空調(diào)系統(tǒng)相比,具有顯著的節(jié)電和節(jié)省運(yùn)行費(fèi)用優(yōu)勢(shì)。實(shí)際應(yīng)用中需綜合考慮蓄冰容量的限制以最終確定系統(tǒng)運(yùn)行方案。
溫濕度獨(dú)立控制;冰蓄冷;輻射空調(diào);運(yùn)行優(yōu)化;能耗
冰蓄冷利用水/冰相變潛熱蓄存能量,夜間用低谷廉價(jià)電制冰蓄冷,白天用電高峰融冰釋冷,滿足建筑物空調(diào)或生產(chǎn)工藝的用冷需求。用冷地區(qū)的峰谷電價(jià)差使其運(yùn)行費(fèi)用顯著減少,已成為區(qū)域供冷冷源選擇的有力競(jìng)爭(zhēng)者[1-5]。而最近研究較多的頂板輻射和置換通風(fēng)相結(jié)合的輻射空調(diào)系統(tǒng)[6-8],其基于溫濕度獨(dú)立控制思想與傳統(tǒng)的熱濕耦合處理空氣的空調(diào)相比,室內(nèi)環(huán)境的熱舒適性、能耗、運(yùn)行費(fèi)用、吹風(fēng)感和噪聲等得以顯著改善。目前,從國(guó)內(nèi)外關(guān)于冰蓄冷和輻射空調(diào)的公開文獻(xiàn)看[9-12],兩個(gè)系統(tǒng)基本都是互相獨(dú)立,很少有將兩種空調(diào)技術(shù)結(jié)合的研究。因此,開發(fā)一種冰蓄冷輻射空調(diào)系統(tǒng),既利用冰蓄冷空調(diào)的電力調(diào)峰潛力,又使輻射空調(diào)發(fā)揮其高舒適性、低能耗的特性,具有重要的研究?jī)r(jià)值和現(xiàn)實(shí)意義。
日本的中野幸夫等[13]研究了一種蓄冰和輻射供冷結(jié)合的空調(diào)系統(tǒng),該系統(tǒng)供冷時(shí)取冰作為蓄冷介質(zhì),從蓄冰槽中制取3 ℃~4 ℃的冷水供給空氣處理機(jī)組,機(jī)組產(chǎn)生的低溫低濕空氣送入安裝在房間頂板與吊頂上的金屬輻射板之間,并通過對(duì)流方式將輻射板冷卻,降溫后的輻射板則通過輻射方式帶走室內(nèi)熱負(fù)荷;冷空氣經(jīng)對(duì)流換熱后升溫,從窗側(cè)的吊頂末端送入房間,用于消除室內(nèi)全部余濕和部分余熱。但由于該系統(tǒng)送入輻射頂板的空氣溫度較低(10 ℃~11 ℃),易于在金屬輻射板表面結(jié)露,影響其使用壽命;此外,該系統(tǒng)中新風(fēng)實(shí)際上承擔(dān)了室內(nèi)全部負(fù)荷,其節(jié)能的特征并不顯著。為此,本文基于溫濕度獨(dú)立控制方法提出了一種新型冰蓄冷輻射空調(diào)系統(tǒng)。
冰蓄冷輻射空調(diào)系統(tǒng)是將冰蓄冷系統(tǒng)和毛細(xì)管輻射空調(diào)系統(tǒng)相結(jié)合。首先,冰蓄冷系統(tǒng)可以在晚間電力低谷時(shí)間段利用雙工況機(jī)組進(jìn)行蓄冰;而在白天電力高峰時(shí)段,再將冷量釋放出來制取冷凍水,來承擔(dān)新風(fēng)除濕所需冷量和一部分毛細(xì)管輻射供冷所需的冷量;獨(dú)立新風(fēng)系統(tǒng)采用置換通風(fēng),使室內(nèi)能保持較高的空氣質(zhì)量;此外,系統(tǒng)還設(shè)置有地埋管構(gòu)成的土壤換熱器,利用地下土壤這一免費(fèi)冷源來為毛細(xì)管輻射供冷提供主要的冷量。圖1為冰蓄冷輻射空調(diào)系統(tǒng)原理圖,主要由三個(gè)子系統(tǒng)構(gòu)成:1)冰蓄冷系統(tǒng);2)獨(dú)立全新風(fēng)除濕系統(tǒng);3)毛細(xì)管輻射供冷系統(tǒng)。
冰蓄冷系統(tǒng)主要由雙工況機(jī)組、蓄冰槽和乙二醇泵,以及末端和下一級(jí)系統(tǒng)相連的板式換熱器組成,其中雙工況機(jī)組的冷凝側(cè)設(shè)有冷卻塔和熱回收裝置,熱回收裝置可以回收一部分冷凝熱來再熱冷凍除濕后的低溫新風(fēng)。系統(tǒng)采用主機(jī)上游串聯(lián)、主機(jī)優(yōu)先工作模式。這樣可使溫度較高的回水首先通過制冷機(jī)組,使機(jī)組保持較高的運(yùn)行效率,節(jié)能性好。同時(shí),可以降低蓄冰槽容量、減少初投資。工作時(shí)可以有四種運(yùn)行模式,分別為:蓄冰模式、融冰供冷模式、制冷機(jī)單獨(dú)供冷和制冷機(jī)融冰聯(lián)合供冷模式。各模式下閥門及泵的開啟狀態(tài),如表1所示。在非標(biāo)準(zhǔn)設(shè)計(jì)日內(nèi),空調(diào)冷負(fù)荷的變化較大,需通過優(yōu)化控制蓄冰槽的有效融冰量以滿足冷負(fù)荷需求。因此,采用制冷機(jī)組優(yōu)先的策略,制冷機(jī)組滿負(fù)荷運(yùn)行,投入的臺(tái)數(shù)及開啟時(shí)間則要由系統(tǒng)整體優(yōu)化控制方案來確定。融冰供冷量為總冷負(fù)荷與制冷機(jī)供冷量的差值,其供冷量隨逐時(shí)負(fù)荷和制冷機(jī)投入臺(tái)數(shù)的變化而變化。
圖1 冰蓄冷輻射空調(diào)系統(tǒng)原理圖
表1 各工作模式下設(shè)備開啟表
室內(nèi)送風(fēng)方式采用置換通風(fēng),即下送風(fēng)上排風(fēng)的方式。置換送風(fēng)有較高的舒適性,由于送風(fēng)速度較小,避免了傳統(tǒng)空調(diào)的風(fēng)感問題。雖然室內(nèi)經(jīng)排風(fēng)口排出的污濁氣體的溫度濕度有所升高,但對(duì)于室外空氣來說,其溫度還是比較低的,所以為避免直接排到室外導(dǎo)致冷量的浪費(fèi),可設(shè)置全熱回收裝置來回收這部分冷量,對(duì)新風(fēng)進(jìn)行預(yù)處理。但由于新風(fēng)承擔(dān)室內(nèi)全部潛熱負(fù)荷和小部分顯熱負(fù)荷,房間的送風(fēng)量很小,為此新風(fēng)的含濕量應(yīng)處理到較小值,才能滿足室內(nèi)除濕要求。一般處理至低于室內(nèi)含濕量4 g/kg以上,由此帶來了新風(fēng)再熱所造成的能源浪費(fèi)問題,可采取以下措施解決:一方面系統(tǒng)利用機(jī)組冷凝熱產(chǎn)生的熱水對(duì)新風(fēng)進(jìn)行再熱升溫處理,機(jī)組停機(jī)時(shí)可以用電加熱來輔助加熱;另一方面在送入房間前的主風(fēng)道上設(shè)置一段室外空氣再熱風(fēng)道,通過調(diào)節(jié)風(fēng)門開度控制進(jìn)入再熱風(fēng)道的風(fēng)量從而調(diào)節(jié)送風(fēng)溫度。圖2為新風(fēng)處理過程示意圖。由系統(tǒng)原理圖可知,在新風(fēng)處理過程中引入雙級(jí)全熱回收空氣處理機(jī)用于新、排風(fēng)的處理。新風(fēng)首先在一級(jí)熱回收裝置中回收室內(nèi)排風(fēng)的冷量,狀態(tài)由W變化到W',由W'變化到W''的冷量則來自于混合空氣由C變化到O'所失去的冷量,新風(fēng)狀態(tài)由 W''經(jīng)過冷凍除濕變化到機(jī)器露點(diǎn) L。隨后通過新、排風(fēng)的混合,新風(fēng)中引入室內(nèi)排風(fēng),達(dá)到送風(fēng)狀態(tài)C。此時(shí),為了實(shí)現(xiàn)室內(nèi)送風(fēng)狀態(tài)點(diǎn)O,需從空氣狀態(tài)O'開始加熱,分別通過預(yù)冷新風(fēng)和引入機(jī)組冷凝熱的方法實(shí)現(xiàn)O'到O的狀態(tài)變化。該系統(tǒng)的空氣處理過程充分利用了室內(nèi)排風(fēng)和機(jī)組冷凝熱,有利于解決新風(fēng)預(yù)冷、再熱過程的輔助能耗問題。
輻射供冷系統(tǒng)主要由毛細(xì)管網(wǎng)輻射末端、分水器、集水器、高溫冷水泵、地埋管環(huán)路以及與上一級(jí)系統(tǒng)相連的板式換熱器組成,輔助設(shè)備主要有膨脹補(bǔ)水箱等。毛細(xì)管網(wǎng)輻射供冷系統(tǒng)是把毛細(xì)管網(wǎng)安裝在室內(nèi)頂棚上,以17 ℃~19 ℃的高溫冷水作為冷源,通過輻射方式進(jìn)行換熱,輻射供冷僅影響室內(nèi)顯熱冷負(fù)荷。
圖2 空氣處理過程示意圖
冰蓄冷系統(tǒng)的運(yùn)行方案需根據(jù)建筑逐時(shí)負(fù)荷特征和設(shè)備配置容量、機(jī)組臺(tái)數(shù)以及當(dāng)?shù)胤骞入妰r(jià)情況等綜合設(shè)計(jì)分析,選取最優(yōu)的運(yùn)行方案。本文分別以100%負(fù)荷、75%負(fù)荷、50%負(fù)荷和25%負(fù)荷為典型負(fù)荷,建立基于全天運(yùn)行費(fèi)用最低的優(yōu)化模型,以獲取該系統(tǒng)在不同負(fù)荷條件下的運(yùn)行方案。計(jì)算模型如下:
設(shè)用戶k時(shí)刻的負(fù)荷為qk,其中制冷機(jī)負(fù)擔(dān)qrk,蓄冰罐負(fù)擔(dān)qik,制冷機(jī)的費(fèi)用為R(qrk),蓄冷罐的費(fèi)用為I(qik)。全天的運(yùn)行費(fèi)M為:
優(yōu)化的約束是:
qrkmax為制冷機(jī)k時(shí)刻的最大制冷能力;qikmax為蓄冰罐k時(shí)刻的最大融冰供冷能力,由蓄冰罐Calmac 1190 A在回水溫度10 ℃、供水溫度6.7 ℃下的融冰曲線用最小二乘法擬合確定:
融冰供冷量:
其中:
x——已融冰供冷量,kW;
t——時(shí)刻,h;
y——各時(shí)刻的最大融冰供冷量,kW·h。
最終得:
由此,優(yōu)化問題轉(zhuǎn)化為取M的最小值:
約束條件:
ak為制冷機(jī)單位供冷負(fù)荷費(fèi)用,等于(94/246)×Ek,Ek為平時(shí)電價(jià)。bk為冰罐負(fù)擔(dān)單位冷負(fù)荷的費(fèi)用,等于(94/162)×E低谷,E低谷為低谷電價(jià),杭州地區(qū)分時(shí)電價(jià)表如表2所示。
表2 杭州地區(qū)分時(shí)電價(jià)
基于杭州地區(qū)某辦公樓進(jìn)行冰蓄冷輻射空調(diào)系統(tǒng)運(yùn)行方案分析。該工程地上6層,高25.1 m。建筑內(nèi)主要有辦公室、會(huì)所、值班室、接待大廳、開敞辦公室、餐廳和休息室等房間。工程位于杭州市,東經(jīng)120°26′、北緯30°10′,建筑面積6,800 m2,總空調(diào)面積為6,150 m2。要求能夠?qū)崿F(xiàn)夏季供冷并能滿足人體的舒適性及衛(wèi)生要求。室內(nèi)外主要設(shè)計(jì)參數(shù)為:夏季空調(diào)室外干球溫度為 35.70 ℃,相對(duì)濕度為80%;夏季室內(nèi)空氣溫度為26 ℃,相對(duì)濕度為 60%。假設(shè)該建筑為辦公室,空調(diào)供冷時(shí)間07:00~18:00,利用清華大學(xué)開發(fā)的建筑環(huán)境設(shè)計(jì)模擬分析軟件DeST (Designer’s SimulationToolkits),計(jì)算得到設(shè)計(jì)日逐時(shí)負(fù)荷圖,如圖3所示。
圖3 設(shè)計(jì)日逐時(shí)負(fù)荷圖
通過上述模型計(jì)算得出了在 100%負(fù)荷、75%負(fù)荷、50%負(fù)荷和25%負(fù)荷情況下不同運(yùn)行方案的日總蓄冰量和日運(yùn)行費(fèi)用,如表3所示。由表3可知在不同負(fù)荷日,當(dāng)總蓄冰量越大時(shí),日運(yùn)行費(fèi)用越小。由于本工程蓄冰量已定,可得出冰蓄冷輻射空調(diào)系統(tǒng)在不同負(fù)荷下的運(yùn)行策略。但由于所建立的模型是一個(gè)線性優(yōu)化模型,在實(shí)際求解冰蓄冷系統(tǒng)運(yùn)行優(yōu)化問題方面,難以僅從模型中所獲得最低運(yùn)行成本角度進(jìn)行選取。
同時(shí),文章還給出了典型負(fù)荷情況下,如100%負(fù)荷、75%負(fù)荷、50%負(fù)荷和 25%負(fù)荷,所優(yōu)化出的運(yùn)行方案,如圖4所示。根據(jù)上述所建立的運(yùn)行優(yōu)化數(shù)學(xué)模型,計(jì)算出冰蓄冷輻射空調(diào)系統(tǒng)在不同負(fù)荷、不同主機(jī)開啟臺(tái)數(shù)情況下的日總蓄冰量及日運(yùn)行費(fèi)用。當(dāng)日運(yùn)行費(fèi)用最少時(shí),總蓄冰量最大,理論上此時(shí)運(yùn)行策略最優(yōu)。但在本工程中蓄冰設(shè)備、蓄冰容量(2,592 kW)已定,使得在日負(fù)荷較大時(shí),不能選取運(yùn)行費(fèi)用最少的運(yùn)行策略。在冷負(fù)荷非常小時(shí),控制策略充分發(fā)揮蓄冰罐的潛力,冷機(jī)完全不運(yùn)行。在冷負(fù)荷較小時(shí),控制策略在滿足用冷負(fù)荷的條件下,盡量利用蓄冰罐融冰來滿足用戶冷負(fù)荷需求。當(dāng)冷負(fù)荷接近典型設(shè)計(jì)日負(fù)荷時(shí),為了滿足用戶冷負(fù)荷的需求,必須控制各時(shí)刻的蓄冰罐的融冰量。
表3 不同負(fù)荷情況下日總蓄冰量和日運(yùn)行費(fèi)用
圖4 典型負(fù)荷運(yùn)行策略圖
冰蓄冷輻射空調(diào)系統(tǒng)與常規(guī)空調(diào)系統(tǒng)配置如表4所示。對(duì)兩個(gè)系統(tǒng)的能耗進(jìn)行對(duì)比,即計(jì)算年耗電量,需明確空調(diào)系統(tǒng)在不同負(fù)荷下的運(yùn)行天數(shù)。故此,本文對(duì)空調(diào)系統(tǒng)的年運(yùn)行情況進(jìn)行必要的假設(shè)。
杭州市地處長(zhǎng)江流域以南,夏季炎熱,供冷時(shí)間約從5月中旬至10月上旬,建筑全年的供冷時(shí)間約為140天,其中100%負(fù)荷日為30天,75%負(fù)荷日為40天,50%負(fù)荷日為45天,25%負(fù)荷日為25天。分別以100%負(fù)荷、75%負(fù)荷、50%負(fù)荷、25%負(fù)荷為基礎(chǔ)進(jìn)行比較,得出夏季本系統(tǒng)與常規(guī)空調(diào)系統(tǒng)日運(yùn)行能耗運(yùn)行耗及整個(gè)供冷季節(jié)運(yùn)行能耗,如表5所示。
表4 冰蓄冷輻射空調(diào)與常規(guī)空調(diào)主要設(shè)備表
表5 冰蓄冷輻射空調(diào)與常規(guī)空調(diào)耗電量及費(fèi)用比較
由表5可知,冰蓄冷輻射空調(diào)系統(tǒng)相對(duì)于常規(guī)空調(diào)系統(tǒng)具有明顯的節(jié)電優(yōu)勢(shì),供冷季節(jié)的能耗僅為常規(guī)空調(diào)系統(tǒng)的47.53%;在經(jīng)濟(jì)上的優(yōu)勢(shì)更加明顯,其在供冷季節(jié)的運(yùn)行費(fèi)僅為常規(guī)空調(diào)系統(tǒng)的39.10%,整個(gè)供冷季節(jié)節(jié)省了17.65萬元的運(yùn)行費(fèi)。但需注意的是,冰蓄冷輻射空調(diào)系統(tǒng)受制于城市峰谷電價(jià)差的影響,與常規(guī)空調(diào)系統(tǒng)的用電時(shí)段及用電量顯著不同。
但是,冰蓄冷輻射空調(diào)的初投資比常規(guī)空調(diào)系統(tǒng)大,主要體現(xiàn)在蓄冷設(shè)備、毛細(xì)管網(wǎng)和毛細(xì)管網(wǎng)的安裝以及地埋管換熱器的制造費(fèi)用上。如能合理設(shè)計(jì)冰蓄冷輻射空調(diào)系統(tǒng),那么影響該空調(diào)系統(tǒng)投資回收期的主要因素就是空調(diào)使用地域的峰谷電價(jià)差;從上面的分析可以看出不同的峰谷電價(jià)差對(duì)投資回收期的影響是非常顯著的,地方政府在其電價(jià)的審批上應(yīng)給予強(qiáng)有力的支持。
1) 基于溫濕度獨(dú)立控制所提出的冰蓄冷輻射空調(diào)系統(tǒng),由冰蓄冷系統(tǒng)、獨(dú)立新風(fēng)除濕系統(tǒng)和毛細(xì)管輻射供冷系統(tǒng)三個(gè)子系統(tǒng)構(gòu)成;該系統(tǒng)將具有電力調(diào)峰能力的冰蓄冷系統(tǒng)和高舒適性、低能耗的輻射空調(diào)系統(tǒng)耦合,并利用地埋管為毛細(xì)管輻射供冷提供主要的高溫冷源。
2) 基于所建立的全天運(yùn)行費(fèi)用最低的優(yōu)化模型,獲取不同負(fù)荷條件下的系統(tǒng)運(yùn)行方案;分別計(jì)算了不同負(fù)荷下的耗電量和運(yùn)行費(fèi)用,與常規(guī)空調(diào)系統(tǒng)相比,文中所提出的冰蓄冷輻射空調(diào)系統(tǒng)具有顯著的節(jié)電和省運(yùn)行費(fèi)用優(yōu)勢(shì)。
3) 在實(shí)際工程中采取的運(yùn)行方案還需考慮蓄冰容量的大小,根據(jù)文中所建立的模型得出的日最少運(yùn)行費(fèi)用的運(yùn)行方案并非實(shí)際運(yùn)行方案,運(yùn)行方案的具體選取還需考慮實(shí)際工程的蓄冰容量限制。
參考文獻(xiàn):
[1]于震,曲凱陽,江億.二次冷媒式過冷水動(dòng)態(tài)蓄冰空調(diào)系統(tǒng)的實(shí)驗(yàn)研究[J].暖通空調(diào),2003,33(4): 44-47.
[2]ZHANG P,MA Z W.An overview of fundamental studies and applications of phase change material slurries to secondary loop refrigeration and air conditioning systems[J].Renewable and Sustainable Energy Reviews,2012,16 (7): 5021-5058.
[3]方貴銀.蓄能空調(diào)技術(shù)[M].北京: 機(jī)械工業(yè)出版社,2006.
[4]WANG H,HE G G,TIAN Q.Experimental study of the supercooling heat exchanger coated with fluorocarbon coating[J].Energy and Buildings,2012,55: 526-532.
[5]張學(xué)軍,田新建,鄭克晴,等.氣體直接接觸式制取冰漿實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào),2010,31(12): 1997-2000.
[6]江億.溫濕度獨(dú)立控制空調(diào)系統(tǒng)[M].北京: 中國(guó)建筑工業(yè)出版社,2006.
[7]張海強(qiáng),劉曉華,江億.溫濕度獨(dú)立控制空調(diào)系統(tǒng)和常規(guī)空調(diào)系統(tǒng)的性能比較[J].暖通空調(diào),2011,41(1): 48-52.
[8]侯東明,王聰.雙冷源溫度濕度獨(dú)立控制空調(diào)系統(tǒng)及其設(shè)計(jì)方法[J].暖通空調(diào),2012,42(11): 100-104.
[9]LI X W,ZHANG X S,SHUO Q.Evaporative supercooling method for ice production[J].Applied Thermal Engineering,2012,37(1): 120-128.
[10]PENG Z B,YUAN Z L,LIANG K F,et al.Ice slurry formation in a cocurrent liquid-liquid flow[J].Chinese Journal of Chemical Engineering,2008,16(4): 552-557.
[11]ZHANG X L,HAN Z,LI Z W.Analysis on IPF influencing factors for vacuum binary ice making method[J].International Journal of Thermal Sciences,2013,67: 210-216.
[12]徐征,劉筱萍,何海亮.溫濕度獨(dú)立控制空調(diào)系統(tǒng)節(jié)能性實(shí)例分析[J].暖通空調(diào),2007,37(6): 129-132.
[13]中野幸夫,宮永俊之,岡建雄,等.冰蓄冷輻射供冷系統(tǒng)的綜合評(píng)價(jià)[J].暖通空調(diào),2002,32(4): 48-52.
Operation Optimization and Energy Consumption Analysis of Radiation Air Conditioning System with Ice Storage
LIANG Kun-feng*,REN Xian-le,JIA Xue-ying,WANG Lin,RUAN Chun-lei,ZHANG Lin-quan
(Institute of Refrigeration and Air Conditioning Technology,Henan University of Science and Technology,Luoyang,Henan 471003,China)
Based on the idea of temperature and humidity independent control,a radiation air conditioning system with ice storage was presented,which combined ice storage system,capillary radiation air conditioning,and ground-source heat exchanger;the new system may realize the power load peak load shifting,the reduction of energy consumption and good thermal comfort.Operation optimization model is established for the ice storage air conditioning system.Taking an office building in Hangzhou area as an example,the load calculation,operation strategy optimization and energy consumption analysis have been conducted.The results show that,the ice storage air conditioning system has the remarkable advantages for saving electricity and operation cost than the conventional air conditioning system.In the practical application,ice storage capacity constraints are considered in order to determine the system operation scheme.
Temperature and humidity independent control;Ice storage;Radiation air conditioning system;Operation optimization;Energy consumption
10.3969/j.issn.2095-4468.2014.06.208
*梁坤峰(1975-),男,博士,副教授。研究方向:蓄能空調(diào)技術(shù)與低品位能源利用。聯(lián)系地址:河南省洛陽市西苑路48號(hào)86#,郵編:471003。聯(lián)系電話:13949230093。E-mail:liangkunf@163.com。
國(guó)家自然科學(xué)基金(U1304521),河南省科技攻關(guān)項(xiàng)目(102102210162)