• <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)境因子分析及其建模的可靠性??梢?,探討不同關(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中的默認值;最后將整體環(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.

    (編輯 胡玲)

    十八禁网站免费在线| 99九九在线精品视频| 777久久人妻少妇嫩草av网站| 午夜视频精品福利| 啦啦啦 在线观看视频| 亚洲精华国产精华精| 免费高清在线观看视频在线观看| 日韩制服骚丝袜av| 美女福利国产在线| 久久久久久久大尺度免费视频| 纯流量卡能插随身wifi吗| 亚洲精品一区蜜桃| 欧美变态另类bdsm刘玥| 精品久久久久久久毛片微露脸 | 亚洲 欧美一区二区三区| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品成人在线| 亚洲国产精品一区三区| 精品少妇一区二区三区视频日本电影| 免费黄频网站在线观看国产| 香蕉丝袜av| 久久精品久久久久久噜噜老黄| 国产在线免费精品| 黑丝袜美女国产一区| 久9热在线精品视频| 日韩中文字幕视频在线看片| 啦啦啦中文免费视频观看日本| 亚洲国产精品一区三区| 欧美乱码精品一区二区三区| 不卡一级毛片| 精品一品国产午夜福利视频| 国产精品久久久久久人妻精品电影 | 久热爱精品视频在线9| 韩国精品一区二区三区| 亚洲视频免费观看视频| 亚洲九九香蕉| 桃花免费在线播放| 日日夜夜操网爽| 一本大道久久a久久精品| 成年av动漫网址| 午夜福利一区二区在线看| 亚洲精品美女久久av网站| 狠狠婷婷综合久久久久久88av| 午夜老司机福利片| 欧美大码av| 午夜91福利影院| 女人被躁到高潮嗷嗷叫费观| 久久天堂一区二区三区四区| 亚洲一码二码三码区别大吗| 熟女少妇亚洲综合色aaa.| 精品一区二区三区av网在线观看 | 女警被强在线播放| 午夜福利视频精品| 国产精品久久久久久精品古装| 精品国产乱子伦一区二区三区 | 国内毛片毛片毛片毛片毛片| 久久久久久亚洲精品国产蜜桃av| 久久久久久久久免费视频了| 热99国产精品久久久久久7| 男女午夜视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 无限看片的www在线观看| a级毛片黄视频| 蜜桃国产av成人99| 久久久水蜜桃国产精品网| 国产成人精品久久二区二区91| 成年人黄色毛片网站| 国产精品 欧美亚洲| 91精品国产国语对白视频| 亚洲欧美精品综合一区二区三区| 男人舔女人的私密视频| 在线精品无人区一区二区三| 亚洲精品日韩在线中文字幕| 久久综合国产亚洲精品| 一区二区三区激情视频| 久久久水蜜桃国产精品网| 亚洲欧美精品自产自拍| 久久久国产成人免费| 日本av手机在线免费观看| 五月天丁香电影| 久久久精品区二区三区| 亚洲国产日韩一区二区| 91国产中文字幕| 久久影院123| 我要看黄色一级片免费的| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| 肉色欧美久久久久久久蜜桃| 国产欧美日韩综合在线一区二区| 一本久久精品| 人人妻,人人澡人人爽秒播| 久久久国产精品麻豆| 日韩人妻精品一区2区三区| 一进一出抽搐动态| 国产精品久久久久成人av| 欧美日韩黄片免| 日本欧美视频一区| 国产黄色免费在线视频| 亚洲avbb在线观看| 国产成人精品久久二区二区免费| 999精品在线视频| 久久精品熟女亚洲av麻豆精品| 大型av网站在线播放| tocl精华| 熟女少妇亚洲综合色aaa.| 国产片内射在线| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲中文av在线| 交换朋友夫妻互换小说| 久久国产精品人妻蜜桃| 99热全是精品| 亚洲第一av免费看| 国产伦人伦偷精品视频| 亚洲国产av新网站| 成年人午夜在线观看视频| 日本五十路高清| 日韩电影二区| 久久精品成人免费网站| 国产男女内射视频| 不卡一级毛片| 欧美激情 高清一区二区三区| 国产一区二区三区在线臀色熟女 | 水蜜桃什么品种好| 欧美97在线视频| 国产三级黄色录像| 欧美一级毛片孕妇| 欧美久久黑人一区二区| 精品一区二区三区av网在线观看 | 黑丝袜美女国产一区| 老熟妇仑乱视频hdxx| 嫩草影视91久久| 免费久久久久久久精品成人欧美视频| 成年动漫av网址| 日本av手机在线免费观看| 日本av免费视频播放| 99精品欧美一区二区三区四区| 另类亚洲欧美激情| 热99久久久久精品小说推荐| www.999成人在线观看| a级片在线免费高清观看视频| 美女高潮喷水抽搐中文字幕| 国产精品一区二区在线观看99| 日韩欧美一区视频在线观看| 久久 成人 亚洲| 狠狠狠狠99中文字幕| 国产成人系列免费观看| 国产精品久久久人人做人人爽| 中文字幕人妻丝袜制服| 色婷婷av一区二区三区视频| www日本在线高清视频| av线在线观看网站| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| 亚洲精品国产av成人精品| 国产无遮挡羞羞视频在线观看| 国产91精品成人一区二区三区 | 男女下面插进去视频免费观看| 国产极品粉嫩免费观看在线| 丝袜脚勾引网站| 在线观看人妻少妇| www.av在线官网国产| 最黄视频免费看| 男女国产视频网站| 水蜜桃什么品种好| 日本精品一区二区三区蜜桃| 制服诱惑二区| 一本大道久久a久久精品| 三上悠亚av全集在线观看| 国产又色又爽无遮挡免| 热99re8久久精品国产| 亚洲av电影在线观看一区二区三区| 咕卡用的链子| 久久亚洲精品不卡| 电影成人av| 日韩制服骚丝袜av| 不卡av一区二区三区| 自线自在国产av| 热99久久久久精品小说推荐| 日韩一卡2卡3卡4卡2021年| 亚洲欧美日韩高清在线视频 | 亚洲第一av免费看| 老鸭窝网址在线观看| 狂野欧美激情性xxxx| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 久久女婷五月综合色啪小说| 精品熟女少妇八av免费久了| 精品国产国语对白av| 久久国产精品男人的天堂亚洲| 国产亚洲欧美在线一区二区| 熟女少妇亚洲综合色aaa.| 色婷婷久久久亚洲欧美| 精品人妻熟女毛片av久久网站| 女人久久www免费人成看片| 80岁老熟妇乱子伦牲交| 老司机深夜福利视频在线观看 | 国产精品秋霞免费鲁丝片| 另类亚洲欧美激情| 在线亚洲精品国产二区图片欧美| 在线观看人妻少妇| 久久久精品免费免费高清| 国产又色又爽无遮挡免| 最近最新中文字幕大全免费视频| 久久免费观看电影| 成年人黄色毛片网站| 18禁裸乳无遮挡动漫免费视频| 老司机福利观看| 国产免费一区二区三区四区乱码| 精品少妇一区二区三区视频日本电影| 国产深夜福利视频在线观看| 国产又色又爽无遮挡免| 久久中文字幕一级| 伦理电影免费视频| 丝袜在线中文字幕| 久热这里只有精品99| 嫁个100分男人电影在线观看| 成人影院久久| 少妇 在线观看| 久久精品aⅴ一区二区三区四区| cao死你这个sao货| 美女福利国产在线| 涩涩av久久男人的天堂| 一级,二级,三级黄色视频| 黑人猛操日本美女一级片| 999久久久精品免费观看国产| 岛国毛片在线播放| 久久这里只有精品19| 91九色精品人成在线观看| 国产色视频综合| 韩国精品一区二区三区| 久久国产精品男人的天堂亚洲| 青春草视频在线免费观看| 80岁老熟妇乱子伦牲交| 亚洲,欧美精品.| 欧美精品人与动牲交sv欧美| 满18在线观看网站| 麻豆av在线久日| 夫妻午夜视频| 亚洲精品中文字幕在线视频| 国产精品一区二区在线观看99| 亚洲av男天堂| 丰满少妇做爰视频| a级毛片黄视频| 国产真人三级小视频在线观看| 中文欧美无线码| 国产成+人综合+亚洲专区| 精品国产一区二区三区四区第35| 中文欧美无线码| 亚洲中文日韩欧美视频| 久久久久久人人人人人| 欧美黄色片欧美黄色片| 久久这里只有精品19| av天堂久久9| 国产在线一区二区三区精| 国产男女超爽视频在线观看| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 国产一区二区激情短视频 | 亚洲一区中文字幕在线| 啦啦啦免费观看视频1| 最近最新免费中文字幕在线| 中文欧美无线码| 亚洲国产av新网站| 久久99热这里只频精品6学生| 黑人操中国人逼视频| 在线 av 中文字幕| 亚洲中文字幕日韩| 婷婷成人精品国产| 啦啦啦免费观看视频1| 亚洲第一青青草原| 女人精品久久久久毛片| 水蜜桃什么品种好| 精品人妻一区二区三区麻豆| 美女主播在线视频| 99九九在线精品视频| 99re6热这里在线精品视频| 99精品久久久久人妻精品| 欧美日韩精品网址| 免费高清在线观看日韩| 成人三级做爰电影| 国产精品麻豆人妻色哟哟久久| 满18在线观看网站| 亚洲欧美一区二区三区久久| 久久亚洲精品不卡| 成人国语在线视频| 国产精品一二三区在线看| 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 一区二区三区乱码不卡18| 亚洲男人天堂网一区| 午夜精品久久久久久毛片777| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美少妇被猛烈插入视频| 欧美成人午夜精品| 国产伦理片在线播放av一区| 777米奇影视久久| 午夜激情av网站| 俄罗斯特黄特色一大片| 国产国语露脸激情在线看| 男人爽女人下面视频在线观看| 一本一本久久a久久精品综合妖精| 中文字幕人妻熟女乱码| 丰满饥渴人妻一区二区三| 99精国产麻豆久久婷婷| 真人做人爱边吃奶动态| h视频一区二区三区| 性少妇av在线| 亚洲五月色婷婷综合| 免费观看av网站的网址| 久久综合国产亚洲精品| 免费在线观看完整版高清| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 国产伦人伦偷精品视频| 亚洲欧美清纯卡通| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| 亚洲全国av大片| 国产成人影院久久av| 亚洲成人手机| 19禁男女啪啪无遮挡网站| 亚洲成国产人片在线观看| 999久久久精品免费观看国产| 青春草视频在线免费观看| 亚洲人成77777在线视频| 婷婷丁香在线五月| 国产av国产精品国产| 啦啦啦啦在线视频资源| videosex国产| videosex国产| 如日韩欧美国产精品一区二区三区| 黄片大片在线免费观看| 国产成人av激情在线播放| 国产无遮挡羞羞视频在线观看| 19禁男女啪啪无遮挡网站| 久久国产精品人妻蜜桃| avwww免费| av又黄又爽大尺度在线免费看| 一级a爱视频在线免费观看| 高清黄色对白视频在线免费看| 精品国产乱子伦一区二区三区 | 性少妇av在线| 韩国精品一区二区三区| 久久免费观看电影| 又紧又爽又黄一区二区| 黄色视频不卡| 国产免费视频播放在线视频| 久久久国产欧美日韩av| 狂野欧美激情性xxxx| 国产无遮挡羞羞视频在线观看| 汤姆久久久久久久影院中文字幕| 国产免费一区二区三区四区乱码| 午夜福利在线免费观看网站| 国产日韩一区二区三区精品不卡| 少妇被粗大的猛进出69影院| 亚洲五月婷婷丁香| 99久久99久久久精品蜜桃| 老熟妇仑乱视频hdxx| 深夜精品福利| 色婷婷av一区二区三区视频| 精品久久久精品久久久| 飞空精品影院首页| 最近最新中文字幕大全免费视频| 国产亚洲欧美精品永久| 99国产精品一区二区三区| 欧美黑人精品巨大| 午夜福利在线观看吧| 桃花免费在线播放| 色婷婷久久久亚洲欧美| h视频一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 大片免费播放器 马上看| 亚洲精品久久成人aⅴ小说| 国产成人精品无人区| 一本—道久久a久久精品蜜桃钙片| 777久久人妻少妇嫩草av网站| 国产欧美日韩一区二区三 | 欧美久久黑人一区二区| 啦啦啦啦在线视频资源| 国产免费福利视频在线观看| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av| 69av精品久久久久久 | 欧美一级毛片孕妇| 亚洲伊人色综图| 成在线人永久免费视频| 久久午夜综合久久蜜桃| 欧美黑人欧美精品刺激| 青春草亚洲视频在线观看| 久久亚洲精品不卡| 中国国产av一级| 天天影视国产精品| av网站在线播放免费| 飞空精品影院首页| 日本一区二区免费在线视频| 色老头精品视频在线观看| 美女高潮到喷水免费观看| 精品高清国产在线一区| 18禁国产床啪视频网站| 成年女人毛片免费观看观看9 | 久久久久国内视频| 欧美中文综合在线视频| 天天添夜夜摸| 精品久久久久久久毛片微露脸 | 蜜桃在线观看..| 久久久精品免费免费高清| 久久人妻熟女aⅴ| 亚洲精品乱久久久久久| 高清欧美精品videossex| 乱人伦中国视频| 美女大奶头黄色视频| 亚洲精品av麻豆狂野| 99久久人妻综合| 99精国产麻豆久久婷婷| 久久人人爽av亚洲精品天堂| 国产欧美日韩一区二区三区在线| 久久狼人影院| 狠狠狠狠99中文字幕| 黑人猛操日本美女一级片| 国产精品久久久久久精品电影小说| 日韩熟女老妇一区二区性免费视频| 亚洲av男天堂| 欧美人与性动交α欧美软件| 亚洲成人国产一区在线观看| 国产精品成人在线| 久久久久久久精品精品| 三上悠亚av全集在线观看| 少妇裸体淫交视频免费看高清 | 99国产综合亚洲精品| 别揉我奶头~嗯~啊~动态视频 | 久久精品久久久久久噜噜老黄| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 黄片小视频在线播放| 精品国产乱码久久久久久小说| 看免费av毛片| 国产不卡av网站在线观看| 欧美日韩视频精品一区| 国产亚洲精品一区二区www | 欧美激情久久久久久爽电影 | 天堂中文最新版在线下载| 精品第一国产精品| 精品人妻熟女毛片av久久网站| 国产日韩一区二区三区精品不卡| 精品国产超薄肉色丝袜足j| 美女主播在线视频| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 午夜福利视频精品| 99精品欧美一区二区三区四区| 日日夜夜操网爽| 日韩欧美一区二区三区在线观看 | 亚洲欧美精品综合一区二区三区| 丝袜美腿诱惑在线| 精品国内亚洲2022精品成人 | 少妇人妻久久综合中文| 蜜桃国产av成人99| 国产一区二区在线观看av| 欧美日韩黄片免| 国产亚洲欧美在线一区二区| 人妻久久中文字幕网| 丝袜脚勾引网站| 免费日韩欧美在线观看| 97人妻天天添夜夜摸| 欧美97在线视频| 男女边摸边吃奶| 国产成人av激情在线播放| av线在线观看网站| 99国产精品一区二区蜜桃av | 成人免费观看视频高清| 免费看十八禁软件| av网站免费在线观看视频| 大码成人一级视频| 欧美日韩精品网址| 少妇被粗大的猛进出69影院| 精品国产乱码久久久久久小说| 亚洲天堂av无毛| 真人做人爱边吃奶动态| 国产淫语在线视频| 青草久久国产| 免费在线观看日本一区| 久久久久久人人人人人| 亚洲精品日韩在线中文字幕| www.精华液| 久久精品亚洲熟妇少妇任你| 不卡一级毛片| 欧美日韩亚洲综合一区二区三区_| 日本vs欧美在线观看视频| 精品高清国产在线一区| 啦啦啦中文免费视频观看日本| 久久国产精品男人的天堂亚洲| 免费在线观看黄色视频的| 91九色精品人成在线观看| www.999成人在线观看| 久久精品久久久久久噜噜老黄| 男男h啪啪无遮挡| 国产成人欧美| 亚洲第一青青草原| 亚洲精品久久久久久婷婷小说| av天堂久久9| 丰满饥渴人妻一区二区三| 久久精品aⅴ一区二区三区四区| 高清黄色对白视频在线免费看| 中文字幕高清在线视频| 三上悠亚av全集在线观看| 国产xxxxx性猛交| 母亲3免费完整高清在线观看| 国产亚洲午夜精品一区二区久久| av电影中文网址| 丁香六月天网| 亚洲国产欧美一区二区综合| 丰满迷人的少妇在线观看| 国产精品影院久久| 日韩,欧美,国产一区二区三区| 欧美日韩福利视频一区二区| 国产成+人综合+亚洲专区| 欧美黑人精品巨大| 麻豆国产av国片精品| 精品国内亚洲2022精品成人 | 亚洲av电影在线进入| 精品久久蜜臀av无| 夜夜骑夜夜射夜夜干| 精品久久久久久久毛片微露脸 | 精品国产一区二区三区四区第35| 满18在线观看网站| 久久久久久久精品精品| 精品亚洲乱码少妇综合久久| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人免费观看mmmm| 一本大道久久a久久精品| 国产野战对白在线观看| 91国产中文字幕| 国产亚洲精品一区二区www | 欧美国产精品va在线观看不卡| 国产淫语在线视频| 国产在线免费精品| 免费在线观看日本一区| 国产xxxxx性猛交| 香蕉国产在线看| 亚洲自偷自拍图片 自拍| 老汉色av国产亚洲站长工具| 亚洲国产精品一区三区| 天天躁夜夜躁狠狠躁躁| 亚洲精华国产精华精| 女人爽到高潮嗷嗷叫在线视频| 男女床上黄色一级片免费看| 日韩制服骚丝袜av| 国产精品一区二区精品视频观看| 19禁男女啪啪无遮挡网站| 两人在一起打扑克的视频| 亚洲专区国产一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 精品一区二区三卡| 色精品久久人妻99蜜桃| 天堂8中文在线网| 精品少妇一区二区三区视频日本电影| 一个人免费在线观看的高清视频 | 老鸭窝网址在线观看| 色视频在线一区二区三区| 嫩草影视91久久| 亚洲欧美一区二区三区黑人| 免费av中文字幕在线| 深夜精品福利| 精品熟女少妇八av免费久了| 又黄又粗又硬又大视频| 国产精品.久久久| av网站在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品自产自拍| 精品国内亚洲2022精品成人 | 18禁观看日本| 免费高清在线观看日韩| 成人18禁高潮啪啪吃奶动态图| 欧美精品高潮呻吟av久久| 亚洲精品日韩在线中文字幕| 国产99久久九九免费精品| 91成人精品电影| 亚洲美女黄色视频免费看| 五月天丁香电影| 69av精品久久久久久 | 国产真人三级小视频在线观看| 肉色欧美久久久久久久蜜桃| 欧美黄色淫秽网站| 亚洲五月色婷婷综合| 欧美老熟妇乱子伦牲交| 黄片大片在线免费观看| www.熟女人妻精品国产| 欧美精品亚洲一区二区| 999精品在线视频| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 国精品久久久久久国模美| 水蜜桃什么品种好| 青春草视频在线免费观看| 亚洲av日韩在线播放| 2018国产大陆天天弄谢| 美女国产高潮福利片在线看| 啦啦啦免费观看视频1| 成年美女黄网站色视频大全免费| 操美女的视频在线观看| 国产老妇伦熟女老妇高清| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 99国产综合亚洲精品| 免费在线观看日本一区| 欧美黑人欧美精品刺激| 日日夜夜操网爽| 99精国产麻豆久久婷婷| 免费日韩欧美在线观看| 19禁男女啪啪无遮挡网站| 男男h啪啪无遮挡| 99精国产麻豆久久婷婷| 久久国产精品影院| 日韩有码中文字幕| 国产精品香港三级国产av潘金莲|