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

    烏梁素海水體砷存在形態(tài)模擬及影響因素分析

    2015-12-06 06:45:44史小紅趙勝男李暢游吳用
    生態(tài)環(huán)境學(xué)報 2015年3期
    關(guān)鍵詞:烏梁素海還原性湖泊

    史小紅,趙勝男,李暢游,吳用

    內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古 呼和浩特 010018

    烏梁素海水體砷存在形態(tài)模擬及影響因素分析

    史小紅,趙勝男*,李暢游,吳用

    內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古 呼和浩特 010018

    為闡明湖泊水體中砷(As)的存在形態(tài)及水環(huán)境條件對重金屬形態(tài)變化的影響規(guī)律,對湖泊重金屬污染修復(fù)治理方面提供科學(xué)依據(jù),采集烏梁素海水樣,對湖水中As含量與存在形態(tài)進行分析,在考慮As污染程度的基礎(chǔ)上,利用PHREEQC研究湖泊水體中As的存在形態(tài),并分析溫度、酸堿度及氧化還原條件對As存在形態(tài)的影響。結(jié)果表明,水樣中As含量范圍在2.27~11.98 μg·L-1之間,平均含量為6.67 μg·L-1(n=15,std=3.41),其所有監(jiān)測點的含量都未超出地表水標準和國家漁業(yè)用水標準. 水體中 As的主要存在形態(tài)為五價態(tài)的 HAsO42-;N13采樣點 As的主要存在形態(tài)較為特殊,主要為三價態(tài)H3AsO3,含量達到 97.59%。溫度對 As各種存在形態(tài)最終分布無大影響。As存在形態(tài)受 pH影響較大,水體 pH較低時(pH=6~7.5),As主要以配合反應(yīng)為主;隨著pH的升高(pH=7.5~11),As逐漸發(fā)生水解反應(yīng)。水體的氧化還原性對As影響較明顯,當(dāng)水體處于還原狀態(tài)時,As主要以H3AsO3和H2AsO3-占主導(dǎo)地位;處于氧化性狀態(tài)時,HAsO42-和H2AsO4-占主導(dǎo)地位。因此對于寒旱區(qū)湖泊,應(yīng)對冰封期湖泊環(huán)境與水體砷含量予以重視。

    烏梁素海;影響因素;存在形態(tài);PHREEQC

    砷(As)是毒性極強的元素,其污染具有持久性、易遷移性、高度生物富集性、強毒性等特性,可通過急性或慢性暴露,對人的皮膚、呼吸系統(tǒng)、神經(jīng)及造血系統(tǒng)產(chǎn)生危害,對人體健康造成極大的傷害(馮衛(wèi)衛(wèi)等,2011;陳保衛(wèi)等,2009;Zhou等,2008;石榮等,2007)。As及其可溶性化合物均有劇毒,已被美國環(huán)保署(EPA)與國際癌癥研究機構(gòu)(IARC)歸類為一級致癌物(Straif等,2009)。由于自然地質(zhì)和人為因素的雙重作用,導(dǎo)致全球數(shù)億人的飲用水、糧食安全面臨 As污染的威脅,因此環(huán)境中 As污染成為全球關(guān)注的重大環(huán)境與健康問題(Consumer Reports,2012;Tao等,2006)。

    目前我國在湖泊 As污染研究方面已經(jīng)開展了一些工作,主要集中于水體、沉積物和水生生物等介質(zhì)中 As的含量、污染狀態(tài)、空間分布特征、生物富集與毒性等方面(張玉寶等,2011;張玉璽等,2012;李波等,2008;徐澤新等,2013;張大文等,2014)。需要注意的是,這些研究均是關(guān)于非寒旱區(qū)的南方湖泊,針對一般外界環(huán)境條件下的 As污染問題,而目前針對寒旱區(qū)湖泊 As污染研究開展的工作十分有限。烏梁素海流域地處中國北方寒旱區(qū),屬于典型灌區(qū)湖泊,主要污染源來自于流域內(nèi)的農(nóng)業(yè)面源污染,化肥流失所攜帶著 As的流入占湖泊 As來源的較大部分。流域內(nèi)礦產(chǎn)資源豐富,每年大約有2×108m3的工業(yè)廢水的排入也會造成烏梁素海水環(huán)境 As的富集與污染。目前,烏梁素海關(guān)于 As的研究集中于湖泊沉積物方面,烏梁素海沉積物中 As的含量在大部分采樣點超標,均值為27.16 mg·kg-1,遠大于背景值9.68 mg·kg-1(Zhao等,2014),按照生態(tài)危害臨界值(TEL)標準,As監(jiān)測值大約是臨界值2~6倍(Zhao等,2014)。沉積物中As的主要存在形態(tài)為殘渣態(tài)(王爽等,2012)?;趯趿核睾3练e物 As的研究可知,烏梁素海面臨潛在As污染風(fēng)險,因此,其水環(huán)境中As研究勢在必行。此外,烏梁素海是內(nèi)蒙古自治區(qū)的主要漁業(yè)基地,每年產(chǎn)魚量約在700多萬噸(Zhao等,2012),其湖水終將退水黃河,每年在枯水期向黃河補水近 2×108m3(趙勝男等,2014),由于 As的特殊物理化學(xué)性質(zhì)和強毒性,因此其水質(zhì)會影響湖區(qū)漁業(yè)與黃河水體安全。

    水體中 As元素的毒性只是由其部分結(jié)合形態(tài)所造成,不同形態(tài)其毒性高低不同。生物對 As的富集程度也與組分存在形態(tài)有關(guān),同時 As的不同形態(tài)對于動植物具有不同的生物有效性和毒性,因此需要在研究As總量的基礎(chǔ)上,對As在湖泊水體中的存在形態(tài)進行研究。水環(huán)境中砷的形態(tài)很多,主要包括砷酸鹽[As(Ⅴ)]、亞砷酸鹽[As(Ⅲ)],還有少量的一甲基胂酸鹽(MMA)和二甲基胂酸鹽(DMA)等(黃清輝等,2008)。不同形態(tài)As生物毒性差異很大,As的低氧化態(tài)比高氧化態(tài)的毒性更大,對于人體而言,亞砷酸比砷酸的毒性大60倍(趙娟等,2013)。因此不能單純依靠總 As含量反映其污染狀況和潛在風(fēng)險,而有必要對水體中 As形態(tài)特征進行研究。此外,水體氧化還原條件、溫度、酸堿度、有機質(zhì)等是對水體 As形態(tài)分布具有較大的影響(Senn等,2002;Audry等,2011;Couture等,2010)。水體的氧化還原性、酸堿性、溫度、離子強度、各種有機及無機組分的種類和濃度的變化會影響到水環(huán)境中 As元素的絡(luò)合與離解、吸附與解吸、沉淀與溶解、氧化與還原等過程,最終導(dǎo)致 As形態(tài)發(fā)生變化。烏梁素海地處中國北方的寒旱區(qū),多年平均氣溫為5.8 ℃,最冷時間出現(xiàn)在1月份,平均為-11~15 ℃,極端低溫為-41 ℃,湖泊大約每年11月下旬開始結(jié)冰,翌年3月開始消融,其季節(jié)特征較為明顯,湖泊在一年中大約有 4─7個月處于冰封狀態(tài),且冰層較厚,體積上占湖泊水體比例較大(Zhang等,2012),因此,湖泊外界環(huán)境條件包括(溫度、氧化還原性)等必然會影響到 As在湖泊中的賦存形態(tài),因此分析水體溫度、氧化還原性、酸堿性的改變對As存在形態(tài)的影響,對提出內(nèi)陸寒區(qū)湖泊具有地區(qū)特點的砷污染治理新思路具有重要的科學(xué)指導(dǎo)意義。基于烏梁素海生態(tài)功能定位、地區(qū)特殊性及砷污染的危險性,研究烏梁素海水環(huán)境中砷含量、分布特征、賦存形態(tài)及影響因素具有重大意義。

    1 材料與方法

    1.1 研究區(qū)概況

    烏梁素海位于內(nèi)蒙古自治區(qū)巴彥淖爾市烏拉特前旗境內(nèi),其地理坐標介于北緯40°36′~41°03′,東經(jīng)108°43′~108°57′之間,現(xiàn)有水域面積293 km2,其中蘆葦區(qū)面積為 118.97 km2,明水區(qū)面積為111.13 km2,明水區(qū)中85.7 km2為沉水植物密集區(qū)。湖水于每年11月初結(jié)冰,翌年3月末到4月初開始融化,冰封期約為5個月。氣溫變化較大,多年平均氣溫為7.3 ℃,全年日照時數(shù)為3185.5 h,多年平均降雨量為224 mm,蒸發(fā)量為1502 mm(趙勝男等,2013;Ren等,2008)。由于河套灌區(qū)化肥和農(nóng)藥的用量在不斷加大,與此同時化肥利用率僅為30%左右,上游工業(yè)廢水、生活污水伴同大量流失化肥的農(nóng)田排水經(jīng)不同的排水溝進入烏梁素海,致使烏梁素海水生環(huán)境日益惡化。

    1.2 取樣與分析

    根據(jù)烏梁素海面積大小、水體流向等情況,共設(shè)置15個水體采樣點,具體布點見圖1。取樣點使用全球GPS定位,破冰取水。

    圖1 水體采樣點Fig. 1 Sampling sites

    在進行水質(zhì)采樣前,采集水樣儀器與盛取水樣的聚乙烯取樣瓶必須使用專用洗液去除油污等,并用超純水洗凈,再用濃度為10%的硝酸浸泡24 h以上,取出用超純水沖洗干凈,取樣過程中,防止污染。利用自制采水器(專利號:ZL201020680966.1)取得水樣后,將水樣放于事先準備好的取樣瓶中。為保持樣品處于低溫狀態(tài),將水樣放置于低溫冷藏箱中,使溫度保持在-4~4 ℃,樣品運回實驗室后,立即進行元素 As分析測定、實驗所用玻璃器皿均用1∶3(體積比)硝酸浸泡24 h以上,并用超純水反復(fù)清洗,備用。相對標準偏差(RSD)小于3.0%。

    水體中 As與其他指標的測定按照《水和廢水監(jiān)測分析方法》(水和廢水監(jiān)測分析方法編委會,2002),采用雙道–原子熒光光度計進行測定(中國林業(yè)科學(xué)研究院分析中心,1993)。水體溫度、pH值、氧化還原電位采用便攜式 pH計進行測定。水體陰陽離子含量利用IS-90離子色譜儀進行測定。

    1.3 研究方法

    PHREEQC軟件是以熱力學(xué)平衡常數(shù)方法為基礎(chǔ)的水文地球化學(xué)模擬軟件程序。主要有4部分組成,包括數(shù)據(jù)庫、輸入文件、標準輸出文件和選擇性輸出文件。其中數(shù)據(jù)庫文件給出了主要離子、礦物質(zhì)、吸附交換、動態(tài)和平衡化學(xué)反應(yīng)等的表達式和常數(shù)。輸入文件是需要用戶編寫的文本文件,文件給出命令(反應(yīng)模式)供模型讀入并進行模擬,也可以在此文件中對數(shù)據(jù)庫進行修改和特別選擇計算輸出結(jié)果。標準輸出文件是PHREEQC在模擬運算過程中輸出結(jié)果,選擇性輸出文件是根據(jù)用戶需要選擇性輸出計算結(jié)果,包括表格輸出和圖形制作兩部分。PHREEQC軟件能夠?qū)ξ镔|(zhì)組分、液相組分形態(tài)、溶解、相、交換作用和表面絡(luò)合作用進行模擬與計算(David等,1999)。目前關(guān)于水文地球化學(xué)模擬PHREEQC軟件應(yīng)用關(guān)于采礦水、地下水中的物質(zhì)存在形態(tài)研究較多,方法技術(shù)較為成熟,取得了較多的研究成果,包括采礦區(qū)地下水中鈾的存在形態(tài)模擬等(Merkel等,2000;Karmegam,2011);海水中痕量元素組分存在形態(tài)及其生物有效性的研究等(魯靜等,2006)。

    HSC Chemistry是集成熱力學(xué)數(shù)據(jù)庫軟件,擁有一個超過20000種無機物詳細熱力學(xué)性質(zhì)的數(shù)據(jù)庫以及針對不同應(yīng)用而設(shè)計的22個計算模塊。

    本文利用PHREEQC研究烏梁素海湖水水體中As的存在形態(tài)及影響因素。首先,測定烏梁素海水體中 As的總量,水體溫度,水體酸堿度,水體氧化還原電位,水體中陰陽離子含量,通過上述數(shù)據(jù)定義水體的化學(xué)組分與環(huán)境特性。將每個水樣監(jiān)測點的實測數(shù)據(jù)按要求設(shè)置成PHREEQC軟件所要求的輸入格式,編制輸入文件。第二步,通過資料查詢的不同 As存在形態(tài),獲取對應(yīng)的水溶態(tài)物種和相關(guān)的反應(yīng)方程式,對模擬數(shù)據(jù)庫進行參數(shù)補充。第三,在獲取反應(yīng)方程式與輸入文件后,運行PHREEQC軟件,得到As在不同狀態(tài)下的存在形態(tài)及濃度的輸出文件。第四,通過改變初始模擬環(huán)境,模擬不同溫度、不同pH、不同氧化還原電位狀態(tài)下,As離子的不同存在形態(tài)及比例。最終,根據(jù)輸出文件,分析外界環(huán)境變化對湖泊水體中離子存在形態(tài)的影響。

    數(shù)據(jù)庫的建立:烏梁素海湖泊水體中存在著大量的無機陰離子,有機配合物等,本文僅考慮單一離子和絡(luò)合陰離子所表示的存在形態(tài)。As在水體中以離子或無機絡(luò)離子形式所示(Merkel,2005),砷(As)的主要存在形態(tài)為 H3AsO3、H2AsO3-、HAsO32-、AsO33-、H4AsO3+、H2AsO4-、H3AsO4、HAsO42-、AsO43-。本文利用HSC Chemistry v5.0軟件以及物理化學(xué)理論,新建了 As元素的存在形態(tài)模擬數(shù)據(jù)庫,結(jié)果如表1所示:

    表1 As的反應(yīng)方程式及參數(shù)Table 1 Chemical reaction equations and parameters of Arsenic

    輸入?yún)?shù)設(shè)置:溫度 T(攝氏度);pH值;pe值,是平衡狀態(tài)下的電子活度,在某種意義上相當(dāng)于氧化還原電位,其值可以相互推導(dǎo),pe值與氧化還原電位 Eh的關(guān)系為:pe=[(n×F)/(2.303×R× T)]×Eh,pe值的獲取可以由實測的氧化還原電位值Eh獲??;陰陽離子值:本文涉及到Ca,Mg,K,Na,Cl以及SO42-為2011年1月實測數(shù)據(jù)。As元素的總濃度為2011年1月同期實測數(shù)據(jù)。由于部分監(jiān)測點的陰陽離子數(shù)據(jù)缺失,選擇I12、K12、L15、M12、N13、P9、Q8、R7及W2等9個點進行砷形態(tài)的分析。

    2 結(jié)果與討論

    2.1 烏梁素海水環(huán)境特征分析

    烏梁素海位于中國北方寒旱區(qū),屬于灌區(qū)草型湖泊,是寒旱區(qū)典型氣候條件和農(nóng)牧交錯區(qū)典型下墊面條件下形成的湖泊,其湖泊水環(huán)境具有一定的特殊性。作者所在團隊為內(nèi)蒙古農(nóng)業(yè)大學(xué)水環(huán)境科研團隊,在烏梁素海從2003年開始進行了連續(xù)12年的每月1次的連續(xù)監(jiān)測,為了描述烏梁素海四季水溫變化特征,選擇2007─2012年的約130個水樣,以1、5、7、10月作為冬、春、夏、秋季的代表月份進行分析。

    湖泊水體溫度四季變化較大(圖2)。烏梁素海湖泊水體冬季水溫最低,平均值為0.36 ℃,最高值1.40 ℃,最低值 0.20 ℃;春季水溫逐步回升,平均值為18.4 ℃,最高值25.40 ℃;夏季水溫達到最高值,平均值為26.3 ℃,最高值33.50 ℃,最低值22.20 ℃;秋水溫逐步降低,平均值為10.4 ℃,最高值19.0 ℃,最低值6.80 ℃;烏梁素海湖泊水體溫度四季變化明顯,夏冬溫度差達到30 ℃。

    烏梁素海湖泊水體 pH值四季變化范圍較小,春季 pH在 6.86~9.65之間變化,夏季 pH在7.35~10.66之間變化;、秋季pH在7.6~10.96之間變化;冬季pH在6.86~9.65之間變化。說明烏梁素海的水環(huán)境呈中性—堿性。

    圖2 不同季節(jié)水體水溫、pH、ORP變化箱圖Fig. 2 Box-plot of water temperature, pH, ORP change In different seasons

    烏梁素海湖泊水體氧化還原電位值隨著季節(jié)變化而發(fā)生的波動較大。冬季氧化還原電位值最低,平均值為70 mV,最低值-367 mV,說明冬季湖泊水體處于強還原性;春季氧化還原電位值逐步回升,平均值為102 mV,最高值208 mV,最低值-146 mV;夏季氧化還原電位值最高,平均值為194 mV,最高值370 mV,最低值27 mV,說明夏季湖泊水體處于強氧化性;秋季氧化還原電位平均值為177 mV,最高值304 mV,最低值11 mV。說明冬季湖泊水體處于強還原性,主要是由于烏梁素海冰封期約為5─7個月,湖泊表面的冰蓋形成使得湖泊水體、沉積物與大氣隔絕,處于厭氧環(huán)境,造成氧化還原電位降低,從而形成強還原條件。

    2.2 烏梁素海水體中As總量評估

    單因子指數(shù)評價法是將各指標濃度代表值與評價標準逐項對比,以單項評價最差項目的類別作為水質(zhì)類別,是中國地表水環(huán)境質(zhì)量標準(GB3838-2002)規(guī)定的評價方法。

    烏梁素海表層水體中As含量范圍在2.27~11.98 μg·L-1之間,平均含量為6.67 μg·L-1(表2),其所有監(jiān)測點的含量都未超出地表水標準和國家漁業(yè)用水標準。以漁業(yè)用水標準作為評價標準,利用單因子指數(shù)評價法對As含量進行評價,其值在0.04~0.23之間,均值為0.13,均小于 1,說明烏梁素海水體還未受到 As污染。但由于諸多水砷污染引發(fā)的砷中毒事件的發(fā)生,國際上通常將飲用水中砷含量的最大允許值設(shè)定為50 μg·L-1,世界衛(wèi)生組織和環(huán)境保護組織從保護人體健康出發(fā),將這一值設(shè)定為10 μg·L-1(Environmental protection agency,1998;Who,2001),將烏梁素海水體與其比較,則部分采樣點稍高于這一標準值,由于烏梁素海的退水將匯入黃河,烏梁素海退水水質(zhì)將直接影響包頭取水口引用水安全,加之As的毒性系數(shù)是10(Li等,2012),基于烏梁素海的生態(tài)功能定位,要對 As的危害與危險程度予以重視。

    表2 烏梁素海水體中As含量統(tǒng)計分析Table 2 Statistical description of the As concentrations in water of Lake Ulansuhai μg·L-1

    2.3 烏梁素海水體中As存在形態(tài)模擬

    關(guān)于湖泊水體 As的研究主要集中總量分析,總量雖然能夠?qū)λw污染狀況做出簡單判斷,但是不能精確的分析出 As的生物有效性,存在一定的局限性。水體中 As的存在形態(tài)隨著其化合價的不同,毒性也發(fā)生改變,三價 As的毒性要遠遠高于五價As,三價As的不同存在形態(tài)的毒性也有所差異,毒性大小為三氫化砷>亞砷酸鹽>砷酸鹽>砷的有機化合物。因此,只有摸清 As的存在形態(tài),才能進一步探求生物地球化學(xué)循環(huán)過程。本文利用HSC Chemistry v5.0軟件以及物理化學(xué)理論,新建了As元素的存在形態(tài)模擬數(shù)據(jù)庫,利用PHREEQC模型進行存在形態(tài)模擬,在PHREEQC的輸出文件中得到各種物質(zhì)的濃度是摩爾濃度與離子活度,由于計算結(jié)果冗長,為了便于結(jié)果的簡化,更直觀的得出元素的存在形態(tài),故本文將各形態(tài)的含量以百分比形式列出,模擬結(jié)果見圖3。

    烏梁素海大部分水體中的 As的主要形式為五價態(tài)的HAsO42-;從圖3中可以看出,所有水樣中HAsO42-占絕對優(yōu)勢,含量在62.58%~98.18%;N13采樣點 As的主要形式較為特殊,主要為三價態(tài)H3AsO3,含量達到97.59%;L15點的As的主要形式有差異,五價態(tài)的HAsO42-含量比三價態(tài)H3AsO3的含量相對多一些,HAsO42-達到62.58%,H3AsO3含量達33.80%。大部分水體中五價態(tài)的H2AsO4-的含量在 1.70%~8.16%,占到一定的比例。水體中H2AsO3-,HAsO32-,AsO33-,H4AsO3+,H3AsO4,AsO43-幾種存在形態(tài)的總量僅占0.0088%~2.39%。

    圖3 烏梁素海水體中As的主要存在形態(tài)及含量百分比Fig. 3 The Arsenic speciation and content percentage in lake water

    天然水體中砷的形態(tài)主要受氧化還原電位(Eh)和酸堿性(pH)控制,在氧化環(huán)境中,HAsO42-和H2AsO4-占優(yōu)勢,pH小于6.9時,H2AsO4-占優(yōu)勢,當(dāng)pH較高時,HAsO42-占優(yōu)勢;在還原環(huán)境中,pH小于9.2時,H3AsO3占優(yōu)勢。結(jié)合烏梁素海水體的氧化還原性分析,采樣點N13,L15的氧化還原電位處于負值,屬于還原性,因此,水體中 H3AsO3含量居多,其它采樣點的水體的 Eh較高,處于氧化性,水體pH值處于8.0附近,因此,水體中HAsO42-占優(yōu)勢,符合正常規(guī)律。

    在烏梁素海四季中,氧化還原電位值多數(shù)在70~200 mV之間,大部分采樣點的砷形態(tài)為五價態(tài)的 HAsO42-,僅在冬季個別極端強還原區(qū)域(氧化還原電位值為-367 mV),以低價態(tài)存在,如N13點的砷形態(tài)為H3AsO3。因此,總體上,烏梁素海大部分水體砷元素的毒性較弱,但由于烏梁素海水體的氧化還原性變化較大,甚至在冬季的冰封狀態(tài)下會出現(xiàn)極強的還原性,因此,必須高度重視As污染,結(jié)合As污染來源,從源頭上進行治理與預(yù)防。

    2.4 烏梁素海水體中As存在形態(tài)影響因素分析

    在湖泊水體中,元素中不同形態(tài)間的轉(zhuǎn)化過程主要受湖泊水體酸堿性、氧化還原性、溫度以及水體中陰陽離子含量的影響。而 pH值(氫離子濃度的負對數(shù))能直接的表達水環(huán)境的酸堿程度,氧化還原電位ORP表示水環(huán)境的氧化—還原性。因此,水體的水質(zhì)指標pH和ORP在一起能更好地反映水體氧化還原反應(yīng)的實際情況。另外,在任何化學(xué)反應(yīng)過程中都伴隨著熱量的轉(zhuǎn)化,因而水溫對水中物質(zhì)的存在形態(tài)及含量也有一定的影響。

    鑒于水體中 As屬于典型的兩性元素,烏梁素海水體的氧化還原性不同,因此選擇兩個代表性采樣點(K12與N13)進行影響因素分析,K12點處于氧化性環(huán)境,N13點處于還原性環(huán)境。

    2.4.1 環(huán)境溫度對砷存在形態(tài)的影響

    模擬過程:水環(huán)境的其他條件為:pH=8,pe=1.6,保持pH與pe不發(fā)生變化,只改變溫度值,根據(jù)對烏梁素海水溫的分析,選取不同的代表溫度(t=0、10、20、30 ℃),對不同溫度下水體中 As的存在形態(tài)進行模擬。分析結(jié)果如圖3所示,圖3中以As形態(tài)為橫坐標,以不同 As形態(tài)的活度的負對數(shù)為縱坐標。

    圖4 水體中As的主要存在形態(tài)含量隨溫度變化圖Fig. 4 The changes of different heavy metals in different temperature

    從圖 4(a)和(b)中可以看出,隨著溫度的改變,As的各種存在形態(tài)發(fā)生微小的改變。結(jié)合前面分析,As的主要存在形態(tài)為五價態(tài)的HAsO42-。隨著溫度的升高,各主要存在形態(tài)含量變化不是較顯著。因此,四季中溫度對其形態(tài)變化影響小。

    2.4.2 環(huán)境酸堿性對As存在形態(tài)的影響

    模擬過程:水環(huán)境的t=0.1 ℃,pe=1.60,保持pe與溫度不發(fā)生變化,只改變pH值,根據(jù)對烏梁素海水體 pH值的分析,選取不同的代表值選取不同的代表值(pH=6、7、8、9、10、11),對水體中的存在形態(tài)進行模擬。分析結(jié)果如圖5所示。

    圖5 水體中As的主要存在形態(tài)含量隨pH變化圖Fig. 5 The arsenic speciation change with pH change

    從整體水環(huán)境層面分析,當(dāng)烏梁素海水體處于弱酸環(huán)境時,水體中As主要以三價態(tài)的H3AsO3存在,當(dāng)水體處于弱酸-中性以及中性-弱堿環(huán)境時,處于氧化性環(huán)境的水體中砷的主要存在形態(tài)時H2AsO4-和HAsO42-;處于還原性環(huán)境的水體中砷的主要存在形態(tài)為H3AsO3和HAsO42-;當(dāng)水體處于強堿環(huán)境時處于氧化性環(huán)境和還原性環(huán)境的水體中的水體中砷的主要存在形態(tài)都為HAsO42-。

    2.4.3 環(huán)境氧化還原性對As存在形態(tài)的影響

    模擬過程:水環(huán)境的pH=8,t=0.1 ℃,保持pH與溫度不發(fā)生變化,只改變pe值,選取不同的代表值(pe=-6、-5、-4、-3、-2、-1、1、3、4、5、6),對水體中As的存在形態(tài)進行模擬。分析結(jié)果如圖6所示。

    水體中 As屬于典型的兩性元素,因此存在形態(tài)變化受pe值的影響較為顯著。從圖6中可以看出,當(dāng)水體處于強還原性時,水體中基本以三價態(tài)砷存在,主要存在形式是H3AsO3,含量在97.8%左右,其他存在形態(tài)為H2AsO3-,含量較低,約在2.3%左右,二者幾乎為砷在強還原性的全部存在形態(tài)。隨著pe值的逐漸增大,H3AsO3和H2AsO3-的含量逐漸減少,HAsO42-和H2AsO4-的含量開始出現(xiàn)增大的趨勢,當(dāng)pe=-1時,三價態(tài)砷的H3AsO3含量在62.3%左右,五價態(tài)的HAsO42-含量增大到35%。當(dāng)pe值開始出現(xiàn)正值,處于氧化性狀態(tài)時,H3AsO3含量迅速降低,當(dāng)pe=1時,三價態(tài)H3AsO3含量幾乎為零,五價態(tài)的 HAsO42-和 H2AsO4-的含量迅速增加,HAsO42-達95%,H2AsO4-達5.5%,隨著pe值繼續(xù)增大,HAsO42-和H2AsO4-含量基本保持不變,達到近乎平衡狀態(tài)。因此,水體處于還原性狀態(tài),砷主要以H3AsO3和H2AsO3-占主導(dǎo)地位;當(dāng)水體處于氧化性狀態(tài),HAsO42-和H2AsO4-占主導(dǎo)地位。

    3 結(jié)論

    烏梁素海表層水體 As含量相對較低,未超出地表水標準和和國家漁業(yè)用水標準,但部分監(jiān)測點超出世界衛(wèi)生組織和環(huán)境保護組織設(shè)定的飲用水標準值。水體中As的主要形式為五價態(tài)的HAsO42-;N13采樣點 As的主要形式較為特殊,主要為三價態(tài)H3AsO3,含量達到97.59%,這主要是受水環(huán)境氧化還原電位的影響。通過模擬分析可知,水溫對水中物質(zhì)的存在形態(tài)及含量也有一定的影響,但并不顯著。水體中As存在形態(tài)受pH影響較大。氧化還原環(huán)境對 As影響較為顯著,水體處于還原性狀態(tài),砷主要以H3AsO3和H2AsO3-占主導(dǎo)地位;當(dāng)水體處于氧化性狀態(tài),HAsO42-和H2AsO4-占主導(dǎo)地位。由于As3+比As5+毒性大,雖然,大部分水體砷元素的毒性都較弱,但由于烏梁素海水體的氧化還原性變化較大,甚至?xí)霈F(xiàn)極強的還原性。因此,必須高度重視As,結(jié)合As的污染來源,從源頭上進行治理與預(yù)防。特別是對于蒙新高原湖區(qū)的其他寒旱區(qū)湖泊,由于其存在較長冰封期,冰層較厚,由于冰層的覆蓋,湖泊水體的氧化還原性變化必然會影響到As在湖泊中的賦存形態(tài)。

    圖6 水體中As的主要存在形態(tài)含量隨pe值變化圖Fig. 6 The arsenic speciation change with pe change

    AUDRY S, POKROVSKY O S, SHIROKOVA L S. 2011. Organic matter mineralization and trace element post-depositional redistribution in western Siberia themokarst lake sediments [J]. Biogeosciences, 8: 3341-3358.

    CONSUMER Reports. 2012. Aresnic in your food [J]. Consumer Reports: 22-27.

    COUTURE R M, WALLSCHLAGER D, ROSE J. 2013. Arsenic binding to organic and inorganic sulfur species during microbial sulfate reduction: a sediment flow-through reactor experiment [J]. Environmental Chemistry, 10(4): 285-294.

    DAVID L, PARKHURST C. 1999. Use Guide to phreeqc(version 2)- a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical [M]. U.S. Geological Survey: 1-309.

    ENVIRONMENTAL PROTECTION AGENCY. 1998. Integrated Information System (IRIS), Arsenic, Inorganic. ASRNU.S. Cincinnati, OH.

    KARMEGAM U, CHIDAMBARAM S, PRASANNA M V. 2011. A study on the mixing proportion in groundwater samples by using Piper diagram and Phreeqc model[J]. Chinese Journal of Geochemistry, 30(4): 490-495.

    LI C Y, ZHAO S N, ZHANG S, et al. 2012. Analysis of heavy metals sources and assessment in sediments of Ulansuhai Lake, China. Fresenius Environmental Bulletin, 21(12): 3666-3674.

    MOHAN D, PITTMAN J C U. 2007. Arsenic removal from water/wastewater using adsorbents-A critical review [J]. Journal of Hazardous Materials, 142(1-2): 1-53.

    REN C T, Li C Y, JIA K L, et al. 2008. Water quality assessment for Ulansuhai Lake using fuzzy clustering and pattern recognition [J]. Chinese Journal of Oceanology and Limnology, 26(3): 339-344.

    SENN D B., HEMOND H F. 2002. Nitrate controls on iron and arsenic in an urban lake [J]. Science, 296(5577): 2373-2376.

    STRAIF K, BENBRAHIM T L, BAAN R. 2009. A review of human carcinogens part C: metals, arsenic, dusts, and fibres [J]. The Lancet Oncology 10(5): 453-454.

    TAO Y, ZHANG S, JIAN W, et al. 2006. Effects of oxalate and phosphate on the release of arsenic from contaminated soils and arsenic accumulation in wheat [J]. Chemosphere, 65(8): 1281-1287.

    WHO. 2001. Arsenic Com pounds, Environmental Health Criteria second ed. World Health Organization, Geneva.

    ZHANG Y, Li C Y, Shi X H. 2012. The migration of total dissolved solids during natural freezing process in Ulansuhai Lake [J]. Journal of Arid Land, 4(1): 85-94.

    ZHAO S N, SHI X H, Li C Y, et al. 2014. Seasonal Variation of Heavy Metals in Sediment of Ulansuhai Lake, China [J]. Chemistry and Ecology, 30(1): 1-14.

    ZHAO S N, SHI X H, ZHANG S, et al. 2012. Analysis of heavy metals sources and assessment in sediments of Ulansuhai Lake, China. Fresenius Environmental Bulletin, 21(12): 3666-3674.

    ZHOU J, MA D, PAN J, et al. 2008. Application of multivariate statistical approach to identify heavy metal sources in sediment and waters: A case study in Yang Hong, China [J]. Environmental Geology, 54(2): 373-380.

    陳保衛(wèi), 那仁滿都拉, 呂美玲. 2009. 砷的代謝機制、毒性和生物監(jiān)測[J].化學(xué)進展, 21(2/3): 474-482.

    馮衛(wèi)衛(wèi), 羅錫明, 劉丹丹, 等. 2011. 寨上金礦礦區(qū)河流沉積物中砷的形態(tài)分析[J]. 生態(tài)環(huán)境學(xué)報, 20(4): 659-662.

    國家環(huán)境保護總局《水和廢水監(jiān)測分析方法》編委會. 2002. 水和廢水監(jiān)測分析方法[M]: 第4版. 北京: 國家環(huán)境科學(xué)出版社.

    黃清輝, 馬志瑋, 李建華, 等. 2008. 2006年春季長江口砷形態(tài)分析及其生物有效性[J]. 環(huán)境科學(xué), 29(8): 2131-2136.

    李波, 劉婭, 方芳, 等. 2008. 吉林省查干湖水中金屬元素與稀土元素分析[J]. 吉林大學(xué)學(xué)報(醫(yī)學(xué)版), 34(5): 919-922.

    魯靜, 李廣玉, 葉思源. 2006. Pb、Zn、Cd和Cr等元素的生物有效性研究——以膠州灣水生系統(tǒng)為例[J]. 地球?qū)W報, 27(6): 603-608.

    MERKEL B, FRIEDRICH B P著. 朱義年, 王焰新譯. 2005. 地下水地球化學(xué)模擬的原理及應(yīng)用[M]. 北京: 中國地質(zhì)大學(xué)出版社.

    MERKEL B, 付素蓉. 2000. 地下水中鈾的反應(yīng)運移模擬[J]. 地球科學(xué), 25(5): 451-455.

    石榮, 賈永鋒, 王承智, 等. 2007. 土壤礦物質(zhì)吸附砷的研究進展[J]. 土壤通報, 38(3): 584-589.

    王爽, 李暢游, 史小紅, 等. 2012. 烏梁素海沉積物中重金屬形態(tài)分布特征及污染狀況評價[J]. 環(huán)境化學(xué), 31(10): 7-15.

    徐澤新, 張敏. 2013. 太湖流域湖蕩濕地沉積物砷汞的空間分布及污染評價[J]. 長江流域資源與環(huán)境, 22(5): 626-632.

    張大文, 羅林廣, 張莉, 等. 2014. 鄱陽湖表層沉積物中砷及重金屬賦存形態(tài)及其潛在生態(tài)風(fēng)險[J]. 長江流域資源與環(huán)境, 23(8): 1132-1138.

    張玉寶, 徐穎, 儲昭升, 等. 2011. 洞庭湖平原中小型湖群沉積物中砷污染特征與評價[J]. 湖泊科學(xué), 23(5): 695-700.

    張玉璽, 向小平, 張英. 2012. 云南陽宗海砷的分布與來源[J]. 環(huán)境科學(xué), 33(11): 3768-3778.

    趙娟, 李育松, 卞建民, 等. 2013. 吉林西部地區(qū)高砷地下水砷的閾值分析及風(fēng)險評價[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 43(1): 251-258.

    趙勝男, 李暢游, 史小紅, 等. 2013. 烏梁素海沉積物重金屬生物活性及環(huán)境污染評估[J]. 生態(tài)環(huán)境學(xué)報, 22(3): 481-489.

    趙勝男, 史小紅, 李暢游, 等. 2014. 烏梁素海水體汞的分布特征及污染風(fēng)險評估[J]. 湖泊科學(xué), 26(2): 221-227.

    中國林業(yè)科學(xué)研究院分析中心. 1993. 現(xiàn)代實用儀器分析方法[M]. 北京:中國林業(yè)出版社.

    Species Distribution and Influence Factors for Arsenic in Wuliangsuhai Lake Water

    SHI Xiaohong, ZHAO Shengnan, LI Changyou, WU Yong
    College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, Huhhot, 010018, China

    To understand the heavy metals Arsenic speciation and influence factors in the water of Wuliangsuhai Lake, the samples were collected. The total Arsenic concentration was analyzed. As different chemical species in water have different toxity to organism, the speciation of the studied metals in WL was analyzed with the PHREEQC model. The impact of temperature, redox and pH on the heavy metal species in WL was simulated with the PHREEQC. The results showed that the average Arsenic concentrations in lake water were 6.67 μg·L-1. All sampling contents were not beyond the surface water standards and national fisheries standards. The mainly chemical species of Arsenic in water solution were HAsO42, however, the N13 water sites were H3AsO3which content is 97.59%. The change of temperature have minor impact on the speciation of Arsenic. pH influences greatly the speciation of Arsenic. In the condition of low pH (pH=6~7.5), the counter hydrolysis reaction of heavy metals is dominant. In the condition of high pH (pH=7.5~11), the hydrolysis reaction of heavy metals is dominant. Water redox condition in the lake has different effects on Arsenic speciation. In the reduction condition, H3AsO3, and H2AsO3-are the dominant species. In the oxidation condition, HAsO42-, and H2AsO4-are dominant species. So, for the cold lake, more attention should be paid to the icebound lake environment and water arsenic content.

    Wuliangsuhai Lake; Influence factors; species; PHREEQC

    10.16258/j.cnki.1674-5906.2015.03.011

    X132

    A

    1674-5906(2015)03-0444-08

    史小紅,趙勝男,李暢游,吳用. 烏梁素海水體砷存在形態(tài)模擬及影響因素分析[J]. 生態(tài)環(huán)境學(xué)報, 2015, 24(3): 444-451.

    SHI Xiaohong, ZHAO Shengnan, LI Changyou, WU Yong. Species Distribution and Influence Factors for Arsenic in Wuliangsuhai Lake Water [J]. Ecology and Environmental Sciences, 2015, 24(3): 444-451.

    國家自然科學(xué)基金項目(51339002;51269016;51269017);內(nèi)蒙古自治區(qū)草原英才創(chuàng)新團隊

    史小紅(1977年生),女,副教授,主要從事水環(huán)境保護與修復(fù)研究。E-mail: iamushixiaohong@163.com *通信作者:E-mail: zhaoshengnan205@163.com

    2015-01-19

    猜你喜歡
    烏梁素海還原性湖泊
    物質(zhì)氧化性與還原性的影響因素分析
    你相信嗎?湖泊也可以“生死輪回”
    “害羞”的湖泊
    巧用化學(xué)中的“經(jīng)驗規(guī)律”化繁為簡
    十字路口的烏梁素海
    中國三峽(2017年4期)2017-06-06 10:44:22
    奇異的湖泊
    烏梁素海地區(qū)近50年氣候變化特征分析
    氧在離子液體中電化學(xué)還原性能研究
    非負載Pt納米顆粒催化劑的電催化氧還原性能研究
    淺談烏梁素海的環(huán)境保護與水土保持工作措施
    中國水利(2015年16期)2015-02-28 15:14:47
    精品一区二区免费观看| 两人在一起打扑克的视频| 国产 一区 欧美 日韩| 国产伦在线观看视频一区| 我的老师免费观看完整版| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| 亚洲最大成人手机在线| 精品人妻一区二区三区麻豆 | 日本 欧美在线| 色视频www国产| 国产色婷婷99| 国内少妇人妻偷人精品xxx网站| 国产精品一及| 嫁个100分男人电影在线观看| 亚洲国产精品久久男人天堂| 欧美区成人在线视频| www日本黄色视频网| 18禁在线播放成人免费| 中文字幕精品亚洲无线码一区| 91麻豆av在线| 色综合欧美亚洲国产小说| 99久久成人亚洲精品观看| 亚洲不卡免费看| 国产成人av教育| 亚洲一区高清亚洲精品| 免费黄网站久久成人精品 | 婷婷色综合大香蕉| 搡女人真爽免费视频火全软件 | 免费看日本二区| 精品久久久久久久久久久久久| 精品久久久久久成人av| 十八禁国产超污无遮挡网站| 日日摸夜夜添夜夜添小说| 天天一区二区日本电影三级| 亚洲电影在线观看av| 能在线免费观看的黄片| 亚洲成a人片在线一区二区| 精品午夜福利视频在线观看一区| 国产综合懂色| 窝窝影院91人妻| 少妇熟女aⅴ在线视频| 在线观看免费视频日本深夜| 嫩草影视91久久| 国产不卡一卡二| 国产欧美日韩一区二区三| 噜噜噜噜噜久久久久久91| 亚洲人与动物交配视频| 波多野结衣巨乳人妻| 国产爱豆传媒在线观看| 国产一区二区激情短视频| 一级作爱视频免费观看| 99久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲熟妇熟女久久| av女优亚洲男人天堂| 亚洲三级黄色毛片| 国产高潮美女av| 国产一区二区激情短视频| bbb黄色大片| 真人做人爱边吃奶动态| 亚洲精品亚洲一区二区| 国产一区二区三区在线臀色熟女| 男人狂女人下面高潮的视频| 亚洲精品一卡2卡三卡4卡5卡| 久久99热这里只有精品18| 欧美xxxx黑人xx丫x性爽| 色尼玛亚洲综合影院| 国产精品98久久久久久宅男小说| 天堂动漫精品| 久久久久性生活片| 久久久久免费精品人妻一区二区| 最新中文字幕久久久久| 午夜福利在线观看吧| 国产精品永久免费网站| 在线天堂最新版资源| 丰满人妻一区二区三区视频av| 成人无遮挡网站| www.熟女人妻精品国产| 在线天堂最新版资源| 精品午夜福利在线看| 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 老熟妇仑乱视频hdxx| 国产亚洲欧美在线一区二区| 国产视频内射| 欧美最新免费一区二区三区 | 99久国产av精品| 成人高潮视频无遮挡免费网站| 美女免费视频网站| 久久久久久久精品吃奶| 级片在线观看| 成人国产一区最新在线观看| 最近中文字幕高清免费大全6 | 欧美最黄视频在线播放免费| 91久久精品电影网| 露出奶头的视频| 夜夜夜夜夜久久久久| 特级一级黄色大片| 日本 av在线| 国模一区二区三区四区视频| 日本一本二区三区精品| 每晚都被弄得嗷嗷叫到高潮| 日本免费a在线| 人妻久久中文字幕网| 国语自产精品视频在线第100页| 久久99热这里只有精品18| 禁无遮挡网站| 亚洲午夜理论影院| 亚洲国产精品999在线| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 亚洲一区二区三区色噜噜| 国产一区二区激情短视频| xxxwww97欧美| 日本成人三级电影网站| 深夜精品福利| 国产69精品久久久久777片| 亚洲无线观看免费| 在线观看一区二区三区| 嫁个100分男人电影在线观看| 久久久久久久亚洲中文字幕 | 人妻久久中文字幕网| av天堂中文字幕网| www.999成人在线观看| 亚洲成人久久爱视频| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 久久久国产成人免费| 赤兔流量卡办理| 亚洲欧美日韩无卡精品| 天天躁日日操中文字幕| 国产大屁股一区二区在线视频| 精品一区二区三区人妻视频| 91在线精品国自产拍蜜月| 高清在线国产一区| 成人av一区二区三区在线看| av视频在线观看入口| av黄色大香蕉| 国产野战对白在线观看| 色精品久久人妻99蜜桃| 日韩大尺度精品在线看网址| 18禁黄网站禁片免费观看直播| 男女视频在线观看网站免费| 怎么达到女性高潮| 人人妻,人人澡人人爽秒播| 狠狠狠狠99中文字幕| 亚洲天堂国产精品一区在线| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 九色国产91popny在线| 日韩人妻高清精品专区| 欧美zozozo另类| 99久久成人亚洲精品观看| 国产精品久久久久久久电影| 亚洲av电影在线进入| 国产69精品久久久久777片| 久久精品国产亚洲av涩爱 | 国产精华一区二区三区| 中文字幕av在线有码专区| av天堂在线播放| 日本免费a在线| 欧美高清性xxxxhd video| 久久国产精品影院| 国产三级黄色录像| 国产精品伦人一区二区| 在线a可以看的网站| 日本精品一区二区三区蜜桃| 欧美黄色淫秽网站| 国产精品1区2区在线观看.| 成年女人看的毛片在线观看| 国产成人福利小说| xxxwww97欧美| 全区人妻精品视频| 成年女人永久免费观看视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av成人不卡在线观看播放网| 啪啪无遮挡十八禁网站| 国产精品电影一区二区三区| 黄色女人牲交| 国产成人欧美在线观看| 亚洲一区二区三区色噜噜| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 欧美+亚洲+日韩+国产| 丝袜美腿在线中文| 日日夜夜操网爽| 亚洲国产精品999在线| 桃红色精品国产亚洲av| 亚洲无线观看免费| 久久久久久久久久黄片| 女人被狂操c到高潮| 在线天堂最新版资源| 国产黄a三级三级三级人| 蜜桃亚洲精品一区二区三区| 日韩人妻高清精品专区| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| 美女大奶头视频| 久久九九热精品免费| 在线观看舔阴道视频| 久久久成人免费电影| 无人区码免费观看不卡| 国产精品1区2区在线观看.| 欧美日韩中文字幕国产精品一区二区三区| 国产综合懂色| 午夜a级毛片| 熟女人妻精品中文字幕| 欧美高清成人免费视频www| 国产精品综合久久久久久久免费| 欧美精品国产亚洲| 亚洲五月婷婷丁香| 怎么达到女性高潮| 波多野结衣高清无吗| 国产一区二区三区在线臀色熟女| 中文字幕精品亚洲无线码一区| 中文字幕av在线有码专区| 好男人电影高清在线观看| 精品久久国产蜜桃| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| 国产午夜福利久久久久久| 少妇丰满av| 有码 亚洲区| 好男人在线观看高清免费视频| 精品99又大又爽又粗少妇毛片 | 三级男女做爰猛烈吃奶摸视频| 亚洲美女视频黄频| 男女之事视频高清在线观看| 国产精品伦人一区二区| 91在线观看av| 国产男靠女视频免费网站| 听说在线观看完整版免费高清| 成熟少妇高潮喷水视频| 亚洲 国产 在线| 精品久久久久久久末码| 欧美乱色亚洲激情| 国产国拍精品亚洲av在线观看| 人妻久久中文字幕网| 国产乱人视频| 久久久色成人| 在线天堂最新版资源| 一个人免费在线观看电影| 免费在线观看日本一区| 色综合亚洲欧美另类图片| 最新在线观看一区二区三区| 嫩草影视91久久| 午夜免费激情av| 日本熟妇午夜| 成年女人永久免费观看视频| 国内精品久久久久精免费| 少妇丰满av| 99热6这里只有精品| 精品不卡国产一区二区三区| 丝袜美腿在线中文| 99久久99久久久精品蜜桃| 亚州av有码| 校园春色视频在线观看| 99国产极品粉嫩在线观看| 欧美乱妇无乱码| 欧美激情久久久久久爽电影| 欧美性感艳星| 午夜免费男女啪啪视频观看 | 一区福利在线观看| 首页视频小说图片口味搜索| 亚洲精品亚洲一区二区| 日韩亚洲欧美综合| 女同久久另类99精品国产91| 久久亚洲真实| 成人国产综合亚洲| 一级毛片久久久久久久久女| 免费av观看视频| 亚洲av美国av| 中文在线观看免费www的网站| 一级a爱片免费观看的视频| 欧美日韩福利视频一区二区| 国内精品美女久久久久久| 欧美黑人巨大hd| 成人永久免费在线观看视频| av在线天堂中文字幕| 可以在线观看毛片的网站| 亚洲av.av天堂| 久久精品国产99精品国产亚洲性色| 日本免费一区二区三区高清不卡| 久9热在线精品视频| 日本成人三级电影网站| 国产精品久久久久久人妻精品电影| 午夜日韩欧美国产| 国产成人福利小说| 欧美一级a爱片免费观看看| 草草在线视频免费看| 久久99热6这里只有精品| 一区二区三区高清视频在线| 久久性视频一级片| 久久国产精品影院| 性插视频无遮挡在线免费观看| 国产午夜精品论理片| 琪琪午夜伦伦电影理论片6080| 99热这里只有是精品50| 91麻豆精品激情在线观看国产| 成人亚洲精品av一区二区| 亚洲真实伦在线观看| 久99久视频精品免费| 一区二区三区激情视频| 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 国产精品自产拍在线观看55亚洲| 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| 啦啦啦观看免费观看视频高清| www.色视频.com| 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大猛烈的视频| 老司机午夜福利在线观看视频| 午夜福利在线在线| 三级毛片av免费| 91久久精品电影网| 国产淫片久久久久久久久 | 欧美激情在线99| 国产伦在线观看视频一区| 久久久久久大精品| a级一级毛片免费在线观看| 十八禁网站免费在线| 国产又黄又爽又无遮挡在线| 黄片小视频在线播放| 国产成人啪精品午夜网站| eeuss影院久久| 99热6这里只有精品| 成人毛片a级毛片在线播放| 免费搜索国产男女视频| 99热精品在线国产| 久久中文看片网| 欧美成人免费av一区二区三区| 男人舔女人下体高潮全视频| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久亚洲av鲁大| 免费在线观看亚洲国产| 欧美又色又爽又黄视频| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 亚洲人成电影免费在线| 国产精品综合久久久久久久免费| 久99久视频精品免费| 人妻夜夜爽99麻豆av| 亚洲av免费高清在线观看| 日韩中文字幕欧美一区二区| 午夜福利视频1000在线观看| 亚洲 欧美 日韩 在线 免费| 毛片女人毛片| 成人av一区二区三区在线看| 国产麻豆成人av免费视频| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 亚洲 国产 在线| 一进一出抽搐动态| 亚洲精品亚洲一区二区| 亚洲人成伊人成综合网2020| 无遮挡黄片免费观看| 亚洲一区高清亚洲精品| 成人三级黄色视频| 亚洲成a人片在线一区二区| 精品人妻1区二区| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 麻豆国产av国片精品| 三级国产精品欧美在线观看| 亚洲国产色片| 97碰自拍视频| 性欧美人与动物交配| 日本与韩国留学比较| 最新在线观看一区二区三区| 亚洲成人中文字幕在线播放| 一卡2卡三卡四卡精品乱码亚洲| 18禁黄网站禁片免费观看直播| 亚洲第一欧美日韩一区二区三区| 色噜噜av男人的天堂激情| 欧美另类亚洲清纯唯美| 国产免费av片在线观看野外av| 亚洲第一电影网av| 免费电影在线观看免费观看| 高潮久久久久久久久久久不卡| 长腿黑丝高跟| 18+在线观看网站| 成年女人永久免费观看视频| 欧美日本视频| 最近中文字幕高清免费大全6 | 国内精品美女久久久久久| 男人和女人高潮做爰伦理| 亚洲第一区二区三区不卡| 久久久色成人| 此物有八面人人有两片| 日本黄色片子视频| 午夜福利欧美成人| 国产av一区在线观看免费| 国产精品一及| 99精品在免费线老司机午夜| 国产精品爽爽va在线观看网站| 韩国av一区二区三区四区| 黄色视频,在线免费观看| 欧美+日韩+精品| 搡老妇女老女人老熟妇| 亚洲av二区三区四区| 欧美成人a在线观看| 欧美最新免费一区二区三区 | 伦理电影大哥的女人| 日本 欧美在线| 国产老妇女一区| 国产精华一区二区三区| 午夜福利高清视频| 老司机午夜十八禁免费视频| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 日韩欧美三级三区| a级一级毛片免费在线观看| 一个人免费在线观看的高清视频| 老司机深夜福利视频在线观看| 嫩草影院新地址| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一小说| 免费在线观看影片大全网站| aaaaa片日本免费| 色在线成人网| 日本在线视频免费播放| www.999成人在线观看| 亚洲av中文字字幕乱码综合| 俄罗斯特黄特色一大片| 成人国产综合亚洲| 国产精华一区二区三区| 国产三级中文精品| 99久国产av精品| 久久99热这里只有精品18| 又黄又爽又刺激的免费视频.| 亚洲熟妇熟女久久| 亚洲一区二区三区不卡视频| 一个人观看的视频www高清免费观看| 在线观看免费视频日本深夜| 成人性生交大片免费视频hd| 国产av在哪里看| 特大巨黑吊av在线直播| 丰满的人妻完整版| 日本一本二区三区精品| 一个人免费在线观看的高清视频| 国产精品国产高清国产av| 日韩高清综合在线| 中文字幕av成人在线电影| 亚洲一区高清亚洲精品| 99久久精品热视频| 看免费av毛片| 综合色av麻豆| 又黄又爽又免费观看的视频| 午夜福利视频1000在线观看| 天堂动漫精品| 美女被艹到高潮喷水动态| 乱人视频在线观看| 亚洲av二区三区四区| 嫩草影院入口| 精品不卡国产一区二区三区| 午夜福利在线观看吧| 亚洲一区二区三区不卡视频| 国产精品亚洲一级av第二区| 91av网一区二区| 老司机福利观看| 亚洲精品日韩av片在线观看| 国产麻豆成人av免费视频| 国产亚洲精品久久久久久毛片| 亚洲五月天丁香| av福利片在线观看| 精品久久久久久久末码| 国产老妇女一区| 亚洲最大成人手机在线| 日韩欧美 国产精品| 中文字幕精品亚洲无线码一区| 最好的美女福利视频网| 国产成人欧美在线观看| 少妇丰满av| 中文字幕人妻熟人妻熟丝袜美| 丰满人妻一区二区三区视频av| 成人美女网站在线观看视频| 亚洲精品粉嫩美女一区| 久久久久久久久中文| 一级毛片久久久久久久久女| 亚洲国产精品合色在线| 色精品久久人妻99蜜桃| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 午夜福利高清视频| 亚洲欧美激情综合另类| 亚洲自拍偷在线| 老司机福利观看| 丰满的人妻完整版| 亚洲在线观看片| 亚洲av日韩精品久久久久久密| 国产熟女xx| 亚洲欧美日韩卡通动漫| 又紧又爽又黄一区二区| 久久午夜亚洲精品久久| ponron亚洲| 麻豆av噜噜一区二区三区| 国产一区二区在线观看日韩| 久久精品国产亚洲av天美| 久久精品国产清高在天天线| 精品久久久久久久久久久久久| 国产高清视频在线观看网站| 波多野结衣巨乳人妻| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 久久久久国产精品人妻aⅴ院| 成人特级av手机在线观看| 看免费av毛片| 丰满乱子伦码专区| 久久久久久久久久成人| 夜夜爽天天搞| 国产午夜精品论理片| 岛国在线免费视频观看| 麻豆成人av在线观看| 免费观看的影片在线观看| 国产精品一区二区性色av| 日韩 亚洲 欧美在线| 色在线成人网| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 亚洲欧美日韩高清在线视频| 国产一区二区激情短视频| 中文字幕av在线有码专区| 精品免费久久久久久久清纯| 91九色精品人成在线观看| 国内毛片毛片毛片毛片毛片| 欧美激情在线99| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| a级毛片免费高清观看在线播放| netflix在线观看网站| 我要看日韩黄色一级片| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 人人妻人人看人人澡| 少妇的逼好多水| 国产欧美日韩精品一区二区| 亚洲熟妇熟女久久| 国产精品不卡视频一区二区 | www.熟女人妻精品国产| 亚洲国产欧洲综合997久久,| 国产伦精品一区二区三区四那| 色综合站精品国产| 国产亚洲精品综合一区在线观看| 精品乱码久久久久久99久播| 国产 一区 欧美 日韩| 成人国产综合亚洲| 午夜福利在线观看吧| 午夜两性在线视频| 欧美色视频一区免费| 91午夜精品亚洲一区二区三区 | 丰满人妻熟妇乱又伦精品不卡| 国产91精品成人一区二区三区| 男女那种视频在线观看| 一进一出好大好爽视频| 美女被艹到高潮喷水动态| 亚洲性夜色夜夜综合| 好看av亚洲va欧美ⅴa在| 最近中文字幕高清免费大全6 | 成人美女网站在线观看视频| 哪里可以看免费的av片| 午夜免费男女啪啪视频观看 | 欧美日韩瑟瑟在线播放| 51午夜福利影视在线观看| 九九在线视频观看精品| 一区二区三区免费毛片| 亚洲中文字幕日韩| 午夜精品在线福利| 在线观看av片永久免费下载| aaaaa片日本免费| eeuss影院久久| 中文字幕人妻熟人妻熟丝袜美| 九色成人免费人妻av| 18+在线观看网站| 欧美激情在线99| 夜夜看夜夜爽夜夜摸| 国内精品久久久久久久电影| 国内揄拍国产精品人妻在线| 亚洲三级黄色毛片| 一夜夜www| 成年女人看的毛片在线观看| 熟女人妻精品中文字幕| 亚洲国产精品sss在线观看| 99久国产av精品| 精品人妻视频免费看| 99热这里只有精品一区| 少妇丰满av| 日韩欧美国产在线观看| 亚洲av免费在线观看| 亚洲欧美日韩高清在线视频| 亚洲最大成人av| 少妇高潮的动态图| 嫁个100分男人电影在线观看| 日韩人妻高清精品专区| 99热6这里只有精品| 婷婷色综合大香蕉| 天堂影院成人在线观看| 女生性感内裤真人,穿戴方法视频| 国产伦在线观看视频一区| 亚洲精品影视一区二区三区av| 每晚都被弄得嗷嗷叫到高潮| 在线看三级毛片| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久久免费视频| or卡值多少钱| 婷婷亚洲欧美| 在现免费观看毛片| 亚洲 国产 在线| 欧美成人免费av一区二区三区| 国产91精品成人一区二区三区| 亚洲成av人片在线播放无| 免费搜索国产男女视频| 免费在线观看影片大全网站|