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

    供水不足地區(qū)園藝植物類溫室自動灌溉系統(tǒng)的設計與試驗

    2018-12-19 08:52:06國俊保高愛玲
    江蘇農(nóng)業(yè)科學 2018年22期
    關鍵詞:蒸騰速率冠層溫室

    國俊保, 高愛玲

    (濰坊科技學院,山東壽光 62700)

    隨著全球氣候變暖的加劇,水資源日益減少,農(nóng)業(yè)生產(chǎn)過程中須要消耗大量的水資源,因此,提高農(nóng)業(yè)生產(chǎn)過程中水資源利用率成為了當前農(nóng)業(yè)領域的研究重點[1]。在我國西北地區(qū),大量區(qū)域目前仍處于水資源匱乏狀態(tài),所以適用于水資源緊缺情況的溫室灌溉系統(tǒng)是一個有實際意義的研究內(nèi)容[2]。并且我國西南部區(qū)域時而發(fā)生旱災,如果能夠設計出水資源緊缺情況的溫室灌溉系統(tǒng),則能夠幫助旱災地區(qū)減少經(jīng)濟損失[3]。

    許多研究人員采用人工智能技術與機器學習技術對溫室的歷史數(shù)據(jù)進行學習,使用訓練獲得的模型預測溫室內(nèi)植物未來的變化趨勢[4-5]。在溫室供水充足的情況下,溫室內(nèi)微氣候的變化趨勢較為穩(wěn)定,所以人工智能與機器學習訓練的模型能夠準確地預測溫室內(nèi)植物參數(shù)的變化趨勢。但是如果發(fā)生突發(fā)性旱災,采用人工智能或者機器學習技術則難以獲得理想的預測與控制效果。

    計算流體力學能夠有效地模擬各種流場的過程[6],溫室中土壤-植物-空氣的水分交換過程是一種典型的流場過程[7],并且計算流體力學不須要通過對歷史數(shù)據(jù)的學習即可建立模型,因此計算流體力學建立的預測模型實時性較高,能夠隨著問題的實時變化而自適應地改變模型的參數(shù),從而能夠應對一些突發(fā)性狀況的發(fā)生。

    為了在不影響植物蒸騰速率的前提下降低灌溉的水量,采用計算流體力學對溫室微氣候、 植物活動參數(shù)與土壤基質參數(shù)均進行建模與預測, 該模型考慮了土壤-植物-空氣水

    分交換過程的各個環(huán)節(jié),并且能夠仿真出溫室內(nèi)蒸騰速率與微氣候的變化過程。

    1 材料與方法

    1.1 試驗環(huán)境與方法

    試驗于2017年6月在山東省濰坊地區(qū)10 m×10 m的Venlo型玻璃溫室進行。溫室由4 mm厚的園景玻璃構成,在屋脊兩側裝備了連續(xù)的屋頂通風口,通風口在室外溫度超過20 ℃的時候打開,一般通風口在白天打開,在夜晚關閉。為了防止溫室溫度過高,采用1個遮陽幕,遮陽幕的成分為鋁(16%)與聚乙烯(84%),導熱率為0.15 W/mK,透射率為50%。溫室地面為混凝土,地面放置4個架子(3 m×1.5 m),架子高度為0.8 m,水仙花花盆均勻地放置于4個架子上,花盆容積為0.74 L,高度為8.7 cm,花盆內(nèi)填充均勻分布的泥炭,密度為0.12 g/cm3。

    架子上植物的冠層面積約為18 m2,葉面積指數(shù)(leaf area index,LAI)為2.36,植物密度為15盆/m2,植物平均高度為24 cm。每天使用滴灌設備為每盆植物澆灌2次完全營養(yǎng)液,其中3個架子的植物進行供水充足的試驗,這些花盆中泥炭的水勢保持在-2 kPa以上;第4個架子的植物進行供水不足的試驗,周期性地停止灌溉,使其處于供水不足的狀態(tài)。

    1.2 溫室微氣候的數(shù)據(jù)采集

    使用傳感器采集溫室內(nèi)的微氣候數(shù)據(jù)。試驗溫室的總體結構如圖1所示,使用溫度傳感器(T)與濕度傳感器(R)測量空氣的溫度與相對濕度,傳感器分別放于3個位置:植物上方15 cm處(T1,RH1),植物內(nèi)部(T2,RH2),架子下方(T3,RH3)。此外測量葉片的溫度(TI),在地面放置PT100鉑探針,測量架子下方的地面溫度。使用CNR4型凈全輻射傳感器采集輻射的長波、短波數(shù)據(jù),在1.29 m的高度放置1個熱球式電風速計,測量冠層的氣流速度。每隔3 s采集1次上述所有參數(shù)數(shù)據(jù),每隔10 min將數(shù)據(jù)保存至數(shù)據(jù)庫中。

    1.3 基質栽培的參數(shù)

    在6個花盆中埋了負壓計與含水量傳感器,分別測量花盆的泥炭基質勢ψ(kPa)與含水量θ(體積分數(shù)),其中3個花盆為供水充足的花盆,另外3個為供水不足的花盆。每隔3 s采集1次上述所有參數(shù)數(shù)據(jù),每隔10 min將數(shù)據(jù)保存至數(shù)據(jù)庫中。

    1.4 植物活動的數(shù)據(jù)

    通過有孔石膏在6張植物葉片下面固定Cu-Cs熱電偶,每隔3 s測量1次6張葉片的溫度,將10 min總數(shù)據(jù)的平均值作為結果保存至數(shù)據(jù)庫中。使用氣孔計測量葉片的氣孔阻力。在冠層的中間放置1個高分辨率Melter-Toledo電子秤,測量冠層蒸騰所流失的水量,電子秤的托盤(56 cm×56 cm)保持水平,每隔1 h記錄花盆流水量。

    2 建立計算流體力學(CFD)模型

    為了建立計算流體力學模型,在感興趣區(qū)域周圍定義1個計算域,然后將計算域離散化為表征體元,最終求解計算域中質量(空氣、水分)與能量的守恒方程量。根據(jù)文獻分析[8],溫室的湍流一般位于溫室的高處,為了記錄湍流的運動變化,將變量分為1個平均分量與1個波動分量,并且將等式平均化處理,獲得RANS等式(雷諾平均等式)。

    使用ANSYS Fluent 15.0流體計算軟件仿真CFD模型,該軟件使用有限體積法求解二維非穩(wěn)態(tài)對流擴散方程。

    2.1 模型的基礎方程

    因為溫室中連續(xù)通風口的主風向一般垂直于通風口,則氣流的二維平面垂直于屋脊,因此模型僅須要考慮二維的情況。二維湍流動能輸運方程的一般形式為

    (1)

    式中:Φ表示非量綱輸運量的集合,即動量、質量與能量;U與V是速度向量的元素;Γ是擴散系數(shù);SΦ是輸運方程的源項(source term)。將k-ε湍流模型設為多組分湍流封閉模型,因為空氣密度(ρ)不僅依賴溫度還依賴空氣的含水量,所以無法應用Boussinesq模型。將空氣密度定義為一個關于溫度與質量的函數(shù):

    (2)

    式中:T是局部溫度,K;R是理想氣體常數(shù),R=8.31 J/(mol·K);Pop為工作壓力,Pa;ωi是組分i的質量值;Mi為組分i的分子量,kg/mol。

    2.2 輻射子模型

    為了區(qū)分輻射中短波(0.1~3 μm)與長波(3~100 μm)的貢獻度,設計了一個輻射子模型。對于有限數(shù)量的離散立體角,可用式(3)求解單色亮度方程,每個離散立體角對應的向量方向為

    (3)

    聚集整個空間的每個單色亮度,即可獲得全局的輻射通量qs,表示為

    (4)

    熱平衡中每個單位體積的輻射貢獻度是能量守恒方程的一個源項(source term),表示了輻射qs(s)的散度。

    (5)

    在當前的研究中,短波長與長波長之間存在差異[9]。與短波輻射相比,冠層吸收長波輻射的量可忽略不計,在白天冠層吸收的長波輻射小于短波輻射的4%。

    2.3 植物子模型

    2.3.1 多孔介質 將植物抽象為均勻分布的多孔介質,該多孔介質對氣流具有阻力(對應植物的氣孔阻力),表示為能量守恒方程中的一個源項SΦ。根據(jù)Darcy-Forchheimer方程,植物每個單位體積的拉力可表示為黏性丟失項與內(nèi)部丟失項2個項的總和,如下式所示。

    (6)

    式中:K是多孔介質的滲透性,%;CF是非線性動量的衰減系數(shù);ρa是空氣密度,kg/m3;μ是空氣的動態(tài)黏度,(N·s)/m2;V是空氣的速度值,m/s。根據(jù)文獻[10]分析,實際的空氣動態(tài)黏度非常低,是10-5(N·s)/m2數(shù)量級,所以可以忽略式(6)的第1項,可得:

    (7)

    式中:LAD是葉密度,定義為LAI/冠層高度(H);CD是拉力系數(shù)。

    從式(7)可以看出,僅須要獲得CD值即可求解植物每個單位體積的拉力,鳳仙花的CD值約為0.32。

    2.3.2 能量平衡模型 冠層的太陽輻射為Rg,植物通過蒸騰作用以及與周圍空氣的熱交換進行能量傳遞,假設蒸騰作用的熱通量密度為Trd,周圍空氣的熱通量密度為Qs。最終,葉片的能量平衡公式寫為

    Rgabs-Trd-Qs=0。

    (8)

    可根據(jù)Beer法則[式(9)]直接推導出冠層每個細胞吸收的輻射量Rgabs,W/m2:

    Rg(y)=Rg0exp[-KcLAD(H-y)];

    (9)

    Rgabs=Rg(yi)-Rg(yi+1)。

    (10)

    式中:yi與yi+1分別是細胞i的上一坐標與下一坐標。Rg0是冠層上的全局輻射,W/m2;Kc是輻射的消光系數(shù);供水充足的情況下該值為0.95,供水不足的情況下該值為0.64,這2個值是使用式(9)計算冠層上面與下面的短波輻射而來。這個差值的原因是在植物缺水的情況下,植物的葉片會自動地旋轉,這個生理機制能夠降低植物吸收的太陽輻射量,因此降低了蒸騰作用。

    冠層的可感知熱量定義為

    Qs=2LADρaCp[(TI-Ta)/Ra]。

    (11)

    式中:ρa是空氣密度,kg/m3;Cp是指定壓強對應的空氣熱量,J/(kg·K);TI與Ta分別是葉片與空氣的溫度。

    一個給定區(qū)域的潛熱密度定義為

    (12)

    式中:γ是溫濕常數(shù);VPDa是空氣的飽和水汽壓差,Pa;Δ是飽和水蒸氣壓曲線的斜率;根據(jù)文獻[11],將氣動阻力Ra設為一個常量,Ra=271 s/m。Rs是氣孔阻力,s/m,采用FFD(全因子試驗設計)[12]方法可推導出Rs,表示為

    Rs=(-115×rg-139×rh-39×τ+139×rg×rh+43×rg×τ+11×rh×τ+661×rg2-368×rg3)×[1+(ψ-11.41)1.05]。

    (13)

    式中:rg=(Rg-75)/75,rh=(RH-65)/10,τ=(T-20.5)/5.5,rg是輻射量,rh是濕度值,τ是溫度因子,ψw是土壤基質勢(kPa)的減少量,RH為相對濕度。將式(8)進行轉換,可計算出細胞的葉片溫度(TI):

    (14)

    2.3.3 水平衡模型 根據(jù)文獻[13]模型,可通過下式計算基質水勢的減少量ψw:

    (15)

    式中:θres與θsat分別是土壤剩余含水量與土壤飽和含水量;α(kPa)與N是2個待校準的參數(shù),這2個參數(shù)來自于拉力計與含水量傳感器的測量數(shù)據(jù)。根據(jù)多線性擬合的結果,可獲得以下參數(shù)值:θsat=0.887 cm3/cm3,θres=0.1 cm3/cm3,α=0.134 kPa-1,N=1.469。通過式(16)計算每個時間點i的含水量,式(16)定義了蒸騰作用與初始化含水量(θ0)的關系。

    (16)

    式中:θi是當前時間點的含水量i,cm3/cm3;θ0是初始含水量,cm3/cm3;ts是時段長度,本研究ts等于3 600 s;ss是架子的總面積;Travg,j是在時間點j時的冠層平均蒸騰速率;λ是水的潛熱量,J/kg;m是植物總量(m=264 kg);ρw是水密度,kg/m3;Vp是花盆的容積,m3。

    2.3.4 模擬程序 對于動量與熱量的輸運方程采用二階迎風差分格式[14]處理,提高輸運方程的計算準確率。對于壓力連接等式采用半隱式格式[15]來求解壓力-動量方程。所有變量的收斂條件是10-6數(shù)量級。

    3 模型驗證

    采用2個準確率性能指標評估本模型的預測準確率,分別為RMSE(均方根)與R2(確定系數(shù)),RMSE也稱為預測系統(tǒng)的擬合標準差,值越小表示預測越精確;R2是通過數(shù)據(jù)的變化來表征一個擬合的好壞,R2取值范圍為[0,1],越接近1,擬合效果越好。

    3.1 試驗結果與分析

    首先對供水充足與供水不足2種情況分別進行分析與研究。將二維非穩(wěn)態(tài)仿真的時間步長設為2 h,周期為從6月14日01:00—23:00,從11:00—20:00打開通風口,6月14日對供水不足的花盆停止灌溉。試驗中測量的數(shù)據(jù)作為預測模型的輸入數(shù)據(jù),測量的結果如圖2所示。從圖2可以看出,太陽輻射的變化曲線與氣溫、地面溫度均具有相關性,在13:00達到最大值。

    在6月14日01:00使用6月14日00:00的數(shù)據(jù)進行一次初步仿真,對于供水不足的情況,初步仿真結果中的初始水量為0.652 cm3/cm3,使用基于水平衡的子模型計算每個時間點的土壤基質水勢。對于供水充足的情況,將土壤基質水勢設為恒定值,該值固定為-1 kPa。

    3.2 空氣溫度與濕度

    2種供水條件下3個測量位置的氣溫隨著時間的變化曲線如圖3與圖4所示。氣溫與葉片溫度的測量值與預測值表現(xiàn)出相同的趨勢,2種供水條件下氣溫的預測值與測量值能夠較好地擬合。由表1可以看出,對于供水充足情況,所有的結果中R2結果均大于或等于0.97,RMSE均小于或等于1.14;對于供水不足情況,所有的結果中R2結果均大于或等于0.94,RMSE均小于或等于 1.62。因為2種供水條件的溫度差異較小,所以供水是否充足對觀察周圍溫度的影響較小。2種供水條件的差異主要由初始的環(huán)境參數(shù)引起,供水不足情況的葉片溫度比供水充足情況約高1K,其原因是供水不足情況的蒸騰作用減弱,而蒸騰作用能夠降低葉片的溫度。2種供水條件的葉片底部溫度幾乎相同,因為葉片底部的蒸騰作用均較弱。

    此外,研究了2種供水條件下3個測量位置的相對濕度,分別表示為RH1、RH2與RH3,相對濕度隨著時間的變化曲線如圖5所示。對于供水充足的情況,結果顯示冠層內(nèi)部相對濕度的預測值與測量值高于冠層的上面與下面。對于供水不足的情況,3個位置相對濕度的測量值與預測值較為接近??傮w而言,2種供水條件的預測值與測量值實現(xiàn)了較好的擬合。從表1可以看出,對于供水充足的情況,R2大于0.88,RMSE小于8.66;對于供水不足的情況,R2大于0.78,RMSE小于6.69。

    3.3 泥炭土的基質勢

    在停止灌溉5 d后,統(tǒng)計了1 d內(nèi)測量與預測的泥炭土的基質勢值(圖6)。因為白天植物的視覺特征較為明顯,因此試驗主要關注白天。因為停止了灌溉,土壤內(nèi)的水量隨著植物的蒸騰作用而減少。從圖6可以看出,泥炭土基質勢的預測值較好地擬合了測量值,并且獲得了較高的預測準確率,預測性能指標分別為r2=0.99,RMSE=-3.38。

    表1 溫室內(nèi)微氣候各個參數(shù)的測量結果

    3.4 觀賞植物的相關參數(shù)

    3.4.1 植物的氣孔阻力 為了分析光照對植物氣孔阻力的影響,對陰影區(qū)域葉片與光照區(qū)域葉片均進行了試驗分析,供水充足與供水不足2種情況的氣孔阻力變化曲線如圖7所示。從圖7-b可以看出,從大約10:00開始,氣孔阻力開始上升,因為此時土壤基質勢開始降低,植物通過增加氣孔阻力來限制植物的蒸騰作用。從圖7-a可以看出,供水充足情況的氣孔阻力也從10:00開始上升,此時陽光輻射與空氣氣壓差均開始上升,所以植物通過增加氣孔阻力來限制植物的蒸騰作用,但其上升幅度遠小于供水不足的情況。

    CFD模型能夠有效地預測出陰影區(qū)域葉片與光照區(qū)域葉片的氣孔阻力變化曲線,供水充足情況的預測性能指標為r2>0.4,RMSE<168;供水不足情況的預測性能指標為r2>0.51,RMSE<527。

    3.4.2 植物的蒸騰速率 供水充足與供水不足2種情況的植物蒸騰速率變化曲線如圖8所示。供水不足情況的蒸騰速率小于供水充足情況。供水不足情況的蒸騰速率最小測量值為1.43 g/(m2·h),低于供水充足情況的最小蒸騰速率??傮w而言,本模型能夠有效地預測供水充足與供水不足2種情況的植物蒸騰速率,供水充足情況的預測性能指標為r2=0.98,RMSE=5.87;供水不足情況的預測性能指標為r2=0.97,RMSE=5.64。此外,計算了24 h蒸騰速率的總預測誤差,供水充足情況的總預測誤差為3.32%,供水不足情況的總預測誤差為4.82%。

    4 結論

    目前溫室自動灌溉系統(tǒng)大多針對供水充足的穩(wěn)定環(huán)境,為了提高供水不足情況的灌溉效率,設計了一種采用物聯(lián)網(wǎng)的園藝植物類溫室自動灌溉系統(tǒng)。為了在不影響植物蒸騰速率的前提下降低灌溉的水量,采用計算流體力學對溫室微氣候、植物活動參數(shù)與土壤環(huán)境參數(shù)均進行建模與預測,該模型考慮了土壤-植物-空氣水分交換過程的各個環(huán)節(jié),并且能夠仿真出溫室內(nèi)蒸騰速率與微氣候的變化過程。在山東濰坊地區(qū)采用水仙花盆景進行了試驗,對供水充足與供水不足2種情況均進行了試驗,該模型對溫室微氣候、植物活動參數(shù)與泥炭土基質勢均獲得了準確的預測效果。

    猜你喜歡
    蒸騰速率冠層溫室
    基于低空遙感的果樹冠層信息提取方法研究
    現(xiàn)代溫室羊肚菌栽培技術
    基于激光雷達的樹形靶標冠層葉面積探測模型研究
    蒼松溫室 蒼松灌溉
    蒼松溫室 蒼松灌溉
    土壤水分狀況對獼猴桃葉片蒸騰速率影響研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    大氣科學(2021年1期)2021-04-16 07:34:18
    黃土丘陵區(qū)山杏人工林蒸騰速率與環(huán)境因子的關系
    可以避免一個溫室化的地球嗎?
    英語文摘(2019年2期)2019-03-30 01:48:28
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    超色免费av| 国产午夜精品久久久久久| 在线观看免费视频日本深夜| 又黄又爽又免费观看的视频| 免费在线观看视频国产中文字幕亚洲| xxx96com| 亚洲情色 制服丝袜| 成人亚洲精品av一区二区 | 亚洲欧美日韩另类电影网站| 老熟妇乱子伦视频在线观看| 国产精品永久免费网站| 高清在线国产一区| 精品第一国产精品| 人妻久久中文字幕网| 亚洲欧美精品综合久久99| 国产麻豆69| 午夜老司机福利片| 日本撒尿小便嘘嘘汇集6| 久9热在线精品视频| 一级作爱视频免费观看| 91国产中文字幕| 香蕉国产在线看| 国产欧美日韩一区二区精品| 中文字幕人妻熟女乱码| 99热国产这里只有精品6| 香蕉国产在线看| 99在线视频只有这里精品首页| 母亲3免费完整高清在线观看| а√天堂www在线а√下载| 99国产精品一区二区蜜桃av| 午夜亚洲福利在线播放| 黄频高清免费视频| 大型av网站在线播放| 免费搜索国产男女视频| 欧美黑人欧美精品刺激| 一a级毛片在线观看| 久久久久亚洲av毛片大全| 成年版毛片免费区| 欧美乱码精品一区二区三区| 99热国产这里只有精品6| 免费人成视频x8x8入口观看| 国产成人精品久久二区二区91| 交换朋友夫妻互换小说| a级毛片在线看网站| 国产成人影院久久av| 久9热在线精品视频| 后天国语完整版免费观看| 成人手机av| 99精品久久久久人妻精品| 桃色一区二区三区在线观看| 色婷婷av一区二区三区视频| av片东京热男人的天堂| 欧美丝袜亚洲另类 | 制服人妻中文乱码| 久久人妻av系列| 可以在线观看毛片的网站| 黑人操中国人逼视频| 久久精品亚洲精品国产色婷小说| 中文字幕高清在线视频| 美国免费a级毛片| 国产精品乱码一区二三区的特点 | 国产精品影院久久| 日韩免费av在线播放| 国产亚洲精品第一综合不卡| 一级a爱片免费观看的视频| 99re在线观看精品视频| 另类亚洲欧美激情| 国产高清激情床上av| 一区二区日韩欧美中文字幕| 午夜福利,免费看| 午夜福利免费观看在线| 免费观看人在逋| 高清黄色对白视频在线免费看| 欧美成人性av电影在线观看| 电影成人av| 黄片大片在线免费观看| 国产熟女午夜一区二区三区| 精品无人区乱码1区二区| 成年女人毛片免费观看观看9| 国产精品久久电影中文字幕| 久久午夜综合久久蜜桃| xxx96com| 91大片在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美日韩高清在线视频| 亚洲精品av麻豆狂野| 国产精品野战在线观看 | 国产成年人精品一区二区 | 亚洲欧洲精品一区二区精品久久久| 淫妇啪啪啪对白视频| 免费在线观看日本一区| 91九色精品人成在线观看| 老司机亚洲免费影院| 精品久久久精品久久久| 黑人巨大精品欧美一区二区mp4| 999久久久精品免费观看国产| 三上悠亚av全集在线观看| 丁香欧美五月| 日韩欧美一区视频在线观看| 国产精品国产av在线观看| 免费高清视频大片| 国产午夜精品久久久久久| av在线天堂中文字幕 | 中出人妻视频一区二区| 免费av中文字幕在线| av天堂久久9| 国产在线精品亚洲第一网站| av国产精品久久久久影院| 国产精品久久久久成人av| 日韩精品中文字幕看吧| 夜夜看夜夜爽夜夜摸 | 久久久国产欧美日韩av| 久久国产精品男人的天堂亚洲| 免费av毛片视频| 十八禁网站免费在线| 久久久精品欧美日韩精品| 一个人观看的视频www高清免费观看 | 欧美+亚洲+日韩+国产| 亚洲精品一区av在线观看| 在线观看免费视频网站a站| 丰满的人妻完整版| 久久天堂一区二区三区四区| 欧美日本中文国产一区发布| 午夜福利免费观看在线| 免费女性裸体啪啪无遮挡网站| 国产精品偷伦视频观看了| 中文字幕人妻熟女乱码| 最近最新免费中文字幕在线| 在线观看午夜福利视频| 中文字幕人妻熟女乱码| 人人妻人人澡人人看| 精品久久久久久电影网| 97超级碰碰碰精品色视频在线观看| 精品午夜福利视频在线观看一区| 可以免费在线观看a视频的电影网站| 精品午夜福利视频在线观看一区| 久久热在线av| 在线观看舔阴道视频| 国产精品98久久久久久宅男小说| 身体一侧抽搐| av中文乱码字幕在线| 国产不卡一卡二| 日韩欧美一区二区三区在线观看| 欧美激情 高清一区二区三区| 美女高潮到喷水免费观看| 久久人妻熟女aⅴ| 香蕉国产在线看| av片东京热男人的天堂| 色尼玛亚洲综合影院| 夜夜躁狠狠躁天天躁| 无限看片的www在线观看| 丰满饥渴人妻一区二区三| 纯流量卡能插随身wifi吗| 精品一品国产午夜福利视频| 十八禁网站免费在线| 久久人妻熟女aⅴ| a在线观看视频网站| 久久久精品国产亚洲av高清涩受| 欧美日韩亚洲国产一区二区在线观看| 最近最新中文字幕大全电影3 | 精品福利观看| 麻豆久久精品国产亚洲av | 女警被强在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲 欧美一区二区三区| 欧美日韩av久久| 麻豆久久精品国产亚洲av | 精品欧美一区二区三区在线| xxx96com| 亚洲va日本ⅴa欧美va伊人久久| 在线观看www视频免费| 久久久国产欧美日韩av| 亚洲熟妇中文字幕五十中出 | 麻豆一二三区av精品| 一二三四在线观看免费中文在| 日韩欧美国产一区二区入口| 亚洲视频免费观看视频| 久久中文字幕一级| 黄片大片在线免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 我的亚洲天堂| 激情在线观看视频在线高清| 久久久久久人人人人人| 日韩 欧美 亚洲 中文字幕| 91精品三级在线观看| 亚洲黑人精品在线| 国产区一区二久久| av天堂久久9| 欧美性长视频在线观看| 国产精品综合久久久久久久免费 | 搡老熟女国产l中国老女人| 国产亚洲av高清不卡| 丁香欧美五月| 中文字幕人妻丝袜一区二区| 久久精品人人爽人人爽视色| 成人免费观看视频高清| 又大又爽又粗| 婷婷丁香在线五月| 色综合婷婷激情| 色哟哟哟哟哟哟| 久久精品aⅴ一区二区三区四区| 亚洲狠狠婷婷综合久久图片| 桃红色精品国产亚洲av| 亚洲第一欧美日韩一区二区三区| 最新美女视频免费是黄的| 一级毛片高清免费大全| 午夜激情av网站| 亚洲第一av免费看| 国产精品爽爽va在线观看网站 | 久久国产亚洲av麻豆专区| 一边摸一边做爽爽视频免费| 亚洲精品一区av在线观看| 黄色女人牲交| 久久99一区二区三区| 极品人妻少妇av视频| 两个人看的免费小视频| 激情在线观看视频在线高清| 国产精品综合久久久久久久免费 | 成年版毛片免费区| 新久久久久国产一级毛片| 亚洲成人国产一区在线观看| 亚洲国产精品一区二区三区在线| 精品一区二区三区四区五区乱码| 91老司机精品| 亚洲五月色婷婷综合| av视频免费观看在线观看| 亚洲av熟女| 99久久精品国产亚洲精品| xxx96com| 亚洲三区欧美一区| 男女下面进入的视频免费午夜 | 黄色丝袜av网址大全| 热99国产精品久久久久久7| 国产不卡一卡二| 天堂动漫精品| 午夜亚洲福利在线播放| 国产精品秋霞免费鲁丝片| 一级毛片女人18水好多| 夜夜看夜夜爽夜夜摸 | 大陆偷拍与自拍| 久久久久九九精品影院| 午夜久久久在线观看| 亚洲午夜理论影院| 午夜视频精品福利| 亚洲美女黄片视频| av在线播放免费不卡| 久久亚洲精品不卡| 日韩欧美一区视频在线观看| 国产精品久久久av美女十八| 亚洲美女黄片视频| 99久久人妻综合| 亚洲伊人色综图| 欧美+亚洲+日韩+国产| 亚洲专区中文字幕在线| 亚洲精品在线美女| 女同久久另类99精品国产91| 女同久久另类99精品国产91| 国产精品综合久久久久久久免费 | 久久午夜亚洲精品久久| 亚洲人成电影观看| 99国产精品一区二区蜜桃av| 精品久久久久久久毛片微露脸| 亚洲av熟女| 大码成人一级视频| 久久久久久久久中文| 亚洲av成人一区二区三| 18禁黄网站禁片午夜丰满| 精品国产美女av久久久久小说| 欧美不卡视频在线免费观看 | 亚洲第一av免费看| 亚洲片人在线观看| 大型av网站在线播放| 久久精品国产99精品国产亚洲性色 | 免费高清在线观看日韩| 亚洲一区高清亚洲精品| 国产免费男女视频| 一级a爱片免费观看的视频| 免费高清在线观看日韩| 国产麻豆69| 亚洲午夜理论影院| 大型黄色视频在线免费观看| 欧美黑人欧美精品刺激| 国产精品成人在线| 黄色a级毛片大全视频| 悠悠久久av| 亚洲精品美女久久av网站| 操出白浆在线播放| 窝窝影院91人妻| 亚洲少妇的诱惑av| 久久中文字幕一级| e午夜精品久久久久久久| 亚洲专区字幕在线| 无遮挡黄片免费观看| 久久中文字幕一级| 亚洲欧美一区二区三区久久| 伊人久久大香线蕉亚洲五| 成人亚洲精品一区在线观看| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 最好的美女福利视频网| 国内毛片毛片毛片毛片毛片| 性欧美人与动物交配| 久久精品亚洲熟妇少妇任你| 国产成人系列免费观看| 亚洲七黄色美女视频| 中文字幕色久视频| av电影中文网址| 免费少妇av软件| 欧美成人性av电影在线观看| 精品欧美一区二区三区在线| 亚洲国产精品合色在线| 精品国产一区二区三区四区第35| 亚洲aⅴ乱码一区二区在线播放 | √禁漫天堂资源中文www| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 国产成人影院久久av| 我的亚洲天堂| 日韩免费av在线播放| 日韩欧美免费精品| 成人影院久久| 久久婷婷成人综合色麻豆| 亚洲av日韩精品久久久久久密| 久99久视频精品免费| 男男h啪啪无遮挡| 99久久人妻综合| 黄色女人牲交| 五月开心婷婷网| 一级片免费观看大全| 久久九九热精品免费| 夜夜夜夜夜久久久久| 免费久久久久久久精品成人欧美视频| 黄片播放在线免费| 国产伦一二天堂av在线观看| 国产高清视频在线播放一区| 国产免费av片在线观看野外av| 亚洲av成人不卡在线观看播放网| 天天影视国产精品| 国产成年人精品一区二区 | 最近最新中文字幕大全电影3 | 9191精品国产免费久久| 国产亚洲av高清不卡| 国产91精品成人一区二区三区| 别揉我奶头~嗯~啊~动态视频| 韩国av一区二区三区四区| 精品国产一区二区三区四区第35| 母亲3免费完整高清在线观看| 9色porny在线观看| 日本免费一区二区三区高清不卡 | 久久亚洲真实| 高潮久久久久久久久久久不卡| www.自偷自拍.com| 欧美激情高清一区二区三区| 国产精品 欧美亚洲| 日本黄色视频三级网站网址| 国产一区二区三区在线臀色熟女 | 亚洲国产精品999在线| 亚洲国产欧美一区二区综合| 国产高清国产精品国产三级| 国产成人欧美| 大陆偷拍与自拍| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 嫩草影院精品99| 熟女少妇亚洲综合色aaa.| 精品国产亚洲在线| 国产精品爽爽va在线观看网站 | 中文字幕人妻丝袜制服| 美女扒开内裤让男人捅视频| 亚洲欧美激情在线| 久久香蕉国产精品| 精品久久久久久成人av| 免费在线观看影片大全网站| 欧美日韩乱码在线| 亚洲欧洲精品一区二区精品久久久| 日本黄色视频三级网站网址| 亚洲人成77777在线视频| 老司机福利观看| 国产一区二区三区视频了| 日韩视频一区二区在线观看| 伦理电影免费视频| 好男人电影高清在线观看| 老鸭窝网址在线观看| 在线观看免费高清a一片| 亚洲第一av免费看| 麻豆国产av国片精品| 亚洲一区二区三区不卡视频| 午夜免费观看网址| 免费在线观看亚洲国产| 女同久久另类99精品国产91| 国产精品自产拍在线观看55亚洲| svipshipincom国产片| 日本免费一区二区三区高清不卡 | 交换朋友夫妻互换小说| 少妇被粗大的猛进出69影院| 后天国语完整版免费观看| 老司机深夜福利视频在线观看| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | 亚洲av第一区精品v没综合| 丰满迷人的少妇在线观看| 日韩精品青青久久久久久| 国产精品久久视频播放| 美国免费a级毛片| 88av欧美| 国产成人欧美在线观看| 一区二区三区激情视频| 丝袜美腿诱惑在线| 色尼玛亚洲综合影院| 免费看a级黄色片| 嫁个100分男人电影在线观看| 美国免费a级毛片| 男女床上黄色一级片免费看| 久久精品国产99精品国产亚洲性色 | 又大又爽又粗| 精品久久久精品久久久| 免费在线观看完整版高清| 大码成人一级视频| 18禁观看日本| 亚洲五月天丁香| 91在线观看av| 岛国视频午夜一区免费看| 亚洲第一欧美日韩一区二区三区| 日本黄色视频三级网站网址| 免费女性裸体啪啪无遮挡网站| 一级毛片女人18水好多| 日本 av在线| 亚洲欧美激情综合另类| 亚洲精品成人av观看孕妇| 日韩欧美免费精品| 久久青草综合色| 日韩精品免费视频一区二区三区| 免费在线观看黄色视频的| 侵犯人妻中文字幕一二三四区| 国产精品九九99| 级片在线观看| 午夜福利欧美成人| 国产成人一区二区三区免费视频网站| 国产精品亚洲一级av第二区| 久久中文字幕一级| xxxhd国产人妻xxx| 男女下面进入的视频免费午夜 | 在线观看66精品国产| 欧美日韩亚洲综合一区二区三区_| 亚洲五月婷婷丁香| 俄罗斯特黄特色一大片| 国产aⅴ精品一区二区三区波| 色在线成人网| 午夜免费激情av| 国产免费男女视频| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 久久精品91无色码中文字幕| 午夜日韩欧美国产| 宅男免费午夜| 99久久综合精品五月天人人| 12—13女人毛片做爰片一| 日韩人妻精品一区2区三区| 在线观看免费高清a一片| 国产精品国产av在线观看| 国产亚洲欧美98| 老汉色∧v一级毛片| 久久久久久久久久久久大奶| 丁香六月欧美| 一进一出好大好爽视频| 日本撒尿小便嘘嘘汇集6| 妹子高潮喷水视频| 交换朋友夫妻互换小说| 香蕉久久夜色| 成年人免费黄色播放视频| 亚洲男人的天堂狠狠| av天堂在线播放| 别揉我奶头~嗯~啊~动态视频| 国产乱人伦免费视频| 精品国产乱码久久久久久男人| 精品国内亚洲2022精品成人| 久久精品aⅴ一区二区三区四区| 色婷婷av一区二区三区视频| 国产99白浆流出| 99国产精品一区二区蜜桃av| 女人被狂操c到高潮| 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 怎么达到女性高潮| 日韩欧美一区视频在线观看| 美女 人体艺术 gogo| 999久久久精品免费观看国产| 国产成人精品久久二区二区91| 亚洲精品av麻豆狂野| 99热只有精品国产| 国产精品永久免费网站| 美女大奶头视频| 日本wwww免费看| 久久久久久免费高清国产稀缺| 一级片免费观看大全| 一本大道久久a久久精品| 亚洲av日韩精品久久久久久密| 我的亚洲天堂| 91麻豆精品激情在线观看国产 | 自线自在国产av| 欧美日韩视频精品一区| 欧美成人午夜精品| 老鸭窝网址在线观看| 女性被躁到高潮视频| 久久人妻熟女aⅴ| 国产熟女xx| 国产一区二区三区在线臀色熟女 | 国产91精品成人一区二区三区| 夫妻午夜视频| 精品国产一区二区三区四区第35| 性少妇av在线| 成人国产一区最新在线观看| 欧美色视频一区免费| 动漫黄色视频在线观看| 久久精品91无色码中文字幕| 精品福利观看| 一个人观看的视频www高清免费观看 | 国产精品综合久久久久久久免费 | 国产一区二区在线av高清观看| 黄片播放在线免费| 欧美激情 高清一区二区三区| 国产人伦9x9x在线观看| 国产av一区二区精品久久| 午夜久久久在线观看| 欧美黑人精品巨大| 欧美另类亚洲清纯唯美| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 成人三级黄色视频| 久久久久久亚洲精品国产蜜桃av| 日韩欧美国产一区二区入口| 久久这里只有精品19| 自线自在国产av| 一进一出抽搐动态| 国产深夜福利视频在线观看| 无人区码免费观看不卡| 久久久久亚洲av毛片大全| 国产成人av教育| 日韩精品青青久久久久久| 69精品国产乱码久久久| 悠悠久久av| 国产极品粉嫩免费观看在线| 成人av一区二区三区在线看| 久久久国产成人免费| 国产精品综合久久久久久久免费 | 亚洲九九香蕉| 老司机深夜福利视频在线观看| 69av精品久久久久久| 97人妻天天添夜夜摸| 中亚洲国语对白在线视频| 国产精品1区2区在线观看.| 99riav亚洲国产免费| 在线观看www视频免费| 久久九九热精品免费| 悠悠久久av| 三级毛片av免费| 黑人巨大精品欧美一区二区蜜桃| 精品久久久久久电影网| 亚洲自拍偷在线| 一级片'在线观看视频| av欧美777| 欧美另类亚洲清纯唯美| 国产高清激情床上av| 中文字幕最新亚洲高清| 岛国视频午夜一区免费看| xxxhd国产人妻xxx| 国产一区二区激情短视频| 色播在线永久视频| 男女之事视频高清在线观看| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 亚洲欧美日韩另类电影网站| 国产精华一区二区三区| 热re99久久精品国产66热6| 青草久久国产| 午夜视频精品福利| 无人区码免费观看不卡| 91字幕亚洲| 国产单亲对白刺激| 亚洲片人在线观看| 午夜免费鲁丝| 精品久久久久久久毛片微露脸| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 亚洲三区欧美一区| 午夜精品久久久久久毛片777| 国产欧美日韩一区二区三区在线| 精品一区二区三区四区五区乱码| 久久久国产欧美日韩av| 成人三级做爰电影| 日韩欧美一区视频在线观看| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 免费在线观看日本一区| 国产一区二区激情短视频| 女人爽到高潮嗷嗷叫在线视频| 在线看a的网站| 亚洲欧美精品综合一区二区三区| 国产精华一区二区三区| 欧美一区二区精品小视频在线| 三级毛片av免费| 久久人人精品亚洲av| 大码成人一级视频| a级片在线免费高清观看视频| 免费观看精品视频网站| 成人免费观看视频高清| 国产熟女xx| 99精国产麻豆久久婷婷| 丁香六月欧美| 19禁男女啪啪无遮挡网站| 国产亚洲精品一区二区www| 精品福利永久在线观看| 亚洲五月色婷婷综合| 中文字幕人妻熟女乱码|