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

?

基于水下自主航行器(AUV)的神狐峽谷谷底塊體搬運(yùn)沉積特征及其對(duì)深水峽谷物質(zhì)輸運(yùn)過程的指示

2021-04-16 05:42:14劉錚陳端新朱友生張廣旭董冬冬
關(guān)鍵詞:濁流聲吶塊體

劉錚,陳端新,朱友生,張廣旭,董冬冬

1.中海石油(中國)有限公司深圳分公司,深圳 518054

2.中國科學(xué)院海洋研究所海洋地質(zhì)與環(huán)境重點(diǎn)實(shí)驗(yàn)室,青島 266061

3.青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋地質(zhì)過程與環(huán)境功能實(shí)驗(yàn)室,青島 266071

4.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(廣州),廣州 511458

5.中國科學(xué)院海洋大科學(xué)研究中心,青島 266071

6.中海油田服務(wù)股份有限公司,天津 300459

7.青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室,青島 266061

海底峽谷通常發(fā)育在陸架陸坡區(qū),是陸源物質(zhì)向深海運(yùn)移的重要通道[1]。海底峽谷既可發(fā)育于主動(dòng)大陸邊緣,又可發(fā)育于被動(dòng)大陸邊緣[2]。濁流通常被認(rèn)為是峽谷內(nèi)部物質(zhì)輸運(yùn)的主要營(yíng)力,尤其對(duì)于那些與陸地河流連接的峽谷。Xu等監(jiān)測(cè)到蒙特利峽谷最大的濁流速度可達(dá)2.60 m/s[3];Paull等推算的濁流的前鋒可以達(dá)到7.2 m/s,并且濁流能夠以至少4 m/s的速度搬運(yùn) 800 kg的物體[4]。2006年臺(tái)灣西南恒春地震引起了多處的滑坡和濁流,有序切斷了1 500~4 000 m 水深的 14 條電纜,其中 6 條在峽谷外,估算的流速為3.7~20 m/s[5]。同樣的事情發(fā)生在2009年,莫拉克臺(tái)風(fēng)引起的高密度流/濁流有序切斷了多條海底通訊電纜,計(jì)算的最大速度達(dá)16.6 m/s[6]。

南海神狐峽谷群是垂直陸坡方向發(fā)育的多條近似平行的限制型海底峽谷(圖1a、b)。神狐峽谷群遠(yuǎn)離陸地,更新世以來沉積速率較高(20~34.16 cm/ka)[7],海底發(fā)育了總面積超過 1 000 km2的滑動(dòng)、滑塌、碎屑流等塊體搬運(yùn)沉積(MTDs, Mass Transport Deposits)[8-9]。與同處南海北部陸緣且有大量碎屑物質(zhì)輸入的臺(tái)灣高屏海底峽谷相比,陸坡限定性峽谷內(nèi)濁流發(fā)生的頻率相對(duì)較小[5]。

以往對(duì)深水區(qū)地層結(jié)構(gòu)的探測(cè)主要依賴二維、三維(2D/3D)反射地震資料和淺地層剖面;前者探測(cè)深度幾千米,分辨率較低,在數(shù)米至數(shù)十米之間;后者探測(cè)深度較淺,幾十米到幾百米,容易受到能量吸收和復(fù)雜作業(yè)環(huán)境干擾,設(shè)備分辨率較高,但實(shí)際獲取的數(shù)據(jù)質(zhì)量較差。即便是能量和探測(cè)精度更高的電火花震源,其分辨率也只有約2 m。這對(duì)于研究深水區(qū)海底沉積層的精細(xì)結(jié)構(gòu)是不夠的。AUV探測(cè)是加載各種探測(cè)設(shè)備的自主水下航行器保持與海底幾十米的距離并按照設(shè)定的路由線路進(jìn)行數(shù)據(jù)采集,具有很好的橫向和縱向的分辨率,例如中海油服3 000 m級(jí)AUV攜帶的EdgeTech2200-M淺剖儀的垂向分辨率可達(dá)6~10 cm[10-11]。國際上AUV探測(cè)已經(jīng)越來越多地應(yīng)用到深海沉積研究中[12-14]。國內(nèi)也研發(fā)出多套AUV設(shè)備,但應(yīng)用范圍有限[15]。

圖1 研究區(qū)位置a.神狐峽谷群的位置,b.神狐峽谷群海底地形圖,c.基于三維地震資料的研究區(qū)海底地形圖,d.基于AUV采集的多波束的研究區(qū)海底地形圖。Fig.1 Locations of the study area a.The location of the studied slope confined canyons, b.the bathymetric map of canyons, c.the bathymetric map of the studied canyon segment based on 3D seismic reflection data, d.the bathymetric map of the studied canyon segment based on multi-beam data acquired by AUV.

對(duì)于神狐峽谷群,以往的研究主要集中在滑坡廣泛發(fā)育的崎嶇的峽谷脊部,但對(duì)較為平坦的峽谷谷底研究較少。本文利用AUV獲取的地球物理資料刻畫神狐峽谷谷底沉積特征,對(duì)研究陸坡限定性峽谷的物質(zhì)輸運(yùn)過程具有重要意義。

1 地質(zhì)概況

神狐峽谷群分布在珠江口外海陸坡區(qū),具體位于珠江口盆地的白云凹陷。珠江口盆地所在的南海北部經(jīng)歷古近紀(jì)的裂陷期和新近紀(jì)―第四紀(jì)的裂后期[16]。在約23.8 Ma的漸新世末期,受南海擴(kuò)張中心南向躍遷的影響,陸架坡折帶從白云凹陷的南部退移到目前的位置[17]。神狐峽谷群開始于13.8 Ma的中中新世,經(jīng)歷了4個(gè)期次的北東向遷移發(fā)育,峽谷覆蓋面積從Ⅰ期(13.8~12.5 Ma)、Ⅱ期(12.5~10.5 Ma)到Ⅲ期(10.5~5.5 Ma)逐漸增大,再到Ⅳ期(1.8 Ma~)縮小[18-19]。濁流、底流及兩者之間的相互作用被認(rèn)為影響了峽谷的遷移發(fā)育[20-22]。

2 數(shù)據(jù)和方法

研究使用了三維地震數(shù)據(jù)和AUV資料進(jìn)行對(duì)比分析。三維地震數(shù)據(jù)的時(shí)間采樣間隔2~4 ms,空間采樣間隔 6.25 m×12.5 m,上部 1 500 m 地層主頻約為40 Hz。三維地震數(shù)據(jù)中的海底反射時(shí)間經(jīng)時(shí)深轉(zhuǎn)換生成水深圖(海水聲速取值 1 500 m/s)。AUV 資料使用 Echo Surveyor III(Kongsberg Hugin 1000 AUV)采集,主要加載了多波束、旁掃聲吶和淺地層剖面儀等設(shè)備。多波束為Kongsberg EM2000,聲脈沖頻率平均 2 Hz,掃描寬度為 240 m,水平分辨率可以達(dá)到0.6 m;旁掃聲吶用來反映地形變化和底質(zhì)類型,采用Edgetech全譜旁掃聲吶,頻率105/410 kHz,聲脈沖頻率 3 Hz,脈沖長(zhǎng)度 9 ms/2 ms,掃描范圍221 m/100 m;淺地層剖面儀為 Edgetech 全譜線性調(diào)頻剖面儀,頻率范圍2~16 kHz,實(shí)際工作頻率2~10 kHz,地層分辨率可達(dá) 3~4 cm,聲脈沖頻率3 Hz,記錄長(zhǎng)度 143 ms。作業(yè)時(shí) AUV 在距海底 35 m的水深處以3~4節(jié)的速度航行。2010年在水深1 302.21 m的塊體搬運(yùn)沉積體上實(shí)施了重力柱狀取樣(GC-3),取樣長(zhǎng)度3.6 m?;贏UV獲取的旁掃聲吶、淺地層剖面和多波束等高分辨率地球物理數(shù)據(jù),對(duì)第14條峽谷進(jìn)行了海底地貌和部分谷底的淺地層分析。

3 結(jié)果

3.1 峽谷地形地貌特征

神狐峽谷群水深 200~2 000 m,由 19 條近似平行的峽谷組成。峽谷起源于陸坡的上部,向深水逐漸加寬,最下端匯聚到珠江大峽谷。峽谷長(zhǎng)3.6~36 km,寬 1~5 km,深 100~400 m。峽谷脊部地形崎嶇不平,溝壑陡崖普遍發(fā)育,峽谷谷底較為平坦(圖1)?;谌S反射地震可以在峽谷谷底觀測(cè)到巨大的陡崖和大型的溝槽,可以大體分辨邊界明顯的MTDs(圖1c)?;贏UV獲取的旁掃聲吶圖則更加清晰地展示了所選取峽谷谷底的精細(xì)地貌特征(圖1b,圖2a)。MTDs在谷底廣泛分布,總面積近 17 km2。MTDs在峽谷頭部的弧形斜坡和谷底都有分布,MTDs呈不規(guī)則的圓形展布,其表面起伏不平(圖2b);在峽谷中下游MTDs沿著谷底狹長(zhǎng)條帶狀分布(圖2c-e)?;贏UV的多波束數(shù)據(jù)(圖1d)和旁掃聲吶數(shù)據(jù)(圖2d),選取了峽谷下游水深1 280~1 360 m、距離峽谷出口 10~15 km 的 MTDs典型發(fā)育區(qū)進(jìn)行了詳細(xì)分析。3個(gè)MTDs分別位于不同臺(tái)階上(圖3,圖4)。MTDs的上面發(fā)育了數(shù)量眾多的寬 3~10 m、長(zhǎng) 20~700 m、深 5~20 cm 的沖蝕溝槽(圖3)。MTDs的厚度都在8.4 m及以下,這在常規(guī)2D/3D地震資料上是難以分辨的(圖4)。

3.2 MTD1

MTD1位于谷底的最西邊,是峽谷谷底主水道的位置,中間寬、上下兩端窄,在下部受到突出的水道堤岸的阻擋而改變方向且寬度變窄。MTD1的上部可以識(shí)別出狹長(zhǎng)的渠道狀的物質(zhì)輸運(yùn)后的殘留物(圖3),而在MTD1鄰近的峽谷翼部沒有識(shí)別出明顯的滑坡體或滑坡遺跡。MTD1長(zhǎng)5.00 km,頭尾兩端厚度小,中間厚度大,厚度最小1.30 m,最大8.40 m,面積 0.90 km2,體積約 4.37 km3(表1)。MTD1可以清晰識(shí)別出連續(xù)的底邊界反射指示底部剪切面,其反射強(qiáng)度低于海底反射但明顯高于其他正常地層反射。MTD1內(nèi)部為振幅較弱的雜亂反射(圖5),僅在底部的局部位置有少量的長(zhǎng)度有限的連續(xù)反射軸(圖5c)。沉積體上發(fā)育了大量平直的沖蝕溝槽,尤其在主水道上更為低洼的地方(圖3)。基于AUV的旁掃聲吶圖(淺色)指示整條峽谷谷底淺表層主要由軟的沉積物組成(圖2);重力柱狀取樣結(jié)果證實(shí)MTD1沉積體由很軟的高可塑性的灰綠色粉砂組成,含非常少的砂(圖2f——g)。

圖2 峽谷的旁掃聲吶圖和巖心圖a.研究區(qū)所在峽谷的旁掃聲吶圖,b——e.局部放大的旁掃聲吶圖,f——g.MTD1上GC-3站位的部分重力柱狀巖心樣品。虛線多邊形指示MTDs的位置。Fig.2 The side-scan sonar map of the canyon floor and sections of the gravity core acquired from the canyon floor a.The side-scan sonar map of the studied canyon, b-e.amplified side-scan sonar maps of the study area, f-g.sediment samples from the gravity core GC-3 over MTD1 at 1.25 m and 3.25 m, respectively.The dashed polygons denote locations of MTDs on the canyon floor.

圖3 研究區(qū)3個(gè)典型MTDs的地形立體圖及穿過MTDs的地形變化曲線Fig.3 Stereo views of three MTDs and bathymetric curves crossing these MTDs

圖4 過同一峽谷位置的三維反射地震剖面和AUV淺地層剖面a.反射地震剖面,b.AUV淺地層剖面。Fig.4 The seismic profile exttracted from the 3D seismic data and sub-bottom profiles acquired by AUV crossing the same section of the canyon a.The seismic profile, b.The AUV based sub-bottom profile.

表1 峽谷谷底MTDs的幾何參數(shù)Table 1 Geometric parameters of MTDs on the canyon floor

3.3 MTD2

MTD2位于峽谷谷底中部的臺(tái)階上,平面上呈不規(guī)則梯形展布(圖1d,圖2d,圖3)。MTD2并非直接出露在海底,上部覆蓋了一層厚約0.6 m的沉積層。MTD2寬度變化不大,為260.0~350.0 m;長(zhǎng)750.0 m;厚 0.90~3.20 m;面積較小,為 0.28 km2;體積約 0.58 km3(表1)。MTD2 沉積體上面也發(fā)育有沖蝕溝槽,此外在沉積體頭部發(fā)育了更多的沖溝。MTD2同樣可以識(shí)別出強(qiáng)反射的底邊界和弱的頂界面。底邊界并非平直光滑,而是出現(xiàn)了大量的凹槽,這些凹槽寬 5.0~20.0 m,高 0.375~0.75 m,其寬度明顯大于沉積體上面的沖蝕溝槽,推測(cè)是由塊體搬運(yùn)沉積體在運(yùn)動(dòng)過程中侵蝕海底形成的侵蝕溝。不同于MTD1兩端較小的厚度,MTD2從頭部到尾部厚度逐漸減?。▓D6)。

3.4 MTD3

MTD3位于谷底最東側(cè)的臺(tái)階上。頭部寬度小,中部和尾部寬度大。MTD3 長(zhǎng) 2.0 km,寬 70~600 m,厚 1.2~3.0 m,面積 1.10 km2,體積約 2.31 km3。垂直于峽谷走向,沉積體呈楔形展布,整個(gè)沉積體分布在中間低兩側(cè)高的大型凹槽中(圖3,圖7a——b)。MTD3在大型凹槽東側(cè)寬度大于其西側(cè)寬度,且東側(cè)厚度逐漸變薄。沿著峽谷走向,沉積體厚度在中部最大,頭部稍小,尾部最薄。MTD3底部反射較為平滑,反射強(qiáng)度較強(qiáng),但明顯小于MTD1和MTD2的底部反射強(qiáng)度。

圖5 過 MTD1 的 AUV 淺地層剖面圖a.上部橫剖面,b.中部橫剖面,c.下部橫剖面。Fig.5 AUV based sub-bottom profiles across MTD1 a.The upper profile perpendicular to the strike, b.the middle profile perpendicular to the strike, c.the lower profile perpendicular to the strike.

圖6 過 MTD2 的 AUV 淺地層剖面圖a.橫剖面,b.縱剖面。Fig.6 AUV based sub-bottom profiles across MTD2 a.The profile perpendicular to the strike, b.the profile parallel to the strike.

圖7 過 MTD3 的 AUV 淺地層剖面圖a.上部橫剖面,b.下部橫剖面,c.縱剖面。Fig.7 AUV based sub-bottom profiles across MTD3 a.The upper profile perpendicular to the strike, b.the lower profile perpendicular to the strike, c.the profile parallel to the strike.

4 討論

4.1 峽谷谷底MTDs的形成

通過對(duì)MTDs埋深的分析,推斷MTDs形成的時(shí)間距今不遠(yuǎn),因?yàn)樵贛TD1和MTD3上面沒有識(shí)別出正常沉積的地層。MTD2的上部發(fā)育厚約0.60 m的沉積地層(圖6),明顯晚于MTD1和MTD2。

基于AUV的地球物理資料對(duì)MTDs特征的分析結(jié)果,我們推測(cè)峽谷中下游大多數(shù)MTDs并不是直接從鄰近的峽谷脊部上搬運(yùn)下來,而是從上游通過滑塌-碎屑流的形式運(yùn)移下來。首先,除了峽谷頭部弧形斜坡分布有大量不規(guī)則圓形展布的MTDs,峽谷中下游的MTDs大都沿峽谷走向狹長(zhǎng)展布,并不是垂直峽谷走向展布(圖2)。雖然脊部沉積物滑移后可以繼續(xù)通過谷底向下繼續(xù)運(yùn)移,但在谷底兩側(cè)靠近坡腳的位置并沒有發(fā)現(xiàn)較厚塊體搬運(yùn)沉積體的存在,在鄰近的脊部也沒有發(fā)現(xiàn)物質(zhì)滑移的證據(jù),相反,我們?cè)贛TD1上部發(fā)現(xiàn)了明顯的物質(zhì)輸運(yùn)后的殘留遺跡(圖5a)。MTD2具有較小的長(zhǎng)寬比以及逐漸減小的厚度,說明該沉積體搬運(yùn)的距離并不遠(yuǎn)(表1,圖6),推測(cè)其是從鄰近的谷底陡坡上搬運(yùn)下來。MTD3 沉積體在大型凹槽東側(cè)寬度要大于西側(cè),且東側(cè)的沉積厚度向東逐漸變薄,說明MTD3不是直接來源于峽谷東側(cè)的脊部,這表明峽谷上游比中下游發(fā)生滑坡的概率更大。三維反射地震資料的解釋也表明MTDs在峽谷上游的數(shù)量更多,分布面積更廣(圖2a)[8-9]。雖然三維反射地震資料也顯示峽谷中下游的脊部有大量MTDs,但這些MTDs厚度較大,在幾十米以上,很多MTDs的內(nèi)部地層變形不嚴(yán)重,反映了塊體滑移或程度較弱的滑塌,但不能反映當(dāng)前較小時(shí)間尺度內(nèi)峽谷谷底的沉積過程。相比較而言,基于AUV的高分辨地球物理資料更真實(shí)反映了峽谷內(nèi)部的沉積特征和過程,即當(dāng)前或較短時(shí)間內(nèi)峽谷谷底分布著大量小規(guī)模的沿峽谷走向呈狹長(zhǎng)展布的MTDs。

研究區(qū)聲學(xué)剖面上的塊體搬運(yùn)沉積體可能是一次塊體流沉積的結(jié)果,也可能是多次塊體流沉積的結(jié)果,但塊體流間隔的時(shí)間較短,沒有形成聲學(xué)可識(shí)別的正常沉積地層。MTD1所處的谷底主水道位置更容易匯聚長(zhǎng)距離搬運(yùn)下來的碎屑流沉積。

4.2 峽谷谷底的物質(zhì)輸運(yùn)

神狐峽谷群距離陸地近250 km,河流高懸浮物注入引起的濁流對(duì)峽谷影響很小,在MTD1上GC-3站位細(xì)粒粉砂沉積物也指示缺少陸源粗粒沉積物的輸入。前人的研究表明峽谷坡度較陡的脊部發(fā)育了大量的MTDs,形成了起伏不平的地貌[23],谷底地形相對(duì)比較平坦。通過以上對(duì)AUV獲取的高分辨率地形和淺部地層資料的解析,我們認(rèn)識(shí)到峽谷谷底并不像它的地表那么簡(jiǎn)單,而是在平坦的地表下發(fā)育了MTDs。這些MTDs的面積和厚度遠(yuǎn)沒有峽谷脊部的大,但在峽谷谷底大量分布。峽谷中下游谷底的大多數(shù)MTDs并不是直接來源于峽谷脊部,而是來自MTDs的上部。這指示了峽谷谷底的物質(zhì)輸送的一種重要途徑可能是通過不斷重復(fù)的發(fā)生塊體搬運(yùn)沉積的形式進(jìn)行的?;虏⒁饾崃鞅徽J(rèn)為是陸坡限定性峽谷物質(zhì)輸運(yùn)的主要營(yíng)力,這種情況下濁流的產(chǎn)生通常是伴隨著滑塌-碎屑流的產(chǎn)生而產(chǎn)生。在高海平面的情況下,以紊流為支撐力的濁流會(huì)侵蝕并攜帶部分沉積物流向下坡方向輸運(yùn),但很難將垂直峽谷走向的MTDs完全改造成沿峽谷走向展布而不留下痕跡,除非發(fā)生大規(guī)模的滑坡事件。但高分辨率AUV資料表明目前的峽谷中下游的谷底只是分布著大量的小規(guī)模MTDs,而不是大量的濁流沉積體。據(jù)此,我們認(rèn)為AUV可辨識(shí)的較小的時(shí)間尺度范圍內(nèi),在特定沉積環(huán)境下(比如高海平面時(shí)期),神狐峽谷內(nèi)物質(zhì)輸運(yùn)的主要營(yíng)力不只是濁流,還應(yīng)該考慮峽谷谷底不斷重復(fù)進(jìn)行的塊體搬運(yùn)沉積過程。AUV資料和三維反射地震資料的解釋結(jié)果并不沖突,只是存在空間尺度和時(shí)間尺度的不同。新形成的MTDs的表面通常會(huì)有沉積物堆積造成的凹凸不平,但谷底MTDs的表面只發(fā)育了大量小型沖蝕溝槽,這很可能是峽谷內(nèi)部的較強(qiáng)的水動(dòng)力對(duì)MTDs表面進(jìn)行了改造。前人研究表明,神狐峽谷內(nèi)部具有復(fù)雜的海洋水動(dòng)力環(huán)境,其中內(nèi)潮普遍發(fā)育且能量最大,峽谷內(nèi)部發(fā)育的內(nèi)潮最大流速可達(dá)50 cm/s,沿著峽谷軸向往復(fù)運(yùn)動(dòng)[24]。MTDs沉積體上面的沖溝是由潮流、余流還是其他類型的水體運(yùn)動(dòng)造成的還不明確。

5 結(jié)論

(1)基于AUV的高分辨率多波束,旁掃聲吶和淺地層剖面數(shù)據(jù)對(duì)峽谷整體和局部的MTDs進(jìn)行了精細(xì)刻畫,在峽谷谷底識(shí)別了常規(guī)地球物理資料不能辨識(shí)或不能清晰辨識(shí)的大量MTDs的分布。峽谷上游弧形斜坡和谷底分布著大量不規(guī)則圓形展布的MTDs,峽谷中下游的MTDs多在谷底呈狹長(zhǎng)展布。

(2)對(duì)MTDs典型發(fā)育區(qū)的研究表明,MTDs沉積體的厚度在8.4 m及以下,沉積體內(nèi)部為反射強(qiáng)度較弱的雜亂反射,推斷為經(jīng)過一定搬運(yùn)距離而充分混合的碎屑流。

(3)通過分析MTDs沉積體的形態(tài)、沉積厚度變化并結(jié)合谷底兩側(cè)峽谷脊部的地層反射特征,認(rèn)為谷底分布的MTDs主要來源于其上部,而不是鄰近的峽谷脊部。在特定沉積環(huán)境和較小時(shí)間尺度范圍內(nèi),除了濁流外,從峽谷上游到下游的塊體搬運(yùn)沉積過程的重復(fù)發(fā)生很可能是峽谷物質(zhì)輸運(yùn)的另一種重要形式。

猜你喜歡
濁流聲吶塊體
探索大洋的“千里眼”——聲吶
海底峽谷內(nèi)濁流流動(dòng)與沉積特征數(shù)值模擬研究
一種新型單層人工塊體Crablock 的工程應(yīng)用
一種便攜式側(cè)掃聲吶舷側(cè)支架的設(shè)計(jì)及實(shí)現(xiàn)
聲吶
“問走”霧霾濁流 “詢來”藍(lán)天碧水
一種Zr 基塊體金屬玻璃的納米壓入蠕變行為研究
上海金屬(2015年3期)2015-11-29 01:09:58
塊體非晶合金及其應(yīng)用
波浪作用下斜坡上護(hù)面塊體斷裂破壞的數(shù)值模擬
水道港口(2015年1期)2015-02-06 01:25:35
COTS技術(shù)在聲吶裝備中的應(yīng)用
舞钢市| 纳雍县| 手机| 札达县| 道孚县| 锦州市| 邹平县| 峨眉山市| 麦盖提县| 杭锦后旗| 南汇区| 大余县| 万载县| 徐闻县| 晋江市| 潢川县| 连州市| 南澳县| 珠海市| 阳春市| 息烽县| 旺苍县| 宜章县| 广西| 甘孜县| 将乐县| 安康市| 伊川县| 晴隆县| 卓资县| 河南省| 延安市| 九龙坡区| 阜城县| 温州市| 灌阳县| 那坡县| 陆丰市| 闽侯县| 海盐县| 藁城市|