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

    松嫩平原含水層沉積物特征及其對砷賦存態(tài)分布的影響

    2015-04-09 23:36:20李付蘭等
    地球科學與環(huán)境學報 2015年1期
    關(guān)鍵詞:氧化物水力沉積物

    李付蘭等

    摘 要: 松嫩平原西部地區(qū)為中國砷中毒區(qū)域,含水層沉積物特征及砷賦存態(tài)對地下水砷遷移、轉(zhuǎn)化的影響并不清楚。在對研究區(qū)高砷區(qū)與低砷區(qū)鉆孔沉積物進行取樣的基礎(chǔ)上,運用激光粒度分析、Tessier五步提取法和Keon分步提取法分別確定了不同深度沉積物的粒度分布和不同賦存態(tài)砷含量,系統(tǒng)研究了沉積物的水力傳導(dǎo)系數(shù)、不同賦存態(tài)砷的分布及其與高砷地下水形成的關(guān)系。結(jié)果表明:沉積物5種賦存態(tài)砷(包括可交換態(tài)、碳酸鹽吸附態(tài)、鐵/錳氧化物吸附態(tài)、有機物或硫化物吸附態(tài)和基質(zhì)態(tài))中,鐵/錳氧化物吸附態(tài)砷含量最高,并且砷含量與鐵、錳含量呈明顯正相關(guān)性; 高砷區(qū)鉆孔沉積物鐵/錳氧化物吸附態(tài)砷質(zhì)量分數(shù)為(1.18~6.11)×10-6,低砷區(qū)為(0.67~2.60)× 10-6;高砷區(qū)鉆孔沉積物以粉砂、黏土為主,水力傳導(dǎo)系數(shù)介于0.004~0.027 m·d-1,低砷區(qū)鉆孔沉積物以細砂、中砂為主,水力傳導(dǎo)系數(shù)介于0.005~0.2 m·d-1;沉積物中粉砂或黏土含量越多,則沉積物顆粒越小,其水力傳導(dǎo)系數(shù)越低,各賦存態(tài)砷含量越高,因此,高砷區(qū)低水力傳導(dǎo)系數(shù)有利于還原環(huán)境的形成。通過地下水的循環(huán),高砷含水層沉積物的淋濾作用和鐵/錳氧化物的還原性溶解形成了研究區(qū)高砷地下水。

    關(guān)鍵詞: 地下水;砷;沉積物;賦存態(tài);粒度;連續(xù)提??;水力傳導(dǎo)系數(shù);松嫩平原

    中圖分類號: P641.3 文獻標志碼: A

    Abstract: The western Songnen Basin is one of typical arsenic poisoning regions in China. However, the effects of sediment characteristics and arsenic occurrence modes on groundwater arsenic migration are still unknown. Based on sediment sampling from boreholes with high or low contents of arsenic, particle size distributions and contents of different arsenic occurrence modes of sediments in different depths were analyzed by laser particle size analyzer, Tessier five-step extraction method and Keon step-by-step extraction method. The roles of hydraulic conductivity coefficient and distribution of different arsenic occurrence modes of sediments in forming groundwater with high content of arsenic were investigated. The results show that content of Fe/Mn oxide absorbed arsenic, which has significantly positive correlation with contents of Fe and Mn, is the highest among the five arsenic occurrence modes including exchangeable, carbonate absorbed, Fe/Mn oxide absorbed, organic or sulfide absorbed and metrical modes; mass fractions of Fe/Mn oxide absorbed arsenics are (1.18-6.11)×10-6 and (0.67-2.60)×10-6 for the sediments from boreholes in the area with high or low contents of arsenic; the sediments from boreholes in the area with high content of arsenic are mainly silt and clay with the hydraulic conductivity coefficient of 0.004-0.027 m·d-1, and the sediments in the area with low content of arsenic are fine and medium sands with the hydraulic conductivity coefficient of 0.005-0.2 m·d-1; sediments with higher contents of silt or clay have lower hydraulic conductivity coefficient and higher contents of different arsenic occurrence modes, therefore, low hydraulic conductivity coefficient is in favour of the formation of reducing environment in the area with high content of arsenic. Because of the groundwater circulation, leaching of sediments and reductive dissolution of Fe/Mn oxides in aquifer with high content of arsenic contribute to the formation of groundwater with high content of arsenic.

    Key words: groundwater; arsenic; sediment; occurrence mode; particle size; sequential extraction; hydraulic conductivity coefficient; Songnen Plain

    0 引 言

    砷是自然界中一種有毒元素,長期飲用高砷地下水會導(dǎo)致人體出現(xiàn)皮膚色素異常、角質(zhì)化、皮膚癌、內(nèi)臟癌癥等慢性砷中毒。 當?shù)叵滤械纳橘|(zhì)量濃度超過世界衛(wèi)生組織(WHO)規(guī)定的飲用水標準10 μg·L-1時, 便可認為是高砷地下水。巖石和沉積物中的砷在生物地球化學和水文地質(zhì)作用下會進入地下水中,因此,加強沉積物及地下水中砷分布的研究已成為科學界非常重要的任務(wù)之一[1-2]。McArthur等認為,沉積盆地地下水中砷分布的差異由含水層上覆地層的滲透性以及上覆沉積物中有機質(zhì)含量決定[3]。Saha等發(fā)現(xiàn),印度東部甘加平原淺層高砷含水層系統(tǒng)中,上覆沉積物巖性以黏土、粉砂為主,且?guī)r層的水力傳導(dǎo)系數(shù)較低[4]。這說明地層沉積物的水力傳導(dǎo)系數(shù)對地下水中砷質(zhì)量濃度有一定的影響。此外,Guo等發(fā)現(xiàn),高砷地下水的形成與還原環(huán)境密切相關(guān),地下水氧化-還原條件的變化也伴隨著砷質(zhì)量濃度的變化[5-6]。鐵的還原性釋放機制被普遍認為是地下水砷富集的重要途徑[7-9]。

    松嫩平原西部地區(qū)是砷中毒區(qū),砷質(zhì)量濃度平均為39.3 μg·L-1,在一些地區(qū)砷質(zhì)量濃度為100~150 μg·L-1 [10]。在垂直剖面上,砷主要富集在4~20 m的潛水含水層和30~50 m的承壓含水層中[11]。其中,所取水樣中,65.5%的淺層地下水和96%的深層地下水砷質(zhì)量濃度超過10 μg·L-1 [12]。高砷地下水的分布不僅與沉積物水力學特征有關(guān),還受其賦存態(tài)砷分布的影響,因此,沉積物粒度分布及其砷的賦存態(tài)對當?shù)馗呱榈叵滤纬捎兄陵P(guān)重要的影響。盡管李月芬等對吉林西部地區(qū)土壤中砷的賦存態(tài)進行了研究[13],但并沒有與高砷地下水的形成結(jié)合起來。

    本文就沉積物中砷的賦存態(tài)及其水力傳導(dǎo)系數(shù)之間的聯(lián)系開展研究。在查清地下水中砷質(zhì)量濃度分布規(guī)律的基礎(chǔ)上,分別選擇高砷區(qū)和低砷區(qū)為鉆孔取樣區(qū),通過采集不同深度沉積物樣品,探討了沉積物中不同賦存態(tài)砷分布特征,并研究了高砷區(qū)鉆孔和低砷區(qū)鉆孔沉積物的粒度分布,討論了含水層水力傳導(dǎo)系數(shù)對不同賦存態(tài)砷分布的影響,揭示了松嫩平原西部地區(qū)地下水循環(huán)對含水層高砷地下水形成的影響。

    1 研究區(qū)概況

    研究區(qū)位于吉林省西北部,經(jīng)緯度為44°13′5″N~ 45°16′27″N,122°2′13″E~123°30′37″E,南北長約116 km,東西寬約112 km,面積8 468 km2(圖1)。多年平均降水量為402 mm,其區(qū)域分布呈從東向西略微退減趨勢。

    本區(qū)域地下水形成分布規(guī)律主要受氣候、水文、巖石地層結(jié)構(gòu)、地質(zhì)構(gòu)造及地貌等多種因素的綜合影響。該區(qū)域大部分地區(qū)埋藏有多個含水層,包括第四系孔隙淺層地下水和承壓水含水層,新近系大安組、泰康組孔隙裂隙含水層和白堊系下統(tǒng)及上統(tǒng)裂隙孔隙含水層。在天然情況下,各含水層之間在平面或剖面上存在直接或間接的水力聯(lián)系,其共同邊界由盆地周邊白堊系之前的各種弱滲透性地層和阻水斷層組成,構(gòu)成一個較完整的含水層系統(tǒng)。其中,第四系孔隙含水層系統(tǒng)在區(qū)內(nèi)分布范圍最廣,含水層為第四系各種成因的松散沉積物,以粒間孔隙為儲水空間。新近系含水層巖石呈微弱膠結(jié),以孔隙-裂隙為儲水空間,基本上成層分布,普遍具承壓性,屬于孔隙-裂隙承壓水含水層。地下水的補給來源以降水入滲占主導(dǎo)地位,以潛水蒸發(fā)和人工開采為主要排泄方式。由于平原區(qū)地形平緩,區(qū)域地下水力坡度小,地下水徑流較為緩慢,所以地下水運動以垂直運動為主,水平運動占次要地位,主要是由盆地四周向中心流動[14],由東部高平原和西部山前傾斜平原向中部水平運動。地下水類型主要為Na-Mg/Ca-HCO3。研究區(qū)的水文地質(zhì)單元可分為3個區(qū),具體可見圖1。

    2 樣品采集與測試

    2.1 樣品采集

    地下水采樣分析發(fā)現(xiàn),高砷地下水具有明顯的空間分布特征[12,14-18],采樣深度范圍為13~120 m,其中淺層地下水為13~32 m。Guo等通過分析所采集的地下水樣品發(fā)現(xiàn),在平面上,從向海水庫到吉林省通榆縣四井子鄉(xiāng),地下水中的砷質(zhì)量濃度逐漸增高,淺層地下水中砷質(zhì)量濃度超過100 μg·L-1的點全部分布在四井子鄉(xiāng)[12];在垂向上,淺層高砷地下水主要分布在10~25 m。圖1顯示四井子鄉(xiāng)為研究區(qū)淺層高砷地下水分布區(qū),向海水庫附近為淺層低砷地下水分布區(qū)。

    根據(jù)圖1中地下水砷的分布以及相關(guān)病情資料[19],利用手動采樣鉆(荷蘭Eijkelkam公司),分別于高砷區(qū)和低砷區(qū)進行鉆探,并采集沉積物樣品。高砷區(qū)鉆孔(ZK1)經(jīng)緯度為44°54′0.28″N,122°48′25.4″E, 位于四井子鄉(xiāng);低砷區(qū)鉆孔(ZK2)經(jīng)緯度為44°58′10.7″N,122°33′36.02″E, 位于興隆漁場(圖1)。 2個鉆孔均位于中部平原潛水承壓水區(qū),地下水水位埋深較淺(淺于8 m),鉆孔深度小于10 m, 共采集沉積物樣品61件。在采樣時,先刮去表層的沉積物,避免各層位沉積物的相互污染。樣品采集后,裝入封口袋中密閉保存。沉積物的采集深度、巖性特征見圖2。

    2.2 沉積物分析測試

    2.2.1 沉積物粒度分布測定

    將沉積物樣品冷凍干燥后,過孔徑為1 mm篩子,后取0.2~0.5 g樣品(黏土少取一些,砂多取一點),放入100 mL燒杯中;向燒杯中加入30%雙氧水約20 mL,在熱板上加熱(100 ℃),消解至無泡沫;待泡沫消失后,加幾滴10%鹽酸,微沸0.5 h,當泡沫消失后取懸浮液,利用激光粒度分析儀(型號為Mastersizer 2000)進行粒度分布測定。

    2.2.2 沉積物的全量分析

    采用X射線熒光光譜分析儀(PANalytical公司Epsilon5型)測定沉積物樣品中SiO2、Al2O3、Fe2O3、MgO、CaO、Na2O、K2O等含量(質(zhì)量分數(shù),下同),測試精度為5%。全消解-ICP-MS標準方法(Thermal公司IRISIntrepid Ⅱ型)測定沉積物中的As、Ba、Cd、Cr、Cu、F、Li、Mn、Mo、Ni、P、Pb、Sb、S、Se、Sr、U、V、Ti、Br、B、Zn、Si等含量,測試精度為2%,其中As含量檢測限為0.2×10-6。

    2.2.3 沉積物中不同賦存態(tài)砷的提取

    根據(jù)深度和巖性變化的特點,選取30件具有代表性的沉積物樣品。在分步提取試驗之前,對樣品進行冷凍干燥并研磨至粒徑小于125 μm。采用Tessier五步提取法[20]分步提取沉積物樣品中不同賦存態(tài)砷的含量。按順序依次采用不同的提取液,以一定的固液比均勻混合后,分別提取不同賦存態(tài)砷、鐵、錳,包括可交換態(tài)、碳酸鹽吸附態(tài)、鐵/錳氧化物吸附態(tài)、有機物或硫化物吸附態(tài)和基質(zhì)態(tài)。具體提取步驟如下:

    (1)在室溫下稱取2.50 g沉積物樣品,用20 mL 0.50 mol·L-1的硝酸鎂溶液萃取30 min,用1∶1(體積比)硝酸調(diào)節(jié)萃取體系pH值至7.0,攪拌5 min(不斷攪拌)得到可交換態(tài)。

    (2)取經(jīng)步驟(1)萃取的離心殘渣,用20 mL 1.0 mol·L-1的醋酸鈉溶液萃取5 h(應(yīng)不間斷攪拌);在通風廚下用乙酸調(diào)節(jié)酸堿度,使萃取液pH值為5.0,攪拌5 min得到碳酸鹽吸附態(tài)。

    (3)取經(jīng)步驟(2)后的離心殘渣,用恒溫加熱器使溫度保持在96 ℃,用50 mL 0.08 mol·L-1鹽酸羥胺的乙酸溶液(體積比為1∶4)提取6 h,攪拌5 min得到鐵/錳氧化物吸附態(tài)。

    (4)取經(jīng)步驟(3)后的離心殘渣,在85 ℃下用7.5 mL 0.02 mol·L-1的硝酸萃取2 h;再加入12.5 mL 300 mol·L-1的過氧化氫(用硝酸調(diào)節(jié)pH值至2.0),在85 ℃下繼續(xù)萃取3 h,冷卻;加入12.5 mL 3.2 mol·L-1醋酸銨的硝酸溶液(體積比為1∶50),用蒸餾水稀釋到50 mL,持續(xù)攪拌30 min得到有機物或硫化物吸附態(tài)。

    (5)取經(jīng)步驟(4)后的離心殘渣,在105 ℃用50 mL濃硝酸提取,直到只剩12.5 mL溶液,得到基質(zhì)態(tài)。

    將每一步的提取液定容后裝瓶,待用電感耦合等離子體光譜法(ICP-AES)測定提取液中砷、鐵、錳含量,然后計算出不同賦存態(tài)砷、鐵、錳的含量。

    此外,利用Keon分步提取法[21-22]對全部沉積物樣品分析提取弱吸附態(tài)砷與強吸附態(tài)砷。具體步驟如下:

    (1)在室溫下,將0.4 g沉積物放入50 mL離心管中,加入20 mL 1.0 mol·L-1 氯化鎂溶液(調(diào)節(jié)pH值至8.0),萃取2 h,不間斷攪拌,重復(fù)2次,將2次的提取液裝入同一瓶中,得到弱吸附態(tài)砷。

    (2) 在室溫下,取經(jīng)步驟(1)萃取的離心殘渣, 加入20 mL 1.0 mol·L-1 磷酸二氫鈉溶液(調(diào)節(jié)pH值至5.0),重復(fù)2次,2次萃取時間分別為16 h和24 h,將2次的提取液轉(zhuǎn)入同一瓶中,得到強吸附態(tài)砷。

    沉積物提取試驗在厭氧操作箱中進行,所有提取液都用超純水制備,超純水要用氮氣吹2 h脫氧。

    3 結(jié)果分析

    3.1 含水層的水力學特征

    沉積物粒度分布見圖3。利用Beyer提出的公式計算相應(yīng)沉積物的水力傳導(dǎo)系數(shù)[23]。Beyer提出的公式為

    從圖3可以看出:鉆孔ZK1中粉砂或黏土(粒徑小于0.006 3 mm)體積分數(shù)平均為54.8%,細砂(粒徑介于0.063~0.2 mm)與粉砂或黏土所占體積分數(shù)平均為95.5%;鉆孔ZK2中粉砂或黏土所占體積分數(shù)平均為31.9%,細砂與中砂(粒徑大于0.2 mm)所占體積分數(shù)平均為68.1%。整體來看,鉆孔ZK1沉積物的粒徑小于鉆孔ZK2沉積物。鉆孔ZK1中,在2~5 m深度的粉砂或黏土與細砂2個粒級體積分數(shù)之和基本為100%,在5~9.25 m深度的細砂和中砂體積分數(shù)增加,粉砂或黏土體積分數(shù)有所減少。鉆孔ZK2中,在1~3.5 m深度的粉砂或黏土和中砂粒級體積分數(shù)較高,平均為89.1%, 在5~5.5 m 深度夾有一層薄的黏土層使其粒度分布在此處有突變。

    從圖4(a)可以看出:鉆孔ZK1沉積物的水力傳導(dǎo)系數(shù)為0.004~0.027 m·d-1;鉆孔ZK2沉積物的水力傳導(dǎo)系數(shù)為0.005~0.2 m·d-1,整體較鉆孔ZK1高。對比圖3、4可知,隨著埋深的增加,顆粒越小,沉積物的水力傳導(dǎo)系數(shù)也越小。

    3.2 沉積物中砷含量

    沉積物總砷的分布見圖5。從圖5可以發(fā)現(xiàn),高砷區(qū)鉆孔ZK1沉積物中砷含量明顯大于低砷區(qū)鉆孔ZK2沉積物。其中,鉆孔ZK1沉積物中砷含量(質(zhì)量分數(shù),下同)介于(2.23~102)×10-6,鉆孔ZK2介于(1.48~7.58)×10-6。在鉆孔ZK1中2.50 m深度處,樣品ZK1-12中砷含量達到最大值,為褐色黏土質(zhì)粉細砂。在鉆孔ZK2中3.5 m深度處,樣品ZK2-13中砷含量達到最大值,為灰黑色粉細砂。在垂直剖面上,鉆孔ZK1沉積物中砷含量分布不均勻,總體呈先增大后減小的趨勢,鉆孔ZK2沉積物中砷含量整體分布較均勻,在整個剖面上都較小。另外,Tessier五步提取法中5種不同賦存形態(tài)砷的總含量與總砷基本吻合。

    3.3 沉積物中不同賦存態(tài)砷

    圖6顯示2個鉆孔沉積物中不同賦存態(tài)砷含量的分布?;|(zhì)態(tài)砷和鐵/錳氧化物吸附態(tài)砷是主要賦存態(tài),而可交換態(tài)砷、碳酸鹽吸附態(tài)砷、有機物或硫化物吸附態(tài)砷所占比例較小。總的來說,各賦存態(tài)砷含量從大到小依次為:基質(zhì)態(tài)砷、鐵/錳氧化物吸附態(tài)砷、碳酸鹽吸附態(tài)砷、有機物或硫化物吸附態(tài)砷、可交換態(tài)砷。鉆孔ZK1沉積物中各賦存態(tài)含量均大于鉆孔ZK2,因此,可推斷沉積物中砷含量的差異可能是地下水高砷分布的原因之一。

    可交換態(tài)砷是指通過擴散作用和外層絡(luò)合作用非專屬性地吸附在土、腐殖質(zhì)及其他成分上的砷,對環(huán)境變化敏感, 易于遷移轉(zhuǎn)化,能被植物吸收。通過離子交換可將可交換態(tài)砷從沉積物樣品中釋放出來。由圖6可見,2個鉆孔沉積物中可交換態(tài)砷含量均較低。 鉆孔ZK1沉積物中可交換態(tài)砷含量介于(0.058~0.22)×10-6,鉆孔ZK2沉積物中可交換態(tài)砷含量低于0.08× 10-6,由此可知,鉆孔ZK1沉積物中可交換態(tài)砷含量大于鉆孔ZK2。

    碳酸鹽吸附態(tài)砷是指吸附在碳酸鹽礦物上的砷,對環(huán)境條件(特別是pH值)較敏感。當pH值下降時,碳酸鹽吸附態(tài)砷易釋放出來并進入環(huán)境;當pH值升高時,其與可交換態(tài)砷相似。2個鉆孔沉積物中碳酸鹽吸附態(tài)砷含量較低。鉆孔ZK1沉積物中碳酸鹽吸附態(tài)砷含量介于 (0.058~0.22)×10-6;鉆孔ZK2沉積物中碳酸鹽吸附態(tài)砷含量介于(0.2~ 0.6)×10-6,個別黏土樣品高于1.0×10-6。

    對于鐵/錳氧化物吸附態(tài)砷,一般情況下,鐵/錳氧化物比表面積大,吸附或共沉淀砷的能力較強。2個鉆孔沉積物中,鐵/錳氧化物吸附態(tài)砷所占平均比例分別為21.2%和39.5%(圖6)。鉆孔ZK1中鐵/錳氧化物吸附態(tài)砷含量最高可達6.0×10-6;鉆孔ZK2中最高為1.8×10-6。這表明鐵/錳氧化物對砷有較強的結(jié)合能力。

    有機物或硫化物吸附態(tài)砷以不同形式進入或包裹在有機質(zhì)顆粒上,與有機質(zhì)螯合或以硫化物形式存在。在氧化條件下,有機質(zhì)可以被降解,并釋放砷。在鉆孔ZK1中,除了個別黏土樣品的有機物或硫化物吸附態(tài)砷含量較高(大于1.0×10-6),其他樣品平均含量介于(0.2~0.6)×10-6。 鉆孔ZK2中有機物或硫化物吸附態(tài)砷分布較為均勻,基本上介于(0.2~0.4)×10-6, 部分樣品大于鉆孔ZK1,部分樣品與鉆孔ZK1基本在一個數(shù)量級范圍之內(nèi)。

    基質(zhì)態(tài)砷是原生礦物中砷最主要的組成部分,一般存在于硅酸鹽、原生礦物和次生礦物等的晶格中,在自然條件下不易釋放,能長期穩(wěn)定在沉積物中,生物可利用性較差。2個鉆孔沉積物中基質(zhì)態(tài)砷平均含量分別為8.85×10-6和1.30×10-6,所占平均比例分別為69.0%和40.3%。

    弱吸附態(tài)砷是指被土壤膠體吸附或與鐵、鋁等結(jié)合的吸附態(tài)砷。從圖7可以看出,鉆孔ZK1、ZK2沉積物中弱吸附態(tài)砷含量均較低。鉆孔ZK1中弱吸附態(tài)砷含量介于(0.15~1.70)×10-6,鉆孔ZK2中低于0.35×10-6。總體來說,鉆孔ZK1弱吸附態(tài)砷含量高于鉆孔ZK2。

    強吸附態(tài)砷主要是以陰離子形式與土壤中帶正電荷的質(zhì)點相互作用結(jié)合。從圖7可以看出:鉆孔ZK1沉積物中強吸附態(tài)砷分布不均勻,含量介于(0.70~25.5)×10-6,在整個剖面1.7~5.5 m深度處較高;在鉆孔ZK2沉積物中強吸附態(tài)砷分布較為均勻,含量最高為0.80×10-6,整體顯著低于鉆孔ZK1沉積物。

    從圖6、7可以看出,鉆孔ZK1沉積物中可交換態(tài)砷含量介于(0.058~0.22)×10-6,明顯小于弱吸附態(tài)砷與強吸附態(tài)砷含量總和。交換態(tài)砷含量以及弱吸附態(tài)砷與強吸附態(tài)砷的含量總和在垂直剖面上的變化趨勢基本一致,說明Keon分步提取法提取可交換態(tài)砷的效果較為顯著。從圖4可以看出,強、弱吸附態(tài)砷含量在水力傳導(dǎo)系數(shù)較小時相對較大,在1.8~3.5 m深度處強、弱吸附態(tài)砷含量在2個鉆孔中相對較高,而水力傳導(dǎo)系數(shù)則相對較低,因此,水力傳導(dǎo)系數(shù)越小,強、弱吸附態(tài)砷含量越高。

    Tessier五步提取法得到的不同賦存態(tài)砷含量與鐵、錳含量存在一定的相關(guān)關(guān)系??山粨Q態(tài)砷、碳酸鹽吸附態(tài)砷和基質(zhì)態(tài)砷含量與鐵、錳含量的相關(guān)性較差(判定系數(shù)小于0.1),而鐵/錳氧化物吸附態(tài)砷和有機物或硫化物吸附態(tài)砷含量與鐵、錳含量成一定的正相關(guān)關(guān)系(判定系數(shù)大于0.52)。李月芬等認為吉林省通榆縣表層土壤中不同賦存態(tài)砷含量從多到少依次為:殘渣態(tài)、腐植酸結(jié)合態(tài)、鐵/錳氧化物吸附態(tài)、碳酸鹽吸附態(tài)、水溶態(tài)、離子交換態(tài)、強有機結(jié)合態(tài)[13]。Guo等提出鐵/錳氧化物礦物在適當?shù)倪€原條件下可以釋放砷進入地下水中[9]。根據(jù)野外數(shù)據(jù)可知,研究區(qū)富砷地下水氧化-還原電位小于0 mV,地下水處于還原環(huán)境,有利于含水層中砷的釋放,因此,在還原環(huán)境中,鐵/錳氧化物礦物的還原性溶解是地下水中砷富集的一個原因。

    4 討 論

    4.1 沉積物的粒徑分布對砷賦存態(tài)的影響

    沉積物粒徑越小,所含強、弱吸附態(tài)砷含量之和越高[圖8(a)、(b)],鐵/錳氧化物吸附態(tài)砷含量也越高[圖8(c)、(d)]。由圖3(a)可知,鉆孔ZK1中的沉積物顆粒粒級較鉆孔ZK2小,基本沒有中砂。因此,在鉆孔ZK1整個剖面上,鐵/錳氧化物吸附態(tài)砷含量遠遠大于鉆孔ZK2。沉積物中弱吸附態(tài)砷、強吸附態(tài)砷、鐵/錳氧化物吸附態(tài)砷含量在粒徑較小的沉積物中較大。一般而言,沉積物顆粒越小,所含黏土礦物、鐵/錳氧化物礦物含量越高。這些礦物能夠有效吸附砷,一方面使強、弱吸附態(tài)砷含量增加, 另一方面使沉積物中鐵/錳氧化物吸附態(tài)砷等含量升高。Kumar等在印度恒河平原發(fā)現(xiàn),埋深為4.57~ 7.62 m之間的沉積物主要為黏土和粉砂,且含有較多的鐵/錳氧化物礦物和黑云母物礦物,沉積物中砷含量也較高[24]。

    4.2 水動力條件對沉積物中砷賦存態(tài)的影響

    由圖9可以看出:在水力傳導(dǎo)系數(shù)較小時,鐵/錳氧化物吸附態(tài)砷含量,強、弱吸附態(tài)砷含量之和及各賦存態(tài)砷總含量都明顯升高;當水力傳導(dǎo)系數(shù)較大時,各吸附態(tài)砷含量明顯降低。由于鉆孔ZK1中水力傳導(dǎo)系數(shù)整體較小,且分布均勻,所以鉆孔ZK1沉積物中鐵/錳氧化物吸附態(tài)砷含量以及強、弱吸附態(tài)砷含量之和都明顯高于鉆孔ZK2。總體而言,水力傳導(dǎo)系數(shù)越高,各吸附態(tài)砷含量越低。 在對孟加拉國的研究中發(fā)現(xiàn),水動力條件越好,地下水對含水層的沖刷作用越強烈[23-24],含水層沉積物中吸附態(tài)砷含量越低。當含水層水力傳導(dǎo)系數(shù)較低時,沖刷作用較弱,沉積物中不同吸附態(tài)砷含量就越高。

    根據(jù)野外所測數(shù)據(jù)可知,四井子鄉(xiāng)附近地下水水位埋深較淺(淺于3 m)。該地區(qū)采用地下水灌溉,地下水水位波動加速了對沉積物的淋濾作用[25-28],而在上覆2~5.5 m深度的沉積物中,強、弱吸附態(tài)砷含量最高。因此,通過充分的水-巖相互作用,沉積物中的強、弱吸附態(tài)砷易釋放進入地下水,這是地下水中富集砷的一個重要途徑。此外,Guo等提出,補給區(qū)地表水和大氣降水快速的補給可以為含水層補充氧氣[29-30]。鉆孔ZK1沉積物的水力傳導(dǎo)系數(shù)整體較鉆孔ZK2要小得多。鉆孔ZK1含水層的水動力條件較差,導(dǎo)致在鉆孔ZK1附近地下水循環(huán)受阻,不利于氧氣進入含水層,使得含水層中形成一個相對還原的環(huán)境。在野外采集地下水樣品時,所測地下水氧化-還原電位低于0 mV也可證實這一點。此外,鉆孔ZK1中鐵/錳氧化物吸附態(tài)砷含量較鉆孔ZK2高,且為沉積物中砷的主要賦存態(tài)。這些砷的賦存態(tài)特征為沉積物中砷的釋放提供了物質(zhì)基礎(chǔ)和環(huán)境條件。在還原環(huán)境中,這些鐵/錳氧化物吸附態(tài)砷易釋放出來,進入地下水。因此,沉積物中鐵/錳氧化物的還原性溶解是地下水中砷來源的一個重要途徑,對四井子鄉(xiāng)淺層地下水中砷的富集貢獻較大。

    5 結(jié) 語

    (1)松嫩平原西部地區(qū)沉積物中砷賦存態(tài)由多到少分別為基質(zhì)態(tài)、鐵/錳氧化物吸附態(tài)、碳酸鹽吸附態(tài)、有機物或硫化物吸附態(tài)、可交換態(tài)。鐵/錳氧化物吸附態(tài)砷與基質(zhì)態(tài)砷含量之和約占總砷80%,鐵/錳氧化物吸附態(tài)是除基質(zhì)態(tài)外的主要賦存態(tài)。高砷區(qū)沉積物的各賦存態(tài)砷含量整體高于低砷區(qū)。

    (2)沉積物中砷含量與鐵、錳含量呈顯著正相關(guān)關(guān)系。此外,鐵/錳氧化物吸附態(tài)砷含量與鐵、錳含量呈顯著正相關(guān)性,說明沉積物中的砷主要與鐵/錳氧化物結(jié)合。

    (3)高砷區(qū)沉積物的顆粒粒徑和水力傳導(dǎo)系數(shù)整體較低砷區(qū)小。沉積物中粉砂或黏土含量越大,沉積物顆粒越小,其水力傳導(dǎo)系數(shù)越低,各賦存態(tài)砷含量越高。在高砷區(qū)鉆孔中,埋深為2.5 m處沉積物中粉砂和黏土體積分數(shù)占86%,其賦存態(tài)砷總含量達到最高(63.1×10-6),水力傳導(dǎo)系數(shù)較低(0.006 03 m·d-1)。

    (4)沉積物粒度分布和含水層的水動力條件通過影響沉積物中砷的賦存態(tài)分布,控制高砷地下水的形成。在低水力傳導(dǎo)系數(shù)、高砷含量的含水層中,砷主要通過2個途徑富集于地下水中:一是沉積物中強、弱吸附態(tài)砷在充分的淋濾作用下,通過水-巖相互作用進入到地下水中;二是高砷區(qū)水動力條件較差,地下水與地表水的循環(huán)被阻隔,含水層中形成了一個相對還原的環(huán)境,使得鐵/錳氧化物礦物發(fā)生還原性溶解。

    參考文獻:

    References:

    [1] 賈永鋒,郭華明.高砷地下水研究的熱點及發(fā)展趨勢[J].地球科學進展,2013,28(1):51-61.

    JIA Yong-feng,GUO Hua-ming.Hot Topics and Trends in the Study of High Arsenic Groundwater[J].Advances in Earth Science,2013,28(1):51-61.

    [2] 郭華明,郭 琦,賈永鋒,等.中國不同區(qū)域高砷地下水化學特征及形成過程[J].地球科學與環(huán)境學報,2013,35(3):83-96.

    GUO Hua-ming,GUO Qi,JIA Yong-feng,et al.Chemical Characteristics and Geochemical Processes of High Arsenic Groundwater in Different Regions of China[J].Journal of Earth Sciences and Environment,2013,35(3):83-96.

    [3] MCARTHUR J M,BANERJEE D M,HUDSON-EDWARDS K A,et al.Natural Organic Matter in Sedimentary Basins and Its Relation to Arsenic in Anoxic Ground Water:The Example of West Bengal and Its Worldwide Implications[J].Applied Geochemistry,2004,19(8):1255-1293.

    [4] SAHA D,SAHU S,CHANDRA P C.Arsenic-safe Alternate Aquifers and their Hydraulic Characteristics in Contaminated Areas of Middle Ganga Plain,Eastern India[J].Environmental Monitoring and Assessment, 2011,175(1/2/3/4):331-348.

    [5] GUO H M,ZHANG B,WANG G C,et al.Geochemical Controls on Arsenic and Rare Earth Elements Approximately Along a Groundwater Flow Path in the Shallow Aquifer of the Hetao Basin,Inner Mongolia[J].Chemical Geology,2010,270(1/2/3/4):117-125.

    [6] 高存榮,劉文波,劉 濱,等.河套平原第四紀沉積物中砷的賦存形態(tài)分析[J].中國地質(zhì),2010,37(3):760-770.

    GAO Cun-rong,LIU Wen-bo,LIU Bin,et al.Modes of Occurrence of Arsenic in Quaternary Sediments of the Hetao Plain[J].Geology in China,2010,37(3):760-770.

    [7] NICKSON R,MCARTHUR J,BURGESS W,et al.Arsenic Poisoning of Bangladesh Groundwater[J].Nature,1998,395:338.

    [8] MCARTHUR J M,RAVENSCROFT P,SAFIULLA S,et al.Arsenic in Groundwater:Testing Pollution Mechanisms for Sedimentary Aquifers in Bangladesh[J].Water Resources Research,2001,37(1):109-117.

    [9] GUO H M,YANG S Z,TANG X H,et al.Groundwater Geochemistry and Its Implications for Arsenic Mobilization in Shallow Aquifers of the Hetao Basin,Inner Mongolia[J].Science of the Total Environment, 2008,393(1):131-144.

    [10] TANG J,BIAN J M,LI Z Y,et al.A Study on Arsenic Distribution Characteristics in Groundwater of Western Jilin Province,P.R.China[C]∥IEEE.The 4th International Conference on Bioinformatics and Biomedical Engineering(iCBBE).Chengdu:IEEE,2010:90-95.

    [11] BIAN J M,TANG J,ZHANG L S,et al.Arsenic Distribution and Geological Factors in the Western Jilin Province,China[J].Journal of Geochemical Exploration, 2012,112:347-356.

    [12] GUO H M,ZHANG D,WEN D G,et al.Arsenic Mobilization in Aquifers of the Southwest Songnen Basin, P.R.China:Evidences from Chemical and Istopic Characteristics[J].Science of the Total Environment,2014,490:590-602.

    [13] 李月芬,王冬艷,湯 潔,等.吉林西部土壤砷的形態(tài)分布及其與土壤性質(zhì)的關(guān)系研究[J].農(nóng)業(yè)環(huán)境科學學報,2012,31(3):516-522.

    LI Yue-fen,WANG Dong-yan,TANG Jie,et al.Speciation of Soil Arsenic and Its Correlation with Soil Properties in Western Jilin Province,China[J].Journal of Agro-environment Science,2012,31(3):516-522.

    [14] 邴智武.松嫩平原地下水氟、砷的富集規(guī)律及影響因素研究[D].長春:吉林大學,2009.

    BING Zhi-wu.Study on the Enrichment Laws and Influencing Factors of Arsenic and Fluoride in Groundwater in the Songnen Plain[D].Changchun:Jilin University,2009.

    [15] LI Z Y,YANG Z M,BIAN J M.Safety Assessment on Drinking Water Based on Database System in Endemic Area in West of Jilin Province[C]∥IEEE.The 4th International Conference on Bioinformatics and Biomedical Engineering(iCBBE).Chengdu:IEEE,2010:105-109.

    [16] 卞建民,湯 潔,封 靈,等.吉林西部砷中毒區(qū)水文地球化學特征[J].水文地質(zhì)工程地質(zhì),2009,36(5):80-83.

    BIAN Jian-min,TANG Jie,F(xiàn)ENG Ling,et al.Hydrogeochemical Characteristics in the Arsenic Poisoning Area in Western Jilin Province[J].Hydrogeology and Engineering Geology,2009,36(5):80-83.

    [17] 封 靈.吉林西部砷中毒病區(qū)地下水砷的賦存及富集規(guī)律研究[D].長春:吉林大學,2007.

    FENG ling.Research on the Regulation of Arsenic Storage and Enrichment in the Groundwater in the Arsenic Poisoning Area in Western Jilin Province[D].Changchun:Jilin University,2007.

    [18] 卞建民,查恩爽,湯 潔,等.吉林西部砷中毒區(qū)高砷地下水反向地球化學模擬[J].吉林大學學報:地球科學版,2010,40(5):1098-1103.

    BIAN Jian-min,ZHA En-shuang,TANG Jie,et al.Inverse Geochemical Modeling of Arsenic Groundwater at Arseniasis Area in the Western of Jilin Province[J].Journal of Jilin University:Earth Science Edition, 2010,40(5):1098-1103.

    [19] 佟建冬,盧振明,張秀麗,等.吉林省飲水高砷區(qū)成人砷中毒病情調(diào)查[J].中國地方病學雜志,2006,25(5):529-531.

    TONG Jian-dong,LU Zhen-ming,ZHANG Xiu-li,et al. Investigation on Arsenic Poisoning in High Arsenic Areas with Drinking Water in Jilin Province[J].Chinese Jouranl of Endemiology,2006,25(5):529-531.

    [20] TESSIER A,CAMPBELL P G C,BISSON M.Sequential Extraction Procedure for the Speciation of Particulate Trace Metals[J].Analytical Chemistry,1979,51(7):844-851.

    [21] KEON N E,SWARTZ C H,BRABANDER D J,et al. Validation of an Arsenic Sequential Extraction Method for Evaluating Mobility in Sediments[J].Environmental Science and Technology,2001,35(13):2778-2784.

    [22] [JP2]VAN GEEN A,ZHENG Y,GOODBRED S,et al.Flushing History as a Hydrogeological Control on the Regional Distribution of Arsenic in Shallow Groundwater of the Bengal Basin[J].Environmental Science and Technology,2008,42(7):2283-2288.

    [23] STUTE M,ZHENG Y,SCHLOSSER P,et al.Hydrological Control of As Concentrations in Bangladesh Groundwater[J].Water Resources Research,2007,43(9):W09417.

    [24] KUMAR P,KUMAR M,RAMANATHAN A L,et al. Tracing the Factors Responsible for Arsenic Enrichment in Groundwater of the Middle Gangetic Plain,India:A Source Identification Perspective[J].Environmental Geochemistry and Health,2010,32(2):129-146.

    [25] WEINMAN B,GOODBRED J S L,ZHENG Y,et al. Contributions of Floodplain Stratigraphy and Evolution to the Spatial Patterns of Groundwater Arsenic in Araihazar,Bangladesh[J]. Geological Society of America Bulletin,2008,120(11/12):1567-1580.

    [26] NEIDHARDT H,BERNER Z,F(xiàn)REIKOWSKI D,et al. Influences of Groundwater Extraction on the Distribution of Dissolved As in Shallow Aquifers of West Bengal,India[J].Journal of Hazardous Materials,2013,262:941-950.

    [27] FENDORF S,MICHAEL H A,VAN GEEN A.Spatial and Temporal Variations of Groundwater Arsenic in South and Southeast Asia[J].Science,2010,328:1123-1127.

    [28] NEUMANN R B,ASHFAQUE K N,BADRUZZAMAN A B M,et al.Anthropogenic Influences on Groundwater Arsenic Concentrations in Bangladesh[J].Nature Geoscience,2010,3(1):46-52.

    [29] GUO H M,ZHANG B,LI Y,et al.Hydrogeological and Biogeochemical Constrains of Arsenic Mobilization in Shallow Aquifers from the Hetao Basin,Inner Mongolia[J].Environmental Pollution,2011,159(4):876-883.

    [30] HARVEY C F,ASHFAQUE K N,YU W,et al. Groundwater Dynamics and Arsenic Contamination in Bangladesh[J].Chemical Geology,2006,228(1/2/3):112-136.

    猜你喜歡
    氧化物水力沉積物
    水力全開
    晚更新世以來南黃海陸架沉積物源分析
    海洋通報(2022年2期)2022-06-30 06:07:04
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對磷的吸收及釋放研究進展
    相轉(zhuǎn)化法在固體氧化物燃料電池中的應(yīng)用
    陶瓷學報(2020年6期)2021-01-26 00:37:56
    細說『碳和碳的氧化物』
    氧化物的分類及其中的“不一定”
    球墨鑄鐵管的水力計算
    水力噴射壓裂中環(huán)空水力封隔全尺寸實驗
    討論用ICP-AES測定土壤和沉積物時鈦對鈷的干擾
    老司机福利观看| 女人精品久久久久毛片| 桃红色精品国产亚洲av| 一夜夜www| а√天堂www在线а√下载| 中文字幕人妻丝袜制服| 日韩av在线大香蕉| 十八禁人妻一区二区| 男人的好看免费观看在线视频 | 青草久久国产| 久久精品成人免费网站| 免费人成视频x8x8入口观看| 欧美成人性av电影在线观看| 国产成人影院久久av| 欧美日韩视频精品一区| 午夜福利在线观看吧| 波多野结衣一区麻豆| 女警被强在线播放| 三级毛片av免费| 国产单亲对白刺激| 最新美女视频免费是黄的| 国产成人精品久久二区二区免费| 黄色成人免费大全| 亚洲三区欧美一区| 日韩免费av在线播放| 久久国产乱子伦精品免费另类| 亚洲第一av免费看| 国产精品二区激情视频| 亚洲色图av天堂| av免费在线观看网站| 麻豆久久精品国产亚洲av | 色精品久久人妻99蜜桃| 亚洲五月婷婷丁香| 视频区图区小说| 亚洲一区高清亚洲精品| 中文字幕最新亚洲高清| 在线国产一区二区在线| 美女国产高潮福利片在线看| 黄色成人免费大全| 三上悠亚av全集在线观看| 男女午夜视频在线观看| 宅男免费午夜| 成年人免费黄色播放视频| 亚洲性夜色夜夜综合| 日日夜夜操网爽| 国产成人免费无遮挡视频| 无限看片的www在线观看| 国产一区二区三区视频了| 女人高潮潮喷娇喘18禁视频| 精品人妻在线不人妻| 久久精品影院6| 欧美日韩中文字幕国产精品一区二区三区 | 人妻久久中文字幕网| 亚洲精品久久午夜乱码| 99国产精品99久久久久| 在线十欧美十亚洲十日本专区| 色哟哟哟哟哟哟| 无遮挡黄片免费观看| 一区二区日韩欧美中文字幕| 久久天堂一区二区三区四区| 亚洲人成77777在线视频| cao死你这个sao货| 波多野结衣一区麻豆| 日韩一卡2卡3卡4卡2021年| 免费久久久久久久精品成人欧美视频| 免费在线观看日本一区| 精品国产国语对白av| 国产视频一区二区在线看| 亚洲精品一区av在线观看| 国产高清激情床上av| 亚洲精品粉嫩美女一区| 淫秽高清视频在线观看| 又大又爽又粗| 久久久国产成人精品二区 | 88av欧美| 亚洲国产精品一区二区三区在线| 免费看十八禁软件| 久久精品aⅴ一区二区三区四区| 中文字幕人妻丝袜一区二区| 免费观看人在逋| 久久午夜亚洲精品久久| 精品久久蜜臀av无| 日日摸夜夜添夜夜添小说| 黄色成人免费大全| 精品久久久久久,| 啦啦啦 在线观看视频| cao死你这个sao货| 国产亚洲欧美98| 夜夜躁狠狠躁天天躁| 亚洲欧美一区二区三区久久| 久久香蕉精品热| 在线观看一区二区三区激情| 成年人免费黄色播放视频| 嫩草影院精品99| 18美女黄网站色大片免费观看| 精品国产一区二区久久| 欧美另类亚洲清纯唯美| 黄片播放在线免费| 欧美 亚洲 国产 日韩一| 免费在线观看影片大全网站| 两个人免费观看高清视频| av在线天堂中文字幕 | e午夜精品久久久久久久| 国产日韩一区二区三区精品不卡| 国产免费男女视频| 亚洲美女黄片视频| 国产黄a三级三级三级人| www.www免费av| 国产亚洲欧美98| 女人被躁到高潮嗷嗷叫费观| 人妻丰满熟妇av一区二区三区| 最新美女视频免费是黄的| 自拍欧美九色日韩亚洲蝌蚪91| 正在播放国产对白刺激| 久久热在线av| 久热爱精品视频在线9| 免费av毛片视频| 天堂中文最新版在线下载| 午夜福利在线观看吧| 91麻豆av在线| 欧美国产精品va在线观看不卡| 最好的美女福利视频网| 精品久久久久久久毛片微露脸| 久久香蕉激情| 国产成人啪精品午夜网站| 精品国产乱子伦一区二区三区| 一本大道久久a久久精品| 嫩草影视91久久| 女生性感内裤真人,穿戴方法视频| 嫁个100分男人电影在线观看| √禁漫天堂资源中文www| 午夜福利影视在线免费观看| 国产野战对白在线观看| 日本免费a在线| 国产亚洲欧美98| 99国产精品99久久久久| 波多野结衣高清无吗| 国产一区二区在线av高清观看| av超薄肉色丝袜交足视频| 最新在线观看一区二区三区| 99国产精品99久久久久| 91字幕亚洲| 免费高清在线观看日韩| 女人精品久久久久毛片| 咕卡用的链子| 91成人精品电影| 老汉色av国产亚洲站长工具| 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 在线观看免费视频日本深夜| 99热只有精品国产| 一区二区三区精品91| 成年女人毛片免费观看观看9| xxxhd国产人妻xxx| 精品一区二区三区视频在线观看免费 | 国产成人精品久久二区二区免费| 国产成+人综合+亚洲专区| 色哟哟哟哟哟哟| 神马国产精品三级电影在线观看 | 身体一侧抽搐| 在线永久观看黄色视频| 国产精品久久久久久人妻精品电影| 女人精品久久久久毛片| 国产免费男女视频| 99国产极品粉嫩在线观看| 一二三四社区在线视频社区8| 丝袜人妻中文字幕| 精品欧美一区二区三区在线| av电影中文网址| av天堂在线播放| 国产精品日韩av在线免费观看 | 亚洲成a人片在线一区二区| 欧美最黄视频在线播放免费 | 免费高清视频大片| 91av网站免费观看| a在线观看视频网站| 亚洲成a人片在线一区二区| 日韩国内少妇激情av| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 欧美最黄视频在线播放免费 | 久久中文字幕人妻熟女| 在线av久久热| 国产精品久久久久成人av| 精品一品国产午夜福利视频| 色精品久久人妻99蜜桃| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 长腿黑丝高跟| 99在线视频只有这里精品首页| 欧美中文综合在线视频| 看免费av毛片| av免费在线观看网站| 18禁观看日本| 午夜日韩欧美国产| 久久久久久久久中文| 大香蕉久久成人网| 最新在线观看一区二区三区| 99久久99久久久精品蜜桃| 99久久人妻综合| 国产成人精品在线电影| 久久久久久亚洲精品国产蜜桃av| 亚洲自偷自拍图片 自拍| 美女午夜性视频免费| av超薄肉色丝袜交足视频| 国产精品一区二区在线不卡| 超碰97精品在线观看| 88av欧美| 欧美日韩视频精品一区| xxxhd国产人妻xxx| 伊人久久大香线蕉亚洲五| 免费在线观看完整版高清| 日韩大码丰满熟妇| 久久影院123| 久久精品亚洲熟妇少妇任你| 久久久久国产精品人妻aⅴ院| 久久欧美精品欧美久久欧美| 国产一区二区三区在线臀色熟女 | 亚洲情色 制服丝袜| 岛国在线观看网站| av免费在线观看网站| 日韩欧美在线二视频| 91在线观看av| 亚洲成人精品中文字幕电影 | 亚洲人成伊人成综合网2020| 欧美性长视频在线观看| 一a级毛片在线观看| 成人亚洲精品一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品野战在线观看 | av片东京热男人的天堂| 国产精品乱码一区二三区的特点 | 99热只有精品国产| 夜夜爽天天搞| 很黄的视频免费| 99精品在免费线老司机午夜| 欧美成人午夜精品| 国产精品一区二区三区四区久久 | 97碰自拍视频| 欧美激情高清一区二区三区| 国产精品久久电影中文字幕| 午夜两性在线视频| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 在线国产一区二区在线| 国产成人系列免费观看| 69av精品久久久久久| 黑人猛操日本美女一级片| 久久精品国产综合久久久| 18禁裸乳无遮挡免费网站照片 | 日韩精品中文字幕看吧| 亚洲精品一区av在线观看| 亚洲欧美一区二区三区黑人| 中文字幕最新亚洲高清| 波多野结衣高清无吗| 日本a在线网址| 国产97色在线日韩免费| 国产精品久久久久成人av| 国产不卡一卡二| 桃红色精品国产亚洲av| 久久久久久久久中文| 巨乳人妻的诱惑在线观看| 最近最新免费中文字幕在线| 91麻豆av在线| 国产精品秋霞免费鲁丝片| 欧美日韩乱码在线| a在线观看视频网站| 丝袜美足系列| 国产又色又爽无遮挡免费看| 一区在线观看完整版| 又黄又粗又硬又大视频| 岛国视频午夜一区免费看| 老司机午夜十八禁免费视频| 亚洲一区二区三区不卡视频| 国产野战对白在线观看| 丝袜人妻中文字幕| 国产99久久九九免费精品| 色哟哟哟哟哟哟| 亚洲精品美女久久av网站| 国产乱人伦免费视频| 精品一区二区三区视频在线观看免费 | 曰老女人黄片| 日本欧美视频一区| 日本黄色日本黄色录像| 午夜亚洲福利在线播放| 香蕉国产在线看| 99久久综合精品五月天人人| 激情在线观看视频在线高清| 美女大奶头视频| 俄罗斯特黄特色一大片| 亚洲avbb在线观看| 一区二区三区激情视频| 亚洲熟妇熟女久久| 成人国语在线视频| 日韩精品中文字幕看吧| 黄色女人牲交| 久久性视频一级片| 国产日韩一区二区三区精品不卡| 999久久久精品免费观看国产| 天堂√8在线中文| 日韩一卡2卡3卡4卡2021年| 色综合欧美亚洲国产小说| 欧美人与性动交α欧美精品济南到| 99精品在免费线老司机午夜| 乱人伦中国视频| 在线观看免费视频日本深夜| 嫩草影视91久久| 亚洲精品中文字幕一二三四区| 久久久久久久久免费视频了| 伊人久久大香线蕉亚洲五| 成人免费观看视频高清| 国产精品秋霞免费鲁丝片| 亚洲国产欧美网| 亚洲专区中文字幕在线| 天堂动漫精品| 制服诱惑二区| xxx96com| 亚洲 欧美 日韩 在线 免费| 国产一卡二卡三卡精品| 99riav亚洲国产免费| 国产在线精品亚洲第一网站| 亚洲 国产 在线| 欧美成人性av电影在线观看| 国产aⅴ精品一区二区三区波| 伊人久久大香线蕉亚洲五| 亚洲精品一二三| 亚洲伊人色综图| 日韩视频一区二区在线观看| 亚洲免费av在线视频| 国产精品国产高清国产av| 黄片小视频在线播放| 亚洲一码二码三码区别大吗| 国产1区2区3区精品| 欧美激情 高清一区二区三区| 欧美成人免费av一区二区三区| 12—13女人毛片做爰片一| 在线看a的网站| 免费一级毛片在线播放高清视频 | 欧美日韩精品网址| 啦啦啦免费观看视频1| 母亲3免费完整高清在线观看| 两个人看的免费小视频| 纯流量卡能插随身wifi吗| 99在线人妻在线中文字幕| www国产在线视频色| 亚洲专区国产一区二区| 亚洲国产精品sss在线观看 | 黄色女人牲交| 一进一出抽搐动态| 变态另类成人亚洲欧美熟女 | av有码第一页| 大码成人一级视频| 欧美亚洲日本最大视频资源| 成年人免费黄色播放视频| 成熟少妇高潮喷水视频| 12—13女人毛片做爰片一| 老汉色av国产亚洲站长工具| 在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲精品国产一区二区精华液| 欧美日韩精品网址| 亚洲av电影在线进入| 久久人妻av系列| 欧美在线黄色| bbb黄色大片| 中文亚洲av片在线观看爽| 日韩av在线大香蕉| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av香蕉五月| 久久久国产成人免费| 女性被躁到高潮视频| 国产成年人精品一区二区 | 亚洲精华国产精华精| 看免费av毛片| 亚洲自拍偷在线| 1024视频免费在线观看| 精品高清国产在线一区| 日韩成人在线观看一区二区三区| 久久香蕉激情| 国产精品一区二区在线不卡| 欧美丝袜亚洲另类 | 18禁国产床啪视频网站| 午夜免费激情av| 欧美性长视频在线观看| 国产一区二区激情短视频| 伦理电影免费视频| 婷婷精品国产亚洲av在线| 桃红色精品国产亚洲av| 久久草成人影院| 午夜精品在线福利| 亚洲成a人片在线一区二区| 好看av亚洲va欧美ⅴa在| 亚洲久久久国产精品| 日韩一卡2卡3卡4卡2021年| 亚洲欧美激情在线| 麻豆av在线久日| 欧美不卡视频在线免费观看 | 视频在线观看一区二区三区| 国产精品久久久久成人av| 亚洲精品国产精品久久久不卡| 亚洲欧美精品综合久久99| 亚洲精品中文字幕在线视频| 久久久久九九精品影院| 欧美日韩瑟瑟在线播放| av在线天堂中文字幕 | 纯流量卡能插随身wifi吗| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 两个人免费观看高清视频| 久久香蕉激情| 丁香六月欧美| av视频免费观看在线观看| 国产97色在线日韩免费| 侵犯人妻中文字幕一二三四区| 免费高清在线观看日韩| 高潮久久久久久久久久久不卡| 成人三级黄色视频| 久99久视频精品免费| 侵犯人妻中文字幕一二三四区| av在线天堂中文字幕 | 国产精品亚洲一级av第二区| 日韩视频一区二区在线观看| 人人澡人人妻人| 婷婷六月久久综合丁香| 亚洲欧美激情综合另类| 精品国产超薄肉色丝袜足j| 最近最新中文字幕大全免费视频| 在线观看免费高清a一片| 国产av一区二区精品久久| 桃色一区二区三区在线观看| 亚洲精品在线美女| 国产精品久久久久久人妻精品电影| 高清欧美精品videossex| 国产成人精品在线电影| 免费观看精品视频网站| 亚洲成人国产一区在线观看| √禁漫天堂资源中文www| 国产97色在线日韩免费| 国产精品亚洲一级av第二区| 中文字幕人妻丝袜制服| 亚洲成人国产一区在线观看| 999久久久国产精品视频| 一二三四社区在线视频社区8| 免费在线观看完整版高清| 很黄的视频免费| 国产高清视频在线播放一区| 麻豆av在线久日| 色综合婷婷激情| 国产精品1区2区在线观看.| 午夜福利一区二区在线看| 欧美乱妇无乱码| 777久久人妻少妇嫩草av网站| 欧美日本中文国产一区发布| 成人亚洲精品av一区二区 | 91av网站免费观看| 黄色毛片三级朝国网站| 99精品欧美一区二区三区四区| 亚洲精品国产一区二区精华液| 男女午夜视频在线观看| 一进一出抽搐动态| 色哟哟哟哟哟哟| 一区二区三区精品91| 黄色成人免费大全| 狠狠狠狠99中文字幕| 韩国av一区二区三区四区| 美女高潮到喷水免费观看| 91老司机精品| 亚洲激情在线av| www.精华液| 午夜a级毛片| 99久久99久久久精品蜜桃| 极品教师在线免费播放| 一级a爱视频在线免费观看| 青草久久国产| 最好的美女福利视频网| 黄色丝袜av网址大全| 国产黄a三级三级三级人| 无限看片的www在线观看| 午夜福利在线观看吧| 少妇粗大呻吟视频| 精品欧美一区二区三区在线| 亚洲,欧美精品.| 久久 成人 亚洲| 亚洲精品中文字幕一二三四区| 国产视频一区二区在线看| 女警被强在线播放| 国产精品香港三级国产av潘金莲| 精品国产乱子伦一区二区三区| 99精国产麻豆久久婷婷| 变态另类成人亚洲欧美熟女 | 亚洲欧美日韩高清在线视频| 久久久精品欧美日韩精品| 老司机在亚洲福利影院| 在线观看日韩欧美| 天天躁狠狠躁夜夜躁狠狠躁| 久久热在线av| 国产有黄有色有爽视频| 午夜两性在线视频| 99久久精品国产亚洲精品| 两个人免费观看高清视频| 啦啦啦在线免费观看视频4| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸 | 女同久久另类99精品国产91| 人人妻人人澡人人看| 一a级毛片在线观看| 亚洲精品久久午夜乱码| 中文字幕另类日韩欧美亚洲嫩草| 亚洲专区国产一区二区| 精品国产一区二区久久| 黄色女人牲交| 午夜日韩欧美国产| 18美女黄网站色大片免费观看| 亚洲全国av大片| 精品久久久久久成人av| 国产三级黄色录像| 美女国产高潮福利片在线看| av福利片在线| 精品一区二区三卡| 亚洲色图综合在线观看| 高清黄色对白视频在线免费看| 亚洲精品中文字幕在线视频| 久久久久国产一级毛片高清牌| 国产成年人精品一区二区 | 国产黄色免费在线视频| 人人妻,人人澡人人爽秒播| 午夜久久久在线观看| 亚洲精品美女久久av网站| 色婷婷久久久亚洲欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 99re在线观看精品视频| 中文字幕av电影在线播放| 人成视频在线观看免费观看| 18禁观看日本| 免费av中文字幕在线| 午夜福利免费观看在线| 精品少妇一区二区三区视频日本电影| 日韩一卡2卡3卡4卡2021年| 欧美午夜高清在线| 精品久久久精品久久久| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 人人澡人人妻人| 一夜夜www| 麻豆久久精品国产亚洲av | 中文字幕人妻丝袜一区二区| 99国产精品一区二区三区| 午夜精品在线福利| 精品一区二区三卡| 一级毛片精品| 一级,二级,三级黄色视频| 亚洲伊人色综图| 91大片在线观看| 欧美激情高清一区二区三区| 不卡一级毛片| av网站在线播放免费| 亚洲精品国产精品久久久不卡| 精品福利观看| 99久久国产精品久久久| 久久天堂一区二区三区四区| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| avwww免费| 美女福利国产在线| www.www免费av| 亚洲成人精品中文字幕电影 | 另类亚洲欧美激情| 黑人操中国人逼视频| 久久精品亚洲av国产电影网| 天天影视国产精品| 国产99久久九九免费精品| 夜夜爽天天搞| 日本五十路高清| 精品第一国产精品| 在线看a的网站| 自线自在国产av| 制服人妻中文乱码| 亚洲精品中文字幕一二三四区| 久久久久久人人人人人| www.熟女人妻精品国产| 天天添夜夜摸| 亚洲一区二区三区色噜噜 | 精品久久久久久,| 精品国内亚洲2022精品成人| 日韩人妻精品一区2区三区| 国产av在哪里看| 50天的宝宝边吃奶边哭怎么回事| cao死你这个sao货| 99国产极品粉嫩在线观看| 亚洲情色 制服丝袜| 亚洲av熟女| 波多野结衣高清无吗| 亚洲情色 制服丝袜| 男女下面插进去视频免费观看| 韩国精品一区二区三区| 欧美久久黑人一区二区| 99国产精品一区二区蜜桃av| 乱人伦中国视频| 淫妇啪啪啪对白视频| 一个人观看的视频www高清免费观看 | 色婷婷久久久亚洲欧美| 欧美在线一区亚洲| 久久久久久免费高清国产稀缺| 一进一出抽搐动态| 涩涩av久久男人的天堂| 久久久国产一区二区| 美女高潮到喷水免费观看| 一边摸一边做爽爽视频免费| 中出人妻视频一区二区| 成人亚洲精品一区在线观看| 亚洲精品在线美女| 欧洲精品卡2卡3卡4卡5卡区| 国产有黄有色有爽视频| 国产欧美日韩一区二区三区在线| 91大片在线观看|