• <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电影在线播放| www日本在线高清视频| 精品少妇一区二区三区视频日本电影| 婷婷精品国产亚洲av在线| 满18在线观看网站| 久久久精品国产亚洲av高清涩受| 老熟妇仑乱视频hdxx| 无遮挡黄片免费观看| 日韩三级视频一区二区三区| 亚洲人成电影免费在线| 国产激情欧美一区二区| 三上悠亚av全集在线观看| 亚洲男人天堂网一区| 免费看十八禁软件| 亚洲情色 制服丝袜| 免费观看人在逋| 热99re8久久精品国产| 国产片内射在线| 9热在线视频观看99| 女人高潮潮喷娇喘18禁视频| 91成年电影在线观看| 午夜精品久久久久久毛片777| 精品第一国产精品| 夜夜看夜夜爽夜夜摸 | av片东京热男人的天堂| 国产成人av激情在线播放| 变态另类成人亚洲欧美熟女 | 最近最新中文字幕大全免费视频| 国产精品久久电影中文字幕| 亚洲五月色婷婷综合| 亚洲精品久久成人aⅴ小说| 桃色一区二区三区在线观看| 又黄又粗又硬又大视频| 巨乳人妻的诱惑在线观看| 国产有黄有色有爽视频| 成人三级做爰电影| 亚洲欧美精品综合一区二区三区| 国产91精品成人一区二区三区| 一级,二级,三级黄色视频| 午夜久久久在线观看| 欧美激情久久久久久爽电影 | 久久亚洲真实| 欧美大码av| 久热爱精品视频在线9| 一区福利在线观看| 丝袜人妻中文字幕| 精品国产乱子伦一区二区三区| 精品国产亚洲在线| 99精品欧美一区二区三区四区| 日韩大码丰满熟妇| 国产av一区在线观看免费| 水蜜桃什么品种好| 国产欧美日韩一区二区三| 午夜福利在线观看吧| 午夜成年电影在线免费观看| 老司机在亚洲福利影院| 久久中文看片网| 免费搜索国产男女视频| 中文字幕色久视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利在线免费观看网站| 性欧美人与动物交配| 韩国av一区二区三区四区| 乱人伦中国视频| 久久久久久久午夜电影 | 久久国产精品男人的天堂亚洲| 露出奶头的视频| 日本vs欧美在线观看视频| 热re99久久国产66热| 亚洲专区中文字幕在线| 香蕉丝袜av| 久久精品国产99精品国产亚洲性色 | 韩国精品一区二区三区| 天堂√8在线中文| 久久香蕉激情| 国产伦一二天堂av在线观看| 欧美黑人精品巨大| 高清毛片免费观看视频网站 | 国产又色又爽无遮挡免费看| 18禁裸乳无遮挡免费网站照片 | 一区二区三区精品91| 最近最新中文字幕大全电影3 | 久久午夜亚洲精品久久| 美女福利国产在线| 欧美av亚洲av综合av国产av| 男女之事视频高清在线观看| 久久伊人香网站| 亚洲专区国产一区二区| 老汉色av国产亚洲站长工具| 免费不卡黄色视频| 热99国产精品久久久久久7| 无人区码免费观看不卡| 久久人人爽av亚洲精品天堂| 999精品在线视频| 欧美日本中文国产一区发布| 亚洲免费av在线视频| 日韩国内少妇激情av| 国产成人系列免费观看| 十八禁人妻一区二区| 国产在线精品亚洲第一网站| 乱人伦中国视频| 日本免费a在线| 人妻久久中文字幕网| 欧美精品亚洲一区二区| 久久香蕉国产精品| 身体一侧抽搐| 色综合站精品国产| 欧美另类亚洲清纯唯美| 老司机靠b影院| 在线天堂中文资源库| 免费在线观看亚洲国产| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 电影成人av| 国内久久婷婷六月综合欲色啪| 亚洲五月色婷婷综合| 91精品三级在线观看| 99在线视频只有这里精品首页| 麻豆一二三区av精品| 男女午夜视频在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美+亚洲+日韩+国产| 淫秽高清视频在线观看| av网站在线播放免费| 亚洲五月色婷婷综合| 老汉色av国产亚洲站长工具| 欧美一级毛片孕妇| 亚洲三区欧美一区| 国产精品av久久久久免费| 中出人妻视频一区二区| 身体一侧抽搐| 亚洲精品中文字幕一二三四区| 精品久久久久久成人av| 黄网站色视频无遮挡免费观看| 在线天堂中文资源库| 女人精品久久久久毛片| 精品国产超薄肉色丝袜足j| 天堂动漫精品| 国产成人欧美在线观看| 97人妻天天添夜夜摸| 午夜免费观看网址| 午夜精品久久久久久毛片777| 精品卡一卡二卡四卡免费| 国产一区在线观看成人免费| 欧美日韩视频精品一区| 99久久99久久久精品蜜桃| 巨乳人妻的诱惑在线观看| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 日韩欧美国产一区二区入口| а√天堂www在线а√下载| 久久伊人香网站| 精品国产超薄肉色丝袜足j| 久久久久久大精品| 国产欧美日韩一区二区三区在线| av在线天堂中文字幕 | 男人舔女人的私密视频| 夜夜爽天天搞| 一级毛片女人18水好多| 国产男靠女视频免费网站| 99国产精品一区二区蜜桃av| 黑人欧美特级aaaaaa片| 久久国产精品男人的天堂亚洲| 超色免费av| 中文字幕高清在线视频| 免费av毛片视频| 一级作爱视频免费观看| 97人妻天天添夜夜摸| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| 亚洲成人久久性| 黄频高清免费视频| 国产99久久九九免费精品| 免费看a级黄色片| cao死你这个sao货| 国产成人免费无遮挡视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国内视频| av电影中文网址| 午夜免费观看网址| 免费在线观看亚洲国产| 亚洲精品在线美女| 人成视频在线观看免费观看| 一区二区三区国产精品乱码| 亚洲成人国产一区在线观看| 国产麻豆69| 色老头精品视频在线观看| 亚洲第一青青草原| www国产在线视频色| 日本黄色视频三级网站网址| 最近最新中文字幕大全电影3 | 超碰成人久久| 一夜夜www| 精品午夜福利视频在线观看一区| 精品少妇一区二区三区视频日本电影| 免费女性裸体啪啪无遮挡网站| 中文字幕最新亚洲高清| 久久国产乱子伦精品免费另类| 日韩大码丰满熟妇| 欧美精品啪啪一区二区三区| 久久久久久久午夜电影 | 日本黄色日本黄色录像| 黄色丝袜av网址大全| 人人妻,人人澡人人爽秒播| 伦理电影免费视频| 国产又爽黄色视频| 国产色视频综合| 欧美乱码精品一区二区三区| 亚洲久久久国产精品| 夜夜夜夜夜久久久久| 亚洲欧美激情在线| 亚洲七黄色美女视频| 黄片播放在线免费| 久久亚洲精品不卡| 少妇的丰满在线观看| 999精品在线视频| 99国产精品一区二区蜜桃av| 国产精品电影一区二区三区| 亚洲激情在线av| 久久这里只有精品19| 欧美日韩乱码在线| 日日摸夜夜添夜夜添小说| 婷婷六月久久综合丁香| 老司机午夜十八禁免费视频| 亚洲精品在线美女| 黄频高清免费视频| 色播在线永久视频| 久久天堂一区二区三区四区| aaaaa片日本免费| 亚洲精品av麻豆狂野| 国产亚洲精品第一综合不卡| 香蕉国产在线看| 久久精品亚洲熟妇少妇任你| 黑丝袜美女国产一区| 成人影院久久| 精品国产一区二区三区四区第35| x7x7x7水蜜桃| 母亲3免费完整高清在线观看| 精品国内亚洲2022精品成人| 国产成人系列免费观看| 亚洲欧美激情综合另类| 免费在线观看影片大全网站| 视频在线观看一区二区三区| 高清欧美精品videossex| 中文字幕人妻丝袜一区二区| 欧美成人免费av一区二区三区| 亚洲国产欧美日韩在线播放| 这个男人来自地球电影免费观看| 欧美国产精品va在线观看不卡| 一区二区日韩欧美中文字幕| 精品国产一区二区久久| 国产99白浆流出| 亚洲精品中文字幕一二三四区| 日日摸夜夜添夜夜添小说| av视频免费观看在线观看| 老鸭窝网址在线观看| 香蕉国产在线看| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 我的亚洲天堂| a级片在线免费高清观看视频| 国产单亲对白刺激| 最好的美女福利视频网| 99热国产这里只有精品6| 欧美成人午夜精品| 中国美女看黄片| 超碰成人久久| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| av视频免费观看在线观看| 首页视频小说图片口味搜索| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩一级在线毛片| 精品一区二区三区视频在线观看免费 | 午夜视频精品福利| 一区二区三区精品91| 美女国产高潮福利片在线看| 性色av乱码一区二区三区2| 精品高清国产在线一区| 欧美黑人欧美精品刺激| 多毛熟女@视频| 亚洲自拍偷在线| 亚洲人成电影免费在线| 在线观看www视频免费| 啪啪无遮挡十八禁网站| 又黄又粗又硬又大视频| 日韩欧美免费精品| 午夜91福利影院| 免费高清视频大片| 亚洲第一欧美日韩一区二区三区| 99久久久亚洲精品蜜臀av| 色精品久久人妻99蜜桃| av超薄肉色丝袜交足视频| 一本综合久久免费| 久久中文字幕人妻熟女| 91在线观看av| 亚洲av电影在线进入| 日韩大尺度精品在线看网址 | 又大又爽又粗| 亚洲五月色婷婷综合| 亚洲一码二码三码区别大吗| 欧美 亚洲 国产 日韩一| 桃红色精品国产亚洲av| 午夜两性在线视频| 十八禁人妻一区二区| 老司机在亚洲福利影院| 啪啪无遮挡十八禁网站| 老熟妇乱子伦视频在线观看| 亚洲一码二码三码区别大吗| 久久久久国产精品人妻aⅴ院| 亚洲性夜色夜夜综合| 国产99久久九九免费精品| 日韩欧美三级三区| 国产精品久久电影中文字幕| 欧美黑人精品巨大| 免费av中文字幕在线| 精品一区二区三卡| 色综合欧美亚洲国产小说| av天堂在线播放| 18禁黄网站禁片午夜丰满| 动漫黄色视频在线观看| 国产色视频综合| 国产精品亚洲一级av第二区| 国产黄色免费在线视频| 午夜两性在线视频| a在线观看视频网站| 九色亚洲精品在线播放| 麻豆av在线久日| 亚洲av五月六月丁香网| 男女下面插进去视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 女生性感内裤真人,穿戴方法视频| 夜夜躁狠狠躁天天躁| 丰满迷人的少妇在线观看| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 神马国产精品三级电影在线观看 | 午夜a级毛片| 女人被狂操c到高潮| 日韩精品青青久久久久久| 中文字幕人妻熟女乱码| 999久久久精品免费观看国产| 麻豆av在线久日| 国产精品一区二区精品视频观看| 国产精品久久视频播放| 女同久久另类99精品国产91| 男人的好看免费观看在线视频 | bbb黄色大片| 不卡一级毛片| 成人黄色视频免费在线看| 激情视频va一区二区三区| 久久久久九九精品影院| 亚洲精品在线美女| 真人一进一出gif抽搐免费| 久久香蕉精品热| 中亚洲国语对白在线视频| 欧美日韩av久久| 久久99一区二区三区| 中文字幕精品免费在线观看视频| 日日干狠狠操夜夜爽| 国产精品久久久久成人av| 亚洲性夜色夜夜综合| 在线观看舔阴道视频| 国产免费现黄频在线看| 国产亚洲精品综合一区在线观看 | 少妇 在线观看| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 精品福利永久在线观看| 麻豆成人av在线观看| 777久久人妻少妇嫩草av网站| 亚洲成人免费电影在线观看| 成年人黄色毛片网站| 不卡一级毛片| 一进一出抽搐动态| 午夜免费成人在线视频| 精品欧美一区二区三区在线| 久热这里只有精品99| 黑人猛操日本美女一级片| 另类亚洲欧美激情| 淫妇啪啪啪对白视频| www.999成人在线观看| 欧美最黄视频在线播放免费 | 精品第一国产精品| 91精品三级在线观看| 欧美最黄视频在线播放免费 | 黑人欧美特级aaaaaa片| 19禁男女啪啪无遮挡网站| 日日爽夜夜爽网站| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| 日韩人妻精品一区2区三区| 欧美精品亚洲一区二区| 亚洲精品国产区一区二| 精品国产一区二区久久| 无遮挡黄片免费观看| 成人国语在线视频| 美女大奶头视频| 久久久国产一区二区| 婷婷六月久久综合丁香| 国产在线精品亚洲第一网站| 五月开心婷婷网| 午夜日韩欧美国产| 18禁美女被吸乳视频| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 一边摸一边做爽爽视频免费| 别揉我奶头~嗯~啊~动态视频| 91字幕亚洲| 极品人妻少妇av视频| 无遮挡黄片免费观看| 成人三级做爰电影| 18禁国产床啪视频网站| 可以在线观看毛片的网站| av有码第一页| 成人永久免费在线观看视频| 欧美乱码精品一区二区三区| 国产成人精品久久二区二区91| 国产精品1区2区在线观看.| 午夜免费激情av| 亚洲精品在线美女| 国产精品秋霞免费鲁丝片| 亚洲美女黄片视频| 欧美中文日本在线观看视频| 在线观看一区二区三区| 丰满迷人的少妇在线观看| 久久精品影院6| 亚洲中文日韩欧美视频| 免费观看人在逋| 激情在线观看视频在线高清| 青草久久国产| 999久久久国产精品视频| 精品无人区乱码1区二区| 亚洲美女黄片视频| 夫妻午夜视频| 国产一区二区三区综合在线观看| 搡老熟女国产l中国老女人| 久久人妻熟女aⅴ| 老司机在亚洲福利影院| 真人一进一出gif抽搐免费| 97碰自拍视频| 精品日产1卡2卡| 精品人妻在线不人妻| 亚洲情色 制服丝袜| 琪琪午夜伦伦电影理论片6080| 亚洲欧美一区二区三区久久| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| tocl精华| 大型黄色视频在线免费观看| 久久久久久亚洲精品国产蜜桃av| 一区福利在线观看| 精品国产乱子伦一区二区三区| videosex国产| 久久精品亚洲熟妇少妇任你| 嫁个100分男人电影在线观看| 欧美日韩亚洲高清精品| 两人在一起打扑克的视频| 国产精品香港三级国产av潘金莲| 一区福利在线观看| 精品一区二区三区av网在线观看| 午夜福利在线免费观看网站| 欧美日韩av久久| 久久午夜亚洲精品久久| 国产三级在线视频| 美女扒开内裤让男人捅视频| 91精品三级在线观看| 亚洲午夜理论影院| av网站免费在线观看视频| 激情视频va一区二区三区| 午夜老司机福利片| 久久国产精品男人的天堂亚洲| 国产精品一区二区在线不卡| 午夜福利免费观看在线| 亚洲av五月六月丁香网| 成人黄色视频免费在线看| 最好的美女福利视频网| 国产片内射在线| 亚洲精品国产精品久久久不卡| 国产伦一二天堂av在线观看| www.www免费av| 91字幕亚洲| 成人亚洲精品av一区二区 | 在线观看舔阴道视频| 精品福利永久在线观看| 黄色 视频免费看| 亚洲一区高清亚洲精品| 18禁国产床啪视频网站| 黑人操中国人逼视频| 久久久久久久久久久久大奶| 国产高清视频在线播放一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲激情在线av| 中文字幕高清在线视频| av超薄肉色丝袜交足视频| 欧美日韩亚洲高清精品| 欧美成人免费av一区二区三区| 99热只有精品国产| 精品国产国语对白av| 夜夜爽天天搞| 日本免费一区二区三区高清不卡 | 国产乱人伦免费视频| 久9热在线精品视频| 天天躁狠狠躁夜夜躁狠狠躁| 9色porny在线观看| 中文字幕另类日韩欧美亚洲嫩草| 一区二区三区激情视频| 97超级碰碰碰精品色视频在线观看| 久久精品国产亚洲av香蕉五月| 免费观看精品视频网站| 久久99一区二区三区| 色哟哟哟哟哟哟| 国产精品亚洲av一区麻豆| 在线永久观看黄色视频| 日韩三级视频一区二区三区| 亚洲国产精品合色在线| 麻豆国产av国片精品| 丁香欧美五月| 黄色a级毛片大全视频| 久久中文字幕一级| 亚洲 欧美 日韩 在线 免费| 女人精品久久久久毛片| 精品乱码久久久久久99久播| 亚洲一区二区三区不卡视频| 视频在线观看一区二区三区| 国产av在哪里看| 中文字幕最新亚洲高清| 日韩视频一区二区在线观看| 一级a爱视频在线免费观看| 99香蕉大伊视频| 99re在线观看精品视频| 日日爽夜夜爽网站| 老鸭窝网址在线观看| 自线自在国产av| √禁漫天堂资源中文www| 国产精品九九99| 国产乱人伦免费视频| 18禁黄网站禁片午夜丰满| 国产高清激情床上av| 日日爽夜夜爽网站| 电影成人av| 亚洲av成人不卡在线观看播放网| 国内久久婷婷六月综合欲色啪| 国产成人精品在线电影| 亚洲五月色婷婷综合| 色婷婷av一区二区三区视频| 在线免费观看的www视频| 久久精品亚洲精品国产色婷小说| 如日韩欧美国产精品一区二区三区| 久久久国产欧美日韩av| 高潮久久久久久久久久久不卡| 亚洲欧美一区二区三区黑人| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区 | 色老头精品视频在线观看| 欧美日韩一级在线毛片| 男女做爰动态图高潮gif福利片 | 激情在线观看视频在线高清| 人人妻人人澡人人看| 免费av毛片视频| 午夜福利影视在线免费观看| 男人操女人黄网站| 久久狼人影院| 久久 成人 亚洲| 1024视频免费在线观看| 免费看a级黄色片| 久久久久久久久免费视频了| 久久精品国产亚洲av高清一级| 欧美一区二区精品小视频在线| 日韩欧美在线二视频| 免费高清在线观看日韩| 亚洲精品在线美女| 成人亚洲精品av一区二区 | 91av网站免费观看| 999久久久精品免费观看国产| 精品日产1卡2卡| 琪琪午夜伦伦电影理论片6080| 久久国产精品影院| 欧美丝袜亚洲另类 | 在线观看日韩欧美| 80岁老熟妇乱子伦牲交| 美国免费a级毛片| 99在线视频只有这里精品首页| 亚洲精品久久成人aⅴ小说| 精品国产国语对白av| 欧美精品一区二区免费开放| 男女下面进入的视频免费午夜 | 日本撒尿小便嘘嘘汇集6| 亚洲第一青青草原| 欧美日韩瑟瑟在线播放| 美女 人体艺术 gogo| 精品国产美女av久久久久小说| 成人永久免费在线观看视频| 别揉我奶头~嗯~啊~动态视频| 亚洲avbb在线观看| 国产精品日韩av在线免费观看 | 亚洲第一青青草原| www日本在线高清视频| 美女 人体艺术 gogo| 深夜精品福利| 色老头精品视频在线观看| 国产麻豆69| 99国产精品免费福利视频| 激情在线观看视频在线高清| 日本vs欧美在线观看视频| 亚洲熟妇中文字幕五十中出 | 18禁黄网站禁片午夜丰满| 欧美老熟妇乱子伦牲交| 757午夜福利合集在线观看| 久久久精品国产亚洲av高清涩受| 深夜精品福利|