• <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)用與展望
    精品欧美国产一区二区三| 麻豆av噜噜一区二区三区| 麻豆成人午夜福利视频| 久久久久久久久大av| 小说图片视频综合网站| 欧美国产日韩亚洲一区| 91在线精品国自产拍蜜月| 99riav亚洲国产免费| 真人一进一出gif抽搐免费| 亚洲久久久久久中文字幕| 人人妻人人看人人澡| 99视频精品全部免费 在线| 国产一区二区三区av在线 | 亚洲最大成人中文| 国产亚洲欧美98| 18禁在线播放成人免费| 亚洲18禁久久av| 91麻豆精品激情在线观看国产| 男女边吃奶边做爰视频| 在现免费观看毛片| 欧美成人免费av一区二区三区| 免费大片18禁| 少妇被粗大猛烈的视频| 99热只有精品国产| 国产一区二区三区在线臀色熟女| 国产精品久久久久久av不卡| 久9热在线精品视频| 免费看日本二区| 日本在线视频免费播放| 看黄色毛片网站| 国语自产精品视频在线第100页| 国内精品美女久久久久久| 亚洲美女黄片视频| 日韩欧美三级三区| 少妇裸体淫交视频免费看高清| 综合色av麻豆| 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 特大巨黑吊av在线直播| 黄色配什么色好看| 一区福利在线观看| 亚洲七黄色美女视频| 嫩草影院入口| 久久精品国产清高在天天线| 国产又黄又爽又无遮挡在线| a级毛片a级免费在线| 久久久久久久久大av| 亚洲精品日韩av片在线观看| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 久久国产精品人妻蜜桃| 欧美日韩黄片免| 久久久久久久久大av| 久久午夜福利片| 九九爱精品视频在线观看| 欧美精品啪啪一区二区三区| 男人舔女人下体高潮全视频| 观看美女的网站| 桃色一区二区三区在线观看| 美女被艹到高潮喷水动态| 久久精品国产99精品国产亚洲性色| 黄色视频,在线免费观看| 婷婷色综合大香蕉| 国产高潮美女av| 一个人看视频在线观看www免费| 国产精品99久久久久久久久| 精品人妻视频免费看| 99九九线精品视频在线观看视频| 美女高潮的动态| 尾随美女入室| 日韩中字成人| 在线观看午夜福利视频| 人妻丰满熟妇av一区二区三区| 观看免费一级毛片| 国产综合懂色| 欧美色欧美亚洲另类二区| 午夜激情欧美在线| 久久亚洲真实| 亚洲熟妇中文字幕五十中出| 此物有八面人人有两片| 国产精品一区二区三区四区免费观看 | 国产高清视频在线观看网站| 久久久色成人| 国产精品不卡视频一区二区| 3wmmmm亚洲av在线观看| 精品久久久久久久末码| 中文字幕久久专区| 日本五十路高清| 露出奶头的视频| 午夜激情福利司机影院| 欧美+日韩+精品| 12—13女人毛片做爰片一| 欧美不卡视频在线免费观看| 国产精品一区二区性色av| 国产精品人妻久久久影院| 性插视频无遮挡在线免费观看| а√天堂www在线а√下载| 一进一出抽搐动态| 国产精品乱码一区二三区的特点| 亚洲综合色惰| 久9热在线精品视频| 麻豆av噜噜一区二区三区| 赤兔流量卡办理| 久久欧美精品欧美久久欧美| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 久久99热6这里只有精品| 91狼人影院| 亚洲精品日韩av片在线观看| 九色成人免费人妻av| 亚洲精品影视一区二区三区av| 国产亚洲欧美98| 天天一区二区日本电影三级| 黄色女人牲交| 国产黄色小视频在线观看| 一本久久中文字幕| 18+在线观看网站| 国产男靠女视频免费网站| 在线观看美女被高潮喷水网站| 欧美成人免费av一区二区三区| .国产精品久久| 国产精品av视频在线免费观看| 在线观看美女被高潮喷水网站| a级毛片免费高清观看在线播放| 我要看日韩黄色一级片| 国产欧美日韩一区二区精品| 欧美日韩国产亚洲二区| 免费av毛片视频| 亚洲欧美日韩高清专用| 欧美+日韩+精品| 久久精品夜夜夜夜夜久久蜜豆| 日本 欧美在线| 尤物成人国产欧美一区二区三区| 亚洲男人的天堂狠狠| 久久久国产成人精品二区| 51国产日韩欧美| 国产精品永久免费网站| 成年女人看的毛片在线观看| 嫩草影院精品99| 极品教师在线视频| 免费在线观看影片大全网站| 99国产精品一区二区蜜桃av| 成人美女网站在线观看视频| www日本黄色视频网| 日韩欧美精品免费久久| 美女被艹到高潮喷水动态| 亚洲美女黄片视频| 精品久久久久久久末码| 真人一进一出gif抽搐免费| 国产精品久久久久久久久免| 全区人妻精品视频| 高清日韩中文字幕在线| 我要搜黄色片| 成人亚洲精品av一区二区| 国产视频一区二区在线看| 高清日韩中文字幕在线| 嫩草影院新地址| 美女高潮的动态| 成人高潮视频无遮挡免费网站| 亚洲美女视频黄频| 18禁在线播放成人免费| 亚洲欧美日韩高清专用| 成人高潮视频无遮挡免费网站| 久久久成人免费电影| 18禁裸乳无遮挡免费网站照片| 淫秽高清视频在线观看| 午夜爱爱视频在线播放| 免费在线观看日本一区| 色综合色国产| 国产精品一区二区三区四区久久| 成年女人毛片免费观看观看9| 久久久久国产精品人妻aⅴ院| 国产av在哪里看| 国国产精品蜜臀av免费| 亚洲精品乱码久久久v下载方式| av.在线天堂| 亚洲自偷自拍三级| 毛片一级片免费看久久久久 | 国产高潮美女av| 亚洲无线在线观看| 在线免费十八禁| 国产精品电影一区二区三区| 又紧又爽又黄一区二区| 国产精品久久久久久久电影| 麻豆久久精品国产亚洲av| 村上凉子中文字幕在线| 啪啪无遮挡十八禁网站| 色综合婷婷激情| 搡女人真爽免费视频火全软件 | 能在线免费观看的黄片| av天堂中文字幕网| 深爱激情五月婷婷| 精品久久久噜噜| 国产一区二区三区av在线 | 午夜免费男女啪啪视频观看 | 亚洲va日本ⅴa欧美va伊人久久| 成人三级黄色视频| 春色校园在线视频观看| 欧美国产日韩亚洲一区| 国产在线男女| 2021天堂中文幕一二区在线观| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 久久久久久久久大av| 久久精品国产亚洲网站| 国产一级毛片七仙女欲春2| 日本精品一区二区三区蜜桃| 午夜亚洲福利在线播放| 一夜夜www| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品久久久com| 亚洲在线自拍视频| avwww免费| 成人av在线播放网站| 国产精品亚洲美女久久久| 国产精品免费一区二区三区在线| 国产精品伦人一区二区| 22中文网久久字幕| 日本 av在线| 久久人人精品亚洲av| 能在线免费观看的黄片| 联通29元200g的流量卡| 国产乱人视频| 超碰av人人做人人爽久久| 极品教师在线视频| 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 一本一本综合久久| 精品午夜福利视频在线观看一区| 又黄又爽又免费观看的视频| 黄片wwwwww| 两个人视频免费观看高清| 老司机福利观看| 老熟妇乱子伦视频在线观看| 久久亚洲精品不卡| 制服丝袜大香蕉在线| 午夜精品久久久久久毛片777| 99热这里只有是精品50| 亚洲人成网站在线播放欧美日韩| 久久精品国产自在天天线| 成人午夜高清在线视频| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 3wmmmm亚洲av在线观看| а√天堂www在线а√下载| 又爽又黄无遮挡网站| 亚洲欧美日韩东京热| 国产 一区精品| 日本a在线网址| 亚洲av电影不卡..在线观看| 欧美精品啪啪一区二区三区| 男人舔女人下体高潮全视频| 啪啪无遮挡十八禁网站| 精品99又大又爽又粗少妇毛片 | 国产私拍福利视频在线观看| 一区二区三区免费毛片| 国产久久久一区二区三区| 网址你懂的国产日韩在线| 国产高清激情床上av| 少妇裸体淫交视频免费看高清| 最新在线观看一区二区三区| 欧美日本视频| 99国产极品粉嫩在线观看| 听说在线观看完整版免费高清| 一进一出抽搐动态| 国产综合懂色| 久久精品91蜜桃| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| 少妇人妻精品综合一区二区 | 亚洲欧美日韩无卡精品| 国国产精品蜜臀av免费| 午夜精品久久久久久毛片777| 真实男女啪啪啪动态图| 亚洲经典国产精华液单| 深夜a级毛片| 日本免费a在线| 两人在一起打扑克的视频| 18+在线观看网站| 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验 | 蜜桃亚洲精品一区二区三区| 日韩国内少妇激情av| 欧美zozozo另类| 亚洲一区高清亚洲精品| 午夜免费男女啪啪视频观看 | h日本视频在线播放| 午夜福利在线在线| 亚洲国产欧美人成| 午夜福利在线观看吧| 桃色一区二区三区在线观看| 国产色爽女视频免费观看| 婷婷亚洲欧美| 2021天堂中文幕一二区在线观| 国产男人的电影天堂91| 能在线免费观看的黄片| 欧美潮喷喷水| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 夜夜夜夜夜久久久久| 可以在线观看毛片的网站| 身体一侧抽搐| 亚洲一级一片aⅴ在线观看| 中亚洲国语对白在线视频| 色尼玛亚洲综合影院| 黄色一级大片看看| 亚洲自拍偷在线| 成年女人永久免费观看视频| 最后的刺客免费高清国语| 午夜久久久久精精品| 免费无遮挡裸体视频| 搡老岳熟女国产| 欧美性猛交╳xxx乱大交人| av女优亚洲男人天堂| 国内精品久久久久精免费| 永久网站在线| 少妇的逼好多水| bbb黄色大片| 免费av观看视频| 亚洲av中文av极速乱 | 乱码一卡2卡4卡精品| 内射极品少妇av片p| 精品久久久久久,| 中文字幕熟女人妻在线| 久久99热6这里只有精品| 日韩国内少妇激情av| 欧美日本亚洲视频在线播放| 乱人视频在线观看| x7x7x7水蜜桃| 国产亚洲精品av在线| 免费看a级黄色片| 有码 亚洲区| 51国产日韩欧美| 一个人看视频在线观看www免费| av在线观看视频网站免费| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区精品| 国产探花在线观看一区二区| 国产一区二区三区视频了| 国产精品一区二区三区四区久久| av女优亚洲男人天堂| netflix在线观看网站| 久久人人精品亚洲av| 国产精品女同一区二区软件 | 亚洲国产高清在线一区二区三| 精品一区二区免费观看| 久久久色成人| 大又大粗又爽又黄少妇毛片口| 我的老师免费观看完整版| 长腿黑丝高跟| 成人性生交大片免费视频hd| 国产伦在线观看视频一区| 91av网一区二区| 国产精品久久久久久av不卡| av黄色大香蕉| 国产女主播在线喷水免费视频网站 | 午夜免费激情av| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 久久人人爽人人爽人人片va| 99久久精品国产国产毛片| 啦啦啦韩国在线观看视频| 欧美日韩综合久久久久久 | 97热精品久久久久久| 黄色欧美视频在线观看| 日本五十路高清| 久久6这里有精品| 久久人妻av系列| 男女啪啪激烈高潮av片| 国产精品爽爽va在线观看网站| 91麻豆av在线| 国产熟女欧美一区二区| 长腿黑丝高跟| 国产视频内射| 又紧又爽又黄一区二区| 国产欧美日韩精品亚洲av| 国内精品久久久久久久电影| 国产精品一区二区三区四区免费观看 | 亚洲av五月六月丁香网| 久久久成人免费电影| 国产真实乱freesex| 99热精品在线国产| 人人妻,人人澡人人爽秒播| 一区二区三区四区激情视频 | 国产男人的电影天堂91| 亚洲人成伊人成综合网2020| 国产免费男女视频| 国产精品免费一区二区三区在线| bbb黄色大片| 一个人看的www免费观看视频| 免费不卡的大黄色大毛片视频在线观看 | 国产成人一区二区在线| 亚洲精品成人久久久久久| 亚洲精品久久国产高清桃花| 成人一区二区视频在线观看| 无遮挡黄片免费观看| 国产精华一区二区三区| 日日摸夜夜添夜夜添小说| 人人妻,人人澡人人爽秒播| 在线观看美女被高潮喷水网站| 久久草成人影院| 少妇熟女aⅴ在线视频| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 久久久久久久亚洲中文字幕| 一个人免费在线观看电影| 大又大粗又爽又黄少妇毛片口| 此物有八面人人有两片| 一个人观看的视频www高清免费观看| 级片在线观看| 色哟哟哟哟哟哟| 综合色av麻豆| 国内久久婷婷六月综合欲色啪| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 性色avwww在线观看| 亚洲,欧美,日韩| 精品乱码久久久久久99久播| 给我免费播放毛片高清在线观看| 欧美潮喷喷水| 精品午夜福利在线看| 亚洲人成网站在线播| 22中文网久久字幕| 美女xxoo啪啪120秒动态图| 欧美中文日本在线观看视频| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 免费搜索国产男女视频| 免费人成在线观看视频色| 又粗又爽又猛毛片免费看| 精品免费久久久久久久清纯| 国产高清三级在线| 亚洲自拍偷在线| 一边摸一边抽搐一进一小说| 日日夜夜操网爽| 欧美bdsm另类| 18禁黄网站禁片午夜丰满| av在线亚洲专区| 国产精品,欧美在线| 麻豆av噜噜一区二区三区| 黄色日韩在线| 男女边吃奶边做爰视频| 亚洲精品456在线播放app | 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 亚洲欧美日韩无卡精品| 男插女下体视频免费在线播放| 久久久久国产精品人妻aⅴ院| 色吧在线观看| 色av中文字幕| 最近最新中文字幕大全电影3| 国产亚洲精品av在线| av中文乱码字幕在线| 亚洲最大成人av| 亚洲一区高清亚洲精品| 国产精品一区二区性色av| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 日韩强制内射视频| 久久久久久久久久久丰满 | 丰满的人妻完整版| 久久99热6这里只有精品| 老司机深夜福利视频在线观看| 88av欧美| 欧美中文日本在线观看视频| 中文亚洲av片在线观看爽| 国产精品精品国产色婷婷| 午夜久久久久精精品| 国产男靠女视频免费网站| 国产成人影院久久av| 搡老熟女国产l中国老女人| 丰满的人妻完整版| 中文字幕精品亚洲无线码一区| 国产高清不卡午夜福利| 国产一区二区三区在线臀色熟女| 精品久久久久久成人av| 悠悠久久av| 亚洲国产精品成人综合色| 免费一级毛片在线播放高清视频| 97人妻精品一区二区三区麻豆| 禁无遮挡网站| 亚洲国产精品合色在线| 我的女老师完整版在线观看| 夜夜夜夜夜久久久久| 免费观看在线日韩| 亚洲av熟女| 麻豆一二三区av精品| 精品国产三级普通话版| 51国产日韩欧美| 精品福利观看| or卡值多少钱| 精品久久久久久久久久久久久| 国产精品福利在线免费观看| 成年人黄色毛片网站| 欧美激情久久久久久爽电影| 久久人人爽人人爽人人片va| 亚洲av电影不卡..在线观看| 亚洲 国产 在线| 欧美成人一区二区免费高清观看| 欧美高清成人免费视频www| 国产私拍福利视频在线观看| 一级av片app| av中文乱码字幕在线| 欧美成人一区二区免费高清观看| 91狼人影院| 亚洲自偷自拍三级| 别揉我奶头 嗯啊视频| 国产 一区精品| 午夜激情福利司机影院| 精品一区二区三区视频在线观看免费| 日日干狠狠操夜夜爽| 国产亚洲精品av在线| 国产黄a三级三级三级人| 国产av一区在线观看免费| 免费看a级黄色片| 免费观看的影片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 床上黄色一级片| 尤物成人国产欧美一区二区三区| 搡老岳熟女国产| 国产视频内射| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频| 成人高潮视频无遮挡免费网站| 亚洲最大成人av| 99热精品在线国产| 亚洲精品色激情综合| 中文字幕人妻熟人妻熟丝袜美| 色在线成人网| 亚洲国产色片| 欧美日韩黄片免| 亚洲va在线va天堂va国产| 黄色一级大片看看| 超碰av人人做人人爽久久| 色综合婷婷激情| 少妇的逼好多水| 99热这里只有是精品50| 99热网站在线观看| 欧美另类亚洲清纯唯美| 午夜影院日韩av| 亚洲成a人片在线一区二区| 乱系列少妇在线播放| 精品一区二区三区人妻视频| 亚洲av免费高清在线观看| 一区二区三区激情视频| 日韩高清综合在线| 欧美精品啪啪一区二区三区| or卡值多少钱| 在线国产一区二区在线| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 国产精品一及| 国产久久久一区二区三区| 久久精品久久久久久噜噜老黄 | 免费看av在线观看网站| 免费看a级黄色片| 久久精品国产亚洲网站| 免费看a级黄色片| 午夜激情福利司机影院| 午夜福利高清视频| 人妻少妇偷人精品九色| 人妻夜夜爽99麻豆av| 在现免费观看毛片| 永久网站在线| 亚洲最大成人手机在线| 国产极品精品免费视频能看的| 18禁裸乳无遮挡免费网站照片| 又紧又爽又黄一区二区| 成人无遮挡网站| 亚洲av二区三区四区| 啪啪无遮挡十八禁网站| 亚洲av二区三区四区| 99久久无色码亚洲精品果冻| 成年版毛片免费区| 欧美在线一区亚洲| 亚洲av五月六月丁香网| 99久久无色码亚洲精品果冻| 精品福利观看| 一区二区三区免费毛片| 小说图片视频综合网站| 免费人成视频x8x8入口观看| 白带黄色成豆腐渣| 亚洲精品影视一区二区三区av| 三级男女做爰猛烈吃奶摸视频| 亚洲精品久久国产高清桃花| 日韩中文字幕欧美一区二区| 日韩欧美精品免费久久| 亚洲国产高清在线一区二区三| 给我免费播放毛片高清在线观看| 国产午夜福利久久久久久| 亚洲av免费高清在线观看| 我的老师免费观看完整版| 麻豆成人午夜福利视频| 一进一出抽搐gif免费好疼| 简卡轻食公司| 22中文网久久字幕| 免费在线观看成人毛片| 午夜福利18| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄 | 91精品国产九色| 男人舔女人下体高潮全视频| 内地一区二区视频在线| 国产伦精品一区二区三区视频9| 级片在线观看| 亚洲18禁久久av| 大又大粗又爽又黄少妇毛片口| 女生性感内裤真人,穿戴方法视频| 成人精品一区二区免费| 深爱激情五月婷婷| 久久久久国产精品人妻aⅴ院|