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

    基于多目標(biāo)模糊規(guī)劃的灌區(qū)多水源優(yōu)化配置

    2017-07-31 20:54:59肖圓圓李天霄
    關(guān)鍵詞:生育期水量流量

    付 強(qiáng) 肖圓圓 崔 嵩 劉 東 李天霄

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

    基于多目標(biāo)模糊規(guī)劃的灌區(qū)多水源優(yōu)化配置

    付 強(qiáng) 肖圓圓 崔 嵩 劉 東 李天霄

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

    以黑龍江省和平灌區(qū)為例,構(gòu)建基于多目標(biāo)模糊規(guī)劃的灌區(qū)多水源優(yōu)化配置模型,該模型能夠在提高農(nóng)業(yè)灌溉用水凈效益的同時有效減少農(nóng)業(yè)灌溉水量,促進(jìn)和平灌區(qū)多水源高效配置。采用具有非線性隸屬度函數(shù)的模糊多目標(biāo)規(guī)劃求解模型,得到不同流量不同水源下的最優(yōu)配水方案。結(jié)果表明:不同流量水平下水稻不同生育階段均存在缺水現(xiàn)象,低流量下需從柳河水庫引入外調(diào)水才能保證水稻的最小需水量。為保證灌區(qū)整體效益,按照引水工程、提水工程、井灌工程的先后順序進(jìn)行配水,并得到多目標(biāo)配水模型在不同情景下的運(yùn)行穩(wěn)定情況。該模型可以高效地進(jìn)行灌區(qū)多水源在作物各生育階段的優(yōu)化配置。

    多目標(biāo)模糊規(guī)劃;多水源;優(yōu)化配置;灌區(qū)

    引言

    如何協(xié)調(diào)好灌區(qū)水資源的供需關(guān)系,確保水資源合理高效利用,是水資源領(lǐng)域研究的熱點(diǎn)問題[1]。灌區(qū)多水源優(yōu)化配置可以理解為基于作物需水量和不同水源的可供水量,通過調(diào)整不同水源的配水量,以實(shí)現(xiàn)農(nóng)業(yè)灌溉用水與農(nóng)業(yè)收益統(tǒng)籌提高的目的[2]。黑龍江省是我國重要的商品糧基地,糧食作物的播種面積逐年增長,糧食生產(chǎn)消耗的灌溉水量也逐年增加[3-4],如何科學(xué)地優(yōu)化且合理地配置有限的灌溉水資源對提高黑龍江省灌區(qū)水資源利用效率和保障糧食安全具有十分重要的意義[5]。

    近年來,國內(nèi)外關(guān)于灌區(qū)水資源優(yōu)化配置的研究不斷涌現(xiàn),其中灌區(qū)水資源多目標(biāo)配置由于能夠?qū)Χ鄠€相互矛盾的目標(biāo)進(jìn)行科學(xué)、合理的優(yōu)選而備受關(guān)注[6-9]。常見的多目標(biāo)優(yōu)化方法有目標(biāo)法[10]、加權(quán)和法[11]、理想點(diǎn)法[12]、模糊優(yōu)選理論[13]、模糊多目標(biāo)規(guī)劃法[14]、遺傳算法[15-18]、人工魚群算法[19]、粒子群算法[20-21]等。模糊多目標(biāo)規(guī)劃法由于較強(qiáng)的操作性和靈活性至今仍被廣泛應(yīng)用于各領(lǐng)域[22-23],但多數(shù)研究只探討線性模糊隸屬度函數(shù)的解法,忽略了隸屬度函數(shù)的非線性特點(diǎn)對優(yōu)化結(jié)果的影響。本文在前人研究的基礎(chǔ)上,構(gòu)建同時考慮灌溉效益和水資源利用率的灌區(qū)多水源優(yōu)化配置模型,并采用基于非線性隸屬度函數(shù)的模糊多目標(biāo)規(guī)劃方法求解模型,以期實(shí)現(xiàn)灌區(qū)多流量多水源在作物不同生育階段內(nèi)的高效配水。

    1 原理介紹

    1.1 模型建立

    灌區(qū)水資源優(yōu)化配置中,管理者需要將灌區(qū)不同水源有限的水資源量在作物各個生育階段內(nèi)進(jìn)行優(yōu)化配置,從而達(dá)到用盡可能少的水量獲得較高的產(chǎn)量或收益的目的。在優(yōu)化配置中,需考慮不同水源的可供給量、水量平衡、供水目標(biāo)、作物需水量等約束,為此,需構(gòu)建基于多目標(biāo)模糊規(guī)劃的灌區(qū)多水源優(yōu)化配置模型。所構(gòu)建模型包括2個目標(biāo)函數(shù),目標(biāo)函數(shù)1為灌區(qū)作物凈效益最大,目標(biāo)函數(shù)2為灌溉用水量最小,目標(biāo)函數(shù)1和目標(biāo)函數(shù)2的結(jié)合可促進(jìn)灌區(qū)用水效率的提升。模型表達(dá)式為:

    約束條件:

    式中 P——作物市場單價,元/kg

    WP——灌溉水分生產(chǎn)率,kg/m3

    I——供水工程類別總數(shù),分別代表引水工程、提水工程、井灌工程,取3

    T——配水周期(即作物生育期)類別總數(shù),分別代表分蘗期、拔節(jié)期、抽穗期、乳熟期,取4

    Xit——i供水工程在t生育期的配水量(決策變量),m3

    ui——i供水工程的供水成本,元/m3

    Qit——i供水工程在t生育期的可用水量,m3

    Rt-1、Rt-2——t-1、t-2生育期的余水量,m3

    Wit——i供水工程在t生育期的供水目標(biāo),m3

    WL——灌溉需水量下限值,m3

    WU——灌溉需水量上限值,m3

    1.2 模型求解

    采用模糊多目標(biāo)規(guī)劃對上述多目標(biāo)模型進(jìn)行求解。然而,在實(shí)際情況中,很難確定模糊多目標(biāo)規(guī)劃中隸屬度函數(shù)的性質(zhì)。不失一般性,采用非線性隸屬度函數(shù)對所構(gòu)建模型的2個目標(biāo)函數(shù)進(jìn)行表示(圖1)。

    令f表示任意的目標(biāo)函數(shù),則相應(yīng)的非線性隸屬度函數(shù)可表示為:

    越大越優(yōu)型

    越小越優(yōu)型

    圖1 非線性隸屬度函數(shù)Fig.1 Nonlinearmembership functions

    式中 μf(x)——f的隸屬度函數(shù)

    fU、fL——目標(biāo)函數(shù)的最大、最小值

    β1、β2——非線性隸屬度函數(shù)的形狀系數(shù)

    β1、β2=1表示線性,β1、β2>1和 0<β1、β2<1表示非線性[24]。引入變量λ,則原模型的目標(biāo)函數(shù)可通過隸屬度函數(shù)的形式,即(μfn(X))βn≥λ(0≤λ≤1)轉(zhuǎn)換成約束條件,與原有約束條件一起構(gòu)成一個目標(biāo)函數(shù)為maxλ的單目標(biāo)模型。對于本文所構(gòu)建模型,目標(biāo)函數(shù)1為經(jīng)濟(jì)效益最大,屬于越大越優(yōu)型,其隸屬度函數(shù)采用式(8)計(jì)算,目標(biāo)函數(shù)2為灌溉用水量最小,屬于越小越優(yōu)型,其隸屬度函數(shù)采用式(9)計(jì)算。根據(jù)上述原理,所構(gòu)建多目標(biāo)模型可轉(zhuǎn)換成如下單目標(biāo)模型形式:

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

    約束條件

    式中 F1max、F1min——目標(biāo)函數(shù)1的單目標(biāo)模型目標(biāo)函數(shù)的最大值和最小值

    F2max、F2min——目標(biāo)函數(shù)2的單目標(biāo)模型目標(biāo)函數(shù)的最大值和最小值

    2 實(shí)例研究

    2.1 研究區(qū)概況

    和平灌區(qū)位于黑龍江省綏化市慶安縣中部,呼蘭河左岸的干支流河漫灘及一級階地上,灌區(qū)范圍由東向西呈帶狀分布。該灌區(qū)主要水源有地表水與地下水,取水方式又可分為3種供水工程,即引水工程、提水工程和井灌工程,井灌工程由地下水提供,其他2個供水工程由地表水提供,柳河水庫為外調(diào)水源。水稻是和平灌區(qū)的主要作物[25],其生育期主要包括返青期、分蘗期、拔節(jié)期、抽穗期、乳熟期、黃熟期,其中返青期和黃熟期的騰發(fā)量都較小,因此,本文根據(jù)水稻生育期需水量需求,將3種供水工程的水量優(yōu)化分配到水稻的分蘗期、拔節(jié)期、抽穗期、乳熟期4個生育期。

    2.2 基礎(chǔ)數(shù)據(jù)確定

    根據(jù)《呼蘭河灌區(qū)工程初期設(shè)計(jì)報(bào)告》以及當(dāng)?shù)厮畡?wù)局提供的調(diào)研數(shù)據(jù),對水稻不同生育階段需水量以及各供水工程灌溉控制面積進(jìn)行分析,得到各供水工程的供水目標(biāo)以及水稻充分灌溉條件下的需水量上、下限值;綜合分析灌區(qū)內(nèi)多年降水量和徑流統(tǒng)計(jì)資料,獲得不同流量水平下水稻各生育期各供水工程的可用水量,基礎(chǔ)數(shù)據(jù)見表1和表2[26]。

    表1 不同生育階段各供水工程的配水目標(biāo)與需水量臨界值Tab.1 Water supply targets and water demand thresholds for water supply projects at different grow th stages

    表2 不同生育階段各供水工程的可用水量Tab.2 Available water of water supp ly projects at different grow th stages

    3 模型結(jié)果與分析

    3.1 隸屬度為線性情況

    根據(jù)上述原理,運(yùn)用LINGO11對模型進(jìn)行求解。得到不同β1和β2組合下不同流量3種供水工程在不同生育期的最優(yōu)配水量,以β1=1、β2=1,即隸屬度函數(shù)為線性情況為例,模型運(yùn)行結(jié)果見表3。從表中可以看出,3個水源中,引水工程和提水工程在4個生育期內(nèi)都有配水,而井灌工程則為部分生育期不予配水,高流量水平下,只有在拔節(jié)期井灌工程才有配水。這是由于井灌工程的供水費(fèi)用在3個供水工程中是最高的,在能夠保證作物各生育期最低需水量的前提下,水量分配順序?yàn)閮?yōu)先采用引水工程,其次采用提水工程,最后采用井灌工程。另外,模型運(yùn)行結(jié)果是作物凈效益最大和灌溉用水量最小2個目標(biāo)綜合作用的結(jié)果,以尋求用水效益和用水量之間的平衡,達(dá)到綜合效益最優(yōu)。以中流量為例,若單純追求作物凈效益最大,則3個供水工程在4個生育期均有分配,若只考慮灌溉用水量最小,則提水工程和井灌工程在抽穗期和乳熟期均未配水,這2個生育期的最小需水量僅由引水工程來滿足,綜合考慮2個目標(biāo)函數(shù),提水工程在穗期和乳熟期有配水,但配水量較單純考慮經(jīng)濟(jì)效益最大目標(biāo)下有所調(diào)整,井灌工程在這2個生育期則不配水。高流量、中流量的總配水量分別為 1 860萬 m3、1 730萬m3,低流量下由于3個水源的總可供水量滿足不了水稻的最小需水量要求,需要從柳河水庫調(diào)水190萬m3,調(diào)水后,低流量的最優(yōu)總配水量為1 643萬m3。圖2(圖中S1、S2、S3、S4分別代表分蘗期、拔節(jié)期、抽穗期、乳熟期)為各流量水平下作物各生育期的缺水量,即平均需水量與最優(yōu)配水量之間的差值。3個流量水平下均存在缺水現(xiàn)象,高流量、中流量和低流量下總?cè)彼糠謩e為315萬m3、445萬m3和532萬m3。

    表3 β1=1和β2=1組合情況下模型最優(yōu)配水結(jié)果Tab.3 Results of optimalwater distribution under β1=1 andβ2=1

    圖2 不同流量水平下不同生育階段的配水量與平均需水量Fig.2 Water allocation and average water demand under different flow levels and growth stages

    3.2 隸屬度為非線性情況

    上述分析了β1=1和β2=1,即模糊多目標(biāo)規(guī)劃中隸屬度函數(shù)為線性的配水結(jié)果,實(shí)際上,目標(biāo)函數(shù)非線性隸屬度函數(shù)也是常見的,如圖1所示。本文將β1和β2進(jìn)行如下組合來分析不同情景下配水方案的變化。情景1:令β1=1,β2=0.3,0.5,0.7,1,3,7,10依次變化;情景2:令β2=1,β1=0.3,0.5,0.7,1,3,7,10依次變化。圖3(圖中P1、P2、P3分別代表引水工程、提水工程和井灌工程)為2種情景下3種水源在不同流量下的最優(yōu)配水量,圖3a為情景1的結(jié)果,圖3b為情景2的結(jié)果。對于情景1,高流量水平下的引水工程配水量保持不變,提水工程配水量在β1為0.3~7之間有略微下降趨勢,井灌工程在這期間配水量均為零。但在β2=10時發(fā)生顯著變化,提水工程和井灌工程配水量均有大幅度提升,表明對于高流量水平,當(dāng)β1=1、β2=10情景下模型穩(wěn)定性較差,決策者在進(jìn)行決策時,應(yīng)避免這種情景。中流量和低流量下的引水工程和提水工程配水量均維持穩(wěn)定狀態(tài),井灌工程配水量有略微調(diào)整,但變化不顯著。對于情景2,在β2=1下,從β1=7開始3個流量水平下的3個供水工程的配水量均發(fā)生顯著變化,表明在這種情景下,β1>7時模型穩(wěn)定性發(fā)生變化。在β1<7的各組合情景中,高流量的引水工程和井灌工程的配水量保持不變,提水工程配水量隨β1的增大有增大趨勢;中流量下引水工程配水量保持不變,提水工程和井灌工程配水量隨β1的增大有增大趨勢;低流量下引水工程和提水工程配水量保持不變,井灌工程的配水量隨β1的增大有增大趨勢,上述結(jié)果表明,引水工程的配水量是最穩(wěn)定的,井灌工程的配水量最不穩(wěn)定,尤其是在低流量水平下具有較大的敏感性。以中流量為例,圖4為不同情景組合不同流量水平下不同生育期的配水量變化情況,β1=1時,拔節(jié)期在任意β2情況下的配水量均保持不變,分蘗期、抽穗期和乳熟期在β2≤1的情況下略微變化,在β2>1的情況下趨于穩(wěn)定。β2=1時,分蘗期和拔節(jié)期的配水量均保持穩(wěn)定,抽穗期和乳熟期在β1<7時略微變化,β1>7時趨于穩(wěn)定。不同的情景組合會導(dǎo)致不同的配水量結(jié)果,決策者可根據(jù)實(shí)際情況結(jié)合個人偏好做出最終決策。

    圖3 不同情景組合下3種供水水源在不同流量水平下的配水量變化情況Fig.3 Water allocation changes under differentwater sources,flow levels and scenarios

    圖4 不同情景組合不同流量水平下不同生育期的配水量變化情況Fig.4 Water allocation changes under different growth stages and scenarios

    4 結(jié)論

    (1)根據(jù)和平灌區(qū)水源特點(diǎn),構(gòu)建多水源多目標(biāo)配水模型,將不同水源(包括引水工程、提水工程和井灌工程)的可供水量在水稻各生育期內(nèi)進(jìn)行高效配置,以達(dá)到用盡可能少的灌溉水資源量獲得盡可能大的經(jīng)濟(jì)效益的目的。

    (2)采用模糊多目標(biāo)規(guī)劃方法對模型進(jìn)行求解,結(jié)果顯示灌區(qū)按照引水工程、提水工程、井灌工程的先后順序進(jìn)行配水。模型得到了灌區(qū)用水效益和用水量之間的平衡,在綜合效益最優(yōu)情況下,3個流量水平下均存在缺水現(xiàn)象,低流量情況下需從柳河水庫調(diào)入190萬m3水才能保證水稻的最小需水量。

    (3)重點(diǎn)探討了多目標(biāo)模型目標(biāo)函數(shù)的模糊不確定性,用非線性隸屬度函數(shù)進(jìn)行量化,通過模型目標(biāo)函數(shù)不同非線性隸屬度函數(shù)的情景組合,討論了模型在不同情景下的運(yùn)行穩(wěn)定情況。引水工程的配水量相對穩(wěn)定,而井灌工程在低流量水平下對目標(biāo)函數(shù)隸屬度函數(shù)的非線性程度具有較大的敏感性。

    參 考 文 獻(xiàn)

    1 姜秋香,付強(qiáng),朱長虹,等.基于多維臨界調(diào)控理論的水資源優(yōu)化配置研究進(jìn)展[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2015,46(7):103-108.JIANG Qiuxiang,F(xiàn)U Qiang,ZHU Changhong,et al.Research progress of water resources optimal allocation based on theory of multidimensional critical regulation and control[J].Journal of Northeast Agricultural University,2015,46(7):103-108.(in Chinese)

    2 李晨洋,張志鑫.基于區(qū)間兩階段模糊隨機(jī)模型的灌區(qū)多水源優(yōu)化配置[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(12):107-114.LIChenyang,ZHANG Zhixin.Multi-water conjunctive optimal allocation based on interval-parameter two-stage fuzzy-stochastic programming[J].Transactions of the CSAE,2016,32(12):107-114.(in Chinese)

    3 付強(qiáng),劉巍,劉東,等.黑龍江省灌溉用水效率指標(biāo)體系空間格局研究[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(12):127-132.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.a(chǎn)spx?flag=1&file_no=20151218&journal_id=jcsam.DOI:10.6041/ j.issn.1000-1298.2015.12.018.FU Qiang,LIUWei,LIU Dong,etal.Spatial distribution of irrigation water use efficiency index system in Heilongjiang Province[J/ OL].Transactions of the Chinese Society for Agricultural Machinery,2015,46(12):127-132.(in Chinese)

    4 付強(qiáng),劉燁,李天霄,等.水足跡視角下的黑龍江省糧食生產(chǎn)用水分析[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(6):184-192.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.a(chǎn)spx?flag=1&file_no=20170624&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2017.06.024.FU Qiang,LIU Ye,LI Tianxiao,et al.Analysis of water utilization in grain production from water footprint perspective in Heilongjiang Province[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2017,48(6):184-192.(in Chinese)

    5 李茉,姜瑤,郭萍,等.考慮不同層次利益主體的灌溉水資源優(yōu)化配置[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(5):199-207.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.a(chǎn)spx?flag=1&file_no=20170225&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2017.05.025.LIMo,JIANG Yao,GUO Ping,et al.Irrigation water optimal allocation considering stakeholders of different levels[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2017,48(5):199-207.(in Chinese)

    6 張智韜,劉俊民,陳俊英,等.基于RS、GIS和蟻群算法的多目標(biāo)渠系配水優(yōu)化[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,41(11):72-78.ZHANG Zhitao,LIU Junmin,CHEN Junying,et al.Water resources allocation of canal system based on multi-objective about RS,GISand ant colony algorithm[J].Transactions of the Chinese Society for Agricultural Machinery,2010,41(11):72-78.(in Chinese)

    7 聶相田,邱林,周波,等.井渠結(jié)合灌區(qū)水資源多目標(biāo)優(yōu)化配置模型與應(yīng)用[J].節(jié)水灌溉,2006(4):26-31.NIE Xiangtian,QIU Lin,ZHOU Bo,et al.A multi-objective optimal allocation model for water resources in well-canal combined irrigated area and its application[J].Water Saving Irrigation,2006(4):26-31.(in Chinese)

    8 趙恩龍.基于多目標(biāo)遺傳算法的灌區(qū)水資源優(yōu)化調(diào)度研究[D].武漢:長江科學(xué)院,2013.ZHAO Enlong.Study of optimalwater resources allocation of irrigation district based on multi-objective genetic algorithm evolution[D].Wuhan:Changjiang River Scientific Research Institute,2013.(in Chinese)

    9 LIM,GUO P,SING V P.Biobjective optimization for efficient irrigation under fuzzy uncertainty[J].Journal of Irrigation and Drainage Engineering,2016,142(8):05016003.

    10 吳慧,呂一兵.一類非線性二層多目標(biāo)規(guī)劃問題的主要目標(biāo)法[J].西南大學(xué)學(xué)報(bào):自然科學(xué)版,2016,38(1):109-115.WU Hui,LYibing.Themain targetmethod for solving a class of nonlinear bilevelmulti-objective programming problems[J].Journal of Southwest University:Natural Science Edition,2016,38(1):109-115.(in Chinese)

    11 喬辰,張國立.幾何加權(quán)法求解多目標(biāo)規(guī)劃問題[J].華北電力大學(xué)學(xué)報(bào):自然科學(xué)版,2011,38(6):107-110.QIAO Chen,ZHANG Guoli.Geometric weighting method for solving multi-objective programming problems[J].Journal of North China Electric Power University:Natural Science Edition,2011,38(6):107-110.(in Chinese)

    12 龔懿,張禮華.基于理想點(diǎn)法的灌區(qū)多目標(biāo)種植結(jié)構(gòu)優(yōu)化模式研究[J].現(xiàn)代農(nóng)業(yè)科技,2014(13):211-218.GONG Yi,ZHANG Lihua.Research on multi-objective mode optimization of planting structure in irrigation area based on ideal pointmethod[J].Modern Agricultural Science and Technology,2014(13):211-218.(in Chinese)

    13 黎坤,陳曉宏,江濤,等.多目標(biāo)系統(tǒng)模糊優(yōu)選理論在城市飲用水源地選址中的應(yīng)用[J].中山大學(xué)學(xué)報(bào):自然科學(xué)版,2005,44(4):120-123.LIKun,CHEN Xiaohong,JIANG Tao,etal.Application ofmulti-objective system fuzzy optimization theory to selecting table-water resources in city[J].Acta Scientlarum Naturalium Universitatis Sunyatseni,2005,44(4):120-123.(in Chinese)

    14 鄒進(jìn),張勇傳.三峽梯級電站短期優(yōu)化調(diào)度的模糊多目標(biāo)動態(tài)規(guī)劃[J].水利學(xué)報(bào),2005,36(8):925-931.ZOU Jin,ZHANG Yongchuan.Short-term optimal scheduling of cascade hydropower stations using fuzzy multi-objective dynamic programming[J].Journal of Hydraulic Engineering,2005,36(8):925-931.(in Chinese)

    15 孫月峰,張勝紅,王曉玲,等.基于混合遺傳算法的區(qū)域大系統(tǒng)多目標(biāo)水資源優(yōu)化配置模型[J].系統(tǒng)工程理論與踐,2009,29(1):139-144.SUN Yuefeng,ZHANG Shenghong,WANG Xiaoling,et al.Multi-objective optimization of regionalwater resources based on mixed genetic algorithm[J].Systems Engineering—Theory&Practice,2009,29(1):139-144.(in Chinese)

    16 游進(jìn)軍,紀(jì)昌明,付湘.基于遺傳算法的多目標(biāo)問題求解方法[J].水利學(xué)報(bào),2003,34(7):64-69.YOU Jinjun,JIChangming,F(xiàn)U Xiang.New method for solvingmulti-objective problem based on genetic algorithm[J].Journal of Hydraulic Engineering,2003,34(7):64-69.(in Chinese)17 陳南祥,李躍鵬,徐晨光.基于多目標(biāo)遺傳算法的水資源優(yōu)化配置[J].水利學(xué)報(bào),2006,37(3):308-313.CHEN Nanxiang,LIYuepeng,XU Chenguang.Optimal deployment of water resources based on multi-objective genetic algorithm[J].Journal of Hydraulic Engineering,2006,37(3):308-313.(in Chinese)

    18 王鵬.基于Pareto Front的多目標(biāo)遺傳算法在灌區(qū)水資源配置中的應(yīng)用[J].節(jié)水灌溉,2005(6):29-32.

    19 高玉芳,張展羽.混沌人工魚群算法及其在灌區(qū)優(yōu)化配水中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(6):7-11.GAO Yufang,ZHANG Zhanyu.Chaotic artificial fish-swarm algorithm and its application in water use optimization irrigated areas[J].Transactions of the CSAE,2007,23(6):7-11.(in Chinese)

    20 張文明,董增川,朱成濤,等.基于粒子群算法的水文模型參數(shù)多目標(biāo)優(yōu)化研究[J].水利學(xué)報(bào),2008,39(5):528-534.ZHANGWenming,DONG Zengchuan,ZHU Chengtao,et al.Automatic calibration of hydrologic model based on multi-objective particle swarm optimization method[J].Journal of Hydraulic Engineering,2008,39(5):528-534.(in Chinese)

    21 周建中,李英海,肖舸,等.基于混合粒子群算法的梯級水電站多目標(biāo)優(yōu)化調(diào)度[J].水利學(xué)報(bào),2010,41(10):1212-1219.ZHOU Jianzhong,LIYinghai,XIAO Ge,et al.Multi-objective optimal dispatch of cascade hydropower stations based on shuffled particle swarm operation algorithm[J].Journal of Hydraulic Engineering,2010,41(10):1212-1219.(in Chinese)

    22 MORANKAR D V,RAJU K S,KUMAR D N.Integrated sustainable irrigation planning with multiobjective fuzzy optimization approach[J].Water Resource Management,2013,27(11):3981-4004.

    23 TEEGAVARAPU R SV,F(xiàn)ERREIRA A R,SIMONOVIC SP.Fuzzymultiobjectivemodels for optimal operation of a hydropower system[J].Water Resource Research,2013,49(6):3180-3193.

    24 SASIKUMAR K,MUJUMDAR PP.Fuzzy optimization model forwater qualitymanagement of a river system[J].Journal ofWater Resources Planning and Management,1998,124(2):79-88.

    25 李丹,李德麗,顏華.用水農(nóng)戶參與水利灌溉管理的沖量過程分析[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2013,44(8):116-121.LIDan,LIDeli,YAN Hua.Impulse process analysis ofwater users participate in irrigationmanagement[J].Journal of Northeast Agricultural University,2013,44(8):116-121.(in Chinese)

    26 付強(qiáng),劉銀鳳,劉東,等.基于區(qū)間多階段隨機(jī)規(guī)劃模型的灌區(qū)多水源優(yōu)化配置[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(1):132-139.FU Qiang,LIU Yinfeng,LIU Dong,et al.Optimal allocation ofmulti-water resources in irrigation area based on interval-parameter multi-stage stochastic programmingmodel[J].Transactions of the CSAE,2016,32(1):132-139.(in Chinese)

    Multi-water Resources Optimal Allocation of Irrigation District Based on Fuzzy Multi-objective Programm ing

    FU Qiang XIAO Yuanyuan CUISong LIU Dong LITianxiao
    (School ofWater Conservancy and Civil Engineering,Northeast Agricultural University,Harbin 150030,China)

    Many problems,such as the shortage of water resources and the low operating efficiency of water project,are ubiquitous in an irrigation district with multi-water resources.Therefore,it is significant and necessary to implement joint dispatching of multi-water resources to alleviate water shortage and promote the synergy ascension of water resources productivity.Taking Heping irrigation district in Heilongjiang Province as an example,a fuzzy multi-objective programming for multi-water resources optimal allocation model was constructed.The proposed model was capable of increasing water net benefit and decreasingwater-use amount simultaneously,which would improve the effective allocation ofmulti-water resources of Heping irrigation district.Fuzzy multi-objective programming with nonlinear membership functions was adopted to solve the model,and the optimal water allocation schemes under different flow levels and water sourceswere obtained.Results showed thatwater shortage existed in each growth period of rice under each flow level,and more water should be transferred from Liuhe reservoir under low flow level to guarantee the minimum water requirement.Water were allocated to rice in the order ofwater diversion project,water lifting project and well irrigation project to ensure the total benefit of Heping irrigation district.The stability of themulti-objectivemodel operation under different scenarios was obtained and analyzed in detail.Water allocation of water diversion project was relatively stable,while water allocation of well irrigation project had higher sensitivity for nonlinearity of membership function of the objective function under low flow level.The developed model can effectively allocate multi-water resources to the different growth stages of rice,which would provide reasonable decision schemes for the corresponding administrative departments.

    fuzzy multi-objective programming;multi-water resources;optimal allocation;irrigation district

    TV213.4;O221.6

    A

    1000-1298(2017)07-0222-06

    2016-12-14

    2017-01-17

    國家自然科學(xué)基金項(xiàng)目(51479032、51279031、51609039)、黑龍江省高校長江學(xué)者后備支持計(jì)劃項(xiàng)目、黑龍江省水利科技項(xiàng)目(201318、201503)和國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFC0406002)

    付強(qiáng)(1973—),男,教授,博士生導(dǎo)師,主要從事農(nóng)業(yè)水土資源優(yōu)化利用與管理研究,E-mail:fuqiang0629@126.com

    10.6041/j.issn.1000-1298.2017.07.028

    猜你喜歡
    生育期水量流量
    冰墩墩背后的流量密碼
    玩具世界(2022年2期)2022-06-15 07:35:36
    小水量超純水制備系統(tǒng)的最佳工藝選擇
    張曉明:流量決定勝負(fù)!三大流量高地裂變無限可能!
    尋找書業(yè)新流量
    出版人(2020年4期)2020-11-14 08:34:26
    基于水力壓裂鉆孔的注水量及壓裂半徑的應(yīng)用研究
    基于作物生育期的潛在蒸散的時空演變特征及R/S 分析
    西藏科技(2015年5期)2015-09-26 11:55:25
    分散藥包千噸注水量的水壓爆破
    五位一體流量平穩(wěn)控制系統(tǒng)
    壓縮機(jī)級間冷卻冷凝析水量計(jì)算
    2013-2014年度二二二團(tuán)冬小麥各生育期氣象條件分析
    一级毛片女人18水好多| 久久久久网色| 18禁观看日本| 亚洲精品粉嫩美女一区| 成人黄色视频免费在线看| 男女高潮啪啪啪动态图| 欧美精品亚洲一区二区| 成年动漫av网址| 国产伦人伦偷精品视频| 久久精品国产99精品国产亚洲性色 | 脱女人内裤的视频| 国产精品亚洲一级av第二区| 一区二区日韩欧美中文字幕| 母亲3免费完整高清在线观看| 中文字幕色久视频| www日本在线高清视频| 国产在线一区二区三区精| 捣出白浆h1v1| tocl精华| 99国产精品99久久久久| 精品人妻1区二区| 国产亚洲欧美在线一区二区| 99热网站在线观看| 国产亚洲精品一区二区www | 精品国产一区二区三区久久久樱花| 中文字幕人妻熟女乱码| 免费在线观看日本一区| 国产精品久久久久久精品古装| 久久天堂一区二区三区四区| 岛国毛片在线播放| 天堂动漫精品| 国产欧美日韩一区二区三| 少妇精品久久久久久久| 精品亚洲成a人片在线观看| 午夜视频精品福利| 十八禁网站网址无遮挡| 99久久国产精品久久久| 亚洲欧美日韩高清在线视频 | 99香蕉大伊视频| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| 他把我摸到了高潮在线观看 | 一边摸一边做爽爽视频免费| 国产一卡二卡三卡精品| 两性夫妻黄色片| 欧美日韩成人在线一区二区| 亚洲人成77777在线视频| 757午夜福利合集在线观看| 成年女人毛片免费观看观看9 | 免费人妻精品一区二区三区视频| 夜夜夜夜夜久久久久| 免费在线观看日本一区| 岛国在线观看网站| 无遮挡黄片免费观看| 欧美激情 高清一区二区三区| 国产福利在线免费观看视频| 黑人猛操日本美女一级片| 无限看片的www在线观看| 嫩草影视91久久| 中文字幕人妻丝袜一区二区| 香蕉国产在线看| 国产片内射在线| 法律面前人人平等表现在哪些方面| 久久久久久久久免费视频了| 男女之事视频高清在线观看| 91国产中文字幕| 亚洲欧洲精品一区二区精品久久久| 成年人午夜在线观看视频| 国产野战对白在线观看| 国产成人精品无人区| 日韩大片免费观看网站| 久久久久网色| 麻豆成人av在线观看| 免费久久久久久久精品成人欧美视频| 国产精品亚洲一级av第二区| 国产激情久久老熟女| 亚洲精华国产精华精| 欧美在线黄色| 亚洲精品在线美女| 麻豆乱淫一区二区| 怎么达到女性高潮| 99精品久久久久人妻精品| 久久人妻熟女aⅴ| 1024香蕉在线观看| 亚洲成a人片在线一区二区| 黄片大片在线免费观看| 亚洲九九香蕉| 免费少妇av软件| 国产精品久久久av美女十八| 国产日韩欧美亚洲二区| 在线观看www视频免费| 黄色视频在线播放观看不卡| 精品国产国语对白av| 一区二区三区国产精品乱码| 人人妻人人添人人爽欧美一区卜| 精品国内亚洲2022精品成人 | 狠狠精品人妻久久久久久综合| 窝窝影院91人妻| 国产黄色免费在线视频| 蜜桃在线观看..| 自线自在国产av| 亚洲色图av天堂| 久久青草综合色| 亚洲国产中文字幕在线视频| 久久精品亚洲av国产电影网| 性色av乱码一区二区三区2| 亚洲第一青青草原| 午夜福利免费观看在线| 精品久久蜜臀av无| 精品少妇一区二区三区视频日本电影| 精品人妻熟女毛片av久久网站| 成年人黄色毛片网站| 怎么达到女性高潮| 无遮挡黄片免费观看| 亚洲av日韩在线播放| 2018国产大陆天天弄谢| 麻豆乱淫一区二区| 久久精品人人爽人人爽视色| 婷婷丁香在线五月| 黑人巨大精品欧美一区二区mp4| 蜜桃在线观看..| 久久九九热精品免费| 淫妇啪啪啪对白视频| 国产无遮挡羞羞视频在线观看| 天堂动漫精品| 19禁男女啪啪无遮挡网站| 久久婷婷成人综合色麻豆| 99精品在免费线老司机午夜| 搡老乐熟女国产| 国产在线精品亚洲第一网站| 亚洲av片天天在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区在线不卡| 日韩免费av在线播放| 美女高潮到喷水免费观看| 黄色片一级片一级黄色片| 亚洲成人手机| 欧美变态另类bdsm刘玥| 欧美日韩成人在线一区二区| 精品国产国语对白av| 国产日韩欧美在线精品| 欧美成狂野欧美在线观看| 久久青草综合色| 国产在线一区二区三区精| 香蕉丝袜av| 欧美另类亚洲清纯唯美| 午夜视频精品福利| 18禁美女被吸乳视频| 免费久久久久久久精品成人欧美视频| 午夜视频精品福利| 满18在线观看网站| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三区在线| 久久久久久久精品吃奶| 色视频在线一区二区三区| 国产高清视频在线播放一区| 天堂中文最新版在线下载| 色尼玛亚洲综合影院| 天堂俺去俺来也www色官网| 99在线人妻在线中文字幕 | 国产精品 欧美亚洲| 欧美成人午夜精品| 国产免费现黄频在线看| 国产亚洲午夜精品一区二区久久| 女人被躁到高潮嗷嗷叫费观| 亚洲av日韩在线播放| 中文字幕最新亚洲高清| 涩涩av久久男人的天堂| 嫁个100分男人电影在线观看| 99精品在免费线老司机午夜| 国产精品香港三级国产av潘金莲| 免费观看a级毛片全部| 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| 女人精品久久久久毛片| 亚洲三区欧美一区| 乱人伦中国视频| 一夜夜www| 法律面前人人平等表现在哪些方面| 亚洲精品在线美女| 51午夜福利影视在线观看| 黑人猛操日本美女一级片| 少妇粗大呻吟视频| 男女下面插进去视频免费观看| 国产人伦9x9x在线观看| 精品人妻1区二区| 一级片免费观看大全| 国产成人av激情在线播放| 欧美 日韩 精品 国产| 国产精品自产拍在线观看55亚洲 | 老司机亚洲免费影院| 久久毛片免费看一区二区三区| 国产欧美日韩一区二区精品| 女性被躁到高潮视频| 日韩欧美三级三区| 亚洲国产欧美在线一区| 国产高清视频在线播放一区| 久久精品91无色码中文字幕| 亚洲av片天天在线观看| 午夜福利,免费看| 日韩一卡2卡3卡4卡2021年| 成年人午夜在线观看视频| 美女高潮喷水抽搐中文字幕| 精品久久蜜臀av无| 久久人人97超碰香蕉20202| 国产精品一区二区在线观看99| av网站在线播放免费| 精品国产亚洲在线| 天堂中文最新版在线下载| 久久人妻av系列| 亚洲 国产 在线| 亚洲国产精品一区二区三区在线| 99国产精品99久久久久| 国产精品国产av在线观看| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 欧美激情极品国产一区二区三区| 成人18禁在线播放| 丝袜喷水一区| 日本一区二区免费在线视频| 亚洲欧美激情在线| 久久人妻熟女aⅴ| 天堂8中文在线网| 丁香六月天网| 国产熟女午夜一区二区三区| 欧美日韩精品网址| 叶爱在线成人免费视频播放| 亚洲 欧美一区二区三区| 日韩视频在线欧美| 色婷婷av一区二区三区视频| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 久久久久久久久免费视频了| 国产av一区二区精品久久| 国产麻豆69| 女同久久另类99精品国产91| 久久人妻熟女aⅴ| 人人澡人人妻人| 搡老岳熟女国产| 一区在线观看完整版| 99国产精品一区二区三区| 精品少妇黑人巨大在线播放| 国产成人啪精品午夜网站| 老司机靠b影院| cao死你这个sao货| 丰满少妇做爰视频| 国产一区二区激情短视频| 免费久久久久久久精品成人欧美视频| 色在线成人网| 日本wwww免费看| 一级a爱视频在线免费观看| 视频区欧美日本亚洲| 色尼玛亚洲综合影院| av天堂在线播放| 嫁个100分男人电影在线观看| 亚洲专区字幕在线| 国产精品.久久久| 法律面前人人平等表现在哪些方面| 成人国产一区最新在线观看| 成人黄色视频免费在线看| 一边摸一边抽搐一进一出视频| 黄色视频不卡| 淫妇啪啪啪对白视频| av在线播放免费不卡| 日韩有码中文字幕| 久久久久久免费高清国产稀缺| 免费女性裸体啪啪无遮挡网站| 啦啦啦中文免费视频观看日本| 亚洲国产看品久久| 女人精品久久久久毛片| 亚洲国产成人一精品久久久| 两性夫妻黄色片| 脱女人内裤的视频| 新久久久久国产一级毛片| 嫁个100分男人电影在线观看| 国产成人av激情在线播放| 女人被躁到高潮嗷嗷叫费观| 日本av免费视频播放| 老熟妇乱子伦视频在线观看| 少妇裸体淫交视频免费看高清 | 国产精品秋霞免费鲁丝片| 亚洲精品在线美女| 桃红色精品国产亚洲av| 亚洲欧美精品综合一区二区三区| 91av网站免费观看| 女人爽到高潮嗷嗷叫在线视频| 精品少妇黑人巨大在线播放| 精品一区二区三区av网在线观看 | 久久精品国产99精品国产亚洲性色 | 亚洲成人手机| 亚洲av日韩精品久久久久久密| 国产高清国产精品国产三级| 国产高清激情床上av| av片东京热男人的天堂| 日本五十路高清| 亚洲熟女精品中文字幕| 亚洲精品中文字幕在线视频| tocl精华| 在线观看免费视频日本深夜| 日韩人妻精品一区2区三区| 一本久久精品| 免费高清在线观看日韩| 99久久99久久久精品蜜桃| 日本黄色日本黄色录像| 久久青草综合色| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 日韩欧美一区二区三区在线观看 | 亚洲精品久久午夜乱码| 日韩欧美一区二区三区在线观看 | 黄色怎么调成土黄色| 在线永久观看黄色视频| 多毛熟女@视频| 国产xxxxx性猛交| 老司机午夜十八禁免费视频| 十八禁高潮呻吟视频| 亚洲一卡2卡3卡4卡5卡精品中文| 12—13女人毛片做爰片一| 又大又爽又粗| 亚洲熟女毛片儿| 亚洲精华国产精华精| www日本在线高清视频| 成在线人永久免费视频| 亚洲成人手机| av有码第一页| 一级毛片精品| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 五月开心婷婷网| 精品国内亚洲2022精品成人 | 国产有黄有色有爽视频| 91九色精品人成在线观看| 久久av网站| 国产激情久久老熟女| 日韩视频在线欧美| 美女扒开内裤让男人捅视频| 波多野结衣av一区二区av| 日本撒尿小便嘘嘘汇集6| 老熟妇仑乱视频hdxx| 久久精品亚洲av国产电影网| 国产精品久久久久久精品古装| 这个男人来自地球电影免费观看| 99国产精品一区二区三区| 国产av精品麻豆| 色播在线永久视频| 叶爱在线成人免费视频播放| 亚洲午夜理论影院| 国内毛片毛片毛片毛片毛片| 欧美 亚洲 国产 日韩一| 日韩欧美免费精品| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 国产99久久九九免费精品| 国产精品二区激情视频| 日日夜夜操网爽| 国产一区二区三区综合在线观看| 啦啦啦中文免费视频观看日本| 一二三四社区在线视频社区8| 青青草视频在线视频观看| 日韩欧美三级三区| 怎么达到女性高潮| 十八禁网站网址无遮挡| 亚洲色图综合在线观看| 美女高潮到喷水免费观看| 啦啦啦视频在线资源免费观看| 99国产精品免费福利视频| 国产成人av激情在线播放| 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 国产精品.久久久| 午夜福利视频精品| 午夜福利在线免费观看网站| 又黄又粗又硬又大视频| cao死你这个sao货| 一级黄色大片毛片| 一区二区三区精品91| 激情视频va一区二区三区| 高清毛片免费观看视频网站 | 精品少妇久久久久久888优播| 一本一本久久a久久精品综合妖精| 搡老乐熟女国产| 欧美黄色片欧美黄色片| 久久热在线av| 亚洲熟妇熟女久久| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 丝袜美足系列| 性色av乱码一区二区三区2| 99国产精品一区二区三区| 男女免费视频国产| 黄色视频在线播放观看不卡| 嫁个100分男人电影在线观看| 黄网站色视频无遮挡免费观看| 两性夫妻黄色片| 黄色 视频免费看| 久9热在线精品视频| 久久婷婷成人综合色麻豆| av福利片在线| 乱人伦中国视频| 首页视频小说图片口味搜索| 精品一品国产午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩一区二区三| 成人影院久久| 久久久水蜜桃国产精品网| 在线观看免费高清a一片| 老司机影院毛片| 午夜91福利影院| 男女之事视频高清在线观看| 一二三四在线观看免费中文在| 高清毛片免费观看视频网站 | 99国产精品99久久久久| 精品国产乱码久久久久久男人| 99国产精品一区二区三区| 欧美精品一区二区大全| 黑人巨大精品欧美一区二区mp4| 在线观看免费日韩欧美大片| 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| 在线观看一区二区三区激情| av一本久久久久| 美女高潮喷水抽搐中文字幕| 精品高清国产在线一区| 亚洲精品中文字幕一二三四区 | 日韩有码中文字幕| 国精品久久久久久国模美| 欧美大码av| 亚洲七黄色美女视频| 久久青草综合色| 欧美午夜高清在线| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久| e午夜精品久久久久久久| 三上悠亚av全集在线观看| 亚洲色图 男人天堂 中文字幕| 一边摸一边做爽爽视频免费| 十八禁网站免费在线| 国产精品美女特级片免费视频播放器 | 黄频高清免费视频| 欧美精品av麻豆av| 精品人妻1区二区| 两性夫妻黄色片| 九色亚洲精品在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 精品人妻在线不人妻| 一本色道久久久久久精品综合| 日韩精品免费视频一区二区三区| 青青草视频在线视频观看| 五月开心婷婷网| 一本综合久久免费| 久久精品91无色码中文字幕| 欧美日韩亚洲高清精品| 日本撒尿小便嘘嘘汇集6| 美女午夜性视频免费| 免费av中文字幕在线| 中文字幕人妻丝袜一区二区| 欧美日韩国产mv在线观看视频| 免费女性裸体啪啪无遮挡网站| av免费在线观看网站| 大香蕉久久成人网| 亚洲全国av大片| 亚洲成国产人片在线观看| 丁香欧美五月| 制服人妻中文乱码| 大香蕉久久网| 欧美 亚洲 国产 日韩一| 1024视频免费在线观看| 狠狠精品人妻久久久久久综合| 日韩大片免费观看网站| 久久久欧美国产精品| 制服人妻中文乱码| 国产成人精品在线电影| 欧美日韩av久久| 久久精品亚洲av国产电影网| 乱人伦中国视频| 2018国产大陆天天弄谢| 精品少妇黑人巨大在线播放| 成人永久免费在线观看视频 | 999精品在线视频| 国产精品秋霞免费鲁丝片| 美女午夜性视频免费| 久久人人爽av亚洲精品天堂| 一个人免费在线观看的高清视频| 9热在线视频观看99| 日本av免费视频播放| 97在线人人人人妻| 久久久精品94久久精品| 国产精品.久久久| 欧美精品啪啪一区二区三区| 香蕉国产在线看| 久久久久国内视频| 三上悠亚av全集在线观看| 成年女人毛片免费观看观看9 | 国产男女超爽视频在线观看| 国产精品一区二区精品视频观看| 亚洲色图 男人天堂 中文字幕| 免费看a级黄色片| 免费av中文字幕在线| 国产欧美亚洲国产| 在线 av 中文字幕| 久久久国产成人免费| tube8黄色片| 亚洲精品国产色婷婷电影| 日韩熟女老妇一区二区性免费视频| 在线 av 中文字幕| 成人黄色视频免费在线看| 午夜福利视频在线观看免费| 亚洲欧洲精品一区二区精品久久久| 美女高潮喷水抽搐中文字幕| 午夜福利影视在线免费观看| 欧美亚洲 丝袜 人妻 在线| 99国产极品粉嫩在线观看| 亚洲avbb在线观看| 亚洲精品粉嫩美女一区| 欧美国产精品va在线观看不卡| 亚洲人成电影观看| 极品少妇高潮喷水抽搐| 国产aⅴ精品一区二区三区波| 日韩熟女老妇一区二区性免费视频| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕| 日韩成人在线观看一区二区三区| 搡老乐熟女国产| 国产男女超爽视频在线观看| 黑人巨大精品欧美一区二区mp4| 热99re8久久精品国产| 日韩欧美一区视频在线观看| 亚洲精品av麻豆狂野| 大型av网站在线播放| 少妇粗大呻吟视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久av美女十八| 我的亚洲天堂| 最近最新中文字幕大全免费视频| 妹子高潮喷水视频| 亚洲精品乱久久久久久| 中文字幕人妻丝袜一区二区| 热99国产精品久久久久久7| 欧美日韩视频精品一区| 日韩视频在线欧美| 日日摸夜夜添夜夜添小说| 亚洲av日韩在线播放| 12—13女人毛片做爰片一| 欧美成人免费av一区二区三区 | 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 国产精品久久电影中文字幕 | 交换朋友夫妻互换小说| 悠悠久久av| 美女扒开内裤让男人捅视频| 视频区欧美日本亚洲| 久热爱精品视频在线9| 老司机午夜十八禁免费视频| 国产精品免费一区二区三区在线 | 成人特级黄色片久久久久久久 | 日韩免费av在线播放| 多毛熟女@视频| 精品久久久久久久毛片微露脸| 黑人巨大精品欧美一区二区mp4| e午夜精品久久久久久久| 精品乱码久久久久久99久播| 一本大道久久a久久精品| 国产黄频视频在线观看| 一个人免费看片子| 青青草视频在线视频观看| 在线看a的网站| 狠狠精品人妻久久久久久综合| svipshipincom国产片| 久久国产精品大桥未久av| 建设人人有责人人尽责人人享有的| 桃花免费在线播放| 亚洲男人天堂网一区| 757午夜福利合集在线观看| 99热国产这里只有精品6| 欧美黑人精品巨大| 中文字幕高清在线视频| 亚洲人成电影观看| 天天影视国产精品| 黄色毛片三级朝国网站| 这个男人来自地球电影免费观看| 一个人免费看片子| 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| 嫩草影视91久久| 久久影院123| 亚洲av日韩在线播放| 伊人久久大香线蕉亚洲五| www.熟女人妻精品国产| 大陆偷拍与自拍| 日本黄色视频三级网站网址 | 777久久人妻少妇嫩草av网站| 亚洲伊人久久精品综合| 天堂动漫精品| 水蜜桃什么品种好| 欧美黄色淫秽网站| 日本精品一区二区三区蜜桃| 下体分泌物呈黄色| 两人在一起打扑克的视频| 精品少妇黑人巨大在线播放| 国产精品99久久99久久久不卡| av片东京热男人的天堂| 亚洲精品美女久久av网站| 一级a爱视频在线免费观看| 成人国语在线视频| 国产在线免费精品| 黄网站色视频无遮挡免费观看| 老司机福利观看| 视频区欧美日本亚洲| 精品国产亚洲在线| 国产激情久久老熟女| 叶爱在线成人免费视频播放| 欧美 日韩 精品 国产| 伦理电影免费视频| 久久精品国产亚洲av高清一级| 国产深夜福利视频在线观看|