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

    基于多目標(biāo)遺傳算法的灌排兩用渠道輸水優(yōu)化調(diào)度

    2020-12-17 08:48:44徐淑琴王雅君喬厚清郭曉婷李仲裕徐恩典
    關(guān)鍵詞:灌排支渠兩用

    徐淑琴,王雅君,喬厚清,郭曉婷,李仲裕,徐恩典

    (東北農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,哈爾濱 150030)

    黑龍江省三江平原地區(qū)為黑龍江省糧食主產(chǎn)區(qū),也是我國(guó)糧食戰(zhàn)略生產(chǎn)和儲(chǔ)備重要基地,糧食商品率一直保持80%以上,多年來(lái)為保證國(guó)家糧食安全貢獻(xiàn)顯著[1]。近年來(lái),黑龍江省水田面積發(fā)展較快,水源以地下水灌溉為主,灌溉用水保障程度低,發(fā)展以地表水源為主的大、中型灌區(qū),將現(xiàn)有排水溝道改造成為灌溉與排水兩用渠道是保證農(nóng)業(yè)可持續(xù)發(fā)展的重要前提[2]。但目前缺乏水資源統(tǒng)一規(guī)劃,尤其是灌區(qū)內(nèi)灌溉用水調(diào)度研究進(jìn)展緩慢,豐富水資源未得到充分利用,需科學(xué)決策并優(yōu)化調(diào)度灌排兩用渠道。因此灌溉、排水系統(tǒng)合一,需制定嚴(yán)謹(jǐn)排澇與灌溉調(diào)度運(yùn)用方案,同時(shí)滿(mǎn)足灌溉與排澇功能,提高灌溉水利用率。趙恩龍利用多目標(biāo)遺傳算法,在水量平衡、庫(kù)容約束基礎(chǔ)上,運(yùn)用Matlab求解模型分別計(jì)算在75%保證率下不同單目標(biāo)優(yōu)化結(jié)果以及多目標(biāo)優(yōu)化結(jié)果,并作優(yōu)化對(duì)比分析[3]。李茉應(yīng)用線性規(guī)劃優(yōu)化灌區(qū)干支渠輸水時(shí)間流量,在此基礎(chǔ)上運(yùn)用0-1線性規(guī)劃法優(yōu)化灌區(qū)渠系輪灌組劃分,實(shí)現(xiàn)灌區(qū)渠系水土資源優(yōu)化配置[4]。梁貞堂等研究黑龍江灌區(qū)灌排兩用渠道斷面設(shè)計(jì)參數(shù),適當(dāng)調(diào)度后灌排兩用渠道可滿(mǎn)足灌溉與排水需要[5]。因此,為提高渠系水灌溉利用率,實(shí)現(xiàn)水資源利用最大化,保持農(nóng)業(yè)灌溉平衡穩(wěn)定發(fā)展,對(duì)灌排兩用渠道作輸水優(yōu)化調(diào)度尤為重要。

    1 研究區(qū)域概況

    1.1 灌區(qū)概況

    以黑龍江省楊木灌區(qū)為例,位于慶豐農(nóng)場(chǎng)中北部,主要以水田灌溉為主,灌區(qū)海拔高度60~65 m。本區(qū)年均降雨量為560 mm,年內(nèi)降水分配不勻,大多集中在7~9月,占全年降雨量一半以上。年均氣溫為2.6℃,年蒸發(fā)量為1 240 mm,積溫為2 700℃,全年日照約為2 500 h,無(wú)霜期為135~156 d,冰凍期一般在11月5~10日。

    1.2 渠道布置概況

    由于地表水置換地下水灌溉需要,楊木灌區(qū)主要以興凱湖水為水源,楊木灌區(qū)總干渠1條,長(zhǎng)17.27 km,分干渠7條,在大穆棱河上設(shè)倒虹吸一座,楊木總干渠以新河總干末端為起點(diǎn),水源穿過(guò)大穆棱河倒虹吸進(jìn)入主干渠,主灌渠由南向北從灌區(qū)中部穿過(guò),通過(guò)橋、涵閘等控制性建筑物進(jìn)入東西兩側(cè)7條分干渠。全區(qū)骨干渠道8條,長(zhǎng)52.46 km,骨干建筑物15座。

    本文研究區(qū)域?yàn)闂钅竟鄥^(qū)四分干渠0+000~8+150段灌排兩用渠道,四分干下布設(shè)3個(gè)支渠、3個(gè)排干,灌溉與排水流向同向。灌區(qū)渠道布置情況如圖1所示。

    圖1 各級(jí)渠道布置Fig.1 Canal layout at all the levels

    2 研究區(qū)域渠道流量及水位分析

    2.1 流量分析

    渠道灌溉時(shí),骨干渠道配水方式采用續(xù)灌,根據(jù)各支渠設(shè)計(jì)流量推求干渠各段設(shè)計(jì)流量,渠道排水時(shí),根據(jù)《興凱湖灌區(qū)初步設(shè)計(jì)》,農(nóng)田排水標(biāo)準(zhǔn)為五年一遇,即灌區(qū)3 d暴雨,4 d排到水田耐淹深度,水田排水模數(shù)為P20%=0.12 m3·km-2。見(jiàn)表1。

    2.2 水位分析

    灌排兩用渠道合理運(yùn)用必須考慮上下級(jí)水位關(guān)系,既滿(mǎn)足上下級(jí)排水水位要求,又滿(mǎn)足上下級(jí)灌溉水位要求[6]。根據(jù)公式(1)計(jì)算溝道設(shè)計(jì)排水時(shí)要求水位,再與實(shí)際灌排兩用渠道(以四分干為例)水位對(duì)比分析,設(shè)計(jì)水位抬高時(shí)對(duì)排澇效益的影響。灌排兩用渠道水位見(jiàn)表2。

    式中,H-溝道要求水位(m);H0-地面參考點(diǎn)高程(m);∑Li-各級(jí)溝道沿程水頭損失之和(m);∑σ-流過(guò)各級(jí)渠道建筑物水頭損失之和;h-低于田塊地面水深(m)。

    2.3 排澇效益分析

    灌排兩用渠道既滿(mǎn)足灌溉要求又滿(mǎn)足除澇要求。灌區(qū)澇害主要為暴雨所致,農(nóng)作物受淹時(shí)間過(guò)長(zhǎng)或深度過(guò)大均引起作物異常生長(zhǎng),輕者減產(chǎn),重者死亡[7]。

    在實(shí)際應(yīng)用中,作物生育期灌水期間,灌溉水位一般高于排水水位,突發(fā)暴雨時(shí)若不及時(shí)開(kāi)啟干渠排水口節(jié)制閘,田間受澇時(shí)部分水量無(wú)法排出,降水滯留使作物受澇時(shí)間過(guò)長(zhǎng),排澇效益折減。文章基于排澇要求水位以及渠道設(shè)計(jì)水位,分析灌溉水位升高對(duì)排澇面積的影響。見(jiàn)圖2。

    3 灌排兩用渠道輸水優(yōu)化調(diào)度及求解算法

    3.1 建模思想

    灌排兩用渠道優(yōu)化輸水實(shí)質(zhì)為優(yōu)化渠道輸水流量、輸水時(shí)間、灌水時(shí)閘門(mén)開(kāi)啟高度以及排水時(shí)閘門(mén)開(kāi)起高度。通過(guò)優(yōu)化輸水方案,實(shí)現(xiàn)灌區(qū)合理灌水和排水優(yōu)化目標(biāo),提高灌區(qū)水資源利用效率及灌區(qū)管理效益[8]。

    表2 灌排兩用渠道各項(xiàng)水位Table 2 Irrigation and drainage canal water level

    圖2 基于水位的排澇面積變化Fig.2 Drainage area variation based on water level

    以楊木灌區(qū)灌排兩用渠道為研究對(duì)象,以輸水流量、輸水時(shí)間以及灌排凈效益為優(yōu)化目標(biāo)。干渠、支渠為續(xù)灌渠道且模型中水量分布應(yīng)滿(mǎn)足以下條件:

    ①實(shí)際輸水應(yīng)滿(mǎn)足灌溉水量需求;②輸水時(shí)間需在灌水輪期內(nèi)變化;③干渠灌溉水位應(yīng)在排水要求水位和滿(mǎn)足自流灌溉水位之間;④渠道輸水時(shí),下級(jí)渠道總流量不大于上級(jí)渠道流量,防止因輸水過(guò)多造成渠道損壞。

    3.2 目標(biāo)函數(shù)及約束條件

    建立在輸水時(shí)間段T內(nèi)渠道水量損失最小、在不同時(shí)刻分干渠凈配水流量變化性最小、灌排凈效益最大目標(biāo)函數(shù)如下:

    目標(biāo)函數(shù)一:

    目標(biāo)函數(shù)二:

    目標(biāo)函數(shù)三:

    式中,Wd、We分別為分干渠、支渠輸水損失量(萬(wàn)m3);N為渠道條數(shù);j為輸水渠道編號(hào);lj為第j渠道長(zhǎng)度(m);qj為第j渠道實(shí)際配水流量(m3·s-1);tj為第j渠道配水時(shí)間(d),且tj1=tj2-tj;A、m分別為渠床土壤透水系數(shù)、渠床土壤透水指數(shù),參照文獻(xiàn)[7]上述參數(shù)取值分別為2.65、0.45;i為渠系輸水不同時(shí)段;T為時(shí)間,d;Q1j代表干渠不同時(shí)刻渠系輸水流量(m3·s-1),Q1j=qif(t),f(t)為0-1變量,代表灌區(qū)內(nèi)渠道輸水狀態(tài),渠道輸水時(shí),f(t)=1,反之則為0;Q1代表上級(jí)渠道各時(shí)段平均流量(m3·s-1)。C1、C2、C3-排澇效益、灌溉效益、支出費(fèi)用;Pw-水泵軸功率(w)且,H為揚(yáng)程(m),η、η01為水泵提水效率、電機(jī)工作效率,取η0=81.2%、η1=91%;T1為工作時(shí)間(d);B為電費(fèi)單價(jià)0.51(元·(kW·h)-1);α1為排澇面積折減系數(shù),取α1=1-(0.75Hs-46.11);k1、k2、ε1分別為工程分?jǐn)傁禂?shù)、治澇投資比例,取k1=14.5%、k2=10%、ε1=10%;S1、S2分別為排澇區(qū)水稻種植面積、水稻灌溉種植面積(hm2);X1、X2水稻單價(jià)為2.8(元·kg-1);A1、A2水稻單產(chǎn)為8696(kg·hm-2)。

    具體約束條件如下:

    流量約束:渠道輸水過(guò)程中,各渠道輸送流量需接近設(shè)計(jì)流量,且滿(mǎn)足輸水流量在渠道最小流量和最大流量之間。渠道流量最小系數(shù)αm、αd加大系數(shù)選取參照《農(nóng)田水利學(xué)》[9],Qd為渠道設(shè)計(jì)流量。即:

    時(shí)間約束:各級(jí)渠道配水開(kāi)始時(shí)間和結(jié)束時(shí)間應(yīng)在輪期內(nèi),且開(kāi)始時(shí)間應(yīng)大于0。即:

    水量平衡約束:下級(jí)渠道在任意時(shí)刻輸水流量之和均等于上級(jí)渠道輸水流量。即:

    水量約束:各級(jí)渠道配水流量與引水時(shí)間之積應(yīng)等于該渠道需配水量。即:

    灌溉水位約束:實(shí)際灌溉水位應(yīng)在排水水位和自流灌溉要求水位之間。即:

    3.3 求解算法

    多目標(biāo)傳統(tǒng)解法包括加權(quán)法、目標(biāo)規(guī)劃法、約束法等多種解法,但以往解法無(wú)法一次運(yùn)行得出問(wèn)題最優(yōu)解級(jí),也無(wú)法滿(mǎn)足實(shí)際應(yīng)用多種需求,且參數(shù)值設(shè)定具有主觀性,導(dǎo)致優(yōu)化結(jié)果較差。隨著以Pareto最優(yōu)理論為基礎(chǔ)多目標(biāo)進(jìn)化算法提出,該問(wèn)題得到較好解決[10]。

    遺傳算法(GA)為一種借鑒生物界自然選擇和群體遺傳學(xué)機(jī)制的隨機(jī)搜索和優(yōu)化算法,具有全局搜索能力,善于搜索復(fù)雜問(wèn)題和非線性問(wèn)題[11]。從代表問(wèn)題可能潛在解集的一個(gè)種群開(kāi)始,初代種群產(chǎn)生后,確定適應(yīng)度函數(shù),并借助自然遺傳學(xué)中遺傳算子選擇、交叉和變異,產(chǎn)生代表新解集種群,得出最優(yōu)解[12]。

    本文選用帶精英策略的非支配排序遺傳算法(Non-dominated sorting genetic algorithms-Ⅱ,NSGA-Ⅱ),該算法通過(guò)非支配排序、擁擠度與擁擠算子計(jì)算以及精英策略等降低運(yùn)算復(fù)雜度,提高運(yùn)算速度,應(yīng)用廣泛[13-14]。多目標(biāo)優(yōu)化問(wèn)題為目標(biāo)直接存在沖突無(wú)法比較,所有目標(biāo)函數(shù)難以同時(shí)得到最優(yōu)解集,且各解之間無(wú)法比較優(yōu)劣[15]。

    運(yùn)用該算法求解步驟如下:①首先,對(duì)模型構(gòu)建中決策變量作實(shí)數(shù)編碼,其中包括輸水流量、輸水時(shí)間、灌溉水位等。②隨機(jī)產(chǎn)生N個(gè)個(gè)體初始種群P0,非支配排序后,通過(guò)選擇、交叉、變異得到第一代子群。③第二代開(kāi)始將父代與子代結(jié)合,作快速非支配排序和擁擠度計(jì)算,再選取合適個(gè)體形成新父代種群。④運(yùn)用遺傳算法基本操作產(chǎn)生新一代子群,依次類(lèi)推,直到滿(mǎn)足條件為止。

    3.4 求解參數(shù)

    文章根據(jù)楊木灌區(qū)2008年蒸發(fā)資料制定灌溉制度,設(shè)計(jì)保證率為80%。研究區(qū)域灌溉面積1 767.55 hm2,排澇面積2 961.48 hm2,作物在移植返青、分蘗、拔節(jié)孕穗以及抽穗開(kāi)花階段總需水量為839.77萬(wàn)m3,灌區(qū)水源為地表水,充分灌溉。各生育階段灌水天數(shù)及灌水次數(shù)見(jiàn)表3,將以上數(shù)據(jù)代入NSGA-Ⅱ模型求解。

    4 實(shí)例結(jié)果分析

    據(jù)文章第3部分所構(gòu)建目標(biāo)函數(shù),運(yùn)用NSGA-Ⅱ求解,綜合考慮灌區(qū)排澇及灌溉效益,作物生育期灌水時(shí),選取最優(yōu)灌溉水位灌溉。灌排兩用渠道進(jìn)口灌溉水位優(yōu)化結(jié)果為61.67 m。

    作物各生育階段一支渠、二支渠、三支渠輸水開(kāi)始時(shí)間與結(jié)束時(shí)間見(jiàn)圖3,可知四分干渠在移植返青期、分蘗期、拔節(jié)孕穗期、抽穗開(kāi)花期4個(gè)階段輸水時(shí)間分別為4~19、8~29、1~13、1~9 d,且優(yōu)化后輸水時(shí)間結(jié)果與優(yōu)化前相比減少,因此輸水時(shí)各級(jí)渠道輸水平均流量接近實(shí)際輸水流量,減少渠道輸水損失。

    表3 作物各生育期灌水參數(shù)Table 3 Irrigation parameters of crops at different growth period

    圖3 生育期灌水持續(xù)時(shí)間Fig.3 Duration of irrigation during the growth period

    優(yōu)化后各級(jí)渠道輸水流量及流速見(jiàn)圖4,由此可見(jiàn),各級(jí)渠道輸水流量波動(dòng)性較小,且流速滿(mǎn)足渠道不沖、不淤流速(V不沖=0.9 m·s-1、V不淤=0.21 m·s-1)。通過(guò)遺傳算法優(yōu)勝劣汰處理,灌區(qū)一、二支渠控制灌溉面積大于三支渠控制灌溉面積,且一、二支渠輸送水量也大于三支渠輸送水量。

    為實(shí)現(xiàn)灌區(qū)灌排兩用渠道輸水調(diào)度自動(dòng)化處理,計(jì)算灌區(qū)在排澇和灌溉期間閘門(mén)開(kāi)啟高度。為控制渠道水位以滿(mǎn)足灌溉排水要求,楊木灌區(qū)四分干設(shè)有節(jié)制閘6座,其中用于灌溉3座,排澇3座。閘門(mén)均為單孔閘門(mén)且閘孔自由出流,閘門(mén)寬1.5 m,高1 m,根據(jù)2008年排澇標(biāo)準(zhǔn)以及優(yōu)化后求解得到各渠道輸水流量計(jì)算閘門(mén)開(kāi)度,得出灌區(qū)灌溉期間各生育階段以及排澇期各閘門(mén)開(kāi)啟高度見(jiàn)圖5,渠道排水及灌溉時(shí)閘門(mén)開(kāi)度均小于渠道水深,且灌排兩用渠道在實(shí)際運(yùn)行中各分水口水位保持穩(wěn)定,滿(mǎn)足下級(jí)渠道輸水流量。優(yōu)化結(jié)果匯總見(jiàn)表4。

    圖4 生育期渠道輸水流量和流速Fig.4 Canal water flow and velocity during the growth period

    圖5 各時(shí)期渠道閘門(mén)開(kāi)啟高度Fig.5 Opening height of channel gate at different periods

    表4 生育期渠道輸水流量和流速Fig.4 Canal water flow and velocity during the growth period

    續(xù)表

    5 結(jié) 論

    a.本文以興凱湖楊木灌區(qū)灌排兩用渠道為例,充分考慮灌溉水位對(duì)排澇、灌溉效益影響以及作物灌水時(shí)渠道輸水損失大、輸水時(shí)間長(zhǎng)等特點(diǎn),建立基于多目標(biāo)灌排兩用渠道輸水調(diào)度模型,運(yùn)用帶精英策略的非支配排序遺傳算法(NSGA-Ⅱ)以及Matlab編程對(duì)目標(biāo)函數(shù)求解。

    b.據(jù)文章第4部分優(yōu)化求解得出灌溉、排水時(shí)對(duì)應(yīng)閘門(mén)開(kāi)啟高度,可控制各分水口水位以及輸水流量。優(yōu)化后結(jié)果與優(yōu)化前相比,作物灌溉期間渠道輸水時(shí)間縮短,輸水過(guò)程中輸水流量均勻穩(wěn)定;作物在排澇期間,閘門(mén)可及時(shí)啟閉,保證渠道在作物內(nèi)將多余水量排出。

    c.文章運(yùn)用模型和算法在滿(mǎn)足灌溉和排水要求時(shí),確定灌排兩用渠道實(shí)際輸水方案,調(diào)節(jié)各渠道灌溉水位、輸水流量、輸水時(shí)間以及閘門(mén)開(kāi)啟高度,灌排兩用渠道在灌溉和排水時(shí)輸水過(guò)程效率更高,為具有灌排兩用渠道灌區(qū)輸水優(yōu)化調(diào)度問(wèn)題提供依據(jù)。由于時(shí)間問(wèn)題,研究使用已有排澇標(biāo)準(zhǔn)替代降雨預(yù)測(cè)結(jié)果計(jì)算排澇,未來(lái)應(yīng)針對(duì)不同灌區(qū)降雨預(yù)測(cè)結(jié)果,實(shí)現(xiàn)實(shí)時(shí)優(yōu)化調(diào)度灌溉排水兩用渠道,完善灌溉及排澇管理控制方案。

    猜你喜歡
    灌排支渠兩用
    晴雨天兩用鞋
    洱海西部灌排溝渠水質(zhì)特征及土地利用的影響
    現(xiàn)代化大型灌排泵站建設(shè)研究
    藥食兩用說(shuō)食姜
    高標(biāo)準(zhǔn)農(nóng)田建設(shè)中灌排溝渠施工關(guān)鍵技術(shù)探討
    以“肥”殺蟲(chóng) 藥肥兩用 一舉多得
    青州張莊灌區(qū)渠道設(shè)計(jì)流量計(jì)算分析
    山東水利(2018年9期)2018-10-12 03:23:14
    渠道流量及流量分段的計(jì)算原則和方法
    墨玉縣普恰克其鄉(xiāng)家依閘口東支渠抗凍脹設(shè)計(jì)
    安陽(yáng)市灌區(qū)支斗渠管道化改造可行性分析
    国产极品粉嫩免费观看在线| 国产日韩一区二区三区精品不卡| 老司机深夜福利视频在线观看| 黑人巨大精品欧美一区二区mp4| 97人妻精品一区二区三区麻豆 | 成人国语在线视频| 18禁黄网站禁片午夜丰满| 日本vs欧美在线观看视频| 三级毛片av免费| 美女 人体艺术 gogo| 欧美另类亚洲清纯唯美| 欧美国产精品va在线观看不卡| 亚洲精品在线观看二区| 欧美日本中文国产一区发布| 久久精品人人爽人人爽视色| 国产精品av久久久久免费| 91成年电影在线观看| 黄色视频,在线免费观看| 一个人免费在线观看的高清视频| 国产一区二区三区综合在线观看| 19禁男女啪啪无遮挡网站| 成在线人永久免费视频| 免费女性裸体啪啪无遮挡网站| 日本免费a在线| 国产成人精品久久二区二区91| 久久久久国内视频| 99香蕉大伊视频| www国产在线视频色| 级片在线观看| 国产精品野战在线观看| 欧美乱色亚洲激情| 国产精品一区二区精品视频观看| 黄色女人牲交| 在线观看免费日韩欧美大片| 国产精品电影一区二区三区| 97碰自拍视频| 成人亚洲精品一区在线观看| 国产精品亚洲av一区麻豆| 国产亚洲精品综合一区在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 女人高潮潮喷娇喘18禁视频| or卡值多少钱| 欧美黄色淫秽网站| 一边摸一边抽搐一进一出视频| 亚洲国产精品sss在线观看| videosex国产| 在线av久久热| www日本在线高清视频| 色在线成人网| 国产亚洲精品久久久久5区| 欧美色欧美亚洲另类二区 | 国产精品久久久久久人妻精品电影| 可以免费在线观看a视频的电影网站| 搡老熟女国产l中国老女人| 乱人伦中国视频| 美女高潮到喷水免费观看| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| 激情视频va一区二区三区| 1024香蕉在线观看| 久久香蕉激情| 美国免费a级毛片| 国产野战对白在线观看| 丝袜美腿诱惑在线| 人妻久久中文字幕网| 亚洲最大成人中文| 婷婷六月久久综合丁香| 久久精品成人免费网站| 高潮久久久久久久久久久不卡| 精品国产亚洲在线| 国产精品九九99| 日韩有码中文字幕| 麻豆一二三区av精品| 天天躁夜夜躁狠狠躁躁| 日韩视频一区二区在线观看| 国产亚洲欧美精品永久| 亚洲电影在线观看av| 国产精品二区激情视频| 久久久精品国产亚洲av高清涩受| 欧美性长视频在线观看| 午夜免费观看网址| 在线观看免费午夜福利视频| 色综合婷婷激情| av中文乱码字幕在线| 免费看美女性在线毛片视频| 日本在线视频免费播放| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 人人妻人人澡人人看| 色播在线永久视频| av在线播放免费不卡| 久久久久久久久久久久大奶| 一级毛片女人18水好多| 村上凉子中文字幕在线| 99国产精品一区二区蜜桃av| 最新美女视频免费是黄的| 亚洲国产高清在线一区二区三 | 亚洲性夜色夜夜综合| 97超级碰碰碰精品色视频在线观看| 美国免费a级毛片| 日本欧美视频一区| 69av精品久久久久久| 国产av一区在线观看免费| 久久香蕉精品热| 国产欧美日韩精品亚洲av| 操出白浆在线播放| 757午夜福利合集在线观看| 日本 av在线| 久久精品成人免费网站| 99精品在免费线老司机午夜| 丁香欧美五月| 亚洲国产精品成人综合色| 国产精品 欧美亚洲| 一级黄色大片毛片| 亚洲中文av在线| 午夜老司机福利片| 中文字幕精品免费在线观看视频| 久久天堂一区二区三区四区| 国产成人av教育| 亚洲一区二区三区不卡视频| 精品一品国产午夜福利视频| 日韩大码丰满熟妇| 日日夜夜操网爽| av中文乱码字幕在线| 欧美成人性av电影在线观看| 久久久久久久久免费视频了| 亚洲美女黄片视频| 欧美成人一区二区免费高清观看 | 90打野战视频偷拍视频| 日韩欧美在线二视频| 亚洲欧美激情在线| 国产成人精品在线电影| 欧美国产日韩亚洲一区| 制服人妻中文乱码| 久久人人97超碰香蕉20202| 俄罗斯特黄特色一大片| 国产日韩一区二区三区精品不卡| 欧美+亚洲+日韩+国产| 免费观看人在逋| www.精华液| 日本黄色视频三级网站网址| 欧美日韩乱码在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一av免费看| 波多野结衣一区麻豆| 午夜日韩欧美国产| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 色综合婷婷激情| 精品国产一区二区三区四区第35| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 久久精品国产清高在天天线| 亚洲一区高清亚洲精品| 在线观看www视频免费| 女性生殖器流出的白浆| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 最近最新中文字幕大全免费视频| 国产精品一区二区免费欧美| 一级a爱片免费观看的视频| 午夜福利在线观看吧| 制服人妻中文乱码| 国产成人系列免费观看| 好男人在线观看高清免费视频 | 人人妻人人澡欧美一区二区 | 悠悠久久av| 久久人妻福利社区极品人妻图片| 精品久久久久久久久久免费视频| 搡老岳熟女国产| 十八禁人妻一区二区| 久久国产亚洲av麻豆专区| 国产在线观看jvid| 国产精品美女特级片免费视频播放器 | 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| 女性被躁到高潮视频| 日韩欧美三级三区| 中出人妻视频一区二区| 国产欧美日韩一区二区三区在线| 免费看十八禁软件| 午夜福利影视在线免费观看| 精品人妻在线不人妻| 久99久视频精品免费| 欧美绝顶高潮抽搐喷水| 久热这里只有精品99| 真人一进一出gif抽搐免费| av电影中文网址| 他把我摸到了高潮在线观看| 首页视频小说图片口味搜索| 老司机靠b影院| 免费观看人在逋| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 亚洲色图 男人天堂 中文字幕| 亚洲国产日韩欧美精品在线观看 | 中国美女看黄片| 一级毛片精品| 国产主播在线观看一区二区| 两人在一起打扑克的视频| 变态另类成人亚洲欧美熟女 | 亚洲国产精品合色在线| 91成年电影在线观看| 日本三级黄在线观看| 国产精品综合久久久久久久免费 | 亚洲午夜理论影院| 国产免费av片在线观看野外av| 国产国语露脸激情在线看| 亚洲午夜精品一区,二区,三区| 国产一区二区三区综合在线观看| 亚洲第一电影网av| 国产成人系列免费观看| 国产成人av教育| 国产精品免费视频内射| 18禁美女被吸乳视频| 极品教师在线免费播放| av片东京热男人的天堂| 老汉色∧v一级毛片| 亚洲七黄色美女视频| 正在播放国产对白刺激| 欧美在线黄色| 免费不卡黄色视频| 午夜视频精品福利| 99riav亚洲国产免费| 精品久久久精品久久久| 女生性感内裤真人,穿戴方法视频| 精品久久久久久成人av| 亚洲精品一区av在线观看| 久久久久国内视频| 亚洲午夜精品一区,二区,三区| 欧美激情极品国产一区二区三区| 国产精品一区二区在线不卡| 欧美在线黄色| av片东京热男人的天堂| 久久久久久免费高清国产稀缺| 制服丝袜大香蕉在线| 久久久久久久久免费视频了| 久久亚洲精品不卡| 欧美中文日本在线观看视频| 此物有八面人人有两片| 一区二区三区精品91| 超碰成人久久| 露出奶头的视频| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 精品一区二区三区四区五区乱码| 免费在线观看影片大全网站| 一本久久中文字幕| 天天一区二区日本电影三级 | 免费在线观看视频国产中文字幕亚洲| 国产一区二区激情短视频| 日本在线视频免费播放| 中文字幕人成人乱码亚洲影| 国产亚洲欧美98| 妹子高潮喷水视频| www.www免费av| 一本大道久久a久久精品| 岛国视频午夜一区免费看| 婷婷丁香在线五月| 欧美日本亚洲视频在线播放| bbb黄色大片| 一边摸一边做爽爽视频免费| 母亲3免费完整高清在线观看| 可以在线观看毛片的网站| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 亚洲人成电影观看| 九色亚洲精品在线播放| 久久精品国产亚洲av高清一级| 国产成人免费无遮挡视频| 国产蜜桃级精品一区二区三区| av片东京热男人的天堂| 黑人巨大精品欧美一区二区蜜桃| 麻豆国产av国片精品| 亚洲欧美精品综合一区二区三区| 久久久久久久精品吃奶| 搞女人的毛片| 久久中文字幕人妻熟女| 国产亚洲精品综合一区在线观看 | 日韩 欧美 亚洲 中文字幕| 1024视频免费在线观看| 国产99久久九九免费精品| 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女| 国产精品电影一区二区三区| 99热只有精品国产| 欧美大码av| 国产又色又爽无遮挡免费看| 大码成人一级视频| 亚洲欧美日韩无卡精品| 日本撒尿小便嘘嘘汇集6| 欧美亚洲日本最大视频资源| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 日韩高清综合在线| 精品少妇一区二区三区视频日本电影| 国产精品 国内视频| 国产亚洲精品一区二区www| 99精品欧美一区二区三区四区| 色av中文字幕| 久久国产精品男人的天堂亚洲| 免费在线观看亚洲国产| 日韩欧美国产一区二区入口| bbb黄色大片| 国产精品99久久99久久久不卡| 国产亚洲精品久久久久5区| 91精品国产国语对白视频| 宅男免费午夜| 国产av一区在线观看免费| 午夜福利欧美成人| 精品日产1卡2卡| 午夜福利18| 久久这里只有精品19| 欧美日韩黄片免| 中国美女看黄片| 又黄又爽又免费观看的视频| 欧美成人性av电影在线观看| 欧美精品啪啪一区二区三区| 满18在线观看网站| 亚洲一区中文字幕在线| 亚洲av熟女| 久久精品亚洲精品国产色婷小说| 中文字幕另类日韩欧美亚洲嫩草| 精品久久久精品久久久| 香蕉丝袜av| 丝袜美足系列| 国产成人精品在线电影| 亚洲午夜理论影院| 黑丝袜美女国产一区| 国产激情欧美一区二区| 国产黄a三级三级三级人| 人人澡人人妻人| 国产91精品成人一区二区三区| 精品久久久久久成人av| 久热这里只有精品99| 欧美中文日本在线观看视频| 丁香欧美五月| 人妻久久中文字幕网| 日韩成人在线观看一区二区三区| 狂野欧美激情性xxxx| 亚洲精品一区av在线观看| 美女 人体艺术 gogo| 999久久久国产精品视频| 午夜精品在线福利| 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| 国产视频一区二区在线看| 中亚洲国语对白在线视频| 成人国产综合亚洲| 久久久久久久久久久久大奶| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 免费无遮挡裸体视频| 日日干狠狠操夜夜爽| 又黄又粗又硬又大视频| 18禁美女被吸乳视频| 熟妇人妻久久中文字幕3abv| 51午夜福利影视在线观看| 欧美乱色亚洲激情| 最近最新中文字幕大全免费视频| 老熟妇乱子伦视频在线观看| 久久久久久久久久久久大奶| 麻豆成人av在线观看| 午夜福利在线观看吧| 国产精品免费视频内射| 久久久久九九精品影院| 久久久久久免费高清国产稀缺| 在线观看免费日韩欧美大片| 麻豆久久精品国产亚洲av| tocl精华| 午夜视频精品福利| 久久久精品国产亚洲av高清涩受| 国产亚洲精品第一综合不卡| 性欧美人与动物交配| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 亚洲国产精品sss在线观看| 欧美中文日本在线观看视频| 在线观看一区二区三区| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 色av中文字幕| 一进一出好大好爽视频| 淫妇啪啪啪对白视频| 91精品三级在线观看| 怎么达到女性高潮| 国产精品日韩av在线免费观看 | 日韩av在线大香蕉| 国产精品99久久99久久久不卡| 老鸭窝网址在线观看| 亚洲一区中文字幕在线| 一边摸一边抽搐一进一小说| 国产一区二区三区视频了| 波多野结衣一区麻豆| 少妇粗大呻吟视频| 老司机午夜十八禁免费视频| 午夜福利18| 久久精品国产99精品国产亚洲性色 | 日韩视频一区二区在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲黑人精品在线| 亚洲欧美日韩另类电影网站| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美日韩在线播放| 免费在线观看亚洲国产| 搡老岳熟女国产| 精品久久久久久成人av| 宅男免费午夜| 国产99久久九九免费精品| 亚洲中文av在线| 亚洲午夜理论影院| 满18在线观看网站| 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| ponron亚洲| 热99re8久久精品国产| 国产精品影院久久| 91麻豆av在线| 国产欧美日韩一区二区三| 久久人人爽av亚洲精品天堂| 国产亚洲欧美在线一区二区| 一本综合久久免费| 日韩av在线大香蕉| 一区二区三区激情视频| 亚洲人成电影免费在线| 免费看a级黄色片| www国产在线视频色| 午夜a级毛片| 日韩欧美国产在线观看| 亚洲精品美女久久av网站| 国产成人免费无遮挡视频| 成在线人永久免费视频| 亚洲午夜精品一区,二区,三区| 午夜福利在线观看吧| 久久精品91蜜桃| 久久精品91无色码中文字幕| 身体一侧抽搐| 非洲黑人性xxxx精品又粗又长| 色在线成人网| 国产精品 欧美亚洲| 国产精品免费一区二区三区在线| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 欧美不卡视频在线免费观看 | 亚洲av熟女| 亚洲五月色婷婷综合| 国产av又大| 伊人久久大香线蕉亚洲五| 在线av久久热| 他把我摸到了高潮在线观看| 国产精品电影一区二区三区| 国内毛片毛片毛片毛片毛片| 老熟妇乱子伦视频在线观看| 午夜精品在线福利| 午夜免费成人在线视频| 熟女少妇亚洲综合色aaa.| 久久中文看片网| 丁香六月欧美| а√天堂www在线а√下载| 一级a爱视频在线免费观看| 国产色视频综合| 日韩精品青青久久久久久| 成年女人毛片免费观看观看9| 国产主播在线观看一区二区| 一级毛片女人18水好多| 丝袜美腿诱惑在线| 变态另类成人亚洲欧美熟女 | 成人三级黄色视频| 亚洲天堂国产精品一区在线| 国产亚洲欧美在线一区二区| 欧美av亚洲av综合av国产av| 午夜久久久在线观看| 美女高潮喷水抽搐中文字幕| 精品卡一卡二卡四卡免费| 国产成人精品久久二区二区91| 亚洲熟女毛片儿| 黄色片一级片一级黄色片| 97人妻精品一区二区三区麻豆 | 国产成人免费无遮挡视频| 欧美日韩乱码在线| 午夜两性在线视频| x7x7x7水蜜桃| 一区二区日韩欧美中文字幕| 免费观看精品视频网站| 日韩欧美免费精品| 精品不卡国产一区二区三区| 国产极品粉嫩免费观看在线| 日韩 欧美 亚洲 中文字幕| 国产在线精品亚洲第一网站| 精品久久久久久久毛片微露脸| 一边摸一边抽搐一进一出视频| 黄片大片在线免费观看| av天堂在线播放| √禁漫天堂资源中文www| 天堂√8在线中文| 亚洲一区二区三区不卡视频| www.www免费av| 自线自在国产av| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品美女久久久久99蜜臀| 国产亚洲精品综合一区在线观看 | 91av网站免费观看| 在线观看免费视频日本深夜| 国产精品秋霞免费鲁丝片| 国产又色又爽无遮挡免费看| 免费在线观看日本一区| 51午夜福利影视在线观看| 激情视频va一区二区三区| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av香蕉五月| 人妻丰满熟妇av一区二区三区| 俄罗斯特黄特色一大片| 国产视频一区二区在线看| 非洲黑人性xxxx精品又粗又长| 午夜a级毛片| 亚洲 欧美一区二区三区| 18禁美女被吸乳视频| 精品一区二区三区视频在线观看免费| 久久久久久亚洲精品国产蜜桃av| 色综合亚洲欧美另类图片| 一级毛片精品| av网站免费在线观看视频| 亚洲国产欧美日韩在线播放| 婷婷丁香在线五月| 国产一区二区三区综合在线观看| 欧美日本视频| 男人的好看免费观看在线视频 | 女人被躁到高潮嗷嗷叫费观| 桃色一区二区三区在线观看| 欧美老熟妇乱子伦牲交| 91在线观看av| 久久久久九九精品影院| 久久人人精品亚洲av| 久久影院123| 叶爱在线成人免费视频播放| 在线观看免费视频日本深夜| 两人在一起打扑克的视频| 国产91精品成人一区二区三区| 久久精品国产综合久久久| 久久精品91无色码中文字幕| 久久久久亚洲av毛片大全| 伦理电影免费视频| 亚洲性夜色夜夜综合| 一边摸一边做爽爽视频免费| 夜夜夜夜夜久久久久| 国产91精品成人一区二区三区| 精品久久久精品久久久| 久久性视频一级片| 丁香欧美五月| 窝窝影院91人妻| 久久久久久久精品吃奶| 亚洲 国产 在线| 亚洲国产精品久久男人天堂| 露出奶头的视频| 女同久久另类99精品国产91| 十分钟在线观看高清视频www| 亚洲激情在线av| 国产亚洲精品一区二区www| 禁无遮挡网站| 巨乳人妻的诱惑在线观看| 啦啦啦韩国在线观看视频| 色播亚洲综合网| 日韩精品青青久久久久久| 级片在线观看| 日本免费a在线| 悠悠久久av| 午夜久久久在线观看| 淫秽高清视频在线观看| 777久久人妻少妇嫩草av网站| 国产麻豆成人av免费视频| av视频免费观看在线观看| 9191精品国产免费久久| 午夜福利视频1000在线观看 | xxx96com| 日日干狠狠操夜夜爽| 日本免费一区二区三区高清不卡 | 啦啦啦 在线观看视频| 午夜激情av网站| aaaaa片日本免费| 亚洲成人国产一区在线观看| 99国产精品一区二区三区| 精品熟女少妇八av免费久了| 久热这里只有精品99| 韩国av一区二区三区四区| 国产亚洲av嫩草精品影院| 国产一区在线观看成人免费| 波多野结衣巨乳人妻| 国产成人影院久久av| 亚洲国产精品sss在线观看| 成人欧美大片| 午夜福利在线观看吧| 久久国产乱子伦精品免费另类| 怎么达到女性高潮| 一级片免费观看大全| 国产精品影院久久| 91老司机精品| 亚洲黑人精品在线| 女人高潮潮喷娇喘18禁视频| 亚洲欧美日韩高清在线视频| 大型av网站在线播放| 亚洲精品在线观看二区| 亚洲在线自拍视频| 啪啪无遮挡十八禁网站| ponron亚洲| 久久午夜综合久久蜜桃| а√天堂www在线а√下载| 男人舔女人的私密视频| 最近最新中文字幕大全电影3 | 亚洲五月婷婷丁香| 中文字幕人成人乱码亚洲影| 亚洲午夜理论影院| 真人一进一出gif抽搐免费| 午夜福利一区二区在线看| 又黄又粗又硬又大视频|