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

    基于證據(jù)權(quán)和卡方自動交互檢測決策樹的滑坡易發(fā)性預(yù)測

    2022-06-22 08:43:39黃發(fā)明石雨歐陽慰平洪安宇曾子強徐富剛
    土木建筑與環(huán)境工程 2022年5期

    黃發(fā)明 石雨 歐陽慰平 洪安宇 曾子強 徐富剛

    摘 要:滑坡與其環(huán)境因子間的非線性關(guān)聯(lián)計算影響滑坡易發(fā)性預(yù)測建模的不確定性。為研究不確定性因素下易發(fā)性建模規(guī)律,以中國延長縣為例,獲取82處滑坡和14種環(huán)境因子,通過頻率比(Frequency Ratio,F(xiàn)R)和證據(jù)權(quán)(Weight of Evidence,WOE)等關(guān)聯(lián)法與卡方自動交互檢測(Chi-squared Automatic Interaction Detector,CHAID)決策樹相耦合進行建模,并用原始環(huán)境因子(稱為“原始因子數(shù)據(jù)”)作為輸入變量的單獨CHAID決策樹進行對比。使用精度、易發(fā)性指數(shù)均值、標(biāo)準(zhǔn)差和平均秩等評價易發(fā)性建模的不確定性。結(jié)果表明:WOE-CHAID模型預(yù)測的滑坡易發(fā)性不確定性低于FR-CHAID模型,可見WOE具有較優(yōu)秀的非線性關(guān)聯(lián)性能;單獨CHAID決策樹預(yù)測的易發(fā)性精度整體略低于WOE-CHAID和FR-CHAID模型,但其建模效率較高;在體現(xiàn)滑坡與其環(huán)境因子空間關(guān)聯(lián)性方面,考慮FR和WOE關(guān)聯(lián)法的CHAID決策樹模型優(yōu)勢顯著。WOE是更優(yōu)秀的關(guān)聯(lián)分析法,CHAID決策樹預(yù)測性能好且預(yù)測效率高,WOE-CHAID決策樹模型的易發(fā)性預(yù)測不確定性較低且更符合實際滑坡概率分布特征。

    關(guān)鍵詞:滑坡易發(fā)性預(yù)測;關(guān)聯(lián)分析;CHAID決策樹;證據(jù)權(quán);環(huán)境因子

    中圖分類號:P642.22 文獻標(biāo)志碼:A 文章編號:2096-6717(2022)05-0016-13

    收稿日期:2021-09-02

    基金項目:國家自然科學(xué)基金(52109089、41807285)

    作者簡介:黃發(fā)明(1988- ),男,博士,副教授,主要從事地質(zhì)災(zāi)害風(fēng)險預(yù)警研究,E-mail:faminghuang@ncu.edu.cn。

    洪安宇(通信作者),女,博士,E-mail:honganyu@ncu.edu.cn。

    Received:2021-09-02

    Foundation items:National Natural Science Foundation of China (No. 52109089, 41807285)

    Author brief:Huang Faming (1988- ), PhD, associate professor, main research interest: geological disaster risk warning, E-mail: faminghuang@ncu.edu.cn.

    HONG Anyu (corresponding author), PhD, E-mail: honganyu@ncu.edu.cn.

    Landslide susceptibility prediction modeling based on weight of evidence and chi-square automatic interactive detection decision tree

    HUANG Faming, SHI Yu, OUYANG Weiping, HONG Anyu, ZENG Ziqiang, XU Fugang

    (School of Civil Engineering and Architecture, Nanchang University, Nanchang 330031, P. R. China)

    Abstract:

    The calculation of the non-linear correlation between the landslide inventories and their environmental factors is an important factor that affects the uncertainty of the landslide susceptibility prediction (LSP) modeling. In order to study the changing patterns of LSP under the influence of the uncertain factors, taking Yanchang County of China as example, 82 landslides and 14 environmental factors are obtained, and the frequency ratio (FR) and weight of evidence (WOE) connection methods are coupled with the chi-squared automatic interaction detector (CHAID) decision tree model to carry out LSP. Then the original environmental factors data (hereinafter referred to as "original data") is used as the input variable to compare the individual CHAID decision tree model to realize the analysis of LSP modeling pattern. ROC accuracy, mean, standard deviation, and average rank are adopted to analyze the uncertainty characteristics in the LSP modeling process. Results show that: 1) LSP uncertainty of the WOE-CHAID model is lower than that of the FR-CHAID model, and WOE has relatively excellent nonlinear correlation performance. 2) The prediction accuracy of individual CHAID decision tree model is slightly lower than that of the WOE-CHAID and FR-CHAID models, but it has higher modeling efficiency. 3) In terms of reflecting the spatial correlation between landslides and its environmental factors, the CHAID decision tree model coupled with FR and WOE connection methods have significant advantages. Generally, WOE is a better connection method and CHAID decision tree model has good prediction performance and high prediction efficiency. Susceptibility prediction by the WOE-CHAID decision tree model is less uncertain and more in line with the actual landslide probability distribution characteristics.

    Keywords:landslide susceptibility prediction; connection method; CHAID decision tree; weight of evidence; environmental factor

    如何有效開展滑坡易發(fā)性預(yù)測制圖是現(xiàn)階段全世界范圍內(nèi)區(qū)域滑坡研究的重點和難點。通過將GIS與數(shù)據(jù)驅(qū)動模型相結(jié)合,以圖像和數(shù)字的方式可構(gòu)建出更高效準(zhǔn)確的易發(fā)性預(yù)測模型。該易發(fā)性制圖的思路對滑坡高發(fā)地區(qū)的防災(zāi)減災(zāi)規(guī)劃具有重要意義。

    滑坡易發(fā)性可定義為特定地點在環(huán)境因子非線性耦合作用下發(fā)生滑坡的空間概率。基于地理相似性規(guī)律,即“地理環(huán)境越相似,地理特征越相近”可知,通過已經(jīng)發(fā)生滑坡的環(huán)境因子來建立預(yù)測模型,則潛在滑坡的空間位置有可能被預(yù)測。很明顯,從滑坡樣本點中確定滑坡易發(fā)性與其環(huán)境因子的關(guān)系式是易發(fā)性預(yù)測的關(guān)鍵所在,因此,選擇用以獲取輸入變量的滑坡環(huán)境因子關(guān)聯(lián)分析法非常重要。隨著遙感和GIS等基礎(chǔ)數(shù)據(jù)源獲取技術(shù)的進步,易發(fā)性建模的空間數(shù)據(jù)源及其質(zhì)量有了較大提升。一般而言,具體研究區(qū)內(nèi)的滑坡環(huán)境因子類型可通過相關(guān)文獻綜述和研究區(qū)的自然地理和地質(zhì)條件確定。筆者重點關(guān)注滑坡易發(fā)性建模過程中滑坡與其環(huán)境因子的非線性關(guān)聯(lián)分析這一不確定性因素,并進一步研究其對滑坡易發(fā)性建模的影響。

    啟發(fā)式模型、數(shù)理統(tǒng)計模型和機器學(xué)習(xí)模型是易發(fā)性預(yù)測過程中常用的3種類型。啟發(fā)式模型和數(shù)理統(tǒng)計模型被大量使用,主要有確定性因子(Certainty Factors,CF)、層次分析法和多元線性回歸等;機器學(xué)習(xí)相關(guān)模型包括邏輯回歸(Logistic Regression,LR)、C5.0決策樹、人工神經(jīng)網(wǎng)絡(luò)、隨機森林(Random Forest,RF)、支持向量機(Support Vector Machines,SVM)、卡方自動交互檢測(Chi-squared Automatic Interaction Detector,CHAID)決策樹和貝葉斯網(wǎng)絡(luò)等。對于哪種類型的模型最適合易發(fā)性預(yù)測,現(xiàn)階段還沒有一致的意見,但優(yōu)秀的機器學(xué)習(xí)模型能夠提高滑坡易發(fā)性預(yù)測精度,對滑坡易發(fā)性區(qū)間劃分有著顯著影響,并可能進一步改變滑坡易發(fā)性級別的劃分。筆者擬用CHAID決策樹這一被廣泛應(yīng)用的典型機器學(xué)習(xí)方法構(gòu)建滑坡易發(fā)性模型并探索建模不確定性特征。

    在將建模預(yù)測出的滑坡易發(fā)性指數(shù)(Landslide Susceptibility Index,LSIs)與各類環(huán)境因子開展聯(lián)系時,需開展滑坡與其基礎(chǔ)環(huán)境因子(不考慮誘發(fā)因子)之間的非線性關(guān)聯(lián)分析,其關(guān)聯(lián)值可直接作為易發(fā)性模型的輸入變量。目前,常用的關(guān)聯(lián)分析法包括確定系數(shù)、頻率比(Frequency Ratio,F(xiàn)R)、熵指數(shù)(Index of Entropy,IOE)和證據(jù)權(quán)重(Weight of Evidence,WOE)等。不同關(guān)聯(lián)分析法的內(nèi)部計算思路具有較大的差異性,導(dǎo)致各方法下的易發(fā)性建模存在不確定性。關(guān)聯(lián)分析法太粗糙會導(dǎo)致部分信息丟失,降低模型預(yù)測精度;優(yōu)秀的關(guān)聯(lián)分析法能獲取較準(zhǔn)確的環(huán)境因子影響滑坡發(fā)育的信息,進一步提高滑坡環(huán)境因子分析及其建模的可靠性??梢姡接懖煌P(guān)聯(lián)分析法對易發(fā)性預(yù)測建模的影響規(guī)律具有重要意義。

    學(xué)者們采用不同關(guān)聯(lián)分析法和模型開展易發(fā)性預(yù)測建模,例如:Zhang等應(yīng)用IOE模型、LR-IOE和SVM-IOE模型獲得了中國陜西省府谷縣滑坡易發(fā)性圖,結(jié)果表明,LR-IOE模型的準(zhǔn)確率最高,其次是IOE模型和SVM-IOE模型。李文彬等深入探討滑坡與其環(huán)境因子間的非線性聯(lián)接以及不同數(shù)據(jù)驅(qū)動模型對滑坡易發(fā)性預(yù)測建模不確定性的影響規(guī)律,結(jié)果表明,RF模型預(yù)測性能最優(yōu),WOE-RF模型預(yù)測的滑坡易發(fā)性不確定性較低。張鐘遠等基于地理信息系統(tǒng)平臺構(gòu)建了云南省鎮(zhèn)康縣滑坡易發(fā)性預(yù)測指標(biāo)體系,結(jié)果顯示,頻率比耦合LR模型具有更高的成功率和預(yù)測率。但大多數(shù)情況下,現(xiàn)有研究使用特定的關(guān)聯(lián)分析法開展易發(fā)性預(yù)測建模,而較少提供可信的依據(jù)和合理的解釋,并且較少深入探討這種不確定性因素對易發(fā)性預(yù)測建模的影響。通過探討關(guān)聯(lián)分析法耦合模型下的滑坡易發(fā)性結(jié)果的不確定性,更能深入理解易發(fā)性預(yù)測的可靠性和可行性,可降低關(guān)聯(lián)分析法不確定性因素帶來的影響。

    筆者采用FR和WOE兩種非線性關(guān)聯(lián)分析法的計算數(shù)據(jù)值與原始環(huán)境因子數(shù)據(jù)(以下簡稱“原始因子數(shù)據(jù)”)作為CHAID決策樹模型的輸入變量,以陜西省延長縣為例,開展滑坡易發(fā)性預(yù)測建模的不確定性分析,包括精度評價、LSIs分布規(guī)律和平均秩等。

    1 滑坡易發(fā)性建模分析

    FR和WOE兩種關(guān)聯(lián)法耦合CHAID決策樹模型時的易發(fā)性預(yù)測建模流程(圖1)如下:

    1)獲取研究區(qū)滑坡編錄及相關(guān)環(huán)境因子數(shù)據(jù)源以便構(gòu)建易發(fā)性建模的空間數(shù)據(jù)集;

    2)將FR、WOE和原始因子數(shù)據(jù)作為CHAID決策樹的輸入變量,形成3種耦合模型;

    3)分別對3種耦合模型開展易發(fā)性預(yù)測建模,然后在GIS中繪制滑坡易發(fā)性圖并劃分易發(fā)性等級;

    4)通過ROC精度、均值、標(biāo)準(zhǔn)差和平均秩等對易發(fā)性預(yù)測結(jié)果進行不確定分析;

    5)通過對比分析找到最佳關(guān)聯(lián)分析法,為易發(fā)性建模提供指導(dǎo)。

    1.1 滑坡與環(huán)境因子的關(guān)聯(lián)分析法

    1.1.1 頻率比

    頻率比(Frequency Ratio,F(xiàn)R)反映了滑坡在各環(huán)境因子類別的分布狀況,闡述環(huán)境因子各屬性區(qū)間對滑坡的相對影響度,并且能夠很好地解釋滑坡與各因子之間的內(nèi)在聯(lián)系。FR>1代表在對應(yīng)的環(huán)境因子條件下利于滑坡事件的發(fā)生;FR<1表明該環(huán)境因子區(qū)間的屬性與滑坡的發(fā)展關(guān)系較弱。利用環(huán)境因子的FR值作為各模型的輸入變量之一,其計算公式如式(1)。

    FR=N/NS/S(1)

    式中:N為環(huán)境因子某區(qū)間中出現(xiàn)的滑坡柵格數(shù);N是全區(qū)已知滑坡所分布柵格的總數(shù);S是環(huán)境因子的單元數(shù);S是全區(qū)柵格總數(shù)。

    1.1.2 證據(jù)權(quán)

    證據(jù)權(quán)(Weight of Evidence,WOE)法在貝葉斯準(zhǔn)則基礎(chǔ)上綜合各類證據(jù)層來實現(xiàn)定量計算某事件的發(fā)生概率。WOE法通過將滑坡編錄和各類環(huán)境因子層進行空間關(guān)聯(lián),從而得到滑坡處各環(huán)境因子的詳細分布特征權(quán)重因子W+和W-,其在每個環(huán)境因子分級中的計算如式(2)、式(3)所示。

    W+=lnB/(B+B)D/(D+D)(2)

    W-=lnB/(B+B)D/(D+D)(3)

    式中:W+和W-分別為環(huán)境因子存在區(qū)和不存在區(qū)的權(quán)重值,對于原始因子數(shù)據(jù)缺失的區(qū)域其權(quán)重值為0;B和D分別為環(huán)境因子存在區(qū)的滑坡和非滑坡單元數(shù),B和D分別是環(huán)境因子不存在區(qū)的滑坡和非滑坡單元數(shù)。證據(jù)層和滑坡點正相關(guān)時,W+>0和W-<0,負相關(guān)時,W+<0和W->0,在不相關(guān)或數(shù)據(jù)缺失時,權(quán)重為0。進一步利用相對系數(shù)C=W+-W-實現(xiàn)證據(jù)層和滑坡間相關(guān)性的衡量。

    1.2 卡方自動交互檢測決策樹

    CHAID決策樹以卡方統(tǒng)計量為基礎(chǔ)實現(xiàn)最優(yōu)決策樹構(gòu)建,也就是通過自變量和因變量間的解釋性來實現(xiàn)因變量的自動判別。CHAID決策樹具有強大的非線性擬合預(yù)測性能,能容忍樣本數(shù)據(jù)缺失及樣本量不足等缺陷。CHAID模型設(shè)定樹生長的層數(shù)、分裂及聚合閾值等停止標(biāo)準(zhǔn)來構(gòu)建準(zhǔn)確高效的預(yù)測或分類模型,同時,為防止過擬合現(xiàn)象而用隨機分成的訓(xùn)練樣本構(gòu)建模型;最后再利用隨機分成的測試樣本對CHAID進行逐步檢驗,以修正模型參數(shù)。

    1.3 不確定性分析方法

    1.3.1 ROC 曲線精度分析

    采用受試者工作特征(Receiver Operating Characteristic,ROC)曲線下面積(Area Under ROC,AUC)值作為一種量化指標(biāo)來整體評估建模性能。ROC曲線對測試集中各樣例進行排序并依序選擇各截斷點,再逐個把樣例作為正例來進行計算,依據(jù)當(dāng)前分類器的“真陽率”和“假陽率”進行ROC曲線的繪制,相關(guān)評價指標(biāo)如表1所示。AUC值等于隨機挑選的正樣本的排名高于隨機挑選的負樣本的概率,AUC值越大,則易發(fā)性模型預(yù)測性能越好。

    1.3.2 易發(fā)性指數(shù)統(tǒng)計規(guī)律分析

    均值(Mean)是集中趨勢的測量,計算如式(4)所示(式中:X為第n個柵格單元的滑坡易發(fā)性指數(shù)值),其量化了研究區(qū)LSIs分布的整體偏向趨勢,反映了LSIs分布的平均水平。標(biāo)準(zhǔn)差(Standard Deviation)是對圍繞平均值的離差的測量,計算如式(5)所示(式中:μ為滑坡易發(fā)性指數(shù)均值;X為第i個柵格單元的滑坡易發(fā)性指數(shù)值),量化了LSIs分布的離散程度,標(biāo)準(zhǔn)差越小,說明LSIs越接近平均值,反之,則說明其與平均值的差異越大。采用均值和標(biāo)準(zhǔn)差從整體上分析LSIs的分布特征,揭示不同關(guān)聯(lián)分析法和模型耦合模型下的預(yù)測性能,為滑坡易發(fā)性研究提供理論指導(dǎo)。

    Mean=X+X+…+Xn(4)

    Standard Deviation=1N∑Ni=1X-μ(5)

    1.3.3 易發(fā)性指數(shù)的差異顯著性

    采用顯著性差異水平進一步分析各耦合模型下易發(fā)性建模的不確定性。具體采用Kendall協(xié)同系數(shù)檢驗法,對任意兩組不同耦合模型下預(yù)測出的LSIs進行差異顯著性檢驗。若Kendall秩相關(guān)系數(shù)W小于1及檢驗結(jié)果的顯著性小于0.05,說明這兩組耦合模型下LSIs的差異是顯著的,拒絕原假設(shè)。本文通過成對因子顯著性檢驗發(fā)現(xiàn),W值為0.139,小于1,且P值均小于0.05,可見,各耦合模型下的LSIs間差異顯著。

    2 延長縣簡介及環(huán)境因子分析

    2.1 延長縣簡介及滑坡編錄

    延長縣位于陜西東部,面積約2 368.7 km,地勢從西北向東南方向傾斜??h境內(nèi)屬黃土高原丘陵溝壑區(qū)(河谷階地、黃土溝谷區(qū)、黃土溝間區(qū)和巖質(zhì)丘陵區(qū)),出露三疊系中上統(tǒng)內(nèi)陸湖相碎屑沉積巖和第四系風(fēng)積、沖洪積和堆積黃土等地層,新近系砂礫巖在研究區(qū)出露較少(圖2)。另外,縣境內(nèi)地質(zhì)構(gòu)造活動強度低,屬于暖溫帶干旱大陸性季風(fēng)氣候,年均降雨量約564 mm且集中在7、8、9月份。

    根據(jù)已有的滑坡野外調(diào)查資料和數(shù)據(jù)庫可知,延長縣共發(fā)生滑坡82處,主要類型為小型淺層覆蓋滑坡,主要運動方式為牽引式(59%)和推移式滑動(41%);縣境內(nèi)的小型滑坡45處(占比54.8%),中型滑坡36處(占比43.9%),大型滑坡只有1處。延長縣滑坡分布位置如圖2所示,滑坡主要分布在縣域西部及周邊地區(qū),東部和中部較少;大部分發(fā)生滑坡的位置地勢較高,距離河流水系也較近。延長縣滑坡的發(fā)生與地層巖性和工程活動密切相關(guān)。

    2.2 環(huán)境因子分析

    2.2.1 環(huán)境因子介紹

    根據(jù)延長境內(nèi)滑坡的特征及相關(guān)參考文獻的介紹,利用遙感影像和GIS軟件系統(tǒng)從數(shù)據(jù)源中提取14類滑坡環(huán)境因子,包括地形、水文、地表覆被和基礎(chǔ)地質(zhì)等。其中,高程、NDVI、NDBI和MNDWI等12個因子為連續(xù)型數(shù)據(jù),而距河流距離和地層巖性2個因子為離散型數(shù)據(jù)(表2)。對于連續(xù)型環(huán)境因子,先通過小間隔對該因子進行等分,再依據(jù)FR和WOE值將數(shù)值相近的區(qū)間合并成一個類別。對于離散型數(shù)據(jù)類型的環(huán)境因子,采用固有的自然分組來進行分級:距河流距離因子按照距河流距離100、300、400、500、800、900、1 000 m和大于1 000 m進行分類;地層巖性因子為三疊系砂巖夾砂質(zhì)泥巖和油頁巖(Tt)、三疊系厚層砂巖夾泥巖(Th)、三疊系細砂層粉砂巖夾與泥巖互層(Ty)、三疊系厚層狀長石石英砂巖(Tw)和第四系更新統(tǒng)風(fēng)積和洪積黃土(Q)。另外,在使用原始因子數(shù)據(jù)作為CHAID決策樹模型的輸入變量時,將距河流的距離和地層巖性兩種離散型數(shù)據(jù)類型的環(huán)境因子進行了“啞變量”處理。

    2.2.2 地形地貌因子

    高程、坡度、坡向、剖面曲率、平面曲率、地形起伏度、地形粗糙度、地形切割深度和地形濕度指數(shù)等環(huán)境因子均從DEM中提?。▓D3) 。以地形起伏度為例,分析其8個等級區(qū)間內(nèi)的FR和WOE值(表2),發(fā)現(xiàn)滑坡發(fā)生概率與研究區(qū)的地形起伏度大小成正比。在20~4區(qū)間內(nèi)發(fā)生滑坡的概率最大,為78.34%;其中,F(xiàn)R值均大于1,WOE值均為正值,35~40區(qū)域內(nèi)FR和WOE值最大,分別為2.843和1.148。FR和WOE值都顯示出地形起伏度大小與滑坡發(fā)生有著較強的正向相關(guān)性,可見關(guān)聯(lián)分析法在表達滑坡與地形起伏度的非線性關(guān)聯(lián)性時具有較為一致的趨勢和計算效果。

    2.2.3 水文環(huán)境因子

    由于河流對邊坡的浸潤和侵蝕作用,越靠近河流的邊坡土壤含水量可能越高,導(dǎo)致斜坡體失穩(wěn)的可能性更高。利用距河流距離和MNDWI來表征水文環(huán)境對滑坡發(fā)育的影響。以距河流的距離因子為例(表2),當(dāng)距河流距離小于400 m時,滑坡發(fā)育的可能性更高(達74.41%),其中,F(xiàn)R值均大于1,WOE值均為正值;在100~300 m區(qū)域內(nèi),F(xiàn)R和WOE值最大,分別為1.873和0.992。

    2.2.4 地表覆被因子

    NDBI和NDVI分別反映了研究區(qū)域內(nèi)的建筑分布和自然植被對滑坡地質(zhì)災(zāi)害發(fā)育的影響。從表2可知,當(dāng)NDVI在0.121~0.424范圍內(nèi)時,其與滑坡有較強的關(guān)系,該區(qū)間包括了研究區(qū)內(nèi)近年來所有的已發(fā)生的滑坡;其中,在0.121~0.182范圍內(nèi),F(xiàn)R值大于1且WOE值為正數(shù)。NDBI能較好地反映研究區(qū)域內(nèi)建筑的分布情況,當(dāng)NBVI在0.730~0.949范圍內(nèi)時幾乎囊括了近年來研究區(qū)內(nèi)所有的滑坡,間接反映了人類工程建設(shè)對滑坡發(fā)育的影響。

    2.2.5 基礎(chǔ)地質(zhì)因子

    巖土類型表征滑坡體的物質(zhì)基礎(chǔ),分析表2可知,Th和Ty巖性區(qū)域面積僅占延長縣面積的10.6%,而區(qū)域內(nèi)滑坡發(fā)生的概率高達23.2%,且FR值均大于1、WOE值均為正值,說明Th和Ty巖性區(qū)域內(nèi)滑坡發(fā)生的頻率較高;在Q巖性條件下,滑坡發(fā)生概率高達76.8%;在Tt巖性區(qū)域內(nèi),無滑坡分布;Tw巖性區(qū)域在研究區(qū)內(nèi)占比比較小,結(jié)果不具有研究意義。

    3 延長縣滑坡易發(fā)性預(yù)測建模

    3.1 數(shù)據(jù)準(zhǔn)備

    30 m分辨率的柵格被廣泛用作滑坡易發(fā)性的制圖單元,基于30 m分辨率,整個延長縣被劃分為2 622 482個柵格,已發(fā)生的82處滑坡被劃分為3 403個滑坡柵格。通過FR和WOE兩種關(guān)聯(lián)法對14個環(huán)境因子各屬性區(qū)間進行重新賦值,作為CHAID決策樹開展易發(fā)性建模的輸入變量;同時,也以原始因子數(shù)據(jù)作為輸入變量開展單獨CHAID決策樹的滑坡易發(fā)性建模。通過SPSS modeler 18.0軟件把3 403個滑坡柵格單元賦值為1,同時隨機挑選與滑坡單元相同數(shù)量的非滑坡單元,并將其易發(fā)性賦值為0,作為模型輸出變量;然后按7∶3隨機劃分滑坡和非滑坡柵格單元(6 806個)及其相關(guān)屬性值,得到模型訓(xùn)練集和測試集。最后將整個研究區(qū)柵格單元的FR和WOE關(guān)聯(lián)分析值以及原始因子數(shù)據(jù)代入訓(xùn)練好的模型中,預(yù)測延長縣LSIs,并將其按照自然間斷點法劃分為5個易發(fā)性級別。

    3.2 延長滑坡易發(fā)性預(yù)測結(jié)果

    在SPSS modeler軟件中進行CHAID決策樹建模。以WOE樣本數(shù)據(jù)為例,首先需從外部源中讀取源節(jié)點,將6 806個滑坡非滑坡樣本數(shù)據(jù)導(dǎo)入SPSS modeler軟件中;接著對字段屬性、測量級別及各字段在建模中的角色進行選擇或修改;再經(jīng)由分區(qū)選擇將樣本數(shù)據(jù)分為訓(xùn)練集(70%)和測試集(30%);然后在CHAID建模節(jié)點字段選項卡中使用預(yù)定義角色,應(yīng)用boosting算法創(chuàng)建一個整體,由其生成模型序列以增強模型預(yù)測的準(zhǔn)確度;選擇CHAID樹生長算法并定制樹的最大深度值為5、父節(jié)點的最小記錄數(shù)為75、子節(jié)點的最小記錄數(shù)為15,以此來限制決策樹的增長;CHAID決策樹的其他參數(shù)使用SPSS modeler中的默認(rèn)值;最后將整體環(huán)境因子的WOE帶入訓(xùn)練好的CHAID決策樹模型中,實現(xiàn)延長縣滑坡LSIs的準(zhǔn)確預(yù)測。FR-CHAID和單獨CHAID決策樹模型的建模步驟和參數(shù)設(shè)置與WOE-CHAID決策樹模型基本一致。

    3.3 滑坡易發(fā)性制圖表達

    分兩步開展滑坡易發(fā)性制圖,首先將3種耦合模型預(yù)測出的LSIs導(dǎo)入GIS軟件中,然后依據(jù)自然間斷點法將延長縣滑坡易發(fā)性劃分為極高、高、中等、低和極低5類等級區(qū)間。WOE-CHAID、FR-CHAID和單獨CHAID決策樹模型下的滑坡易發(fā)性結(jié)果如圖4所示。延長縣大部分地區(qū)屬于低和極低易發(fā)區(qū),滑坡高和極高易發(fā)區(qū)主要位于坡度和高程中等且距離河流較近的山地丘陵地區(qū)。但3種耦合模型下得到的滑坡易發(fā)性級別存在顯著差異,圖4中延長縣內(nèi)已發(fā)生的82處滑坡幾乎都落在WOE-CHAID和FR-CHAID決策樹模型預(yù)測的極高與高易發(fā)性等級區(qū)域內(nèi),而單獨CHAID決策樹模型預(yù)測的極高與高易發(fā)性等級區(qū)域與82處滑坡位置存在些許偏差。

    4 滑坡易發(fā)性預(yù)測不確定性分析

    4.1 ROC精度評價

    采用測試集AUC值作為具體指標(biāo)量化不同耦合模型的預(yù)測性能,AUC值越大,表明耦合模型預(yù)測性能越優(yōu)。WOE-CHAID、FR-CHAID和單獨CHAID決策樹模型的滑坡易發(fā)性結(jié)果ROC曲線如圖5所示。從圖5中可知,3種耦合模型下的結(jié)果均較好且相對穩(wěn)定,表現(xiàn)出良好的滑坡易發(fā)性性能。AUC精度從大到小依次為:AUC(WOE-CHAID)>AUC(FR-CHAID)>AUC(單獨的CHAID),說明FR和WOE兩種關(guān)聯(lián)分析法在CHAID決策樹模型中具有比原始因子數(shù)據(jù)更穩(wěn)定的易發(fā)性預(yù)測性能。WOE耦合CHAID決策樹模型的易發(fā)性預(yù)測效果最好且預(yù)測效率最高,AUC精度較FR提高了2.1%,較原始因子數(shù)據(jù)提高了3.1%。

    4.2 滑坡易發(fā)性指數(shù)分布規(guī)律

    采用均值和標(biāo)準(zhǔn)差分別反映LSIs分布的平均水平和離散程度,并以此分析耦合模型下的易發(fā)性預(yù)測不確定性。WOE-CHAID、FR-CHAID和單獨CHAID決策樹模型預(yù)測的LSIs分布不確定性規(guī)律較為一致,在極低和低易發(fā)區(qū)分布較集中而在高和極高易發(fā)區(qū)分布逐漸減少。LSIs平均值從小到大排名為:單獨的CHAID (0.364)<FR-CHAID(0.385)<WOE-CHAID(0.399);標(biāo)準(zhǔn)差從小到大排名為:FR-CHAID(0.178)<WOE-CHAID(0.219)<單獨的CHAID(0.228)(圖6)。其中,WOE-CHAID決策樹模型預(yù)測的LSIs平均值較小,標(biāo)準(zhǔn)差較大,表明其對研究區(qū)內(nèi)的滑坡均有較好

    的識別能力,另外,對滑坡易發(fā)性的區(qū)分度也較好,能很好地反映出不同柵格單元內(nèi)LSIs的差異,用較少的高LSIs反映盡可能多的滑坡編錄信息。結(jié)合AUC精度發(fā)現(xiàn),WOE-CHAID決策樹模型預(yù)測效果略好于FR-CHAID決策樹模型,單獨的CHAID決策樹模型預(yù)測效果最差。

    4.3 耦合模型預(yù)測易發(fā)性指數(shù)的差異性分析

    采用顯著性差異水平來進一步分析各耦合模型下易發(fā)性建模的不確定性,通過該試驗計算各耦合模型下預(yù)測的LSIs的平均秩,以便對易發(fā)性模型性能排序。平均秩越小則模型性能越好,最終模型比較結(jié)果為:WOE-CHAID決策樹模型預(yù)測LSIs的平均秩(值為1.85)最小,其次是FR-CHAID(值為2.06) 和單獨的CHAID決策樹(值為2.09)模型。顯著性差異水平和平均秩顯示出各耦合模型的易發(fā)性建模存在不確定性,如何規(guī)避這些不確定性是獲得可靠的易發(fā)性模型的重要研究內(nèi)容。

    4.4 滑坡環(huán)境因子重要性分析

    滑坡環(huán)境因子的重要性反映了已發(fā)生的滑坡事件受該環(huán)境因子影響程度的大小。由于原始因子數(shù)據(jù)和不同的關(guān)聯(lián)分析值在易發(fā)性預(yù)測建模中有著不同的表現(xiàn),基于CHAID決策樹模型中自帶的分類器屬性來評估在原始因子數(shù)據(jù)、FR和WOE等輸入變量下各個環(huán)境因子的重要性。另外,易發(fā)性建模中共使用14個環(huán)境因子(原始因子數(shù)據(jù)含“啞變量”類型,共23個環(huán)境因子),排名10名之后的環(huán)境因子重要性均小于0.04,因此僅展示重要性排名前10的環(huán)境因子。從圖7可知,坡度、地形起伏度、距河流的距離(原始因子數(shù)據(jù)中為100~300 m和500~800 m的兩個“啞變量”因子)、地形切割深度和地形粗糙度等5個環(huán)境因子在單獨CHAID、FR-CHAID和WOE-CHAID決策樹易發(fā)性預(yù)測中有著較大的貢獻,占據(jù)重要性排名均在前5位,重要性均大于0.08。其次,平面曲率和地形濕度指數(shù)在所有決策樹模型中也發(fā)揮著相對重要的作用,重要性均大于0.04。

    4.5 各關(guān)聯(lián)分析法的性能分析

    關(guān)聯(lián)分析法通過定量統(tǒng)計可直觀表現(xiàn)各環(huán)境因

    子不同屬性區(qū)間對滑坡易發(fā)性空間的影響性。Li等、Saha等對上述部分關(guān)聯(lián)分析法反映滑坡與其環(huán)境因子空間關(guān)聯(lián)的性能進行了對比分析,所得結(jié)果與筆者研究基本一致。由上述分析可知,環(huán)境因子與滑坡間的空間信息的關(guān)聯(lián)性表達越充分,則LSIs的區(qū)分度越大,進一步的易發(fā)性預(yù)測效果就越佳。在FR和WOE關(guān)聯(lián)分析法的環(huán)境因子分級中,WOE更能反映環(huán)境因子內(nèi)部影響滑坡發(fā)育的空間信息的差異,具有更優(yōu)的預(yù)測精度(AUC=86.3%);FR相較于WOE法更加簡潔高效,在保證易發(fā)性精度的同時能有效避免太復(fù)雜的統(tǒng)計分析;基于原始因子數(shù)據(jù)進行的單獨CHAID決策樹模型易發(fā)性預(yù)測精度略小于FR-CHAID和WOE-CHAID決策樹模型。此外,單獨的CHAID、FR-CHAID和WOE-CHAID決策樹模型預(yù)測的LSIs平均值逐漸減小而標(biāo)準(zhǔn)差逐漸增大,且平均秩也逐漸減小??梢婈P(guān)聯(lián)分析法的易發(fā)性預(yù)測建模效果較好,WOE優(yōu)于FR,而原始因子數(shù)據(jù)的易發(fā)性建模效果較差。

    由文獻[27,42]可知,滑坡與環(huán)境因子(不考慮誘發(fā)因子)之間的非線性關(guān)聯(lián)分析法種類繁多。筆者僅使用FR和WOE兩種關(guān)聯(lián)分析法耦合CHAID決策樹模型進行滑坡易發(fā)性的不確定性對比分析而并未考慮其他關(guān)聯(lián)分析法,在下一步研究中可以考慮使用概率法、信息量、確定性系數(shù)和熵指數(shù)等其他關(guān)聯(lián)分析法,耦合多種不同類型的模型開展更加全面的易發(fā)性預(yù)測不確定性分析。

    5 結(jié)論

    1)WOE-CHAID決策樹模型易發(fā)性預(yù)測的AUC精度最高,且均值和平均秩較小,標(biāo)準(zhǔn)差較大;FR-CHAID決策樹的AUC精度略低于WOE-CHAID,可見WOE具有更優(yōu)秀的非線性關(guān)聯(lián)性能。

    2)將原始因子直接用作輸入變量的單獨CHAID決策樹模型的易發(fā)性預(yù)測精度整體略低于關(guān)聯(lián)分析法的耦合模型。為了提高滑坡易發(fā)性建模效率,可直接使用單獨CHAID決策樹模型,但要體現(xiàn)滑坡與其環(huán)境因子的空間關(guān)聯(lián)性或分析環(huán)境因子各子區(qū)間對滑坡發(fā)育的影響規(guī)律,則使用關(guān)聯(lián)分析法和CHAID決策樹模型耦合建模的優(yōu)勢顯著。

    3)總體來說,WOE-CHAID決策樹模型的易發(fā)性預(yù)測結(jié)果可靠性最高,預(yù)測出的LSIs與實際的滑坡概率分布特征更加相符。

    參考文獻:

    [1] 黃發(fā)明, 殷坤龍, 蔣水華, 等. 基于聚類分析和支持向量機的滑坡易發(fā)性評價[J]. 巖石力學(xué)與工程學(xué)報, 2018, 37(1): 156-167.

    HUANG F M, YIN K L, JIANG S H, et al. Landslide susceptibility assessment based on clustering analysis and support vector machine [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 156-167. (in Chinese)

    [2] 黃發(fā)明, 葉舟, 姚池, 等. 滑坡易發(fā)性預(yù)測不確定性:環(huán)境因子不同屬性區(qū)間劃分和不同數(shù)據(jù)驅(qū)動模型的影響[J]. 地球科學(xué), 2020, 45(12): 4535-4549.

    HUANG F M, YE Z, YAO C, et al. Uncertainties of landslide susceptibility prediction: Different attribute interval divisions of environmental factors and different data-based models [J]. Earth Science, 2020, 45(12): 4535-4549. (in Chinese)

    [3] 李利峰, 張曉虎, 鄧慧琳, 等. 基于熵指數(shù)與邏輯回歸耦合模型的滑坡災(zāi)害易發(fā)性評價: 以藍田縣為例[J]. 科學(xué)技術(shù)與工程, 2020, 20(14): 5536-5543.

    LI L F, ZHANG X H, DENG H L, et al. Assessment of landslide susceptibility based on coupling model of index of entropy and logistic regression: A case study of Lantian County [J]. Science Technology and Engineering, 2020, 20(14): 5536-5543. (in Chinese)

    [4] HUANG F M, CAO Z S, JIANG S H, et al. Landslide susceptibility prediction based on a semi-supervised multiple-layer perceptron model [J]. Landslides, 2020, 17(12): 2919-2930.

    [5] HUANG F M, CAO Z S, GUO J F, et al. Comparisons of heuristic, general statistical and machine learning models for landslide susceptibility prediction and mapping [J]. CATENA, 2020, 191: 104580.

    [6] 楊光, 徐佩華, 曹琛, 等. 基于確定性系數(shù)組合模型的區(qū)域滑坡敏感性評價[J]. 工程地質(zhì)學(xué)報, 2019, 27(5): 1153-1163.

    YANG G, XU P H, CAO C, et al. Assessment of regional landslide susceptibility based on combined model of certainty factor method [J]. Journal of Engineering Geology, 2019, 27(5): 1153-1163. (in Chinese)

    [7] 許嘉慧, 孫德亮, 王月, 等. 基于GIS與改進層次分析法的奉節(jié)縣滑坡易發(fā)性區(qū)劃[J]. 重慶師范大學(xué)學(xué)報(自然科學(xué)版), 2020, 37(2): 36-44, 2, 142.

    XU J H, SUN D L, WANG Y, et al. Landslide susceptibility mapping of Fengjie County based on GIS and improved analytic hierarchy process [J]. Journal of Chongqing Normal University (Natural Science), 2020, 37(2): 36-44, 2, 142. (in Chinese)

    [8] 黃發(fā)明, 陳佳武, 唐志鵬, 等. 不同空間分辨率和訓(xùn)練測試集比例下的滑坡易發(fā)性預(yù)測不確定性[J]. 巖石力學(xué)與工程學(xué)報, 2021, 40(6): 1155-1169.

    HUANG F M, CHEN J W, TANG Z P, et al. Uncertainties of landslide susceptibility prediction due to different spatial resolutions and different proportions of training and testing datasets [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1155-1169. (in Chinese)

    [9] 楊永剛, 殷坤龍, 趙海燕, 等. 基于C5.0決策樹-快速聚類模型的萬州區(qū)庫岸段鄉(xiāng)鎮(zhèn)滑坡易發(fā)性區(qū)劃[J]. 地質(zhì)科技情報, 2019, 38(6): 189-197.

    YANG Y G, YIN K L, ZHAO H Y, et al. Landslide susceptibility evaluation for township units of bank section in Wanzhou district based on C5.0 decision tree and K-means cluster model [J]. Geological Science and Technology Information, 2019, 38(6): 189-197. (in Chinese)

    [10] 田乃滿, 蘭恒星, 伍宇明, 等. 人工神經(jīng)網(wǎng)絡(luò)和決策樹模型在滑坡易發(fā)性分析中的性能對比[J]. 地球信息科學(xué)學(xué)報, 2020, 22(12): 2304-2316.

    TIAN N M, LAN H X, WU Y M, et al. Performance comparison of BP artificial neural network and CART decision tree model in landslide susceptibility prediction [J]. Journal of Geo-Information Science, 2020, 22(12): 2304-2316. (in Chinese)

    [11] 喬世范, 王超. 基于遺傳模擬退火算法的滑坡位移預(yù)測方法 [J]. 土木與環(huán)境工程學(xué)報(中英文), 2021, 43(1): 25-35.

    QIAO S F, WANG C. Landslide displacement prediction based on the Genetic Simulated Annealing algorithm[J]. Journal of Civil and Environmental Engineering, 2021, 43(1): 25-35. (in Chinese)

    [12] 吳潤澤, 胡旭東, 梅紅波, 等. 基于隨機森林的滑坡空間易發(fā)性評價:以三峽庫區(qū)湖北段為例[J]. 地球科學(xué), 2021, 46(1): 321-330.

    WU R Z, HU X D, MEI H B, et al. Spatial susceptibility assessment of landslides based on random forest: A case study from Hubei section in the Three Gorges reservoir area [J]. Earth Science, 2021, 46(1): 321-330. (in Chinese)

    [13] 牛瑞卿, 彭令, 葉潤青, 等. 基于粗糙集的支持向量機滑坡易發(fā)性評價[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2012, 42(2): 430-439.

    NIU R Q, PENG L, YE R Q, et al. Landslide susceptibility assessment based on rough sets and support vector machine [J]. Journal of Jilin University (Earth Science Edition), 2012, 42(2): 430-439. (in Chinese)

    [14] ALTHUWAYNEE O F, PRADHAN B, AHMAD N. Landslide susceptibility mapping using decision-tree based Chi-squared automatic interaction detection (CHAID) and Logistic regression (LR) integration [J]. IOP Conference Series: Earth and Environmental Science, 2014, 20: 012032.

    [15] CHEN W W, ZHANG S. GIS-based comparative study of Bayes network, Hoeffding tree and logistic model tree for landslide susceptibility modeling [J]. CATENA, 2021, 203: 105344.

    [16] DEMIR G. GIS-based landslide susceptibility mapping for a part of the North Anatolian Fault Zone between Re?adiye and Koyulhisar (Turkey) [J]. CATENA, 2019, 183: 104211.

    [17] 張玘愷, 凌斯祥, 李曉寧, 等. 九寨溝縣滑坡災(zāi)害易發(fā)性快速評估模型對比研究[J]. 巖石力學(xué)與工程學(xué)報, 2020, 39(8): 1595-1610.

    ZHANG Q K, LING S X, LI X N, et al. Comparison of landslide susceptibility mapping rapid assessment models in Jiuzhaigou County, Sichuan Province, China [J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(8): 1595-1610. (in Chinese)

    [18] HUANG F M, ZHANG J, ZHOU C B, et al. A deep learning algorithm using a fully connected sparse autoencoder neural network for landslide susceptibility prediction [J]. Landslides, 2020, 17(1): 217-229.

    [19] 劉璐瑤, 高惠瑛. 基于證據(jù)權(quán)與Logistic回歸模型耦合的滑坡易發(fā)性評價[J/OL]. 工程地質(zhì)學(xué)報. https://doi.org/10.13544/j.cnki.jeg.2020-482.

    LIU L Y, GAO H Y. Landslide susceptibility assessment based on coupling of WOE model and Logistic regression model [J/OL]. Journal of Engineering Geology. https://doi.org/10.13544/j.cnki.jeg.2020-482. (in Chinese)

    [20] 徐勝華, 劉紀(jì)平, 王想紅, 等. 熵指數(shù)融入支持向量機的滑坡災(zāi)害易發(fā)性評價方法: 以陜西省為例[J]. 武漢大學(xué)學(xué)報·信息科學(xué)版, 2020, 45(8): 1214-1222.

    XU S H, LIU J P, WANG X H, et al. Landslide susceptibility assessment method incorporating index of entropy based on support vector machine: A case study of Shaanxi Province [J]. Geomatics and Information Science of Wuhan University, 2020, 45(8): 1214-1222. (in Chinese)

    [21] 楊華陽, 許向?qū)帲?楊鴻發(fā). 基于證據(jù)權(quán)法的九寨溝地震滑坡危險性評價[J]. 中國地質(zhì)災(zāi)害與防治學(xué)報, 2020, 31(3): 20-29.

    YANG H Y, XU X N, YANG H F. The Jiuzhaigou co-seismic landslide hazard assessment based on weight of evidence method [J]. The Chinese Journal of Geological Hazard and Control, 2020, 31(3): 20-29. (in Chinese)

    [22] ZHANG T Y, HAN L, CHEN W, et al. Hybrid integration approach of entropy with logistic regression and support vector machine for landslide susceptibility modeling [J]. Entropy, 2018, 20(11): 884.

    [23] 李文彬, 范宣梅, 黃發(fā)明, 等. 不同環(huán)境因子聯(lián)接方法和數(shù)據(jù)驅(qū)動模型對滑坡易發(fā)性預(yù)測建模的影響規(guī)律[J/OL]. 地球科學(xué). https://kns.cnki.net/kcms/detail/42.1874.P.20210506.1457.004.html.

    LI W B, FAN X M, HUANG F M, et al. Influence law of different environmental factor connection methods and data-based models on landslide susceptibility prediction modeling [J/OL]. Earth Science. https://kns.cnki.net/kcms/detail/42.1874.P.20210506.1457.004.html. (in Chinese)

    [24] 張鐘遠, 鄧明國, 徐世光, 等. 鎮(zhèn)康縣滑坡易發(fā)性評價模型對比研究 [J]. 巖石力學(xué)與工程學(xué)報, 2022, 41(1):157-171.

    ZHANG Z Y, DENG M G, XU S G, et al. Comparison of landslide susceptibility assessment models in Zhenkang County, Yunnan Province, China [J/OL]. Chinese Journal of Rock Mechanics and Engineering. 2022, 41(1):157-171. (in Chinese)

    [25] 郭子正, 殷坤龍, 黃發(fā)明, 等. 基于滑坡分類和加權(quán)頻率比模型的滑坡易發(fā)性評價[J]. 巖石力學(xué)與工程學(xué)報, 2019, 38(2): 287-300.

    GUO Z Z, YIN K L, HUANG F M, et al. Evaluation of landslide susceptibility based on landslide classification and weighted frequency ratio model [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(2): 287-300. (in Chinese)

    [26] KONTOES C, LOUPASAKIS C, PAPOUTSIS I, et al. Landslide susceptibility mapping of central and western Greece, combining NGI and WoE methods, with remote sensing and ground truth data [J]. Land, 2021, 10(4): 402.

    [27] LI W B, FAN X M, HUANG F M, et al. Uncertainties analysis of collapse susceptibility prediction based on remote sensing and GIS: Influences of different data-based models and connections between collapses and environmental factors [J]. Remote Sensing, 2020, 12(24): 4134.

    [28] 車文超, 秦勝伍, 苗強, 等. 滑坡敏感性評價中因子分類方法的研究 [J]. 工程地質(zhì)學(xué)報, 2020, 28(Sup1): 116-124.

    CHE W C, QIN S W, MIAO Q, et al. Research on factor classification method of landslide susceptibility mapping [J]. [J].Journal of Engineering Geology, 2020, 28(Sup1): 116-124. (in Chinese)

    [29] 羅路廣, 裴向軍, 崔圣華, 等. 九寨溝地震滑坡易發(fā)性評價因子組合選取研究[J]. 巖石力學(xué)與工程學(xué)報, 2021, 40(11): 2306-2319.

    LUO L G, PEI X J, CUI S H, et al. Combined selection of susceptibility assessment factors for Jiuzhaigou earthquake-induced landslides [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(11): 2306-2319. (in Chinese)

    [30] GUO Z Z, SHI Y, HUANG F M, et al. Landslide susceptibility zonation method based on C5.0 decision tree and K-means cluster algorithms to improve the efficiency of risk management [J]. Geoscience Frontiers, 2021, 12(6): 101249.

    [31] 郭天頌, 張菊清, 韓煜, 等. 基于粒子群優(yōu)化支持向量機的延長縣滑坡易發(fā)性評價[J]. 地質(zhì)科技情報, 2019, 38(3): 236-243.

    GUO T S, ZHANG J Q, HAN Y, et al. Evaluation of landslide susceptibility in Yanchang County based on particle swarm optimization-based support vector machine [J]. Geological Science and Technology Information, 2019, 38(3): 236-243. (in Chinese)

    [32] HUANG F M, TAO S Y, CHANG Z L, et al. Efficient and automatic extraction of slope units based on multi-scale segmentation method for landslide assessments [J]. Landslides, 2021, 18(11): 3715-3731.

    [33] 沈玲玲, 劉連友, 許沖, 等. 基于多模型的滑坡易發(fā)性評價: 以甘肅岷縣地震滑坡為例[J]. 工程地質(zhì)學(xué)報, 2016, 24(1): 19-28.

    SHEN L L, LIU L Y, XU C, et al. Multi-models based landslide susceptibility evaluation: Illustrated with landslides triggered by Minxian earthquake [J]. Journal of Engineering Geology, 2016, 24(1): 19-28. (in Chinese)

    [34] 楊世豪, 蘇立君, 張崇磊, 等. 強降雨作用下昔格達邊坡滲流特性及穩(wěn)定性分析[J]. 土木與環(huán)境工程學(xué)報(中英文), 2020, 42(4): 19-27.

    YANG S H, SU L J, ZHANG C L, et al. Analysis of seepage characteristics and stability of Xigeda formation slope under heavy rainfall [J]. Journal of Civil and Environmental Engineering, 2020, 42(4): 19-27. (in Chinese)

    [35] 于憲煜, 胡友健, 牛瑞卿. 基于RS-SVM模型的滑坡易發(fā)性評價因子選擇方法研究[J]. 地理與地理信息科學(xué), 2016, 32(3): 23-28, 2.

    YU X Y, HU Y J, NIU R Q. Research on the method to select landslide susceptibility evaluation factors based on RS-SVM model [J]. Geography and Geo-Information Science, 2016, 32(3): 23-28, 2. (in Chinese)

    [36] 吳雨辰, 周晗旭, 車愛蘭. 基于粗糙集神經(jīng)網(wǎng)絡(luò)的IBURI地震滑坡易發(fā)性研究[J]. 巖石力學(xué)與工程學(xué)報, 2021, 40(6): 1226-1235.

    WU Y C, ZHOU H X, CHE A L. Susceptibility of landslides caused by IBURI earthquake based on rough set-neural network [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1226-1235. (in Chinese)

    [37] 黃發(fā)明, 汪洋, 董志良, 等. 基于灰色關(guān)聯(lián)度模型的區(qū)域滑坡敏感性評價[J]. 地球科學(xué), 2019(2): 664-676.

    HUANG F M, WANG Y, DONG Z L, et al. Regional landslide susceptibility mapping based on grey relational degree model [J]. Earth Science, 2019(2): 664-676. (in Chinese)

    [38] 許英姿, 盧玉南, 李東陽, 等. 基于GIS和信息量模型的廣西花崗巖分布區(qū)滑坡易發(fā)性評價[J]. 工程地質(zhì)學(xué)報, 2016, 24(4): 693-703.

    XU Y Z, LU Y N, LI D Y, et al. GIS and information model based landslide susceptibility assessment in granite area of Guangxi Province [J]. Journal of Engineering Geology, 2016, 24(4): 693-703. (in Chinese)

    [39] 解明禮, 巨能攀, 趙建軍, 等. 區(qū)域地質(zhì)災(zāi)害易發(fā)性分級方法對比分析研究[J]. 武漢大學(xué)學(xué)報·信息科學(xué)版, 2021, 46(7): 1003-1014.

    XIE M L, JU N P, ZHAO J J, et al. Comparative analysis on classification methods of geological disaster susceptibility assessment [J]. Geomatics and Information Science of Wuhan University, 2021, 46(7): 1003-1014. (in Chinese)

    [40] ZHU L, WANG G J, HUANG F M, et al. Landslide susceptibility prediction using sparse feature extraction and machine learning models based on GIS and remote sensing [J]. IEEE Geoscience and Remote Sensing Letters, 4029, PP(99): 1-5.

    [41] SAHA A, SAHA S. Comparing the efficiency of weight of evidence, support vector machine and their ensemble approaches in landslide susceptibility modelling: A study on Kurseong region of Darjeeling Himalaya, India [J]. Remote Sensing Applications: Society and Environment, 2020, 19: 100323.

    [42] SANDRIC I, IONITA C, CHITU Z, et al. Using CUDA to accelerate uncertainty propagation modelling for landslide susceptibility assessment [J]. Environmental Modelling & Software, 2019, 115: 176-186.

    (編輯 胡玲)

    亚洲七黄色美女视频| 精品国产一区二区三区四区第35| 亚洲成av片中文字幕在线观看| 男女免费视频国产| 欧美日韩视频精品一区| 校园春色视频在线观看| 热99国产精品久久久久久7| 国产精品综合久久久久久久免费 | 国产单亲对白刺激| 欧美精品av麻豆av| 国内久久婷婷六月综合欲色啪| 久久国产亚洲av麻豆专区| 少妇 在线观看| av欧美777| 午夜亚洲福利在线播放| 很黄的视频免费| 久久 成人 亚洲| 精品国产亚洲在线| 国产蜜桃级精品一区二区三区 | 黄网站色视频无遮挡免费观看| 欧美av亚洲av综合av国产av| 人人妻人人澡人人爽人人夜夜| 欧美日韩福利视频一区二区| 久久 成人 亚洲| 悠悠久久av| 一本一本久久a久久精品综合妖精| av超薄肉色丝袜交足视频| 国产精品亚洲av一区麻豆| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 亚洲久久久国产精品| 18禁裸乳无遮挡免费网站照片 | 精品亚洲成a人片在线观看| 90打野战视频偷拍视频| 国产97色在线日韩免费| 久久久久国产精品人妻aⅴ院 | 久久午夜综合久久蜜桃| 亚洲精品国产区一区二| 亚洲久久久国产精品| 91国产中文字幕| av天堂久久9| 一级a爱视频在线免费观看| 精品国产一区二区三区四区第35| 免费看十八禁软件| 国产片内射在线| 亚洲精品一卡2卡三卡4卡5卡| av欧美777| cao死你这个sao货| 高清av免费在线| 国产精品乱码一区二三区的特点 | 老熟女久久久| 免费高清在线观看日韩| 久久久久久亚洲精品国产蜜桃av| netflix在线观看网站| 9热在线视频观看99| 变态另类成人亚洲欧美熟女 | 身体一侧抽搐| 亚洲自偷自拍图片 自拍| 午夜成年电影在线免费观看| 午夜精品在线福利| 国产av精品麻豆| 亚洲伊人色综图| 又黄又粗又硬又大视频| 老司机午夜十八禁免费视频| 99国产精品免费福利视频| av中文乱码字幕在线| 日韩成人在线观看一区二区三区| 两个人看的免费小视频| 母亲3免费完整高清在线观看| 成年动漫av网址| 国产亚洲精品久久久久久毛片 | 午夜福利视频在线观看免费| 看免费av毛片| 国产熟女午夜一区二区三区| 国产精华一区二区三区| 丝袜在线中文字幕| 国产亚洲精品第一综合不卡| 欧美黑人欧美精品刺激| 麻豆av在线久日| 操出白浆在线播放| 国产成人免费观看mmmm| 热re99久久国产66热| 天堂俺去俺来也www色官网| 天天影视国产精品| 90打野战视频偷拍视频| 欧美日韩福利视频一区二区| 男女高潮啪啪啪动态图| 国产男女内射视频| 美女扒开内裤让男人捅视频| 淫妇啪啪啪对白视频| 一二三四在线观看免费中文在| 首页视频小说图片口味搜索| 精品午夜福利视频在线观看一区| 亚洲avbb在线观看| 国产男女内射视频| 欧美乱码精品一区二区三区| 老司机深夜福利视频在线观看| 欧美国产精品va在线观看不卡| 91九色精品人成在线观看| www日本在线高清视频| 熟女少妇亚洲综合色aaa.| 亚洲精华国产精华精| 国产午夜精品久久久久久| 男女免费视频国产| 国产极品粉嫩免费观看在线| 搡老乐熟女国产| 女性被躁到高潮视频| 嫩草影视91久久| 亚洲视频免费观看视频| 亚洲av第一区精品v没综合| 国产精品美女特级片免费视频播放器 | av视频免费观看在线观看| 亚洲片人在线观看| 在线永久观看黄色视频| 女人爽到高潮嗷嗷叫在线视频| 免费一级毛片在线播放高清视频 | 又紧又爽又黄一区二区| 午夜影院日韩av| 久久中文看片网| 日韩视频一区二区在线观看| 大陆偷拍与自拍| 老司机靠b影院| 啦啦啦免费观看视频1| 人人妻人人爽人人添夜夜欢视频| 成人亚洲精品一区在线观看| 别揉我奶头~嗯~啊~动态视频| 黄色a级毛片大全视频| 国产乱人伦免费视频| 热99久久久久精品小说推荐| 一区二区三区精品91| 久久国产精品人妻蜜桃| 麻豆成人av在线观看| 精品人妻在线不人妻| 人成视频在线观看免费观看| 国产精品久久电影中文字幕 | 侵犯人妻中文字幕一二三四区| 亚洲全国av大片| 日韩三级视频一区二区三区| 一级片免费观看大全| 国产亚洲欧美精品永久| 香蕉久久夜色| 亚洲成人免费av在线播放| 人人妻人人爽人人添夜夜欢视频| 欧美日韩乱码在线| 亚洲五月婷婷丁香| 高清av免费在线| av电影中文网址| 亚洲国产欧美日韩在线播放| 欧美在线黄色| 国产精品免费一区二区三区在线 | 久久久久精品国产欧美久久久| 国产亚洲精品一区二区www | 在线观看免费午夜福利视频| 黄色丝袜av网址大全| www日本在线高清视频| 午夜影院日韩av| 在线免费观看的www视频| 国产成人精品久久二区二区91| 中文字幕色久视频| 黄片播放在线免费| 亚洲欧美精品综合一区二区三区| 伊人久久大香线蕉亚洲五| 在线观看免费视频网站a站| 91成年电影在线观看| 丝袜人妻中文字幕| 少妇的丰满在线观看| 99国产精品99久久久久| 日韩免费高清中文字幕av| 亚洲欧美精品综合一区二区三区| 日韩 欧美 亚洲 中文字幕| 日韩视频一区二区在线观看| 嫁个100分男人电影在线观看| 午夜精品久久久久久毛片777| 精品欧美一区二区三区在线| 国产人伦9x9x在线观看| 美女福利国产在线| www日本在线高清视频| 亚洲国产欧美日韩在线播放| 免费少妇av软件| 人妻一区二区av| 免费观看人在逋| 丝袜人妻中文字幕| 国产亚洲欧美98| 久久ye,这里只有精品| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品国产av在线观看| 亚洲熟妇中文字幕五十中出 | 国产精品.久久久| 亚洲一码二码三码区别大吗| tocl精华| 操美女的视频在线观看| 亚洲精品国产区一区二| 少妇被粗大的猛进出69影院| 国产一区二区三区综合在线观看| 90打野战视频偷拍视频| 99国产精品一区二区蜜桃av | 久久久精品区二区三区| 宅男免费午夜| 亚洲五月婷婷丁香| 久久影院123| 建设人人有责人人尽责人人享有的| 国产一区二区三区综合在线观看| 两人在一起打扑克的视频| 色婷婷久久久亚洲欧美| 18禁国产床啪视频网站| 热99re8久久精品国产| 老司机亚洲免费影院| 国产欧美日韩一区二区三| 国产欧美日韩一区二区精品| 午夜福利影视在线免费观看| 国产高清激情床上av| 波多野结衣一区麻豆| 国产精品 欧美亚洲| 丰满饥渴人妻一区二区三| 久久香蕉国产精品| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区mp4| 成年人免费黄色播放视频| 老司机午夜十八禁免费视频| 777久久人妻少妇嫩草av网站| 精品视频人人做人人爽| 午夜精品在线福利| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 无人区码免费观看不卡| 亚洲一区中文字幕在线| 天堂√8在线中文| 亚洲一区高清亚洲精品| 国产亚洲精品久久久久久毛片 | www.精华液| 亚洲成人免费电影在线观看| 欧美激情 高清一区二区三区| 欧美日韩视频精品一区| 91在线观看av| 亚洲成人免费av在线播放| 啪啪无遮挡十八禁网站| 国产在线精品亚洲第一网站| 亚洲专区中文字幕在线| 久久中文字幕一级| 法律面前人人平等表现在哪些方面| 国产亚洲精品久久久久久毛片 | 看免费av毛片| 日韩欧美国产一区二区入口| 亚洲av第一区精品v没综合| av一本久久久久| 久久久久久人人人人人| 国产一区二区三区在线臀色熟女 | 精品久久久久久久久久免费视频 | 97人妻天天添夜夜摸| 人妻一区二区av| 亚洲人成伊人成综合网2020| 一边摸一边抽搐一进一出视频| 亚洲精品久久成人aⅴ小说| 久久性视频一级片| 色综合婷婷激情| 国产淫语在线视频| 国产精品偷伦视频观看了| 夜夜躁狠狠躁天天躁| 少妇 在线观看| 亚洲熟妇熟女久久| 精品欧美一区二区三区在线| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| 国产色视频综合| 黄色怎么调成土黄色| 身体一侧抽搐| 视频区欧美日本亚洲| 91国产中文字幕| 国产精品一区二区免费欧美| 热re99久久精品国产66热6| 丁香欧美五月| a级毛片在线看网站| 免费观看精品视频网站| 国产在线一区二区三区精| 欧美激情高清一区二区三区| 久久久久国产一级毛片高清牌| 少妇 在线观看| 久久天堂一区二区三区四区| 国产精品成人在线| 亚洲色图综合在线观看| 精品电影一区二区在线| 性色av乱码一区二区三区2| 法律面前人人平等表现在哪些方面| 97人妻天天添夜夜摸| av视频免费观看在线观看| 这个男人来自地球电影免费观看| 成年女人毛片免费观看观看9 | 国产成人系列免费观看| 天天操日日干夜夜撸| 欧美日韩亚洲国产一区二区在线观看 | 韩国av一区二区三区四区| 亚洲国产精品sss在线观看 | 久久久久久久午夜电影 | 亚洲精品自拍成人| 激情视频va一区二区三区| a级片在线免费高清观看视频| 男人的好看免费观看在线视频 | 桃红色精品国产亚洲av| 国产精品偷伦视频观看了| 啦啦啦视频在线资源免费观看| 久久久精品免费免费高清| 国产精品秋霞免费鲁丝片| 97人妻天天添夜夜摸| 日韩熟女老妇一区二区性免费视频| 老司机福利观看| 国产亚洲一区二区精品| 国产男女超爽视频在线观看| 黑人猛操日本美女一级片| 三上悠亚av全集在线观看| 人成视频在线观看免费观看| 色播在线永久视频| 在线永久观看黄色视频| 在线av久久热| 免费在线观看影片大全网站| 欧美日韩乱码在线| 我的亚洲天堂| 极品少妇高潮喷水抽搐| 亚洲午夜精品一区,二区,三区| av天堂在线播放| 欧美日韩成人在线一区二区| 在线观看免费高清a一片| 99re在线观看精品视频| 国产精品成人在线| 色播在线永久视频| 在线观看免费视频网站a站| 国产极品粉嫩免费观看在线| 亚洲精品一卡2卡三卡4卡5卡| 国产男女超爽视频在线观看| 91大片在线观看| 精品亚洲成a人片在线观看| 最近最新免费中文字幕在线| 精品久久久久久久久久免费视频 | 成年人午夜在线观看视频| 午夜福利,免费看| 高清av免费在线| 乱人伦中国视频| 午夜福利在线观看吧| 欧美日韩亚洲国产一区二区在线观看 | 成人黄色视频免费在线看| 美女福利国产在线| 精品国产一区二区三区四区第35| 午夜福利在线观看吧| 久久热在线av| 怎么达到女性高潮| 男女高潮啪啪啪动态图| 欧美成狂野欧美在线观看| 十八禁人妻一区二区| 亚洲熟女精品中文字幕| 久久精品国产99精品国产亚洲性色 | 精品福利观看| 久久久久久久国产电影| av有码第一页| 日本一区二区免费在线视频| 国产aⅴ精品一区二区三区波| 中文字幕精品免费在线观看视频| 少妇 在线观看| 国产又爽黄色视频| 久久久国产成人免费| 亚洲一区中文字幕在线| 国产欧美亚洲国产| √禁漫天堂资源中文www| 欧美日韩黄片免| 一个人免费在线观看的高清视频| 777米奇影视久久| 久久久久久久国产电影| 久久天躁狠狠躁夜夜2o2o| 老司机午夜十八禁免费视频| 国产一区二区三区视频了| 人成视频在线观看免费观看| 91在线观看av| 一二三四社区在线视频社区8| 亚洲avbb在线观看| 精品福利观看| 国产成人精品在线电影| 十分钟在线观看高清视频www| 少妇粗大呻吟视频| 欧美在线一区亚洲| 窝窝影院91人妻| 久久国产乱子伦精品免费另类| 欧美乱妇无乱码| 精品久久久精品久久久| 1024香蕉在线观看| 久久性视频一级片| 久久精品国产综合久久久| 亚洲国产精品sss在线观看 | 精品久久久久久久毛片微露脸| 黑人猛操日本美女一级片| 午夜免费成人在线视频| 大香蕉久久成人网| 一本大道久久a久久精品| 亚洲,欧美精品.| 操美女的视频在线观看| 一区二区三区激情视频| 黄色视频,在线免费观看| 欧美久久黑人一区二区| 午夜影院日韩av| 成人手机av| 搡老熟女国产l中国老女人| 久久中文看片网| 久久久国产成人精品二区 | 超碰成人久久| 一a级毛片在线观看| 香蕉国产在线看| 午夜日韩欧美国产| 亚洲第一青青草原| 免费在线观看亚洲国产| 成年人午夜在线观看视频| 99热只有精品国产| 99热网站在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品免费一区二区三区在线 | 夜夜爽天天搞| 黄色丝袜av网址大全| 757午夜福利合集在线观看| 亚洲国产欧美日韩在线播放| 他把我摸到了高潮在线观看| 黄色视频,在线免费观看| 日本vs欧美在线观看视频| 亚洲熟女毛片儿| 日本a在线网址| 欧美乱码精品一区二区三区| 天堂动漫精品| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| a在线观看视频网站| 国产精品免费视频内射| 夜夜爽天天搞| 色婷婷久久久亚洲欧美| 露出奶头的视频| 国产精品一区二区在线不卡| 国产淫语在线视频| 日韩有码中文字幕| 777久久人妻少妇嫩草av网站| 国产精品.久久久| 国产人伦9x9x在线观看| 亚洲片人在线观看| 国产区一区二久久| 国产精品1区2区在线观看. | 嫁个100分男人电影在线观看| 97人妻天天添夜夜摸| 大型av网站在线播放| 欧美激情久久久久久爽电影 | 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 9热在线视频观看99| 亚洲中文日韩欧美视频| 国产真人三级小视频在线观看| 桃红色精品国产亚洲av| 国产亚洲精品一区二区www | 自拍欧美九色日韩亚洲蝌蚪91| 动漫黄色视频在线观看| 极品教师在线免费播放| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 欧美丝袜亚洲另类 | 黄色 视频免费看| 在线看a的网站| 精品高清国产在线一区| 亚洲精品一二三| 精品久久蜜臀av无| 精品电影一区二区在线| 视频在线观看一区二区三区| 成人永久免费在线观看视频| 亚洲精品美女久久av网站| 高清av免费在线| 亚洲av成人不卡在线观看播放网| 亚洲aⅴ乱码一区二区在线播放 | 18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 久久精品aⅴ一区二区三区四区| 免费在线观看亚洲国产| 欧美成人午夜精品| 高清黄色对白视频在线免费看| 亚洲中文日韩欧美视频| 精品国产美女av久久久久小说| 国产亚洲欧美98| 婷婷成人精品国产| 一级片免费观看大全| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 好男人电影高清在线观看| 色精品久久人妻99蜜桃| 91九色精品人成在线观看| www.999成人在线观看| 法律面前人人平等表现在哪些方面| 波多野结衣av一区二区av| 国产成人欧美在线观看 | 国产单亲对白刺激| 在线观看一区二区三区激情| 婷婷丁香在线五月| 精品视频人人做人人爽| 国内久久婷婷六月综合欲色啪| 国产蜜桃级精品一区二区三区 | 一a级毛片在线观看| 久久久久久亚洲精品国产蜜桃av| av中文乱码字幕在线| 操美女的视频在线观看| tube8黄色片| 欧美日韩亚洲国产一区二区在线观看 | 黑人巨大精品欧美一区二区蜜桃| 日本vs欧美在线观看视频| 村上凉子中文字幕在线| 极品人妻少妇av视频| 乱人伦中国视频| 亚洲免费av在线视频| 国产99白浆流出| 嫁个100分男人电影在线观看| 国产精品免费视频内射| 中国美女看黄片| 国产精品亚洲一级av第二区| 国产精品成人在线| 妹子高潮喷水视频| 久久久国产成人免费| 中文字幕制服av| 国产极品粉嫩免费观看在线| 嫩草影视91久久| 操美女的视频在线观看| 在线视频色国产色| 色综合欧美亚洲国产小说| 美女午夜性视频免费| 国产区一区二久久| 久久精品熟女亚洲av麻豆精品| www.熟女人妻精品国产| 国产不卡av网站在线观看| 久久人人爽av亚洲精品天堂| 99精品欧美一区二区三区四区| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| 怎么达到女性高潮| 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 在线观看免费视频网站a站| 午夜久久久在线观看| 他把我摸到了高潮在线观看| 极品人妻少妇av视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲欧美精品永久| 精品人妻在线不人妻| 搡老乐熟女国产| 亚洲av成人av| 亚洲成人手机| 国产一区二区激情短视频| 王馨瑶露胸无遮挡在线观看| 欧美激情高清一区二区三区| 一边摸一边抽搐一进一出视频| 交换朋友夫妻互换小说| 国产精品1区2区在线观看. | 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 黄色a级毛片大全视频| 国产亚洲欧美在线一区二区| 成人永久免费在线观看视频| 交换朋友夫妻互换小说| 亚洲五月色婷婷综合| 青草久久国产| 亚洲精品国产区一区二| 在线观看免费日韩欧美大片| 国产高清视频在线播放一区| 免费在线观看亚洲国产| 国产亚洲精品第一综合不卡| 亚洲免费av在线视频| 搡老岳熟女国产| 校园春色视频在线观看| 在线永久观看黄色视频| 人人澡人人妻人| 亚洲久久久国产精品| 久久青草综合色| 亚洲av成人av| 一区二区三区精品91| 国产亚洲av高清不卡| 亚洲精华国产精华精| 在线观看66精品国产| a在线观看视频网站| 精品国产超薄肉色丝袜足j| av福利片在线| 老熟女久久久| 欧美黄色淫秽网站| 99re在线观看精品视频| 亚洲avbb在线观看| 两个人免费观看高清视频| 高清在线国产一区| 他把我摸到了高潮在线观看| 久久午夜综合久久蜜桃| av欧美777| 9色porny在线观看| 亚洲中文av在线| 国产亚洲精品一区二区www | 欧美日韩瑟瑟在线播放| 天天操日日干夜夜撸| 在线观看免费视频日本深夜| 一区在线观看完整版| 乱人伦中国视频| 久久久久视频综合| 亚洲精品久久成人aⅴ小说| 伦理电影免费视频| 高清视频免费观看一区二区| 日日夜夜操网爽| 一个人免费在线观看的高清视频| 亚洲一区二区三区不卡视频| 男女免费视频国产| 国产精品免费视频内射| a在线观看视频网站| 动漫黄色视频在线观看| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| x7x7x7水蜜桃| 在线观看一区二区三区激情| 夫妻午夜视频| 午夜免费观看网址| 高潮久久久久久久久久久不卡| 国产三级黄色录像| 亚洲专区字幕在线| 交换朋友夫妻互换小说| 亚洲成人免费电影在线观看|