• <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ù)矩陣的增減性及中值定理
    欧美激情极品国产一区二区三区| 国产伦在线观看视频一区| 国产亚洲欧美精品永久| 久久天躁狠狠躁夜夜2o2o| 夜夜爽天天搞| 日韩大码丰满熟妇| 国产高清有码在线观看视频 | 亚洲精品国产区一区二| 黄色成人免费大全| 亚洲av五月六月丁香网| 美女国产高潮福利片在线看| 18禁观看日本| 视频区欧美日本亚洲| 日韩大码丰满熟妇| 久久午夜综合久久蜜桃| 国产视频一区二区在线看| 18禁黄网站禁片免费观看直播| 免费高清在线观看日韩| e午夜精品久久久久久久| 国产精品野战在线观看| 伊人久久大香线蕉亚洲五| 啦啦啦 在线观看视频| 国产黄色小视频在线观看| 亚洲人成伊人成综合网2020| 久久久久免费精品人妻一区二区 | 中亚洲国语对白在线视频| 亚洲第一欧美日韩一区二区三区| 岛国在线观看网站| 一区二区日韩欧美中文字幕| 一边摸一边抽搐一进一小说| 男女那种视频在线观看| 少妇粗大呻吟视频| www.999成人在线观看| 俺也久久电影网| √禁漫天堂资源中文www| www日本在线高清视频| 人人妻人人澡欧美一区二区| 免费女性裸体啪啪无遮挡网站| 亚洲欧美激情综合另类| 国产精品99久久99久久久不卡| 色在线成人网| 欧美在线一区亚洲| 久久久久精品国产欧美久久久| 在线十欧美十亚洲十日本专区| 村上凉子中文字幕在线| av在线天堂中文字幕| 久久久精品国产亚洲av高清涩受| 一本大道久久a久久精品| 精品电影一区二区在线| 嫩草影院精品99| 国产成人一区二区三区免费视频网站| 久久 成人 亚洲| 亚洲 国产 在线| 国产三级黄色录像| av片东京热男人的天堂| 久久精品国产综合久久久| 男人舔女人的私密视频| 免费高清在线观看日韩| av片东京热男人的天堂| 成年版毛片免费区| 黄片小视频在线播放| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费视频内射| 国产精品久久久人人做人人爽| 欧美黑人巨大hd| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 18禁国产床啪视频网站| 麻豆国产av国片精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲 国产 在线| 91成年电影在线观看| 国产麻豆成人av免费视频| 色播在线永久视频| 午夜福利18| 性色av乱码一区二区三区2| 亚洲成国产人片在线观看| 亚洲自偷自拍图片 自拍| 亚洲美女黄片视频| 老鸭窝网址在线观看| 午夜成年电影在线免费观看| 亚洲人成网站在线播放欧美日韩| 日韩三级视频一区二区三区| 欧美日韩福利视频一区二区| 久久久久久久久中文| 亚洲真实伦在线观看| 女性生殖器流出的白浆| 18美女黄网站色大片免费观看| 国产一卡二卡三卡精品| 成年版毛片免费区| 婷婷精品国产亚洲av| 在线观看免费视频日本深夜| 亚洲电影在线观看av| 久久久久久久久久黄片| 香蕉国产在线看| 啦啦啦观看免费观看视频高清| 色哟哟哟哟哟哟| 天堂√8在线中文| 91国产中文字幕| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 国产成人啪精品午夜网站| xxx96com| 少妇 在线观看| 日本黄色视频三级网站网址| 超碰成人久久| 白带黄色成豆腐渣| 精品福利观看| 老熟妇仑乱视频hdxx| 亚洲欧美日韩无卡精品| 国产av一区二区精品久久| 免费在线观看完整版高清| videosex国产| 免费在线观看完整版高清| 哪里可以看免费的av片| 欧美激情高清一区二区三区| 国产成人精品久久二区二区91| 一本久久中文字幕| 亚洲一区高清亚洲精品| 久久精品国产综合久久久| 日本五十路高清| 国产高清激情床上av| 久久久久久人人人人人| 亚洲第一欧美日韩一区二区三区| 国产精品 欧美亚洲| 国产高清videossex| 亚洲专区字幕在线| 变态另类成人亚洲欧美熟女| 在线观看免费视频日本深夜| 色在线成人网| 黄片小视频在线播放| 99国产精品一区二区蜜桃av| 麻豆成人午夜福利视频| 亚洲男人的天堂狠狠| 中文字幕高清在线视频| 一边摸一边做爽爽视频免费| 人人妻人人澡欧美一区二区| 国产欧美日韩精品亚洲av| 一级片免费观看大全| www.熟女人妻精品国产| 18禁观看日本| 久久天躁狠狠躁夜夜2o2o| 夜夜躁狠狠躁天天躁| 亚洲全国av大片| 99国产综合亚洲精品| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 色播亚洲综合网| 91大片在线观看| 久久婷婷成人综合色麻豆| 午夜福利免费观看在线| 国产激情久久老熟女| 亚洲熟女毛片儿| 亚洲精品粉嫩美女一区| 国产亚洲欧美精品永久| 国内少妇人妻偷人精品xxx网站 | 丁香欧美五月| 成年人黄色毛片网站| 欧美一级a爱片免费观看看 | 免费人成视频x8x8入口观看| 人人妻人人澡人人看| 可以在线观看毛片的网站| 日本五十路高清| 老鸭窝网址在线观看| 高潮久久久久久久久久久不卡| 嫩草影视91久久| 欧美zozozo另类| 最新美女视频免费是黄的| 成人国产一区最新在线观看| 又大又爽又粗| 中文字幕高清在线视频| 国产欧美日韩一区二区三| 欧美日韩一级在线毛片| 美国免费a级毛片| 日韩欧美一区视频在线观看| 日韩有码中文字幕| 国产乱人伦免费视频| 久久久水蜜桃国产精品网| 两人在一起打扑克的视频| 婷婷精品国产亚洲av| 国产av不卡久久| 国产精品亚洲av一区麻豆| 亚洲精品一卡2卡三卡4卡5卡| www.999成人在线观看| 亚洲一区高清亚洲精品| 国产高清视频在线播放一区| 亚洲精华国产精华精| 欧美最黄视频在线播放免费| av欧美777| 真人一进一出gif抽搐免费| 亚洲av成人一区二区三| 亚洲久久久国产精品| 亚洲成a人片在线一区二区| 欧美日韩精品网址| 国产精品爽爽va在线观看网站 | 2021天堂中文幕一二区在线观 | ponron亚洲| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 亚洲激情在线av| 三级毛片av免费| 亚洲欧美精品综合一区二区三区| av在线天堂中文字幕| 中文在线观看免费www的网站 | 极品教师在线免费播放| 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| tocl精华| 午夜激情av网站| 久久久久久大精品| 精品第一国产精品| 日本 av在线| 亚洲成av人片免费观看| 一二三四在线观看免费中文在| 国产亚洲欧美在线一区二区| 制服人妻中文乱码| 一a级毛片在线观看| 欧美日韩一级在线毛片| 成人特级黄色片久久久久久久| 99国产综合亚洲精品| 久久热在线av| 亚洲片人在线观看| 免费在线观看成人毛片| 成人av一区二区三区在线看| 18美女黄网站色大片免费观看| 亚洲午夜理论影院| 精品国产美女av久久久久小说| 午夜老司机福利片| 黄色成人免费大全| 国产高清有码在线观看视频 | 成人国产综合亚洲| av有码第一页| 两个人看的免费小视频| 国产蜜桃级精品一区二区三区| 女警被强在线播放| 亚洲国产高清在线一区二区三 | 国产单亲对白刺激| 国产精品二区激情视频| 一个人观看的视频www高清免费观看 | 成年版毛片免费区| 亚洲欧美日韩无卡精品| av在线天堂中文字幕| 久久香蕉激情| 看片在线看免费视频| 久久香蕉精品热| 亚洲国产精品999在线| 亚洲av第一区精品v没综合| 成年免费大片在线观看| 午夜久久久久精精品| 窝窝影院91人妻| 少妇粗大呻吟视频| 日韩三级视频一区二区三区| 欧美日韩福利视频一区二区| 91九色精品人成在线观看| 午夜福利18| 免费搜索国产男女视频| 精品久久蜜臀av无| av在线天堂中文字幕| 人成视频在线观看免费观看| 在线播放国产精品三级| 午夜激情av网站| 精品乱码久久久久久99久播| 欧美日韩福利视频一区二区| 成人永久免费在线观看视频| 国产精品影院久久| 免费看美女性在线毛片视频| 日本一区二区免费在线视频| 夜夜躁狠狠躁天天躁| 黄色成人免费大全| 日韩欧美三级三区| 国产色视频综合| 久久伊人香网站| 欧美日本视频| 免费看a级黄色片| av在线天堂中文字幕| 99在线视频只有这里精品首页| 午夜久久久久精精品| 在线十欧美十亚洲十日本专区| 亚洲精品国产区一区二| 日本精品一区二区三区蜜桃| 久久人妻福利社区极品人妻图片| 欧美丝袜亚洲另类 | 日日夜夜操网爽| 久久99热这里只有精品18| 少妇的丰满在线观看| 在线观看免费视频日本深夜| 1024视频免费在线观看| 777久久人妻少妇嫩草av网站| 天天添夜夜摸| 黑人操中国人逼视频| 99在线人妻在线中文字幕| 日韩欧美国产在线观看| av欧美777| 性欧美人与动物交配| 成人国产一区最新在线观看| 日本一区二区免费在线视频| 精品电影一区二区在线| 久久九九热精品免费| 国产精品乱码一区二三区的特点| 制服丝袜大香蕉在线| 亚洲av片天天在线观看| 琪琪午夜伦伦电影理论片6080| 男女那种视频在线观看| 久久久久久久久久黄片| 久久伊人香网站| 美女午夜性视频免费| 无人区码免费观看不卡| 久久久久久九九精品二区国产 | 亚洲在线自拍视频| 久久精品国产清高在天天线| 亚洲第一电影网av| 国产精品久久久久久精品电影 | 老汉色av国产亚洲站长工具| 日韩大尺度精品在线看网址| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 99精品欧美一区二区三区四区| 日韩欧美在线二视频| www日本黄色视频网| 日韩一卡2卡3卡4卡2021年| ponron亚洲| 日韩大码丰满熟妇| 国产精品 欧美亚洲| av福利片在线| 亚洲国产精品久久男人天堂| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区字幕在线| 18禁黄网站禁片免费观看直播| 免费在线观看亚洲国产| 最好的美女福利视频网| 午夜老司机福利片| 日本五十路高清| 黄色片一级片一级黄色片| 成人av一区二区三区在线看| 可以在线观看的亚洲视频| 亚洲精品在线美女| 成年人黄色毛片网站| 亚洲三区欧美一区| 男人的好看免费观看在线视频 | 亚洲全国av大片| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 午夜激情福利司机影院| 日韩欧美国产一区二区入口| av片东京热男人的天堂| 国产成人系列免费观看| 少妇的丰满在线观看| 亚洲自拍偷在线| 国内精品久久久久久久电影| 久久精品人妻少妇| 免费在线观看亚洲国产| 波多野结衣av一区二区av| 国产日本99.免费观看| 久久久久亚洲av毛片大全| 久久久精品国产亚洲av高清涩受| 免费看a级黄色片| 国产亚洲欧美98| 午夜免费成人在线视频| 桃色一区二区三区在线观看| 国产亚洲精品第一综合不卡| 久久精品国产99精品国产亚洲性色| 99国产极品粉嫩在线观看| 99热6这里只有精品| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 男人舔女人的私密视频| xxxwww97欧美| 无遮挡黄片免费观看| e午夜精品久久久久久久| 免费看a级黄色片| 成年免费大片在线观看| 美女午夜性视频免费| 欧美日韩亚洲国产一区二区在线观看| 黄色女人牲交| 精品久久久久久久久久免费视频| 亚洲av成人av| 一级毛片女人18水好多| 国产精品九九99| 久久国产亚洲av麻豆专区| 久久久久久国产a免费观看| 亚洲精品av麻豆狂野| 色综合站精品国产| 黑人操中国人逼视频| 久久中文字幕一级| 激情在线观看视频在线高清| 最近最新免费中文字幕在线| 国产亚洲精品久久久久久毛片| 在线观看免费日韩欧美大片| 一区二区三区高清视频在线| 久久这里只有精品19| 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 日韩成人在线观看一区二区三区| 国产精品久久电影中文字幕| 国产精品,欧美在线| 成人18禁高潮啪啪吃奶动态图| 在线免费观看的www视频| 欧美一级a爱片免费观看看 | 国产区一区二久久| 亚洲国产看品久久| 国产激情久久老熟女| 国产成人欧美| 热re99久久国产66热| 久久精品人妻少妇| av在线天堂中文字幕| 91老司机精品| 中文字幕人成人乱码亚洲影| 久99久视频精品免费| 身体一侧抽搐| 18禁裸乳无遮挡免费网站照片 | 老司机靠b影院| 国产精品 国内视频| 99国产精品一区二区三区| 国产在线观看jvid| 禁无遮挡网站| 中文字幕av电影在线播放| 此物有八面人人有两片| 亚洲aⅴ乱码一区二区在线播放 | 一级毛片精品| 一进一出抽搐gif免费好疼| 久久久久久久午夜电影| e午夜精品久久久久久久| 亚洲精品在线美女| 成在线人永久免费视频| x7x7x7水蜜桃| 国产99白浆流出| 日本a在线网址| 久久久久久久久中文| 久久久久久久久久黄片| 色av中文字幕| 黄色女人牲交| 啦啦啦韩国在线观看视频| 国产亚洲欧美在线一区二区| 亚洲自偷自拍图片 自拍| 男女视频在线观看网站免费 | 精品乱码久久久久久99久播| 亚洲av中文字字幕乱码综合 | 日韩中文字幕欧美一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久亚洲av毛片大全| 久久草成人影院| av在线播放免费不卡| 在线播放国产精品三级| 国产黄a三级三级三级人| 在线国产一区二区在线| 99国产精品一区二区蜜桃av| 51午夜福利影视在线观看| 岛国在线观看网站| 青草久久国产| 欧美+亚洲+日韩+国产| 国产单亲对白刺激| 国产私拍福利视频在线观看| 成人欧美大片| 久久中文字幕一级| 久久国产精品影院| 久久香蕉激情| 成人三级黄色视频| 久久久久精品国产欧美久久久| 亚洲精品美女久久久久99蜜臀| 中文字幕av电影在线播放| 欧美日韩亚洲综合一区二区三区_| 在线十欧美十亚洲十日本专区| 亚洲精品久久国产高清桃花| 一区二区三区高清视频在线| 免费看日本二区| 男女做爰动态图高潮gif福利片| 亚洲在线自拍视频| 久久国产精品影院| 一级作爱视频免费观看| 国产成人系列免费观看| 最新美女视频免费是黄的| 琪琪午夜伦伦电影理论片6080| 美国免费a级毛片| 日本 欧美在线| 国产一级毛片七仙女欲春2 | 久久狼人影院| 精品久久久久久久久久免费视频| av在线天堂中文字幕| 777久久人妻少妇嫩草av网站| tocl精华| 韩国精品一区二区三区| 国产激情久久老熟女| 免费高清视频大片| 美国免费a级毛片| 中国美女看黄片| 精品电影一区二区在线| 午夜福利成人在线免费观看| www.999成人在线观看| 国产1区2区3区精品| 亚洲av第一区精品v没综合| 夜夜看夜夜爽夜夜摸| 亚洲人成伊人成综合网2020| 午夜福利一区二区在线看| 一级毛片精品| 日韩av在线大香蕉| 欧美亚洲日本最大视频资源| 亚洲免费av在线视频| 久热这里只有精品99| 窝窝影院91人妻| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看 | 中文资源天堂在线| 亚洲精品久久成人aⅴ小说| 女人被狂操c到高潮| 婷婷精品国产亚洲av在线| 国产麻豆成人av免费视频| 亚洲全国av大片| 三级毛片av免费| 欧美又色又爽又黄视频| 99热6这里只有精品| 97碰自拍视频| 免费在线观看成人毛片| av中文乱码字幕在线| 日韩一卡2卡3卡4卡2021年| 国产av在哪里看| 日韩一卡2卡3卡4卡2021年| 国产男靠女视频免费网站| 很黄的视频免费| 91成人精品电影| 精品久久蜜臀av无| 精品国产一区二区三区四区第35| 免费看日本二区| 午夜a级毛片| 黄色丝袜av网址大全| 18禁国产床啪视频网站| 久久久久久久久免费视频了| 久久久久久国产a免费观看| 欧美不卡视频在线免费观看 | 国产1区2区3区精品| 免费搜索国产男女视频| 久久精品人妻少妇| 国产av一区在线观看免费| 亚洲精品国产精品久久久不卡| 亚洲国产精品合色在线| 亚洲色图 男人天堂 中文字幕| 丁香欧美五月| 成人一区二区视频在线观看| av视频在线观看入口| 成人av一区二区三区在线看| 精品久久久久久成人av| 欧美性猛交黑人性爽| 侵犯人妻中文字幕一二三四区| 黄色片一级片一级黄色片| 国产蜜桃级精品一区二区三区| 人人妻人人看人人澡| av超薄肉色丝袜交足视频| 日本免费一区二区三区高清不卡| 国产精品 欧美亚洲| 精品久久久久久久毛片微露脸| 麻豆成人av在线观看| 中文字幕精品亚洲无线码一区 | 99久久综合精品五月天人人| 日韩三级视频一区二区三区| 欧美精品亚洲一区二区| 18禁黄网站禁片午夜丰满| 日韩欧美一区二区三区在线观看| 韩国av一区二区三区四区| 国产精华一区二区三区| 巨乳人妻的诱惑在线观看| 啦啦啦观看免费观看视频高清| 国产成人精品无人区| 成人av一区二区三区在线看| 国产久久久一区二区三区| 亚洲av日韩精品久久久久久密| 国产精品98久久久久久宅男小说| 日韩一卡2卡3卡4卡2021年| 欧美黄色片欧美黄色片| 757午夜福利合集在线观看| 国产精品美女特级片免费视频播放器 | 久久国产亚洲av麻豆专区| 久久精品91无色码中文字幕| 午夜影院日韩av| 精品国产乱子伦一区二区三区| 亚洲欧美一区二区三区黑人| 免费在线观看影片大全网站| √禁漫天堂资源中文www| 成在线人永久免费视频| 亚洲一区二区三区色噜噜| 精品午夜福利视频在线观看一区| 国产精品日韩av在线免费观看| 精品少妇一区二区三区视频日本电影| 最新在线观看一区二区三区| 亚洲成av片中文字幕在线观看| 亚洲一区中文字幕在线| 精品一区二区三区四区五区乱码| 美女高潮喷水抽搐中文字幕| 亚洲最大成人中文| 女警被强在线播放| 免费电影在线观看免费观看| 中文字幕精品免费在线观看视频| 欧美日韩福利视频一区二区| 亚洲激情在线av| 日本熟妇午夜| 日韩有码中文字幕| 宅男免费午夜| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成国产人片在线观看| 国产黄a三级三级三级人| 又紧又爽又黄一区二区| 成年女人毛片免费观看观看9| 亚洲中文av在线| 女同久久另类99精品国产91| 亚洲第一欧美日韩一区二区三区| 国产精品日韩av在线免费观看| 中文亚洲av片在线观看爽| 老司机福利观看| 国产精品自产拍在线观看55亚洲| 999久久久精品免费观看国产| 法律面前人人平等表现在哪些方面| 丝袜美腿诱惑在线| 亚洲成人久久性| 国内久久婷婷六月综合欲色啪|