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

    基于APSO-ELM和模糊邏輯的灌溉時(shí)間決策

    2020-06-12 02:33:20田強(qiáng)明溫宗周李麗敏張順鋒
    中國農(nóng)村水利水電 2020年4期
    關(guān)鍵詞:土壤濕度灌溉作物

    田強(qiáng)明,溫宗周,李麗敏,張順鋒

    (西安工程大學(xué) 電子信息學(xué)院,西安 710048)

    目前對(duì)灌溉控制方式的研究主要有3種。第一種是基于土壤水分和作物水分的方法,該方法一方面需要大量的傳感器,成本過高,另一方面需要前期大量的試驗(yàn)數(shù)據(jù)作灌溉決策支撐;第二種是基于模型灌溉決策的方法,GU Zhe等人利用農(nóng)業(yè)系統(tǒng)模型,對(duì)作物水分脅迫進(jìn)行預(yù)測(cè)和評(píng)價(jià),從而確定灌溉的時(shí)間。經(jīng)過實(shí)驗(yàn)驗(yàn)證,該模型在保證作物產(chǎn)量的基礎(chǔ)上,節(jié)水率達(dá)到了30%,但要通過大量的田間試驗(yàn)數(shù)據(jù)來對(duì)模型進(jìn)行不斷的校準(zhǔn)[1-3];第三種是基于實(shí)際作物蒸發(fā)蒸騰量和水量平衡的方法,這也是目前較為普遍并且容易執(zhí)行的方法,BARTLETT A C等人利用土壤水分平衡方法和農(nóng)業(yè)氣象網(wǎng)絡(luò)數(shù)據(jù)的實(shí)時(shí)查詢,創(chuàng)建了一種基于蒸發(fā)蒸騰量的在線灌溉調(diào)度工具,使得作物可以根據(jù)蒸發(fā)蒸騰量及時(shí)進(jìn)行灌溉[4],這種方式只考慮了蒸發(fā)蒸騰量對(duì)灌溉的影響,忽略了土壤和作物本身因素的影響。PEREA R G等人根據(jù)土壤數(shù)據(jù)和與液壓系統(tǒng)有關(guān)的信息設(shè)計(jì)了一套應(yīng)用灌溉系統(tǒng),向工作人員提供有關(guān)所需灌溉時(shí)間的每日信息,該應(yīng)用已在幾個(gè)商業(yè)農(nóng)場(chǎng)得到應(yīng)用,節(jié)水幅度從11%到33%不等,驗(yàn)證了該系統(tǒng)在實(shí)際應(yīng)用中是可行的[5],但這種方式忽略了環(huán)境因子對(duì)作物實(shí)際需水量的影響。

    在當(dāng)前農(nóng)業(yè)精準(zhǔn)灌溉的需求下,參考作物蒸發(fā)蒸騰量是必須考慮的問題之一。作物蒸發(fā)蒸騰量的計(jì)算大致分為3種。第一種是物理測(cè)定法,包括蒸發(fā)器觀測(cè)、梯度測(cè)量、質(zhì)量平衡等方法,但這些方法操作繁瑣,多用于科學(xué)研究;第二種是基于經(jīng)驗(yàn)?zāi)P头?。該方法主要是利用環(huán)境因子通過數(shù)理統(tǒng)計(jì)的方法建立起來的經(jīng)驗(yàn)?zāi)P?,王娟等人以氣象?shù)據(jù)為基礎(chǔ)對(duì)修正后的Rana和Katerji模型對(duì)華北地區(qū)的作物蒸發(fā)蒸騰量進(jìn)行預(yù)測(cè)[6];閆世程等人根據(jù)土壤水分?jǐn)?shù)據(jù)對(duì)SIMDualKc模型參數(shù)進(jìn)行修正,并證明該模型能夠較好的預(yù)測(cè)西北半干旱地區(qū)滴灌玉米的作物蒸發(fā)蒸騰量[7];田瑋瑋等人通過情景分析法,利用CERES-Wheat模型分析了晉中地區(qū)小麥需水量、小麥蒸發(fā)蒸騰量與產(chǎn)量的關(guān)系,并利用該模型準(zhǔn)確預(yù)測(cè)了小麥的作物蒸發(fā)蒸騰量[8]。這些模型具有較強(qiáng)的可操作性,但針對(duì)性較強(qiáng),具有較差的可移植性,并不適合大范圍推廣;第三種是數(shù)值模擬法,該方法具有效率高,精度高,并且具有較強(qiáng)的泛化性能,是當(dāng)前參考作物蒸發(fā)蒸騰量預(yù)測(cè)研究的熱點(diǎn)問題。

    因標(biāo)準(zhǔn)參考作物蒸發(fā)蒸騰量計(jì)算模型參數(shù)過多,容易造成參數(shù)獲取成本過高,并且有部分參數(shù)不易獲取,所以李家明等人利用通徑分析法研究出漳河灌區(qū)參考作物蒸發(fā)蒸騰量的主要環(huán)境因子[9]。張志政等人在BP神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上引入了PSO,通過組合不同的參考作物蒸發(fā)蒸騰量的影響因子確定了輸入量,最終對(duì)參考作物蒸發(fā)蒸騰量的預(yù)測(cè)精度達(dá)到了0.96[10],但因BP的正向傳播輸出層得不到預(yù)期的效果時(shí),則反向傳播需要不斷地對(duì)網(wǎng)絡(luò)各連接權(quán)值進(jìn)行修正,則會(huì)導(dǎo)致網(wǎng)絡(luò)訓(xùn)練的速度下降。因ELM模型具有較強(qiáng)的泛化性能,所以崔寧博等人利用該模型對(duì)西北地區(qū)的參考作物蒸發(fā)蒸騰量進(jìn)行了預(yù)測(cè),預(yù)測(cè)結(jié)果達(dá)到了0.96[11,12],相較于BP,精度達(dá)到了良好的效果,但卻忽略了訓(xùn)練速度的影響。

    模糊控制在灌溉時(shí)間控制中的研究十分廣泛。吳興利等人利用土壤濕度偏差和濕度偏差變化率作為模糊控制的輸入,以電磁閥的控制量作為輸出,間接實(shí)現(xiàn)對(duì)灌溉時(shí)間的控制[13]。 X PENG等人利用土壤濕度和空氣溫度作為模糊控制的輸入[14],L XIAO等人利用最佳土壤濕度和當(dāng)前土壤濕度的差作為輸入[15],S FAZACKERLEY等人利用土壤濕度、水流失率的數(shù)學(xué)函數(shù)來計(jì)算灌溉時(shí)間[16],劉東等人利用土壤溫度和土壤濕度下降率作為模糊控制的輸入[17]。上述控制的輸出都是灌溉時(shí)間,但考慮的全是土壤因素對(duì)灌溉時(shí)間的影響,忽略了作物本身因素帶來的影響。

    因此,本文提出了一種基于APSO-ELM和模糊邏輯的灌溉時(shí)間決策的控制方法。采用基于實(shí)際作物蒸發(fā)蒸騰量和水量平衡的控制方法,先用APSO-ELM網(wǎng)絡(luò)預(yù)測(cè)參考作物蒸發(fā)蒸騰量,再結(jié)合作物各生長期的作物系數(shù),計(jì)算出實(shí)際作物蒸發(fā)蒸騰量。利用實(shí)際作物蒸發(fā)蒸騰量和土壤濕度下降率通過模糊系統(tǒng),最終輸出作物需要的灌溉時(shí)間。

    1 APSO-ELM算法

    APSO是在PSO的基礎(chǔ)上改進(jìn)而來的,粒子群優(yōu)化算法(PSO)是一種群體智能算法[18]?,F(xiàn)假設(shè)粒子的初始狀態(tài)設(shè)置屬性如下式所示。

    (1)

    (2)

    標(biāo)準(zhǔn)PSO速度和位置更新公式如下所示。

    (3)

    (4)

    由于標(biāo)準(zhǔn)PSO具有全局搜索能力差,尋優(yōu)速度慢等問題,所以采用如下改進(jìn)參數(shù)因子。w大小可以影響算法的全局搜索能力,所以選用非線性動(dòng)態(tài)慣性權(quán)重策略,公式為:

    (5)

    式中:w(t)表示粒子在t時(shí)刻的慣性權(quán)重;k為控制因子,用于調(diào)節(jié)w和t變化曲線的平滑度,經(jīng)過多次試驗(yàn),其k值范圍可取[3,4];wstart為慣性權(quán)重的初始值;wend為慣性權(quán)重終值;tmax為最大迭代次數(shù);t為當(dāng)前迭代次數(shù),本次取值范圍[0.9,0.4]。

    為了避免把局部最優(yōu)解誤作為全局最優(yōu)解以及尋優(yōu)速度過慢的問題,本系統(tǒng)采用學(xué)習(xí)因子自調(diào)整策略,其調(diào)整公式如下。

    (6)

    極限學(xué)習(xí)機(jī)(ELM)是基于前饋神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)算法[19],其主要特點(diǎn)是隱含層節(jié)點(diǎn)參數(shù)可以是隨機(jī)或人為給定的且不需要調(diào)整,學(xué)習(xí)過程僅需計(jì)算輸出權(quán)重。ELM具有學(xué)習(xí)效率高和泛化能力強(qiáng)的優(yōu)點(diǎn),其結(jié)構(gòu)圖如圖1所示。其中,xj表示神經(jīng)元的輸入,N表示輸入層的個(gè)數(shù),L表示隱含層的神經(jīng)元個(gè)數(shù),bi表示每個(gè)神經(jīng)元的偏置,βL表示輸出權(quán)重,oj表示網(wǎng)絡(luò)的輸出。

    圖1 ELM算法結(jié)構(gòu)圖Fig.1 ELM algorithm structure diagram

    2 灌溉時(shí)間決策系統(tǒng)設(shè)計(jì)

    土壤濕度下降率是直接反映作物水分虧缺程度的關(guān)鍵因素,參考作物蒸發(fā)蒸騰量是反映維持作物基本生理活動(dòng)的有效參數(shù)。因此,灌溉量的控制有必要將兩者相結(jié)合,進(jìn)而決定對(duì)灌溉時(shí)間的控制。

    2.1 標(biāo)準(zhǔn)參考作物蒸發(fā)蒸騰量計(jì)算模型

    因?yàn)樵跍厥抑锌梢院雎燥L(fēng)速的影響,所以不適宜采用聯(lián)合國糧農(nóng)組織(FAO)推薦使用的Penman-Monteith (PM)公式計(jì)算參考作物蒸發(fā)蒸騰量的模型,而采用PM在溫室中的修正公式[20]如下:

    (7)

    式中:ET0-PMF56指修正后的參考作物蒸發(fā)蒸騰量,mm/d;T為平均氣溫,℃;es為地表2 m高度處平均飽和蒸氣壓,kPa;ea為地表2 m高度處實(shí)際飽和蒸氣壓,kPa;Rn為凈輻射,MJ/(m2·d);γ為濕度表常數(shù),kPa/℃;G為土壤熱通量,MJ/(m2·d);Δ為溫度~飽和水汽壓關(guān)系曲線在T處的切線斜率。

    通過文獻(xiàn)[21]方法計(jì)算出Δ的值,如下式。

    (8)

    式中:U為溫室內(nèi)空氣的相對(duì)濕度,%。

    實(shí)際作物蒸發(fā)蒸騰量ET0如式(9)所示。

    ET0=ksET0-PMF56

    (9)

    式中:ks為作物生長系數(shù)。

    2.2 APSO-ELM對(duì)參考作物蒸發(fā)蒸騰量預(yù)測(cè)流程

    為了避免由于ELM輸入層與隱含層隨機(jī)生成的權(quán)值和閾值給網(wǎng)絡(luò)預(yù)測(cè)精度帶來不穩(wěn)定性,引入APSO使之可以避免陷入局部最優(yōu),優(yōu)化網(wǎng)絡(luò)帶來的不穩(wěn)定性。整個(gè)樣本分為測(cè)試集和訓(xùn)練集,首先需要初始化ELM和粒子群的基本參數(shù),然后通過APSO不斷的迭代更新,得到ELM網(wǎng)絡(luò)的最優(yōu)權(quán)值和閾值,最終得到網(wǎng)絡(luò)的輸出。其算法流程如圖2所示。

    圖2 APSO-ELM算法流程Fig.2 APSO-ELM algorithm process

    2.3 灌溉時(shí)間模糊邏輯系統(tǒng)流程

    設(shè)置好數(shù)據(jù)采集的時(shí)間間隔t以后,根據(jù)專家經(jīng)驗(yàn)和實(shí)際需求設(shè)置模糊系統(tǒng)規(guī)則庫,利用獲取到的實(shí)際作物蒸發(fā)蒸騰量ET0和土壤濕度下降率ΔE,通過模糊邏輯系統(tǒng),得到灌溉時(shí)間T。其系統(tǒng)流程設(shè)計(jì)如圖3所示。

    圖3 灌溉時(shí)間模糊邏輯系統(tǒng)流程Fig.3 Fuzzy logic system flow of irrigation time

    2.4 隸屬度函數(shù)確定

    系統(tǒng)的輸入量確定為兩個(gè),第一個(gè)輸入量是土壤濕度的下降率ΔE,下降率可以反映變化的快慢;第二個(gè)輸入量是實(shí)際作物蒸發(fā)蒸騰量ET0,ET0可以反映出維持作物生理活動(dòng)最基本需水量。輸出量為灌溉時(shí)間T。其中ΔE語言值設(shè)置為4個(gè),ET0和T語言值設(shè)置為7個(gè)。土壤水分下降率變化偏差ΔE、實(shí)際作物蒸發(fā)蒸騰量ET0、灌溉時(shí)間T的隸屬度函數(shù)都選用三角形隸屬度。輸入輸出隸屬度函數(shù)定義如圖4所示。

    圖4 輸入輸出隸屬度函數(shù)Fig.4 Input/output membership function

    其中輸入變量ΔE的語言變量為:NB、NM、NS、ZO,分別表示土壤水分下降率快速減少、中速減少、低速減少、不變,論域?yàn)閇-3,0];ET0的語言變量為:NB、NM、NS、ZO、PS、PM、PB,分別為零、非常少、少、偏中、中、偏多、多,論域?yàn)閇-3,3];T的語言變量為:ZO、PS、PS+、PM、PM+、PB、PB+分別表示不灌,灌溉時(shí)間很短、灌溉時(shí)間較短、灌溉時(shí)間中等、灌溉時(shí)間較長、灌溉時(shí)間很長、灌溉時(shí)間最長,論域?yàn)閇0,6]。

    2.5 模糊規(guī)則確定

    模糊規(guī)則設(shè)計(jì)的原則要求在土壤濕度下降率以及實(shí)際作物蒸發(fā)蒸騰量變化時(shí),輸出的灌溉時(shí)間要能根據(jù)其變化滿足作物的水分需求。本文灌溉控制模型輸入ΔE有4個(gè)語言變量,ET0語言變量有7個(gè),因此灌溉控制規(guī)則總共有28條。模糊控制器控制規(guī)則如表1所示。例如,當(dāng)土壤濕度下降率ΔE是NB時(shí),即快速減少,實(shí)際蒸發(fā)蒸騰量ET0是PB時(shí),即ET0變化多時(shí),則輸出的灌溉時(shí)間為PB+,即這時(shí)的灌溉時(shí)間最長。經(jīng)過解模糊以后得到實(shí)際的灌溉時(shí)間。

    表1 模糊控制規(guī)則表Tab.1 Fuzzy control rule table

    3 系統(tǒng)性能驗(yàn)證及分析

    本文以番茄種植為例。番茄生長周期經(jīng)歷育苗期、定植期、開花期和結(jié)果期。通過陜西某地2017-2018的番茄生長氣象環(huán)境資料和各個(gè)時(shí)期的作物系數(shù),建立參考作物蒸發(fā)蒸騰量預(yù)測(cè)模型,進(jìn)而得到實(shí)際作物蒸發(fā)蒸騰量。作物系數(shù)采用聯(lián)合國糧農(nóng)組織推薦使用的數(shù)據(jù)得到番茄的各個(gè)時(shí)期的作物生長系數(shù)ks,如表2所示。

    表2 番茄生長系數(shù)Tab.2 Tomato growth coefficient

    3.1 參考作物蒸發(fā)蒸騰量仿真結(jié)果分析

    本文選取參考作物蒸發(fā)蒸騰量中較易獲取的日平均氣溫、日照長度、平均相對(duì)濕度、大氣壓4個(gè)環(huán)境參數(shù)作為模型的輸入。實(shí)驗(yàn)數(shù)據(jù)選取陜西某地2017-2018年的氣象觀測(cè)資料。選取120組數(shù)據(jù),其中103組數(shù)據(jù)用于模型訓(xùn)練,17組數(shù)據(jù)用于模型測(cè)試。圖5顯示了模型仿真均方誤差MSE的變化趨勢(shì),從圖5可以看出模型迭代到45次左右的時(shí)候誤差趨于穩(wěn)定,MSE約等于0.0115。

    圖5 均方誤差MSE變化趨勢(shì)Fig.5 Trend of mean square error(MSE)

    圖6顯示了標(biāo)準(zhǔn)ET0-PMF56、APSO-ELM和ELM之間的關(guān)系。從圖6可以看出ELM模型預(yù)測(cè)出的參考作物蒸發(fā)蒸騰量決定系數(shù)R2只有0.949,這與標(biāo)準(zhǔn)ET0-PMF56的誤差比較大。而APSO-ELM模型預(yù)測(cè)的結(jié)果達(dá)到了0.981,這與ET0-PMF56具有很高的擬合度。

    圖6 ET0-PMF56、APSO-ELM和ELM結(jié)果對(duì)比Fig.6 ET0-PMF56, comparison of APSO-ELM and ELM results

    3.2 系統(tǒng)有效性分析

    根據(jù)上述實(shí)驗(yàn),番茄在其生長期內(nèi),實(shí)際作物蒸發(fā)蒸騰量ET0范圍可達(dá)到6 mm/d。現(xiàn)環(huán)境信息采集時(shí)間間隔為1 h,相鄰時(shí)間間隔內(nèi),土壤水分下降率變化的偏差在0%~15%,滴頭流量為1.6 L/h,灌溉時(shí)間最大可達(dá)到35 min。根據(jù)模糊規(guī)則設(shè)置的參數(shù)條件,得到灌溉時(shí)間T、土壤濕度下降率ΔE、實(shí)際作物蒸發(fā)蒸騰量ET0之間的變化關(guān)系。如圖7顯示了不同實(shí)際作物蒸發(fā)蒸騰量ET0在土壤濕度下降率ΔE相同的情況下帶來灌溉時(shí)間T的不同。從圖7可以看出灌溉時(shí)間T隨著土壤濕度下降率ΔE的提高在增加,在土壤濕度下降率ΔE較小的時(shí)候,灌溉時(shí)間T隨著不同實(shí)際作物蒸發(fā)蒸騰量ET0的不同有明顯的差異;在土壤濕度下降率ΔE較大的時(shí)候,灌溉時(shí)間T差異明顯變小,這時(shí)再加上實(shí)際作物蒸發(fā)蒸騰量ET0比較大,所以灌溉時(shí)間T就會(huì)增加;這時(shí)實(shí)際作物蒸發(fā)蒸騰量ET0比較小的時(shí)候,還是需要比較長的灌溉時(shí)間T,是因?yàn)槌薊T0以外,還有地表滲漏,地表徑流等耗水量因素造成土壤濕度下降率ΔE快速下降,所以需要較長的灌溉時(shí)間,這與表1的規(guī)則相符。

    圖7 灌溉時(shí)間T隨ΔE變化情況Fig.7 Irrigation T varied with ΔE

    圖8顯示了灌溉時(shí)間T隨實(shí)際作物需水量ET0的變化情況。從圖8可以得到土壤濕度下降率ΔE變化緩慢,因?yàn)镋T0是一個(gè)滯后的過程,作物蒸散的水分需要作物根系從周圍土壤中吸收,而吸收需要一個(gè)過程,所以造成土壤濕度的下降率ΔE的變化不是很明顯。另外,ET0變化比較大的時(shí)候,不管ΔE的變化多少,都需要比較長的灌溉時(shí)間T。例如,當(dāng)ET0達(dá)到5 mm/d的時(shí)候,灌溉時(shí)間T達(dá)到了25 min,這與表1的規(guī)則相符。

    圖8 灌溉時(shí)間T隨ET0變化情況Fig.8 Irrigation T changes with ET0

    4 結(jié) 語

    本文通過驗(yàn)證APSO-ELM模型計(jì)算實(shí)際作物蒸發(fā)蒸騰量ET0與ET0-PMF56計(jì)算的需水量具有很高的準(zhǔn)確度,減少了數(shù)據(jù)獲取成本過高以及計(jì)算復(fù)雜的問題。同時(shí)通過實(shí)驗(yàn)數(shù)據(jù)表明,ET0和土壤濕度下降率ΔE作為模糊邏輯系統(tǒng)的輸入,能夠準(zhǔn)確輸出灌溉時(shí)間T,該系統(tǒng)不僅考慮了土壤因素,而且考慮了作物本身生理活動(dòng)造成的水分流失,可以針對(duì)作物各個(gè)時(shí)期進(jìn)行按需灌溉。

    猜你喜歡
    土壤濕度灌溉作物
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    作物遭受霜凍該如何補(bǔ)救
    四種作物 北方種植有前景
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    基于51單片機(jī)控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    無人機(jī)遙感在作物監(jiān)測(cè)中的應(yīng)用與展望
    亚洲精品乱码久久久v下载方式| 国产免费一级a男人的天堂| 久久国产乱子免费精品| 国产视频首页在线观看| 亚洲精品第二区| 极品少妇高潮喷水抽搐| 全区人妻精品视频| 在线观看一区二区三区| 国产av国产精品国产| 麻豆精品久久久久久蜜桃| 一级毛片aaaaaa免费看小| 欧美日韩一区二区视频在线观看视频在线 | 日本一二三区视频观看| 国产熟女欧美一区二区| 亚洲av中文字字幕乱码综合| 青春草视频在线免费观看| 国产在线一区二区三区精| 有码 亚洲区| 黄片wwwwww| 日日啪夜夜撸| 免费高清在线观看视频在线观看| 男女啪啪激烈高潮av片| 永久免费av网站大全| 亚洲欧美精品自产自拍| av在线观看视频网站免费| 一区二区三区精品91| 欧美另类一区| 久久人人爽人人片av| 免费少妇av软件| 插逼视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| av在线播放精品| 亚洲精品成人av观看孕妇| 免费在线观看成人毛片| 亚洲人与动物交配视频| 赤兔流量卡办理| 久久久久久九九精品二区国产| 五月玫瑰六月丁香| 亚洲av免费高清在线观看| 国产 一区 欧美 日韩| 国产精品一区二区在线观看99| 欧美bdsm另类| 18禁在线无遮挡免费观看视频| 18禁动态无遮挡网站| 内射极品少妇av片p| 久久99蜜桃精品久久| 一级二级三级毛片免费看| 最近的中文字幕免费完整| 99久久精品一区二区三区| 国产午夜福利久久久久久| 寂寞人妻少妇视频99o| 高清欧美精品videossex| 免费观看无遮挡的男女| 嫩草影院新地址| 99九九线精品视频在线观看视频| 欧美一区二区亚洲| 成人特级av手机在线观看| 97在线视频观看| 欧美性感艳星| 国产精品熟女久久久久浪| 久久久色成人| 少妇 在线观看| 美女视频免费永久观看网站| 99久久人妻综合| 久久久久久久亚洲中文字幕| 久热这里只有精品99| 国产成人一区二区在线| 免费大片黄手机在线观看| 22中文网久久字幕| 色哟哟·www| 精品一区二区三区视频在线| 18+在线观看网站| 男女啪啪激烈高潮av片| 天美传媒精品一区二区| 国产精品伦人一区二区| 国产色爽女视频免费观看| 亚洲av二区三区四区| 日韩精品有码人妻一区| 亚洲国产精品成人综合色| av卡一久久| 91在线精品国自产拍蜜月| av女优亚洲男人天堂| 日韩亚洲欧美综合| a级毛色黄片| 看黄色毛片网站| 黄色日韩在线| 精品国产露脸久久av麻豆| av网站免费在线观看视频| av在线观看视频网站免费| 亚洲自拍偷在线| 18禁在线播放成人免费| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看| 精品久久久久久电影网| 国产亚洲精品久久久com| tube8黄色片| 国产精品久久久久久久电影| 午夜激情久久久久久久| 乱系列少妇在线播放| 国产午夜精品一二区理论片| 五月伊人婷婷丁香| 国产免费一区二区三区四区乱码| 黄色欧美视频在线观看| 国产精品偷伦视频观看了| 超碰av人人做人人爽久久| 亚洲精品中文字幕在线视频 | 18禁在线无遮挡免费观看视频| 国产在线一区二区三区精| 在线天堂最新版资源| 五月玫瑰六月丁香| 最近2019中文字幕mv第一页| www.av在线官网国产| 九九久久精品国产亚洲av麻豆| 国产免费一区二区三区四区乱码| 午夜免费观看性视频| 婷婷色麻豆天堂久久| 99九九线精品视频在线观看视频| 国产久久久一区二区三区| 色网站视频免费| 爱豆传媒免费全集在线观看| 欧美xxxx性猛交bbbb| 亚洲无线观看免费| 久久久成人免费电影| 国产精品久久久久久av不卡| 国产欧美日韩精品一区二区| 久久久午夜欧美精品| 欧美性感艳星| 99久久人妻综合| 亚洲第一区二区三区不卡| 国产成人精品婷婷| 国产精品久久久久久精品古装| 精品酒店卫生间| 亚洲国产av新网站| 熟女人妻精品中文字幕| 国产成人午夜福利电影在线观看| 最新中文字幕久久久久| 亚洲欧美中文字幕日韩二区| 亚洲av中文av极速乱| 国产精品爽爽va在线观看网站| 一级av片app| 网址你懂的国产日韩在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产伦在线观看视频一区| 婷婷色综合大香蕉| 美女内射精品一级片tv| 亚洲人成网站高清观看| 午夜福利高清视频| 欧美zozozo另类| videos熟女内射| 麻豆久久精品国产亚洲av| 精品一区在线观看国产| 高清毛片免费看| 欧美三级亚洲精品| 一边亲一边摸免费视频| 日韩一区二区三区影片| 免费大片黄手机在线观看| 日韩精品有码人妻一区| 免费观看的影片在线观看| 亚洲av中文字字幕乱码综合| 少妇人妻一区二区三区视频| 99久久精品国产国产毛片| 精品久久久噜噜| 亚洲一区二区三区欧美精品 | 人体艺术视频欧美日本| 偷拍熟女少妇极品色| 成年av动漫网址| 欧美三级亚洲精品| 久久久久久久精品精品| 国产日韩欧美在线精品| 久久国产乱子免费精品| 久久久久精品性色| 白带黄色成豆腐渣| 久久精品综合一区二区三区| 国产欧美日韩一区二区三区在线 | 国产人妻一区二区三区在| 亚洲精品一二三| 精品人妻视频免费看| 中国国产av一级| 国产爽快片一区二区三区| 肉色欧美久久久久久久蜜桃 | 国产成人精品福利久久| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 中文天堂在线官网| 夫妻午夜视频| 尤物成人国产欧美一区二区三区| 91午夜精品亚洲一区二区三区| 国产高清有码在线观看视频| 免费黄频网站在线观看国产| 深夜a级毛片| 日韩免费高清中文字幕av| 久热这里只有精品99| 久久国产乱子免费精品| 精品一区二区三区视频在线| 久久鲁丝午夜福利片| 中国三级夫妇交换| 在线亚洲精品国产二区图片欧美 | 国产淫片久久久久久久久| 精品久久久久久久人妻蜜臀av| 精品国产乱码久久久久久小说| 99久久精品国产国产毛片| 国产精品不卡视频一区二区| 日韩不卡一区二区三区视频在线| 韩国高清视频一区二区三区| 在线a可以看的网站| 97超视频在线观看视频| 国产精品久久久久久精品古装| 欧美日韩一区二区视频在线观看视频在线 | 高清欧美精品videossex| 亚洲精华国产精华液的使用体验| 免费电影在线观看免费观看| 干丝袜人妻中文字幕| 黄色一级大片看看| 国产黄a三级三级三级人| 又爽又黄无遮挡网站| 国产有黄有色有爽视频| 日本爱情动作片www.在线观看| 老女人水多毛片| 国产黄色免费在线视频| 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 一区二区三区精品91| 免费观看性生交大片5| 欧美老熟妇乱子伦牲交| av女优亚洲男人天堂| 别揉我奶头 嗯啊视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲内射少妇av| 亚洲国产最新在线播放| 国产爽快片一区二区三区| 99re6热这里在线精品视频| 国产永久视频网站| 国产一区亚洲一区在线观看| 丰满少妇做爰视频| 欧美日韩亚洲高清精品| 黑人高潮一二区| 精品人妻一区二区三区麻豆| av.在线天堂| 美女国产视频在线观看| 日韩免费高清中文字幕av| 别揉我奶头 嗯啊视频| 一个人看的www免费观看视频| 国产一区二区三区av在线| 在线观看美女被高潮喷水网站| 免费在线观看成人毛片| 中国美白少妇内射xxxbb| 亚洲精品日韩在线中文字幕| 天美传媒精品一区二区| 精品亚洲乱码少妇综合久久| 欧美xxⅹ黑人| 久久影院123| 一级黄片播放器| 亚洲国产日韩一区二区| 久久综合国产亚洲精品| kizo精华| 麻豆国产97在线/欧美| 晚上一个人看的免费电影| 永久免费av网站大全| 国产成人a∨麻豆精品| 亚洲精华国产精华液的使用体验| 精品久久国产蜜桃| 黄色欧美视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 蜜桃亚洲精品一区二区三区| 91久久精品电影网| 精品久久久精品久久久| 看十八女毛片水多多多| 久久ye,这里只有精品| 亚洲图色成人| 亚洲精品成人久久久久久| 国产精品一区www在线观看| 久久女婷五月综合色啪小说 | 观看免费一级毛片| 搡女人真爽免费视频火全软件| 国产乱来视频区| 永久免费av网站大全| videossex国产| 久久精品国产亚洲av涩爱| 国产黄a三级三级三级人| 国产在线一区二区三区精| 精品人妻视频免费看| 亚洲人与动物交配视频| 男的添女的下面高潮视频| 国产免费一区二区三区四区乱码| 国产精品.久久久| 狂野欧美激情性xxxx在线观看| 日本一本二区三区精品| 少妇人妻 视频| 美女高潮的动态| 成人毛片60女人毛片免费| 精品久久久久久久末码| 99久久九九国产精品国产免费| 狂野欧美激情性xxxx在线观看| 夜夜爽夜夜爽视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产免费福利视频在线观看| 禁无遮挡网站| 亚洲电影在线观看av| 亚洲aⅴ乱码一区二区在线播放| 中文资源天堂在线| 精品一区在线观看国产| 小蜜桃在线观看免费完整版高清| 夫妻午夜视频| 午夜视频国产福利| 精品久久国产蜜桃| 人妻 亚洲 视频| 亚洲国产欧美在线一区| 亚洲人成网站高清观看| 久久久久久伊人网av| 国产91av在线免费观看| 黄片无遮挡物在线观看| 国产视频内射| 国产精品成人在线| a级毛色黄片| 男女啪啪激烈高潮av片| 国产亚洲精品久久久com| 国产又色又爽无遮挡免| 久久精品国产自在天天线| 亚洲国产高清在线一区二区三| av在线亚洲专区| 国产成人免费观看mmmm| 午夜激情福利司机影院| 91狼人影院| 亚洲av免费在线观看| 欧美日韩亚洲高清精品| 亚洲精品亚洲一区二区| 亚洲精品国产色婷婷电影| 亚洲国产精品成人久久小说| 尾随美女入室| av免费在线看不卡| 亚洲精品中文字幕在线视频 | 成人午夜精彩视频在线观看| 国产精品一区二区性色av| 久久久久久国产a免费观看| 成年免费大片在线观看| 国产高清不卡午夜福利| a级毛片免费高清观看在线播放| 蜜桃久久精品国产亚洲av| 国产精品一区二区性色av| 噜噜噜噜噜久久久久久91| 啦啦啦啦在线视频资源| 欧美+日韩+精品| 国产一区二区在线观看日韩| 国产亚洲av嫩草精品影院| 久久久欧美国产精品| 黄片wwwwww| 中文乱码字字幕精品一区二区三区| 久久久精品免费免费高清| 成人黄色视频免费在线看| 欧美日韩国产mv在线观看视频 | 免费av观看视频| 少妇被粗大猛烈的视频| 亚洲最大成人手机在线| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 男插女下体视频免费在线播放| 久久久久久久久久成人| 国产精品女同一区二区软件| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品乱码久久久久久按摩| av女优亚洲男人天堂| 精品熟女少妇av免费看| 亚洲四区av| 全区人妻精品视频| 国产精品嫩草影院av在线观看| 国产午夜福利久久久久久| 最近手机中文字幕大全| 欧美xxⅹ黑人| 午夜爱爱视频在线播放| 国产乱人视频| 一级黄片播放器| 日韩人妻高清精品专区| 在线免费观看不下载黄p国产| 我的老师免费观看完整版| 日本欧美国产在线视频| 麻豆久久精品国产亚洲av| 老司机影院成人| 麻豆久久精品国产亚洲av| 国产精品嫩草影院av在线观看| 大话2 男鬼变身卡| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 狂野欧美白嫩少妇大欣赏| 男女边摸边吃奶| 日韩 亚洲 欧美在线| 亚洲精品亚洲一区二区| 性色av一级| 自拍偷自拍亚洲精品老妇| 日本黄色片子视频| eeuss影院久久| 欧美极品一区二区三区四区| 亚洲伊人久久精品综合| 国产大屁股一区二区在线视频| 我的女老师完整版在线观看| 亚洲av成人精品一区久久| 97超视频在线观看视频| 成人一区二区视频在线观看| 听说在线观看完整版免费高清| 国产精品人妻久久久久久| a级毛片免费高清观看在线播放| 久久久欧美国产精品| 亚洲精品国产成人久久av| 成人一区二区视频在线观看| 国产高清国产精品国产三级 | 国产女主播在线喷水免费视频网站| 在线 av 中文字幕| 午夜日本视频在线| 国产一级毛片在线| 草草在线视频免费看| 久久精品国产a三级三级三级| 亚洲综合色惰| 男人添女人高潮全过程视频| 久久久久久久午夜电影| 青春草视频在线免费观看| 高清av免费在线| 波多野结衣巨乳人妻| 国产精品爽爽va在线观看网站| 中文资源天堂在线| 日本爱情动作片www.在线观看| 精品酒店卫生间| 欧美xxxx黑人xx丫x性爽| 亚洲国产高清在线一区二区三| 男人和女人高潮做爰伦理| 1000部很黄的大片| a级毛色黄片| 哪个播放器可以免费观看大片| 国产精品蜜桃在线观看| 精品一区二区免费观看| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 免费人成在线观看视频色| 男女国产视频网站| 国产一区二区三区综合在线观看 | 热re99久久精品国产66热6| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av天美| 日韩大片免费观看网站| 丝袜喷水一区| 国产精品三级大全| 亚洲国产精品999| 色视频在线一区二区三区| 国产精品一区二区三区四区免费观看| 日韩免费高清中文字幕av| 国产高潮美女av| 天天一区二区日本电影三级| 亚洲经典国产精华液单| 国精品久久久久久国模美| 亚洲性久久影院| 男女国产视频网站| 日韩一区二区视频免费看| 久久久成人免费电影| 在线亚洲精品国产二区图片欧美 | 一边亲一边摸免费视频| 涩涩av久久男人的天堂| 免费电影在线观看免费观看| 老女人水多毛片| 最近最新中文字幕免费大全7| 亚洲精品国产av蜜桃| 午夜福利在线在线| 最近2019中文字幕mv第一页| 国产永久视频网站| 老司机影院成人| 成年女人看的毛片在线观看| 一本久久精品| 夫妻性生交免费视频一级片| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 一二三四中文在线观看免费高清| 日韩国内少妇激情av| 天美传媒精品一区二区| 不卡视频在线观看欧美| 一级二级三级毛片免费看| 91精品一卡2卡3卡4卡| 中国国产av一级| 国产伦理片在线播放av一区| 一级a做视频免费观看| 夜夜爽夜夜爽视频| 国产成年人精品一区二区| 亚洲av不卡在线观看| 国产永久视频网站| 久久久久九九精品影院| 特级一级黄色大片| xxx大片免费视频| 男人添女人高潮全过程视频| 国产一区亚洲一区在线观看| 亚洲av在线观看美女高潮| 男女下面进入的视频免费午夜| 超碰97精品在线观看| 中文欧美无线码| 亚洲国产精品国产精品| 欧美少妇被猛烈插入视频| 亚洲在久久综合| 精品久久久精品久久久| 亚洲av中文av极速乱| 中文字幕制服av| 国产成年人精品一区二区| 日韩在线高清观看一区二区三区| 亚洲精品影视一区二区三区av| 1000部很黄的大片| 蜜桃亚洲精品一区二区三区| 亚洲精品色激情综合| 日韩成人伦理影院| 一级毛片我不卡| 午夜福利高清视频| 亚洲精品国产av蜜桃| 中文精品一卡2卡3卡4更新| 人妻一区二区av| 一区二区三区免费毛片| 欧美xxⅹ黑人| 国产片特级美女逼逼视频| 欧美三级亚洲精品| 久久国内精品自在自线图片| 插逼视频在线观看| 国产免费一区二区三区四区乱码| h日本视频在线播放| 最近手机中文字幕大全| 一区二区三区精品91| 日日啪夜夜爽| 激情 狠狠 欧美| .国产精品久久| 日韩 亚洲 欧美在线| 26uuu在线亚洲综合色| 国产成人a区在线观看| 97在线视频观看| 少妇高潮的动态图| 国产 精品1| 看非洲黑人一级黄片| 大又大粗又爽又黄少妇毛片口| 免费大片18禁| 免费观看在线日韩| 婷婷色av中文字幕| 亚洲国产欧美在线一区| 伊人久久精品亚洲午夜| 高清欧美精品videossex| 国产精品福利在线免费观看| av在线老鸭窝| 天美传媒精品一区二区| 国产免费福利视频在线观看| 国产v大片淫在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 在线精品无人区一区二区三 | 亚洲精品一区蜜桃| 国产高清国产精品国产三级 | 国产成人免费观看mmmm| 久久亚洲国产成人精品v| 插阴视频在线观看视频| 自拍偷自拍亚洲精品老妇| 国产亚洲av片在线观看秒播厂| 日韩强制内射视频| 色网站视频免费| 99热这里只有是精品在线观看| 国产成人91sexporn| 国产精品不卡视频一区二区| 久久久久九九精品影院| 黄片无遮挡物在线观看| 欧美一级a爱片免费观看看| 少妇裸体淫交视频免费看高清| 免费观看在线日韩| 国产一级毛片在线| 国产男女内射视频| 亚洲成色77777| 午夜激情福利司机影院| 亚洲av.av天堂| 亚洲性久久影院| 在现免费观看毛片| 男女下面进入的视频免费午夜| 人妻 亚洲 视频| videossex国产| 国产精品偷伦视频观看了| 国产一区有黄有色的免费视频| 久久久国产一区二区| 国产精品不卡视频一区二区| 最近中文字幕高清免费大全6| 久久女婷五月综合色啪小说 | 久久这里有精品视频免费| 国产精品精品国产色婷婷| 欧美激情久久久久久爽电影| 18禁裸乳无遮挡免费网站照片| 成人毛片60女人毛片免费| 久久精品国产亚洲网站| 久久精品国产自在天天线| 亚洲美女搞黄在线观看| 五月玫瑰六月丁香| 91精品一卡2卡3卡4卡| 日日摸夜夜添夜夜爱| av.在线天堂| 91精品国产九色| 3wmmmm亚洲av在线观看| 人妻系列 视频| 成人毛片a级毛片在线播放| 性色avwww在线观看| 免费看光身美女| 免费大片18禁| 日日摸夜夜添夜夜添av毛片| 欧美激情在线99| 国产午夜福利久久久久久| 99热全是精品| 视频区图区小说| 欧美一级a爱片免费观看看| 亚洲av电影在线观看一区二区三区 | 三级国产精品欧美在线观看| 麻豆乱淫一区二区| 搡女人真爽免费视频火全软件| 国产色婷婷99| 嫩草影院精品99| 久久久久久久久久久免费av| 精品一区二区三区视频在线| 国内精品美女久久久久久| 三级国产精品片| 国产国拍精品亚洲av在线观看| 国产高清不卡午夜福利| 免费看日本二区| 精品一区在线观看国产| 国产免费福利视频在线观看| 别揉我奶头 嗯啊视频|