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

    基于NSGA-II算法的灌區(qū)配水渠道輪灌分組優(yōu)化調(diào)控研究

    2023-03-11 05:40:28樊煜高占義陳皓銳王云輝房本巖徐寧
    灌溉排水學(xué)報(bào) 2023年2期
    關(guān)鍵詞:斗渠分組灌溉

    樊煜,高占義,陳皓銳*,王云輝,房本巖,徐寧

    ?灌溉水源與輸配水系統(tǒng)?

    基于NSGA-II算法的灌區(qū)配水渠道輪灌分組優(yōu)化調(diào)控研究

    樊煜1,2,高占義1,2,陳皓銳1,2*,王云輝3,房本巖3,徐寧3

    (1.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100038;2.國(guó)家節(jié)水灌溉北京工程技術(shù)研究中心,北京 100048;3.山東省濱州市引黃灌溉服務(wù)中心,山東 濱州 256600)

    【目的】基于NSGA-II算法對(duì)灌區(qū)配水渠道輪灌分組進(jìn)行優(yōu)化調(diào)控,實(shí)現(xiàn)灌區(qū)高效配水?!痉椒ā恳韵录?jí)渠道流量和輪灌分組數(shù)為決策變量,以總配水時(shí)間短、同一輪灌組內(nèi)最上游渠道與最下游渠道間距小、同一輪灌組內(nèi)各渠道配水時(shí)間差最小為目標(biāo),考慮輪期約束、出水口狀態(tài)約束、渠道過(guò)流能力約束、水量平衡約束,建立多目標(biāo)配水優(yōu)化模型,利用NSGA-II算法對(duì)模型進(jìn)行求解?!窘Y(jié)果】對(duì)于簸箕李灌區(qū)的石皮東支渠,模型將斗渠劃分為9個(gè)輪灌組,模型的配水流量和配水總時(shí)長(zhǎng)均滿足約束要求;模型盡可能地將相鄰的渠道劃分到一個(gè)輪灌組并同時(shí)啟閉,方便管理人員操作閘門(mén)?!窘Y(jié)論】本文建立的灌區(qū)配水模型能夠?qū)崿F(xiàn)流量大、短歷時(shí)、配水工作量少,可為灌區(qū)配水工作提供技術(shù)支撐。

    多目標(biāo)配水;渠道配水;輪灌;優(yōu)化模型

    0 引言

    【研究意義】灌溉渠道水量?jī)?yōu)化分配是指在滿足作物灌溉需求的前提下,基于灌區(qū)各級(jí)渠道的現(xiàn)狀輸水能力,通過(guò)各種技術(shù)手段,對(duì)配水渠道的流量和配水時(shí)間進(jìn)行優(yōu)化,以達(dá)到設(shè)定目標(biāo)[1]。在灌區(qū)輸配水過(guò)程中,科學(xué)合理的渠系配水可以提高供水服務(wù)水平、減少無(wú)效棄水、提高灌溉水利用率,從而促進(jìn)灌區(qū)水資源高效利用、有效緩解水資源供需矛盾、實(shí)現(xiàn)灌區(qū)農(nóng)業(yè)生產(chǎn)優(yōu)質(zhì)高效發(fā)展,為灌區(qū)生產(chǎn)帶來(lái)良好的經(jīng)濟(jì)與生態(tài)環(huán)境效益[2-5]。

    【研究進(jìn)展】在灌溉渠系配水過(guò)程中,上級(jí)渠道往往不能滿足所有下級(jí)渠道的灌溉,灌區(qū)常常采用輪灌配水方式。Suryavanshi等[6]提出了0-1線性規(guī)劃模型,假定上級(jí)配水渠道是由一組等流量的流管組成,且下級(jí)渠道流量相等,通過(guò)尋找下級(jí)渠道的最優(yōu)組合方式來(lái)減少渠道工程投資,獲得灌溉配水的最優(yōu)調(diào)度,這種輪灌分組方式已得到了廣泛應(yīng)用。王智等[7]提出了在一定輪灌期和來(lái)水流量限制條件下實(shí)現(xiàn)下級(jí)配水渠道最優(yōu)組合及配水時(shí)間最優(yōu)的配水模型;李彬等[8]將輪灌最優(yōu)組合方案與計(jì)算機(jī)控制系統(tǒng)進(jìn)行結(jié)合,提高了灌區(qū)灌溉管理的自動(dòng)化水平;駱莉等[9]以輪灌組數(shù)最少為目標(biāo),建立了配水渠道輪灌組合優(yōu)化模型,并用遺傳算法計(jì)算模型,得到了最優(yōu)灌溉組合;劉照等[10]基于RS、GIS技術(shù),選擇總配水時(shí)間最短和輪灌組間引水持續(xù)時(shí)間差最小為目標(biāo),運(yùn)用粒子群算法和蟻群算法求解模型,解決了渠系配水問(wèn)題;Anwar等[11]在單階段模型中考慮了每個(gè)用水戶的預(yù)定時(shí)間與實(shí)際配水時(shí)間的差值,在二階段模型中使第一階段處于劣勢(shì)的用戶優(yōu)先配水,建立了輪灌配水模型,可用于有時(shí)間需求的灌溉計(jì)劃。以上研究均假設(shè)下級(jí)渠道的引水流量相同,而灌區(qū)實(shí)際配水是一個(gè)復(fù)雜的變化過(guò)程,下級(jí)渠道的配水流量通常不可能完全相同。周美林等[12]以引水流量和引水時(shí)間差最小為目標(biāo),建立了下級(jí)渠道流量不等條件下的配水模型,確定了灌區(qū)最優(yōu)輪灌組合;Tonny等[13]通過(guò)簡(jiǎn)單多機(jī)調(diào)度和復(fù)雜多機(jī)調(diào)度,分別確定了在各用水戶流量相同和不相同2種情況下的配水時(shí)間表??偨Y(jié)以往研究結(jié)果可發(fā)現(xiàn),同一輪灌組內(nèi)渠道的配水時(shí)間通常存在差異,給灌區(qū)管理帶來(lái)不便。呂宏興等[14]按照“定流量、變歷時(shí)”的輪灌方式,建立了0-1線性整數(shù)規(guī)劃模型,并將各輪灌組的引水時(shí)間均一化處理,使各輪灌組在同一時(shí)間結(jié)束灌溉;王慶杰等[15]將“組間續(xù)灌、組內(nèi)輪灌”的渠道配水方式轉(zhuǎn)化為“組間輪灌、組內(nèi)續(xù)灌”,在尋找最優(yōu)配水方案的同時(shí),對(duì)各輪灌組內(nèi)渠道的配水時(shí)間進(jìn)行均一化處理,使各輪灌組內(nèi)渠道可以同時(shí)開(kāi)啟和關(guān)閉。

    【切入點(diǎn)】已有的輪灌配水模型為灌區(qū)配水工作提供了一定的技術(shù)參考,但目前的輪灌模型通常將上下游距離較遠(yuǎn)的渠道分到同一組內(nèi),在配水時(shí),灌區(qū)管理人員需要頻繁往返于上下游之間啟閉閘門(mén),給管理工作帶來(lái)較大不便?!緮M解決的關(guān)鍵問(wèn)題】鑒于此,本研究采取“組間輪灌、組內(nèi)續(xù)灌”的灌溉方式,考慮工作人員的工作量,建立輪灌配水模型,并采用NSGA-II對(duì)模型進(jìn)行求解,使模型更加符合實(shí)際的配水情況,為灌區(qū)管理工作提供便利。

    1 模型構(gòu)建與求解

    1.1 灌區(qū)輪灌配水模型的構(gòu)建

    假設(shè)某上級(jí)渠道有條下級(jí)渠道,輪灌組數(shù)的劃分按照式(1)計(jì)算:

    式中:為輪灌組數(shù);d為上級(jí)渠道的設(shè)計(jì)流量(m3/s);=1, 2, …,為下級(jí)渠道的序號(hào);qd為各下級(jí)渠道的設(shè)計(jì)流量(m3/s);ceil為向上取整;floor為向下取整。

    1.1.1 決策變量

    模型決策變量為輪灌編組和下級(jí)渠道的配水流量。決策變量X={0,1}表示第輪灌組第出水口的開(kāi)關(guān)狀態(tài),X=0表示出水口關(guān)閉(1, 2, … ,),X=1表示出水口開(kāi)啟;決策變量q表示下級(jí)渠道的配水流量(m3/s)。

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

    模型的目標(biāo)為配水時(shí)間短和灌區(qū)管理方便。配水時(shí)間短的函數(shù)建立見(jiàn)式(2);為了使管理人員操作方便,建立同一分組內(nèi)各渠道配水時(shí)間差最小的函數(shù),見(jiàn)式(3);為了縮短管理人員的工作距離,同一分組內(nèi),盡可能地使互相靠近的渠道在一個(gè)分組內(nèi),即同一組內(nèi)最上游渠道和最下游渠道的距離最小,見(jiàn)式(4)。

    式中:t為第條下級(jí)渠道中每條渠道的配水時(shí)間(h);w為每條下級(jí)渠道的需水量(m3);ab為同一輪灌組內(nèi)任意2條下級(jí)渠道的配水時(shí)間(h);Δ為同一輪灌組內(nèi)任意2條下級(jí)渠道的配水時(shí)間差(h);M為第個(gè)輪灌組內(nèi)最上游渠道與最下游渠道之間的距離(km)。

    1.1.3 約束條件

    渠道一次性引水約束:任意1條下級(jí)渠道在輪期內(nèi)只能開(kāi)啟1次,即:

    下級(jí)渠道流量約束:任意1條下級(jí)渠道的配水流量應(yīng)在其設(shè)計(jì)流量的0.6~1.0倍以內(nèi),即:

    水量平衡約束:任意時(shí)刻,同時(shí)配水的下級(jí)渠道流量之和應(yīng)等于上級(jí)渠道流量,即:

    式中:Q為上級(jí)渠道時(shí)刻的流量。

    上級(jí)渠道流量約束:任意時(shí)刻,上級(jí)渠道流量應(yīng)在其設(shè)計(jì)流量的0.6~1.0倍之間,即:

    時(shí)間約束:每個(gè)灌溉組的配水時(shí)間不超過(guò)配水渠道最大允許輸水時(shí)間,即:

    1.2 模型求解算法

    本研究建立的模型為非線性約束條件下含有整數(shù)規(guī)劃的模型,并存在多個(gè)約束條件和目標(biāo)函數(shù),模型求解困難。遺傳算法能夠解決多約束多目標(biāo)條件下的大搜索空間優(yōu)化問(wèn)題,相比其他智能算法更具有優(yōu)越性[16-17]。本文采用精英策略的非支配排序遺傳算法(NSGA-II)求解模型。NSGA-II由Deb等[18]提出,是傳統(tǒng)遺傳算法的改進(jìn),降低了計(jì)算復(fù)雜度,避免了較優(yōu)解的丟失。本研究的編碼方式采用二進(jìn)制編碼,交叉算子為兩點(diǎn)交叉,變異算子為均勻變異算子,編程工具為python,算法流程如圖1所示[19]。NSGA-II算法的步驟為:隨機(jī)產(chǎn)生初始種群0;然后根據(jù)目標(biāo)函數(shù)值進(jìn)行非劣排序,再通過(guò)選擇和遺傳操作(交叉、變異)可得新種群Q,精英保留策略有R=RR,經(jīng)過(guò)對(duì)R進(jìn)行非劣排序得到種群P+1;然后迭代運(yùn)行,當(dāng)達(dá)到最大進(jìn)化代數(shù)時(shí)停止運(yùn)行,最終得到一系列非劣解集。

    2 模型應(yīng)用

    2.1 應(yīng)用區(qū)域概況

    簸箕李灌區(qū)位于山東省濱州市最西部,涉及惠民縣、陽(yáng)信縣、無(wú)棣縣,地理坐標(biāo)為東經(jīng)117°14′37″—117°58′44″、北緯37°07′41″—38°14′57″,控制土地面積為22.4萬(wàn)hm2,設(shè)計(jì)灌溉面積為6萬(wàn)hm2。灌區(qū)多年平均氣溫為12~14 ℃,多年平均降水量為570.6 mm,多年平均蒸發(fā)量為1 195.1 mm。本文選取簸箕李灌區(qū)的石皮東支渠及其斗渠為研究對(duì)象(圖2),石皮東支渠位于惠民縣,灌溉面積為1 853 hm2,其控制范圍大部分在惠民縣,少部分在陽(yáng)信縣。灌區(qū)主要種植作物為冬小麥、玉米和棉花,3種作物的灌溉制度見(jiàn)表1。本文以小麥返青期的優(yōu)化配水為研究目標(biāo)。小麥返青期灌水定額為540 m3/hm2,輪灌期為10 d,渠道相關(guān)參數(shù)見(jiàn)表2。

    2.2 模型參數(shù)

    模型決策變量的數(shù)量為60個(gè)。遺傳算法種群規(guī)模為150,交叉概率為0.8,變異概率為0.05,最大遺傳代數(shù)為500。

    圖2 渠系位置及分布

    表1 不同作物灌溉制度

    2.3 結(jié)果與分析

    表3為模型的輪灌組劃分情況和每條斗渠的配水時(shí)間。模型將30條斗渠劃分為9個(gè)組,9個(gè)輪灌組之間依次進(jìn)行灌溉,每個(gè)分組內(nèi)的斗渠同時(shí)進(jìn)行灌溉。灌溉開(kāi)始時(shí),第一分組內(nèi)的1、2、3、4、6斗同時(shí)開(kāi)啟,當(dāng)?shù)谝环纸M內(nèi)的斗渠全部灌溉結(jié)束時(shí),第二分組的5、10、12斗開(kāi)始灌溉,以此類推,直至第九組灌溉結(jié)束。每個(gè)分組的配水時(shí)長(zhǎng)為該分組內(nèi)斗渠配水時(shí)間的最大值,9個(gè)組的配水時(shí)長(zhǎng)分別為15.64、22.9、15.83、16.99、21.51、38.21、20.37、21.06、28.49 h,總配水時(shí)長(zhǎng)為9個(gè)組的配水時(shí)長(zhǎng)之和,為201.54 h,與配水輪期240 h相比,本文提出的配水模型縮短了總配水時(shí)間。在輪灌組內(nèi),各斗渠之間的配水時(shí)間也較為接近,斗渠之間配水時(shí)間的最大差值為第一輪灌組內(nèi)1斗渠和3斗渠的配水時(shí)間差值,為15 min。

    表2 石皮東支渠及斗渠參數(shù)

    表3 斗渠分組與配水時(shí)間

    圖3為斗渠配水流量與設(shè)計(jì)流量的對(duì)比,所有斗渠配水流量均在設(shè)計(jì)流量的0.60~1倍之間,滿足斗渠過(guò)流能力要求。

    圖3 斗渠配水流量與設(shè)計(jì)流量之間的對(duì)比

    圖4反映了斗渠啟閉順序和相鄰啟閉斗渠的間距。模型考慮了同一輪灌內(nèi)最上游斗渠與最下游斗渠間距最小、同一輪灌組內(nèi)斗渠灌溉時(shí)間差距最小2個(gè)目標(biāo)。模型最后的分組結(jié)果顯示,同一組內(nèi)的斗渠基本上是相鄰斗渠,且同一輪灌組內(nèi)的斗渠灌溉時(shí)間最大差值在15 min內(nèi),可近似認(rèn)為同一輪灌組之間的斗渠是可以統(tǒng)一關(guān)閉的,所以開(kāi)啟和關(guān)閉斗渠的順序是相同的,斗渠啟閉順序?yàn)?→2→3→4→6→5→10→12→7→8→9→11→13→14→15→16→17→18→19→23→24→20→21→22→25→26→30→27→28→29。在斗渠啟閉過(guò)程中,盡可能按照渠道上游到下游的順序來(lái)啟閉,以此來(lái)減少工作人員的工作距離,管理員的工作總距離為相鄰啟閉的斗渠間距之和,即30.35 km;管理員逆向工作的距離為從下游到上游的總距離,即7.21 km,模型實(shí)現(xiàn)了最大程度降低管理員的工作距離。

    圖4 模型斗渠啟閉順序和相鄰啟閉斗渠的間距(km)

    3 討論

    本文在建立模型時(shí),流量下限設(shè)置為設(shè)計(jì)流量的0.6倍,模型最后的配水結(jié)果也顯示所有斗渠的流量均滿足此要求,實(shí)現(xiàn)了大流量配水。圖5為石皮東支渠配水流量與設(shè)計(jì)流量的對(duì)比,支渠配水流量在1.88~2 m3/s之間,是設(shè)計(jì)流量的0.94~1倍,支渠同樣也實(shí)現(xiàn)了大流量運(yùn)行,充分發(fā)揮了渠道工程能力。

    圖5 支渠配水流量與設(shè)計(jì)流量之間的對(duì)比

    以往配水模型的配水方式大多是“組間輪灌,組內(nèi)續(xù)灌”,按照這種方式進(jìn)行灌溉時(shí),灌區(qū)管理人員需要頻繁往返于上、下游渠道之間啟閉閘門(mén),給管理帶來(lái)巨大不便。本文將配水方式調(diào)整為“組間續(xù)灌,組內(nèi)輪灌”,在這種灌溉方式下,建立了2個(gè)目標(biāo)函數(shù)。第一個(gè)目標(biāo)是同一組內(nèi)渠道配水時(shí)間差最小,實(shí)現(xiàn)同一組內(nèi)的渠道同時(shí)開(kāi)啟或關(guān)閉,方便灌區(qū)管理人員管理,這在具有遠(yuǎn)程控制系統(tǒng)的灌區(qū)更為實(shí)用;第二個(gè)目標(biāo)是同一分組內(nèi)的渠道應(yīng)盡可能地相互靠近,該目標(biāo)是為了使灌區(qū)管理人員在啟閉同一組之間的閘門(mén)時(shí),能夠減少其工作距離,不需要頻繁往返于上下游渠道之間。結(jié)果表明,每個(gè)分組內(nèi)的渠道基本相互臨近,上下游相距較遠(yuǎn)的渠道不會(huì)被分到一個(gè)組內(nèi);每個(gè)分組內(nèi),斗渠配水時(shí)間差值最大為15 min,可近似認(rèn)為灌溉在同一時(shí)間結(jié)束,可同時(shí)關(guān)閉閘門(mén)。本文沒(méi)有考慮多級(jí)渠道,存在一定局限性,未來(lái)仍需進(jìn)一步完善。

    4 結(jié)論

    以各斗渠流量和輪灌分組情況為決策變量,以配水時(shí)間最短、同一輪灌內(nèi)最上游渠道與最下游渠道間距最小、同一分組內(nèi)的渠道灌溉時(shí)間差最小為目標(biāo),建立配水模型。所建立的配水模型在斗渠配水流量和支渠配水流量方面均滿足過(guò)流能力要求,并實(shí)現(xiàn)了大流量配水;模型將渠道分為9組,總配水時(shí)長(zhǎng)為201.54 h,滿足輪期要求,縮短了灌溉時(shí)間;根據(jù)本文模型進(jìn)行渠道配水時(shí),灌區(qū)管理人員不需頻繁往返于上下游之間,同一分組內(nèi)的渠道可同時(shí)啟閉,大大降低了管理人員的工作量。

    [1] 張智韜, 李援農(nóng), 陳俊英, 等. 基于3S技術(shù)和蟻群算法的灌區(qū)渠系優(yōu)化配水[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào), 2010, 38(7): 221-226.

    ZHANG Zhitao, LI Yuannong, CHEN Junying, et al. Optimization of water allocation management for irrigation canals based on 3S technology and colony algorithm[J]. Journal of Northwest A&F University, 2010, 38(7): 221-226.

    [2] 李彤姝, 黃睿, 孫志鵬, 等. 基于多目標(biāo)粒子群算法的渠系優(yōu)化配水研究[J]. 灌溉排水學(xué)報(bào), 2020, 39(9): 95-100, 125.

    LI Tongshu, HUANG Rui, SUN Zhipeng, et al. Optimizing water distribution in canal networks using multi-objective particle swarm optimization method[J]. Journal of Irrigation and Drainage, 2020, 39(9): 95-100, 125.

    [3] 徐淑琴, 高凱茹, 樂(lè)靜, 等. 基于NSGA-II灌區(qū)兩級(jí)渠道輸配水優(yōu)化調(diào)度[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2020, 51(3): 71-78.

    XU Shuqin, GAO Kairu, YUE Jing, et al. Based on the NSGA-Ⅱ irrigation channels of two-stage water delivery and distribution optimization scheduling[J]. Journal of Northeast Agricultural University, 2020, 51(3): 71-78.

    [4] 沈來(lái)銀, 胡鐵松, 周姍, 等. 基于SHAW模型的河套灌區(qū)秋澆渠系優(yōu)化配水模型研究[J]. 水利學(xué)報(bào), 2020, 51(4): 458-467.

    SHEN Laiyin, HU Tiesong, ZHOU Shan, et al. Study on optimal water distribution model of canal system for autumn irrigation in Hetao Irrigation Area based on SHAW model[J]. Journal of Hydraulic Engineering, 2020, 51(4): 458-467.

    [5] 賀建文, 何英. 基于NSGA-II遺傳算法的豐收灌區(qū)優(yōu)化配水研究[J]. 人民黃河, 2020, 42(2): 276-278, 286.

    [6] SURYAVANSHI A R, REDDY J M. Optimal operation schedule of irrigation distribution systems[J]. Agricultural Water Management, 1986, 11: 23-30.

    [7] 王智, 朱鳳書(shū). 灌溉渠系流量最優(yōu)調(diào)配的0-1規(guī)劃模型[J]. 灌溉排水, 1992, 11(3): 8-13.

    WANG Zhi, ZHU Fengshu. 0-1 Programming model for optimal flow regulation in irrigation canal systems[J]. Irrigation and Drainage, 1992, 11(3): 8-13.

    [8] 李彬, 曾慶祝. 灌區(qū)計(jì)算機(jī)監(jiān)控系統(tǒng)與灌溉優(yōu)化調(diào)度模型集成應(yīng)用研究[J]. 排灌機(jī)械, 2005, 23(4): 25-28.

    LI Bin, ZENG Qingzhu. Study on integration between computer monitoring system and optimization model in irrigation channel distribution[J]. Journal of Drainage and Irrigation Machinery Engineering, 2005, 23(4): 25-28.

    [9] 駱莉, 袁宏源. 配水渠道輪灌組合優(yōu)化模型及其遺傳算法[J]. 中國(guó)農(nóng)村水利水電, 2000, 10: 11-13.

    [10] 劉照, 華慶偉, 張成才, 等. 基于RS、GIS及智能算法的渠系優(yōu)化配水[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào), 2017, 45(4): 213-222, 229.

    LIU Zhao, HUA Qingwei, ZHANG Chengcai, et al. Optimal irrigation water distribution based on RS, GIS and intelligent algorithms[J]. Journal of Northwest A & F University, 2017, 45(4): 213-222, 229.

    [11] ANWAR A A, CLARKE, D. Irrigation scheduling using mixed-integer linear programming[J]. Journal of Irrigation and Drainage Engineering, 2001, 127: 63-69.

    [12] 周美林, 呂宏興, 韓文霆. 渠系配水優(yōu)化模型和多目標(biāo)遺傳算法研究[J]. 中國(guó)農(nóng)村水利水電, 2014(9): 5-7.

    ZHOU Meilin, LYU Hongxing, HAN Wenting. Optimal water scheduling model and multi-objective genetic algorithm[J]. China Rural Water and Hydropowe, 2014(9): 5-7.

    [13] TONNY TV, ANWAR A. Irrigation scheduling using complex machine scheduling[J]. Journal of Irrigation and Drainage Engineering, 2014, 4(14 065): 1-8.

    [14] 呂宏興, 熊運(yùn)章, 汪志農(nóng). 灌溉渠道支斗渠輪灌配水與引水時(shí)間優(yōu)化模型[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2000, 16(6): 43-46.

    LYU Hongxing, XIONG Yunzhang, WANG Zhinong. Optimal model of rotation irrigation distribution channel and branch canal and delivery time[J]. Transactions of the Chinese Society of Agricultural Engineering, 2000, 16(6): 43-46.

    [15] 王慶杰, 岳春芳, 李藝珍, 等. 基于改進(jìn)粒子群算法的兩級(jí)渠道水資源優(yōu)化配置[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2019, 37(4): 26-33.

    WANG Qingjie, YUE Chunfang, LI Yizhen, et al. Optimal allocation of water resources with two-level channel based on improved particle swarm optimization algorithm[J]. Agricultural Research in the Arid Areas, 2019, 37(4): 26-33.

    [16] 趙恩龍. 基于多目標(biāo)遺傳算法的灌區(qū)水資源優(yōu)化調(diào)度研究[D]. 武漢:長(zhǎng)江科學(xué)院, 2013.

    ZHAO Enlong. Study of optimal water resources allocation of irrigation district based on multi-objective genetic algorithm evolution[D]. Wuhan: Changjiang River Scientific Research Institute, 2013.

    [17] 趙文舉, 馬孝義, 劉哲, 等. 多級(jí)灌溉渠系配水優(yōu)化編組模型與算法研究[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(2): 11-16.

    ZHAO Wenju, MA Xiaoyi, LIU Zhe, et al. Optimal water allocation marshalling model of multilevel irrigation canal system and model Solution[J]. Transactions of the Chinese Society of Agricultural Engineering, 2008, 24(2): 11-16.

    [18] KALYANMOY Deb, AMRIT Pratap, SAMEER Agarwal, et al. A fast and elitist multi-objective genetic algorithm[J]. IEEE Transactions on Evolutionary Computation, 2002, 6(2): 182-197.

    [19] 王瑞峰, 陶榮杰. 基于NSGA2算法的ZPW2000A軌道電路維修策略優(yōu)化研究王瑞峰[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2018, 15(9): 2 394-2 400.

    WANG Ruifeng, TAO Rongjie. Research on optimization of ZPW2000A track circuit maintenance strategy based on NSGA2 algorithm[J]. Journal of Railway Science and Engineering, 2018, 15(9): 2 394-2 400.

    Optimizing Water Distribution in Irrigation Canals Using the NSGA-II Algorithm

    FAN Yu1,2, GAO Zhanyi1,2, CHEN Haorui1,2*, WANG Yunhui3, FANG Benyan3, XU Ning3

    (1. State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China;2. National Center of Efficient Irrigation Engineering and Technology Research, Beijing 100048, China; 3. Yellow River Irrigation Service Center, Binzhou Municipality, Binzhou 256600, China)

    【Objective】Optimization of water distribution in irrigation districts is a way to improve water use efficiency and water resource management. This paper proposes a method to optimize water distribution in irrigation canals. 【Method】The flow rate in the canals at the lower end and the number of rotational groups are taken as the decision variables. The objectives of the optimization are to minimize the total water distribution duration, the distance between the upstream and downstream canals in the same group, and the difference in water distribution duration in each canal in the same rotational group. The optimization considers the constraints of time, outlet state, canal overflow capacity, and water balance. The model is applied to the Shipi East branch canal in Bojili irrigation regions, and the optimization is solved using the NSGA-II algorithm by dividing the lower-end canals into nine rotation irrigation groups. 【Result】Dividing the adjacent canals into one round grouting group as possible as it could, and opening and closing them simultaneously to facilitate the management, the water distribution flow rate and total water distribution time calculated by the model meet the operational requirements.【Conclusion】The optimization model proposed in this paper can distribute water between canals to achieve high flow rate, short duration and less workload. It can be used to help improve water management in irrigation districts.

    multi-objective water distribution; canals water distribution; rotational irrigation; optimization model

    樊煜, 高占義, 陳皓銳, 等. 基于NSGA-II算法的灌區(qū)配水渠道輪灌分組優(yōu)化調(diào)控研究[J]. 灌溉排水學(xué)報(bào), 2023, 42(2): 130-135.

    FAN Yu, GAO Zhanyi, CHEN Haorui, et al.Optimizing Water Distribution in Irrigation Canals Using the NSGA-II Algorithm[J]. Journal of Irrigation and Drainage, 2023, 42(2): 130-135.

    1672 - 3317(2023)02 - 0130 - 06

    S274.3

    A

    10.13522/j.cnki.ggps.2022346

    2022-06-22

    “科技興蒙”行動(dòng)重點(diǎn)專項(xiàng)(2021EEDSCXSFQZD010);中國(guó)水利水電科學(xué)研究院技術(shù)創(chuàng)新團(tuán)隊(duì)項(xiàng)目(ID0145B022021)

    樊煜(1993-),男。博士研究生,主要從事灌區(qū)渠系配水及調(diào)控研究。E-mail: 673000238@qq.com

    陳皓銳(1984-),男。教授級(jí)高級(jí)工程師,主要從事節(jié)水灌溉理論與技術(shù)要求。E-mail:chenhr@iwhr.com

    責(zé)任編輯:韓 洋

    猜你喜歡
    斗渠分組灌溉
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    基于灌溉渠水質(zhì)保護(hù)的都江堰灌區(qū)斗渠廊道LID生態(tài)措施研究
    分組搭配
    怎么分組
    分組
    某縣新增千億斤糧食產(chǎn)能規(guī)劃渠道工程典型設(shè)計(jì)
    治淮(2016年7期)2016-09-02 02:55:41
    安陽(yáng)市灌區(qū)支斗渠管道化改造可行性分析
    极品人妻少妇av视频| 国产精品综合久久久久久久免费 | 中文字幕制服av| 久久精品国产a三级三级三级| 日韩欧美一区二区三区在线观看 | 成人精品一区二区免费| 老汉色∧v一级毛片| 久久国产精品男人的天堂亚洲| 老汉色av国产亚洲站长工具| 最新的欧美精品一区二区| 久久久久国产一级毛片高清牌| 脱女人内裤的视频| 国产高清视频在线播放一区| 久久香蕉国产精品| 久久久久国内视频| 精品第一国产精品| 女人被躁到高潮嗷嗷叫费观| 亚洲中文av在线| 久久午夜亚洲精品久久| 怎么达到女性高潮| 久久久久久久国产电影| 免费高清在线观看日韩| 国产成人av教育| 久久青草综合色| 精品国产亚洲在线| 日韩三级视频一区二区三区| 国产欧美日韩一区二区精品| 真人做人爱边吃奶动态| 免费看a级黄色片| 久久久久久久午夜电影 | 18禁黄网站禁片午夜丰满| 91老司机精品| 老汉色∧v一级毛片| 欧美黄色片欧美黄色片| 又大又爽又粗| av天堂久久9| av免费在线观看网站| 日韩中文字幕欧美一区二区| 国产成人av激情在线播放| 亚洲人成伊人成综合网2020| 国内毛片毛片毛片毛片毛片| 叶爱在线成人免费视频播放| 欧美在线一区亚洲| 高清毛片免费观看视频网站 | 国产一区二区三区在线臀色熟女 | 最近最新免费中文字幕在线| 亚洲人成电影免费在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩一区二区三区精品不卡| 黑人操中国人逼视频| 黄色a级毛片大全视频| x7x7x7水蜜桃| 岛国在线观看网站| 亚洲成人国产一区在线观看| 免费观看a级毛片全部| 淫妇啪啪啪对白视频| 天堂俺去俺来也www色官网| 两个人免费观看高清视频| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 欧美亚洲日本最大视频资源| 岛国毛片在线播放| 国产97色在线日韩免费| 无限看片的www在线观看| 新久久久久国产一级毛片| 女性生殖器流出的白浆| 三级毛片av免费| 一级毛片高清免费大全| 一级a爱片免费观看的视频| 亚洲欧美一区二区三区黑人| 多毛熟女@视频| 精品福利观看| 一区福利在线观看| 母亲3免费完整高清在线观看| 老司机福利观看| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕| 黄片大片在线免费观看| 一级毛片女人18水好多| 老汉色∧v一级毛片| 十八禁网站免费在线| 波多野结衣一区麻豆| 人人澡人人妻人| 国产成人欧美| 久久国产精品大桥未久av| 久久久精品国产亚洲av高清涩受| 午夜福利影视在线免费观看| 日韩精品免费视频一区二区三区| 欧美亚洲 丝袜 人妻 在线| 国产日韩一区二区三区精品不卡| 91字幕亚洲| 午夜福利在线免费观看网站| 成人三级做爰电影| 亚洲精品在线美女| 国产一区二区三区在线臀色熟女 | 午夜福利影视在线免费观看| 国产蜜桃级精品一区二区三区 | 韩国精品一区二区三区| 精品福利永久在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品在线观看二区| 久久久久久久久久久久大奶| 黄色怎么调成土黄色| 女人被躁到高潮嗷嗷叫费观| 狠狠婷婷综合久久久久久88av| 国产97色在线日韩免费| 欧美不卡视频在线免费观看 | 久久精品成人免费网站| 黄色毛片三级朝国网站| 午夜福利视频在线观看免费| 久久精品成人免费网站| 成人精品一区二区免费| avwww免费| 中国美女看黄片| 亚洲在线自拍视频| 国产亚洲精品第一综合不卡| 高清黄色对白视频在线免费看| 免费久久久久久久精品成人欧美视频| 久久婷婷成人综合色麻豆| 19禁男女啪啪无遮挡网站| 国产麻豆69| 欧美 亚洲 国产 日韩一| 午夜福利乱码中文字幕| 大香蕉久久成人网| 在线观看舔阴道视频| 丁香欧美五月| 成人黄色视频免费在线看| 免费一级毛片在线播放高清视频 | 好男人电影高清在线观看| 日韩视频一区二区在线观看| 国产精品永久免费网站| 大型黄色视频在线免费观看| av网站免费在线观看视频| 久久影院123| 久久久久久人人人人人| 欧美激情高清一区二区三区| 最新美女视频免费是黄的| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 久久久久久亚洲精品国产蜜桃av| 女同久久另类99精品国产91| 午夜福利在线免费观看网站| 亚洲国产中文字幕在线视频| 老司机靠b影院| 色94色欧美一区二区| 久久精品国产99精品国产亚洲性色 | 少妇粗大呻吟视频| 亚洲av日韩在线播放| 丰满人妻熟妇乱又伦精品不卡| 天堂俺去俺来也www色官网| 成人免费观看视频高清| 熟女少妇亚洲综合色aaa.| 午夜两性在线视频| 日韩欧美免费精品| 亚洲在线自拍视频| 国产一区在线观看成人免费| 亚洲熟女精品中文字幕| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 精品无人区乱码1区二区| 欧美精品啪啪一区二区三区| 欧美日本中文国产一区发布| 人妻久久中文字幕网| 日本黄色视频三级网站网址 | 久久久国产成人精品二区 | 两性夫妻黄色片| 欧美日韩福利视频一区二区| 人妻久久中文字幕网| 韩国av一区二区三区四区| 999精品在线视频| 午夜福利一区二区在线看| 性色av乱码一区二区三区2| 免费在线观看黄色视频的| 国产高清videossex| 久久精品亚洲熟妇少妇任你| 三上悠亚av全集在线观看| 欧美老熟妇乱子伦牲交| 人成视频在线观看免费观看| 老司机午夜十八禁免费视频| 国产99久久九九免费精品| 精品国产亚洲在线| 午夜两性在线视频| 亚洲成人国产一区在线观看| 久久久国产成人精品二区 | 国产欧美日韩一区二区三| 国产欧美日韩综合在线一区二区| 老司机在亚洲福利影院| 国产精品一区二区在线观看99| 色综合欧美亚洲国产小说| 超碰97精品在线观看| 欧美日韩乱码在线| 久久青草综合色| 天天影视国产精品| 免费在线观看黄色视频的| 亚洲国产精品合色在线| 亚洲av电影在线进入| 国产精品久久电影中文字幕 | 欧美亚洲 丝袜 人妻 在线| 王馨瑶露胸无遮挡在线观看| 欧美黑人精品巨大| 操美女的视频在线观看| 婷婷成人精品国产| 亚洲精品中文字幕一二三四区| 亚洲在线自拍视频| av在线播放免费不卡| 99精品久久久久人妻精品| 色婷婷久久久亚洲欧美| 国产一区二区三区在线臀色熟女 | 国产亚洲精品久久久久久毛片 | 久热爱精品视频在线9| 国产一区二区激情短视频| 色综合欧美亚洲国产小说| 欧美日韩亚洲高清精品| 又大又爽又粗| a级毛片黄视频| 国产一区二区三区综合在线观看| av网站在线播放免费| ponron亚洲| 欧美在线一区亚洲| 欧美乱色亚洲激情| 亚洲av日韩精品久久久久久密| 免费久久久久久久精品成人欧美视频| 国产日韩欧美亚洲二区| 久久香蕉激情| 亚洲一区中文字幕在线| 国产精品九九99| 欧洲精品卡2卡3卡4卡5卡区| 日本五十路高清| 99riav亚洲国产免费| 黑人操中国人逼视频| 在线免费观看的www视频| 99热网站在线观看| 老汉色∧v一级毛片| 大型av网站在线播放| 一级毛片精品| 天堂动漫精品| www.熟女人妻精品国产| 一级a爱片免费观看的视频| 免费观看a级毛片全部| 亚洲色图av天堂| 国产亚洲欧美在线一区二区| 黑人猛操日本美女一级片| 国产男靠女视频免费网站| 超色免费av| 日韩免费高清中文字幕av| 人人妻人人爽人人添夜夜欢视频| 可以免费在线观看a视频的电影网站| 国产成人欧美| 久久精品亚洲av国产电影网| 变态另类成人亚洲欧美熟女 | 捣出白浆h1v1| 精品第一国产精品| 国产精品久久久人人做人人爽| 操美女的视频在线观看| 成人国产一区最新在线观看| 久久午夜亚洲精品久久| 亚洲熟女毛片儿| а√天堂www在线а√下载 | 一级,二级,三级黄色视频| 亚洲av片天天在线观看| 亚洲精品粉嫩美女一区| 精品国产乱子伦一区二区三区| 一区福利在线观看| 老司机福利观看| 亚洲av成人不卡在线观看播放网| 久久久国产成人免费| 一二三四在线观看免费中文在| 国产精品自产拍在线观看55亚洲 | 午夜亚洲福利在线播放| 色综合婷婷激情| 日韩欧美三级三区| 欧美另类亚洲清纯唯美| 亚洲国产毛片av蜜桃av| 操美女的视频在线观看| 日韩欧美国产一区二区入口| 69精品国产乱码久久久| 手机成人av网站| 亚洲九九香蕉| 99久久99久久久精品蜜桃| 免费在线观看完整版高清| 欧美一级毛片孕妇| 男女午夜视频在线观看| 国产欧美日韩一区二区精品| 999久久久精品免费观看国产| 99国产精品一区二区蜜桃av | 国产精品 国内视频| 久久精品亚洲精品国产色婷小说| 99精国产麻豆久久婷婷| 少妇猛男粗大的猛烈进出视频| 亚洲av美国av| av天堂久久9| 国产精品香港三级国产av潘金莲| 美女 人体艺术 gogo| 精品国产亚洲在线| 久久性视频一级片| 变态另类成人亚洲欧美熟女 | 久久中文字幕人妻熟女| 亚洲av成人一区二区三| 欧美日韩中文字幕国产精品一区二区三区 | 捣出白浆h1v1| xxxhd国产人妻xxx| 精品乱码久久久久久99久播| 国产片内射在线| 国产精品久久久av美女十八| 久久久国产成人精品二区 | 99国产综合亚洲精品| 黄色毛片三级朝国网站| 国产三级黄色录像| 欧美最黄视频在线播放免费 | 午夜福利在线免费观看网站| 色综合婷婷激情| 日韩欧美三级三区| 日韩有码中文字幕| 91精品国产国语对白视频| 亚洲欧美日韩另类电影网站| 亚洲第一av免费看| 91成年电影在线观看| 欧美成狂野欧美在线观看| 久久久久久久久免费视频了| 婷婷成人精品国产| 俄罗斯特黄特色一大片| 欧美一级毛片孕妇| 亚洲专区国产一区二区| 99在线人妻在线中文字幕 | 女人被躁到高潮嗷嗷叫费观| 女人久久www免费人成看片| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜精品一区,二区,三区| 日韩欧美三级三区| 精品一区二区三区四区五区乱码| 无限看片的www在线观看| 午夜免费观看网址| 美女 人体艺术 gogo| 最新美女视频免费是黄的| 亚洲国产毛片av蜜桃av| 中出人妻视频一区二区| 精品一区二区三卡| 成人国语在线视频| 日韩成人在线观看一区二区三区| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| 成人国语在线视频| 免费少妇av软件| 国产成人欧美在线观看 | 黑人猛操日本美女一级片| 日韩精品免费视频一区二区三区| 欧美中文综合在线视频| 国产av一区二区精品久久| 国产男女内射视频| 两个人免费观看高清视频| 欧美中文综合在线视频| 国产精品国产av在线观看| 国产精品欧美亚洲77777| 成年动漫av网址| 国产亚洲欧美精品永久| 久久久精品免费免费高清| 免费久久久久久久精品成人欧美视频| 80岁老熟妇乱子伦牲交| 久久影院123| 91麻豆精品激情在线观看国产 | 一边摸一边做爽爽视频免费| 国产精品偷伦视频观看了| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 亚洲成人免费电影在线观看| tocl精华| 久久久久久久国产电影| 性少妇av在线| 亚洲国产毛片av蜜桃av| 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| svipshipincom国产片| av欧美777| 久久久国产一区二区| 国产真人三级小视频在线观看| 精品午夜福利视频在线观看一区| 久热这里只有精品99| 国产精品欧美亚洲77777| 日韩欧美一区视频在线观看| 国产aⅴ精品一区二区三区波| 国产91精品成人一区二区三区| 伦理电影免费视频| 亚洲成av片中文字幕在线观看| xxx96com| 岛国毛片在线播放| 91九色精品人成在线观看| 日韩视频一区二区在线观看| 不卡av一区二区三区| 久久精品国产亚洲av高清一级| 757午夜福利合集在线观看| 亚洲av熟女| 国产精品国产高清国产av | 三上悠亚av全集在线观看| 人人妻人人澡人人看| 亚洲 国产 在线| x7x7x7水蜜桃| 亚洲第一av免费看| 高清在线国产一区| 一二三四社区在线视频社区8| 精品国内亚洲2022精品成人 | 最新美女视频免费是黄的| 欧美黄色片欧美黄色片| 国产99久久九九免费精品| 久久久久久久午夜电影 | 久久久国产欧美日韩av| 色婷婷av一区二区三区视频| 99re在线观看精品视频| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 中文字幕最新亚洲高清| 日本黄色日本黄色录像| 久久精品国产综合久久久| 午夜亚洲福利在线播放| 九色亚洲精品在线播放| 欧美乱妇无乱码| 悠悠久久av| 国产成人一区二区三区免费视频网站| 波多野结衣av一区二区av| 老司机福利观看| av视频免费观看在线观看| 国产精品一区二区在线不卡| 99精品欧美一区二区三区四区| av国产精品久久久久影院| 在线天堂中文资源库| 成人精品一区二区免费| 18禁美女被吸乳视频| 水蜜桃什么品种好| 国产精华一区二区三区| 久久婷婷成人综合色麻豆| 三级毛片av免费| 波多野结衣av一区二区av| 国产亚洲欧美精品永久| 99热只有精品国产| 看黄色毛片网站| 欧美激情高清一区二区三区| 精品亚洲成a人片在线观看| 丝袜人妻中文字幕| 女人被躁到高潮嗷嗷叫费观| 波多野结衣一区麻豆| 亚洲精品成人av观看孕妇| 国产成人影院久久av| 国产99白浆流出| 在线观看免费午夜福利视频| 日本五十路高清| 自拍欧美九色日韩亚洲蝌蚪91| 丰满的人妻完整版| 免费一级毛片在线播放高清视频 | 国产免费av片在线观看野外av| 成人影院久久| 日本精品一区二区三区蜜桃| 一级毛片高清免费大全| 亚洲国产精品一区二区三区在线| 亚洲欧美日韩另类电影网站| 免费在线观看黄色视频的| 村上凉子中文字幕在线| 高清在线国产一区| 国产亚洲欧美在线一区二区| 老司机亚洲免费影院| 麻豆av在线久日| 久久久久久免费高清国产稀缺| 欧美日韩一级在线毛片| 新久久久久国产一级毛片| 777米奇影视久久| 国产av一区二区精品久久| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区在线观看99| 欧美精品啪啪一区二区三区| 久久精品亚洲av国产电影网| av天堂在线播放| 欧美成人免费av一区二区三区 | 婷婷丁香在线五月| 18禁国产床啪视频网站| 亚洲成人手机| 欧美精品人与动牲交sv欧美| 悠悠久久av| 欧美黄色淫秽网站| 一区二区日韩欧美中文字幕| 99精品在免费线老司机午夜| 18禁美女被吸乳视频| 在线观看免费视频日本深夜| bbb黄色大片| 国产一区在线观看成人免费| 亚洲在线自拍视频| 精品国产乱码久久久久久男人| 无限看片的www在线观看| 精品国产一区二区三区久久久樱花| 精品国产乱子伦一区二区三区| 国产不卡一卡二| 操美女的视频在线观看| 亚洲精品在线观看二区| av片东京热男人的天堂| 国产精品久久久人人做人人爽| 日日夜夜操网爽| 18禁美女被吸乳视频| 亚洲成人免费电影在线观看| 日本撒尿小便嘘嘘汇集6| 一区福利在线观看| 亚洲精品av麻豆狂野| 国产av精品麻豆| 亚洲欧美精品综合一区二区三区| 99精品久久久久人妻精品| 国产在线观看jvid| 精品久久久久久久久久免费视频 | 亚洲精品自拍成人| 久久久久久久国产电影| 亚洲avbb在线观看| 国产一区二区激情短视频| 欧美日韩瑟瑟在线播放| 国产熟女午夜一区二区三区| 久久国产精品男人的天堂亚洲| 久久人人爽av亚洲精品天堂| 丝袜美腿诱惑在线| 午夜两性在线视频| 国产97色在线日韩免费| 精品乱码久久久久久99久播| 午夜福利影视在线免费观看| 桃红色精品国产亚洲av| 国产野战对白在线观看| 久久ye,这里只有精品| 国产精品一区二区在线观看99| tube8黄色片| 伦理电影免费视频| 老司机影院毛片| 侵犯人妻中文字幕一二三四区| 免费观看a级毛片全部| 搡老岳熟女国产| 丰满迷人的少妇在线观看| 极品教师在线免费播放| 亚洲三区欧美一区| 国产精品亚洲av一区麻豆| 91麻豆精品激情在线观看国产 | 亚洲美女黄片视频| 9色porny在线观看| 久久人妻熟女aⅴ| 成人特级黄色片久久久久久久| 亚洲人成伊人成综合网2020| 久久精品国产清高在天天线| 女人爽到高潮嗷嗷叫在线视频| 亚洲一区二区三区欧美精品| 亚洲av成人不卡在线观看播放网| 国产片内射在线| 亚洲性夜色夜夜综合| 亚洲欧美激情在线| 欧美黄色淫秽网站| 91精品三级在线观看| 欧美精品啪啪一区二区三区| 在线观看日韩欧美| 精品一区二区三卡| 真人做人爱边吃奶动态| 视频在线观看一区二区三区| 免费看a级黄色片| 9色porny在线观看| 亚洲七黄色美女视频| 女人被躁到高潮嗷嗷叫费观| 法律面前人人平等表现在哪些方面| 国产精品.久久久| 在线观看66精品国产| 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 亚洲第一av免费看| 国产激情欧美一区二区| 国产精品久久久久久人妻精品电影| 999久久久国产精品视频| svipshipincom国产片| 亚洲人成77777在线视频| 久久香蕉精品热| 国产单亲对白刺激| 亚洲精品美女久久久久99蜜臀| 亚洲一区中文字幕在线| 无人区码免费观看不卡| 另类亚洲欧美激情| 香蕉久久夜色| 欧美色视频一区免费| av欧美777| 亚洲国产欧美网| 18在线观看网站| videos熟女内射| 精品少妇久久久久久888优播| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区蜜桃| 久久午夜综合久久蜜桃| 精品免费久久久久久久清纯 | 成熟少妇高潮喷水视频| 女人久久www免费人成看片| 又紧又爽又黄一区二区| 亚洲精品乱久久久久久| 国产精品九九99| 国产精品亚洲一级av第二区| 狂野欧美激情性xxxx| 亚洲精品美女久久av网站| 老熟妇乱子伦视频在线观看| 极品人妻少妇av视频| 天堂动漫精品| 欧美精品亚洲一区二区| 亚洲免费av在线视频| 在线天堂中文资源库| 一级黄色大片毛片| 可以免费在线观看a视频的电影网站| 19禁男女啪啪无遮挡网站| 一级黄色大片毛片| 少妇的丰满在线观看| 捣出白浆h1v1| 一级黄色大片毛片| 久久久久视频综合| 五月开心婷婷网| 日本黄色视频三级网站网址 | 国产一区二区三区综合在线观看| 国产精品1区2区在线观看. | 老熟妇仑乱视频hdxx| 丰满人妻熟妇乱又伦精品不卡| 免费观看人在逋| 男人舔女人的私密视频| 丰满迷人的少妇在线观看| 交换朋友夫妻互换小说| 国产野战对白在线观看|