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

    遺傳算法在節(jié)水高效灌溉制度優(yōu)化中的應(yīng)用
    ----以冬小麥為例

    2016-03-22 03:03:28
    節(jié)水灌溉 2016年12期
    關(guān)鍵詞:冬小麥生育期算子

    王 福 增

    (河南省人民勝利渠管理局,河南 新鄉(xiāng) 453003)

    0 引 言

    農(nóng)業(yè)作為水資源的用水大戶,其用水量約占世界水資源總量的70%,占我國水資源總量的80%以上,農(nóng)業(yè)用水中有90%以上都用于灌溉[1]。然而,農(nóng)業(yè)用水的有效利用率僅為30%~40%左右[2]。而且,在社會飛速發(fā)展和人民生活水平不斷提升的今天,水資源的供需矛盾日益突出[3],季節(jié)性水資源短缺問題也越來越嚴(yán)重,農(nóng)業(yè)節(jié)水灌溉越來越受到人們的關(guān)注。非充分灌溉正是在這種灌溉水利用效率不高和水資源短缺的現(xiàn)狀下提出的[3]。

    在非充分灌溉條件下進(jìn)行最優(yōu)灌溉制度的設(shè)計就是在作物各個生育階段對一定的灌溉水量進(jìn)行合理的分配,簡單地說就是明確各個階段的配水量,使得在供水不充足的情況下能收獲的產(chǎn)量最高。目前設(shè)計的冬小麥灌溉制度大多采用的是隨機(jī)動態(tài)規(guī)劃和動態(tài)規(guī)劃的方法[4],在這些方法中大多都把作物的實際騰發(fā)量看作決策變量,而事實上,作物的實際騰發(fā)量是受土壤含水量影響的。兩者在動態(tài)規(guī)劃計算中需要進(jìn)行大量的計算,目前還有“維數(shù)災(zāi)”的問題沒有得到很好的處理。

    所謂的遺傳算法是在生物優(yōu)勝劣汰進(jìn)化機(jī)制的基礎(chǔ)上形成的一種方法,這種方法可以對全局進(jìn)行搜索來得到最優(yōu)解,并且對于難度較大的高度非線性化函數(shù)也可進(jìn)行優(yōu)化。目前,該方法在數(shù)值優(yōu)化、工程設(shè)計以及作業(yè)調(diào)度等方面已經(jīng)得到了應(yīng)用[5]。張兵[6]等在靜態(tài)的作物水分生產(chǎn)函數(shù)Jensen模型的基礎(chǔ)上建立了使產(chǎn)量最大的有限灌水量優(yōu)化灌溉模型,并利用遺傳算法對其進(jìn)行求解,結(jié)果證明在多種不同的情況下該算法都能找到模型的最優(yōu)解,在作物生育期內(nèi)對有限的灌溉水量進(jìn)行合理分配。路振廣[7]等利用遺傳算法對灌區(qū)的水資源優(yōu)化配置模型進(jìn)行求解,分析冬小麥不同水文年的最優(yōu)灌溉制度,發(fā)現(xiàn)依據(jù)優(yōu)化的灌溉制度進(jìn)行灌溉可以在冬小麥的生長期內(nèi)對水資源進(jìn)行優(yōu)化配置,提高水資源的利用效率。上述模型中沒有充分考慮到灌溉對于農(nóng)田土壤水分動態(tài)變化和田間騰發(fā)的影響,這可能會對灌溉制度的確定產(chǎn)生一定影響。

    本文提出采用遺傳算法來求解非充分灌溉條件下冬小麥的最優(yōu)灌溉制度設(shè)計,同時也給出了與之對應(yīng)的約束條件及處理方法,最后通實例計算來說明該算法及處理方法的可行性,對冬小麥整個生育期內(nèi)的灌溉水量進(jìn)行合理的分配,在冬小麥的缺水敏感期內(nèi)進(jìn)行灌溉,節(jié)約灌溉用水,提高水分利用效率。

    1 數(shù)學(xué)模型的構(gòu)建

    1.1 目標(biāo)函數(shù)的確定

    把規(guī)劃期產(chǎn)量最大作為目標(biāo),對冬小麥各個生育期的灌水量進(jìn)行優(yōu)化設(shè)計。使用的作物水分生產(chǎn)函數(shù)是在供水不足的情況下,具有代表性的作物實際產(chǎn)量和水量的連乘模型,目標(biāo)函數(shù)是在單位面積下實際產(chǎn)量和最高產(chǎn)量的比值的最大值[2],可用下式表示:

    (1)

    式中:F表示冬小麥整個生育期的產(chǎn)量;Ya、Ym分別表示單位面積的實際產(chǎn)量和最高產(chǎn)量;ETi、ETmi分別表示第i個生育期的實際騰發(fā)量和最大騰發(fā)量,mm;λi表示第i個生育期作物在缺水條件下的敏感指數(shù)。具體數(shù)據(jù)可以由試驗相關(guān)資料得到。

    1.2 邊界條件的確立

    1.2.1 供水量條件約束

    整個生育期的可灌水量應(yīng)等于各個生育階段的供水量之和,可用下式表示:

    (2)

    式中:Ii表示第i個生育期分配的水量,mm;Q表示作物整個生育期可分配的水量,mm。

    1.2.2 水量平衡方程[9]

    Wi+1=Wi+Pi+Ii-ETi-Qi

    (3)

    式中:Wi、Wi+1分別表示第i個生育期初、時段末根系層的土壤含水量,mm;Pi表示第i個生育期的有效降雨量,mm;Ii表示第i個生育期的灌溉水量(即分配水量),mm;ETi表示第i個生育期內(nèi)的騰發(fā)量,mm;Qi表示第i個生育期內(nèi)的根系層土壤下邊界水分通量,mm。

    1.2.3 上下邊界約束

    Ii≥0

    (4)

    Wi,min≤Wi≤Wi,max

    (5)

    0≤ETi≤ETmi

    (6)

    式中:Wi,max、Wi,min分別表示第i個生育期內(nèi)土壤根系層土壤含水量的上下限,mm。

    1.3 冬小麥騰發(fā)量、有效降雨量以及土壤下界面水分通量的計算

    1.3.1 實際騰發(fā)量的計算[8]

    冬小麥第i個生育期的實際騰發(fā)量與這個時期的實際灌水量是相關(guān)的,通??刹捎靡韵鹿竭M(jìn)行計算:

    ETi=KsiETmiETmi=KciET0

    (7)

    式中:ET0為參考騰發(fā)量,mm,本文選用FAO推薦的彭曼-蒙特斯公式結(jié)合實測的氣象數(shù)據(jù)來逐日計算作物各個生育期的參考騰發(fā)量;Ksi表示土壤水分脅迫系數(shù),是一個與實時土壤含水量相關(guān)的一個無量綱數(shù)值,可由以下經(jīng)驗公式進(jìn)行計算:

    (8)

    式中:Wi表示第i個生育期內(nèi)作物根系層土壤含水量;Wp表示作物根系層達(dá)到凋萎點時的土壤含水量;Wj表示不產(chǎn)生水分脅迫時的臨界含水量;n表示一無量綱的經(jīng)驗指數(shù)。

    Kci為作物系數(shù),是一個無量綱的數(shù),受作物種類、品種以及生育期影響,可用下式進(jìn)行表示:

    (9)

    式中:Kcm表示作物整個生育期內(nèi)作物系數(shù)的最大值;t表示從作物播種到計算點的天數(shù),d;tm表示從作物播種到作物系數(shù)達(dá)到最大時的天數(shù),d;c表示形狀參數(shù),無量綱。

    1.3.2 土壤下界面水分通量的計算

    第i個生育期內(nèi)作物根系層土壤下界面的水分通量可以用以下經(jīng)驗公式進(jìn)行計算:

    (10)

    式中:Wf表示作物根系層田間持水量;Wc表示作物根系層底部進(jìn)行水分交換的臨界含水量;α、d表示無量綱的經(jīng)驗系數(shù)。

    2 求解方法

    2.1 混合遺傳算法過程

    混合遺傳算法[10]是根據(jù)進(jìn)化過程中各代中最差的個體中所包含的需要進(jìn)行優(yōu)化的信息,把下降搜索算子套嵌在遺傳算法中,以此來提高該算法進(jìn)行局部尋找的能力。具體的過程如下所示:

    (1)初始化設(shè)置。進(jìn)化代數(shù)計數(shù)器設(shè)置為t←0;對最大進(jìn)化代數(shù)T進(jìn)行設(shè)置;把隨機(jī)生成的M個個體當(dāng)作初始的群體P(t):t←1。

    (2)對群體中的個體進(jìn)行評估。計算群體中的每個個體的適應(yīng)度,找到最優(yōu)的個體(即種群中適應(yīng)度最大的個體)。

    (3)選擇算子,并對算子進(jìn)行交叉、變異。先選擇算子,并對所選擇的算子進(jìn)行交叉、變異,然后作用于原始群體,可得到一個新的群體P′(t)={x′k(t)|k=0,1,2,…,M},該新的種群中個體總數(shù)為M+1個。

    (4)對新種群中個體進(jìn)行評估。計算新群體P′(t)中各個個體的適應(yīng)度,找到新群體中最差的個體。

    (5)下降算子。 除去最差的個體,沿x′k(t)→x′worst(t)方向?qū)ζ渌麄€體x′k(t)進(jìn)行線性搜索,這樣可以得到第三個群體P*(t)={x*k(t)|k=0,1,2,…,M-1},并且有f(x*k(t))=minf(x*k(t)+λ[x′k(t)-x′worst(t)]),k=0,1,2,…,M-1。將P*(t)看作是下一代新的群體,依次類推進(jìn)行計算。

    (6)終止條件。如果t

    利用這種混合遺傳算法可以使局部搜索的能力得到提升,同時也使全局搜索能力有所增強(qiáng)。

    2.2 約束條件的處理

    在遺傳算法中[11-12],往往只能得到變量的上下限約束問題,而對于非界限的約束,目前常用 的方法是罰函數(shù)法。可是這種方法并不能很好地計算線性約束,得不到很好的計算效果。

    本文中冬小麥的最優(yōu)灌溉制度模型是一個有約束的非線性規(guī)劃模型,以各個生育期的灌溉水量作為決策變量,在確定了灌溉水量后,在計算的基礎(chǔ)上可以確定其他變量。所以該模型中需通過構(gòu)建初始種群變異算子實現(xiàn)對非界限約束式(2)進(jìn)行處理。

    在進(jìn)行處理時,首先應(yīng)該在I1的取值范圍內(nèi)生成I1(0),接著再生成I2(0),這樣類推進(jìn)行,可以生成In-1(0)。對于最后In(0)不是通過隨機(jī)生成得到,而是通過進(jìn)行計算得到,只有這樣才能保證式(2)能夠總是成立。

    得到各個生育期灌水量后可以通過計算得到其他一些變量的值,使每一個個體都有可行解,加速了算法的收斂速度。

    2.3 選擇遺傳算子

    該遺傳算法中有3個基本的算子,分別是選擇算子、交叉算子和變異算子。對于前兩種算子,可以用聯(lián)賽方式來進(jìn)行選擇算子;算術(shù)交叉的方式進(jìn)行交叉算子,在經(jīng)過算術(shù)交叉后得到的新群體仍能滿足變量上下限約束和線性約束。對于變異算子,本文沒有使用均勻變異,而是按照群體的產(chǎn)生方式,改變其中的一個變量會對其后的所有變量產(chǎn)生影響。所以,如果生成的隨機(jī)數(shù)比變異概率大,就可隨機(jī)選擇一個個體看作是變異算子,繼續(xù)采用2.2中所用的方法來生成新個體。

    3 實例驗證

    3.1 試驗區(qū)概況

    試驗于2013-2014年在河南省人民勝利渠灌區(qū)進(jìn)行,試驗區(qū)是典型的平原氣候,多年平均風(fēng)速為2.5 m/s,多年平均降水量570 mm,土壤干容重1.47 g/cm3,孔隙率42.89%,土壤田間持水量為22.32%(占干土重)。地下水埋深大于8 m。冬小麥整個生育期都沒有地下水的補(bǔ)給,即公式(3)中的Qi是0。

    3.2 基本數(shù)據(jù)資料

    根據(jù)彭曼-蒙特斯公式結(jié)合實測的氣象數(shù)據(jù)可得到逐日的參考作物需水量,變化過程如圖1所示。

    圖1 冬小麥全生育期潛在作物需水量變化過程

    依據(jù)冬小麥吸收水分抑制程度并參考當(dāng)?shù)氐臍庀蟛糠侄←湼珊抵笜?biāo)確定冬小麥非充分灌溉的上限為90%的田間持水量。作物根系層達(dá)到凋萎點時的土壤貯水量Wp為97 mm;不產(chǎn)生水分脅迫時的臨界含水量Wj為194 mm;n取1。試驗結(jié)果表明,充分灌溉條件下小麥的最高產(chǎn)量是9 500 kg/hm2。冬小麥的作物系數(shù)、缺水敏感指數(shù)等指標(biāo)見表1。

    表1 計算需要的資料

    3.3 實例計算

    利用遺傳算法,結(jié)合獲得的資料可以分別計算出灌溉定額為90、135、180和225 mm時的水量分配情況。

    由表2可知,本文的遺傳算法得到的灌水日期和最大相對產(chǎn)量與動態(tài)規(guī)劃得到的結(jié)果基本是一致的,大部分的灌水都分布在4月到5月之間,此時正值冬小麥的拔節(jié)期和抽穗期,這之間冬小麥的敏感指數(shù)明顯高于其他各個時期。所以,利用該方法得到的計算結(jié)果與將有限的水量灌溉在缺水敏感期的原則是相符合的。

    4 結(jié) 論

    文中主要對作物最優(yōu)的灌溉制度進(jìn)行設(shè)計,在利用遺傳算法進(jìn)行計算的同時,也提出了一種加快收斂且保證解可行性的偽隨機(jī)方式。最后通過實例證明,遺傳算法和利用偽隨機(jī)方式處理約束條件的策略是可行的,利用該方法對冬小麥各個生育期的灌溉水量進(jìn)行優(yōu)化,可以把有限的水灌溉到作物產(chǎn)量缺水最敏感的時期,節(jié)約灌溉水量,同時也可提高水分利用效率。

    表2 優(yōu)化灌溉制度計算成果表

    [1] 趙 永,蔡煥杰,張朝勇. 非充分灌溉研究現(xiàn)狀及存在問題[J]. 中國農(nóng)村水利水電,2004,(4):1-4.

    [2] 上官周平,邵明安,薛增召. 旱地作物需水量預(yù)報決策輔助系統(tǒng)[J]. 農(nóng)業(yè)工程學(xué)報,2001,(3):42-46.

    [3] 楊 靜,王玉萍,王 群,等. 非充分灌溉的研究進(jìn)展及展望[J]. 安徽農(nóng)業(yè)科學(xué),2008,36(8):3 301-3 303.

    [4] 崔遠(yuǎn)來. 非充分灌溉優(yōu)化配水技術(shù)研究綜述[J]. 灌溉排水學(xué)報,2000,19(2):66-70.

    [5] 宋朝紅,崔遠(yuǎn)來,羅 強(qiáng). 基于遺傳算法的非充分灌溉下最優(yōu)灌溉制度設(shè)計[J]. 灌溉排水學(xué)報,2005,24(6):45-48.

    [6] 張 兵,袁壽其,李 紅,等. 基于遺傳算法求解的冬小麥優(yōu)化灌溉產(chǎn)量模型研究[J]. 農(nóng)業(yè)工程學(xué)報,2006,22(8):12-15.

    [7] 路振廣,孟春紅,王曉麗. 基于遺傳算法的冬小麥優(yōu)化灌溉制度研究[J].人民黃河,2009,31(11):76-78.

    [8] 汪志農(nóng). 灌溉排水工程學(xué)[M]. 2版.北京:中國農(nóng)業(yè)出版社,2013:35,50.

    [9] 張志宇,郄志紅,吳鑫淼. 冬小麥-夏玉米輪作體系灌溉制度多目標(biāo)優(yōu)化模型[J]. 農(nóng)業(yè)工程學(xué)報,2013,29(16):102-111.

    [10] 宋朝紅,羅 強(qiáng),紀(jì)昌明. 基于下降搜索的混合遺傳算法[J]. 武漢大學(xué)學(xué)報(工學(xué)版),2003,36,(1):109-122.

    [11] 遲道才,曲 霞,崔 磊,等. 基于遺傳算法的投影尋蹤回歸模型在參考作物滕發(fā)量預(yù)測中的應(yīng)用[J]. 節(jié)水灌溉,2011,(2):5-7.

    [12] 張 兵,張建生,黃文生. 基于遺傳算法的西北旱區(qū)灌溉效益決策模型研究[J]. 節(jié)水灌溉,2013,(1):61-63.

    猜你喜歡
    冬小麥生育期算子
    擬微分算子在Hp(ω)上的有界性
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    甘肅冬小麥田
    Roper-Suffridge延拓算子與Loewner鏈
    基于作物生育期的潛在蒸散的時空演變特征及R/S 分析
    西藏科技(2015年5期)2015-09-26 11:55:25
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    2013-2014年度二二二團(tuán)冬小麥各生育期氣象條件分析
    阿拉爾市2012年棉花生育期氣象條件分析
    冬小麥——新冬18號
    国产精品久久久久久精品电影| 夜夜爽天天搞| 深夜精品福利| 一二三四在线观看免费中文在| 一本综合久久免费| 日日摸夜夜添夜夜添小说| 久久久精品大字幕| 啦啦啦韩国在线观看视频| 免费观看精品视频网站| 两个人看的免费小视频| 久久伊人香网站| 国产又黄又爽又无遮挡在线| 少妇的丰满在线观看| 精品欧美国产一区二区三| 国产成人精品无人区| 51午夜福利影视在线观看| 老汉色av国产亚洲站长工具| 黄色视频,在线免费观看| 男女下面进入的视频免费午夜| 久久精品国产清高在天天线| 九九热线精品视视频播放| av视频在线观看入口| 两性夫妻黄色片| 非洲黑人性xxxx精品又粗又长| 午夜福利成人在线免费观看| 欧美一级a爱片免费观看看| 国产亚洲精品av在线| 一二三四在线观看免费中文在| 国语自产精品视频在线第100页| 全区人妻精品视频| 巨乳人妻的诱惑在线观看| 国产黄a三级三级三级人| 黄色视频,在线免费观看| 99久久无色码亚洲精品果冻| 亚洲成av人片免费观看| av片东京热男人的天堂| 日韩欧美国产在线观看| 91九色精品人成在线观看| 美女扒开内裤让男人捅视频| 亚洲av成人精品一区久久| 日韩欧美在线乱码| 精品国产亚洲在线| 亚洲精品一区av在线观看| 丰满的人妻完整版| 国产精品女同一区二区软件 | 变态另类丝袜制服| 亚洲18禁久久av| 午夜免费观看网址| 国产精品久久久久久精品电影| 午夜免费观看网址| 视频区欧美日本亚洲| 免费电影在线观看免费观看| 舔av片在线| 国产精品日韩av在线免费观看| 日日干狠狠操夜夜爽| 青草久久国产| 国产成人av教育| 美女高潮的动态| 欧美极品一区二区三区四区| 最近最新中文字幕大全电影3| 国产av麻豆久久久久久久| 成人国产综合亚洲| 久久精品人妻少妇| 色在线成人网| 宅男免费午夜| 俄罗斯特黄特色一大片| 亚洲黑人精品在线| 一进一出抽搐动态| 国产高清videossex| 欧美日本视频| 国产亚洲精品一区二区www| 真人做人爱边吃奶动态| 小蜜桃在线观看免费完整版高清| 国产成人精品久久二区二区免费| 午夜日韩欧美国产| 国产毛片a区久久久久| 欧美不卡视频在线免费观看| 国内精品一区二区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品午夜福利视频在线观看一区| 国产美女午夜福利| 国内精品久久久久久久电影| 久久草成人影院| 欧美一级a爱片免费观看看| 怎么达到女性高潮| 天堂av国产一区二区熟女人妻| 亚洲一区二区三区不卡视频| 国产精品亚洲美女久久久| 亚洲国产日韩欧美精品在线观看 | 黄色片一级片一级黄色片| a在线观看视频网站| 欧美日韩乱码在线| 精品久久久久久久末码| 欧美激情在线99| 亚洲av第一区精品v没综合| 欧美激情在线99| 国产精品久久视频播放| 一边摸一边抽搐一进一小说| 久久精品91无色码中文字幕| 午夜免费激情av| 99riav亚洲国产免费| 国产欧美日韩一区二区精品| 午夜精品久久久久久毛片777| 国产淫片久久久久久久久 | 欧美3d第一页| 久久草成人影院| 级片在线观看| 国产aⅴ精品一区二区三区波| 亚洲欧洲精品一区二区精品久久久| 在线播放国产精品三级| 人妻丰满熟妇av一区二区三区| 男女那种视频在线观看| 国产久久久一区二区三区| 日本黄色片子视频| 成熟少妇高潮喷水视频| 成人18禁在线播放| 啦啦啦观看免费观看视频高清| 久久久久久久久中文| 黑人操中国人逼视频| 老司机深夜福利视频在线观看| 国产精品亚洲美女久久久| 老熟妇乱子伦视频在线观看| 老汉色av国产亚洲站长工具| 久久中文看片网| 琪琪午夜伦伦电影理论片6080| 色综合站精品国产| 最近视频中文字幕2019在线8| 亚洲精品国产精品久久久不卡| 免费大片18禁| 日日干狠狠操夜夜爽| 夜夜爽天天搞| a级毛片在线看网站| 国产精华一区二区三区| 久久久国产成人免费| 首页视频小说图片口味搜索| 日韩欧美在线二视频| 免费观看精品视频网站| 一级毛片精品| 国产免费男女视频| av福利片在线观看| 国内精品美女久久久久久| 好男人电影高清在线观看| 日韩免费av在线播放| 女人高潮潮喷娇喘18禁视频| 99热只有精品国产| 一进一出好大好爽视频| 亚洲激情在线av| 国产欧美日韩精品亚洲av| 99久久精品热视频| 哪里可以看免费的av片| 国产一区二区在线av高清观看| 女生性感内裤真人,穿戴方法视频| 国产三级黄色录像| 少妇裸体淫交视频免费看高清| 亚洲精品中文字幕一二三四区| 丝袜人妻中文字幕| 成人性生交大片免费视频hd| 香蕉国产在线看| 午夜激情欧美在线| 伦理电影免费视频| 欧美乱色亚洲激情| 婷婷精品国产亚洲av| 三级国产精品欧美在线观看 | 偷拍熟女少妇极品色| 国产成人精品无人区| 一边摸一边抽搐一进一小说| 久久久色成人| 欧美中文日本在线观看视频| 亚洲精品456在线播放app | 成年免费大片在线观看| 成年免费大片在线观看| 国产在线精品亚洲第一网站| 色综合婷婷激情| 亚洲欧美日韩东京热| 国语自产精品视频在线第100页| 18禁观看日本| 少妇丰满av| 亚洲精品国产精品久久久不卡| 少妇裸体淫交视频免费看高清| 国模一区二区三区四区视频 | 国产毛片a区久久久久| 黄色日韩在线| 伊人久久大香线蕉亚洲五| 午夜a级毛片| 久久久久久大精品| 精品久久久久久久末码| 一级作爱视频免费观看| 成年女人毛片免费观看观看9| 国产综合懂色| 久久久久久久久中文| 欧美不卡视频在线免费观看| 久久久久久人人人人人| av中文乱码字幕在线| 久久精品国产清高在天天线| 亚洲美女视频黄频| 欧美日韩中文字幕国产精品一区二区三区| or卡值多少钱| 男插女下体视频免费在线播放| 欧美日韩亚洲国产一区二区在线观看| 国产免费男女视频| 亚洲国产精品sss在线观看| 19禁男女啪啪无遮挡网站| 国产精品精品国产色婷婷| 国产成人影院久久av| 老司机午夜福利在线观看视频| 制服人妻中文乱码| 亚洲va日本ⅴa欧美va伊人久久| 国产高潮美女av| 一a级毛片在线观看| 国产三级黄色录像| 男插女下体视频免费在线播放| av片东京热男人的天堂| 午夜福利成人在线免费观看| 黑人欧美特级aaaaaa片| 99久久99久久久精品蜜桃| 精品国产超薄肉色丝袜足j| 亚洲成人中文字幕在线播放| 日韩欧美三级三区| 听说在线观看完整版免费高清| 国模一区二区三区四区视频 | 精品不卡国产一区二区三区| 亚洲精品在线观看二区| 久久天堂一区二区三区四区| 亚洲欧美日韩卡通动漫| 午夜免费激情av| 国产亚洲欧美98| 国产成人av激情在线播放| 麻豆成人午夜福利视频| 日本黄色片子视频| 天堂影院成人在线观看| 国产亚洲精品av在线| 在线观看免费视频日本深夜| 丁香欧美五月| 亚洲九九香蕉| av女优亚洲男人天堂 | 国产爱豆传媒在线观看| 精品国产美女av久久久久小说| av在线天堂中文字幕| 可以在线观看的亚洲视频| 午夜福利成人在线免费观看| 一级毛片精品| 国产精品一区二区精品视频观看| 1000部很黄的大片| 欧美在线黄色| 舔av片在线| 黄片大片在线免费观看| 久久精品国产99精品国产亚洲性色| 91老司机精品| 亚洲欧美日韩东京热| 久久精品人妻少妇| 国产高清激情床上av| 国产三级黄色录像| 天天躁日日操中文字幕| 99久久精品国产亚洲精品| 亚洲欧美日韩高清在线视频| 国产精品99久久久久久久久| 久久精品91蜜桃| 草草在线视频免费看| 欧美日韩乱码在线| 国产成人av教育| 精品久久久久久久久久久久久| 欧美激情久久久久久爽电影| 神马国产精品三级电影在线观看| 国产精品久久久久久精品电影| 日本a在线网址| 非洲黑人性xxxx精品又粗又长| 久久热在线av| 又紧又爽又黄一区二区| 亚洲片人在线观看| 久久久精品欧美日韩精品| tocl精华| 亚洲avbb在线观看| 韩国av一区二区三区四区| 99视频精品全部免费 在线 | 91久久精品国产一区二区成人 | 精品乱码久久久久久99久播| 亚洲黑人精品在线| 国产成人欧美在线观看| 国产激情欧美一区二区| 成人特级黄色片久久久久久久| 一区福利在线观看| а√天堂www在线а√下载| aaaaa片日本免费| 国产探花在线观看一区二区| 亚洲专区字幕在线| 国产精品久久久久久亚洲av鲁大| 亚洲精品色激情综合| 亚洲国产欧美一区二区综合| 免费看日本二区| 日韩国内少妇激情av| 国产高潮美女av| 变态另类丝袜制服| 毛片女人毛片| 成人18禁在线播放| 色老头精品视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 综合色av麻豆| 午夜精品久久久久久毛片777| 欧美一级a爱片免费观看看| 一个人看视频在线观看www免费 | 欧美日韩瑟瑟在线播放| 神马国产精品三级电影在线观看| 深夜精品福利| 两性夫妻黄色片| 两个人看的免费小视频| 岛国在线观看网站| avwww免费| 欧美日韩乱码在线| 女人被狂操c到高潮| 亚洲成av人片在线播放无| 此物有八面人人有两片| 国产黄a三级三级三级人| 久久久久国内视频| 国产在线精品亚洲第一网站| av天堂中文字幕网| 美女午夜性视频免费| 一个人免费在线观看电影 | 黑人巨大精品欧美一区二区mp4| 18禁黄网站禁片免费观看直播| 免费大片18禁| 可以在线观看的亚洲视频| 国产精品野战在线观看| 三级国产精品欧美在线观看 | 夜夜夜夜夜久久久久| 成人欧美大片| 999久久久精品免费观看国产| avwww免费| 免费人成视频x8x8入口观看| 色尼玛亚洲综合影院| 小说图片视频综合网站| 日本与韩国留学比较| 俄罗斯特黄特色一大片| e午夜精品久久久久久久| 脱女人内裤的视频| 在线观看免费视频日本深夜| 97人妻精品一区二区三区麻豆| 香蕉丝袜av| 性欧美人与动物交配| 又黄又粗又硬又大视频| 午夜福利18| or卡值多少钱| 在线看三级毛片| 一本综合久久免费| 日韩人妻高清精品专区| 国产精品免费一区二区三区在线| 丝袜人妻中文字幕| 婷婷丁香在线五月| 日本熟妇午夜| 国产免费av片在线观看野外av| 色精品久久人妻99蜜桃| 国内久久婷婷六月综合欲色啪| 757午夜福利合集在线观看| 精华霜和精华液先用哪个| 男女之事视频高清在线观看| 蜜桃久久精品国产亚洲av| 黄色成人免费大全| 欧美黑人欧美精品刺激| 欧美色视频一区免费| 午夜日韩欧美国产| 在线视频色国产色| 夜夜爽天天搞| 久久午夜综合久久蜜桃| 夜夜看夜夜爽夜夜摸| 亚洲中文av在线| 性欧美人与动物交配| 手机成人av网站| 美女黄网站色视频| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 国产三级在线视频| 亚洲aⅴ乱码一区二区在线播放| 精品国产乱子伦一区二区三区| 日韩人妻高清精品专区| 国产精品自产拍在线观看55亚洲| 成人亚洲精品av一区二区| 国产精品98久久久久久宅男小说| 亚洲18禁久久av| 欧美乱妇无乱码| 亚洲中文日韩欧美视频| 亚洲 国产 在线| av在线天堂中文字幕| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 国产亚洲欧美在线一区二区| 精品日产1卡2卡| 99视频精品全部免费 在线 | 一个人免费在线观看电影 | 日韩精品中文字幕看吧| 国产一区二区激情短视频| 成人一区二区视频在线观看| 床上黄色一级片| 欧美最黄视频在线播放免费| 亚洲无线在线观看| 看黄色毛片网站| 成人特级av手机在线观看| 99国产精品99久久久久| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| 国产不卡一卡二| 91麻豆av在线| 国产一区在线观看成人免费| 99热6这里只有精品| 一区二区三区激情视频| 麻豆成人午夜福利视频| 国产精品影院久久| 久久草成人影院| 国产真人三级小视频在线观看| 午夜福利成人在线免费观看| 国产免费av片在线观看野外av| 国产高清videossex| 国产精品,欧美在线| 91老司机精品| 久久久久久久午夜电影| 成人永久免费在线观看视频| cao死你这个sao货| 亚洲人成网站在线播放欧美日韩| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国内精品久久久久久久电影| 18禁黄网站禁片免费观看直播| 午夜福利成人在线免费观看| 国产精品久久久久久精品电影| 国产aⅴ精品一区二区三区波| 亚洲av免费在线观看| 啦啦啦观看免费观看视频高清| 麻豆国产97在线/欧美| 久久久久久久久中文| 俄罗斯特黄特色一大片| 手机成人av网站| www.www免费av| 亚洲一区二区三区色噜噜| 亚洲人成电影免费在线| 亚洲乱码一区二区免费版| 91麻豆av在线| 欧美午夜高清在线| 99精品在免费线老司机午夜| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 欧美黑人巨大hd| 中文亚洲av片在线观看爽| 亚洲色图 男人天堂 中文字幕| 日本免费一区二区三区高清不卡| 国产69精品久久久久777片 | 两个人看的免费小视频| 天堂√8在线中文| 日韩欧美在线乱码| 成人午夜高清在线视频| 全区人妻精品视频| bbb黄色大片| 99国产极品粉嫩在线观看| 久久草成人影院| 久久精品国产亚洲av香蕉五月| 两个人看的免费小视频| 90打野战视频偷拍视频| 又黄又爽又免费观看的视频| 嫁个100分男人电影在线观看| 久久欧美精品欧美久久欧美| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 久久草成人影院| 久久精品夜夜夜夜夜久久蜜豆| 久久久久精品国产欧美久久久| 91麻豆精品激情在线观看国产| 免费观看精品视频网站| 国产精品av久久久久免费| 91九色精品人成在线观看| 欧美激情在线99| 亚洲中文字幕日韩| 在线免费观看不下载黄p国产 | 国产久久久一区二区三区| 欧美黄色片欧美黄色片| 一级毛片高清免费大全| 欧美成人性av电影在线观看| 一区二区三区激情视频| 99热6这里只有精品| 巨乳人妻的诱惑在线观看| 中文字幕最新亚洲高清| 国产亚洲欧美在线一区二区| 欧美精品啪啪一区二区三区| а√天堂www在线а√下载| 亚洲欧美精品综合一区二区三区| 久久久国产成人免费| 国产野战对白在线观看| 伦理电影免费视频| 琪琪午夜伦伦电影理论片6080| 一个人看视频在线观看www免费 | 国产在线精品亚洲第一网站| 久久久国产精品麻豆| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 1000部很黄的大片| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 中文资源天堂在线| 男女视频在线观看网站免费| 一进一出抽搐gif免费好疼| 欧美中文日本在线观看视频| 国产精品,欧美在线| 久久久久久久久中文| cao死你这个sao货| 不卡一级毛片| 桃色一区二区三区在线观看| 丰满的人妻完整版| 国产成人精品久久二区二区免费| 国产高清有码在线观看视频| 亚洲中文字幕一区二区三区有码在线看 | xxx96com| 精品久久久久久久毛片微露脸| www日本在线高清视频| 欧美xxxx黑人xx丫x性爽| 精品电影一区二区在线| 熟女少妇亚洲综合色aaa.| 亚洲精品美女久久久久99蜜臀| 一本综合久久免费| 午夜免费激情av| 午夜福利在线在线| 高清在线国产一区| 国产亚洲欧美在线一区二区| 白带黄色成豆腐渣| 亚洲人成电影免费在线| 一级a爱片免费观看的视频| 久久国产精品人妻蜜桃| 国产精品av久久久久免费| 在线观看免费视频日本深夜| 免费在线观看成人毛片| 国产欧美日韩一区二区三| 人人妻人人看人人澡| 天天躁日日操中文字幕| 好男人电影高清在线观看| 国产激情久久老熟女| 亚洲精品粉嫩美女一区| 欧美不卡视频在线免费观看| 在线观看日韩欧美| 亚洲片人在线观看| 波多野结衣巨乳人妻| 男女视频在线观看网站免费| 91字幕亚洲| 国产精品 欧美亚洲| 琪琪午夜伦伦电影理论片6080| 老熟妇乱子伦视频在线观看| 女警被强在线播放| 久久婷婷人人爽人人干人人爱| 精品久久久久久久人妻蜜臀av| 免费搜索国产男女视频| 伊人久久大香线蕉亚洲五| 欧美中文综合在线视频| 欧美一区二区国产精品久久精品| 白带黄色成豆腐渣| 国产成人av教育| 久久婷婷人人爽人人干人人爱| 中文字幕高清在线视频| 亚洲精品久久国产高清桃花| 亚洲精品一区av在线观看| 最近最新中文字幕大全电影3| 偷拍熟女少妇极品色| 亚洲 欧美 日韩 在线 免费| 久久久国产成人精品二区| 99精品欧美一区二区三区四区| 欧美绝顶高潮抽搐喷水| 1024手机看黄色片| 最新中文字幕久久久久 | 久久中文看片网| 一本综合久久免费| 后天国语完整版免费观看| 99re在线观看精品视频| 丁香六月欧美| 午夜福利在线观看吧| 女人高潮潮喷娇喘18禁视频| 国产高清三级在线| 亚洲七黄色美女视频| 一区二区三区高清视频在线| 亚洲国产中文字幕在线视频| 九九热线精品视视频播放| 日本成人三级电影网站| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久| 岛国在线观看网站| 国产成人av教育| 午夜免费观看网址| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美日韩高清专用| 一区福利在线观看| 9191精品国产免费久久| 国产亚洲av嫩草精品影院| 巨乳人妻的诱惑在线观看| 国产精品野战在线观看| 欧美成人一区二区免费高清观看 | 欧美国产日韩亚洲一区| 综合色av麻豆| 日本黄色视频三级网站网址| 婷婷丁香在线五月| 视频区欧美日本亚洲| 亚洲熟女毛片儿| 亚洲欧美精品综合久久99| 国产三级黄色录像| 久久香蕉精品热| 欧美一级毛片孕妇| 亚洲av日韩精品久久久久久密| 99视频精品全部免费 在线 | 亚洲精品色激情综合| 国产毛片a区久久久久| 99热6这里只有精品| 亚洲av美国av| 亚洲av五月六月丁香网| 99久久精品国产亚洲精品| 亚洲国产欧洲综合997久久,| 精品福利观看| 日本黄色片子视频| 床上黄色一级片| 黑人巨大精品欧美一区二区mp4| 女人被狂操c到高潮| 国产三级黄色录像| 日本成人三级电影网站| 久久草成人影院|