• <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)用
    亚洲va在线va天堂va国产| 校园人妻丝袜中文字幕| 汤姆久久久久久久影院中文字幕 | 人妻夜夜爽99麻豆av| 中国美白少妇内射xxxbb| 免费观看精品视频网站| 真实男女啪啪啪动态图| 99久久人妻综合| 蜜桃久久精品国产亚洲av| 乱系列少妇在线播放| 亚洲最大成人手机在线| 国模一区二区三区四区视频| 一区二区三区免费毛片| 国产高清有码在线观看视频| 久久久久久国产a免费观看| 亚洲欧美成人综合另类久久久 | 亚洲欧美日韩卡通动漫| 日韩强制内射视频| 日本猛色少妇xxxxx猛交久久| 国产一区二区在线观看日韩| 91aial.com中文字幕在线观看| 汤姆久久久久久久影院中文字幕 | 久久久久久久亚洲中文字幕| 亚洲人成网站在线播| 晚上一个人看的免费电影| 精品无人区乱码1区二区| 国产美女午夜福利| 性色avwww在线观看| 性插视频无遮挡在线免费观看| 国产精品精品国产色婷婷| 国产一区二区亚洲精品在线观看| 精品久久久久久久久av| 变态另类丝袜制服| 在线免费观看的www视频| 亚洲av男天堂| 色综合亚洲欧美另类图片| 婷婷六月久久综合丁香| 亚洲欧美成人综合另类久久久 | 亚洲国产精品成人久久小说| 久久国内精品自在自线图片| 中文乱码字字幕精品一区二区三区 | 久久精品国产鲁丝片午夜精品| 免费黄网站久久成人精品| 日本三级黄在线观看| 国产成人精品一,二区| 成人综合一区亚洲| 午夜久久久久精精品| 最近最新中文字幕大全电影3| 国产又黄又爽又无遮挡在线| 国产大屁股一区二区在线视频| 国产91av在线免费观看| 色视频www国产| 18禁在线无遮挡免费观看视频| 中文精品一卡2卡3卡4更新| 亚洲av福利一区| 91aial.com中文字幕在线观看| 精华霜和精华液先用哪个| 精品国产一区二区三区久久久樱花 | 亚洲一级一片aⅴ在线观看| 国产高清视频在线观看网站| 日本猛色少妇xxxxx猛交久久| 麻豆国产97在线/欧美| 激情 狠狠 欧美| 成人漫画全彩无遮挡| 国产免费视频播放在线视频 | 好男人视频免费观看在线| 久久久久九九精品影院| 亚洲最大成人手机在线| 水蜜桃什么品种好| 午夜激情欧美在线| 国产乱来视频区| 国产精品伦人一区二区| 精品人妻一区二区三区麻豆| 成人综合一区亚洲| 欧美高清成人免费视频www| 午夜免费男女啪啪视频观看| 国产伦理片在线播放av一区| 麻豆成人av视频| 亚洲经典国产精华液单| 国产极品精品免费视频能看的| 婷婷色麻豆天堂久久 | 免费在线观看成人毛片| 晚上一个人看的免费电影| 亚洲精品乱码久久久v下载方式| 国产亚洲av片在线观看秒播厂 | 国产中年淑女户外野战色| 精品99又大又爽又粗少妇毛片| 男的添女的下面高潮视频| 亚洲成色77777| 欧美精品一区二区大全| 亚洲美女视频黄频| 国产黄a三级三级三级人| 天美传媒精品一区二区| 春色校园在线视频观看| 午夜爱爱视频在线播放| 蜜桃亚洲精品一区二区三区| 久久久久久久久久久丰满| 亚洲在线观看片| 97超碰精品成人国产| 亚洲在久久综合| 九九在线视频观看精品| 男的添女的下面高潮视频| 精品人妻偷拍中文字幕| 自拍偷自拍亚洲精品老妇| 久久亚洲精品不卡| 日韩三级伦理在线观看| 国产亚洲午夜精品一区二区久久 | 亚洲在线自拍视频| 午夜视频国产福利| 日本熟妇午夜| 精品国产一区二区三区久久久樱花 | 亚洲人与动物交配视频| a级毛色黄片| 又爽又黄a免费视频| 日本一本二区三区精品| 精品久久久噜噜| 亚洲成av人片在线播放无| 亚洲国产精品合色在线| 国产免费福利视频在线观看| 亚洲精品亚洲一区二区| 国产成人91sexporn| 一级二级三级毛片免费看| 一个人观看的视频www高清免费观看| 日韩亚洲欧美综合| 成人鲁丝片一二三区免费| 欧美高清性xxxxhd video| 国产探花在线观看一区二区| 精品久久久久久久人妻蜜臀av| 国产精品99久久久久久久久| 欧美区成人在线视频| 国产亚洲5aaaaa淫片| www日本黄色视频网| 国产精品一区二区三区四区久久| 久久人人爽人人爽人人片va| 一个人免费在线观看电影| 午夜福利视频1000在线观看| 久久精品国产亚洲网站| 高清毛片免费看| 免费人成在线观看视频色| 在线观看66精品国产| 22中文网久久字幕| 村上凉子中文字幕在线| 亚洲av男天堂| 激情 狠狠 欧美| 91精品伊人久久大香线蕉| 中文乱码字字幕精品一区二区三区 | 亚洲av男天堂| 精品久久国产蜜桃| av视频在线观看入口| 国产精华一区二区三区| 国产亚洲一区二区精品| 丝袜喷水一区| 久久亚洲国产成人精品v| 插逼视频在线观看| 久久久午夜欧美精品| 国产伦精品一区二区三区视频9| 精品久久久久久久久av| 精品免费久久久久久久清纯| 丰满少妇做爰视频| 老司机影院成人| 99热这里只有精品一区| 可以在线观看毛片的网站| 大香蕉久久网| 成人毛片a级毛片在线播放| eeuss影院久久| 国内揄拍国产精品人妻在线| 狂野欧美激情性xxxx在线观看| 亚洲婷婷狠狠爱综合网| 国产麻豆成人av免费视频| 又爽又黄a免费视频| 国产精品不卡视频一区二区| 日本免费在线观看一区| 国产av在哪里看| 国产精品人妻久久久影院| 我要搜黄色片| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩高清专用| 国产一级毛片在线| 欧美潮喷喷水| 波野结衣二区三区在线| 日韩精品青青久久久久久| 七月丁香在线播放| 91久久精品国产一区二区三区| 好男人视频免费观看在线| 秋霞在线观看毛片| 久久久精品大字幕| 色综合色国产| 国产极品天堂在线| 色吧在线观看| 久久午夜福利片| 国语对白做爰xxxⅹ性视频网站| 美女内射精品一级片tv| 少妇熟女欧美另类| 国产免费又黄又爽又色| 伦理电影大哥的女人| 国产一区二区三区av在线| 天堂√8在线中文| 亚洲精品影视一区二区三区av| 干丝袜人妻中文字幕| 联通29元200g的流量卡| 日本黄大片高清| 国产成人免费观看mmmm| 天天躁日日操中文字幕| 国产亚洲一区二区精品| 天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 看非洲黑人一级黄片| 亚洲最大成人av| 中文乱码字字幕精品一区二区三区 | 一级毛片久久久久久久久女| 午夜激情福利司机影院| 欧美97在线视频| 插逼视频在线观看| 日本wwww免费看| 校园人妻丝袜中文字幕| 色视频www国产| 欧美潮喷喷水| 国产成人精品久久久久久| 搡女人真爽免费视频火全软件| 最近中文字幕高清免费大全6| 三级国产精品片| 亚洲国产日韩欧美精品在线观看| 91久久精品电影网| 色吧在线观看| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 一个人看视频在线观看www免费| 日韩av不卡免费在线播放| 黑人高潮一二区| 亚洲五月天丁香| 女的被弄到高潮叫床怎么办| 欧美人与善性xxx| 只有这里有精品99| a级毛片免费高清观看在线播放| 老司机影院成人| 十八禁国产超污无遮挡网站| 在线观看66精品国产| eeuss影院久久| 狠狠狠狠99中文字幕| 建设人人有责人人尽责人人享有的 | 天美传媒精品一区二区| 国产精品国产三级国产专区5o | 在线免费观看的www视频| 国产精品久久久久久av不卡| 91精品国产九色| av在线亚洲专区| 日韩一本色道免费dvd| 天堂√8在线中文| 美女cb高潮喷水在线观看| 亚洲久久久久久中文字幕| 看免费成人av毛片| 国内精品美女久久久久久| 亚洲国产精品专区欧美| videossex国产| 国产大屁股一区二区在线视频| 少妇熟女aⅴ在线视频| 五月伊人婷婷丁香| 少妇裸体淫交视频免费看高清| 成人性生交大片免费视频hd| 啦啦啦观看免费观看视频高清| 久久久久久久久久久免费av| 日韩中字成人| 免费观看性生交大片5| 又爽又黄无遮挡网站| 亚洲欧美成人精品一区二区| 国产爱豆传媒在线观看| 国产精品人妻久久久影院| 亚洲av免费高清在线观看| 国产在视频线在精品| 三级经典国产精品| 91av网一区二区| 视频中文字幕在线观看| 99热全是精品| 国产激情偷乱视频一区二区| 国产成人精品久久久久久| 五月玫瑰六月丁香| 99久久精品国产国产毛片| 亚洲精品一区蜜桃| 国语自产精品视频在线第100页| 少妇丰满av| 久久6这里有精品| 一卡2卡三卡四卡精品乱码亚洲| 欧美+日韩+精品| 黄片无遮挡物在线观看| 久久综合国产亚洲精品| 在线免费观看不下载黄p国产| 久久这里有精品视频免费| 日本-黄色视频高清免费观看| 97超视频在线观看视频| 禁无遮挡网站| 国产高清三级在线| 国产av码专区亚洲av| 成人国产麻豆网| 国产精品久久视频播放| 亚洲国产精品久久男人天堂| 一级黄色大片毛片| 免费av毛片视频| 亚洲国产最新在线播放| 国产在视频线在精品| av.在线天堂| 国产精品三级大全| 国产精品久久视频播放| 国产精品久久电影中文字幕| 亚洲国产精品合色在线| 男人舔女人下体高潮全视频| 一个人观看的视频www高清免费观看| 亚洲国产精品国产精品| 美女脱内裤让男人舔精品视频| 99久久精品热视频| 最近手机中文字幕大全| 大又大粗又爽又黄少妇毛片口| 人妻夜夜爽99麻豆av| 国国产精品蜜臀av免费| 中文字幕精品亚洲无线码一区| 久久亚洲国产成人精品v| 欧美又色又爽又黄视频| 国产精品精品国产色婷婷| 蜜桃亚洲精品一区二区三区| 国产精品国产高清国产av| 搡女人真爽免费视频火全软件| 18禁裸乳无遮挡免费网站照片| 高清在线视频一区二区三区 | 国内揄拍国产精品人妻在线| 国产成人a∨麻豆精品| 高清在线视频一区二区三区 | 国产成人精品婷婷| 国产真实乱freesex| 91精品一卡2卡3卡4卡| 亚洲国产精品sss在线观看| 免费观看性生交大片5| 久久久久久久国产电影| 夜夜看夜夜爽夜夜摸| av又黄又爽大尺度在线免费看 | 99热精品在线国产| 男女啪啪激烈高潮av片| 国产精品爽爽va在线观看网站| 人妻系列 视频| 亚洲在久久综合| 日韩大片免费观看网站 | 亚洲天堂国产精品一区在线| 成年女人看的毛片在线观看| 国产成人精品婷婷| 非洲黑人性xxxx精品又粗又长| 精品少妇黑人巨大在线播放 | 99久久无色码亚洲精品果冻| 丝袜喷水一区| 精品国产露脸久久av麻豆 | 国产人妻一区二区三区在| 夜夜爽夜夜爽视频| 人妻夜夜爽99麻豆av| 亚洲婷婷狠狠爱综合网| 啦啦啦韩国在线观看视频| 亚洲精品成人久久久久久| 亚洲成人久久爱视频| 可以在线观看毛片的网站| 日韩欧美三级三区| 色尼玛亚洲综合影院| 午夜日本视频在线| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 舔av片在线| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| 1024手机看黄色片| 内射极品少妇av片p| 精品一区二区免费观看| 国产精品久久久久久精品电影小说 | 午夜福利视频1000在线观看| 国产成人午夜福利电影在线观看| 久久精品久久久久久噜噜老黄 | 黄色一级大片看看| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 岛国毛片在线播放| 国产av不卡久久| 18禁动态无遮挡网站| 精品熟女少妇av免费看| 欧美性感艳星| av黄色大香蕉| 日韩精品青青久久久久久| 国产免费视频播放在线视频 | 日本三级黄在线观看| 精品人妻偷拍中文字幕| 99国产精品一区二区蜜桃av| 青春草视频在线免费观看| 亚洲欧美日韩无卡精品| 一边摸一边抽搐一进一小说| 秋霞伦理黄片| 久久99热这里只有精品18| 三级国产精品片| 插逼视频在线观看| 国产视频首页在线观看| 内地一区二区视频在线| 午夜福利高清视频| 亚洲av熟女| 国产成人a区在线观看| av在线蜜桃| 久久韩国三级中文字幕| 干丝袜人妻中文字幕| 国产午夜福利久久久久久| 观看美女的网站| 麻豆一二三区av精品| 我的女老师完整版在线观看| 亚洲精品日韩av片在线观看| 国产成人福利小说| 午夜激情欧美在线| 99视频精品全部免费 在线| 美女xxoo啪啪120秒动态图| 国产精品久久电影中文字幕| 国产免费一级a男人的天堂| 狠狠狠狠99中文字幕| 日韩视频在线欧美| 3wmmmm亚洲av在线观看| 成人性生交大片免费视频hd| 女人久久www免费人成看片 | 免费在线观看成人毛片| 老司机影院毛片| 国产亚洲最大av| 九九爱精品视频在线观看| 亚洲精品乱码久久久久久按摩| 精品久久久久久久末码| 免费av观看视频| 久久99热这里只有精品18| 亚洲一区高清亚洲精品| 精品一区二区三区人妻视频| 国产高清有码在线观看视频| 成人午夜精彩视频在线观看| 亚洲综合色惰| 亚洲av电影不卡..在线观看| 亚洲最大成人中文| 国语自产精品视频在线第100页| 变态另类丝袜制服| 精品久久久久久久久久久久久| 国产不卡一卡二| 欧美丝袜亚洲另类| 中文天堂在线官网| 麻豆成人午夜福利视频| 全区人妻精品视频| 视频中文字幕在线观看| 天堂网av新在线| 欧美日韩精品成人综合77777| 日韩 亚洲 欧美在线| 国产激情偷乱视频一区二区| 免费看光身美女| 精品久久国产蜜桃| 亚洲不卡免费看| 日产精品乱码卡一卡2卡三| 日韩欧美精品免费久久| 国产一区二区在线av高清观看| 精品久久久久久久末码| 麻豆成人av视频| 日本猛色少妇xxxxx猛交久久| 尤物成人国产欧美一区二区三区| 午夜日本视频在线| 国产亚洲精品久久久com| 三级经典国产精品| 18禁在线播放成人免费| 亚洲三级黄色毛片| 久久这里只有精品中国| 国产69精品久久久久777片| 国内少妇人妻偷人精品xxx网站| 一本久久精品| 国产成人freesex在线| 成年女人永久免费观看视频| 国产一区有黄有色的免费视频 | 日韩欧美国产在线观看| 最新中文字幕久久久久| 三级毛片av免费| 亚洲第一区二区三区不卡| 亚洲婷婷狠狠爱综合网| 国产成人freesex在线| 欧美性感艳星| 亚洲丝袜综合中文字幕| 日韩欧美三级三区| 亚洲精品乱码久久久v下载方式| 97人妻精品一区二区三区麻豆| 99热6这里只有精品| 国产免费福利视频在线观看| 亚洲高清免费不卡视频| 亚洲国产精品合色在线| 我要搜黄色片| 18+在线观看网站| 有码 亚洲区| 中文字幕久久专区| 国产私拍福利视频在线观看| 亚洲国产精品久久男人天堂| av国产免费在线观看| 少妇熟女欧美另类| 欧美精品一区二区大全| 99久久精品国产国产毛片| 波野结衣二区三区在线| 卡戴珊不雅视频在线播放| 内射极品少妇av片p| 嫩草影院精品99| 小蜜桃在线观看免费完整版高清| 亚洲国产欧洲综合997久久,| 国产大屁股一区二区在线视频| 亚洲av男天堂| 国产久久久一区二区三区| 欧美日韩精品成人综合77777| 成年免费大片在线观看| 日韩av不卡免费在线播放| 亚洲精品国产av成人精品| .国产精品久久| 国产精品女同一区二区软件| 天堂av国产一区二区熟女人妻| 黑人高潮一二区| 国产亚洲av片在线观看秒播厂 | 亚洲中文字幕一区二区三区有码在线看| 卡戴珊不雅视频在线播放| 成人高潮视频无遮挡免费网站| av在线播放精品| 亚洲在久久综合| 国产av不卡久久| 久久久精品大字幕| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 亚洲国产精品sss在线观看| 麻豆av噜噜一区二区三区| 亚洲av中文av极速乱| 国产精品美女特级片免费视频播放器| 成年女人看的毛片在线观看| 一区二区三区乱码不卡18| 国产高清不卡午夜福利| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 亚洲av福利一区| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 国产一区有黄有色的免费视频 | 1000部很黄的大片| 欧美成人一区二区免费高清观看| 网址你懂的国产日韩在线| 欧美zozozo另类| 国产探花极品一区二区| 91久久精品国产一区二区三区| 大又大粗又爽又黄少妇毛片口| 亚洲一级一片aⅴ在线观看| av在线老鸭窝| 美女黄网站色视频| 日韩欧美国产在线观看| 精品久久久久久电影网 | 亚洲国产精品sss在线观看| 欧美成人免费av一区二区三区| 免费观看a级毛片全部| 午夜激情福利司机影院| 又爽又黄a免费视频| 哪个播放器可以免费观看大片| 欧美成人精品欧美一级黄| 亚洲精品日韩av片在线观看| 晚上一个人看的免费电影| 亚洲怡红院男人天堂| 真实男女啪啪啪动态图| 少妇的逼好多水| 国产精品不卡视频一区二区| 精品免费久久久久久久清纯| 亚洲成色77777| 亚洲精品国产av成人精品| 两个人的视频大全免费| 最近的中文字幕免费完整| 久久99精品国语久久久| 美女高潮的动态| 黄片wwwwww| 色尼玛亚洲综合影院| 欧美日本视频| 日日撸夜夜添| 久久精品国产鲁丝片午夜精品| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片| 男人舔女人下体高潮全视频| 日本三级黄在线观看| 亚洲精品乱久久久久久| 国产精品麻豆人妻色哟哟久久 | 欧美日本亚洲视频在线播放| 国产视频内射| 国产av在哪里看| 国产精品蜜桃在线观看| 精品国产三级普通话版| 蜜桃久久精品国产亚洲av| 国产成人91sexporn| 亚洲av成人精品一区久久| 高清av免费在线| 国产又色又爽无遮挡免| 我要看日韩黄色一级片| 天天一区二区日本电影三级| 毛片一级片免费看久久久久| av天堂中文字幕网| 18禁裸乳无遮挡免费网站照片| 亚洲中文字幕日韩| 在线免费十八禁| 久久久午夜欧美精品| 麻豆av噜噜一区二区三区| 美女大奶头视频| 亚洲国产精品久久男人天堂| 免费av毛片视频| 国产探花在线观看一区二区| 成人特级av手机在线观看| 亚洲无线观看免费| 岛国在线免费视频观看| 网址你懂的国产日韩在线| 亚洲国产最新在线播放| 全区人妻精品视频| 午夜亚洲福利在线播放| 特级一级黄色大片| 免费观看人在逋| 久久久久久久久久成人| 九九热线精品视视频播放| 日本一本二区三区精品| 久久这里只有精品中国| 亚洲国产成人一精品久久久| 亚洲精品,欧美精品| 最近的中文字幕免费完整| 国内少妇人妻偷人精品xxx网站| 美女cb高潮喷水在线观看|