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

    基于不同評價因子組合的鐵路沿線滑坡危險性評價

    2022-11-01 11:21:22郭榮昌余嶺燕
    蘭州交通大學(xué)學(xué)報 2022年5期
    關(guān)鍵詞:評價模型

    賀 攀,郭榮昌,張 蕊,余嶺燕

    (蘭州交通大學(xué)自動化與電氣工程學(xué)院,蘭州 730070)

    滑坡是土壤、泥漿、巖屑或巖石的坡面運動,是陸地斜坡環(huán)境中普遍存在的一種地質(zhì)災(zāi)害[1].鐵路沿線常遭受滑坡危害,滑坡對乘客生命財產(chǎn)和鐵路線路造成了嚴(yán)重的威脅,據(jù)初步統(tǒng)計,成都至白玉段鐵路沿線發(fā)育有滑坡126處[2].滑坡危險性(發(fā)生的概率大?。┰u價可為滑坡防治提供技術(shù)支撐,滑坡評價因子合理選取是危險性評價的重要部分.

    已有研究常依據(jù)專家經(jīng)驗以及分析滑坡的觸發(fā)機制、孕災(zāi)環(huán)境來進行滑坡評價因子選取,并采用不同的算法計算評價因子對于滑坡的貢獻率,進行因子篩選.針對滑坡危險性評價因子選取和篩選,文獻[3]通過計算因子間協(xié)方差和相關(guān)矩陣,選取滿足相互獨立的評價因子用于滑坡易發(fā)性評價;文獻[4]從16個評價因子中選取30組評價因子組合,對比不同評價因子組合對應(yīng)的模型精度,得到了最佳評價因子組合,也解決了因子共線性問題;文獻[5]選用主成分分析和特征遞歸消除,從92個特征中分別篩選出15個和30個特征用于滑坡災(zāi)害空間預(yù)測,實驗表明特征遞歸消除法篩選得到的評價因子組合所對應(yīng)的模型預(yù)測精度更高;文獻[6]使用因子相關(guān)分析、主成分分析和粗糙集對滑坡評價因子進行篩選和約簡,經(jīng)實驗驗證,篩選和約簡后評價因子可提高評價結(jié)果的精度和準(zhǔn)確性;文獻[7]提出Apriori算法對滑坡評價因子進行強關(guān)聯(lián)分析,用于因子篩選,實驗結(jié)果表明采用Apriori算法從預(yù)選因子中選取導(dǎo)致滑坡發(fā)生可能性較大的因子進行組合,所得到評價結(jié)果與實際滑坡分布更加吻合;文獻[8]采用遺傳算法優(yōu)化后的粗糙集對滑坡評價因子進行約簡,表明選用約簡后的評價因子所對應(yīng)的評價模型精度更高.以上研究,多采用單個因子篩選法生成評價因子組合,較少進行因子篩選方法對比,因此,有必要采用不同因子篩選法生成多個評價因子組合,對比評價因子組合的合理性和準(zhǔn)確性,選取最適合研究區(qū)域的評價因子組合,用于滑坡危險性評價.

    綜合已有研究,按照以下三點進行評價因子選取和篩選:1)綜合考慮專家經(jīng)驗與實地滑坡調(diào)查,預(yù)選滑坡評價因子;2)選取的評價因子間應(yīng)相對獨立,降低評價結(jié)果過擬合現(xiàn)象[9];3)在評價因子篩選的過程中,要考慮各評價因子的權(quán)重大小以及相關(guān)性,應(yīng)選取滑坡發(fā)育的主要影響因素作為評價因子,舍去次要和冗余因子,以達到簡化評價模型的目的.因此,首先,通過遙感解譯影像、收集的滑坡資料以及野外調(diào)查,預(yù)選滑坡評價因子,建立滑坡災(zāi)害空間數(shù)據(jù)庫,并對預(yù)選因子進行共線性檢驗,確保因子間相互獨立;然后,選用主成分分析法、粗糙集和灰色關(guān)聯(lián)3種因子篩選方法,生成不同的評價因子組合;最后,將不同的評價因子組合輸入到滑坡危險性評價性能較好的隨機森林[10]、支持向量機[11]、邏輯回歸[12]3種評價模型,通過對比模型精度,選擇最優(yōu)的評價因子組合和模型進行危險性評價,使之達到最佳的鐵路沿線滑坡危險性評價效果,為鐵路的規(guī)劃建設(shè)和后期列車安全運營、防災(zāi)減害提供參考依據(jù).

    1 研究區(qū)域與數(shù)據(jù)來源

    1.1 研究區(qū)域

    雅安至巴塘段鐵路位于四川省西部,全長約520 km,途徑雅安、天全、瀘定、康定、雅江、理塘和巴塘等站.區(qū)域內(nèi)的鐵路沿線區(qū)域氣候復(fù)雜多變,斷裂構(gòu)造發(fā)育,地質(zhì)災(zāi)害點多面廣、險情重、危害大.該區(qū)域內(nèi)滑坡多為降雨型滑坡,年降水量在500~1 600 mm之間,降水主要集中在5~9月.

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

    通過對研究區(qū)域滑坡點的勘察和四川省自然資源廳對研究區(qū)域滑坡災(zāi)害的信息公布,分析滑坡的成災(zāi)原因和孕災(zāi)環(huán)境,并結(jié)合前人的研究[13-15],選取剖面曲率D1、平面曲率D2、坡向D3、高程D4、坡度D5、地層巖性 D6、距斷層距離 D7、距河流距離 D8、植被覆蓋率D9、降雨量D10、土地利用類型D11和人類工程活動D12共12個因子作為預(yù)選評價因子,這些因子可以綜合呈現(xiàn)研究區(qū)域的地質(zhì)條件、地貌條件以及滑坡災(zāi)害的觸發(fā)條件等情況.其中:選取當(dāng)次降雨過程中最大24 h降雨量表示D10;用距道路距離表示D12.各評價因子數(shù)據(jù)來源見表1.

    表1 評價因子數(shù)據(jù)來源Tab.1 Data sources of evaluation factors

    該區(qū)域歷史滑坡數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心,包括滑坡位置信息、災(zāi)害險情等級、滑坡傷亡人數(shù)、財產(chǎn)損失、滑坡構(gòu)造和防治建議等信息,鐵路沿線滑坡分布情況如圖1所示.提取研究區(qū)域內(nèi)160處歷史滑坡災(zāi)害和160處非滑坡建立滑坡災(zāi)害空間數(shù)據(jù)庫,其中:非滑坡為歷史滑坡點附近500 m以內(nèi)隨機非滑坡點.

    圖1 研究區(qū)域鐵路沿線滑坡分布圖Fig.1 Distribution map of landslides along the railway in the study area

    2 基于不同評價因子組合的鐵路沿線滑坡危險性評價模型

    基于不同評價因子組合的鐵路沿線滑坡危險性評價流程如圖2所示,主要分為4個步驟:

    圖2 基于不同評價因子組合的鐵路沿線滑坡危險性評價流程圖Fig.2 Flow chart of landslide risk assessment along railway based on different combination of evaluation factors

    Step1:對預(yù)選的12個滑坡評價因子進行共線性檢驗;

    Step2:采用主成分分析法(principal component analysis,PCA)、粗糙集(rough set,RS)和灰色關(guān)聯(lián)分析法(grey relation analysis,GRA)3種因子篩選法生成不同評價因子組合;

    Step3:將評價因子組合輸入到隨機森林、支持向量機、邏輯回歸3種評價模型,并計算模型精度;

    Step4:對比模型精度大小,選擇最佳的評價因子組合和評價模型用于滑坡危險性評價.

    2.1 因子共線性檢驗

    預(yù)選的滑坡評價因子間可能不是完全獨立的,存在因子共線性問題,導(dǎo)致危險性評價結(jié)果的準(zhǔn)確性降低,因此,在因子篩選和構(gòu)建滑坡危險性評價模型之前,需要對預(yù)選因子進行多重共線性檢驗.當(dāng)因子存在多重共線性時,因子重要性就會被影響,導(dǎo)致無法正確解釋和理解特征;當(dāng)因子不存在多重共線性時,方可用于因子篩選和模型訓(xùn)練.通過計算容忍度(T)和方差膨脹因子(VIF)的方式進行因子多重共線性檢驗,計算公式為

    其中:A2為評價因子間的方差.

    2.2 因子篩選

    主成分分析法通過對數(shù)據(jù)進行降維,構(gòu)造原始因素之間的線性組合,減少數(shù)據(jù)的維度,降維后的數(shù)據(jù)產(chǎn)生的主成分之間互不相關(guān),且主成分能夠反映原始因素所提供的大部分信息[16].計算各評價因子的特征值和貢獻率,利用各主成分累計貢獻率(一般達到85%以上)確定主成分的個數(shù)[17],從而確定選取和剔除的因子.

    粗糙集是一種能夠處理不完整、不確定信息的歸納、分類方法[18],其主要用于約簡初始的評價指標(biāo),使得在保持屬性分類性能基本不變的情況下,剔除冗余因子,從而保留核因子.

    灰色關(guān)聯(lián)分析法是根據(jù)指標(biāo)之間的相似程度來衡量樣本與參考樣本之間關(guān)聯(lián)程度的一種方法[19].首先,選取評價因子的最優(yōu)值或最劣值組成參考樣本x0;接著,采用均值化處理評價因子,得到n組樣本m個指標(biāo)的樣本矩陣;然后,計算比較樣本xi與參考樣本的絕對差值 |x0(k)-xi(k)|(k=1,2,…,m;i=1,2,…,n)以及絕對差值的最大值與最小值;最后,通過式(2)~(3)計算關(guān)聯(lián)系數(shù)ξi和關(guān)聯(lián)度ri.關(guān)聯(lián)度為各因子對于滑坡事件發(fā)生的貢獻率,按照因子的關(guān)聯(lián)度由大到小進行累加,選取滿足累計貢獻率要求的因子,從而剔除其余因子.參考主成分分析法的累計貢獻率,設(shè)置灰色關(guān)聯(lián)分析法累計貢獻率為85%.

    其中:ρ為分辨系數(shù),一般取值為 0.5[19].

    預(yù)選評價因子中可能存在貢獻率較弱的因子,這些因子在造成復(fù)雜計算程度的同時,還將影響計算的效率以及準(zhǔn)確度.使用PCA,RS,GRA對預(yù)選因子進行篩選或約簡,根據(jù)各自的量化方法,選取貢獻率較大的滑坡評價因子或者核因子,生成3種評價因子組合.

    2.3 評價模型

    隨機森林(random forest,RF)是機器學(xué)習(xí)算法中的一種組合分類器,從N個樣本中有放回抽取n個樣本(n<N),從K個總屬性中選取k個屬性(k<K),基于基尼指數(shù)選擇最佳分割屬性創(chuàng)建決策樹,通過Bagging算法可集成多個決策樹,從而組成隨機森林[20].隨機森林訓(xùn)練速度快,不易“過擬合”,對噪聲和異常值的容忍度好[21].

    支持向量機(support vector machine,SVM)在處理二分類問題中尋找最優(yōu)的超平面,以最大間隔區(qū)分它們[22].對于線性可分數(shù)據(jù),尋找最大間距來完成分類;處理非線性可分數(shù)據(jù)時,通過核函數(shù)將其映射到高維空間,完成線性可分,引入松弛變量來控制分類誤差.

    邏輯回歸(logistic regression,LR)原理是將線性回歸的結(jié)果映射到0~1之間,適用于因變量為二分類的問題,且對數(shù)據(jù)要求較為寬松[23].用于滑坡危險性評價時,自變量為滑坡評價因子,因變量為滑坡是否發(fā)生,其中:0表示滑坡未發(fā)生;1表示滑坡發(fā)生.通過建立最優(yōu)的擬合函數(shù)揭示自變量和因變量之間關(guān)系.

    對隨機森林、支持向量機進行參數(shù)調(diào)整,用于優(yōu)化模型精度;然后,將篩選后的評價因子組合分別輸入到隨機森林、支持向量機、邏輯回歸3種評價模型,并計算3種評價因子組合分別輸入到3種評價模型的預(yù)測精度.

    2.4 滑坡危險性評價

    選用預(yù)測精度最佳的評價因子組合和模型生成研究區(qū)域滑坡危險性分布圖,根據(jù)危險性大小劃分風(fēng)險等級.統(tǒng)計分析鐵路線路危險區(qū)域分布情況,并對滑坡隱患點進行驗證,完成鐵路沿線滑坡危險性評價.

    3 實驗分析

    3.1 因子共線性檢驗

    通過計算評價因子容忍度和方差膨脹因子的方式進行因子多重共線性檢驗,檢驗結(jié)果見表2.當(dāng)容忍度小于0.1或方差膨脹因子大于10,認為因子存在共線性[24],表2中各因子容忍度和方差膨脹因子不符合共線性條件,因此,認為所選因子之間相關(guān)性弱或不相關(guān),預(yù)選的12個滑坡評價因子滿足相互獨立的要求,可直接用于因子篩選和參與模型訓(xùn)練,并將此12個評價因子作為評價因子組合一.

    表2 評價因子共線性檢驗表Tab.2 Evaluation factor collinearity test table

    3.2 評價因子篩選

    預(yù)選的12個滑坡評價因子中可能包括相關(guān)性較差的評價因子,因此,選用SPSS25軟件對評價因子進行主成分分析,計算各評價因子特征值和累計貢獻率,計算結(jié)果見表3.當(dāng)累計貢獻率達到86.94%時,剔除未參與計算累計貢獻率的評價因子,即剔除剖面曲率、平面曲率、到斷層距離和土地利用類型,保留降雨、坡度、到河流距離、植被覆蓋率、人類工程活動、地層巖性、坡向和高程,并將此8個評價因子作為評價因子組合二.

    表3 評價因子的特征值和累計貢獻率Tab.3 Eigenvalue and cum ulative contribution rate of evaluation factors

    用粗糙集進行因子篩選,將預(yù)選的12個評價因子作為條件屬性,選擇滑坡發(fā)生與不發(fā)生作為決策屬性,構(gòu)建初始決策表,對初始決策表進行屬性約簡,找出最小的條件屬性集和滑坡核屬性.使用Rosetta軟件對條件屬性和相應(yīng)的決策屬性進行屬性約簡,最終得到8個核因子,分別是降雨、地層巖性、到斷層距離、坡度、坡向、到河流距離、植被覆蓋率和人類工程活動,將此8個核因子作為評價因子組合三,剔除剖面曲率、平面曲率、土地利用類型和高程.

    對預(yù)選的12個評價因子進行灰色關(guān)聯(lián)分析,確定參考樣本后生成樣本矩陣,通過計算比較樣本和參考樣本的差值,求得關(guān)聯(lián)系數(shù)以及關(guān)聯(lián)度.各因子的關(guān)聯(lián)度為因子對于滑坡事件發(fā)生的貢獻率大小,關(guān)聯(lián)度和累計貢獻率計算結(jié)果見表4.當(dāng)累計貢獻率達到85.1%時,剔除未參與計算累計貢獻率的評價因子,即剔除到河流距離、坡向、平面曲率和剖面曲率,保留降雨、坡度、地層巖性、人類工程活動、植被覆蓋率、土地利用類型、到斷層距離和高程,將此8個評價因子作為評價因子組合四.

    表4 評價因子的灰色關(guān)聯(lián)度和累計貢獻率Tab.4 Grey relation degree and cumulative contribution rate of evaluation factors

    3.3 評價模型優(yōu)化與精度對比

    將320組數(shù)據(jù)按照7:3劃分訓(xùn)練集和測試集,即224組數(shù)據(jù)用于訓(xùn)練,96組數(shù)據(jù)用于測試,將滑坡數(shù)據(jù)標(biāo)記為1,非滑坡數(shù)據(jù)標(biāo)記為0.為了使評價模型滿足精度要求,對隨機森林模型各參數(shù)進行尋優(yōu),再利用網(wǎng)格搜索和交叉驗證得到各參數(shù)最優(yōu)值為:n-estimators=69,max-features=3,min-samplesleaf=3,max-depth=12;對支持向量機模型進行參數(shù)調(diào)整,選擇核函數(shù)為徑向基核函數(shù),通過交叉驗證得到核函數(shù)中的參數(shù)γ為0.8,以及支持向量機模型中懲罰系數(shù)c為10.

    受試者工作特征曲線(receiver operating characteristic curve,ROC曲線)已廣泛應(yīng)用于地質(zhì)災(zāi)害危險性評價結(jié)果精度分析中[25].ROC曲線使用曲線下面積(area under curve,AUC)值衡量模型的預(yù)測精度,AUC值為0至1,越接近于1,說明模型精度越好.將評價因子組合一、二、三、四分別輸入到RF,SVM和LR 3種評價模型,生成ROC曲線并計算模型AUC值,最后共得到12種不同組合的ROC曲線,實驗結(jié)果如圖3所示.圖 3(a)、(b)、(c)分別為4種評價因子組合輸入到RF,SVM和LR模型中的ROC曲線,12種組合的模型AUC值見表5.由表5可知:評價因子組合四輸入隨機森林模型中所得到的預(yù)測精度最佳,且AUC值為0.886 1.對比12種組合模型的AUC值可知:PCA,RS和GRA 3種篩選因子的方法可提高評價模型的預(yù)測精度,且AUC值平均提高0.034 6,0.044 8和0.054 5.對比 RF,SVM和LR 3種評價模型AUC值可知,隨機森林整體預(yù)測精度較高.

    表5 12種組合模型的AUC值Tab.5 AUC values of 12 combined models

    圖3 12種組合模型的ROC曲線Fig.3 ROC curves of 12 combined models

    綜上,在雅安至巴塘段鐵路沿線滑坡危險性評價中,評價因子組合四輸入隨機森林模型中所得到的預(yù)測精度最佳,因此,選用評價因子組合四和隨機森林模型開展研究區(qū)滑坡危險性評價.

    3.4 滑坡危險性評價

    選用灰色關(guān)聯(lián)分析法篩選得到的評價因子組合四輸入到隨機森林模型,進行滑坡危險性評價,計算滑坡危險性大小,并生成滑坡危險性分布圖,如圖4所示.基于自然斷點法將滑坡危險性評價結(jié)果劃分為低、較低、中、較高和高5個風(fēng)險等級,普遍認為較高和高風(fēng)險地區(qū)為易受滑坡威脅的區(qū)域.選取103個滑坡隱患點用于檢驗滑坡危險性評價結(jié)果,可得86.4%的滑坡隱患點位于較高和高風(fēng)險地區(qū),即認為驗證精度為86.4%.對圖4滑坡危險性分布圖進行統(tǒng)計分析,得到各風(fēng)險度級別中鐵路長度及比例(見表6),其中:位于高和較高風(fēng)險地區(qū)的鐵路線路長度分別為10.9 km(2.1%)和48.7 km(9.4%),且較多分布在雅安地區(qū).評價結(jié)果與現(xiàn)場調(diào)查結(jié)果基本一致,經(jīng)過現(xiàn)場調(diào)查分析,雅安地區(qū)鐵路沿線滑坡危險性較高的主要原因為雅安降雨量較大,造成坡體的抗剪強度降低,加速坡體解體破壞,其中:雅安市近60年平均年降雨量為1 664.3 mm,為研究區(qū)域內(nèi)最高,而巴塘縣近60年平均年降雨量為497.4 mm,為研究區(qū)域最低,同時雅安人類工程活動較強烈,道路和居民地分布較為密集,為滑坡發(fā)生提供了良好的孕災(zāi)環(huán)境.在線路建設(shè)時,使用科學(xué)手段監(jiān)測危險性較大的坡體,提前做好防御滑坡的措施,以防止滑坡對鐵路線路造成巨大的損害;在后期鐵路運營時,將滑坡災(zāi)害危險性評價結(jié)果與鐵路信號預(yù)警系統(tǒng)緊密結(jié)合,保證列車安全運行,為鐵路運營提供安全保障.

    表6 各風(fēng)險級別中鐵路長度及比例Tab.6 Length and proportion of railway in each risk level

    圖4 滑坡危險性分布圖Fig.4 Landslide risk distribution map

    4 結(jié)論

    以雅安至巴塘段鐵路為研究區(qū)域,基于PCA,RS和GRA 3種因子篩選方法與RF,SVM和LR 3種評價模型組合的方法,選用預(yù)測精度最優(yōu)的方法進行滑坡危險性評價,得到以下結(jié)論:

    1)3種因子篩選方法使用不同的量化方法對評價因子的相關(guān)性進行分析,剔除相關(guān)性較差的因子,保留主要因子.篩選結(jié)果表明:降雨、坡度和地層巖性為導(dǎo)致滑坡事件發(fā)生的重要因子,平面曲率和剖面曲率與滑坡發(fā)生相關(guān)性較弱,分析結(jié)果與實際滑坡災(zāi)害經(jīng)驗分析相符.

    2)12種組合模型中,灰色關(guān)聯(lián)分析法篩選出的評價因子組合輸入隨機森林模型,其對應(yīng)的模型預(yù)測精度最佳,AUC值為0.886 1.通過灰色關(guān)聯(lián)分析法進行因子篩選提高了模型預(yù)測精度,AUC面積增大了0.071 9.

    3)基于不同評價因子組合的鐵路沿線滑坡危險性評價,選取最優(yōu)的評價因子組合和評價模型,由其得到的滑坡危險性評價結(jié)果與滑坡災(zāi)害隱患相比,吻合程度達86.4%,證明此研究方法可靠、準(zhǔn)確,具有一定的實用價值,可為雅安至巴塘段鐵路的中期建設(shè)和后期運營中的風(fēng)險管理提供一定的參考和技術(shù)支撐.

    猜你喜歡
    評價模型
    一半模型
    SBR改性瀝青的穩(wěn)定性評價
    石油瀝青(2021年4期)2021-10-14 08:50:44
    中藥治療室性早搏系統(tǒng)評價再評價
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    基于Moodle的學(xué)習(xí)評價
    關(guān)于項目后評價中“專項”后評價的探討
    保加利亞轉(zhuǎn)軌20年評價
    成人综合一区亚洲| 国产亚洲精品久久久com| .国产精品久久| 女人被狂操c到高潮| 免费人成视频x8x8入口观看| 美女国产视频在线观看| 中文字幕av成人在线电影| 欧美激情国产日韩精品一区| 亚洲人成网站在线播| 国产亚洲av嫩草精品影院| 国产高清有码在线观看视频| 亚州av有码| 国产精品av视频在线免费观看| 国产乱人视频| 国产女主播在线喷水免费视频网站 | 亚洲第一区二区三区不卡| 久久这里有精品视频免费| 久久久久性生活片| 日韩欧美精品免费久久| 午夜亚洲福利在线播放| 欧美性猛交黑人性爽| 国产白丝娇喘喷水9色精品| 精品久久久久久久久久久久久| 美女 人体艺术 gogo| 美女内射精品一级片tv| 黄色视频,在线免费观看| 色哟哟·www| av卡一久久| 久久久久性生活片| 久久久久免费精品人妻一区二区| 亚洲真实伦在线观看| АⅤ资源中文在线天堂| 一区二区三区四区激情视频 | 91午夜精品亚洲一区二区三区| 亚州av有码| 男女视频在线观看网站免费| 亚洲精品乱码久久久v下载方式| 特级一级黄色大片| 黄色欧美视频在线观看| 51国产日韩欧美| 亚洲自偷自拍三级| 91久久精品电影网| 国产伦一二天堂av在线观看| 女的被弄到高潮叫床怎么办| 一进一出抽搐动态| 丝袜美腿在线中文| 国产一区二区亚洲精品在线观看| 美女黄网站色视频| 久久人人精品亚洲av| 91久久精品电影网| 99久久人妻综合| 两个人视频免费观看高清| АⅤ资源中文在线天堂| 久久久久久久久久黄片| 亚洲国产高清在线一区二区三| 伦精品一区二区三区| 欧美一级a爱片免费观看看| 亚洲美女搞黄在线观看| 成人av在线播放网站| 一边摸一边抽搐一进一小说| 亚洲高清免费不卡视频| 一个人观看的视频www高清免费观看| or卡值多少钱| 久久精品综合一区二区三区| 久久精品国产自在天天线| 一区二区三区高清视频在线| 岛国毛片在线播放| 成年版毛片免费区| 国语自产精品视频在线第100页| 久久久a久久爽久久v久久| 国产美女午夜福利| 欧美成人一区二区免费高清观看| 成人午夜高清在线视频| 日本av手机在线免费观看| 狂野欧美白嫩少妇大欣赏| 联通29元200g的流量卡| 久久精品国产鲁丝片午夜精品| 国产精品不卡视频一区二区| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 成年免费大片在线观看| 亚洲欧美精品综合久久99| 中文字幕制服av| 国产片特级美女逼逼视频| 精品免费久久久久久久清纯| 97在线视频观看| 美女xxoo啪啪120秒动态图| 久久人人爽人人片av| 小蜜桃在线观看免费完整版高清| 亚洲精品日韩av片在线观看| 美女内射精品一级片tv| av黄色大香蕉| h日本视频在线播放| 全区人妻精品视频| av福利片在线观看| 天堂网av新在线| 老女人水多毛片| 美女被艹到高潮喷水动态| 成人永久免费在线观看视频| av天堂中文字幕网| 国产亚洲91精品色在线| a级毛色黄片| 亚洲在久久综合| av国产免费在线观看| 亚洲无线在线观看| 插阴视频在线观看视频| 欧美色欧美亚洲另类二区| 99久国产av精品国产电影| 成人性生交大片免费视频hd| 乱系列少妇在线播放| 晚上一个人看的免费电影| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| 国产精品女同一区二区软件| 欧美成人免费av一区二区三区| 国产精品.久久久| 国产私拍福利视频在线观看| 欧美xxxx性猛交bbbb| 91精品一卡2卡3卡4卡| 亚洲欧美日韩无卡精品| 性色avwww在线观看| 国产成人a区在线观看| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 麻豆乱淫一区二区| 婷婷亚洲欧美| 国产高清视频在线观看网站| 国产日本99.免费观看| 不卡视频在线观看欧美| 在线免费观看的www视频| 国产片特级美女逼逼视频| 国产色婷婷99| 欧美日本亚洲视频在线播放| 两个人视频免费观看高清| 日韩欧美国产在线观看| av天堂在线播放| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 此物有八面人人有两片| 六月丁香七月| 看黄色毛片网站| 婷婷色综合大香蕉| 国产亚洲av嫩草精品影院| 欧美变态另类bdsm刘玥| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 黄色配什么色好看| 只有这里有精品99| 午夜福利在线观看免费完整高清在 | 九色成人免费人妻av| 免费电影在线观看免费观看| 国产精品av视频在线免费观看| 激情 狠狠 欧美| 欧美激情在线99| 男人的好看免费观看在线视频| 欧美xxxx性猛交bbbb| 看片在线看免费视频| 久久99蜜桃精品久久| 久久6这里有精品| 亚洲天堂国产精品一区在线| 悠悠久久av| 99九九线精品视频在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 九九在线视频观看精品| 青春草亚洲视频在线观看| 成人亚洲精品av一区二区| 全区人妻精品视频| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 三级国产精品欧美在线观看| 久久久精品94久久精品| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 97在线视频观看| 中文字幕制服av| 国产一级毛片七仙女欲春2| 欧美性感艳星| 99久久成人亚洲精品观看| 中文字幕av成人在线电影| 亚洲欧美日韩卡通动漫| 国产色婷婷99| 男女做爰动态图高潮gif福利片| 亚洲国产色片| .国产精品久久| 欧美一区二区国产精品久久精品| 国产精品久久电影中文字幕| 久久精品人妻少妇| 国模一区二区三区四区视频| 看黄色毛片网站| 草草在线视频免费看| 国产淫片久久久久久久久| 亚洲性久久影院| 热99在线观看视频| 51国产日韩欧美| 一个人看的www免费观看视频| 国产精品,欧美在线| 亚洲最大成人手机在线| 看免费成人av毛片| 色综合站精品国产| 乱码一卡2卡4卡精品| 日本五十路高清| 免费观看人在逋| 久久草成人影院| 日韩欧美在线乱码| 国产免费一级a男人的天堂| 我的女老师完整版在线观看| 悠悠久久av| 高清日韩中文字幕在线| 国产一区二区在线av高清观看| 国产成人a∨麻豆精品| 国产精品一区www在线观看| 国产精品麻豆人妻色哟哟久久 | 97热精品久久久久久| 干丝袜人妻中文字幕| 久久久欧美国产精品| 在线观看av片永久免费下载| 国产av不卡久久| 国产精品久久久久久精品电影小说 | 日产精品乱码卡一卡2卡三| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 三级毛片av免费| 午夜精品一区二区三区免费看| 国产精品一区二区性色av| 国产精品一及| videossex国产| 亚洲欧美精品综合久久99| 日韩大尺度精品在线看网址| 午夜福利视频1000在线观看| 一进一出抽搐动态| 免费人成在线观看视频色| 伦理电影大哥的女人| 91aial.com中文字幕在线观看| 少妇猛男粗大的猛烈进出视频 | 三级经典国产精品| 九九久久精品国产亚洲av麻豆| 18禁裸乳无遮挡免费网站照片| 成人一区二区视频在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲精品国产成人久久av| 久久综合国产亚洲精品| 99热这里只有是精品在线观看| 国产精品不卡视频一区二区| 欧美又色又爽又黄视频| 哪里可以看免费的av片| 日本免费a在线| 欧美色欧美亚洲另类二区| 国产伦在线观看视频一区| 成人三级黄色视频| 91久久精品国产一区二区成人| 99热这里只有精品一区| 欧美最新免费一区二区三区| 亚洲av.av天堂| 午夜福利视频1000在线观看| 99精品在免费线老司机午夜| 欧美+亚洲+日韩+国产| 日日干狠狠操夜夜爽| 亚洲av电影不卡..在线观看| 国产成人一区二区在线| 久久99蜜桃精品久久| 中文字幕熟女人妻在线| 日韩一区二区三区影片| 在线国产一区二区在线| 2021天堂中文幕一二区在线观| 亚洲av第一区精品v没综合| 黄色欧美视频在线观看| 美女黄网站色视频| 又爽又黄a免费视频| 久久国产乱子免费精品| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 99国产精品一区二区蜜桃av| 日本熟妇午夜| 久久久精品欧美日韩精品| 国产 一区 欧美 日韩| 国产精品久久久久久亚洲av鲁大| 天美传媒精品一区二区| 精品久久久久久久久久久久久| 欧美日韩综合久久久久久| 在线免费观看的www视频| 精品人妻熟女av久视频| 不卡视频在线观看欧美| 亚洲av男天堂| 精品人妻熟女av久视频| 伦精品一区二区三区| 亚洲av男天堂| 精品一区二区三区视频在线| 91狼人影院| 99久久精品热视频| 国产片特级美女逼逼视频| 欧美人与善性xxx| 91麻豆精品激情在线观看国产| 国产精品一区www在线观看| 亚洲在久久综合| 美女黄网站色视频| 国产中年淑女户外野战色| 日韩,欧美,国产一区二区三区 | av免费观看日本| 成人无遮挡网站| 免费一级毛片在线播放高清视频| 村上凉子中文字幕在线| 欧美另类亚洲清纯唯美| 成人欧美大片| 国产真实伦视频高清在线观看| 高清毛片免费观看视频网站| a级一级毛片免费在线观看| www日本黄色视频网| 日本欧美国产在线视频| 亚洲成av人片在线播放无| 男女做爰动态图高潮gif福利片| 高清毛片免费观看视频网站| 一级毛片aaaaaa免费看小| 精品人妻偷拍中文字幕| 国产男人的电影天堂91| 非洲黑人性xxxx精品又粗又长| 欧美最黄视频在线播放免费| 男女啪啪激烈高潮av片| a级毛片免费高清观看在线播放| 人妻系列 视频| 亚洲色图av天堂| 丰满乱子伦码专区| 国产伦精品一区二区三区四那| av又黄又爽大尺度在线免费看 | 精品久久久噜噜| 91麻豆精品激情在线观看国产| 九九爱精品视频在线观看| 亚洲精品成人久久久久久| 国产亚洲欧美98| 2022亚洲国产成人精品| 热99re8久久精品国产| 亚洲第一区二区三区不卡| 一级毛片我不卡| 免费无遮挡裸体视频| а√天堂www在线а√下载| 成人永久免费在线观看视频| 日本av手机在线免费观看| 色5月婷婷丁香| 一级毛片电影观看 | 搞女人的毛片| 我要搜黄色片| 日本熟妇午夜| 一个人观看的视频www高清免费观看| 看免费成人av毛片| 久久韩国三级中文字幕| av天堂中文字幕网| 亚洲精品国产av成人精品| 女同久久另类99精品国产91| 禁无遮挡网站| 国产熟女欧美一区二区| 99久国产av精品国产电影| 精品熟女少妇av免费看| 日本一二三区视频观看| 亚洲第一区二区三区不卡| 禁无遮挡网站| 99久久久亚洲精品蜜臀av| 99在线视频只有这里精品首页| 永久网站在线| 国产片特级美女逼逼视频| 身体一侧抽搐| 国产黄色视频一区二区在线观看 | 国国产精品蜜臀av免费| 免费av观看视频| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 日本av手机在线免费观看| 一个人看视频在线观看www免费| 国产极品天堂在线| 韩国av在线不卡| 网址你懂的国产日韩在线| 国产av不卡久久| 国内久久婷婷六月综合欲色啪| 久久精品人妻少妇| 免费搜索国产男女视频| 啦啦啦韩国在线观看视频| 搡女人真爽免费视频火全软件| 亚洲欧美中文字幕日韩二区| 国产成人freesex在线| 国产精品免费一区二区三区在线| 韩国av在线不卡| 国产精品1区2区在线观看.| 成年女人永久免费观看视频| 亚洲av电影不卡..在线观看| 亚洲欧美精品自产自拍| 精品久久国产蜜桃| 中出人妻视频一区二区| 国内精品美女久久久久久| 成人无遮挡网站| 熟女电影av网| 国产淫片久久久久久久久| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 国产精品福利在线免费观看| 欧美极品一区二区三区四区| 91麻豆精品激情在线观看国产| 国产高清三级在线| 亚洲欧洲日产国产| 欧美日韩国产亚洲二区| 日韩av在线大香蕉| 熟女电影av网| 国产人妻一区二区三区在| 男女啪啪激烈高潮av片| 男人狂女人下面高潮的视频| 婷婷六月久久综合丁香| 乱码一卡2卡4卡精品| 亚洲自偷自拍三级| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| 日韩在线高清观看一区二区三区| 91在线精品国自产拍蜜月| 三级经典国产精品| 国产亚洲精品久久久com| 12—13女人毛片做爰片一| 能在线免费观看的黄片| 国产极品天堂在线| 中文字幕制服av| 久久99热这里只有精品18| 直男gayav资源| 日韩三级伦理在线观看| 看黄色毛片网站| 少妇裸体淫交视频免费看高清| 国内精品一区二区在线观看| 性色avwww在线观看| av在线蜜桃| 国产精品美女特级片免费视频播放器| 国产亚洲5aaaaa淫片| 久久草成人影院| 国产高清激情床上av| 观看免费一级毛片| 亚洲精品粉嫩美女一区| 国产精品不卡视频一区二区| 亚洲aⅴ乱码一区二区在线播放| 人人妻人人澡人人爽人人夜夜 | 日韩一本色道免费dvd| 色视频www国产| eeuss影院久久| av黄色大香蕉| 精品久久国产蜜桃| 精品少妇黑人巨大在线播放 | 在线播放国产精品三级| 国产成人精品婷婷| 亚洲欧美成人精品一区二区| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 能在线免费观看的黄片| 美女大奶头视频| 亚洲国产精品国产精品| 麻豆国产97在线/欧美| 免费av观看视频| 国产高清有码在线观看视频| 国产精品久久久久久久电影| 国产美女午夜福利| 日韩中字成人| 精品久久久久久久久久免费视频| 久久精品综合一区二区三区| 久久鲁丝午夜福利片| 一级毛片我不卡| 欧美一区二区精品小视频在线| 亚洲自拍偷在线| 一级av片app| 国产高清激情床上av| 亚洲在线自拍视频| 18禁在线播放成人免费| 国语自产精品视频在线第100页| 国产亚洲av嫩草精品影院| 亚洲精品乱码久久久久久按摩| 蜜桃亚洲精品一区二区三区| 久久久成人免费电影| 波多野结衣高清无吗| 亚洲人与动物交配视频| 91av网一区二区| 亚洲图色成人| 亚洲欧洲日产国产| 国内精品一区二区在线观看| 国产亚洲精品久久久com| 波多野结衣高清无吗| 欧美极品一区二区三区四区| 国产男人的电影天堂91| 欧美激情国产日韩精品一区| 插阴视频在线观看视频| 国产精品久久久久久亚洲av鲁大| 免费黄网站久久成人精品| 亚洲国产欧美在线一区| 欧美+日韩+精品| 91久久精品国产一区二区成人| 天堂中文最新版在线下载 | 国产午夜福利久久久久久| 精品无人区乱码1区二区| 少妇丰满av| 午夜精品一区二区三区免费看| 中国国产av一级| 黄色一级大片看看| 99在线视频只有这里精品首页| 中文欧美无线码| 日韩欧美国产在线观看| 国产精品爽爽va在线观看网站| 少妇熟女aⅴ在线视频| 日韩高清综合在线| 久久精品人妻少妇| 九九爱精品视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲无线观看免费| 一本久久精品| 国产三级在线视频| 国产又黄又爽又无遮挡在线| a级毛片a级免费在线| 国产毛片a区久久久久| 国产精品嫩草影院av在线观看| 亚洲乱码一区二区免费版| 天堂影院成人在线观看| 舔av片在线| 99热精品在线国产| 久久综合国产亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 级片在线观看| 亚洲最大成人av| 三级男女做爰猛烈吃奶摸视频| a级毛片免费高清观看在线播放| 成人高潮视频无遮挡免费网站| 男插女下体视频免费在线播放| 蜜臀久久99精品久久宅男| 国产av一区在线观看免费| 91麻豆精品激情在线观看国产| 久久久a久久爽久久v久久| 免费看a级黄色片| 黑人高潮一二区| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 免费黄网站久久成人精品| 日韩中字成人| 色尼玛亚洲综合影院| 丰满人妻一区二区三区视频av| 夜夜夜夜夜久久久久| 免费观看a级毛片全部| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 久久久久网色| 只有这里有精品99| 精品人妻熟女av久视频| 国产精品久久久久久亚洲av鲁大| 欧美成人精品欧美一级黄| 热99在线观看视频| 国产色爽女视频免费观看| 亚洲av免费在线观看| 亚洲第一电影网av| 菩萨蛮人人尽说江南好唐韦庄 | 偷拍熟女少妇极品色| 日韩一区二区视频免费看| 国产一区二区在线av高清观看| 99久久精品一区二区三区| 中国国产av一级| 天天一区二区日本电影三级| 国产精品久久久久久av不卡| 在线观看66精品国产| 男人和女人高潮做爰伦理| 日本免费一区二区三区高清不卡| 一级黄片播放器| 久久人人精品亚洲av| 亚洲av不卡在线观看| 亚洲综合色惰| 日韩大尺度精品在线看网址| 成人特级av手机在线观看| 国产免费男女视频| 免费搜索国产男女视频| 在线观看免费视频日本深夜| 好男人在线观看高清免费视频| 又粗又爽又猛毛片免费看| 国产日本99.免费观看| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av天美| 爱豆传媒免费全集在线观看| 免费看日本二区| 精华霜和精华液先用哪个| 国内久久婷婷六月综合欲色啪| 亚洲自偷自拍三级| 国产在线男女| 女人十人毛片免费观看3o分钟| 一区二区三区高清视频在线| 国产精品永久免费网站| 亚洲最大成人av| 一个人看的www免费观看视频| 免费观看的影片在线观看| 一区二区三区高清视频在线| 免费看日本二区| 美女黄网站色视频| 日韩欧美在线乱码| 丝袜美腿在线中文| 国产高清不卡午夜福利| 99在线视频只有这里精品首页| 日本五十路高清| 校园春色视频在线观看| 国产精品,欧美在线| 联通29元200g的流量卡| 久久综合国产亚洲精品| 一个人看的www免费观看视频| 日本欧美国产在线视频| 亚洲av中文av极速乱| 精品日产1卡2卡| 97超视频在线观看视频| 波多野结衣巨乳人妻| 嫩草影院新地址| 国产精品久久久久久久电影| 亚洲人与动物交配视频| 亚洲精品成人久久久久久| 国产 一区 欧美 日韩| a级毛片免费高清观看在线播放| 久久鲁丝午夜福利片| 禁无遮挡网站| 国产极品天堂在线| 精品久久久久久成人av| 黑人高潮一二区| 尾随美女入室| videossex国产| 久久精品人妻少妇| 中文字幕制服av| 亚洲国产欧洲综合997久久,|