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

    土壤中重金屬累積特征及生態(tài)風險評價①
    ——以烏拉特后旗有色金屬冶煉企業(yè)集中區(qū)為例

    2016-07-19 07:33:10榮,史銳*,張
    土壤 2016年2期
    關鍵詞:土壤重金屬成因

    岳 榮,史 銳*,張 紅

    (1 巴彥淖爾市環(huán)境科學研究所,內(nèi)蒙古臨河 015000;2 山西大學環(huán)境與資源學院,太原 030006)

    ?

    土壤中重金屬累積特征及生態(tài)風險評價①
    ——以烏拉特后旗有色金屬冶煉企業(yè)集中區(qū)為例

    岳榮1,史銳1*,張紅2

    (1 巴彥淖爾市環(huán)境科學研究所,內(nèi)蒙古臨河015000;2 山西大學環(huán)境與資源學院,太原030006)

    摘要:在內(nèi)蒙古烏拉特后旗有色金屬冶煉企業(yè)集中區(qū)域內(nèi)土壤中測定了 7種重金屬元素的含量,利用地統(tǒng)計學方法、地累積指數(shù)法及Hankanson指數(shù)法研究了土壤中重金屬元素的空間變異特征及成因,綜合評價了重金屬環(huán)境風險。結(jié)果表明:研究區(qū)內(nèi)重金屬元素含量已經(jīng)很大程度上受到了工業(yè)、農(nóng)業(yè)生產(chǎn)活動的影響,結(jié)構(gòu)性因素的影響已經(jīng)不顯著;重金屬元素的地累積指數(shù)均值排列順序為 Cd>Zn>Pb>Cu>As>Hg>Cr,100% 的土壤樣本表現(xiàn)出極強的 Cd污染程度,7種重金屬的潛在生態(tài)風險由強至弱依次為:Cd>Cu>Pb>As>Zn>Hg>Cr,研究區(qū)土壤總體處于輕微風險程度;研究區(qū)重金屬含量的累積成因各有不同,總體來說由區(qū)域的有色金屬冶煉、城鎮(zhèn)交通運輸、燃煤和耕地區(qū)域內(nèi)的污水灌溉及農(nóng)業(yè)施肥等導致。

    關鍵詞:土壤重金屬;空間變異;生態(tài)風險評價;成因

    土壤質(zhì)量變化是土壤學家和環(huán)境學家十分關注的問題,隨著城市化、工業(yè)化和農(nóng)業(yè)集約化的快速發(fā)展, 我國土壤環(huán)境和由此導致的糧食產(chǎn)量、食品安全問題越來越突出,這不僅制約了我國經(jīng)濟的快速發(fā)展,還嚴重影響著人類健康,因此,高強度經(jīng)濟開發(fā)所伴隨的土壤環(huán)境質(zhì)量變化問題成為急需研究的課題。重金屬是典型的土壤污染物,具有隱蔽性、難降解、移動性差和易被富集等特點[1-4],不但可影響土壤微生物區(qū)系和微生物過程,從而影響生態(tài)系統(tǒng)的結(jié)構(gòu)與功能,還可以在土壤溶液的作用下發(fā)生吸附、絮凝,并易在環(huán)境條件發(fā)生變化時懸浮和溶出形成二次污染,造成在土壤中濃度水平偏高,最終通過食物鏈在人類體內(nèi)富集。目前國內(nèi)外對土壤中重金屬的研究主要為土壤重金屬累積特征及其環(huán)境行為等[5-8],但多數(shù)研究是基于單個樣點的點過程研究,直接從宏觀角度利用空間采樣和空間模型研究較大尺度重金屬污染的面過程很少。而且,現(xiàn)行的土壤重金屬綜合污染評價方法主要為內(nèi)梅羅污染指數(shù)評價法和綜合污染指數(shù)評價法,它們用比較明確的界線對土壤重金屬污染程度加以區(qū)分和量化,而實際上土壤重金屬的污染狀況是漸變的,因此,在研究區(qū)域內(nèi)網(wǎng)格化布點可以直觀地描述土壤重金屬污染狀況的漸變性,使評價結(jié)果更加準確可靠。本研究采用網(wǎng)格化布點來分析描述土壤重金屬含量的空間特征,評價其污染程度及生態(tài)風險,這對預測重金屬污染風險,探明重金屬污染場地的污染來源及制定防治措施具有十分重要的意義[9-10]。

    內(nèi)蒙古烏拉特后旗是內(nèi)蒙古自治區(qū)兩個較大的有色金屬采選冶煉基地之一,有著悠久的有色金屬開采冶煉歷史。研究區(qū)域內(nèi)交通發(fā)達,人口密集,農(nóng)田較為集中,城鎮(zhèn)集中區(qū)及大量的有色金屬開采和冶煉企業(yè)分布于陰山南側(cè),長期的開采和冶煉產(chǎn)生大量的廢水、廢氣和廢渣,對烏拉特后旗內(nèi)陰山一線區(qū)域土壤質(zhì)量造成了巨大的影響。目前,文獻報道主要針對有色金屬較為集中的平原區(qū)域,針對沿山區(qū)域罕有報道。

    1 研究區(qū)概況及布點采樣

    研究區(qū)位于內(nèi)蒙古自治區(qū)烏拉特后旗陰山南側(cè)青山工業(yè)園區(qū)及周邊(圖 1),涵蓋了城鎮(zhèn)建成區(qū)、農(nóng)業(yè)耕作區(qū)、工業(yè)園區(qū)(主要為有色金屬冶煉區(qū))等多種景觀格局。區(qū)域內(nèi)有色金屬冶煉企業(yè)均有 10年以上的采冶歷史,是典型的重金屬污染場地。根據(jù)當?shù)氐耐恋乩妙愋?,按照多樣化土壤性質(zhì)類別與多樣化立地植被原則選取樣點,在研究區(qū)域按照網(wǎng)格布點原則共布設 61 個采樣點,樣點的分布見圖 2。2014 年 10月在研究區(qū)域內(nèi),采用 GPS 定位采集地表 0 ~ 20 cm土層的土壤樣品,每個樣品均由鄰近相同類型的 5 個點的土壤樣品均勻混合而成,土壤經(jīng)登記編號后,用聚氯乙烯塑料袋封裝,風干后剔除植物殘體和石塊,磨碎、過 200目篩,混勻、保存?zhèn)溆谩?/p>

    圖1 研究區(qū)地理位置及采樣、工礦企業(yè)布點圖Fig. 1 The sampling sites, mining enterprises and location of research area

    2 研究方法

    2.1測試方法

    土壤重金屬元素 As 與 Hg 采用北京海光儀器廠AFS-3100 原子熒光儀測定,Cu、Cd 與 Zn 采用德國耶拿公司 NovAA400P 原子吸收分光光度計測定,Cr 及 Pb 采用石墨爐原子吸收分光光度法測定。首先,將粒徑 ≤74 μm 的樣品在 105℃ 條件下烘干 2 h,然后稱取 4.0 g 樣品裝入模具中,采用低壓聚乙烯粉鑲邊墊底,在 30 t 的壓力下制成直徑 32 mm、鑲邊外徑 40 mm 的樣片,放入干燥器中待測。采用GBW07423-GBW07430 (土壤)國家一級標準物質(zhì)作為標準樣品。方法檢出限采用校準曲線法計算,首先測量多個低含量試樣并校準曲線,計算出含量為零時的射線熒光強度,并計算各組分的檢出限及分析結(jié)果,結(jié)果見表 1,由表 1 可知,本研究的重金屬分析方法的檢出限和準確度能滿足土壤樣品的分析要求。

    表1 土壤重金屬元素含量分析方法標樣比對分析結(jié)果Table 1 Detection limits and precision of the test methods of heavy metals in the soil of research area

    2.2空間變異特性研究方法

    半變異函數(shù)是分析變量空間結(jié)構(gòu)的主要工具,用以反映變量的空間自相關性[11]。 根據(jù)公式(1)計算出半變異函數(shù) ()hγ的散點圖, 并用不同類型的模型進行擬合, 計算其參數(shù)及擬合度最好的模型類型,計算方法如下:

    式中: Z( Xi)為i點位的重金屬含量觀測值,N( h)為由距離h分割觀測點的對數(shù),h為步長。

    2.3綜合評價方法

    采用地累積指數(shù)法[12](式2)和Hankanson潛在生態(tài)風險指數(shù)法[13](式3),評價研究區(qū)重金屬污染風險。計算公式分別為:

    式中:nC為樣品中n元素含量的實測值;nB為n元素的地球化學背景值,nR為n元素的參比值,nE為n元素的潛在生態(tài)風險,RI為總體潛在生態(tài)風險。為保證兩種評價法能夠獲得準確的土壤中重金屬潛在生態(tài)風險,本文取內(nèi)蒙古地區(qū) A層土壤中重金屬含量為背景值[14],標準化的重金屬毒性響應系數(shù)nT參照劉勇等[6]的研究取值。

    2.4遙感影像解譯及數(shù)據(jù)處理

    選用內(nèi)蒙古自治區(qū)烏拉特后旗 2014年 10月 Landsat8 TM 影像(http://www.gscloud.cn/),運用遙感圖像處理軟件 ERDAS IMAGINE 8.6 進行監(jiān)督分類,參考中國常用土地分類標準和實際情況,將研究區(qū)土地分為農(nóng)田、城鎮(zhèn)及工礦用地、沙地、林地、水域和未利用地 6 類,參考內(nèi)蒙古烏拉特后旗2014年度 1︰100 000 土地利用現(xiàn)狀圖,對文件中的判讀錯誤區(qū)域進行糾正,得到研究區(qū)2014年土地利用圖,用來分析判別土壤重金屬來源。常規(guī)數(shù)理統(tǒng)計分析采用 SPSS20.0軟件完成, 地統(tǒng)計分析、地累積指數(shù)及潛在生態(tài)風險分布制圖采用 GS+Version7軟件。

    3 結(jié)果與分析

    3.1重金屬的描述性統(tǒng)計結(jié)果

    研究區(qū)土壤重金屬含量的描述性統(tǒng)計結(jié)果見表2,從表 2 中可以看出,Zn、Pb、Cr、Cu、As、Hg 及 Cd 含量的平均值分別為 226.80、39.18、0.08、33.04、2.58、0.03、89.44 mg/kg,與我國土壤環(huán)境質(zhì)量一級標準(自然背景值)相比,Zn 和 Pb 已明顯構(gòu)成污染,均值分別超過標準值的 2.27 與 1.12 倍,并且于排放源附近樣點檢出的最大值超標 7.1 倍與6.59 倍,與內(nèi)蒙古土壤環(huán)境背景值相比,高出近 4倍與 2.3 倍。此外,Cd 與 Cu 含量的平均值已經(jīng)接近標準值,Cd、Hg、Cu 均發(fā)現(xiàn)超標樣點,存在污染風險,As 與 Hg 含量較低,檢測平均值相當于國家級標準的 0.17 倍和 0.20 倍,相比內(nèi)蒙古城市土壤環(huán)境背景也較低。

    表2 重金屬含量的描述性統(tǒng)計結(jié)果Table 2 Descriptive statistics of the heavy metal contents in soils of research area

    3.2重金屬的空間分異特征

    研究區(qū)土壤重金屬半變異函數(shù)擬合見圖 2,半變異函數(shù)的主要特征參數(shù)見表 3。由表 3 與圖 2 可知,Pb、Cr、Hg 和 Cd 的半變異函數(shù)最優(yōu)擬合模型為高斯模型(Gaussian),Cu 則為球狀模型(Spherical),As 為指數(shù)模型(Exponential),各種重金屬元素半變異函數(shù)模型的擬合度均在 0.75 以上,模型的選取基本符合要求[15-16]。

    塊金值反映的是最小抽樣尺度以下變量的變異性及測量誤差,通常認為由實驗誤差和小于實驗取樣尺度引起的變異,基臺值是半方差值隨步長增加到一個相對穩(wěn)定的水平上時對應的值,是區(qū)域變量總體特征的體現(xiàn),數(shù)值上等于空間結(jié)構(gòu)值和塊金值之和。塊金值與基臺值的比值為空間相關度,表明系統(tǒng)變量的空間相關性的程度[17],變程是指變異函數(shù)在有限步長上達到基臺值時對應的值。由表 3 可以看出,研究區(qū)內(nèi)重金屬元素一致表現(xiàn)為基臺值>塊金值,說明研究區(qū)重金屬空間變異的隨機性較大,這種隨機性主要源于低于采樣間距的小尺度下重金屬含量的空間變化。Zn、Pb 與 As 的塊金值與基臺值比值在 0.25 ~ 0.75 ,表明這 3 種重金屬含量具有中等的空間相關性,其中 As 與 Zn 的變程為 98.5 km 與 97.6 km,已接近變程最大理論值(100 km),表明 As 與 Zn 元素在大尺度上呈均一化分布形式。Cd、Cr、Cu 與 Hg的塊金值/基臺值>0.75,表明相關性很弱,其中,Cu 與 Cd 塊金值/基臺值均大于 0.9,基本不具有空間相關性,這與研究區(qū)重金屬排放源的地理位置、排放量不均勻時空分布有關。綜上,基本可以推斷研究區(qū)內(nèi)土壤中重金屬含量已經(jīng)很大程度上受到了人類活動的影響,結(jié)構(gòu)性因素如氣候、母質(zhì)、地形、土壤類型等因素的影響已經(jīng)不顯著。

    表3 重金屬含量的半變異函數(shù)及主要特征參數(shù)Table 3 Theoretical models and relevant parameters of soil heavy metal contents

    圖2 研究區(qū)土壤重金屬的半變異函數(shù)圖Fig. 2 Semi-variogram of the heavy metals and spherical model fitness

    3.3重金屬污染的綜合評價

    3.3.1地累積指數(shù)評價及來源分析由地積累指數(shù)(Igeo)計算結(jié)果可知(表 4),除 Cr 外,研究區(qū)土壤均受到其他重金屬元素的污染,受到中度以上 Zn污染的土壤樣本占比達 67% 以上,其中 5% 污染程度達到強度;84% 的土壤樣本受到 Pb 污染,其中 70% 的樣本處于無-中度污染級別,2% 受到強度污染;54%、44% 的樣本受到 Cu 及 As 中度以下程度污染,47% 樣本受到 Hg 污染,均為中度-強度程度以下,100% 的樣本受到了 Cd 污染,程度均為極強。各采樣點重金屬元素的 Igeo均值排列順序為Cd(9.07)>Zn(1.67)>Pb(0.52)>Cu(0.12)>As(-0.06)>Hg(-0.24)>Cr(-10.19)。相比其他研究結(jié)果[18],本研究區(qū)污染樣本占比率較高,強度相對較高,土壤環(huán)境整體上受到污染。表 3 顯示 Zn、As 的變程均大于研究區(qū)域長度(45 km),其原因有待進一步分析。從其余元素的評價結(jié)果來看,Cr 基本未受到污染,濃度水平維持在本底值附近。研究區(qū)所有采樣點的 Cd均處于極強的污染程度,且累積強度較大。有研究表明,Cd 在土壤中可隨水滲漏遷移到 40 cm 以下的犁底層,并且容易被土豆、小麥等作物吸收進入食物鏈[19],因此,有必要開展針對 Cd 的污染控制與土壤修復,保護當?shù)剞r(nóng)產(chǎn)品食品安全。

    表4 研究區(qū)土壤重金屬污染地積累指數(shù)分級評價Table 4 The classification of heavy metal pollution based on the index of geo-accumulation in the soil of research area

    3.3.2潛在生態(tài)風險評價表 5 為潛在生態(tài)風險指數(shù)分級,研究區(qū)域土壤 7 種重金屬的平均潛在生態(tài)風險指數(shù)和潛在生態(tài)風險綜合指數(shù)樣本占比如表 6 所示,綜合表 5 與表 6 可知,Zn、Cr、Cu、As、Hg 元素在所有樣點均處于輕微潛在生態(tài)風險程度內(nèi),2% 的樣本受 Pb 影響有中等潛在生態(tài)風險,55%的樣本受 Cd 影響有中等潛在生態(tài)風險,其中,12%的樣點達到了強潛在生態(tài)風險級別。從各元素 En 平均值來看來看,Cd 元素處于中等危害,其他元素均處于輕微危害,7 種重金屬的潛在生態(tài)風險由強至弱依次為:Cd>Cu>Pb>As>Zn>Hg>Cr,其中,Cd 的貢獻最突出。

    表5 生態(tài)危害系數(shù)、指數(shù)與危害程度分級Table 5 Ecological risk coefficient, risk index and classification of risk intensity

    表6 研究區(qū)土壤重金屬生態(tài)危害評價指數(shù)Table 6 Risk index of heavy metals in the soil of research area

    對照研究區(qū)土地利用與潛在生態(tài)風險指數(shù)(RI)空間分布來看(圖3),研究區(qū)土壤中重金屬的RI在西部城鎮(zhèn)用地附近的未利用地處、中部的城鎮(zhèn)工礦用地處及其東北側(cè)的農(nóng)田區(qū)域偏高,在其他區(qū)域內(nèi)普遍偏低。在實地調(diào)查中發(fā)現(xiàn),研究區(qū)西部未利用地原為某大型采掘企業(yè)礦石堆場,再加上研究區(qū)地勢西高東低,中部城鎮(zhèn)區(qū)生活污水常年隨地表徑流進入東北側(cè)的農(nóng)田區(qū)域中,導致該區(qū)域生態(tài)風險增加,此外,農(nóng)田區(qū)域較其他未利用地生態(tài)風險較高可能與農(nóng)業(yè)施肥有關。對比采樣分布圖中有色金屬采掘企業(yè)分布來看,RI基本上呈現(xiàn)以排放源為中心擴散式分布,表明其分布基本上服從于冶煉企業(yè)的分布,揭示了高環(huán)境污染區(qū)將伴隨高的生態(tài)風險,同時也說明人口密集、工業(yè)化程度較高的城鎮(zhèn)區(qū)域更容易造成生態(tài)危害,這將為土壤整治采取有效措施、控制污染源頭提供技術(shù)支持。綜合來看,研究區(qū)域土壤重金屬的潛在生態(tài)風險可能由礦山采冶、生活源排放和耕地區(qū)域內(nèi)的農(nóng)業(yè)施肥等導致[18]。基于各元素平均含量計算的研究區(qū)土壤的RI值為85.61,表明研究區(qū)土壤總體處于輕微風險程度。

    3.3.3討論對比可知,研究區(qū)內(nèi)土壤中重金屬的地積累指數(shù)和潛在生態(tài)風險指數(shù)存在差異,主要區(qū)別在Cu、As 和 Zn,這主要由重金屬毒性系數(shù)的影響導致。As 的毒性系數(shù)最高(Tn= 10),但濃度上與內(nèi)蒙古地區(qū)背景值比較相差不大,導致其在兩種評價方法中的排位基本不變,Cu 的毒性系數(shù)(Tn=5)>Zn(Tn=1),因此,Cu 元素的潛在生態(tài)風險排位超過 Zn。

    地積累指數(shù)法側(cè)重于重金屬含量與背景值的對比評價,主要表征外源重金屬的富集程度,利用地累積指數(shù)可以很好地推斷外源的位置與貢獻量,本研究的地累積指數(shù)與空間變異特征結(jié)果分析基本一致,在一定程度上表明研究區(qū)土壤中 Cd 含量明顯受到外源的影響。Zn、Pb 及 As 呈現(xiàn)中度空間相關性,地累積指數(shù)卻顯示土壤受到污染的現(xiàn)象,這是由本底值與實測值差距較大及受到外源影響共同造成的,兩者結(jié)果并不矛盾。潛在生態(tài)風險指數(shù)評價方法相比地累積指數(shù)法增加考慮了重金屬的生物毒性影響,更注重于重金屬對生態(tài)系統(tǒng)的危害。對本研究區(qū)來說,兩種方法計算的結(jié)果均表明研究區(qū)土壤受到污染,并且具有一定程度的潛在生態(tài)風險。

    3.4重金屬可能的來源分析

    結(jié)合研究區(qū)土壤重金屬地累積指數(shù)空間分布圖(圖 4)分析可知,As、Cu、Hg 元素基本維持在無-中度以下,研究區(qū) Cr 元素地累積指數(shù)均小于 0,表明研究區(qū)土壤基本未受到 Cr 的累積污染。Pb 高值區(qū)均主要分布在城鎮(zhèn)及工礦用地,Zn 高值區(qū)主要分布在中部城鎮(zhèn)區(qū)及周邊農(nóng)田。Cd 整體指數(shù)較高,且分布較凌亂,但大致呈現(xiàn)從中心向兩側(cè)遞減的趨勢,主要高值區(qū)集中在中部農(nóng)田及中部城鎮(zhèn)區(qū)東側(cè)。

    Pb、Zn 及 Cd 常被作為交通污染源的標識元素[19],汽車輪胎與地面摩擦可產(chǎn)生含 Pb、Zn、Cd 較高的顆粒物[20],使用的含 Pb 汽油及拉運的礦石細粉在風的擴散作用下,在某種程度上也可以加重了上述元素的富集。Zn 和 Cd 的高值區(qū)均有出現(xiàn)在耕地區(qū)域內(nèi)的現(xiàn)象,可能表明農(nóng)業(yè)源是其主要污染源,含磷復合肥及以城市垃圾、污泥為原料的肥料中均混雜有Zn 和 Cd[21]。在實際調(diào)查中發(fā)現(xiàn),當?shù)剞r(nóng)業(yè)耕作中大量使用含磷化肥,一般過磷酸鹽中含有較多的 Cd 與 Zn,連續(xù)施用磷肥可引起表土中這些元素含量的顯著增加,從而導致土壤中有累積現(xiàn)象發(fā)生。此外,研究區(qū)域地勢西高東低,常年主導風向為西北偏西,研究區(qū)中部的城鎮(zhèn)尚未解決集中供暖,因此,推斷研究區(qū) Pb、Zn、Cd 的累積也可能與冬季采暖小鍋爐排放的煙塵有關。對比土樣中重金屬元素半變異函數(shù)的變程來看,As 與 Zn 的變程均接近理論最大值,這是由于該元素在土壤中不易遷移所致。Cu 和 Pb地累積指數(shù)的變化基本呈現(xiàn)以采冶企業(yè)為中心遞減的規(guī)律,表明有色金屬采冶是其主要的影響因素。Cd 高值區(qū)域主要集中在城鎮(zhèn)工礦用地地表水流向的下游方向,這與城鎮(zhèn)生活污水散排易于空間擴散的特征相符[22]。綜合以上分析可知,研究區(qū)土壤重金屬的來源受多種因素影響,包括城鎮(zhèn)工礦區(qū)域和耕地區(qū)域內(nèi)的日常生活和生產(chǎn)活動,比如有色金屬采冶、交通運輸、燃煤釋放、農(nóng)業(yè)施肥和污水灌溉等。

    4 結(jié)論

    1) 研究區(qū) Cd、Cu 與 Hg 空間相關性很弱,Cr 與 Cd 基本不具有空間相關性,推斷研究區(qū)內(nèi)重金屬含量已經(jīng)很大程度上受到了人類活動的影響,工農(nóng)業(yè)生產(chǎn),特別是工業(yè)排放的時空分布格局使土壤重金屬的空間相關性減弱,結(jié)構(gòu)性因素如氣候、母質(zhì)、地形、土壤類型等因素的影響已經(jīng)不顯著。

    2) 研究區(qū)各采樣點重金屬元素的地累積指數(shù)均值排列順序為 Cd>Zn>Pb>Cu>As>Hg>Cr,除 Cr外,各重金屬元素在研究區(qū)內(nèi)有不同程度累積,其中,100% 的土壤樣本表現(xiàn)出極強的 Cd 污染程度。

    3) 研究區(qū)Cd 元素處于中等危害,其他元素均處于輕微危害,7 種重金屬元素的潛在生態(tài)風險由強至弱依次為:Cd>Cu>Pb>As>Zn>Hg>Cr,其中,Cd 的貢獻最突出,研究區(qū)土壤總體處于輕微生態(tài)風險程度。4) 研究區(qū)重金屬含量的累積主要由區(qū)域內(nèi)的有色金屬冶煉、城鎮(zhèn)交通運輸、燃煤、耕地農(nóng)業(yè)施肥及污水灌溉等因素共同導致。

    圖3 研究區(qū)域土地利用及潛在生態(tài)危險指數(shù)分布Fig. 3 The distribution of potential ecological risk index as related with land use types in research area

    圖4 研究區(qū)域土壤重金屬污染地累積指數(shù)空間分布Fig. 4 The spatial variation of heavy metals based on the index of geoaccumulation in research area

    參考文獻:

    [1]Hu X F, Wu H X, Hu X, et al. Impact of urbanization on Shanghai's soil environmental quality[J]. Pedosphere, 2004,14(2): 151-158

    [2] Li X P, Huang C C. Environment impact of heavy metals on urban soil in the vicinity of industrial area of Baoji city, P.R.China[J]. Environmental Geology, 2007,52(8): 1 631-1 637

    [3] 李曉燕, 陳同斌, 雷梅, 等. 不同土地利用方式下北京城區(qū)土壤的重金屬累積特征[J] . 環(huán)境科學學報, 2010,30(11): 2 285-2 293

    [4] 汪權(quán)方, 陳百明, 李家永, 等. 城市土壤研究進展與中國城市土壤生態(tài)保護研究[J]. 水土保持學報, 2003, 17(4):142-145

    [5] 吳新民, 李戀卿, 潘根興, 等. 南京市不同功能城區(qū)土壤中重金屬Cu、Zn、Pb 和Cd 的污染特征[J]. 環(huán)境科學, 2003, 24(3): 105-111

    [6] 劉勇, 岳玲玲, 李晉昌. 太原市土壤重金屬污染及其潛在生態(tài)風險評價[J]. 環(huán)境科學學報, 2011, 31(6): 1 285-1 293

    [7] 柳云龍, 章立佳, 韓曉非, 等. 上海城市樣帶土壤重金屬空間變異特征及污染評價[J]. 環(huán)境科學, 2012, 22(2):599-605

    [8] 師榮光, 趙玉杰, 周啟星, 等. 蘇北優(yōu)勢農(nóng)業(yè)區(qū)土壤砷含量空間變異性研究[J]. 農(nóng)業(yè)工程學報, 2008, 24(1):80-84

    [9]周靜, 崔紅標, 梁家妮, 等. 重金屬污染土壤修復技術(shù)的選擇和面臨的問題-以江銅貴冶九牛崗土壤修復示范工程項目為例[J]. 土壤, 2015, 47(2): 283-288

    [10] 曹雪瑩, 張莎娜, 譚長銀, 等. 中南大型有色金屬冶煉廠周邊農(nóng)田土壤重金屬污染特征研究[J]. 土壤, 2015,47(1): 94-99

    [11] Isaaks E H, Srivastava R M. An introduction to applied geo statistics[M]. New York:Oxford University Press, 1989

    [12] 許振成, 楊曉云, 溫勇, 等. 北江中上游底泥重金屬污染及其潛在生態(tài)危害評價[J]. 環(huán)境科學, 2009, 30(11):3 262-3 268

    [13] Hakanson L. An ecological risk index for aquatic pollution control-A sediment logical approach[J]. Water Research,1980, 14(8): 975-1 001

    [14] 高紅霞, 王喜寬, 張青, 等. 內(nèi)蒙古河套地區(qū)土壤背景值特征[J]. 地質(zhì)與資源, 2007, 16(3): 209-212

    [15] 侯景儒, 黃競先, 地質(zhì)統(tǒng)計學的理論與方法[M]. 北京:地質(zhì)出版社, 1990

    [16] 候景儒, 尹鎮(zhèn)南, 李維明. 實用地質(zhì)統(tǒng)計學[M]. 北京:地質(zhì)出版社, 1998

    [17] 鐘曉蘭, 周生路, 李江濤, 等. 長江三角洲地區(qū)土壤重金屬污染的空間變異特征——以江蘇省太倉市為例[J].土壤學報, 2007, 44(1): 33-40

    [18] 陳玉東, 王火焰, 周健民, 等. 黑龍江海倫市農(nóng)田土壤重金屬與磷素含量的特征研究[J]. 土壤,2012, 44(4):613-620

    [19] 肖鵬飛, 李法云, 付寶榮, 等. 土壤重金屬污染及其植物修復研究[J]. 遼寧大學學報( 自然科學版), 2004,31(3): 279-283

    [20] Friedlander S. Chemical element balances and identification of air pollution sources[J]. Environmental Science &Technology, 1973. 7(3): 235-240

    [21] 王軍, 陳振樓, 王初, 等. 上海崇明島蔬菜地土壤重金屬含量與生態(tài)風險預警評估[J]. 環(huán)境科學, 2007, 28(3):647-653

    [22] 吳紹華, 周生路, 楊得志, 等. 宜興市近郊土壤重金屬來源與空間分布研究[J]. 科學通報, 2008, 53(增刊):162-170

    Evaluation on Heavy Metal Contamination and Its Potential Ecological Risk in Soil: A Case Study of Nonferrous Metal Smelting Zone in Wulatehouqi

    YUE Rong1, SHI Rui1*, ZHANG Hong2
    (1 Institute of Environmental Science of Bayannur, Inner Mongolia Autonomous Region, Linhe, Inner Mongolia015000, China;2 College of Environmental Science and Resources, Shanxi University, Taiyuan030006, China)

    Abstract:The concentrations of 7 heavy metals in soil in the nonferrous smelting area in Wulatehouqi, Inner Mongolia Autonomous Region were determined. The spatial variability, potential sources, and ecological risk of heavy metals were analyzed using geostatisticals, geoaccumulation index, and potential ecological risk index. The results indicated that: 1) the concentration and distribution of 7 heavy metals in soil of this area were mainly influenced by human activities, and the effects of soil forming factors were not significant; 2) the geoaccumulation index for these metals followed the order: Cd>Zn>Pb>Cu>As>Hg>Cr, and the potential ecological risk index followed the order: Cd>Cu>Pb>As>Zn>Hg>Cr. All soil samples were significantly polluted by Cd; 3) the comprehensive ecological risk was at a relatively low level. The potential pollution sources were nonferrous metal mining and smelting, transportation, coal combustion, and irrigation with sewage.

    Key words:Soil heavy metals; Spatial variability; Ecological risk assessment; Potential pollution sources

    中圖分類號:X 82

    DOI:10.13758/j.cnki.tr.2016.02.016

    基金項目:①巴彥淖爾市重金屬專項研究項目資助。

    * 通訊作者(84650100@qq.com)

    作者簡介:岳榮(1968—),男,內(nèi)蒙古巴彥淖爾市人,高級工程師,主要從事環(huán)境評估與生態(tài)恢復研究。E-mail: yr2652@163.com

    猜你喜歡
    土壤重金屬成因
    說說兩種『白氣』的成因
    色散的成因和應用
    談學困生的成因及轉(zhuǎn)化策略
    甘肅教育(2020年12期)2020-04-13 06:24:44
    暈紋石成因解讀(上)
    寶藏(2017年7期)2017-08-09 08:15:19
    有機肥對土壤重金屬生物有效性影響研究進展
    基于Kriging插值的礦區(qū)周邊土壤重金屬空間分布規(guī)律研究
    科技視界(2016年13期)2016-06-13 08:06:36
    翻譯實踐問題及成因
    唐山文學(2016年11期)2016-03-20 15:25:54
    貴州擺紀磷石膏尾礦庫周邊土壤重金屬污染評價
    某地區(qū)表層土壤重金屬分布特征研究
    土壤重金屬的污染特點及治理策略研究
    科學與財富(2014年9期)2014-09-12 06:41:53
    tube8黄色片| 日韩制服丝袜自拍偷拍| 男女边吃奶边做爰视频| 亚洲色图 男人天堂 中文字幕| 精品免费久久久久久久清纯 | 亚洲国产av影院在线观看| 母亲3免费完整高清在线观看| 亚洲久久久国产精品| 国产野战对白在线观看| 久久人人爽人人片av| 国产日韩欧美视频二区| 一级黄片播放器| 少妇 在线观看| 久久久久精品性色| 一边摸一边抽搐一进一出视频| 精品少妇内射三级| 成人18禁高潮啪啪吃奶动态图| 大片电影免费在线观看免费| 色婷婷av一区二区三区视频| 婷婷色综合大香蕉| 久久精品亚洲av国产电影网| 亚洲精品自拍成人| 精品福利永久在线观看| 欧美日韩综合久久久久久| 深夜精品福利| 性高湖久久久久久久久免费观看| 欧美日韩亚洲综合一区二区三区_| 国产av码专区亚洲av| 一本一本久久a久久精品综合妖精| 精品人妻在线不人妻| 日韩电影二区| 国产淫语在线视频| 日韩视频在线欧美| 国产一区有黄有色的免费视频| 免费黄网站久久成人精品| 妹子高潮喷水视频| 国产乱来视频区| 国产精品久久久久成人av| 亚洲精品国产av成人精品| 丁香六月欧美| 亚洲 欧美一区二区三区| 亚洲 欧美一区二区三区| 久久久久久久精品精品| 99久国产av精品国产电影| 在线免费观看不下载黄p国产| 九色亚洲精品在线播放| 老司机在亚洲福利影院| 99久国产av精品国产电影| 人妻人人澡人人爽人人| 一本久久精品| 韩国高清视频一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 一本一本久久a久久精品综合妖精| 香蕉丝袜av| 超色免费av| 日本wwww免费看| 99久久精品国产亚洲精品| 久久久精品国产亚洲av高清涩受| 久久精品久久久久久噜噜老黄| 久久精品久久久久久噜噜老黄| 精品一区二区免费观看| 日韩不卡一区二区三区视频在线| 精品国产一区二区久久| 国产 一区精品| 香蕉国产在线看| 如日韩欧美国产精品一区二区三区| 男的添女的下面高潮视频| 久久久亚洲精品成人影院| 欧美在线一区亚洲| 亚洲精品aⅴ在线观看| 国产极品粉嫩免费观看在线| 一级毛片电影观看| 国产视频首页在线观看| 国产亚洲最大av| 国产精品女同一区二区软件| 青春草亚洲视频在线观看| 成人18禁高潮啪啪吃奶动态图| 久久久久人妻精品一区果冻| 综合色丁香网| xxxhd国产人妻xxx| 亚洲精品av麻豆狂野| 卡戴珊不雅视频在线播放| 9191精品国产免费久久| 午夜福利视频精品| 丝袜美足系列| 欧美 亚洲 国产 日韩一| 国产成人欧美| 国产午夜精品一二区理论片| 国产免费现黄频在线看| 波野结衣二区三区在线| 成人免费观看视频高清| 99精国产麻豆久久婷婷| 一级毛片黄色毛片免费观看视频| 满18在线观看网站| 超色免费av| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级专区第一集| 街头女战士在线观看网站| 国产黄频视频在线观看| 精品国产一区二区久久| av又黄又爽大尺度在线免费看| 日韩大片免费观看网站| 精品久久蜜臀av无| 在线观看免费高清a一片| 亚洲精品中文字幕在线视频| 18禁观看日本| 十八禁高潮呻吟视频| 巨乳人妻的诱惑在线观看| 欧美在线黄色| 国产一卡二卡三卡精品 | 日韩人妻精品一区2区三区| 街头女战士在线观看网站| 精品国产乱码久久久久久小说| 国产成人免费无遮挡视频| 中文字幕高清在线视频| 国产片特级美女逼逼视频| 晚上一个人看的免费电影| 韩国精品一区二区三区| 少妇精品久久久久久久| 久久国产精品男人的天堂亚洲| 宅男免费午夜| 国产在线免费精品| 天堂中文最新版在线下载| 美女大奶头黄色视频| 亚洲精品日韩在线中文字幕| 成人黄色视频免费在线看| e午夜精品久久久久久久| www.熟女人妻精品国产| 中国国产av一级| av在线播放精品| av电影中文网址| 另类精品久久| av视频免费观看在线观看| 亚洲国产日韩一区二区| 午夜福利在线免费观看网站| 天堂中文最新版在线下载| 久久天躁狠狠躁夜夜2o2o | 中文字幕另类日韩欧美亚洲嫩草| 免费观看av网站的网址| 婷婷色麻豆天堂久久| 国产极品粉嫩免费观看在线| 亚洲精品国产一区二区精华液| 国产亚洲最大av| 亚洲,一卡二卡三卡| 女性被躁到高潮视频| 久久久久精品性色| 日韩中文字幕欧美一区二区 | 人妻 亚洲 视频| 国产av码专区亚洲av| 久久精品亚洲av国产电影网| 狠狠精品人妻久久久久久综合| 欧美日韩av久久| 九色亚洲精品在线播放| 国产精品久久久久久久久免| 丁香六月天网| 国产成人一区二区在线| 午夜福利免费观看在线| 黄片无遮挡物在线观看| 国产精品 欧美亚洲| 一本大道久久a久久精品| 丝袜脚勾引网站| 18禁动态无遮挡网站| 国产成人系列免费观看| 午夜免费男女啪啪视频观看| 久久久国产精品麻豆| 麻豆av在线久日| 精品久久久久久电影网| 伦理电影免费视频| 欧美人与性动交α欧美软件| 日韩中文字幕欧美一区二区 | 亚洲国产看品久久| 成人亚洲欧美一区二区av| 91aial.com中文字幕在线观看| 一级毛片 在线播放| 青春草亚洲视频在线观看| 老司机影院成人| 在线观看www视频免费| 少妇精品久久久久久久| 岛国毛片在线播放| 亚洲美女黄色视频免费看| 亚洲一级一片aⅴ在线观看| 国产一区二区 视频在线| 午夜91福利影院| 欧美精品亚洲一区二区| 亚洲美女搞黄在线观看| 欧美激情极品国产一区二区三区| av在线观看视频网站免费| 伦理电影免费视频| a级毛片黄视频| 亚洲美女搞黄在线观看| kizo精华| 精品免费久久久久久久清纯 | 亚洲av日韩精品久久久久久密 | 在现免费观看毛片| 亚洲国产精品国产精品| 欧美激情高清一区二区三区 | 国产一级毛片在线| 国产精品香港三级国产av潘金莲 | 日韩熟女老妇一区二区性免费视频| 免费观看av网站的网址| 人成视频在线观看免费观看| 欧美人与性动交α欧美精品济南到| 亚洲综合色网址| 美女福利国产在线| 少妇人妻精品综合一区二区| 午夜福利视频在线观看免费| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区91 | 777久久人妻少妇嫩草av网站| 久久精品aⅴ一区二区三区四区| 黄色视频不卡| 国产成人精品无人区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产一区二区| 久久97久久精品| 九色亚洲精品在线播放| 91成人精品电影| 亚洲成av片中文字幕在线观看| 亚洲av欧美aⅴ国产| 精品国产一区二区久久| 国产精品女同一区二区软件| 我要看黄色一级片免费的| 妹子高潮喷水视频| 欧美黑人欧美精品刺激| 男人操女人黄网站| 亚洲情色 制服丝袜| 欧美日韩av久久| 亚洲国产精品国产精品| 欧美激情极品国产一区二区三区| 在线看a的网站| 国产精品秋霞免费鲁丝片| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 精品亚洲成a人片在线观看| 啦啦啦 在线观看视频| 久久久久久久久久久久大奶| 国产欧美日韩综合在线一区二区| 亚洲av电影在线观看一区二区三区| 曰老女人黄片| 80岁老熟妇乱子伦牲交| 亚洲国产欧美一区二区综合| 少妇人妻久久综合中文| 看免费av毛片| 国产亚洲av高清不卡| 国产日韩欧美在线精品| 成人亚洲精品一区在线观看| 在线观看免费日韩欧美大片| 日韩av免费高清视频| 日本欧美国产在线视频| 久久人人97超碰香蕉20202| 80岁老熟妇乱子伦牲交| 久久99一区二区三区| 男人操女人黄网站| 伦理电影免费视频| 欧美人与性动交α欧美软件| 九九爱精品视频在线观看| 在线天堂中文资源库| 婷婷色麻豆天堂久久| 美女扒开内裤让男人捅视频| 90打野战视频偷拍视频| 99精品久久久久人妻精品| 国产成人av激情在线播放| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 久久精品人人爽人人爽视色| 制服诱惑二区| 国产日韩欧美亚洲二区| h视频一区二区三区| 亚洲七黄色美女视频| 日日撸夜夜添| 国产99久久九九免费精品| 亚洲国产成人一精品久久久| 国产精品一区二区在线观看99| 桃花免费在线播放| 国产成人a∨麻豆精品| 天天躁夜夜躁狠狠躁躁| 国产精品国产av在线观看| 大话2 男鬼变身卡| 桃花免费在线播放| 丁香六月欧美| 晚上一个人看的免费电影| 免费在线观看完整版高清| 男的添女的下面高潮视频| xxx大片免费视频| netflix在线观看网站| 久久精品国产a三级三级三级| 欧美黑人欧美精品刺激| 一区在线观看完整版| 啦啦啦在线免费观看视频4| www.自偷自拍.com| 波野结衣二区三区在线| 亚洲伊人久久精品综合| 日日爽夜夜爽网站| 丝袜在线中文字幕| 久久精品国产亚洲av涩爱| 汤姆久久久久久久影院中文字幕| 国产亚洲午夜精品一区二区久久| 高清欧美精品videossex| 91成人精品电影| 欧美日韩亚洲综合一区二区三区_| 一区二区av电影网| 一区二区三区四区激情视频| 久久久精品区二区三区| 午夜福利,免费看| 亚洲国产av影院在线观看| 欧美国产精品va在线观看不卡| 亚洲第一青青草原| 丝袜人妻中文字幕| 久久久久精品国产欧美久久久 | 日韩欧美一区视频在线观看| 在线看a的网站| 大码成人一级视频| 男人爽女人下面视频在线观看| 国产一区有黄有色的免费视频| 黄频高清免费视频| 最新在线观看一区二区三区 | 免费人妻精品一区二区三区视频| 亚洲欧洲精品一区二区精品久久久 | 中文字幕亚洲精品专区| 宅男免费午夜| 91国产中文字幕| 欧美日韩视频高清一区二区三区二| 侵犯人妻中文字幕一二三四区| 午夜福利视频在线观看免费| 老司机影院毛片| 国产成人免费观看mmmm| 亚洲国产最新在线播放| 亚洲伊人色综图| 美女高潮到喷水免费观看| 亚洲色图 男人天堂 中文字幕| 久久 成人 亚洲| 国产精品香港三级国产av潘金莲 | 久久人妻熟女aⅴ| 国产亚洲午夜精品一区二区久久| 欧美精品亚洲一区二区| 国产成人一区二区在线| 麻豆av在线久日| 国产视频首页在线观看| 99热全是精品| 精品人妻熟女毛片av久久网站| 如日韩欧美国产精品一区二区三区| 久久99热这里只频精品6学生| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 免费久久久久久久精品成人欧美视频| 亚洲三区欧美一区| 国产欧美亚洲国产| 日本一区二区免费在线视频| 男男h啪啪无遮挡| 亚洲av在线观看美女高潮| 欧美日韩一区二区视频在线观看视频在线| 成人国产av品久久久| 色吧在线观看| 在线观看www视频免费| 女人久久www免费人成看片| 性少妇av在线| 五月天丁香电影| 欧美老熟妇乱子伦牲交| 两性夫妻黄色片| 高清不卡的av网站| 亚洲美女搞黄在线观看| 免费观看人在逋| 一边亲一边摸免费视频| 久热爱精品视频在线9| 亚洲,欧美,日韩| 制服诱惑二区| 精品少妇内射三级| 91成人精品电影| 亚洲精品视频女| 又黄又粗又硬又大视频| 少妇被粗大的猛进出69影院| 日韩大码丰满熟妇| 天美传媒精品一区二区| 韩国精品一区二区三区| 国产免费福利视频在线观看| 熟女少妇亚洲综合色aaa.| 欧美日韩精品网址| av国产久精品久网站免费入址| 欧美在线黄色| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人a∨麻豆精品| 搡老岳熟女国产| 1024视频免费在线观看| 搡老乐熟女国产| 国产精品久久久久久久久免| 中文字幕人妻丝袜一区二区 | 久久久久国产精品人妻一区二区| 亚洲成人国产一区在线观看 | 别揉我奶头~嗯~啊~动态视频 | 电影成人av| 欧美久久黑人一区二区| 中文字幕制服av| 久久国产精品大桥未久av| 一级片'在线观看视频| 水蜜桃什么品种好| 精品人妻熟女毛片av久久网站| 日韩一卡2卡3卡4卡2021年| 国产老妇伦熟女老妇高清| avwww免费| av免费观看日本| 国产老妇伦熟女老妇高清| 宅男免费午夜| 久久久国产一区二区| 在线观看www视频免费| 色94色欧美一区二区| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 欧美精品高潮呻吟av久久| 精品国产乱码久久久久久男人| 日韩av免费高清视频| 99热全是精品| 男女床上黄色一级片免费看| 街头女战士在线观看网站| 免费日韩欧美在线观看| 国产在线视频一区二区| 亚洲av男天堂| 9热在线视频观看99| 高清视频免费观看一区二区| 精品一区二区三区四区五区乱码 | 中国三级夫妇交换| 最近的中文字幕免费完整| 人人妻人人爽人人添夜夜欢视频| 国产成人91sexporn| 超碰97精品在线观看| 在线天堂最新版资源| 成人国语在线视频| 热re99久久精品国产66热6| 看十八女毛片水多多多| 国产亚洲欧美精品永久| 国产亚洲一区二区精品| 免费高清在线观看日韩| 久久午夜综合久久蜜桃| 精品一区二区三卡| 亚洲国产精品国产精品| 久久人妻熟女aⅴ| 国产av国产精品国产| 男女午夜视频在线观看| 日韩人妻精品一区2区三区| 校园人妻丝袜中文字幕| 91成人精品电影| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 日韩视频在线欧美| 欧美人与性动交α欧美精品济南到| a级毛片黄视频| 久久精品人人爽人人爽视色| 亚洲欧洲日产国产| 纯流量卡能插随身wifi吗| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 午夜福利乱码中文字幕| 女人精品久久久久毛片| 热re99久久精品国产66热6| 国产不卡av网站在线观看| 免费在线观看黄色视频的| 最近2019中文字幕mv第一页| 亚洲av电影在线进入| 成人漫画全彩无遮挡| 国产色婷婷99| 日韩伦理黄色片| 女性被躁到高潮视频| 女性生殖器流出的白浆| 丰满乱子伦码专区| 99久久人妻综合| 超色免费av| 色婷婷久久久亚洲欧美| 亚洲国产精品一区二区三区在线| 日本爱情动作片www.在线观看| 久久精品熟女亚洲av麻豆精品| 欧美人与性动交α欧美精品济南到| 男人爽女人下面视频在线观看| 欧美激情极品国产一区二区三区| 国产精品三级大全| 亚洲精品一区蜜桃| 1024视频免费在线观看| 午夜福利影视在线免费观看| 老汉色av国产亚洲站长工具| e午夜精品久久久久久久| 亚洲精品在线美女| 久久免费观看电影| 我要看黄色一级片免费的| 亚洲久久久国产精品| 欧美日韩亚洲国产一区二区在线观看 | 免费黄网站久久成人精品| 超色免费av| 久久久久网色| 欧美最新免费一区二区三区| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 电影成人av| 欧美日韩精品网址| 中文字幕另类日韩欧美亚洲嫩草| 观看美女的网站| 久久久精品免费免费高清| 最新在线观看一区二区三区 | 国产成人av激情在线播放| 男人添女人高潮全过程视频| 日韩制服骚丝袜av| 久久精品国产综合久久久| 国产精品久久久人人做人人爽| 国产亚洲欧美精品永久| 午夜福利视频在线观看免费| 少妇被粗大猛烈的视频| 麻豆精品久久久久久蜜桃| 男女下面插进去视频免费观看| 国产免费一区二区三区四区乱码| 亚洲在久久综合| av天堂久久9| 久久婷婷青草| 精品久久蜜臀av无| 亚洲av在线观看美女高潮| 成人黄色视频免费在线看| 麻豆精品久久久久久蜜桃| 欧美精品一区二区免费开放| 国产亚洲av高清不卡| 久久久亚洲精品成人影院| 我要看黄色一级片免费的| 国产精品一区二区在线观看99| 青春草亚洲视频在线观看| 在线观看免费视频网站a站| 少妇被粗大猛烈的视频| 亚洲欧美激情在线| 午夜福利网站1000一区二区三区| netflix在线观看网站| 不卡视频在线观看欧美| 欧美精品av麻豆av| 久久久精品免费免费高清| 午夜久久久在线观看| 成人午夜精彩视频在线观看| 色综合欧美亚洲国产小说| 一级爰片在线观看| av片东京热男人的天堂| 老熟女久久久| 国产精品无大码| 国产乱来视频区| 国产女主播在线喷水免费视频网站| 欧美人与性动交α欧美软件| 久久女婷五月综合色啪小说| 亚洲 欧美一区二区三区| 男女床上黄色一级片免费看| 久久久精品94久久精品| 国产又爽黄色视频| 2018国产大陆天天弄谢| 国产xxxxx性猛交| 国产激情久久老熟女| 国产色婷婷99| 亚洲国产欧美在线一区| 搡老岳熟女国产| 亚洲一级一片aⅴ在线观看| 日韩制服丝袜自拍偷拍| 操出白浆在线播放| 国产一级毛片在线| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 亚洲av欧美aⅴ国产| 曰老女人黄片| 一级片免费观看大全| a级毛片黄视频| 色精品久久人妻99蜜桃| 亚洲欧美激情在线| 欧美精品一区二区免费开放| 性色av一级| 亚洲精品,欧美精品| 亚洲五月色婷婷综合| 美女大奶头黄色视频| 中文字幕精品免费在线观看视频| 欧美日韩综合久久久久久| 丝袜美足系列| 亚洲av日韩精品久久久久久密 | 少妇被粗大猛烈的视频| 三上悠亚av全集在线观看| 少妇人妻久久综合中文| 18禁观看日本| 日本vs欧美在线观看视频| 日韩制服丝袜自拍偷拍| 亚洲欧美成人综合另类久久久| 国产欧美日韩一区二区三区在线| av福利片在线| av线在线观看网站| 亚洲欧洲国产日韩| 18禁观看日本| 国产成人91sexporn| 国产乱人偷精品视频| 国产精品亚洲av一区麻豆 | 99香蕉大伊视频| 中文字幕制服av| 又大又黄又爽视频免费| 欧美亚洲日本最大视频资源| 国产麻豆69| 青春草视频在线免费观看| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 人人澡人人妻人| 女人久久www免费人成看片| 在线观看免费高清a一片| 看十八女毛片水多多多| 精品国产国语对白av| 精品国产乱码久久久久久男人| 免费日韩欧美在线观看| 成人毛片60女人毛片免费| 国产精品一区二区精品视频观看| 色精品久久人妻99蜜桃| 最近的中文字幕免费完整| 国产黄色免费在线视频| 在线观看www视频免费| 日本午夜av视频| 日日啪夜夜爽| 18禁裸乳无遮挡动漫免费视频| 欧美在线黄色| 国产人伦9x9x在线观看| 九草在线视频观看| 五月天丁香电影| 两性夫妻黄色片| 国产深夜福利视频在线观看| 七月丁香在线播放|