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

    農(nóng)用地土壤重金屬污染調(diào)查最優(yōu)網(wǎng)格尺度及布點優(yōu)化方法

    2021-09-19 10:35:46秦旭芝鄭涵文何瑞成丁文娟張軍崔冠男史沛麗謝云峰
    關(guān)鍵詞:布點插值尺度

    秦旭芝,鄭涵文,何瑞成,丁文娟,張軍,崔冠男,史沛麗,謝云峰*

    1.廣西壯族自治區(qū)生態(tài)環(huán)境監(jiān)測中心 2.桂林理工大學(xué)環(huán)境科學(xué)與工程學(xué)院 3.生態(tài)環(huán)境部土壤與農(nóng)業(yè)農(nóng)村生態(tài)環(huán)境監(jiān)管技術(shù)中心

    土壤重金屬污染來源多樣,由于重金屬具有難降解、難遷移等特點,一旦形成污染會在土壤中長期存在,并且在空間分布上呈現(xiàn)出明顯的異質(zhì)性和不連續(xù)性,加大準確判斷土壤中污染物空間分布的難度[1]?,F(xiàn)有技術(shù)條件下缺乏直接確定土壤污染物空間分布的方法,只能基于采樣調(diào)查數(shù)據(jù)通過數(shù)學(xué)模型進行空間插值[2-4],而這種研究方法的準確性受前期采樣調(diào)查結(jié)果的影響很大。大量研究表明,對于空間分布變異性很大的污染物,布點方案是影響污染調(diào)查結(jié)果準確性的重要因素,如Jenkins等[5]對土壤中三硝基甲苯空間分布變異性的研究發(fā)現(xiàn),95%的變異性來源于采樣位置;Argyraki等[6]發(fā)現(xiàn)由采樣位置造成的不確定性占整個污染物含量估計的不確定性的50%。結(jié)合調(diào)查區(qū)域的歷史背景、地理環(huán)境,科學(xué)有效地制定布點方案能大幅提高土壤污染調(diào)查的準確性[7-9]。目前常用的布點方法主要包括兩大類[10]:1)基于模型的布點方案,如網(wǎng)格布點法、隨機布點法等[11];2)基于設(shè)計的布點方案,研究者基于調(diào)查區(qū)域的環(huán)境背景及土地利用歷史等對可能存在污染的區(qū)域進行加密布點[12]。Brus等[13-15]對2種布點方案進行了詳細論述和比較。

    實際案例中多種布點方案常協(xié)同使用[16-20],但如何確定最優(yōu)采樣密度是土壤污染調(diào)查的重點和難點問題[21-24]。謝志宜等[25]在開展珠三角區(qū)域耕地土壤質(zhì)量研究時發(fā)現(xiàn),布點密度越大,調(diào)查區(qū)域土壤質(zhì)量刻畫越細致,但調(diào)查成本也越高,8 km×8 km的網(wǎng)格,能較好地保證調(diào)查結(jié)果的準確性,同時也能夠有效降低調(diào)查成本。工業(yè)場地土壤通常具有重金屬污染嚴重、污染重金屬單一且分布較為均勻等特點[26-29],重金屬污染調(diào)查尺度較小,調(diào)查時通常先采用系統(tǒng)布點法,排查整個調(diào)查區(qū)域內(nèi)可能存在的土壤環(huán)境風(fēng)險,然后結(jié)合分區(qū)布點法和專業(yè)判斷法對可能存在環(huán)境風(fēng)險的區(qū)域進行重點調(diào)查。蘭鵬鵬等[30]以40 m×40 m網(wǎng)格對湖南某冶煉場地進行分區(qū)布點,結(jié)果發(fā)現(xiàn)場地污染主要在1.0~2.5 m深層空間內(nèi)。

    為確保調(diào)查區(qū)域內(nèi)污染物空間分布的準確性,常需采用分階段的方式對研究區(qū)域進行調(diào)查,主要包括基礎(chǔ)信息調(diào)查、初步調(diào)查和加密調(diào)查等環(huán)節(jié)。場地污染物空間分布調(diào)查的誤差主要來源于污染物空間分布插值結(jié)果不確定性較大的區(qū)域[31-34],污染程度高估會導(dǎo)致調(diào)查費用和修復(fù)費用的增加,污染程度低估會導(dǎo)致環(huán)境風(fēng)險無法得到有效控制。點位優(yōu)化的目的在于在尋求降低調(diào)查結(jié)果不確定性的最佳樣本量,同時有效降低調(diào)查成本。本研究旨在探索農(nóng)用地土壤重金屬污染調(diào)查最優(yōu)布點密度及布點優(yōu)化,以提高土壤污染調(diào)查方法對污染區(qū)范圍和污染程度的估計精度,并為土壤污染調(diào)查提供方法學(xué)支持。

    1 材料與方法

    1.1 研究區(qū)域

    研究區(qū)域位于中南某老工業(yè)園區(qū)周邊的農(nóng)用地,地處111°58′E~113°05′E,27°20′55″N~28°05′40″N,調(diào)查區(qū)域面積約75萬m2,所在區(qū)域總體地貌輪廓北、西、南地勢高,中部、東部地勢低平,但地勢起伏較為和緩,高程為52.17~79.69 m,高差不大,地貌類型多樣,山地、丘陵、崗地、平原、水面俱備(圖1)。區(qū)域四季分明,降水充沛,盛夏高溫,冬季寒冷。年均降水量為 1 425 mm,4—7月降水較集中,期間多有洪水發(fā)生,年均氣溫為17.5 ℃。冬季多西北風(fēng),夏季多東南風(fēng)。研究區(qū)域周邊的工業(yè)園區(qū)具有50多年的化工生產(chǎn)歷史,工業(yè)企業(yè)數(shù)量多,“三廢”排放量大,污染物種類多,對周邊農(nóng)用地土壤環(huán)境影響較大。

    注: 紅色表示高程較高;黃色表示高程中等;綠色表示高程較低。圖1 研究區(qū)域內(nèi)地形地勢示意Fig.1 Map of terrain within the scope of study

    1.2 土壤采集與分析

    按照70 m×70 m網(wǎng)格設(shè)計要求,共布設(shè)168個表層土壤(0~20 cm)采樣點,現(xiàn)場利用GPS獲取采樣點經(jīng)緯度,使用木鏟進行表層土壤采樣,低溫保存于密封袋中帶回實驗室。土壤樣品經(jīng)風(fēng)干去除雜物后,研磨過100目尼龍篩。采用HNO3-H2O2消解法進行樣品預(yù)處理,采用石墨爐原子吸收光譜法測定土壤Cd濃度。每個樣品設(shè)3組平行樣,選10%的樣品做重復(fù)測試,并進行空白試驗,測定過程中添加國家標(biāo)準物質(zhì)(GSS-6)進行質(zhì)量控制。同時測定土壤理化性質(zhì)(表1)。另外還檢測了其他5種重金屬(Cu、Ni、Pb、Hg、Zn)濃度。檢測結(jié)果顯示,除Ni外其他重金屬均存在超標(biāo)情況,其中Cd超標(biāo)點位最多,其他重金屬僅個別點位超標(biāo)。

    表1 土壤理化性質(zhì)參數(shù)Table 1 Physical and chemical properties of soil

    1.3 子樣本的抽取

    將調(diào)查區(qū)域分為70 m×70 m、100 m×100 m、160 m×160 m和200 m×200 m 4種網(wǎng)格尺度。先隨機抽取18個點位作為獨立驗證集〔圖2(a)〕,構(gòu)建除70 m×70 m尺度外的其他3種尺度網(wǎng)格,將距離網(wǎng)格中心點最近的點位作為該尺度下的采樣點,最終得到100 m×100 m尺度網(wǎng)格點位數(shù)87個,160 m×160 m 尺度網(wǎng)格點位數(shù)37個,200 m×200 m尺度網(wǎng)格點位數(shù)24個(圖2)。

    圖2 抽取樣本子集及其驗證點位分布Fig.2 Sampling subset and distribution of verification points

    1.4 土壤Cd濃度空間插值及其精度評價

    利用ArcGIS 10.5軟件對不同調(diào)查尺度下土壤Cd濃度進行空間分析,利用Matheron矩估計法,以權(quán)重最小二次方擬合理論半方差函數(shù)模型。采用克里格(Kriging)法對4種采樣尺度測定結(jié)果進行插值,分析不同采樣尺度檢測結(jié)果的精確度。由于4種采樣尺度土壤中Cd濃度均不符合正態(tài)分布,使用BOX-COX變換進行數(shù)據(jù)轉(zhuǎn)換,經(jīng)K-S 檢驗后符合正態(tài)分布。評價結(jié)果采用半方差函數(shù)及交叉驗證確定最佳精確度。結(jié)合調(diào)查區(qū)域土壤pH,根據(jù)GB 15618—2018《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險管控標(biāo)準(試行)》,采用0.6 mg/kg作為區(qū)域內(nèi)Cd污染評價標(biāo)準。

    2 結(jié)果與評價

    2.1 不同采樣尺度下土壤Cd濃度變異性

    對4種采樣尺度下土壤Cd進行半變異分析,結(jié)果如圖3和表2所示。從圖3看出,4種采樣尺度的采樣點間距離均小于空間相關(guān)距離,說明4種采樣尺度均適用于空間變異分析,但不同尺度下的空間相關(guān)距離存在差異。從圖3還可看出,土壤Cd濃度均表現(xiàn)出了一定的線性趨勢,說明存在較好的半方差結(jié)構(gòu),且采樣尺度越小線性趨勢越明顯。塊金值通常被認為是代表隨機變異的量,而基臺值代表變量空間變異的結(jié)構(gòu)性方差,二者比值為塊金系數(shù)。當(dāng)塊金系數(shù)小于25%時,表示變量具有強烈的空間相關(guān)性;塊金系數(shù)為25%~50%,表示變量具有明顯的空間自相關(guān);塊金系數(shù)為50%~75%,表示變量具有中等的空間自相關(guān);塊金系數(shù)大于75%時,表示變量空間自相關(guān)性較小。從表2可知,70 m×70 m、100 m×100 m、160 m×160 m、200 m×200 m 4種采樣尺度下土壤Cd濃度的塊金系數(shù)分別為65.82%、46.40%、58.65%和45.37%,均在25%~75%之間,表明研究區(qū)土壤Cd濃度具有中等空間相關(guān)性,受結(jié)構(gòu)性因素和隨機因素影響。

    圖3 不同采樣尺度下土壤Cd濃度半方差分布Fig.3 Semi-variance distribution of Cd content in soils at different sampling scales

    表2 不同采樣尺度下土壤Cd半變異函數(shù)參數(shù)Table 2 Semi-variant function parameters of Cd in soils at different sampling scales

    2.2 土壤Cd濃度空間預(yù)測精度比較

    將4種網(wǎng)格尺度下土壤Cd濃度克里格插值進度交叉驗證(Cross Validation),結(jié)果見表3。由表3可見,總體而言4種采樣尺度下采樣網(wǎng)格越小,平均誤差(ME)絕對值越小。70 m×70 m和100 m×100 m 采樣尺度的均方根誤差(RMSE)和平均標(biāo)準誤差(ASE)較小且接近,表明70 m×70 m和100 m×100 m網(wǎng)格采樣精度接近,且能較好地滿足預(yù)測要求,預(yù)測值與真實值相近。

    表3 不同采樣尺度下克里格法插值精度估計Table 3 Estimation of Kriging interpolation accuracy at different sampling scales

    圖4為不同采樣網(wǎng)格尺度下克里格插值方法土壤Cd濃度的獨立驗證結(jié)果。由圖4可知,網(wǎng)格密度越小,克里格插值預(yù)測精度越高,且4種采樣密度下|RMSE-ASE|隨采樣網(wǎng)格尺度的增大而增大,說明采樣網(wǎng)格越小,插值預(yù)測的結(jié)果越接近真實值。200 m×200 m網(wǎng)格尺度下精度遠小于其他較小的網(wǎng)格尺度??傮w而言,采樣網(wǎng)格越大,土壤Cd的預(yù)測越平滑,對細節(jié)反映能力越弱,結(jié)合半方差函數(shù)分析及ME、RMSE、ASE 3項指標(biāo),100 m×100 m網(wǎng)格能較為準確地反映場地Cd濃度分布狀況。

    圖4 不同采樣網(wǎng)格尺度下土壤Cd濃度克里格 插值的獨立驗證結(jié)果Fig.4 Independent validation results of Kriging interpolation of Cd content in soils at different sampling grid scales

    2.3 污染范圍預(yù)測不確定性分析

    污染范圍界定得準確與否以及對預(yù)測不確定性區(qū)域的量化對場地的修復(fù)治理邊界以及相關(guān)決策的制定具有重要影響,污染范圍界定過大,將實際未污染的區(qū)域插值預(yù)測為污染區(qū)域,會過多地耗費人力、物力,加大修復(fù)成本,造成一定的經(jīng)濟負擔(dān);而污染范圍界定過小,將實際污染的區(qū)域插值預(yù)測為非污染區(qū)域,導(dǎo)致場地污染風(fēng)險不能有效降低或消除。

    利用污染預(yù)測標(biāo)準誤差圖,量化不同空間位置預(yù)測結(jié)果的不確定性。結(jié)果顯示,采樣點周邊位置的誤差通常很小,預(yù)測標(biāo)準誤差大的區(qū)域主要集中在樣點較為稀疏且有高值點的區(qū)域。為進一步明確預(yù)測結(jié)果的不確定性區(qū)域,將污染范圍預(yù)測圖和預(yù)測標(biāo)準誤差圖進行空間疊加計算,假定將預(yù)測值與誤差值之差大于0.6 mg/kg的范圍劃分為污染區(qū)域,二者之和小于0.6 mg/kg的范圍劃為未污染區(qū)域,其他為不確定性區(qū)域,劃定結(jié)果見圖5。由圖5可見,不確定區(qū)域主要分布在西南側(cè)和北側(cè)區(qū)域,這2個區(qū)域從20世紀60年代至21世紀初先后存在多家從事化工生產(chǎn)活動的企業(yè),因生產(chǎn)條件有限、排污過程的跑冒滴漏會導(dǎo)致地塊土壤受到污染,污染物在地下水水力作用的影響下遷移擴散,使得地塊周邊農(nóng)田土壤成片受到污染,且西南側(cè)和北側(cè)區(qū)域內(nèi)還存在池塘等使前期該區(qū)域采樣點較少,即不確定區(qū)域主要存在樣點稀疏且高值區(qū)域。

    圖5 土壤Cd污染調(diào)查點位分布Fig.5 Distribution of samples for investigation of Cd pollution in soils

    2.4 不確定性區(qū)域加密驗證

    通過預(yù)測標(biāo)準誤差圖劃定不確定性區(qū)域,按初步調(diào)查超標(biāo)點位鄰近原則,在不確定區(qū)域內(nèi)的初步調(diào)查超標(biāo)點位周邊選取最近點位作為詳細調(diào)查點位(圖5),對初步調(diào)查和詳細調(diào)查Cd濃度結(jié)果進行統(tǒng)計分析,結(jié)果見表4。從表4可知,100 m×100 m網(wǎng)格初步調(diào)查結(jié)果空間位置精確度為78.89%,對場地污染范圍預(yù)測不確定性影響較大的主要是高估區(qū)域,占比為14.23%,在不確定性區(qū)域加密布點以后,空間位置精確度提高到86.25%,提高了7.36個百分點,且污染高估區(qū)域占比明顯下降。對污染物濃度統(tǒng)計分析發(fā)現(xiàn),100 m×100 m網(wǎng)格初步調(diào)查與詳細調(diào)查結(jié)果差距不大,與所有樣點進行對比發(fā)現(xiàn),100 m×100 m網(wǎng)格Cd平均濃度略大于70 m×70 m網(wǎng)格,但二者誤差滿足HJ/T 166—2004《土壤環(huán)境監(jiān)測技術(shù)規(guī)范》所規(guī)定的25%范圍。圖6為土壤Cd污染范圍預(yù)測結(jié)果。由圖6可見,Cd污染主要集中在場地西側(cè),100 m×100 m網(wǎng)格加密前和加密后的差異主要集中在不確定性區(qū)域范圍內(nèi)。100 m×100 m 網(wǎng)格精確度有明顯提升,同時調(diào)查點位數(shù)量下降了35%。

    表4 土壤Cd統(tǒng)計特征Table 4 Statistical characteristics of Cd in soil heavy metals

    圖6 土壤Cd污染范圍預(yù)測結(jié)果Fig.6 Prediction results of Cd pollution range in soils

    本研究旨在探索科學(xué)合理的農(nóng)用地土壤污染調(diào)查最佳布點密度及點位優(yōu)化方法,在分析土壤污染物空間自相關(guān)性的基礎(chǔ)上,利用地統(tǒng)計學(xué)方法進行空間分布插值,基于預(yù)測誤差和污染評價標(biāo)準,劃定調(diào)查結(jié)果不確定性區(qū)域,結(jié)合初步調(diào)查超標(biāo)點位和地形、地貌等條件,在不確定性區(qū)域中加密布點,從而提高調(diào)查精度。該方法依據(jù)地統(tǒng)計學(xué)基本假設(shè),適用于空間分布具有中等相關(guān)性及以上的污染物,胡克林[35-37]等研究表明,表現(xiàn)出較明顯空間相關(guān)性,但對于空間自相關(guān)性較弱,僅在局部區(qū)域個別點位超標(biāo)的重金屬、多環(huán)芳烴等污染物,該方法的適用性需進一步驗證,這與謝云峰等[1]的研究結(jié)果類似。本研究采樣布點方法與目前概率加密布點法和模擬退火算法的各類加密方法相似[38-39],基于初步調(diào)查點位數(shù)據(jù)進行地統(tǒng)計學(xué)分析,判斷污染物是否具有空間自相關(guān)性,對于具有空間自相關(guān)性的污染地塊再進行參數(shù)化計算,從而優(yōu)化加密布點區(qū)域,做到既有效降低調(diào)查成本,又能提高調(diào)查精度。

    3 結(jié)論

    (1)對于采用基于模型化采樣布點的污染場地調(diào)查而言,通過獨立驗證和交叉驗證模型方法,能夠?qū)Σ蓸泳W(wǎng)格尺度進行有效優(yōu)化,本研究中100 m×100 m網(wǎng)格能夠在準確反映場地污染程度及范圍的前提下,降低采樣及分析成本。

    (2)基于空間插值預(yù)測標(biāo)準誤差與預(yù)測值進行運算,能合理確定場地污染的不確定性范圍,通過對不確定性區(qū)域進行針對性加密布點,能顯著提高預(yù)測精確度。

    猜你喜歡
    布點插值尺度
    2024年新增專業(yè)及布點高校名單
    財產(chǎn)的五大尺度和五重應(yīng)對
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    甘肅高校商科專業(yè)布點問題研究
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    江西省綠色通道車輛貨物檢測點布點方案探討
    9
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    久久精品久久久久久久性| 国产一区二区三区av在线| 亚洲国产精品专区欧美| 啦啦啦视频在线资源免费观看| 乱系列少妇在线播放| 熟妇人妻不卡中文字幕| 大香蕉久久网| av播播在线观看一区| 国产精品欧美亚洲77777| 日本欧美国产在线视频| 亚洲国产毛片av蜜桃av| 午夜福利在线观看免费完整高清在| 日韩伦理黄色片| 色视频在线一区二区三区| 99热全是精品| 熟女av电影| 中文字幕av成人在线电影| av又黄又爽大尺度在线免费看| 久久 成人 亚洲| 干丝袜人妻中文字幕| kizo精华| 黄片无遮挡物在线观看| 成人18禁高潮啪啪吃奶动态图 | 日本-黄色视频高清免费观看| 草草在线视频免费看| 久久99热这里只频精品6学生| 99热这里只有是精品在线观看| a级一级毛片免费在线观看| 国产女主播在线喷水免费视频网站| 各种免费的搞黄视频| 天天躁日日操中文字幕| 日韩不卡一区二区三区视频在线| 免费不卡的大黄色大毛片视频在线观看| 六月丁香七月| 国产色爽女视频免费观看| .国产精品久久| 在线观看一区二区三区| 午夜福利视频精品| 亚洲精品一区蜜桃| 日本-黄色视频高清免费观看| 成人国产av品久久久| www.色视频.com| 91久久精品国产一区二区成人| 久热久热在线精品观看| 女人久久www免费人成看片| 欧美bdsm另类| 能在线免费看毛片的网站| av线在线观看网站| 免费观看在线日韩| 97在线视频观看| 国产 精品1| av专区在线播放| 两个人的视频大全免费| 97在线视频观看| 狂野欧美激情性bbbbbb| 尾随美女入室| 在线观看一区二区三区| 亚洲精品国产av成人精品| 直男gayav资源| 欧美最新免费一区二区三区| 亚洲综合色惰| 国产精品嫩草影院av在线观看| 亚洲av成人精品一区久久| 免费观看在线日韩| 亚洲精品国产成人久久av| 最近2019中文字幕mv第一页| 蜜桃在线观看..| 日日啪夜夜撸| 久久6这里有精品| 久热久热在线精品观看| 99视频精品全部免费 在线| 美女国产视频在线观看| 成年人午夜在线观看视频| 男人添女人高潮全过程视频| 国产精品国产三级专区第一集| 国产亚洲最大av| 日韩 亚洲 欧美在线| h视频一区二区三区| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 国产精品久久久久久久久免| av天堂中文字幕网| 精品午夜福利在线看| 久久久色成人| 在线播放无遮挡| 国产一区二区三区av在线| 久久青草综合色| 纯流量卡能插随身wifi吗| 国产在线免费精品| 黑人高潮一二区| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 国产一区二区在线观看日韩| 中文字幕精品免费在线观看视频 | 2018国产大陆天天弄谢| 韩国高清视频一区二区三区| 国产黄片美女视频| 高清av免费在线| 熟女人妻精品中文字幕| 国产成人freesex在线| 美女主播在线视频| 99久国产av精品国产电影| 夜夜爽夜夜爽视频| 亚洲精品,欧美精品| 综合色丁香网| 中文在线观看免费www的网站| 全区人妻精品视频| 国产亚洲最大av| 亚洲av综合色区一区| 成人高潮视频无遮挡免费网站| 日韩不卡一区二区三区视频在线| 麻豆成人午夜福利视频| 国产伦理片在线播放av一区| 亚洲伊人久久精品综合| 乱系列少妇在线播放| 亚洲精华国产精华液的使用体验| 日韩制服骚丝袜av| 丝袜脚勾引网站| kizo精华| 高清日韩中文字幕在线| 午夜激情福利司机影院| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 国产精品久久久久成人av| 亚洲丝袜综合中文字幕| 一区在线观看完整版| av在线播放精品| 亚洲图色成人| 免费观看av网站的网址| 亚洲怡红院男人天堂| 97热精品久久久久久| 少妇熟女欧美另类| 国产伦在线观看视频一区| 免费播放大片免费观看视频在线观看| 简卡轻食公司| 女性生殖器流出的白浆| 黑人高潮一二区| 狂野欧美激情性bbbbbb| 三级经典国产精品| 欧美国产精品一级二级三级 | 在线精品无人区一区二区三 | 亚洲三级黄色毛片| 一区在线观看完整版| 丰满少妇做爰视频| 国产精品人妻久久久久久| 亚洲精品日韩av片在线观看| 久久久久网色| 亚洲最大成人中文| 国产亚洲欧美精品永久| 亚洲av综合色区一区| av天堂中文字幕网| 青春草亚洲视频在线观看| 如何舔出高潮| 成人毛片60女人毛片免费| 狂野欧美激情性xxxx在线观看| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看 | 亚洲精品成人av观看孕妇| 欧美人与善性xxx| 亚洲国产精品专区欧美| 美女高潮的动态| 亚洲成色77777| 在线观看三级黄色| 18禁动态无遮挡网站| 全区人妻精品视频| 久久青草综合色| 亚洲成色77777| 免费人妻精品一区二区三区视频| 最后的刺客免费高清国语| 久热这里只有精品99| 日本与韩国留学比较| 久久韩国三级中文字幕| 亚洲第一av免费看| 久久精品久久久久久噜噜老黄| videossex国产| 人妻 亚洲 视频| av免费观看日本| 亚洲欧美清纯卡通| 久久久久久久久久成人| 国产男女内射视频| 亚洲av在线观看美女高潮| 国产女主播在线喷水免费视频网站| 久久人人爽av亚洲精品天堂 | 色哟哟·www| 欧美日韩亚洲高清精品| 亚洲精品视频女| 美女福利国产在线 | 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 少妇高潮的动态图| 我要看日韩黄色一级片| 极品少妇高潮喷水抽搐| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| 亚洲丝袜综合中文字幕| 少妇人妻久久综合中文| 亚洲av综合色区一区| 亚洲国产高清在线一区二区三| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 精华霜和精华液先用哪个| 国产亚洲精品久久久com| 国产成人一区二区在线| 日本vs欧美在线观看视频 | 午夜免费男女啪啪视频观看| 最近最新中文字幕免费大全7| 99热全是精品| 身体一侧抽搐| 少妇猛男粗大的猛烈进出视频| 国产色婷婷99| 国产精品av视频在线免费观看| 欧美精品亚洲一区二区| 九九在线视频观看精品| 一级二级三级毛片免费看| 日日啪夜夜爽| 91aial.com中文字幕在线观看| 性色avwww在线观看| 欧美变态另类bdsm刘玥| 亚洲av成人精品一二三区| 99久国产av精品国产电影| 伊人久久精品亚洲午夜| 97在线视频观看| av线在线观看网站| 精华霜和精华液先用哪个| 国产老妇伦熟女老妇高清| 亚洲最大成人中文| 亚洲av电影在线观看一区二区三区| 国产深夜福利视频在线观看| 国产 一区精品| 在现免费观看毛片| 成人亚洲欧美一区二区av| 丰满迷人的少妇在线观看| 新久久久久国产一级毛片| 汤姆久久久久久久影院中文字幕| 麻豆成人午夜福利视频| 成人亚洲欧美一区二区av| 99九九线精品视频在线观看视频| 久久国产精品男人的天堂亚洲 | 精品亚洲成a人片在线观看 | 国产日韩欧美在线精品| 亚洲国产av新网站| 18禁动态无遮挡网站| 中文字幕av成人在线电影| 国产免费一区二区三区四区乱码| 欧美老熟妇乱子伦牲交| 精品久久国产蜜桃| 夜夜爽夜夜爽视频| 亚洲一级一片aⅴ在线观看| 免费看av在线观看网站| 亚洲精品视频女| 日本一二三区视频观看| 三级经典国产精品| 久久97久久精品| 少妇的逼水好多| 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美 | 青春草国产在线视频| 国产色爽女视频免费观看| 国产爱豆传媒在线观看| 我要看黄色一级片免费的| 国产欧美亚洲国产| 超碰av人人做人人爽久久| 国产淫语在线视频| 高清日韩中文字幕在线| 成年av动漫网址| 日日啪夜夜撸| 狂野欧美白嫩少妇大欣赏| 国产有黄有色有爽视频| 美女内射精品一级片tv| www.av在线官网国产| 精品一区在线观看国产| 大又大粗又爽又黄少妇毛片口| 最新中文字幕久久久久| 久久久久国产精品人妻一区二区| 国产精品一区二区性色av| 亚洲成人手机| 人妻一区二区av| 日本猛色少妇xxxxx猛交久久| 又爽又黄a免费视频| av播播在线观看一区| 下体分泌物呈黄色| 我的老师免费观看完整版| 91久久精品国产一区二区成人| 日韩在线高清观看一区二区三区| 一级a做视频免费观看| 新久久久久国产一级毛片| 精华霜和精华液先用哪个| 在线免费观看不下载黄p国产| 永久网站在线| 99久久精品国产国产毛片| www.av在线官网国产| 国产精品精品国产色婷婷| 国产成人精品一,二区| 亚洲av成人精品一二三区| 免费看av在线观看网站| 日本与韩国留学比较| 水蜜桃什么品种好| 高清黄色对白视频在线免费看 | 尤物成人国产欧美一区二区三区| 99热国产这里只有精品6| 久久精品久久久久久久性| 成人毛片a级毛片在线播放| 成人无遮挡网站| 亚洲久久久国产精品| 麻豆国产97在线/欧美| 国产 一区 欧美 日韩| 视频区图区小说| 免费看日本二区| 这个男人来自地球电影免费观看 | 卡戴珊不雅视频在线播放| 亚洲性久久影院| 精品久久久久久久久av| 大又大粗又爽又黄少妇毛片口| av免费在线看不卡| 丝瓜视频免费看黄片| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 在线观看免费日韩欧美大片 | 我要看黄色一级片免费的| 亚洲经典国产精华液单| 欧美zozozo另类| 男人添女人高潮全过程视频| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 自拍欧美九色日韩亚洲蝌蚪91 | 国产有黄有色有爽视频| 一本久久精品| 熟女电影av网| av国产久精品久网站免费入址| 大香蕉久久网| 国产成人a∨麻豆精品| 内射极品少妇av片p| 精品久久久久久久久av| 人妻 亚洲 视频| 少妇的逼好多水| 亚洲人成网站高清观看| av又黄又爽大尺度在线免费看| videos熟女内射| 精品一区二区免费观看| 精品亚洲成a人片在线观看 | 国产精品不卡视频一区二区| 成人特级av手机在线观看| 在线观看免费视频网站a站| 午夜激情久久久久久久| 国产高清三级在线| 大片电影免费在线观看免费| 成人综合一区亚洲| 国产精品三级大全| 少妇的逼水好多| 国产精品欧美亚洲77777| 国产欧美日韩精品一区二区| 成人特级av手机在线观看| 久久这里有精品视频免费| 久久国内精品自在自线图片| 国产毛片在线视频| 亚洲国产精品一区三区| 91狼人影院| 亚洲国产精品一区三区| 久久精品国产亚洲av涩爱| 综合色丁香网| 少妇高潮的动态图| 欧美成人一区二区免费高清观看| 在线天堂最新版资源| 日韩欧美精品免费久久| 欧美三级亚洲精品| 亚洲国产色片| 男女边摸边吃奶| 51国产日韩欧美| 欧美精品人与动牲交sv欧美| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 超碰97精品在线观看| 深夜a级毛片| 国产一级毛片在线| 日韩中字成人| 亚洲精品一区蜜桃| 亚洲国产精品专区欧美| 亚洲av成人精品一区久久| 久久久久视频综合| 中文字幕久久专区| 精品一区在线观看国产| 亚洲,欧美,日韩| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 偷拍熟女少妇极品色| 91狼人影院| 久久久久久久精品精品| h日本视频在线播放| 亚洲婷婷狠狠爱综合网| 国产有黄有色有爽视频| av女优亚洲男人天堂| 亚洲av欧美aⅴ国产| 在线免费十八禁| 久久热精品热| 色视频www国产| 国产真实伦视频高清在线观看| 国产精品久久久久久久久免| 日本wwww免费看| 高清在线视频一区二区三区| 一本色道久久久久久精品综合| 国产精品麻豆人妻色哟哟久久| 亚州av有码| 成人综合一区亚洲| 搡老乐熟女国产| 高清视频免费观看一区二区| av在线播放精品| 国产精品.久久久| 一区二区三区免费毛片| 大话2 男鬼变身卡| 国产成人aa在线观看| 男女边摸边吃奶| 天堂俺去俺来也www色官网| 久久国产乱子免费精品| 九九爱精品视频在线观看| 国产成人91sexporn| 一个人免费看片子| 一边亲一边摸免费视频| 亚洲精品一二三| 一级毛片久久久久久久久女| 大香蕉久久网| 日本一二三区视频观看| 天天躁夜夜躁狠狠久久av| 日韩一区二区三区影片| 国产一区二区在线观看日韩| 国产亚洲av片在线观看秒播厂| 国产男人的电影天堂91| 国产成人午夜福利电影在线观看| 午夜免费鲁丝| av天堂中文字幕网| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠久久av| 麻豆精品久久久久久蜜桃| 日本色播在线视频| 精品亚洲成a人片在线观看 | 亚洲国产日韩一区二区| 一区二区三区乱码不卡18| 制服丝袜香蕉在线| 久久精品国产亚洲网站| 在线观看一区二区三区| 久久久久久久亚洲中文字幕| 欧美国产精品一级二级三级 | 欧美日韩在线观看h| 欧美 日韩 精品 国产| 日韩不卡一区二区三区视频在线| 国产日韩欧美在线精品| 我要看日韩黄色一级片| 女性生殖器流出的白浆| 国内精品宾馆在线| 亚洲真实伦在线观看| 中文字幕制服av| 国产伦精品一区二区三区视频9| 老女人水多毛片| 亚洲国产欧美在线一区| 日韩成人伦理影院| 午夜福利视频精品| 人妻夜夜爽99麻豆av| 久久精品夜色国产| 日韩伦理黄色片| 久久久a久久爽久久v久久| 国产亚洲午夜精品一区二区久久| 欧美+日韩+精品| 嫩草影院入口| 交换朋友夫妻互换小说| 美女视频免费永久观看网站| 成人国产av品久久久| 久久99热6这里只有精品| 美女主播在线视频| 亚洲国产毛片av蜜桃av| 久久久久人妻精品一区果冻| 国产成人精品福利久久| 在现免费观看毛片| 欧美精品一区二区免费开放| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 久久av网站| 一二三四中文在线观看免费高清| av.在线天堂| 欧美xxxx性猛交bbbb| 观看av在线不卡| 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 十八禁网站网址无遮挡 | 尤物成人国产欧美一区二区三区| 色哟哟·www| 伊人久久精品亚洲午夜| 久久国内精品自在自线图片| 国产精品熟女久久久久浪| 精品人妻偷拍中文字幕| 欧美日韩亚洲高清精品| 在线观看三级黄色| 欧美激情国产日韩精品一区| 欧美日韩国产mv在线观看视频 | 国产一级毛片在线| 各种免费的搞黄视频| 国产成人午夜福利电影在线观看| 一二三四中文在线观看免费高清| 成年美女黄网站色视频大全免费 | 国产精品久久久久成人av| freevideosex欧美| 久久精品国产a三级三级三级| 九色成人免费人妻av| 伊人久久精品亚洲午夜| 亚洲av免费高清在线观看| 一区二区三区精品91| 蜜桃久久精品国产亚洲av| 亚洲av福利一区| 欧美+日韩+精品| 插阴视频在线观看视频| 欧美亚洲 丝袜 人妻 在线| 国产深夜福利视频在线观看| 国产精品三级大全| 女性生殖器流出的白浆| 乱系列少妇在线播放| 国产一区二区在线观看日韩| 国产精品免费大片| 色综合色国产| 亚洲欧美中文字幕日韩二区| 婷婷色综合大香蕉| 美女中出高潮动态图| 亚洲精品久久午夜乱码| 亚洲第一区二区三区不卡| 麻豆精品久久久久久蜜桃| 日本vs欧美在线观看视频 | 高清黄色对白视频在线免费看 | 免费大片黄手机在线观看| 国产精品嫩草影院av在线观看| 欧美日韩在线观看h| 又黄又爽又刺激的免费视频.| 亚洲国产精品一区三区| 熟女电影av网| 国产69精品久久久久777片| 日韩大片免费观看网站| 人妻一区二区av| 亚洲精品国产av成人精品| 久久99精品国语久久久| 亚洲国产精品国产精品| 狂野欧美激情性xxxx在线观看| av天堂中文字幕网| 麻豆精品久久久久久蜜桃| 七月丁香在线播放| 国产成人一区二区在线| 亚洲四区av| 黄片无遮挡物在线观看| 在线观看一区二区三区激情| 人人妻人人添人人爽欧美一区卜 | 观看美女的网站| 成年美女黄网站色视频大全免费 | 国产 精品1| 成人国产麻豆网| 国产在线男女| 成人无遮挡网站| 亚洲真实伦在线观看| 老女人水多毛片| 少妇裸体淫交视频免费看高清| 亚洲电影在线观看av| 免费看av在线观看网站| 国产亚洲91精品色在线| 亚洲最大成人中文| 丝袜脚勾引网站| 搡女人真爽免费视频火全软件| 午夜免费观看性视频| 嫩草影院入口| 这个男人来自地球电影免费观看 | 干丝袜人妻中文字幕| 精品午夜福利在线看| 九草在线视频观看| 狂野欧美激情性bbbbbb| 亚洲av.av天堂| 欧美+日韩+精品| 久久久久网色| 一区二区av电影网| 午夜免费男女啪啪视频观看| 一本久久精品| 亚洲天堂av无毛| 精品视频人人做人人爽| 久久韩国三级中文字幕| 汤姆久久久久久久影院中文字幕| 最近最新中文字幕大全电影3| 天美传媒精品一区二区| 纵有疾风起免费观看全集完整版| 久久久久久久国产电影| 日韩强制内射视频| 国产有黄有色有爽视频| 国产v大片淫在线免费观看| 国产亚洲午夜精品一区二区久久| 久久久久久久久久人人人人人人| 美女高潮的动态| 久久久国产一区二区| 亚洲成色77777| 国产亚洲一区二区精品| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 青春草国产在线视频| 在线观看三级黄色| 一个人看的www免费观看视频| 麻豆精品久久久久久蜜桃| 777米奇影视久久| 成人亚洲精品一区在线观看 | 国产精品一区二区性色av| 免费黄网站久久成人精品| 搡老乐熟女国产| 大话2 男鬼变身卡| 国产久久久一区二区三区| 青春草亚洲视频在线观看| 国产精品一区二区性色av| 亚洲图色成人| 少妇裸体淫交视频免费看高清| 男的添女的下面高潮视频| 国产日韩欧美亚洲二区| 免费观看性生交大片5| 亚洲精品日韩在线中文字幕| 少妇精品久久久久久久| 卡戴珊不雅视频在线播放| 亚洲精品,欧美精品| 久久精品熟女亚洲av麻豆精品|