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

?

多層螺旋霧幕技術(shù)的噴霧模擬及降塵試驗(yàn)研究*

2019-09-06 07:00荊德吉葛少成孟祥曦
關(guān)鍵詞:水霧除塵粉塵

荊德吉,徐 放,葛少成,張 天,孟祥曦

(1.遼寧工程技術(shù)大學(xué) 安全科學(xué)與工程學(xué)院,遼寧 阜新 123000; 2.遼寧工程技術(shù)大學(xué) 安全科學(xué)與工程研究院,遼寧 阜新 123000; 3.礦山熱動(dòng)力災(zāi)害與防治教育部重點(diǎn)實(shí)驗(yàn)室(遼寧工程技術(shù)大學(xué)),遼寧 阜新 123000; 4.太原理工大學(xué) 安全與應(yīng)急管理工程學(xué)院,山西 太原 030024)

0 引言

掘進(jìn)工作面是煤礦的主要產(chǎn)塵區(qū)域,隨著采掘技術(shù)的進(jìn)步,礦井掘進(jìn)的效率逐步提高,然而高效率的掘進(jìn)作業(yè)隨之帶來(lái)的是產(chǎn)塵量的增加,導(dǎo)致更嚴(yán)重的粉塵污染[1-3]。

目前對(duì)于煤礦的開(kāi)放性塵源所采用的除塵方法多為噴霧除塵[4],然而傳統(tǒng)的內(nèi)外噴霧有很多弊端,如水霧的擴(kuò)散速度較慢,水霧范圍較小,不能有效阻止粉塵的擴(kuò)散,導(dǎo)致除塵效率不高[5-6]。多項(xiàng)模擬和實(shí)驗(yàn)表明,利用霧幕相比傳統(tǒng)噴霧控塵效果更好,且風(fēng)幕對(duì)集塵與除塵也有良好的效果。劉寶明等對(duì)掘進(jìn)機(jī)旋轉(zhuǎn)噴霧降塵進(jìn)行了應(yīng)用研究[7];張建卓等基于Fluent模擬論證提出了掘進(jìn)面集塵風(fēng)幕除塵方案[8];李雨成等設(shè)計(jì)風(fēng)幕控塵裝置并通過(guò)試驗(yàn)驗(yàn)證裝置的優(yōu)越性[2];Nie等采用數(shù)值模擬與現(xiàn)場(chǎng)測(cè)試的方法對(duì)旋轉(zhuǎn)風(fēng)幕對(duì)掘進(jìn)工作面的抑塵效果展開(kāi)研究[9];錢(qián)杰等在現(xiàn)有旋轉(zhuǎn)風(fēng)幕除塵的基礎(chǔ)上加入了噴霧,對(duì)旋轉(zhuǎn)水霧除塵進(jìn)行了研究[10]。

為了進(jìn)一步提高掘進(jìn)工作面噴霧除塵效率,以旋風(fēng)霧幕除塵為依據(jù),結(jié)合已知除塵方法進(jìn)行改進(jìn),提出了多層螺旋霧幕除塵方法。與風(fēng)機(jī)提供風(fēng)力不同,該方法以空壓機(jī)提供氣壓,通過(guò)噴嘴射出氣流。相比傳統(tǒng)的單層霧幕,噴嘴的多層螺旋排列方式也有助于形成多層且范圍更廣的螺旋霧幕。利用氣水兩相流的方式噴霧,研究這種方式下噴射出的氣流對(duì)水霧的運(yùn)動(dòng)軌跡及結(jié)構(gòu)的影響以及對(duì)除塵效果的影響。通過(guò)Comsol模擬軟件的模擬仿真得出多層霧幕裝置工作時(shí)的風(fēng)流場(chǎng)和霧滴運(yùn)動(dòng)軌跡并分析其對(duì)除塵效果的影響,為噴霧除塵試驗(yàn)提供理論依據(jù)。再通過(guò)搭建試驗(yàn)平臺(tái)來(lái)測(cè)試自然沉降、傳統(tǒng)噴霧和多層螺旋霧幕3種除塵方法的除塵效果,確定多層螺旋霧幕除塵方法的高效性。

1 多層螺旋噴霧模擬計(jì)算分析

利用Comsol軟件模擬計(jì)算多層螺旋噴霧的氣流流場(chǎng)、風(fēng)速、粒子軌跡,以此說(shuō)明該方法對(duì)氣流及水霧流向的影響,為下一步的試驗(yàn)提供保證和可靠依據(jù)[11-12]。

1.1 裝置設(shè)計(jì)與模型建立

該方法主要應(yīng)用于掘進(jìn)工作面,通過(guò)在掘進(jìn)機(jī)上安裝多層螺旋霧幕裝置形成多層霧幕來(lái)提高控塵效果。在眾多的可形成多層霧幕的裝置構(gòu)造方案中,將噴嘴呈螺旋排布的方法結(jié)構(gòu)簡(jiǎn)單容易實(shí)現(xiàn),且該構(gòu)造可使射出的水霧有向前的分速度,有利于水霧的向前推進(jìn)。如圖1所示,該裝置設(shè)計(jì)為:在掘進(jìn)機(jī)的截割臂上安裝管路支架,用以固定管路和噴嘴;供水管與供氣管一端分別與水泵、空壓機(jī)相連,另一端管體呈螺旋狀盤(pán)繞在管路支架上;支架固定的供水管與供氣管之間連接若干噴嘴,由水泵和空壓機(jī)提供噴霧動(dòng)力,由若干噴嘴噴出水霧形成多層霧幕,且有效除塵區(qū)域可隨著掘進(jìn)機(jī)頭位置的變化而變化,進(jìn)一步提高掘進(jìn)工作面的除塵效果。

1.噴嘴;2.供水管;3.供氣管;4.水泵;5.空壓機(jī)。圖1 多層霧幕除塵裝置設(shè)計(jì)示意Fig.1 Schematic diagram of multi-layer fog curtain dust removal device

利用Comsol數(shù)值模擬研究該裝置應(yīng)用在掘進(jìn)工作面時(shí)的工作效果,分析該方法的特點(diǎn)及可行性。以2層霧幕為例,如圖2所示,建立寬度和高度均為4 m的巷道模型,以巷道中心軸為準(zhǔn)建立10只以2層螺旋排布的空氣霧化噴嘴,噴嘴的氣流出口與噴霧出口的方向均為沿螺旋線向前的切線方向,螺旋的外環(huán)半徑為0.15 m,軸向截距為0.35 m,徑向截距為0.25 m。

圖2 模型建立及網(wǎng)格劃分Fig.2 Model building and meshing

1.2 湍流模型

風(fēng)流場(chǎng)和風(fēng)速的計(jì)算采用標(biāo)準(zhǔn)k-ω湍流模型,該模型對(duì)逆壓梯度流場(chǎng)的計(jì)算有較高的精度,可應(yīng)用于墻壁束縛流動(dòng),即可應(yīng)用于對(duì)掘進(jìn)巷道的氣流流向、流速的計(jì)算[13]。

(1)

設(shè)巷道內(nèi)空氣不可壓縮,湍流動(dòng)能k的輸運(yùn)方程為:

ρ(u·)k=

(2)

湍流耗散率ω的輸運(yùn)方程為:

ρ(u·)ω=·[(μ+μTσω)
ω=om

(3)

式中:ρ為密度,kg/m3;μT為湍流黏性系數(shù);u為風(fēng)速,m/s;σk和σω分別為湍流動(dòng)能k和耗散率ω的普朗系數(shù);-pl為雷諾應(yīng)力,N;F為其他外力,N。

1.3 曳力模型

水霧粒子軌跡的可利用流體曳力模型計(jì)算。只要物體在流體中存在相對(duì)速度,這個(gè)物體就會(huì)受到流體給它的作用力,即為流體曳力。

流體曳力服從Schiller-Naumann定律[14]:

(4)

式中:FD為曳力,N;τp為剪切應(yīng)力,N;CD為曳力系數(shù);Rer為雷諾數(shù);mp為液滴質(zhì)量,kg;u,v分別為風(fēng)速和液滴運(yùn)動(dòng)速度,m/s;ρp為液滴密度,kg/m3;dp為液滴直徑,μm;μ為氣體動(dòng)力黏度,Pa·s。

1.4 液滴破碎模型

K-H破碎模型用來(lái)解釋2種運(yùn)動(dòng)流體產(chǎn)生的不穩(wěn)定性,當(dāng)2種流體做平行運(yùn)動(dòng)且達(dá)到一定的相對(duì)速度,就會(huì)過(guò)渡至不穩(wěn)定狀態(tài)而破碎成更小的液滴[15]?;谠撛砜梢阅M水霧的破碎情況,計(jì)算出霧滴粒徑。

基于液體射流的線性穩(wěn)定分析,液滴破碎時(shí)半徑變化率為:

(5)

式中,rch為破碎后液滴半徑,μm;t為時(shí)間,s;rKH為破碎前液滴半徑,μm;τKH為破碎時(shí)間,s。

(6)

(7)

(8)

(9)

式中:ΛKH為不穩(wěn)定波波長(zhǎng);ΩKH為增長(zhǎng)最快的不穩(wěn)定波頻率;τKH為破碎時(shí)間;Z為液體Ohnesorge數(shù);T為泰勒數(shù);Wel為液體韋伯?dāng)?shù);Weg為氣體韋伯?dāng)?shù);Rel為液體雷諾數(shù);Urel為氣液速度差,m/s;rp為初始的液滴半徑,μm;ρ為氣體密度,kg/m3;σp為液體表面張力N/m;μp為液體動(dòng)力黏度,Pa·s;BKH為經(jīng)驗(yàn)?zāi)P蛥?shù),取值與噴孔尺寸設(shè)計(jì)和噴霧狀態(tài)有關(guān),通常取10。

破碎后形成的液滴半徑:

當(dāng)B0ΛKH≤rp時(shí),

rch=B0ΛKH

(10)

當(dāng)B0ΛKH≥rp時(shí),

(11)

式中:B0為模型經(jīng)驗(yàn)參數(shù),用來(lái)控制破碎時(shí)間。

1.5 模型條件參數(shù)設(shè)定

以下參數(shù)設(shè)置均以環(huán)境溫度為20 ℃時(shí)的工作狀態(tài)為依據(jù)設(shè)定。環(huán)境和噴霧參數(shù)設(shè)定如表1所示。

表1 參數(shù)設(shè)定Table 1 Parameter setting

1.6 模擬結(jié)果分析

計(jì)算待風(fēng)流穩(wěn)定后的穩(wěn)態(tài)結(jié)果,風(fēng)流走向如圖3所示。

圖3 風(fēng)流走向Fig.3 Wind flow direction

從圖3可以看出,在各個(gè)噴嘴射出的風(fēng)流之間相互擾動(dòng)與巷道墻壁對(duì)風(fēng)流的順向引導(dǎo)的作用下,形成了旋風(fēng),且旋風(fēng)布滿(mǎn)整個(gè)巷道。根據(jù)箭頭的走向也可直觀地看出其風(fēng)向是旋轉(zhuǎn)的。

將該巷道的速度場(chǎng)沿y軸方向分為4層切面,觀察并分析各層之間風(fēng)速分布的規(guī)律(見(jiàn)圖4)。

圖4 風(fēng)速分布Fig.4 Wind speed distribution map

由圖4可以看出,y軸正方向0~1 m處的最大風(fēng)速在9~10 m/s左右,y軸正方向2~3 m處的最大風(fēng)速在7~9 m/s左右,風(fēng)速緩慢減小,各層最低風(fēng)速在2~3 m/s左右。風(fēng)速的分布規(guī)律為:巷道中間風(fēng)速較小,各層風(fēng)速大于7 m/s的區(qū)域?yàn)榄h(huán)狀分布,且向前環(huán)狀區(qū)域的范圍逐漸變大。

將曳力模型的計(jì)算與湍流模型的計(jì)算結(jié)果耦合進(jìn)行水霧粒子軌跡的計(jì)算,輔以K-H破碎模型計(jì)算霧滴的粒徑變化[13]。設(shè)置每只噴嘴在1 s內(nèi)噴出霧滴的數(shù)量為50,霧化模擬結(jié)果顯示了霧滴分布和運(yùn)動(dòng)軌跡,分析水霧分布的特點(diǎn)及對(duì)除塵的影響。

圖5的模擬結(jié)果顯示出了噴霧第2 s時(shí)水霧的運(yùn)動(dòng)情況及形成的結(jié)構(gòu)。噴嘴附近的水霧剛剛由噴嘴射出而沒(méi)有大量擴(kuò)散,因此霧滴分布密集,且噴嘴呈2層螺旋排列,在旋風(fēng)的作用下會(huì)形成2層明顯的螺旋霧幕并覆蓋巷道截面,2層霧幕相比單層霧幕可進(jìn)一步阻止粉塵的擴(kuò)散,提高除塵效果。水霧剛剛噴出時(shí)粒徑分布約在10 μm左右,而水霧也會(huì)受到風(fēng)流的影響破碎成為粒徑更小的霧滴。

圖5 噴霧2 s時(shí)水霧分布Fig.5 Water fog distribution at 2 s after spray

圖6為噴霧第5 s時(shí)粒子軌跡,此時(shí)2層霧幕仍然可見(jiàn),然而由于旋風(fēng)的作用,2層霧幕中的水霧很快擴(kuò)散并充斥整個(gè)巷道,并根據(jù)圖4可知風(fēng)速分布對(duì)于粒子軌跡的影響。風(fēng)速的不同決定了水霧疏密程度的不同,根據(jù)伯努利能量守恒原理,風(fēng)速大的區(qū)域氣壓相對(duì)較低,氣壓低的區(qū)域更易聚集水霧與粉塵,由此可見(jiàn)噴霧5 s后水霧分布與圖4的風(fēng)速分布相近,內(nèi)部風(fēng)速較快,水霧的流動(dòng)也較快,水霧分布也會(huì)比外部更密集,這便形成了喇叭狀的旋轉(zhuǎn)水霧。粉塵由于氣壓差的作用更易卷入水霧密集區(qū)域,在風(fēng)流作用下快速與霧滴結(jié)合,實(shí)現(xiàn)控塵效果。巷道墻壁附近也有旋轉(zhuǎn)水霧的分布,可對(duì)擴(kuò)散至低風(fēng)速取區(qū)的粉塵進(jìn)行控制,且由于風(fēng)速較小,擴(kuò)散至低風(fēng)速區(qū)水霧的粒徑變化率較小,有利于水霧擴(kuò)散至更遠(yuǎn)的區(qū)域,使除塵范圍更廣。

圖6 噴霧5 s時(shí)水霧分布Fig.6 Water fog distribution at 5 s after spray

2 試驗(yàn)?zāi)P图霸O(shè)備

2.1 掘進(jìn)工作面模型

試驗(yàn)將仿真掘進(jìn)工作面在掘進(jìn)過(guò)程中產(chǎn)塵和除塵過(guò)程,目的是將傳統(tǒng)噴霧和多層環(huán)繞霧幕的除塵效果進(jìn)行對(duì)比。本試驗(yàn)以長(zhǎng)、寬、高分別為4,4,4 m的試驗(yàn)棚作為掘進(jìn)巷道,前端模擬掘進(jìn)過(guò)程中的產(chǎn)塵過(guò)程,后端模擬噴霧過(guò)程。多層旋渦霧幕裝置設(shè)計(jì)如圖7所示,用以固定噴嘴的管路支架結(jié)構(gòu)與數(shù)值模擬中的噴嘴排列結(jié)構(gòu)相同。

圖7 多層旋渦霧幕試驗(yàn)?zāi)P驮O(shè)計(jì)Fig.7 Design of experimental model for multi-layer vortex fog curtain

2.2 主要試驗(yàn)設(shè)備和材料

如圖8所示,主要試驗(yàn)設(shè)備有:水泵、空氣壓縮機(jī)、連接10只SV882氣動(dòng)噴嘴的供水管與供氣管、螺旋狀管路支架、由鼓風(fēng)機(jī)和風(fēng)筒組成的發(fā)塵裝置、煤粉、粉塵濃度測(cè)試儀等。

圖8 試驗(yàn)平臺(tái)場(chǎng)景Fig.8 Scene of experimental platform

2.3 試驗(yàn)過(guò)程及方法

該試驗(yàn)在室內(nèi)進(jìn)行,試驗(yàn)棚的一端設(shè)置發(fā)塵裝置,由風(fēng)筒以及后端連接的鼓風(fēng)機(jī)組成,可將煤粉吹入棚內(nèi)。另一端模擬掘進(jìn)機(jī)的多層螺旋霧幕除塵裝置,由螺旋管和10只噴嘴組成,空壓機(jī)與水泵共同提供噴霧動(dòng)力,噴霧裝置距塵源2.5 m,且正對(duì)塵源。由于噴霧除塵時(shí)工作人員距產(chǎn)塵點(diǎn)較遠(yuǎn),且后端角落處極易聚集粉塵,所以測(cè)點(diǎn)設(shè)置在實(shí)驗(yàn)棚最左側(cè)后端的一角。噴霧時(shí)水泵的水壓設(shè)為0.12 MPa,空壓機(jī)氣壓設(shè)為0.52 MPa。

2.4 3種不同除塵方式試驗(yàn)對(duì)比

用秤分別量取3組1.5 kg的煤粉用來(lái)分別進(jìn)行3組試驗(yàn),第1組用自然沉降的方法,第2組用傳統(tǒng)噴霧除塵方法,第3組用多層螺旋霧幕除塵的方法。發(fā)塵裝置啟動(dòng)后大約45 s發(fā)塵完畢,測(cè)量3組試驗(yàn)的粉塵濃度變化情況,粉塵濃度測(cè)試儀放在測(cè)點(diǎn)處。

第1組:取1.5 kg的煤粉發(fā)塵,發(fā)塵完畢使煤塵布滿(mǎn)整個(gè)空間后,立即開(kāi)啟粉塵濃度測(cè)定儀進(jìn)行濃度測(cè)量,發(fā)塵完畢以后的粉塵濃度C0及自然沉降1,2,3 min時(shí)粉塵濃度C1,C2,C3。

第2組:將噴嘴排列成環(huán)形,噴射方向調(diào)至向前并向外輻射,噴射角度60°,進(jìn)行傳統(tǒng)噴霧除塵。將1.5 kg煤粉裝入發(fā)塵裝置發(fā)塵,發(fā)塵完畢后啟動(dòng)水泵和空壓機(jī),水霧從噴嘴噴出,與此同時(shí)測(cè)量在發(fā)塵完畢后的粉塵濃度C0及噴霧1,2,3 min時(shí)粉塵濃度C1,C2,C3。

第3組:將噴嘴固定在管路支架上,噴射方向?yàn)檠刂Ъ苈菪那芯€方向。取1.5 kg的煤粉發(fā)塵,發(fā)塵完畢后立即啟動(dòng)水泵和空壓機(jī),同時(shí)測(cè)出發(fā)塵完畢后的粉塵濃度C0及噴霧1,2,3 min時(shí)粉塵濃度C1,C2,C3。

記錄3組試驗(yàn)各個(gè)時(shí)間段的粉塵濃度,根據(jù)該濃度數(shù)據(jù),計(jì)算各組試驗(yàn)在除塵進(jìn)行3 min時(shí)的除塵速率vn,用公式(12)計(jì)算:

(12)

式中:vn為第n階段的除塵速率;C0為初始粉塵濃度,Cn為第nmin的粉塵濃度;h為除塵時(shí)間,min。

試驗(yàn)結(jié)果匯總見(jiàn)表2。

2.5 試驗(yàn)結(jié)果分析

由表2數(shù)據(jù)可知,1.5 kg煤塵發(fā)塵完畢后,試驗(yàn)棚內(nèi)測(cè)點(diǎn)處的初始粉塵濃度始終保持在470 mg/m3左右,浮動(dòng)微弱,可近似看做濃度相同。通過(guò)3組數(shù)據(jù)可知,傳統(tǒng)噴霧的降塵效果好于自然降塵,多層螺旋霧幕的除塵效果好于傳統(tǒng)噴霧除塵,其中僅用多層螺旋霧幕除塵的方法,在3 min后濃度降為3.68 mg/m3,達(dá)到了國(guó)家標(biāo)準(zhǔn)的4 mg/m3以下,呼塵的濃度也降到了1 mg/m3以下。在進(jìn)行雙層螺旋霧幕試驗(yàn)時(shí),可清晰觀察到2層密集水霧(拍照時(shí),照片顯示不明顯,故未附圖),霧滴粒徑更加細(xì)小,這與數(shù)值模擬結(jié)果一致,并且粉塵濃度變化結(jié)果也驗(yàn)證了該方案的正確性,說(shuō)明多層螺旋霧幕除塵可更加有效地降低粉塵濃度。

表2 3種除塵方式下的各個(gè)時(shí)間的粉塵濃度變化及除塵速率對(duì)比Table 2 Comparison of dust concentration variation and dust removal efficiency at each time by three dust removal modes

多層螺旋霧幕除塵相比傳統(tǒng)噴霧,有3點(diǎn)重要的優(yōu)勢(shì):

1)傳統(tǒng)噴霧除塵方法射出的氣流無(wú)法相互作用形成連續(xù)的風(fēng)流,容易造成氣流紊亂并停滯。且水霧擴(kuò)散速度較慢,水霧在巷道角落區(qū)域也僅處于緩慢沉降或布朗運(yùn)動(dòng)狀態(tài),無(wú)法使其與粉塵更快碰撞結(jié)合,造成除塵速率較低。而多層螺旋霧幕除塵方法由于可形成多層霧幕,旋風(fēng)的作用使多層霧幕中的水霧很快擴(kuò)散形成旋轉(zhuǎn)水霧,水霧的覆蓋范圍更廣且運(yùn)動(dòng)速度較快,可使霧滴與粉塵有更多的機(jī)會(huì)接觸。

2)風(fēng)流會(huì)使粒徑較大的霧滴破碎成多個(gè)粒徑較小的霧滴,粒徑較小的霧滴更易與粉塵碰撞結(jié)合,且霧滴的數(shù)量增多也可增加霧滴與粉塵結(jié)合的可能性,提高控塵效果。

3)高速旋風(fēng)形成負(fù)壓區(qū),由于氣壓的作用可使粉塵被吸入旋風(fēng)流中,使粉塵有更多機(jī)會(huì)與水霧結(jié)合。

3 結(jié)論

1)通過(guò)湍流模型與流體曳力模型耦合計(jì)算表明,采用多層螺旋霧幕的方法,將噴嘴以2層螺旋排列,可產(chǎn)生2層較大范圍的螺旋霧幕。且旋風(fēng)的作用可使霧幕中的水霧很快擴(kuò)散形成旋轉(zhuǎn)水霧,快速流動(dòng)的水霧可使粉塵迅速與霧滴結(jié)合,同時(shí)風(fēng)流也會(huì)對(duì)水霧起到一定的破碎作用,使霧滴與粉塵的碰撞幾率增加,提高控塵效果。

2)通過(guò)試驗(yàn)將自然降塵、傳統(tǒng)噴霧除塵、多層螺旋霧幕除塵的除塵速率對(duì)比,多層螺旋霧幕對(duì)全塵與呼塵的除塵效果明顯好于自然降塵與傳統(tǒng)噴霧除塵。且僅使用多層螺旋霧幕除塵3 min后可使近470 mg/m3的全塵濃度降至國(guó)家標(biāo)準(zhǔn)4 mg/m3以下,并可將呼塵的濃度降到1 mg/m3以下。多層螺旋霧幕的試驗(yàn)也驗(yàn)證了數(shù)值模擬的正確性。

猜你喜歡
水霧除塵粉塵
專(zhuān)利名稱(chēng):一種用于再生混凝土混料的除塵裝置及其除塵方法
細(xì)水霧滅火系統(tǒng)對(duì)油池火抑制效果研究
礦井膠帶蔓延火災(zāi)的細(xì)水霧抑制影響因素研究*
高壓粉塵防爆電機(jī)過(guò)厚粉塵層的形成與分析
自然之道
如何鑒別及選購(gòu)一臺(tái)好的石材除塵設(shè)備
高爐出鐵場(chǎng)粉塵高效捕集系統(tǒng)優(yōu)化及實(shí)踐
水鹽體系法生產(chǎn)硫酸鉀除塵方式選擇
粉塵大戰(zhàn)
大佛除塵
延长县| 澳门| 舟曲县| 行唐县| 个旧市| 胶南市| 定结县| 浦东新区| 万州区| 兰州市| 星座| 巴里| 临沂市| 利川市| 凤冈县| 桐梓县| 祁阳县| 永修县| 颍上县| 福安市| 杭州市| 聂拉木县| 南澳县| 津南区| 潮州市| 临安市| 永福县| 绥中县| 贵定县| 阳城县| 深圳市| 昭通市| 通化市| 云梦县| 石阡县| 东乌| 栾川县| 新兴县| 南京市| 故城县| 南郑县|