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

    錫林郭勒露天煤礦礦區(qū)草原土壤重金屬分布特征

    2016-08-10 09:15:53楊勇劉愛軍朝魯孟其其格單玉梅烏尼圖陳海軍王明玖1內蒙古農業(yè)大學生態(tài)環(huán)境學院內蒙古呼和浩特010019內蒙古自治區(qū)草原勘察規(guī)劃院內蒙古呼和浩特010051內蒙古自治區(qū)農牧業(yè)科學院內蒙古呼和浩特010031中國科學院內蒙古草業(yè)研究中心內蒙古呼和浩特010031內蒙古自治區(qū)生物技術研究院內蒙古呼和浩特010010
    生態(tài)環(huán)境學報 2016年5期
    關鍵詞:土壤重金屬錫林郭勒分布特征

    楊勇,劉愛軍,朝魯孟其其格,單玉梅,烏尼圖,陳海軍,王明玖1. 內蒙古農業(yè)大學生態(tài)環(huán)境學院,內蒙古 呼和浩特 010019;. 內蒙古自治區(qū)草原勘察規(guī)劃院,內蒙古 呼和浩特 010051;3. 內蒙古自治區(qū)農牧業(yè)科學院,內蒙古 呼和浩特 010031;4. 中國科學院內蒙古草業(yè)研究中心,內蒙古 呼和浩特 010031;5. 內蒙古自治區(qū)生物技術研究院,內蒙古 呼和浩特 010010

    ?

    錫林郭勒露天煤礦礦區(qū)草原土壤重金屬分布特征

    楊勇1, 2,劉愛軍2,朝魯孟其其格2,單玉梅3, 4,烏尼圖2,陳海軍5*,王明玖1*
    1. 內蒙古農業(yè)大學生態(tài)環(huán)境學院,內蒙古 呼和浩特 010019;2. 內蒙古自治區(qū)草原勘察規(guī)劃院,內蒙古 呼和浩特 010051;3. 內蒙古自治區(qū)農牧業(yè)科學院,內蒙古 呼和浩特 010031;4. 中國科學院內蒙古草業(yè)研究中心,內蒙古 呼和浩特 010031;5. 內蒙古自治區(qū)生物技術研究院,內蒙古 呼和浩特 010010

    摘要:在內蒙古典型草原上,以錫林郭勒勝利煤田西一號和西二號露天礦為研究對象,采用野外樣帶取樣和室內分析相結合的方法,研究了露天煤礦區(qū)土壤重金屬含量(Cr、Cu、Mn、Ni、Zn和Pb元素)的空間分布特征,并探討土壤重金屬來源,旨在為典型草原露天礦區(qū)土壤重金屬防治和草原生態(tài)系統(tǒng)可持續(xù)發(fā)展提供科學理論依據(jù)。結果表明:研究區(qū)土壤重金屬Cr、Cu、Mn、Ni、Zn和Pb的平均質量分數(shù)分別為27.53、16.17、353.49、14.78、44.65和3.44 mg·kg-1,低于國家土壤環(huán)境質量二級標準,僅Cu含量高于內蒙古背景值;礦區(qū)東北方向樣帶,所有采樣點的土壤重金屬含量高于內蒙古背景值;露天礦開采對周圍土壤重金屬空間分布有顯著影響,土壤重金屬元素含量在礦區(qū)中心處最高,并向四周逐漸降低;礦區(qū)周圍0.5 km范圍內土壤重金屬含量均超過內蒙古自治區(qū)背景值。統(tǒng)計分析(相關分析、主成分分析和聚類分析)結果顯示,土壤重金屬Cr、Cu、Mn、Ni和Zn來源可劃歸為一類,Pb單獨為一類。

    關鍵詞:錫林郭勒;露天煤礦;典型草原;土壤重金屬;分布特征

    引用格式:楊勇, 劉愛軍, 朝魯孟其其格, 單玉梅, 烏尼圖, 陳海軍, 王明玖. 錫林郭勒露天煤礦礦區(qū)草原土壤重金屬分布特征[J]. 生態(tài)環(huán)境學報, 2016, 25(5): 885-892.

    YANG Yong, LIU Aijun, CHAO Lumengqiqige, SHAN Yumei, WU Nitu, CHEN Haijun, WANG Mingjiu. Spatial Distribution of Soil Heavy Metals of Opencut Coal Mining in Inner Mongolia Xilingol Typical Steppe [J]. Ecology and Environmental Sciences,2016, 25(5): 885-892.

    煤炭資源是社會發(fā)展和國民經濟的重要物質基礎,在社會經濟發(fā)展中占有非常重要的戰(zhàn)略地位,具有不可替代的作用(白中科等,2006),煤炭開采有力地推動了地區(qū)經濟與社會的發(fā)展,同時也帶來了一系列嚴重的生態(tài)環(huán)境問題(Ferandez-manso et al.,2012)。我國的大型露天開采煤礦大多分布在生態(tài)脆弱的干旱、半干旱草原地區(qū),這些地區(qū)植被覆蓋率較低,水土流失和土地荒漠化十分嚴重。露天煤礦開采對原本脆弱的草原生態(tài)造成嚴重破壞,導致了當?shù)夭菰鷳B(tài)系統(tǒng)的進一步退化,已引起人們的高度關注。礦區(qū)土壤環(huán)境問題已成為我國目前環(huán)境污染領域的研究熱點(白中科等,2006;楊勝香等,2012;李國平等,2013;張廣勝等,2015)。

    國內外學者針對礦區(qū)土壤重金屬進行了大量研究。Panov et al.(1999)對俄羅斯煤田周圍土壤重金屬研究表明,土壤中Hg、As、Pb、Zn和Cd含量超標。Teixeira et al.(2001)對巴西某煤礦區(qū)內河流底部沉積物中的重金屬研究表明,礦區(qū)受到了Pb、Cu、Ni、Fe和Zn等重金屬污染。Benhaddya et al.(2014)研究了阿爾及利亞表層土壤重金屬污染空間分布特征,為研究區(qū)重金屬防治提供了數(shù)據(jù)支持。我國崔龍鵬等(2004)對淮南礦區(qū)土壤重金屬進行了研究,表明采礦活動是礦區(qū)土壤受到重金屬污染的主要原因。楊建等(2008)對河南焦作礦區(qū)演馬礦煤矸石堆周圍土壤中重金屬進行調查研究,發(fā)現(xiàn)研究區(qū)地形地貌特征和風向會影響重金屬的遷移,煤矸石堆周圍土壤重金屬含量超標出現(xiàn)在距離較近處。張連科等(2016)分析了包頭某鋁廠周邊土壤重金屬的空間分布,發(fā)現(xiàn)周邊土壤中重金屬含量高于內蒙古背景值,并具有明顯的空間分布特征。劉碩等(2016)對龍口煤礦區(qū)土壤重金屬污染與空間分布特征進行了研究,結果表明研究區(qū)土壤重金屬綜合污染程度較高。

    目前,我國開展的礦區(qū)重金屬研究主要集中于南方地區(qū),在內蒙古露天煤礦區(qū)僅涉及礦區(qū)植被恢復模式或排土場水土流失的研究(耿海清等,2010;郭建英等,2015),鮮有針對內蒙古典型草原露天煤礦區(qū)開展的較大范圍的土壤重金屬研究。在此背景下,本研究以內蒙古錫林郭勒典型草原露天煤礦區(qū)為研究對象,應用多變量統(tǒng)計方法分析土壤重金屬污染的來源,探討典型草原露天煤礦區(qū)土壤重金屬含量的分布特征,旨在為典型草原露天煤礦區(qū)土壤重金屬污染防治提供基礎數(shù)據(jù)和科學依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    在內蒙古錫林郭勒典型草原上,選取錫林郭勒盟勝利煤田西一號露天煤礦(原烏蘭圖嘎露天煤礦)和西二號露天煤礦(原宏文露天煤礦)為研究對象,分析露天煤礦開采對周邊草原土壤重金屬分布特征的影響。勝利煤田位于內蒙古自治區(qū)錫林郭勒盟錫林浩特市境內,南東邊界距錫林浩特市區(qū)西界約4 km,地理位置為N43°57′~44°14′,E115°30′~116°26′(圖1A)。該區(qū)域氣候屬于溫帶半干旱草原氣候。氣溫日較差和年較差大,冬季寒冷干燥,夏季溫暖濕潤。3—5月常有大風,月平均風速達4.9 m·s-1。根據(jù)當?shù)囟嗄隁庀笥^測資料,多年平均氣溫為0.3 ℃,月均溫度最低為-21.6 ℃(1月),最高為 19.0 ℃(7月)。年均降水量為 336.9 mm,且60%~80%降水集中在生長季的5—8月。無霜期為4月末至 10月初,大約為 150 d。年蒸發(fā)量為1600~1800 mm,相當于降雨量的4~5倍。

    1.2 樣地設置

    為研究露天煤礦開采對周圍土壤重金屬分布特征的影響,首先依據(jù)Google earth遙感影像,通過人工目視解譯確定 2014年勝利煤田西一號露天煤礦和西二號露天煤礦邊界限,并將礦區(qū)及距礦區(qū)邊界8 km范圍內作為具體的重點研究區(qū)域。以礦區(qū)為中心,向礦區(qū)四周7個不同方向輻射設置調查樣線。根據(jù)遙感解譯結果和礦區(qū)實際情況,由于礦區(qū)東邊界距離錫林浩特市4.5 km,在東向設置了長度為4 km的調查樣線;礦區(qū)西北方向有處大型風力發(fā)電場,為避免風力發(fā)電場的干擾,未在礦區(qū)西北方向設置調查樣線。

    圖1 研究區(qū)地理位置(A)和采樣點分布圖(B)Fig. 1 Location map of the study area (A) and sampling locations (B)

    以礦區(qū)為中心,向東(E)、南(S)、西(W)和北(N),以及東北(NE)、東南(SE)和西南(SW)7個輻射方向各設置1條研究樣線,并根據(jù)礦區(qū)實際情況,在礦區(qū)東向樣線上,在距礦區(qū)外圍邊界0、0.5、1、2、4 km處各設置1個調查樣地,在其余6個方向樣帶上,在距礦區(qū)外邊界0、0.5、1、2、4、6、8 km處各設置1個調查樣地(圖1B)。用GPS詳細記錄每個樣地的地理信息。

    1.3 土壤樣品的采集

    2014年7月至 9月中旬,在預設樣地上選擇100 m×100 m區(qū)域進行取樣,并在其對角線上設置1條100 m樣線,在樣線上等距布設6個1 m×1 m調查樣方,分別編號為C1~C6。在C1、C3和C5 3個樣方內采集試驗土樣。取樣時,將土壤表層的殘留物和雜質清理干凈,用直徑5cm土鉆采集表層土樣(0~10 cm)。每個樣方內鉆取5次,均勻混合后裝入塑封袋中,作好標記帶回實驗室。將每份樣品采用四分法取樣,經自然風干后,去除非目標物質(沙礫和植物殘體等),過100目篩子后待分析。

    1.4 土壤重金屬含量測定

    在對煤矸石中重金屬含量分析的基礎上,選擇Cr、Cu、Mn、Ni、Zn和Pb 等6種重金屬元素作為測定對象。全量測定采用HNO3-HCl- HClO4開放式消煮法,用電感耦合等離子體發(fā)射光譜儀( Inductively Coupled Plasma Optical Emission Spectrometer,Optima 2000 DV,PerkinElmer,USA)測定。為保證分析質量,用國家地球化學標準樣進行質量控制,同時進行空白和試劑的校正試驗(鮑士旦,2010)。

    1.5 數(shù)據(jù)分析與處理

    采用Excel 2003對所有數(shù)據(jù)進行錄入及基本統(tǒng)計分析,包括所有數(shù)據(jù)的平均值、最小值、最大值、標準差和變異系數(shù);采用SPSS 16.0進行相關分析、主成分分析和聚類分析;采用SigmaPlot 12.0(Systat Software Inc.)進行作圖。

    2 結果與分析

    2.1 土壤重金屬含量統(tǒng)計分析

    露天煤礦區(qū)土壤Cr、Cu、Mn、Ni、Zn和Pb 6種重金屬含量描述統(tǒng)計分析表明,Mn的平均值最大,其次為Zn,Ni和Pb的平均值最小。6種重金屬元素含量均低于國家土壤環(huán)境質量二級標準(GB 15618—1995)。與內蒙古土壤背景相比,除 Cu高于內蒙古背景值外,其余5種重金屬含量均值均低于內蒙古背景值,但有部分樣點的含量高于內蒙古背景值,說明土壤在一定程度上受到外界干擾。Pb在不同方向變化最大,變異系數(shù)最大,Cr和Mn次之,Cu變異系數(shù)最小,為15.13%??傮w上,礦區(qū)6種土壤重金屬元素的變異系數(shù)在15.13%~31.27%之間,屬于中等變異性,說明這6種重金屬受人類活動影響強度較大(表1)。

    表1 土壤重金屬統(tǒng)計分析Table 1 The statistical feature of soil heavy metal contents

    2.2 不同方向土壤重金屬含量分析

    對露天煤礦區(qū)東、南、西、北、東南、東北和西南7個方向土壤重金屬含量統(tǒng)計分析表明,土壤重金屬在不同方向上含量相差較大。Cr含量在7個方向大小順序依次為:北向>東北>西南>東>西>東南>南。Cu在礦區(qū)東向和東北方向的土壤中平均含量較大,含量最低為礦區(qū)西向,為14.47 mg·kg-1,與最大值相差4.05 mg·kg-1。Mn在礦區(qū)東向和北向的土壤中平均含量較大,在礦區(qū)西向含量最低,為307.01 mg·kg-1,與最大值相差96.27 mg·kg-1;其余方向上平均含量在307.17~386.04 mg·kg-1之間。Ni在礦區(qū)東北方向和東向的土壤中平均含量較大,分別為17.30和17.04 mg·kg-1,最低含量位于為礦區(qū)南向,與最大值相差4.26 mg·kg-1。Zn在礦區(qū)東向和北向的土壤中平均含量最大,分別為 52.03和47.69 mg·kg-1,平均最低含量為礦區(qū)東南向土壤,為40.63 mg·kg-1,與最大值相差11.40 mg·kg-1。Pb在礦區(qū)西向和東北方向平均含量較大,平均最低含量為礦區(qū)東向土壤,為2.78 mg·kg-1,與最大值相差1.60 mg·kg-1(表2)。

    表2 不同方向土壤重金屬含量Table 2 Soil heavy metals contents in different sampling directions               mg·kg-1

    2.3 不同采樣距離土壤重金屬含量變化

    為了分析露天煤礦對周圍土壤重金屬空間分布的影響,計算7個方向上,土壤Cr、Cu、Mn、Ni、Zn和Pb 6種重金屬含量隨距礦區(qū)距離增加的變化特征。結果表明,總體上Cr、Cu、Mn、Ni和Zn 5種重金屬元素隨著距礦區(qū)距離的增加,其含量逐漸減小,尤其0~2 km范圍內,下降較快。從2 km后,基本無明顯變化。Pb隨著距礦區(qū)距離的增加,無明顯變化規(guī)律,基本呈上下波動變化。但對每種重金屬元素,其含量在不同方向上隨著采樣點距礦區(qū)距離的增加,變化也不相同(圖2)。

    Cr元素含量在礦區(qū)東向、東北方向上隨著距礦區(qū)距離的增加,呈現(xiàn)先減小后增大的趨勢。在東向2 km和東北方向1 km處,Cr含量達到最低值,隨后 Cr含量隨著距礦區(qū)距離的增加而增加。在礦區(qū)南向、西向、北向和西南,Cr含量基本表現(xiàn)為在0~4 km范圍內逐漸下降,下降幅度介于19.26%~35.15%之間。在礦區(qū)東南方向上,變化較復雜(圖2A)。

    Cu元素含量在不同方向都為距礦區(qū)0 km處最高。在礦區(qū)東向、東南和東北方向上隨著距礦區(qū)距離的增加,呈現(xiàn)先減小后增大的趨勢,在東向和東北方向2 km處達到最低值,在東南4 km處達到最低值。在礦區(qū)南向、西向、北向和西南4 km處,較0 km下降了23.51%~36.73%,隨后出現(xiàn)上下波動,變化幅度較?。▓D2B)。

    Mn元素含量在礦區(qū)東向、東北方向上隨著距礦區(qū)距離的增加,呈現(xiàn)先減小后增大的趨勢,在距礦區(qū)2 km處達到最低值。在南向、北向、東南和西南0~2 km范圍內呈直線下降趨勢,隨后在波動變化中逐漸下降。在礦區(qū)西向,在距礦區(qū)0.5 km處出現(xiàn)下降,隨后逐漸增加,在2 km處達到第2個峰值(414.22 mg·kg-1),隨后在4 km處出現(xiàn)大幅度下降,下降幅度達54.58%,在6~8 km范圍內變化幅度較小(圖2C)。

    Ni元素含量在東向呈現(xiàn)先減小后增大的趨勢,在2 km處降到最低。在南向和東北向呈先降低后升高的變化趨勢,在南向2 km處和東北向處最低,之后有小幅升高趨勢。在西向和北向0~4 km范圍內下降幅度較大,之后變化較小。在東南和西南呈不斷下降趨勢,在0~2 km下降較快,之后趨于平緩(圖2D)。

    圖2 土壤重金屬含量隨采樣距離的變化Fig. 2 Variation of soil heavy metal contents with sampling distance

    Zn元素含量在礦區(qū)東向、東南、東北和西南方向上隨著距礦區(qū)距離的增加,呈現(xiàn)先減小后增大的趨勢,在2 km處達到最低值,隨后在東向和東北方向上呈大幅度升高趨勢,在東南和西南方向上有小幅度升高。在南向呈上下波動變化。在西向最大值出現(xiàn)在0.5 km處,在4 km處下降了41.35%,隨后出現(xiàn)升高—降低—升高的波動變化。在北向 0~1 km,出現(xiàn)較大幅度降低,隨后呈增大—減小的變化趨勢,但變化范圍較?。▓D2E)。

    Pb含量在東向隨著距離的增加無顯著變化規(guī)律。在南向呈波動式逐漸降低的變化趨勢。在西向呈現(xiàn)先增大后減小的拋物線趨勢。北向呈線性減小趨勢,減小速率較低。在西南0~1 km范圍內逐漸增大,隨后開始降低,在4 km處降到最低,4~8 km有小幅度增大趨勢。在礦區(qū)東北方向上呈現(xiàn)先減小后增大的趨勢,在1 km處達到最低值。在東南向呈上下波動變化,并在4~8 km范圍內達到穩(wěn)定(圖2F)。

    2.4 土壤重金屬來源解析

    2.4.1 土壤重金屬元素間的相關分析

    礦區(qū)土壤中不同重金屬元素含量之間的相關性分析表明,Pb與Cr、Cu、Mn、Ni、Zn間均無顯著相關性(P>0.05)。Cr、Cu、Mn、Ni和Zn兩兩間呈極顯著正相關性(P<0.01),其中 Cr和 Zn相關系數(shù)最低,Cr和Ni相關性較高,相關系數(shù)達0.855,相關系數(shù)數(shù)值在0.583~0.855之間(表3)??梢猿醪酵茢?,土壤中Cr、Cu、Mn、Ni和Zn這5種元素的地球化學性質相近,在相同的外界環(huán)境條件下其變化趨勢基本一致,具有相同的來源,極可能來源于露天煤礦。而Pb與以上元素不同,為單獨來源。

    表3 土壤重金屬含量的相關系數(shù)Table 3 Pearson′s correlation coefficients (r) among soil heavy metal contents

    2.4.2 土壤重金屬含量特征的主成分分析

    運用最大方差旋轉方法對土壤中6種重金屬含量進行主成分分析(Principal component analysis,PCA),由初始特征值可以看出,前 2個主成分的累積方差貢獻率已經達到 79.41%,可以解釋接近80%的總方差。旋轉前提取的第1個主成分的特征根大于1,旋轉后提取的2個主成分的特征根都大于1(表4)。因此,土壤中所有重金屬元素可以提取為2個主成分,這2個主成分的累積方差貢獻率達79.41%。

    根據(jù)各主成分旋轉載荷矩陣可以得出,第1主成分主要反映Cr、Cu、Mn、Ni和Zn的組成信息,貢獻率為60.15%;第2主成分主要反映Pb的富集信息,貢獻率為19.26%(表5)。

    表4 土壤重金屬含量主成分方差貢獻Table 4 The contribution of PCA variance for soil heavy metal contents

    表5 土壤重金屬含量主成分分析Table 5 The analysis of PCA of soil heavy metal contents

    2.4.3 土壤重金屬含量的聚類分析

    聚類分析指將物理或抽象對象的集合分組為由類似的對象組成的多個類的分析過程,是目前被廣泛應用于環(huán)境領域的一種多元統(tǒng)計方法。系統(tǒng)聚類法(分層聚類)是目前應用較多的一種聚類分析方法。通過分層聚類樹狀圖可以直觀地反映出土壤元素間的相關性或親屬關系,可有效揭示土壤重金屬污染物的來源。本文以相關系數(shù)為距離測度方法,采用組間連接方法,對土壤中6種重金屬進行聚類分析,其聚類樹狀圖形象、直觀地反映出了 6種重金屬元素的親屬關系。土壤重金屬可分為二聚類,第一聚類包括Cr、Cu、Mn、Ni和Zn,第二聚類為Pb(圖 3)。聚類分析的結果和主成分分析的結果保持一致,可以認為Cr、Cu、Mn、Ni和 Zn來源相同,Pb則有單獨來源。

    圖3 土壤重金屬元素聚類分析Fig. 3 Hierarchical cluster analysis of soil heavy metal concentrations

    3 討論

    土壤重金屬空間分布的影響因素較為復雜,除受土壤母質影響外,還受外來人為因素的干擾(Spijker et al.,2005)。礦業(yè)活動是造成礦區(qū)周邊土壤中重金屬含量增加的重要因素,但由于擴散條件不同,會導致土壤中重金屬在污染源不同方向上含量相差較大。孫葉芳等(2005)以浙江某鉛鋅礦區(qū)為研究對象,對污染土壤中的重金屬(Cu、Zn、Pb、Cd)進行測定,并采用美國最新的重金屬污染評價方 法TCLP( Toxicity characteristic leaching procedure)對土壤重金屬污染狀況進行評價,結果表明,礦區(qū)附近土壤不同程度地受到Cu、Zn、Pb、Cd的污染,且土壤重金屬在不同方向上(東、南、西、北、東南、東北和西南)含量相差較大,處于煤矸石堆下風向和地下水下游方向,土壤中重金屬元素含量較高。

    本研究表明,露天煤礦區(qū)土壤重金屬含量在不同方向相差較大,基本上表現(xiàn)出東、東北、北3個方向土壤重金屬含量較高,這與研究區(qū)主導風向有關。礦區(qū)周圍無明顯地面河流,錫林河在礦區(qū)東向已斷流,污染物主要以大氣及粉塵擴散為主。根據(jù)錫林浩特市多年氣象資料顯示,礦區(qū)年主導風向為西南風,春季和冬季主導風向為西南風,夏季主導風向為東南風,秋季主導風向為西南風。因此造成處于下風向的礦區(qū)北向和東北向土壤重金屬含量較大。同時本研究在礦區(qū)東向設置的采樣點,除受到主導風向的影響外,也受到錫林浩特市的影響。東向E4采樣點,已達錫林浩特市近郊,土壤中Cr、Cu、Mn、Ni、Zn和Pb 6種重金屬含量都有升高趨勢,可能與該樣點距離錫林浩特市較近有關。因為礦業(yè)和城市活動為主的人為干擾是影響土壤重金屬空間分布的重要原因(李艷霞等,2007)。人類活動較為密集的地區(qū)土壤重金屬富集程度也較高(付傳城等,2014)。城市發(fā)展會對城市及周邊土壤造成重金屬污染(柳云龍等,2012),隨著城市人口的不斷增加、城鎮(zhèn)工業(yè)化的不斷發(fā)展、城市機動車保有輛的不斷增加,城市帶來的土壤重金屬污染將更加嚴峻(陳秀端等,2011)。

    露天煤礦開采會向外排放大量的煤矸石,煤矸石中的重金屬在長期堆積、風化和淋溶作用下向土壤遷移,并出現(xiàn)一定的累積效應(劉玉榮等,2003)。煤矸石對土壤造成的重金屬污染主要有兩種途徑:一是經風蝕后的矸石揚塵懸浮于大氣中,隨風降落于矸石堆周圍土壤;二是矸石受大氣降水的沖刷和淋溶作用而使重金屬隨地表徑流進入土壤(陳峰等,2006)。因此,其污染范圍與遷移途徑和遷移條件有關。王心義等(2006)對河南焦作礦區(qū)演馬礦煤矸石堆周圍土壤中重金屬的研究發(fā)現(xiàn),煤矸石含Zn 123.89 mg·kg-1,但在距離煤矸石山340 m內土壤含量 Zn超過 2900 mg·kg-1,最高達 9400 mg·kg-1。研究區(qū)地形地貌特征和風向會影響重金屬的遷移。在本研究中,土壤重金屬含量在礦區(qū)中心較高,高于內蒙古土壤背景值的樣點分布于距礦區(qū)0.5 km范圍內,且隨著距礦區(qū)距離的增加而逐漸減小,基本與上述結果相一致。同類其他研究也表明煤矸石山周邊土壤重金屬含量隨著距離煤矸石山距離的增加而減?。ㄌ锊氏嫉?,2007;魏忠義等,2008)。在本研究中,受到主導風向和城市建設的影響,Ni和Zn在礦區(qū)東向4 km處和東北向6 km、8 km處存在土壤重金屬含量高于內蒙古土壤背景值的情況。

    礦區(qū)土壤重金屬的來源分析,可為防控礦區(qū)土壤重金屬污染和周邊草地系統(tǒng)的有效管理提供科學信息。目前,相關分析、主成分分析和聚類分析被廣泛應用于土壤重金屬來源分析的研究中(呂建樹等,2012;謝小進等,2010;Mostert et al.,2010;Chabukdhara et al.,2012)。一般情況下,如果土壤中不同重金屬含量之間呈顯著相關性,可以說明金屬間具有同源關系或是復合污染(柴世偉等,2004;Yan et al.,2010)。在本研究中,Pb與Cr無顯著相關性(P>0.05),與Cu、Mn、Zn和Ni呈顯著負相關性(P<0.05),Cr、Cu、Mn、Ni和Zn兩兩之間呈極顯著正相關性(P<0.01)。通過主成分分析和聚類分析也表明,土壤重金屬Cr、Cu、Mn、Ni和Zn可歸為一類,另一類為 Pb,礦區(qū)周圍土壤重金屬Cr、Cu、Mn、Ni和Zn來源相同,可能來源于礦區(qū)開采,Pb則有單獨來源,可能主要來源于土壤背景或道路車輛,包括汽車尾氣排放、工業(yè)煙塵、Pb和As殺蟲氣和污泥施用等(王軍等,2007;Besnard et al.,2001)。謝小進等(2010)對上海農田土壤重金屬分布的研究中發(fā)現(xiàn),Pb的積聚與交通業(yè)有密切的關系,高值區(qū)集中分布在路網密集的區(qū)域。

    4 結論

    (1)研究區(qū)6種土壤重金屬平均含量低于國家土壤環(huán)境質量二級標準,僅有Cu含量高于內蒙古自治區(qū)背景值。礦區(qū)東北方向樣帶,所有采樣點的土壤重金屬含量高于內蒙古背景值。

    (2)土壤重金屬元素含量在礦區(qū)中心處最高,并向四周逐漸降低。礦區(qū)周圍0.5 km范圍內土壤重金屬含量均超過于內蒙古自治區(qū)背景值。

    (3)土壤重金屬Cr、Cu、Mn、Ni和Zn元素來源可歸為一類,Pb單獨為一類。

    參考文獻:

    BENHADDYA M L, HADJEL M. 2014. Spatial distribution and contamination assessment of heavy metals in surface soils of Hassi Messaoud, Algeria [J]. Environmental Earth Sciences, 71(3): 1473-1486.

    BESNARD E, CHENU C, ROBERT M. 2001. Influence of organic amendments on copper distribution among particle-size and density fractions in Champagne vineyard soils [J]. Environmental Pollution,112(3): 329-337.

    CHABUKDHARA M, NEMA A K. 2012. Assessment of heavy metal contamination in Hindon River sediments: A chemometric and geochemical approach [J]. Chemosphere, 87(8): 945-953.

    FERANDEZ-MANSO A, QUINTANO C, ROBERTS D A. 2012. Evaluation of potential of multiple endmember spectral mixture analysis (MESMA) for surface coal mining affected area mapping in different world forest ecosystems [J]. Remote Sensing of Environment,127(12): 181-193.

    MOSTERT M R, AYOKO G A, KOKOT S. 2010. Application of chemometrics to analysis of soil pollutants [J]. Trac Trends in Analytical Chemistry, 29(5): 430-445.

    PANOV B S, DUDIK A M, SHEVCHENKO O A, et al. 1999. On pollution of the biosphere in industrial areas: the example of the donets coal basin [J]. International Journal of Coal Geology, 40(2-3): 199-210.

    SPIJKER J, VRIEND S P, GAANS P F M V. 2005. Natural and anthropogenic patterns of covariance and spatial variability of minor and trace elements in agricultural topsoil [J]. Geoderma, 127(1-2): 24-35.

    TEIXEIRA E, ORTIZL L, ALVES M, et al. 2001. Distribution of selected heavy metals in fluvial sediments of coal mining region of Baixo Jacui,RS, Brazil [J]. Environmental Geology, 41(1): 145-154.

    YAN C Z, LI Q Z, ZHANG X, et al. 2010. Mobility and ecological risk assessment of heavy metals in surface sediments of Xia-men Bay and its adjacent areas, China [J]. Environmental Earth Sciences, 60(7): 1469-1479.

    GB 15618—1995, 土壤環(huán)境質量標準[S].

    白中科, 付梅臣, 趙中秋. 2006. 論礦區(qū)土壤環(huán)境問題[J]. 生態(tài)環(huán)境,15(5): 1122-1125.

    白中科, 耿海青, 郭二民, 等. 2006. 關于煤炭開發(fā)生態(tài)補償?shù)娜舾梢庖姡跩]. 環(huán)境保護, (9): 46-48.

    鮑士旦. 2010. 土壤農化分析[M]. 北京: 中國農業(yè)出版社.

    柴世偉, 溫琰茂, 韋獻革, 等. 2004. 珠江三角洲主要城市郊區(qū)農業(yè)土壤的重金屬含量特征[J]. 中山大學學報(自然科學版), 43(4): 90-94.

    陳峰, 胡振琪, 柏玉, 等. 2006. 矸石山周圍土壤重金屬污染的生態(tài)風險評價[J]. 農業(yè)環(huán)境科學學報, 25(增刊): 575-578.

    陳秀端, 盧新衛(wèi), 趙彩鳳, 等. 2011. 西安市二環(huán)內表層土壤重金屬空間分布特征[J]. 地理學報, 66(9): 1281-1288.

    崔龍鵬, 白建峰, 史永紅, 等. 2004. 采礦活動對煤礦區(qū)土壤中重金屬污染研究[J]. 土壤學報, 41(6): 896-904.

    付傳城, 王文勇, 潘劍君, 等. 2014. 城鄉(xiāng)結合帶土壤重金屬時空變異特征與源解析-以南京市柘塘鎮(zhèn)為例[J]. 土壤學報, 51(5): 1066-1077.

    耿海清, 陳帆, 劉杰, 等. 2010. 煤炭富集區(qū)開發(fā)模式解析—以錫林郭勒盟為例[J]. 地域研究與開發(fā), 29(4): 32-37.

    郭建英, 何京麗, 李錦榮, 等. 2015. 典型草原大型露天煤礦排土場邊坡水蝕控制效果[J]. 農業(yè)工程學報, 31(3): 296-303.

    李國平, 郭江. 2013. 能源資源富集區(qū)生態(tài)環(huán)境治理問題研究[J]. 中國人口·資源與環(huán)境, 23(7): 42-48.

    李艷霞, 徐理超, 熊雄, 等. 2007. 典型礦業(yè)城市農田土壤重金屬含量的空間結構特征-以遼寧省阜新市為例[J]. 環(huán)境科學學報, 27(4): 679-687.

    劉碩, 吳泉源, 曹學江. 2016. 龍口煤礦區(qū)土壤重金屬污染評價與空間分布特征[J]. 環(huán)境科學, 37(1): 270-279.

    劉玉榮, 黨志, 尚愛安. 2003. 煤矸石風化土壤中重金屬的環(huán)境效應研究[J]. 農業(yè)環(huán)境科學學報, 22(1): 64-66.

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

    呂建樹, 張祖陸, 劉洋, 等. 2012. 日照市土壤重金屬來源解析及環(huán)境風險評價[J]. 地理學報, 67(7): 971-984.

    孫葉芳, 謝正苗, 徐建明, 等. 2005. TCLP法評價礦區(qū)土壤重金屬的生態(tài)環(huán)境風險[J]. 環(huán)境科學, 26(3): 153-156.

    田彩霞, 郭保華. 2007. 煤矸石堆對周圍土壤中重金屬元素影響分析[J].礦業(yè)安全與環(huán)保, 34(3): 23-25.

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

    王心義, 楊建, 郭慧霞. 2006. 礦區(qū)煤矸石堆放引起土壤重金屬污染研究[J]. 煤炭學報, 31(6): 808-812.

    魏忠義, 陸亮, 王秋兵. 2008. 撫順西露天礦大型煤矸石山及其周邊土壤重金屬污染研究[J]. 土壤通報, 39(4): 946-949.

    謝小進, 康建成, 李衛(wèi)江, 等. 2010. 上海寶山區(qū)農用土壤重金屬分布與來源分析[J]. 環(huán)境科學, 31(3): 768-774.

    謝小進, 康建成, 閆國東, 等. 2010. 黃浦江中上游地區(qū)農用土壤重金屬含量特征分析[J]. 中國環(huán)境科學, 30(8): 1110-1117.

    楊建, 陳家軍, 王心義. 2008. 煤矸石堆周圍土壤重金屬污染空間分布及評價[J]. 農業(yè)環(huán)境科學學報, 27(3): 873-878.

    楊勝香, 袁志忠, 李朝陽, 等. 2012. 湘西花垣礦區(qū)土壤重金屬污染及其生物有效性[J]. 環(huán)境科學, 33(5): 1718-1724.

    張廣勝, 徐文彬, 李俊翔, 等. 2015. 一個未開采的鉛鋅礦周邊土壤重金屬含量及生態(tài)安全評價[J]. 生態(tài)環(huán)境學報, 24(3): 522-528.

    張連科, 李海鵬, 黃學敏, 等. 2016. 包頭某鋁廠周邊土壤重金屬的空間分布及來源解析[J]. 環(huán)境科學, 37(3): 1139-1146.

    DOI:10.16258/j.cnki.1674-5906.2016.05.023

    中圖分類號:X53

    文獻標志碼:A

    文章編號:1674-5906(2016)05-0885-08

    基金項目:內蒙古科技創(chuàng)新引導獎勵資金計劃“大型露天煤礦生態(tài)修復技術集成與示范”(20140701);內蒙古自治區(qū)科技計劃項目“不同生態(tài)類型區(qū)域生態(tài)評估及風險預測技術研究”(2013年立項);內蒙古自治區(qū)自然科學基金項目(2015MS3084)

    作者簡介:楊勇(1984年生),男,助理研究員,博士研究生,研究方向為草地生態(tài)系統(tǒng)結構與功能。E-mail: yangyong606@Gmail.com

    *通信作者。E-mail: wangmj_0540@163.com; chenhaijun2004@163.com

    收稿日期:2016-04-15

    Spatial Distribution of Soil Heavy Metals of Opencut Coal Mining in Inner Mongolia Xilingol Typical Steppe

    YANG Yong1, 2, LIU Aijun2, CHAO Lumengqiqige2, SHAN Yumei3, 4, WU Nitu2,CHEN Haijun3*, WANG Mingjiu1*
    1. College of Ecology and Environmental Science, Inner Mongolia Agricultural University, Hohhot 010019, China;
    2. Inner Mongolia Institute of Grassland Survey and Planning, Hohhot 010051, China;
    3. Inner Mongolia Academy of Agriculture & Animal Husbandry Science, Hohhot 010031, China;
    4. Inner Mongolia Research Center for Prataculture, Chinese Academy of Science, Hohhot 010031, China;
    5. Inner Mongolia Institute of Biotechnology, Hohhot 010010, China

    Abstract:In this study, the Spatial distribution of Soil Heavy Metals of Opencut Coal Mining was analized in the two Shengli Coal Field in Xilingol typical steppe by the approach of combining field survey and lobarotory analyses. The sources of heavy metals in the mining area was discussed to provide evidence for preventing soil heavy metal pollution and for grassland ecosystem safety. The results indicated that the average concentrations of Cr, Cu, Mn, Ni, Zn, and Pb in the soil of the study area was 27.53, 16.17, 353.49,14.78, 44.65 and 3.44 mg·kg-1, respectively, being lower than the second class level of the National Soil Environment Safety Standard. But the average concentration of the Cu was higher than the background value of soil in the Inner Mongolia. The open-cut coal mining exerted a significant impact on the spatial distribution of heavy metals in the soil. The content of soil heavy metal elements was the highest in the center of the mining area, and gradually decreased to the surrounding areas, and was higher in the northeastern of the mining area. Within 0.5 km from the mining area, the content of soil heavy metals exceeded soil background values in the Inner Mongolia. The results of correlation analysis, principal component analysis and cluster analysis showed that Cr,Cu, Mn, Ni and Zn in the mining area come from the same source, while Pb came from the other source.

    Key words:Xilingol; opencut coal mining; typical steppe; soil heavy metal; distribution characteristics

    猜你喜歡
    土壤重金屬錫林郭勒分布特征
    內蒙古錫林郭勒那達慕大會 草原夏日
    天籟之音 生命之歌——錫林郭勒長調民歌編譯
    草原歌聲(2020年1期)2020-07-25 01:45:18
    天堂草原錫林郭勒
    草原歌聲(2017年3期)2017-04-23 05:13:46
    嶺南地區(qū)冠心病中醫(yī)證型及證素分布的地域性特征研究
    四川宜賓瀘州長江河谷地帶雷電監(jiān)測數(shù)據(jù)特征分析
    基于Kriging插值的礦區(qū)周邊土壤重金屬空間分布規(guī)律研究
    科技視界(2016年13期)2016-06-13 08:06:36
    貴州擺紀磷石膏尾礦庫周邊土壤重金屬污染評價
    五壘島灣海域石油烴分布特征
    某地區(qū)表層土壤重金屬分布特征研究
    池塘淤泥氮營養(yǎng)物質分布特征研究
    免费一级毛片在线播放高清视频| 国产私拍福利视频在线观看| 国产精品伦人一区二区| 亚洲成人久久性| 婷婷精品国产亚洲av| 午夜爱爱视频在线播放| 亚洲精品国产成人久久av| 天天躁夜夜躁狠狠久久av| 亚洲,欧美,日韩| 高清在线视频一区二区三区 | 亚洲中文字幕日韩| 亚洲三级黄色毛片| 亚洲欧美精品综合久久99| 一本一本综合久久| 欧美zozozo另类| 国产黄a三级三级三级人| av国产免费在线观看| 国产精品综合久久久久久久免费| 婷婷亚洲欧美| 中文资源天堂在线| 老司机福利观看| 午夜激情欧美在线| 天天躁日日操中文字幕| 久久久久久久久大av| 欧美又色又爽又黄视频| 99视频精品全部免费 在线| 欧美人与善性xxx| 级片在线观看| 性欧美人与动物交配| 青青草视频在线视频观看| 精品午夜福利在线看| 少妇被粗大猛烈的视频| 午夜福利在线观看免费完整高清在 | 国产av一区在线观看免费| 中文字幕精品亚洲无线码一区| 成人无遮挡网站| 一个人免费在线观看电影| 亚洲一级一片aⅴ在线观看| 日本与韩国留学比较| 国产精品久久久久久久电影| 午夜精品一区二区三区免费看| 国产乱人偷精品视频| 欧美另类亚洲清纯唯美| 日本黄色片子视频| 亚洲一区高清亚洲精品| 老司机影院成人| 在线观看午夜福利视频| 一边摸一边抽搐一进一小说| 中国美女看黄片| av在线播放精品| 国产成人福利小说| 亚洲av.av天堂| av.在线天堂| 久久久久网色| 久久久久久大精品| 99热只有精品国产| 欧美bdsm另类| av在线观看视频网站免费| 日本撒尿小便嘘嘘汇集6| 久久久久久久久大av| 亚洲最大成人中文| 99热精品在线国产| 别揉我奶头 嗯啊视频| 亚洲精品456在线播放app| 秋霞在线观看毛片| 美女 人体艺术 gogo| 一本一本综合久久| 亚洲国产色片| 国产视频内射| 麻豆乱淫一区二区| 在线免费观看不下载黄p国产| 久久久色成人| 插阴视频在线观看视频| 一本久久精品| 亚洲欧美日韩卡通动漫| 舔av片在线| 亚洲无线在线观看| 99热只有精品国产| 久久精品夜夜夜夜夜久久蜜豆| 久久99蜜桃精品久久| 精品人妻偷拍中文字幕| 免费观看在线日韩| 99久久九九国产精品国产免费| 国产伦在线观看视频一区| 国产白丝娇喘喷水9色精品| 国产伦在线观看视频一区| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人精品欧美一级黄| 亚洲人与动物交配视频| 午夜精品在线福利| 日韩欧美精品免费久久| 欧美成人a在线观看| 搡老妇女老女人老熟妇| 深爱激情五月婷婷| 淫秽高清视频在线观看| 长腿黑丝高跟| 日韩三级伦理在线观看| 丰满人妻一区二区三区视频av| 欧美日本亚洲视频在线播放| 身体一侧抽搐| 亚洲国产高清在线一区二区三| 免费av不卡在线播放| 国产熟女欧美一区二区| 黄色一级大片看看| 国产免费男女视频| 久久久久国产网址| 18+在线观看网站| 99热精品在线国产| 日韩大尺度精品在线看网址| 午夜久久久久精精品| 天堂网av新在线| 又爽又黄a免费视频| 国产精品久久久久久久久免| 久久草成人影院| 成人二区视频| 久久午夜亚洲精品久久| 97热精品久久久久久| videossex国产| 久久九九热精品免费| 国产精品久久久久久精品电影小说 | 精品久久久久久久久av| АⅤ资源中文在线天堂| 91麻豆精品激情在线观看国产| 久久99热这里只有精品18| 日韩中字成人| 男女下面进入的视频免费午夜| 亚洲欧美成人综合另类久久久 | av.在线天堂| 悠悠久久av| 欧美一区二区精品小视频在线| 国内少妇人妻偷人精品xxx网站| 成人毛片a级毛片在线播放| 我的老师免费观看完整版| 日本免费一区二区三区高清不卡| 波多野结衣高清无吗| 久久精品国产亚洲av香蕉五月| 久久中文看片网| 久久久久久久午夜电影| 日韩av在线大香蕉| 国产视频内射| 精品国内亚洲2022精品成人| 成人高潮视频无遮挡免费网站| 中国美白少妇内射xxxbb| 啦啦啦韩国在线观看视频| 国产伦在线观看视频一区| 免费观看在线日韩| 波多野结衣高清无吗| 亚洲美女视频黄频| 国产视频首页在线观看| 嫩草影院入口| 非洲黑人性xxxx精品又粗又长| 婷婷色av中文字幕| 99热这里只有是精品50| 亚洲欧美成人精品一区二区| 变态另类成人亚洲欧美熟女| 不卡一级毛片| 日产精品乱码卡一卡2卡三| 一区福利在线观看| 人人妻人人澡欧美一区二区| 国产精品久久久久久亚洲av鲁大| 看非洲黑人一级黄片| 高清日韩中文字幕在线| 亚洲国产欧洲综合997久久,| 我要看日韩黄色一级片| av卡一久久| 国产麻豆成人av免费视频| 国产精品久久久久久久久免| 欧美潮喷喷水| 三级经典国产精品| 日韩欧美 国产精品| 欧美bdsm另类| 麻豆国产97在线/欧美| 三级国产精品欧美在线观看| 亚洲人成网站在线播放欧美日韩| 欧美日本亚洲视频在线播放| 99久久中文字幕三级久久日本| 一级毛片久久久久久久久女| 日韩亚洲欧美综合| 国产精品蜜桃在线观看 | 人人妻人人看人人澡| 国产精品一区二区性色av| 男女做爰动态图高潮gif福利片| 国产高潮美女av| 亚洲在线自拍视频| 秋霞在线观看毛片| 国产综合懂色| 国产午夜精品一二区理论片| 全区人妻精品视频| av国产免费在线观看| 国产精品三级大全| 六月丁香七月| 一级黄片播放器| 最近最新中文字幕大全电影3| 国内精品一区二区在线观看| 国产精品综合久久久久久久免费| 三级男女做爰猛烈吃奶摸视频| 国产乱人视频| 看非洲黑人一级黄片| 热99在线观看视频| 深夜精品福利| 亚洲精品粉嫩美女一区| 日本av手机在线免费观看| 久久久色成人| 人妻久久中文字幕网| 91狼人影院| 我的老师免费观看完整版| 99热网站在线观看| 美女xxoo啪啪120秒动态图| 日本欧美国产在线视频| 亚洲国产精品久久男人天堂| 在线播放无遮挡| 天堂中文最新版在线下载 | 久久精品国产亚洲av涩爱 | 免费观看的影片在线观看| 国产黄a三级三级三级人| 亚洲va在线va天堂va国产| 亚洲av电影不卡..在线观看| 亚洲成av人片在线播放无| 麻豆久久精品国产亚洲av| 国产伦精品一区二区三区视频9| 男女边吃奶边做爰视频| 亚洲欧美日韩高清专用| 国产精品麻豆人妻色哟哟久久 | 日韩av在线大香蕉| 麻豆av噜噜一区二区三区| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区免费观看| 亚洲精品亚洲一区二区| 丰满的人妻完整版| av在线播放精品| 干丝袜人妻中文字幕| 欧美变态另类bdsm刘玥| 国产女主播在线喷水免费视频网站 | 18禁在线播放成人免费| 好男人在线观看高清免费视频| 亚洲天堂国产精品一区在线| 少妇高潮的动态图| 91午夜精品亚洲一区二区三区| 男女做爰动态图高潮gif福利片| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 麻豆乱淫一区二区| 亚洲在线自拍视频| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在 | 丰满人妻一区二区三区视频av| 69av精品久久久久久| 国产片特级美女逼逼视频| 国产麻豆成人av免费视频| 午夜精品一区二区三区免费看| 免费av不卡在线播放| 亚洲自拍偷在线| 日韩成人av中文字幕在线观看| 韩国av在线不卡| 日韩视频在线欧美| 永久网站在线| 亚洲综合色惰| 69av精品久久久久久| 少妇人妻精品综合一区二区 | 亚洲经典国产精华液单| 欧美日本亚洲视频在线播放| 亚洲精品色激情综合| 亚洲av熟女| 夜夜夜夜夜久久久久| 国产成人精品一,二区 | 亚洲无线在线观看| 国模一区二区三区四区视频| 亚洲欧美精品综合久久99| 波多野结衣巨乳人妻| 人人妻人人澡人人爽人人夜夜 | 日韩欧美精品免费久久| 欧美高清性xxxxhd video| 亚洲,欧美,日韩| 亚洲精品日韩av片在线观看| 欧美性猛交╳xxx乱大交人| 乱系列少妇在线播放| 日本一本二区三区精品| 国产极品精品免费视频能看的| 国产精品久久久久久精品电影| 免费搜索国产男女视频| 国产av一区在线观看免费| 亚洲欧美日韩东京热| 日产精品乱码卡一卡2卡三| 日韩三级伦理在线观看| 99久久无色码亚洲精品果冻| 免费黄网站久久成人精品| 国产免费男女视频| 国产三级在线视频| 国产一区二区三区在线臀色熟女| 亚洲欧美中文字幕日韩二区| 精品人妻一区二区三区麻豆| 婷婷亚洲欧美| 久久久成人免费电影| 在线观看午夜福利视频| 成人鲁丝片一二三区免费| 久久精品国产亚洲av涩爱 | 一本一本综合久久| 国产探花在线观看一区二区| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产| 日韩欧美精品免费久久| 久久人妻av系列| av.在线天堂| 国产精品国产三级国产av玫瑰| 97热精品久久久久久| 直男gayav资源| 亚洲经典国产精华液单| 男人舔女人下体高潮全视频| 色播亚洲综合网| 高清毛片免费看| 国产极品精品免费视频能看的| 九色成人免费人妻av| 欧美色视频一区免费| 精品久久久久久久久久免费视频| 亚洲人成网站在线播放欧美日韩| 亚洲av男天堂| 一本久久中文字幕| av在线老鸭窝| 国产精品综合久久久久久久免费| 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 国产伦在线观看视频一区| 少妇被粗大猛烈的视频| 岛国在线免费视频观看| 成年av动漫网址| 变态另类成人亚洲欧美熟女| 美女黄网站色视频| 给我免费播放毛片高清在线观看| 国产又黄又爽又无遮挡在线| 99视频精品全部免费 在线| 97超视频在线观看视频| 日韩av在线大香蕉| 国产精品永久免费网站| 亚洲经典国产精华液单| 日日撸夜夜添| 日韩欧美精品免费久久| 国产大屁股一区二区在线视频| 永久网站在线| 床上黄色一级片| 久久人人爽人人片av| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 亚洲精品影视一区二区三区av| 丝袜喷水一区| 欧美精品一区二区大全| a级毛片a级免费在线| a级毛色黄片| 大香蕉久久网| 日韩欧美精品v在线| 成人特级av手机在线观看| 中文字幕精品亚洲无线码一区| 日日摸夜夜添夜夜添av毛片| 少妇人妻一区二区三区视频| 两个人视频免费观看高清| av黄色大香蕉| 国产精品综合久久久久久久免费| 18禁在线播放成人免费| 人妻久久中文字幕网| av天堂在线播放| 狂野欧美激情性xxxx在线观看| 亚洲va在线va天堂va国产| 99久久人妻综合| 插逼视频在线观看| 我的女老师完整版在线观看| 日韩亚洲欧美综合| 国产av麻豆久久久久久久| 日产精品乱码卡一卡2卡三| 两个人视频免费观看高清| 99热这里只有精品一区| 九九在线视频观看精品| 免费无遮挡裸体视频| 久久婷婷人人爽人人干人人爱| 日本一二三区视频观看| 国产 一区 欧美 日韩| 国产精品久久久久久久电影| 国产av不卡久久| 美女cb高潮喷水在线观看| 免费一级毛片在线播放高清视频| 日本成人三级电影网站| 黑人高潮一二区| 三级毛片av免费| 欧美人与善性xxx| 国产精品福利在线免费观看| 免费搜索国产男女视频| 国产色爽女视频免费观看| 日韩欧美 国产精品| 舔av片在线| 嘟嘟电影网在线观看| 人妻制服诱惑在线中文字幕| 高清日韩中文字幕在线| 毛片女人毛片| 久久久欧美国产精品| 深夜a级毛片| 亚洲aⅴ乱码一区二区在线播放| 国产在视频线在精品| 三级国产精品欧美在线观看| 久久综合国产亚洲精品| 国内久久婷婷六月综合欲色啪| 日本黄大片高清| 久久久精品94久久精品| 人妻久久中文字幕网| 国产精品1区2区在线观看.| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 亚洲经典国产精华液单| 少妇的逼好多水| 成年免费大片在线观看| 少妇熟女aⅴ在线视频| 国产蜜桃级精品一区二区三区| 国产精华一区二区三区| 午夜视频国产福利| 国产一区二区在线观看日韩| 日本一二三区视频观看| 国产伦理片在线播放av一区 | 国产黄a三级三级三级人| 桃色一区二区三区在线观看| 久久久久国产网址| 又粗又硬又长又爽又黄的视频 | 婷婷亚洲欧美| 亚州av有码| 日韩一区二区三区影片| 伊人久久精品亚洲午夜| 国产一区亚洲一区在线观看| 亚洲成a人片在线一区二区| 成年版毛片免费区| 色综合站精品国产| 国产69精品久久久久777片| 国产白丝娇喘喷水9色精品| 欧美一区二区国产精品久久精品| 亚洲精品456在线播放app| av.在线天堂| 在线观看美女被高潮喷水网站| 成人亚洲精品av一区二区| 国产极品精品免费视频能看的| 一个人免费在线观看电影| 中文字幕久久专区| 精品人妻视频免费看| 午夜精品在线福利| 欧美不卡视频在线免费观看| 久久精品影院6| 精品日产1卡2卡| 亚洲欧美精品自产自拍| 嫩草影院新地址| 国产成人精品一,二区 | 最近2019中文字幕mv第一页| 国产毛片a区久久久久| 一边亲一边摸免费视频| 中文字幕免费在线视频6| 欧美一区二区精品小视频在线| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜添av毛片| 黑人高潮一二区| 男人狂女人下面高潮的视频| 午夜福利在线观看免费完整高清在 | 日韩,欧美,国产一区二区三区 | 欧美一区二区国产精品久久精品| 九九爱精品视频在线观看| 全区人妻精品视频| 狠狠狠狠99中文字幕| a级一级毛片免费在线观看| 最新中文字幕久久久久| 热99re8久久精品国产| 白带黄色成豆腐渣| 变态另类丝袜制服| 日本在线视频免费播放| 蜜臀久久99精品久久宅男| 国模一区二区三区四区视频| 亚洲成av人片在线播放无| 在线免费观看的www视频| 久久人人爽人人片av| 夜夜看夜夜爽夜夜摸| 嫩草影院入口| 久久久久免费精品人妻一区二区| 国产伦精品一区二区三区四那| 美女大奶头视频| 成人性生交大片免费视频hd| 全区人妻精品视频| 在线观看午夜福利视频| 成人高潮视频无遮挡免费网站| 久久久精品94久久精品| 最近最新中文字幕大全电影3| av国产免费在线观看| 国产精品美女特级片免费视频播放器| 免费观看在线日韩| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 亚洲欧洲国产日韩| 亚洲欧美精品综合久久99| 日本五十路高清| 久久久久网色| 一级黄色大片毛片| 嫩草影院新地址| 欧美日本视频| 九色成人免费人妻av| 国产av一区在线观看免费| www日本黄色视频网| 精品一区二区三区视频在线| 亚洲av二区三区四区| 亚洲av一区综合| 99久久精品一区二区三区| 国产精品一及| 女的被弄到高潮叫床怎么办| 波多野结衣高清无吗| 亚洲四区av| 永久网站在线| 六月丁香七月| 97超碰精品成人国产| 变态另类丝袜制服| 尤物成人国产欧美一区二区三区| 人人妻人人澡欧美一区二区| 日本一二三区视频观看| 免费av不卡在线播放| 欧美激情在线99| 日韩国内少妇激情av| 嫩草影院精品99| 99久久精品一区二区三区| 五月玫瑰六月丁香| 亚洲第一电影网av| eeuss影院久久| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 国产淫片久久久久久久久| 夫妻性生交免费视频一级片| 一夜夜www| 国产午夜精品论理片| 日韩欧美三级三区| 亚洲精品乱码久久久v下载方式| 一区福利在线观看| 亚洲av不卡在线观看| 91av网一区二区| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 男人和女人高潮做爰伦理| 你懂的网址亚洲精品在线观看 | 丝袜美腿在线中文| 久久久久久久久久成人| 欧美性猛交黑人性爽| 日韩欧美在线乱码| 久久久精品94久久精品| 啦啦啦观看免费观看视频高清| 国产精品,欧美在线| 精品国内亚洲2022精品成人| 国产极品精品免费视频能看的| 高清日韩中文字幕在线| 午夜精品国产一区二区电影 | 国产伦一二天堂av在线观看| 黑人高潮一二区| 极品教师在线视频| 国产精品一区二区性色av| 偷拍熟女少妇极品色| 日韩大尺度精品在线看网址| 亚洲18禁久久av| 国产日韩欧美在线精品| 色吧在线观看| 99热这里只有精品一区| 哪个播放器可以免费观看大片| 麻豆精品久久久久久蜜桃| av天堂在线播放| 人妻系列 视频| 2021天堂中文幕一二区在线观| 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 99久久精品一区二区三区| 日韩制服骚丝袜av| 蜜桃久久精品国产亚洲av| av免费观看日本| 精品日产1卡2卡| av国产免费在线观看| 午夜福利成人在线免费观看| 看非洲黑人一级黄片| 极品教师在线视频| 成人国产麻豆网| 老师上课跳d突然被开到最大视频| 国产国拍精品亚洲av在线观看| 国产午夜精品一二区理论片| 少妇裸体淫交视频免费看高清| 国产精品.久久久| 嫩草影院入口| 亚洲成人av在线免费| 麻豆国产av国片精品| 蜜臀久久99精品久久宅男| 国产v大片淫在线免费观看| 哪个播放器可以免费观看大片| 国产精品久久久久久久久免| 精品国产三级普通话版| 九九爱精品视频在线观看| 亚洲av.av天堂| 啦啦啦观看免费观看视频高清| 国产精品三级大全| 国产精品不卡视频一区二区| av在线播放精品| 成人特级黄色片久久久久久久| 免费观看的影片在线观看| 久久精品人妻少妇| 秋霞在线观看毛片| 日本色播在线视频| 色尼玛亚洲综合影院| 免费不卡的大黄色大毛片视频在线观看 | 在线a可以看的网站| 亚洲欧美清纯卡通| 99久久久亚洲精品蜜臀av| 国产美女午夜福利| 日产精品乱码卡一卡2卡三| 性欧美人与动物交配| 26uuu在线亚洲综合色| 欧美另类亚洲清纯唯美| 久久久久久久久大av| 又黄又爽又刺激的免费视频.| 有码 亚洲区| 91精品一卡2卡3卡4卡| 搡老妇女老女人老熟妇| 亚洲经典国产精华液单| 免费观看的影片在线观看| 日韩一本色道免费dvd|