• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于隨機(jī)森林的川藏鐵路沿線縣域地質(zhì)環(huán)境承載力評價(jià)

    2022-09-29 08:04:44支澤民陳瓊周強(qiáng)夏興生劉峰貴
    科學(xué)技術(shù)與工程 2022年24期
    關(guān)鍵詞:鐵路沿線凍土縣域

    支澤民,陳瓊,周強(qiáng),夏興生,劉峰貴,2*

    (1.青海師范大學(xué)地理科學(xué)學(xué)院,西寧 810008;2.高原科學(xué)與可持續(xù)發(fā)展研究院,西寧 810008)

    隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展,人類活動(dòng)對于自然環(huán)境的改造作用逐漸加劇,對于地表形態(tài)的改造趨勢增強(qiáng),社會(huì)發(fā)展與自然環(huán)境之間的矛盾日益突出,逐漸成為區(qū)域社會(huì)發(fā)展的重要制約因素[1]。為科學(xué)合理對區(qū)域開發(fā)進(jìn)行指導(dǎo),對區(qū)域內(nèi)資源環(huán)境承載力的合理評估就變得極為必要,而資源環(huán)境承載力包括資源承載力與環(huán)境承載力,其中環(huán)境承載力又包括水環(huán)境承載力[2-3]、土壤環(huán)境承載力[4]以及地質(zhì)環(huán)境承載力[5-6]等,其中又以地質(zhì)環(huán)境承載力為承載力之根本。作為人類社會(huì)經(jīng)濟(jì)活動(dòng)的基礎(chǔ),地質(zhì)環(huán)境的好壞直接關(guān)系到人類社會(huì)經(jīng)濟(jì)的發(fā)展方向與發(fā)展速度。

    目前,關(guān)于地質(zhì)環(huán)境承載力方面的研究較為成熟,主要包括區(qū)域地質(zhì)環(huán)境承載力[7-8]、城市地質(zhì)環(huán)境承載力[9-10]以及礦山地質(zhì)環(huán)境承載力[11-12]等,多通過構(gòu)建地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系,對指標(biāo)進(jìn)行合理量化,從而實(shí)現(xiàn)地質(zhì)環(huán)境承載力的評價(jià)。汪宙峰等[8]選取都汶公路沿線區(qū)域,通過層次分析法確定指標(biāo)權(quán)重,對公路沿線區(qū)生態(tài)地質(zhì)環(huán)境承載力進(jìn)行評價(jià),揭示了區(qū)域生態(tài)地質(zhì)環(huán)境承載力的影響因素與空間分布;楊盼盼等[9]選取2010年、2013年以及2017年臨潼區(qū)23個(gè)街道的面板數(shù)據(jù)構(gòu)建了行政區(qū)劃級別的城鎮(zhèn)化—地質(zhì)環(huán)境耦合協(xié)調(diào)度概念模型,對地質(zhì)環(huán)境的時(shí)空演化規(guī)律進(jìn)行了分析;何芳等[12]基于不同時(shí)期的高分影像從木里礦區(qū)的地質(zhì)環(huán)境背景、地質(zhì)環(huán)境破壞以及地質(zhì)環(huán)境恢復(fù)等多個(gè)層面對礦區(qū)地質(zhì)環(huán)境承載力進(jìn)行評價(jià),揭示了礦區(qū)地形地貌以及土地資源等破壞狀況,為礦區(qū)修復(fù)治理提供了可靠的基礎(chǔ)資料。盡管目前對于地質(zhì)環(huán)境承載力的研究較多,但多集中于中東部區(qū)域,針對西部高原地區(qū)的相關(guān)研究較少。且受高原地區(qū)凍土分布的影響,地區(qū)地質(zhì)環(huán)境承載力評價(jià)體系有待修正,很多研究在進(jìn)行地質(zhì)環(huán)境承載力時(shí)未能將其考慮在內(nèi),使得高原地區(qū)地質(zhì)環(huán)境承載力的評價(jià)存在缺陷。為此,在綜合其他區(qū)域地質(zhì)環(huán)境承載力指標(biāo)的同時(shí),加入凍土分布這一極具區(qū)域特色的評價(jià)指標(biāo)對川藏鐵路沿線縣域這一重要區(qū)域進(jìn)行評價(jià)。

    川藏鐵路的建成對于西藏地區(qū)社會(huì)經(jīng)濟(jì)的崛起具有極其重要的戰(zhàn)略地位[13],其沿線地質(zhì)地貌環(huán)境復(fù)雜,人類活動(dòng)日漸加劇,而目前關(guān)于川藏鐵路沿線地質(zhì)環(huán)境承載力的研究較少。因此,對于川藏鐵路沿線地質(zhì)環(huán)境承載力的評估就變得及其重要。選取川藏鐵路沿線縣作為研究區(qū),基于隨機(jī)森林模型,構(gòu)建包含地貌環(huán)境、地質(zhì)環(huán)境、生態(tài)環(huán)境以及社會(huì)環(huán)境4個(gè)準(zhǔn)則層12個(gè)指標(biāo)的川藏鐵路沿線縣域地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系,受該區(qū)域主要位于青藏高原這一特殊地理位置的影響,加入凍土的空間分布作為評價(jià)指標(biāo),該指標(biāo)的選取可很好地契合本區(qū)域?qū)嶋H地質(zhì)環(huán)境。同時(shí),結(jié)合地理信息系統(tǒng)(geographic information system,GIS)空間分析對川藏鐵路沿線縣域地質(zhì)環(huán)境承載力進(jìn)行評價(jià),以期獲得川藏鐵路沿線縣域的地質(zhì)環(huán)境承載力空間分布狀況,明析其地貌、地質(zhì)以及社會(huì)等多環(huán)境質(zhì)量空間分布,從而為川藏鐵路沿線周圍縣域的社會(huì)經(jīng)濟(jì)發(fā)展與災(zāi)害防治提供參考依據(jù)。

    1 研究區(qū)

    川藏鐵路自東由四川省成都起,向西經(jīng)過雅安、康定、昌都、林芝、山南,至西藏拉薩結(jié)束,全長1 543 km[13],是一條由橫斷山高山峽谷區(qū)向青藏高原內(nèi)部延伸的鐵路干線,橫跨揚(yáng)子板塊、川滇地塊、羌塘地塊和拉薩地塊等大地構(gòu)造單元[14],翻越4 000 m海拔以上高山8座,跨過7條大江大河,分布有龍門山斷裂、鮮水河斷裂、玉農(nóng)希斷裂等眾多區(qū)域性活動(dòng)斷裂[15],地質(zhì)、地貌環(huán)境復(fù)雜(圖1),活動(dòng)斷裂密集,地層活動(dòng)較為活躍,地震烈度較為強(qiáng)烈,沿途自然災(zāi)害頻發(fā),尤其是地質(zhì)災(zāi)害[14,16-18],且隨著鐵路線自東向西向高原內(nèi)部的延伸,凍融侵蝕風(fēng)險(xiǎn)極大,以未完工林芝—雅安段較為嚴(yán)重[19]。同時(shí),隨著中國西部大開發(fā)戰(zhàn)略的推進(jìn),川藏鐵路沿線的人類活動(dòng)強(qiáng)度更是逐年增大,城鎮(zhèn)用地?cái)U(kuò)張顯著[20]。在自然環(huán)境與人類活動(dòng)的共同作用下,川藏鐵路沿線脆弱的地質(zhì)環(huán)境極易遭到不可逆的破壞。因此,對于川藏鐵路沿線縣域地質(zhì)環(huán)境承載力的研究就顯得極為重要。選取自雅安市出發(fā)的川藏鐵路經(jīng)過的33個(gè)縣區(qū)作為研究區(qū),對其沿線縣域的地質(zhì)環(huán)境承載力進(jìn)行評價(jià)。

    圖1 研究區(qū)概況圖Fig.1 Overview of the study area

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

    2.1 數(shù)據(jù)來源

    數(shù)據(jù)主要以地質(zhì)、地貌、凍土分布、降水、地質(zhì)災(zāi)害數(shù)據(jù)等數(shù)據(jù)為主,以基礎(chǔ)地理信息數(shù)據(jù)為輔,通過野外調(diào)查與網(wǎng)絡(luò)收集等多渠道、多方式對數(shù)據(jù)進(jìn)行獲取。其中,地質(zhì)災(zāi)害點(diǎn)數(shù)據(jù)來自野外調(diào)查與網(wǎng)上資料收集,同時(shí),筆者參加第二次青藏高原綜合科學(xué)考察項(xiàng)目對調(diào)查數(shù)據(jù)進(jìn)行部分點(diǎn)的驗(yàn)證,主要包括八宿縣、卡若區(qū)、江達(dá)縣、類烏齊縣等地區(qū)大型災(zāi)害點(diǎn);工程巖組數(shù)據(jù)、斷層數(shù)據(jù)來源于中國地質(zhì)調(diào)查局1∶2 500 000中國地質(zhì)圖(https://www.cgs.gov.cn/),矢量數(shù)據(jù),為使巖性數(shù)據(jù)可帶入指標(biāo)體系中,參考已有相關(guān)研究[21-23]對工程巖組進(jìn)行賦值劃分(表1);海拔高程數(shù)據(jù)分辨率為30 m,來源于地理空間數(shù)據(jù)云平臺(tái)(http://www.gscloud.cn/home)的ASTER GDEM 30 m數(shù)據(jù),坡度與地形起伏度數(shù)據(jù)基于高程數(shù)據(jù)通過ArcGIS10.7平臺(tái)提取獲得;年降水量數(shù)據(jù)為2017—2019年全國降水量數(shù)據(jù)所取平均值,柵格數(shù)據(jù),分辨率為1 km,來源于國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心(http://www.geodata.cn/);凍土及凍土類型分布數(shù)據(jù)來源于青藏高原國家科學(xué)數(shù)據(jù)中心(http://data.tpdc.ac.cn/zh-hans/),柵格數(shù)據(jù),分辨率為1 km;人口密度數(shù)據(jù)來源世界人口數(shù)據(jù)網(wǎng)站world pop(https://www.worldpop.org/),柵格數(shù)據(jù),分辨率為100 m;歸一化植被指數(shù)(normalized difference vegetation index,NDVI)數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(https://www.resdc.cn/),柵格數(shù)據(jù),分辨率為1 km;各級行政區(qū)劃、河流等數(shù)據(jù)以及道路密度與居民點(diǎn)密度數(shù)據(jù)來源于全國地理信息資源服務(wù)系統(tǒng)(https://www.webmap.cn/main.do?method=index),矢量數(shù)據(jù),通過ArcGIS10.7進(jìn)行密度分析獲取。其中,地質(zhì)地貌數(shù)據(jù)、降水?dāng)?shù)據(jù)、凍土類型分布數(shù)據(jù)以及人口、道路等數(shù)據(jù)均為構(gòu)建川藏鐵路沿線地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系的指標(biāo)數(shù)據(jù),行政區(qū)劃數(shù)據(jù)用以確定研究區(qū)范圍。

    表1 巖性賦值表[21-23]Table 1 The assignment Table of lithology

    2.2 研究方法

    主要以隨機(jī)森林模型(random forest,RF)對指標(biāo)層數(shù)據(jù)進(jìn)行權(quán)重確定,屬于組合式的深度學(xué)習(xí)方法[24],作為一種較為常用的機(jī)器學(xué)習(xí)方法,是一種分類回歸樹的數(shù)據(jù)挖掘方法[25],該模型結(jié)合了回歸樹、隨機(jī)選擇和裝袋的算法思想[26],在解決非線性問題上具有較大的優(yōu)勢[27]。同時(shí),該方法對于運(yùn)算數(shù)據(jù)的要求較低,對于多要素類的問題研究可定量提供解釋變量對于被解釋變量的貢獻(xiàn)率與重要性程度,即可確定其權(quán)重。因此,對于構(gòu)建指標(biāo)體系的相關(guān)研究來說,較為合適,且對于被解釋變量的影響因素的分析也具有一定的便利性,準(zhǔn)確性相對更高更具有客觀性[28]。與其他機(jī)器學(xué)習(xí)算法相比,隨機(jī)森林模型對于模型運(yùn)算中產(chǎn)生的噪聲與異常值具有較高的容忍度[24,29],應(yīng)用十分廣泛,包括生態(tài)風(fēng)險(xiǎn)[27]、災(zāi)害易發(fā)性[28]、金融相關(guān)[30]、災(zāi)害風(fēng)險(xiǎn)[31]等研究。主要使用R 4.0.4環(huán)境下的R studio編譯器調(diào)用RF軟件包進(jìn)行迭代運(yùn)算,其原理如圖2[28]所示。

    CART(classification and regression tree)為分類回歸樹;OOB(out of bag)為袋外誤差圖2 隨機(jī)森林模型原理圖[28]Fig.2 Schematic diagram of random forest model[28]

    3 指標(biāo)體系構(gòu)建

    3.1 評價(jià)指標(biāo)構(gòu)建

    地質(zhì)環(huán)境承載力評價(jià)指標(biāo)的選取對于研究結(jié)果的準(zhǔn)確度至關(guān)重要,其目標(biāo)是協(xié)調(diào)人地之間的關(guān)系,合理的指標(biāo)選取可使評價(jià)結(jié)果更接近實(shí)際狀況。參考文獻(xiàn)[32-35],根據(jù)研究區(qū)實(shí)際區(qū)域環(huán)境特征,綜合考慮數(shù)據(jù)的可獲取性,構(gòu)建川藏鐵路沿線地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系。從川藏鐵路沿線各縣的地形地貌條件、地質(zhì)環(huán)境條件、社會(huì)經(jīng)濟(jì)條件以及生態(tài)環(huán)境4個(gè)準(zhǔn)則層分別選取海拔、坡度、地形起伏度、斷層密度、工程巖組、凍土分布、歷史災(zāi)害點(diǎn)密度、植被覆蓋指數(shù)、年均降水量、人口密度、居民點(diǎn)密度、道路密度作為評價(jià)指標(biāo),構(gòu)建川藏鐵路沿線各縣地質(zhì)環(huán)境質(zhì)量評價(jià)指標(biāo)體系(圖3)。

    圖3 地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系Fig.3 Geo-environmental carrying capacity evaluation system

    作為制約區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的重要因素,地質(zhì)環(huán)境承載力的大小受到多種要素的影響。文中地貌環(huán)境利用海拔、坡度以及地形起伏度表征,其與地質(zhì)環(huán)境承載力的大小呈負(fù)相關(guān)關(guān)系。地質(zhì)環(huán)境包括斷層密度、工程巖組、凍土的分布以及歷史災(zāi)害點(diǎn)密度,其中歷史地質(zhì)災(zāi)害點(diǎn)對于地質(zhì)環(huán)境質(zhì)量、地質(zhì)環(huán)境承載力具有極好的表示意義,當(dāng)某一地區(qū)地質(zhì)災(zāi)害發(fā)生較多時(shí),則表明該區(qū)域地質(zhì)環(huán)境質(zhì)量較差,地質(zhì)環(huán)境承載力較低。生態(tài)環(huán)境利用歸一化植被覆蓋指數(shù)NDVI和年均降水量表示,川藏鐵路沿線區(qū)域海拔較高,除林芝、波密等縣外,其余地區(qū)均降水量均較少,因此,降水量成為區(qū)域生態(tài)環(huán)境變化的限制因素。社會(huì)環(huán)境指標(biāo)包括人口密度、居民點(diǎn)密度和道路密度,隨著區(qū)域內(nèi)社會(huì)經(jīng)濟(jì)的快速發(fā)展,人類活動(dòng)對于區(qū)域地質(zhì)環(huán)境承載力的占用逐漸增大,而人口、建筑以及道路等指標(biāo)可較好地的表征人類生產(chǎn)活動(dòng)的強(qiáng)度(圖4)。

    圖4 地質(zhì)環(huán)境承載力評價(jià)指標(biāo)空間分布Fig.4 Spatial distribution of Geo-environmental carrying capacity evaluation indicators

    3.2 指標(biāo)分級

    不同指標(biāo)對于地質(zhì)環(huán)境承載力的貢獻(xiàn)存在差異,評價(jià)指標(biāo)中既存在數(shù)值量,又存在類型量,各數(shù)值指標(biāo)之間存在較大差異,為去除數(shù)值量與類型量指標(biāo)之間的差異,且為避免對凍土分布和工程巖組巖性等類型量指標(biāo)的人為賦值所造成的實(shí)驗(yàn)誤差,體現(xiàn)計(jì)算的科學(xué)性,在進(jìn)行計(jì)算前對數(shù)據(jù)進(jìn)行了預(yù)處理,進(jìn)行等級劃分。同時(shí)為體現(xiàn)指標(biāo)分級的客觀性,數(shù)值量指標(biāo)采用“自然裂點(diǎn)法”將各個(gè)指標(biāo)劃分為5級,分別為 “差”“較差”“中等”“較好”以及“好”,等級“差”賦值為1,“較差”賦值為2,“中等”賦值為3,“較好”賦值為4,“好”賦值為5,以便進(jìn)行下一步計(jì)算。類型量指標(biāo)主要包括工程巖組硬度和凍土的分布,其中工程巖組硬度參照文獻(xiàn)[21-23]關(guān)于工程巖組巖性硬度的劃分將其劃分為5級(表1),堅(jiān)硬巖類、較堅(jiān)硬巖類、較堅(jiān)硬-較軟弱巖類、較軟巖類、軟巖類。凍土分布數(shù)據(jù)則依據(jù)凍土類型特點(diǎn)將未分布凍土等級劃分為“好”,多年凍土劃分為“中等”,季節(jié)性凍土劃分為“差”3個(gè)等級。其中,海拔、坡度、地形起伏度、斷層密度、歷史災(zāi)害點(diǎn)密度、年均降水量、人口密度、聚居民點(diǎn)密度、道路密度為負(fù)向指標(biāo),年均降水量與NDVI為正向指標(biāo),據(jù)此,等級劃分如表2所示。

    3.3 指標(biāo)權(quán)重確定及評價(jià)

    首先,對各個(gè)指標(biāo)數(shù)據(jù)進(jìn)行處理,將不同分辨率的柵格數(shù)據(jù)統(tǒng)一重采樣為1 km分辨率,將矢量數(shù)據(jù)通過“密度分析”工具轉(zhuǎn)化為1 km分辨率柵格數(shù)據(jù),后依據(jù)評價(jià)指標(biāo)分級標(biāo)準(zhǔn)(表2)對各指標(biāo)數(shù)據(jù)進(jìn)行重分類,并分別賦值。同時(shí),考慮到研究區(qū)域面積較大,選取常用的正方形格網(wǎng)法進(jìn)行地質(zhì)環(huán)境承載力評價(jià)。創(chuàng)建分辨率為1 km的格網(wǎng)對研究區(qū)地質(zhì)環(huán)境承載力進(jìn)行評價(jià),將其劃分為177 107個(gè)評價(jià)單元,并對每個(gè)格網(wǎng)對應(yīng)各指標(biāo)的類型賦值進(jìn)行提取。利用ArcGIS“加權(quán)求和”工具將各準(zhǔn)則層包含指標(biāo)進(jìn)行加權(quán)求和作為隨機(jī)森林模型中的目標(biāo)變量,將各指標(biāo)作為預(yù)測變量。隨機(jī)選取70%的評價(jià)單元代入隨機(jī)森林模型進(jìn)行運(yùn)算,30%的評價(jià)單元進(jìn)行模型驗(yàn)證。隨機(jī)森林模型的構(gòu)建中,回歸樹的數(shù)量ntree與樹節(jié)點(diǎn)mtry作為隨機(jī)森林的重要參數(shù),通過設(shè)置不同的ntree與mtry值,從而計(jì)算出不同隨機(jī)特征數(shù)下的隨機(jī)森林的袋外誤差OOB[28],OOB袋外誤差越小,模型預(yù)測精度越高[27],而后通過%IncMSE值進(jìn)行重要性分析,基于其重要性程度確定指標(biāo)層權(quán)重(表3)。經(jīng)過多次重復(fù)調(diào)試,地貌環(huán)境和地質(zhì)環(huán)境準(zhǔn)則層設(shè)置ntree=800,mtry=5,生態(tài)環(huán)境與社會(huì)經(jīng)濟(jì)環(huán)境準(zhǔn)則層ntree=900,mtry=5,模型整體方差解釋率R2均大于0.80,誤差Error均小于0.01。

    表2 評價(jià)指標(biāo)分級標(biāo)準(zhǔn)Table 2 Evaluation index grading criteria

    表3 評價(jià)指標(biāo)權(quán)重Table 3 Weights of evaluation index

    確定各指標(biāo)及準(zhǔn)則層權(quán)重后,將各個(gè)指標(biāo)基于ArcGIS“柵格計(jì)算器”進(jìn)行疊加得到川藏鐵路沿線縣域各準(zhǔn)則層,將各準(zhǔn)則層進(jìn)行疊加計(jì)算后得到川藏鐵路沿線縣域地質(zhì)環(huán)境承載力分布,可表示為

    (1)

    式(1)中:Q為地質(zhì)環(huán)境承載力;Ri為影響地質(zhì)環(huán)境質(zhì)量的準(zhǔn)則i;wi為準(zhǔn)則層i的權(quán)重;n為準(zhǔn)則層個(gè)數(shù)。

    4 結(jié)果分析

    4.1 地質(zhì)環(huán)境承載力評價(jià)

    根據(jù)式(1)對各指標(biāo)層數(shù)據(jù)與準(zhǔn)則層數(shù)據(jù)進(jìn)行計(jì)算,分別得到川藏線縣域地貌環(huán)境、地質(zhì)環(huán)境、生態(tài)環(huán)境、社會(huì)經(jīng)濟(jì)環(huán)境以及地質(zhì)環(huán)境承載力空間分布(圖5)。

    圖5 評價(jià)結(jié)果空間分布Fig.5 Spatial distribution of evaluation results

    結(jié)果表明,川藏鐵路沿線縣域地貌環(huán)境總體上表現(xiàn)較好,以中等地貌環(huán)境及較好地貌環(huán)境分布最多范圍最廣,面積最大,占比達(dá)84.46%;地貌環(huán)境較差區(qū)域主要分布于研究區(qū)西部波密縣、加查縣區(qū)域;等級為好主要分布于研究區(qū)東部雅安市城區(qū)、蘆山縣、滎經(jīng)等縣以及研究區(qū)內(nèi)大江大河流經(jīng)所形成河谷地區(qū),區(qū)內(nèi)坡度較小,海拔較低,地形起伏度不大,相對較為平坦,有利于人類活動(dòng)展開,從而使得該類型區(qū)域地貌環(huán)境總體表現(xiàn)較好[圖5(a)]。

    川藏鐵路沿線地質(zhì)環(huán)境與地貌環(huán)境空間分布不同,地質(zhì)環(huán)境較好的區(qū)域主要集中于研究區(qū)西部波密縣、米林縣、加查縣以及拉薩市區(qū)周邊區(qū)域;中等及以上地質(zhì)環(huán)境區(qū)域面積占比為75.76%;地質(zhì)環(huán)境差、較差區(qū)域面積占比24.24%,主要分布于研究區(qū)中東部八宿縣北部、卡若區(qū)、察雅縣以及康定市、天全縣、寶興縣等區(qū)域,該區(qū)域工程巖組巖性硬度整體偏軟,活動(dòng)斷層分布密集,地質(zhì)活動(dòng)頻繁,且由于該區(qū)域主要分布季節(jié)性凍土,伴隨著全球氣候變化的加劇,冰川及凍土消融整體呈現(xiàn)加劇的趨勢,區(qū)域地質(zhì)環(huán)境整體偏差[圖5(b)]。

    川藏鐵路沿線縣域生態(tài)環(huán)境整體偏好,中等及以上生態(tài)環(huán)境占比84.15%,生態(tài)環(huán)境分布趨勢自西向東整體呈現(xiàn)“好-差-好”,中部區(qū)域八宿縣及波密縣北部區(qū)域生態(tài)環(huán)境表現(xiàn)較差。通過分析我們可以發(fā)現(xiàn),波密縣北部區(qū)域坐落于念青唐古拉東段南麓,海拔整體偏高,冰川發(fā)育完全,造成縣域北部區(qū)域植被覆蓋相對較差,生態(tài)環(huán)境相對較差,而八宿縣屬于高原溫帶半干旱季風(fēng)性氣候,日照足,蒸發(fā)量大,同時(shí)受到怒江干熱河谷影響,年降水量較小,植被發(fā)育差。因此,該區(qū)域生態(tài)環(huán)境也相對較差[圖5(c)]。

    社會(huì)環(huán)境作為地質(zhì)環(huán)境承載力評價(jià)的負(fù)向因子,即社會(huì)環(huán)境越好,表明該區(qū)域交通條件較差,人類活動(dòng)強(qiáng)度較小,對地質(zhì)環(huán)境承載力影響較小,從而表示該區(qū)域的開發(fā)潛力較大[33]。川藏線縣域社會(huì)環(huán)境“差”“較差”占比僅為1.45%,且在空間分布上主要分布于雅安市周圍以及拉薩市周邊,該區(qū)域作為川藏鐵路沿線周邊發(fā)展較快的城市,人口分布較多,城市化程度較高,導(dǎo)致其社會(huì)環(huán)境較差,在保證地質(zhì)環(huán)境質(zhì)量的前提下的發(fā)展?jié)摿^??;中等以上社會(huì)環(huán)境占比達(dá)到98.55%,較好、好的社會(huì)環(huán)境面積占比更是達(dá)到81.24%,從理論角度來講,該區(qū)域社會(huì)環(huán)境對于地質(zhì)環(huán)境質(zhì)量的影響是極小的,區(qū)域發(fā)展?jié)摿O大,但由于川藏鐵路沿線自然環(huán)境條件的限制,僅有部分河谷區(qū)域適合發(fā)展[圖5(d)]。

    將地貌環(huán)境、地質(zhì)環(huán)境、生態(tài)環(huán)境以及社會(huì)環(huán)境按照準(zhǔn)則層權(quán)重進(jìn)行疊加分析得到川藏鐵路沿線地質(zhì)環(huán)境承載力分布結(jié)果如圖5(e)所示??傮w上,東部地質(zhì)環(huán)境承載力較差,西部地質(zhì)環(huán)境承載力較好的特征。從面積分布來看,川藏鐵路沿線地質(zhì)環(huán)境承載力以中等最多,占比達(dá)到66.69%,從其空間分布看,中等質(zhì)量在整個(gè)研究區(qū)域都有分布,較集中分布于曲水縣、扎囊縣、貢嘎縣等山南市大部以及巴塘縣、理塘縣和雅江縣等川西高原區(qū)域,其余呈現(xiàn)出零星分布的特征;其次,地質(zhì)環(huán)境承載力較差面積占比達(dá)18.40%。在其空間分布上,較差質(zhì)量主要分布于研究區(qū)中部八宿縣、察雅縣和卡若區(qū)以及寶興縣、堆龍德慶區(qū)、林周縣區(qū)域;較好質(zhì)量面積占比為14.90%,主要分布于米林縣、巴宜區(qū)以及雅江縣中南部[圖5(e)]。

    4.2 影響因素分析

    根據(jù)上文提到的隨機(jī)森林模型的優(yōu)點(diǎn)之一,由模型可以得出指標(biāo)因子的重要性程度,通常將重要性程度作為計(jì)算的權(quán)重,基于這一特點(diǎn),本研究對川藏鐵路沿線縣域地質(zhì)環(huán)境承載力影響因素進(jìn)行分析。將各準(zhǔn)則層與指標(biāo)層重要性程度進(jìn)行統(tǒng)一標(biāo)準(zhǔn)化處理,可以得到川藏鐵路沿線縣域地質(zhì)環(huán)境承載力影響因子重要性程度(圖6)。

    圖6 川藏鐵路沿線地質(zhì)環(huán)境承載力指標(biāo)重要性程度Fig.6 Importance of geo-environmental carrying capacity indicators along the Sichuan-Tibet Railway

    工程巖組作為地質(zhì)環(huán)境質(zhì)量重要的指標(biāo),巖土體類型不同導(dǎo)致巖性硬度及其應(yīng)力分布就存在不同,對地質(zhì)災(zāi)害的發(fā)生與可承受人類活動(dòng)強(qiáng)度產(chǎn)生影響。川藏鐵路沿線工程巖組主要表現(xiàn)為西部地區(qū)巖性硬度較硬,卡若區(qū)、八宿縣至東部康定市這一區(qū)域巖性硬度較軟。由圖6可知,工程巖組對于地質(zhì)環(huán)境承載力的影響為所有指標(biāo)占比最大,為重要性程度最高的指標(biāo)。由于該區(qū)域位于青藏高原東緣,地殼運(yùn)動(dòng)較為活躍,活動(dòng)斷層對于地質(zhì)環(huán)境承載力的影響也相對較大,在地質(zhì)環(huán)境承載力表現(xiàn)較差的八宿縣、察雅縣一帶,班公錯(cuò)-怒江斷裂、日土-丁青主斷裂橫橫亙近東西,從而導(dǎo)致該區(qū)域地質(zhì)環(huán)境質(zhì)量較差。NDVI作為常用表征植被覆蓋指數(shù)的指標(biāo),數(shù)值越大,表明該區(qū)域植被覆蓋度越好,則植被對于地表松散物質(zhì)的固著能力越強(qiáng)。同時(shí),其涵養(yǎng)水分能力也較高,從而對于地質(zhì)環(huán)境承載力的提高具有顯著正向作用;根據(jù)圖6知,NDVI作為影響川藏鐵路沿線地質(zhì)環(huán)境承載力的第二大影響因素,對于沿線地質(zhì)環(huán)境承載力的影響較大。其次,川藏鐵路沿線近年來,隨著區(qū)域社會(huì)經(jīng)濟(jì)的發(fā)展,基礎(chǔ)設(shè)施建設(shè)逐漸完善,在川藏鐵路建設(shè)完成前,川藏公路作為內(nèi)地進(jìn)入西藏的主要途徑,道路的修建往往會(huì)破壞區(qū)域原有的地貌形態(tài),從而使其應(yīng)力結(jié)構(gòu)發(fā)生變化。又道路的修建必然加大車輛的流動(dòng),車輛擾動(dòng)的影響增大,從而導(dǎo)致道路密度成為影響地質(zhì)環(huán)境承載力的重要指標(biāo)。其他影響因素重要性排序依次為海拔、人口密度、坡度、地形起伏度、年均降水量、凍土類型的分布以及居民點(diǎn)密度與災(zāi)地質(zhì)災(zāi)害點(diǎn)密度。

    4.3 評價(jià)結(jié)果驗(yàn)證

    根據(jù)地質(zhì)環(huán)境承載力的指標(biāo)構(gòu)成以及其對于人類社會(huì)發(fā)展的影響,當(dāng)?shù)刭|(zhì)環(huán)境承載力較好時(shí),地貌環(huán)境、地質(zhì)環(huán)境、生態(tài)環(huán)境以及社會(huì)環(huán)境均相對較好。在此背景下,區(qū)域內(nèi)地質(zhì)災(zāi)害的發(fā)生于區(qū)域地質(zhì)環(huán)境承載力的大小呈明顯的負(fù)相關(guān)關(guān)系。因此,本文選取地質(zhì)災(zāi)害的發(fā)生數(shù)量在地質(zhì)環(huán)境承載力等級的占比對地質(zhì)環(huán)境承載力評價(jià)結(jié)果進(jìn)行驗(yàn)證。通過統(tǒng)計(jì)分析,川藏線沿線縣域共發(fā)育地質(zhì)災(zāi)害點(diǎn)3 863處,主要發(fā)育于地質(zhì)環(huán)境承載力“較差”“中等”“較好”3個(gè)等級區(qū)域,其占比分別為17.32%、72.17%以及10.51%。基于此,認(rèn)為該研究地質(zhì)環(huán)境承載力評價(jià)結(jié)果合理,具有一定的科學(xué)性,符合實(shí)際。

    5 結(jié)論

    從川藏鐵路沿線縣域?qū)嶋H出發(fā),從地貌環(huán)境、地質(zhì)環(huán)境、生態(tài)環(huán)境和社會(huì)環(huán)境四個(gè)準(zhǔn)則層選取12個(gè)評價(jià)指標(biāo),同時(shí),考慮到青藏高原獨(dú)特得到凍土分布狀況,選取凍土分布作為地質(zhì)環(huán)境承載力評價(jià)的指標(biāo)因子,構(gòu)建了川藏鐵路沿線縣域地質(zhì)環(huán)境承載力評價(jià)指標(biāo)體系,通過隨機(jī)森林模型對各個(gè)指標(biāo)進(jìn)行權(quán)重賦值,利用GIS空間分析對川藏鐵路沿線縣域地質(zhì)環(huán)境承載力進(jìn)行評價(jià),得出結(jié)論。

    (1)川藏鐵路沿線地質(zhì)環(huán)境承載力受區(qū)域地質(zhì)背景即工程巖組與活動(dòng)斷層的影響較大,其次,以NDVI為代表的生態(tài)環(huán)境對川藏鐵路沿線的地質(zhì)環(huán)境承載力影響也比較大,因此,在改善區(qū)域地質(zhì)環(huán)境質(zhì)量時(shí)應(yīng)該著重生態(tài)修復(fù)。

    (2)川藏鐵路沿線縣域地質(zhì)環(huán)境承載力總體表現(xiàn)良好,中等質(zhì)量以上區(qū)域面積占比超過了60%,且質(zhì)量較差區(qū)域主要集中于八宿縣、察雅縣一帶,面積占比相對較小,除已經(jīng)建好的“成雅段”與“拉林段”鐵路外,雅安至林芝路段鐵路修建中應(yīng)該尤其注意八宿縣、察雅縣一帶。同時(shí),在后期川藏鐵路的維護(hù)過程中還應(yīng)注重加強(qiáng)地質(zhì)災(zāi)害監(jiān)測與預(yù)警,提高區(qū)域防災(zāi)減災(zāi)能力,為區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展保駕護(hù)航。

    猜你喜歡
    鐵路沿線凍土縣域
    縣域消防專項(xiàng)規(guī)劃研究
    北極凍土在求救
    基于分布式光纖傳感技術(shù)的鐵路沿線外部入侵監(jiān)測系統(tǒng)
    杭黃鐵路沿線車站探營
    ——千島湖站
    法人(2019年1期)2019-01-10 01:30:38
    凍土下的猛犸墳場
    山東縣域GDP排名出爐
    商周刊(2017年6期)2017-08-22 03:42:46
    鐵路沿線光電纜溫度智能監(jiān)測系統(tǒng)
    鐵路沿線景觀錯(cuò)亂區(qū)的問題解析及優(yōu)化策略——以哈牡線牡丹江段為例
    住宅科技(2015年11期)2015-12-08 08:20:02
    縣域就診率為何差了40%
    26
    精品久久蜜臀av无| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美在线精品| 亚洲欧美一区二区三区久久| 亚洲国产欧美在线一区| 最新中文字幕久久久久| 久久久久久久久久久久大奶| 久久99蜜桃精品久久| 91aial.com中文字幕在线观看| 亚洲欧美色中文字幕在线| 日产精品乱码卡一卡2卡三| 母亲3免费完整高清在线观看 | 亚洲欧美精品自产自拍| 国产毛片在线视频| 国产毛片在线视频| 伦理电影大哥的女人| 免费观看av网站的网址| 国产一区亚洲一区在线观看| 亚洲精品日本国产第一区| 91成人精品电影| 欧美精品高潮呻吟av久久| 丁香六月天网| 国产极品粉嫩免费观看在线| 久久99精品国语久久久| av免费观看日本| 久久99蜜桃精品久久| 亚洲色图 男人天堂 中文字幕| 国产欧美日韩综合在线一区二区| 色吧在线观看| 国产激情久久老熟女| 国产亚洲最大av| 欧美人与性动交α欧美软件| 亚洲国产看品久久| 久久精品国产亚洲av高清一级| 日日撸夜夜添| av电影中文网址| 国产亚洲精品第一综合不卡| 亚洲精品aⅴ在线观看| 性少妇av在线| 夫妻午夜视频| 久久久精品94久久精品| 18+在线观看网站| 亚洲一区中文字幕在线| 中国国产av一级| 国产精品一区二区在线不卡| 老司机亚洲免费影院| 三级国产精品片| 久久久国产欧美日韩av| 午夜免费男女啪啪视频观看| 99精国产麻豆久久婷婷| 一级,二级,三级黄色视频| 亚洲av成人精品一二三区| 桃花免费在线播放| 免费观看性生交大片5| 久久久久国产精品人妻一区二区| 国产野战对白在线观看| 亚洲第一青青草原| 在线观看免费视频网站a站| 久久久久国产网址| 午夜福利网站1000一区二区三区| 人人妻人人澡人人看| 黄片播放在线免费| 亚洲欧美中文字幕日韩二区| 伊人久久大香线蕉亚洲五| 国产精品不卡视频一区二区| 国产在线视频一区二区| 久久久久久久久久人人人人人人| 永久网站在线| 少妇被粗大猛烈的视频| 中文字幕人妻熟女乱码| 各种免费的搞黄视频| 最近手机中文字幕大全| 国产免费又黄又爽又色| 国产极品天堂在线| 黄色毛片三级朝国网站| 99久久人妻综合| 黑丝袜美女国产一区| 少妇人妻久久综合中文| av国产久精品久网站免费入址| 纯流量卡能插随身wifi吗| 免费大片黄手机在线观看| 日韩中字成人| 丝袜美足系列| 一级爰片在线观看| 大香蕉久久成人网| 91在线精品国自产拍蜜月| 性色av一级| 国产精品蜜桃在线观看| 大片电影免费在线观看免费| 免费av中文字幕在线| www.熟女人妻精品国产| 一本久久精品| 午夜精品国产一区二区电影| 免费看不卡的av| 肉色欧美久久久久久久蜜桃| 亚洲国产av新网站| 青青草视频在线视频观看| 国产免费一区二区三区四区乱码| 最近手机中文字幕大全| 韩国高清视频一区二区三区| 国产一区有黄有色的免费视频| 久久毛片免费看一区二区三区| 90打野战视频偷拍视频| 只有这里有精品99| 国产成人a∨麻豆精品| 精品亚洲乱码少妇综合久久| 男女啪啪激烈高潮av片| 美女xxoo啪啪120秒动态图| 成人二区视频| av.在线天堂| 久久久久国产一级毛片高清牌| 亚洲国产精品成人久久小说| 久久久欧美国产精品| 丝袜人妻中文字幕| av卡一久久| 精品亚洲成国产av| 国产黄色免费在线视频| 亚洲国产欧美网| 久久久国产欧美日韩av| 成人手机av| 又粗又硬又长又爽又黄的视频| 黑人猛操日本美女一级片| 纵有疾风起免费观看全集完整版| 亚洲欧洲国产日韩| 亚洲成人手机| 亚洲美女搞黄在线观看| 人人妻人人爽人人添夜夜欢视频| 天天操日日干夜夜撸| 69精品国产乱码久久久| 久久精品久久久久久噜噜老黄| 欧美精品av麻豆av| 99久久综合免费| 一区二区三区激情视频| 国产淫语在线视频| 69精品国产乱码久久久| 久久久久精品性色| 亚洲中文av在线| 国产97色在线日韩免费| 久久久精品区二区三区| 日本-黄色视频高清免费观看| 亚洲内射少妇av| 免费黄色在线免费观看| 大陆偷拍与自拍| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| 成人免费观看视频高清| av又黄又爽大尺度在线免费看| 最近的中文字幕免费完整| 伦精品一区二区三区| 久久这里只有精品19| 亚洲四区av| 免费看不卡的av| 26uuu在线亚洲综合色| 最近手机中文字幕大全| 高清视频免费观看一区二区| 久久久久国产一级毛片高清牌| 韩国av在线不卡| 天天躁狠狠躁夜夜躁狠狠躁| 深夜精品福利| 欧美精品人与动牲交sv欧美| 久久久国产一区二区| 中文字幕精品免费在线观看视频| 91成人精品电影| 午夜福利网站1000一区二区三区| 97在线视频观看| 少妇 在线观看| 久久久欧美国产精品| av在线观看视频网站免费| 男女边吃奶边做爰视频| 国产黄色视频一区二区在线观看| 蜜桃在线观看..| 成人手机av| xxxhd国产人妻xxx| 久久久久久人妻| 免费高清在线观看日韩| av片东京热男人的天堂| 丝袜美足系列| 老女人水多毛片| 免费高清在线观看日韩| 性少妇av在线| 久久精品国产综合久久久| 亚洲av综合色区一区| 国产精品国产av在线观看| 你懂的网址亚洲精品在线观看| 少妇人妻 视频| 日本爱情动作片www.在线观看| 一边摸一边做爽爽视频免费| 国产又色又爽无遮挡免| 亚洲欧洲国产日韩| 人妻系列 视频| 久久久久久久久久久免费av| 精品午夜福利在线看| 亚洲国产色片| 国产成人精品一,二区| 国产白丝娇喘喷水9色精品| 街头女战士在线观看网站| 欧美另类一区| 亚洲欧美清纯卡通| 欧美成人精品欧美一级黄| 日日摸夜夜添夜夜爱| 少妇人妻精品综合一区二区| 欧美97在线视频| 男人操女人黄网站| a 毛片基地| 伦理电影大哥的女人| √禁漫天堂资源中文www| 亚洲精品国产av蜜桃| 亚洲av日韩在线播放| 综合色丁香网| 赤兔流量卡办理| 日韩制服骚丝袜av| 蜜桃国产av成人99| 男女边摸边吃奶| 国产免费一区二区三区四区乱码| 不卡视频在线观看欧美| 国产男人的电影天堂91| 日日啪夜夜爽| 欧美日韩成人在线一区二区| 欧美日韩视频高清一区二区三区二| 日韩欧美精品免费久久| 亚洲国产欧美在线一区| 免费在线观看视频国产中文字幕亚洲 | 在现免费观看毛片| 天天躁夜夜躁狠狠躁躁| 五月开心婷婷网| 成人国产麻豆网| 色播在线永久视频| 国产日韩欧美亚洲二区| 亚洲成色77777| 久久久久久伊人网av| 日韩欧美精品免费久久| 久久久精品区二区三区| 久久韩国三级中文字幕| 蜜桃国产av成人99| 免费在线观看完整版高清| 国产福利在线免费观看视频| 色婷婷av一区二区三区视频| 青春草亚洲视频在线观看| 建设人人有责人人尽责人人享有的| 另类亚洲欧美激情| freevideosex欧美| 女性被躁到高潮视频| 久久鲁丝午夜福利片| 高清av免费在线| 亚洲图色成人| 如日韩欧美国产精品一区二区三区| 在线 av 中文字幕| 99国产综合亚洲精品| 黄片播放在线免费| 国产精品久久久av美女十八| 亚洲国产av新网站| 狠狠精品人妻久久久久久综合| 国产成人精品久久二区二区91 | 日本爱情动作片www.在线观看| 成年女人在线观看亚洲视频| 少妇的丰满在线观看| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 最近的中文字幕免费完整| 香蕉丝袜av| 亚洲精品久久午夜乱码| 久久亚洲国产成人精品v| 在线观看三级黄色| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久伊人网av| 国产成人欧美| 中文欧美无线码| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 亚洲 欧美一区二区三区| 80岁老熟妇乱子伦牲交| 国产精品免费大片| 午夜影院在线不卡| 一本色道久久久久久精品综合| 国产精品久久久久久久久免| 欧美精品一区二区免费开放| 日韩中文字幕欧美一区二区 | 午夜福利在线免费观看网站| 久久97久久精品| 18+在线观看网站| 久久久久国产一级毛片高清牌| 欧美成人午夜免费资源| 中文字幕人妻熟女乱码| 晚上一个人看的免费电影| 这个男人来自地球电影免费观看 | 精品一区二区三卡| 各种免费的搞黄视频| 久久人人97超碰香蕉20202| 欧美日韩成人在线一区二区| 韩国av在线不卡| 永久免费av网站大全| av国产久精品久网站免费入址| 老汉色av国产亚洲站长工具| 国产成人av激情在线播放| 午夜免费观看性视频| 成人国语在线视频| 精品亚洲成国产av| 七月丁香在线播放| 亚洲少妇的诱惑av| 日日撸夜夜添| 极品少妇高潮喷水抽搐| av片东京热男人的天堂| av在线app专区| 久久精品国产自在天天线| 黄网站色视频无遮挡免费观看| 母亲3免费完整高清在线观看 | 亚洲成色77777| 我的亚洲天堂| 人人妻人人爽人人添夜夜欢视频| 9色porny在线观看| 免费少妇av软件| 哪个播放器可以免费观看大片| 亚洲色图综合在线观看| 国产成人免费无遮挡视频| 一级爰片在线观看| 大香蕉久久成人网| 欧美人与性动交α欧美软件| 亚洲国产日韩一区二区| 国产亚洲午夜精品一区二区久久| 99久久人妻综合| 欧美亚洲 丝袜 人妻 在线| 王馨瑶露胸无遮挡在线观看| 伦精品一区二区三区| 大码成人一级视频| 一区二区三区激情视频| 久久久久久久久久久免费av| 久久久国产一区二区| 亚洲成人一二三区av| 午夜福利在线观看免费完整高清在| 啦啦啦在线免费观看视频4| 日韩一区二区视频免费看| av免费观看日本| 性色avwww在线观看| 久久这里有精品视频免费| 只有这里有精品99| 王馨瑶露胸无遮挡在线观看| 成年动漫av网址| 久久久久精品久久久久真实原创| 美女国产高潮福利片在线看| 日韩,欧美,国产一区二区三区| 亚洲成国产人片在线观看| 国产老妇伦熟女老妇高清| 99re6热这里在线精品视频| 中文字幕色久视频| 在线免费观看不下载黄p国产| 夫妻性生交免费视频一级片| 亚洲情色 制服丝袜| 欧美精品高潮呻吟av久久| av福利片在线| 丰满饥渴人妻一区二区三| 午夜福利,免费看| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 日本欧美国产在线视频| 国产在线免费精品| 两个人看的免费小视频| 国产日韩欧美在线精品| 久久精品国产自在天天线| 狂野欧美激情性bbbbbb| 久久久国产精品麻豆| 伦理电影免费视频| 国产亚洲av片在线观看秒播厂| av福利片在线| 少妇 在线观看| 免费大片黄手机在线观看| 国产精品亚洲av一区麻豆 | 亚洲精品av麻豆狂野| 日日啪夜夜爽| 日产精品乱码卡一卡2卡三| xxxhd国产人妻xxx| videos熟女内射| 视频在线观看一区二区三区| 亚洲婷婷狠狠爱综合网| 18+在线观看网站| 日韩三级伦理在线观看| 日韩电影二区| 国产片内射在线| 久久久欧美国产精品| 成人影院久久| 中文字幕人妻丝袜一区二区 | 欧美成人精品欧美一级黄| 国产熟女午夜一区二区三区| 韩国精品一区二区三区| 七月丁香在线播放| 日韩视频在线欧美| 丰满乱子伦码专区| 人妻一区二区av| 日韩成人av中文字幕在线观看| 久久精品国产亚洲av涩爱| 啦啦啦中文免费视频观看日本| www.av在线官网国产| 在线观看免费高清a一片| 国产免费现黄频在线看| 亚洲综合精品二区| 交换朋友夫妻互换小说| 侵犯人妻中文字幕一二三四区| 亚洲精品自拍成人| 免费看av在线观看网站| 人人妻人人澡人人看| 伊人久久国产一区二区| av又黄又爽大尺度在线免费看| 精品一区二区三卡| 国产亚洲最大av| 亚洲第一青青草原| 免费观看无遮挡的男女| 国产av码专区亚洲av| av不卡在线播放| 大话2 男鬼变身卡| 精品少妇一区二区三区视频日本电影 | 精品人妻熟女毛片av久久网站| 我的亚洲天堂| 亚洲成色77777| 美国免费a级毛片| 男女免费视频国产| 亚洲五月色婷婷综合| 两性夫妻黄色片| 亚洲国产成人一精品久久久| 久久久久久伊人网av| 亚洲国产毛片av蜜桃av| 美女主播在线视频| 99热全是精品| 亚洲欧美色中文字幕在线| 欧美日本中文国产一区发布| 亚洲欧洲日产国产| 久久精品国产鲁丝片午夜精品| 久久精品国产自在天天线| 黄色 视频免费看| 亚洲人成电影观看| 亚洲av免费高清在线观看| 成年人免费黄色播放视频| 欧美少妇被猛烈插入视频| 亚洲欧美日韩另类电影网站| 男女下面插进去视频免费观看| 欧美人与性动交α欧美软件| 欧美xxⅹ黑人| 美女国产高潮福利片在线看| 九草在线视频观看| av国产精品久久久久影院| 十八禁网站网址无遮挡| 日本vs欧美在线观看视频| 久久精品久久久久久久性| 狠狠婷婷综合久久久久久88av| 免费观看性生交大片5| 久久这里有精品视频免费| 十分钟在线观看高清视频www| 国产爽快片一区二区三区| 国产精品国产三级专区第一集| 丝袜脚勾引网站| 黄频高清免费视频| 丰满少妇做爰视频| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线免费观看视频4| 亚洲精品美女久久av网站| 久久综合国产亚洲精品| 国产又爽黄色视频| 国产精品免费大片| 国产成人一区二区在线| 国产免费一区二区三区四区乱码| 边亲边吃奶的免费视频| 亚洲精品中文字幕在线视频| 亚洲av男天堂| 欧美变态另类bdsm刘玥| 日韩成人av中文字幕在线观看| 午夜福利在线免费观看网站| 下体分泌物呈黄色| 波多野结衣一区麻豆| 女人高潮潮喷娇喘18禁视频| 黑人猛操日本美女一级片| 亚洲精品美女久久久久99蜜臀 | av免费在线看不卡| 777米奇影视久久| 亚洲综合色惰| 精品少妇内射三级| 18禁观看日本| 久久毛片免费看一区二区三区| 又粗又硬又长又爽又黄的视频| 男人添女人高潮全过程视频| 97在线视频观看| 亚洲精品一区蜜桃| 日韩一卡2卡3卡4卡2021年| 伊人久久国产一区二区| 最新中文字幕久久久久| av线在线观看网站| av免费在线看不卡| 久久精品国产亚洲av天美| 国产精品秋霞免费鲁丝片| 国产深夜福利视频在线观看| 国产一区二区在线观看av| 99久久人妻综合| 亚洲国产av新网站| 亚洲图色成人| 精品久久久久久电影网| 精品亚洲乱码少妇综合久久| 亚洲国产av影院在线观看| 亚洲精品国产av蜜桃| 欧美bdsm另类| 狠狠精品人妻久久久久久综合| 久久精品国产自在天天线| 18禁观看日本| videossex国产| 国产免费现黄频在线看| 国产精品久久久久久av不卡| 国产欧美亚洲国产| 久久青草综合色| 久久鲁丝午夜福利片| 国精品久久久久久国模美| 国产成人精品福利久久| 亚洲,一卡二卡三卡| 五月开心婷婷网| 最近中文字幕高清免费大全6| 精品少妇黑人巨大在线播放| 在线观看美女被高潮喷水网站| 在线精品无人区一区二区三| 综合色丁香网| 亚洲国产日韩一区二区| 国产免费视频播放在线视频| 精品酒店卫生间| 亚洲国产欧美网| 熟女电影av网| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 中文乱码字字幕精品一区二区三区| 满18在线观看网站| 欧美日韩av久久| 亚洲激情五月婷婷啪啪| 秋霞在线观看毛片| 国产精品成人在线| 999精品在线视频| 最新中文字幕久久久久| 高清在线视频一区二区三区| 国产97色在线日韩免费| 亚洲欧美色中文字幕在线| 高清av免费在线| av视频免费观看在线观看| 制服人妻中文乱码| 91精品三级在线观看| 亚洲,一卡二卡三卡| 国产成人精品在线电影| av电影中文网址| 亚洲一码二码三码区别大吗| 国产成人aa在线观看| 免费高清在线观看视频在线观看| 热re99久久精品国产66热6| av福利片在线| 十八禁高潮呻吟视频| 精品久久蜜臀av无| 国产淫语在线视频| 在线观看国产h片| 中文字幕人妻丝袜一区二区 | 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| av国产久精品久网站免费入址| 大话2 男鬼变身卡| 一级黄片播放器| 美女主播在线视频| 丝瓜视频免费看黄片| 精品一区二区免费观看| 久久久国产欧美日韩av| 国产亚洲欧美精品永久| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 欧美激情极品国产一区二区三区| 久久午夜福利片| 久久精品国产亚洲av高清一级| 水蜜桃什么品种好| 在线观看免费日韩欧美大片| 亚洲国产成人一精品久久久| 亚洲av男天堂| 成年动漫av网址| 亚洲情色 制服丝袜| 亚洲中文av在线| 亚洲欧美清纯卡通| 男的添女的下面高潮视频| 建设人人有责人人尽责人人享有的| 欧美日韩一级在线毛片| 韩国精品一区二区三区| 精品午夜福利在线看| 日韩在线高清观看一区二区三区| 黄色毛片三级朝国网站| 99精国产麻豆久久婷婷| 中文乱码字字幕精品一区二区三区| 丝袜人妻中文字幕| 飞空精品影院首页| 少妇人妻久久综合中文| 国产色婷婷99| 丰满乱子伦码专区| 99国产精品免费福利视频| 精品卡一卡二卡四卡免费| 国产成人精品久久久久久| 在线免费观看不下载黄p国产| 国产 一区精品| 黑人巨大精品欧美一区二区蜜桃| 国产男人的电影天堂91| 久久亚洲国产成人精品v| 卡戴珊不雅视频在线播放| 永久网站在线| 久久精品国产鲁丝片午夜精品| 久久久久视频综合| 少妇的逼水好多| 在线观看免费高清a一片| 欧美日韩精品网址| 在线观看免费日韩欧美大片| 亚洲av在线观看美女高潮| 黑人猛操日本美女一级片| 国产成人精品福利久久| 9热在线视频观看99| 99香蕉大伊视频| 午夜影院在线不卡| 久久人人爽av亚洲精品天堂| 成人亚洲欧美一区二区av| 一边摸一边做爽爽视频免费| 性少妇av在线| 午夜免费观看性视频| 啦啦啦在线免费观看视频4|