并聯(lián)型有源電力濾波器的復(fù)合電流控制方法

田飛燕，高云廣，2，宋建成，鄭麗君，楊健康，康勤根

（1.煤礦電氣設(shè)備與智能控制山西省重點(diǎn)實(shí)驗(yàn)室（太原理工大學(xué)），山西 太原 030024；2.太原科技大學(xué)電子與信息工程學(xué)院，山西 太原 030024；3.山西晉城無(wú)煙煤礦業(yè)集團(tuán)有限責(zé)任公司，山西 晉城 048000）

并聯(lián)型有源電力濾波器的復(fù)合電流控制方法

田飛燕1，高云廣1，2，宋建成1，鄭麗君1，楊健康3，康勤根3

（1.煤礦電氣設(shè)備與智能控制山西省重點(diǎn)實(shí)驗(yàn)室（太原理工大學(xué)），山西 太原 030024；2.太原科技大學(xué)電子與信息工程學(xué)院，山西 太原 030024；3.山西晉城無(wú)煙煤礦業(yè)集團(tuán)有限責(zé)任公司，山西 晉城 048000）

針對(duì)非線性負(fù)荷對(duì)電網(wǎng)電能質(zhì)量的影響，研究一種并聯(lián)型有源電力濾波器拓?fù)浣Y(jié)構(gòu)。根據(jù)此拓?fù)浣Y(jié)構(gòu)，建立其在三相靜止坐標(biāo)系下和兩相同步旋轉(zhuǎn)坐標(biāo)系下的數(shù)學(xué)模型，推導(dǎo)模型的傳遞函數(shù)，設(shè)計(jì)基于PΙ控制的電流控制器。通過(guò)分析PΙ控制器在有源電力濾波器電流跟蹤控制中存在的局限性，提出將PΙ控制與重復(fù)控制并聯(lián)的復(fù)合控制策略，以改善系統(tǒng)的穩(wěn)態(tài)性能。通過(guò)實(shí)驗(yàn)對(duì)所提出的控制方法進(jìn)行研究，結(jié)果表明該復(fù)合控制方法具有很高的穩(wěn)態(tài)補(bǔ)償精度，能夠有效地抑制由非線性負(fù)載引起的諧波電流，大大改善電網(wǎng)電流波形。

諧波治理；有源電力濾波器；PΙ控制；重復(fù)控制；復(fù)合控制

隨著電力電子裝置的普遍應(yīng)用，大量的諧波和無(wú)功電流注入到電網(wǎng)中［1］，導(dǎo)致各種電力設(shè)備無(wú)法正常工作，甚至?xí)鹪O(shè)備損壞和事故發(fā)生。因此諧波治理技術(shù)的研究被廣泛關(guān)注。

目前，有源電力濾波技術(shù)受到廣泛的研究。有源電力濾波器（active power filter，APF）能夠抵消電網(wǎng)中由非線性負(fù)載引起的諧波和無(wú)功電流，改善電網(wǎng)電流波形。與無(wú)源濾波器相比，APF占用空間小，不受電網(wǎng)阻抗的影響，能夠有效跟蹤頻率和幅值均變化的諧波電流，并且還可以消除濾波器阻抗和網(wǎng)絡(luò)阻抗之間可能存在的諧振［2-4］。

有源電力濾波器是通過(guò)控制逆變器的輸出電流來(lái)快速地跟蹤由非線性負(fù)載引起的諧波電流，實(shí)現(xiàn)對(duì)諧波電流的有效補(bǔ)償。在此過(guò)程中，要使逆變器的輸出電流快速精確地跟蹤上諧波電流的變化，關(guān)鍵在于電流跟蹤控制算法的實(shí)時(shí)性和有效性。電流跟蹤控制方法主要有：滯環(huán)比較控制、三角載波控制、無(wú)差拍控制等。其中三角載波控制較為簡(jiǎn)單易行，且動(dòng)態(tài)響應(yīng)好、開(kāi)關(guān)頻率固定，其將電流的偏差經(jīng)放大器后再與三角波比較，比較所得的矩形脈沖作為控制信號(hào)來(lái)控制逆變器的開(kāi)關(guān)元件。放大器通常采用PΙ控制器，而PΙ控制只有當(dāng)給定是直流量時(shí)才能實(shí)現(xiàn)系統(tǒng)的無(wú)靜差跟蹤［5-7］。因?yàn)橛性措娏V波器控制的給定值是含有多次諧波的交流量，因此，單獨(dú)PΙ控制無(wú)法對(duì)諧波信號(hào)實(shí)現(xiàn)無(wú)靜差跟蹤。

本文在傳統(tǒng)PΙ控制的基礎(chǔ)上，結(jié)合重復(fù)控制思想，提出將PΙ控制和重復(fù)控制相結(jié)合的電流控制方法，以提高APF的補(bǔ)償性能，降低電網(wǎng)電流的總諧波失真（total harmonics distortion，THD）。

1 APF數(shù)學(xué)模型及控制系統(tǒng)

系統(tǒng)主要有3部分：電源、非線性負(fù)載和APF。其中APF包括：電壓源逆變器、直流母線電容和三相輸出電感。圖1所示為三相并聯(lián)型有源電力濾波器的結(jié)構(gòu)原理圖，主電路為三相電壓型PWM變換器，非線性負(fù)載為帶阻性負(fù)載的三相不控整流橋。

圖1 SAPF系統(tǒng)原理圖Fig.1 The system schematic of three-phase SAPF

1.1 APF數(shù)學(xué)模型

由圖1可得，三相有源電力濾波器在三相坐標(biāo)系下的數(shù)學(xué)模型為

式中：usa,usb,usc為電網(wǎng)電壓；ura,urb,urc為PWM變換器輸出電壓；ica,icb,icc為補(bǔ)償電流；r為輸出電感及開(kāi)關(guān)損耗的等效電阻。

為了便于對(duì)有功和無(wú)功分量進(jìn)行獨(dú)立控制，將式（1）通過(guò)abc/dq坐標(biāo)變換，轉(zhuǎn)換到d-q坐標(biāo)系中，下式即為APF在d-q坐標(biāo)系中的數(shù)學(xué)模型：

根據(jù)此模型即可以進(jìn)行PΙ控制器的設(shè)計(jì)。

1.2 APF的控制系統(tǒng)

在有源濾波器的控制中，不僅需要對(duì)電流進(jìn)行快速地跟蹤和補(bǔ)償，而且還需對(duì)PWM變換器直流側(cè)電壓進(jìn)行穩(wěn)定控制，以維持系統(tǒng)特有的穩(wěn)定性。根據(jù)補(bǔ)償電流快速性要優(yōu)于直流側(cè)電壓控制的特點(diǎn)，采用電壓外環(huán)電流內(nèi)環(huán)的控制方法。

圖2是SAPF控制結(jié)構(gòu)圖，電壓環(huán)采用PΙ控制，其作用是保證APF的直流側(cè)電壓穩(wěn)定在給定值，使電流內(nèi)環(huán)能夠有效地補(bǔ)償諧波［8-11］。本文所研究的電流環(huán)控制器為PΙ控制器和復(fù)合控制器，電流環(huán)的作用是控制APF輸出與給定諧波電流大小相等方向相反的電流。

圖2 SAPF控制結(jié)構(gòu)圖Fig.2 Control structure diagram of SAPF

2 傳統(tǒng)PI控制方法

由式（2）可知，d，q軸變量之間相互耦合，為了簡(jiǎn)化控制器設(shè)計(jì)，可采用帶有電壓前饋和輸出電流交叉解耦項(xiàng)的PΙ控制器。urd,urq的控制方程如下：

將式（3）帶入式（2）可得：

式中：KPi為電流內(nèi)環(huán)比例系數(shù)；KIi為積分系數(shù)；i*d,i*q分別為id,iq的諧波電流給定值。

從式（4）可以看出，通過(guò)控制i*d-id和i*q-iq可以對(duì)id,iq進(jìn)行獨(dú)立控制。由式（4）可知，控制對(duì)象實(shí)際上是一階慣性環(huán)節(jié)。圖3為采用PΙ控制器對(duì)電流進(jìn)行調(diào)節(jié)時(shí)，電流環(huán)的簡(jiǎn)化控制框圖。

圖3 d-q坐標(biāo)系下的電流環(huán)簡(jiǎn)化控制圖Fig.3 Simplified current control block diagram of APF in d-q frame

根據(jù)圖3可知，電流閉環(huán)控制系統(tǒng)傳遞函數(shù)為

采用零極點(diǎn)對(duì)消法，即令

則式（5）可以簡(jiǎn)化為

式中：Td為d，q軸電流閉環(huán)系統(tǒng)電流響應(yīng)延時(shí)時(shí)間常數(shù)。

由式（8）可知，Td與電流控制器的比例系數(shù)KPi成反比，與輸出電感L成正比。當(dāng)電感一定的情況下，若增大KPi,則延時(shí)時(shí)間會(huì)減小。但是KPi不能無(wú)限制地增大，過(guò)大的KPi會(huì)導(dǎo)致系統(tǒng)振蕩。單獨(dú)的PΙ控制在控制諧波電流時(shí)受到帶寬范圍的限制，難以消除指令延遲的影響，補(bǔ)償效果不理想。

3 APF復(fù)合控制器設(shè)計(jì)

針對(duì)有源電力濾波器電流環(huán)單獨(dú)PΙ控制補(bǔ)償性能有限，本文提出基于PΙ控制和重復(fù)控制的復(fù)合控制系統(tǒng)來(lái)滿(mǎn)足系統(tǒng)對(duì)穩(wěn)態(tài)性能的要求。重復(fù)控制的基本思想源于控制理論中的內(nèi)模原理，即把系統(tǒng)外部信號(hào)的動(dòng)態(tài)模型植入控制器內(nèi)，在穩(wěn)定的閉環(huán)系統(tǒng)中包含外部輸入信號(hào)的數(shù)學(xué)模型，以構(gòu)成高精度的反饋控制系統(tǒng)［11-13］。

在有源電力濾波器中，給定信號(hào)為檢測(cè)到的諧波交流信號(hào)，其由許多不同頻率的交流量疊加而成。若要實(shí)現(xiàn)無(wú)差跟蹤，則要求控制模型中的內(nèi)模能夠構(gòu)建出每一種交流信號(hào)。雖然檢測(cè)到的諧波信號(hào)含有很多不同頻率的交流量，但是在每一個(gè)基波周期諧波信號(hào)的波形都是重復(fù)出現(xiàn)的。因此，可以將基波周期作為這些諧波信號(hào)的重復(fù)周期。

3.1 APF重復(fù)控制器的內(nèi)模

圖4為APF重復(fù)器內(nèi)模，其中N=T/Ts，T為基波周期，Ts為控制器離散過(guò)程中的采樣周期，z-N為純延時(shí)環(huán)節(jié)。

圖4 重復(fù)信號(hào)發(fā)生器內(nèi)模模型Fig.4 Internal model of repetitive signal generator

從圖4可以得到y(tǒng)(z)和e(z)的關(guān)系，如下式：

由式（9）可以看出，每隔1個(gè)周期（N步），則將上周期輸出量的削弱值與當(dāng)前的輸入量累加，該累加值作為當(dāng)前周期的輸出。根據(jù)諧波的周期性特點(diǎn)，采用重復(fù)控制可以提高跟蹤信號(hào)的準(zhǔn)確性和穩(wěn)定性。

3.2 APF重復(fù)控制器結(jié)構(gòu)

圖5是重復(fù)控制器結(jié)構(gòu)，其中P(z)為控制對(duì)象，要求其本身是穩(wěn)定的；Q(z)和C(z)分別為重復(fù)控制輔助補(bǔ)償器和針對(duì)受控對(duì)象的補(bǔ)償器，C(z)=KrzkS(z)，Kr為重要控制器增益，k為相位補(bǔ)償系數(shù)；r為需要跟蹤的重復(fù)性指令。重復(fù)控制器檢測(cè)指令r與實(shí)際輸出y之間的誤差e，內(nèi)模對(duì)誤差e進(jìn)行周期性積分，以便在誤差消失或變得較小時(shí)控制器仍能輸出相應(yīng)的控制量。在周期延遲環(huán)節(jié)和補(bǔ)償器的作用下，根據(jù)內(nèi)模積分的結(jié)果，重復(fù)控制器將在下一周期的適當(dāng)時(shí)刻輸出控制量ur，以減小誤差e。

圖5 重復(fù)控制器結(jié)構(gòu)Fig.5 The structure of the repetitive controller

3.3 APF復(fù)合控制策略

重復(fù)控制器對(duì)誤差信號(hào)是以基波周期為步長(zhǎng)進(jìn)行的累加，它無(wú)法將動(dòng)態(tài)響應(yīng)時(shí)間減小到1個(gè)基波周期內(nèi)。而PΙ控制策略動(dòng)態(tài)響應(yīng)速度很快，因此將傳統(tǒng)PΙ控制策略和重復(fù)控制策略并聯(lián)使用，構(gòu)成新的控制系統(tǒng)，同時(shí)兼顧系統(tǒng)的穩(wěn)態(tài)性能和動(dòng)態(tài)性能。復(fù)合控制框圖如圖6所示。

圖6中，GPI(z)是PΙ控制器，通過(guò)對(duì)控制對(duì)象進(jìn)行頻域分析并校正，可得到PΙ控制下穩(wěn)定運(yùn)行參數(shù)。在PΙ控制的穩(wěn)定系統(tǒng)下，重復(fù)控制器的參數(shù)具體計(jì)算步驟如下：

圖6 PI控制和重復(fù)控制的并聯(lián)式復(fù)合控制框圖Fig.6 Parallel compound control block diagram of the PI control and repetitive control

1）計(jì)算1個(gè)周期的采樣次數(shù)N；

2）選取Q(z)。Q(z)=1時(shí)，重復(fù)控制系統(tǒng)的穩(wěn)態(tài)誤差可達(dá)到零，但是要使系統(tǒng)保持穩(wěn)定，則必須有Q(z)＜1?？紤]穩(wěn)定性和穩(wěn)態(tài)誤差兩方面的要求，Q(z)取常數(shù)為0.98；

3）S(z)由控制對(duì)象的幅頻特性選擇，一般補(bǔ)償器S(z)可選二階低通濾波器。當(dāng)補(bǔ)償50次以?xún)?nèi)的諧波電流時(shí)，截止頻率設(shè)為fc=2 800 Hz，則所設(shè)計(jì)的二階低通濾波器：

4）為了簡(jiǎn)化，設(shè)計(jì)k估算值為3即z3，主要作用是補(bǔ)償S(z)和控制對(duì)象在中低頻段的相位滯后；

5）Kr越小系統(tǒng)穩(wěn)定性越好，但收斂速度變慢且穩(wěn)態(tài)誤差上升，Kr通過(guò)實(shí)驗(yàn)選取為0.7。

4 實(shí)驗(yàn)結(jié)果

為了驗(yàn)證上述控制策略的有效性，在實(shí)驗(yàn)室搭建了APF實(shí)驗(yàn)平臺(tái)進(jìn)行特性實(shí)驗(yàn)。實(shí)驗(yàn)系統(tǒng)主要包括3部分：電源系統(tǒng)、三相二極管整流裝置和并聯(lián)型有源電力濾波器。不控橋負(fù)載電阻為20 Ω，并聯(lián)型有源電力濾波器的系統(tǒng)組成框圖如圖7所示。

圖7 并聯(lián)型有源電力濾波器系統(tǒng)組成框圖Fig.7 Diagram of shunt active power filter system

主電路ΙGBT參數(shù)為100 A/1 200 V，控制器DSP采用TΙ公司的TMS320F28335，諧波測(cè)量?jī)x器采用Fluke公司的Fluke434電能質(zhì)量分析儀，HDO6104示波器以及Tektronix P5210電壓測(cè)量探頭和HOKΙ 3274電流測(cè)量探頭等。APF實(shí)驗(yàn)主要參數(shù)為：電源相電壓110 V，電網(wǎng)頻率50 Hz，APF輸出電感2 mH，APF直流側(cè)電容3 300 μF，APF直流側(cè)電壓110 V，開(kāi)關(guān)頻率9 kHz。

4.1 實(shí)驗(yàn)波形

圖8為負(fù)載（帶阻性負(fù)載的三相不控整流橋）電流波形頻譜分布情況，單獨(dú)采用PΙ控制器時(shí)得到的電網(wǎng)電流頻譜和波形如圖9和圖10所示。

圖8 負(fù)載電流波形頻譜圖Fig.8 Spectrum waveform of load current

圖9 PI控制補(bǔ)償后電網(wǎng)電流頻譜圖Fig.9 Spectrum waveform of grid current with PI controller

圖10 單獨(dú)PI控制時(shí)電流波形Fig.10 Current waveforms with PI controller

采用PΙ與重復(fù)控制相結(jié)合的復(fù)合控制策略對(duì)電流進(jìn)行跟蹤控制時(shí)，電網(wǎng)電流諧波分布及頻譜如表1和圖11所示，電流波形如圖12所示。

表1 復(fù)合控制補(bǔ)償后電網(wǎng)電流諧波分布Tab.1 The harmonic distribution of grid current with compound controller

圖11 復(fù)合控制補(bǔ)償后電網(wǎng)電流頻譜Fig.11 Spectrum waveform of grid current with compound control

圖12 復(fù)合控制時(shí)電流波形Fig.12 Current waveforms with compound control

4.2 實(shí)驗(yàn)結(jié)果分析

從圖9可看到，采用PΙ控制進(jìn)行補(bǔ)償后，電網(wǎng)電流THD從28.6%降低到了14.0%，補(bǔ)償效果不理想。這是由于在系統(tǒng)穩(wěn)定的前提下PΙ控制的補(bǔ)償能力有限，即PΙ控制器無(wú)法對(duì)交流信號(hào)實(shí)現(xiàn)無(wú)差調(diào)節(jié)。圖10中補(bǔ)償后的網(wǎng)側(cè)電流波形有周期性尖脈沖出現(xiàn)，其原因是三相不控整流在換相時(shí)，一相電流瞬時(shí)發(fā)生跳變，在陡變的地方，電流頻譜較為復(fù)雜，含有不確定的高頻分量，導(dǎo)致檢測(cè)到的電流指令超出電流環(huán)PΙ調(diào)節(jié)器能夠控制的帶寬范圍。此外，數(shù)字控制中總有一拍延遲，APF無(wú)法實(shí)時(shí)地輸出需要的諧波以抵消這些陡變量，造成實(shí)際系統(tǒng)電流波形上出現(xiàn)了這些尖脈沖。從圖11和圖12可以看出，采用PΙ+重復(fù)控制算法進(jìn)行補(bǔ)償后，網(wǎng)側(cè)電流波形有明顯改善，電流THD值降低到了3.3%。從圖12可以看出，加入重復(fù)控制器后，電網(wǎng)電流波形的正弦度大大改善，在固定位置上的尖沖也得到了有效抑制，這是因?yàn)橹貜?fù)控制器實(shí)際上是一種以基波周期為步長(zhǎng)對(duì)誤差信號(hào)進(jìn)行累加的積分控制，通過(guò)對(duì)波形誤差進(jìn)行逐周期的補(bǔ)償，即可達(dá)到很高的穩(wěn)態(tài)精度。

5 結(jié)論

本文針對(duì)非線性負(fù)荷對(duì)電網(wǎng)質(zhì)量的影響，對(duì)非線性負(fù)荷特性、并聯(lián)型有源濾波器的拓?fù)浣Y(jié)構(gòu)、數(shù)學(xué)模型、控制策略和特性試驗(yàn)進(jìn)行了系統(tǒng)研究，研究結(jié)論如下：1）提出的將PΙ控制與重復(fù)控制并聯(lián)使用的復(fù)合控制策略，有效解決了單獨(dú)采用PΙ控制時(shí)補(bǔ)償精度不高的問(wèn)題，提高了系統(tǒng)穩(wěn)態(tài)精度，實(shí)現(xiàn)了實(shí)時(shí)準(zhǔn)確跟蹤補(bǔ)償電流的目的；2）設(shè)計(jì)了基于復(fù)合控制的并聯(lián)型有源電力濾波器，有效抑制了由非線性負(fù)載引起的諧波電流。采用APF進(jìn)行補(bǔ)償后，電網(wǎng)電流波形正弦度大大提高，近似于正弦波，電流THD減小到5%以下，達(dá)到了諧波抑制的國(guó)家標(biāo)準(zhǔn)。

［1］王兆安，楊君，劉進(jìn)軍.諧波抑制和無(wú)功功率補(bǔ)償［M］.北京：機(jī)械工業(yè)出版社，2006.

［2］魏學(xué)良，戴珂，方昕，等.三相并聯(lián)型有源電力濾波器補(bǔ)償電流性能分析與改進(jìn)［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2007，27（28）：113-119.

［3］Chen B，Zha X，Gong J，et al.Realization and Improvement of Repetitive Control in Rotating Frame for Active Power Filter System［J］.Applied Power Electronics and Exposion，2010：887-894.

［4］楊新華，謝斌.基于電流環(huán)復(fù)合控制的有源電力濾波器研究［J］.電氣傳動(dòng)，2014，44（2）：55-58.

［5］Rshmani S，Hamadi A，AL-Haddad K.A Lyapunov-func?tion-based Control for a Three-phase Shunt Hybrid Active Fil?ter［J］.IEEE Transactions on Industrial Electronics，2012，59（3）：1418-1429.

［6］郭偉峰，武健，徐殿國(guó)，等.新型滑模控制的并聯(lián)混合有源電力濾波器［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2009，29（27）：29-35.

［7］魏學(xué)良.三相三線并聯(lián)型APF電流環(huán)數(shù)字化控制研究［D］.武漢：華中科技大學(xué)，2007.

［8］于晶榮，粟梅，孫堯.有源電力濾波器的改進(jìn)重復(fù)控制及其優(yōu)化設(shè)計(jì)［J］.電工技術(shù)學(xué)報(bào)，2012，27（2）：235-242.

［9］Asiminoaei L，Aeloiza E，Enjeti P N，et al.Shunt Active-pow?er-filter Topology Based on Parallel Interleaved Inverters［J］.IEEE Transactions on Industrial Electronics，2008，55（3）：1175-1189.

［10］Xie C，Wang Y，Zhong X，et al.A Novel Active Damping Method for LCL-filter-based Shunt Active Power Filter［C］∕∕IEEE In?ternational Symposium on Industrial Elcetronics，2012：64-69.

［11］宮金武，查曉明，陳佰鋒.一種快速重復(fù)控制策略在APF中的實(shí)現(xiàn)和分析［J］.電工技術(shù)學(xué)報(bào)，2011，26（10）：110-117.

［12］許明夏，林平，張濤，等.有源電力濾波器重復(fù)控制方法的設(shè)計(jì)［J］.電源學(xué)報(bào)，2012（2）：16-20.

［13］劉飛，宮金武，彭光強(qiáng)，等.一種復(fù)合式重復(fù)控制在并聯(lián)型有源電力濾波器設(shè)計(jì)中的應(yīng)用［J］.電工技術(shù)學(xué)報(bào)，2012，27（12）：138-145.

［14］Zou Z，Wang Z，Cheng M，et al.Active Power Filter for Harmo?nie Compensation Using a Digital Dual-mode-structure Repeti?tive Control Approach［C］∕∕IEEE International Symposium on Power Electronics for Distributed Generation Systems,2012：161-166.

Compound Current Control Method for Shunt Active Power Filters

TIAN Feiyan1，GAO Yunguang1，2，SONG Jiancheng1，ZHENG Lijun1，YANG Jiankang3，KANG Qingen3
（1.Shanxi Key Laboratory of Coal Mine Electrical Equipment and Intelligent Control，Taiyuan University of Technology，Taiyuan 030024，Shanxi,China；2.School of Electronics and Information Engineering，Taiyuan University of Science and Technology，Taiyuan 030024，Shanxi,China；3.Shanxi Jincheng Anthracite Coal Mining Group Co.，Ltd.，Jincheng 048000，Shanxi,China）

Aiming at the influence of nonlinear load on the quality of power grid，a shunt active power filter（SAPF）topology was researched.According to this topology，the mathematical model of the SAPF was established in the three-phase static coordinate system and the two-phase d-q synchronous rotating coordinate system.The transfer function of the model was deduced，and the current controller based on PI control was designed.By analyzing the limitation of PI controllers in current tracking control for active power filters，a compound control strategy combining repetitive control with PI control was proposed to improve steady state performance of the system.Experimental results demonstrate that the current control method combining PI control method with repetitive control method has a high steady-state compensation accuracy.It is able to effectively suppress harmonic currents caused by nonlinear loads and greatly improve grid current waveform.

harmonics suppression；active power filter；PI control；repetitive control；compound control

TM72

A

10.19457∕j.1001-2095.20170713

2016-06-15

修改稿日期：2016-08-20

晉煤集團(tuán)科技攻關(guān)項(xiàng)目（JSYJ-JSKF-2014-0050）；山西省重點(diǎn)實(shí)驗(yàn)室開(kāi)放課題（MEI-201501）

田飛燕（1991-），女，碩士，Email：feiy0218@163.com