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

?

天然氣水合物鉆采一體化模擬實驗系統(tǒng)及降壓法開采初步實驗

2019-09-03 06:45劉昌嶺李彥龍劉樂樂胡高偉吳能友孟慶國
天然氣工業(yè) 2019年6期
關鍵詞:模擬實驗產(chǎn)氣水合物

劉昌嶺 李彥龍 劉樂樂 胡高偉 陳 強 吳能友 孟慶國

1. 自然資源部天然氣水合物重點實驗室·青島海洋地質(zhì)研究所

2. 青島海洋科學與技術試點國家實驗室海洋礦產(chǎn)資源評價與探測技術功能實驗室

0 引言

海域天然氣水合物(以下簡稱水合物)試采成本昂貴,水合物開采技術研究仍以室內(nèi)模擬實驗為主?,F(xiàn)有的水合物開采技術實驗研究通常在較小尺寸的模擬實驗裝置中進行,并多關注較小尺度下水合物分解產(chǎn)氣和產(chǎn)水規(guī)律[1-6]。由于反應釜樣品尺寸小,導致明顯的邊界效應,且主導物理效應與現(xiàn)場實際情況不同,實驗結果難以在現(xiàn)場中得到應用。因此,研發(fā)大尺寸的水合物綜合開采實驗系統(tǒng)是非常有必要的。目前,已有多個國家研制了一批較大尺寸的水合物開采模擬實驗系統(tǒng),主要包括日本產(chǎn)業(yè)技術綜合研究所研制的大型實驗裝置[7]、德國波茨坦亥姆霍茲中心建立的實驗系統(tǒng)[8]、中國科學院廣州能源所研制的實驗系統(tǒng)[9]和中國石油大學(北京)研制的實驗系統(tǒng)[10],這些實驗系統(tǒng)重點關注水合物藏的形成機理及其分解產(chǎn)氣規(guī)律研究。近期西南石油大學依據(jù)水合物固態(tài)流化開采法的工藝流程,建立了大型海洋水合物固態(tài)流化開采實驗室[11],并初步開展水合物固態(tài)流化開采仿真[12]。

現(xiàn)有的水合物開采技術研究通常關注水合物分解產(chǎn)氣和產(chǎn)水規(guī)律、溫度、壓力和水合物飽和度關系研究等,而對水合物開采的產(chǎn)砂預測與防治等問題的關注較少[13-15]。為此,青島海洋地質(zhì)研究所研發(fā)了一套大型水合物鉆、采一體化模擬實驗系統(tǒng),可模擬實際地質(zhì)條件制備海洋水合物樣品,通過電阻層析成像技術實時探測水合物成藏與分布情況,研發(fā)了鉆井實驗模塊,可模擬鉆井、降壓開采工藝與過程,實時測量開采過程中出砂與管道流動等過程中多個物理參數(shù)變化,實現(xiàn)試采全過程仿真實驗模擬。筆者將重點介紹該系統(tǒng)的研發(fā)思路和主要構成,以及基于該實驗系統(tǒng)的初步實驗結果,以期為后續(xù)開展南海多類型水合物開采工藝實驗研究奠定基礎。

1 模擬實驗系統(tǒng)的研發(fā)

1.1 主要功能與結構組成

水合物鉆采一體化模擬實驗系統(tǒng)具備以下功能:①制備符合自然條件的含水合物沉積層;②實時監(jiān)測水合物沉積層整體和局部的水合物飽和度;③模擬井筒鉆進樣品過程中的儲存動態(tài)響應規(guī)律;④實施多級分步降壓開采過程;⑤實現(xiàn)產(chǎn)氣、產(chǎn)水、產(chǎn)砂過程的實時監(jiān)測與分離處理。

該系統(tǒng)采用模塊化設計,主要由主體高壓裝置、鉆采一體化、氣液供給、圍壓加載、回壓控制、氣液固分離及在線監(jiān)測、溫度控制、開采工作制度控制、數(shù)據(jù)測控與后處理等模塊組成(圖1、2)。

1.2 各子系統(tǒng)結構組成及主要技術參數(shù)

1.2.1 主體高壓裝置模塊

主體高壓裝置模塊是整個系統(tǒng)的核心,主要由高壓反應釜、圍壓膠套、壓力測柱、鎧裝溫度測柱、電阻層析成像測點等部分組成。

其中,高壓反應釜內(nèi)部尺寸為?750 mm×1 180 mm,容積為521 L,采用整體鍛造工藝加工。反應釜設計壽命為30年,工作壓力30 MPa,工作溫度范圍-20~50℃。高壓反應釜配置專用工作平臺,方便各測點的布置、填充物及釜蓋的裝卸等;高壓反應釜自帶制冷系統(tǒng),2~3 h的時間可使高壓反應器從室溫降低至-20℃;高壓反應器溫度控制精度為±0.1℃。

為滿足不同地層圍壓條件的模擬,反應釜內(nèi)筒采用氯丁橡膠制成的圍壓膠套,內(nèi)部尺寸?600 mm×1 145 mm,厚度10 mm,膠套上安裝320個電阻率層析成像測點。膠套內(nèi)置含水合物沉積物的樣品尺寸為?600 mm×1 000 mm,容積為282.6 L。圍壓膠套組件能夠承受井筒鉆進和開采模擬過程中產(chǎn)生的扭剪力,保證圍壓膠桶與傳感器以及上下端蓋連接處的密封性。

壓力測柱、鎧裝溫度測柱、電阻層析成像測點是主體高壓裝置模塊的主要測試單元。其中壓力測柱為測管式組件,分別插入沉積物4個層面,4個層面距離相隔均為200 mm;平面布置以軸線為中心,90°對稱距離分別為280 mm和440 mm。鎧裝溫度測柱由5個PT 100A鉑電阻裝入l根316 L的?10 mm鋼管中,測量柱表面粗糙化處理,防止氣液沿壁串流,精度為0.1 K。系統(tǒng)內(nèi)置6根鎧裝溫度測柱以獲取水合物成藏—開采過程中的溫度場分布。

電阻層析成像(ERT)技術可以利用不同介質(zhì)(鹽水、氣體、沉積物、水合物)的電導率差異,識別處于敏感場的電導率分布,并進行三維成像,從而識別模擬儲層中的水合物分布規(guī)律[16-18]。該系統(tǒng)選用四點法電阻層析成像測量系統(tǒng)[19],測量電極安裝布局在圍壓膠套上,電極沿樣品軸向間隔50 mm布置20層,每層電極個數(shù)為16個,共320個電極,每層電極均布在膠桶圓周上。

1.2.2 鉆采一體化模塊

鉆采一體化模塊是該套模擬實驗系統(tǒng)最具特色的部分,在國際天然氣水合物模擬裝置上尚屬首創(chuàng)。該模塊可實現(xiàn)以下功能:①模擬鉆具在高壓狀態(tài)下鉆入水合物儲層的現(xiàn)場鉆采工藝過程;②模擬鉆井液循環(huán)過程中水合物儲層的物化參數(shù)響應規(guī)律;③實現(xiàn)高壓狀態(tài)下鉆進的同時完成井下防砂工具安裝;④模擬不同工況下儲層水合物分解產(chǎn)出過程中的井筒攜砂、攜液流動規(guī)律。

圖1 天然氣水合物鉆采一體化模擬實驗系統(tǒng)組成示意圖

圖2 天然氣水合物鉆采一體化模擬實驗系統(tǒng)圖

為了滿足上述功能需求,鉆采一體化模塊設計了相應的鉆機總成、循環(huán)注入泵、流量計、壓力計、溫度計、背壓調(diào)節(jié)閥、截止閥、固液分離器、泥漿罐、攪拌器及鉆井液制冷機組,工藝流程如圖3所示,鉆機動力鉆具與主體高壓裝置模塊滿足動密封需求。

鉆采一體化模塊的基本設計工藝參數(shù)為:鉆機旋轉轉速50~400 r/min(無極調(diào)速),動密封最大允許鉆速6 m/min,工作壓力30 MPa;循環(huán)管路流量6 L/min,循環(huán)注入泵選用三缸柱塞泵,流量400 L/h、電機功率15 kW;泥漿攪拌速度400 r/min、攪拌功率0.5 kW;制冷速率:泥漿罐裝滿額定泥漿0.5~1.0 h溫度從室溫降至-20 ℃,與高壓反應釜水夾套制冷機組共用。

圖3 鉆采一體化模塊循環(huán)流程圖

模擬動力鉆具與模擬開采井管同心安裝,鉆頭帶動模擬動力鉆具和模擬開采井管鉆入模擬水合物儲層,鉆達預定深度后導通外部流程,開采過程中鉆桿充當內(nèi)部油管,外部模擬開采井筒上根據(jù)實際的地層需求安裝篩網(wǎng)或礫石包裹層,以模擬開采過程中的出砂防控情況。為了達到模擬實際試采現(xiàn)場的目的,井管表面噴涂絕熱絕緣涂層并作表面粗糙化處理,防止氣液沿壁串流、熱量流失,特別是防止對層析成像電阻率測試場的干擾等。

1.2.3 氣液供給模塊

氣液供給模塊用于實驗過程中天然氣和孔隙水的供給回收,由氣體供給子模塊和液相供給子模塊構成。液相供給子模塊包括液體容器、高壓注入泵、流量計及相應的高壓管閥等,供液速率介于0~500 mL/min,計量精度為±0.5 mL/min。為了實驗模擬降壓、注熱、注化學劑等多種開采過程,液相供給子模塊的管閥均經(jīng)過耐化學劑涂層處理,并配備相應的蒸汽發(fā)生器。

氣體供給子模塊由儲氣罐組、空壓機、增壓泵、壓力調(diào)節(jié)器、氣體流量計及管閥件組成,用于水合物生成時,向反應釜提供穩(wěn)定壓力與流量的天然氣;并滿足CO2置換開采過程中CO2的定壓、定速率注入。

為滿足持續(xù)成藏—開采過程的實驗,氣液供給模塊與產(chǎn)出分離計量模塊之間設計循環(huán)回路。

1.2.4 圍壓加載模塊

圍壓加載模塊用于提供水合物成藏—開采過程模擬中的儲層圍壓,模塊設計壓力高于系統(tǒng)主體高壓裝置模塊設計壓力5 MPa。系統(tǒng)設計定壓跟蹤和定壓差跟蹤2種模式,根據(jù)設定值自動跟蹤。主要由液體容器、高壓注入泵、圍壓背壓閥及相應的高壓管件組成。

1.2.5 開采工作制度控制模塊

水合物開采工作制度的2個關鍵仿真參數(shù)分別是模擬油嘴開度和井口油壓。通過油嘴開度和井口油壓的控制實現(xiàn)水合物降壓開采過程中井底流壓、生產(chǎn)壓降幅度、壓降速率的模擬。為此,采用帶計量刻度的針形閥安裝在氣液分離器之后的氣路流通通道,實現(xiàn)模擬油嘴的主要功能。

井口油壓控制主要通過回壓控制閥和回壓緩沖容器等設備的聯(lián)合實現(xiàn),壓力控制精度不大于0.1 MPa。通過合理設計,解決溫度變化的高壓實驗系統(tǒng)出口流動狀態(tài)不穩(wěn)定的問題,采用全自動回壓控制閥,使降壓開采時回壓控制可調(diào)。其基本工作原理為:計算機設置壓力值,控制傳感器反饋管道中的壓力信號,調(diào)節(jié)器和氣動調(diào)節(jié)氣源配合使用調(diào)節(jié)管道壓力;當管道壓力高于或低于設定壓力值時,通過傳感器反饋給計算機,計算機自動調(diào)節(jié)壓力值,使兩端壓力值達到期望水平,從而確保出口流動狀態(tài)相對穩(wěn)定。

產(chǎn)出分離計量模塊主要用于水合物開采過程中氣—液—固三相的分離及實時計量,主要由兩級氣液固三相分離器、背壓式減壓閥、Y形過濾器、光柵液位計、在線激光粒度儀、電子天平、緩沖氣瓶、氣體流量計、截止閥及管線等組成。

分離后的氣體通過氣體流量計持續(xù)計量,光柵液位計實時記錄水砂混合液體積、估算混合體系中的固相含量;采用一定時間間隔分離取樣的手段進行出砂量的計量。在線激光粒度儀可以實時在線記錄產(chǎn)出地層砂粒徑隨時間的變化規(guī)律。

1.2.7 溫度控制模塊

溫度控制模塊是水合物成藏、開采過程中儲層溫度模擬的必要模塊,包括高壓反應釜冷卻水夾套、制冷機組、電加熱器、循環(huán)泵等設備以及步進式低溫恒溫室,高壓反應器溫度控制精度為±0.1 ℃;供液模塊溫度控制精度為±1 ℃;步進式低溫恒溫室的尺寸為6 m×6 m×4 m(長×寬×高),步進式低溫恒溫室主要放置主體高壓裝置,室內(nèi)工作溫度介于-20 ℃~室溫,溫度控制精度為±2 ℃。

1.2.8 數(shù)據(jù)測控與后處理模塊

數(shù)據(jù)測控與后處理模塊主要由測控硬件及水合物開采室內(nèi)測控軟件2個部分構成。測控硬件主要包括溫度傳感器、壓力傳感器、流量控制器、電阻層析成像測量儀、氣液分離計量系統(tǒng)等數(shù)據(jù)測量設備和計算機、打印機、A/D采集卡、I/O控制板、軟件等數(shù)據(jù)采集設備。高壓反應釜溫度測量精度為±0.1℃,壓力測量精度為±0.1 MPa,出口氣體測量精度為±1.0 L/min,出口液體測量精度為±0.5 mL/min。測量參數(shù)包括出砂速率、出砂量、出砂粒徑、產(chǎn)水速率和產(chǎn)氣速率以及反應體系的溫度與孔隙壓力。

水合物開采室內(nèi)測控軟件(著作權登記號:2018SR625845)平臺是配套研發(fā)的,其主要功能模塊包括:①水合物開采過程仿真控制模塊;②水合物儲層多物理場在線監(jiān)測與反演模塊;③水合物出砂模擬與預測模塊;④水合物儲層產(chǎn)氣、產(chǎn)水規(guī)律測量模塊;⑤水合物開采工作制度調(diào)控模塊;⑥水合物開采實際工況仿真模塊等。該平臺專門用于水合物開采過程大尺寸室內(nèi)模擬仿真的測控,具備數(shù)據(jù)采集、仿真數(shù)據(jù)庫管理、數(shù)據(jù)粗處理、二維/三維數(shù)據(jù)呈現(xiàn)、開采降壓程序控制、圍壓自動跟蹤控制等功能。

2 初步模擬實驗與主要認識

CO2水合物具有與甲烷水合物相似的成藏、分解規(guī)律,且其生成條件相對溫和,成藏效率高,安全性強。因此,為了驗證實驗系統(tǒng)的整體可靠性,筆者以CO2水合物為例,初步展示冰點附近CO2水合物的降壓分解規(guī)律,分析篩管預埋工況下的氣—水—砂產(chǎn)出規(guī)律及儲層溫度場、壓力場的演化規(guī)律。

雙乙烯酮是一種極為活潑的乙?;瘎?,也是有機化學合成中一種非常重要的中間體,目前已廣泛地應用于醫(yī)藥、食品、農(nóng)藥等領域,由于其分子結構的特殊性,已被化學合成界默認為萬能的碳鏈延長劑。雙乙烯酮生產(chǎn)過程控制要求嚴格,嚴禁與水、無機酸、堿性物質(zhì)及含鹵素的化合物相接觸,由于其自身結構特殊且非常不穩(wěn)定,所以生產(chǎn)過程中的雙乙烯酮多呈深黑色,并且含有多種雜質(zhì)和乙烯酮的多聚物。另外,由于雙乙烯酮生產(chǎn)以醋酸為原料,生產(chǎn)過程中會副產(chǎn)醋酐等,副產(chǎn)物的性質(zhì)和雙乙烯酮極為相似,從而給中控分析帶來一定的因難。

2.1 水合物儲層的制備

在研發(fā)的水合物鉆采一體化模擬實驗系統(tǒng)中制備CO2水合物儲層。實驗用砂為利用石英砂與北京昌平河砂按照一定比例配比而成的混合砂,其粒度分布曲線如圖4所示。實驗用砂最小砂粒徑0.5 μm,最大砂粒徑450 μm,粒度中值為230 μm,累積填砂量約為500 kg。從整體來看,離真實模擬南海儲層還很遠,但已經(jīng)基本接近日本Nankai Trough的儲層粒度特征。

采用篩管預埋方式開展模擬仿真,根據(jù)圖4所示的模擬沉積物粒度分布曲線,選用標稱擋砂精度為90 um的精密篩管作為井筒控砂介質(zhì),篩管由東營市瑞豐石油技術發(fā)展有限責任公司提供。

CO2水合物成藏過程中維持實驗溫度保持在0~2 ℃,CO2合成過程持續(xù)30 d,通過氣體緩沖腔持續(xù)供氣并維持沉積物內(nèi)部平均壓力在2.5~3.0 MPa,滿足CO2水合物生成的基本溫壓條件,合成結束后根據(jù)耗氣量估算沉積物中的平均水合物飽和度約為40%。詳細的水合物成藏過程及其儲層物性演化規(guī)律將在另文中敘述。

圖4 實驗所用的模擬地層砂粒直徑分布曲線圖

2.2 降壓開采工作制度控制

實際開采過程中通過電潛泵泵頻、井口油嘴聯(lián)合調(diào)節(jié)實現(xiàn)對井底流壓和產(chǎn)氣速率的控制。本實驗中由于不存在電潛泵,因此采用氣路回壓和氣嘴開度相結合的手段,以達到控制模擬井筒內(nèi)部流壓和產(chǎn)氣速率的目的。水合物降壓開采過程中,降壓速率和降壓幅度是影響儲層穩(wěn)定性、出砂規(guī)律的主要因素,進而影響正常的開采產(chǎn)氣過程。筆者首先將?8 mm氣嘴調(diào)至最小開度(10%)以1.1 MPa/h的降壓速率將氣路回壓從2.2 MPa降低值0.2 MPa,此過程中控制儲層圍壓降低至2.4 MPa,一方面始終保持儲層承受的有效圍壓為0.9~1.0 MPa,另一方面防止有效圍壓過大擠毀膠套。

實驗過程中氣路回壓、氣嘴開度及絕對圍壓控制參數(shù)變化曲線如圖5所示。氣路回壓降低到設定值后(0.2 MPa、2 h),調(diào)節(jié)氣嘴開度,觀察冰點附近CO2水合物的分解產(chǎn)出規(guī)律。開采過程中根據(jù)沉積物內(nèi)部外圍孔隙壓力實時調(diào)節(jié)絕對圍壓值,以保證沉積物有效圍壓始終為0.9~1.0 MPa。因此整個開采過程的工作制度控制主要分為“緩慢降壓→穩(wěn)壓產(chǎn)氣”2個階段,以下主要依據(jù)這兩個階段分析冰點附近CO2水合物的分解產(chǎn)出規(guī)律及儲層溫壓演化規(guī)律。

圖5 CO2水合物圍壓控制參數(shù)變化圖

2.3 氣液產(chǎn)出規(guī)律

開采過程中,經(jīng)過兩級三相分離系統(tǒng),產(chǎn)出的氣體通過氣路產(chǎn)出并利用氣體質(zhì)量流量計實時記錄,產(chǎn)出的水砂混合物通過非連續(xù)人工取樣方式采集,然后進行沉淀分離。整個開采過程中的產(chǎn)氣、產(chǎn)水變化規(guī)律分別如圖6、7所示,產(chǎn)砂規(guī)律在此不再贅述。

圖6 冰點附近CO2水合物分解產(chǎn)氣規(guī)律圖

圖7 冰點附近CO2水合物分解產(chǎn)水及水氣比變化規(guī)律圖

由圖7可知,在緩慢降壓階段,產(chǎn)氣過程連續(xù),并伴隨一定的產(chǎn)水量。根據(jù)CO2水合物相平衡條件判斷,該區(qū)間產(chǎn)出的氣體大部分為孔隙中未合成水合物的游離氣,大量游離氣產(chǎn)出過程中攜帶井筒積液產(chǎn)出,該階段的平均水氣比維持在70~120。

在穩(wěn)壓產(chǎn)期階段,CO2產(chǎn)氣速率表現(xiàn)出明顯的波動,CO2水合物的分解產(chǎn)出表現(xiàn)為明顯的非連續(xù)性,大量產(chǎn)氣過程主要集中在3個典型時間區(qū)間(2~8 h、25~33 h、45~50 h),而這3個時間區(qū)間的平均產(chǎn)水速率也明顯增大。特別是在第一大量產(chǎn)氣時間區(qū)間(2~8 h)的前半程(2.0~4.5 h),水氣比急劇上升,然后維持在190~220左右直到第一產(chǎn)氣區(qū)間結束。在第二、第三典型產(chǎn)氣區(qū)間,平均水氣比相比于產(chǎn)氣低谷區(qū)間均表現(xiàn)出明顯的上升,整個降壓開采仿真過程中的最高氣水比大250。在水合物分解低谷區(qū),氣水比下降,最低水氣比約為120,主要是該階段產(chǎn)氣速率較低,氣體攜水產(chǎn)出量也相應降低,但總體而言,穩(wěn)壓產(chǎn)氣階段的綜合氣水比均高于降壓階段,主要是水合物分解水對產(chǎn)水量的貢獻。

2.4 降壓開采過程儲層參數(shù)演化

為了展示水合物開采過程中儲層的整體溫度場、壓力場演化規(guī)律,從而為水合物分解陣面演化規(guī)律提供基礎,選取典型時刻節(jié)點沉積物內(nèi)部的溫度場、壓力場分布規(guī)律分別如圖8、9所示。從圖中可以看出,整個開采過程中儲層壓力呈明顯的非均勻下降趨勢;起始條件下由于成藏模擬過程水合物合成不均勻、外界溫度不均勻等條件的影響,儲層溫度場表現(xiàn)出明顯的“上高下低”,隨著水合物分解過程的持續(xù),儲層整體溫度降低且不同位置的溫度趨于一致,這是由于水合物開采過程中的吸熱效應和儲層傳熱作用共同作用的結果,水合物分解過程儲層最大的溫度降幅為5 ℃。

結合圖6、7的產(chǎn)氣、產(chǎn)水規(guī)律分析,產(chǎn)氣/產(chǎn)水速率越快,水合物分解越大,導致地層溫度下降越快,因此儲層溫度、壓力波動最顯著的3個階段正好是氣水集中產(chǎn)出的3個區(qū)間(2~8 h、25~33 h、45~ 50 h)。

3 結論與展望

1)基于天然氣水合物降壓法開采思路和工藝流程,采用模塊化設計研發(fā)了一套海洋天然氣水合物鉆采一體化模擬實驗系統(tǒng),可模擬實際地質(zhì)條件制備接近海洋水合物儲層的樣品,通過電阻層析成像技術實時探測水合物成藏與分布情況,通過水合物鉆采實驗技術模擬鉆井、降壓開采工藝與過程,實時測量開采過程中出砂與管道流動等過程中多個物理參數(shù)變化,實現(xiàn)了試采全過程實驗模擬。

2)采用該裝置進行了冰點附近CO2水合物的初步開采模擬實驗,基于實驗結果建立了數(shù)據(jù)獲取、數(shù)據(jù)分析的基本流程,通過對開采過程參數(shù)分析,初步獲得了CO2水合物開采過程中儲層溫、壓場變化以及產(chǎn)氣、產(chǎn)水規(guī)律。隨著技術的進步和該模擬實驗系統(tǒng)的不斷完善,其實驗結果將為我國南海天然氣水合物開采提供重要依據(jù)。

圖8 典型時刻儲層壓力場分布規(guī)律圖

圖9 典型時刻儲層溫度場分布規(guī)律圖

猜你喜歡
模擬實驗產(chǎn)氣水合物
天然氣水合物儲運技術研究進展
基于分子模擬的氣體水合物結構特征及儲氣特性研究
濕垃圾與病死豬混合厭氧消化產(chǎn)氣性能研究
海域天然氣水合物三維地震處理關鍵技術應用
氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
羊毒素型產(chǎn)氣莢膜梭菌臨床癥狀及病理變化
加溫加堿預處理對污泥厭氧消化產(chǎn)氣量影響研究
Meso-mechanical model of concrete under a penetration load
模擬實驗在中職學校會計教育中的重要性研究
彈道修正模擬實驗裝置的研究
永胜县| 手机| 昌宁县| 桃源县| 武平县| 德令哈市| 郸城县| 临夏市| 太原市| 宝清县| 郎溪县| 北碚区| 黄龙县| 壶关县| 石楼县| 兰州市| 叙永县| 阆中市| 钟祥市| 泰来县| 隆尧县| 桐梓县| 尤溪县| 定兴县| 峨眉山市| 桂林市| 扬州市| 崇文区| 梓潼县| 阿克| 栾川县| 新营市| 宁安市| 汝南县| 衡东县| 溆浦县| 泸州市| 临桂县| 东方市| 肥东县| 江源县|