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

    農(nóng)業(yè)旱情評(píng)估方法研究
    ——基于農(nóng)田水分循環(huán)模擬與遙感影像信息同化

    2016-03-23 03:53:52倪深海閆娜娜水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室南京210029南京水利科學(xué)研究院南京210029中國科學(xué)院遙感與數(shù)字地球研究所北京10009山西省水文水資源勘測(cè)局太原00001
    中國農(nóng)村水利水電 2016年1期
    關(guān)鍵詞:旱情農(nóng)田作物

    倪深海,顧 穎,閆娜娜,常 勝,趙 凱,申 瑜(1.水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,南京 210029;2.南京水利科學(xué)研究院,南京 210029;.中國科學(xué)院遙感與數(shù)字地球研究所,北京 10009;.山西省水文水資源勘測(cè)局,太原 00001)

    0 引 言

    農(nóng)業(yè)旱情是以土壤中實(shí)際貯存的、可供作物利用的水量多少為依據(jù)判別作物的受旱情況。無論是降水還是灌溉,水分都要首先貯存在土壤中,然后被作物生長逐步吸收利用;另外,區(qū)域的氣象、水文、灌溉等條件以及灌溉用水管理狀況,也都能夠很好地通過農(nóng)業(yè)旱情信息得以體現(xiàn)[1-3]。因此,可通過土壤墑情信息、農(nóng)作物缺水信息及遙感影像識(shí)別農(nóng)業(yè)旱情。如何判別農(nóng)業(yè)旱情是旱情預(yù)測(cè)預(yù)警的關(guān)鍵,為表征農(nóng)業(yè)旱情的發(fā)生與發(fā)展過程,可采用多種判別方法。采用仿真技術(shù)模擬作物生長過程,識(shí)別農(nóng)作物需水、缺水信息,并利用遙感信息進(jìn)行同化是判別農(nóng)業(yè)旱情行之有效的方法[4,5],同化后的模型模擬農(nóng)業(yè)旱情與實(shí)測(cè)結(jié)果較好的擬合度,實(shí)現(xiàn)對(duì)旱情及趨勢(shì)的全面監(jiān)視和評(píng)估的目的。

    信息同化技術(shù)應(yīng)用到旱情監(jiān)測(cè)和評(píng)估中,可以克服單一信息評(píng)估旱情的片面性。研究表明,基于農(nóng)田水分循環(huán)模擬與遙感信息同化的農(nóng)業(yè)旱情評(píng)估模型是可行的,為及時(shí)、準(zhǔn)確識(shí)別評(píng)估農(nóng)業(yè)旱情提供一種新的方法,對(duì)我國旱情監(jiān)測(cè)評(píng)估方面有著重要應(yīng)用和推廣價(jià)值。

    1 信息同化技術(shù)與方法

    1.1 信息同化含義

    信息同化是將各種不同來源,不同時(shí)空、不同觀測(cè)手段獲得的數(shù)據(jù)和資料與數(shù)學(xué)模型有機(jī)結(jié)合,不斷更新系統(tǒng)狀態(tài)與參數(shù),建立數(shù)據(jù)與模型相互協(xié)調(diào)的優(yōu)化關(guān)系,提高物理過程模擬或預(yù)報(bào)精度的技術(shù)。

    早在20世紀(jì)50年代信息同化就被成功應(yīng)用于數(shù)值天氣預(yù)報(bào),之后在海洋預(yù)測(cè)系統(tǒng)中也得到廣泛應(yīng)用,直到20世紀(jì)90年代才被用來研究陸面過程. 陸面信息同化是在大氣和海洋數(shù)據(jù)同化研究的基礎(chǔ)上迅速成長起來的一個(gè)嶄新的領(lǐng)域,其核心思想是在陸面過程模型的動(dòng)力框架內(nèi),融合不同來源和不同分辨率的直接與間接觀測(cè),將陸面過程模型和各種觀測(cè)算子集成為不斷地依靠觀測(cè)而自動(dòng)調(diào)整型軌跡,并且減小誤差的預(yù)報(bào)系統(tǒng)[6,7]。旱情信息的同化正是陸面過程同化的一個(gè)重要內(nèi)容,利用遙感監(jiān)測(cè)旱情信息與農(nóng)田水分循環(huán)模擬數(shù)據(jù)進(jìn)行同化,以不斷更新農(nóng)田水分循環(huán)模擬模型狀態(tài)變量與參數(shù),提高農(nóng)田水分循環(huán)過程模擬與預(yù)報(bào)精度。

    1.2 信息同化算法

    信息同化算法是連接遙感影像信息與模型模擬旱情信息的關(guān)鍵核心部分[8],本文采用逐步訂正法進(jìn)行信息同化。逐步訂正法采用模擬旱情信息作為初估場(chǎng),不斷插入遙感監(jiān)測(cè)旱情信息,使用這些遙感監(jiān)測(cè)值以及給定的影響半徑對(duì)初估場(chǎng)進(jìn)行一次訂正,然后使用訂正后的分析場(chǎng)作為下一次訂正的初估場(chǎng),同時(shí)縮小影響半徑,進(jìn)行下一次訂正。這樣的循環(huán)過程就構(gòu)成四維數(shù)據(jù)同化。由于模擬旱情信息和遙感監(jiān)測(cè)旱情信息相近,并且分析場(chǎng)空間連續(xù)性較好,同時(shí)又考慮了時(shí)間的連續(xù)性,因此,在這樣的基礎(chǔ)上進(jìn)行連續(xù)訂正,并在訂正過程中縮小影響半徑,首先可以去掉初估場(chǎng)的大尺度誤差,然后使模擬的旱情信息分析場(chǎng)越來越逼近實(shí)際旱情信息場(chǎng),就比直接插值容易得到較好的效果。

    (1)

    式中:fbi是背景場(chǎng)在第i個(gè)計(jì)算單元上的值;f0i是計(jì)算單元i上零次迭代估計(jì)值(f0i=fbi);fni是格點(diǎn)i上第n次迭代估計(jì)值;f0k是計(jì)算單元i周圍的第k個(gè)觀測(cè);ε2是觀測(cè)誤差方差與背景場(chǎng)誤差方差的比率的估計(jì);Kni是距離格點(diǎn)i為Rn內(nèi)的觀測(cè)值總數(shù);wnik為權(quán)重。

    逐步迭代法的權(quán)重:

    (2)

    式中:r2ik為觀測(cè)點(diǎn)rk和格點(diǎn)ri之間距離的平方。

    逐步訂正法得到的最后分析值實(shí)際上是各種可利用的信息加權(quán)平均。逐步訂正法簡單經(jīng)濟(jì),能產(chǎn)生合理的分析。

    2 農(nóng)業(yè)旱情信息同化數(shù)據(jù)源

    農(nóng)業(yè)旱情信息來源主要來自于農(nóng)田水分循環(huán)過程模擬結(jié)果、遙感監(jiān)測(cè)信息以及對(duì)農(nóng)作物外在的觀察和描述。

    2.1 農(nóng)田水分循環(huán)模型模擬的旱情信息

    農(nóng)田水分循環(huán)模型可連續(xù)模擬逐日的土壤含水量、作物缺水量,可建立基于逐旬時(shí)間尺度的旱情評(píng)估指數(shù),比以往基于月尺度的干旱指數(shù)更能準(zhǔn)確描述干旱的發(fā)生、結(jié)束和程度,進(jìn)而為旱情監(jiān)測(cè)預(yù)測(cè)及抗旱決策提供更加及時(shí)、可靠的支持,達(dá)到有效減災(zāi)和防災(zāi)的目的。基于農(nóng)田水分循環(huán)過程模擬的農(nóng)業(yè)旱情,考慮降水、蒸發(fā)、植被和土壤特性對(duì)干旱的綜合影響,能夠更加真實(shí)地反映實(shí)際發(fā)生的旱情。

    從農(nóng)田水量平衡原理出發(fā),以農(nóng)田表層為原型,建立農(nóng)田水分循環(huán)模型,模擬農(nóng)作物生長期農(nóng)田水分循環(huán)過程及逐日農(nóng)作物需水、缺水信息。根據(jù)區(qū)域氣候和下墊面條件,可確定作物不同生長階段的逐旬農(nóng)業(yè)旱情等級(jí)[9]。

    2.2 遙感監(jiān)測(cè)的旱情信息

    遙感技術(shù)可以獲取時(shí)空連續(xù)的地表參量信息,反映時(shí)空尺度上地表狀態(tài)的變化。借助更高分辨率的MODIS遙感數(shù)據(jù),構(gòu)建作物生長季的葉面缺水指數(shù)(NDWI),可表征作物受旱狀況,為旱情監(jiān)測(cè)與驗(yàn)證提供重要信息源。

    獲取并處理與模擬分析同期的MODIS源數(shù)據(jù),MODIS1B數(shù)據(jù)處理過程包括輻射校正、幾何校正、大氣校正和云檢測(cè)等。綠色植物的反射光譜的光譜特征在0.9~2.5 μm區(qū)域是由液態(tài)水吸收控制的,同時(shí)也受一些其他的生物化學(xué)成分吸收的微弱影響。在近紅外波段,植被液態(tài)水的吸收可以忽略不計(jì),在短波紅外波段,水的吸收很弱,從而利用植被在紅外波段(NIR)與短波紅外波段的光譜可以很靈敏的反應(yīng)植被冠層水的含量。而葉面含水的多寡與水分脅迫有直接的關(guān)系。因此可以用葉面缺水指數(shù)(NDWI)進(jìn)行旱情監(jiān)測(cè)。根據(jù)監(jiān)測(cè)時(shí)段逐日NDWI指數(shù)數(shù)據(jù)集產(chǎn)品,采用平均值合成法得到旬NDWI指數(shù),分析NDWI指數(shù)與旱情等級(jí)關(guān)系,可得逐旬遙感監(jiān)測(cè)的旱情信息。

    3 農(nóng)業(yè)旱情信息同化過程與驗(yàn)證

    3.1 農(nóng)業(yè)旱情信息同化過程

    由于土壤結(jié)構(gòu)和質(zhì)地的復(fù)雜性及其空間不一致性,并受限于土壤參數(shù)及其資料,精確地模擬空間分布的農(nóng)業(yè)旱情是有困難的。由于模型計(jì)算的作物需水量和作物缺水量代表著模型計(jì)算單元一天24 h的平均值,而遙感觀測(cè)值僅代表計(jì)算單元瞬時(shí)的葉面缺水指數(shù)。考慮到農(nóng)田水分循環(huán)模型模擬得到的作物旱情等級(jí)與遙感監(jiān)測(cè)旱情信息的時(shí)空尺度差異,采用信息同化的逐步訂正法進(jìn)行模型參數(shù)調(diào)試,可提高農(nóng)業(yè)旱情模擬進(jìn)度。因此,經(jīng)同化的農(nóng)田水分循環(huán)模擬模型可用于農(nóng)業(yè)旱情評(píng)估與預(yù)測(cè)。農(nóng)業(yè)旱情信息同化逐步訂正法分析步驟如下。

    (1)采用農(nóng)田水分循環(huán)模擬模型計(jì)算逐日作物需水量和作物缺水量值,產(chǎn)生時(shí)間分辨率為1旬,空間分辨率為計(jì)算單元(縣級(jí)行政單元)的農(nóng)業(yè)旱情信息作為初估場(chǎng)。

    (2)采集MODIS遙感數(shù)據(jù),采用平均值合成法計(jì)算旬葉面缺水指數(shù)(NDWI),可得逐旬遙感監(jiān)測(cè)的旱情信息。

    (3)比對(duì)各計(jì)算單元模擬農(nóng)業(yè)旱情信息和遙感監(jiān)測(cè)旱情信息并計(jì)算誤差協(xié)方差,依次使用逐個(gè)計(jì)算單元遙感監(jiān)測(cè)值調(diào)整農(nóng)田水分循環(huán)模擬模型參數(shù),模擬出時(shí)間步長推進(jìn)1旬的農(nóng)業(yè)旱情信息,作為下一次訂正的初估場(chǎng),逐旬進(jìn)行下一次訂正。

    (4)從農(nóng)田水分循環(huán)模擬值和遙感監(jiān)測(cè)值的計(jì)算得到誤差場(chǎng),對(duì)誤差場(chǎng)的樣本進(jìn)行統(tǒng)計(jì)計(jì)算,農(nóng)業(yè)旱情信息同化直至誤差系列期望值為0、協(xié)方差為給定的先驗(yàn)誤差為止。

    3.2 同化后的模擬模型結(jié)果驗(yàn)證

    選擇山西的太原、晉中和呂梁三個(gè)地級(jí)市為研究示范區(qū),屬于黃河流域的汾河水系和東川河水系,區(qū)域總面積44 213 km2,結(jié)合水資源分區(qū)和行政區(qū)劃,將示范區(qū)分為24個(gè)計(jì)算單元(縣級(jí)行政單元),分別進(jìn)行農(nóng)田水分循環(huán)過程模擬和同期遙感影像采集。示范區(qū)主要種植小麥和玉米,2009年農(nóng)作物總播種面積86.34萬hm2。示范區(qū)多年平均降水量430~525 mm,多年平均蒸發(fā)量1 710~1 780 mm,以春夏連旱為主。根據(jù)1970-2009年各計(jì)算單元的逐日降水、蒸發(fā)資料,逐日模擬小麥玉米生長期需水、缺水過程,所選資料具有可靠性和代表性。

    在空間尺度上,利用仿真模擬的方法計(jì)算得到逐月的作物受旱情況與遙感影像監(jiān)測(cè)的農(nóng)作物受旱狀況有差異,以旱情較重的2001年為例(見圖1),運(yùn)用逐步訂正法,根據(jù)遙感監(jiān)測(cè)旱情信息進(jìn)行同化,調(diào)整優(yōu)化農(nóng)田水分循環(huán)模型狀態(tài)變量與參數(shù)。

    農(nóng)業(yè)旱情模擬結(jié)果的可靠性取決于基本資料參數(shù)的可靠程度和農(nóng)業(yè)干旱等級(jí)的劃分標(biāo)準(zhǔn)是否合適??捎玫湫透珊的陮?shí)測(cè)旱情對(duì)同化后的模擬模型進(jìn)行驗(yàn)證。山西示范區(qū)2001年3-7月實(shí)際旱情比較嚴(yán)重,將農(nóng)田水分循環(huán)模型模擬旱情信息與實(shí)際旱情進(jìn)行比較,24個(gè)單元中80%以上的單元旱情等級(jí)吻合較好(旱情等級(jí)為1、2、3、4分別表示輕旱、中旱、重旱和特旱),農(nóng)田水分循環(huán)模型模擬旱情等級(jí)驗(yàn)證準(zhǔn)確度表見表1,模擬旱情與同期實(shí)際旱情比對(duì)圖見圖2。研究結(jié)果表明:同化后的農(nóng)田水分循環(huán)模型模擬旱情等級(jí)與實(shí)測(cè)旱情結(jié)果具有較好的擬合度。

    表1 農(nóng)田水分循環(huán)模型模擬旱情等級(jí)驗(yàn)證準(zhǔn)確度Tab.1 Verify the accuracy of farmland water cycle model to simulate the drought level

    圖1 2001年4-6月山西示范區(qū)仿真模擬結(jié)果(左邊)及對(duì)應(yīng)時(shí)間的遙感監(jiān)測(cè)結(jié)果(右邊)Fig.1 Simulation result (left) and remote sensing image (right) of the model district of Shanxi among April to June in 2001

    圖2 山西示范區(qū)2001年3-7月模擬旱情與同期實(shí)際旱情比對(duì)圖Fig.2 Compares to the chart between simulation drought to actual occurrence one of the model district of Shanxi among March to July in 2001

    4 結(jié) 語

    (1)提出的農(nóng)田水分循環(huán)模型沒有考慮旱作區(qū)地表以下水平方向的水量交換。在模型應(yīng)用時(shí),對(duì)作物實(shí)際受旱情況需進(jìn)行監(jiān)視和跟蹤;運(yùn)用同化技術(shù),根據(jù)遙感觀測(cè)信息對(duì)模型參數(shù)進(jìn)行實(shí)時(shí)修正,仿真模擬旱情等級(jí)與實(shí)際旱情等級(jí)進(jìn)行了比對(duì),80%以上單元旱情等級(jí)相吻合,所以,可用同化后的農(nóng)田水分循環(huán)模擬模型準(zhǔn)確預(yù)測(cè)作物缺水率,預(yù)估由缺水造成的作物減產(chǎn)量。

    (2)仿真模擬技術(shù)的應(yīng)用是研究區(qū)域作物旱情規(guī)律的一種行之有效途徑。利用水分模擬結(jié)果,掌握作物旱情發(fā)生、發(fā)展和緩解的過程,分析研究區(qū)域作物旱情時(shí)空分布規(guī)律,可為農(nóng)業(yè)旱情預(yù)測(cè)預(yù)警提供技術(shù)支撐,為抗旱決策提供科學(xué)依據(jù)。

    [1] 康紹忠,劉曉明,熊運(yùn)章,等.土壤-植物-大氣連續(xù)體水分傳輸理論及其應(yīng)用[M].北京:水利出版社,1992.

    [2] 吳志勇,陸桂華,張建云,等.基于VIC模型的逐日土壤含水量模擬[J].地理科學(xué),2007,27(3):359-364

    [3] 張文勝,劉錦娣,王建朋,等.復(fù)雜多因素影響下區(qū)域旱情評(píng)價(jià)體系研究[J].人民黃河,2013,35(5):41-43.

    [4] Launay M, Guerif M. Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications[J]. Agriculture Ecosystems and Environment, 2005, 111(1/2/3/4):321-339.

    [5] Vazlkdoust M,van Dam J C,Bamlaanssen W, et al. Assimilation of satellite data into agrohydrological modes to improve crop yield forecasts[J]. International Journal of Remote Sensing, 2009,30(1):2 523-2 545.

    [6] 柳欽火,仲 波,吳紀(jì)桃,等.環(huán)境遙感定量反演與同化[M].北京:科學(xué)出版社,2011.

    [7] 王 文,寇小華.水文數(shù)據(jù)同化方法及遙感數(shù)據(jù)在水文數(shù)據(jù)同化中的應(yīng)用進(jìn)展[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,37(5):556-562.

    [8] 梁順林,李 新.陸面觀測(cè)、模擬與數(shù)據(jù)同化[M].北京:高等教育出版社,2013.

    [9] 倪深海,顧 穎,戚建國,等.農(nóng)業(yè)干旱動(dòng)態(tài)模擬技術(shù)研究與應(yīng)用[J].中國農(nóng)村水利水電,2012,(8):79-82.

    猜你喜歡
    旱情農(nóng)田作物
    人工增雨解旱情
    達(dá)爾頓老伯的農(nóng)田
    基于不同旱情指數(shù)的石羊河流域春旱監(jiān)測(cè)研究
    作物遭受霜凍該如何補(bǔ)救
    四種作物 北方種植有前景
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    農(nóng)田創(chuàng)意秀
    無人機(jī)遙感在作物監(jiān)測(cè)中的應(yīng)用與展望
    農(nóng)田搞養(yǎng)殖需辦哪些證
    農(nóng)田制作所
    亚洲一区二区三区不卡视频| 久久精品国产亚洲av涩爱 | 日韩欧美国产一区二区入口| 日韩人妻高清精品专区| 久久国产乱子免费精品| 露出奶头的视频| 岛国在线免费视频观看| 午夜精品久久久久久毛片777| 日韩亚洲欧美综合| 欧美中文日本在线观看视频| av黄色大香蕉| 亚洲男人的天堂狠狠| 狠狠狠狠99中文字幕| eeuss影院久久| 人人妻人人看人人澡| 特大巨黑吊av在线直播| 两人在一起打扑克的视频| 最后的刺客免费高清国语| 97热精品久久久久久| 精品久久久久久久末码| 十八禁人妻一区二区| 最近在线观看免费完整版| 1024手机看黄色片| av女优亚洲男人天堂| 免费人成视频x8x8入口观看| 国产一区二区亚洲精品在线观看| 久久久国产成人免费| 亚洲,欧美精品.| 丁香六月欧美| 国产成人aa在线观看| 精品人妻熟女av久视频| 久久久久久久精品吃奶| 久久久久久久午夜电影| 欧美黑人欧美精品刺激| 宅男免费午夜| 国产精品av视频在线免费观看| 久久精品国产自在天天线| 嫩草影视91久久| 国产精品一区二区性色av| 久久国产精品影院| 少妇高潮的动态图| 中文字幕高清在线视频| 又黄又爽又刺激的免费视频.| 国产中年淑女户外野战色| 啪啪无遮挡十八禁网站| 国产亚洲欧美在线一区二区| 91av网一区二区| 老熟妇仑乱视频hdxx| 91狼人影院| 久久久精品大字幕| 97超级碰碰碰精品色视频在线观看| av在线天堂中文字幕| 久久久久亚洲av毛片大全| 欧美不卡视频在线免费观看| 麻豆一二三区av精品| 天堂影院成人在线观看| 亚洲欧美清纯卡通| 免费高清视频大片| 国产亚洲精品久久久久久毛片| 在线观看舔阴道视频| 国产在线男女| 夜夜夜夜夜久久久久| 好男人在线观看高清免费视频| 国内毛片毛片毛片毛片毛片| 日韩成人在线观看一区二区三区| 在线观看66精品国产| 永久网站在线| 久久欧美精品欧美久久欧美| 婷婷精品国产亚洲av| 搡老妇女老女人老熟妇| 88av欧美| 国产av不卡久久| 欧美黑人欧美精品刺激| 99久久精品热视频| 色吧在线观看| 最近视频中文字幕2019在线8| 国产黄a三级三级三级人| 九九久久精品国产亚洲av麻豆| 国产蜜桃级精品一区二区三区| 色吧在线观看| 亚洲av美国av| 99国产综合亚洲精品| 久久午夜福利片| 桃色一区二区三区在线观看| 久久国产精品人妻蜜桃| 亚洲男人的天堂狠狠| 国产单亲对白刺激| 99久久无色码亚洲精品果冻| 精品午夜福利视频在线观看一区| 精品久久久久久久久av| 国产精品一区二区性色av| 日韩欧美在线二视频| 非洲黑人性xxxx精品又粗又长| 99在线人妻在线中文字幕| 少妇丰满av| 亚洲真实伦在线观看| 看十八女毛片水多多多| 欧美日韩瑟瑟在线播放| 亚洲三级黄色毛片| 色尼玛亚洲综合影院| 99在线视频只有这里精品首页| 一夜夜www| 精品乱码久久久久久99久播| 欧美bdsm另类| 日韩免费av在线播放| 亚洲av中文字字幕乱码综合| 十八禁国产超污无遮挡网站| 日本撒尿小便嘘嘘汇集6| 欧美成人a在线观看| 长腿黑丝高跟| 免费大片18禁| 国产日本99.免费观看| 真实男女啪啪啪动态图| 国产大屁股一区二区在线视频| 亚洲不卡免费看| 99久久成人亚洲精品观看| 精品午夜福利视频在线观看一区| 日本黄色片子视频| 日本撒尿小便嘘嘘汇集6| 免费在线观看亚洲国产| 久久久久国产精品人妻aⅴ院| 久久精品国产99精品国产亚洲性色| 免费观看的影片在线观看| 桃红色精品国产亚洲av| 91久久精品国产一区二区成人| 欧美性感艳星| 黄色视频,在线免费观看| 免费高清视频大片| 波多野结衣高清作品| 国产免费av片在线观看野外av| 精品久久久久久久久亚洲 | 欧美激情国产日韩精品一区| 91狼人影院| 亚洲成av人片在线播放无| 久久久精品大字幕| 免费观看精品视频网站| 日韩欧美在线乱码| 午夜日韩欧美国产| 一本综合久久免费| 欧美乱色亚洲激情| 少妇熟女aⅴ在线视频| 人妻制服诱惑在线中文字幕| 国产伦一二天堂av在线观看| 91狼人影院| 国产精品一区二区三区四区免费观看 | 蜜桃久久精品国产亚洲av| 99热这里只有精品一区| www.熟女人妻精品国产| 精品久久久久久久久久免费视频| 免费人成在线观看视频色| 亚洲午夜理论影院| 少妇熟女aⅴ在线视频| 麻豆成人av在线观看| 熟女电影av网| 999久久久精品免费观看国产| 亚洲中文字幕日韩| 91字幕亚洲| 欧美在线黄色| 搡女人真爽免费视频火全软件 | 少妇的逼好多水| 亚洲熟妇中文字幕五十中出| a级毛片a级免费在线| 麻豆国产97在线/欧美| 亚洲国产高清在线一区二区三| 99久久成人亚洲精品观看| 极品教师在线免费播放| 久久国产乱子免费精品| 久99久视频精品免费| 搡老岳熟女国产| 欧美最新免费一区二区三区 | 亚洲成人精品中文字幕电影| 91av网一区二区| 搡女人真爽免费视频火全软件 | 最后的刺客免费高清国语| 毛片女人毛片| 国产精品美女特级片免费视频播放器| 极品教师在线视频| 亚洲avbb在线观看| 亚洲av成人精品一区久久| 午夜免费男女啪啪视频观看 | 成人三级黄色视频| 久久国产乱子伦精品免费另类| av欧美777| 一进一出好大好爽视频| 亚洲18禁久久av| 亚洲国产欧美人成| 国产伦人伦偷精品视频| 我要搜黄色片| 欧美激情在线99| 午夜老司机福利剧场| 十八禁人妻一区二区| 精品欧美国产一区二区三| 桃色一区二区三区在线观看| 中文字幕免费在线视频6| 真实男女啪啪啪动态图| 看黄色毛片网站| 国产毛片a区久久久久| 12—13女人毛片做爰片一| 成年女人看的毛片在线观看| av天堂中文字幕网| 变态另类成人亚洲欧美熟女| 激情在线观看视频在线高清| 欧美激情国产日韩精品一区| 在线播放无遮挡| 亚洲成av人片在线播放无| 国产亚洲精品综合一区在线观看| 特级一级黄色大片| 又黄又爽又免费观看的视频| 欧美高清成人免费视频www| 久9热在线精品视频| 97超级碰碰碰精品色视频在线观看| 夜夜夜夜夜久久久久| 免费在线观看成人毛片| 国产日本99.免费观看| 久久99热6这里只有精品| 一级av片app| 日韩欧美在线二视频| 久久伊人香网站| 又爽又黄无遮挡网站| 中出人妻视频一区二区| 国产单亲对白刺激| 国产精品亚洲av一区麻豆| 成人鲁丝片一二三区免费| 久久久国产成人精品二区| 久久久国产成人精品二区| 欧美丝袜亚洲另类 | 久久九九热精品免费| 欧美日韩福利视频一区二区| 久99久视频精品免费| 成年女人毛片免费观看观看9| 日韩欧美在线二视频| 啦啦啦观看免费观看视频高清| 亚洲,欧美,日韩| 精品一区二区三区视频在线观看免费| 日韩欧美国产一区二区入口| 色在线成人网| 美女高潮喷水抽搐中文字幕| 哪里可以看免费的av片| 日韩欧美免费精品| 久久久久九九精品影院| 欧美色视频一区免费| 欧美xxxx性猛交bbbb| 免费黄网站久久成人精品 | 搡老岳熟女国产| 色吧在线观看| 我要看日韩黄色一级片| 免费搜索国产男女视频| 一本久久中文字幕| 桃色一区二区三区在线观看| 日韩欧美国产在线观看| 丁香六月欧美| av福利片在线观看| 久久精品人妻少妇| 男女下面进入的视频免费午夜| 亚洲18禁久久av| 亚洲黑人精品在线| 一进一出抽搐动态| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久久久久免费视频| 搞女人的毛片| 亚洲国产高清在线一区二区三| 国产三级中文精品| 91av网一区二区| 宅男免费午夜| 久久久久久大精品| 色播亚洲综合网| 亚洲精品粉嫩美女一区| 国产精品,欧美在线| 亚洲av美国av| 五月伊人婷婷丁香| 人人妻人人澡欧美一区二区| 亚洲 欧美 日韩 在线 免费| 国产一级毛片七仙女欲春2| 日韩人妻高清精品专区| 亚洲精品一区av在线观看| 精品久久久久久久久亚洲 | 久久亚洲精品不卡| 99热只有精品国产| 噜噜噜噜噜久久久久久91| 嫩草影院精品99| 欧美黑人巨大hd| 麻豆成人午夜福利视频| 欧美成人免费av一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产探花在线观看一区二区| 欧美+日韩+精品| 国产乱人视频| 国产蜜桃级精品一区二区三区| 欧美性猛交╳xxx乱大交人| 99riav亚洲国产免费| 夜夜夜夜夜久久久久| av在线老鸭窝| 国产亚洲精品综合一区在线观看| 国产黄片美女视频| 麻豆久久精品国产亚洲av| 中文字幕人成人乱码亚洲影| 琪琪午夜伦伦电影理论片6080| 男人的好看免费观看在线视频| 免费看日本二区| 成人美女网站在线观看视频| 免费看光身美女| 永久网站在线| 精品一区二区免费观看| 最近视频中文字幕2019在线8| 成年人黄色毛片网站| 国产精品亚洲av一区麻豆| 99久久精品一区二区三区| 国产精品1区2区在线观看.| 97碰自拍视频| 丰满的人妻完整版| 69人妻影院| 国产日本99.免费观看| 99热这里只有是精品50| 韩国av一区二区三区四区| 国产免费av片在线观看野外av| 亚洲内射少妇av| 国产爱豆传媒在线观看| 亚洲片人在线观看| 午夜福利免费观看在线| 国产成人影院久久av| 很黄的视频免费| 色尼玛亚洲综合影院| 伊人久久精品亚洲午夜| 欧美日本亚洲视频在线播放| 精品国产三级普通话版| 一个人免费在线观看的高清视频| 青草久久国产| 日本在线视频免费播放| 国产成+人综合+亚洲专区| 亚洲一区二区三区色噜噜| 亚洲片人在线观看| 欧美日韩乱码在线| 老司机午夜十八禁免费视频| 国产精品永久免费网站| 亚洲欧美日韩东京热| 此物有八面人人有两片| 变态另类丝袜制服| 18美女黄网站色大片免费观看| 美女高潮的动态| 级片在线观看| 好看av亚洲va欧美ⅴa在| 国产一级毛片七仙女欲春2| 亚洲国产精品sss在线观看| or卡值多少钱| 好男人在线观看高清免费视频| 亚洲美女搞黄在线观看 | 亚洲av日韩精品久久久久久密| 午夜免费激情av| 婷婷丁香在线五月| 国产伦一二天堂av在线观看| 国产老妇女一区| 国产91精品成人一区二区三区| 我要看日韩黄色一级片| 一级作爱视频免费观看| 色综合欧美亚洲国产小说| 赤兔流量卡办理| 日韩精品青青久久久久久| 九九久久精品国产亚洲av麻豆| 国产精品精品国产色婷婷| 欧美黑人巨大hd| 色视频www国产| 国产久久久一区二区三区| 国产av麻豆久久久久久久| 亚洲av二区三区四区| 久久欧美精品欧美久久欧美| 宅男免费午夜| 国产高清有码在线观看视频| 欧美色欧美亚洲另类二区| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 国产一区二区在线观看日韩| 床上黄色一级片| 嫩草影院入口| 精品午夜福利视频在线观看一区| 免费一级毛片在线播放高清视频| 亚洲经典国产精华液单 | 国产精品嫩草影院av在线观看 | 亚洲av电影在线进入| 毛片女人毛片| 欧美极品一区二区三区四区| 十八禁国产超污无遮挡网站| 黄色一级大片看看| 婷婷色综合大香蕉| 偷拍熟女少妇极品色| 一本一本综合久久| 欧美高清性xxxxhd video| 伊人久久精品亚洲午夜| 久久国产乱子伦精品免费另类| 舔av片在线| 18禁黄网站禁片午夜丰满| 成人欧美大片| 国产精品久久久久久亚洲av鲁大| 女生性感内裤真人,穿戴方法视频| 日韩有码中文字幕| 中文字幕高清在线视频| 嫩草影视91久久| 噜噜噜噜噜久久久久久91| 日本撒尿小便嘘嘘汇集6| 亚洲午夜理论影院| 日韩中字成人| 久久天躁狠狠躁夜夜2o2o| 免费看日本二区| 伦理电影大哥的女人| avwww免费| 久久久久国内视频| 国产色爽女视频免费观看| 欧美三级亚洲精品| 深夜精品福利| 亚洲av免费高清在线观看| 高潮久久久久久久久久久不卡| 一进一出抽搐gif免费好疼| 男女床上黄色一级片免费看| 午夜精品久久久久久毛片777| 日日摸夜夜添夜夜添小说| 男人舔奶头视频| 欧美zozozo另类| 十八禁国产超污无遮挡网站| 国内少妇人妻偷人精品xxx网站| 日韩欧美免费精品| 久久精品国产亚洲av天美| 亚洲人成网站在线播放欧美日韩| 日本 av在线| 最近在线观看免费完整版| 看黄色毛片网站| 神马国产精品三级电影在线观看| 国产精品,欧美在线| 香蕉av资源在线| 亚洲成人久久性| 国产在线男女| 好男人在线观看高清免费视频| 国产精品,欧美在线| 亚洲性夜色夜夜综合| 日韩 亚洲 欧美在线| 韩国av一区二区三区四区| 在线观看一区二区三区| 亚洲av成人av| 国产成人aa在线观看| 一本综合久久免费| 神马国产精品三级电影在线观看| 国产精品自产拍在线观看55亚洲| 少妇人妻精品综合一区二区 | 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 男女那种视频在线观看| 欧美一区二区精品小视频在线| 日韩欧美一区二区三区在线观看| 国产亚洲精品久久久com| 一卡2卡三卡四卡精品乱码亚洲| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久久亚洲 | 亚洲七黄色美女视频| av黄色大香蕉| 亚洲精品在线观看二区| 18美女黄网站色大片免费观看| 99精品在免费线老司机午夜| 日韩欧美精品v在线| 性色avwww在线观看| 美女被艹到高潮喷水动态| 91av网一区二区| 最后的刺客免费高清国语| 国产激情偷乱视频一区二区| 18禁黄网站禁片免费观看直播| www.www免费av| 成人精品一区二区免费| 3wmmmm亚洲av在线观看| 亚洲精品久久国产高清桃花| 99热精品在线国产| 精品久久久久久久久亚洲 | 国产精品久久久久久人妻精品电影| 三级毛片av免费| x7x7x7水蜜桃| 757午夜福利合集在线观看| 亚洲内射少妇av| 成人鲁丝片一二三区免费| 免费av不卡在线播放| 国产三级黄色录像| 免费黄网站久久成人精品 | 美女cb高潮喷水在线观看| 国产午夜精品论理片| 国产成年人精品一区二区| 国产精品久久久久久久久免 | 宅男免费午夜| 成年免费大片在线观看| 午夜福利在线观看免费完整高清在 | 国产熟女xx| 自拍偷自拍亚洲精品老妇| 亚洲五月婷婷丁香| 韩国av一区二区三区四区| 久久热精品热| 90打野战视频偷拍视频| 麻豆成人午夜福利视频| 免费在线观看日本一区| 日本黄大片高清| АⅤ资源中文在线天堂| 久9热在线精品视频| 久久久精品大字幕| 日韩精品青青久久久久久| 露出奶头的视频| 亚洲成a人片在线一区二区| 色5月婷婷丁香| 一a级毛片在线观看| 一边摸一边抽搐一进一小说| av黄色大香蕉| 波野结衣二区三区在线| 99久久精品热视频| 最后的刺客免费高清国语| 国产精品一区二区免费欧美| 国产午夜福利久久久久久| 免费高清视频大片| 久久久国产成人免费| 欧美最黄视频在线播放免费| 1000部很黄的大片| 精品久久国产蜜桃| 国产日本99.免费观看| 国产精品99久久久久久久久| 搡老岳熟女国产| 久久久久久久精品吃奶| 99久久成人亚洲精品观看| 好男人电影高清在线观看| 十八禁人妻一区二区| 国产精品一区二区三区四区免费观看 | 亚洲五月婷婷丁香| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 真实男女啪啪啪动态图| 久久国产乱子伦精品免费另类| 国产人妻一区二区三区在| 亚洲中文字幕日韩| 欧美黑人巨大hd| 国产午夜精品论理片| 长腿黑丝高跟| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 成人鲁丝片一二三区免费| 好男人电影高清在线观看| 青草久久国产| av福利片在线观看| 色播亚洲综合网| 国产探花极品一区二区| 俄罗斯特黄特色一大片| 91九色精品人成在线观看| 亚洲内射少妇av| 日本一本二区三区精品| 国产免费av片在线观看野外av| 首页视频小说图片口味搜索| 91午夜精品亚洲一区二区三区 | 亚洲七黄色美女视频| 色5月婷婷丁香| 精品久久久久久久久av| 尤物成人国产欧美一区二区三区| 国内精品一区二区在线观看| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3| 久久精品夜夜夜夜夜久久蜜豆| 夜夜夜夜夜久久久久| 国产高清视频在线播放一区| 性色avwww在线观看| 亚洲色图av天堂| av在线天堂中文字幕| 日日摸夜夜添夜夜添小说| 亚洲av成人av| 日本一本二区三区精品| 亚洲综合色惰| 欧美日韩亚洲国产一区二区在线观看| 国产午夜精品久久久久久一区二区三区 | 制服丝袜大香蕉在线| 看片在线看免费视频| 波多野结衣高清作品| 99热这里只有是精品50| 色综合亚洲欧美另类图片| 色综合欧美亚洲国产小说| 丁香欧美五月| 免费看日本二区| 18禁黄网站禁片免费观看直播| 麻豆av噜噜一区二区三区| 精品人妻1区二区| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 好男人在线观看高清免费视频| 国产视频一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 黄色日韩在线| 啪啪无遮挡十八禁网站| 搞女人的毛片| 国内毛片毛片毛片毛片毛片| 国产在线男女| 亚洲最大成人av| 人人妻人人看人人澡| 国产亚洲欧美在线一区二区| 免费一级毛片在线播放高清视频| 女生性感内裤真人,穿戴方法视频| 亚洲av成人精品一区久久| 久久中文看片网| 久久国产精品人妻蜜桃| 美女高潮喷水抽搐中文字幕| 十八禁国产超污无遮挡网站| АⅤ资源中文在线天堂| 精品久久久久久久人妻蜜臀av| 国产精品爽爽va在线观看网站| 搡老熟女国产l中国老女人| 观看美女的网站| 757午夜福利合集在线观看| 成人午夜高清在线视频| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 欧美不卡视频在线免费观看| 亚洲五月天丁香| 国产精品亚洲一级av第二区| 给我免费播放毛片高清在线观看| 日韩中文字幕欧美一区二区| 久久性视频一级片| 我的老师免费观看完整版| 欧美xxxx性猛交bbbb|