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

    基于EEMD-SVM模型的邊坡變形預(yù)測(cè)

    2023-11-17 01:36:50李建新
    廣東水利水電 2023年10期
    關(guān)鍵詞:原始數(shù)據(jù)分量模態(tài)

    李建新,肖 鋼,唐 鈴

    (1.東莞市地理信息與規(guī)劃編制研究中心,廣東 東莞 523000;2.贛州市國(guó)土空間調(diào)查規(guī)劃研究中心,江西 贛州 341000)

    1 概述

    邊坡是自然或人工開(kāi)挖形成的,也是工程建設(shè)中最常見(jiàn)的工程形式,而邊坡失穩(wěn)造成的滑坡時(shí)有發(fā)生,造成嚴(yán)重的經(jīng)濟(jì)損失與人員傷亡[1]。實(shí)際上,滑坡災(zāi)害的發(fā)生并不是一蹴而就,而是邊坡變形達(dá)到一定程度的結(jié)果。如果及時(shí)建立了基于邊坡變形監(jiān)測(cè)數(shù)據(jù)的高精度預(yù)測(cè)模型,準(zhǔn)確識(shí)別邊坡變形趨勢(shì),一定程度上可以避免災(zāi)害的發(fā)生。因此,建立高精度的邊坡變形預(yù)測(cè)模型是必要的,具有重要的社會(huì)價(jià)值和經(jīng)濟(jì)價(jià)值。

    邊坡變形過(guò)程復(fù)雜,不僅受到地質(zhì)條件的控制,還受到其他因素的影響,包括地下水、降雨量和人工活動(dòng)等,難以通過(guò)傳統(tǒng)方法建立高精度的邊坡變形預(yù)測(cè)模型[2]。隨著近些年來(lái)機(jī)器學(xué)習(xí)的不斷發(fā)展,BP神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)(SVM)和極限學(xué)習(xí)機(jī)等方法廣泛應(yīng)用于邊坡變形預(yù)測(cè)[3-8]。SVM是Vapnik基于統(tǒng)計(jì)學(xué)習(xí)理論提出的一種機(jī)器學(xué)習(xí)方法[9],對(duì)于非線性、小樣本等問(wèn)題具有較好的預(yù)測(cè)精度。研究表明SVM模型應(yīng)用于邊坡或滑坡其他領(lǐng)域同樣具有較高的精度,具有很好的適應(yīng)性[10-11]。

    邊坡變形監(jiān)測(cè)過(guò)程中受到多種因素的干擾,導(dǎo)致監(jiān)測(cè)數(shù)據(jù)中存在噪聲,進(jìn)而影響預(yù)測(cè)精度。近些年來(lái)基于集成經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)的去噪方法逐漸應(yīng)用在各個(gè)領(lǐng)域,結(jié)果表明EEMD能夠作為一種有效的去噪方法[12-13]。EEMD 方法的最大優(yōu)點(diǎn)在于其能夠以自適應(yīng)方式提取信號(hào)的各個(gè)分量和變化趨勢(shì)。同時(shí),相較于經(jīng)驗(yàn)?zāi)B(tài)分解方法,EEMD 方法顯著減少了模態(tài)混疊現(xiàn)象,并且也克服了小波變換容易造成許多虛假諧波。簡(jiǎn)而言之,EEMD 方法是一種高效的信號(hào)處理方法,能夠更加準(zhǔn)確地提取出信號(hào)中的各種信息。因此通過(guò)EEMD方法去除邊坡變形數(shù)據(jù)中的噪聲,使得SVM模型能夠更精確的識(shí)別邊坡的變形趨勢(shì),提高邊坡變形的預(yù)測(cè)精度。

    本文提出一種由EEMD去噪和SVM組合而成的預(yù)測(cè)模型,以提高邊坡變形預(yù)測(cè)的精度。模型首先通過(guò)集成經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)算法將邊坡變形時(shí)間序列進(jìn)行分解,基于相關(guān)系數(shù)的計(jì)算去除噪聲后重構(gòu)邊坡變形時(shí)間序列,然后將重構(gòu)后的邊坡變形時(shí)間序列作為SVM模型的輸入,從而建立EEMD-SVM模型,既而預(yù)測(cè)邊坡變形值。

    2 模型理論

    2.1 集成經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)

    經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)是一種能夠?qū)⒎瞧椒€(wěn)和非線性時(shí)間序列自適應(yīng)的分解為IMFs分量的方法,而不需像小波分解、傅里葉變換等方法那樣需要提供先驗(yàn)知識(shí)。但是傳統(tǒng)EMD方法存在模態(tài)混疊的問(wèn)題,為了克服這個(gè)問(wèn)題,Wu等[14]在EMD的基礎(chǔ)上加入噪聲輔助,稱為集成經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)。EEMD先將白噪聲添加到原始數(shù)據(jù)信號(hào)中,為時(shí)頻空間提供一個(gè)統(tǒng)一的參考框架,再對(duì)添加了白噪聲的數(shù)據(jù)進(jìn)行EMD分解,重復(fù)足夠數(shù)量的實(shí)驗(yàn),最后將多次分解得到的同頻率的IMF分量求平均作為最終的IMF分量。由于添加的白噪聲零均值的特性,經(jīng)過(guò)整體平均計(jì)算后能夠消除白噪聲對(duì)最終IMFs分量的影響。

    EEMD算法具體計(jì)算步驟如下。

    1)假設(shè)原始時(shí)間序列數(shù)據(jù)為x(t),將白噪聲nj(t)添加到原始時(shí)間序列當(dāng)中:

    xj(t)=x(t)+nj(t)j=1,2,…,M

    (1)

    式中:

    xj(t)——第j次實(shí)驗(yàn)已經(jīng)添加白噪聲的數(shù)據(jù);

    M——最大實(shí)驗(yàn)次數(shù)。

    2)在第j次實(shí)驗(yàn)中,通過(guò)EMD將添加了白噪聲的數(shù)據(jù)xj(t)分解為若干個(gè)IMF:

    (2)

    式中:

    ci,j(t)——第j次實(shí)驗(yàn)分解出的第i個(gè)IMF分量;

    rj(t)——第j次實(shí)驗(yàn)分解得到的殘差分量;

    N——每次實(shí)驗(yàn)分解得到的IMF數(shù)量。

    3)重復(fù)步驟(1)和(2),直至達(dá)到最大實(shí)驗(yàn)次數(shù),最終分解結(jié)果通過(guò)式(3)和(4)作整體平均所得:

    (3)

    (4)

    2.2 支持向量機(jī)(SVM)

    SVM通過(guò)引入核函數(shù)將輸入數(shù)據(jù)映射到高維特征空間,從而將非線性問(wèn)題轉(zhuǎn)化為線性問(wèn)題,并且通過(guò)核函數(shù)避免了高維特征空間中的點(diǎn)積運(yùn)算。SVM的目標(biāo)函數(shù)為:

    y=w·φ(x)+b

    (5)

    式中:

    x——SVM的輸入,可以是多維數(shù)據(jù);

    y——SVM的輸出;

    w——權(quán)重向量;

    b——偏置常數(shù);

    φ(x)——非線性映射函數(shù)。

    (6)

    式中:

    ε——訓(xùn)練時(shí)回歸函數(shù)允許誤差;

    C——懲罰系數(shù),C值越大對(duì)訓(xùn)練誤差大于ε的樣本的懲罰越大。

    二次優(yōu)化問(wèn)題通常建立拉格朗日方程,對(duì)方程求參數(shù)偏導(dǎo),且根據(jù)KKT條件偏導(dǎo)數(shù)為0,代回原方程可得對(duì)偶優(yōu)化問(wèn)題,通過(guò)SMO算法或者其他二次規(guī)劃算法可得支持向量回歸機(jī)預(yù)測(cè)模型:

    (7)

    式中:

    n——樣本個(gè)數(shù);

    K(xi,x)——核函數(shù),常用核函數(shù)為徑向基核函數(shù)。

    采用核函數(shù)進(jìn)行預(yù)測(cè)時(shí),核參數(shù)和懲罰參數(shù)對(duì)算法的表現(xiàn)起著非常重要的作用,核參數(shù)的值太小(0.1~1)或懲罰參數(shù)的C值太大會(huì)對(duì)訓(xùn)練集造成過(guò)學(xué)習(xí)現(xiàn)象,核參數(shù)的值太大(100~100 000)或c值太小(0.1~10)則會(huì)對(duì)訓(xùn)練集造成欠學(xué)習(xí)現(xiàn)象。本文通過(guò)粒子群尋優(yōu)算法對(duì)核參數(shù)和懲罰參數(shù)進(jìn)行尋優(yōu)。

    2.3 模型評(píng)價(jià)指標(biāo)

    本文引用3個(gè)統(tǒng)計(jì)量來(lái)評(píng)估模型的預(yù)測(cè)精度,分別為平均絕對(duì)誤差(MAE)、平均絕對(duì)百分比誤差(MAPE)和均方根誤差(RMSE),定義如下:

    (8)

    (9)

    (10)

    式中:

    n——測(cè)試集樣本數(shù)目。

    3 邊坡變形預(yù)測(cè)模型構(gòu)建

    本文首先通過(guò)EEMD算法去除邊坡監(jiān)測(cè)數(shù)據(jù)中的噪聲,將去噪之后的變形數(shù)據(jù)作為SVM模型的輸入,建立EEMD-SVM模型。在EEMD去噪過(guò)程中,EEMD分解會(huì)產(chǎn)生若干分量和殘差序列,計(jì)算各個(gè)分量與原始數(shù)據(jù)的相關(guān)系數(shù),將相關(guān)系數(shù)較小的分量當(dāng)作噪聲直接去除,重構(gòu)剩余分量生成無(wú)噪聲的邊坡變形數(shù)據(jù)。對(duì)于EEMD-SVM模型,核函數(shù)采用徑向基核函數(shù),為了提高模型的訓(xùn)練速度,將去噪后的邊坡變形數(shù)據(jù)歸一化到[0,1]區(qū)間:

    (11)

    式中:

    xi——去噪后的邊坡變形值;

    xmax——去噪后的邊坡變形最大值;

    xmin——去噪后的邊坡變形最小值;

    圖1給出了構(gòu)建邊坡變形預(yù)測(cè)模型的流程,詳細(xì)步驟如下:

    圖1 EEMD-SVM模型流程示意

    1)將邊坡變形時(shí)間序列劃分為訓(xùn)練樣本集和測(cè)試樣本集,對(duì)訓(xùn)練樣本集進(jìn)行EEMD分解,得到若干分量和殘差序列。根據(jù)Wu等[12]的推薦,加入噪聲的標(biāo)準(zhǔn)差應(yīng)為原始數(shù)據(jù)標(biāo)準(zhǔn)的0.2倍,最大實(shí)驗(yàn)次數(shù)為100。

    2)分別計(jì)算各分量和殘差序列與邊坡變形時(shí)間序列的相關(guān)系數(shù)。

    3)將相關(guān)系數(shù)小于0.2的分量作為噪聲直接去除,重構(gòu)剩余分量得到無(wú)噪聲邊坡變形時(shí)間序列。

    4)將通過(guò)去噪后重構(gòu)的邊坡變形時(shí)間序列作為SVM模型的輸入,并對(duì)SVM模型進(jìn)行訓(xùn)練,并且預(yù)測(cè)未來(lái)邊坡變形值。

    5)根據(jù)測(cè)試樣本集的邊坡變形實(shí)測(cè)值,通過(guò)式(8)~(10)評(píng)價(jià)EEMD-SVM模型的預(yù)測(cè)精度。

    4 工程實(shí)例

    溪洛渡水電站是金沙江下游河段梯級(jí)開(kāi)發(fā)規(guī)劃的第3個(gè)梯級(jí),是一座以發(fā)電為主,兼有防洪、攔沙和改善下游航運(yùn)等綜合利用效益的特大型水利水電樞紐工程。工程樞紐位于四川省雷波縣和云南省永善縣接壤的金沙江溪洛渡峽谷。水庫(kù)正常蓄水位為600 m,死水位為540 m,汛期限制水位為560 m,總庫(kù)容為126.7億m3,調(diào)節(jié)庫(kù)容為64.6億m3。水庫(kù)在蓄水過(guò)程及運(yùn)行過(guò)程中,由于庫(kù)水位的動(dòng)態(tài)變化使庫(kù)岸邊坡地質(zhì)結(jié)構(gòu)遭受強(qiáng)烈改造,從而可能導(dǎo)致滑坡災(zāi)害,帶來(lái)嚴(yán)重后果。

    為了保障水電站的正常運(yùn)行,建立了溪洛渡水電站水庫(kù)安全監(jiān)測(cè)系統(tǒng),系統(tǒng)通過(guò)GNSS技術(shù)監(jiān)測(cè)了水庫(kù)長(zhǎng)達(dá)199 km的庫(kù)岸邊坡。本文以花坪子邊坡的TP02-HPZ監(jiān)測(cè)點(diǎn)為例,收集到2018年1月1日至2018年4月30日的邊坡變形數(shù)據(jù),總共120期數(shù)據(jù),前110期為訓(xùn)練樣本,后10期為測(cè)試樣本。

    將訓(xùn)練樣本通過(guò)EEMD算法進(jìn)行分解,加入的噪聲標(biāo)準(zhǔn)差為原始邊坡序列的0.2倍,最大實(shí)驗(yàn)次數(shù)設(shè)置為100,圖2給出了訓(xùn)練樣本的EEMD分解結(jié)果。其中,第1行至第5行分別為分解后的分量IMF1~I(xiàn)MF5,反映了從高頻到低頻不同時(shí)間尺度的邊坡變形波動(dòng)特征,并且相對(duì)于原始的時(shí)間數(shù)據(jù)序列變化較為平穩(wěn),第7行為剩余分量。表1也給出了各分量與原始數(shù)據(jù)的Pearson相關(guān)系數(shù),由表1可知,IMF1、IMF2和 IMF3分量與原始數(shù)據(jù)相關(guān)系數(shù)均低于0.2,并且不具有顯著性,說(shuō)明這3個(gè)分量的變化趨勢(shì)與原始數(shù)據(jù)存在較大差異,可以將其當(dāng)作噪聲直接剔除。而IMF4、IMF5和Residual分量與原始數(shù)據(jù)都存在顯著相關(guān)性,其中IMF5和Residual分量相關(guān)系數(shù)均在0.7以上,說(shuō)明原始邊坡變形趨勢(shì)主要是由這3個(gè)分量組成的。圖3給出了剩下的IMF4、IMF5和Residual分量重構(gòu)得到的去噪后的邊坡變形時(shí)間序列,可以看出去噪后的邊坡變形時(shí)間序列更加平滑,能夠更加準(zhǔn)確的描述邊坡的真實(shí)變形趨勢(shì),為SVM模型準(zhǔn)確識(shí)別邊坡變形趨勢(shì)奠定了基礎(chǔ)。

    圖2 訓(xùn)練樣本EEMD分解結(jié)果示意

    表1 各分量與原始數(shù)據(jù)的Pearson相關(guān)系數(shù)

    圖3 去噪后的邊坡變形時(shí)間序列示意

    為了證明EEMD-SVM模型的有效性,引入SVM、BP神經(jīng)網(wǎng)絡(luò)(BP)和EEMD-BP模型作為對(duì)比。其中SVM和BP模型將原始邊坡變形時(shí)間序列作為輸入,而EEMD-SVM和EEMD-BP模型均采用去噪后的邊坡變形時(shí)間序列作為輸入。圖4給出了4種模型的預(yù)測(cè)結(jié)果,為了更加直觀詳細(xì)了解每個(gè)時(shí)刻的預(yù)測(cè)誤差情況,表2給出了4種模型的預(yù)測(cè)值。由表2可以看出,EEMD-SVM模型預(yù)測(cè)結(jié)果更加接近真實(shí)值,相較于其他3種模型,預(yù)測(cè)誤差主要集中在0.2以下。

    圖4 模型預(yù)測(cè)結(jié)果示意

    表2 模型預(yù)測(cè)值

    表3給出了4種模型的預(yù)測(cè)精度。預(yù)測(cè)精度選用平均絕對(duì)誤差(MAE)、平均絕對(duì)百分比誤差(MAPE)和均方根誤差(RMSE)這3種評(píng)價(jià)指標(biāo)對(duì)預(yù)測(cè)精度進(jìn)行分析。由表3可以看出,SVM和EEMD-SVM模型的預(yù)測(cè)精度分別優(yōu)于BP和EEMD-BP模型,說(shuō)明SVM模型相較于BP模型更適合邊坡變形預(yù)測(cè)。而EEMD-SVM和EEMD-BP模型的預(yù)測(cè)精度分別優(yōu)于SVM和BP模型,表明EEMD去噪算法可以消除邊坡變形數(shù)據(jù)的噪聲,提取出更準(zhǔn)確的邊坡變形值。其中EEMD-SVM模型的3個(gè)精度評(píng)價(jià)指標(biāo)均是最小的,相較于傳統(tǒng)SVM模型,MAE、MAPE和RMSE由0.371、0.694和0.414分別降到了0.263、0.491和0.347,MAE降低了29.1%,模型預(yù)測(cè)精度有著明顯提升。

    表3 模型預(yù)測(cè)精度

    4 結(jié)語(yǔ)

    1)邊坡的變形預(yù)測(cè)對(duì)于及時(shí)掌握邊坡的變形趨勢(shì)是至關(guān)重要的,本文通過(guò)將支持向量機(jī)與集成經(jīng)驗(yàn)?zāi)B(tài)分解方法結(jié)合,建立了可用于邊坡變形預(yù)測(cè)的EEMD-SVM模型。工程實(shí)例表明,EEMD-SVM模型具有較高的預(yù)測(cè)精度,可以為類似工程提供一種邊坡變形預(yù)測(cè)的參考方法。

    2)在邊坡變形數(shù)據(jù)采集過(guò)程中,誤差是不可避免的,EEMD算法能夠?qū)⑦吰伦冃螘r(shí)間序列進(jìn)行分解,消除無(wú)關(guān)分量,得到的去噪邊坡變形時(shí)間序列更加平穩(wěn),變化趨勢(shì)更加清晰,能夠有效應(yīng)用于邊坡變形預(yù)測(cè)。

    3)傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)模型在小樣本預(yù)測(cè)時(shí)能力有限,而支持向量機(jī)模型預(yù)測(cè)精度更高,模型運(yùn)行也更簡(jiǎn)單。

    猜你喜歡
    原始數(shù)據(jù)分量模態(tài)
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    帽子的分量
    受特定變化趨勢(shì)限制的傳感器數(shù)據(jù)處理方法研究
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    論《哈姆雷特》中良心的分量
    全新Mentor DRS360 平臺(tái)借助集中式原始數(shù)據(jù)融合及直接實(shí)時(shí)傳感技術(shù)實(shí)現(xiàn)5 級(jí)自動(dòng)駕駛
    汽車零部件(2017年4期)2017-07-12 17:05:53
    分量
    國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    由單個(gè)模態(tài)構(gòu)造對(duì)稱簡(jiǎn)支梁的抗彎剛度
    久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 蜜桃亚洲精品一区二区三区| 国内精品久久久久精免费| 国产精品国产高清国产av| 中文字幕av成人在线电影| 亚洲av美国av| а√天堂www在线а√下载| 午夜精品久久久久久毛片777| 精品人妻一区二区三区麻豆 | 欧美区成人在线视频| 91av网一区二区| 国产精品爽爽va在线观看网站| 久9热在线精品视频| 亚洲天堂国产精品一区在线| 国产精品98久久久久久宅男小说| 在线观看一区二区三区| 极品教师在线免费播放| 色视频www国产| 小说图片视频综合网站| 欧美日本视频| 男插女下体视频免费在线播放| 久9热在线精品视频| 97超级碰碰碰精品色视频在线观看| 日本三级黄在线观看| 亚洲三级黄色毛片| 大又大粗又爽又黄少妇毛片口| 免费观看精品视频网站| 中文字幕高清在线视频| 国产欧美日韩精品一区二区| 免费看光身美女| 日本五十路高清| 国产精品,欧美在线| 国产一区二区在线av高清观看| 制服丝袜大香蕉在线| 国产亚洲av嫩草精品影院| 真人一进一出gif抽搐免费| 别揉我奶头~嗯~啊~动态视频| 欧美性猛交╳xxx乱大交人| 听说在线观看完整版免费高清| 亚洲七黄色美女视频| 亚洲va在线va天堂va国产| 亚洲成人久久爱视频| 中文字幕免费在线视频6| 国产探花在线观看一区二区| 免费看光身美女| 色哟哟哟哟哟哟| 午夜福利18| 免费不卡的大黄色大毛片视频在线观看 | 99热这里只有精品一区| 国产色婷婷99| 大型黄色视频在线免费观看| 男女边吃奶边做爰视频| 亚洲成a人片在线一区二区| 成人美女网站在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 日本黄色视频三级网站网址| 亚洲av中文av极速乱 | 午夜福利视频1000在线观看| 欧美在线一区亚洲| 午夜精品久久久久久毛片777| 国产亚洲欧美98| 久久久久性生活片| 老司机福利观看| 国产精品一区二区免费欧美| 小蜜桃在线观看免费完整版高清| 婷婷亚洲欧美| 悠悠久久av| 91av网一区二区| 国产成人福利小说| 女生性感内裤真人,穿戴方法视频| 国产精品精品国产色婷婷| 国产黄色小视频在线观看| 在线观看美女被高潮喷水网站| 一边摸一边抽搐一进一小说| 禁无遮挡网站| 久久午夜亚洲精品久久| 午夜激情欧美在线| 别揉我奶头 嗯啊视频| 国产av麻豆久久久久久久| 老熟妇乱子伦视频在线观看| 一区二区三区高清视频在线| 在线观看一区二区三区| 有码 亚洲区| 欧美成人a在线观看| 中文字幕精品亚洲无线码一区| 一卡2卡三卡四卡精品乱码亚洲| 联通29元200g的流量卡| 99久久九九国产精品国产免费| 91午夜精品亚洲一区二区三区 | 蜜桃亚洲精品一区二区三区| 中文资源天堂在线| 国内揄拍国产精品人妻在线| 在线观看午夜福利视频| 窝窝影院91人妻| 男女啪啪激烈高潮av片| av视频在线观看入口| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av| 久久99热这里只有精品18| 一边摸一边抽搐一进一小说| 成人精品一区二区免费| 2021天堂中文幕一二区在线观| or卡值多少钱| 蜜桃久久精品国产亚洲av| 精品不卡国产一区二区三区| 久久精品综合一区二区三区| 国产v大片淫在线免费观看| 黄色一级大片看看| 88av欧美| 99热这里只有是精品在线观看| 国产探花极品一区二区| 在线免费观看的www视频| 国产欧美日韩一区二区精品| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 日日干狠狠操夜夜爽| 国产久久久一区二区三区| 全区人妻精品视频| 久久久久久大精品| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区人妻视频| 免费大片18禁| 黄色丝袜av网址大全| 国内精品美女久久久久久| 麻豆国产97在线/欧美| 一级黄片播放器| 啦啦啦韩国在线观看视频| 国产成人一区二区在线| 乱系列少妇在线播放| 真人一进一出gif抽搐免费| 日韩欧美国产在线观看| 日日摸夜夜添夜夜添av毛片 | www.色视频.com| 国产精品一及| 亚洲精品国产成人久久av| 一区二区三区四区激情视频 | 成年女人看的毛片在线观看| 婷婷色综合大香蕉| 国产黄色小视频在线观看| 精品人妻熟女av久视频| 国产综合懂色| 精品国内亚洲2022精品成人| 精品一区二区三区视频在线观看免费| 国产主播在线观看一区二区| 色哟哟·www| 亚洲欧美日韩高清专用| 天堂网av新在线| 看免费成人av毛片| 欧美区成人在线视频| 婷婷六月久久综合丁香| 免费人成视频x8x8入口观看| 亚洲不卡免费看| 成人永久免费在线观看视频| 成年女人看的毛片在线观看| 日韩欧美一区二区三区在线观看| 夜夜爽天天搞| 国产综合懂色| 国产高清不卡午夜福利| 一进一出抽搐gif免费好疼| 精品久久久久久久久av| 成年人黄色毛片网站| 99热精品在线国产| 波多野结衣高清无吗| 中文字幕熟女人妻在线| 国产视频内射| 中文字幕av在线有码专区| 91久久精品电影网| АⅤ资源中文在线天堂| 国产美女午夜福利| 尾随美女入室| 亚洲精品一区av在线观看| 国产伦精品一区二区三区视频9| 美女cb高潮喷水在线观看| 亚洲图色成人| .国产精品久久| 成人三级黄色视频| 黄色视频,在线免费观看| 成人性生交大片免费视频hd| 国产大屁股一区二区在线视频| 久久香蕉精品热| 波多野结衣高清作品| 男插女下体视频免费在线播放| 十八禁国产超污无遮挡网站| 国产极品精品免费视频能看的| 成人国产麻豆网| 日日干狠狠操夜夜爽| 亚洲男人的天堂狠狠| 一区二区三区激情视频| 免费不卡的大黄色大毛片视频在线观看 | 成人永久免费在线观看视频| 美女免费视频网站| 成人午夜高清在线视频| 特大巨黑吊av在线直播| 欧美潮喷喷水| 亚洲成人免费电影在线观看| 日本 av在线| 精品免费久久久久久久清纯| 亚洲av不卡在线观看| 亚洲欧美日韩东京热| 精品国产三级普通话版| 男女视频在线观看网站免费| 精品人妻一区二区三区麻豆 | 中文字幕高清在线视频| 久久精品人妻少妇| 麻豆国产97在线/欧美| 国产精品一区www在线观看 | 女生性感内裤真人,穿戴方法视频| 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 日本黄色视频三级网站网址| 99热只有精品国产| 亚洲最大成人手机在线| 久久久久久久亚洲中文字幕| 12—13女人毛片做爰片一| 男女边吃奶边做爰视频| 又黄又爽又刺激的免费视频.| 2021天堂中文幕一二区在线观| 嫁个100分男人电影在线观看| 国产伦一二天堂av在线观看| 国产国拍精品亚洲av在线观看| 老司机福利观看| 国产爱豆传媒在线观看| 高清毛片免费观看视频网站| 亚洲色图av天堂| 亚洲av日韩精品久久久久久密| 日本三级黄在线观看| 人人妻,人人澡人人爽秒播| www.www免费av| 日日摸夜夜添夜夜添av毛片 | 欧美日本视频| 有码 亚洲区| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 精品欧美国产一区二区三| 看免费成人av毛片| 国内毛片毛片毛片毛片毛片| 国产 一区 欧美 日韩| 国内精品美女久久久久久| 亚洲专区国产一区二区| 18+在线观看网站| 久久香蕉精品热| 久久久成人免费电影| а√天堂www在线а√下载| 中出人妻视频一区二区| 成人高潮视频无遮挡免费网站| avwww免费| 国产精品爽爽va在线观看网站| 又爽又黄无遮挡网站| 少妇熟女aⅴ在线视频| 欧美日韩精品成人综合77777| 一个人观看的视频www高清免费观看| 国产精品永久免费网站| av天堂中文字幕网| 一进一出好大好爽视频| 中亚洲国语对白在线视频| 欧美另类亚洲清纯唯美| 亚洲av二区三区四区| 免费av毛片视频| 午夜精品在线福利| 联通29元200g的流量卡| 91午夜精品亚洲一区二区三区 | 欧美另类亚洲清纯唯美| 国产高清不卡午夜福利| 美女高潮喷水抽搐中文字幕| 能在线免费观看的黄片| 免费一级毛片在线播放高清视频| 久久精品人妻少妇| 美女xxoo啪啪120秒动态图| 亚洲av电影不卡..在线观看| 又粗又爽又猛毛片免费看| www.www免费av| 日韩,欧美,国产一区二区三区 | 国产精品一区二区三区四区免费观看 | 黄色视频,在线免费观看| 欧美一级a爱片免费观看看| 男女做爰动态图高潮gif福利片| 久久99热6这里只有精品| 国产高潮美女av| 免费搜索国产男女视频| 内射极品少妇av片p| 国产黄a三级三级三级人| 亚洲图色成人| 亚洲美女视频黄频| 成熟少妇高潮喷水视频| 97人妻精品一区二区三区麻豆| 亚洲成人久久爱视频| 日本黄色视频三级网站网址| 高清毛片免费观看视频网站| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 亚洲午夜理论影院| 久久久久久国产a免费观看| 日本一本二区三区精品| 色吧在线观看| 国产aⅴ精品一区二区三区波| eeuss影院久久| 桃色一区二区三区在线观看| 久久欧美精品欧美久久欧美| 热99re8久久精品国产| av在线蜜桃| 亚洲欧美日韩卡通动漫| 亚洲精华国产精华液的使用体验 | 婷婷六月久久综合丁香| 欧洲精品卡2卡3卡4卡5卡区| 色噜噜av男人的天堂激情| 亚洲一区高清亚洲精品| 大型黄色视频在线免费观看| 国产人妻一区二区三区在| 国产高清不卡午夜福利| 亚洲aⅴ乱码一区二区在线播放| 免费av观看视频| 十八禁网站免费在线| 国产在视频线在精品| 精品久久久久久久久亚洲 | 精品99又大又爽又粗少妇毛片 | 尤物成人国产欧美一区二区三区| 黄色女人牲交| 国产精品永久免费网站| 亚洲综合色惰| 亚洲第一电影网av| 日韩精品有码人妻一区| 日韩欧美在线二视频| 婷婷色综合大香蕉| 中文字幕免费在线视频6| 久久人人精品亚洲av| 美女 人体艺术 gogo| 一区二区三区高清视频在线| 熟女人妻精品中文字幕| 久久热精品热| 非洲黑人性xxxx精品又粗又长| 国产爱豆传媒在线观看| 一进一出抽搐动态| 亚洲专区中文字幕在线| x7x7x7水蜜桃| av中文乱码字幕在线| 乱系列少妇在线播放| 我的老师免费观看完整版| 给我免费播放毛片高清在线观看| av在线蜜桃| 成人无遮挡网站| 欧美人与善性xxx| 九九久久精品国产亚洲av麻豆| 在线免费十八禁| 国产精品乱码一区二三区的特点| 国产精品人妻久久久影院| 日本黄色片子视频| 亚洲人成网站在线播放欧美日韩| 亚洲最大成人av| 国产在线精品亚洲第一网站| 日韩精品有码人妻一区| 欧美中文日本在线观看视频| 免费在线观看影片大全网站| 日本爱情动作片www.在线观看 | 国产亚洲精品av在线| 亚洲一级一片aⅴ在线观看| 桃色一区二区三区在线观看| 全区人妻精品视频| 日韩欧美免费精品| 欧美一级a爱片免费观看看| 97人妻精品一区二区三区麻豆| 久久久国产成人免费| 亚洲国产精品合色在线| 少妇的逼好多水| 中文字幕久久专区| 有码 亚洲区| 99久久精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 欧美最黄视频在线播放免费| 色哟哟·www| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 亚洲av中文字字幕乱码综合| 国产精品精品国产色婷婷| 搡老熟女国产l中国老女人| 极品教师在线免费播放| 亚洲最大成人中文| 一个人观看的视频www高清免费观看| 国产精品一区二区三区四区免费观看 | 亚洲va在线va天堂va国产| 在线a可以看的网站| 成人无遮挡网站| 熟女电影av网| 真人一进一出gif抽搐免费| 99久久久亚洲精品蜜臀av| 欧美黑人巨大hd| 国产午夜精品论理片| 国产精品人妻久久久影院| 中文字幕免费在线视频6| 中文字幕av成人在线电影| 成人无遮挡网站| 人人妻人人看人人澡| 婷婷色综合大香蕉| 人妻少妇偷人精品九色| 日韩精品中文字幕看吧| 国产欧美日韩精品亚洲av| 亚洲av电影不卡..在线观看| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 一个人看视频在线观看www免费| 久久人人爽人人爽人人片va| 免费黄网站久久成人精品| 亚洲精品在线观看二区| 给我免费播放毛片高清在线观看| 欧美激情久久久久久爽电影| 欧美一区二区国产精品久久精品| av黄色大香蕉| 97人妻精品一区二区三区麻豆| bbb黄色大片| 亚洲第一区二区三区不卡| 夜夜爽天天搞| 亚洲性久久影院| 天堂√8在线中文| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 日日夜夜操网爽| 色综合站精品国产| 久久精品国产亚洲av涩爱 | 国产精品爽爽va在线观看网站| 亚洲无线观看免费| 啪啪无遮挡十八禁网站| 亚洲电影在线观看av| 俺也久久电影网| 1000部很黄的大片| 两人在一起打扑克的视频| 国产真实伦视频高清在线观看 | 久久人人爽人人爽人人片va| netflix在线观看网站| 亚洲av一区综合| 色综合色国产| 男插女下体视频免费在线播放| 两个人的视频大全免费| 久久99热6这里只有精品| 日韩中字成人| 一本精品99久久精品77| 少妇的逼好多水| 亚洲午夜理论影院| 91久久精品国产一区二区成人| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美清纯卡通| 22中文网久久字幕| 老司机午夜福利在线观看视频| 少妇丰满av| 老师上课跳d突然被开到最大视频| 国产成年人精品一区二区| 午夜免费成人在线视频| 精品午夜福利视频在线观看一区| 免费av毛片视频| 中文字幕久久专区| 偷拍熟女少妇极品色| 老熟妇乱子伦视频在线观看| 国内精品宾馆在线| 99国产极品粉嫩在线观看| 91麻豆av在线| 美女高潮喷水抽搐中文字幕| 看片在线看免费视频| 一个人免费在线观看电影| 精品福利观看| 毛片女人毛片| 婷婷精品国产亚洲av在线| 性欧美人与动物交配| 久久亚洲精品不卡| 免费观看精品视频网站| 日韩欧美免费精品| 女同久久另类99精品国产91| 18禁黄网站禁片午夜丰满| 亚洲国产精品sss在线观看| 亚洲第一电影网av| 国产黄a三级三级三级人| 国产精品亚洲美女久久久| 丝袜美腿在线中文| 色哟哟哟哟哟哟| 一级黄片播放器| 少妇丰满av| 亚洲在线观看片| 最近视频中文字幕2019在线8| 又黄又爽又免费观看的视频| 亚洲av免费高清在线观看| 亚洲精品在线观看二区| 欧美国产日韩亚洲一区| 国产一区二区激情短视频| 国产探花极品一区二区| 禁无遮挡网站| 欧美日本亚洲视频在线播放| 又爽又黄a免费视频| 波多野结衣高清作品| 亚洲无线观看免费| 午夜老司机福利剧场| 中文字幕av成人在线电影| 99久国产av精品| 亚洲一级一片aⅴ在线观看| 少妇的逼水好多| 亚洲aⅴ乱码一区二区在线播放| 成年女人看的毛片在线观看| 最新在线观看一区二区三区| 最新中文字幕久久久久| 国产在线精品亚洲第一网站| ponron亚洲| 桃色一区二区三区在线观看| 国产免费av片在线观看野外av| 人妻丰满熟妇av一区二区三区| 日韩强制内射视频| 成人国产综合亚洲| 欧美色视频一区免费| 午夜免费成人在线视频| 日本三级黄在线观看| АⅤ资源中文在线天堂| 一个人看视频在线观看www免费| 日韩欧美精品免费久久| 精品人妻1区二区| 日本撒尿小便嘘嘘汇集6| 日本与韩国留学比较| 91麻豆av在线| 一区福利在线观看| 国产久久久一区二区三区| 欧美bdsm另类| 如何舔出高潮| 精品久久久久久,| 赤兔流量卡办理| 日本在线视频免费播放| 国产伦精品一区二区三区视频9| 久久精品国产亚洲av天美| 久久精品国产亚洲av香蕉五月| 色尼玛亚洲综合影院| 久久精品国产清高在天天线| 女的被弄到高潮叫床怎么办 | 两人在一起打扑克的视频| 色在线成人网| 精品午夜福利视频在线观看一区| 国产不卡一卡二| 内地一区二区视频在线| 99久久精品热视频| 真实男女啪啪啪动态图| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 舔av片在线| 特级一级黄色大片| 日日干狠狠操夜夜爽| 亚洲经典国产精华液单| 国产精品久久视频播放| 夜夜看夜夜爽夜夜摸| 别揉我奶头~嗯~啊~动态视频| 久久精品人妻少妇| 日本黄色片子视频| 少妇丰满av| 国产又黄又爽又无遮挡在线| 波多野结衣高清作品| 中文字幕av在线有码专区| 岛国在线免费视频观看| 国产老妇女一区| 日韩 亚洲 欧美在线| 午夜影院日韩av| 精品99又大又爽又粗少妇毛片 | 精品久久久噜噜| 亚洲七黄色美女视频| 国产女主播在线喷水免费视频网站 | 又爽又黄a免费视频| 嫩草影院精品99| 一级黄色大片毛片| 国产亚洲欧美98| 亚洲18禁久久av| 亚洲va日本ⅴa欧美va伊人久久| 1024手机看黄色片| 久久精品夜夜夜夜夜久久蜜豆| 窝窝影院91人妻| 色精品久久人妻99蜜桃| 国产男靠女视频免费网站| 一本精品99久久精品77| 观看美女的网站| 99久久久亚洲精品蜜臀av| 蜜桃久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 亚洲成人免费电影在线观看| 亚洲精品456在线播放app | 搞女人的毛片| 在现免费观看毛片| 国内揄拍国产精品人妻在线| 国产精品人妻久久久久久| 桃红色精品国产亚洲av| 亚洲精品成人久久久久久| 国产激情偷乱视频一区二区| 国模一区二区三区四区视频| 色在线成人网| 可以在线观看毛片的网站| 男女视频在线观看网站免费| 日本五十路高清| av黄色大香蕉| 搡老妇女老女人老熟妇| 一进一出抽搐gif免费好疼| 国内久久婷婷六月综合欲色啪| 深爱激情五月婷婷| 12—13女人毛片做爰片一| 国产精品永久免费网站| 日韩欧美免费精品| 日韩欧美 国产精品| 色尼玛亚洲综合影院| 一级毛片久久久久久久久女| 一区二区三区高清视频在线| 国产av一区在线观看免费| 俄罗斯特黄特色一大片| 亚洲一区二区三区色噜噜| 一卡2卡三卡四卡精品乱码亚洲| or卡值多少钱| 色哟哟哟哟哟哟| 色尼玛亚洲综合影院| 免费观看人在逋| 91久久精品电影网| 级片在线观看| 狠狠狠狠99中文字幕| 此物有八面人人有两片| 男女那种视频在线观看| 一个人观看的视频www高清免费观看| 亚洲精品一区av在线观看| 国产探花在线观看一区二区| 精品久久久久久久末码|