• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)及模型試驗研究

      2022-11-08 10:38:10關(guān)喜彬
      隧道建設(中英文) 2022年10期
      關(guān)鍵詞:寒區(qū)外界射流

      關(guān)喜彬

      (中鐵十九局集團第六工程有限公司,江蘇 無錫 214028)

      0 引言

      在高緯度、高海拔地區(qū)修建隧道的數(shù)量逐年增多,如吉圖暉鐵路沿線的草木溝隧道、五峰山隧道和哈爾巴嶺2號隧道等。大量寒區(qū)隧道建設和運營情況表明,低溫引起的凍害普遍存在于既有隧道中,這些凍害既對隧道結(jié)構(gòu)造成破壞,同時也給鐵路運營埋下安全隱患[1-3]。

      近年來,學者們采用多種方法來減少凍害對寒區(qū)隧道結(jié)構(gòu)的破壞和影響,并采取多種防凍措施來預防凍害發(fā)生[4]。例如:保溫層法是指在隧道內(nèi)鋪設輕質(zhì)隔熱材料,利用其導熱系數(shù)低的特點,減少洞內(nèi)冷空氣與圍巖之間的熱傳導,阻止寒冷氣流入侵圍巖內(nèi)部;但保溫層法不能隨著外界環(huán)境氣溫變化而變化,屬于被動的保溫方法,當外界環(huán)境溫度過低時,其保溫效果并不能令人滿意[5]。防寒門法是指在隧道洞口建造防寒門來阻隔寒冷氣流入侵隧道內(nèi)部,其保溫效果較好,但僅適應于交通流量小的隧道[6]。電伴熱法是指采用加熱電纜對隧道直接供暖,保溫效果較好,但后期維護能耗較大[7]。地熱源泵法是指采用地面換熱器吸收地熱能對隧道直接供暖,其保溫效果較好,適用于地熱能資源豐富的地區(qū)[8]。因此,針對現(xiàn)有保溫技術(shù)的局限性和目前寒區(qū)隧道凍害形勢的嚴峻性,有必要探索一種新型的保溫方法來解決隧道凍害問題。

      空氣幕裝置因其良好的空氣阻隔效果,最初被用于隔斷室內(nèi)外氣流的熱交換和控制粉塵污染的擴散[9]。20世紀30年代,前蘇聯(lián)科學家開展了大量空氣幕的空氣流場模型試驗,提出了空氣幕運動軌跡方程和阻隔效率的計算方法,為空氣幕的設計奠定了理論基礎(chǔ)[10]。2019年,Yang等[11]探究了不同風速、風向條件下空氣幕性能的變化規(guī)律。目前,空氣幕技術(shù)廣泛應用于冷庫、礦用巷道和超市等領(lǐng)域,但在寒區(qū)隧道保溫領(lǐng)域尚未有應用空氣幕的先例。

      目前,國內(nèi)有不少學者采用模型試驗的方法來探究如何解決寒區(qū)隧道凍害問題。夏才初等[12]開展?jié)u凍隧道演化模擬試驗,揭示寒區(qū)隧道凍害的發(fā)展規(guī)律。郭瑞等[13]開展寒區(qū)隧道縱向溫度場模型試驗,探究設防長度與洞內(nèi)風速的關(guān)系。高焱等[14]開展寒區(qū)高速鐵路隧道溫度場模型試驗,揭示了列車運行條件下隧道洞內(nèi)溫度場的變化規(guī)律。渠孟飛等[15]開展寒區(qū)隧道襯砌凍脹力的模型試驗,得到巴郎山隧道襯砌凍脹力的變化規(guī)律。但目前尚未發(fā)現(xiàn)關(guān)于空氣幕作用下保溫室內(nèi)空氣流場和溫度場方面的報道,因此開展相關(guān)的模型試驗具有一定的創(chuàng)新意義和工程應用價值。

      為預防凍害的發(fā)生,本文提出一種新型的寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng),設計并研制不同外界風速、風幕機射流風速和射流角度等條件下的試驗系統(tǒng),探究保溫裝置內(nèi)空氣流場和溫度場的變化規(guī)律。

      1 隧道側(cè)吹式空氣幕保溫系統(tǒng)及控制方程

      1.1 寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)

      天平鐵路關(guān)山隧道全長15 634 m,位于甘肅省天水市與平?jīng)鍪薪唤绲囟?。關(guān)山隧道極端最低氣溫為-25.5 ℃,年平均降水量為568.5 mm,年最大降水量為824.4 mm,最大季節(jié)凍結(jié)深度為79 cm。為預防關(guān)山隧道凍害的發(fā)生,提出一種新型的寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)。該保溫系統(tǒng)由保溫裝置、側(cè)吹式風幕機、溫度和風速風向傳感器及PLC智能控制器4部分組成。當隧道洞口的溫度和風速檢測傳感器檢測到外界氣溫和風速達到設定值時,PLC智能控制器控制風幕機的工作狀態(tài),調(diào)整空氣幕的射流角度、射流風速、射流溫度等參數(shù),通過風幕機吹出的豎向強風,阻隔寒冷氣流入侵隧道內(nèi)部,確保隧道洞口的溫度達到0 ℃以上,預防隧道凍害的發(fā)生。當隧道采用縱向式通風時,隧道中的射流風機從洞口引入新鮮空氣,通過另外一個洞口排出。此時,空氣幕保溫系統(tǒng)中的風速風向傳感器監(jiān)測到風向由隧道內(nèi)部向外部流動,關(guān)閉空氣幕,使隧道順利完成通風。寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)如圖1所示。

      圖1 寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)

      1.2 隧道空氣幕阻隔效率控制方程

      隧道空氣幕阻隔效率控制方程的計算模型見圖2。

      圖2 空氣幕阻隔效率計算模型

      如圖2所示,設定空氣幕的氣流為f,隧道外流入隧道的氣流為f1,空氣幕傾斜平射氣流為f2[10],則:

      f=f1+f2。

      (1)

      (2)

      (3)

      式(2)—(3)中:ω為隧道外自然風速,m/s;L為隧道的最大跨徑,m;H為隧道平均高度,m;ω0為射流風速,m/s;B1為空氣幕噴口長度,m;ε為湍流系數(shù);α為射流角度,(°);b0為空氣幕射流厚度,m;x、y為距離坐標原點的距離。

      當由隧道外流入隧道的氣流與空氣幕傾斜平射氣流合成的空氣幕氣流垂直于x軸時,則在空氣幕作用下,外界氣流與空氣幕的合成氣流

      (4)

      外界氣流與空氣幕的合成氣流fG同樣可表現(xiàn)為空氣幕射流量Q0和空氣幕工作時未能阻擋的空氣量Qm的總和:

      fG=Q0+Qm。

      (5)

      空氣幕射流量

      Q0=ω0B1b0。

      (6)

      空氣幕工作時未能阻擋的空氣量

      (7)

      空氣幕不工作時,侵入隧道的總風量Q等同于隧道外流入隧道的氣流f1,表示為:

      (8)

      得到空氣幕阻隔效率

      (9)

      保溫系統(tǒng)參數(shù)標定的依據(jù)是以阻隔效率最大為標準,為進一步研究空氣幕射流角度、射流風速和射流厚度等參數(shù)與阻隔效率之間的最優(yōu)關(guān)系,通過最優(yōu)關(guān)系找出最優(yōu)設計方案,開展不同空氣幕參數(shù)與阻隔效率之間的研究。計算參數(shù)如表1所示。

      表1 計算參數(shù)

      1.3 射流角度與阻隔效率之間的變化規(guī)律

      空氣幕射流角度變化范圍為0°~90°,假設洞外風速為1.5 m/s,射流風速為6 m/s,射流厚度為0.05 m,將以上計算參數(shù)代入式(9),計算射流角度與阻隔效率的關(guān)系,再采用Matlab軟件繪制射流角度與阻隔效率之間的變化規(guī)律,見圖3。

      圖3 射流角度與阻隔效率之間的變化規(guī)律曲線

      如圖3所示,阻隔效率隨射流角度的變化趨勢為先遞增、后遞減,其最大值為1.091,此時對應風幕機最佳的射流角度為77.9°。當空氣幕阻隔效率大于等于1時,說明空氣幕對空氣流動的阻隔效率為100%,也說明此時空氣幕的開啟功率過大。

      1.4 外界風速與阻隔效率之間的變化規(guī)律

      外界風速為0~12 m/s,假設射流角度為77.6°,射流風速為6 m/s,射流厚度為0.5 m,探究并采用Origin數(shù)據(jù)分析軟件繪制外界風速與阻隔效率之間的變化規(guī)律,見圖4。

      圖4 外界風速與阻隔效率之間的變化規(guī)律曲線

      由圖4可知,外界風速與阻隔效率之間呈反比例關(guān)系;當外界風速為1.5 m/s,射流風速為6 m/s時,阻隔效率為1.091,說明射流風速為4倍及以上的外界風速時,阻隔效果好。

      1.5 射流風速與阻隔效率之間的變化規(guī)律

      假設射流風速為0~20 m/s,外界風速為1.5 m/s,射流角度為77.6°,射流厚度為0.5 m時,探究并采用Origin數(shù)據(jù)分析軟件繪制射流風速與阻隔效率之間的變化規(guī)律,見圖5。

      圖5 射流風速與阻隔效率之間的變化規(guī)律曲線

      由圖5可知,射流風速與阻隔效率之間呈正比例關(guān)系;當外界風速為1.5 m/s,射流風速為5.5 m/s時,此時阻隔效率為0.99,說明此種情況可以有效阻隔外界寒冷氣流入侵隧道內(nèi)部。

      2 試驗系統(tǒng)

      隧道內(nèi)的空氣流場和溫度場是非常復雜的對流換熱和圍巖熱傳導耦合問題,考慮空氣幕作用后使得問題更為復雜。因此,為了簡化,作出如下假定:1)氣流為不可壓縮的理想流體;2)氣體流動為湍流;3)氣流、圍巖均為連續(xù)介質(zhì)。

      2.1 相似比設計

      開展模型試驗研究時,需要滿足幾何相似、運動相似和熱力學相似3個相似準則,綜合考慮試驗精度和可實施性[16-17],該試驗的幾何相似比確定為lm∶lp=1∶20(下角標m和p分別表示模型和原型)。

      在確保模型幾何相似的基礎(chǔ)上,達到與原型的運動相似,模型與原型的阿基米德數(shù)、雷諾數(shù)這2個準則數(shù)應該相等;同時還需要滿足熱力學相似。即在幾何和運動相似的基礎(chǔ)上,模型與原型在空氣熱傳導、熱對流和熱輻射3種熱傳遞方式上相同,因此模型與原型的貝克萊準則應該相等。

      雷諾數(shù)是用于判斷流體是層流還是湍流的相似準則。當雷諾數(shù)較小時,黏滯力對流場的影響大于慣性,流體為層流狀態(tài);隨著雷諾數(shù)的增加,慣性對流場的影響大于黏滯力,流體逐漸由層流狀態(tài)轉(zhuǎn)變成湍流狀態(tài);當雷諾數(shù)繼續(xù)增大時,流體湍流的激烈程度幾乎不再變化,其值只與相對粗糙度有關(guān),處于相同自模區(qū)內(nèi)流體的流動狀態(tài)相似。保溫室原型和模型內(nèi)壁面材質(zhì)分別為混凝土和保溫板,原型的相對粗糙度為0.002,模型的相對粗糙度為0.03。由莫迪圖可知,當相對粗糙度分別為0.002和0.03時,雷諾數(shù)Re分別大于600 000和30 000時,即流體進入第2自模區(qū)的湍流穩(wěn)定狀態(tài)。因此,針對本次試驗,模型流場的相似比主要由阿基米德數(shù)準則決定,即:

      (Ar)p=(Ar)m。

      (10)

      (11)

      (12)

      式(10)—(12)中:(Ar)p、(Ar)m分別為原型和模型的阿基米德數(shù);g為重力加速度;lp、lm分別為原型與模型的保溫室長度;Tp、Tm分別為原型與模型的保溫室溫度;ΔTp、ΔTm分別為原型與模型的保溫室的溫差;vp、vm分別為原型與模型的風速;C1為試驗的幾何相似比。

      2.2 試驗系統(tǒng)

      試驗系統(tǒng)由外界環(huán)境控制系統(tǒng)、隧道模型、隧道空氣幕保溫裝置和數(shù)據(jù)監(jiān)測系統(tǒng)4部分組成。試驗系統(tǒng)原理如圖6所示。

      圖6 試驗系統(tǒng)原理圖

      空氣幕裝置安裝在保溫室前端兩側(cè),空氣幕的射流風速為1~20 m/s,射流溫度為30~55 ℃,通過抽取外界氣流,經(jīng)空氣幕加熱后,以一定的角度噴射,阻隔外界冷空氣的侵入。風速傳感器布置在保溫室內(nèi),共設置6個測點,用于實時測量侵入隧道的風速,風速傳感器的測量范圍為0.3~30 m/s,精度為0.1 m/s。溫度傳感器布置在保溫室后端,共設置6個測點,用于實時測量溫度的變化情況,溫度傳感器的測量范圍為-55~125 ℃,精度為0.1 ℃。溫度傳感器和風速傳感器測點布置如圖7所示。

      圖7 溫度傳感器和風速傳感器測點布置圖

      試驗系統(tǒng)的搭建主要有以下步驟:1)搭建外界環(huán)境控制系統(tǒng),其溫度調(diào)控范圍為0~-30 ℃;2)澆筑隧道襯砌結(jié)構(gòu),按1∶20的幾何相似比例制作隧道襯砌;3)在隧道前端安裝空氣幕保溫室;4)將風速傳感器和溫度傳感器安裝在保溫室內(nèi)。隧道模型及空氣幕保溫室見圖8和圖9。

      圖8 隧道模型

      圖9 空氣幕保溫室

      3 試驗結(jié)果

      3.1 空氣幕最優(yōu)射流角度

      選取測點5和測點6為研究對象,當外界風速分別為4 m/s和5 m/s、空氣幕的射流厚度為5 cm、射流風速為10 m/s時,改變空氣幕的射流角度,不同射流角度下保溫室內(nèi)風速的變化規(guī)律如圖10和圖11所示。

      圖10 不同外界風速情況下測點5的實測數(shù)據(jù)

      圖11 不同外界風速情況下測點6的實測數(shù)據(jù)

      由試驗測試結(jié)果可知,測點5和測點6的最優(yōu)射流角度為75°,采用式(9)計算得到空氣幕最優(yōu)射流角度為77.9°,兩者得到的最優(yōu)射流角度結(jié)果相近。由此可知,隧道空氣幕阻隔效率控制方程可靠可信,可以用于寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)的參數(shù)優(yōu)化設計。

      3.2 外界風速與阻隔效率之間的關(guān)系

      選取測點1—6為研究對象,當外界風速分別為1.5、3、4、5 m/s,空氣幕的射流厚度為5 cm,射流風速為10 m/s,射流角度為77.9°時,得出不同外界風速情況下保溫室內(nèi)風速的變化規(guī)律,如圖12所示。

      由試驗測試結(jié)果可知,空氣幕能夠有效阻隔外界氣流的侵入。隨著外界風速的增加,保溫室內(nèi)風速也隨之增加,即空氣幕的阻隔效率逐漸降低;由測點1—6的實測數(shù)據(jù)可知,保溫室底部的風速最大,中下位置的風速最小。因此,在實際工程應用時,該系統(tǒng)的風速監(jiān)測器應安裝于隧道底部。

      3.3 射流風速與阻隔效率之間的關(guān)系

      采用控制變量法,當空氣幕的射流厚度b0為5 cm,外界風速ω為1.5、3、5 m/s,射流角度為0°~90°時,分別研究射流風速分別為10 m/s和12 m/s條件下保溫室內(nèi)風速變化規(guī)律,如圖13和圖14所示。

      由圖13和圖14試驗結(jié)果對比分析可知,當射流風速增大時,同一測點位置的風速均有所下降,阻隔效率有所增加;這與圖5理論計算中阻隔效率隨著射流風速的增加而增大的規(guī)律相同。相同射流厚度和外界風速條件下,射流風速對射流角度的最優(yōu)值影響較小。當外界風速為1.5 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為55°~75°;當外界風速為3 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為65°~75°;當外界風速為5 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為75°左右。由此可見,不同的外界風速越小,最優(yōu)射流角度越小,式(9)所計算出的最優(yōu)射流角度77.9°可用于寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)的參數(shù)優(yōu)化設計。工程實際應用中,可根據(jù)外界情況和洞內(nèi)風速傳感器適當降低射流角度。

      (a)b0=5 cm, ω0=10 m/s, ω=1.5 m/s

      (a)b0=5 cm, ω0=12 m/s, ω=1.5 m/s

      3.4 側(cè)吹式射流保溫裝置保溫效果分析

      以關(guān)山隧道2021年最冷月平均氣溫-19.2 ℃為環(huán)境氣溫,選取溫度傳感器測點1—6為研究對象,當外界風速為1.5 m/s、空氣幕的射流角度為75°、射流厚度為5 cm、射流風速為10 m/s、射流溫度為45 ℃時,空氣幕開啟后保溫室內(nèi)溫度場的變化情況如圖15所示。

      圖15 保溫室內(nèi)溫度傳感器測點1—6的實測數(shù)據(jù)

      當隧道外界的環(huán)境氣溫為-19.2 ℃時,保溫室內(nèi)的溫度低于0 ℃,此種條件下常溫空氣幕無法滿足隧道防寒保溫的需求。當熱溫空氣幕開啟8 min后,保溫室內(nèi)空氣除底部測點3和測點6外,其余部分均開始出現(xiàn)正溫,此種條件下進入隧道洞口的氣溫已經(jīng)達到0 ℃以上。由此可見,熱溫空氣幕具有較好的保溫效果,有助于預防關(guān)山隧道凍害的發(fā)生。

      4 經(jīng)濟性對比

      以關(guān)山隧道為例,對比分析保溫層法與空氣幕保溫系統(tǒng)100年使用年限內(nèi)的經(jīng)濟費用。

      4.1 保溫層法總費用

      保溫層法的總成本主要包括保溫材料成本、人工成本和維護成本。關(guān)山隧道現(xiàn)場冬季最冷月平均氣溫為-25.5 ℃,文獻[18]給出了保溫層鋪設長度的計算式,即:

      y=-0.8128x2-53.448x+150.73。

      (13)

      式中:y為隧道入口處保溫層鋪設長度,m;x是最冷月份的平均溫度,℃。

      由式(13)計算得出關(guān)山隧道保溫層鋪設長度為985.13 m。假設隧道入口和出口處保溫層的鋪設長度相等,則隧道所需保溫層的總長度為1 970.26 m;關(guān)山隧道每m鋪設0.05 m厚隔熱層所需材料約為0.86 m3;保溫層材料選用硬質(zhì)聚氨酯,市場價格約1 300元/m3。根據(jù)計算,保溫材料成本為220.25萬元;鋪設保溫層的人工成本為4 500 元/m,因此,施工人工成本為886.62萬元。

      硬質(zhì)聚氨酯保溫材料的使用壽命為25年,在隧道100年的使用期內(nèi)需要更換3次保溫材料。人工拆除保溫層的費用為167元/m3,計算得出拆除費用為32.9萬元。

      因此,隧道運營100年,保溫層法總成本為4 526.18萬元。

      4.2 空氣幕系統(tǒng)總費用

      空氣幕系統(tǒng)主要由保溫裝置、側(cè)吹式風幕機、溫度和風速風向傳感器和PLC智能控制器4部分組成。保溫裝置主體結(jié)構(gòu)采用鋼結(jié)構(gòu)保溫框架,建設費用約為10.6萬元,使用年限為50年,保溫裝置主體結(jié)構(gòu)全周期總費用共計21.2萬元。側(cè)吹式風幕機采用定制型可調(diào)熱風幕機,單臺價格約為3.9萬元,2臺共7.8萬元,使用年限為10年,空氣幕全周期總費用共計78萬元。PLC智能控制器費用為3.6萬元,使用年限為10年,在隧道100年全周期總費用為36萬元。

      空氣幕系統(tǒng)年耗電量主要包括空氣幕和PLC智能控制器。PLC智能控制器耗電量較低,因此只計算空氣幕在凍結(jié)期150 d內(nèi)全功率運行所需的電量和電費。單臺風幕機功率為9.9 kW,2臺風幕機功率為19.8 kW。年凍結(jié)期內(nèi)全功率運行的功率為71 280 kW,按照中國工業(yè)用電0.86元/kW·h,計算可得風幕機年用電費用為61 300.8元,即風幕機在隧道100年全周期中用電總費用為613萬元。

      因此,隧道運營100年,隧道出入口處空氣幕系統(tǒng)總成本為1 496.4萬元。

      對比可知,在隧道100年全周期中,保溫層法總成本為4 526.18萬元,空氣幕系統(tǒng)總成本為1 496.4萬元。根據(jù)本文計算的費用,在隧道內(nèi)采用空氣幕系統(tǒng)進行防寒保溫具有很大的經(jīng)濟效益。但本文僅粗略估算了空氣幕系統(tǒng)的成本,具體施工及用電成本需結(jié)合現(xiàn)場實際情況進行詳細計算。

      5 結(jié)論與討論

      1)根據(jù)流函數(shù)疊加原理推導出寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)空氣幕阻隔效率的控制方程,并在模型試驗的基礎(chǔ)上驗證了該控制方程的正確性。該控制方程可用于優(yōu)化寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng)的參數(shù)設計。

      2)為預防和控制關(guān)山隧道凍害的發(fā)生,提出寒區(qū)隧道側(cè)吹式空氣幕保溫系統(tǒng),設計并研制不同外界風速、風幕機射流風速和射流角度等條件下的試驗系統(tǒng),為解決寒區(qū)隧道凍害問題提供新方法和試驗數(shù)據(jù)支撐。

      3)模型試驗結(jié)果表明,當外界風速為1.5 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為55°~75°;當外界風速為3 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為65°~75°;當外界風速為5 m/s時,不同射流風速下空氣幕的最優(yōu)射流角度為75°左右。

      4)模型試驗結(jié)果表明,當外界環(huán)境溫度為-19.2 ℃時,常溫空氣幕難以滿足隧道保溫需求。熱空氣幕開啟8 min后,保溫室內(nèi)溫度均達到0 ℃以上,能夠有效預防隧道凍害的發(fā)生。

      猜你喜歡
      寒區(qū)外界射流
      修德箴言
      深海逃逸艙射流注水均壓過程仿真分析
      低壓天然氣泄漏射流擴散特性研究
      煤氣與熱力(2022年4期)2022-05-23 12:45:00
      外界誘惑大 讓你轉(zhuǎn)錢必有詐
      云南畫報(2022年4期)2022-05-05 06:00:26
      不同施肥處理對寒區(qū)苜蓿草地產(chǎn)量的影響
      人體對外界環(huán)境的感知
      寒區(qū)某水電站大管徑PCCP管施工過程
      設計的外界包容性分析
      人間(2015年20期)2016-01-04 12:47:20
      高原寒區(qū)野外駐訓衛(wèi)勤保障回顧性探討
      駐寒區(qū)某部士兵常見病調(diào)查及對策建議
      新巴尔虎左旗| 兰西县| 鹤峰县| 高陵县| 壶关县| 金寨县| 施甸县| 石河子市| 和田县| 博乐市| 松原市| 凌源市| 寿光市| 巫山县| 甘孜县| 平果县| 新兴县| 仙桃市| 襄城县| 高安市| 浙江省| 通化县| 镇安县| 赞皇县| 邛崃市| 新郑市| 彰武县| 绥滨县| 安吉县| 安阳市| 张家港市| 郓城县| 沙河市| 九龙县| 清河县| 施秉县| 龙泉市| 桃园市| 云和县| 颍上县| 武宣县|