• <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在线观看| 91午夜精品亚洲一区二区三区| 久久99蜜桃精品久久| 国产男人的电影天堂91| 最后的刺客免费高清国语| 国产单亲对白刺激| 国产精品野战在线观看| 国产精品综合久久久久久久免费| 日韩大片免费观看网站 | 草草在线视频免费看| 丰满人妻一区二区三区视频av| 欧美激情在线99| 99热全是精品| 秋霞伦理黄片| 日韩一本色道免费dvd| 国产精品国产三级国产av玫瑰| 99久久无色码亚洲精品果冻| 麻豆av噜噜一区二区三区| 夜夜爽夜夜爽视频| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 亚洲欧美日韩卡通动漫| 欧美另类亚洲清纯唯美| 网址你懂的国产日韩在线| 亚洲精品日韩av片在线观看| 亚洲av熟女| 免费看av在线观看网站| 内射极品少妇av片p| 日本五十路高清| 国国产精品蜜臀av免费| 51国产日韩欧美| 成人综合一区亚洲| 国产爱豆传媒在线观看| 亚洲国产欧美人成| 国产乱来视频区| a级一级毛片免费在线观看| 91久久精品电影网| 国产视频首页在线观看| 麻豆精品久久久久久蜜桃| 啦啦啦韩国在线观看视频| av在线亚洲专区| 欧美日韩综合久久久久久| 男人舔奶头视频| 五月伊人婷婷丁香| av视频在线观看入口| 婷婷色av中文字幕| 欧美激情久久久久久爽电影| 日本wwww免费看| 一二三四中文在线观看免费高清| 亚洲图色成人| 亚洲电影在线观看av| av免费在线看不卡| 国产精品嫩草影院av在线观看| 久久久a久久爽久久v久久| 精品不卡国产一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | av在线天堂中文字幕| 最近的中文字幕免费完整| 一级毛片aaaaaa免费看小| av在线蜜桃| 精品午夜福利在线看| 国产亚洲精品久久久com| 日韩欧美国产在线观看| 精品不卡国产一区二区三区| 精品久久国产蜜桃| 成年免费大片在线观看| 网址你懂的国产日韩在线| 草草在线视频免费看| 好男人视频免费观看在线| 精品人妻一区二区三区麻豆| 国产av码专区亚洲av| 七月丁香在线播放| 国产精品熟女久久久久浪| 国产一区有黄有色的免费视频 | 蜜臀久久99精品久久宅男| 高清av免费在线| 亚洲国产色片| 亚洲精华国产精华液的使用体验| 中文欧美无线码| 超碰av人人做人人爽久久| 男人狂女人下面高潮的视频| 网址你懂的国产日韩在线| 亚洲国产精品专区欧美| 午夜亚洲福利在线播放| 日韩亚洲欧美综合| 狂野欧美白嫩少妇大欣赏| 建设人人有责人人尽责人人享有的 | 丝袜美腿在线中文| 欧美色视频一区免费| 国产精品熟女久久久久浪| 国内精品宾馆在线| 男插女下体视频免费在线播放| 成人无遮挡网站| 成人亚洲精品av一区二区| 免费av观看视频| 边亲边吃奶的免费视频| 亚洲欧美成人综合另类久久久 | 男人狂女人下面高潮的视频| 亚洲人成网站高清观看| 高清午夜精品一区二区三区| 男人舔女人下体高潮全视频| 国产av不卡久久| 一级毛片电影观看 | 美女cb高潮喷水在线观看| 91精品伊人久久大香线蕉| 免费搜索国产男女视频| 禁无遮挡网站| 免费看av在线观看网站| 国产高清三级在线| 国产精品99久久久久久久久| 亚洲第一区二区三区不卡| 精华霜和精华液先用哪个| 欧美精品国产亚洲| 成年免费大片在线观看| 一个人看的www免费观看视频| 搡女人真爽免费视频火全软件| 国产成人a∨麻豆精品| 国产大屁股一区二区在线视频| 精品熟女少妇av免费看| 中文字幕熟女人妻在线| 亚洲国产精品久久男人天堂| 午夜福利网站1000一区二区三区| 嫩草影院新地址| 蜜桃久久精品国产亚洲av| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看 | 日韩 亚洲 欧美在线| 免费av观看视频| 欧美区成人在线视频| 免费不卡的大黄色大毛片视频在线观看 | 国产视频内射| 亚洲成人av在线免费| 亚洲精品色激情综合| 久99久视频精品免费| 国产亚洲精品av在线| 看非洲黑人一级黄片| 欧美成人一区二区免费高清观看| 亚洲av中文字字幕乱码综合| 国产精品国产三级专区第一集| 欧美性感艳星| 亚洲人成网站高清观看| 国产一区二区三区av在线| 亚洲五月天丁香| 十八禁国产超污无遮挡网站| 村上凉子中文字幕在线| 亚洲最大成人手机在线| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 久久精品国产亚洲网站| 小说图片视频综合网站| 久久久精品94久久精品| 国产精品一二三区在线看| 免费无遮挡裸体视频| 久久99热这里只有精品18| 午夜精品在线福利| 一个人看视频在线观看www免费| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 亚洲欧美成人精品一区二区| 七月丁香在线播放| 又爽又黄无遮挡网站| 国产白丝娇喘喷水9色精品| 97超视频在线观看视频| 干丝袜人妻中文字幕| 国语自产精品视频在线第100页| 日韩三级伦理在线观看| 免费看av在线观看网站| 欧美精品一区二区大全| 91在线精品国自产拍蜜月| videos熟女内射| 日本免费一区二区三区高清不卡| 99久久无色码亚洲精品果冻| 国产激情偷乱视频一区二区| 只有这里有精品99| 午夜免费激情av| 国产伦精品一区二区三区四那| 日本色播在线视频| 在现免费观看毛片| 九草在线视频观看| 97在线视频观看| 蜜桃亚洲精品一区二区三区| 人妻夜夜爽99麻豆av| 少妇裸体淫交视频免费看高清| 色尼玛亚洲综合影院| 国产高潮美女av| 视频中文字幕在线观看| 欧美一级a爱片免费观看看| 亚洲自拍偷在线| 亚洲精品影视一区二区三区av| 午夜精品一区二区三区免费看| 蜜臀久久99精品久久宅男| 美女内射精品一级片tv| 在线播放无遮挡| 热99在线观看视频| 久久久久性生活片| 波野结衣二区三区在线| 一级二级三级毛片免费看| 国产精品永久免费网站| 只有这里有精品99| 性色avwww在线观看| 赤兔流量卡办理| 日本wwww免费看| 最近的中文字幕免费完整| 一个人看的www免费观看视频| 成年女人永久免费观看视频| 亚洲久久久久久中文字幕| 日本免费一区二区三区高清不卡| 久久久久九九精品影院| 国产伦一二天堂av在线观看| 美女黄网站色视频| 老司机福利观看| 国内精品美女久久久久久| 免费看美女性在线毛片视频| 夜夜爽夜夜爽视频| 久久99热这里只有精品18| 热99re8久久精品国产| av福利片在线观看| 国产高清国产精品国产三级 | or卡值多少钱| 国产日韩欧美在线精品| 黄片wwwwww| 国产片特级美女逼逼视频| 欧美人与善性xxx| 亚洲国产欧美在线一区| 精品久久久久久久久久久久久| 边亲边吃奶的免费视频| 99在线人妻在线中文字幕| 亚洲五月天丁香| 午夜福利成人在线免费观看| 精品久久久噜噜| 内射极品少妇av片p| 久久婷婷人人爽人人干人人爱| 最近中文字幕高清免费大全6| 简卡轻食公司| 国产在线一区二区三区精 | 欧美日韩综合久久久久久| 久久精品人妻少妇| 七月丁香在线播放| 少妇猛男粗大的猛烈进出视频 | 美女大奶头视频| 毛片女人毛片| 成人亚洲欧美一区二区av| 18禁在线播放成人免费| 国产午夜精品一二区理论片| 亚洲美女视频黄频| 极品教师在线视频| 国产av在哪里看| 日本欧美国产在线视频| 国产午夜精品论理片| 免费不卡的大黄色大毛片视频在线观看 | 国产淫片久久久久久久久| 高清在线视频一区二区三区 | 亚洲成人久久爱视频| 伦精品一区二区三区| 免费电影在线观看免费观看| 国产精品久久久久久久久免| 村上凉子中文字幕在线| 国产乱来视频区| 日本黄大片高清| 久久韩国三级中文字幕| av在线亚洲专区| 变态另类丝袜制服| 丝袜喷水一区| 亚洲av男天堂| 中文字幕制服av| 成人毛片60女人毛片免费| 精品久久久久久久末码| 黄色配什么色好看| 国产精品精品国产色婷婷| 精品人妻一区二区三区麻豆| 成人国产麻豆网| av女优亚洲男人天堂| 黄色日韩在线| 少妇裸体淫交视频免费看高清| 综合色丁香网| 日韩一区二区视频免费看| 亚洲一级一片aⅴ在线观看| 婷婷色av中文字幕| 蜜臀久久99精品久久宅男| 99在线视频只有这里精品首页| 亚洲成色77777| 久久久精品欧美日韩精品| 大又大粗又爽又黄少妇毛片口| 少妇人妻一区二区三区视频| 亚洲四区av| 亚洲最大成人中文| 成人漫画全彩无遮挡| 久99久视频精品免费| 九九在线视频观看精品| 亚洲欧美日韩卡通动漫| 亚洲国产精品合色在线| 亚洲av福利一区| 在线播放国产精品三级| 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 久久精品91蜜桃| 少妇猛男粗大的猛烈进出视频 | 精品久久久噜噜| 亚洲国产精品专区欧美| 国产精品永久免费网站| 九九久久精品国产亚洲av麻豆| 成年免费大片在线观看| 国产爱豆传媒在线观看| 亚洲欧美中文字幕日韩二区| 永久网站在线| 午夜免费激情av| 成年版毛片免费区| 午夜亚洲福利在线播放| 欧美+日韩+精品| 99在线人妻在线中文字幕| 久久人妻av系列| 国产不卡一卡二| 青春草国产在线视频| 午夜福利在线观看免费完整高清在| 亚洲av熟女| 伦精品一区二区三区| 激情 狠狠 欧美| 精品国内亚洲2022精品成人| 久久精品影院6| 啦啦啦啦在线视频资源| 成年版毛片免费区| 天天躁日日操中文字幕| 成人高潮视频无遮挡免费网站| 深爱激情五月婷婷| 天堂av国产一区二区熟女人妻| 久久99热这里只有精品18| 女的被弄到高潮叫床怎么办| 看非洲黑人一级黄片| 亚洲精品456在线播放app| 久久久亚洲精品成人影院| 日韩欧美精品免费久久| 欧美最新免费一区二区三区| 最后的刺客免费高清国语| 少妇高潮的动态图| 免费搜索国产男女视频| 老司机影院毛片| 亚洲国产精品成人综合色| 搞女人的毛片| 搡老妇女老女人老熟妇| 搞女人的毛片| 最近中文字幕高清免费大全6| 高清视频免费观看一区二区 | 99热全是精品| 国产亚洲5aaaaa淫片| 非洲黑人性xxxx精品又粗又长| 日韩人妻高清精品专区| 1024手机看黄色片| 日韩欧美 国产精品| 美女cb高潮喷水在线观看| 一级黄色大片毛片| 美女内射精品一级片tv| 99视频精品全部免费 在线| 一个人观看的视频www高清免费观看| 日本熟妇午夜| 国产成人a∨麻豆精品| 国产毛片a区久久久久| 一级毛片久久久久久久久女| 天天一区二区日本电影三级| 天美传媒精品一区二区| 伦精品一区二区三区| 偷拍熟女少妇极品色| kizo精华| 亚洲av电影在线观看一区二区三区 | 久99久视频精品免费| 边亲边吃奶的免费视频| 狂野欧美白嫩少妇大欣赏| 国产乱人视频| 1024手机看黄色片| 少妇高潮的动态图| 亚洲国产欧美人成| 免费黄网站久久成人精品| 午夜福利网站1000一区二区三区| 欧美一区二区亚洲| 欧美一级a爱片免费观看看| 久久国产乱子免费精品| 亚洲高清免费不卡视频| 可以在线观看毛片的网站| 少妇丰满av| 国产黄色视频一区二区在线观看 | 欧美性猛交╳xxx乱大交人| 欧美高清性xxxxhd video| 午夜福利网站1000一区二区三区| 一二三四中文在线观看免费高清| 美女被艹到高潮喷水动态| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 国产精品国产高清国产av| 中文精品一卡2卡3卡4更新| 国产精品伦人一区二区| 日本免费在线观看一区| 最近中文字幕高清免费大全6| 国产探花在线观看一区二区| 久久这里只有精品中国| 久久精品国产亚洲av天美| 免费av毛片视频| 蜜桃久久精品国产亚洲av| 国产伦精品一区二区三区视频9| 99热这里只有精品一区| 亚洲综合色惰| 精品久久久噜噜| 天堂√8在线中文| 午夜视频国产福利| 国产视频首页在线观看| 欧美激情国产日韩精品一区| 毛片一级片免费看久久久久| 精品酒店卫生间| 欧美性感艳星| 一级毛片aaaaaa免费看小| av国产免费在线观看| 欧美xxxx性猛交bbbb| 免费无遮挡裸体视频| 国产黄a三级三级三级人| 中文字幕免费在线视频6| 成年女人看的毛片在线观看| 国产黄色小视频在线观看| 国产 一区 欧美 日韩| 久久欧美精品欧美久久欧美| 一区二区三区免费毛片| 联通29元200g的流量卡| 日韩欧美精品v在线| 亚洲最大成人手机在线| 中文资源天堂在线| 久久久久久久久大av| 高清毛片免费看| 少妇丰满av| 午夜福利在线观看吧| 国产成人午夜福利电影在线观看| 久久久成人免费电影| 99热全是精品| 熟女电影av网| 日韩成人av中文字幕在线观看| 久久久精品大字幕| 99久久精品热视频| 久久精品国产99精品国产亚洲性色| 久久热精品热| 禁无遮挡网站| 国产精品一区二区在线观看99 | 日韩精品有码人妻一区| 色视频www国产| 国产高潮美女av| 一级毛片我不卡| 国产毛片a区久久久久| 久久久精品欧美日韩精品| 欧美97在线视频| 久久久成人免费电影| 乱人视频在线观看| 欧美+日韩+精品| 一二三四中文在线观看免费高清| 高清午夜精品一区二区三区| 男女啪啪激烈高潮av片| 中文乱码字字幕精品一区二区三区 | 成年女人看的毛片在线观看| 97热精品久久久久久| 欧美日韩综合久久久久久| 国产亚洲精品久久久com| 国产精品1区2区在线观看.| 校园人妻丝袜中文字幕| 色网站视频免费| 精品欧美国产一区二区三| 久久久午夜欧美精品| 久热久热在线精品观看| 99久久精品热视频| 美女被艹到高潮喷水动态| 美女xxoo啪啪120秒动态图| 三级经典国产精品| 老司机影院成人| 国产三级中文精品| 国产精品蜜桃在线观看| 丰满人妻一区二区三区视频av| 一级毛片aaaaaa免费看小| 久久综合国产亚洲精品| 亚洲成人av在线免费| 在现免费观看毛片| 亚洲成人精品中文字幕电影| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 国产精品人妻久久久久久| 狠狠狠狠99中文字幕| 在线观看66精品国产| 亚洲欧美中文字幕日韩二区| 国产一区二区三区av在线| 国产乱人偷精品视频| 天堂√8在线中文| 欧美日韩国产亚洲二区| 亚洲五月天丁香| 欧美一区二区精品小视频在线| 美女高潮的动态| 久久久久久伊人网av| 免费一级毛片在线播放高清视频| 国产高清有码在线观看视频| 日本免费在线观看一区| 久久久久性生活片| 久久婷婷人人爽人人干人人爱| 黄色欧美视频在线观看| 国产黄色小视频在线观看| 最近中文字幕2019免费版| 熟女人妻精品中文字幕| 少妇丰满av| 国产探花在线观看一区二区| 观看美女的网站| 夜夜看夜夜爽夜夜摸| 在线免费观看不下载黄p国产| 欧美bdsm另类| 成人美女网站在线观看视频| 国产精品女同一区二区软件| 全区人妻精品视频| 99热6这里只有精品| 午夜a级毛片| 中文字幕久久专区| 免费看美女性在线毛片视频| 亚洲中文字幕日韩| 国产69精品久久久久777片| 免费av毛片视频| 18+在线观看网站| 成人午夜高清在线视频| 国产一级毛片七仙女欲春2| 色视频www国产| 久久久久久伊人网av| 国产成人午夜福利电影在线观看| 性插视频无遮挡在线免费观看| 边亲边吃奶的免费视频| 欧美另类亚洲清纯唯美| 床上黄色一级片| 欧美+日韩+精品| 99久久精品国产国产毛片| av在线观看视频网站免费| 久久久久久久久久久丰满| 色哟哟·www| 亚洲欧美日韩东京热| 国产亚洲av嫩草精品影院| 国产极品天堂在线| 亚洲国产精品久久男人天堂| 久久久久久大精品| 国产高清不卡午夜福利| 建设人人有责人人尽责人人享有的 | 亚洲aⅴ乱码一区二区在线播放| 亚洲av日韩在线播放| 欧美丝袜亚洲另类| 国产伦精品一区二区三区视频9| 麻豆一二三区av精品| 欧美成人a在线观看| 成人毛片a级毛片在线播放| 日本黄色片子视频| 免费不卡的大黄色大毛片视频在线观看 | 大香蕉久久网| 日韩亚洲欧美综合| 久久精品国产自在天天线| 国产淫语在线视频| 卡戴珊不雅视频在线播放| 日韩一本色道免费dvd| 成年av动漫网址| 一级毛片电影观看 | 一区二区三区四区激情视频| 久久精品91蜜桃| 久热久热在线精品观看| 天堂√8在线中文| 国产精品.久久久| 人妻少妇偷人精品九色| 草草在线视频免费看| 亚洲成av人片在线播放无| 午夜老司机福利剧场| 村上凉子中文字幕在线| 久久久精品大字幕| 日日摸夜夜添夜夜爱| 欧美一区二区精品小视频在线| 精品久久久久久成人av| 国产亚洲最大av| 免费av不卡在线播放| 国产亚洲一区二区精品| 亚洲国产欧美在线一区| 人妻制服诱惑在线中文字幕| 水蜜桃什么品种好| 亚洲av男天堂| 免费看光身美女| 国产精品嫩草影院av在线观看| 男人舔女人下体高潮全视频| 小蜜桃在线观看免费完整版高清| 插阴视频在线观看视频| 成人三级黄色视频| 国产激情偷乱视频一区二区| 日韩欧美精品免费久久| 深夜a级毛片| 国产精品电影一区二区三区| 国产精品久久视频播放| 国产亚洲5aaaaa淫片| 国产一区有黄有色的免费视频 | 99在线视频只有这里精品首页| 男人舔女人下体高潮全视频| 日韩欧美在线乱码| 国内精品宾馆在线| 免费播放大片免费观看视频在线观看 | 日韩人妻高清精品专区| 国产在视频线在精品| 18禁在线播放成人免费| 边亲边吃奶的免费视频| 久久精品国产鲁丝片午夜精品| 青青草视频在线视频观看| 日本黄大片高清| 亚洲精品乱码久久久v下载方式| 最近视频中文字幕2019在线8| 精品国产露脸久久av麻豆 | 国产成人精品一,二区| 丝袜美腿在线中文| 2022亚洲国产成人精品| 国产日韩欧美在线精品| 男女视频在线观看网站免费| 久久久午夜欧美精品| 亚洲在线观看片| 免费搜索国产男女视频| 自拍偷自拍亚洲精品老妇| 亚洲美女搞黄在线观看| 精品午夜福利在线看| 尾随美女入室| 卡戴珊不雅视频在线播放| 欧美三级亚洲精品|