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

?

基于貝葉斯網(wǎng)絡(luò)對全國PM2.5濃度影響因素分析

2020-07-26 14:19:51李超群
科技創(chuàng)新與應(yīng)用 2020年22期
關(guān)鍵詞:網(wǎng)絡(luò)結(jié)構(gòu)貝葉斯節(jié)點

李超群

摘? 要:為能更有針對性的控制PM2.5濃度,對2000-2017年間全國31個省市PM2.5濃度數(shù)值和由專家先驗得出的影響PM2.5的六種人為因素分別建立了基于BIC評分函數(shù)、K2評分函數(shù)進(jìn)行結(jié)構(gòu)學(xué)習(xí)的兩種貝葉斯網(wǎng)絡(luò)模型、支持向量機(jī)模型、K-近鄰模型進(jìn)行分析。利用5-折交叉驗證對四種模型進(jìn)行評估。發(fā)現(xiàn)在樣本量不太大的情況下,貝葉斯網(wǎng)絡(luò)表現(xiàn)出更好地穩(wěn)健性與優(yōu)越性,而基于K2評分函數(shù)進(jìn)行結(jié)構(gòu)學(xué)習(xí)的貝葉斯網(wǎng)絡(luò)模型具有更好地分類性能。為政府相關(guān)部門對我國PM2.5濃度更加有效的控制,以及采取更加具有針對性的治理方案提供了思路與方案。

關(guān)鍵詞:PM2.5;貝葉斯網(wǎng)絡(luò);結(jié)構(gòu)學(xué)習(xí);參數(shù)學(xué)習(xí);模型評估

中圖分類號:TP18? ? ? ? ? 文獻(xiàn)標(biāo)志碼:A? ? ? ? ?文章編號:2095-2945(2020)22-0001-06

Abstract: For controlling the concentration of the PM2.5, with the data of PM2.5 concentration values of 31 provinces in China during theyears 2000 to 2017, and six influencing factors of PM2.5 concentration obtained by experts' prior experience. Two Bayesian Network Models based on BIC score function and K2 score function respectively, Support Vector Machine, K-Nearest Neighbor Model are analyzed. Using the method of 5-fold cross validation, the four models are trained and tested. In conclusion, in the case of small sample size, the Bayesian Network shows better robustness and superiority, while the Bayesian Network Model based on K2 score function for structural learning has better performance. It is meaningful to provide a new idea and scheme for the control of PM2.5 concentration in China.

Keywords: PM2.5; bayesian network; structure learning; parameter learning; assessment of models

引言

近幾年,隨著我國不斷發(fā)展,城市化的程度不斷提高,在有了越來越多的機(jī)遇同時,也面臨著環(huán)境的巨大挑戰(zhàn)。我國PM2.5數(shù)值頻繁“爆表”,長期的霧霾天不僅嚴(yán)重影響了我國正常的生產(chǎn)進(jìn)程,還影響了人們正常的生活。因此,對PM2.5的影響因素進(jìn)行分析,就顯得十分迫切,但是PM2.5濃度的變化是一個不確定的、復(fù)雜的過程,這也增加了建模過程以及分析的難度。近年來,許多專家學(xué)者致力于對PM2.5影響因素進(jìn)行分析。國內(nèi)學(xué)者解蕾等[1]使用二元定距變量的相關(guān)分析,分析了兩兩變量之間的關(guān)聯(lián)程度,最終得出結(jié)論PM2.5與SO2,PM2.5與CO,PM2.5與PM10之間都存在著較強(qiáng)的相關(guān)性,之后建立PM2.5與PM10的線性回歸方程,進(jìn)一步分析二者之間的線性關(guān)系。學(xué)者盧德彬等[2]采用Theil-Sen median趨勢分析,Mann-Kendall、R/S和相關(guān)分析方法,分析了PM2.5的時空格局、空間變化特征、污染來源。學(xué)者羅毅等[3]將BP神經(jīng)網(wǎng)絡(luò)與支持向量機(jī)相結(jié)合,構(gòu)建了PM2.5濃度組合的模擬模型,分析了PM2.5污染規(guī)律和趨勢。學(xué)者王娟[4]利用灰色關(guān)聯(lián)度分析,找到了主要影響因素,建立主要影響因素與PM2.5的二次多元回歸方程,綜合考慮了各因素的影響。雖然上述模型均取得了不錯的結(jié)果,然而,利用上述方法在對PM2.5進(jìn)行分析時仍存在著一些不足:(1)將所有可能影響因素拆分開來單獨分析,僅考慮單因素的影響,而忽略了多因素之間的交互作用;(2)相關(guān)分析只能反映出兩個或者多個因素之間的線性相關(guān)程度,而不能得出具體的因果關(guān)系;(3)由于PM2.5濃度變化的復(fù)雜性,各因素與其之間的關(guān)系是不確定的,多元回歸模型可能并不能反映出其變化的真實規(guī)律;(4)大多文獻(xiàn)

側(cè)重于對影響PM2.5濃度的自然因素進(jìn)行分析,而簡化了對人為源的分析?;诖耍疚睦靡延械膶<蚁闰炛R,找出影響我國PM2.5濃度的6個主要人為因素,分別建立基于K2評分進(jìn)行結(jié)構(gòu)學(xué)習(xí)的貝葉斯網(wǎng)絡(luò)以及根據(jù)BIC評分進(jìn)行結(jié)構(gòu)學(xué)習(xí)的兩種貝葉斯網(wǎng)絡(luò)模型,對在各因素共同作用下,PM2.5數(shù)值是否能達(dá)到我國二級標(biāo)準(zhǔn)進(jìn)行判斷,同時建立支持向量機(jī)(SVM)、K-近鄰(KNN)模型,利用5-折交叉驗證對四種模型的分類精度進(jìn)行對比,考察四種模型在樣本數(shù)目不太大的情況下的分類精度,找出更為有效的單一分類模型。通過分析影響PM2.5濃度的直接因素、間接因素或無關(guān)因素,據(jù)此可以為政府相關(guān)部門控制PM2.5濃度提供更加有針對性,更為有效的建議方法。

貝葉斯網(wǎng)絡(luò)是由節(jié)點以及有向邊組成的概率圖模型,是一種在不確定、不完整的因素影響下,進(jìn)行推理的一種有效的工具。貝葉斯網(wǎng)絡(luò)提供了一種展示變量之間的因果關(guān)系的框架結(jié)構(gòu)。[5]目前貝葉斯網(wǎng)絡(luò)被廣泛應(yīng)用于故障分析以及相關(guān)領(lǐng)域。例如,在故障分析領(lǐng)域,李爽等[6]基于BN-ELM方法對煤礦瓦斯安全態(tài)勢進(jìn)行研究;熊宇峰等[7]借助樹形貝葉斯網(wǎng)絡(luò),實現(xiàn)對配電網(wǎng)故障的快速災(zāi)情推斷,Hu[8]等借助貝葉斯網(wǎng)絡(luò)對地震液化勢進(jìn)行預(yù)測。在金融領(lǐng)域,嚴(yán)冠等[9]等利用貝葉斯方法建立銀行同業(yè)借貸網(wǎng)絡(luò),對其中的系統(tǒng)風(fēng)險進(jìn)行分析研究。

1 研究區(qū)域與數(shù)據(jù)來源

本文面向我國31個省市,收集了2000-2017年相關(guān)數(shù)據(jù)558條。其中,我國2000-2017年的PM2.5值年度數(shù)據(jù),由全國城市空氣質(zhì)量實時發(fā)布平臺獲得。根據(jù)專家知識,得到影響PM2.5濃度的6個人為因素,包括:城市化水平(Urbanize)、產(chǎn)業(yè)結(jié)構(gòu)(IS)、能源結(jié)構(gòu)(ES)、外商直接投資額(Open)、人均生產(chǎn)總值(GRP)、政府治理(PT)。本文使用城鎮(zhèn)居民的消費水平與居民消費水平的比值表示城市化水平;產(chǎn)業(yè)結(jié)構(gòu)用當(dāng)?shù)氐诙a(chǎn)業(yè)增加值占該地區(qū)生產(chǎn)總值的比重來表示;用該地區(qū)煤炭消耗量與該地區(qū)生產(chǎn)總值的比值反映能源結(jié)構(gòu),外商直接投資額以外商投資企業(yè)進(jìn)出口總額與地區(qū)生產(chǎn)總值的比值來衡量;由政府治理廢氣項目完成投資占工業(yè)污染治理完成投資的比重來衡量政府治理的力度。相關(guān)數(shù)據(jù)來自《中國統(tǒng)計年鑒》。對數(shù)據(jù)的統(tǒng)計描述如圖1。

從圖中可以看出,地區(qū)人均生產(chǎn)總值存在較多離群點,且各點之間差距較大,這也從側(cè)面反映出我國的經(jīng)濟(jì)發(fā)展仍然是不均衡的,地區(qū)之間的經(jīng)濟(jì)水平差異較大。此外,可以看到外商直接投資也存在較多離群點,根據(jù)我們的認(rèn)識,經(jīng)濟(jì)水平更高、地理位置更優(yōu)越的地區(qū)可以吸引更多的外商投資,因此造成了該樣本值分布不均、差異較大。為了更加真實的反應(yīng)各因素對PM2.5濃度的影響,本文將離群點全部保留。

2 研究方法

本文根據(jù)歷史數(shù)據(jù)以及專家知識,利用5-折交叉檢驗方法,對兩種貝葉斯網(wǎng)絡(luò)、SVM模型、K-近鄰,四種分類模型進(jìn)行訓(xùn)練與檢驗,并對四種模型在樣本數(shù)目不太大情況下的分類的準(zhǔn)確性進(jìn)行比較,得出更有效的分類模型。

2.1 貝葉斯網(wǎng)絡(luò)

一個完整的貝葉斯網(wǎng)絡(luò)模型由貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)、條件概率表兩部分組成。其網(wǎng)絡(luò)結(jié)構(gòu)是一個有向無環(huán)圖,[10]其不同節(jié)點代表了不同的變量,節(jié)點之間的有向弧由父節(jié)點指向子節(jié)點,表示子節(jié)點的取值會受到父節(jié)點取值的直接影響。

2.1.1 概率推理

與傳統(tǒng)的概率派直接根據(jù)樣本數(shù)據(jù)建模分析不同,貝葉斯派為了更加充分利用先驗信息,首先引入?yún)?shù)先驗分布,再根據(jù)樣本數(shù)據(jù),計算得到參數(shù)的后驗概率,從而加深對先驗信息的認(rèn)識,在不斷迭代,重復(fù)上述方法的過程中,不斷得到新的先驗概率,再將更新修正后的先驗信息繼續(xù)用于求參數(shù)的后驗概率,從而實現(xiàn)更加準(zhǔn)確的推斷。通過樣本數(shù)據(jù),可以對一個給出節(jié)點的位置與其條件概率分布的貝葉斯網(wǎng)絡(luò)重新估計其變量的概率分布,其公式基礎(chǔ)便是如下貝葉斯公式與乘法公式:

其中,P(X)為先驗概率,P(Y)為邊際分布,P(Y|X)是由樣本數(shù)據(jù)得到的似然概率。由于貝葉斯網(wǎng)絡(luò)的聯(lián)合分布具有馬爾可夫性,即變量只與其父節(jié)點取值有關(guān),因此公式(2)可以寫成(3)的形式,π(xi)是變量xi的父節(jié)點概率。

2.1.2 數(shù)據(jù)處理

由于貝葉斯網(wǎng)絡(luò)對于處理離散數(shù)據(jù)具有更強(qiáng)的穩(wěn)健性,因此先對數(shù)據(jù)進(jìn)行離散化處理。按照PM2.5的濃度是否達(dá)到我國二級標(biāo)準(zhǔn),即年均值是否小于35μg/m3。將均值小于等于35μg/m3賦值1,否則賦值為0。由于外商直接投資額、政府治理兩個因素含有大量缺失數(shù)據(jù),因此本文將這兩個因素中的缺失數(shù)據(jù)賦值為0,非缺失數(shù)據(jù)賦值為1,以考察這兩個因素的有無對于PM2.5濃度能否達(dá)標(biāo)的影響。對于其他變量的缺失值,本文用插補(bǔ)法進(jìn)行插補(bǔ),之后利用k均值的方法將變量人均生產(chǎn)總值、能源結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、城市化水平進(jìn)行聚類,其聚類中心如表1:

表1 變量聚類中心

離散化后的結(jié)果如表2所示。

2.1.3 結(jié)構(gòu)學(xué)習(xí)

貝葉斯網(wǎng)絡(luò)進(jìn)行結(jié)構(gòu)學(xué)習(xí)的方法主要分為兩大類:(1)基于約束的結(jié)構(gòu)學(xué)習(xí)法。這類方法的核心思想是:首先對訓(xùn)練數(shù)據(jù)集進(jìn)行條件獨立性測試等統(tǒng)計測試,[11]從而找出變量間的條件獨立關(guān)系。DAG圖就是按照盡可能多地覆蓋條件獨立性關(guān)系的原則所構(gòu)造出來的。[12](2)基于評分搜索的結(jié)構(gòu)學(xué)習(xí)方法。一種直接的想法是,希望能夠遍歷DAG結(jié)構(gòu)中的每一種情況,根據(jù)某種評價準(zhǔn)則,去對每一種可能情況進(jìn)行評分,選擇該標(biāo)準(zhǔn)下的最優(yōu)結(jié)構(gòu)模型。然而隨著貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)模型中節(jié)點數(shù)的增加,DAG結(jié)構(gòu)的數(shù)量將會大幅增加,Robinson曾給出一個計算DAG結(jié)構(gòu)數(shù)量的計算公式:

由公式可以看出,隨著節(jié)點數(shù)增加,DAG結(jié)構(gòu)的數(shù)目的增加速度是指數(shù)級別的,因此,在節(jié)點數(shù)比較多的情況下,想要遍歷每一種結(jié)構(gòu)形式在有限的時間內(nèi)是很難完成的。為了解決這一問題,許多學(xué)者便提出了利用評分函數(shù)的方法,例如爬山算法,面向完全數(shù)據(jù)模型的K2算法、以及面向含有少量缺失數(shù)據(jù)的EM算法以及馬爾可夫-蒙特卡洛(MCMC)算法,用評分函數(shù)作為準(zhǔn)則,對各結(jié)構(gòu)進(jìn)行打分,進(jìn)而選擇出得分最高的結(jié)構(gòu),該結(jié)構(gòu)即為該準(zhǔn)則下最優(yōu)。評分函數(shù)主要分為兩大類,一類是包括BDeu評分、K2評分等的貝葉斯評分函數(shù)。假設(shè)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)G的先驗概率為P(G) ,針對給定樣本集D,根據(jù)貝葉斯公式,網(wǎng)絡(luò)結(jié)構(gòu)G的后驗概率可以表示為:

由于P(D)為樣本的函數(shù),與拓?fù)浣Y(jié)構(gòu)無關(guān),因此對式(5)求極值等價于對其分子P(G)P(D|G)求極值,使得分子達(dá)到最大的G,即為所求的,具有最大后驗概率的網(wǎng)絡(luò)結(jié)構(gòu)。為了計算的方便,定義

logP(G,D)=log(P(G)P(D|G))=logP(G)+logP(D|G)(6)

為網(wǎng)絡(luò)結(jié)構(gòu)的貝葉斯評分,即為MAP測度。[13]另一類是基于信息論的評分函數(shù),BIC評分、MDL評分以及AIC評分等都是使用頻率較高的評分函數(shù)。爬山算法是貪婪算法的一種,通過不斷迭代最大程度提高分?jǐn)?shù),一旦找到局部最大值,搜索將終止,并返回相應(yīng)的局部最優(yōu)的結(jié)果。本文以爬山算法作為搜索策略,分別選擇貝葉斯評分函數(shù)中的K2評分函數(shù),基于信息論的BIC評分函數(shù)對模型進(jìn)行打分,結(jié)合專家先驗知識確定出兩個局部最優(yōu)的DAG網(wǎng)絡(luò)結(jié)構(gòu)。為了表述的方便,下文分別將基于K2評分函數(shù)的貝葉斯網(wǎng)絡(luò)模型與基于BIC評分的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)模型記為BN-K2模型與BN-BIC模型。兩種貝葉斯網(wǎng)絡(luò)模型的結(jié)構(gòu)圖如下圖2:

其中,P表示精確率,R表示召回率。

AUC(Area Under Curve)指ROC曲線下的面積,取值介于0.1到1之間,并且越接近于1表示模型的性能越好?;诨煜仃嚨母髦笜?biāo)的計算結(jié)果匯總?cè)绫?。

由表7可以看出,BN-K2模型優(yōu)勢比較明顯,各項得分均為最高。根據(jù)F1得分可以看出,BN-K2模型得分最高,為60.366%。根據(jù)AUC的值表現(xiàn)出的結(jié)果與F1得分類似,BN-K2模型取值為0.546,是表現(xiàn)最好的模型。支持向量機(jī)模型的AUC值最低,只有0.509。

3.2 討論

綜合考慮5-折交叉驗證得分與基于混淆矩陣的評分指標(biāo),可以得出結(jié)論,兩種貝葉斯網(wǎng)絡(luò)模型是四種模型中最為穩(wěn)健的,并且模型對于先驗知識的利用,也使其在樣本量不大的情況下,有著更為優(yōu)秀的表現(xiàn)。對樣本量不大的問題仍有著較為不錯的表現(xiàn)也是貝葉斯網(wǎng)絡(luò)模型的優(yōu)勢之一。

由于樣本量以及先驗信息的限制,模型能夠從數(shù)據(jù)中學(xué)習(xí)到的信息并不夠多,因此模型的預(yù)測精度并不高,還需要進(jìn)一步的提高。接下來的工作中,需要改進(jìn)提升的地方還有很多,例如:(1)增大樣本量,使模型能夠從中得到充分的信息,以增加模型的精度;(2)在對數(shù)據(jù)缺失值處理的時候,可以采用更加有效的處理方式;(3)收集更多的專家先驗知識,對模型的結(jié)構(gòu)與參數(shù)進(jìn)行調(diào)整與優(yōu)化;(4)能夠使用混合模型對影響因素進(jìn)行更加深入的分析。

4 結(jié)論與建議

本文對我國2000-2017年全國各年平均地PM2.5濃度值,以及根據(jù)專家先驗得出的六種PM2.5濃度值的影響因素分別建立了BIC-BN模型、K2-BN模型、SVM模型、K-近鄰模型四種模型。根據(jù)BIC-BN模型的DAG結(jié)構(gòu)圖(圖2),我們可以看到,地區(qū)產(chǎn)業(yè)結(jié)構(gòu)是影響PM2.5濃度的直接因素,城市化水平是影響其濃度的間接因素;K2-BN模型得出的DAG結(jié)構(gòu)圖復(fù)雜得多,揭示了更多的直接影響因素。我們可以發(fā)現(xiàn),地區(qū)的人均生產(chǎn)總值、能源結(jié)構(gòu)城市化水平、政府治理與外商直接投資額都會對PM2.5濃度產(chǎn)生直接影響。借助5-折交叉驗證與AUC對模型進(jìn)行評估,得出結(jié)論,即在樣本量不太大的情況下,兩種貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)模型穩(wěn)健性更好,模型的性能也較好,而在兩種貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)模型中,基于K2評分的貝葉斯網(wǎng)絡(luò)模型性能更佳。

最后,根據(jù)本文分析的結(jié)果,綜合考慮兩種模型。各個政府要實現(xiàn)對PM2.5濃度的控制,實現(xiàn)達(dá)到國家二級標(biāo)準(zhǔn)的目標(biāo),需要在大力發(fā)展地區(qū)經(jīng)濟(jì)水平的同時,要重視調(diào)整產(chǎn)業(yè)結(jié)構(gòu)與能源結(jié)構(gòu)使其變得更加合理,例如提高第三產(chǎn)業(yè)的占比,減少煤炭的使用量,增加清潔能源的用量,加大環(huán)境治理的力度。根據(jù)貝葉斯網(wǎng)絡(luò)模型得出的有關(guān)PM2.5濃度影響因素的結(jié)論,為地方政府采取措施提供了更加有效的方向,地方政府可以采取更加有針對性的治理措施,更加有效的解決大氣污染問題。

參考文獻(xiàn):

[1]解蕾,狄光智.基于R語言的城市PM2.5影響因素分析[J].軟件工程,2019,22(05):15-17+8.

[2]盧德彬,毛婉柳,楊東陽,等.基于多源遙感數(shù)據(jù)的中國PM2.5變化趨勢與影響因素分析[J].長江流域資源與環(huán)境,2019,28(03):651-660.

[3]羅毅,鄧瓊飛,楊昆,等.近20年來中國典型區(qū)域PM2.5時空演變過程[J].環(huán)境科學(xué),2018,39(07):3003-3013.

[4]王娟.基于多元回歸分析的PM2.5預(yù)測研究[J].微型電腦應(yīng)用,2020,36(03):48-51.

[5]丁艷麗,楊敏,楊殿微.型號研制可靠性工作項目轉(zhuǎn)階段風(fēng)險預(yù)警[J].項目管理技術(shù),2009,7(07):37-40.

[6]李爽,李丁煒,猶夢潔,等.基于BN-ELM的煤礦瓦斯安全態(tài)勢預(yù)測方法研究[J].系統(tǒng)工程,2020,38(03):132-140.

[7]熊宇峰,周剛,陳穎,等.基于樹形貝葉斯網(wǎng)絡(luò)的配電網(wǎng)快速災(zāi)情推斷[J].電網(wǎng)技術(shù),2020,44(06):2222-2230.

[8]Ji-Lei Hu,Xiao-Wei Tang,Jiang-Nan Qiu. Assessment of seismic liquefaction potential based on Bayesian network constructed from domain knowledge and history data[J].Soil Dynamics and Earthquake Engineering,2016,89.

[9]嚴(yán)冠,劉志東.基于貝葉斯方法的中國商業(yè)銀行同業(yè)借貸網(wǎng)絡(luò)中系統(tǒng)風(fēng)險研究[J].中國管理科學(xué),2020,28(04):14-26.

[10]解晶,何楨,馮楠.基于貝葉斯網(wǎng)絡(luò)的計算機(jī)硬件保修成本分析[J].計算機(jī)工程與應(yīng)用,2007(24):104-106.

[11]晏文娟.模糊貝葉斯網(wǎng)絡(luò)的研究及其在電子商務(wù)發(fā)展水平中的應(yīng)用[D].華南理工大學(xué),2009.

[12]李瑋瑋.貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)方法的研究[D].南京航空航天大學(xué),2009.

[13]胡春玲.貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)及其應(yīng)用研究[D].合肥工業(yè)大學(xué),2011.

猜你喜歡
網(wǎng)絡(luò)結(jié)構(gòu)貝葉斯節(jié)點
CM節(jié)點控制在船舶上的應(yīng)用
Analysis of the characteristics of electronic equipment usage distance for common users
基于AutoCAD的門窗節(jié)點圖快速構(gòu)建
貝葉斯公式及其應(yīng)用
基于貝葉斯估計的軌道占用識別方法
基于互信息的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)
知識網(wǎng)絡(luò)結(jié)構(gòu)維對于創(chuàng)新績效的作用機(jī)制——遠(yuǎn)程創(chuàng)新搜尋的中介作用
滬港通下A+ H股票網(wǎng)絡(luò)結(jié)構(gòu)演化的實證分析
一種基于貝葉斯壓縮感知的說話人識別方法
電子器件(2015年5期)2015-12-29 08:43:15
復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)比對算法研究進(jìn)展
鲁甸县| 麻栗坡县| 前郭尔| 奇台县| 沐川县| 广宗县| 白水县| 上虞市| 定日县| 嘉定区| 平江县| 石景山区| 凤冈县| 日土县| 正定县| 凉山| 霍邱县| 丰顺县| 江永县| 新郑市| 太仓市| 广饶县| 五大连池市| 镇安县| 昌乐县| 交城县| 合作市| 方城县| 大同市| 东源县| 酉阳| 崇左市| 三江| 朝阳县| 梁河县| 木里| 通化市| 崇明县| 佛学| 扬中市| 同江市|