一種新型風(fēng)力發(fā)電機與飛輪儲能聯(lián)合系統(tǒng)的并網(wǎng)運行控制

劉世林 孫海順 顧明磊 文勁宇

（華中科技大學(xué)強電磁工程與新技術(shù)國家重點實驗室 武漢 430074）

1 引言

近年來，隨著能源和環(huán)境的壓力增大，風(fēng)力發(fā)電作為最成熟的可再生能源發(fā)電技術(shù)得到迅猛發(fā)展。由于風(fēng)能具有隨機性，而且通常對風(fēng)力發(fā)電機進(jìn)行最大功率追蹤控制，導(dǎo)致了風(fēng)電功率的隨機波動很大。隨著風(fēng)電機組單機容量和風(fēng)電場裝機容量的不斷增大，大型風(fēng)電場直接并入輸電系統(tǒng)，風(fēng)力發(fā)電對電力系統(tǒng)的安全穩(wěn)定運行帶來越來越大的影響。為了減小風(fēng)力發(fā)電對電網(wǎng)穩(wěn)定性的影響，電網(wǎng)公司對并網(wǎng)運行風(fēng)電場的最大功率波動進(jìn)行了規(guī)定，且要求風(fēng)電機組采用電壓控制方式。

目前，抑制風(fēng)力發(fā)電有功功率波動通常采用降額發(fā)電的方法，但這種方法直接影響了風(fēng)能的利用效率，且調(diào)節(jié)能力有限[1,2]。而風(fēng)電的無功功率波動通常采用 STATCOM 等 FACTS裝置進(jìn)行調(diào)節(jié)[3]，但無功補償裝置無法抑制有功功率波動。利用電力儲能技術(shù)，在風(fēng)力發(fā)電系統(tǒng)中引入儲能元件，既能在較寬范圍內(nèi)調(diào)節(jié)有功功率，又能調(diào)節(jié)無功功率，是解決風(fēng)力發(fā)電功率波動問題的一種有效途徑[4-8]。常用的儲能方式包括蓄電池、飛輪、超導(dǎo)磁儲能和超級電容器等，各種方式具有各自的特點和適用場合

[9]。在用于改善風(fēng)力發(fā)電機輸出功率特性方面，飛輪儲能在儲能時間、響應(yīng)速度和使用壽命等方面具有突出的優(yōu)勢[8]，而且，在最大儲能時間約為 10min的應(yīng)用場合，飛輪儲能具有最好的經(jīng)濟性[10]，因而飛輪儲能在風(fēng)力發(fā)電中的應(yīng)用成為當(dāng)前研究的熱點。

文獻(xiàn)[7，8]分別研究了兩種基于感應(yīng)電機（IM）的飛輪儲能裝置與風(fēng)力發(fā)電機聯(lián)合運行的方法，但這兩種方法中的飛輪儲能裝置都只能進(jìn)行有功調(diào)節(jié)，不能參與無功功率調(diào)節(jié)。多功能柔性功率調(diào)節(jié)器（Flexible Power Conditioner，F(xiàn)PC）是一種基于雙饋感應(yīng)電機（DFIG）的飛輪儲能裝置，其同時具有有功和無功的獨立調(diào)節(jié)能力。相對基于 IM的飛輪儲能裝置，F(xiàn)PC主要通過定子側(cè)進(jìn)行功率交換，在儲能容量相同的情況下，只需要較小容量的電力電子變換器，從而可以一定程度上減少硬件成本[11-15]。因而，本文研究將 FPC應(yīng)用于風(fēng)力發(fā)電系統(tǒng)中，用于改善風(fēng)力發(fā)電機并網(wǎng)運行特性。

本文根據(jù)變速恒頻風(fēng)力發(fā)電機和FPC的結(jié)構(gòu)特點，提出一種新型風(fēng)力發(fā)電機與 FPC的聯(lián)合系統(tǒng)（WG-FPC），并根據(jù)WG-FPC的結(jié)構(gòu)特點和控制目標(biāo)，給出相應(yīng)的并網(wǎng)運行控制策略。建立了WG-FPC的數(shù)學(xué)建模，在Matlab/Simulink環(huán)境下進(jìn)行了仿真實驗，驗證了所提方法和控制策略的可行性與有效性。

2 數(shù)學(xué)建模

2.1 系統(tǒng)結(jié)構(gòu)

WG-FPC的結(jié)構(gòu)和功率關(guān)系如圖1所示。在該結(jié)構(gòu)中，省掉了FPC的網(wǎng)側(cè)變換器，而將轉(zhuǎn)子側(cè)變換器（RSC）接在風(fēng)力發(fā)電機的全功率變換器的中間直流母線上。FPC用于抑制風(fēng)力發(fā)電機的輸出功率波動和提供無功補償以維持機端電壓穩(wěn)定。

圖1 WG-FPC系統(tǒng)結(jié)構(gòu)與功率關(guān)系Fig.1 Configuration and power relation of WG-FPC

系統(tǒng)中各變換器的功能定義如下：風(fēng)力發(fā)電機側(cè)變換器（WSC），控制發(fā)電機工作在最優(yōu)運行點，實現(xiàn)最大功率追蹤，減小發(fā)電機側(cè)諧波電流；網(wǎng)側(cè)變換器（GSC），采用單位功率因數(shù)控制并維持直流母線的電容電壓（VC）恒定；轉(zhuǎn)子側(cè)變換器（RSC）為FPC提供勵磁，控制FPC的有功功率和無功功率。本文不具體討論各變換器的開關(guān)器件控制，將其視為理想變換器。

2.2 DFIG數(shù)學(xué)模型

由于本文主要關(guān)注接入FPC 對風(fēng)力發(fā)電機并網(wǎng)運行特性的改善，不考慮風(fēng)力發(fā)電機本身的控制過程，認(rèn)為風(fēng)力發(fā)電機按最大功率追蹤控制，僅對其輸出功率進(jìn)行描述，不對其數(shù)學(xué)模型進(jìn)行討論。

本節(jié)討論的DFIG 模型，定子側(cè)參考方向采用發(fā)電機慣例，電流、功率流出定子時為正向；轉(zhuǎn)子側(cè)參考方向采用電動機慣例，電流、功率流入轉(zhuǎn)子時為正向[13,14]。同時假設(shè)：①忽略定子磁鏈的暫態(tài)過程，并將定子磁鏈sψ˙定義在同步坐標(biāo)系的d 軸上；②忽略定子電阻Rs。設(shè)ψs、Vs分別為定子磁鏈和電壓的幅值，則有

基于上述假設(shè)，可推導(dǎo)出此同步旋轉(zhuǎn)坐標(biāo)系下DFIG的轉(zhuǎn)子電壓方程[14]

進(jìn)一步推導(dǎo)出轉(zhuǎn)子轉(zhuǎn)速ωr、定子側(cè)有功功率Ps和無功功率Qs的表達(dá)式為

式中，np為電機極對數(shù)；ωs為同步電速度；J為電機轉(zhuǎn)動慣量；M為與轉(zhuǎn)速成正比的轉(zhuǎn)矩阻尼系數(shù)。式（3）～式（5）即構(gòu)成了完整的DFIG數(shù)學(xué)模型。

3 控制策略

3.1 WG-FPC輸出功率參考值的確定

文獻(xiàn)[16]提出采用提前 15min的短期預(yù)測風(fēng)速計算出的預(yù)測功率作為參考值，這種方法實現(xiàn)的前提是飛輪儲能容量和功率足夠大，且需要由上位機實時向風(fēng)電機組發(fā)出功率指令，因此在實際應(yīng)用中很難實現(xiàn)。為了抑制有功功率波動，同時維持飛輪轉(zhuǎn)速在允許的范圍內(nèi)，本文采用模糊推理的方法確定WG-FPC輸出功率的參考值[8]。

模糊推理的輸入為風(fēng)力發(fā)電機的實際輸出功率Pwgf和飛輪的轉(zhuǎn)速Ω，輸出為 WG-FPC的輸出功率指令值 Preg，原理如圖 2所示。輸入變量的隸屬度函數(shù)如圖3所示，S、M、B分別表示小、中、大三個模糊子集；輸出變量的隸屬度函數(shù)如圖4所示，VS、S、SM、M、BM、B、VB分別表示很小、小、小中、中、大中、大、很大6個模糊子集。

圖2 模糊推理原理圖Fig.2 Block diagram of the fuzzy logic inference

圖3 輸入變量的隸屬函數(shù)Fig.3 Membership functions of input variables

圖4 輸出變量的隸屬函數(shù)Fig.4 Membership function of output variable

模糊推理的結(jié)果決定了 FPC的儲能或發(fā)電狀態(tài)，推理的基本原則如下：如果飛輪轉(zhuǎn)速太低，則優(yōu)先讓FPC儲能，風(fēng)力發(fā)電機的大部分輸出功率用于對FPC充電，少量注入電網(wǎng)；如果轉(zhuǎn)速太高，則FPC優(yōu)先放電，風(fēng)力發(fā)電機的大部分輸出功率注入電網(wǎng)，少量用于對FPC充電；如果飛輪轉(zhuǎn)速適中，則 FPC按正常速度運轉(zhuǎn)，系統(tǒng)注入電網(wǎng)的功率為Pwgf。根據(jù)上述基本原則，模糊推理規(guī)則見表1，推理結(jié)果如圖5所示。

表1 模糊推理表Tab.1 Fuzzy logic inference table

圖5 模糊推理結(jié)果Fig.5 Result of fuzzy logic inference

3.2 GSC控制

GSC的有功功率控制是通過對直流母線電容電壓的控制實現(xiàn)的。根據(jù)圖1中假設(shè)的電流和功率方向，從功率的角度看，對于直流母線電容電壓VC，有如下等式

式中，Pwg、Pr、Pf分別為風(fēng)力發(fā)電機輸出的有功功率、FPC轉(zhuǎn)子側(cè)有功功率和GSC的有功功率；C為直流母線電容。

忽略濾波器電阻Rf，并根據(jù)式（1）和式（2）的假設(shè)，則有[17]

式中，vfd、vfq分別是GSC交流側(cè)電壓的d、q分量；ifd、ifq是濾波器電流的d、q分量。進(jìn)一步可推導(dǎo)出GSC與電網(wǎng)交換的有功功率 Pf和無功功率 Qf分別為

由式（8）可知，可以分別通過對vfq和vfd的控制，實現(xiàn)對 Pf和 Qf的解耦控制。結(jié)合式（6），直流母線電容電壓VC可以由vfd控制，控制框圖如圖6所示。

圖6 直流母線電壓控制框圖Fig.6 Block diagram of DC bus voltage controller

無功功率Qf的控制框圖如圖7所示?？紤]GSC的容量問題，采用單位功率因數(shù)控制，即給定值Qf_ref為 0。

圖7 GSC無功功率控制框圖Fig.7 Block diagram of GSC reactive power controller

3.3 FPC控制

FPC按功率、電流雙閉環(huán)控制。本文采用基于定子磁鏈定向的PWM矢量控制方法，實現(xiàn)FPC的有功功率和無功功率的解耦控制。具體過程已有文獻(xiàn)討論[12-14]，此處不再贅述。下面對FPC的有功、無功功率的給定值確定方法進(jìn)行討論。

根據(jù)圖1所示的功率關(guān)系和3.1節(jié)的討論易知，F(xiàn)PC的有功功率給定值Pref為

為了確定FPC的無功功率給定值，首先分析機端電壓穩(wěn)定性的影響因素。由圖1可得

式中，Pg、Qg分別為WG-FPC注入電網(wǎng)的有功和無功功率；Rl、Ll為等值線路的阻抗參數(shù)；是電網(wǎng)電壓。假設(shè)是恒定的，且相角為零（即以為參考相量），則有

式中，Vg是電網(wǎng)電壓的幅值。從而，機端電壓 Vs可表示為

一般情況下，有

因此，機端電壓表達(dá)式可簡化為

由式（14）可知，在假設(shè)無窮大系統(tǒng)電壓 Vg恒定，且等值線路的參數(shù)Rl、Ll確定時，機端電壓水平完全由WG-FPC系統(tǒng)的有功、無功功率決定。在有功功率波動的情況下，若要保持機端電壓穩(wěn)定，必須從電網(wǎng)吸收無功功率進(jìn)行補償，控制原理如圖8所示。為了減少無功功率的傳輸，由FPC提供需要的無功功率，F(xiàn)PC無功功率參考值Qref數(shù)值上等于用于機端電壓控制的無功功率期望值

圖8 機端電壓控制框圖Fig.8 Block diagram of machine terminal voltage controller

4 FPC儲能容量討論

本文中的FPC主要用于抑制風(fēng)力發(fā)電機輸出功率波動，儲能容量的大小直接影響到輸出功率特性。對于風(fēng)電功率中的低頻分量，由于其波動比較緩慢，功率變化率較小，注入電網(wǎng)時電力系統(tǒng)有充足的時間進(jìn)行響應(yīng)，然而當(dāng)高頻分量與其疊加后，導(dǎo)致功率變化率變大[18]，短時間內(nèi)對電網(wǎng)造成嚴(yán)重的沖擊，給電力系統(tǒng)安全運行帶來隱患。因此，抑制風(fēng)電功率波動的目標(biāo)是要抑制風(fēng)電功率中的高頻成分，減小風(fēng)電功率的變化率，為電力系統(tǒng)提供較為穩(wěn)定的功率輸出。文獻(xiàn)[19]研究表明，風(fēng)力發(fā)電的0.01～1Hz頻段內(nèi)的波動功率對電能質(zhì)量的影響最大，而1Hz以上成分主要被風(fēng)力渦輪機的慣性所吸收。故若不考慮FPC的功率限制，認(rèn)為Pwg中的0.01Hz以上的波動功率全部被FPC吸收，則FPC的瞬時功率為

式中，F(xiàn)(s)為一階高通濾波器，其穿越頻率fc=0.01Hz，時間常數(shù)τ =16s。FPC吸收的能量Efpc可表示為功率Pfpc的時間積分，則有

利用卷積方法，可將式（16）變換到時域范圍內(nèi)，則有

因為風(fēng)力發(fā)電機的輸出功率Pwg介于0到額定功率Pwgr之間，且，則有

相關(guān)研究表明，在風(fēng)電機組正常運行情況下，使用的儲能容量還不到上限τPwgr的一半[19]。故可取FPC的儲能容量為

式中，cs的值介于0和1之間，當(dāng)τ =16s時，cs可約取0.3。為了保證極端運行條件下，F(xiàn)PC不會超出其儲能容量，可采用模糊推理的方法得到FPC的瞬時功率（見3.1節(jié)所述）。另外，F(xiàn)PC的額定功率也與濾波時間常數(shù)有關(guān)，假設(shè)FPC按額定功率Pfpcr充電，其轉(zhuǎn)速從額定轉(zhuǎn)速Ωn升高到最高轉(zhuǎn)速Ωmax時，則有成立[20]。

5 仿真實驗及結(jié)果分析

5.1 WG-FPC仿真模型

圖9 WG-FPC運行控制原理框圖Fig.9 Block diagram of operation control for WG-FPC

根據(jù)上述分析結(jié)果，在Matlab/Simulink環(huán)境下建立WG-FPC運行控制仿真模型，如圖9所示。仿真中使用的FPC參數(shù)來自實驗室研制的樣機[14]；風(fēng)力發(fā)電機額定功率為100kW；其他相關(guān)參數(shù)：電容C=2 200×5μF，直流母線電壓VC=800V，濾波器電阻Rf=0、電感Lf=0.01H，等值線路電阻Rl=0.682Ω、電感Ll=0.019H，線端電壓Vs=Vg=690V。需要指出的是，在有功功率控制中，Pref是 FPC的總有功功率給定值，由于轉(zhuǎn)子通過RSC與其他部分存在有功功率交換，因此在有功指令中加入了轉(zhuǎn)子有功補償[13]。

5.2 仿真結(jié)果分析

本文所有仿真結(jié)果是在圖10所示的實際測量風(fēng)速（來自于文獻(xiàn)[8]）下進(jìn)行的，采用文獻(xiàn)[1]中所述的風(fēng)速-功率特性曲線求取風(fēng)力發(fā)電機實時輸出功率，仿真周期為250s，仿真結(jié)果如圖11所示。

圖10 實測風(fēng)速曲線Fig.10 Field-measured wind speed profile

圖11 實測風(fēng)速下WG-FPC的動態(tài)響應(yīng)曲線Fig.11 Dynamic responses of WG-FPC under field-measured wind speed condition

圖11a表明，通過直流母線電壓控制，能很好地將直流母線電壓穩(wěn)定在給定值上。圖11b是飛輪轉(zhuǎn)速曲線，表明FPC在抑制風(fēng)力發(fā)電機功率波動過程中的充、放電情況，且與圖10對照容易看出：當(dāng)風(fēng)速較高時，飛輪充電，轉(zhuǎn)速升高，當(dāng)風(fēng)速較低時，飛輪放電，轉(zhuǎn)速降低，從而達(dá)到平抑功率波動的目的，而且飛輪轉(zhuǎn)速始終維持在限定范圍內(nèi)。

圖11c表明在沒有FPC的情況下，按照最大功率追蹤控制的風(fēng)力發(fā)電機向電網(wǎng)注入的功率隨風(fēng)速的波動而發(fā)生劇烈波動，而接入FPC可以較好地抑制功率的波動，改善了風(fēng)力發(fā)電機輸出功率特性。FPC輸出的調(diào)節(jié)功率如圖11d所示。

圖11e是機端電壓曲線。由圖可以看出，沒有接入FPC時，機端電壓隨風(fēng)力發(fā)電機的輸出功率波動而發(fā)生波動；在接入FPC，但沒有對FPC的無功功率進(jìn)行控制的情況下，由于有功的波動得到抑制，機端電壓波動相對變小，但仍有一定程度的波動；而在FPC對無功進(jìn)行補償?shù)那闆r下，機端電壓可以很好地穩(wěn)定在給定值上，從而很大程度地改善了風(fēng)電系統(tǒng)的電壓穩(wěn)定性。FPC輸出的補償無功功率如圖11f所示。以上仿真結(jié)果表明，應(yīng)用FPC與風(fēng)力發(fā)電機聯(lián)合運行，可以很好地改善風(fēng)力發(fā)電機的輸出功率特性。

6 結(jié)論

本文提出的風(fēng)力發(fā)電機與FPC的聯(lián)合運行系統(tǒng)，能夠根據(jù)風(fēng)力發(fā)電機的輸出功率和飛輪轉(zhuǎn)速自動調(diào)節(jié)注入電網(wǎng)的功率，具有很強的適應(yīng)性。系統(tǒng)中的FPC可以直接參與風(fēng)力發(fā)電機的有功和無功調(diào)節(jié)，在平抑風(fēng)力發(fā)電機輸出功率波動的同時，增強了其參與電力系統(tǒng)的調(diào)壓、調(diào)頻控制能力，從而使隨機性很強的風(fēng)電場具有類似于傳統(tǒng)電源的可調(diào)度性，從根本上解決風(fēng)力發(fā)電并網(wǎng)困難問題。由于篇幅有限，WG-FPC參與系統(tǒng)調(diào)頻控制、變換器GSC與FPC的無功協(xié)調(diào)控制以及儲能容量對并網(wǎng)運行特性的影響分析等問題將另文討論。

[1]韓濤, 盧繼平, 喬梁, 等.大型并網(wǎng)風(fēng)電場儲能容量優(yōu)化方案[J]. 電網(wǎng)技術(shù), 2010, 34(1): 169-173.Han Tao, Lu Jiping, Qiao Liang, et al.Optimized scheme of energy storage capacity for grid connected large scale wind farm[J]. Power System Technology,2010, 34(1): 169-173.

[2]Luo C, Banakarh, Shen B, et al. Strategies to smooth wind power fluctuations of wind turbine generator[J].IEEE Transactions on Energy Conversion, 2007, 22(2): 341-349.

[3]Chong H, Huang A Q, Baran M E, et al. STATCOM impact study on the integration of a large wind farm into a weak loop power system[J]. IEEE Transactions on Energy Conversion, 2008, 23(1): 226-233.

[4]Lu M S, Chang C L, Lee W R, et al. Combining the wind power generation system with energy storage equipment[J]. IEEE Transactions on Industry Applications, 2009, 45(6): 2109-2115.

[5]Chad A, Kai S, Geza J.A knowledge-based approach for control of two-level energy storage for wind energy system[J]. IEEE Transactions on Energy Conversion, 2009, 24(2): 539-547.

[6]李霄, 胡長生, 劉昌金, 等. 基于超級電容儲能的風(fēng)電場功率調(diào)節(jié)系統(tǒng)建模與控制[J]. 電力系統(tǒng)自動化, 2009, 33(9): 86-90.Li Xiao, Hu Changsheng, Liu Changjin, et al. Design of SMES control system for smoothing power fluctuations in wind farms[J]. Automation of Electric Power Systems, 2009, 33(9): 86-90.

[7]Lilia J, Lotfi K, Abderrazak O. A fuzzy logic supervisor for active and reactive power control of a variable speed wind energy conversion system associated to a flywheel storage system[J]. Electric Power Systems Research, 2009, 79: 919-925.

[8]Gabriel O C, Christophe S, Benoit R, et al. Control and performance evaluation of a flywheel energystorage system associated to a variable-speed wind generator[J]. IEEE Transactions on Industrial Electronics, 2006, 53(4): 1074-1085.

[9]張文亮, 丘明,來小康.儲能技術(shù)在電力系統(tǒng)中的應(yīng)用[J]. 電網(wǎng)技術(shù), 2008, 32(7): 1-9.Zhang Wenliang, Qiu Ming, Lai Xiaokang.Application of energy storage technologies in power grids[J]. Power System Technology, 2008, 32(7): 1-9.

[10]John P B, David G I. Energy storage and its use with intermittent renewable energy[J]. IEEE Transactions on Energy Conversion, 2004, 19(2): 441-448.

[11]文勁宇, 程時杰, 李剛, 等. 一種增強電力系統(tǒng)穩(wěn)定性的多功能柔性功率調(diào)節(jié)器[J]. 中國電機工程學(xué)報, 2005, 25(25): 6-11.Wen Jinyu, Cheng Shijie, Li Gang, et al. A multifunctional flexible power conditioner for power system stabilities enhancement[J]. Proceedings of the CSEE, 2005, 25(25): 6-11.

[12]鄒旭東, 劉新民, 段善旭, 等.儲能調(diào)相功率調(diào)制系統(tǒng)柔性功率調(diào)節(jié)器[J]. 電工技術(shù)學(xué)報, 2009, 24(6): 146-153.Zou Xudong, Liu Xinmin, Duan Shanxu, et al.Flexible power conditioner of energy storage power modulation system[J]. Transactions of China Electrotechnical Society, 2009, 24(6): 146-153.

[13]李剛, 程時杰, 文勁宇, 等.利用柔性功率調(diào)節(jié)器提高電力系統(tǒng)穩(wěn)定性[J]. 中國電機工程學(xué)報, 2006, 26(23): 1-6.Li Gang, Cheng Shijie, Wen Jinyu, et al. Power system stability enhancement by the flexible power conditioner[J]. Proceedings of the CSEE, 2006, 26(23): 1-6.

[14]楊浩, 文勁宇, 李剛, 等. 多功能柔性功率調(diào)節(jié)器運行特性的仿真研究[J]. 中國電機工程學(xué)報, 2006,26(2): 19-24.Yang Hao, Wen Jinyu, Li Gang, et al. Investigation on operation characteristics of multi-functional flexible power conditioner[J]. Proceedings of the CSEE, 2006, 26(2):19-24.

[15]Hirofumi A, Hikaru S.Control and performance of a doubly-fed induction machine intended for a flywheel energy storage system[J]. IEEE Transactions on Power Electronics, 2002, 17(1): 109-116.

[16]孫春順, 王耀南, 李欣然. 飛輪輔助的風(fēng)力發(fā)電系統(tǒng)功率和頻率綜合控制[J]. 中國電機工程學(xué)報,2008, 28(29): 111-116.Sun Chunshun, Wang Yaonan, Li Xinran. Synthesized power and frequency control of wind power generation system assisted through flywheels[J].Proceedings of the CSEE, 2008, 28(29): 111-116.

[17]陳堅. 電力電子學(xué)—電力電子變換和控制技術(shù)[M].北京: 高等教育出版社, 2004.

[18]John P B, David G I. Energy storage and its use with intermittent renewable energy[J]. IEEE Transactions on Energy Conversion, 2004, 19(2): 441-448.

[19]Li W, Joos G, Abbey C. Wind power impact on system frequency deviation and an ESS based power filtering algorithm solution[C]. Proceedings of IEEE PES Power Systems Conference and Exposition,Atlanta, GA, USA, 2006: 2077-2084.

[20]Gabriel C, Stefan B, Miecea M R, et al. Design and control strategies of an induction-machine-based flywheel energy storage system associated to a variable-speed wind generator[J]. IEEE Transactions on Energy Conversion, 2010, 25(2): 526-534.