精品人妻系列无码人妻漫画,久久精品国产一区二区三区,国产精品无码专区,无码人妻少妇伦在线电影,亚洲人妻熟人中文字幕一区二区,jiujiuav在线,日韩高清久久AV

中國(guó)儲(chǔ)能網(wǎng)歡迎您!
當(dāng)前位置: 首頁(yè) >數(shù)智化儲(chǔ)能系統(tǒng)>產(chǎn)品檢測(cè)與認(rèn)證 返回

基于深度學(xué)習(xí)的直流微電網(wǎng)虛假數(shù)據(jù)注入攻擊二階段檢測(cè)方法

作者:陶磊 羅萍萍 林濟(jì)鏗 來(lái)源:中國(guó)電力 發(fā)布時(shí)間:2024-10-09 瀏覽:

中國(guó)儲(chǔ)能網(wǎng)訊:

     摘要

  直流微電網(wǎng)是一個(gè)網(wǎng)絡(luò)物理信息系統(tǒng),在信息傳遞的過(guò)程中容易遭受網(wǎng)絡(luò)攻擊的影響。虛假數(shù)據(jù)注入信息通道會(huì)影響微電網(wǎng)的系統(tǒng)安全。檢測(cè)并修正虛假數(shù)據(jù)注入攻擊,能夠提升微電網(wǎng)系統(tǒng)運(yùn)行的安全性。針對(duì)這一問(wèn)題,提出了一種基于卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network,CNN)和長(zhǎng)短期記憶網(wǎng)絡(luò)(long short-term memory,LSTM)聯(lián)合最大互信息系數(shù)(maximum information coefficient,MIC)的二階段虛假數(shù)據(jù)注入攻擊檢測(cè)方法。首先,使用CNN從直流微電網(wǎng)運(yùn)行的時(shí)序數(shù)列中提取時(shí)序特征,運(yùn)用LSTM模型結(jié)合CNN提取的時(shí)序特征運(yùn)行得到直流微電網(wǎng)運(yùn)行狀態(tài)預(yù)測(cè)值,與直流微電網(wǎng)運(yùn)行的實(shí)際值對(duì)比,初步判斷系統(tǒng)中是否存在虛假數(shù)據(jù);其次,考慮到CNN-LSTM模型存在一定的誤報(bào)率,構(gòu)建MIC校驗(yàn)器,進(jìn)一步判斷系統(tǒng)中是否存在虛假數(shù)據(jù)并恢復(fù);最后,通過(guò)直流微電網(wǎng)Matlab仿真分析,驗(yàn)證了所提方法的合理性和可行性。

  1 直流微電網(wǎng)控制模型

  1.1 直流微電網(wǎng)系統(tǒng)描述

  本文考慮的直流微電網(wǎng)包含m個(gè)分布式發(fā)電單元(DG)和n個(gè)負(fù)載,它們通過(guò)并聯(lián)的方式連接在直流母線上,并且通過(guò)下垂控制和二次控制維持系統(tǒng)的穩(wěn)定。圖1展示了本文研究的直流微電網(wǎng)的一般網(wǎng)絡(luò)物理信息系統(tǒng)。

圖1 直流微電網(wǎng)控制模型

Fig.1 Control model of DC microgrid

  1.2 直流微電網(wǎng)控制模型

  為確保直流微電網(wǎng)的穩(wěn)定運(yùn)行,通常采用下垂控制實(shí)現(xiàn)各個(gè)DG間的電流分配。由下垂控制方程得到參考電壓Uref為

  1.3 虛假數(shù)據(jù)對(duì)直流微電網(wǎng)的影響

  針對(duì)直流微電網(wǎng)的攻擊通常是DoS攻擊和FDIA,F(xiàn)DIA相較于DoS攻擊而言有著更高的隱蔽性和破壞性。本文對(duì)FDIA進(jìn)行更深入的討論。在微電網(wǎng)中,當(dāng)攻擊者通過(guò)通信通道對(duì)系統(tǒng)注入虛假數(shù)據(jù)時(shí),會(huì)導(dǎo)致信息層發(fā)出錯(cuò)誤指令,造成控制系統(tǒng)紊亂,對(duì)整個(gè)微電網(wǎng)網(wǎng)絡(luò)的穩(wěn)定造成嚴(yán)重影響。

  虛假數(shù)據(jù)信號(hào)注入表現(xiàn)形式以疊加的方式注入鄰居間的信息傳輸通道中,因此控制器遭受FDIA時(shí),微電網(wǎng)受攻擊模型可以表示為

  2 虛假數(shù)據(jù)注入攻擊檢測(cè)模型

  2.1 CNN-LSTM檢測(cè)模型

  CNN是一種深度學(xué)習(xí)模型,采用局部連接和共享權(quán)值的方式,直接從原始數(shù)據(jù)中通過(guò)卷積層和池化層交替使用來(lái)獲取有效表征,自動(dòng)提取數(shù)據(jù)的局部特征,并建立特征向量。本文使用CNN模型提取數(shù)據(jù)特征值。

  LSTM在隱藏層中加入記憶單元控制時(shí)間序列數(shù)據(jù)的記憶信息。信息在隱藏層的不同細(xì)胞之間通過(guò)幾個(gè)可控的門(mén)(遺忘門(mén)、輸入門(mén)、輸出門(mén))傳遞,從而可以控制前一信息和當(dāng)前信息的記憶和遺忘程度。LSTM具有長(zhǎng)期記憶功能,可以避免梯度消失問(wèn)題。

  本文檢測(cè)模型所使用的數(shù)據(jù)從仿真的微電網(wǎng)運(yùn)行中獲得,其中負(fù)荷在不斷變化,此時(shí)收集數(shù)據(jù)的時(shí)間為15 s,采樣的時(shí)間為5 ms,可以得到3000組數(shù)據(jù),其中80%為訓(xùn)練集,20%為測(cè)試集,所要預(yù)測(cè)的目標(biāo)為下一時(shí)刻的各DG單元的電流電壓值。本文的CNN-LSTM輸入數(shù)據(jù)為多元時(shí)間序列,采用滑動(dòng)窗口的方式將多元時(shí)間序列依次輸入,滑動(dòng)窗口的寬度設(shè)置為20條記錄,步長(zhǎng)設(shè)置為1,輸入的特征如圖2所示。圖2中t為時(shí)間尺度,T為每5 ms中的采樣時(shí)刻,T+n為T(mén)時(shí)刻后第n個(gè)5 ms的時(shí)刻。

圖2 CNN-LSTM網(wǎng)絡(luò)模型輸入數(shù)據(jù)結(jié)構(gòu)

Fig.2 Input data structure of CNN-LSTM network model

  本文所提出的CNN-LSTM網(wǎng)絡(luò)模型結(jié)構(gòu)如圖3所示。2個(gè)卷積核的數(shù)目分別為16和32,池化層池的大小為2×1,dropout大小為0.2。2層LSTM神經(jīng)元個(gè)數(shù)為16、32。最后通過(guò)全連接層輸出指定格式向量,即下一時(shí)刻的電流和電流值。

圖3 CNN-LSTM檢測(cè)模型流程

Fig.3 Flow chart of CNN-LSTM detection model

  2.2 MIC校驗(yàn)器

  直流微電網(wǎng)網(wǎng)絡(luò)層遭受虛假數(shù)據(jù)注入攻擊時(shí),改變了直流微電網(wǎng)電流分配比例和電壓的穩(wěn)定,而直流微電網(wǎng)負(fù)荷變化也會(huì)同時(shí)改變直流微電網(wǎng)電壓和電流的變化。如何區(qū)別是由正常系統(tǒng)變化引起的還是由FDIA引起的電壓電流變化是檢測(cè)模塊所需要完成的任務(wù)。MIC是衡量2個(gè)變量間關(guān)系相關(guān)性的一種方法。當(dāng)系統(tǒng)中存在虛假數(shù)據(jù)時(shí),系統(tǒng)中的相關(guān)變量之間的相關(guān)性會(huì)在一定程度上解耦,因此可用MIC方法判斷系統(tǒng)中是否存在虛假數(shù)據(jù)。

  2.3 二階段FDIA檢測(cè)模型

  基于CNN-LSTM聯(lián)合MIC的二階段檢測(cè)機(jī)制如圖4所示。離線訓(xùn)練時(shí),選取正常運(yùn)行狀態(tài)下的微電網(wǎng)歷史運(yùn)行數(shù)據(jù)集x?Dtrain,CNN-LSTM檢測(cè)模型通過(guò)訓(xùn)練擬合得到準(zhǔn)確的預(yù)測(cè)值;與此同時(shí),該歷史數(shù)據(jù)集計(jì)算MIC校驗(yàn)器正常運(yùn)行狀態(tài)下電流與電壓誤差之間的相關(guān)性,得到相關(guān)性數(shù)值,用以后續(xù)FDIA檢驗(yàn)。測(cè)試時(shí),首先將直流微電網(wǎng)運(yùn)行時(shí)的量測(cè)值傳輸至CNN-LSTM檢測(cè)模型中,當(dāng)實(shí)際運(yùn)行值與預(yù)測(cè)值相匹配時(shí),檢測(cè)系統(tǒng)正常運(yùn)轉(zhuǎn)不動(dòng)作,輸出為0代表正常狀態(tài);在二者間存在誤差大于設(shè)定閾值時(shí),輸出為1,代表受到網(wǎng)絡(luò)攻擊狀態(tài)。由于直流微電網(wǎng)負(fù)荷波動(dòng)會(huì)導(dǎo)致系統(tǒng)運(yùn)行狀態(tài)發(fā)生改變,導(dǎo)致CNN-LSTM檢測(cè)模型將其錯(cuò)誤判別為異常值,需要繼續(xù)使用MIC校驗(yàn)器檢驗(yàn)是否存在FDIA,此時(shí)若存在FDIA,線路電流和母線電壓誤差之間的相關(guān)度會(huì)大幅降低,這是因?yàn)镕DIA篡改了信息傳輸過(guò)程中的電壓或者電流,但由于數(shù)據(jù)不是自然變化,此時(shí)二者之間的關(guān)系受到破壞,導(dǎo)致系統(tǒng)控制紊亂,MIC對(duì)于此類(lèi)數(shù)據(jù)變化極為敏感。

圖4 基于CNN-LSTM聯(lián)合MIC二階段FDIA檢測(cè)機(jī)制

Fig.4 Two-stage FDIA detection mechanism based on CNN-LSTM combined with MIC

  2.4 FDIA檢測(cè)與修正

  針對(duì)直流微電網(wǎng)受到FDIA時(shí),基于神經(jīng)網(wǎng)絡(luò)模型,將注入的虛假數(shù)據(jù)值λj剔除,使用正確的傳輸值uj(t)作為二次控制信號(hào)的輸入。基于以上策略,本文首先通過(guò)神經(jīng)網(wǎng)絡(luò)估計(jì)ui(t)的值,并將所得到的估計(jì)值作為二次控制收到信號(hào)的參考值。修正模型由神經(jīng)網(wǎng)絡(luò)、PI控制器和輸出組成。圖5展現(xiàn)了本文的神經(jīng)網(wǎng)絡(luò)檢測(cè)模型的構(gòu)架。在圖5中,使用神經(jīng)網(wǎng)絡(luò)來(lái)估計(jì)uj的實(shí)際值,使用PI控制器將輸入收斂到神經(jīng)網(wǎng)絡(luò)的輸出。若PI控制器的輸出命名為μ(t),則PI控制器的輸入為

圖5 虛假數(shù)據(jù)修正模型

Fig.5 False data correction model

  采用神經(jīng)網(wǎng)絡(luò)對(duì)直流母線的電壓值進(jìn)行估計(jì)。當(dāng)檢測(cè)模型認(rèn)為通信通道中存在虛假數(shù)據(jù),為了消除虛假數(shù)據(jù)的影響,使用神經(jīng)網(wǎng)絡(luò)的估計(jì)值代替Ui,直到檢測(cè)模型確認(rèn)虛假數(shù)據(jù)已被系統(tǒng)剔除,微電網(wǎng)恢復(fù)正常運(yùn)行?;谒岢龅奶摷贁?shù)據(jù)檢測(cè)模型,作為估計(jì)器的神經(jīng)網(wǎng)絡(luò)模型有重要的作用,需要估計(jì)值與系統(tǒng)當(dāng)前運(yùn)行狀態(tài)下的Ui誤差很小。

  3 算例分析

  基于Matlab/Simulink平臺(tái)構(gòu)建了圖1所示的直流微電網(wǎng)模型,模型由4個(gè)DG組成,每個(gè)DG都為一個(gè)300 V的直流電壓源,直流母線的參考電壓為370 V,模型參數(shù)為R1=0.4Ω,R2=0.6Ω,R3=0.3Ω,R4=0.5Ω。

  本文在模擬微電網(wǎng)遭遇FDIA時(shí),假設(shè)攻擊者可以操控傳輸數(shù)據(jù)對(duì)微電網(wǎng)發(fā)動(dòng)攻擊。本文考慮3種測(cè)試場(chǎng)景。

  3.1 場(chǎng)景1:正常運(yùn)行狀態(tài)

  在場(chǎng)景1中,改變直流微電網(wǎng)的負(fù)載,通過(guò)比較直流微電網(wǎng)真實(shí)運(yùn)行狀態(tài)和檢測(cè)模型估計(jì)值,驗(yàn)證負(fù)荷不斷變化情況下,模型預(yù)測(cè)的準(zhǔn)確性。直流微電網(wǎng)在正常狀態(tài)下運(yùn)行,在t=0s時(shí)處于穩(wěn)態(tài),在t=1.5s時(shí)接入負(fù)載,仿真結(jié)果如圖6所示。由圖6可以看出,在直流微電網(wǎng)正常運(yùn)行時(shí),預(yù)測(cè)值和真實(shí)值之間的誤差為0,當(dāng)微電網(wǎng)中存在負(fù)荷變化時(shí),出現(xiàn)了極小的估計(jì)誤差,此時(shí)估計(jì)誤差的值小于0.5且經(jīng)過(guò)短暫的時(shí)間后,誤差重新收斂為0,結(jié)果表明,文中所提出的深度學(xué)習(xí)檢測(cè)模型能夠準(zhǔn)確預(yù)測(cè)出直流微電網(wǎng)正常運(yùn)行時(shí)的電流、電壓值。

圖6 添加負(fù)載時(shí)的直流微電網(wǎng)電流電壓值

Fig.6 The current and voltage value of DC microgrid with load added

  3.2 場(chǎng)景2:虛假數(shù)據(jù)注入攻擊

  在場(chǎng)景2中,對(duì)直流微電網(wǎng)信息傳輸通道加入固定值的虛假數(shù)據(jù)和時(shí)變的虛假數(shù)據(jù),對(duì)比加入檢測(cè)模型的直流微電網(wǎng)和不加入模型的直流微電網(wǎng)驗(yàn)證檢測(cè)模型的有效性。直流微電網(wǎng)在正常的狀態(tài)下運(yùn)行,在t=0s時(shí)處于穩(wěn)態(tài),在t=1s時(shí)對(duì)直流微電網(wǎng)DG2注入λ=?20的虛假數(shù)據(jù),在t=2s時(shí)在第一次攻擊的基礎(chǔ)上在DG4注入攻擊為λ=20sin?2πt,結(jié)果如圖7所示。由圖7 a)可以看出,在加入虛假數(shù)據(jù)后,母線電壓在二次控制的作用下發(fā)生了改變,但是在很短的時(shí)間內(nèi)恢復(fù)平穩(wěn)運(yùn)行。由圖7 b)可以看出,t=2s時(shí),只有DG2受到攻擊時(shí),可以消除攻擊信號(hào)對(duì)系統(tǒng)的不利影響,但是當(dāng)攻擊信號(hào)的幅度和數(shù)量增加時(shí),這些攻擊信號(hào)不能被完全抑制,導(dǎo)致系統(tǒng)中出現(xiàn)較大的波動(dòng)。圖7 c)中的母線電壓和電流由于虛假數(shù)據(jù)的存在,控制系統(tǒng)受到虛假數(shù)據(jù)的影響最終導(dǎo)致失控。結(jié)果表明當(dāng)直流微電網(wǎng)受攻擊后,檢測(cè)模型能夠檢測(cè)出虛假數(shù)據(jù)的存在且能夠快速響應(yīng)消除存在的虛假數(shù)據(jù)。

圖7 加入虛假數(shù)據(jù)后的直流微電網(wǎng)電壓電流值

Fig.7 The voltage and current value of DC microgrid after adding false data

  3.3 場(chǎng)景3:負(fù)荷變化和虛假數(shù)據(jù)注入攻擊

  場(chǎng)景3是為了證明在負(fù)荷不斷變化的情況下,注入虛假數(shù)據(jù)后,檢測(cè)模型能夠快速準(zhǔn)確檢測(cè)出虛假數(shù)據(jù)并排除虛假數(shù)據(jù)的影響。系統(tǒng)開(kāi)始時(shí)穩(wěn)定運(yùn)行,在t=1s時(shí)添加負(fù)荷,在t=2.5s時(shí)對(duì)DG3傳輸通道加入λ=15的虛假數(shù)據(jù),在t=4s時(shí)將虛假數(shù)據(jù)改為λ=10sin?2πt?10并減少負(fù)荷,結(jié)果如表1和圖8所示。

表1 攻擊前后的MIC值

Table 1 The MIC value before and after the attack

圖8 直流微電網(wǎng)電壓電流值

Fig.8 The voltage and current value of DC microgrid

  如表1所示,當(dāng)系統(tǒng)正常運(yùn)行和負(fù)荷變化時(shí),MIC值在0.8~1的區(qū)間內(nèi),當(dāng)系統(tǒng)中存在虛假數(shù)據(jù)時(shí),MIC的值低于0.5。可以看出負(fù)荷變化對(duì)于MIC校驗(yàn)器影響很小,而存在虛假數(shù)據(jù)時(shí),變量間的關(guān)系受到破壞,MIC校驗(yàn)器對(duì)于變量間微小的變化極為敏感。

  由圖8可以看出,當(dāng)發(fā)生負(fù)荷變化時(shí),直流微電網(wǎng)正常運(yùn)行,檢測(cè)模型不動(dòng)作;當(dāng)注入了虛假數(shù)據(jù)后,檢測(cè)模型可以迅速響應(yīng),使微電網(wǎng)運(yùn)行不受到虛假數(shù)據(jù)的影響。

  3.4 MIC校驗(yàn)性能分析

  使用MIC方法分別測(cè)試在負(fù)荷變化時(shí)和虛假數(shù)據(jù)注入攻擊時(shí),線路電流和母線電壓之間的相關(guān)性,通過(guò)實(shí)驗(yàn)證明MIC模型的鑒別能力。

  首先將負(fù)荷設(shè)置為不斷發(fā)生改變的狀態(tài),在t=1s時(shí)對(duì)DG3傳輸通道加入λ=10的虛假數(shù)據(jù)。DG3受到攻擊前后1 s的電流值和母線電壓誤差值如圖9所示。圖9 a)中通過(guò)式(15)可以計(jì)算得到MIC的值為0.85,圖9 b)中通過(guò)式(15)可以計(jì)算得到MIC的值為0.08。可見(jiàn)MIC模型對(duì)于系統(tǒng)中存在的虛假數(shù)據(jù)有著很高的靈敏度,由此可得,當(dāng)系統(tǒng)中存在虛假數(shù)據(jù)時(shí),線路電流和母線電壓誤差之間的相關(guān)性會(huì)受到破壞。

圖9 攻擊前后線路電流和母線電壓誤差

Fig.9 Line current and bus voltage errors before and after the attack

  3.5 檢測(cè)效果對(duì)比分析

  使用4個(gè)指標(biāo)來(lái)評(píng)價(jià)檢測(cè)模型的效果:準(zhǔn)確率A表示被正確分類(lèi)的比例;精確率P表示被預(yù)測(cè)為異常的樣本中預(yù)測(cè)正確的比例;召回率R表示異常樣本中預(yù)測(cè)正確的比例;F值表示精確率與召回率的調(diào)和平均值。

  式中:σTP、σTN分別表示正確分類(lèi)的正、負(fù)樣本;σFN表示錯(cuò)誤分類(lèi)為負(fù)的正樣本;σFP表示錯(cuò)誤分類(lèi)為正的負(fù)樣本。

  本文使用文獻(xiàn)[23]中的虛假數(shù)據(jù)生成方法,模擬了直流微電網(wǎng)不同類(lèi)型的虛假數(shù)據(jù)注入攻擊和正常運(yùn)行狀態(tài),使用自編碼器生成類(lèi)似的場(chǎng)景,通過(guò)Bootstrap模型分類(lèi)不同類(lèi)型的樣本集,并從樣本集中隨機(jī)采樣實(shí)例來(lái)完成本文的檢測(cè)實(shí)驗(yàn)。

  為了清晰反映本文所使用的檢測(cè)模型的性能,測(cè)試集中存在1343個(gè)虛假數(shù)據(jù)樣本和1078個(gè)正常運(yùn)行樣本,根據(jù)評(píng)價(jià)指標(biāo),不同模型的結(jié)果如表2所示。可以看出本文所使用的檢測(cè)模型對(duì)于正常樣本的誤報(bào)率最低,存在虛假數(shù)據(jù)時(shí),檢測(cè)率更高。

表2 不同模型的檢測(cè)結(jié)果

Table 2 Test results of different models

  4 結(jié)語(yǔ)

  直流微電網(wǎng)使用信息網(wǎng)絡(luò)傳遞控制信號(hào)時(shí)容易遭受網(wǎng)絡(luò)攻擊,為此,本文提出了一種檢測(cè)和修正直流微電網(wǎng)FDIA的防御策略。該方法利用CNN-LSTM網(wǎng)絡(luò)估計(jì)直流微電網(wǎng)中所有發(fā)電單元的直流電壓和直流輸出電流。當(dāng)直流微電網(wǎng)中不存在FDIA時(shí),檢測(cè)模型不動(dòng)作,但當(dāng)攻擊者試圖向直流微電網(wǎng)系統(tǒng)中注入虛假數(shù)據(jù)時(shí),估計(jì)器所得到的估計(jì)值和系統(tǒng)運(yùn)行的真實(shí)值會(huì)產(chǎn)生較大的偏差,此時(shí)判定系統(tǒng)中可能存在虛假數(shù)據(jù),經(jīng)過(guò)MIC校驗(yàn)器確認(rèn)FDIA的存在。本文通過(guò)離線數(shù)字時(shí)域仿真有效地驗(yàn)證了所提出的檢測(cè)方法可以成功識(shí)別出FDIA且能恢復(fù)系統(tǒng)的真實(shí)運(yùn)行狀態(tài)。本文的方法可以區(qū)分微電網(wǎng)正常運(yùn)行下的負(fù)荷變化和FDIA,最大限度降低FDIA的影響。

  本文所設(shè)計(jì)的基于直流微電網(wǎng)的虛假數(shù)據(jù)檢測(cè)模型能夠有效抵御直流電網(wǎng)中FDIA,但研究?jī)?nèi)容尚處于起步階段,應(yīng)用于未來(lái)更復(fù)雜的直流微電網(wǎng)中還需要更深入的研究。

  注:本文內(nèi)容呈現(xiàn)略有調(diào)整,如需要請(qǐng)查看原文。


分享到:

關(guān)鍵字:微電網(wǎng)

中國(guó)儲(chǔ)能網(wǎng)版權(quán)說(shuō)明:

1、凡注明來(lái)源為“中國(guó)儲(chǔ)能網(wǎng):xxx(署名)”,除與中國(guó)儲(chǔ)能網(wǎng)簽署內(nèi)容授權(quán)協(xié)議的網(wǎng)站外,未經(jīng)本網(wǎng)授權(quán),任何單位及個(gè)人不得轉(zhuǎn)載、摘編或以其它方式使用上述作品。

2、凡本網(wǎng)注明“來(lái)源:xxx(非中國(guó)儲(chǔ)能網(wǎng))”的作品,均轉(zhuǎn)載與其他媒體,目的在于傳播更多信息,但并不代表中國(guó)儲(chǔ)能網(wǎng)贊同其觀點(diǎn)、立場(chǎng)或證實(shí)其描述。其他媒體如需轉(zhuǎn)載,請(qǐng)與稿件來(lái)源方聯(lián)系,如產(chǎn)生任何版權(quán)問(wèn)題與本網(wǎng)無(wú)關(guān)。

3、如因作品內(nèi)容、版權(quán)以及引用的圖片(或配圖)內(nèi)容僅供參考,如有涉及版權(quán)問(wèn)題,可聯(lián)系我們直接刪除處理。請(qǐng)?jiān)?0日內(nèi)進(jìn)行。

4、有關(guān)作品版權(quán)事宜請(qǐng)聯(lián)系:13661266197、 郵箱:[email protected]

衡东县| 大洼县| 醴陵市| 铜陵市| 视频| 泰安市| 腾冲县| 隆德县| 高台县| 金湖县| 临沭县| 威远县| 上蔡县| 古田县| 渭源县| 乌兰察布市| 锡林浩特市| 通江县| 南和县| 保靖县| 青岛市| 于田县| 旬阳县| 石景山区| 乌海市| 凤翔县| 喀喇| 湟源县| 平凉市| 普定县| 桃园市| 泊头市| 桐柏县| 吉林市| 建湖县| 宿州市| 独山县| 甘泉县| 邛崃市| 年辖:市辖区| 上虞市|