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

    基于高光譜的干旱區(qū)鹽漬化土壤鹽分含量估算

    2018-12-19 08:52:06如則麥麥提米吉提買(mǎi)買(mǎi)提沙吾提麥爾耶姆亞森馬春玥
    江蘇農(nóng)業(yè)科學(xué) 2018年22期
    關(guān)鍵詞:鹽漬化含鹽量綠洲

    如則麥麥提·米吉提, 買(mǎi)買(mǎi)提·沙吾提, 麥爾耶姆·亞森, 馬春玥

    (新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院/新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,新疆烏魯木齊 830046)

    土壤鹽漬化是土地退化的主要表現(xiàn),灌溉引起的土壤次生鹽漬化成為限制干旱區(qū)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的因素之一,也是影響綠洲生態(tài)環(huán)境穩(wěn)定與安全的重要因素[1]。鹽漬土是重要的耕地后備資源,及時(shí)獲取土壤鹽漬化信息,對(duì)于我國(guó)西北部干旱區(qū)的可持續(xù)發(fā)展至關(guān)重要[2-3]。遙感技術(shù)能宏觀(guān)、快速地提取地表時(shí)空信息,已成為動(dòng)態(tài)監(jiān)測(cè)土壤鹽漬化的新方法[4]。鹽漬化土壤的光譜反射特征是由土壤理化性質(zhì)決定的,傳統(tǒng)遙感數(shù)據(jù)并不能完全反映土壤的光譜特征,而高光譜數(shù)據(jù)能夠提供鹽漬化土壤細(xì)微的光譜特征。光譜分辨率高、波段具有連續(xù)性強(qiáng)、光譜信息量大的優(yōu)點(diǎn)。闡明土壤屬性與其高光譜的關(guān)系,成為高光譜遙感技術(shù)定量監(jiān)測(cè)和提取土壤信息的重要方法,為實(shí)現(xiàn)大尺度、實(shí)時(shí)監(jiān)測(cè)土壤提供了新途徑,使定量反演土壤鹽漬化狀況成為可能[5-12]。目前,國(guó)內(nèi)外己有眾多學(xué)者對(duì)土壤鹽漬化的高光譜特征及定量、半定量做了大量的研究工作,Dehaan等觀(guān)測(cè)澳大利亞Murray-Darling盆地的鹽殼及重度、中度、輕度鹽漬土的光譜,發(fā)現(xiàn)這4類(lèi)土壤在不同光譜波段有明顯的吸收作用[13];Leone等應(yīng)用實(shí)驗(yàn)室光譜進(jìn)行多元統(tǒng)計(jì)分析以評(píng)價(jià)土壤開(kāi)發(fā)和土壤退化程度并建立預(yù)測(cè)模型[14-16]。劉慶生等分析了鹽分含量與上覆植被光譜的關(guān)系[17]。馬諾等通過(guò)野外調(diào)查以及實(shí)地測(cè)點(diǎn),利用便攜式光譜儀測(cè)量研究區(qū)域內(nèi)不同鹽漬程度的土壤光譜,分析土壤光譜曲線(xiàn)與土壤鹽漬程度之間的關(guān)系[18]。劉亞秋等利用熱紅外光譜儀建立了鹽分與鹽漬化土壤發(fā)射率一階導(dǎo)數(shù)模型[19]。彭杰等對(duì)野外高光譜數(shù)據(jù)與土壤含鹽量進(jìn)行耦合分析并建立多元線(xiàn)性回歸建模,決定系數(shù)達(dá)到0.80[20]。丁建麗等利用高光譜數(shù)據(jù)建立的監(jiān)測(cè)模型為提取鹽漬化土壤信息提供了依據(jù)[21]。

    基于此,本研究以新疆維吾爾自治區(qū)渭干河-庫(kù)車(chē)河三角洲綠洲為研究區(qū),利用57個(gè)表層土壤樣品的室內(nèi)高光譜數(shù)據(jù),并結(jié)合實(shí)測(cè)土壤含鹽量,研究15種數(shù)學(xué)變換處理在高光譜數(shù)據(jù)估算荒漠土壤含鹽量的可能性,并應(yīng)用偏最小二乘回歸、主成分回歸、多元逐步線(xiàn)性回歸建立土壤含鹽量估算模型,以期為干旱區(qū)綠洲土壤鹽漬化的高精度遙感監(jiān)測(cè)提供一種更為有效的手段。

    1 材料與方法

    1.1 研究區(qū)概況

    胃干河-庫(kù)車(chē)河流域簡(jiǎn)稱(chēng)胃-庫(kù)綠洲,是位于我國(guó)西北部的典型干旱區(qū)綠洲。本地區(qū)行政上轄阿克蘇地區(qū)的庫(kù)車(chē)縣、沙雅縣、新河縣,地理坐標(biāo)分別為80°37′~83°59′E,41°06′~41°40′N(xiāo)。綠洲北部跟天山南部接壤,南部延伸到塔克拉瑪干沙漠北緣。胃-庫(kù)綠洲是胃干河-庫(kù)車(chē)河長(zhǎng)期沉淀形成的山前沖洪積扇平原,是塔里木盆地的主要綠洲及我國(guó)重要的產(chǎn)棉區(qū)之一。該綠洲位于中緯度地區(qū),遠(yuǎn)離海洋,氣候干燥,風(fēng)沙頻繁,降水稀少,蒸發(fā)量大,晝夜溫差大,夏季干熱,冬季干冷。年均氣溫10.5~11.4 ℃,年降水50.5~66.5 mm,年均蒸發(fā)量2 000.7~2 092.0 mm,年均無(wú)霜期 240 d。土壤以輕壤和沙壤為最多,中壤、重壤、黏土有少量分布,土壤構(gòu)成物顆粒細(xì),透水性差。

    1.2 數(shù)據(jù)采集

    本研究于2016年4月中旬,在渭干河-庫(kù)車(chē)河綠洲進(jìn)行野外調(diào)查,采集57個(gè)采樣點(diǎn)的土壤樣品,采樣點(diǎn)分布見(jiàn)圖1。每個(gè)點(diǎn)分3層取樣(0~10 cm、10~30 cm、30~50 cm),共采集171個(gè)土樣。樣品在室內(nèi)條件下自然風(fēng)干,磨碎后過(guò)2 mm篩子。其中,每個(gè)樣品分為2份,1份樣品用于測(cè)定土壤,鹽分、溶解性固體(TDS)、pH值等指標(biāo),另1份用于土壤光譜的測(cè)定。在土水體積比為1 ∶5的土壤懸濁液中,使用便攜式多參數(shù)分析儀Multi 3420 SETB(wissens chaftlich technische werkst?tten,德國(guó))進(jìn)行指標(biāo)測(cè)定。

    測(cè)出的土壤參數(shù)值統(tǒng)計(jì)特征見(jiàn)表1。本研究區(qū)含鹽量最大值為69.8 g/kg,最小值為0 g/kg,平均值為20.19 g/kg,變異系數(shù)為99.15%,屬?gòu)?qiáng)變異,表明土壤鹽分分布很不均勻。土壤樣本的pH值變異系數(shù)為8.26%,屬低變異。根據(jù)土壤的分級(jí)標(biāo)準(zhǔn),57個(gè)樣本中,非鹽漬土樣本為13個(gè),輕、中、重度鹽漬化樣本數(shù)分別為11、10、23個(gè)。非鹽漬化土壤樣本主要以綠洲內(nèi)部為主,內(nèi)部到綠洲與荒漠邊界土壤含鹽量逐漸增加。中、重度鹽漬化土壤分布在研究區(qū)的東南部區(qū)域。

    表1 土壤指標(biāo)統(tǒng)計(jì)特征

    注:樣本數(shù)=171個(gè)。

    表2 土壤鹽漬化程度分級(jí)標(biāo)準(zhǔn)[22]

    1.3 土壤樣本光譜測(cè)定

    使用FieldSpec3型光譜儀(analytical spectral devices),波段范圍為350~2 500 nm,光譜采樣間隔分別為 1.4 nm(波段350~1 000 nm區(qū)間)和2 nm(波段1 000~2 500 nm區(qū)間),采樣間隔為1 nm,輸出波段數(shù)為2 151個(gè)。土壤樣本光譜測(cè)定在能控制光照條件的暗室內(nèi)進(jìn)行,測(cè)試之前都要進(jìn)行白板校正,儀器光源為50 W鹵素?zé)?,天頂角?5°,8°視場(chǎng)角探頭垂直土樣表面,距離為10 cm。每個(gè)土樣分別測(cè)10次,并取其算術(shù)平均值作為該樣品的實(shí)際光譜反射率。

    1.4 數(shù)據(jù)預(yù)處理與變換

    將57個(gè)采樣點(diǎn)的表層土壤含鹽量從高到低進(jìn)行排序,等間隔選取37個(gè)作為建模集,20個(gè)作為驗(yàn)證集,分別用于模型的建立以及精度驗(yàn)證。將通過(guò)顯著性檢驗(yàn)(P<0.01)的波段作為自變量,利用偏最小二乘回歸法(Partial Least-Squares Regression,PLSR)、主成分回歸模型(Principle Component Regression model,PCM)和多元逐步回歸模型(Multiple Stepwise Regression Model,MSRM)建立土壤含鹽量的估算模型,通過(guò)對(duì)比各模型的決定系數(shù)(Determination of Coefficients,R2)、均方根誤差(Root Mean Square Error of Prediction,RMSE)、相對(duì)分析誤差(Relative Prediction Deviation,RPD)篩選出最優(yōu)模型。R2用于判定模型的穩(wěn)定程度,取值范圍為 0~1,其值越大說(shuō)明該模型的穩(wěn)定性越好;RMSE用于表征模型的準(zhǔn)確性,其值越小,表明模型的精度越高。另外,當(dāng)RPD的值小于1.4時(shí),模型基本上可用;當(dāng)RPD的值在1.4~2.0之間時(shí),模型估算精度一般,RPD值越大于2.0,表明模型的定量預(yù)測(cè)能力越好[25]。

    2 結(jié)果與分析

    2.1 土壤樣本反射光譜曲線(xiàn)特征

    將57個(gè)土壤表層樣本反射率以土壤鹽漬化程度分級(jí)標(biāo)準(zhǔn)進(jìn)行分類(lèi)再求平均,作為不同程度鹽漬化土壤樣本的光譜反射率曲線(xiàn)。從圖2可以看出,4條曲線(xiàn)變化趨勢(shì)大致相似,土壤光譜從350~1 300 nm內(nèi)波段范圍持續(xù)上升,至 1 400 nm 處出現(xiàn)一個(gè)淺水分吸收谷,之后從1 400~1 850 nm波段范圍內(nèi)繼續(xù)升高,但增加幅度較為平緩,至1 900 nm處再次出現(xiàn)一個(gè)深的水分吸收谷。土壤成分含量值在一定程度上決定了土壤光譜曲線(xiàn)反射率值的大小和曲線(xiàn)整體的走勢(shì),含鹽量越高,反射率也高。

    2.2 土壤含鹽量與光譜相關(guān)性分析

    將土壤含鹽量與光譜反射率及其15種變換形式做相關(guān)系數(shù)在0.01水平上的顯著性檢驗(yàn)及相關(guān)性分析。從圖3可以看出,從土壤含鹽量與光譜的相關(guān)性分析來(lái)看,反射率一階微分(圖3-b)、反射率二階微分(圖3-c)、對(duì)數(shù)相反數(shù)一階微分(圖3-i)、連續(xù)統(tǒng)去除一階微分(圖3-o)、平方根一階微分(圖3-m)、立方根一階微分(圖3-k)與土壤含鹽量數(shù)據(jù)的整體相關(guān)系數(shù)較好,并通過(guò)顯著性檢驗(yàn)(P<0.01)。土壤含鹽量與R、1/R、(1/R)′、lgR、1/lgR、連續(xù)統(tǒng)去除法、反射率平方根、立方根的相關(guān)性較小,沒(méi)有通過(guò)0.01水平顯著性檢驗(yàn)的波段。因此,從光譜反射率及其15種變換形式中,遴選出以上5種,只對(duì)以上5種光譜反射率變換結(jié)果分別與土壤含鹽量進(jìn)行估算模型的建立與驗(yàn)證。

    結(jié)果表明,單個(gè)波段所建立的模型很難預(yù)測(cè)土壤含鹽量,微分變換較好地體現(xiàn)了原始數(shù)據(jù)中隱晦的光譜信息,相關(guān)性得以提升,據(jù)此運(yùn)用光低階微分處理方法,可以篩選出對(duì)土壤含鹽量敏感的特征光譜波段[26]。光譜反射率不同變換形式與土壤含鹽量特征波段的選取是選擇相關(guān)系數(shù)的峰值所對(duì)應(yīng)的光譜波段。

    2.3 土壤含鹽量建模

    根據(jù)以上相關(guān)分析結(jié)果得到的特征光譜波段,針對(duì)本研究的高光譜指數(shù)分別建立土壤含鹽量的偏最小二乘回歸模型,主成分回歸模型和多元逐步線(xiàn)性回歸模型。DPS數(shù)據(jù)處理軟件將57個(gè)土壤樣本中隨機(jī)選取的37個(gè)樣本用于建立模型,以土壤含鹽量作為因變量,以特征波段的光譜反射率作為自變量,分別建立土壤含鹽量的偏最小二乘回歸模型、主成分回歸模型和多元逐步線(xiàn)性回歸模型。運(yùn)用偏最小二乘回歸模型見(jiàn)表3,主成分回歸模型和多元逐步線(xiàn)性回歸方法的敏感光譜波段與土壤含鹽量之間的預(yù)測(cè)模型。從表3可以看出,以原始光譜與土壤含鹽量相關(guān)性較差,相關(guān)性最大的不到 0.2,沒(méi)有超過(guò)0.01顯著性檢驗(yàn)線(xiàn),因此建立的模型擬合精度很差。經(jīng)過(guò)微分變換后,土壤含鹽量與高光譜的相關(guān)性有所提高,最大可達(dá)到0.47,為挖掘敏感波段呈現(xiàn)良好的效果,建立的模型估算精度較好。所選擇的5種數(shù)學(xué)變換形式當(dāng)中,基于連續(xù)統(tǒng)去除一階微分模型的決定系數(shù)最差,立方根一階微分的最好。根據(jù)決定系數(shù)R2和相對(duì)分析誤差RPD值最大、均方根誤差最小的原則,對(duì)于各種變換形式,3種預(yù)測(cè)模型都是土壤含鹽量立方根一階微分光譜特征波段為自變量時(shí)擬合最優(yōu),R2值分別為0.68、0.65、0.85。利用多元逐步線(xiàn)性回歸方法建立的模型的精度最高,模型的R2、RPD和RMSE分別達(dá)到0.85、2.36、1.8。可以使用這種預(yù)測(cè)模型來(lái)定量反演土壤含鹽量,提取土壤鹽漬化信息。

    表3 含鹽量回歸模型

    2.4 模型精度檢驗(yàn)

    根據(jù)模型的精度參數(shù),對(duì)建立的模型進(jìn)行篩選后發(fā)現(xiàn),RPD≥2的模型只有2個(gè),是基于平方根一階微分和立方根一階微分的多元逐步線(xiàn)性回歸模型建立的模型,其他模型的RPD屬于小于或等于1.5。對(duì)模型的R2、RMSE、RPD以及波段數(shù)量進(jìn)行綜合對(duì)比發(fā)現(xiàn),基于平方根一階微分和立方根一階微分建立多元逐步回歸模型RPD、RMSE、R2最優(yōu)且參與建模的波段數(shù)量較多(19個(gè)和16個(gè))。一階和二階微分模型對(duì)驗(yàn)證樣品的含鹽量決定系數(shù)分別為0.34和 0.57,偏離1 ∶1線(xiàn),模型估算精度不高。連續(xù)統(tǒng)去除一階微分模型對(duì)驗(yàn)證樣品的含鹽量決定系數(shù)最差,只達(dá)到0.08。對(duì)3種建模方法而言,基于多元逐步線(xiàn)性回歸方法的模型精度高于其他2個(gè)方法的建模精度;利用主成分回歸方法建立的模型精度最低,R2值最高到0.68。對(duì)原數(shù)據(jù)的數(shù)學(xué)變換后整數(shù)階微分而言,一階微分處理增強(qiáng)了原數(shù)據(jù)的敏感性及建模能力。說(shuō)明原始數(shù)據(jù)的平方根一階微分和立方根一階微分一定程度上可以增強(qiáng)光譜對(duì)含鹽量的敏感程度,并利用多元逐步線(xiàn)性回歸方法對(duì)估算土壤鹽含量具有較好的定量反演能力,且計(jì)算量相對(duì)較小,因而確定該模型為最優(yōu)模型(圖4)。

    3 結(jié)論

    在本研究中,通過(guò)對(duì)土壤含鹽量與高光譜指數(shù)之間的相關(guān)性分析、3種預(yù)測(cè)模型的建立以及模型精度驗(yàn)證,可以得出以下結(jié)論。

    在不同鹽漬化水平下,土壤的光譜特征在形態(tài)上趨于一致;在整個(gè)波譜區(qū)間,土壤的光譜反射率隨鹽漬化程度的加重呈現(xiàn)增加趨勢(shì)。

    對(duì)原始高光譜數(shù)據(jù)的15種數(shù)學(xué)變換中,利用微分處理是體現(xiàn)光譜數(shù)據(jù)敏感波段效率好,經(jīng)微分變換后的高光譜指數(shù)與土壤含鹽量進(jìn)行相關(guān)性分析,篩選出對(duì)土壤含鹽量變化響應(yīng)的敏感波段,構(gòu)建基于偏最小二乘回歸、主成分回歸和多元逐步線(xiàn)性回歸的土壤鹽分監(jiān)測(cè)模型。結(jié)果表明,立方根一階微分變換形式為構(gòu)建預(yù)測(cè)模型的最佳光譜指標(biāo)。

    對(duì)土壤含鹽量的3種回歸模型預(yù)測(cè)效果進(jìn)行檢驗(yàn),比較實(shí)測(cè)值與預(yù)測(cè)值的相關(guān)關(guān)系。結(jié)果表明,多元逐步線(xiàn)性回歸模型的穩(wěn)定性和預(yù)測(cè)精度高于其他2種回歸模型,證明了該方法的普適性。

    本研究運(yùn)用偏最小二乘回歸、主成分回歸和多元逐步線(xiàn)性回歸模型在實(shí)驗(yàn)室內(nèi)對(duì)干旱區(qū)土壤鹽分含量的測(cè)定,對(duì)比3種回歸方法所建立模型的精度,發(fā)現(xiàn)多元逐步線(xiàn)性回歸模型的預(yù)測(cè)精度最高。但影響土壤光譜反射率是土壤的粒徑、粗糙度和水分含量等多種因素的綜合反映。在今后的研究中要考慮這些因素的影響。

    猜你喜歡
    鹽漬化含鹽量綠洲
    蔬菜大棚土壤鹽漬化成因及防治措施
    含鹽量及含水率對(duì)鹽漬土凍脹規(guī)律影響試驗(yàn)研究*
    土地質(zhì)量地球化學(xué)調(diào)查成果在判定土壤鹽漬化、沙化中的應(yīng)用
    綠洲里的老先生
    黃河三角洲鹽漬土有機(jī)氮組成及氮有效性對(duì)土壤含鹽量的響應(yīng)*
    甘肅蘇干湖濕地土壤鹽漬化、地下水位埋深及其對(duì)生態(tài)環(huán)境的影響
    沙漠綠洲
    什么是水的含鹽量?
    瑪納斯河流域土壤鹽漬化時(shí)空動(dòng)態(tài)變化
    秦陵陪葬坑土遺址安全含鹽量探究
    精品国内亚洲2022精品成人| 69人妻影院| 在现免费观看毛片| 又爽又黄a免费视频| 久久精品国产亚洲av涩爱 | 99久久精品国产亚洲精品| 亚洲欧美日韩高清在线视频| 亚洲第一区二区三区不卡| 草草在线视频免费看| 国产伦人伦偷精品视频| 色av中文字幕| 欧美精品啪啪一区二区三区| 国产真实伦视频高清在线观看 | 三级毛片av免费| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 婷婷丁香在线五月| 国产精品一区二区三区四区免费观看 | 高清在线国产一区| 亚洲中文字幕日韩| 99热只有精品国产| 精品一区二区三区视频在线| av专区在线播放| 一区二区三区激情视频| 国产国拍精品亚洲av在线观看| 日韩高清综合在线| 久久久精品大字幕| 中文字幕av成人在线电影| 国产国拍精品亚洲av在线观看| 亚洲最大成人手机在线| 亚洲三级黄色毛片| 欧美日韩综合久久久久久 | 国产欧美日韩精品亚洲av| 男插女下体视频免费在线播放| 丰满人妻熟妇乱又伦精品不卡| 免费无遮挡裸体视频| 搡老岳熟女国产| 国产黄色小视频在线观看| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| 久久久久久大精品| 精品熟女少妇八av免费久了| 男人狂女人下面高潮的视频| 天堂影院成人在线观看| 999久久久精品免费观看国产| 国产91精品成人一区二区三区| 午夜福利在线观看吧| 日本在线视频免费播放| 免费无遮挡裸体视频| 一区福利在线观看| 97热精品久久久久久| 直男gayav资源| 一个人看的www免费观看视频| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| 69人妻影院| 老司机福利观看| 亚洲最大成人中文| 国产精品av视频在线免费观看| 一区二区三区四区激情视频 | 亚洲,欧美精品.| 一边摸一边抽搐一进一小说| 久久人妻av系列| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 能在线免费观看的黄片| 色综合婷婷激情| 久久久久亚洲av毛片大全| 国产精品爽爽va在线观看网站| 国产日本99.免费观看| 欧美高清成人免费视频www| 神马国产精品三级电影在线观看| 搡老妇女老女人老熟妇| 国产v大片淫在线免费观看| 日本熟妇午夜| 亚洲激情在线av| 中亚洲国语对白在线视频| 亚洲精品乱码久久久v下载方式| 网址你懂的国产日韩在线| 草草在线视频免费看| 在线观看午夜福利视频| 99热这里只有是精品50| 亚洲精品亚洲一区二区| av在线天堂中文字幕| 国产精品亚洲av一区麻豆| 亚洲欧美日韩东京热| 中文字幕人成人乱码亚洲影| 午夜免费激情av| 亚洲在线自拍视频| 日韩免费av在线播放| 欧美日韩中文字幕国产精品一区二区三区| 又爽又黄无遮挡网站| 国产精品乱码一区二三区的特点| 国产欧美日韩一区二区三| 亚洲欧美清纯卡通| 成人精品一区二区免费| 精华霜和精华液先用哪个| 女人十人毛片免费观看3o分钟| 色av中文字幕| 三级国产精品欧美在线观看| 俄罗斯特黄特色一大片| 黄色丝袜av网址大全| 午夜两性在线视频| 国产黄a三级三级三级人| 少妇人妻精品综合一区二区 | 琪琪午夜伦伦电影理论片6080| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 亚洲精品色激情综合| 一级作爱视频免费观看| 色综合婷婷激情| 亚洲国产精品sss在线观看| 一区二区三区四区激情视频 | 久久久久久大精品| 亚洲美女搞黄在线观看 | 国产白丝娇喘喷水9色精品| 丰满的人妻完整版| 国产成人aa在线观看| 在线观看舔阴道视频| 五月伊人婷婷丁香| a在线观看视频网站| 脱女人内裤的视频| 亚洲乱码一区二区免费版| 亚洲七黄色美女视频| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 少妇的逼好多水| x7x7x7水蜜桃| 亚洲经典国产精华液单 | 免费av观看视频| 美女大奶头视频| 久久天躁狠狠躁夜夜2o2o| 99热这里只有是精品50| 亚洲精品在线观看二区| 国产精品久久视频播放| 久久久久久久久中文| 亚洲乱码一区二区免费版| 夜夜爽天天搞| 中文字幕高清在线视频| 国产高清激情床上av| 青草久久国产| 亚洲无线观看免费| 深夜精品福利| 亚洲精品456在线播放app | 色综合亚洲欧美另类图片| 国产亚洲av嫩草精品影院| 美女被艹到高潮喷水动态| 最好的美女福利视频网| 老鸭窝网址在线观看| 午夜精品一区二区三区免费看| 婷婷精品国产亚洲av在线| 欧美zozozo另类| 国产色爽女视频免费观看| 亚洲精品亚洲一区二区| 国产成人影院久久av| 精品久久久久久久人妻蜜臀av| 精品人妻视频免费看| 国产高清有码在线观看视频| 99热精品在线国产| 亚洲片人在线观看| 国产人妻一区二区三区在| 国产免费男女视频| 一二三四社区在线视频社区8| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 深夜精品福利| 校园春色视频在线观看| 1000部很黄的大片| 中亚洲国语对白在线视频| 国产欧美日韩精品一区二区| 国产主播在线观看一区二区| 中文字幕精品亚洲无线码一区| 伊人久久精品亚洲午夜| 床上黄色一级片| 国产成人欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 夜夜爽天天搞| 赤兔流量卡办理| 男人舔女人下体高潮全视频| 少妇人妻精品综合一区二区 | 51午夜福利影视在线观看| 尤物成人国产欧美一区二区三区| 最好的美女福利视频网| 最近在线观看免费完整版| 内射极品少妇av片p| 国产精品一及| 亚洲第一欧美日韩一区二区三区| 一区二区三区高清视频在线| 天天躁日日操中文字幕| 色综合婷婷激情| 草草在线视频免费看| 熟女人妻精品中文字幕| 久久国产精品影院| 我要搜黄色片| 丁香六月欧美| 午夜a级毛片| 久久久国产成人免费| 久久精品国产亚洲av涩爱 | 在线观看舔阴道视频| 亚洲综合色惰| 美女高潮的动态| 国产精品亚洲一级av第二区| 一级毛片久久久久久久久女| 脱女人内裤的视频| 久久久久久久久大av| 亚洲人成伊人成综合网2020| 精品一区二区三区人妻视频| 69av精品久久久久久| av专区在线播放| avwww免费| 精品熟女少妇八av免费久了| 全区人妻精品视频| 日韩av在线大香蕉| 久久久久久久久久成人| 99热只有精品国产| 欧美高清性xxxxhd video| ponron亚洲| 久久久久久久精品吃奶| 在线观看舔阴道视频| 国产精品,欧美在线| 色视频www国产| 成人特级黄色片久久久久久久| 18禁在线播放成人免费| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 亚洲久久久久久中文字幕| 亚洲专区中文字幕在线| 高清毛片免费观看视频网站| 欧美激情久久久久久爽电影| 有码 亚洲区| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 麻豆成人午夜福利视频| 久久精品91蜜桃| 精品免费久久久久久久清纯| 中文字幕高清在线视频| 免费高清视频大片| 在线播放国产精品三级| 精品久久久久久久久亚洲 | 小蜜桃在线观看免费完整版高清| 在线免费观看的www视频| 欧美3d第一页| 欧美日韩黄片免| 久久6这里有精品| 亚洲av中文字字幕乱码综合| 亚州av有码| 成人亚洲精品av一区二区| 91字幕亚洲| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 99国产综合亚洲精品| 国产成人a区在线观看| 国产一区二区在线观看日韩| 亚洲一区高清亚洲精品| 久久久久久九九精品二区国产| 亚洲精品影视一区二区三区av| 欧美日本视频| 久久久久久久午夜电影| 久久久精品大字幕| 亚洲,欧美精品.| 午夜福利成人在线免费观看| av在线观看视频网站免费| 亚洲成av人片免费观看| 一级黄色大片毛片| 我要搜黄色片| 久久久久免费精品人妻一区二区| 90打野战视频偷拍视频| 窝窝影院91人妻| 免费电影在线观看免费观看| 国产久久久一区二区三区| av在线观看视频网站免费| 99在线视频只有这里精品首页| 成年版毛片免费区| 99国产精品一区二区蜜桃av| 久久久久久久久久成人| 免费看美女性在线毛片视频| 国产亚洲精品综合一区在线观看| 国产亚洲欧美在线一区二区| 国内精品久久久久精免费| 国产主播在线观看一区二区| 精品99又大又爽又粗少妇毛片 | 欧美日韩乱码在线| 欧美日韩国产亚洲二区| 亚洲av电影在线进入| 精品一区二区三区av网在线观看| 一区福利在线观看| 中文字幕免费在线视频6| 国产精品影院久久| 国产麻豆成人av免费视频| 婷婷精品国产亚洲av| 欧美一区二区亚洲| 日日摸夜夜添夜夜添小说| 丰满人妻一区二区三区视频av| 极品教师在线视频| 亚洲精品日韩av片在线观看| 国产精品亚洲av一区麻豆| 嫩草影视91久久| 男人和女人高潮做爰伦理| 女同久久另类99精品国产91| 一区二区三区免费毛片| 久久国产精品影院| 一边摸一边抽搐一进一小说| 麻豆成人av在线观看| а√天堂www在线а√下载| 在线播放无遮挡| АⅤ资源中文在线天堂| 国产色婷婷99| 国产精品一区二区三区四区久久| 国产在视频线在精品| 少妇被粗大猛烈的视频| 成熟少妇高潮喷水视频| 五月伊人婷婷丁香| 色视频www国产| 嫩草影视91久久| 别揉我奶头~嗯~啊~动态视频| 久久精品91蜜桃| 18美女黄网站色大片免费观看| 国产精品乱码一区二三区的特点| 亚洲七黄色美女视频| 国产午夜精品论理片| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 久久久久久久久久黄片| av在线老鸭窝| 欧美黑人欧美精品刺激| 淫妇啪啪啪对白视频| 他把我摸到了高潮在线观看| 黄色丝袜av网址大全| 久久精品夜夜夜夜夜久久蜜豆| 97超级碰碰碰精品色视频在线观看| 中文字幕av成人在线电影| 久久久国产成人免费| 精品国产亚洲在线| 真人一进一出gif抽搐免费| 午夜两性在线视频| 久久九九热精品免费| 在线a可以看的网站| 人妻夜夜爽99麻豆av| 成人午夜高清在线视频| 国产免费男女视频| 观看美女的网站| 国产熟女xx| 俄罗斯特黄特色一大片| 熟妇人妻久久中文字幕3abv| 性色av乱码一区二区三区2| 欧美日韩国产亚洲二区| 精品免费久久久久久久清纯| 国产精品国产高清国产av| 成人午夜高清在线视频| 亚洲av二区三区四区| 男女视频在线观看网站免费| 亚洲经典国产精华液单 | 久久久久久久久久成人| 两个人视频免费观看高清| 免费看光身美女| 欧美日韩乱码在线| 国产精品一区二区三区四区免费观看 | 午夜福利18| 免费在线观看亚洲国产| 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 99国产精品一区二区三区| 天美传媒精品一区二区| 亚洲av五月六月丁香网| 男女视频在线观看网站免费| 99视频精品全部免费 在线| 性色avwww在线观看| 在线天堂最新版资源| 国产成人影院久久av| 亚洲中文字幕日韩| 国产淫片久久久久久久久 | 十八禁人妻一区二区| 九九在线视频观看精品| 国产亚洲精品久久久com| 在线天堂最新版资源| 亚洲五月婷婷丁香| av在线蜜桃| 99在线视频只有这里精品首页| 1000部很黄的大片| 国产中年淑女户外野战色| 嫩草影视91久久| 国产日本99.免费观看| 日本与韩国留学比较| 最新在线观看一区二区三区| 亚洲精品成人久久久久久| 精品日产1卡2卡| 18禁黄网站禁片免费观看直播| 一级av片app| 亚洲精品日韩av片在线观看| 亚洲不卡免费看| 欧美午夜高清在线| 国产色婷婷99| 一个人免费在线观看的高清视频| 亚洲,欧美,日韩| 欧美成人性av电影在线观看| 三级男女做爰猛烈吃奶摸视频| 日韩高清综合在线| 12—13女人毛片做爰片一| 能在线免费观看的黄片| 国产三级在线视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产自在天天线| 99热这里只有精品一区| 美女黄网站色视频| 免费一级毛片在线播放高清视频| 久久久久久国产a免费观看| 99久久99久久久精品蜜桃| 亚洲av一区综合| 欧美精品啪啪一区二区三区| www.www免费av| 欧美日韩黄片免| 国产精品美女特级片免费视频播放器| 日本与韩国留学比较| 91麻豆av在线| 高清日韩中文字幕在线| 欧美最黄视频在线播放免费| 舔av片在线| 日韩国内少妇激情av| 国产精品久久久久久人妻精品电影| 国产黄片美女视频| 毛片女人毛片| 91字幕亚洲| 美女高潮的动态| 在线观看一区二区三区| 观看免费一级毛片| 成年女人永久免费观看视频| 国产伦在线观看视频一区| 亚洲av第一区精品v没综合| 亚洲av五月六月丁香网| 欧美成狂野欧美在线观看| 国产爱豆传媒在线观看| 国产69精品久久久久777片| 国产一区二区在线观看日韩| 九色国产91popny在线| 麻豆国产97在线/欧美| 在线观看舔阴道视频| 日日干狠狠操夜夜爽| 黄色日韩在线| 最近在线观看免费完整版| 久久久久亚洲av毛片大全| 免费看日本二区| 久久久成人免费电影| 日韩大尺度精品在线看网址| 亚洲狠狠婷婷综合久久图片| 国产精品国产高清国产av| 十八禁人妻一区二区| 欧美一区二区精品小视频在线| 欧美色欧美亚洲另类二区| 国产一区二区在线观看日韩| av在线老鸭窝| 我要看日韩黄色一级片| 国产乱人伦免费视频| 免费av毛片视频| 中文在线观看免费www的网站| 欧美精品国产亚洲| 亚洲午夜理论影院| 天天躁日日操中文字幕| 99久久久亚洲精品蜜臀av| 亚洲欧美激情综合另类| 午夜福利欧美成人| 国产在线精品亚洲第一网站| 最近视频中文字幕2019在线8| 日韩国内少妇激情av| 久久久久国内视频| 成人毛片a级毛片在线播放| 欧美3d第一页| 宅男免费午夜| 91在线观看av| 99久久无色码亚洲精品果冻| 每晚都被弄得嗷嗷叫到高潮| 一级av片app| 国产精品一及| 亚洲午夜理论影院| 嫩草影院入口| 国产精品综合久久久久久久免费| 性欧美人与动物交配| 国产精品久久久久久精品电影| 免费人成在线观看视频色| 日本免费a在线| 久久国产精品影院| 三级毛片av免费| 色噜噜av男人的天堂激情| 亚洲av成人精品一区久久| 舔av片在线| 国产精品1区2区在线观看.| 国产欧美日韩一区二区精品| 日韩成人在线观看一区二区三区| 国产 一区 欧美 日韩| 色哟哟·www| 亚洲va日本ⅴa欧美va伊人久久| 一本久久中文字幕| 亚洲 国产 在线| 一卡2卡三卡四卡精品乱码亚洲| 日本精品一区二区三区蜜桃| 久久久久免费精品人妻一区二区| 日韩欧美精品免费久久 | 成熟少妇高潮喷水视频| 欧美精品国产亚洲| 国产精品乱码一区二三区的特点| 男人和女人高潮做爰伦理| 国产午夜精品久久久久久一区二区三区 | 国产一区二区在线av高清观看| 国产伦在线观看视频一区| 在线观看午夜福利视频| 亚洲精品在线观看二区| 尤物成人国产欧美一区二区三区| 综合色av麻豆| 又紧又爽又黄一区二区| 欧美xxxx性猛交bbbb| 少妇被粗大猛烈的视频| 免费人成在线观看视频色| 69人妻影院| 色哟哟哟哟哟哟| 久久草成人影院| 精品乱码久久久久久99久播| 亚洲一区二区三区不卡视频| 国产成人影院久久av| av天堂在线播放| 尤物成人国产欧美一区二区三区| 噜噜噜噜噜久久久久久91| 亚洲精品一卡2卡三卡4卡5卡| 精品熟女少妇八av免费久了| 国内揄拍国产精品人妻在线| 丁香六月欧美| 亚洲七黄色美女视频| 色哟哟哟哟哟哟| 如何舔出高潮| 麻豆国产97在线/欧美| 伦理电影大哥的女人| 亚洲成人久久爱视频| 热99在线观看视频| 看黄色毛片网站| 男女视频在线观看网站免费| 97碰自拍视频| 午夜久久久久精精品| 乱人视频在线观看| 久久精品91蜜桃| 国产成人欧美在线观看| 亚洲av成人av| 欧美精品啪啪一区二区三区| 精品熟女少妇八av免费久了| 亚洲美女视频黄频| 亚洲成av人片在线播放无| 欧美成人一区二区免费高清观看| 久久欧美精品欧美久久欧美| 中文在线观看免费www的网站| 欧美性感艳星| 久久精品91蜜桃| 久久精品国产99精品国产亚洲性色| 亚洲av免费在线观看| 伊人久久精品亚洲午夜| 成人av在线播放网站| 久久精品国产自在天天线| 国产激情偷乱视频一区二区| 18美女黄网站色大片免费观看| av中文乱码字幕在线| 国产精品一区二区免费欧美| 国产伦人伦偷精品视频| 一区二区三区高清视频在线| 日本三级黄在线观看| 一个人免费在线观看的高清视频| 中文亚洲av片在线观看爽| 首页视频小说图片口味搜索| 亚洲自偷自拍三级| 欧美精品啪啪一区二区三区| 国产精品美女特级片免费视频播放器| 久久久成人免费电影| 亚洲一区二区三区不卡视频| 欧美色视频一区免费| 男人舔女人下体高潮全视频| 日韩欧美在线乱码| 久久九九热精品免费| 亚洲精品日韩av片在线观看| 99久久成人亚洲精品观看| 男女床上黄色一级片免费看| 免费一级毛片在线播放高清视频| 美女xxoo啪啪120秒动态图 | 十八禁人妻一区二区| 大型黄色视频在线免费观看| 国产精品久久久久久久久免 | 十八禁人妻一区二区| 国产高潮美女av| 国产视频内射| 亚洲乱码一区二区免费版| 国产精品久久电影中文字幕| 国产精品,欧美在线| 国产精品野战在线观看| 欧美日韩乱码在线| 免费av不卡在线播放| 一级作爱视频免费观看| 午夜影院日韩av| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 亚洲成av人片在线播放无| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 狂野欧美白嫩少妇大欣赏| av国产免费在线观看| 最近最新中文字幕大全电影3| 色综合站精品国产| 欧美绝顶高潮抽搐喷水| 日韩欧美精品免费久久 | 可以在线观看的亚洲视频| 亚洲自拍偷在线| 午夜久久久久精精品| 综合色av麻豆| 亚洲色图av天堂| 麻豆国产av国片精品| 午夜福利18| 色av中文字幕| 精品久久久久久久久av| 欧美区成人在线视频|