• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      永磁同步電機(jī)的改進(jìn)無(wú)差拍預(yù)測(cè)抗擾前饋控制

      2022-07-19 02:57:02靳東松劉凌
      關(guān)鍵詞:無(wú)差觀測(cè)器滑模

      靳東松,劉凌

      (1.西安交通大學(xué)電力設(shè)備電氣絕緣國(guó)家重點(diǎn)實(shí)驗(yàn)室,710049,西安;2.西安交通大學(xué)電氣工程學(xué)院,710049,西安)

      由于永磁同步電機(jī)(PMSM)數(shù)字控制系統(tǒng)具有高精度、高效率等卓越的控制性能,因此該控制系統(tǒng)廣泛應(yīng)用于數(shù)控系統(tǒng)、電動(dòng)汽車驅(qū)動(dòng)與機(jī)器人等工業(yè)領(lǐng)域[1-3]。在實(shí)際的PMSM控制系統(tǒng)中,電流控制器控制的電流環(huán)時(shí)間常數(shù)短、響應(yīng)速度快,并且輸出的電流將直接作用于電機(jī)來(lái)產(chǎn)生電磁轉(zhuǎn)矩,是影響控制系統(tǒng)性能的關(guān)鍵因素。為實(shí)現(xiàn)PMSM控制系統(tǒng)的快速響應(yīng)及高穩(wěn)態(tài)精度,國(guó)內(nèi)外關(guān)于電流控制的研究主要有預(yù)測(cè)控制[4-13],滯環(huán)控制[14]以及滑模變結(jié)構(gòu)控制[15-21]等。

      預(yù)測(cè)電流控制方法具有高動(dòng)態(tài)性能、低電流諧波分量等優(yōu)點(diǎn),其需要精確的模型參數(shù)來(lái)保證其高控制性能。但是電機(jī)本體在運(yùn)行時(shí)是一個(gè)非線性時(shí)變系統(tǒng),其定子電阻與電感等電參數(shù)隨溫升變化而發(fā)生改變。電機(jī)電參數(shù)的變化將引起電磁轉(zhuǎn)矩脈動(dòng)和轉(zhuǎn)速波動(dòng),限制電流控制性能。所以有效地解決模型參數(shù)不匹配問題是在預(yù)測(cè)控制中實(shí)現(xiàn)高精度電流控制的關(guān)鍵。文獻(xiàn)[10]中在無(wú)差拍控制中引入魯棒電流預(yù)測(cè)算法,提高了系統(tǒng)電流環(huán)的動(dòng)態(tài)性能和穩(wěn)態(tài)精度,但是其中只針對(duì)電感參數(shù)失配。電機(jī)運(yùn)行時(shí)并不只有電感參數(shù)失配,實(shí)際使用時(shí)抗擾能力會(huì)大大受限。文獻(xiàn)[11]設(shè)計(jì)了兩種并行的擾動(dòng)觀測(cè)器來(lái)分別補(bǔ)償模型參數(shù)不匹配和逆變器非線性對(duì)電流控制精度的影響,提高了魯棒性并抵制了電流諧波分量,但是兩種觀測(cè)器分別設(shè)置在兩相坐標(biāo)軸上,總共就是四個(gè)觀測(cè)器,數(shù)據(jù)量與計(jì)算量大,對(duì)控制器要求較高,實(shí)際應(yīng)用將受限。文獻(xiàn)[12]設(shè)計(jì)了定子電流和擾動(dòng)觀測(cè)器來(lái)同時(shí)補(bǔ)償模型參數(shù)失配和一步時(shí)滯,提高了系統(tǒng)的抗擾性能,但是該方法未補(bǔ)償負(fù)載變化對(duì)轉(zhuǎn)速的影響。文獻(xiàn)[13]建立了一種考慮參數(shù)不確定性與外部干擾的PMSM模型,設(shè)計(jì)了一個(gè)高階滑模觀測(cè)器,分別用于速度和電流的前饋補(bǔ)償,提高了速度魯棒性和電流跟蹤精度。

      滯環(huán)控制的優(yōu)勢(shì)在于控制算法與硬件電路較為簡(jiǎn)單,且電流響應(yīng)相對(duì)平穩(wěn)。但是滯環(huán)控制開關(guān)頻率的變化較大,這導(dǎo)致控制器輸出電流中會(huì)有較多的高次諧波,會(huì)使電流畸變甚至失真,這一不足使其難以用于PMSM的高精度控制。

      滑??刂聘鶕?jù)系統(tǒng)狀態(tài)設(shè)計(jì)滑模面,通過切換控制量,迫使系統(tǒng)狀態(tài)沿設(shè)計(jì)的軌跡做高頻率、小幅度的滑動(dòng),增強(qiáng)系統(tǒng)狀態(tài)的抗擾性,該控制響應(yīng)迅速、易于實(shí)現(xiàn),適用于非線性時(shí)變系統(tǒng),如PMSM的速度控制[18]、負(fù)載前饋補(bǔ)償[19-20]及無(wú)位置傳感器控制[21]。

      PMSM驅(qū)動(dòng)系統(tǒng)中無(wú)差拍預(yù)測(cè)電流控制方法依賴于電機(jī)模型參數(shù),而滑??刂茖?duì)PMSM系統(tǒng)參數(shù)不匹配和干擾具有強(qiáng)魯棒性,所以可研究將兩種算法結(jié)合。文獻(xiàn)[22]將滑??刂婆c無(wú)差拍預(yù)測(cè)電流控制結(jié)合,分別將改進(jìn)算法用于電流環(huán)與轉(zhuǎn)速環(huán),提高了系統(tǒng)動(dòng)態(tài)性能,但是計(jì)算量大,對(duì)控制器的要求較高;文獻(xiàn)[23]將積分滑??刂婆c無(wú)差拍預(yù)測(cè)電流控制結(jié)合,消除了傳統(tǒng)滑模控制的到達(dá)階段,增強(qiáng)了系統(tǒng)的動(dòng)態(tài)性能與穩(wěn)定性,超螺旋算法的應(yīng)用抑制了系統(tǒng)穩(wěn)定階段抖振,但是增大了啟動(dòng)階段的抖振。針對(duì)上述工作中未考慮在轉(zhuǎn)速環(huán)中負(fù)載轉(zhuǎn)矩及干擾補(bǔ)償對(duì)無(wú)差拍控制模型失配的影響,本文將進(jìn)一步研究無(wú)差拍控制的滑模擾動(dòng)補(bǔ)償。

      本文提出了面向永磁同步電機(jī)無(wú)差拍控制的滑模擾動(dòng)前饋補(bǔ)償策略,設(shè)計(jì)了滑模負(fù)載觀測(cè)器(SMLTO)并前饋補(bǔ)償?shù)诫娏骺刂破?采用傳統(tǒng)指數(shù)趨近律設(shè)計(jì)了一種自適應(yīng)參數(shù)的改進(jìn)指數(shù)趨近律并用于改進(jìn)滑模負(fù)載觀測(cè)器(MSMLTO)。通過觀測(cè)器對(duì)預(yù)測(cè)電流控制器的前饋補(bǔ)償,改善了電流控制器對(duì)模型參數(shù)不匹配及負(fù)載變化的抗擾能力,并完成了對(duì)穩(wěn)定性的證明。最后,數(shù)值仿真結(jié)果驗(yàn)證了該控制算法的可行性和有效性。

      1 預(yù)測(cè)電流控制

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

      假定磁路不飽和,不計(jì)磁滯與渦流損耗影響,空間磁場(chǎng)呈正弦分布,表貼式永磁同步電機(jī)(SPMSM)在同步旋轉(zhuǎn)坐標(biāo)系下的定子電壓方程為[24]

      (1)

      轉(zhuǎn)矩方程為

      Te=1.5pnψfiq

      (2)

      運(yùn)動(dòng)方程為

      (3)

      式中:ud、uq、id、iq分別為電機(jī)定子電壓、定子電流在d、q軸上的分量;Ld=Lq=Ls為電機(jī)的直軸與交軸同步電感;pn是電機(jī)極對(duì)數(shù);ωe為電角速度;p是微分算子;Rd=Rq=Rs是電機(jī)定子d、q軸電阻;ψf為電機(jī)轉(zhuǎn)子磁通;Te是電機(jī)電磁轉(zhuǎn)矩;Tl是電機(jī)負(fù)載轉(zhuǎn)矩;B是黏滯系數(shù);J是轉(zhuǎn)動(dòng)慣量。

      1.2 無(wú)差拍預(yù)測(cè)電流控制

      當(dāng)控制系統(tǒng)的采樣周期TSC足夠短時(shí),PMSM的離散時(shí)間模型可以用一階泰勒級(jí)數(shù)展開描述,則PMSM的離散電流預(yù)測(cè)模型如下式所示[25]

      (4)

      式中

      其中下標(biāo)k表示離散量。

      傳統(tǒng)的無(wú)差拍預(yù)測(cè)電流控制中,根據(jù)離散預(yù)測(cè)模型(4),使實(shí)際電流矢量經(jīng)過一個(gè)調(diào)制周期后達(dá)到參考電流的電機(jī)定子電壓如下式所示

      (5)

      電機(jī)運(yùn)行時(shí)電機(jī)參數(shù)會(huì)變化,不妨記Rs、Ls、ψf為電機(jī)實(shí)時(shí)參數(shù),Rs0、Ls0、ψf0為電機(jī)初始參數(shù)。由式(4)分析可知,電機(jī)電參數(shù)的變化將會(huì)影響ud(k)與uq(k),進(jìn)而影響電機(jī)的電磁轉(zhuǎn)矩與轉(zhuǎn)速,造成系統(tǒng)不穩(wěn)定。

      1.3 負(fù)載轉(zhuǎn)矩變化分析

      采用id=0的矢量控制策略,則由式(2)和式(3)可推導(dǎo)得

      (6)

      假設(shè)電機(jī)負(fù)載轉(zhuǎn)矩變化量為ΔTl>0,則

      (7)

      KD[eω(k)-eω(k-1)]

      (8)

      2 積分滑模負(fù)載轉(zhuǎn)矩觀測(cè)器

      由式(4)與式(5)分析,Rs與Ls對(duì)電流控制器輸出量的作用范圍較大,影響較大,Rs與Ls的變化將影響i(k+1)的值,影響電流控制器輸出電壓矢量值,進(jìn)而造成電機(jī)電磁轉(zhuǎn)矩與轉(zhuǎn)速的畸變。本文設(shè)計(jì)一種無(wú)需使用Rs與Ls兩電參數(shù)的負(fù)載轉(zhuǎn)矩觀測(cè)器,將觀測(cè)值在轉(zhuǎn)換后前饋補(bǔ)償至iq(k+1)以補(bǔ)償Rs與Ls變化的影響,采用id=0的矢量控制策略,id接近于0,影響較小,故不對(duì)id做出補(bǔ)償。

      2.1 SMLTO的設(shè)計(jì)

      PMSM的數(shù)字控制系統(tǒng)中,數(shù)字控制器的開關(guān)頻率一般較高,且電機(jī)系統(tǒng)的機(jī)械周期往往遠(yuǎn)大于其電氣周期,所以在數(shù)字控制器的一個(gè)采樣周期中,可近似地認(rèn)為負(fù)載轉(zhuǎn)矩沒有變化,即dTl/dt=0。針對(duì)研究的對(duì)象及所采取的控制策略,在公式的基礎(chǔ)上,建立關(guān)于電機(jī)轉(zhuǎn)速和負(fù)載轉(zhuǎn)矩的狀態(tài)方程

      (9)

      將Tl和ω作為觀測(cè)值,建立SMLTO的狀態(tài)方程如下

      (10)

      式(9)減去式(10),可得

      (11)

      選擇積分滑模面作為該觀測(cè)器的滑模面

      (12)

      式中α,β>0,均為滑模面系數(shù)。

      2.2 新型變系數(shù)指數(shù)趨近律設(shè)計(jì)

      傳統(tǒng)的指數(shù)趨近律如下所示

      (13)

      式中k1與λ為趨近律系數(shù)。

      對(duì)式(13)從0至t積分,可得其趨近時(shí)間

      (14)

      式中s0為滑動(dòng)模態(tài)的初始狀態(tài)值。

      系統(tǒng)狀態(tài)s偏離滑模面較遠(yuǎn)時(shí),趨近速率為k1sgn(s)與λs的和,此時(shí)趨近速率較大;系統(tǒng)狀態(tài)s離滑模面較近時(shí),k1sgn(s)在起主要作用,趨近速率基本由k1決定。若想獲得較小的趨近時(shí)間,就需要較大的k1值,但同時(shí)又帶來(lái)較大抖振;反之減小抖振又會(huì)增大趨近時(shí)間,所以需在兩者間取舍。

      新型變系數(shù)指數(shù)趨近律如下所示

      (15)

      (16)

      式中:K、λ、k2、δ、ε、γ均為該趨近率的系數(shù),其中λ>0,k2>0,δ>0,0<γ<1,1<ε<2。

      式(15)中飽和函數(shù)sigmoid(s)的表達(dá)式如下

      (17)

      式中Δ為邊界層厚度。

      分析式(15)與式(16),當(dāng)系統(tǒng)狀態(tài)s離滑模面較遠(yuǎn)時(shí),K值可簡(jiǎn)化為

      (18)

      由式(18)可知,當(dāng)|s|較大時(shí),(2-ε)e-δ|s|趨于0,故此時(shí)K較大。

      若k1=k2,則改進(jìn)趨近律的趨近速率大于式(13)所示的傳統(tǒng)指數(shù)趨近律。

      當(dāng)系統(tǒng)狀態(tài)s離滑模面較近時(shí),K值簡(jiǎn)化為

      (19)

      式(19)的分母中,當(dāng)|e1|趨近于0時(shí),K值也隨之趨近于0,這時(shí)趨近速率逐漸減小,系統(tǒng)抖振也逐漸減小。

      假定s(t)=0,s0為系統(tǒng)初始狀態(tài),對(duì)式(15)從0到t積分,可得趨近時(shí)間

      t=(2-ε+1/|e1|)[|s0|e-δ|s0|+

      (20)

      選擇趨近律的參數(shù)λ與δ使其滿足不等式|s0|(δ2+λδ+1)<1-λ-δ,同時(shí)取k1=k2,又已知ε>0,0

      (21)

      則有t-t0<0,即按系統(tǒng)條件選擇參數(shù)時(shí)該改進(jìn)趨近律的趨近速率大于傳統(tǒng)的指數(shù)趨近律。

      若取t=t0,則由式(14)與式(20)可得

      (22)

      整理式(22)可得k2

      2.3 SMLTO的穩(wěn)定性分析

      整理式(11)、式(12)及式(15)可得控制函數(shù)

      (23)

      當(dāng)系統(tǒng)狀態(tài)s抵達(dá)滑模面時(shí),有

      (24)

      再將式(24)代入式(11)整理得

      (25)

      由式(25)可解得

      (26)

      式中c為常數(shù)。由穩(wěn)定性條件,為使e2在有限時(shí)間內(nèi)收斂至0,反饋增益m需滿足m<0。

      選定李雅普諾夫函數(shù)

      (27)

      求李雅普諾夫函數(shù)的導(dǎo)數(shù),可得

      s(-Ksigmoid(s)-λs)=-sKsigmoid(s)-λs2

      (28)

      由式(16)及各參數(shù)選擇條件可得

      (29)

      綜上所述,改進(jìn)負(fù)載轉(zhuǎn)矩觀測(cè)器是穩(wěn)定的。

      2.4 觀測(cè)負(fù)載轉(zhuǎn)矩的前饋補(bǔ)償

      (30)

      (31)

      補(bǔ)償后的電機(jī)離散電流預(yù)測(cè)方程如下

      (32)

      圖1 基于MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制框圖

      3 數(shù)值仿真分析

      PMSM參數(shù)如表1所示。為驗(yàn)證改進(jìn)算法的有效性,使用Matlab/Simulink進(jìn)行仿真。仿真中的電流采樣周期為10 μs,仿真結(jié)果將與傳統(tǒng)預(yù)測(cè)電流控制、傳統(tǒng)SMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制和MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制比較。因?yàn)殡姍C(jī)在運(yùn)行中會(huì)有散熱系統(tǒng)給轉(zhuǎn)子降溫,而且電機(jī)轉(zhuǎn)子磁鏈的測(cè)量往往是在電機(jī)正常運(yùn)轉(zhuǎn)時(shí)進(jìn)行,所以在額定參數(shù)內(nèi)運(yùn)行時(shí)電機(jī)轉(zhuǎn)子磁鏈的變化一般可忽略不計(jì)。

      表1 永磁同步電機(jī)參數(shù)

      對(duì)同一臺(tái)PMSM使用不同控制策略時(shí)的仿真波形圖如圖2至圖4所示。速度指令為3 000 r/min,空載啟動(dòng),在0.5 s加8 N·m負(fù)載轉(zhuǎn)矩,在0.8 s將負(fù)載轉(zhuǎn)矩減至3 N·m,運(yùn)行至1 s。初始參數(shù)下不同控制方法的仿真結(jié)果如圖2所示。分析圖2中各波形圖可知,此時(shí)在3種控制方法下電機(jī)均可正常運(yùn)行,但是傳統(tǒng)DPCC下加速過程相對(duì)前饋補(bǔ)償下的DPCC加速過程遲約0.09 s,且自圖(a)至圖(c)加載時(shí)觀測(cè)轉(zhuǎn)矩的超調(diào)量依次為18.75%、16.7%與5%,減載時(shí)為26.73%、26.67%與16.67%,抖振依次減小,可知電機(jī)在初始參數(shù)下運(yùn)行前饋補(bǔ)償能有效提高加速性能,改進(jìn)觀測(cè)器能更好地減小超調(diào)量和減輕抖振。

      (a)傳統(tǒng)無(wú)差拍預(yù)測(cè)電流控制

      (a)傳統(tǒng)無(wú)差拍預(yù)測(cè)電流控制

      (a)傳統(tǒng)無(wú)差拍預(yù)測(cè)電流控制

      圖3為L(zhǎng)=0.5L0時(shí),即電機(jī)定子電感衰減為初始值的0.5倍時(shí),不同控制方法下的仿真波形圖。由圖3(a)可知,未加負(fù)載時(shí)電機(jī)運(yùn)行基本正常,加上負(fù)載后電機(jī)運(yùn)行很不平穩(wěn),控制系統(tǒng)受影響較大;而有負(fù)載轉(zhuǎn)矩前饋補(bǔ)償?shù)那闆r下電機(jī)運(yùn)行較穩(wěn)定,將圖3(b)與圖3(c)分別與圖2(b)與圖2(c)對(duì)應(yīng)控制策略比較,可知兩種控制策略下電機(jī)抖振均有增大,但是其中傳統(tǒng)SMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制抖振增大較嚴(yán)重,加載時(shí)超調(diào)量約為13.75%,減載時(shí)超調(diào)量約為26.67%,而MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制則基本未受影響,加載時(shí)超調(diào)量約為8.75%,減載時(shí)超調(diào)量約為16.67%,對(duì)比圖3(a)與圖3(b)、圖3(c),傳統(tǒng)DPCC時(shí)系統(tǒng)反應(yīng)較采用轉(zhuǎn)矩補(bǔ)償策略下的反應(yīng)時(shí)間更長(zhǎng),分析觀測(cè)轉(zhuǎn)矩可知,針對(duì)初始加速指令,圖3(b)中響應(yīng)時(shí)間比圖3(a)短約0.1 s,圖3(c)又比圖3(b)短約0.02 s,而且負(fù)載變化時(shí)電磁轉(zhuǎn)矩恢復(fù)時(shí)間更短,可見所提的補(bǔ)償策略有效提高了系統(tǒng)的響應(yīng)速度。分析可知電感衰減時(shí)前饋補(bǔ)償可使系統(tǒng)正常運(yùn)行,動(dòng)態(tài)性能更好,而經(jīng)MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制則對(duì)電感衰減的抗擾性能更強(qiáng),系統(tǒng)運(yùn)行更穩(wěn)定。

      圖4為L(zhǎng)=0.5L0、R=5R0,即定子電感衰減為初始值的0.5倍、定子電阻增加為初始值的5倍時(shí),不同控制方法下的仿真波形圖。分析圖4可知,電機(jī)在空載時(shí)受影響較小,加負(fù)載后,傳統(tǒng)無(wú)差拍預(yù)測(cè)電流控制下觀測(cè)轉(zhuǎn)矩跟蹤失準(zhǔn),受影響較大,在實(shí)際的PMSM驅(qū)動(dòng)系統(tǒng)中電機(jī)難以運(yùn)行;而有負(fù)載轉(zhuǎn)矩前饋補(bǔ)償?shù)那闆r下,比較圖4(b)與圖4(c),圖4(b)中觀測(cè)轉(zhuǎn)矩抖振較大,加載時(shí)超調(diào)量約為15%,減載時(shí)超調(diào)量約為26.16%,且負(fù)載轉(zhuǎn)矩跟蹤效果較差,而圖4(c)中觀測(cè)轉(zhuǎn)矩基本未受影響,加載時(shí)超調(diào)量約為8.69%,減載時(shí)超調(diào)量約為15.49%,針對(duì)初始加速指令,圖4(b)中響應(yīng)時(shí)間比圖4(a)短約0.11 s,圖4(c)又比圖4(b)短約0.02 s,可知MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制抵抗參數(shù)變化的性能較傳統(tǒng)SMLTO有所改善。

      4 結(jié) 論

      本文根據(jù)永磁同步電機(jī)dq兩相旋轉(zhuǎn)坐標(biāo)系下的定子電壓方程、轉(zhuǎn)矩方程與運(yùn)動(dòng)方程等動(dòng)態(tài)數(shù)學(xué)模型,提出了一種永磁同步電機(jī)的改進(jìn)無(wú)差拍預(yù)測(cè)抗擾前饋控制策略。采用改進(jìn)指數(shù)趨近律與滑模負(fù)載轉(zhuǎn)矩觀測(cè)器結(jié)合,將觀測(cè)轉(zhuǎn)矩轉(zhuǎn)換為電流量后對(duì)傳統(tǒng)無(wú)差拍預(yù)測(cè)電流控制進(jìn)行前饋補(bǔ)償以提高控制器對(duì)模型參數(shù)不匹配及負(fù)載轉(zhuǎn)矩變化的抗擾能力。將改進(jìn)的預(yù)測(cè)控制方法用于永磁同步電機(jī)控制系統(tǒng)進(jìn)行仿真驗(yàn)證,結(jié)果表明,與傳統(tǒng)的無(wú)差拍預(yù)測(cè)電流控制相比,傳統(tǒng)SMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制有更快的加速性能和較強(qiáng)的抵抗模型參數(shù)不匹配和負(fù)載變化的能力,MSMLTO補(bǔ)償?shù)念A(yù)測(cè)電流控制較傳統(tǒng)SMLTO有更好的抵抗模型參數(shù)不匹配和負(fù)載變化的能力,控制性能進(jìn)一步改善。

      猜你喜歡
      無(wú)差觀測(cè)器滑模
      基于無(wú)差拍電流預(yù)測(cè)控制的PMSM電感失配研究
      雙BP神經(jīng)網(wǎng)絡(luò)的永磁同步電機(jī)磁鏈與轉(zhuǎn)矩?zé)o差拍控制
      基于組合滑??刂频慕^對(duì)重力儀兩級(jí)主動(dòng)減振設(shè)計(jì)
      測(cè)控技術(shù)(2018年4期)2018-11-25 09:47:26
      并網(wǎng)逆變器逆系統(tǒng)自學(xué)習(xí)滑??箶_控制
      基于觀測(cè)器的列車網(wǎng)絡(luò)控制
      基于非線性未知輸入觀測(cè)器的航天器故障診斷
      基于電感電流無(wú)差拍控制的三相并網(wǎng)逆變器
      基于干擾觀測(cè)器的PI控制單相逆變器
      基于電流預(yù)測(cè)的有源濾波器無(wú)差拍控制
      北海市| 类乌齐县| 清镇市| 乐至县| 乌兰浩特市| 项城市| 沂水县| 红桥区| 松桃| 伊川县| 广河县| 刚察县| 夏邑县| 通渭县| 章丘市| 康保县| 通山县| 宣威市| 汶川县| 运城市| 横峰县| 虎林市| 贵港市| 赤水市| 竹北市| 龙门县| 东山县| 桂东县| 吕梁市| 山丹县| 宁武县| 铜陵市| 商水县| 自贡市| 榕江县| 汪清县| 徐州市| 会东县| 耒阳市| 东辽县| 汾西县|