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

    基于模型預(yù)測控制的泵站群聯(lián)合調(diào)度研究

    2023-12-06 06:42:44張少愷
    海河水利 2023年11期
    關(guān)鍵詞:泵站時(shí)段水位

    張少愷,龍 巖,管 一,李 樹

    (1.河北工程大學(xué)水利水電學(xué)院,河北 邯鄲 056038;2.河北工程大學(xué)河北省智慧水利重點(diǎn)實(shí)驗(yàn)室,河北 邯鄲 056038;3.通遼市水利事業(yè)發(fā)展中心,內(nèi)蒙古 通遼 028000)

    受全球氣候變暖的影響和城市化的快速發(fā)展,由極端惡劣天氣引發(fā)的城市洪澇災(zāi)害愈發(fā)嚴(yán)重,受災(zāi)范圍愈發(fā)廣泛[1]。泵站群作為福州市排澇的主要工程措施,在汛期發(fā)揮重大作用[2]。本文以福州市汛期承擔(dān)主要排水功能的晉安河流域?yàn)檠芯繉?duì)象,從安全、經(jīng)濟(jì)層面考慮構(gòu)建泵站群聯(lián)合防洪優(yōu)化調(diào)度模型,同時(shí)引入模型預(yù)測控制(MPC:Model Predictive Control)策略,通過優(yōu)化算法求解,提出福州市晉安河流域汛期調(diào)度方案。

    1 研究區(qū)域

    1.1 流域狀況

    晉安河發(fā)源于福州市北峰,橫跨福州市晉安區(qū)南北,是福州市最長的城市內(nèi)河。晉安河全長6.7 km,流域面積69.7 km2。晉安河下游為光明港,最后注入閩江。

    1.2 水利工程現(xiàn)狀

    晉安河上游為琴亭湖,湖水日常水位保持在4.50 m 左右。琴亭湖周邊地面高程均在7.80 m 以上,當(dāng)湖水位在7.80 m 以下時(shí),周邊不會(huì)產(chǎn)生內(nèi)澇。琴亭湖泵站水泵數(shù)量為5 臺(tái),總裝機(jī)3000 kW,總設(shè)計(jì)流量25 m3/s,平均設(shè)計(jì)流量5 m3/s,設(shè)計(jì)揚(yáng)程0~7.0 m。光明港下游的魁岐排澇站承擔(dān)主要排水任務(wù),水泵數(shù)量為3 臺(tái),總裝機(jī)6600 kW,總設(shè)計(jì)流量120 m3/s,平均設(shè)計(jì)流量40 m3/s,設(shè)計(jì)凈揚(yáng)程2.6 m。琴亭湖、魁岐泵站群調(diào)度概化,如圖1所示。

    圖1 琴亭湖、魁岐泵站群調(diào)度概化

    2 模型預(yù)測控制策略

    為了更好地解決模型優(yōu)化過程中的失配、畸變、干擾等因素引起的不確定性[3],引入基于MPC 的滾動(dòng)優(yōu)化策略。MPC 是一種滾動(dòng)優(yōu)化控制算法,主要包括模型預(yù)測、滾動(dòng)優(yōu)化、反饋校正3 個(gè)環(huán)節(jié)。MPC滾動(dòng)優(yōu)化策略,如圖2所示。

    圖2 MPC滾動(dòng)優(yōu)化策略示意

    圖2中,Np為預(yù)測時(shí)域,Nc為控制時(shí)域。在t時(shí),利用預(yù)測時(shí)域內(nèi)的預(yù)測數(shù)據(jù),在控制時(shí)域內(nèi)進(jìn)行多目標(biāo)優(yōu)化,得到控制時(shí)域內(nèi)最優(yōu)的控制序列;但在當(dāng)前t時(shí)段,只執(zhí)行第一個(gè)控制序列的指令。令t=t+1,重復(fù)進(jìn)行上一過程,直至完成整個(gè)時(shí)間軸的優(yōu)化。

    3 河道水動(dòng)力模型構(gòu)建

    為實(shí)現(xiàn)對(duì)晉安河的洪澇過程模擬,需獲得琴亭湖上游來水和琴亭湖至光明港的區(qū)間來水,建立基于SWMM 的晉安河河道水力學(xué)模型。SWMM 模型,如圖3所示。

    圖3 SWMM 模型

    通過查閱相關(guān)文獻(xiàn),以三八泵站和魁岐排澇站為例,模擬2021年福州市的2 次降雨,采用Nash-Sutcliffe 效率系數(shù)法,對(duì)SWMM 模型參數(shù)進(jìn)行率定與驗(yàn)證。模擬水位過程線與實(shí)測水位過程線形狀比較吻合,水位峰值及峰值出現(xiàn)時(shí)間的偏差均較小。三八泵站和魁岐排澇站6月29日降雨事件率定結(jié)果如圖4—5所示,其納什效率系數(shù)分別為0.53、0.92;8月8 日降雨事件驗(yàn)證結(jié)果如圖6—7 所示,其納什效率系數(shù)分別為0.68、0.55。

    圖4 三八泵站率定結(jié)果

    圖5 魁岐排澇站率定結(jié)果

    圖6 三八泵站驗(yàn)證結(jié)果

    圖7 魁岐排澇站驗(yàn)證結(jié)果

    4 基于模型預(yù)測控制的泵站群調(diào)度模型構(gòu)建

    4.1 預(yù)測模型

    本文在泵站群優(yōu)化調(diào)度模型內(nèi)建立基于水量平衡原理的預(yù)測模型,選擇時(shí)間間隔為1 h 進(jìn)行采樣,對(duì)未來8 h的水位進(jìn)行預(yù)測。預(yù)測模型為:

    式中:Vi,t為i水庫在t時(shí)段的庫容值(m3);Vi,t+1為i水庫在t+ 1 時(shí)段的庫容值(m3);Ii,t為i水庫在t時(shí)段的入庫流量(m3/s);Qi,t為i水庫在t時(shí)段的出庫流量(m3/s);Vi,0為i水庫初始庫容值(m3);b0i為i水庫初始水位值(m);Vi,T為i水庫T時(shí)段末庫容值(m3);bT i為i水庫T時(shí)段末水位值(m);t為調(diào)度時(shí)段( )

    h ,t=1,2,…,T;T為調(diào)度總時(shí)段( )h ,取8 h;ΔT表示單位時(shí)間長度為1 h;i=1代表琴亭湖,i=2代表光明港區(qū)域。

    4.2 梯級(jí)泵站優(yōu)化調(diào)度模型

    4.2.1 目標(biāo)函數(shù)

    (1)調(diào)度優(yōu)化目標(biāo)1。琴亭湖超警風(fēng)險(xiǎn)最小。其目標(biāo)函數(shù)為:

    式中:F1為調(diào)度周期內(nèi)琴亭湖各時(shí)段距死水位的最小值(m);Zx為琴亭湖汛限水位(m);Z1,t為t時(shí)段琴亭湖水位(m);Zs為琴亭湖死水位(m);t、T含義及取值同上。

    (2)調(diào)度優(yōu)化目標(biāo)2。光明港峰值水位最低。其目標(biāo)函數(shù)為:

    式中:F2為調(diào)度周期內(nèi)光明港各時(shí)段峰值水位的最小值(m);Z2,t為t時(shí)段光明港的水位(m);t的含義及取值同上。

    (3)調(diào)度優(yōu)化目標(biāo)3。調(diào)度運(yùn)行期間抽水費(fèi)用最少。其目標(biāo)函數(shù)為:

    式中:F3為調(diào)度周期內(nèi)總抽水費(fèi)用(元);k為電費(fèi)單價(jià)(元/kW·h);N為泵站數(shù)量,n=1 代表琴亭湖泵站,n=2 代表魁岐排澇站;ρ為水的密度(g/cm3);g為重力加速度(m/s2),取9.8 m/s2;Hn,t為第n個(gè)泵站在t時(shí)段的揚(yáng)程(m);Qn,t為第n個(gè)泵站在t時(shí)段的過泵流量(m3/s);η為水泵損失效率;ΔT、t、T含義及取值同上。

    (4)調(diào)度優(yōu)化目標(biāo)4。閘門啟閉次數(shù)最少。泵站的檢修次數(shù)和檢修成本與泵站機(jī)組在使用過程中的啟閉有關(guān)。閘門啟閉次數(shù)的目標(biāo)函數(shù)為:

    式中:F4為調(diào)度期間內(nèi)泵站機(jī)組啟閉次數(shù)之和;L(n,t)為n泵站的t時(shí)段與上一時(shí)段泵站機(jī)組相比變化的臺(tái)數(shù);n= 1 代表琴亭湖泵站,n= 2 代表魁岐排澇站;t、T含義及取值同上。

    (5)調(diào)度優(yōu)化目標(biāo)5。琴亭湖調(diào)度時(shí)段末水位最低。其目標(biāo)函數(shù)為:

    式中:F5為琴亭湖調(diào)度周期內(nèi)各時(shí)段末最低水位(m);Z1,t為t時(shí)段琴亭湖水位(m);t的含義及取值同上。

    (6)調(diào)度優(yōu)化目標(biāo)6。光明港調(diào)度時(shí)段末水位最低。其目標(biāo)函數(shù)為:

    式中:F6為光明港調(diào)度周期內(nèi)各時(shí)段末最低水位(m);Z2,t為t時(shí)段光明港水位(m);t的含義及取值同上。

    4.2.2 約束條件

    (1)水位約束。其約束條件公式為:

    式中:Zn,t為第n個(gè)泵站在t時(shí)段的上游水位( m );、分別為第n個(gè)泵站的上游水位上、下限約束(m),為0,為7.0 m,為0,為5.0 m;n= 1 代表琴亭湖泵站,n= 2 代表魁岐排澇站;t的含義及取值同上。

    (2)流量約束。其約束條件公式為:

    (3)水泵揚(yáng)程約束。其約束條件公式為:

    4.2.3 模型求解

    由于琴亭湖、魁岐泵站群聯(lián)合優(yōu)化調(diào)度是多目標(biāo)優(yōu)化問題,所以采用構(gòu)造懲罰函數(shù)的方法求解。將琴亭湖調(diào)度時(shí)段末水位低于目標(biāo)水位、光明港峰值水位最低和調(diào)度時(shí)段末水位低于目標(biāo)水位、泵站啟閉次數(shù)最少作為懲罰函數(shù),融入抽排水費(fèi)用最低目標(biāo)函數(shù)中,基于JAVA 語言通過Eclipse 軟件編寫優(yōu)化調(diào)度程序,運(yùn)用人工蜂群算法計(jì)算。

    4.3 MPC設(shè)計(jì)

    本文將模型預(yù)測應(yīng)用于泵站群優(yōu)化調(diào)度模型中,在已知當(dāng)前時(shí)刻琴亭湖、光明港的水位、水位庫容曲線和未來一段時(shí)間的入庫流量前提下,進(jìn)行晉安河流域泵站群未來最優(yōu)調(diào)控,使得琴亭湖、光明港充分發(fā)揮調(diào)蓄能力。MPC優(yōu)化調(diào)度方案步驟,如圖8所示。

    圖8 MPC優(yōu)化調(diào)度方案步驟

    結(jié)合圖8,將MPC 應(yīng)用于對(duì)未來時(shí)段晉安河防洪調(diào)度的最優(yōu)控制決策,具體步驟如下。

    步驟1:首先進(jìn)行人工蜂群算法的相關(guān)參數(shù)設(shè)置。①初始種群規(guī)模:本文種群規(guī)模取值為200。種群規(guī)模SN包含采蜜蜂和觀察蜂,各取種群數(shù)量的1/2。②一般食物源為種群的1/2。③設(shè)置limit限制,超過此限制沒有更新的采蜜蜂變?yōu)橛^察蜂。④設(shè)置maxcycle為100。

    步驟2:數(shù)據(jù)輸入。設(shè)定采樣時(shí)間為1 h 和預(yù)測步長為8 h,總調(diào)度時(shí)長為24 h。輸入k時(shí)段未來8 h的琴亭湖上游來水、光明港區(qū)間來水。

    步驟3:在已知k時(shí)刻琴亭湖、光明港初始水位并已知未來8 h入庫流量的前提下,將晉安河泵站群起始數(shù)據(jù)代入預(yù)測模型。

    步驟4:利用人工蜂群算法計(jì)算懲罰函數(shù),求解出使性能指標(biāo)最小的最優(yōu)控制序列,輸出最優(yōu)控制序列中的首個(gè)控制指令進(jìn)行晉安河泵站群優(yōu)化調(diào)度。

    步驟5:將MPC 防洪優(yōu)化調(diào)度最優(yōu)控制序列結(jié)果代入河道動(dòng)力模型中,獲取k+ 1 時(shí)刻初始水位,進(jìn)行水位反饋校正。

    步驟6:進(jìn)入下一時(shí)刻最優(yōu)控制序列的求解。

    5 實(shí)例分析

    選取100 a 一遇設(shè)計(jì)降雨序列,作為雨量邊界,輸入河道動(dòng)力模型獲取琴亭湖、光明港入流過程。琴亭湖主要考慮上游來水,光明港主要考慮琴亭湖至光明港的區(qū)間來水以及魁岐排澇站泵、閘排水方式的影響。當(dāng)福州城區(qū)晉安河水位低于閩江潮位時(shí),將沿江的閘門全部關(guān)閉,若超過警戒水位,此時(shí)開啟魁岐排澇站泵站機(jī)組,降低晉安河水位;當(dāng)福州城區(qū)內(nèi)河水位高于閩江潮位時(shí),將魁岐泵站機(jī)組關(guān)閉,僅采用魁岐水閘排水,將閘內(nèi)水位降到警戒水位以下。琴亭湖、光明港區(qū)間來水過程如圖9所示,閩江潮位如圖10所示。

    圖9 琴亭湖、光明港區(qū)間來水

    圖10 閩江潮位過程

    在100 a一遇的設(shè)計(jì)洪水情況下,結(jié)合所建水力學(xué)模型和實(shí)際防洪排澇需求,得出以下結(jié)論。

    (1)從經(jīng)濟(jì)運(yùn)行角度分析,MPC 優(yōu)化調(diào)度運(yùn)行費(fèi)用為2.36萬元,規(guī)則調(diào)度運(yùn)行費(fèi)用為2.48萬元,運(yùn)行費(fèi)下降了4.8%;從設(shè)備安全及維護(hù)角度分析,MPC優(yōu)化調(diào)度的泵站啟閉次數(shù)為43,規(guī)則調(diào)度的泵站啟閉次數(shù)為46。規(guī)則調(diào)度、MPC 優(yōu)化調(diào)度泵站啟閉情況對(duì)比,詳見表1。

    表1 規(guī)則調(diào)度、MPC優(yōu)化調(diào)度泵站啟閉情況對(duì)比

    (2)從防洪效果上看,MPC 優(yōu)化調(diào)度方案琴亭湖、光明港峰值水位分別為6.15、5.621 m,規(guī)則調(diào)度方案琴亭湖、光明港峰值水位分別為7.11、5.732 m,洪峰水位分別降低了13.5%、1.93%。琴亭湖、光明港100 a一遇調(diào)度,如圖11—12所示。

    圖11 琴亭湖100 a一遇調(diào)度

    圖12 光明港100 a一遇調(diào)度

    (3)從實(shí)際工程角度出發(fā),因MPC 優(yōu)化調(diào)度采用滾動(dòng)策略實(shí)時(shí)計(jì)算,可以應(yīng)對(duì)突發(fā)情況,相比其他依靠既定規(guī)則的優(yōu)化調(diào)度,實(shí)用性大幅提升。

    6 結(jié)語

    本文基于實(shí)際工程需求,構(gòu)建了基于模型預(yù)測控制的泵站群聯(lián)合調(diào)度模型,為福州市江北城區(qū)應(yīng)對(duì)100 a一遇的典型洪水場景提供了優(yōu)化調(diào)度方案,可作為福州市汛期實(shí)時(shí)排澇調(diào)度的參考依據(jù)。

    猜你喜歡
    泵站時(shí)段水位
    張家邊涌泵站建設(shè)難點(diǎn)及技術(shù)創(chuàng)新實(shí)踐
    四個(gè)養(yǎng)生黃金時(shí)段,你抓住了嗎
    2016年河南省己建成泵站數(shù)量
    全省已建成泵站數(shù)量
    基于MFAC-PID的核電站蒸汽發(fā)生器水位控制
    河南省2014年已建成泵站數(shù)量
    傍晚是交通事故高發(fā)時(shí)段
    分時(shí)段預(yù)約在PICC門診維護(hù)中的應(yīng)用與探討
    基于PLC的水位控制系統(tǒng)的設(shè)計(jì)與研究
    河南科技(2014年4期)2014-02-27 14:07:11
    分時(shí)段預(yù)約掛號(hào)的實(shí)現(xiàn)與應(yīng)用
    av网站免费在线观看视频 | 亚洲内射少妇av| 免费观看无遮挡的男女| 免费看av在线观看网站| 久久这里有精品视频免费| 91狼人影院| 肉色欧美久久久久久久蜜桃 | 日韩强制内射视频| 久久久久精品性色| 精品一区二区三卡| 在线天堂最新版资源| 久久久久精品久久久久真实原创| 免费看美女性在线毛片视频| 日韩一区二区三区影片| 久久这里有精品视频免费| 三级国产精品欧美在线观看| 2018国产大陆天天弄谢| 午夜福利在线观看吧| 国产免费视频播放在线视频 | 亚洲国产成人一精品久久久| 如何舔出高潮| 国产在视频线在精品| 成年女人在线观看亚洲视频 | 亚洲经典国产精华液单| 精品国产露脸久久av麻豆 | a级毛片免费高清观看在线播放| 少妇裸体淫交视频免费看高清| 能在线免费看毛片的网站| 在线免费十八禁| 97精品久久久久久久久久精品| 老师上课跳d突然被开到最大视频| 麻豆成人av视频| 男女边摸边吃奶| 少妇的逼好多水| 午夜老司机福利剧场| 丝瓜视频免费看黄片| 国产探花极品一区二区| 欧美bdsm另类| 精品久久久久久久末码| 老女人水多毛片| 人妻制服诱惑在线中文字幕| 婷婷色综合大香蕉| 欧美性猛交╳xxx乱大交人| 亚洲人与动物交配视频| 久久久欧美国产精品| 日韩人妻高清精品专区| 国国产精品蜜臀av免费| 久久精品久久久久久久性| 欧美高清成人免费视频www| 日韩人妻高清精品专区| 一个人免费在线观看电影| 成人鲁丝片一二三区免费| 亚洲国产日韩欧美精品在线观看| 美女大奶头视频| 少妇裸体淫交视频免费看高清| 干丝袜人妻中文字幕| 中国美白少妇内射xxxbb| 亚洲av中文av极速乱| 欧美另类一区| 国产单亲对白刺激| 三级经典国产精品| 国产不卡一卡二| 美女黄网站色视频| 只有这里有精品99| 久久久午夜欧美精品| 成人国产麻豆网| 麻豆国产97在线/欧美| 午夜福利视频精品| 日本黄色片子视频| 天堂网av新在线| 伊人久久精品亚洲午夜| 欧美三级亚洲精品| 美女大奶头视频| 麻豆国产97在线/欧美| 卡戴珊不雅视频在线播放| 麻豆成人午夜福利视频| 99热这里只有精品一区| 可以在线观看毛片的网站| 成人亚洲欧美一区二区av| 神马国产精品三级电影在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产乱人偷精品视频| 少妇高潮的动态图| 日韩国内少妇激情av| 日韩av不卡免费在线播放| 一级毛片电影观看| 欧美 日韩 精品 国产| 亚洲真实伦在线观看| av又黄又爽大尺度在线免费看| 超碰av人人做人人爽久久| 五月伊人婷婷丁香| 亚洲美女视频黄频| 欧美精品国产亚洲| 久久久久久国产a免费观看| av卡一久久| 青春草国产在线视频| 久久久a久久爽久久v久久| 亚洲无线观看免费| 午夜日本视频在线| 最近最新中文字幕免费大全7| 最近视频中文字幕2019在线8| 亚洲国产最新在线播放| 成年免费大片在线观看| 国产探花在线观看一区二区| 亚洲成人av在线免费| 91午夜精品亚洲一区二区三区| 成人国产麻豆网| 淫秽高清视频在线观看| 日日干狠狠操夜夜爽| 麻豆乱淫一区二区| 国产伦在线观看视频一区| 日日啪夜夜撸| eeuss影院久久| 日韩成人av中文字幕在线观看| 简卡轻食公司| 亚洲国产成人一精品久久久| 精品国内亚洲2022精品成人| 白带黄色成豆腐渣| 老司机影院毛片| 免费观看在线日韩| av播播在线观看一区| 看非洲黑人一级黄片| 五月玫瑰六月丁香| 日韩成人av中文字幕在线观看| 精品欧美国产一区二区三| 国产亚洲一区二区精品| 国产亚洲av片在线观看秒播厂 | 日本三级黄在线观看| 免费观看精品视频网站| 三级经典国产精品| 一级毛片黄色毛片免费观看视频| 日本黄大片高清| 日本与韩国留学比较| 97超视频在线观看视频| 久久99热这里只频精品6学生| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美清纯卡通| 97热精品久久久久久| 插阴视频在线观看视频| 麻豆av噜噜一区二区三区| 日韩视频在线欧美| 亚洲最大成人中文| 国产永久视频网站| 精品久久久久久电影网| 乱人视频在线观看| 狂野欧美白嫩少妇大欣赏| 黄色欧美视频在线观看| 女人久久www免费人成看片| 国产午夜福利久久久久久| 91久久精品电影网| 国产精品一区二区性色av| 午夜福利在线在线| 免费不卡的大黄色大毛片视频在线观看 | 国产高清有码在线观看视频| 91久久精品国产一区二区成人| 久久人人爽人人片av| 国产综合懂色| 黄色一级大片看看| 国产午夜精品久久久久久一区二区三区| 中文精品一卡2卡3卡4更新| 国产免费视频播放在线视频 | 街头女战士在线观看网站| 成人特级av手机在线观看| 午夜爱爱视频在线播放| 国产日韩欧美在线精品| 亚洲成人av在线免费| 国产精品av视频在线免费观看| 欧美性感艳星| 少妇熟女欧美另类| 啦啦啦韩国在线观看视频| kizo精华| 国产免费福利视频在线观看| 久久这里只有精品中国| 能在线免费看毛片的网站| 日韩欧美一区视频在线观看 | 伊人久久国产一区二区| 亚洲精品,欧美精品| av卡一久久| 99久久精品一区二区三区| 少妇人妻一区二区三区视频| 亚洲四区av| 亚洲欧美成人精品一区二区| 777米奇影视久久| 婷婷色av中文字幕| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 亚洲三级黄色毛片| 热99在线观看视频| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| or卡值多少钱| 中文字幕久久专区| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| 伦精品一区二区三区| 国产成人免费观看mmmm| 国产v大片淫在线免费观看| 亚洲电影在线观看av| 国国产精品蜜臀av免费| 免费黄频网站在线观看国产| eeuss影院久久| 国产精品久久久久久精品电影小说 | 国产69精品久久久久777片| 成人毛片60女人毛片免费| 国产亚洲午夜精品一区二区久久 | 国产白丝娇喘喷水9色精品| 国产亚洲91精品色在线| 高清在线视频一区二区三区| 91aial.com中文字幕在线观看| 婷婷色综合www| 欧美+日韩+精品| 日韩大片免费观看网站| 少妇熟女欧美另类| 尾随美女入室| 91午夜精品亚洲一区二区三区| 日韩三级伦理在线观看| 亚洲综合精品二区| 国产探花极品一区二区| 成人特级av手机在线观看| 久久久久久久久久久丰满| 久久久午夜欧美精品| 一区二区三区乱码不卡18| 91精品一卡2卡3卡4卡| 国产精品一及| 中文字幕免费在线视频6| 我的老师免费观看完整版| 久久99热这里只频精品6学生| 青春草视频在线免费观看| 3wmmmm亚洲av在线观看| 麻豆乱淫一区二区| 国产永久视频网站| 久久久色成人| 亚洲成色77777| 青春草亚洲视频在线观看| 男女那种视频在线观看| 免费播放大片免费观看视频在线观看| 天堂俺去俺来也www色官网 | 性色avwww在线观看| 内射极品少妇av片p| 亚洲国产欧美在线一区| 中文在线观看免费www的网站| 中文欧美无线码| 好男人在线观看高清免费视频| 女人久久www免费人成看片| a级毛片免费高清观看在线播放| 精品人妻一区二区三区麻豆| 国产一区二区在线观看日韩| 午夜激情欧美在线| 你懂的网址亚洲精品在线观看| 免费大片18禁| 亚洲人成网站在线观看播放| 丰满人妻一区二区三区视频av| 极品教师在线视频| 亚洲成人精品中文字幕电影| 在线免费观看不下载黄p国产| 国产有黄有色有爽视频| 91精品伊人久久大香线蕉| 91精品一卡2卡3卡4卡| xxx大片免费视频| 高清毛片免费看| 国产亚洲5aaaaa淫片| 春色校园在线视频观看| 亚洲成人av在线免费| 日本av手机在线免费观看| 大又大粗又爽又黄少妇毛片口| 久久精品夜夜夜夜夜久久蜜豆| 三级毛片av免费| 国产一级毛片七仙女欲春2| 日韩av在线大香蕉| 人妻制服诱惑在线中文字幕| 午夜亚洲福利在线播放| 国产伦理片在线播放av一区| 噜噜噜噜噜久久久久久91| 亚洲人成网站高清观看| 色吧在线观看| av黄色大香蕉| 亚洲欧美成人综合另类久久久| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 午夜免费激情av| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 日韩,欧美,国产一区二区三区| 国产亚洲午夜精品一区二区久久 | 永久网站在线| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 听说在线观看完整版免费高清| 老女人水多毛片| 黄片wwwwww| 99热全是精品| 男女边摸边吃奶| 美女cb高潮喷水在线观看| 岛国毛片在线播放| 久久6这里有精品| 久久久欧美国产精品| 国产成年人精品一区二区| 干丝袜人妻中文字幕| av网站免费在线观看视频 | 黄色一级大片看看| 久久6这里有精品| 国产伦在线观看视频一区| 一级毛片电影观看| 中文字幕av在线有码专区| 国产欧美另类精品又又久久亚洲欧美| 免费黄频网站在线观看国产| 丰满人妻一区二区三区视频av| 免费黄网站久久成人精品| 亚洲国产精品sss在线观看| 亚洲精品成人av观看孕妇| 我的老师免费观看完整版| 春色校园在线视频观看| 欧美最新免费一区二区三区| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 午夜福利在线观看吧| 黄色一级大片看看| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 国产免费视频播放在线视频 | 国产亚洲av嫩草精品影院| 成年版毛片免费区| 亚洲国产精品成人综合色| 国产成人一区二区在线| 天美传媒精品一区二区| 少妇熟女aⅴ在线视频| 久久综合国产亚洲精品| 久久久久免费精品人妻一区二区| 免费看美女性在线毛片视频| 搞女人的毛片| 久久97久久精品| 高清视频免费观看一区二区 | 三级经典国产精品| 精品久久久久久久末码| 精品国产一区二区三区久久久樱花 | 国产又色又爽无遮挡免| 成人特级av手机在线观看| 欧美潮喷喷水| 精品午夜福利在线看| 极品少妇高潮喷水抽搐| 97超视频在线观看视频| 亚洲天堂国产精品一区在线| 国产精品一及| 天堂俺去俺来也www色官网 | 亚洲精品,欧美精品| 一夜夜www| 亚洲美女搞黄在线观看| 亚洲欧美精品自产自拍| 好男人在线观看高清免费视频| 欧美三级亚洲精品| 久久久久网色| 两个人的视频大全免费| 国产精品女同一区二区软件| 欧美极品一区二区三区四区| 夜夜看夜夜爽夜夜摸| 我要看日韩黄色一级片| 高清毛片免费看| 亚洲欧美一区二区三区黑人 | 国产精品嫩草影院av在线观看| 大又大粗又爽又黄少妇毛片口| av播播在线观看一区| 人妻系列 视频| 高清午夜精品一区二区三区| 国产亚洲午夜精品一区二区久久 | 一区二区三区乱码不卡18| 国产精品不卡视频一区二区| 亚洲av一区综合| 久久久久免费精品人妻一区二区| 观看美女的网站| 精品少妇黑人巨大在线播放| 久久6这里有精品| www.av在线官网国产| 亚洲欧美日韩卡通动漫| 夫妻性生交免费视频一级片| 丝袜美腿在线中文| .国产精品久久| av卡一久久| 99久久九九国产精品国产免费| 偷拍熟女少妇极品色| 国产午夜精品一二区理论片| 乱人视频在线观看| 日韩国内少妇激情av| 观看美女的网站| 日本熟妇午夜| 六月丁香七月| 少妇高潮的动态图| 国语对白做爰xxxⅹ性视频网站| 国产熟女欧美一区二区| 国产片特级美女逼逼视频| 大话2 男鬼变身卡| 久99久视频精品免费| 国产高潮美女av| 国产淫片久久久久久久久| ponron亚洲| 男人和女人高潮做爰伦理| 午夜福利网站1000一区二区三区| 少妇人妻精品综合一区二区| 国产av在哪里看| 又爽又黄无遮挡网站| 精品久久久久久久久亚洲| 熟妇人妻不卡中文字幕| 日日撸夜夜添| 啦啦啦韩国在线观看视频| 国产亚洲午夜精品一区二区久久 | 亚洲精品视频女| 一级黄片播放器| 街头女战士在线观看网站| 欧美 日韩 精品 国产| 最新中文字幕久久久久| 又爽又黄a免费视频| 爱豆传媒免费全集在线观看| 成人一区二区视频在线观看| 激情五月婷婷亚洲| 伊人久久精品亚洲午夜| 欧美xxxx黑人xx丫x性爽| 国产精品1区2区在线观看.| av国产久精品久网站免费入址| 亚洲国产高清在线一区二区三| 在现免费观看毛片| 国产极品天堂在线| 女的被弄到高潮叫床怎么办| 高清日韩中文字幕在线| 秋霞伦理黄片| 美女xxoo啪啪120秒动态图| 亚洲最大成人中文| 一级毛片我不卡| 国产免费一级a男人的天堂| 能在线免费观看的黄片| 亚洲成人精品中文字幕电影| 色尼玛亚洲综合影院| 国产成年人精品一区二区| 成人欧美大片| 国产黄色小视频在线观看| 国产白丝娇喘喷水9色精品| 一本久久精品| 91久久精品电影网| 成人鲁丝片一二三区免费| 男女下面进入的视频免费午夜| av线在线观看网站| 精品午夜福利在线看| 赤兔流量卡办理| 午夜亚洲福利在线播放| 色网站视频免费| 国产片特级美女逼逼视频| 2022亚洲国产成人精品| 狠狠精品人妻久久久久久综合| 久久99热这里只有精品18| 亚洲人成网站高清观看| 国产高清国产精品国产三级 | 中文字幕av在线有码专区| 日本免费a在线| 人人妻人人澡人人爽人人夜夜 | 免费av不卡在线播放| 一本一本综合久久| 日日啪夜夜撸| 丝瓜视频免费看黄片| 久久久久久久久久黄片| 亚洲精品成人av观看孕妇| 寂寞人妻少妇视频99o| 精品久久久久久久末码| 人妻一区二区av| 国产黄a三级三级三级人| av在线老鸭窝| 免费播放大片免费观看视频在线观看| 国产视频内射| 22中文网久久字幕| 精品久久久精品久久久| 国产黄片美女视频| 欧美zozozo另类| 成人高潮视频无遮挡免费网站| 有码 亚洲区| 在线观看人妻少妇| 色综合亚洲欧美另类图片| 熟妇人妻不卡中文字幕| 欧美日韩精品成人综合77777| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 在线免费十八禁| 亚洲精品自拍成人| 成年人午夜在线观看视频 | 99久久精品热视频| 中文字幕免费在线视频6| 深夜a级毛片| 免费看av在线观看网站| 国产黄色免费在线视频| 久久99热6这里只有精品| 精品久久久久久久久亚洲| 又黄又爽又刺激的免费视频.| 国产黄a三级三级三级人| 精品少妇黑人巨大在线播放| 美女cb高潮喷水在线观看| 久久久亚洲精品成人影院| 特大巨黑吊av在线直播| 黑人高潮一二区| 中国美白少妇内射xxxbb| 久久韩国三级中文字幕| ponron亚洲| 天堂俺去俺来也www色官网 | 99久国产av精品| 国产精品一及| 中文资源天堂在线| 美女内射精品一级片tv| ponron亚洲| 久热久热在线精品观看| 97超视频在线观看视频| 80岁老熟妇乱子伦牲交| 天天躁日日操中文字幕| 国产免费一级a男人的天堂| 成人二区视频| 亚洲av男天堂| 尾随美女入室| 美女内射精品一级片tv| 男女边摸边吃奶| 国产一级毛片在线| 亚洲av免费在线观看| 久久这里有精品视频免费| 久久6这里有精品| 全区人妻精品视频| 亚洲一区高清亚洲精品| 亚洲国产精品专区欧美| 亚洲精品一二三| 女人十人毛片免费观看3o分钟| 成年女人看的毛片在线观看| 韩国av在线不卡| 美女国产视频在线观看| 亚洲av中文字字幕乱码综合| 校园人妻丝袜中文字幕| 国产熟女欧美一区二区| 美女高潮的动态| 久久99热这里只有精品18| 午夜免费男女啪啪视频观看| 国产黄片视频在线免费观看| 亚洲电影在线观看av| 成人鲁丝片一二三区免费| 久久久久久伊人网av| 中文字幕av在线有码专区| 一个人观看的视频www高清免费观看| ponron亚洲| 免费在线观看成人毛片| 91午夜精品亚洲一区二区三区| 男人和女人高潮做爰伦理| 成人毛片a级毛片在线播放| 久久久久九九精品影院| 全区人妻精品视频| 超碰97精品在线观看| 精品一区二区免费观看| 国产伦精品一区二区三区四那| 国精品久久久久久国模美| 免费av不卡在线播放| 精品久久久久久久末码| av福利片在线观看| 久久6这里有精品| 国产一区二区三区综合在线观看 | 少妇被粗大猛烈的视频| 亚洲精品一区蜜桃| 一级毛片 在线播放| 国产高清三级在线| 国产一区亚洲一区在线观看| 真实男女啪啪啪动态图| 免费电影在线观看免费观看| 亚洲av成人精品一区久久| 国产av不卡久久| 久久精品熟女亚洲av麻豆精品 | 韩国av在线不卡| 69av精品久久久久久| 全区人妻精品视频| 亚洲精品456在线播放app| 最近最新中文字幕免费大全7| 久久亚洲国产成人精品v| 午夜福利在线观看吧| 国产黄频视频在线观看| 国产精品麻豆人妻色哟哟久久 | 免费观看av网站的网址| 秋霞在线观看毛片| 我的女老师完整版在线观看| 91狼人影院| 国产老妇伦熟女老妇高清| 亚洲乱码一区二区免费版| 夜夜爽夜夜爽视频| 久久99热这里只频精品6学生| 日韩欧美 国产精品| 晚上一个人看的免费电影| 床上黄色一级片| 国产成人精品福利久久| 秋霞在线观看毛片| 免费av毛片视频| 搞女人的毛片| 午夜精品在线福利| 亚洲欧洲国产日韩| 精品久久久久久久人妻蜜臀av| 最后的刺客免费高清国语| 亚洲乱码一区二区免费版| 联通29元200g的流量卡| 91精品伊人久久大香线蕉| 亚洲精品日韩在线中文字幕| 国产久久久一区二区三区| 国产视频首页在线观看| 全区人妻精品视频| 91精品国产九色| 男女国产视频网站| 欧美最新免费一区二区三区| 欧美丝袜亚洲另类| 五月伊人婷婷丁香| 午夜久久久久精精品| xxx大片免费视频| 91精品一卡2卡3卡4卡| 亚洲av电影在线观看一区二区三区 | 亚洲成人一二三区av| 欧美bdsm另类| 又黄又爽又刺激的免费视频.| 久久久国产一区二区| 国产成人一区二区在线| 一区二区三区免费毛片| 国产精品一区二区性色av| 韩国高清视频一区二区三区|