国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

基于多種群合作協(xié)同算法的電力系統(tǒng)調(diào)度優(yōu)化

2016-03-24 06:40:32伍大清楊治平
關(guān)鍵詞:適應(yīng)度種群粒子

劉 勇, 伍大清, 2, 3, , 5, 孫 莉, 楊治平

(1. 四川省人工智能重點實驗室, 四川 自貢 643000; 2. 南華大學(xué) 計算機(jī)科學(xué)與技術(shù)學(xué)院, 湖南 衡陽 421001;

3. 安徽大學(xué) 教育部計算智能與信號處理重點實驗室, 安徽 合肥 230039; 4. 東華大學(xué) 旭日工商管理學(xué)院,

上海 200051; 5. 農(nóng)業(yè)部農(nóng)業(yè)信息服務(wù)技術(shù)重點實驗室, 北京 100081)

基于多種群合作協(xié)同算法的電力系統(tǒng)調(diào)度優(yōu)化

劉 勇1, 伍大清1, 2, 3, 4, 5, 孫 莉4, 楊治平4

(1. 四川省人工智能重點實驗室, 四川 自貢 643000; 2. 南華大學(xué) 計算機(jī)科學(xué)與技術(shù)學(xué)院, 湖南 衡陽 421001;

3. 安徽大學(xué) 教育部計算智能與信號處理重點實驗室, 安徽 合肥 230039; 4. 東華大學(xué) 旭日工商管理學(xué)院,

上海 200051; 5. 農(nóng)業(yè)部農(nóng)業(yè)信息服務(wù)技術(shù)重點實驗室, 北京 100081)

為解決電力系統(tǒng)環(huán)境經(jīng)濟(jì)調(diào)度這一復(fù)雜多目標(biāo)約束優(yōu)化問題,提出一種有效多種群合作協(xié)同優(yōu)化算法,采用多個種群對搜索空間進(jìn)行搜索,運(yùn)用新型的速度位移更新方式以及種群周期內(nèi)拓?fù)浣Y(jié)構(gòu)重組策略. 結(jié)果表明,該算法能對解空間進(jìn)行更加全面、充分的探索,可快速找到一組分布具有盡可能好的逼近性、寬廣性和均勻性的最優(yōu)解集合.將該算法應(yīng)用到某電力系統(tǒng)的環(huán)境經(jīng)濟(jì)調(diào)度中, 其仿真計算結(jié)果與其他求解方法結(jié)果的對比分析表明, 該算法可以有效兼顧全局收斂性和Pareto非劣調(diào)度方案的多樣性, 具有較高的效率以及魯棒性.

合作協(xié)同優(yōu)化; 多種群優(yōu)化; 粒子群優(yōu)化; 環(huán)境經(jīng)濟(jì)調(diào)度

近年來,隨著人類對環(huán)境保護(hù)的日益重視,電力系統(tǒng)中傳統(tǒng)的經(jīng)濟(jì)調(diào)度模式已不能滿足環(huán)保的需求,一種綜合考慮發(fā)電成本和污染氣體排放的新型環(huán)境經(jīng)濟(jì)調(diào)度(environmental-economic dispatch, EED)模式受到國內(nèi)外學(xué)者的廣泛關(guān)注.國內(nèi)外學(xué)者通過分析運(yùn)用電力調(diào)度的約束模型,將電力調(diào)度的經(jīng)濟(jì)最小化作為單一目標(biāo)進(jìn)行優(yōu)化,也有采用多目標(biāo)遺傳算法NSGA-Ⅱ[1]、粒子群優(yōu)化算法[2]、進(jìn)化規(guī)劃算法[3]和差分進(jìn)化算法[4]等來求解電力系統(tǒng)環(huán)境經(jīng)濟(jì)調(diào)度問題,并取得了豐富的研究成果,但這些方法仍存在諸如計算速度慢、搜索精度不高、難以處理實際工程問題中復(fù)雜的約束條件等問題.

本文以求解EED的多目標(biāo)約束優(yōu)化問題為背景,將多目標(biāo)分解思想與人類社交行為有機(jī)地結(jié)合起來,提出一種有效的多種群合作協(xié)同微粒群優(yōu)化算法(an effective co-evolutionary multi-swarm particle swarm optimization,簡稱ECMPSO).ECMPSO將多目標(biāo)優(yōu)化問題分解為若干個單目標(biāo)優(yōu)化子問題,每個子問題對應(yīng)一個子種群,不同子種群間通過共享外部存檔集相互合作,有效地近似整個Pareto前沿面,提出了基于人類獎懲行為的獎勵-懲罰學(xué)習(xí)因子的微粒群優(yōu)化算法速度更新方式.該算法的理論分析與試驗結(jié)果表明,ECMPSO在收斂性和多樣性方面均表現(xiàn)出明顯優(yōu)勢,能夠避免陷入局部最優(yōu)解,可高效搜索外部存檔中新粒子,能夠獲得均勻的Pareto前沿.

1 多目標(biāo)環(huán)境經(jīng)濟(jì)調(diào)度問題的數(shù)學(xué)描述

電力系統(tǒng)環(huán)境經(jīng)濟(jì)調(diào)度問題是在已知機(jī)組開停計劃的情況下,以某時間段機(jī)組的有功出力為決策變量,將系統(tǒng)總負(fù)荷合理分配,使整個機(jī)組消耗的燃料以及污染氣體排放量達(dá)到最少化[5].

1.1 約束條件

(1) 發(fā)電機(jī)運(yùn)行容量約束.每個機(jī)組的發(fā)電功率應(yīng)介于其最小輸出功率與最大輸出功率之間,即

(1)

(2) 有功功率平衡約束.負(fù)載總的需求功率與網(wǎng)絡(luò)損耗之和等于各機(jī)組發(fā)電功率之和,即

(2)

式中:Ploss為系統(tǒng)網(wǎng)絡(luò)損耗;PD為系統(tǒng)總負(fù)荷需求;NG為系統(tǒng)內(nèi)發(fā)電機(jī)總數(shù).

采用B系數(shù)法,系統(tǒng)網(wǎng)損與發(fā)電機(jī)有功功率的關(guān)系為

(3)

式中:Bij, Bi0, B00稱為B系數(shù).

1.2 目標(biāo)函數(shù)

(1) 燃料消耗目標(biāo)函數(shù)為

(4)

式中:ai, bi, ci均為系統(tǒng)參數(shù).

(2) 污染氣體排放量函數(shù)為

式中:αi, βi, ξi, γi均為系統(tǒng)參數(shù).

2 基于多種群合作協(xié)同的環(huán)境經(jīng)濟(jì)調(diào)度算法

2.1 算法原理

基本粒子群優(yōu)化(PSO)算法[6]通常用來求解單目標(biāo)優(yōu)化問題,每個粒子視為搜索空間中待優(yōu)化問題的可行解,假設(shè)n個粒子在D維的空間中,粒子之間通過相互競爭與協(xié)作來尋找問題的最優(yōu)解,在每次迭代中,粒子個體搜索到的最好位置稱為個體最優(yōu),用pbest表示,群體中所有粒子搜索到最好位置稱為全局最優(yōu),用gbest表示,通常以尋找最低適應(yīng)度值為最優(yōu).

在基本的PSO算法中,所有的粒子僅僅圍繞個體最優(yōu)和群體最優(yōu)學(xué)習(xí),顯然,這是一個理想社會條件.然而,本文算法引進(jìn)人的社交行為特點,而社交過程中環(huán)境因素給人帶來巨大影響.粒子在搜索過程中也同樣如此,不僅受到周圍好環(huán)境(適應(yīng)度值低的鄰居粒子)的影響,也可能受到差環(huán)境(適應(yīng)度值高的鄰居粒子)的影響.但如果客觀、理性看待這些“壞習(xí)慣/行為”,提前對這些“壞習(xí)慣/行為”的粒子有所警惕,采取措施,那么對搜索是非常有益的.因此,ECMPSO算法中引進(jìn)最大局部鄰居粒子Nmax,即每次迭代,每個小種群都會產(chǎn)生一個局部最優(yōu)粒子,所對應(yīng)的適應(yīng)度值排序,最大的適應(yīng)度粒子被選作Nmax, Nmax(t)=argmax{f(pbest1), f(pbest2), …, f(pbestn)}, f(·)代表每個粒子的適應(yīng)度函數(shù)的評估值.ECMPSO引進(jìn)一個隨機(jī)學(xué)習(xí)因子r3,且r3服從標(biāo)準(zhǔn)正態(tài)N(0, 1)分布的隨機(jī)數(shù).每一代新個體的生成,判斷是否對共享的外部存檔集做出貢獻(xiàn)決定r3的取值.如果產(chǎn)生的新個體能支配外部存檔集中的非支配解,則產(chǎn)生r3>0的獎勵學(xué)習(xí)因子,加快粒子的飛行速度,則有利于提高粒子勘測區(qū)域深度;如果產(chǎn)生的新個體未能支配外部存檔集中的非支配解,則產(chǎn)生r3<0懲罰學(xué)習(xí)因子,將減緩粒子的飛行速度,有利于提高粒子開發(fā)廣度.

基本PSO算法中,粒子在個體最優(yōu)以及全局最優(yōu)的引導(dǎo)下尋找區(qū)域的最優(yōu)值,搜索過程如圖1(a)所示,該搜索方式能促使粒子搜索速度加快,但很容易陷入局部最優(yōu)解.ECMPSO算法的粒子搜索過程如圖1(b)和1(c)所示,從圖1(b)和1(c)能清楚看到,獎勵學(xué)習(xí)項以及懲罰學(xué)習(xí)項能促使粒子探索更廣闊的空間,有助于粒子跳出局部最優(yōu)解,從而提高種群多樣性,使其快速收斂到全局最優(yōu)值.

ECMPSO算法中的子種群采用環(huán)形拓?fù)浣Y(jié)構(gòu),為了粒子向周圍環(huán)境中的不同的粒子學(xué)習(xí),充分保持

(a) 基本PSO搜索過程

(b) ECMPSO獎勵學(xué)習(xí)過程

(c) ECMPSO懲罰學(xué)習(xí)過程

Fig.1 Process of particles searching

2.2 數(shù)學(xué)描述

為了減少參數(shù)靈敏性對算法帶來的影響,本文算法使用隨機(jī)學(xué)習(xí)因子r1以及r2代替基本粒子群中的加速因子,每個粒子通過式(7)和(8)更新速度和位移.

(6)

(7)

2.3 算法描述

ECMPSO算法基于多目標(biāo)技術(shù),種群中的粒子相互協(xié)同合作尋找到分布盡可能好的逼近性、寬廣性和均勻性Pareto最優(yōu)解集合.基于ECMPSO算法的環(huán)境經(jīng)濟(jì)發(fā)電調(diào)度流程如下:

步驟1 獲得具有多臺發(fā)電機(jī)組的電力系統(tǒng)中每臺機(jī)組的出力數(shù)據(jù)上限與下限、燃料消耗函數(shù)的系數(shù)數(shù)據(jù)、有害氣體排放量函數(shù)的系數(shù)數(shù)據(jù)、輸電線路損耗的B系數(shù)數(shù)據(jù)和系統(tǒng)總負(fù)荷數(shù)據(jù).

步驟2 建立電力系統(tǒng)環(huán)境經(jīng)濟(jì)調(diào)度問題的數(shù)學(xué)優(yōu)化模型.

工作經(jīng)驗告訴我,面對幸福的選擇困難情況時,相信自己的直覺,結(jié)果往往美滋滋。這一次,我毫不猶豫地選擇GT,法拉利的GT跑車。

步驟3 初始化種群.隨機(jī)初始化粒子數(shù)為N的群體以及對應(yīng)的適應(yīng)度, 針對多目標(biāo)問題的目標(biāo)個數(shù)N,將種群劃分成N個子種群,各個粒子所對應(yīng)的速度,最大迭代次數(shù),粒子的局部最優(yōu)位置pbest,群體最佳位置gbest以及每個子群體中Nmax,初始化外部存檔集.

步驟4 迭代更新.對于每個種群的每個粒子進(jìn)行如下操作:(1)對種群中所有的粒子采用式(6)和(7)更新粒子的速度和位置;(2)計算新粒子的各個目標(biāo)適應(yīng)度值,判斷是否能支配外部存檔中的非支配解,若是則產(chǎn)生獎勵學(xué)習(xí)因子,若否則產(chǎn)生懲罰學(xué)習(xí)因子,分別更新pbest、gbest以及Nmax.注意, 這里對速度越界的粒子適應(yīng)度值不進(jìn)行更新.

步驟5 更新外部存檔集,判斷外部存檔大小是否超過種群規(guī)模,若是,使用擁擠距離更新外部存檔.

步驟6 迭代計數(shù)器累增1,是否能整除拓?fù)浣Y(jié)構(gòu)重組周期L.若是則執(zhí)行步驟4;若否則執(zhí)行步驟7.

步驟7 迭代次數(shù)累增,判斷是否滿足算法終止條件.若滿足,則執(zhí)行步驟8,否則轉(zhuǎn)步驟4.

步驟8 輸出Pareto最優(yōu)前沿面, 算法結(jié)束.

步驟9 將確定的最終解作為指令,通過自動發(fā)電控制裝置發(fā)送給相關(guān)發(fā)電廠或機(jī)組的自動控制調(diào)節(jié)裝置,實現(xiàn)對機(jī)組發(fā)電功率的控制.

3 仿真試驗結(jié)果與分析

3.1 測試函數(shù)

試驗部分選用了國際測試函數(shù)ZDT系列的4個雙目標(biāo)問題來驗證ECMPSO算法的有效性,這些測試問題能夠在不同的方面對進(jìn)化多目標(biāo)問題進(jìn)行測試.試驗結(jié)果與NSGAII[1], MOABC[8], MOCLPSO[9]進(jìn)行了對比.其中,NSGAII是一種公認(rèn)而有效地進(jìn)化多目標(biāo)優(yōu)化算法,MOABC是最近提出來的蜂群多目標(biāo)優(yōu)化方法, MOCLPSO是一種非常有效的微粒群多目標(biāo)優(yōu)化算法.表1列出了4個測試函數(shù)的表達(dá)式.

表1 測試函數(shù)

綜合評價指標(biāo)IGD(inverted generational distance)[10]被用來評估算法的性能.假定Q*為多目標(biāo)優(yōu)化問題理想PF上的一組均勻采樣,Q為多目標(biāo)優(yōu)化算法求得的一組理想PF的逼近解,則解集Q的IGD指標(biāo)定義如下

(8)

表2給出了4種算法在不同測試函數(shù)上的IGD度量結(jié)果,其中的數(shù)據(jù)為算法運(yùn)行30次得到的平均值和方差.由表2可知,ECMPSO算法在ZDT1、 ZDT2、 ZDT3、 ZDT6收斂性和分布性上均優(yōu)于其他3種算法,在ZDT3函數(shù)中MOCLPSO算法獲得的IGD值與ECMPSO算法差不多,兩者無顯著差異.

表2 4種算法分別求得4個測試函數(shù)IGD評價指標(biāo)

試驗數(shù)據(jù)表明ECMPSO算法具有很強(qiáng)穩(wěn)定性.圖2中分別給出4種算法在ZDT1、 ZDT2測試函數(shù)上的Pareto前沿面.由圖2表明,ECMPSO算法使用獎懲學(xué)習(xí)因子,種群周期內(nèi)拓?fù)浣Y(jié)構(gòu)重組策略在解群分布的均勻性和寬廣性方面要明顯優(yōu)于其他3種算法.

3.2 環(huán)境經(jīng)濟(jì)調(diào)度優(yōu)化仿真試驗

本試驗用標(biāo)準(zhǔn)的IEEE 30節(jié)點六發(fā)電機(jī)測試系統(tǒng),該系統(tǒng)的負(fù)荷需求為2.834 p.u,所有的線路通過41條輸電線路進(jìn)行連接.詳細(xì)的燃料成本數(shù)據(jù)和氮氧化合物的排放系數(shù)見文獻(xiàn)[11].為了清晰表明不同算法的性能,測試時考慮有網(wǎng)損和無網(wǎng)損兩種

圖2 各算法在2 000次迭代后的最優(yōu)Pareto前沿面

情況.兩種情況具體如下:

情況1 忽略6個機(jī)組的系統(tǒng)損耗,只考慮容量約束和平衡約束;

情況2 考慮系統(tǒng)損耗,同時考慮容量約束和平衡約束.

為驗證本文算法的可靠性,在參數(shù)設(shè)置一致的情況下,每組試驗獨立運(yùn)行20次,分別以最低的發(fā)電成本以及相應(yīng)較好的排放、最低的排放量以及相應(yīng)較低的發(fā)電成本兩種方式來評估ECMPSO算法、DE-IMOPSO算法[12]、NSGAII算法[1]以及2LB-MOPSO算法.4種算法在情況1和2條件下的運(yùn)算結(jié)果如表3~6所示.

表3 不同算法在情況1中燃料花費(fèi)最小的運(yùn)算結(jié)果

表4 不同算法在情況1中氣體排放量最少的運(yùn)算結(jié)果

表5 不同算法在情況2中燃料花費(fèi)最小的運(yùn)算結(jié)果

表6 不同算法在情況2中氣體排放量最少的運(yùn)算結(jié)果

由表3和4可知,在情況1中ECMPSO算法求得的燃料花費(fèi)最小值為600.004 $/h和氣體排放量最小值為0.194 02 t/h,優(yōu)于DE-IMOPSO、 NSGA-Ⅱ、 2LB-MOPSO 3種算法的結(jié)果.由表5和6可知,在情況2中ECMPSO算法求得的燃料花費(fèi)最小值為605.810 $/h和氣體排放量最小值為0.189 72 t/h,均優(yōu)于DE-IMOPSO、 NSGA-Ⅱ、 2LB-MOPSO算法的結(jié)果.表3~6結(jié)果表明,ECMPSO算法取得了更好的最優(yōu)解.ECMPSO算法在情況1和2中求得Pareto最優(yōu)前沿如圖3所示.由圖3可以觀察到,ECMPSO算法所得解遍布整個均衡面,在保持多樣性的同時均勻分布性也較好.

(a) 情況1

(b) 情況2

4 結(jié) 語

針對電力系統(tǒng)的新型環(huán)境經(jīng)濟(jì)調(diào)度模式,提出了一種多種群協(xié)同微粒群優(yōu)化(ECMPSO)算法,并將其與其他算法進(jìn)行仿真試驗對比,結(jié)果表明,ECMPSO算法在兼顧調(diào)度方案多樣性的同時具有更高的求解精度,求出的調(diào)度方案分布均勻, 可為決策者提供更多有價值的候選調(diào)度方案.

[1] DEB K, PRATAP A, AGARWAL S, et al. A fast and elitist multi-objective genetic algorithm: NSGA-Ⅱ[J]. IEEE Transactions on Evolutionary Computation, 2002,6(2):182-197.

[2] 伍大清,鄭建國. 基于混合策略自適應(yīng)學(xué)習(xí)的并行粒子群優(yōu)化算法[J].控制與決策,2013,28(7):1087-1093.

[3] 劉貞,張希良, 張繼紅,等.排污總量控制下的電力交易市場與排污權(quán)交易市場分析[J].電力系統(tǒng)保護(hù)與控制,2009,37(22):4-8.

[4] 盧有麟,周建中,覃暉,等.差分進(jìn)化算法在電力系統(tǒng)環(huán)境經(jīng)濟(jì)調(diào)度中的應(yīng)用[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2010(8):121-124.

[5] LIU G, PENG C, XIANG L. Economic-environmental dispatch using improved multi-objective particle swarm optimization [J]. Power System Technology, 2011, 35(7):139-144.

[6] HU W, LI Z S. A simpler and more effective particle swarm optimization algorithm [J]. Journal of Software, 2007, 18(4):861-868.

[7] WU D Q, ZHENG Z G. A dynamic multistage hybrid swarm intelligence optimization algorithm for function optimization[J]. Discrete Dynamics in Nature and Society, 2012, doi:10.1155/2012/578064.

[8] ZOU W P, ZHU Y L, CHEN H N, et al. Solving multi-objective optimization problems using artificial bee colony algorithm [J]. Discrete Dynamics in Nature and Society, 2011, doi:10.1155/2011/569784.

[9] HUANG V L, SUGANTHAN P N, LIANG J. Comprehensive learning particle swarm optimizer for solving multi-objective optimization problems [J]. International Journal of Intelligent Systems, 2006, 21(2):209-226.

[10] 李密青,鄭金華.一種多目標(biāo)進(jìn)化算法解集分布廣度評價方法[J].計算機(jī)學(xué)報,2011,34(4):647-664.

[11] HEMAMALINI S, SIMON P S. Emission constrained economic dispatch with valve-point effect using particle swarm optimization [C]//TENCON 2008 IEEE Region 10 Conference,India,2008.

[12] 徐麗青, 吳亞麗.求解環(huán)境經(jīng)濟(jì)調(diào)度問題的多目標(biāo)差分粒子群優(yōu)化算法[J].西安理工大學(xué)學(xué)報,2011,27(1):62-68.

[13] 瞿博陽,梁靜,SUGANTHAN P N.雙局部粒子群算法解決環(huán)境經(jīng)濟(jì)調(diào)度問題[J].計算機(jī)工程與應(yīng)用,2014,5(11):45-53.

Power System Optimization Based on Multi-swarm Cooperation Collaborative Algorithm

LIUYong1,WUDa-qing1, 2, 3, 4, 5,SUNLi4,YANGZhi-ping4

(1. Artificial Intelligence Key Laboratory of Sichuan Province, Zigong 643000, China; 2. Computer Science and Technology Institute, University of South China, Hengyang 421001, China; 3. Key Laboratory of Intelligent Computing & Signal Processing, Ministry of Education, Anhui University, Hefei 230039, China; 4. Glorious Sun School of Business and Management, Donghua University, Shanghai 200051, China; 5. Key Laboratory of Agri-Informatics, Ministry of Agriculture, Beijing 100081, China)

An effective co-evolutionary multi-swarm particle swarm optimization (ECMPSO) is proposed for environmental-economic dispatch (EED), which is a complicated multi-objective constrained problem. The strategies such as multi-swarm, new particle swarm optimization (PSO) update as well as restructuring topology structure strategy with dynamic population in period are introduced in ECMPSO, so the algorithm has strong global search ability and good robust performance. Some typical multi-objective optimization functions are tested to verify the algorithm, the results show that the proposed algorithm has superior performance of faster convergence speed and strong ability to jump out of local optimum. The proposed algorithm is applied to solve the EED problem of a power system, experimental results demonstrate the successful application of ECMPSO, compared with other methods, ECMPSO outperforms in globe convergence and the diversity of the Pareto set with higher effectiveness and robustness.

cooperation collaborative optimization; multi-swarm optimization; particle swarm optimization; environmental economic dispatch

1671-0444(2016)01-0110-07

2014-09-25

四川省人工智能重點實驗室基金資助項目(2015RYJ01);湖南省科技計劃資助項目(2015JC3089);教育部計算智能與信號處理重點實驗室基金資助項目(2014C16304);上海市自然科學(xué)基金資助項目(15ZR1401600);農(nóng)業(yè)部農(nóng)業(yè)信息服務(wù)技術(shù)重點實驗室基金資助項目(2015-AIST-02);四川省模式識別與智能信息處理重點實驗室基金資助項目

劉 勇(1981—),男,四川自貢人,實驗師,碩士,研究方向為復(fù)雜系統(tǒng)優(yōu)化.E-mail:24297633@qq.com

伍大清(聯(lián)系人),女,講師,E-mail:dqw_1982@126.com

TM 734

A

猜你喜歡
適應(yīng)度種群粒子
改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
山西省發(fā)現(xiàn)刺五加種群分布
基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
中華蜂種群急劇萎縮的生態(tài)人類學(xué)探討
紅土地(2018年7期)2018-09-26 03:07:38
基于粒子群優(yōu)化極點配置的空燃比輸出反饋控制
基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
中國塑料(2016年11期)2016-04-16 05:26:02
基于Matlab的α粒子的散射實驗?zāi)M
物理與工程(2014年4期)2014-02-27 11:23:08
崗更湖鯉魚的種群特征
基于兩粒子糾纏態(tài)隱形傳送四粒子GHZ態(tài)
少數(shù)民族大學(xué)生文化適應(yīng)度調(diào)查
类乌齐县| 林口县| 芷江| 泰兴市| 乌拉特前旗| 双桥区| 青神县| 南阳市| 麦盖提县| 阿拉善左旗| 台湾省| 丹阳市| 华容县| 石河子市| 巫溪县| 延川县| 牙克石市| 随州市| 锦屏县| 湘阴县| 乐都县| 习水县| 仙桃市| 卢湾区| 通州区| 同心县| 淮滨县| 岳西县| 定襄县| 普兰县| 驻马店市| 平陆县| 苍溪县| 崇礼县| 江西省| 三都| 庆云县| 永善县| 大厂| 西和县| 绿春县|