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

    基于MODIS NDVI的冬小麥?zhǔn)斋@指數(shù)遙感提取

    2022-07-29 08:58:00劉劍鋒張喜旺
    江蘇農(nóng)業(yè)科學(xué) 2022年13期
    關(guān)鍵詞:開(kāi)花期冬小麥籽粒

    劉劍鋒, 方 鵬, 陳 琳, 張喜旺

    (1.黃河水利職業(yè)技術(shù)學(xué)院,河南開(kāi)封 475004; 2.河南大學(xué)黃河中下游數(shù)字地理技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,河南開(kāi)封 475004)

    收獲指數(shù)(harvest index)是影響農(nóng)作物產(chǎn)量十分重要的理化參數(shù),也是諸多作物估產(chǎn)模型的關(guān)鍵參數(shù)之一。作物收獲指數(shù)是農(nóng)作物在收獲期的籽??傎|(zhì)量和作物總生物量的比值,反映了作物光合同化產(chǎn)物在籽粒和營(yíng)養(yǎng)器官上的分配比例。根據(jù)作物生理學(xué)理論,生物量表示作物生理潛力,籽粒質(zhì)量由籽粒庫(kù)容限定,所以收獲指數(shù)標(biāo)志著經(jīng)濟(jì)產(chǎn)物與同化產(chǎn)物的效率,也反映了源的生理效能,源、庫(kù)、流在一定程度上會(huì)影響作物的產(chǎn)量?,F(xiàn)代小麥育種不斷關(guān)注小麥?zhǔn)斋@指數(shù)與小麥生物量的關(guān)系。一些研究表明,收獲指數(shù)與產(chǎn)量之間是正相關(guān)的,收獲指數(shù)不斷提高也會(huì)促進(jìn)小麥產(chǎn)量的增長(zhǎng)。因此,關(guān)注和研究冬小麥?zhǔn)斋@指數(shù),對(duì)于準(zhǔn)確理解冬小麥生理機(jī)制和冬小麥增產(chǎn)、估產(chǎn)具有重要意義。

    傳統(tǒng)農(nóng)業(yè)大多通過(guò)實(shí)地測(cè)量方法獲取小麥?zhǔn)斋@指數(shù),這種方法費(fèi)時(shí)、費(fèi)力、效率較低而且難以大范圍獲取小麥?zhǔn)斋@指數(shù)。也有部分學(xué)者通過(guò)氣象資料分析構(gòu)建冬小麥氣象模型來(lái)計(jì)算收獲指數(shù)。這種方法綜合考慮了各個(gè)氣象因素對(duì)收獲指數(shù)的影響力,但是需要大量的年際氣象累積數(shù)據(jù),而且能否適用于小區(qū)域還有待驗(yàn)證。除此之外,更多的農(nóng)業(yè)研究者關(guān)注的是小麥生長(zhǎng)自然環(huán)境、遺傳基因等因素對(duì)小麥?zhǔn)斋@指數(shù)的影響,而鮮見(jiàn)利用遙感技術(shù)提取小麥?zhǔn)斋@指數(shù)。

    遙感技術(shù)的快速發(fā)展為諸多研究領(lǐng)域提供了新的研究手段,但國(guó)內(nèi)外利用遙感數(shù)據(jù)提取冬小麥?zhǔn)斋@指數(shù)的文章偏少,可用模型精度也不高。Moriondo等通過(guò)分析冬小麥開(kāi)花期至成熟期和出苗期到開(kāi)花期的歸一化植被指數(shù)(NDVI)變化關(guān)系,建立了收獲指數(shù)與NDVI的模型,但側(cè)重考慮生長(zhǎng)過(guò)程中溫度和水的負(fù)面影響,因此適應(yīng)水分脅迫較大的地區(qū),對(duì)于其他區(qū)域的效果有待驗(yàn)證。Walter等采用新穎的攝影測(cè)量方法,通過(guò)數(shù)字?jǐn)z影相機(jī)和相關(guān)軟件建立小面積的小麥3D點(diǎn)云圖來(lái)預(yù)測(cè)生物量、植株高度、收獲指數(shù),取得了較好的效果。但攝影相機(jī)成像區(qū)域較小,難以大面積推廣。Richards等通過(guò)對(duì)干旱環(huán)境下冬小麥生物量、葉面積、收獲指數(shù)等關(guān)系的分析,提出了一種非線性的模型來(lái)計(jì)算收獲指數(shù),因此適用于干旱環(huán)境,且具有經(jīng)驗(yàn)偏差。Li等綜合分析了現(xiàn)有的幾類冬小麥?zhǔn)斋@指數(shù)反演模型在不同土壤氮肥條件下的反演精度,指出了利用NDVI反演冬小麥?zhǔn)斋@指數(shù)是最可行的。任建強(qiáng)等基于NDVI時(shí)序數(shù)據(jù),提出了一種冬小麥生殖階段與營(yíng)養(yǎng)階段比值的收獲指數(shù)參數(shù)模型,但該模型忽略了冬小麥生殖生長(zhǎng)階段的籽粒質(zhì)量也是總生物量的一部分,也導(dǎo)致了反演結(jié)果偏低。

    總體而言,國(guó)內(nèi)外利用遙感衛(wèi)星數(shù)據(jù)提取冬小麥?zhǔn)斋@指數(shù)的研究偏少,精度普遍偏低。為了改進(jìn)收獲指數(shù)遙感反演模型,提高利用衛(wèi)星數(shù)據(jù)反演冬小麥?zhǔn)斋@指數(shù)精度,從而更加深入地解析遙感數(shù)據(jù)與收獲指數(shù)的關(guān)聯(lián)機(jī)制,本研究依據(jù)前人研究基礎(chǔ),提出了一種新的構(gòu)建收獲指數(shù)的參數(shù)模型,建立了參數(shù)模型與實(shí)測(cè)數(shù)據(jù)的定量關(guān)系,以期為冬小麥?zhǔn)斋@指數(shù)遙感反演工作提供思路。

    1 研究區(qū)與數(shù)據(jù)

    1.1 研究區(qū)

    館陶縣地處華北平原,與山東省接壤,臨近河南省,地理位置為36°46′35″~36°27′8″N,115°7′12″~115°28′ 5″E,總面積456 km。館陶縣屬于古河道沖積平原,地勢(shì)西南高,東北低,平均海拔約為 40 m??h內(nèi)有衛(wèi)運(yùn)河和章河流經(jīng),衛(wèi)運(yùn)河縣界內(nèi)全長(zhǎng)40.5 km,章河全長(zhǎng)4.2 km,屬季節(jié)性河流。氣候?yàn)榈湫偷呐瘻貛О霛駶?rùn)大陸性季風(fēng)氣候,夏季炎熱干燥,冬季寒冷,日照充足,四季分明,年平均氣溫14 ℃,全年無(wú)霜期200 d,年日照2 557 h。良好的自然條件和氣候條件,使其成為華北平原重要的冬小麥種植區(qū)(圖1)。

    1.2 數(shù)據(jù)

    1.2.1 遙感影像數(shù)據(jù) 本研究所使用的遙感影像數(shù)據(jù)為16 d合成的MODIS13Q1 VI數(shù)據(jù),來(lái)源于美國(guó)地質(zhì)調(diào)查局(USGS,https://earthexplorer.usgs.gov/),選取250 m分辨率的NDVI。該數(shù)據(jù)采用限定視角內(nèi)最大值合成法(CV-MVC),在合成時(shí)間段內(nèi)擁有更大的NDVI像元值、更小的觀測(cè)視角、更接近星下點(diǎn),并且能最大程度地減小雙向反射的影響,使合成的NDVI數(shù)據(jù)精度更高。時(shí)間段分別為2015年10月16日至2016年6月25日、2016年10月16日至2017年6月25日、2017年10月16日至2018年6月25日。每期數(shù)據(jù)16景影像,共計(jì)48景影像。

    1.2.2 實(shí)測(cè)收獲指數(shù) 本研究采用2016年、2017年、2018年連續(xù)3年的地面實(shí)測(cè)收獲指數(shù)數(shù)據(jù)??紤]到MODIS 13Q1 NDVI數(shù)據(jù)的空間分辨率為 250 m,研究區(qū)內(nèi)設(shè)置10個(gè)邊長(zhǎng)為250 m的大樣方,每個(gè)樣方內(nèi)設(shè)置5個(gè)樣點(diǎn)。每個(gè)樣點(diǎn)實(shí)測(cè)種植密度并取10株樣本,烘箱設(shè)置溫度為105 ℃殺青1 h,再用80 ℃烘干24 h至恒質(zhì)量,得到小麥的所有干質(zhì)量生物量。通過(guò)實(shí)測(cè)每株的籽粒質(zhì)量和生物量獲取實(shí)測(cè)收獲指數(shù)。

    1.2.3 冬小麥空間分布 基于外業(yè)調(diào)查資料,利用Landsat遙感影像進(jìn)行冬小麥種植區(qū)識(shí)別,獲得2016—2018年的河北省館陶縣冬小麥種植區(qū)分布數(shù)據(jù)。

    2 研究方法

    2.1 時(shí)序數(shù)據(jù)平滑處理

    MODIS NDVI雖然經(jīng)過(guò)最大值合成,但提取的時(shí)序數(shù)據(jù)噪聲依然很大,呈現(xiàn)明顯的鋸齒狀波動(dòng),從而影響定量反演冬小麥?zhǔn)斋@指數(shù)的精度。針對(duì)這種情況,經(jīng)比較現(xiàn)有方法的優(yōu)缺點(diǎn),本研究選取Savitzky-Golay濾波對(duì)時(shí)序數(shù)據(jù)進(jìn)行平滑處理,以便更清晰地描述冬小麥生長(zhǎng)關(guān)鍵時(shí)期時(shí)間序列特征曲線的變化。處理模型如下:

    (1)

    式中:+是擬合之前的第(+)個(gè)NDVI時(shí)間序列數(shù)據(jù);′是擬合之后的NDVI時(shí)間序列數(shù)據(jù);為滑動(dòng)窗口的寬度;為濾波器長(zhǎng)度;為濾波卷積系數(shù)。

    2.2 模型構(gòu)建

    2.2.1 模型理論基礎(chǔ) 從作物生理學(xué)角度看,作物收獲指數(shù)是收獲期籽粒總質(zhì)量和作物地上總生物量的比值,實(shí)質(zhì)上反映了作物干物質(zhì)在籽粒和營(yíng)養(yǎng)器官上的分配比例。冬小麥在一個(gè)生長(zhǎng)周期內(nèi)要經(jīng)過(guò)漫長(zhǎng)的生長(zhǎng)階段,其中返青期至開(kāi)花期、開(kāi)花期至乳熟期是最重要的2個(gè)生長(zhǎng)階段。返青期至開(kāi)花期階段,小麥營(yíng)養(yǎng)生長(zhǎng)器官發(fā)育完善,奠定了小麥生長(zhǎng)過(guò)程的生物量基礎(chǔ);而在開(kāi)花期至乳熟期階段,小麥營(yíng)養(yǎng)生長(zhǎng)器官發(fā)育達(dá)到極致后,小麥干物質(zhì)逐漸集中到籽粒部分,這是小麥產(chǎn)量形成的關(guān)鍵階段?;谶@2個(gè)關(guān)鍵階段開(kāi)展收獲指數(shù)研究成為探討遙感反演方法的焦點(diǎn)。

    2.2.2 現(xiàn)有模型存在的問(wèn)題 國(guó)內(nèi)外相關(guān)研究不多,本研究對(duì)現(xiàn)有收獲指數(shù)模型進(jìn)行合理改進(jìn)和優(yōu)化,提出了一種新的模型,以期實(shí)現(xiàn)模型在精度方面的提升。

    目前城建模型如公式(2)、公式(3)所示。

    =-×(1-/)。

    (2)

    式中:為收獲指數(shù)模擬參數(shù);為理論最大收獲指數(shù),取值0.48;為收獲指數(shù)最大值和最小值之差,取值0.18;為開(kāi)花期到成熟期NDVI均值;為返青期到開(kāi)花期NDVI均值。該模型考慮到冬小麥生長(zhǎng)過(guò)程水分脅迫對(duì)收獲指數(shù)的影響,因此最大收獲指數(shù)當(dāng)作某一常數(shù)以減輕水分脅迫帶來(lái)的負(fù)面影響。然而,該模型的區(qū)域適應(yīng)性不足,且獲取困難。

    =∑/∑NDVI。

    (3)

    式中:為收獲指數(shù)模擬參數(shù);∑為冬小麥開(kāi)花期到乳熟期的NDVI累計(jì)值;∑為返青期到開(kāi)花期的NDVI累計(jì)值。該模型依據(jù)冬小麥生長(zhǎng)階段的時(shí)序NDVI累計(jì)值計(jì)算收獲指數(shù),在區(qū)域尺度獲取冬小麥?zhǔn)斋@指數(shù)也是可行的。然而,只考慮到生物量與營(yíng)養(yǎng)生長(zhǎng)階段有關(guān),而忽略了生殖生長(zhǎng)階段所產(chǎn)生的籽粒質(zhì)量也是生物量的一部分這一關(guān)鍵問(wèn)題。

    2.2.3模型構(gòu)建 基于上述分析,根據(jù)作物生理學(xué)和農(nóng)學(xué)理論,冬小麥地上總生物量不僅包括營(yíng)養(yǎng)生長(zhǎng)階段的干物質(zhì)在莖桿部分的分配,還包括生殖生長(zhǎng)階段的干物質(zhì)在籽粒部分的分配。因此,利用與+的比值建立的參數(shù)模型,不僅具有嚴(yán)謹(jǐn)?shù)睦碚摶A(chǔ),而且實(shí)際反演冬小麥?zhǔn)斋@指數(shù)時(shí)更加簡(jiǎn)單易行。

    由于NDVI與生物量之間的相關(guān)性,以及現(xiàn)有模型的經(jīng)驗(yàn),本研究依照冬小麥生長(zhǎng)期內(nèi)營(yíng)養(yǎng)生長(zhǎng)階段和生殖生長(zhǎng)階段的時(shí)序NDVI特征曲線的變化,建立參數(shù)模型,公式如下。

    (4)

    式中:為收獲指數(shù)模擬參數(shù);是冬小麥開(kāi)花期至成熟期的NDVI均值,用以表征冬小麥生殖生長(zhǎng)階段生產(chǎn)的籽粒質(zhì)量;是冬小麥返青期到開(kāi)花期的NDVI均值,用以表征冬小麥生長(zhǎng)階段的生產(chǎn)生物量。

    3 結(jié)果與分析

    3.1 提取模型參數(shù)

    利用平滑后的2016年和2017年的NDVI時(shí)間序列,基于收獲指數(shù)參數(shù)模型(2)、(3)、(4)計(jì)算得到冬小麥?zhǔn)斋@指數(shù)模擬參數(shù),如圖2、圖3、圖4所示。

    經(jīng)統(tǒng)計(jì),在上述2期的收獲指數(shù)空間分布圖中,模型(2)所提取的收獲指數(shù)范圍為0.35~0.65,模型(3)提取的參數(shù)范圍為0.36~0.65,模型(4)提取的參數(shù)范圍為0.43~0.66。由于模型都是依據(jù)NDVI時(shí)間序列數(shù)據(jù)變化來(lái)響應(yīng)冬小麥營(yíng)養(yǎng)生長(zhǎng)階段和生殖生長(zhǎng)階段的變化,都具農(nóng)學(xué)理論支撐,所得數(shù)據(jù)基本合理。

    3.2 模擬參數(shù)與實(shí)測(cè)值的擬合

    基于上述計(jì)算結(jié)果,以樣方為基礎(chǔ)與實(shí)測(cè)值進(jìn)行回歸分析,探測(cè)各模型的擬合程度,結(jié)果如圖5、圖6、圖7所示。

    由圖5可知,基于模型(2)的模擬參數(shù)與實(shí)測(cè)收獲指數(shù)的擬合結(jié)果,2016年為0.289 3,2017年為0.302 1??傮w而言,利用模型(2)反演得到的收獲指數(shù)具有一定的精度,但是并不高。該模型最初是由Moriondo等提出并將其應(yīng)用到意大利格羅塞托省和福賈省的小麥?zhǔn)罩笖?shù)研究。該區(qū)域氣候類型為典型的地中海氣候,夏季炎熱干燥,降水較少,小麥?zhǔn)芩窒拗戚^大;冬季溫和多雨,適合小麥生長(zhǎng)。與我國(guó)華北平原相比,該區(qū)域在冬季氣溫偏高,降水較多,具有良好的小麥生長(zhǎng)條件,且該區(qū)域種植的冬小麥沒(méi)有返青期,因此該模型應(yīng)用到華北平原冬小麥種植區(qū)域的收獲指數(shù)反演時(shí)精度偏低。

    圖6、圖7為基于模型(3)和(4)得到的模擬參數(shù)與實(shí)測(cè)收獲指數(shù)的擬合結(jié)果。2016年分別為0.389 2和0.423 0;2017年分別為0.590 4和0.676 6。

    與模型(2)相比,模型(3)和(4)側(cè)重考慮了冬小麥的生長(zhǎng)機(jī)制和農(nóng)學(xué)基礎(chǔ),根據(jù)冬小麥NDVI曲線變化與冬小麥生物量變化的關(guān)系來(lái)構(gòu)建收獲指數(shù)反演模型,具有較強(qiáng)的理論支撐。因此,擬合優(yōu)度有了很大的提高。而且模型(4)在模型(3)的基礎(chǔ)上進(jìn)行了改進(jìn),具有更加嚴(yán)謹(jǐn)?shù)霓r(nóng)學(xué)依據(jù)。

    綜合上述分析可知,本研究構(gòu)建的參數(shù)模型與實(shí)測(cè)收獲指數(shù)數(shù)據(jù)擬合程度最高,其最高達(dá)到0.676 6。進(jìn)而證明了參數(shù)模型反演收獲指數(shù)的優(yōu)越性。

    3.3 預(yù)測(cè)模型建立

    根據(jù)2016—2017年實(shí)測(cè)數(shù)據(jù)以及的計(jì)算結(jié)果,構(gòu)建預(yù)測(cè)冬小麥?zhǔn)斋@指數(shù)的線性預(yù)測(cè)函數(shù),如式(5)所示。

    =0057 2+0511 2。

    (5)

    式中:為預(yù)測(cè)的收獲指數(shù);為從時(shí)序數(shù)據(jù)中提取的參數(shù)。

    3.4 收獲指數(shù)空間分布

    根據(jù)模型(4)計(jì)算得到的2018年研究區(qū)冬小麥模擬參數(shù),帶入公式(5)中,并以冬小麥種植區(qū)為掩模,提取得到2018年館陶縣冬小麥?zhǔn)斋@指數(shù)空間分布,如圖8所示。可以看出,2018年的冬小麥?zhǔn)斋@指數(shù)范圍主要集中在0.53~0.54,與2018年實(shí)測(cè)的館陶縣冬小麥?zhǔn)斋@指數(shù)數(shù)值分布基本一致,反演精度較高。從空間分布情況看,中部較高,四周偏低;北部冬小麥?zhǔn)斋@指數(shù)比南部冬小麥?zhǔn)斋@指數(shù)高。

    將2018年實(shí)測(cè)的冬小麥?zhǔn)斋@指數(shù)與預(yù)測(cè)結(jié)果進(jìn)行分析,結(jié)果如圖9所示??梢钥闯?,預(yù)測(cè)結(jié)果與2018年實(shí)測(cè)的收獲指數(shù)呈現(xiàn)出明顯的線性關(guān)系,為0.632 9。分析結(jié)果表明,通過(guò)參數(shù)模型建立的2018年冬小麥?zhǔn)斋@指數(shù)模型能很好地反映實(shí)際的收獲指數(shù)情況。

    4 結(jié)語(yǔ)

    利用遙感技術(shù)開(kāi)展作物收獲指數(shù)的研究尚不夠成熟,目前的研究精度普遍不高,因此很有進(jìn)行深入研究的必要。本研究依據(jù)農(nóng)學(xué)和作物生理學(xué)原理,利用冬小麥營(yíng)養(yǎng)生長(zhǎng)階段與全生育階段的NDVI均值之比,作為構(gòu)建與收獲指數(shù)相關(guān)的參數(shù)模型,相對(duì)于現(xiàn)有模型,理論上更加合理。從驗(yàn)證結(jié)果看,該模型反演精度也最高,為今后的冬小麥?zhǔn)斋@指數(shù)反演提供了一種可行的方法。

    但是在利用遙感技術(shù)反演冬小麥?zhǔn)斋@指數(shù)時(shí),還應(yīng)注意以下問(wèn)題:(1)目前國(guó)內(nèi)外對(duì)于收獲指數(shù)的研究很少涉及到作物生長(zhǎng)環(huán)境因素對(duì)作物收獲指數(shù)的影響,如氣溫、降水、蒸散發(fā)等環(huán)境因素。在今后的研究中應(yīng)加以重視。(2)MODIS13Q1 NDVI雖然經(jīng)過(guò)最大值合成法處理,但依然存在一定的精度問(wèn)題。通過(guò)Savitzky-Golay濾波對(duì)NDVI時(shí)間序列數(shù)據(jù)進(jìn)行平滑,剔除部分異常值,獲得更接近實(shí)際的冬小麥時(shí)序特征曲線對(duì)于研究收獲指數(shù)具有明顯的益處。最后,本研究提出的研究方案,從理論上來(lái)講,也可用于玉米、水稻、高粱等糧食作物,有待進(jìn)一步研究。

    猜你喜歡
    開(kāi)花期冬小麥籽粒
    高寒草原針茅牧草花期物候變化特征及其影響因子分析
    籽粒莧的飼用價(jià)值和高產(chǎn)栽培技術(shù)
    籽粒莧的特性和種植技術(shù)
    玉米機(jī)械脫粒籽粒含水量與破碎率的相關(guān)研究
    商麥1619 籽粒灌漿的特性
    甘肅冬小麥田
    SOLVABILITY OF A PARABOLIC-HYPERBOLIC TYPE CHEMOTAXIS SYSTEM IN 1-DIMENSIONAL DOMAIN?
    初春氣象條件對(duì)蘋果開(kāi)花期的影響分析
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    冬小麥——新冬18號(hào)
    国产黄色免费在线视频| 最后的刺客免费高清国语| 免费观看a级毛片全部| 人妻制服诱惑在线中文字幕| 极品教师在线视频| 欧美成人a在线观看| 日韩大片免费观看网站| 国产爱豆传媒在线观看| 久久久久久久久久成人| 最近最新中文字幕大全电影3| 少妇熟女欧美另类| 色婷婷久久久亚洲欧美| 欧美激情在线99| 嫩草影院新地址| 亚洲内射少妇av| 搡女人真爽免费视频火全软件| 亚洲成人中文字幕在线播放| 亚洲av日韩在线播放| 欧美+日韩+精品| 80岁老熟妇乱子伦牲交| 内地一区二区视频在线| 超碰97精品在线观看| 97精品久久久久久久久久精品| 五月天丁香电影| 久久久久久久精品精品| 欧美另类一区| 99久久中文字幕三级久久日本| 亚洲真实伦在线观看| 国内揄拍国产精品人妻在线| 欧美97在线视频| 女人被狂操c到高潮| 亚洲va在线va天堂va国产| 超碰av人人做人人爽久久| 国产精品三级大全| 精品人妻偷拍中文字幕| 亚洲经典国产精华液单| 黄色一级大片看看| 亚洲在久久综合| 国产探花在线观看一区二区| 国产淫语在线视频| 在线观看一区二区三区激情| 老师上课跳d突然被开到最大视频| 国产大屁股一区二区在线视频| 久热久热在线精品观看| 亚洲国产色片| 久久这里有精品视频免费| 国产精品伦人一区二区| 亚洲成人久久爱视频| 亚洲色图综合在线观看| 色视频在线一区二区三区| 一级二级三级毛片免费看| 一级毛片电影观看| 高清av免费在线| 精品久久国产蜜桃| 国产精品久久久久久精品电影| 欧美日韩视频精品一区| 日韩一本色道免费dvd| av天堂中文字幕网| 人妻制服诱惑在线中文字幕| 成人亚洲欧美一区二区av| 网址你懂的国产日韩在线| 成人午夜精彩视频在线观看| 国产亚洲av片在线观看秒播厂| 99久久人妻综合| 男人狂女人下面高潮的视频| 国产片特级美女逼逼视频| 99久国产av精品国产电影| 精品久久久久久久末码| 久久ye,这里只有精品| 亚洲内射少妇av| 大码成人一级视频| 中文精品一卡2卡3卡4更新| 少妇人妻久久综合中文| 久久99热这里只有精品18| 国产真实伦视频高清在线观看| 亚洲av日韩在线播放| h日本视频在线播放| 内射极品少妇av片p| 国产成人一区二区在线| 人妻 亚洲 视频| 在线播放无遮挡| 高清毛片免费看| 精品99又大又爽又粗少妇毛片| 国产黄色视频一区二区在线观看| 午夜老司机福利剧场| 免费观看无遮挡的男女| 啦啦啦中文免费视频观看日本| 色视频在线一区二区三区| 国产av国产精品国产| 啦啦啦中文免费视频观看日本| 久久久精品欧美日韩精品| 一边亲一边摸免费视频| 日韩制服骚丝袜av| 亚洲精品,欧美精品| 日日啪夜夜撸| 国产黄频视频在线观看| 日本与韩国留学比较| 日本色播在线视频| 一区二区三区精品91| 亚洲欧美成人综合另类久久久| 国产精品一及| 777米奇影视久久| 菩萨蛮人人尽说江南好唐韦庄| 精品一区二区免费观看| 3wmmmm亚洲av在线观看| 免费高清在线观看视频在线观看| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久久久按摩| 天天躁夜夜躁狠狠久久av| 2022亚洲国产成人精品| 国产成人精品久久久久久| 国产精品一区二区性色av| freevideosex欧美| 日韩av在线免费看完整版不卡| av播播在线观看一区| 国产真实伦视频高清在线观看| 最新中文字幕久久久久| 亚洲av成人精品一二三区| 久久99精品国语久久久| 日韩中字成人| 午夜日本视频在线| 国内精品美女久久久久久| 久久久久久九九精品二区国产| 久久99精品国语久久久| 日韩欧美精品v在线| 少妇人妻一区二区三区视频| 男女下面进入的视频免费午夜| 视频区图区小说| 黑人高潮一二区| 久久ye,这里只有精品| 亚洲欧美日韩卡通动漫| 久久韩国三级中文字幕| videossex国产| 王馨瑶露胸无遮挡在线观看| 黄色视频在线播放观看不卡| 一级黄片播放器| 美女高潮的动态| 国产色爽女视频免费观看| 亚洲人成网站在线观看播放| 国产欧美日韩精品一区二区| 国产在线男女| 亚州av有码| 成人黄色视频免费在线看| 亚洲精品日韩在线中文字幕| 日产精品乱码卡一卡2卡三| 亚洲精品自拍成人| 神马国产精品三级电影在线观看| 免费高清在线观看视频在线观看| 久久精品久久久久久噜噜老黄| 久久精品国产自在天天线| av在线观看视频网站免费| 熟女电影av网| 五月开心婷婷网| 日本黄大片高清| 寂寞人妻少妇视频99o| 2022亚洲国产成人精品| 3wmmmm亚洲av在线观看| 国产精品一及| 我的女老师完整版在线观看| 精品久久久久久久久av| 日韩强制内射视频| 久久综合国产亚洲精品| 免费观看性生交大片5| 性色avwww在线观看| 男人爽女人下面视频在线观看| 免费看日本二区| 一级爰片在线观看| 亚洲精华国产精华液的使用体验| 亚洲国产日韩一区二区| 成人毛片60女人毛片免费| 少妇丰满av| 最近最新中文字幕免费大全7| 国产黄a三级三级三级人| 精品久久久久久久久av| 国国产精品蜜臀av免费| 成人特级av手机在线观看| 一级a做视频免费观看| 91久久精品国产一区二区三区| a级毛片免费高清观看在线播放| 亚洲国产欧美人成| 日韩强制内射视频| 视频中文字幕在线观看| 黄色怎么调成土黄色| 免费观看无遮挡的男女| 18禁裸乳无遮挡免费网站照片| 国产精品嫩草影院av在线观看| 下体分泌物呈黄色| 日韩一区二区三区影片| 久热久热在线精品观看| av线在线观看网站| 国产探花在线观看一区二区| 特级一级黄色大片| 99精国产麻豆久久婷婷| 在线看a的网站| 97人妻精品一区二区三区麻豆| 搡女人真爽免费视频火全软件| 国产av不卡久久| 欧美国产精品一级二级三级 | 菩萨蛮人人尽说江南好唐韦庄| 日韩中字成人| 欧美日韩视频精品一区| 国产真实伦视频高清在线观看| 99热网站在线观看| 青青草视频在线视频观看| 一级毛片我不卡| 狂野欧美白嫩少妇大欣赏| 久久99热6这里只有精品| 国产爽快片一区二区三区| 国产免费又黄又爽又色| 久热这里只有精品99| 伦精品一区二区三区| 日韩国内少妇激情av| 十八禁网站网址无遮挡 | 丰满人妻一区二区三区视频av| 久久人人爽人人爽人人片va| 18禁在线播放成人免费| 禁无遮挡网站| 极品少妇高潮喷水抽搐| 中文字幕人妻熟人妻熟丝袜美| 大又大粗又爽又黄少妇毛片口| 久久久久国产网址| 成人黄色视频免费在线看| 一级a做视频免费观看| 中国国产av一级| 亚洲va在线va天堂va国产| 少妇人妻精品综合一区二区| 欧美日韩亚洲高清精品| 久久久久久久亚洲中文字幕| 久久久精品免费免费高清| 日韩成人伦理影院| 欧美老熟妇乱子伦牲交| 久久鲁丝午夜福利片| 国产精品久久久久久久电影| 日韩av在线免费看完整版不卡| 久久国产乱子免费精品| 我的女老师完整版在线观看| 国产免费视频播放在线视频| 啦啦啦啦在线视频资源| 国产精品国产av在线观看| 亚洲经典国产精华液单| 亚洲精品久久久久久婷婷小说| 一级毛片 在线播放| 午夜亚洲福利在线播放| 国产淫片久久久久久久久| 久热久热在线精品观看| 久久99热这里只频精品6学生| 久久ye,这里只有精品| 国产精品久久久久久av不卡| 少妇熟女欧美另类| av专区在线播放| 欧美日韩综合久久久久久| 久久99精品国语久久久| 国产精品久久久久久精品电影| 亚洲av一区综合| 少妇 在线观看| 黄色欧美视频在线观看| 成人亚洲欧美一区二区av| 美女国产视频在线观看| 18禁在线播放成人免费| 51国产日韩欧美| 日本一本二区三区精品| 人妻制服诱惑在线中文字幕| 国内少妇人妻偷人精品xxx网站| 国内揄拍国产精品人妻在线| 丰满乱子伦码专区| 国产淫片久久久久久久久| av黄色大香蕉| 中国美白少妇内射xxxbb| 熟女电影av网| 少妇丰满av| av一本久久久久| 嫩草影院入口| 久久精品国产亚洲网站| 国产大屁股一区二区在线视频| 日韩成人伦理影院| 日本-黄色视频高清免费观看| 亚洲天堂av无毛| 麻豆国产97在线/欧美| 欧美zozozo另类| 亚洲欧美中文字幕日韩二区| 久久久久久久国产电影| 熟女av电影| 91aial.com中文字幕在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 真实男女啪啪啪动态图| 亚洲国产日韩一区二区| 国产极品天堂在线| 国产一区二区亚洲精品在线观看| 欧美日韩国产mv在线观看视频 | 简卡轻食公司| 一级爰片在线观看| 免费看光身美女| 成人毛片60女人毛片免费| 99热国产这里只有精品6| 久久99热这里只频精品6学生| 婷婷色av中文字幕| 欧美日本视频| 午夜福利在线观看免费完整高清在| 高清日韩中文字幕在线| 国产成人免费无遮挡视频| 成人鲁丝片一二三区免费| www.av在线官网国产| 久久久久久伊人网av| 午夜福利视频1000在线观看| 国产精品久久久久久久久免| 免费av观看视频| 伦理电影大哥的女人| 亚洲欧美日韩无卡精品| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 亚洲婷婷狠狠爱综合网| 在线亚洲精品国产二区图片欧美 | 免费黄网站久久成人精品| 尾随美女入室| 午夜日本视频在线| 国产精品人妻久久久影院| 嫩草影院入口| 婷婷色av中文字幕| 欧美变态另类bdsm刘玥| 亚洲国产高清在线一区二区三| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 午夜免费鲁丝| 最近的中文字幕免费完整| 亚洲av成人精品一区久久| 美女脱内裤让男人舔精品视频| 久久久成人免费电影| 国产精品爽爽va在线观看网站| 久久午夜福利片| 极品少妇高潮喷水抽搐| 亚洲成人av在线免费| 亚洲精品日韩av片在线观看| 亚洲国产欧美在线一区| 能在线免费看毛片的网站| 夜夜爽夜夜爽视频| av在线观看视频网站免费| 少妇猛男粗大的猛烈进出视频 | 插逼视频在线观看| 免费电影在线观看免费观看| 性色avwww在线观看| 国产精品伦人一区二区| 久久精品熟女亚洲av麻豆精品| 国内精品美女久久久久久| 欧美少妇被猛烈插入视频| 亚洲欧美一区二区三区国产| 国产精品麻豆人妻色哟哟久久| 日韩,欧美,国产一区二区三区| 夫妻性生交免费视频一级片| 日韩伦理黄色片| 亚洲高清免费不卡视频| av线在线观看网站| 亚洲电影在线观看av| 精品国产三级普通话版| 亚洲aⅴ乱码一区二区在线播放| 亚洲av在线观看美女高潮| 成人亚洲精品一区在线观看 | 美女被艹到高潮喷水动态| 中文天堂在线官网| 久久久久久国产a免费观看| 国产淫片久久久久久久久| 另类亚洲欧美激情| 日韩精品有码人妻一区| 啦啦啦中文免费视频观看日本| 欧美成人午夜免费资源| 黄色配什么色好看| 最新中文字幕久久久久| 成年人午夜在线观看视频| 亚洲,欧美,日韩| 亚洲av成人精品一二三区| 在现免费观看毛片| 色5月婷婷丁香| 亚洲自拍偷在线| 亚洲av二区三区四区| 精品国产露脸久久av麻豆| 99九九线精品视频在线观看视频| 精品亚洲乱码少妇综合久久| 汤姆久久久久久久影院中文字幕| 亚洲av免费在线观看| 亚洲欧美日韩另类电影网站 | 成年人午夜在线观看视频| 人妻制服诱惑在线中文字幕| 国产av不卡久久| 中文字幕免费在线视频6| 日本免费在线观看一区| 丰满乱子伦码专区| 亚洲精品中文字幕在线视频 | 欧美xxⅹ黑人| 久久久久久九九精品二区国产| 久久影院123| 晚上一个人看的免费电影| 白带黄色成豆腐渣| 18禁裸乳无遮挡动漫免费视频 | 国产高清不卡午夜福利| 男女国产视频网站| 极品少妇高潮喷水抽搐| 边亲边吃奶的免费视频| 一级a做视频免费观看| 色视频www国产| 大香蕉久久网| 日韩大片免费观看网站| 午夜爱爱视频在线播放| av在线老鸭窝| 亚洲国产日韩一区二区| 日韩伦理黄色片| 国产一区二区亚洲精品在线观看| 国产人妻一区二区三区在| 欧美老熟妇乱子伦牲交| 国产 一区精品| 国精品久久久久久国模美| 别揉我奶头 嗯啊视频| 边亲边吃奶的免费视频| 精品少妇久久久久久888优播| 丝瓜视频免费看黄片| av在线蜜桃| 久久精品久久精品一区二区三区| 国产综合懂色| 国产精品国产三级专区第一集| 天美传媒精品一区二区| 免费黄网站久久成人精品| 国产成人精品久久久久久| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 久久热精品热| av在线播放精品| 看黄色毛片网站| 3wmmmm亚洲av在线观看| 国语对白做爰xxxⅹ性视频网站| 可以在线观看毛片的网站| 国产黄色免费在线视频| 嫩草影院精品99| 国产精品蜜桃在线观看| 伦理电影大哥的女人| av又黄又爽大尺度在线免费看| 国产一级毛片在线| 国产有黄有色有爽视频| 91久久精品国产一区二区三区| 九九爱精品视频在线观看| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 欧美一级a爱片免费观看看| 日本与韩国留学比较| 色婷婷久久久亚洲欧美| 成年版毛片免费区| 国产精品人妻久久久影院| 成人亚洲精品av一区二区| 一级毛片黄色毛片免费观看视频| 欧美少妇被猛烈插入视频| 晚上一个人看的免费电影| av天堂中文字幕网| 一级a做视频免费观看| 亚洲成人一二三区av| 99视频精品全部免费 在线| 亚洲内射少妇av| 女的被弄到高潮叫床怎么办| 三级国产精品欧美在线观看| 日韩成人av中文字幕在线观看| 成人美女网站在线观看视频| 自拍偷自拍亚洲精品老妇| 禁无遮挡网站| 亚洲成色77777| 亚洲av日韩在线播放| 偷拍熟女少妇极品色| 国产色爽女视频免费观看| 18+在线观看网站| 日韩不卡一区二区三区视频在线| 97热精品久久久久久| 亚洲国产精品999| 亚洲精品色激情综合| 干丝袜人妻中文字幕| 日韩人妻高清精品专区| 性色avwww在线观看| 国产午夜精品一二区理论片| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 波野结衣二区三区在线| 国产亚洲最大av| 高清毛片免费看| 国内少妇人妻偷人精品xxx网站| 国产成人免费无遮挡视频| 高清在线视频一区二区三区| 亚洲精品亚洲一区二区| 天天躁日日操中文字幕| 日日摸夜夜添夜夜添av毛片| 成人国产麻豆网| 亚洲精品456在线播放app| 国产精品一区二区性色av| 一区二区三区免费毛片| 亚洲无线观看免费| 国产亚洲午夜精品一区二区久久 | 内射极品少妇av片p| 久久人人爽人人爽人人片va| 国产精品伦人一区二区| 我的老师免费观看完整版| av国产精品久久久久影院| kizo精华| 综合色av麻豆| 日韩视频在线欧美| 3wmmmm亚洲av在线观看| 成年免费大片在线观看| 午夜福利高清视频| 男人爽女人下面视频在线观看| 九九在线视频观看精品| 国产成人a∨麻豆精品| 免费大片黄手机在线观看| 欧美一区二区亚洲| av专区在线播放| 高清日韩中文字幕在线| 在线亚洲精品国产二区图片欧美 | 成年女人在线观看亚洲视频 | 久久精品国产亚洲网站| 99视频精品全部免费 在线| 欧美日韩视频高清一区二区三区二| 国产欧美日韩精品一区二区| 亚洲美女视频黄频| 成人二区视频| 亚洲av中文av极速乱| 在线 av 中文字幕| 亚洲av日韩在线播放| 18禁裸乳无遮挡免费网站照片| 日本熟妇午夜| 国产精品三级大全| 99热这里只有是精品50| 久久精品久久久久久噜噜老黄| 建设人人有责人人尽责人人享有的 | av免费观看日本| 国产精品一区www在线观看| 国产男女超爽视频在线观看| 午夜免费鲁丝| 成人午夜精彩视频在线观看| 国产精品一及| 又爽又黄a免费视频| 精品人妻视频免费看| 老司机影院成人| 国产精品一区二区三区四区免费观看| 青青草视频在线视频观看| 尾随美女入室| 青春草亚洲视频在线观看| 日韩一区二区视频免费看| 亚洲精品色激情综合| 乱码一卡2卡4卡精品| 国产精品秋霞免费鲁丝片| 欧美 日韩 精品 国产| 精品人妻视频免费看| 啦啦啦啦在线视频资源| 国产精品偷伦视频观看了| 人妻系列 视频| 亚洲精品亚洲一区二区| 狂野欧美激情性bbbbbb| 22中文网久久字幕| 色婷婷久久久亚洲欧美| 国产有黄有色有爽视频| 能在线免费看毛片的网站| 不卡视频在线观看欧美| 亚洲精品国产成人久久av| 久久久久久久亚洲中文字幕| 久久久久久久久大av| 中文字幕人妻熟人妻熟丝袜美| 日韩不卡一区二区三区视频在线| 久久久久久久大尺度免费视频| 亚洲成人精品中文字幕电影| 国产亚洲一区二区精品| 91aial.com中文字幕在线观看| 99热全是精品| 欧美精品一区二区大全| 国产精品熟女久久久久浪| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 亚洲第一区二区三区不卡| 美女cb高潮喷水在线观看| 人体艺术视频欧美日本| 亚洲va在线va天堂va国产| 大香蕉久久网| 在线观看免费高清a一片| 国产黄片美女视频| 久久韩国三级中文字幕| 国产白丝娇喘喷水9色精品| 欧美激情在线99| 特级一级黄色大片| 在线观看美女被高潮喷水网站| kizo精华| 青青草视频在线视频观看| 人人妻人人爽人人添夜夜欢视频 | 99久久精品国产国产毛片| 欧美三级亚洲精品| 2021少妇久久久久久久久久久| 午夜福利在线观看免费完整高清在| 在线 av 中文字幕| 亚洲av欧美aⅴ国产| 一本一本综合久久| 国产免费福利视频在线观看| 国产精品久久久久久精品电影小说 | 欧美精品人与动牲交sv欧美| 99久久精品国产国产毛片| 全区人妻精品视频| 亚洲综合精品二区| 国产av码专区亚洲av| 插阴视频在线观看视频| 国产高清不卡午夜福利| 性插视频无遮挡在线免费观看| 看免费成人av毛片| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| 欧美成人a在线观看| 免费观看在线日韩| 熟女人妻精品中文字幕| 日日啪夜夜撸| 女人被狂操c到高潮| 亚洲自拍偷在线| 国产精品久久久久久精品古装| 日本与韩国留学比较| 热re99久久精品国产66热6| 街头女战士在线观看网站| 色视频www国产| 精华霜和精华液先用哪个| 国产精品久久久久久久久免|