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

    基于FEEMD-GRU-FC模型的滑坡位移預(yù)測(cè)

    2024-12-31 00:00:00何清李麗琳林子安
    人民長(zhǎng)江 2024年7期
    關(guān)鍵詞:八字降雨滑坡

    摘要:

    水庫(kù)庫(kù)區(qū)滑坡同時(shí)受到降雨和水庫(kù)水位周期性變化的影響,這會(huì)干擾滑坡位移預(yù)測(cè)模型在該類區(qū)域的預(yù)測(cè)性能。為此,提出了一種針對(duì)環(huán)境周期性變化特點(diǎn)的FEEMD-GRU-FC(Fast Ensemble Empirical Mode Decomposition-Gated Recurrent Unit-Fully Connected Layer)滑坡位移預(yù)測(cè)模型。該模型充分考慮位移數(shù)據(jù)不同頻率之間的噪聲影響,將位移分解為不同頻率的子序列分別預(yù)測(cè),采用PCC(Pearson Correlation Coefficient)方法選擇影響較大的外部環(huán)境因素,添加全連接層完善GRU模型提取環(huán)境數(shù)據(jù)中周期性變化特征的能力,模型結(jié)合三峽庫(kù)區(qū)八字門(mén)滑坡的實(shí)際數(shù)據(jù)進(jìn)行模擬仿真,使用不同預(yù)測(cè)模型和評(píng)價(jià)指標(biāo)進(jìn)行多角度對(duì)比驗(yàn)證。實(shí)驗(yàn)結(jié)果表明:所提模型在八字門(mén)滑坡的預(yù)測(cè)效果要優(yōu)于其他3種模型,其擬合度達(dá)到了98.6%,3種不同誤差指標(biāo)均為最低,證明了該模型可以準(zhǔn)確預(yù)測(cè)庫(kù)區(qū)滑坡位移。研究成果可為防災(zāi)減災(zāi)工作提供指導(dǎo)。

    關(guān)" 鍵" 詞:

    滑坡位移預(yù)測(cè); 皮爾遜相關(guān)系數(shù); 門(mén)控循環(huán)單元; 全連接層; 非線性; 八字門(mén)滑坡; 三峽庫(kù)區(qū)

    中圖法分類號(hào): TV697.23; P642.22

    文獻(xiàn)標(biāo)志碼: A

    DOI:10.16232/j.cnki.1001-4179.2024.07.014

    收稿日期:

    2024-01-12;接受日期:

    2024-04-22

    基金項(xiàng)目:

    國(guó)家自然科學(xué)基金項(xiàng)目(U23A20280);廣西科技重大專項(xiàng)(桂科AA23062038);廣西重點(diǎn)研發(fā)項(xiàng)目(桂科AB21196041)

    作者簡(jiǎn)介:

    何" 清,女,工程師,主要從事智能交通和滑坡預(yù)警研究。E-mail:2606261503@qq.com

    通信作者:

    林子安,男,工程師,博士,主要從事人工智能和滑坡預(yù)警研究。E-mail:20031102010@mails.guet.edu.cn

    Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

    文章編號(hào):1001-4179(2024) 07-0108-07

    引用本文:

    何清,李麗琳,林子安.

    基于FEEMD-GRU-FC模型的滑坡位移預(yù)測(cè)

    [J].人民長(zhǎng)江,2024,55(7):108-114.

    0" 引 言

    滑坡是對(duì)人類及其生存環(huán)境造成危害或破壞的復(fù)雜而常見(jiàn)的地質(zhì)自然災(zāi)害[1],因其發(fā)生頻率高、影響范圍廣、持續(xù)時(shí)間長(zhǎng),已成為中國(guó)最嚴(yán)重的地質(zhì)災(zāi)害之一。根據(jù)中國(guó)2021年國(guó)家統(tǒng)計(jì)年鑒,2020年中國(guó)共發(fā)生地質(zhì)災(zāi)害7 840起,其中滑坡4 810起,占中國(guó)地質(zhì)災(zāi)害總數(shù)的61.3%,處理滑坡問(wèn)題迫在眉睫。然而,由于滑坡形成機(jī)理和誘發(fā)因素的復(fù)雜性,對(duì)滑坡的預(yù)測(cè)仍是一項(xiàng)艱巨的任務(wù)。而滑坡位移預(yù)測(cè)的研究是解決滑坡問(wèn)題、更好認(rèn)識(shí)滑坡運(yùn)動(dòng)的有效手段[2]。

    滑坡可以被認(rèn)為是一種非線性的系統(tǒng)。滑坡的發(fā)生可以有多種原因,如地質(zhì)、降雨或者人類活動(dòng)等。關(guān)于滑坡位移預(yù)測(cè)的研究較多,現(xiàn)有滑坡位移預(yù)測(cè)模型可分為基于物理機(jī)理的模型和基于數(shù)據(jù)驅(qū)動(dòng)的模型[3]。但是基于物理原理的模型由于滑坡自身的復(fù)雜性,并沒(méi)有取得很好的預(yù)測(cè)效果[4]。近年來(lái)關(guān)于滑坡位移的預(yù)測(cè)模型大多都是基于數(shù)據(jù)驅(qū)動(dòng)的。例如Wu等提出新的灰色預(yù)測(cè)模型,實(shí)現(xiàn)滑坡的位移預(yù)測(cè)[5]。李麗敏等將滑坡位移當(dāng)作時(shí)間序列數(shù)據(jù)處理,使用原始GRU模型進(jìn)行預(yù)測(cè),但忽略了外部環(huán)境數(shù)據(jù)的周期性變化[6]。袁乾博等基于時(shí)間序列的傳統(tǒng)GRU模型對(duì)三峽庫(kù)區(qū)的八字門(mén)滑坡進(jìn)行預(yù)測(cè),對(duì)比一般的GRNN模型,其預(yù)測(cè)效果有了明顯提高[7]。段功豪等設(shè)計(jì)了一種基于多演化階段插值的滑坡位移預(yù)測(cè),結(jié)果表明該方法得到的誤差更低,相較于傳統(tǒng)方法更具優(yōu)勢(shì)[8]。

    然而,滑坡位移通常隨著時(shí)間變化而逐漸增加[9],傳統(tǒng)的靜態(tài)預(yù)測(cè)模型沒(méi)有考慮到滑坡位移是一種動(dòng)態(tài)的時(shí)間序列數(shù)據(jù)[10],也忽略了滑坡容易受外部環(huán)境周期性影響的特點(diǎn)[11]。針對(duì)以上情況,本文提出一種針對(duì)周期性時(shí)間序列數(shù)據(jù)的FEEMD-GRU-FC滑坡位移預(yù)測(cè)模型。該模型使用FEEMD算法將滑坡位移及環(huán)境變量分解得到不同頻率的子序列,通過(guò)PCC方法選出對(duì)滑坡位移影響較大的環(huán)境因素,根據(jù)GRU模型善于處理滑坡位移與環(huán)境因素之間非線性關(guān)系的優(yōu)點(diǎn),以及針對(duì)滑坡外部環(huán)境數(shù)據(jù)周期性變化特點(diǎn),給GRU模型增加了全連接層以此提高模型提取周期性特征的能力。最后使用三峽庫(kù)區(qū)的八字門(mén)滑坡相關(guān)數(shù)據(jù)進(jìn)行仿真驗(yàn)證。

    1" 模型介紹

    1.1" 模型流程

    FEEMD-GRU-FC滑坡位移的預(yù)測(cè)過(guò)程包括5個(gè)步驟,圖1展示了完整的實(shí)現(xiàn)過(guò)程。

    步驟1:獲取滑坡位移以及滑坡外部的環(huán)境數(shù)據(jù)。

    步驟2:使用FEEMD算法對(duì)滑坡位移和環(huán)境數(shù)據(jù)進(jìn)行分解,得到多個(gè)不同頻率的子序列。

    步驟3:基于每個(gè)滑坡位移子序列與所有的環(huán)境數(shù)據(jù)子序列計(jì)算PCC值,選擇相關(guān)性較大的環(huán)境數(shù)據(jù)子序列作為滑坡位移子序列的輸入變量。

    步驟4:GRU-FC模型對(duì)每個(gè)滑坡位移子序列進(jìn)行分別預(yù)測(cè),所有預(yù)測(cè)結(jié)果相加得到最終的滑坡預(yù)測(cè)位移。

    步驟5:使用MAE、RMSE、MAPE和R2評(píng)估FEEMD-GRU-FC模型的預(yù)測(cè)性能。

    1.2" FEEMD模型

    FEEMD(Fast Ensemble Empirical Mode Decomposition)模型由Wang等[12]在2014年提出,是一種改進(jìn)的EEMD算法,該算法可以有效提高EEMD的實(shí)時(shí)計(jì)算性能[13]。作為一種輔助信號(hào)處理方法,F(xiàn)EEMD將高斯白噪聲引入到原始時(shí)間序列中再應(yīng)用,其非常適合于非平穩(wěn)和非線性信號(hào)的分解。FEEMD方法通過(guò)優(yōu)化停止條件減少篩選次數(shù)來(lái)提升算法的執(zhí)行效率,且能較好地對(duì)信號(hào)進(jìn)行去噪[14],在有效去除噪聲的同時(shí)也保證了有用信號(hào)的完整性[15]。FEEMD的分解過(guò)程遵循以下基本步驟:

    (1) 設(shè)置FEEMD的集合數(shù)I=100;

    (2) 在原始數(shù)據(jù)上添加選定振幅的白噪聲序列,得到生成的序列x(t);

    (3) 利用3次樣條插值將x(t)的局部極大值點(diǎn)和局部極小值點(diǎn)聯(lián)系起來(lái),分別得到上下包絡(luò)Ui(t)和Li(t);

    (4) 計(jì)算上下包絡(luò)Ui(t)和Li(t)的平均值mi(t),并指定xi(t)和mi(t)之間的差值為hi(t);

    hi(t)=xi(t)-mi(t)

    (1)

    如果hi(t)滿足,整個(gè)數(shù)據(jù)集中零點(diǎn)和極值點(diǎn)的個(gè)數(shù)相差1以內(nèi)或者上下包絡(luò)Ui(t)和Li(t)在任意點(diǎn)的平均值為零任意一個(gè)條件,那么算法會(huì)繼續(xù)執(zhí)行,否則會(huì)重復(fù)進(jìn)入上述過(guò)程,直到hi(t)滿足要求。

    得到hi(t)后,計(jì)算出余項(xiàng)數(shù)據(jù)ri(t):

    ri(t)=xi(t)-hi(t)

    (2)

    重復(fù)以上步驟,定義標(biāo)準(zhǔn)偏差SD用作停止條件,當(dāng)SD小于0.2時(shí)會(huì)停止上述步驟并輸出計(jì)算得到的N個(gè)IMF分量[16],SD的公式如下:

    SD=Ni=1|hiM-1(t)-hiM(t)|2hiM-1(t)2

    (3)

    式中:M是篩選次數(shù)。

    xi(t)表示為

    xi(t)=Ni=1IMF(i)+ri(t)

    (4)

    式中:N表述IMF的數(shù)量,r表示余項(xiàng)數(shù)據(jù)。

    1.3" PCC方法

    PCC(Pearson Correlation Coefficient)是測(cè)量變量之間關(guān)系的常用方法之一,目前已被廣泛用于計(jì)算事物之間的相關(guān)性 [17]。通常它是基于協(xié)方差用矩陣的數(shù)據(jù)計(jì)算來(lái)評(píng)價(jià)兩個(gè)向量之間的相互強(qiáng)弱關(guān)系,其相關(guān)系數(shù)為1或者-1時(shí),相關(guān)性越強(qiáng);當(dāng)相關(guān)系數(shù)接近0時(shí)則表示較弱的相關(guān)性。假設(shè)有兩個(gè)變量ɑi和ɑj,則兩者相關(guān)系數(shù)的公式可表示為

    p=Cov(ɑi,ɑj)Var(ɑi)×Var(ɑj)

    (5)

    式中:Cov(ɑi,ɑj)是協(xié)方差,Var(ɑi)是ɑi的方差,Var(ɑj)是ɑj的方差。

    當(dāng)相關(guān)系數(shù)大于1時(shí),一般認(rèn)為p≤0.40為弱相關(guān),p在0.40~0.70之間為中等相關(guān),p在0.70~1.00之間為強(qiáng)相關(guān)或高相關(guān),p≥0.9為極高相關(guān)[18]。本文也根據(jù)這一相關(guān)性劃分選擇中等相關(guān)及以上的環(huán)境因素參與滑坡位移的預(yù)測(cè),認(rèn)為中等相關(guān)及以上的因素對(duì)滑坡位移產(chǎn)生了實(shí)際影響。

    1.4" GRU模型

    GRU(Gated Recurrent Unit)是一種新型的卷積神經(jīng)網(wǎng)絡(luò)模型,是RNN模型和LSTM模型的改進(jìn)模型,優(yōu)化了LSTM模型的門(mén)函數(shù),將神經(jīng)元狀態(tài)和隱藏狀態(tài)混合[19]。GRU模型在保留了RNN處理時(shí)序數(shù)據(jù)優(yōu)勢(shì)的基礎(chǔ)上,解決了長(zhǎng)程依賴問(wèn)題,能夠使模型建立長(zhǎng)距離時(shí)刻的依賴關(guān)系[20],在每個(gè)循環(huán)單元計(jì)算中可以自適應(yīng)捕獲不同時(shí)間尺度的相關(guān)性[21]。GRU模型也善于處理時(shí)間序列數(shù)據(jù)以及數(shù)據(jù)之間的非線性關(guān)系,可以有效緩解RNN模型梯度消失和梯度爆炸的問(wèn)題,并且參數(shù)少于LSTM模型,可以有效減少處理數(shù)據(jù)的時(shí)間,其收斂速度快并且運(yùn)行成本低[22]。GRU模型的結(jié)構(gòu)如圖2所示[23],公式為

    rt=σ(wr·[ht-1,xt])(6)

    zt=σ(wz·[ht-1,xt])(7)

    Ct=(wc·[rt×ht-1,xt])(8)

    ht=(I-zt)×ht-1+zt×Ct(9)

    yt=σ(wo·ht)(10)

    式中:xt是當(dāng)前時(shí)刻的輸入數(shù)據(jù),ht-1是前一個(gè)時(shí)刻的狀態(tài)儲(chǔ)存變量,ht是當(dāng)前時(shí)刻的狀態(tài)儲(chǔ)存變量,rt是更新門(mén)的狀態(tài),zt是重置門(mén)的狀態(tài),Ct是當(dāng)前候選集的狀態(tài),yt是當(dāng)前時(shí)刻的輸出數(shù)據(jù),wr是更新門(mén)的權(quán)重,wz是重置門(mén)的權(quán)重,wc是候選集的權(quán)重,wo是yt和xt與ht-1組成矩陣的權(quán)重,σ是sigmoid函數(shù),

    SymbolFC@ 是tanh函數(shù),[]表示直接連接兩個(gè)矩陣向量,·表示兩個(gè)矩陣向量相乘。

    1.5" GRU-FC模型

    滑坡主要受降雨影響,而每年的降雨通常具有周期性和重復(fù)性。所以為了提高模型處理周期性數(shù)據(jù)的能力,對(duì)GRU模型進(jìn)行改進(jìn),在GRU模型的隱藏層之后增加全連接層(Fully Connected Layer),使得全連接層的每個(gè)節(jié)點(diǎn)在預(yù)測(cè)時(shí)都不會(huì)忽略前面和后面的降雨信息。最終該模型使用輸出單元計(jì)算預(yù)測(cè)的滑坡位移數(shù)據(jù)。目前全連接層的組合模型已應(yīng)用于場(chǎng)景圖像定位[24]、電網(wǎng)異常檢測(cè)[25]等多個(gè)領(lǐng)域。GRU-FC模型的結(jié)構(gòu)如圖3所示。

    1.6" 評(píng)價(jià)指標(biāo)

    本文為了準(zhǔn)確估計(jì)模型的預(yù)測(cè)能力,采用評(píng)價(jià)指標(biāo)用于評(píng)價(jià)模型的性能和精度[2],采用平均絕對(duì)誤差(MAE)、均方根誤差(RMSE)、平均絕對(duì)百分比誤差(MAPE)和決定系數(shù)(R2)4種不同的評(píng)價(jià)指標(biāo)。MAE是一種基本的評(píng)價(jià)指標(biāo),反映了預(yù)測(cè)值和實(shí)際數(shù)據(jù)的整體差距。RMSE對(duì)異常值很敏感,如果預(yù)測(cè)數(shù)據(jù)不合理,RMSE值將受到很大影響。MAPE從整體考慮了預(yù)測(cè)值、實(shí)際數(shù)據(jù)以及誤差之間的關(guān)系。R2則表示預(yù)測(cè)數(shù)據(jù)擬合實(shí)際數(shù)據(jù)的程度。

    MAE=1nni=1|zi-yi|

    (11)

    RMSE=1nni=1(zi-yi)2

    (12)

    MAPE=1nni=1zi-yiyi

    (13)

    R2=1-ni=1(zi-yi)2ni=1(pi-yi)2

    (14)

    式中:n表示數(shù)據(jù)的數(shù)量;z={z1,z2,…,zn},表示預(yù)測(cè)的數(shù)據(jù);y={y1,y2,…,yn},表示實(shí)際的測(cè)量數(shù)據(jù);p={p1,p2,…,pn},表示實(shí)際測(cè)量數(shù)據(jù)的平均值。RMSE、MAE和MAPE這3個(gè)評(píng)價(jià)指標(biāo)都是誤差指標(biāo),其值越小,模型的預(yù)測(cè)性能和準(zhǔn)確性越優(yōu)秀;R2是擬合系數(shù),其值越大,模型的預(yù)測(cè)性能和準(zhǔn)確性越優(yōu)秀。

    2" 案例仿真

    2.1" 數(shù)據(jù)來(lái)源

    八字門(mén)滑坡位于三峽庫(kù)區(qū)長(zhǎng)江北岸支流香溪河右岸,地處湖北省秭歸縣歸州鎮(zhèn),距離三峽大壩31 km。香溪河南北走向,幾乎橫穿長(zhǎng)江,三峽水庫(kù)淹沒(méi)滑坡前緣。滑坡區(qū)內(nèi)共有3個(gè)GPS變形監(jiān)測(cè)點(diǎn)(ZG111、ZG110、ZG112)。本次研究選擇ZG111作為八字門(mén)滑坡的監(jiān)測(cè)點(diǎn),因?yàn)樵撐恢玫默F(xiàn)有監(jiān)測(cè)數(shù)據(jù)序列最長(zhǎng),滑坡變形量最大。八字門(mén)滑坡地形如圖4所示。

    本次研究采用八字門(mén)滑坡60組數(shù)據(jù)。每一組都包含滑坡位移、降雨和水庫(kù)水位數(shù)據(jù)。數(shù)據(jù)開(kāi)始于2008年1月,結(jié)束于2012年12月,采集頻率為每月一次,數(shù)據(jù)如圖5所示。

    2.2" FEEMD算法分解原始數(shù)據(jù)

    由圖5可以看出,滑坡位移在雨季的前幾個(gè)月和水庫(kù)水位下降時(shí)開(kāi)始迅速增加。雨季過(guò)后,水庫(kù)水位穩(wěn)定后,滑坡位移變化基本結(jié)束,所以降雨和水庫(kù)水位的變化是滑坡位移變化的誘發(fā)因素。本文對(duì)滑坡位移、降雨和水庫(kù)水位進(jìn)行FEEMD分解,將各數(shù)據(jù)分解得到不同頻率的4個(gè)子序列,以便在之后預(yù)測(cè)單個(gè)子序列的過(guò)程中減少其他頻率數(shù)據(jù)帶來(lái)的誤差,分解如圖6所示。

    2.3" PCC方法計(jì)算相關(guān)性

    滑坡位移、降雨和水庫(kù)水位經(jīng)過(guò)FEEMD算法分解后,得到了多個(gè)不同頻率的子序列。但并不是所有的降雨子序列和水庫(kù)水位子序列對(duì)滑坡位移都會(huì)產(chǎn)生較大的影響,所以為了預(yù)測(cè)的準(zhǔn)確性,本文采用PCC方法計(jì)算滑坡位移與環(huán)境影響因素的相關(guān)性,分別選出相關(guān)性較大的環(huán)境因素子序列作為各個(gè)滑坡位移子序列的輸入數(shù)據(jù),PCC方法計(jì)算結(jié)果如表1所列。

    可以從表1中看出,降雨相對(duì)水庫(kù)水位而言會(huì)對(duì)滑坡位移產(chǎn)生更大的影響,這可能是因?yàn)榻涤陼?huì)直接沖刷滑坡表面,提高滑坡含水率增加滑坡自身重量。根據(jù)表1的數(shù)據(jù)以及PCC的相關(guān)性等級(jí)劃分,選擇PCC≥0.4的環(huán)境因素作為各滑坡位移子序列預(yù)測(cè)模型的輸入變量,IMF1、IMF2、IMF3和r的輸入變量個(gè)數(shù)分別為6,7,6,8個(gè)。

    2.4" FEEMD-GRU-FC模型預(yù)測(cè)滑坡位移

    將60個(gè)月的數(shù)據(jù)分為兩個(gè)部分,前48個(gè)月數(shù)據(jù)作為訓(xùn)練集訓(xùn)練FEEMD-GRU-FC預(yù)測(cè)模型,后12個(gè)月的數(shù)據(jù)作為驗(yàn)證集對(duì)模型進(jìn)行驗(yàn)證,每個(gè)滑坡位

    移子序列的預(yù)測(cè)過(guò)程完全一致,最后將各滑坡位移子序列的預(yù)測(cè)結(jié)果相加,得到最終的預(yù)測(cè)結(jié)果。模型的各參數(shù)進(jìn)行如下設(shè)置:集合數(shù)I設(shè)置為100[12];偏差SD停止條件設(shè)置為0.2[16];PCC大于0.4為選擇相關(guān)影響因素條件[18];由于八字門(mén)滑坡所處區(qū)域受到季節(jié)性降雨氣候影響,所以將時(shí)間步長(zhǎng)設(shè)置為12個(gè)月;經(jīng)過(guò)多次實(shí)驗(yàn),在仿真過(guò)程中發(fā)現(xiàn)學(xué)習(xí)率與隱藏層節(jié)點(diǎn)數(shù)不宜過(guò)大或過(guò)小,將其分別設(shè)置為0.01與50時(shí)模型效果最佳。預(yù)測(cè)結(jié)果與實(shí)際結(jié)果對(duì)比如圖7~8所示,從圖7~8中可以看出本文提出的模型能較準(zhǔn)確地預(yù)測(cè)八字門(mén)的滑坡位移。

    3" 討 論

    為了證明GRU-FC模型預(yù)測(cè)滑坡位移的優(yōu)越性能,使用同為深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)的RNN-FC模型、LSTM-FC模型和原始GRU模型對(duì)八字門(mén)滑坡位移進(jìn)行預(yù)測(cè),并與GRU-FC模型進(jìn)行對(duì)比。為了便于比較,本文將GRU-FC模型、RNN-FC模型、LSTM-FC模型和GRU模型的時(shí)間步長(zhǎng)統(tǒng)一設(shè)置為12,學(xué)習(xí)率設(shè)置為0.01,隱藏層節(jié)點(diǎn)數(shù)設(shè)置為50。預(yù)測(cè)結(jié)果如圖9所示。

    為了更直觀地對(duì)比這幾種滑坡位移預(yù)測(cè)模型的性能,列出這幾種模型的兩種評(píng)價(jià)指標(biāo)結(jié)果見(jiàn)表2。

    從表2中可以看出,無(wú)論是MAE還是其他3種指標(biāo),GRU-FC模型在這4種模型里面都是最優(yōu)秀的。GRU-FC與RNN-FC和LSTM-FC模型相比,擁有更好的預(yù)測(cè)性能,這是因?yàn)镚RU-FC模型解決了RNN-FC模型存在的梯度消失和梯度爆炸問(wèn)題;而GRU-FC優(yōu)于LSTM-FC模型,則是因?yàn)長(zhǎng)STM比GRU模型的結(jié)構(gòu)更加復(fù)雜,參數(shù)量更多,所以在訓(xùn)練模型的時(shí)候成本更高,在數(shù)據(jù)量不多的情況下GRU-FC比LSTM-FC模型更占優(yōu)勢(shì)。GRU-FC模型的預(yù)測(cè)效果比GRU模型要好,是因?yàn)槿B接層可以同時(shí)收集到滑坡發(fā)展過(guò)程中每個(gè)周期有著近似變化的數(shù)據(jù),提高了GRU-FC提取周期性數(shù)據(jù)特征的能力,在面對(duì)八字門(mén)滑坡周期性的降雨和水庫(kù)水位變化時(shí)會(huì)具備更優(yōu)秀的預(yù)測(cè)能力。

    本文提出的FEEMD-GRU-FC模型雖然可以對(duì)三峽庫(kù)區(qū)的八字門(mén)滑坡位移進(jìn)行有效預(yù)測(cè),但仍有一些待改進(jìn)的地方。模型沒(méi)有經(jīng)過(guò)在其他地區(qū)的滑坡或者其他類型滑坡位移數(shù)據(jù)的仿真驗(yàn)證,所以之后的工作還需要將該模型與其他滑坡數(shù)據(jù)結(jié)合使用,以提高模型的適用性。

    4" 結(jié) 語(yǔ)

    滑坡位移預(yù)測(cè)的研究已經(jīng)進(jìn)行了很長(zhǎng)時(shí)間,但仍然是一個(gè)具有挑戰(zhàn)性的研究,威脅著世界各地人民的生命和財(cái)產(chǎn)安全。針對(duì)這一挑戰(zhàn),本文提出一種針對(duì)非線性周期性時(shí)間序列數(shù)據(jù)的FEEMD-GRU-FC庫(kù)區(qū)滑坡位移預(yù)測(cè)模型。該模型采用FEEMD算法快速將滑坡位移與環(huán)境因素分解成不同頻率的子序列,再結(jié)合GRU模型和全連接層的優(yōu)點(diǎn),有效處理滑坡位移及其外部環(huán)境的非線性關(guān)系及周期性變化,準(zhǔn)確地預(yù)測(cè)滑坡位移。并且與其他3種人工智能深度學(xué)習(xí)模型進(jìn)行對(duì)比,仿真結(jié)果也證明了本文提出的FEEMD-GRU-FC模型能夠有效地預(yù)測(cè)滑坡位移,為同類型非線性周期性變化滑坡的預(yù)防和治理提供了理論依據(jù)。未來(lái)可在實(shí)際項(xiàng)目過(guò)程中對(duì)模型進(jìn)行修改完善,進(jìn)一步提升模型的可行性和實(shí)用性,驗(yàn)證所提模型對(duì)于實(shí)際工程生產(chǎn)項(xiàng)目的幫助。此外,F(xiàn)EEMD-GRU-FC模型還可應(yīng)用于其他預(yù)測(cè)領(lǐng)域,如降雨預(yù)測(cè)、發(fā)電量預(yù)測(cè)等,協(xié)助決策者不斷改進(jìn)合理判斷的過(guò)程。

    參考文獻(xiàn):

    [1]" LIAN C,ZENG Z G,WANG X P,et al.Landslide displacement interval prediction using lower upper bound estimation method with pre-trained random vector functional link network initialization[J].Neural Networks,2020,130:286-296.

    [2]" LIN Z A,SUN X Y,JI Y F.Landslide displacement prediction model using time series analysis method and modified LSTM model[J].Electronics,2022,11:1519.

    [3]" LIU Z Q,GUO D,LACASSE S,et al.Algorithms for intelligent prediction of landslide displacements[J].Journal of Zhejiang University (Science A),2020,21(6):412-429.

    [4]" DU H,SONG D Q,CHEN Z,et al.Prediction model oriented for landslide displacement with step-like curve by applying ensemble empirical mode decomposition and the PSO-ELM method[J].Journal of Cleaner Production,2020,270:122248.

    [5]" WU L Z,LI S H,HUANG R Q,et al.A new grey prediction model and its application to predicting landslide displacement[J].Applied Soft Computing Journal,2020,45:5347-5361.

    [6]" 李麗敏,張明岳,溫宗周,等.基于奇異譜分析法和長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)組合模型的滑坡位移預(yù)測(cè)[J].信息與控制,2021,50(4):459-469,482.

    [7]" 袁乾博,肖詩(shī)榮,李春霞,等.基于GLSSVM模型的三門(mén)洞滑坡變形預(yù)測(cè)研究[J].人民長(zhǎng)江,2020,51(7):130-135.

    [8]" 段功豪,馬迅,付杰.基于多演化階段插值方法的滑坡預(yù)測(cè)研究[J].人民長(zhǎng)江,2023,54(1):106-110,118.

    [9]" ZHANG J R,TANG H M,TANNANT D D,et al.Combined forecasting model with CEEMD-LCSS reconstruction and the ABC-SVR method for landslide displacement prediction[J].Journal of Cleaner Production,2021,293:126205.

    [10]YANG B B,YIN K L,LACASSE S,et al.Time series analysis and long short-term memory neural network to predict landslide displacement[J].Landslides,2019,16:677-694.

    [11]LIN Z A,JI Y F,LIANG W B,et al.Landslide displacement prediction based on time-frequency analysis and LMD-BiLSTM model [J].Mathematics,2022,10(13):22032227.

    [12]WANG Y H,YEH C H,YOUNG V H W,et al.On the computational complexity of the empirical mode decomposition algorithm[J].Physica A Statistical Mechanics amp; Its Applications,2014,400:159-167.

    [13]LIU H,TIAN H Q,LI Y F.Comparison of new hybrid FEEMD-MLP,F(xiàn)EEMD-ANFIS,Wavelet Packet-MLP and Wavelet Packet-ANFIS for wind speed preditions[J].Energy Conversion and Management,2015,89:1-11.

    [14]饒鴻江,何永泰,彭潔,等.基于FEEMD-NTEO的風(fēng)電場(chǎng)送出線路故障定位[J].電機(jī)與控制應(yīng)用,2024,51(1):106-116.

    [15]殷禮勝,魏帥康,孫雙晨,等.基于FEEMD-SAPSO-BiLSTM組合模型的短時(shí)交通流預(yù)測(cè)[J].電子測(cè)量與儀器學(xué)報(bào),2021,35(10):72-81.

    [16]蒲嫻怡,畢貴紅,王凱,等.基于FEEMD-PACF-BP_AdaBoost模型的風(fēng)電功率超短期預(yù)測(cè)[J].計(jì)算機(jī)應(yīng)用與軟件,2021,38(11):91-97.

    [17]劉若男,辛義忠,李巖.基于皮爾遜相關(guān)系數(shù)的動(dòng)態(tài)簽名驗(yàn)證方法[J].儀器儀表學(xué)報(bào),2022,43(7):279-287.

    [18]FU T,TANG X B,CAI Z K,et al.Correlation research of phase angle variation and coating performance by means of Pearson′s correlation coefficient[J].Progress in Organic Coatings,2020,139:1005459.

    [19]趙兵,王增平,紀(jì)維佳,等.基于注意力機(jī)制的 CNN-GRU 短期電力負(fù)荷預(yù)測(cè)方法[J].電網(wǎng)技術(shù),2019,43(12):4370-4376.

    [20]鄭宏舟,趙宇宸,孟飛.基于深度學(xué)習(xí)GRU網(wǎng)絡(luò)的UWB室內(nèi)定位優(yōu)化[J].上海理工大學(xué)學(xué)報(bào),2024,46(1):78-86.

    [21]韋小多,陳艷,敖煜測(cè),等.一種基于STGCN-GRU模型的交通流量預(yù)測(cè)[J].中國(guó)科技信息,2023(23):114-116.

    [22]趙海峰,諸立凱,劉長(zhǎng)松,等.基于注意力機(jī)制的CNN-GRU煤層氣產(chǎn)能預(yù)測(cè)方法研究[J].煤礦安全,2023,54(12):11-17.

    [23]FRANCIS C,HYODONG H,OOK L,et al.A novel hybrid deep learning approach to code generation aimed at mitigating the real-time network attack in the mobile experiment Via GRU-LM and Word2vec[J].Mobile Information Systems,2022,2022(1):1-11.

    [24]薛朝輝,周逸飏,強(qiáng)永剛,等.融合 NetVLAD 和全連接層的三元神經(jīng)網(wǎng)絡(luò)交叉視角場(chǎng)景圖像定位[J].遙感學(xué)報(bào),2021,25(5):1095-1107.

    [25]李清濤,任宇馳,王遠(yuǎn),等.基于人工神經(jīng)網(wǎng)絡(luò)全連接層優(yōu)化的線損異常診斷方法研究[J].電氣應(yīng)用,2020,39(4):82-88.

    (編輯:劉 媛)

    Landslide displacement prediction based on FEEMD-GRU-FC model

    HE Qing1,LI Lilin1,LIN Zian2

    (1.Science and Technology Information Department,Guangxi Transportation Science and Technology Group Co.,Ltd.,Nanning 530000,China;" 2.School of Computer and Information Security,Guilin University of Electronic Technology,Guilin 541000,China)

    Abstract:

    Landslides in reservoir areas are affected by rainfall and periodic changes in reservoir water level,which will interfere with the prediction performance of landslide displacement prediction models.A FEEMD-GRU-FC (fast ensemble empirical mode decomposition-gated recurrent unit-fully connected layer) landslide displacement prediction model based on the periodic environmental changes was proposed.The model fully considered the influence of noise among different frequencies of displacement data,divided the displacement into sub-sequences of different frequencies to predict them respectively,and used the PCC (Pearson correlation coefficient) method to select external environmental factors with large influence.The ability of the GRU model to extract periodic change features from environmental data was improved by adding a full connection layer.The model was applied by combining with the actual data of the Bazimen Landslide in the Three Gorges Reservoir area.Different prediction models and evaluation indicators were used for multi-perspective comparison and verification.The experimental results showed that the model proposed was better than the other three models in predicting Bazimen Landslide displacement.Its fitting degree was 98.6%,and three kinds of error indexes were the lowest,which proved that the model can accurately predict the landslide displacement in the reservoir area.The research results can guide disaster prevention and reduction.

    Key words:

    landslide displacement prediction; Pearson correlation coefficient; gated recurrent unit; fully connected layer; nonlinearity; Bazimen Landslide; Three Gorges Reservoir area

    猜你喜歡
    八字降雨滑坡
    短短八字,道盡一生
    意林(2023年18期)2023-11-05 10:45:10
    滑坡推力隱式解與顯式解對(duì)比分析——以河北某膨脹土滑坡為例
    5類人跟泡腳“八字不合”
    滄州市2016年“7.19~7.22”與“8.24~8.25”降雨對(duì)比研究
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    基于Fluent的滑坡入水過(guò)程數(shù)值模擬
    紅黏土降雨入滲的定量分析
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    南方降雨不斷主因厄爾尼諾
    江埡水庫(kù)降雨徑流相關(guān)圖的建立
    亚洲国产日韩一区二区| 精品一区二区三区视频在线| 久久国内精品自在自线图片| av卡一久久| 久久午夜福利片| 国产 精品1| 一边亲一边摸免费视频| 亚洲精品日韩av片在线观看| 精品久久久久久久久亚洲| 欧美日韩精品成人综合77777| 亚洲综合色惰| 国产成人91sexporn| 国产欧美另类精品又又久久亚洲欧美| 国产精品一及| 一边亲一边摸免费视频| 少妇裸体淫交视频免费看高清| 一区二区三区乱码不卡18| 精品久久久精品久久久| 中文乱码字字幕精品一区二区三区| 国产中年淑女户外野战色| 亚洲欧美日韩另类电影网站 | 亚洲第一av免费看| 国产免费又黄又爽又色| 欧美成人午夜免费资源| 久久国产精品男人的天堂亚洲 | 精品视频人人做人人爽| 久久久a久久爽久久v久久| 直男gayav资源| 国产欧美日韩一区二区三区在线 | 97在线视频观看| 成人影院久久| 国产高清国产精品国产三级 | 久久久久精品久久久久真实原创| 亚洲激情五月婷婷啪啪| 国产午夜精品久久久久久一区二区三区| 国产黄色免费在线视频| 国产极品天堂在线| 亚洲激情五月婷婷啪啪| 少妇人妻一区二区三区视频| 欧美精品一区二区免费开放| 久久99蜜桃精品久久| 日韩中字成人| 一级毛片 在线播放| 欧美精品国产亚洲| 亚洲国产欧美人成| 狂野欧美激情性xxxx在线观看| 在线观看国产h片| 免费观看性生交大片5| 91精品伊人久久大香线蕉| 精品久久久精品久久久| 国产人妻一区二区三区在| 婷婷色综合www| 日本黄色日本黄色录像| 性高湖久久久久久久久免费观看| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 高清黄色对白视频在线免费看 | 高清黄色对白视频在线免费看 | 午夜日本视频在线| 国产精品国产三级专区第一集| 男女下面进入的视频免费午夜| 日本wwww免费看| 一本—道久久a久久精品蜜桃钙片| 色婷婷av一区二区三区视频| 777米奇影视久久| 成人特级av手机在线观看| 欧美精品国产亚洲| 国产成人免费观看mmmm| 久久久国产一区二区| 视频中文字幕在线观看| 视频中文字幕在线观看| 久久 成人 亚洲| 精品久久久精品久久久| 99久国产av精品国产电影| 深夜a级毛片| 国产黄片美女视频| 国产成人精品福利久久| 亚洲成人av在线免费| 午夜福利高清视频| 亚洲av中文av极速乱| 哪个播放器可以免费观看大片| 男女下面进入的视频免费午夜| 国产精品秋霞免费鲁丝片| 亚洲精品乱久久久久久| 亚洲内射少妇av| 啦啦啦视频在线资源免费观看| 国产精品熟女久久久久浪| 国产精品国产三级国产专区5o| 精品人妻一区二区三区麻豆| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| kizo精华| 女人久久www免费人成看片| 免费观看av网站的网址| 男人舔奶头视频| 99热网站在线观看| 久久久精品94久久精品| 午夜福利视频精品| 热99国产精品久久久久久7| 亚洲成人一二三区av| 乱码一卡2卡4卡精品| 日韩在线高清观看一区二区三区| 色综合色国产| 国产爱豆传媒在线观看| 亚洲欧美精品专区久久| 精品国产乱码久久久久久小说| 久久99热这里只有精品18| 精品久久久久久久久av| 在线观看国产h片| 一级毛片我不卡| 国产中年淑女户外野战色| 久久精品国产a三级三级三级| 国产色婷婷99| 久久久久久久久久久免费av| 亚洲一区二区三区欧美精品| 亚洲怡红院男人天堂| 多毛熟女@视频| 王馨瑶露胸无遮挡在线观看| 日韩 亚洲 欧美在线| 99久久人妻综合| 联通29元200g的流量卡| 黄色配什么色好看| 91久久精品电影网| www.色视频.com| 久久青草综合色| 在现免费观看毛片| 中文字幕亚洲精品专区| 免费黄频网站在线观看国产| 国产在线视频一区二区| 99精国产麻豆久久婷婷| 精品国产露脸久久av麻豆| 国产精品福利在线免费观看| 亚洲自偷自拍三级| 91在线精品国自产拍蜜月| 国产人妻一区二区三区在| 欧美极品一区二区三区四区| 国产精品麻豆人妻色哟哟久久| 熟女av电影| 国产在线免费精品| 交换朋友夫妻互换小说| 亚洲美女视频黄频| 国产亚洲午夜精品一区二区久久| 欧美成人a在线观看| 亚洲av免费高清在线观看| 亚洲精品456在线播放app| 联通29元200g的流量卡| 国产欧美亚洲国产| 亚洲欧美一区二区三区黑人 | 十分钟在线观看高清视频www | 精品人妻视频免费看| 国产一区亚洲一区在线观看| 国产综合精华液| 国产精品久久久久成人av| 中文字幕制服av| av视频免费观看在线观看| 肉色欧美久久久久久久蜜桃| 最黄视频免费看| 少妇的逼好多水| 亚洲色图综合在线观看| 久久久久国产精品人妻一区二区| 丝瓜视频免费看黄片| 人人妻人人看人人澡| 国产熟女欧美一区二区| 高清黄色对白视频在线免费看 | 亚洲熟女精品中文字幕| 国产精品人妻久久久影院| 啦啦啦中文免费视频观看日本| 国产 精品1| 国产日韩欧美亚洲二区| 成年人午夜在线观看视频| 国产视频内射| 涩涩av久久男人的天堂| 国产探花极品一区二区| 婷婷色麻豆天堂久久| 日韩大片免费观看网站| 午夜福利视频精品| 热re99久久精品国产66热6| 欧美一级a爱片免费观看看| 国产白丝娇喘喷水9色精品| 午夜福利在线在线| 国产精品免费大片| 亚洲婷婷狠狠爱综合网| 久久久久性生活片| 少妇的逼水好多| 三级经典国产精品| 蜜臀久久99精品久久宅男| 777米奇影视久久| 高清欧美精品videossex| av线在线观看网站| 国产精品99久久久久久久久| 91午夜精品亚洲一区二区三区| 亚洲精品色激情综合| 国产在视频线精品| 99热这里只有精品一区| 免费看光身美女| 亚洲美女黄色视频免费看| 日本爱情动作片www.在线观看| 少妇熟女欧美另类| 亚洲av综合色区一区| 中文在线观看免费www的网站| 国产一区二区三区综合在线观看 | 日日撸夜夜添| 国产高清有码在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产精品麻豆人妻色哟哟久久| av天堂中文字幕网| 婷婷色综合大香蕉| 国产黄片视频在线免费观看| 啦啦啦啦在线视频资源| 一个人看视频在线观看www免费| 久久99蜜桃精品久久| 狂野欧美激情性xxxx在线观看| 麻豆成人午夜福利视频| 丝袜喷水一区| 国产高清不卡午夜福利| 免费高清在线观看视频在线观看| 有码 亚洲区| 久久国内精品自在自线图片| 久久精品国产a三级三级三级| 久久精品国产亚洲av涩爱| 午夜福利在线在线| 亚洲国产欧美人成| 在线免费十八禁| av在线app专区| 看十八女毛片水多多多| 婷婷色综合大香蕉| 校园人妻丝袜中文字幕| 人妻系列 视频| 成年女人在线观看亚洲视频| 免费观看av网站的网址| 一本色道久久久久久精品综合| 午夜福利视频精品| 日日啪夜夜爽| 国产精品一区二区在线观看99| 免费播放大片免费观看视频在线观看| 偷拍熟女少妇极品色| 国内精品宾馆在线| 精品亚洲乱码少妇综合久久| 亚洲经典国产精华液单| 日韩精品有码人妻一区| 亚洲精品久久午夜乱码| 国产一级毛片在线| 日日摸夜夜添夜夜添av毛片| 国内少妇人妻偷人精品xxx网站| 日本猛色少妇xxxxx猛交久久| 亚洲国产最新在线播放| 久久久久久久久久久丰满| 欧美人与善性xxx| av国产免费在线观看| 国精品久久久久久国模美| 中文字幕制服av| 久久久久久久久久人人人人人人| 亚洲天堂av无毛| 综合色丁香网| 美女视频免费永久观看网站| 免费看光身美女| 国产真实伦视频高清在线观看| 国产精品精品国产色婷婷| 久久久久国产精品人妻一区二区| 亚洲av中文字字幕乱码综合| 男女啪啪激烈高潮av片| 一级毛片电影观看| 欧美精品国产亚洲| 亚洲va在线va天堂va国产| 亚洲不卡免费看| 黄色日韩在线| 99久久精品热视频| 最近手机中文字幕大全| 日韩视频在线欧美| 天天躁夜夜躁狠狠久久av| 国产国拍精品亚洲av在线观看| 亚洲av免费高清在线观看| av专区在线播放| 男女边摸边吃奶| 国产精品一区二区在线不卡| 亚洲av综合色区一区| 少妇高潮的动态图| 少妇被粗大猛烈的视频| 久久午夜福利片| 最近最新中文字幕大全电影3| 精品熟女少妇av免费看| 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说 | 成人无遮挡网站| 丝瓜视频免费看黄片| 国产av码专区亚洲av| av免费在线看不卡| 婷婷色综合www| 中文字幕久久专区| 国产有黄有色有爽视频| 午夜激情久久久久久久| 最近2019中文字幕mv第一页| 国产又色又爽无遮挡免| 日本与韩国留学比较| 啦啦啦中文免费视频观看日本| 欧美高清性xxxxhd video| 岛国毛片在线播放| 草草在线视频免费看| 高清欧美精品videossex| 精品视频人人做人人爽| 91午夜精品亚洲一区二区三区| 亚洲国产精品一区三区| 三级国产精品片| 国产又色又爽无遮挡免| 国产黄色免费在线视频| 伦理电影免费视频| 下体分泌物呈黄色| 精品一区二区三卡| 亚洲精品日韩在线中文字幕| 91久久精品国产一区二区三区| 狂野欧美激情性bbbbbb| 国产午夜精品一二区理论片| 国产亚洲5aaaaa淫片| 免费大片18禁| 久久久久久人妻| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍三级| 国产在视频线精品| 国产 一区 欧美 日韩| 精品酒店卫生间| 少妇人妻久久综合中文| 美女内射精品一级片tv| 亚洲精品乱码久久久v下载方式| 国产男女超爽视频在线观看| 亚洲av成人精品一二三区| 大片免费播放器 马上看| 一级av片app| 久久久久久久久久久免费av| 九九久久精品国产亚洲av麻豆| 久久6这里有精品| 晚上一个人看的免费电影| 老司机影院毛片| 99热国产这里只有精品6| 欧美成人a在线观看| 国产久久久一区二区三区| 精品国产一区二区三区久久久樱花 | .国产精品久久| 99久久精品热视频| 日韩免费高清中文字幕av| 午夜视频国产福利| 在线观看av片永久免费下载| 国国产精品蜜臀av免费| 亚洲人与动物交配视频| 久久毛片免费看一区二区三区| 日本av免费视频播放| 国产成人a∨麻豆精品| 国产一区有黄有色的免费视频| 久久女婷五月综合色啪小说| 十分钟在线观看高清视频www | 寂寞人妻少妇视频99o| 日韩强制内射视频| 91精品国产九色| 免费大片黄手机在线观看| 国产精品福利在线免费观看| 91久久精品国产一区二区成人| 又爽又黄a免费视频| 边亲边吃奶的免费视频| 天堂8中文在线网| 国产乱人偷精品视频| 欧美三级亚洲精品| 少妇裸体淫交视频免费看高清| 欧美成人午夜免费资源| 两个人的视频大全免费| 啦啦啦中文免费视频观看日本| 中文在线观看免费www的网站| 午夜福利在线在线| 男男h啪啪无遮挡| 国产大屁股一区二区在线视频| 久久久欧美国产精品| 久久99蜜桃精品久久| 插阴视频在线观看视频| 久久久欧美国产精品| 亚洲伊人久久精品综合| 久久久久网色| 国产欧美亚洲国产| 一级a做视频免费观看| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 成人18禁高潮啪啪吃奶动态图 | av在线蜜桃| 精品少妇久久久久久888优播| 国产深夜福利视频在线观看| 午夜福利高清视频| 日本av手机在线免费观看| 欧美97在线视频| 久久精品国产自在天天线| 乱系列少妇在线播放| 亚洲在久久综合| 免费久久久久久久精品成人欧美视频 | 亚洲av在线观看美女高潮| 在线看a的网站| 国产精品爽爽va在线观看网站| 蜜桃久久精品国产亚洲av| 人妻夜夜爽99麻豆av| 校园人妻丝袜中文字幕| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| 我要看黄色一级片免费的| 国产伦在线观看视频一区| 97超视频在线观看视频| 成人二区视频| 最近中文字幕2019免费版| 丰满迷人的少妇在线观看| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 我要看黄色一级片免费的| 久久鲁丝午夜福利片| 黑人高潮一二区| 免费看日本二区| 国产免费一级a男人的天堂| 久久青草综合色| 久久这里有精品视频免费| 国产精品av视频在线免费观看| 一级黄片播放器| 一级毛片黄色毛片免费观看视频| 天美传媒精品一区二区| 精品视频人人做人人爽| 久久国产乱子免费精品| 在线免费十八禁| 三级经典国产精品| 国产高清国产精品国产三级 | 亚洲美女视频黄频| 免费在线观看成人毛片| 黑丝袜美女国产一区| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 又大又黄又爽视频免费| 国产精品99久久99久久久不卡 | 97超视频在线观看视频| 国产欧美日韩一区二区三区在线 | 国产亚洲av片在线观看秒播厂| 亚洲成人中文字幕在线播放| 国产视频内射| 国产乱人视频| 在线 av 中文字幕| 老司机影院毛片| 日本av手机在线免费观看| 午夜免费观看性视频| 成人国产av品久久久| 青春草国产在线视频| 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| 亚洲精品国产av蜜桃| 日日啪夜夜撸| 我要看日韩黄色一级片| 成年女人在线观看亚洲视频| av国产精品久久久久影院| 极品少妇高潮喷水抽搐| 国产有黄有色有爽视频| 天天躁日日操中文字幕| 日本欧美视频一区| 哪个播放器可以免费观看大片| 黄色怎么调成土黄色| 嫩草影院入口| 一级毛片aaaaaa免费看小| 日韩一本色道免费dvd| 99久久综合免费| 久久av网站| 丰满乱子伦码专区| 男人添女人高潮全过程视频| 国产精品国产三级专区第一集| 天堂8中文在线网| 赤兔流量卡办理| 哪个播放器可以免费观看大片| 午夜福利视频精品| 人妻一区二区av| 亚洲国产欧美人成| 一级爰片在线观看| 精品一区二区免费观看| 欧美成人精品欧美一级黄| 欧美亚洲 丝袜 人妻 在线| 欧美精品亚洲一区二区| av.在线天堂| videossex国产| 十八禁网站网址无遮挡 | 男女啪啪激烈高潮av片| 色婷婷久久久亚洲欧美| 国产精品.久久久| 观看美女的网站| 91精品一卡2卡3卡4卡| 黑人猛操日本美女一级片| 久久久a久久爽久久v久久| 少妇人妻久久综合中文| 国产爽快片一区二区三区| 国产一区二区在线观看日韩| 日本黄色片子视频| 日韩av不卡免费在线播放| 亚洲av免费高清在线观看| 99视频精品全部免费 在线| 免费观看性生交大片5| 亚洲精品一二三| 亚洲精品国产成人久久av| 久久国产精品男人的天堂亚洲 | 日日摸夜夜添夜夜爱| 中国美白少妇内射xxxbb| 一级毛片电影观看| a级一级毛片免费在线观看| 2018国产大陆天天弄谢| 成人美女网站在线观看视频| 欧美另类一区| 精品人妻一区二区三区麻豆| 精华霜和精华液先用哪个| 热re99久久精品国产66热6| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 免费观看在线日韩| av不卡在线播放| 99精国产麻豆久久婷婷| 免费观看无遮挡的男女| 亚洲久久久国产精品| 我要看日韩黄色一级片| 人人妻人人看人人澡| 一边亲一边摸免费视频| 亚洲av成人精品一区久久| 日韩中文字幕视频在线看片 | 日本色播在线视频| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 成人特级av手机在线观看| 成年免费大片在线观看| 久久这里有精品视频免费| 男的添女的下面高潮视频| 国产成人精品婷婷| 国产亚洲精品久久久com| 国产久久久一区二区三区| 岛国毛片在线播放| 成人黄色视频免费在线看| 国国产精品蜜臀av免费| 亚洲国产高清在线一区二区三| 欧美成人a在线观看| 欧美最新免费一区二区三区| 女人十人毛片免费观看3o分钟| 高清视频免费观看一区二区| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 一级毛片久久久久久久久女| 99re6热这里在线精品视频| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 亚洲色图综合在线观看| 欧美高清成人免费视频www| 国产视频首页在线观看| 人妻 亚洲 视频| 久久久精品94久久精品| 国产成人一区二区在线| 国精品久久久久久国模美| 亚洲成人av在线免费| 高清不卡的av网站| 一个人看的www免费观看视频| 亚洲国产最新在线播放| 国产av码专区亚洲av| 亚洲精品国产av成人精品| 国产成人精品久久久久久| a 毛片基地| 91午夜精品亚洲一区二区三区| 久久久久国产精品人妻一区二区| 日日啪夜夜撸| 免费看日本二区| 又大又黄又爽视频免费| 99热这里只有精品一区| 国产成人免费观看mmmm| 涩涩av久久男人的天堂| 午夜激情久久久久久久| 久久精品国产鲁丝片午夜精品| 高清午夜精品一区二区三区| 国产片特级美女逼逼视频| 亚洲精品乱码久久久久久按摩| 韩国高清视频一区二区三区| 亚洲精品日韩av片在线观看| 免费看日本二区| 美女主播在线视频| 亚洲精品乱码久久久v下载方式| 爱豆传媒免费全集在线观看| 一个人免费看片子| 国产精品三级大全| 欧美高清成人免费视频www| 天堂8中文在线网| 亚洲中文av在线| 丰满少妇做爰视频| 亚洲av日韩在线播放| 在线免费观看不下载黄p国产| 国产精品偷伦视频观看了| 国产人妻一区二区三区在| 久久精品人妻少妇| 国产精品国产三级国产专区5o| 美女主播在线视频| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 日本欧美视频一区| 国产精品爽爽va在线观看网站| 在线精品无人区一区二区三 | 日韩在线高清观看一区二区三区| 午夜免费鲁丝| av天堂中文字幕网| 高清欧美精品videossex| 久久精品人妻少妇| 建设人人有责人人尽责人人享有的 | 这个男人来自地球电影免费观看 | 在线播放无遮挡| av在线播放精品| 久久女婷五月综合色啪小说| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 麻豆精品久久久久久蜜桃| 日韩伦理黄色片| 国产成人精品一,二区| 伦理电影免费视频| 国产一区二区三区综合在线观看 | 色视频www国产| 亚洲欧洲国产日韩| 少妇精品久久久久久久| 国产男女内射视频| 不卡视频在线观看欧美|