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

    基于序貫高斯條件模擬的土壤重金屬含量預測與不確定性評價
    ——以宜興市土壤Hg為例

    2018-08-27 03:29:30徐辰星濮勵杰徐彩瑤
    土壤學報 2018年4期
    關鍵詞:克里高斯不確定性

    徐辰星 濮勵杰,2? 朱 明,2 徐彩瑤 張 濛 許 艷

    (1 南京大學地理與海洋科學學院,南京 210023)

    (2 國土資源部海岸帶開發(fā)與保護重點實驗室,南京 210023)

    (3 蘇州科技大學環(huán)境科學與工程學院,江蘇蘇州 215009)

    土壤屬性制圖的主要目的是表達土壤異質(zhì)的空間信息,但由于空間預測模型的限制,由采樣點位置的觀測值來預測土壤性質(zhì)空間分布的過程中包含的不確定性,將通過一系列的傳遞并對最終結果產(chǎn)生深刻影響[1]。克里格估值只能提供簡化變異的空間分布模式,在極值和空間連續(xù)性模式敏感的應用中存在很大的局限性[2-4],比如典型的土壤污染數(shù)據(jù)。隨機模擬方法能克服“平滑效應”,它將數(shù)據(jù)作為一個整體,最大可能地接近真實的空間分布,因此近年來空間隨機模擬方法的應用逐漸成為地統(tǒng)計學研究的主要趨勢之一[5-8]。目前建模土壤性質(zhì)常用的隨機模擬方法包括序貫高斯條件模擬、序貫高斯協(xié)同模擬、序貫高斯指示模擬和序貫高斯指示協(xié)同模擬[8-10]。

    據(jù)中國2014年發(fā)布的《全國土壤污染狀況調(diào)查公報》顯示:全國不同區(qū)域的土壤樣品中總Hg含量為0.017~3.32 mg kg-1,表層0~20 cm的土壤Hg含量最高,點位超標率為1.6%[11]。長江三角洲重金屬污染嚴重,Hg是主要的重金屬污染元素。Hg進入土壤后,95%以上被固定在土壤表層[12],且受人為活動影響較大。近年來蘇南城市化建設與經(jīng)濟發(fā)展如雨后春筍,土壤污染的報道較為常見,本文以宜興市土壤重金屬地質(zhì)調(diào)查數(shù)據(jù)為例,利用序貫高斯條件模擬方法描述土壤Hg的空間分布特征,并利用高斯序貫指示模擬進行土壤Hg超過設定閾值的概率分布評價,探索該區(qū)域土壤重金屬含量空間預測的不確定性評價方法,對控制區(qū)域土壤重金屬風險具有重要的理論和現(xiàn)實意義。

    1 材料與方法

    1.1 研究區(qū)概況

    宜興市地處江蘇省南部(31°07′N~31°37′N,119°31′E~120°03′E),總面積1 996 km2(含太湖水域面積242.3 km2)。宜興市四季分明,全年溫暖濕潤,年平均氣溫15.7℃,1月平均氣溫2.8℃,7月平均氣溫28.3℃。年平均無霜期240多天。年平均雨日136.6 d,年平均降水量約1 177 mm,春夏雨水集中[13]。宜興市自改革開放以來工業(yè)發(fā)展迅速,地形因素(地勢南高北低,南部為丘陵山區(qū),北部為平原區(qū),東部為太湖瀆區(qū),西部為低洼圩區(qū))也使得土壤重金屬分布的異質(zhì)性較強[14]。研究區(qū)的土壤以黃棕壤、水稻土、潮土為主。

    1.2 樣品采集與分析

    數(shù)據(jù)來源于無錫市耕地質(zhì)量生態(tài)地球化學調(diào)查,野外采樣于2009年9月1日—15日完成。樣點布設采用網(wǎng)格法,網(wǎng)格寬度為1 km,實際野外采樣過程中,若樣點落在建筑物或交通道路上,則選擇采樣點周圍具備采集條件的地點采自表層0~20 cm,按照梅花采樣法,采集半徑20 m內(nèi)的4~8個點混合樣,裝入樣品袋帶回實驗室待測。各樣點均用GPS記錄其經(jīng)緯度,并詳細記錄周圍的環(huán)境信息。共采集了1 157個土壤表層樣品(圖1),土樣風干后用瑪瑙研缽研碎過20目尼龍篩,裝袋備用,土壤Hg含量依照GB/T 22105《土壤質(zhì)量總汞、總砷、總鉛的測定 原子熒光法》,采用電感耦合等離子體原子熒光光譜法(ICP-MS)測定。

    圖1 研究區(qū)土壤采樣點分布圖Fig. 1 Sampling sites distribution of study area

    1.3 數(shù)據(jù)處理

    經(jīng)實驗室化學分析后,有兩個樣點的值分別為1.34和0.825 mg kg-1,考慮到重金屬元素分布變異程度較高的固有屬性,參照拉依達準則將誤差標準擴充至平均值與10倍標準差的和,因此將這兩個數(shù)據(jù)視為實驗誤差剔除,實際用于本研究的樣本數(shù)量為1 155個。原始數(shù)據(jù)不符合正態(tài)分布,偏度和峰度分別為1.795和9.495。研究區(qū)內(nèi)有25個樣品高于國家土壤環(huán)境質(zhì)量二級標準(GB15618.1995)0.3 mg kg-1,占樣本總數(shù)的2.16%,高于全國平均水平[11]。本文利用GS+9.0進行半方差函數(shù)模型的擬合,運用ArcGIS10.1的高斯地統(tǒng)計模擬模塊實現(xiàn)序貫高斯條件模擬。為了研究地統(tǒng)計模擬的特點和可靠性,將原始數(shù)據(jù)分為模擬數(shù)據(jù)集和驗證數(shù)據(jù)集,以4 km×4 km網(wǎng)格對原采樣點進行重采樣,得到由309個樣本組成的模擬數(shù)據(jù)集,其余846個樣本則作為驗證數(shù)據(jù)集用于模擬精度的驗證。

    表1 土壤Hg含量描述性統(tǒng)計Table 1 Descriptive statistics for soil Hg content

    本研究使用Box-Cox法對模擬數(shù)據(jù)集進行轉(zhuǎn)換嘗試變化參數(shù),然后檢驗變換后的正態(tài)性是否足夠理想,最終設置參數(shù)為0.202,轉(zhuǎn)換后的偏度和峰度分別為0.0009和3.62,通過SPSS中的單樣本K-S檢驗,說明數(shù)據(jù)符合正態(tài)分布,可用于進一步分析。對模擬數(shù)據(jù)集中的樣點做半方差云圖和趨勢分析,未發(fā)現(xiàn)明顯的各向異性和趨勢。

    1.4 研究方法

    (1)半方差模型和序貫高斯條件模擬。本文采用半方差模型描述連續(xù)的土壤屬性的空間自相關結構,揭示可能影響數(shù)據(jù)分布的因素,并用序貫高斯條件模擬(Sequential Gaussian Conditional Simulation, SGCS)定量刻畫土壤屬性的變異性和不確定性,通過隨機變量的系列結果來統(tǒng)計不確定性。

    (2)精度評價。為了測試半方差模型的擬合效果以及驗證高斯地統(tǒng)計模擬的精確度和可靠性,需利用一種客觀的檢驗方法進行評價。對序貫高斯條件模擬的精度,采用驗證數(shù)據(jù)集中樣點的Hg含量的平均預測誤差(Mean Prediction Error, MPE)和均方根預測誤差(Root Mean Square Prediction Error, RMSPE)來評價[16]。

    (3)不確定性評價。本文采用高斯序貫指示模擬探討土壤Hg的不確定性[10,17]。單一空間位置x′上Hg含量的不確定性可以用該點Hg含量y′大于某個給定閾值 yt的概率來定量表示,空間不確定性(spatial uncertainty)可以用于評價基于概率Pexceed描繪重金屬空間分布的可信度[20]。

    2 結果與討論

    2.1 土壤Hg含量空間變異結構特征

    研究區(qū)土壤表層Hg含量的平均值為0.127 mg kg-1,變異系數(shù)達52.4%。半方差模型參數(shù)見表2,塊金系數(shù)達0.504,屬于中等強度變異,說明土壤Hg受到土壤母質(zhì)、地形等結構因素和耕作、土壤改良等隨機因素的共同影響。本文依據(jù)決定系數(shù)的大小,選擇了擬合效果(R2=0.930)最佳的球狀模型,并用此模型的參數(shù)進行克里格插值。

    克里格插值產(chǎn)生的交叉驗證報告,均方根預測誤差,標準平均值,平均標準誤差已達到最小,并且均方根預測誤差與平均標準誤差非常接近,符合一般情況下模型最優(yōu)的準則。本研究將序貫高斯條件模擬方法的實現(xiàn)次數(shù)設置為100,E-type模擬的土壤Hg空間分布與簡單克里格估值呈現(xiàn)了相似的趨勢(圖2),高值集中在研究區(qū)的東北部太湖瀆區(qū)和西部的洼地圩區(qū)的。

    表2 Hg含量空間變異結構特征Table 2 Structural characteristics of spatial variation of Hg content

    2.2 空間預測精度

    將序貫高斯條件模擬得到E-type型,與抽取的第1、25、50、75、100次實現(xiàn),提取柵格上的數(shù)值到驗證數(shù)據(jù)集中的846個點上,對比模擬值與真實測量值。隨機抽取的5個單次實現(xiàn)的平均預測誤差與均方根預測誤差均較簡單克里格插值大(相差在0.01以內(nèi)),E-type型與簡單克里格的平均預測誤差和均方根誤差相差在0.000 1以內(nèi),可見兩種不同原理的方法在研究區(qū)尺度上運用的精度可以視作無差別。將閾值設定為當?shù)赝寥繦g的背景值0.2 mg kg-1和用于不確定性研究的0.15 mg kg-1的分類錯誤率范圍分別為11.70%~12.88%和27.19%~30.14%。雖然5個單次實現(xiàn)的分類錯誤率均高于E-type,但是簡單克里格估值與E-type型的分類錯誤率不存在顯著差異,隨機模擬與傳統(tǒng)的地統(tǒng)計方法在精度上差異不顯著。

    簡單克里格得到的空間分布更為平滑,E-type型會隨著隨機模擬指定次數(shù)的增加而更趨向于克里格估值(圖2a、圖2b)。平滑效應造成估值的數(shù)值值域范圍較?。?1],克里格估值的取值范圍為0.050~0.265 mg kg-1,E-type型模擬取值范圍為0.048~0.276 mg kg-1。模擬數(shù)據(jù)集的峰度稍大,兩種方法均會造成柵格數(shù)據(jù)向中值靠攏。在這種情形下,克里格插值的平滑效應尤為明顯,數(shù)值集聚在0.121~0.18 mg kg-1。抽取的5個單次實現(xiàn)的模擬值均呈現(xiàn)正偏態(tài),對柵格上的數(shù)值統(tǒng)計:0.091~0.12 mg kg-1的頻數(shù)最大,0.061~0.09 mg kg-1和0.121~0.15 mg kg-1的頻數(shù)次之,與輸入數(shù)據(jù)的描述性統(tǒng)計一致。標準差在距離采樣點較遠的地方以及研究區(qū)邊界有較大的誤差方差,簡單克里格得到的標準差在采樣點半徑500 m內(nèi)較小,E-type型標準差(圖2c、圖2d)只在采樣點局部較小,主要是因為克里格方差僅依賴于采樣點的分布情況。

    序貫高斯條件模擬借助Monto-Carlo方法,依據(jù)原始數(shù)據(jù)遵循的概率密度函數(shù)求解更具優(yōu)勢,模擬值的方差在任何時候均等于原始數(shù)據(jù)所遵循的高斯分布的方差,避免了平滑效應。此外,序貫高斯模擬條件將局部變異性重新添加到了其生成的表面中,克服了克里格方差與實際樣本值無關的缺點[21]??死锔窆乐翟诿總€位置得到的是最優(yōu)估值,不存在任何可能的變化,誤差表征的不確定性表現(xiàn)在不同的預測點位置之間。序貫高斯條件模擬的隨機性服從條件,每次模擬結果再現(xiàn)輸入數(shù)據(jù)的統(tǒng)計特征,即均值、方差、協(xié)方差、變異函數(shù)均與原始數(shù)據(jù)的均值、方差、協(xié)方差、變異函數(shù)保持一致。若干次模擬實現(xiàn)就在每個位置上產(chǎn)生了多個模擬結果,具備了以概率論為基礎的不確定性評價,可以對該位置空間變量可能的取值區(qū)間做數(shù)理統(tǒng)計。

    2.3 空間預測的不確定性評價

    本文以評價標準0.15 mg kg-1作為土壤Hg含量的閾值為例進行序貫指示模擬,模擬次數(shù)設為1 000次,利用局部不確定性評價[17]計算各個點位超過該臨界閾值的概率,對這些區(qū)域做超閾值的聯(lián)合分布概率,可以評價“污染區(qū)”的風險程度。

    本文以1 000次高斯序貫指示模擬為例探討空間多個位置的不確定性,最大的標準差為0.053 8,落在了未采土壤樣本的水域,其次較大的標準差0.045落在臨界概率為55%~65%的區(qū)域內(nèi)。標準差與臨界概率兩者比值的最大值可用作表征相對最大誤差,最大比值為8.18%表明結果可靠。將指示柵格的大小設定為100 m2,用超閾值概率[17]表征不同置信度下的土壤Hg超閾值區(qū)域。按95%、90%、85%、75%的閾值概率分別有206個、307個、511個和1 838個柵格,再對這些區(qū)域進行超閾值的聯(lián)合評價,由表5所示,單點閾值概率達到95%、90%、85%區(qū)域的聯(lián)合概率分別為24.76%、16.61%和10.00%,多點聯(lián)合概率較單點臨界概率的評價標準更加嚴格。當單點的臨界概率達到75%時,采用聯(lián)合概率進行計算只有2.77%,說明盡管臨界概率比較高,但置信度不足以用作劃定依據(jù)。

    圖2 Hg空間預測圖和標準差圖Fig. 2 Spatial prediction map and standard deviation map of soil Hg content

    圖3 序貫高斯條件模擬的單次實現(xiàn)Fig. 3 Values distributions of SGCS

    表層(0~20 cm)土壤Hg含量空間表達波動性大的區(qū)域,主要集中在該市距離太湖約3 km的新莊鎮(zhèn)及和橋鎮(zhèn),同時該區(qū)域也是Hg含量較高的一帶。結合土地利用類型和地圖,新莊鎮(zhèn)有五十多家燈具企業(yè),江蘇省首座含汞燈管處理裝置在和橋鎮(zhèn)投入運行,用于粉碎和酸洗回收金屬汞和硫酸汞。固體廢棄物如燈管的露天堆放及管理不善,易造成包膜破碎和污染物質(zhì)向土壤滲透。污染物可通過大氣沉降、徑流等途徑[22],造成汞化合物對農(nóng)田、水源的嚴重污染。經(jīng)核實地圖變更信息,宜興市曾有六十余家涉鉛、鎘、汞、鉻等重金屬的化工企業(yè),上百家涉汞燈具企業(yè),數(shù)十家鉛蓄電池企業(yè)。以及紫砂壺、玻璃、陶瓷、搪瓷、橡膠使用到的化工無機染料,農(nóng)業(yè)有機肥的大量累積施用,均有可能造成土壤Hg含量的累積。土壤母質(zhì)是土壤汞的最基本來源,但外界干擾是引起土壤Hg含量和空間波動性升高的重要因素[12]。結合高程DEM,土地利用方式和坡度是影響汞流失的重要原因,南部丘陵地區(qū)多個山莊和風景區(qū)的土壤Hg含量低,且空間表達的波動性很小。

    圖4 臨界概率Pc條件下獲取的土壤Hg含量>0.15 mg kg-1的區(qū)域分布Fig. 4 Spatial distribution of Hg content>0.15 mg kg-1 under the conditional of critical probability Pc

    3 結 論

    兩種方法在本文的應用中精度無異;克里格估值具有明顯的“平滑效應”,體現(xiàn)在空間分布的連續(xù)、估值值域范圍的收束與數(shù)值分布向擬合面的集中,而序貫高斯條件模擬的數(shù)值服從概率密度函數(shù),尊重了原始數(shù)據(jù)的分布,較好地保留了數(shù)據(jù)的變異特征;高斯序貫指示模擬可以評價空間任意位置上Hg含量大于該閾值的概率,衍生的不確定性研究結果有助于定量評估環(huán)境風險或健康風險、檢測各種資源環(huán)境應用模型的穩(wěn)健性以及誤差傳播等問題。

    猜你喜歡
    克里高斯不確定性
    小高斯的大發(fā)現(xiàn)
    今晚不能去你家玩啦!
    知識窗(2023年12期)2024-01-03 01:38:55
    我可以咬一口嗎?
    知識窗(2023年2期)2023-03-05 11:28:27
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    你今天真好看
    天才數(shù)學家——高斯
    英鎊或繼續(xù)面臨不確定性風險
    中國外匯(2019年7期)2019-07-13 05:45:04
    你今天真好看
    讀者(2018年24期)2018-12-04 03:01:34
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    有限域上高斯正規(guī)基的一個注記
    黄片大片在线免费观看| av中文乱码字幕在线| 国产黄色免费在线视频| 99热国产这里只有精品6| 村上凉子中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 亚洲av日韩在线播放| 久久狼人影院| 黄色 视频免费看| 亚洲国产欧美网| 视频区图区小说| av网站免费在线观看视频| av视频免费观看在线观看| 亚洲精品国产色婷婷电影| 狠狠狠狠99中文字幕| 国产淫语在线视频| 女人被躁到高潮嗷嗷叫费观| a级毛片黄视频| 日日夜夜操网爽| 一区二区三区精品91| 自线自在国产av| 手机成人av网站| av网站在线播放免费| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 免费久久久久久久精品成人欧美视频| 90打野战视频偷拍视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人啪精品午夜网站| 亚洲九九香蕉| 每晚都被弄得嗷嗷叫到高潮| av一本久久久久| 国产欧美亚洲国产| 国产乱人伦免费视频| 夫妻午夜视频| 欧美日韩成人在线一区二区| 99国产极品粉嫩在线观看| www.熟女人妻精品国产| 久99久视频精品免费| 久久精品成人免费网站| 免费观看人在逋| 日本欧美视频一区| 色老头精品视频在线观看| 成人国语在线视频| 男女之事视频高清在线观看| 丰满的人妻完整版| 国产精品亚洲一级av第二区| 一级作爱视频免费观看| 亚洲精品中文字幕一二三四区| 日韩 欧美 亚洲 中文字幕| 少妇猛男粗大的猛烈进出视频| 人人妻人人澡人人看| 亚洲精品一二三| 国产亚洲一区二区精品| 国产成人啪精品午夜网站| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 欧美中文综合在线视频| 天堂中文最新版在线下载| 亚洲第一av免费看| 日韩制服丝袜自拍偷拍| 精品一区二区三区视频在线观看免费 | 国产免费av片在线观看野外av| 曰老女人黄片| 欧美国产精品一级二级三级| 三级毛片av免费| 亚洲性夜色夜夜综合| 欧美黑人精品巨大| 亚洲国产毛片av蜜桃av| 天堂俺去俺来也www色官网| 国产精品1区2区在线观看. | 在线观看免费午夜福利视频| 露出奶头的视频| videosex国产| 黑丝袜美女国产一区| 欧美日本中文国产一区发布| 极品少妇高潮喷水抽搐| 国产精品99久久99久久久不卡| 黄色毛片三级朝国网站| 丝袜美足系列| netflix在线观看网站| 成人18禁在线播放| 制服诱惑二区| 俄罗斯特黄特色一大片| www.自偷自拍.com| 亚洲欧美精品综合一区二区三区| 亚洲色图 男人天堂 中文字幕| 欧美日韩亚洲综合一区二区三区_| 国产免费现黄频在线看| 亚洲第一av免费看| 欧美精品啪啪一区二区三区| 精品卡一卡二卡四卡免费| 正在播放国产对白刺激| 久9热在线精品视频| 久久亚洲真实| 丝袜人妻中文字幕| 欧美在线一区亚洲| 中文字幕人妻丝袜一区二区| 亚洲国产欧美日韩在线播放| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区视频了| 极品少妇高潮喷水抽搐| 久久精品aⅴ一区二区三区四区| 99re6热这里在线精品视频| 国产高清激情床上av| 12—13女人毛片做爰片一| 99久久99久久久精品蜜桃| 亚洲欧洲精品一区二区精品久久久| 女人被狂操c到高潮| 国产色视频综合| 国产欧美亚洲国产| av有码第一页| 亚洲精品粉嫩美女一区| 亚洲国产毛片av蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 高潮久久久久久久久久久不卡| 搡老乐熟女国产| 看免费av毛片| 国产精品久久久久成人av| 欧美日韩国产mv在线观看视频| 曰老女人黄片| 日韩制服丝袜自拍偷拍| 国产91精品成人一区二区三区| а√天堂www在线а√下载 | 在线观看免费午夜福利视频| 亚洲av熟女| 美女高潮到喷水免费观看| 99在线人妻在线中文字幕 | 久久久久国产一级毛片高清牌| 多毛熟女@视频| 老司机午夜十八禁免费视频| 国产欧美日韩综合在线一区二区| 国产精华一区二区三区| 国产亚洲精品久久久久5区| 99热只有精品国产| 精品一品国产午夜福利视频| 人人妻人人爽人人添夜夜欢视频| 精品乱码久久久久久99久播| 人人妻人人澡人人爽人人夜夜| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产亚洲av香蕉五月 | 午夜福利一区二区在线看| 久久天堂一区二区三区四区| 一区二区日韩欧美中文字幕| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 国产蜜桃级精品一区二区三区 | 久久国产精品人妻蜜桃| 啦啦啦视频在线资源免费观看| 丝袜美足系列| 女人爽到高潮嗷嗷叫在线视频| 一区二区三区国产精品乱码| 正在播放国产对白刺激| 午夜老司机福利片| 亚洲国产精品sss在线观看 | 人妻丰满熟妇av一区二区三区 | 不卡一级毛片| 一边摸一边做爽爽视频免费| 免费在线观看视频国产中文字幕亚洲| 最近最新中文字幕大全电影3 | 国产精品国产av在线观看| 亚洲成人免费av在线播放| 天天躁夜夜躁狠狠躁躁| 亚洲精品国产区一区二| 精品卡一卡二卡四卡免费| 国产精品一区二区免费欧美| 一级毛片精品| 国产无遮挡羞羞视频在线观看| 在线视频色国产色| 啦啦啦 在线观看视频| 99精国产麻豆久久婷婷| 99国产精品一区二区蜜桃av | 国产精品偷伦视频观看了| 在线观看午夜福利视频| 午夜福利免费观看在线| 欧美日韩瑟瑟在线播放| 91老司机精品| 国产精品久久电影中文字幕 | 咕卡用的链子| www日本在线高清视频| 欧美日韩精品网址| 欧美成人午夜精品| 亚洲精品国产区一区二| 成年版毛片免费区| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美激情在线| 亚洲欧美一区二区三区久久| 欧美精品啪啪一区二区三区| 777久久人妻少妇嫩草av网站| 国产乱人伦免费视频| 九色亚洲精品在线播放| 日韩免费高清中文字幕av| 精品久久久精品久久久| 麻豆乱淫一区二区| 亚洲 国产 在线| 美女午夜性视频免费| 搡老岳熟女国产| 亚洲第一av免费看| 免费少妇av软件| 免费在线观看影片大全网站| 老司机在亚洲福利影院| 久热爱精品视频在线9| 日韩精品免费视频一区二区三区| 日韩欧美一区二区三区在线观看 | 啦啦啦免费观看视频1| 午夜免费成人在线视频| 一区二区三区激情视频| 一本综合久久免费| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 波多野结衣av一区二区av| 亚洲国产看品久久| 欧美日韩亚洲高清精品| 国产男靠女视频免费网站| 搡老岳熟女国产| 欧美亚洲日本最大视频资源| av免费在线观看网站| 热99国产精品久久久久久7| 欧美不卡视频在线免费观看 | 国产成人精品久久二区二区免费| 两个人看的免费小视频| 久久久久久久久久久久大奶| 99久久国产精品久久久| 高潮久久久久久久久久久不卡| 亚洲熟女毛片儿| 91麻豆精品激情在线观看国产 | 国产蜜桃级精品一区二区三区 | 国产精品一区二区免费欧美| 精品高清国产在线一区| 欧美精品亚洲一区二区| 亚洲少妇的诱惑av| 两个人免费观看高清视频| 亚洲人成电影观看| 看免费av毛片| 亚洲av成人不卡在线观看播放网| 久久中文字幕一级| 国产精品免费大片| 亚洲欧美精品综合一区二区三区| 99热只有精品国产| 免费在线观看黄色视频的| 色播在线永久视频| 免费少妇av软件| 一边摸一边抽搐一进一小说 | 一本一本久久a久久精品综合妖精| 后天国语完整版免费观看| 精品国产国语对白av| 一级a爱片免费观看的视频| 免费人成视频x8x8入口观看| 亚洲国产中文字幕在线视频| 久久久久久免费高清国产稀缺| 亚洲精品国产精品久久久不卡| 国产精品秋霞免费鲁丝片| 99久久99久久久精品蜜桃| 亚洲国产精品sss在线观看 | 国产1区2区3区精品| 欧美乱码精品一区二区三区| 777久久人妻少妇嫩草av网站| 国产精品香港三级国产av潘金莲| 亚洲av第一区精品v没综合| 91国产中文字幕| 亚洲精品乱久久久久久| 日日摸夜夜添夜夜添小说| 免费在线观看日本一区| 男男h啪啪无遮挡| 午夜福利影视在线免费观看| 亚洲专区国产一区二区| 国产成人欧美| 国产精品国产av在线观看| 欧美日韩亚洲综合一区二区三区_| 电影成人av| 视频区图区小说| 欧美一级毛片孕妇| 最新美女视频免费是黄的| 激情视频va一区二区三区| 久久精品国产亚洲av香蕉五月 | 久久精品成人免费网站| 一级作爱视频免费观看| av国产精品久久久久影院| 纯流量卡能插随身wifi吗| 大香蕉久久网| 一进一出抽搐动态| 国产人伦9x9x在线观看| 久久久久国产精品人妻aⅴ院 | 日本一区二区免费在线视频| 日韩欧美免费精品| 高清毛片免费观看视频网站 | 在线天堂中文资源库| 免费女性裸体啪啪无遮挡网站| 十八禁人妻一区二区| 女警被强在线播放| √禁漫天堂资源中文www| 欧美最黄视频在线播放免费 | 日韩熟女老妇一区二区性免费视频| 女人久久www免费人成看片| 国产成人精品无人区| aaaaa片日本免费| 欧美久久黑人一区二区| 老司机亚洲免费影院| 妹子高潮喷水视频| 久久狼人影院| 国产伦人伦偷精品视频| 午夜精品在线福利| av福利片在线| 精品欧美一区二区三区在线| 国产aⅴ精品一区二区三区波| 欧美 亚洲 国产 日韩一| 天天影视国产精品| 满18在线观看网站| 成人国语在线视频| 97人妻天天添夜夜摸| 久久久精品免费免费高清| 天天添夜夜摸| aaaaa片日本免费| 亚洲国产欧美网| 亚洲国产欧美日韩在线播放| 巨乳人妻的诱惑在线观看| 国产一卡二卡三卡精品| 19禁男女啪啪无遮挡网站| 国产精品偷伦视频观看了| 久久久久精品国产欧美久久久| 欧美不卡视频在线免费观看 | 亚洲av成人av| 69av精品久久久久久| 欧美成人免费av一区二区三区 | 亚洲国产精品合色在线| 亚洲片人在线观看| 国产一区二区激情短视频| 大香蕉久久网| bbb黄色大片| 男男h啪啪无遮挡| 欧美黑人精品巨大| 高清毛片免费观看视频网站 | 亚洲精品乱久久久久久| 又紧又爽又黄一区二区| 国产男靠女视频免费网站| 男女之事视频高清在线观看| 午夜两性在线视频| 老熟妇仑乱视频hdxx| 国产成人精品在线电影| 男人操女人黄网站| 在线播放国产精品三级| 精品福利观看| 日本欧美视频一区| 国产精品二区激情视频| 亚洲七黄色美女视频| 黄色片一级片一级黄色片| 久久人人爽av亚洲精品天堂| 国产精品久久久久久精品古装| av欧美777| 国内毛片毛片毛片毛片毛片| av福利片在线| 精品一区二区三区av网在线观看| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 两人在一起打扑克的视频| 黄色成人免费大全| 欧美黑人精品巨大| cao死你这个sao货| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 成年动漫av网址| 欧美在线一区亚洲| 成年女人毛片免费观看观看9 | 亚洲欧洲精品一区二区精品久久久| 欧美 亚洲 国产 日韩一| 91成人精品电影| 天堂中文最新版在线下载| e午夜精品久久久久久久| 美国免费a级毛片| 男女午夜视频在线观看| 天堂动漫精品| 男女之事视频高清在线观看| 久久久国产成人免费| 正在播放国产对白刺激| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 国产av又大| 狂野欧美激情性xxxx| 黄色毛片三级朝国网站| 91字幕亚洲| 亚洲九九香蕉| 欧美乱妇无乱码| 久久精品亚洲av国产电影网| av网站在线播放免费| 熟女少妇亚洲综合色aaa.| 国产精品电影一区二区三区 | 一二三四社区在线视频社区8| 69精品国产乱码久久久| 看片在线看免费视频| 免费看a级黄色片| 国内毛片毛片毛片毛片毛片| 狂野欧美激情性xxxx| 免费在线观看日本一区| 久久久久久免费高清国产稀缺| 国产91精品成人一区二区三区| 久久九九热精品免费| 亚洲情色 制服丝袜| 亚洲全国av大片| 国产欧美日韩精品亚洲av| 十八禁人妻一区二区| 国产人伦9x9x在线观看| 国产乱人伦免费视频| 在线观看免费视频网站a站| a级毛片在线看网站| 免费观看精品视频网站| 999久久久国产精品视频| 久久国产精品大桥未久av| 国产国语露脸激情在线看| 精品国产国语对白av| 欧美精品高潮呻吟av久久| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色 | 亚洲七黄色美女视频| 777米奇影视久久| 国产精品98久久久久久宅男小说| 久久精品国产亚洲av香蕉五月 | 久久国产亚洲av麻豆专区| 女人被躁到高潮嗷嗷叫费观| 亚洲精品粉嫩美女一区| 日本wwww免费看| 精品卡一卡二卡四卡免费| avwww免费| 久久精品91无色码中文字幕| 丝袜美足系列| 欧美av亚洲av综合av国产av| 色精品久久人妻99蜜桃| 成人精品一区二区免费| 久久天堂一区二区三区四区| 久久性视频一级片| 啦啦啦在线免费观看视频4| 国产日韩欧美亚洲二区| 欧美激情极品国产一区二区三区| 国产欧美日韩综合在线一区二区| a级毛片在线看网站| 欧美在线黄色| av福利片在线| 又紧又爽又黄一区二区| 久久人妻熟女aⅴ| 少妇猛男粗大的猛烈进出视频| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 久久精品亚洲精品国产色婷小说| 丰满的人妻完整版| 欧美色视频一区免费| 夜夜躁狠狠躁天天躁| 国产aⅴ精品一区二区三区波| 黑人巨大精品欧美一区二区mp4| 黄色视频,在线免费观看| 欧美另类亚洲清纯唯美| 美女视频免费永久观看网站| 免费女性裸体啪啪无遮挡网站| 国产精品免费大片| 97人妻天天添夜夜摸| a级片在线免费高清观看视频| 99香蕉大伊视频| 欧美日韩国产mv在线观看视频| 欧美精品高潮呻吟av久久| 久久天躁狠狠躁夜夜2o2o| 99国产精品99久久久久| 国产xxxxx性猛交| 日本一区二区免费在线视频| 在线观看www视频免费| 一边摸一边抽搐一进一小说 | 亚洲成人免费av在线播放| 桃红色精品国产亚洲av| 9191精品国产免费久久| 在线天堂中文资源库| 自拍欧美九色日韩亚洲蝌蚪91| 一级黄色大片毛片| 欧美日韩中文字幕国产精品一区二区三区 | 欧美黄色片欧美黄色片| 99热网站在线观看| 在线观看免费视频日本深夜| 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 黄色成人免费大全| 亚洲熟妇熟女久久| 老司机午夜福利在线观看视频| 国产视频一区二区在线看| 国产精品九九99| 国产高清videossex| 老熟女久久久| 午夜福利在线观看吧| 欧美乱码精品一区二区三区| 捣出白浆h1v1| av网站在线播放免费| 亚洲 国产 在线| 国产精品 国内视频| 黑人欧美特级aaaaaa片| 99久久人妻综合| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 国产淫语在线视频| 国产欧美日韩综合在线一区二区| 岛国在线观看网站| 精品国产一区二区三区久久久樱花| 久久热在线av| 两个人免费观看高清视频| av片东京热男人的天堂| 中文亚洲av片在线观看爽 | 成人18禁高潮啪啪吃奶动态图| 国产野战对白在线观看| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色 | 91老司机精品| 国产高清激情床上av| 日韩一卡2卡3卡4卡2021年| 国产精品一区二区精品视频观看| 久久中文字幕人妻熟女| 精品第一国产精品| 极品少妇高潮喷水抽搐| 这个男人来自地球电影免费观看| 亚洲av成人一区二区三| 在线免费观看的www视频| 欧美另类亚洲清纯唯美| av线在线观看网站| 18禁裸乳无遮挡免费网站照片 | 欧美黄色淫秽网站| 久久国产精品影院| 人人妻人人澡人人爽人人夜夜| 最新美女视频免费是黄的| 免费少妇av软件| 人人妻人人澡人人看| 一区二区三区激情视频| 亚洲一码二码三码区别大吗| 成人精品一区二区免费| 99热只有精品国产| 成人三级做爰电影| 韩国精品一区二区三区| 麻豆国产av国片精品| 在线观看免费视频日本深夜| 国产精品 国内视频| 丁香六月欧美| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看 | 免费观看精品视频网站| 99re在线观看精品视频| 美女 人体艺术 gogo| 国产一区二区三区视频了| 国产精品 国内视频| 美国免费a级毛片| 久久香蕉国产精品| 看免费av毛片| 国产精品九九99| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看| 精品久久久久久久久久免费视频 | 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看| 18禁裸乳无遮挡动漫免费视频| 18禁观看日本| 男人操女人黄网站| 少妇 在线观看| 狠狠狠狠99中文字幕| 色播在线永久视频| 亚洲精品粉嫩美女一区| 黑人欧美特级aaaaaa片| 久热这里只有精品99| 夫妻午夜视频| 丁香欧美五月| 中文字幕色久视频| 国产淫语在线视频| 欧美国产精品一级二级三级| 国产在线观看jvid| 成人国语在线视频| 伊人久久大香线蕉亚洲五| 制服诱惑二区| 女同久久另类99精品国产91| 精品免费久久久久久久清纯 | 欧美黑人精品巨大| 中文字幕色久视频| 国产野战对白在线观看| 久久久精品国产亚洲av高清涩受| 韩国精品一区二区三区| 国产成人免费观看mmmm| 亚洲少妇的诱惑av| 亚洲成国产人片在线观看| 丁香欧美五月| 很黄的视频免费| 18禁观看日本| 精品一区二区三区四区五区乱码| 久久人人97超碰香蕉20202| 老司机午夜十八禁免费视频| 天堂俺去俺来也www色官网| 午夜亚洲福利在线播放| 色播在线永久视频| 亚洲av日韩精品久久久久久密| 少妇猛男粗大的猛烈进出视频| 日韩欧美免费精品| 欧美激情久久久久久爽电影 | 久久久久久久国产电影| 亚洲人成伊人成综合网2020| 三级毛片av免费| 一进一出抽搐gif免费好疼 | 超色免费av| 丰满饥渴人妻一区二区三| a在线观看视频网站| 午夜两性在线视频| 国产精品成人在线| 久久久精品区二区三区| 欧美人与性动交α欧美精品济南到| 日韩免费av在线播放| 在线av久久热| 欧美在线黄色| 久久ye,这里只有精品| 老汉色av国产亚洲站长工具| 国产在线精品亚洲第一网站| videosex国产| 女性被躁到高潮视频| 国产黄色免费在线视频| 国产一区二区三区视频了| 女警被强在线播放| 亚洲少妇的诱惑av|