• <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久久男人的天堂| 中文欧美无线码| 老司机影院成人| freevideosex欧美| 国产视频内射| 日本黄色日本黄色录像| 欧美性感艳星| 亚洲色图av天堂| 欧美zozozo另类| 欧美精品人与动牲交sv欧美| 男人和女人高潮做爰伦理| 在线看a的网站| 夫妻午夜视频| 日韩视频在线欧美| 欧美一区二区亚洲| 免费黄网站久久成人精品| 成人二区视频| 99精国产麻豆久久婷婷| 午夜福利视频精品| 欧美区成人在线视频| 人妻系列 视频| 日韩一本色道免费dvd| 久久精品久久久久久久性| 久久影院123| 亚洲精品亚洲一区二区| av专区在线播放| 国产精品久久久久久精品古装| 日韩不卡一区二区三区视频在线| 成人午夜精彩视频在线观看| 美女福利国产在线 | 久久久久网色| 最黄视频免费看| 国产精品人妻久久久久久| 亚洲精品成人av观看孕妇| 人人妻人人添人人爽欧美一区卜 | 国产av精品麻豆| 18禁动态无遮挡网站| 免费大片18禁| 在线观看免费日韩欧美大片 | 国产成人精品久久久久久| 国产高清三级在线| 亚洲四区av| 精品久久久久久久久av| 国产精品无大码| 九草在线视频观看| 色视频www国产| 亚洲精品久久午夜乱码| 亚洲精品视频女| 少妇的逼水好多| 久久人人爽人人片av| 少妇的逼好多水| 男女边摸边吃奶| 成人一区二区视频在线观看| 一区二区三区四区激情视频| 久久久久久久久久久免费av| av女优亚洲男人天堂| 91久久精品电影网| 青青草视频在线视频观看| 成人毛片60女人毛片免费| 免费黄网站久久成人精品| 成人三级做爰电影| 成人午夜精彩视频在线观看| 我的亚洲天堂| 天堂中文最新版在线下载| 国产精品一区二区在线不卡| 久久精品久久精品一区二区三区| av天堂在线播放| 亚洲欧美精品自产自拍| 亚洲欧美日韩另类电影网站| 在线观看免费午夜福利视频| 午夜福利,免费看| 999久久久国产精品视频| 青青草视频在线视频观看| kizo精华| 中国国产av一级| 最近最新中文字幕大全免费视频 | 成年人午夜在线观看视频| 最新的欧美精品一区二区| 一个人免费看片子| 国产精品三级大全| 黑丝袜美女国产一区| 久久精品久久久久久久性| 操美女的视频在线观看| 久久久国产精品麻豆| 国产黄色视频一区二区在线观看| 一本色道久久久久久精品综合| 亚洲av成人不卡在线观看播放网 | 婷婷色av中文字幕| 日韩中文字幕视频在线看片| 宅男免费午夜| 国产av国产精品国产| 久久人妻福利社区极品人妻图片 | 精品卡一卡二卡四卡免费| 久久天躁狠狠躁夜夜2o2o | 精品一区二区三卡| 亚洲,一卡二卡三卡| 男的添女的下面高潮视频| 美女中出高潮动态图| 亚洲国产欧美网| 一级片'在线观看视频| 一级毛片 在线播放| 久久午夜综合久久蜜桃| 男女午夜视频在线观看| 精品少妇一区二区三区视频日本电影| 国产日韩欧美亚洲二区| 亚洲精品久久成人aⅴ小说| 啦啦啦啦在线视频资源| 久久99热这里只频精品6学生| 成年女人毛片免费观看观看9 | 高潮久久久久久久久久久不卡| 首页视频小说图片口味搜索 | 久久99精品国语久久久| 亚洲少妇的诱惑av| 中文乱码字字幕精品一区二区三区| 又粗又硬又长又爽又黄的视频| 免费高清在线观看日韩| 欧美+亚洲+日韩+国产| 一本—道久久a久久精品蜜桃钙片| 岛国毛片在线播放| 老司机在亚洲福利影院| 成人影院久久| 久久久久精品人妻al黑| 久久 成人 亚洲| 欧美成狂野欧美在线观看| 久久狼人影院| 国产成人a∨麻豆精品| 天堂中文最新版在线下载| 国产片特级美女逼逼视频| 精品熟女少妇八av免费久了| 亚洲精品在线美女| 久久国产亚洲av麻豆专区| 后天国语完整版免费观看| 成人午夜精彩视频在线观看| 精品一区二区三区av网在线观看 | 天天操日日干夜夜撸| 秋霞在线观看毛片| 中文欧美无线码| 在现免费观看毛片| 天天添夜夜摸| 丁香六月欧美| 一级黄片播放器| 人人妻人人澡人人爽人人夜夜| svipshipincom国产片| 欧美日韩黄片免| 男人操女人黄网站| 一二三四在线观看免费中文在| 精品免费久久久久久久清纯 | 狠狠婷婷综合久久久久久88av| 久久久久久久久免费视频了| 国产一区二区三区综合在线观看| 亚洲国产精品一区三区| 日韩大码丰满熟妇| 久久毛片免费看一区二区三区| 欧美变态另类bdsm刘玥| 亚洲国产成人一精品久久久| 美女福利国产在线| 极品人妻少妇av视频| 肉色欧美久久久久久久蜜桃| 91成人精品电影| 亚洲一区中文字幕在线| 搡老乐熟女国产| 中文欧美无线码| 国产午夜精品一二区理论片| 老司机影院成人| 国产视频首页在线观看| 亚洲av成人精品一二三区| 男女免费视频国产| 精品一区二区三卡| 99热全是精品| 可以免费在线观看a视频的电影网站| 美女高潮到喷水免费观看| 手机成人av网站| 天堂8中文在线网| 日韩一本色道免费dvd| 国产精品.久久久| 一边摸一边做爽爽视频免费| 国产免费视频播放在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品国产三级专区第一集| 国产精品.久久久| 日韩欧美一区视频在线观看| 亚洲国产成人一精品久久久| 少妇 在线观看| 亚洲欧美日韩高清在线视频 | 后天国语完整版免费观看| 女人久久www免费人成看片| 久久久欧美国产精品| 国产无遮挡羞羞视频在线观看| 日日夜夜操网爽| 精品一区在线观看国产| 汤姆久久久久久久影院中文字幕| 日韩,欧美,国产一区二区三区| 狠狠精品人妻久久久久久综合| www.熟女人妻精品国产| 中文欧美无线码| 欧美黑人精品巨大| 超碰成人久久| 一本大道久久a久久精品| 久久久国产一区二区| 日韩av在线免费看完整版不卡| 国产在线视频一区二区| av在线老鸭窝| 99热全是精品| 欧美精品啪啪一区二区三区 | 两个人免费观看高清视频| 欧美成人午夜精品| 久久青草综合色| 日韩伦理黄色片| 亚洲七黄色美女视频| 日本五十路高清| 两个人看的免费小视频| 免费观看av网站的网址| 日本a在线网址| 黄色 视频免费看| 一本综合久久免费| 欧美激情 高清一区二区三区| 丰满饥渴人妻一区二区三| 亚洲中文日韩欧美视频| 黑人猛操日本美女一级片| 日韩伦理黄色片| av网站免费在线观看视频| 大香蕉久久网| 成年美女黄网站色视频大全免费| 欧美xxⅹ黑人| 老司机影院毛片| 精品一区二区三卡| 免费在线观看影片大全网站 | 中文字幕高清在线视频| 午夜久久久在线观看| 亚洲国产av影院在线观看| 69精品国产乱码久久久| 日韩制服骚丝袜av| 国产熟女午夜一区二区三区| 久久久久久久久久久久大奶| 国产熟女欧美一区二区| 国产一区二区三区综合在线观看| 久久毛片免费看一区二区三区| 亚洲伊人色综图| 女人被躁到高潮嗷嗷叫费观| 国产日韩欧美视频二区| 夜夜骑夜夜射夜夜干| 日日摸夜夜添夜夜爱| 国产片内射在线| 国产极品粉嫩免费观看在线| a级毛片在线看网站| 久久精品人人爽人人爽视色| 成年动漫av网址| 日韩制服丝袜自拍偷拍| 国产视频一区二区在线看| 午夜福利乱码中文字幕| 国产av国产精品国产| 亚洲欧洲精品一区二区精品久久久| 99久久精品国产亚洲精品| 国产免费现黄频在线看| 国产成人系列免费观看| 亚洲少妇的诱惑av| 另类精品久久| 在线观看免费高清a一片| 丰满迷人的少妇在线观看| 国产一区二区在线观看av| 国产一区二区三区av在线| 汤姆久久久久久久影院中文字幕| 一级,二级,三级黄色视频| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影 | 男女无遮挡免费网站观看| 亚洲精品一区蜜桃| 晚上一个人看的免费电影| 精品久久久精品久久久| 美女福利国产在线| 高潮久久久久久久久久久不卡| av天堂久久9| 高清av免费在线| 首页视频小说图片口味搜索 | 国产精品三级大全| 国产精品人妻久久久影院| 亚洲精品国产色婷婷电影| 脱女人内裤的视频| 久久久久久久久免费视频了| 校园人妻丝袜中文字幕| 欧美国产精品一级二级三级| 精品一区在线观看国产| 搡老乐熟女国产| 国产视频首页在线观看| 50天的宝宝边吃奶边哭怎么回事| 韩国高清视频一区二区三区| 大香蕉久久网| 操美女的视频在线观看| 精品第一国产精品| 中国国产av一级| 欧美日韩视频精品一区| 精品熟女少妇八av免费久了| 免费在线观看黄色视频的| 操出白浆在线播放| 亚洲天堂av无毛| 如日韩欧美国产精品一区二区三区| 超色免费av| 国产精品偷伦视频观看了| 国产精品成人在线| 久久天躁狠狠躁夜夜2o2o | 视频区欧美日本亚洲| 在线天堂中文资源库| 成在线人永久免费视频| 久久狼人影院| 亚洲av美国av| 精品亚洲成国产av| 极品少妇高潮喷水抽搐| 国产无遮挡羞羞视频在线观看| 韩国高清视频一区二区三区| 在线观看www视频免费| 飞空精品影院首页| 伊人久久大香线蕉亚洲五| 丰满饥渴人妻一区二区三| 91成人精品电影| 日韩,欧美,国产一区二区三区| 日韩熟女老妇一区二区性免费视频| 99re6热这里在线精品视频| 欧美激情 高清一区二区三区| 国产免费现黄频在线看| 女人久久www免费人成看片| 欧美激情极品国产一区二区三区| av在线老鸭窝| 日韩 欧美 亚洲 中文字幕| 丝袜喷水一区| 亚洲av片天天在线观看| 91精品伊人久久大香线蕉| 亚洲欧美一区二区三区久久| 日日摸夜夜添夜夜爱| 日日夜夜操网爽| 欧美日韩视频高清一区二区三区二| 国产精品一区二区在线不卡| 国产极品粉嫩免费观看在线| 久9热在线精品视频| av天堂久久9| av福利片在线| 精品久久久久久电影网| 十分钟在线观看高清视频www| 可以免费在线观看a视频的电影网站| e午夜精品久久久久久久| 亚洲av成人精品一二三区| 99热国产这里只有精品6| 国产欧美日韩综合在线一区二区| 香蕉国产在线看| 男人操女人黄网站| bbb黄色大片| 两人在一起打扑克的视频| 国产亚洲午夜精品一区二区久久| 91老司机精品| 欧美日韩国产mv在线观看视频| 欧美国产精品va在线观看不卡| 亚洲欧美精品综合一区二区三区| 999久久久国产精品视频| 一个人免费看片子| 深夜精品福利| 国产精品一区二区在线观看99| 日本a在线网址| 国产成人影院久久av| 国产精品秋霞免费鲁丝片| 成人国产av品久久久| 亚洲精品美女久久久久99蜜臀 | 国产又爽黄色视频| 美女午夜性视频免费| 国产一卡二卡三卡精品| 99精国产麻豆久久婷婷| 汤姆久久久久久久影院中文字幕| 色精品久久人妻99蜜桃| 女性生殖器流出的白浆| 欧美激情 高清一区二区三区| 久久狼人影院| 精品人妻1区二区| 建设人人有责人人尽责人人享有的| 丝袜美足系列| 高潮久久久久久久久久久不卡| 婷婷丁香在线五月| 中文欧美无线码| 久久天躁狠狠躁夜夜2o2o | 中文字幕人妻丝袜制服| 国产一卡二卡三卡精品| av在线app专区| 妹子高潮喷水视频| 黑人欧美特级aaaaaa片| 国产免费视频播放在线视频| 亚洲欧美中文字幕日韩二区| 成人国产一区最新在线观看 | 国产精品一区二区在线观看99| 美女脱内裤让男人舔精品视频| 天天操日日干夜夜撸| 又大又黄又爽视频免费| bbb黄色大片| 男女国产视频网站| 精品国产国语对白av| 久热爱精品视频在线9| 中文字幕亚洲精品专区| av天堂在线播放| 免费在线观看完整版高清| 国产精品久久久久久人妻精品电影 | 日本午夜av视频| 精品一区二区三区av网在线观看 | 国产野战对白在线观看| 亚洲国产日韩一区二区| 人成视频在线观看免费观看| 精品熟女少妇八av免费久了| 在线观看国产h片| 精品卡一卡二卡四卡免费| 欧美成人午夜精品| 国产亚洲精品第一综合不卡| 免费在线观看完整版高清| 1024视频免费在线观看| 成人国语在线视频| 高清不卡的av网站| 欧美亚洲日本最大视频资源| av线在线观看网站| 国产主播在线观看一区二区 | 亚洲人成网站在线观看播放| 成年人免费黄色播放视频| 天堂中文最新版在线下载| 大香蕉久久成人网| 亚洲成国产人片在线观看| 91精品国产国语对白视频| 国产精品99久久99久久久不卡| 性少妇av在线| 久久久久久久久免费视频了| 免费日韩欧美在线观看| 黄色 视频免费看| 丝袜在线中文字幕| 80岁老熟妇乱子伦牲交| 国产又色又爽无遮挡免| 亚洲成av片中文字幕在线观看| 午夜91福利影院| 欧美日韩亚洲综合一区二区三区_| 另类精品久久| 黄色a级毛片大全视频| 亚洲,一卡二卡三卡| 国产亚洲精品第一综合不卡| 亚洲九九香蕉| 欧美人与性动交α欧美精品济南到| 国产精品三级大全| 亚洲av片天天在线观看| 蜜桃国产av成人99| 亚洲精品国产区一区二| 国产无遮挡羞羞视频在线观看| 日韩精品免费视频一区二区三区| 女人久久www免费人成看片| 日韩一本色道免费dvd| 人人妻人人澡人人看| 国产高清videossex| 在线观看免费高清a一片| 亚洲欧美一区二区三区国产| av线在线观看网站| 尾随美女入室| 岛国毛片在线播放| 在线观看免费日韩欧美大片| 中文精品一卡2卡3卡4更新| 91精品国产国语对白视频| 久久人人爽人人片av| 女人高潮潮喷娇喘18禁视频| 性高湖久久久久久久久免费观看| 精品国产一区二区久久| 又大又黄又爽视频免费| 成人影院久久| 中文字幕高清在线视频| 美女国产高潮福利片在线看| 国产在线观看jvid| 久久人妻福利社区极品人妻图片 | 中文欧美无线码| 国产成人欧美在线观看 | 99久久精品国产亚洲精品| 免费av中文字幕在线| 久久青草综合色| 国产亚洲av高清不卡| avwww免费| 老鸭窝网址在线观看| 国产精品亚洲av一区麻豆| 午夜91福利影院| 极品少妇高潮喷水抽搐| 五月开心婷婷网| 美女主播在线视频| 中文乱码字字幕精品一区二区三区| 一级黄色大片毛片| 日韩精品免费视频一区二区三区| 久久精品国产亚洲av涩爱| 成人手机av| 午夜久久久在线观看| 欧美精品啪啪一区二区三区 | 色播在线永久视频| 又黄又粗又硬又大视频| 男女午夜视频在线观看| 叶爱在线成人免费视频播放| 国产成人啪精品午夜网站| 精品欧美一区二区三区在线| 2018国产大陆天天弄谢| 水蜜桃什么品种好| 午夜福利在线免费观看网站| svipshipincom国产片| 国产野战对白在线观看| 嫁个100分男人电影在线观看 | 在线观看人妻少妇| 免费人妻精品一区二区三区视频| 亚洲中文日韩欧美视频| 91麻豆精品激情在线观看国产 | 老汉色∧v一级毛片| 久久精品国产亚洲av高清一级| 中文字幕av电影在线播放| 男女无遮挡免费网站观看| 99精品久久久久人妻精品| 2021少妇久久久久久久久久久| 日本黄色日本黄色录像| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 免费看十八禁软件| 最近手机中文字幕大全| 亚洲精品美女久久av网站| 2018国产大陆天天弄谢| 自线自在国产av| 一区二区三区精品91| 看免费av毛片| 精品视频人人做人人爽| 国产国语露脸激情在线看| 建设人人有责人人尽责人人享有的| 久久人妻熟女aⅴ| 国产成人免费无遮挡视频| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 黄网站色视频无遮挡免费观看| 国产亚洲精品久久久久5区| 9色porny在线观看| 青春草视频在线免费观看| 视频区图区小说| 亚洲熟女毛片儿| 国产成人精品久久二区二区91| 国产成人a∨麻豆精品| 蜜桃在线观看..| 一边亲一边摸免费视频| 亚洲精品一区蜜桃| 久9热在线精品视频| 水蜜桃什么品种好| 国产成人精品久久久久久| 精品高清国产在线一区| 亚洲久久久国产精品| 在线观看国产h片| 久久精品aⅴ一区二区三区四区| 一本色道久久久久久精品综合| 久久久久久免费高清国产稀缺| 一级片免费观看大全| 免费在线观看视频国产中文字幕亚洲 | 午夜福利视频精品| 国产成人一区二区在线| 2018国产大陆天天弄谢| 十八禁高潮呻吟视频| 黄色片一级片一级黄色片| 少妇粗大呻吟视频| 免费高清在线观看视频在线观看| 老司机在亚洲福利影院| 你懂的网址亚洲精品在线观看| 免费在线观看日本一区| 国产伦理片在线播放av一区| 熟女少妇亚洲综合色aaa.| 少妇精品久久久久久久| 一级毛片 在线播放| 中文字幕制服av| 久久性视频一级片| 久热爱精品视频在线9| 纵有疾风起免费观看全集完整版| 一本综合久久免费| 90打野战视频偷拍视频| 免费一级毛片在线播放高清视频 | 亚洲国产看品久久| 欧美黄色淫秽网站| 老熟女久久久| av福利片在线| 欧美成人精品欧美一级黄| 午夜av观看不卡| 99热全是精品| 汤姆久久久久久久影院中文字幕| 国产精品免费视频内射| 男的添女的下面高潮视频| 国产在视频线精品| 免费高清在线观看日韩| 国产成人a∨麻豆精品| 啦啦啦 在线观看视频| xxxhd国产人妻xxx| 老汉色av国产亚洲站长工具| 青青草视频在线视频观看| 亚洲第一青青草原| 欧美精品av麻豆av| 国产一卡二卡三卡精品| 国产片内射在线| 赤兔流量卡办理| 夫妻午夜视频| 考比视频在线观看| 国产精品 欧美亚洲| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 在线观看人妻少妇| 亚洲成国产人片在线观看| 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| 精品亚洲成国产av| 久久久久视频综合| 国产成人系列免费观看| 777米奇影视久久| 精品亚洲成国产av| 国产一区亚洲一区在线观看| 午夜影院在线不卡| 欧美日韩亚洲国产一区二区在线观看 | 黄色片一级片一级黄色片| 好男人电影高清在线观看| 日本一区二区免费在线视频| 午夜精品国产一区二区电影| 国产成人一区二区三区免费视频网站 | 亚洲专区国产一区二区| 日韩av在线免费看完整版不卡| 亚洲欧美日韩另类电影网站|