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

    山區(qū)土壤重金屬空間插值性能的改進(jìn)

    2014-03-27 02:26:04張榆霞李寶磊萬(wàn)國(guó)盛
    關(guān)鍵詞:分位中值精確度

    張榆霞, 李寶磊, 施 擇,萬(wàn)國(guó)盛

    1.云南省環(huán)境監(jiān)測(cè)中心站, 云南 昆明 650034 2.云南大學(xué)信息學(xué)院, 云南 昆明 650091

    土壤重金屬污染給世界各地帶來(lái)嚴(yán)重的環(huán)境和健康問(wèn)題[1]。土壤中重金屬的空間分布圖被廣泛應(yīng)用于識(shí)別污染源、預(yù)測(cè)污染趨勢(shì)以及控制潛在的污染風(fēng)險(xiǎn)。山區(qū)水域、懸崖和山谷等自然條件給土壤均勻采樣帶來(lái)了巨大的困難。根據(jù)客觀實(shí)際條件,利用有限采樣點(diǎn)的樣品測(cè)試數(shù)據(jù),獲得精準(zhǔn)的土壤重金屬空間分布圖,對(duì)于進(jìn)一步進(jìn)行空間評(píng)價(jià)、分析和預(yù)防環(huán)境污染非常有意義。為此,探索采樣稀疏、采樣點(diǎn)分布不規(guī)則的情況下,如何獲得精確穩(wěn)定的插值結(jié)果顯得至關(guān)重要。

    目前廣泛應(yīng)用于土壤性質(zhì)空間插值的方法主要有反距離(IDW)[2]、普通克里格(OK)[3]和徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(RBFANNs)[4]插值法?;趲缀螌W(xué)意義的IDW插值方法忽略了土壤重金屬含量高度復(fù)雜的非線性特征而導(dǎo)致精度不高。 基于統(tǒng)計(jì)學(xué)意義的OK插值法雖然得到了廣泛的應(yīng)用,但其應(yīng)用的假設(shè)條件和采樣要求[5],使其不能很好地描述具有非線性系統(tǒng)特征的土壤性質(zhì)空間分布。近年來(lái),許多學(xué)者將人工神經(jīng)網(wǎng)絡(luò)較強(qiáng)的非線性映射能力應(yīng)用于土壤性質(zhì)的空間分布研究中[6],并取得了較理想的結(jié)果,但由于人工神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過(guò)程具有隨機(jī)性,會(huì)降低其插值結(jié)果的穩(wěn)定性。為此,研究引入了一種集成徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(IRBFANNs)模型,用以提高山區(qū)土壤重金屬濃度插值的精確性和穩(wěn)定性。基于云南省楚雄市南部山區(qū)表層土壤中重金屬錳和釩的樣品測(cè)試數(shù)據(jù),進(jìn)行了3種不同等級(jí)采樣密度下的IDW、OK、RBFANNs和IRBFANNs插值法比較實(shí)驗(yàn)研究。

    1 實(shí)驗(yàn)部分

    1.1 IRBFANNs空間插值模型

    集成技術(shù)的基本思想是通過(guò)訓(xùn)練多個(gè)神經(jīng)網(wǎng)絡(luò),并將其預(yù)測(cè)結(jié)果進(jìn)行平均以期消除誤差,提供更精確、穩(wěn)定的預(yù)測(cè)[4]。IRBFANNs預(yù)測(cè)需要2個(gè)步驟:首先在使用該模型之前用Bagging[7]方法對(duì)總訓(xùn)練數(shù)據(jù)重復(fù)取樣獲得不同的子訓(xùn)練數(shù)據(jù)集,用以訓(xùn)練各個(gè)RBFANNs模型,并通過(guò)式(1)計(jì)算各個(gè)模型的權(quán)重;然后把預(yù)測(cè)點(diǎn)信息Xin(被預(yù)測(cè)點(diǎn)的經(jīng)度、緯度、臨近5個(gè)采樣點(diǎn)采樣值組成的輸入向量)輸入到各個(gè)RBFANNs模型,通過(guò)式(2)計(jì)算IRBFANNs模型的輸出。

    (1)

    (2)

    1.2 實(shí)驗(yàn)數(shù)據(jù)

    以云南省楚雄市南部以及周邊地區(qū)為試驗(yàn)區(qū)。該區(qū)域系云貴高原中部,紅河水系與金沙江水系分水嶺地帶,地跨100°59′E~101°52′E,24°1′N(xiāo)~25°3′N(xiāo),面積達(dá)9 938.641 5 km2,山地是該區(qū)域主要的地貌。地勢(shì)西北高、東南低,海拔為556~3 657 m,海拔落差達(dá)3 101 m。土壤主要為紫色土、水稻土、紅壤和黃棕壤。該區(qū)域河溝縱橫、山地海拔落差大的特點(diǎn)給采樣點(diǎn)布設(shè)帶來(lái)較大難度。以該區(qū)域內(nèi)42個(gè)采樣點(diǎn)的土壤重金屬錳和釩檢測(cè)值為研究數(shù)據(jù)。在實(shí)驗(yàn)之前,所有數(shù)據(jù)都做了歸一化處理,其分布滿(mǎn)足標(biāo)準(zhǔn)正態(tài)分布。

    1.3 實(shí)驗(yàn)設(shè)計(jì)

    為了比較IDW、OK、RBFANNs和IRBFANNs插值法插值的精確性和穩(wěn)定性及其受采樣密度的影響情況,基于42個(gè)采樣點(diǎn)數(shù)據(jù)進(jìn)行了A、B、C 3項(xiàng)實(shí)驗(yàn),分別隨機(jī)選取41、26和16個(gè)測(cè)試樣本子集,以保證研究結(jié)果的廣泛適用性。每項(xiàng)實(shí)驗(yàn)中都進(jìn)行了100次獨(dú)立的隨機(jī)測(cè)試,每次測(cè)試中都記錄預(yù)測(cè)誤差的均方值(RMSE)用于統(tǒng)計(jì)分析。IDW插值法通過(guò)預(yù)測(cè)點(diǎn)周?chē)?0個(gè)采樣點(diǎn)的歐氏距離以及采樣值進(jìn)行預(yù)測(cè);OK插值法通過(guò)使用Matlab工具箱中的dace函數(shù)實(shí)現(xiàn),回歸模型和相關(guān)函數(shù)分別為Regpoly2和Corrgauss, 相關(guān)函數(shù)的初始參數(shù)theta是10;RBFANNs通過(guò)使用Matlab工具箱中的newrb函數(shù)實(shí)現(xiàn),其中采用“試錯(cuò)法”確定的錳和釩的最優(yōu)散布常數(shù)分別為0.3和0.05,其他參數(shù)使用工具箱提供的默認(rèn)參數(shù);在IRBFANNs中,每個(gè)子訓(xùn)練集都是通過(guò)使用Bagging算法從總訓(xùn)練數(shù)據(jù)中隨機(jī)重復(fù)抽取80%數(shù)據(jù)生成。

    插值性能的評(píng)價(jià)指標(biāo)RMSE如式(3)所示:

    (3)

    式中,n為檢驗(yàn)數(shù)據(jù)集中的檢驗(yàn)數(shù)據(jù)的個(gè)數(shù),z*(xi)為插值方法對(duì)采樣點(diǎn)xi處土壤重金屬含量的預(yù)測(cè)值,z(xi)為采樣點(diǎn)xi處采樣獲得的土壤重金屬含量的真實(shí)值。RMSE值越小則預(yù)測(cè)誤差就越小,插值精度就越高;RMSE值的波動(dòng)范圍越小,則插值算法的穩(wěn)定性就越高。

    2 結(jié)果與討論

    2.1 精確度與穩(wěn)定性統(tǒng)計(jì)分析

    圖1展示了100個(gè)測(cè)試實(shí)驗(yàn)中得到的RMSE統(tǒng)計(jì)量。

    圖1 RMSE箱圖圖例

    由圖1可見(jiàn),該箱圖描述了統(tǒng)計(jì)數(shù)據(jù)的中值、上五分位值、下五分位值、最大值以及最小值。中值越小則算法的插值精度就越高,上、下五分位值的差越小,則誤差波動(dòng)范圍就越集中,算法的穩(wěn)定性就越強(qiáng)。

    4種方法在3種采樣密度下對(duì)土壤中錳和釩元素總量的預(yù)測(cè)實(shí)驗(yàn)結(jié)果如圖2和圖3所示。

    圖2 錳RMSE箱圖

    圖3 釩RMSE箱圖

    2.1.1 采樣密度對(duì)插值性能的影響

    由圖2和圖3可見(jiàn):①4種方法中,RMSE指標(biāo)的中值都隨著訓(xùn)練數(shù)據(jù)的減少而增加,說(shuō)明4種方法的精確度隨著采樣密度的降低而降低,并且RMSE指標(biāo)的上下五分之一分位值之間的差值隨著訓(xùn)練數(shù)據(jù)的減少而減少,這是由于訓(xùn)練數(shù)據(jù)減少時(shí)模型的精確度下降,預(yù)測(cè)誤差都偏大,導(dǎo)致RMSE值的范圍減小。②各類(lèi)方法性能下降的程度不一致。相比較而言,IDW和OK方法的RMSE中值隨著訓(xùn)練數(shù)據(jù)的減少具有較明顯的增加,這說(shuō)明IDW和OK方法的精確度易受采樣密度的影響。RBFANNs方法的RMSE中值的增加幅度沒(méi)有IDW和OK方法明顯,但是其上下五分之一分位值之間的差值遠(yuǎn)遠(yuǎn)大于其他方法,說(shuō)明受采樣密度降低的影響,RBFANNs方法的精確度雖然下降得不是那么明顯,但是穩(wěn)定性有明顯的下降。③在所有方法中,隨著采樣密度的降低,IRBFANNs方法的RMSE中值上升幅度最小,上下五分之一分位值之間的差值變化不大,說(shuō)明IRBFANNs精確度受采樣密度的影響較小且能保證插值的穩(wěn)定性。

    總之,雖然隨著訓(xùn)練數(shù)據(jù)的減少,所有方法的預(yù)測(cè)性能都有所下降,但I(xiàn)RBFANNs的RMSE中值增加最少,特別在訓(xùn)練數(shù)據(jù)數(shù)目為16時(shí),IRBFANNs的RMSE中值以及上下五分之一分位值之間的差值都是最小的,這表明IRBFANNs方法應(yīng)用在采樣點(diǎn)密度稀疏的情況時(shí),精確度和穩(wěn)定性?xún)?yōu)于其他方法。

    2.1.2 不同方法插值性能比較

    觀察圖2和圖3,同一個(gè)項(xiàng)目,不同方法之間的實(shí)驗(yàn)結(jié)果表明:

    1)訓(xùn)練數(shù)據(jù)數(shù)目最多(項(xiàng)目A,41點(diǎn))時(shí),同一種元素各個(gè)方法的RMSE中值和上下五分之一分位值之間的差值差異不大,表明訓(xùn)練數(shù)據(jù)數(shù)目最多時(shí),各個(gè)方法的差異不是很明顯。而OK方法獲得最小的RMSE中值,其他方法獲得的RMSE中值略大于OK方法,這說(shuō)明采樣密集時(shí),基于地理統(tǒng)計(jì)學(xué)的OK方法表現(xiàn)出眾。由圖3可見(jiàn),在釩的實(shí)驗(yàn)中,IRBFANNs方法的RMSE中值以及上下五分之一分位值之間的差值明顯小于RBFANNs方法,這說(shuō)明集成有助于提高基于神經(jīng)網(wǎng)絡(luò)方法的插值精度和穩(wěn)定性。

    2)當(dāng)訓(xùn)練數(shù)據(jù)數(shù)目為26點(diǎn)時(shí),各個(gè)方法獲得的中值以及上下五分之一分位值差值差異相對(duì)于項(xiàng)目A趨于明顯。傳統(tǒng)的IDW和OK方法,RMSE中值以及上下五分之一分位值差值最小。而基于RBFANNs方法獲得的中值以及上下五分之一分位值之間的差值都較大,這說(shuō)明傳統(tǒng)神經(jīng)網(wǎng)絡(luò)方法的精確度、穩(wěn)定性不如其他傳統(tǒng)方法。

    3)當(dāng)訓(xùn)練樣本減少到16點(diǎn)時(shí),傳統(tǒng)的IDW和OK方法獲得的中值以及上下五分之一分位值之間的差值均比RBFANNs方法小,但比IRBFANNs方法大。這說(shuō)明采樣密度稀疏時(shí),集成神經(jīng)網(wǎng)絡(luò)方法的精確度和穩(wěn)定性?xún)?yōu)于其他方法。IRBFANNs方法在RMSE中值以及上下五分之一分位值之間的差值都明顯小于RBFANNs方法,這一結(jié)論與項(xiàng)目A、B一致。這進(jìn)一步說(shuō)明集成有助于提高基于神經(jīng)網(wǎng)絡(luò)方法的插值精確度和穩(wěn)定性。所有方法中,IREBANNs方法獲得最小的RMSE中值以及上下五分之一分位值之間的差值,這表明IRBFANNs方法應(yīng)用在采樣點(diǎn)密度稀疏的情況時(shí)具有較好的插值精確度和穩(wěn)定性。

    總之,樣本點(diǎn)最多時(shí),各種方法的插值精確度和穩(wěn)定性差異不大,傳統(tǒng)OK和IDW方法略?xún)?yōu)于其他方法。樣本點(diǎn)適中時(shí),各個(gè)方法的插值精確度和穩(wěn)定性差異趨于明顯,各個(gè)方法對(duì)于不同的元素表現(xiàn)不一致。樣本點(diǎn)最少時(shí),各個(gè)方法的插值精確度和穩(wěn)定性差異明顯,IRBFANNs優(yōu)于其他方法,其插值精確度和穩(wěn)定性都最好。

    綜上所述,訓(xùn)練數(shù)據(jù)減少時(shí),所有模型的預(yù)測(cè)精確度和穩(wěn)定性都有不同程度下降,并且各種方法對(duì)于不同的元素下降程度不同;對(duì)于所有元素樣本點(diǎn)最少時(shí),盡管各種方法的性能都有明顯的下降,但是IRBFANNs方法插值精確度和穩(wěn)定性都優(yōu)于其他方法。

    2.2 空間插值可視化分析

    為了更直觀地觀察各個(gè)方法的插值效果,圖4顯示了錳元素基于4種方法,在不同采樣密度條件下的插值結(jié)果空間分布。

    2.2.1 采樣密度對(duì)插值效果的影響

    由圖4可見(jiàn),同一個(gè)方法,隨著采樣點(diǎn)的減少,各個(gè)方法獲得的空間分布圖的分辨率降低,其所描述的空間分布情況變得模糊。但是各個(gè)方法性能降低的程度有所不同。相比較而言,隨著采樣密度的降低,傳統(tǒng)的IDW和OK方法獲得的插值空間分布圖的分辨率、細(xì)節(jié)描述性有明顯的下降。基于神經(jīng)網(wǎng)絡(luò)的RBFANNs和IRBFANNs方法,獲得的空間插值分布圖所描述的空間分布情況和質(zhì)量,受采樣點(diǎn)減少的影響較小,細(xì)節(jié)保留較為完整。

    2.2.2 不同方法插值效果比較

    由圖4可見(jiàn),樣本點(diǎn)最多時(shí),IDW方法雖然提供了一個(gè)確定的錳元素含量插值表面,但是插值空間分布圖粗糙,所描述的空間分布情況缺乏空間連續(xù)性,這可能是由于樣本點(diǎn)較集中區(qū)域中,離插值點(diǎn)極近的樣本點(diǎn)對(duì)該插值點(diǎn)的估計(jì)值影響特別大,而孤立樣本點(diǎn)對(duì)各個(gè)方向插值點(diǎn)的作用是穩(wěn)定的逐漸衰弱的。

    OK方法獲得的插值空間分布圖多斑點(diǎn),顏色差異性較小,其描述的空間分布情況平滑,缺乏細(xì)節(jié)信息,可見(jiàn)OK方法受樣本點(diǎn)非均勻分布的影響容易產(chǎn)生孤島效應(yīng),其插值的結(jié)果缺乏合理性。

    與OK方法相比,RBFANNs方法獲得的插值空間分布圖描述了更多的空間差異性,并且從圖中還可以看出該元素的空間分布趨勢(shì),這在OK方法獲得的插值空間分布圖中是很難觀察出來(lái)的。RBFANNs方法較IDW和OK方法,在錳元素空間差異細(xì)節(jié)性描述以及空間分布趨勢(shì)方面又有一定程度的提高,但是該方法獲得的插值數(shù)據(jù)的范圍與樣本數(shù)據(jù)范圍(標(biāo)準(zhǔn)正態(tài)分布)不一致,插值數(shù)據(jù)最大值和最小值遠(yuǎn)遠(yuǎn)高于或者低于樣本數(shù)據(jù)的最大值和最小值,這是由于RBFANNs的插值表面極易受邊緣效應(yīng)的影響,研究區(qū)域的邊緣插值結(jié)果往往偏大或者偏小,造成插值結(jié)果不穩(wěn)定。

    IRBFANNs方法獲得的插值空間分布圖不僅平滑并具有較好的連續(xù)性,極少有斑點(diǎn),空間變異明顯,且插值數(shù)據(jù)范圍與樣本數(shù)據(jù)范圍相一致,可見(jiàn)其能夠合理詳細(xì)地描述元素空間分布的局部細(xì)節(jié)以及分布趨勢(shì)。同樣的采樣密度下,樣本點(diǎn)適中和樣本點(diǎn)最少時(shí),通過(guò)觀察所有元素的插值空間分布圖,可以得出與樣本點(diǎn)最多時(shí)相一致的結(jié)論。

    綜上所述,IRBFANNs插值方法能夠有效合理地描述土壤重金屬空間分布的空間變異性的細(xì)節(jié)和空間分布趨勢(shì)。與其他插值方法相比,IRBFANNs在樣本點(diǎn)較少的情況下能夠獲得質(zhì)量相對(duì)較好的土壤重金屬空間分布圖。這與“2.1節(jié)”中的統(tǒng)計(jì)分析結(jié)論相一致。

    3 結(jié)論

    使用IRBFANNs插值方法,可以提高土壤重金屬含量空間插值的性能。通過(guò)誤差統(tǒng)計(jì)和插值可視化分析,與傳統(tǒng)插值方法相比較,該方法在采樣密度稀疏的情況下,產(chǎn)生預(yù)測(cè)誤差的均值、中值以及上下五分位值都最小,表明在樣本點(diǎn)數(shù)量減少時(shí),IRBFANNs算法能夠獲得最好的插值精確度和穩(wěn)定性,據(jù)此獲得了更準(zhǔn)確的區(qū)域重金屬分布及趨勢(shì)圖,從而提高山區(qū)土壤重金屬分布預(yù)測(cè)性能。

    [1] Alloway B J, Ayres D C. Chemical principles of environmental pollution[M]. Florida: CRC Press, 1997.

    [2] Tomczak M. Spatial interpolation and its uncertainty using automated anisotropic inverse distance weighting (IDW) cross validation /jackknife approach[J]. Journal of Geographic Information and Decision Analysis,1998,2(2):18-30.

    [3] Oliver M A, Webster R. Kriging: a method of interpolation for geographical information systems[J]. International Journal of Geographical Information System,1990,4(3):313-332.

    [4] 劉思聰.B/S 結(jié)構(gòu)的云南省土壤重金屬空間插值分析系統(tǒng)[D].昆明:云南大學(xué),2012.

    [5] 王政權(quán).地統(tǒng)計(jì)學(xué)及在生態(tài)學(xué)中的應(yīng)用[M]. 北京:科學(xué)出版社,1999.

    [6] 何勇,張淑娟,方慧.基于人工神經(jīng)網(wǎng)絡(luò)的田間信息插值方法研究[J].農(nóng)業(yè)工程學(xué)報(bào),2004,20(3):120-123.

    [7] Breiman L. Bagging predictors[J]. Machine learning,1996,24(2):123-140.

    猜你喜歡
    分位中值精確度
    基于兩階段Expectile回歸的風(fēng)險(xiǎn)保費(fèi)定價(jià)
    當(dāng)前市場(chǎng)位置和潛在空間
    研究核心素養(yǎng)呈現(xiàn)特征提高復(fù)習(xí)教學(xué)精確度
    “硬核”定位系統(tǒng)入駐兗礦集團(tuán),精確度以厘米計(jì)算
    Lagrange中值定理的巧妙應(yīng)用
    微分中值定理教法研討
    基于“業(yè)績(jī)與薪酬雙對(duì)標(biāo)”的國(guó)有企業(yè)負(fù)責(zé)人薪酬研究
    西部論叢(2017年10期)2017-02-23 06:31:36
    后中值波電流脈沖MIG焊工藝
    一個(gè)有效估計(jì):半?yún)?shù)非時(shí)齊擴(kuò)散模型的局部線性復(fù)合分位回歸估計(jì)
    函數(shù)矩陣的增減性及中值定理
    大香蕉97超碰在线| av在线天堂中文字幕| 91精品伊人久久大香线蕉| 女人被狂操c到高潮| 成人综合一区亚洲| 夫妻性生交免费视频一级片| 欧美日韩精品成人综合77777| 久久人人爽人人爽人人片va| 男人和女人高潮做爰伦理| 狠狠精品人妻久久久久久综合| 91狼人影院| 1000部很黄的大片| 国产 亚洲一区二区三区 | 美女xxoo啪啪120秒动态图| 熟女人妻精品中文字幕| 少妇的逼水好多| 边亲边吃奶的免费视频| 亚洲欧洲日产国产| 久99久视频精品免费| 亚洲av不卡在线观看| 亚洲国产成人一精品久久久| 中文字幕av成人在线电影| 美女被艹到高潮喷水动态| 一本一本综合久久| 精品久久久久久久久av| 久久99热这里只有精品18| 国产极品天堂在线| 80岁老熟妇乱子伦牲交| 一级毛片aaaaaa免费看小| 亚洲人成网站高清观看| 最近最新中文字幕大全电影3| 高清毛片免费看| 亚洲精品成人av观看孕妇| 建设人人有责人人尽责人人享有的 | 国产成人一区二区在线| 高清视频免费观看一区二区 | 丝袜美腿在线中文| 久久97久久精品| 久久精品久久精品一区二区三区| 日日干狠狠操夜夜爽| 丝袜美腿在线中文| av在线播放精品| 亚洲国产精品成人综合色| 欧美一区二区亚洲| 欧美性猛交╳xxx乱大交人| 国产一区二区亚洲精品在线观看| 有码 亚洲区| 爱豆传媒免费全集在线观看| 久久国产乱子免费精品| 亚洲怡红院男人天堂| 免费播放大片免费观看视频在线观看| 国产精品国产三级国产av玫瑰| 亚洲熟妇中文字幕五十中出| 97在线视频观看| 精品国产三级普通话版| 高清av免费在线| av一本久久久久| 美女xxoo啪啪120秒动态图| 久久精品久久精品一区二区三区| 91精品伊人久久大香线蕉| 毛片一级片免费看久久久久| 亚洲国产最新在线播放| 亚洲综合精品二区| 免费看美女性在线毛片视频| 日韩成人av中文字幕在线观看| 综合色av麻豆| 一级二级三级毛片免费看| 美女xxoo啪啪120秒动态图| 国产成人a∨麻豆精品| 午夜福利在线在线| 2022亚洲国产成人精品| 免费看a级黄色片| 国产亚洲5aaaaa淫片| 久久精品人妻少妇| 亚洲精品视频女| 亚洲精品久久久久久婷婷小说| 亚洲精品aⅴ在线观看| 亚洲激情五月婷婷啪啪| 午夜福利在线在线| 国产不卡一卡二| 日韩大片免费观看网站| 一级片'在线观看视频| 有码 亚洲区| 精品一区在线观看国产| 高清av免费在线| a级毛片免费高清观看在线播放| 国产亚洲午夜精品一区二区久久 | 国产精品人妻久久久影院| 一级片'在线观看视频| 亚洲内射少妇av| 激情五月婷婷亚洲| 最近视频中文字幕2019在线8| 最近视频中文字幕2019在线8| 在线观看美女被高潮喷水网站| 一区二区三区免费毛片| 国产淫片久久久久久久久| 免费看av在线观看网站| 女的被弄到高潮叫床怎么办| 老师上课跳d突然被开到最大视频| 久久久欧美国产精品| 国产成人freesex在线| 肉色欧美久久久久久久蜜桃 | 亚洲激情五月婷婷啪啪| 国产精品久久久久久精品电影| 又粗又硬又长又爽又黄的视频| 成人特级av手机在线观看| 最近中文字幕高清免费大全6| 一区二区三区免费毛片| 午夜爱爱视频在线播放| 欧美 日韩 精品 国产| 免费观看精品视频网站| 欧美人与善性xxx| 九九在线视频观看精品| 亚洲乱码一区二区免费版| 亚洲综合色惰| 成人漫画全彩无遮挡| 嫩草影院入口| 好男人在线观看高清免费视频| 国产精品一二三区在线看| 乱码一卡2卡4卡精品| 免费av观看视频| 久久久久久久久久久免费av| 午夜福利高清视频| 欧美日韩在线观看h| 国产精品国产三级国产av玫瑰| 97在线视频观看| 国产精品久久视频播放| 韩国av在线不卡| 亚洲美女搞黄在线观看| 国产淫语在线视频| 黄色配什么色好看| 亚洲欧美一区二区三区黑人 | 免费大片黄手机在线观看| 少妇熟女aⅴ在线视频| 亚洲电影在线观看av| 亚洲精品亚洲一区二区| 九九久久精品国产亚洲av麻豆| 男女边吃奶边做爰视频| 国产免费视频播放在线视频 | 精品国产露脸久久av麻豆 | 国产精品女同一区二区软件| 18禁裸乳无遮挡免费网站照片| 久99久视频精品免费| 成人亚洲精品av一区二区| 99久久精品一区二区三区| 亚洲经典国产精华液单| 亚洲精品日韩在线中文字幕| 成人国产麻豆网| 深爱激情五月婷婷| 六月丁香七月| 国产精品美女特级片免费视频播放器| 日本三级黄在线观看| 久久久精品免费免费高清| 日韩国内少妇激情av| 最近中文字幕2019免费版| 国产黄色免费在线视频| 亚洲美女视频黄频| 超碰av人人做人人爽久久| 久久久久九九精品影院| 成年免费大片在线观看| 又黄又爽又刺激的免费视频.| 日韩在线高清观看一区二区三区| 在线观看美女被高潮喷水网站| 国产精品久久久久久av不卡| 国产中年淑女户外野战色| 欧美xxxx黑人xx丫x性爽| 欧美xxxx性猛交bbbb| 99久久人妻综合| 亚洲av中文字字幕乱码综合| 国产黄片美女视频| 日产精品乱码卡一卡2卡三| 国产精品一及| 欧美一区二区亚洲| 日本欧美国产在线视频| 一级av片app| 一级爰片在线观看| 国产伦一二天堂av在线观看| 精品国产露脸久久av麻豆 | 国产高清国产精品国产三级 | 日本熟妇午夜| 纵有疾风起免费观看全集完整版 | 日韩成人av中文字幕在线观看| 久久久久网色| 国产精品伦人一区二区| 国产日韩欧美在线精品| 久久久久久久久久久免费av| 22中文网久久字幕| 精品午夜福利在线看| 亚洲国产精品成人综合色| 最近的中文字幕免费完整| 高清在线视频一区二区三区| 亚洲,欧美,日韩| 日韩精品青青久久久久久| 精品国产三级普通话版| 偷拍熟女少妇极品色| 老司机影院成人| 高清在线视频一区二区三区| www.av在线官网国产| 国产精品国产三级专区第一集| 最近最新中文字幕免费大全7| 成人午夜高清在线视频| 99热6这里只有精品| 国产成人精品一,二区| 在线a可以看的网站| 亚洲熟女精品中文字幕| 色吧在线观看| 丝袜美腿在线中文| 美女高潮的动态| 午夜免费观看性视频| 国产免费福利视频在线观看| 亚洲精品成人久久久久久| av一本久久久久| 日韩欧美国产在线观看| 尾随美女入室| ponron亚洲| 免费电影在线观看免费观看| 两个人视频免费观看高清| 亚洲欧美中文字幕日韩二区| 高清在线视频一区二区三区| 免费看光身美女| 免费观看在线日韩| 亚洲国产av新网站| 欧美日韩视频高清一区二区三区二| 国产精品美女特级片免费视频播放器| 国产在视频线在精品| 人妻一区二区av| av网站免费在线观看视频 | 色哟哟·www| 18禁在线无遮挡免费观看视频| 成人二区视频| 国产午夜精品一二区理论片| 大话2 男鬼变身卡| 日韩欧美三级三区| 性插视频无遮挡在线免费观看| 国产精品国产三级国产专区5o| 一级爰片在线观看| 精品国产一区二区三区久久久樱花 | 男女啪啪激烈高潮av片| 亚洲精品色激情综合| 午夜福利视频精品| 亚洲国产成人一精品久久久| 亚洲欧美成人精品一区二区| 亚洲电影在线观看av| av免费在线看不卡| 精品久久久久久久久av| 日本免费a在线| 久久国内精品自在自线图片| 简卡轻食公司| 欧美一级a爱片免费观看看| 天堂中文最新版在线下载 | 亚洲av国产av综合av卡| 亚洲内射少妇av| 国产老妇伦熟女老妇高清| 99热这里只有是精品50| 欧美区成人在线视频| 久久99热这里只有精品18| av国产久精品久网站免费入址| 免费看日本二区| 在线播放无遮挡| 深夜a级毛片| 欧美日韩视频高清一区二区三区二| 啦啦啦啦在线视频资源| 国产探花极品一区二区| 亚洲成人精品中文字幕电影| 超碰av人人做人人爽久久| 国产单亲对白刺激| 欧美激情国产日韩精品一区| 边亲边吃奶的免费视频| 久久久久久久久久黄片| 亚洲精品日本国产第一区| 国产探花在线观看一区二区| 国产精品一区二区性色av| 国产一区亚洲一区在线观看| 大香蕉97超碰在线| 男女那种视频在线观看| 久久久久国产网址| 草草在线视频免费看| 亚洲国产最新在线播放| 99久久人妻综合| 亚洲电影在线观看av| 777米奇影视久久| 国产人妻一区二区三区在| 联通29元200g的流量卡| 纵有疾风起免费观看全集完整版 | 国语对白做爰xxxⅹ性视频网站| 亚洲av成人精品一区久久| 亚洲av日韩在线播放| 国产成人福利小说| 三级毛片av免费| 亚洲性久久影院| 欧美日韩精品成人综合77777| 一二三四中文在线观看免费高清| 亚洲精品一区蜜桃| 国产老妇女一区| 九九久久精品国产亚洲av麻豆| 午夜福利在线观看吧| 亚洲人成网站高清观看| 国产爱豆传媒在线观看| 深爱激情五月婷婷| 久久人人爽人人片av| 丝瓜视频免费看黄片| 国产毛片a区久久久久| 国产精品一区二区性色av| 国产精品一区二区三区四区久久| 在线 av 中文字幕| 看黄色毛片网站| 男人舔奶头视频| 久久精品久久精品一区二区三区| 精品亚洲乱码少妇综合久久| 黄色日韩在线| 久久久久精品久久久久真实原创| 好男人视频免费观看在线| 午夜精品一区二区三区免费看| 纵有疾风起免费观看全集完整版 | 一本一本综合久久| 午夜视频国产福利| 深夜a级毛片| 97热精品久久久久久| 99热这里只有是精品50| 亚洲精品aⅴ在线观看| 激情 狠狠 欧美| 中文字幕亚洲精品专区| 国产男女超爽视频在线观看| 国产精品一区二区三区四区久久| 欧美3d第一页| 伦理电影大哥的女人| 成人国产麻豆网| 成人av在线播放网站| 国产一区二区三区综合在线观看 | 亚洲va在线va天堂va国产| a级一级毛片免费在线观看| 亚洲自拍偷在线| 嫩草影院入口| 一级毛片aaaaaa免费看小| 亚洲成人精品中文字幕电影| 成人无遮挡网站| 成年女人在线观看亚洲视频 | 欧美bdsm另类| 国语对白做爰xxxⅹ性视频网站| av专区在线播放| 久久精品国产自在天天线| 亚洲高清免费不卡视频| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 性插视频无遮挡在线免费观看| 亚洲三级黄色毛片| 亚洲精品乱码久久久久久按摩| 嫩草影院新地址| 日日啪夜夜爽| av网站免费在线观看视频 | 国产精品蜜桃在线观看| 全区人妻精品视频| 我要看日韩黄色一级片| av线在线观看网站| 1000部很黄的大片| 一夜夜www| 国产在视频线在精品| 神马国产精品三级电影在线观看| 一二三四中文在线观看免费高清| 天堂网av新在线| 国产一区二区在线观看日韩| 老司机影院成人| 日本熟妇午夜| 我要看日韩黄色一级片| 精品久久久久久成人av| 边亲边吃奶的免费视频| 久久99热这里只频精品6学生| 别揉我奶头 嗯啊视频| 91久久精品电影网| 免费观看性生交大片5| 国产一区有黄有色的免费视频 | 中国国产av一级| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 女人十人毛片免费观看3o分钟| 亚洲三级黄色毛片| 97人妻精品一区二区三区麻豆| 久久国内精品自在自线图片| av在线亚洲专区| 国产精品人妻久久久影院| 又爽又黄无遮挡网站| 亚洲欧美精品自产自拍| 久久综合国产亚洲精品| 精品亚洲乱码少妇综合久久| 中文天堂在线官网| 两个人视频免费观看高清| 欧美一级a爱片免费观看看| 91av网一区二区| 国产精品嫩草影院av在线观看| 国产精品不卡视频一区二区| 日日撸夜夜添| 国产精品一二三区在线看| 国产有黄有色有爽视频| 中国国产av一级| 国产高清三级在线| 亚洲自拍偷在线| 婷婷色综合www| 久久精品国产亚洲av涩爱| av又黄又爽大尺度在线免费看| 麻豆乱淫一区二区| 99久久九九国产精品国产免费| 麻豆国产97在线/欧美| 亚洲国产精品国产精品| 只有这里有精品99| 亚洲va在线va天堂va国产| 国产伦一二天堂av在线观看| 五月玫瑰六月丁香| 嘟嘟电影网在线观看| 日韩av不卡免费在线播放| 国产探花极品一区二区| 国产成人精品婷婷| 中文字幕免费在线视频6| 边亲边吃奶的免费视频| 大片免费播放器 马上看| kizo精华| 国产欧美另类精品又又久久亚洲欧美| 亚洲一区高清亚洲精品| 国产精品久久久久久久电影| 一夜夜www| 天堂av国产一区二区熟女人妻| 欧美日韩精品成人综合77777| 亚洲国产色片| 国产黄片美女视频| 欧美性感艳星| 久久精品熟女亚洲av麻豆精品 | 国精品久久久久久国模美| 国产欧美日韩精品一区二区| 精品一区二区免费观看| 国产精品一区www在线观看| 久久鲁丝午夜福利片| 亚洲精品国产成人久久av| 国产日韩欧美在线精品| 26uuu在线亚洲综合色| 一级黄片播放器| 三级国产精品欧美在线观看| 啦啦啦中文免费视频观看日本| av专区在线播放| 久久久成人免费电影| 免费在线观看成人毛片| 欧美最新免费一区二区三区| av在线天堂中文字幕| 日韩一区二区视频免费看| 99热网站在线观看| 如何舔出高潮| 久久99蜜桃精品久久| 能在线免费观看的黄片| 国产精品1区2区在线观看.| 在线观看美女被高潮喷水网站| 成年女人在线观看亚洲视频 | 亚洲精品乱久久久久久| 国产单亲对白刺激| 久久久亚洲精品成人影院| 天堂√8在线中文| 韩国av在线不卡| 国产高潮美女av| 禁无遮挡网站| 欧美不卡视频在线免费观看| 有码 亚洲区| 丝袜喷水一区| 国产一区二区在线观看日韩| 色视频www国产| 日韩av在线免费看完整版不卡| 99久国产av精品国产电影| 成人亚洲精品一区在线观看 | 秋霞伦理黄片| 久久久久久久久大av| 亚洲成人一二三区av| 日产精品乱码卡一卡2卡三| av女优亚洲男人天堂| 精品国产一区二区三区久久久樱花 | 亚洲精品影视一区二区三区av| 欧美人与善性xxx| 搡老乐熟女国产| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄| 一级黄片播放器| 久久久久久九九精品二区国产| 最后的刺客免费高清国语| www.av在线官网国产| 国产综合精华液| 少妇的逼水好多| 国产av不卡久久| 亚洲人成网站高清观看| 美女国产视频在线观看| 亚洲人成网站在线播| 国产精品国产三级专区第一集| 能在线免费观看的黄片| 国产精品一区二区性色av| 国模一区二区三区四区视频| 黄色配什么色好看| 国产黄片视频在线免费观看| 91狼人影院| 成年人午夜在线观看视频 | 久久这里只有精品中国| 免费看日本二区| 日韩欧美精品v在线| 国产精品一二三区在线看| av网站免费在线观看视频 | 午夜免费观看性视频| 精品久久久久久久久av| 亚洲三级黄色毛片| 亚洲av福利一区| 一级爰片在线观看| 日日撸夜夜添| 狠狠精品人妻久久久久久综合| 亚洲熟女精品中文字幕| 成年免费大片在线观看| 国产三级在线视频| 在线免费观看不下载黄p国产| 成人毛片60女人毛片免费| 日本一二三区视频观看| 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| 精品酒店卫生间| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 日本av手机在线免费观看| 亚洲国产日韩欧美精品在线观看| 一级毛片电影观看| 秋霞伦理黄片| 一级片'在线观看视频| 成人一区二区视频在线观看| 777米奇影视久久| 亚洲在线自拍视频| 亚洲精品色激情综合| 十八禁国产超污无遮挡网站| 色综合亚洲欧美另类图片| 国产在线一区二区三区精| 亚洲美女视频黄频| 日日啪夜夜爽| 日韩大片免费观看网站| 国产精品麻豆人妻色哟哟久久 | 真实男女啪啪啪动态图| 亚洲一级一片aⅴ在线观看| 精品亚洲乱码少妇综合久久| 精品久久久久久电影网| 尤物成人国产欧美一区二区三区| 国产成人免费观看mmmm| 国产在视频线在精品| 亚洲一级一片aⅴ在线观看| 丝袜喷水一区| 国产男人的电影天堂91| 大又大粗又爽又黄少妇毛片口| 国产av不卡久久| 青春草视频在线免费观看| 欧美精品国产亚洲| 精品国内亚洲2022精品成人| 少妇的逼好多水| 国产激情偷乱视频一区二区| 亚洲精品国产成人久久av| 日韩在线高清观看一区二区三区| 乱人视频在线观看| 欧美激情国产日韩精品一区| 亚洲av一区综合| 日本熟妇午夜| 毛片一级片免费看久久久久| 看十八女毛片水多多多| 成人无遮挡网站| 国产成人精品久久久久久| 麻豆成人午夜福利视频| 日韩人妻高清精品专区| 午夜激情福利司机影院| 亚洲乱码一区二区免费版| 欧美 日韩 精品 国产| 热99在线观看视频| 国产精品嫩草影院av在线观看| 人体艺术视频欧美日本| 午夜精品在线福利| 国产爱豆传媒在线观看| 一级黄片播放器| 亚洲欧美精品专区久久| 国产精品蜜桃在线观看| 99久久精品热视频| 真实男女啪啪啪动态图| 99热网站在线观看| 亚洲人成网站在线播| 久久99蜜桃精品久久| 成年人午夜在线观看视频 | 肉色欧美久久久久久久蜜桃 | 七月丁香在线播放| 免费看av在线观看网站| 中文字幕av成人在线电影| 免费观看的影片在线观看| 欧美激情久久久久久爽电影| or卡值多少钱| 国产欧美日韩精品一区二区| 国产黄a三级三级三级人| 免费少妇av软件| 草草在线视频免费看| 欧美 日韩 精品 国产| 亚洲伊人久久精品综合| 久久久a久久爽久久v久久| 久久人人爽人人爽人人片va| 看黄色毛片网站| 日本爱情动作片www.在线观看| 色综合色国产| 午夜老司机福利剧场| 永久免费av网站大全| 秋霞在线观看毛片| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 国产成人免费观看mmmm| 亚洲丝袜综合中文字幕| 亚洲av电影在线观看一区二区三区 | 久久精品久久久久久久性| 久久久久精品久久久久真实原创| 国产在视频线在精品| 亚洲av成人av| 中文资源天堂在线| 女人十人毛片免费观看3o分钟| 欧美日韩综合久久久久久| 男人爽女人下面视频在线观看| 久久久午夜欧美精品| 久久久久久久大尺度免费视频| 国产成人免费观看mmmm|