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

    基于LMBP算法的崇陽(yáng)溪流域降雨徑流預(yù)報(bào)模型研究

    2019-12-21 02:46:42金保明杜倫閱李斐斐
    關(guān)鍵詞:崇陽(yáng)雨量站武夷山

    金保明,王 偉,杜倫閱,李斐斐

    (福州大學(xué)土木工程學(xué)院,福建 福州 350108)

    0 引言

    我國(guó)山川眾多,洪澇災(zāi)害頻發(fā),給人民群眾生命財(cái)產(chǎn)帶來(lái)?yè)p失,影響國(guó)民經(jīng)濟(jì)的發(fā)展.隨著國(guó)內(nèi)防洪工程的建設(shè)和非工程措施系統(tǒng)的日益完善,國(guó)內(nèi)大江大河的防洪能力得到了極大的改善,能夠抵擋較大洪水,遇特大洪水也有一定的應(yīng)急措施[1].但是山區(qū)中小流域由于地形原因,流域降雨強(qiáng)度大,降雨集中,河流洪水具有暴漲暴落的特點(diǎn),往往會(huì)導(dǎo)致洪水災(zāi)害,難于預(yù)報(bào)和預(yù)防,給當(dāng)?shù)貛?lái)巨大的損失[2].洪水預(yù)報(bào)是重要的防洪非工程措施,目前主要水文預(yù)報(bào)模型包括集總式和分布式水文模型.集總式水文模型由于結(jié)構(gòu)簡(jiǎn)單,應(yīng)用廣泛,但是模型沒(méi)有考慮流域內(nèi)的氣候變化、土地利用、植被覆蓋、地形等下墊面因素的變化,是一種概化的數(shù)學(xué)模型.分布式水文模型則考慮了上述因子的變化,能較好模擬實(shí)際的洪水過(guò)程,但它需要大量的水文信息和空間數(shù)據(jù).山區(qū)中小流域一般實(shí)測(cè)資料系列相對(duì)較短,空間數(shù)據(jù)不易獲得,所以該類(lèi)模型在山區(qū)中小流域的應(yīng)用并不能得到很好的預(yù)報(bào)效果,預(yù)報(bào)誤差相對(duì)較大[3].相比較而言,人工神經(jīng)網(wǎng)絡(luò)模型具有很強(qiáng)的非線性映射能力,它不需要空間地形資料、只需要水文資料,可以應(yīng)用于山區(qū)中小流域水文預(yù)報(bào).

    目前大約80%的人工神經(jīng)網(wǎng)絡(luò)模型采取了BP網(wǎng)絡(luò)或BP網(wǎng)絡(luò)的變化形式[4].王建金等[5]選取海南省定安河流域作為研究區(qū)域,采用改進(jìn)的BP神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)報(bào),結(jié)果表明預(yù)報(bào)精度較高.Veintimilla-Reyes等[6]將BP神經(jīng)網(wǎng)絡(luò)輸出層與隱含層權(quán)重更新方法進(jìn)行優(yōu)化,對(duì)位于厄瓜多爾的托梅班巴河小流域的降雨徑流進(jìn)行預(yù)報(bào),結(jié)果表明,通過(guò)優(yōu)化明顯提高了預(yù)報(bào)精度.Raj等[7]運(yùn)用小波算法的BP神經(jīng)網(wǎng)絡(luò)對(duì)印度泰米爾納德邦的卡尼亞庫(kù)馬里地區(qū)的短期降雨量進(jìn)行預(yù)測(cè),預(yù)測(cè)效果良好.傳統(tǒng)的BP神經(jīng)網(wǎng)絡(luò)模型,由于采用最速下降法[8],存在收斂速度慢[9]、容易陷入局部極小值點(diǎn)[10]、隱含層節(jié)點(diǎn)數(shù)難以確定及無(wú)法研究物理機(jī)理[11]等問(wèn)題,因此目前一般采用改進(jìn)的算法.Levenberg-Marquardt算法就是其中的一種改進(jìn)算法.閔祥宇[12]提出了改進(jìn)的GA-LMBP模型,使模型能夠跳出局部極小值,同時(shí)網(wǎng)絡(luò)加速收斂,減少了訓(xùn)練時(shí)間.張素瓊等[13]建立LMBP神經(jīng)網(wǎng)絡(luò)模型,對(duì)月徑流進(jìn)行預(yù)報(bào),預(yù)報(bào)結(jié)果與實(shí)測(cè)值十分接近.劉佩瑤等[14]采用新安江模型和LMBP神經(jīng)網(wǎng)絡(luò)模型對(duì)富屯溪流域洪水進(jìn)行預(yù)報(bào),經(jīng)過(guò)比較,表明LMBP神經(jīng)網(wǎng)絡(luò)模型的精度不僅略優(yōu)于新安江模型,而且相比于傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)具有收斂速度快、時(shí)效性好的優(yōu)點(diǎn).

    因此利用LMBP神經(jīng)網(wǎng)絡(luò)對(duì)山區(qū)中小流域進(jìn)行洪水預(yù)報(bào)是可行的.本研究根據(jù)崇陽(yáng)溪武夷山水文站的前期流量資料以及上游雨量站的時(shí)段雨量資料,建立LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流預(yù)報(bào)模型,進(jìn)行洪水流量預(yù)報(bào).

    1 LMBP算法

    當(dāng)LMBP算法用于多層神經(jīng)網(wǎng)絡(luò)訓(xùn)練問(wèn)題時(shí),其性能優(yōu)化指數(shù)F(x)[15]可表示為

    (1)

    其中:x是權(quán)值和偏置值組成的向量;tq是實(shí)測(cè)輸出;aq是網(wǎng)絡(luò)輸出;v(x)是誤差向量;Q是目標(biāo)輸出的個(gè)數(shù).

    其梯度函數(shù)[15]為

    ▽F(x)=2JT(x)v(x)

    (2)

    其中:J(x)為雅可比矩陣.

    進(jìn)而推導(dǎo)出Levenberg-Marquardt算法的公式[15]為

    xk+1=xk-[JT(xk)J(xk)+μkI]-1JT(xk)v(xk)

    (3)

    其中:J(xk)為雅可比矩陣;v(xk)為誤差向量;I為單位矩陣;μk為一個(gè)大于零的系數(shù)(開(kāi)始計(jì)算時(shí)取μk=0.01);k為算法迭代次數(shù).

    LMBP算法的計(jì)算過(guò)程如下.

    1)將所有輸入數(shù)據(jù)提交網(wǎng)絡(luò),數(shù)據(jù)逐層向前傳播,計(jì)算到網(wǎng)絡(luò)最后一層的輸出值aq.用式(1)計(jì)算所有輸入的平方誤差之和F(x).

    2)計(jì)算網(wǎng)絡(luò)的Marquardt敏感度并進(jìn)行反向傳播,計(jì)算出雅可比矩陣J(x).

    3)通過(guò)式(3)求解Δxk.

    4)用xk+Δxk重復(fù)計(jì)算平方誤差之和F(x).如果F(x)的值減小,則用μk除以一個(gè)大于1的正數(shù)θ,然后更新xk+1轉(zhuǎn)至第1)步;如果F(x)的值未減小,則用μk乘以一個(gè)大于1的正數(shù)θ[15],不更新xk+1并轉(zhuǎn)至第3)步.

    如此反復(fù)訓(xùn)練,當(dāng)平方誤差和F(x)小于某一給定值時(shí)(本次取0.001),算法被認(rèn)為收斂,模型參數(shù)計(jì)算完成.

    2 崇陽(yáng)溪流域LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流預(yù)報(bào)模型

    2.1 流域單元面積劃分

    崇陽(yáng)溪上游流域位于福建省西北部的武夷山市,武夷山站以上流域面積為1 078 km2,流域水系圖如圖1所示.在ArcGIS平臺(tái)上,利用泰森多邊形法將流域劃分為6個(gè)子流域,求出洋莊、吳邊、大安、坑口、嶺陽(yáng)、嵐谷雨量站控制子流域的面積分別為:129、239、191、161、176、182 km2,相應(yīng)的面積權(quán)重分別為:0.12、0.22、0.18、0.15、0.16、0.17.圖2為崇陽(yáng)溪上游流域泰森多邊形劃分圖.通過(guò)資料分析,確定武夷山水文站上游6個(gè)雨量站凈雨形成的徑流到武夷山水文站的匯流時(shí)間分別為1、1.5、2.5、3、4、3 h.

    圖1 崇陽(yáng)溪上游流域水系圖Fig.1 River system of the upper Chongyang River

    圖2 崇陽(yáng)溪上游流域泰森多邊形劃分圖Fig.2 Thiessen polygon in upper reaches ofChongyang River

    2.2 模型建立

    根據(jù)上述6個(gè)雨量站的時(shí)段雨量資料和武夷山水文站的前期流量資料,建立數(shù)學(xué)模型

    Q(t)=F[μiPi(t-ti),Q(t-1)]+ε(t)

    (4)

    式中:Q(t)為預(yù)報(bào)流量(m3·s-1);μi為單元子流域面積權(quán)重,i=1、2、3、4、5、6;Pi為第i個(gè)雨量站的時(shí)段降雨量(mm),i=1、2、3、4、5、6;Q(t-1)為武夷山水文站前1 h流量(m3·s-1);ε(t)為模型誤差(m3·s-1);ti為第i個(gè)雨量站凈雨形成的徑流到武夷山水文站匯流時(shí)間(h).

    將各單元子流域面積權(quán)值和各雨量站洪水傳播時(shí)間代入公式(4)可將數(shù)學(xué)模型具體化

    (5)

    式中:P1、P2、P3、P4、P5、P6分別為洋莊、吳邊、大安、坑口、嶺陽(yáng)、嵐谷站的時(shí)段雨量;Q(t-1)、Q(t)的含義同式(4).

    模型采用3層網(wǎng)絡(luò)結(jié)構(gòu),輸入層與隱含層的傳遞函數(shù)采用Logsig函數(shù),隱含層到輸出層采用線性函數(shù).輸入層包括式(5)中的7個(gè)節(jié)點(diǎn);中間為隱含層,隱含層的節(jié)點(diǎn)個(gè)數(shù)采用試算法最終確定;輸出層為一個(gè)節(jié)點(diǎn),即武夷山水文站的洪水流量Q(t).

    2.3 樣本數(shù)據(jù)的生成和處理

    選取1997—2014年崇陽(yáng)溪上游流域的21場(chǎng)洪水過(guò)程,將其中14場(chǎng)洪水作為訓(xùn)練樣本,7場(chǎng)洪水作為檢驗(yàn)樣本,用于模型的評(píng)價(jià)和最終的檢驗(yàn).

    模型中輸入流量的量級(jí)在[102,103]范圍內(nèi),降雨量的量級(jí)在[0,101]范圍內(nèi).為了使所有元素在神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程中受到同等重視,將輸入數(shù)據(jù)利用下式進(jìn)行歸一化處理,使數(shù)據(jù)變換到[-1,1]范圍內(nèi).

    (6)

    式中:Ym為歸一化處理后的輸入值;Xm為實(shí)測(cè)降雨輸入(mm)或流量輸入(m3·s-1);Xmax為最大實(shí)測(cè)降雨量(mm)或流量值(m3·s-1);Xmin為最小實(shí)測(cè)降雨量(mm)或流量值(m3·s-1).

    網(wǎng)絡(luò)最后的輸出結(jié)果需要經(jīng)過(guò)反歸一化處理才能得到實(shí)際的輸出流量,如下式所示

    (7)

    式中:X為反歸一化后的輸出流量(m3·s-1);Y為L(zhǎng)MBP神經(jīng)網(wǎng)絡(luò)輸出值;Xmax為最大實(shí)測(cè)流量值(m3·s-1);Xmin為最小實(shí)測(cè)流量值(m3·s-1).

    2.4 模型參數(shù)的率定

    模型輸入層單元數(shù)為7,輸出層單元數(shù)為1,而隱含層單元數(shù)的確定,則利用試算法,分別選取n=5、6、7、8、9、10、11單元數(shù)進(jìn)行試算,通過(guò)程序進(jìn)行網(wǎng)絡(luò)訓(xùn)練,最終對(duì)比分析確定最優(yōu)的單元節(jié)點(diǎn)數(shù)為8,得到模型各層權(quán)值和閥值如表1所示.

    表1 LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流預(yù)報(bào)模型各層權(quán)值、閥值Tab.1 LMBP neural network runoff prediction model weight and threshold value of each layer

    2.5 模型預(yù)報(bào)結(jié)果

    選擇“1998·06·14”、“1999·06·17”、“2003·06·25”、“2005·06·19”、“2007·06·14”、“2008·07·19”、“2013·05·27”7場(chǎng)洪水流量過(guò)程作為檢驗(yàn)樣本,對(duì)LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流預(yù)報(bào)模型進(jìn)行檢驗(yàn),得到7場(chǎng)洪水預(yù)報(bào)過(guò)程,詳見(jiàn)圖3~9.

    圖3 武夷山站“1998·06·14”洪水流量預(yù)報(bào)過(guò)程線圖Fig.3 Flood flow hydrograph of Wuyishan hydrologicstation on June 14th,1998

    圖4 武夷山站“1999·06·17”洪水流量預(yù)報(bào)過(guò)程線圖Fig.4 Flood flow hydrograph of Wuyishan hydrologicstation on June 17th,1999

    圖5 武夷山站“2003·06·25”洪水流量預(yù)報(bào)過(guò)程線圖Fig.5 Flood flow hydrograph of Wuyishan hydrologicstation on June 25th,2003

    圖6 武夷山站“2005·06·19”洪水流量預(yù)報(bào)過(guò)程線圖Fig.6 Flood flow hydrograph of Wuyishan hydrologicstation on June 19th,2005

    圖7 武夷山站“2007·06·14”洪水流量預(yù)報(bào)過(guò)程線圖Fig.7 Flood flow hydrograph of Wuyishan hydrologicstation on June 14th,2007

    圖8 武夷山站“2008·07·19”洪水流量預(yù)報(bào)過(guò)程線圖Fig.8 Flood flow hydrograph of Wuyishan hydrologic station on July 19th,2008

    圖9 武夷山站“2013·05·27”洪水流量預(yù)報(bào)過(guò)程線圖Fig.9 Flood flow hydrograph of Wuyishan hydrologicstation on May 27th,2013

    對(duì)每場(chǎng)洪水過(guò)程的預(yù)報(bào)結(jié)果進(jìn)行誤差分析,采用絕對(duì)誤差最大值、相對(duì)誤差最大值、平均誤差、確定性系數(shù)等指標(biāo)來(lái)評(píng)價(jià)LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流預(yù)報(bào)模型的精度,列于表2.同時(shí),對(duì)每場(chǎng)洪水的洪峰流量誤差進(jìn)行統(tǒng)計(jì)分析,列于表3.

    表2 7場(chǎng)洪水流量過(guò)程預(yù)報(bào)誤差分析表 Tab.2 Prediction error analysis table of seven flood discharge processes

    從表2可看出,武夷山水文站7場(chǎng)洪水流量過(guò)程相對(duì)誤差的平均值均小于20%,且確定性系數(shù)均在0.9以上.經(jīng)過(guò)分析,相對(duì)誤差最大值均出現(xiàn)在低水部分,由于低水流量值比較小,導(dǎo)致相對(duì)誤差較大.從表3可看出,各場(chǎng)洪水的實(shí)測(cè)洪峰流量與預(yù)報(bào)洪峰流量的相對(duì)誤差均在10%以內(nèi),預(yù)報(bào)精度符合規(guī)范要求.可以看出該模型對(duì)于洪水過(guò)程高水部分及洪峰流量部分預(yù)報(bào)效果較好.

    表3 洪峰流量誤差分析表Tab.3 Error analysis table of flood peak flow

    3 結(jié)語(yǔ)

    運(yùn)用泰森多邊形法對(duì)流域進(jìn)行單元?jiǎng)澐郑_定流域內(nèi)6個(gè)雨量站和武夷山水文站的控制流域面積權(quán)重,同時(shí)確定上游6個(gè)雨量站的凈雨形成的徑流到武夷山水文站的匯流時(shí)間.選擇崇陽(yáng)溪上游流域6個(gè)雨量站的時(shí)段雨量與其相應(yīng)的控制流域面積權(quán)重的乘積、武夷山水文站前期流量作為L(zhǎng)MBP模型的輸入,武夷山水文站的相應(yīng)流量作為輸出,建立崇陽(yáng)溪上游流域LMBP神經(jīng)網(wǎng)絡(luò)降雨徑流模型.結(jié)果表明,該模型運(yùn)算時(shí)間短,時(shí)效性好.除低水部分個(gè)別數(shù)據(jù)外,其余部分預(yù)報(bào)精度符合規(guī)范要求,特別是中高水及洪峰流量部分預(yù)報(bào)準(zhǔn)確度高.因此該模型可用于山區(qū)中小流域洪水預(yù)警預(yù)報(bào).

    猜你喜歡
    崇陽(yáng)雨量站武夷山
    遼西山洪災(zāi)害重點(diǎn)防治區(qū)遙測(cè)雨量站點(diǎn)優(yōu)化布設(shè)研究
    它們的家園——武夷山
    等天晴的時(shí)候
    信息熵方法在遼寧省不同分區(qū)雨量站網(wǎng)布設(shè)的應(yīng)用研究
    大美武夷山
    崇陽(yáng)《雙合蓮》的自由抗?fàn)幣c詩(shī)性傳統(tǒng)
    歌海(2020年4期)2020-08-25 01:57:49
    風(fēng)雨廊橋——武夷山馀慶橋
    紅土地(2018年12期)2018-04-29 09:16:50
    武夷山
    中小流域雨量站網(wǎng)密度規(guī)劃與研究
    大力實(shí)施商標(biāo)戰(zhàn)略 建設(shè)鄂南品牌強(qiáng)鎮(zhèn)——崇陽(yáng)市天城鎮(zhèn)發(fā)展紀(jì)實(shí)
    亚洲一区二区三区欧美精品| 日本vs欧美在线观看视频| 男人添女人高潮全过程视频| 亚洲精品久久午夜乱码| 伦理电影免费视频| 久久人人爽人人爽人人片va| 18在线观看网站| 欧美日韩精品成人综合77777| 爱豆传媒免费全集在线观看| 日韩中文字幕视频在线看片| 国产女主播在线喷水免费视频网站| 国产午夜精品一二区理论片| 成年人免费黄色播放视频| 国产永久视频网站| 日韩三级伦理在线观看| 欧美性感艳星| 热99国产精品久久久久久7| 国产无遮挡羞羞视频在线观看| 久久久久久久亚洲中文字幕| 一本大道久久a久久精品| 欧美日韩av久久| 国产精品国产三级国产专区5o| 久热久热在线精品观看| 亚洲av不卡在线观看| 91精品国产九色| 视频在线观看一区二区三区| 国产免费现黄频在线看| 岛国毛片在线播放| 极品少妇高潮喷水抽搐| 欧美人与善性xxx| 精品熟女少妇av免费看| 精品国产乱码久久久久久小说| 久久精品国产自在天天线| 日本色播在线视频| 最新中文字幕久久久久| 我的女老师完整版在线观看| 啦啦啦在线观看免费高清www| 久久ye,这里只有精品| 在线观看国产h片| 男男h啪啪无遮挡| 精品卡一卡二卡四卡免费| 最近的中文字幕免费完整| 亚洲精品中文字幕在线视频| 美女视频免费永久观看网站| 极品少妇高潮喷水抽搐| 国产免费又黄又爽又色| a级片在线免费高清观看视频| 在线播放无遮挡| 久久久国产精品麻豆| 夜夜爽夜夜爽视频| 国产欧美亚洲国产| 亚洲精品亚洲一区二区| 欧美日韩成人在线一区二区| 亚洲怡红院男人天堂| 高清av免费在线| 综合色丁香网| 国产精品三级大全| 能在线免费看毛片的网站| 日韩一区二区三区影片| 亚洲欧美色中文字幕在线| 9色porny在线观看| 美女福利国产在线| 久久久久人妻精品一区果冻| 亚洲综合精品二区| 2021少妇久久久久久久久久久| 晚上一个人看的免费电影| 午夜免费观看性视频| 黄色一级大片看看| 波野结衣二区三区在线| 欧美日本中文国产一区发布| 老司机影院毛片| 亚洲成人av在线免费| av网站免费在线观看视频| 免费高清在线观看视频在线观看| 99热国产这里只有精品6| 国产精品人妻久久久影院| 久久影院123| 精品久久久精品久久久| 狂野欧美白嫩少妇大欣赏| 丝袜在线中文字幕| 精品亚洲成a人片在线观看| 热99久久久久精品小说推荐| 亚洲成色77777| 91aial.com中文字幕在线观看| 天堂8中文在线网| 国产无遮挡羞羞视频在线观看| 国产探花极品一区二区| 极品人妻少妇av视频| 精品久久久精品久久久| 午夜视频国产福利| 成人黄色视频免费在线看| 校园人妻丝袜中文字幕| 在线观看人妻少妇| 最新的欧美精品一区二区| 视频中文字幕在线观看| av在线app专区| 免费观看在线日韩| 久久影院123| 丝袜在线中文字幕| 免费观看av网站的网址| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女视频黄频| 美女国产高潮福利片在线看| 97超碰精品成人国产| 在线观看免费日韩欧美大片 | 性高湖久久久久久久久免费观看| 久久久久久久久大av| 国产精品久久久久成人av| 青春草视频在线免费观看| 久久热精品热| 中文字幕av电影在线播放| 久久精品国产鲁丝片午夜精品| 国产熟女欧美一区二区| 夫妻午夜视频| 欧美成人精品欧美一级黄| 亚洲中文av在线| 国产一区二区在线观看av| av不卡在线播放| 久久婷婷青草| 最新中文字幕久久久久| 黄色视频在线播放观看不卡| 亚洲精品美女久久av网站| 看非洲黑人一级黄片| 一级毛片我不卡| 黄色欧美视频在线观看| 国产黄频视频在线观看| 热re99久久国产66热| 一级爰片在线观看| 蜜臀久久99精品久久宅男| 丰满迷人的少妇在线观看| 亚洲国产欧美日韩在线播放| 国产成人精品久久久久久| √禁漫天堂资源中文www| 久久久精品区二区三区| 欧美日韩一区二区视频在线观看视频在线| 国产日韩一区二区三区精品不卡 | 亚洲激情五月婷婷啪啪| 精品少妇内射三级| 考比视频在线观看| 日本午夜av视频| 国模一区二区三区四区视频| 亚洲欧美成人精品一区二区| 日韩欧美一区视频在线观看| 一区二区三区精品91| 亚洲av电影在线观看一区二区三区| 婷婷色av中文字幕| 国产免费福利视频在线观看| 狂野欧美激情性xxxx在线观看| 女的被弄到高潮叫床怎么办| 亚洲国产精品一区二区三区在线| 综合色丁香网| 91成人精品电影| 两个人的视频大全免费| 人人妻人人添人人爽欧美一区卜| 国产色爽女视频免费观看| 满18在线观看网站| 久久人人爽av亚洲精品天堂| 99久久人妻综合| 亚洲av日韩在线播放| 国产片内射在线| 国产精品久久久久久久久免| 永久免费av网站大全| 亚洲av电影在线观看一区二区三区| 熟妇人妻不卡中文字幕| av网站免费在线观看视频| 人人妻人人添人人爽欧美一区卜| 国产一区有黄有色的免费视频| 亚洲一区二区三区欧美精品| 亚洲综合精品二区| 国产在线视频一区二区| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区国产| 日韩熟女老妇一区二区性免费视频| 日韩中文字幕视频在线看片| 久久久久久久大尺度免费视频| 日本免费在线观看一区| 久久99热这里只频精品6学生| 在线观看人妻少妇| 黄片无遮挡物在线观看| 国产精品久久久久久精品电影小说| 亚洲久久久国产精品| 亚洲综合色惰| 国产精品久久久久久精品古装| 成人亚洲欧美一区二区av| 一级a做视频免费观看| 18禁动态无遮挡网站| 欧美激情国产日韩精品一区| 日本色播在线视频| 一级毛片我不卡| 秋霞伦理黄片| 久久99一区二区三区| av福利片在线| 欧美国产精品一级二级三级| 婷婷色综合www| 亚洲av在线观看美女高潮| 中国美白少妇内射xxxbb| 交换朋友夫妻互换小说| 菩萨蛮人人尽说江南好唐韦庄| 国语对白做爰xxxⅹ性视频网站| 国产精品蜜桃在线观看| 国产精品人妻久久久久久| 天天操日日干夜夜撸| 亚洲一级一片aⅴ在线观看| 午夜福利影视在线免费观看| 亚洲图色成人| 黄色一级大片看看| 91aial.com中文字幕在线观看| 色视频在线一区二区三区| 欧美亚洲日本最大视频资源| 两个人的视频大全免费| videossex国产| a级片在线免费高清观看视频| 99热网站在线观看| 嘟嘟电影网在线观看| 国产免费福利视频在线观看| 热re99久久国产66热| 两个人免费观看高清视频| 国产免费一级a男人的天堂| 欧美国产精品一级二级三级| 啦啦啦啦在线视频资源| 国产精品一区二区三区四区免费观看| 午夜久久久在线观看| 18禁观看日本| 黄片无遮挡物在线观看| 狂野欧美激情性bbbbbb| 亚洲欧美精品自产自拍| 午夜91福利影院| a级片在线免费高清观看视频| 亚洲av.av天堂| 韩国高清视频一区二区三区| 国产有黄有色有爽视频| 亚洲成人手机| 十八禁网站网址无遮挡| 搡老乐熟女国产| av在线app专区| 91aial.com中文字幕在线观看| 亚洲av免费高清在线观看| 国产探花极品一区二区| 搡老乐熟女国产| 少妇精品久久久久久久| 亚洲国产av影院在线观看| 狂野欧美激情性bbbbbb| 久久99一区二区三区| 少妇猛男粗大的猛烈进出视频| 一区二区日韩欧美中文字幕 | 日本猛色少妇xxxxx猛交久久| 亚洲av不卡在线观看| 一级,二级,三级黄色视频| 一区二区三区乱码不卡18| 亚洲无线观看免费| 国产精品欧美亚洲77777| 精品亚洲乱码少妇综合久久| 国产极品粉嫩免费观看在线 | 国产女主播在线喷水免费视频网站| 少妇猛男粗大的猛烈进出视频| 国产成人午夜福利电影在线观看| 日韩av免费高清视频| 亚洲精品一区蜜桃| 最黄视频免费看| 久久久久国产精品人妻一区二区| 制服丝袜香蕉在线| 夜夜骑夜夜射夜夜干| 岛国毛片在线播放| 最近中文字幕高清免费大全6| 99久久精品一区二区三区| 亚洲综合色惰| 久久精品国产亚洲av天美| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 日韩人妻高清精品专区| 春色校园在线视频观看| 欧美日本中文国产一区发布| 视频中文字幕在线观看| 欧美变态另类bdsm刘玥| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 又大又黄又爽视频免费| 18禁动态无遮挡网站| 18在线观看网站| 国产老妇伦熟女老妇高清| 免费av不卡在线播放| 天美传媒精品一区二区| 欧美激情国产日韩精品一区| 国产精品女同一区二区软件| 97超碰精品成人国产| 国产亚洲欧美精品永久| 丝袜喷水一区| 美女国产视频在线观看| 自线自在国产av| 91aial.com中文字幕在线观看| 国产成人aa在线观看| 欧美三级亚洲精品| 一本大道久久a久久精品| 日韩熟女老妇一区二区性免费视频| 久久综合国产亚洲精品| 人妻一区二区av| 有码 亚洲区| 精品亚洲成国产av| 欧美bdsm另类| 欧美成人精品欧美一级黄| 少妇人妻精品综合一区二区| 国产免费一级a男人的天堂| 午夜老司机福利剧场| 岛国毛片在线播放| 啦啦啦视频在线资源免费观看| 亚洲性久久影院| 精品一区在线观看国产| 亚洲丝袜综合中文字幕| 国产精品国产三级国产av玫瑰| 精品99又大又爽又粗少妇毛片| 久久99精品国语久久久| 日韩一区二区三区影片| 久久99精品国语久久久| 天堂俺去俺来也www色官网| 久久久久久久久久久久大奶| 午夜激情av网站| 国产免费福利视频在线观看| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 最近手机中文字幕大全| 久久女婷五月综合色啪小说| 精品国产国语对白av| 女人久久www免费人成看片| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| 国产精品欧美亚洲77777| 中文字幕久久专区| 亚洲内射少妇av| 夜夜爽夜夜爽视频| 免费久久久久久久精品成人欧美视频 | 少妇熟女欧美另类| 视频在线观看一区二区三区| 久久精品国产亚洲网站| 蜜臀久久99精品久久宅男| 欧美xxxx性猛交bbbb| av网站免费在线观看视频| 我要看黄色一级片免费的| 波野结衣二区三区在线| 丰满少妇做爰视频| 肉色欧美久久久久久久蜜桃| 99久国产av精品国产电影| 国产有黄有色有爽视频| 国内精品宾馆在线| 精品国产一区二区三区久久久樱花| 亚洲成人手机| 人体艺术视频欧美日本| 99久久精品国产国产毛片| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 九九在线视频观看精品| 中文字幕精品免费在线观看视频 | 九九久久精品国产亚洲av麻豆| 在现免费观看毛片| 国产一区二区三区综合在线观看 | 99热这里只有精品一区| 久久毛片免费看一区二区三区| 亚洲精品,欧美精品| 国产老妇伦熟女老妇高清| 精品国产露脸久久av麻豆| 精品久久国产蜜桃| 国产深夜福利视频在线观看| 男女啪啪激烈高潮av片| 亚洲精品国产色婷婷电影| 晚上一个人看的免费电影| 亚洲av日韩在线播放| 国产精品女同一区二区软件| 乱码一卡2卡4卡精品| 国产精品 国内视频| 一区二区三区免费毛片| 国产午夜精品一二区理论片| 亚洲怡红院男人天堂| 建设人人有责人人尽责人人享有的| 国产av码专区亚洲av| 综合色丁香网| 精品视频人人做人人爽| 午夜av观看不卡| 大香蕉97超碰在线| 欧美人与善性xxx| 桃花免费在线播放| 亚洲av综合色区一区| 亚洲欧美一区二区三区黑人 | 少妇高潮的动态图| 在线观看美女被高潮喷水网站| 黄色视频在线播放观看不卡| 秋霞伦理黄片| 亚洲国产av影院在线观看| 如何舔出高潮| 国产精品国产三级专区第一集| 嫩草影院入口| 亚洲人成网站在线观看播放| 九九久久精品国产亚洲av麻豆| 少妇 在线观看| 精品国产露脸久久av麻豆| 综合色丁香网| 国产又色又爽无遮挡免| 99久久中文字幕三级久久日本| 亚洲欧美日韩另类电影网站| 全区人妻精品视频| 国产亚洲最大av| 精品亚洲成国产av| 国产熟女午夜一区二区三区 | 成人亚洲精品一区在线观看| 亚洲精品国产色婷婷电影| av在线老鸭窝| 免费播放大片免费观看视频在线观看| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 51国产日韩欧美| 国产极品天堂在线| 久久久久久久大尺度免费视频| 一区二区av电影网| 久久国产精品男人的天堂亚洲 | 九色亚洲精品在线播放| 成人手机av| 婷婷色av中文字幕| 日韩制服骚丝袜av| 男人添女人高潮全过程视频| 2018国产大陆天天弄谢| 日韩在线高清观看一区二区三区| 日日爽夜夜爽网站| 久久久久久久久久人人人人人人| 国产片特级美女逼逼视频| 一级,二级,三级黄色视频| 久久精品久久久久久久性| 如日韩欧美国产精品一区二区三区 | 国产精品麻豆人妻色哟哟久久| 免费不卡的大黄色大毛片视频在线观看| 男人添女人高潮全过程视频| 丰满迷人的少妇在线观看| 亚洲精品一二三| 久久精品久久精品一区二区三区| 搡女人真爽免费视频火全软件| 日韩 亚洲 欧美在线| 天天操日日干夜夜撸| 美女国产高潮福利片在线看| av播播在线观看一区| 18禁观看日本| 777米奇影视久久| 午夜免费观看性视频| av天堂久久9| 欧美日本中文国产一区发布| 一本久久精品| 国产女主播在线喷水免费视频网站| 国产黄色视频一区二区在线观看| av播播在线观看一区| 亚洲国产成人一精品久久久| av免费在线看不卡| 另类亚洲欧美激情| 亚洲精品久久午夜乱码| 精品国产露脸久久av麻豆| 曰老女人黄片| 国产极品天堂在线| 国产欧美日韩一区二区三区在线 | 两个人免费观看高清视频| 黑人欧美特级aaaaaa片| 亚洲精品视频女| 最新中文字幕久久久久| 妹子高潮喷水视频| 午夜日本视频在线| 水蜜桃什么品种好| 天天操日日干夜夜撸| 性色avwww在线观看| 老司机影院毛片| 丝袜美足系列| 国产探花极品一区二区| 黑人猛操日本美女一级片| 精品久久久噜噜| 日本午夜av视频| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 丰满少妇做爰视频| 免费av不卡在线播放| 国产成人a∨麻豆精品| 国产成人freesex在线| 欧美性感艳星| 蜜桃在线观看..| 夫妻午夜视频| 99久久精品一区二区三区| 久久人人爽人人爽人人片va| 国产日韩欧美在线精品| 七月丁香在线播放| 极品少妇高潮喷水抽搐| 国产一区二区在线观看av| av线在线观看网站| 国产伦精品一区二区三区视频9| 精品99又大又爽又粗少妇毛片| 国产午夜精品久久久久久一区二区三区| 超碰97精品在线观看| 久久久久久久久大av| 考比视频在线观看| 一二三四中文在线观看免费高清| 亚洲国产精品一区二区三区在线| 啦啦啦啦在线视频资源| 最近最新中文字幕免费大全7| 少妇被粗大的猛进出69影院 | 日本爱情动作片www.在线观看| 熟女电影av网| 99热这里只有是精品在线观看| 亚洲国产av影院在线观看| 秋霞伦理黄片| 在线观看免费高清a一片| 欧美xxⅹ黑人| 欧美亚洲日本最大视频资源| 黑人高潮一二区| 人妻一区二区av| 午夜老司机福利剧场| 亚洲精品自拍成人| 晚上一个人看的免费电影| 一本—道久久a久久精品蜜桃钙片| 亚洲精品一二三| 99九九在线精品视频| 亚洲精品亚洲一区二区| 乱人伦中国视频| 欧美国产精品一级二级三级| videosex国产| 少妇人妻精品综合一区二区| 国产一区二区在线观看av| 欧美变态另类bdsm刘玥| 好男人视频免费观看在线| 亚洲精品自拍成人| 97在线人人人人妻| 制服人妻中文乱码| 国产av精品麻豆| 97精品久久久久久久久久精品| 搡女人真爽免费视频火全软件| 日韩欧美一区视频在线观看| 国产精品久久久久久av不卡| 99视频精品全部免费 在线| 九色亚洲精品在线播放| 九九久久精品国产亚洲av麻豆| 欧美变态另类bdsm刘玥| av国产久精品久网站免费入址| 制服丝袜香蕉在线| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 不卡视频在线观看欧美| av不卡在线播放| 男人操女人黄网站| 午夜影院在线不卡| 亚洲色图综合在线观看| 国产永久视频网站| 久久精品国产鲁丝片午夜精品| 中文字幕最新亚洲高清| 在线天堂最新版资源| 热re99久久精品国产66热6| 国产一区二区在线观看日韩| 三上悠亚av全集在线观看| 日韩av不卡免费在线播放| 我的女老师完整版在线观看| 国产欧美日韩一区二区三区在线 | 大香蕉97超碰在线| 全区人妻精品视频| 美女cb高潮喷水在线观看| 免费高清在线观看日韩| 王馨瑶露胸无遮挡在线观看| 亚洲成人av在线免费| 最近手机中文字幕大全| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| av播播在线观看一区| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 午夜影院在线不卡| 色哟哟·www| 亚洲人成网站在线观看播放| 青青草视频在线视频观看| 国产在线免费精品| 一本久久精品| 久久热精品热| 欧美3d第一页| 国产一区二区在线观看av| 国产日韩欧美在线精品| 欧美少妇被猛烈插入视频| 91精品三级在线观看| 亚洲婷婷狠狠爱综合网| 亚洲第一区二区三区不卡| 97精品久久久久久久久久精品| 街头女战士在线观看网站| 熟妇人妻不卡中文字幕| 国产精品一区二区在线不卡| 日本-黄色视频高清免费观看| 特大巨黑吊av在线直播| 精品久久国产蜜桃| 最近中文字幕2019免费版| 国产精品久久久久久精品古装| 亚洲内射少妇av| 国产高清不卡午夜福利| 亚洲成人一二三区av| 自拍欧美九色日韩亚洲蝌蚪91| 五月开心婷婷网| 国产精品国产三级国产专区5o| 亚洲av福利一区| 亚洲经典国产精华液单| 国产69精品久久久久777片| 一边亲一边摸免费视频| 亚洲成色77777| 妹子高潮喷水视频| 精品久久国产蜜桃| 国产av码专区亚洲av| 久久久久久久大尺度免费视频| 亚洲精品国产色婷婷电影| 免费观看av网站的网址| 成年av动漫网址| 青春草亚洲视频在线观看| 高清黄色对白视频在线免费看| 久久精品国产a三级三级三级| 人人澡人人妻人| 永久免费av网站大全| av黄色大香蕉| 一区二区三区四区激情视频| 这个男人来自地球电影免费观看 | 国国产精品蜜臀av免费| 九九在线视频观看精品| 国产精品蜜桃在线观看| 搡老乐熟女国产|