王石,吳峰(河海大學(xué)能源與電氣學(xué)院,江蘇南京 211100)
?
考慮風(fēng)電置信水平的機(jī)組組合優(yōu)化方法
王石,吳峰
(河海大學(xué)能源與電氣學(xué)院,江蘇南京211100)
摘要:置信水平是將風(fēng)電功率波動(dòng)由概率形式體現(xiàn)的。隨著節(jié)能減排的需求越來越大,考慮置信水平能夠使系統(tǒng)機(jī)組不過多的安排機(jī)組備用,以免增加系統(tǒng)的運(yùn)行成本?;跈C(jī)組組合的方法,提出了一種考慮置信水平下的備用決策方法。通過考慮風(fēng)電功率波動(dòng)的概率分布和負(fù)荷波動(dòng)的概率分布,建立聯(lián)合概率密度函數(shù),然后加入置信水平,得到不同置信水平下的系統(tǒng)基本發(fā)電成本和校正調(diào)度成本。采用IEEERTS算例,得出系統(tǒng)總成本。
關(guān)鍵詞:置信水平;機(jī)組組合;概率分布;基本運(yùn)行成本;校正調(diào)度成本
Project Supported by the National Natural Science Foundation of China(51422701).
由于風(fēng)電的間歇性和隨機(jī)性,其輸出功率隨機(jī)波動(dòng)[1],電網(wǎng)必須維持足夠的旋轉(zhuǎn)備用保證電網(wǎng)的安全穩(wěn)定運(yùn)行。但由于正負(fù)旋轉(zhuǎn)備用都具有一定的成本,過量的旋轉(zhuǎn)備用總量會(huì)降低電網(wǎng)的經(jīng)濟(jì)性。因此,對(duì)備用總量的決策與優(yōu)化是維持電網(wǎng)安全穩(wěn)定運(yùn)行、提高電網(wǎng)運(yùn)行經(jīng)濟(jì)性的關(guān)鍵。
目前,針對(duì)旋轉(zhuǎn)備用的決策與優(yōu)化,相比于日內(nèi)多級(jí)的有功調(diào)度[2]而言,日前機(jī)組組合方法的研究更為成熟[3-4]。為了應(yīng)對(duì)風(fēng)電的隨機(jī)性和間歇性,常常過于頻繁地調(diào)整機(jī)組出力,這對(duì)電網(wǎng)經(jīng)濟(jì)性帶來巨大沖擊。
為了防止不同時(shí)段機(jī)組調(diào)整出力過多,減少系統(tǒng)調(diào)度成本,本文通過考慮風(fēng)力發(fā)電輸出功率的置信水平[5-6],將風(fēng)電功率出力的不確定性以概率形式體現(xiàn),研究在不同置信水平下,系統(tǒng)備用應(yīng)對(duì)風(fēng)電功率波動(dòng)的分配與決策,以提高電網(wǎng)運(yùn)行的經(jīng)濟(jì)性[7-8]。
文獻(xiàn)[9]考慮風(fēng)速波動(dòng)與負(fù)荷波動(dòng)的聯(lián)合概率分布,確定備用決策與分配,但并沒有考慮風(fēng)速預(yù)測(cè)隨時(shí)間越長(zhǎng)越難以預(yù)測(cè)這一特點(diǎn);文獻(xiàn)[10]以發(fā)電成本與停電損失最小來決策出備用決策與安排,但沒有考慮機(jī)組自身的爬坡約束;文獻(xiàn)[11]通過對(duì)負(fù)荷波動(dòng)、風(fēng)電功率波動(dòng)進(jìn)行抽樣統(tǒng)計(jì),目標(biāo)是為了使發(fā)電成本與停電損失最小,同樣沒有考慮機(jī)組爬坡約束和網(wǎng)絡(luò)約束的影響。
本文考慮風(fēng)電預(yù)測(cè)的置信度,采用風(fēng)電功率區(qū)間預(yù)測(cè),風(fēng)電可能的波動(dòng)區(qū)間隨時(shí)間變長(zhǎng)越來越大,同時(shí),在考慮機(jī)組爬坡約束和網(wǎng)絡(luò)安全約束的基礎(chǔ)上,使得電網(wǎng)運(yùn)行經(jīng)濟(jì)性達(dá)到最優(yōu)。
電網(wǎng)運(yùn)行時(shí)留有備用是為了應(yīng)對(duì)系統(tǒng)中的不確定性,進(jìn)而保持系統(tǒng)安全穩(wěn)定運(yùn)行?,F(xiàn)在電網(wǎng)中的不確定性主要包括系統(tǒng)故障的不確定性,風(fēng)電功率預(yù)測(cè)的不確定性和負(fù)荷預(yù)測(cè)的不確定性。本文著重研究后兩者不確定性條件下,電力系統(tǒng)備用的安排和決策。
1.1風(fēng)電功率概率模型
準(zhǔn)確地預(yù)測(cè)風(fēng)電功率,是保證含大規(guī)模風(fēng)力發(fā)電電力系統(tǒng)安全、可靠和經(jīng)濟(jì)運(yùn)行的關(guān)鍵。由經(jīng)驗(yàn)可知,風(fēng)速服從韋伯爾分布,但目前對(duì)風(fēng)電功率并沒有一個(gè)準(zhǔn)確的模型。文獻(xiàn)[12]得出風(fēng)電功率服從正態(tài)分布,PW
為風(fēng)電功率預(yù)測(cè)值,ΔPW為風(fēng)電功率預(yù)
測(cè)誤差,則實(shí)際風(fēng)電功率PW=PW+ΔPW,實(shí)際風(fēng)電功率概率分布函數(shù)為:
式中:PW為實(shí)際風(fēng)電功率;σPW為風(fēng)電功率的標(biāo)準(zhǔn)差。
1.2負(fù)荷預(yù)測(cè)誤差概率模型
目前,負(fù)荷的預(yù)測(cè)已經(jīng)比較準(zhǔn)確,但還是存在一定的偏差,常采用正態(tài)分布來表示負(fù)合的預(yù)測(cè)誤差[11],ΔPL為負(fù)荷預(yù)測(cè)誤差,其概率密度函數(shù)為:
式中:ΔPL為負(fù)荷預(yù)測(cè)誤差;σL為負(fù)荷預(yù)測(cè)誤差的標(biāo)準(zhǔn)差。
1.3風(fēng)電功率和負(fù)荷的聯(lián)合概率密度函數(shù)
本文綜合考慮2種不確定性因素對(duì)電網(wǎng)運(yùn)行的影響,設(shè)Z為系統(tǒng)功率偏差,即:
式中:PGΣ為所有常規(guī)機(jī)組出力之和。當(dāng)Z>0時(shí),表示發(fā)電量過多,已經(jīng)超過了系統(tǒng)需求,需要采取棄風(fēng)的方式;當(dāng)Z<0時(shí),表示發(fā)電量不足以滿足系統(tǒng)需求,需要采用切負(fù)荷的方式。
由于風(fēng)電預(yù)測(cè)功率和負(fù)荷預(yù)測(cè)誤差之間相互獨(dú)立,可對(duì)2部分的概率密度函數(shù)進(jìn)行卷積,可得隨機(jī)變量Z的概率密度函數(shù)為:
式中:p0為風(fēng)電場(chǎng)輸出功率為0時(shí)的概率;pN為風(fēng)電場(chǎng)輸出功率為額定功率PN時(shí)的概率;fPW(pW)為風(fēng)電場(chǎng)輸出功率為(0,PN)的概率密度函數(shù)。風(fēng)電功率在(0,PN)之間波動(dòng),目前廣泛采用變速恒頻風(fēng)機(jī)模型[12],即:
式中:vin、vout和vN為切入、切出風(fēng)速和額定風(fēng)速。單純風(fēng)電功率點(diǎn)預(yù)測(cè)不能體現(xiàn)隨時(shí)間變長(zhǎng)越難以預(yù)測(cè)這一特點(diǎn),通過采用區(qū)間預(yù)測(cè)的方式,不僅能直觀地體現(xiàn)出時(shí)間特性,同時(shí)也能更好地表現(xiàn)出風(fēng)電出力的概率形式。
2.1風(fēng)電功率預(yù)測(cè)區(qū)間下的數(shù)學(xué)模型
伴隨時(shí)間尺度的變大,風(fēng)電功率預(yù)測(cè)區(qū)間也會(huì)變大,功率率波動(dòng)PW在[PW分別為風(fēng)電功率預(yù)測(cè)區(qū)間上、下限。由此,原先由點(diǎn)預(yù)測(cè)產(chǎn)生的隨機(jī)變量Z的概率密度函數(shù)變?yōu)椋海琍W]之間,PW、PW
式中:pW為風(fēng)電功率取PW為風(fēng)電功率取時(shí)的概率;pWPW時(shí)的概率。由概率論[13-15]可知,在連續(xù)取值范圍內(nèi),某一點(diǎn)取值的概率為0。則式(6)可變換成:
2.2考慮風(fēng)電功率置信水平的數(shù)學(xué)模型變換
風(fēng)電場(chǎng)不光提供風(fēng)電功率基準(zhǔn)值和風(fēng)電功率預(yù)測(cè)區(qū)間,一般還會(huì)提供風(fēng)電功率的概率密度函數(shù),本文采用正態(tài)分布函數(shù)進(jìn)行模擬。根據(jù)概率密度函數(shù),取不同的風(fēng)電置信水平要求,可以得到不同的風(fēng)電功率允許區(qū)間。
假設(shè)已知風(fēng)電功率的概率密度函數(shù),根據(jù)文獻(xiàn)[16]可知,通過基于Nataf變換的蒙特卡洛法可得到風(fēng)電功率的累計(jì)概率函數(shù)。由此可得到風(fēng)電功率的允許區(qū)間上下限:
式中:β為置信水平,ξt為風(fēng)電功率在置信水平β下允許區(qū)間的邊界值。由此,可得出24時(shí)段的風(fēng)電功率預(yù)測(cè)區(qū)間,根據(jù)式(7),可將原式改成:
、ξt
當(dāng)正旋轉(zhuǎn)備用全部投入時(shí),功率仍有缺額,即Z+RU<0時(shí),需要切負(fù)荷,可得出失負(fù)荷期望為:
式中:ΔPloss=-(Z+RUt);
當(dāng)負(fù)旋轉(zhuǎn)備用全部投入時(shí),功率仍有富余,即Z-RD>0時(shí),需要棄風(fēng),可得出棄風(fēng)期望為:
式中:ΔWloss=Z-RDt。
3.1備用決策
由于正旋轉(zhuǎn)備用與失負(fù)荷相關(guān),負(fù)旋轉(zhuǎn)備用與棄風(fēng)相關(guān),而正負(fù)旋轉(zhuǎn)備用、失負(fù)荷與棄風(fēng)都存在一定的經(jīng)濟(jì)成本,可建立2個(gè)以正、負(fù)旋轉(zhuǎn)備用為變量的函數(shù):
式中:cU正負(fù)旋轉(zhuǎn)備用成本;cL為失負(fù)荷成本;cW為棄風(fēng)成本。將式(11)和式(12)代入式(13),可得:R、cDR為
基于式(11)進(jìn)行優(yōu)化,即可得到經(jīng)濟(jì)最優(yōu)時(shí)對(duì)應(yīng)的總正負(fù)旋轉(zhuǎn)備用RUt、RDt。
3.2備用分配
采用優(yōu)先順序法[17-18]確定機(jī)組啟停狀態(tài),優(yōu)先順序法的指標(biāo)為機(jī)組的滿負(fù)荷平均煤耗率:
式中:ai、bi和ci為機(jī)組i的發(fā)電成本系數(shù)。根據(jù)滿負(fù)荷平均煤耗率對(duì)機(jī)組進(jìn)行排序,依次從大到小加入機(jī)組,且在各時(shí)段滿足機(jī)組出力約束、機(jī)組爬坡約束、機(jī)組最小啟停約束和網(wǎng)絡(luò)安全約束,通過采用二分法的經(jīng)濟(jì)調(diào)度方式分配機(jī)組出力,流程圖如圖1所示。
在電力系統(tǒng)運(yùn)行中,成本分為基本運(yùn)行成本和調(diào)度成本?;具\(yùn)行成本S0指的就是發(fā)電機(jī)組出力成本之和,其中:
而校正調(diào)度成本SJ為重新調(diào)整機(jī)組出力成本、釋放備用容量成本、再棄風(fēng)和再切負(fù)荷成本成本,其中:
式中:c0為重新調(diào)整單位機(jī)組出力成本;~Pi,t為 第i機(jī)組第t時(shí)刻的出力變化量;cR為釋放單位備用容量成本;Rt為第t時(shí)刻釋放備用容量總和;L′los(st)和W′los(st)為第t時(shí)刻再切負(fù)荷與再棄風(fēng)期望。
由本文目的可知,目標(biāo)函數(shù)為:
式中:S為系統(tǒng)總成本,根據(jù)改變?chǔ)碌娜≈?,找到系統(tǒng)最小總成本與此時(shí)對(duì)應(yīng)β的取值。
本文中采用IEEE24節(jié)點(diǎn)系統(tǒng)負(fù)荷與風(fēng)電預(yù)測(cè)區(qū)間數(shù)據(jù),如圖2所示,該風(fēng)電場(chǎng)接在9號(hào)節(jié)點(diǎn)。由于各時(shí)段預(yù)測(cè)區(qū)間占u±3δ區(qū)間的99.7%[16],由此可得到風(fēng)電風(fēng)電功率的累計(jì)函數(shù)。取置信水平β=0.5,
圖2風(fēng)電功率預(yù)測(cè)區(qū)間Fig. 2 Prediction interval of wind power
再根據(jù)式(25)和式(26),可得到此時(shí)風(fēng)電功率的允許區(qū)間,如圖3所示。由于當(dāng)前對(duì)風(fēng)電功率的概率密度函數(shù)沒有一個(gè)確定的模型,本文主要考慮備用的決策與分配,取正旋轉(zhuǎn)備用成本為20 $/(MW·h),失負(fù)荷成本為1 000 $/(MW·h),負(fù)旋轉(zhuǎn)備用成本為11 $/(MW·h),棄風(fēng)成本為50 $/(MW·h),再根據(jù)式(11)可得出風(fēng)電預(yù)測(cè)區(qū)間下的正旋轉(zhuǎn)備用與切負(fù)荷成本、負(fù)旋轉(zhuǎn)備用與棄風(fēng)成本。根據(jù)計(jì)算可得,隨著正旋轉(zhuǎn)備用的逐漸增加,失負(fù)荷期望逐漸減小,失負(fù)荷和正旋轉(zhuǎn)備用成本先減小后增大,且隨著負(fù)旋轉(zhuǎn)備用的逐漸增加,棄風(fēng)期望逐漸減小,棄風(fēng)和負(fù)旋轉(zhuǎn)備用成本先減小后增大。由此可得到各個(gè)時(shí)段的正負(fù)旋轉(zhuǎn)備用總量,根據(jù)二分調(diào)度發(fā)對(duì)備用進(jìn)行分配,且滿足各約束條件。同理,求出允許區(qū)間下的備用決策與分配。由圖4、圖5可知,不同置信水平?jīng)Q策出的備用總量不一樣,隨著置信水平的增加,所需備用容量在減少,同時(shí)意味著調(diào)度中所需釋放的備用容量在增加。
通過采用蒙特卡羅法,抽取1 000組在預(yù)測(cè)區(qū)間內(nèi)、服從正態(tài)分布的風(fēng)電功率出力曲線,算出各組數(shù)據(jù)下的調(diào)度成本,最后取均值,可得到校正調(diào)度成本。設(shè)重新調(diào)整機(jī)組出力成本為80 $/(MW·h),備用釋放成本為50 $/(MW·h),切負(fù)荷成本為1 000 $/(MW·h)。對(duì)于不同的置信水平,可得到不同的調(diào)度成本,如圖6所示。隨著置信水平的增加,系統(tǒng)校正調(diào)度成本在逐漸增加。
圖3風(fēng)電功率預(yù)測(cè)區(qū)間與允許區(qū)間Fig. 3 Prediction interval and allowable interval of wind power
圖4不同置信水平下的正旋轉(zhuǎn)備用總量Fig. 4 Up spinning reserve amount in different confidence levels
圖5不同置信水平下的負(fù)旋轉(zhuǎn)備用總量Fig. 5 Down spinning reserve amount in different confidence levels
圖6校正調(diào)度成本Fig. 6 Calibration schedule cost
而系統(tǒng)基本運(yùn)行成本在不斷減小。由此可知,系統(tǒng)總成本如圖8所示,在置信水平為0.6時(shí)達(dá)到經(jīng)濟(jì)最優(yōu)。
圖7基本運(yùn)行成本Fig. 7 Basic operating cost
圖8系統(tǒng)總成本Fig. 8 The system’s total cost
本文針對(duì)風(fēng)電功率的預(yù)測(cè)方式不同,由原先的點(diǎn)預(yù)測(cè)變成區(qū)間預(yù)測(cè),可展現(xiàn)出風(fēng)電功率隨時(shí)間的增加,越來越難以預(yù)測(cè)的特點(diǎn)。加入置信水平之后,可通過改變風(fēng)電功率區(qū)間,使波動(dòng)變小。通過比較預(yù)測(cè)區(qū)間下和加入置信水平允許區(qū)間下的機(jī)組出力與備用決策和安排,可以看出,隨著置信水平的增加,風(fēng)電功率波動(dòng)越來越小,但所需要付出的校正調(diào)度成本也在逐漸增加。由此可知,加入置信水平后,可以權(quán)衡系統(tǒng)基本運(yùn)行成本和校正調(diào)度成本,使系統(tǒng)達(dá)到經(jīng)濟(jì)最優(yōu)。
參考文獻(xiàn)
[1] GU Yingzhong,XIE Le.Look-ahead dispatch with forecast uncertainty and infeasibility managemen[C]// IEEE Power Engineering Society General Meeting,July 22-26,2012,San Diego,USA:7p.
[2]張國(guó)強(qiáng),吳文傳,張伯明.考慮風(fēng)電接入的有功備用協(xié)調(diào)優(yōu)化[J].電力系統(tǒng)自動(dòng)化,2012,40(13): 110-114. ZHANG Guoqiang,WU Wenchuan,ZHANG Boming.Optimization of operation reserve coordination considering wind power integration[J].Automation of Electric Power Systems,2012,40(13):110-114(in Chinese).
[3]陳皓勇,王錫凡.機(jī)組組合問題的優(yōu)化方法綜述[J].電力系統(tǒng)自動(dòng)化,1991,23(5):51-56. CHEN Haoyong,WANG Xifan.A survey of optimizationbased methods for unit commitment[J].Automation of Electric Power Systems,1991,23(5):51-56(in Chinese).
[4]楊朋朋.機(jī)組組合理論與算法研究[D].濟(jì)南:山東大學(xué),2008.
[5]崔曉丹,李威,任先成,等.大規(guī)模風(fēng)電接入的輸電網(wǎng)規(guī)劃研究評(píng)述[J].江蘇電機(jī)工程,2012,31(6):1-6.CUI Xiaodan,LI Wei,REN Xiancheng,et al.Review of the research on transmission network planning with large scale wind power connected[J].Jiangsu Electrical Engineering,2012,31(6):1-6.
[6]朱雨彤,靳晶.多時(shí)間尺度相協(xié)調(diào)的電力市場(chǎng)價(jià)值評(píng)估[J].江蘇電機(jī)工程,2015,34(5):47-50.ZHU Yutong,JIN Jing.Evaluation on the value of multiple time-scale electricity market[J].Jiangsu Electrical Engineering,2015,34(5):47-50(in Chinese).
[7]李智,張新松,郭曉麗.大規(guī)模風(fēng)電接入火電系統(tǒng)的最優(yōu)旋轉(zhuǎn)備用容量研究[J].電力系統(tǒng)保護(hù)與控制,2012 (13):110-114.LI Zhi,ZHANG Xinsong,GUO Xiaoli.Study on the optimal spinning reserve in thermal power system with significant wind power penetration[J]. Power System Protection and Control,2012,40(13):110-114(in Chinese).
[8]張新松,袁越,傅質(zhì)馨.基于隱性備用約束的機(jī)組組合模型[J].電力系統(tǒng)保護(hù)與控制,2013(1):136-142.ZHANG Xinsong,YUAN Yue,F(xiàn)U Zhixin.The unit commitment formulation with implicit reserve constraint[J]. Power System Protection and Control,2013,41(1):136-142(in Chinese).
[9]趙晉泉,唐潔,羅衛(wèi)華,等.一種含風(fēng)電電力系統(tǒng)的日前發(fā)電計(jì)劃和旋轉(zhuǎn)備用決策模型[J].電力自動(dòng)化設(shè)備,2014(5):21-27.ZHAO Jinquan,TANG Jie,LUO Weihua,et al.Dayahead generation scheduling and spinning reserve decisionmaking model for power grid containing wind power[J].Electric Power Automation Equipment,2014,34(5):21-27(in Chinese).
[10]李志剛,吳文傳,張伯明.消納大規(guī)模風(fēng)電的魯棒區(qū)間經(jīng)濟(jì)調(diào)度(一)調(diào)度模式與數(shù)學(xué)模型[J].電力系統(tǒng)自動(dòng)化,2014(20):33-39.LI Zhigang,WU Wenchuan,ZHANG Boming.A robust interval economic dispatch method accommodating largescale wind power generation:Part one dispatch scheme and mathematical model[J].Automation of Electric Power Systems,2014,38(20):33-39(in Chinese).
[11]李志剛,吳文傳,張伯明.消納大規(guī)模風(fēng)電的魯棒區(qū)間經(jīng)濟(jì)調(diào)度(二)不確定集合構(gòu)建與保守度調(diào)節(jié)[J].電力系統(tǒng)自動(dòng)化,2014(21):32-38.LI Zhigang,WU Wenchuan,ZHANG Boming.A robust interval economic dispatch method accommodating largescale wind power generation:Part two uncertainty set modeling and conservativeness adjustment[J].Automation of Electric Power Systems,2014,38(21):32-38(in Chinese).
[12]范高峰,王偉勝,劉純,等.基于人工神經(jīng)網(wǎng)絡(luò)的風(fēng)電功率預(yù)測(cè)[J].中國(guó)電機(jī)工程學(xué)報(bào),2008,28(34):118-123.FAN Gaofeng,WANG Weisheng,LIU Chun,et al.Wind power perdiction based on artificial neural network[J].Proceedings of the CSEE,2008,28(34):118-123(in Chinese).
[13] JUAN M M,JUAN P R.Point estimate schemes to solve the probabilistic power flow[J].IEEE Transactions on Power Systems,2007,22(4):1594-1601.
[14] ZHAO W,CHEN Z,Blaabjerg F.Probabilistic capacity of a grid connected wind farm based on optimization method[J].Renewable Energy,2006,31(13):2171-2187.
[15]盛驟,謝式千,潘程毅.概率論與數(shù)理統(tǒng)計(jì)[M].北京:高等教育出版社,2004.
[16] CHEN Y,WEN J,CHENG S.Probabilistic load flow method based on Nataf transformation and Latin hypercube sampling[J].IEEE Trans on Sustainable Energy,2013,4 (2):294-301.
[17]黎靜華.適合于機(jī)組組合問題的擴(kuò)展優(yōu)先順序法[J].電力系統(tǒng)保護(hù)與控制,2010,38(2):1-7.LI Jinghua.Extended priority list method for unit commitment problem[J].Power System Protection and Control[J]. 2010,38(2):1-7(in Chinese).
[18] PUKELSHEIM F.The three sigma rule[J].The American Statistician,1994,48(2):88-91.
王石(1991—),男,碩士研究生,研究方向?yàn)樾履茉唇尤雮溆萌萘績(jī)?yōu)化;
吳峰(1977—),男,博士,教授,博士生導(dǎo)師,研究方向?yàn)榭稍偕茉窗l(fā)電系統(tǒng)的建模與控制。
(編輯董小兵)
Optimization Methods of Unit Commitment Using a Confidence Levels’Method
WANG Shi,WU Feng
(College of Energy & Electrical Engineering,Hohai University,Nanjing 211100,Jiangsu,China)
ABSTRACT:Confidence levels reflect wind power volatility in the way of probability. With the growth of demand of energy conservation and emissions reduction,it is useful for avoiding too much unit reserve to consider confidence levels,so that it can avoid increasing the cost of system operation. This paper is based on unit commitment to come up with a way of considering confidence levels’reserve determination. Considering the probability distributions of the wind power volatility and the load volatility,it can establish the function of joint probability density,and add the confidence levels to the function,so as to calculate the system basic operation cost and calibration schedule cost in different confidence levels. IEEE-RTS is taken as the example to demonstrate the efficiency of the method proposed in this paper.
KEY WORDS:confidence levels;unit commitment;probability distribution;basic operating cost;calibration schedule cost
作者簡(jiǎn)介:
收稿日期:2015-10-17。
基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51422701)。
文章編號(hào):1674- 3814(2016)02- 0116- 06
中圖分類號(hào):TM614
文獻(xiàn)標(biāo)志碼:A