两个人的电影免费视频_国产精品久久久久久久久成人_97视频在线观看播放_久久这里只有精品777_亚洲熟女少妇二三区_4438x8成人网亚洲av_内谢国产内射夫妻免费视频_人妻精品久久久久中国字幕

一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法與流程

文檔序號(hào):11155503閱讀:640來源:國(guó)知局
一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法與制造工藝

本發(fā)明涉及輸變電設(shè)備檢測(cè)領(lǐng)域,具體涉及一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法。



背景技術(shù):

眾所周知,設(shè)備的運(yùn)行狀態(tài)對(duì)于電廠生產(chǎn)影響巨大,設(shè)備在正常工況運(yùn)行過程中體現(xiàn)出來的功能特征、表象特征以及電氣特征與異常狀態(tài)下的特征是有差異的。根據(jù)設(shè)備狀態(tài)預(yù)警,做出檢修計(jì)劃,將節(jié)省電廠巨大的人力、物力、財(cái)力,具有非常大的經(jīng)濟(jì)潛力。

目前,在設(shè)備狀態(tài)預(yù)警方面,設(shè)備生產(chǎn)廠商根據(jù)制造工藝設(shè)定各個(gè)部件不同預(yù)警級(jí)別的閾值方法,現(xiàn)有設(shè)備異常檢測(cè)方法采用靜態(tài)閾值設(shè)置的方法,即以現(xiàn)有數(shù)據(jù)平均值為基準(zhǔn),上下分別以均值的固定百分比(例如150%、300%)作為上下閾值,進(jìn)行數(shù)據(jù)異常檢測(cè)。

上下閾值之間的區(qū)域?yàn)椤罢!眳^(qū)域,上下閾值之外的區(qū)域?yàn)椤爱惓!眳^(qū)域。靜態(tài)閾值設(shè)置無法反映數(shù)據(jù)的動(dòng)態(tài)實(shí)時(shí)特點(diǎn)。閾值上下界設(shè)置采用業(yè)務(wù)專家的主觀經(jīng)驗(yàn),無法反映數(shù)據(jù)客觀規(guī)律。閾值的確定直接影響了對(duì)于變電設(shè)備的檢測(cè)精度以及數(shù)據(jù)處理的效率,然而目前靜態(tài)閾值的方式已經(jīng)無法滿足日趨復(fù)雜的變電設(shè)備檢測(cè)。



技術(shù)實(shí)現(xiàn)要素:

本發(fā)明的目的在于克服現(xiàn)有技術(shù)的不足,提供一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法,解決變電設(shè)備狀態(tài)參數(shù)異常檢測(cè)閾值區(qū)間較多依賴變電設(shè)備專家離線計(jì)算及主觀經(jīng)驗(yàn),而無法根據(jù)設(shè)備運(yùn)行實(shí)時(shí)狀態(tài)自適應(yīng)方式變化的問題,采用了多種歸一化方式進(jìn)行數(shù)據(jù)處理,并且對(duì)采集的實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)首先進(jìn)行多重的去野值判斷,提高了數(shù)據(jù)的有效性和處理效率。

本發(fā)明提供了一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法,依次包括如下步驟:

步驟1:在變電系統(tǒng)數(shù)據(jù)庫里,找到變電設(shè)備的所有歷史指標(biāo)信息,選取一段時(shí)間內(nèi)變電設(shè)備在正常運(yùn)行狀態(tài)下和非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

步驟2:

A.對(duì)選取的變電設(shè)備的正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)按照如下公式進(jìn)行歸一化處理:

其中,hn1為h1歸一化后的值,hmax1和hmin1分別為變電設(shè)備正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)的最大值和最小值,h1為變電設(shè)備正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

B.對(duì)選取的變電設(shè)備的非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)按照如下公式進(jìn)行歸一化處理:

其中,hn2為h2歸一化后的值,hmax2和hmin2分別為變電設(shè)備非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)的最大值和最小值,h2為變電設(shè)備非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

步驟3:

A:將歸一化后的hn1做平均值計(jì)算,得到平均值H1,并且設(shè)置閾值補(bǔ)償參數(shù)L1,設(shè)置正常狀態(tài)下初始野值濾除閾值的下限值為H1-L1,上限值為H1+L1,;

B:將歸一化后的hn2做平均值計(jì)算,得到平均值H2,計(jì)算ΔH=|H1-H2|,并且設(shè)置誤差閾值補(bǔ)償參數(shù)L2,設(shè)置非正常狀態(tài)下初始野值濾除誤差閾值的下限值為ΔH-L2,上限值為ΔH+L2;

步驟4:

A.將正常狀態(tài)下初始野值濾除閾值的下限值為H-L,上限值為H+L分別利用如下公式進(jìn)行還原,得到還原后的初始野值濾除閾值的上、下限值:

B.將非正常狀態(tài)下初始野值濾除誤差閾值的下限值為ΔH-L2,上限值為ΔH+L2分別利用如下公式進(jìn)行還原,得到還原后的初始野值濾除誤差閾值的上、下限值:

步驟5:實(shí)時(shí)獲取變電設(shè)備運(yùn)行狀態(tài)多維度實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),判斷:

A:比較實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)和初始野值濾除閾值的上、下限值的差值,如果差值沒有落入初始野值濾除閾值的上、下限值構(gòu)成的閾值范圍,則剔除此實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),如果都落入則進(jìn)入下一步驟;

B:判斷實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)是否落入初始野值濾除閾值的上、下限值構(gòu)成的閾值范圍,如果未落入則剔除此實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),如果落入則進(jìn)行歸一化處理,獲得變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列;

步驟6:通過對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列應(yīng)用降維操作,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列;

步驟7:應(yīng)用可配置的預(yù)測(cè)器,對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列進(jìn)行預(yù)測(cè)計(jì)算,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的預(yù)測(cè)數(shù)據(jù)序列;

步驟8:應(yīng)用變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)數(shù)據(jù)及預(yù)測(cè)數(shù)據(jù)進(jìn)行分值計(jì)算,得到變電設(shè)備運(yùn)行狀態(tài)分值;

步驟9:確定變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,包括上閾值界、下閾值界;

步驟10:應(yīng)用變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,判斷當(dāng)前時(shí)刻分值是否在閾值區(qū)間范圍,判斷變電設(shè)備的狀態(tài)異常。

其中,步驟5中進(jìn)行歸一化處理,獲得變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列具體為:

多維度實(shí)時(shí)變電設(shè)備參數(shù)xt定義為t時(shí)刻變電設(shè)備數(shù)據(jù)集合,x(i)定義為多維度變電設(shè)備數(shù)據(jù)的第i維,i=1,2,…,n,表示為:

xt=(x(1),x(2),…,x(i),…,x(n))

對(duì)多維度實(shí)時(shí)變電設(shè)備參數(shù)xt進(jìn)行正規(guī)化,實(shí)現(xiàn)歸一化運(yùn)算,

得到t時(shí)刻變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列到nt,其中NOR(xt)是正規(guī)化運(yùn)算,||xt||,是xt的模。

其中,步驟6具體為對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列nt進(jìn)行降維運(yùn)算:

yt=PCA(nt)

得到t時(shí)刻變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列yt,其中PCA(xt)是降維運(yùn)算,通過對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列應(yīng)用降維操作,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列。

其中,步驟7具體為:應(yīng)用變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列yt和變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)數(shù)據(jù)及預(yù)測(cè)數(shù)據(jù)進(jìn)行分值計(jì)算,

得到變電設(shè)備運(yùn)行狀態(tài)分值到zt。

其中,步驟9具體為:

[Q1-1.5IQR,Q3+1.5IQR]

其中,Q1-1.5IQR為上閾值;Q3+1.5IQR為下閾值;四分位數(shù)間距IQR=Q3-Q1。

其中,步驟10具體為:使用變電設(shè)備運(yùn)行狀態(tài)分值到zt與變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間[Q1-1.5IQR,Q3+1.5IQR]進(jìn)行比較,如果zt在該區(qū)間范圍,則變電設(shè)備狀態(tài)正常;否則,變電設(shè)備狀態(tài)異常。

本發(fā)明的多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法,可以實(shí)現(xiàn):

解決變電設(shè)備狀態(tài)參數(shù)異常檢測(cè)閾值區(qū)間較多依賴變電設(shè)備專家離線計(jì)算及主觀經(jīng)驗(yàn),而無法根據(jù)設(shè)備運(yùn)行實(shí)時(shí)狀態(tài)自適應(yīng)方式變化的問題,利用特定順序和特殊的處理方式,并且采用了多種歸一化方式進(jìn)行數(shù)據(jù)處理,并且對(duì)采集的實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)首先進(jìn)行多重的去野值判斷,提高了數(shù)據(jù)的有效性和處理效率,提高了運(yùn)算效率,并且能夠?qū)崟r(shí)快速、自動(dòng)化的檢測(cè)設(shè)備異常情況。

附圖說明

圖1為野值濾除閾值處理流程圖

圖2為變電設(shè)備參數(shù)異常檢測(cè)方法流程圖

具體實(shí)施方式

下面詳細(xì)說明本發(fā)明的具體實(shí)施,有必要在此指出的是,以下實(shí)施只是用于本發(fā)明的進(jìn)一步說明,不能理解為對(duì)本發(fā)明保護(hù)范圍的限制,該領(lǐng)域技術(shù)熟練人員根據(jù)上述本發(fā)明內(nèi)容對(duì)本發(fā)明做出的一些非本質(zhì)的改進(jìn)和調(diào)整,仍然屬于本發(fā)明的保護(hù)范圍。

本發(fā)明提供了一種多維度實(shí)時(shí)變電設(shè)備數(shù)據(jù)的閾值自適應(yīng)設(shè)置異常檢測(cè)方法,如圖1和2所示,具體包括如下的步驟:

步驟1:在變電系統(tǒng)數(shù)據(jù)庫里,找到變電設(shè)備的所有歷史指標(biāo)信息,選取一段時(shí)間內(nèi)變電設(shè)備在正常運(yùn)行狀態(tài)下和非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

步驟2:

A.對(duì)選取的變電設(shè)備的正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)按照如下公式進(jìn)行歸一化處理:

其中,hn1為h1歸一化后的值,hmax1和hmin1分別為變電設(shè)備正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)的最大值和最小值,h1為變電設(shè)備正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

B.對(duì)選取的變電設(shè)備的非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)按照如下公式進(jìn)行歸一化處理:

其中,hn2為h2歸一化后的值,hmax2和hmin2分別為變電設(shè)備非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù)的最大值和最小值,h2為變電設(shè)備非正常運(yùn)行狀態(tài)下的歷史數(shù)據(jù);

步驟3:

A:將歸一化后的hn1做平均值計(jì)算,得到平均值H1,并且設(shè)置閾值補(bǔ)償參數(shù)L1,設(shè)置正常狀態(tài)下初始野值濾除閾值的下限值為H1-L1,上限值為H1+L1,其中L1可以為根據(jù)實(shí)際需要濾除野值而設(shè)置的,也可以根據(jù)長(zhǎng)時(shí)間變電設(shè)備的歷史數(shù)據(jù)做出的經(jīng)驗(yàn)選擇,例如可以是0.1,0.2等;

B:將歸一化后的hn2做平均值計(jì)算,得到平均值H2,計(jì)算ΔH=|H1-H2|,并且設(shè)置誤差閾值補(bǔ)償參數(shù)L2,設(shè)置非正常狀態(tài)下初始野值濾除誤差閾值的下限值為ΔH-L2,上限值為ΔH+L2,其中L2可以為根據(jù)實(shí)際要濾除野值而設(shè)置的,也可以根據(jù)長(zhǎng)時(shí)間變電設(shè)備的歷史數(shù)據(jù)做出的經(jīng)驗(yàn)選擇,例如可以是0.03,0.05等;

步驟4:

A.將正常狀態(tài)下初始野值濾除閾值的下限值為H-L,上限值為H+L分別利用如下公式進(jìn)行還原,得到還原后的初始野值濾除閾值的上、下限值:

B.將非正常狀態(tài)下初始野值濾除誤差閾值的下限值為ΔH-L2,上限值為ΔH+L2分別利用如下公式進(jìn)行還原,得到還原后的初始野值濾除誤差閾值的上、下限值:

步驟5:實(shí)時(shí)獲取變電設(shè)備運(yùn)行狀態(tài)多維度實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),判斷:

A:比較實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)和初始野值濾除閾值的上、下限值的差值,如果差值沒有落入初始野值濾除閾值的上、下限值構(gòu)成的閾值范圍,則剔除此實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),如果都落入則進(jìn)入下一步驟;

B:判斷實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù)是否落入初始野值濾除閾值的上、下限值構(gòu)成的閾值范圍,如果未落入則剔除此實(shí)時(shí)測(cè)點(diǎn)數(shù)據(jù),如果落入則進(jìn)行歸一化處理,獲得變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列;

進(jìn)行歸一化處理,獲得變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列;

其中對(duì)測(cè)點(diǎn)數(shù)據(jù)進(jìn)行歸一化運(yùn)算,具體為;

多維度實(shí)時(shí)變電設(shè)備參數(shù)xt定義為t時(shí)刻變電設(shè)備數(shù)據(jù)集合,x(i)定義為多維度變電設(shè)備數(shù)據(jù)的第i維,i=1,2,…,n,表示為:

xt=(x(1),x(2),…,x(i),…,x(n))

所述歸一化處理,是指對(duì)多維度實(shí)時(shí)變電設(shè)備參數(shù)xt進(jìn)行正規(guī)化,實(shí)現(xiàn)歸一化運(yùn)算,

得到t時(shí)刻變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列到nt,其中NOR(xt)是正規(guī)化運(yùn)算,||xt||,是xt的模。

其中所述變電設(shè)備包括(但不限于)開關(guān)刀閘、變壓器等,所述變電設(shè)備狀態(tài)參數(shù)包括(但不限于)電流、電壓、功率、功率因數(shù)、油溫、油中氣體數(shù)據(jù)等。

步驟6:通過對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列應(yīng)用降維操作,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列;

對(duì)歸一化處理后的數(shù)據(jù)序列進(jìn)行多維度變電設(shè)備數(shù)據(jù)降維,具體為:

所述進(jìn)行降維操作,是指對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列nt進(jìn)行降維運(yùn)算:

yt=PCA(nt)

得到t時(shí)刻變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列yt,其中PCA(xt)是降維運(yùn)算,通過對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)多維度歸一化數(shù)據(jù)序列應(yīng)用降維操作,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列;所述降維運(yùn)算PCA(xt)包括(但不限于)離散余弦變換、主成分分析、離散小波變換等。

步驟7:應(yīng)用可配置的預(yù)測(cè)器,對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列進(jìn)行預(yù)測(cè)計(jì)算,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的預(yù)測(cè)數(shù)據(jù)序列;

對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列進(jìn)行預(yù)測(cè)計(jì)算,產(chǎn)生變電設(shè)備運(yùn)行狀態(tài)的預(yù)測(cè)數(shù)據(jù)序列,具體為:

所述預(yù)測(cè)計(jì)算,是指對(duì)變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列yt進(jìn)行預(yù)測(cè)計(jì)算,

得到t時(shí)刻變電設(shè)備運(yùn)行狀態(tài)的預(yù)測(cè)數(shù)據(jù)序列其中AR(yt)是預(yù)測(cè)器;其中所述預(yù)測(cè)器包括(但不限于):多項(xiàng)式回歸預(yù)測(cè)器、支持向量機(jī)預(yù)測(cè)器、人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)器、卡爾曼預(yù)測(cè)濾波器、決策樹回歸預(yù)測(cè)器。

步驟8:應(yīng)用變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)數(shù)據(jù)及預(yù)測(cè)數(shù)據(jù)進(jìn)行分值計(jì)算,得到變電設(shè)備運(yùn)行狀態(tài)分值;

利用變電設(shè)備實(shí)際數(shù)據(jù)與預(yù)測(cè)數(shù)據(jù)進(jìn)行分值計(jì)算,得到當(dāng)前時(shí)刻電設(shè)備運(yùn)行狀態(tài)分值,具體為:

所述分值計(jì)算,是指應(yīng)用變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)降維數(shù)據(jù)序列yt和變電設(shè)備運(yùn)行狀態(tài)的實(shí)時(shí)數(shù)據(jù)及預(yù)測(cè)數(shù)據(jù)進(jìn)行分值計(jì)算,

得到變電設(shè)備運(yùn)行狀態(tài)分值到zt;

步驟9:確定變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,包括上閾值界、下閾值界;

確定變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,包括上閾值界、下閾值界,計(jì)算正常閾值范圍,具體為:

所述變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,是指

[Q1-1.5IQR,Q3+1.5IQR]

其中,Q1-1.5IQR為上閾值;Q3+1.5IQR為下閾值;四分位數(shù)間距IQR=Q3-Q1。

步驟10:應(yīng)用變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間,判斷當(dāng)前時(shí)刻分值是否在閾值區(qū)間范圍,并做出變電設(shè)備運(yùn)行狀態(tài)是否異常的結(jié)論。

所述判斷當(dāng)前時(shí)刻分值是否在閾值區(qū)間范圍,并做出變電設(shè)備運(yùn)行狀態(tài)是否異常的結(jié)論是指使用變電設(shè)備運(yùn)行狀態(tài)分值到zt與變電設(shè)備運(yùn)行狀態(tài)閾值區(qū)間[Q1-1.5IQR,Q3+1.5IQR]進(jìn)行比較,如果zt在該區(qū)間范圍,則變電設(shè)備狀態(tài)正常;否則,變電設(shè)備狀態(tài)異常。

盡管為了說明的目的,已描述了本發(fā)明的示例性實(shí)施方式,但是本領(lǐng)域的技術(shù)人員將理解,不脫離所附權(quán)利要求中公開的發(fā)明的范圍和精神的情況下,可以在形式和細(xì)節(jié)上進(jìn)行各種修改、添加和替換等的改變,而所有這些改變都應(yīng)屬于本發(fā)明所附權(quán)利要求的保護(hù)范圍,并且本發(fā)明要求保護(hù)的產(chǎn)品各個(gè)部門和方法中的各個(gè)步驟,可以以任意組合的形式組合在一起。因此,對(duì)本發(fā)明中所公開的實(shí)施方式的描述并非為了限制本發(fā)明的范圍,而是用于描述本發(fā)明。相應(yīng)地,本發(fā)明的范圍不受以上實(shí)施方式的限制,而是由權(quán)利要求或其等同物進(jìn)行限定。

當(dāng)前第1頁1 2 3 
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評(píng)論。精彩留言會(huì)獲得點(diǎn)贊!
1
宝坻区| 安顺市| 金山区| 藁城市| 安仁县| 临猗县| 西宁市| 祁东县| 东乌珠穆沁旗| 乌鲁木齐县| 金沙县| 吉水县| 谢通门县| 永春县| 萨嘎县| 乐陵市| 宝清县| 贵定县| 新河县| 滁州市| 清涧县| 固镇县| 浠水县| 旬阳县| 怀来县| 剑阁县| 神池县| 北辰区| 台东市| 罗源县| 阜阳市| 鄱阳县| 崇义县| 西乌珠穆沁旗| 阜阳市| 漳平市| 长顺县| 怀化市| 昌平区| 泗阳县| 临邑县|