• <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| 国产精品麻豆人妻色哟哟久久 | 亚洲高清免费不卡视频| 国产白丝娇喘喷水9色精品| 一卡2卡三卡四卡精品乱码亚洲| 日日摸夜夜添夜夜爱| 亚洲av成人精品一二三区| 免费无遮挡裸体视频| 国产精品永久免费网站| 亚洲国产色片| 欧美一区二区亚洲| 在线观看一区二区三区| 国产黄片美女视频| 国产亚洲精品av在线| 午夜久久久久精精品| 99热全是精品| 一个人观看的视频www高清免费观看| 变态另类丝袜制服| av免费在线看不卡| 中文精品一卡2卡3卡4更新| 国产午夜福利久久久久久| 午夜视频国产福利| 亚洲av不卡在线观看| 亚洲国产精品成人久久小说| 好男人视频免费观看在线| av国产免费在线观看| 人体艺术视频欧美日本| 欧美潮喷喷水| 成人鲁丝片一二三区免费| 毛片一级片免费看久久久久| 亚洲熟妇中文字幕五十中出| 国产亚洲5aaaaa淫片| 国产精品久久久久久久电影| 国产激情偷乱视频一区二区| 极品教师在线视频| 日韩亚洲欧美综合| 国产黄a三级三级三级人| 国产 一区精品| 国产午夜福利久久久久久| 搞女人的毛片| 亚洲三级黄色毛片| 一卡2卡三卡四卡精品乱码亚洲| 禁无遮挡网站| av在线播放精品| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久亚洲| 亚洲av日韩在线播放| 人人妻人人澡人人爽人人夜夜 | 免费看a级黄色片| 网址你懂的国产日韩在线| 久99久视频精品免费| 亚洲国产欧美人成| 日韩欧美 国产精品| 国产精品99久久久久久久久| 欧美97在线视频| 亚洲av不卡在线观看| 国产美女午夜福利| 免费不卡的大黄色大毛片视频在线观看 | 哪个播放器可以免费观看大片| 插逼视频在线观看| 亚洲精华国产精华液的使用体验| 国产av一区在线观看免费| 欧美成人精品欧美一级黄| 国产乱来视频区| 成人av在线播放网站| 插阴视频在线观看视频| 毛片一级片免费看久久久久| 亚洲四区av| 久99久视频精品免费| 久热久热在线精品观看| 能在线免费观看的黄片| 97热精品久久久久久| 插阴视频在线观看视频| 国产精品精品国产色婷婷| 免费不卡的大黄色大毛片视频在线观看 | 麻豆乱淫一区二区| 欧美成人一区二区免费高清观看| 超碰av人人做人人爽久久| 99在线人妻在线中文字幕| 男女下面进入的视频免费午夜| 天堂中文最新版在线下载 | 亚洲成色77777| 国产伦精品一区二区三区四那| 99热这里只有是精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 美女脱内裤让男人舔精品视频| 国产探花在线观看一区二区| 久久99热6这里只有精品| 91精品伊人久久大香线蕉| 亚洲av免费在线观看| 哪个播放器可以免费观看大片| 日本免费在线观看一区| 成人高潮视频无遮挡免费网站| 久久久a久久爽久久v久久| 最近最新中文字幕大全电影3| 欧美成人精品欧美一级黄| 秋霞在线观看毛片| 大香蕉97超碰在线| 日日摸夜夜添夜夜爱| 欧美性猛交╳xxx乱大交人| 青春草亚洲视频在线观看| 日韩精品青青久久久久久| 亚洲经典国产精华液单| 中文资源天堂在线| 国产又色又爽无遮挡免| 亚洲在线观看片| 久久久国产成人精品二区| 久久久久久国产a免费观看| 不卡视频在线观看欧美| 久99久视频精品免费| av免费在线看不卡| 白带黄色成豆腐渣| 国产美女午夜福利| 国产精品人妻久久久影院| 国产黄片美女视频| 高清在线视频一区二区三区 | 亚洲av一区综合| 淫秽高清视频在线观看| 长腿黑丝高跟| 久久精品影院6| 亚洲精品国产成人久久av| 男人和女人高潮做爰伦理| 内地一区二区视频在线| 啦啦啦韩国在线观看视频| kizo精华| 亚洲欧美成人精品一区二区| 免费av观看视频| 亚洲人成网站在线观看播放| 小蜜桃在线观看免费完整版高清| 成人一区二区视频在线观看| 丝袜喷水一区| 国产日韩欧美在线精品| 国产成人精品一,二区| 中文字幕人妻熟人妻熟丝袜美| 高清在线视频一区二区三区 | 欧美变态另类bdsm刘玥| 精品国内亚洲2022精品成人| 亚洲美女搞黄在线观看| 欧美高清性xxxxhd video| 中文字幕亚洲精品专区| 精品一区二区免费观看| av.在线天堂| 小说图片视频综合网站| 伦理电影大哥的女人| 欧美不卡视频在线免费观看| 国产白丝娇喘喷水9色精品| 精品酒店卫生间| www.av在线官网国产| 少妇高潮的动态图| 精品酒店卫生间| 欧美激情国产日韩精品一区| 亚洲欧美精品综合久久99| 三级男女做爰猛烈吃奶摸视频| 亚洲精品久久久久久婷婷小说 | 最近最新中文字幕免费大全7| 三级男女做爰猛烈吃奶摸视频| 亚洲成色77777| 黄色配什么色好看| 伊人久久精品亚洲午夜| 亚洲无线观看免费| 日韩欧美在线乱码| 久久久久免费精品人妻一区二区| 国产精品野战在线观看| 免费观看精品视频网站| 青春草视频在线免费观看| 成人一区二区视频在线观看| 国产成人午夜福利电影在线观看| 国产伦在线观看视频一区| 国产在线男女| 中文资源天堂在线| 麻豆成人av视频| 国产白丝娇喘喷水9色精品| 国产精品一区二区在线观看99 | 国产精品美女特级片免费视频播放器| 精品欧美国产一区二区三| 在现免费观看毛片| 真实男女啪啪啪动态图| 成人亚洲欧美一区二区av| 日本wwww免费看| 成人国产麻豆网| 国产不卡一卡二| 少妇丰满av| 国产黄色小视频在线观看| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| videos熟女内射| 成人美女网站在线观看视频| av在线老鸭窝| 久久久a久久爽久久v久久| 精品国内亚洲2022精品成人| 美女被艹到高潮喷水动态| 国产午夜福利久久久久久| 国产精品无大码| 久久国产乱子免费精品| 最新中文字幕久久久久| 国产精品av视频在线免费观看| www日本黄色视频网| 久久99热这里只有精品18| www.色视频.com| 草草在线视频免费看| 亚洲18禁久久av| 我要看日韩黄色一级片| 狠狠狠狠99中文字幕| 成人无遮挡网站| 春色校园在线视频观看| 国产精品嫩草影院av在线观看| 国内精品宾馆在线| 韩国av在线不卡| 欧美另类亚洲清纯唯美| 国产综合懂色| 少妇的逼水好多| 人妻制服诱惑在线中文字幕| av.在线天堂| 久久欧美精品欧美久久欧美| 又粗又爽又猛毛片免费看| 久久99蜜桃精品久久| 日韩中字成人| 午夜a级毛片| 国产精品一区二区三区四区久久| 超碰av人人做人人爽久久| 亚洲欧美清纯卡通| 中文字幕精品亚洲无线码一区| 极品教师在线视频| 欧美潮喷喷水| 国产精品1区2区在线观看.| 三级国产精品欧美在线观看| 久久精品国产自在天天线| 国产一区亚洲一区在线观看| 一个人免费在线观看电影| 春色校园在线视频观看| 干丝袜人妻中文字幕| 99热这里只有是精品50| 日韩国内少妇激情av| 亚洲电影在线观看av| 国产一区二区三区av在线| 国产成年人精品一区二区| 99九九线精品视频在线观看视频| 国产精品久久视频播放| 精品人妻熟女av久视频| 女的被弄到高潮叫床怎么办| 天堂网av新在线| 伊人久久精品亚洲午夜| 夜夜看夜夜爽夜夜摸| 久久99精品国语久久久| 国产一区亚洲一区在线观看| 国产黄片视频在线免费观看| 国产淫片久久久久久久久| 免费看日本二区| 青春草视频在线免费观看| 国内精品一区二区在线观看| 亚洲第一区二区三区不卡| 亚洲成人中文字幕在线播放| 久久久久久久午夜电影| 一区二区三区乱码不卡18| 久久精品国产99精品国产亚洲性色| 少妇的逼水好多| 亚洲中文字幕一区二区三区有码在线看| 久久精品久久久久久噜噜老黄 | 别揉我奶头 嗯啊视频| 亚洲性久久影院| 又粗又爽又猛毛片免费看| 在线免费观看的www视频| 成人午夜高清在线视频| 在线免费观看不下载黄p国产| 又黄又爽又刺激的免费视频.| 少妇高潮的动态图| 亚洲欧美日韩无卡精品| 国产国拍精品亚洲av在线观看| 毛片一级片免费看久久久久| 女的被弄到高潮叫床怎么办| 日韩成人伦理影院| 寂寞人妻少妇视频99o| 日本wwww免费看| 91久久精品国产一区二区成人| 韩国高清视频一区二区三区| 欧美潮喷喷水| 高清午夜精品一区二区三区| www日本黄色视频网| 超碰av人人做人人爽久久| 狂野欧美白嫩少妇大欣赏| 久久久精品大字幕| 欧美日韩一区二区视频在线观看视频在线 | av播播在线观看一区| 久久精品综合一区二区三区| 久久久色成人| 亚洲av成人精品一区久久| 欧美人与善性xxx| 两个人的视频大全免费| 久久久成人免费电影| av在线播放精品| 国产亚洲一区二区精品| 亚洲人成网站在线观看播放| 日本黄大片高清| 国产成人免费观看mmmm| 亚洲av一区综合| 久久人人爽人人片av| 久久人妻av系列| 天天一区二区日本电影三级| 波多野结衣高清无吗| 国产成年人精品一区二区| 亚洲国产最新在线播放| 国产一级毛片在线| 亚洲国产日韩欧美精品在线观看| 日本黄大片高清| 日韩强制内射视频| 亚洲最大成人手机在线| 免费黄网站久久成人精品| 欧美人与善性xxx| 日本五十路高清| 男女边吃奶边做爰视频| 久久久久网色| 成人美女网站在线观看视频| 在线a可以看的网站| 国内揄拍国产精品人妻在线| 日本av手机在线免费观看| 欧美3d第一页| 久久久久久久久中文| 寂寞人妻少妇视频99o| 国产成人免费观看mmmm| 精品久久久久久久末码| 久久久久久久国产电影| 日本熟妇午夜| 精品人妻偷拍中文字幕| 国产亚洲午夜精品一区二区久久 | 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 欧美性猛交黑人性爽| 日韩亚洲欧美综合| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩东京热| 精品人妻偷拍中文字幕| 高清av免费在线| 中文乱码字字幕精品一区二区三区 | 一级毛片电影观看 | 日韩欧美国产在线观看| 久久久久久久久久久免费av| av免费观看日本| 国产高清有码在线观看视频| 国产麻豆成人av免费视频| 啦啦啦观看免费观看视频高清| 亚洲第一区二区三区不卡| 最近最新中文字幕免费大全7| 能在线免费看毛片的网站| 亚洲精品乱码久久久v下载方式| 成人高潮视频无遮挡免费网站| 国产精品福利在线免费观看| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 久久精品熟女亚洲av麻豆精品 | 久久精品国产鲁丝片午夜精品| 熟妇人妻久久中文字幕3abv| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 精品人妻熟女av久视频| 赤兔流量卡办理| 亚洲真实伦在线观看| 国产精品一区二区在线观看99 | 麻豆一二三区av精品| 1000部很黄的大片| 色哟哟·www| 久久久色成人| 啦啦啦韩国在线观看视频| 国产精品久久久久久久电影| 婷婷六月久久综合丁香| 亚洲欧美日韩卡通动漫| 欧美高清性xxxxhd video| .国产精品久久| av免费在线看不卡| 亚洲最大成人手机在线| 成人综合一区亚洲| 国产乱人视频| 成人漫画全彩无遮挡| 一级av片app| 黄色配什么色好看| 国产成人91sexporn| 亚洲欧美清纯卡通| 午夜精品一区二区三区免费看| 全区人妻精品视频| 国产亚洲91精品色在线| 狂野欧美白嫩少妇大欣赏| 亚洲成av人片在线播放无| 国产 一区精品| 99久久无色码亚洲精品果冻| 高清在线视频一区二区三区 | eeuss影院久久| 久久午夜福利片| 高清在线视频一区二区三区 | 99久久精品热视频| 禁无遮挡网站| 精品久久久久久久人妻蜜臀av| 自拍偷自拍亚洲精品老妇| 精品免费久久久久久久清纯| 日韩欧美国产在线观看| 人妻夜夜爽99麻豆av| 噜噜噜噜噜久久久久久91| 日本黄色片子视频| 人妻制服诱惑在线中文字幕| 国产伦精品一区二区三区视频9| 国模一区二区三区四区视频| 老女人水多毛片| 久久热精品热| 校园人妻丝袜中文字幕| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 综合色av麻豆| 欧美丝袜亚洲另类| 日日摸夜夜添夜夜添av毛片| 美女大奶头视频| 18禁在线无遮挡免费观看视频| 日韩欧美国产在线观看| 91久久精品电影网| 一级毛片久久久久久久久女| 国产91av在线免费观看| av天堂中文字幕网| 国产av不卡久久| 99热这里只有是精品在线观看| 国产精品一区二区三区四区免费观看| 我的女老师完整版在线观看| 人妻少妇偷人精品九色| 国产精品99久久久久久久久| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩欧美精品在线观看| eeuss影院久久| 国产精华一区二区三区| 亚洲人成网站在线观看播放| 免费播放大片免费观看视频在线观看 | 久久精品影院6| 别揉我奶头 嗯啊视频| h日本视频在线播放| 亚洲av.av天堂| 能在线免费观看的黄片| 99久久人妻综合| 国产成人午夜福利电影在线观看| 亚洲在久久综合| av在线亚洲专区| 午夜视频国产福利| 午夜免费男女啪啪视频观看| 秋霞伦理黄片| 亚洲丝袜综合中文字幕| 天堂网av新在线| 极品教师在线视频| 久久久久久久久中文| 九九在线视频观看精品| 国产精品福利在线免费观看| 国国产精品蜜臀av免费| 亚洲一级一片aⅴ在线观看| 校园人妻丝袜中文字幕| 最近最新中文字幕大全电影3| 日韩精品青青久久久久久| 91精品伊人久久大香线蕉| 亚洲丝袜综合中文字幕| 久久99蜜桃精品久久| videossex国产| 99久久精品一区二区三区| av国产免费在线观看| 亚洲av免费高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久大av| 三级毛片av免费| 波野结衣二区三区在线| 91久久精品国产一区二区三区| 欧美极品一区二区三区四区| 亚洲精品aⅴ在线观看| av视频在线观看入口| 中文字幕久久专区| 91aial.com中文字幕在线观看| 国产精华一区二区三区| 精品少妇黑人巨大在线播放 | 亚洲人成网站在线播| 日韩欧美 国产精品| 国产精品久久久久久av不卡| 在线观看美女被高潮喷水网站| 观看免费一级毛片| 2022亚洲国产成人精品| 一本久久精品| 日韩国内少妇激情av| 国产一区亚洲一区在线观看| 婷婷色av中文字幕| 国产一级毛片七仙女欲春2| 搡老妇女老女人老熟妇| 亚洲丝袜综合中文字幕| 精品国产三级普通话版| 热99re8久久精品国产| 色网站视频免费| 中国国产av一级| 欧美97在线视频| 亚洲美女视频黄频| 久久人人爽人人片av| 全区人妻精品视频| 嫩草影院新地址| 国产成人精品一,二区| 久久99蜜桃精品久久| 1024手机看黄色片| 三级经典国产精品| 中文字幕熟女人妻在线| 十八禁国产超污无遮挡网站| 精品久久久久久久久亚洲| 国产精品美女特级片免费视频播放器| 国产精品1区2区在线观看.| 免费大片18禁| 联通29元200g的流量卡| www.色视频.com| 国产精品伦人一区二区| 日韩精品青青久久久久久| 禁无遮挡网站| 国产白丝娇喘喷水9色精品| 综合色丁香网| 欧美潮喷喷水| 亚洲一级一片aⅴ在线观看| 一卡2卡三卡四卡精品乱码亚洲| av在线观看视频网站免费| 亚洲欧美日韩高清专用| 亚洲色图av天堂| 99九九线精品视频在线观看视频| 国产免费又黄又爽又色| 久久这里只有精品中国| 一级av片app| 日产精品乱码卡一卡2卡三| 中文乱码字字幕精品一区二区三区 | 1024手机看黄色片| 三级男女做爰猛烈吃奶摸视频| 岛国在线免费视频观看| 亚洲性久久影院| 我要看日韩黄色一级片| 秋霞伦理黄片| 亚洲自偷自拍三级| 亚洲天堂国产精品一区在线| 3wmmmm亚洲av在线观看| 亚洲一区高清亚洲精品| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 高清午夜精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 激情 狠狠 欧美| 日本一本二区三区精品| 国产精品一及| 极品教师在线视频| 国产精品国产三级专区第一集| 国产麻豆成人av免费视频| 亚洲中文字幕一区二区三区有码在线看| 六月丁香七月| 看免费成人av毛片| 日本免费a在线| 人人妻人人看人人澡| 国产男人的电影天堂91| 男女那种视频在线观看| 最近最新中文字幕免费大全7| 久久精品国产亚洲av天美| 少妇丰满av| 极品教师在线视频| 成人三级黄色视频| 成人二区视频| 村上凉子中文字幕在线| 国产亚洲一区二区精品| 只有这里有精品99| 久久99热6这里只有精品| 久久精品人妻少妇| av女优亚洲男人天堂| 身体一侧抽搐| 好男人在线观看高清免费视频| 永久网站在线| 午夜免费激情av| 边亲边吃奶的免费视频| 国产精品嫩草影院av在线观看| 亚洲电影在线观看av| 久久韩国三级中文字幕| 精品酒店卫生间| 99久久无色码亚洲精品果冻| 一个人看的www免费观看视频| 69av精品久久久久久| 大香蕉久久网| 99久久无色码亚洲精品果冻| 国产黄片美女视频| 六月丁香七月| 欧美区成人在线视频| 午夜老司机福利剧场| 日日摸夜夜添夜夜添av毛片| 嫩草影院精品99| 国产伦一二天堂av在线观看| 久久久久久大精品| 亚洲国产精品合色在线| 男插女下体视频免费在线播放| 成人毛片a级毛片在线播放| 免费在线观看成人毛片| 成人国产麻豆网| 中文字幕久久专区| 亚洲最大成人av| 中文资源天堂在线| 日本免费a在线| 一区二区三区乱码不卡18| 一区二区三区四区激情视频| av天堂中文字幕网| 男人的好看免费观看在线视频| 寂寞人妻少妇视频99o| 国产真实乱freesex| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 欧美日韩一区二区视频在线观看视频在线 | 成人漫画全彩无遮挡| 久久精品国产亚洲av涩爱| 欧美zozozo另类| 国产高清不卡午夜福利| 午夜福利在线观看吧| 亚洲电影在线观看av| 国产黄色视频一区二区在线观看 | 亚洲av不卡在线观看| 国产一区亚洲一区在线观看| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区 | 国产精品无大码| 精品久久久久久成人av|