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

    松嫩平原西部地下水氟和砷的富集機(jī)理與動(dòng)態(tài)變化特征

    2023-10-26 09:40:58鄧遠(yuǎn)東冶雪艷吳亞敏杜新強(qiáng)
    中國(guó)環(huán)境科學(xué) 2023年10期
    關(guān)鍵詞:通榆縣淺層深層

    鄧遠(yuǎn)東,冶雪艷,吳亞敏,杜新強(qiáng)*

    松嫩平原西部地下水氟和砷的富集機(jī)理與動(dòng)態(tài)變化特征

    鄧遠(yuǎn)東1,2,冶雪艷1,2,吳亞敏1,2,杜新強(qiáng)1,2*

    (1.吉林大學(xué)新能源與環(huán)境學(xué)院,長(zhǎng)春 130021;2.吉林大學(xué)地下水資源與環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)春 130021)

    為揭示松嫩平原西部地下水中氟和砷的演變趨勢(shì),以通榆縣為典型研究區(qū),系統(tǒng)開展地下水氟和砷濃度現(xiàn)狀調(diào)查,利用離子比圖、自組織映射網(wǎng)絡(luò)等方法進(jìn)行水質(zhì)數(shù)據(jù)分析.結(jié)果表明:潛水氟含量受含水介質(zhì)沉淀溶解作用的制約,蒸發(fā)濃縮作用對(duì)大部分淺層高、低氟地下水形成都有一定影響,但并非高氟地下水形成的主導(dǎo)作用.水中鐵、錳氧化物的溶解促進(jìn)淺、深層地下水As的解吸釋放.1995年以來(lái)區(qū)內(nèi)中東部淺層地下水氟濃度整體水平有所降低,可能與當(dāng)?shù)卮笠?guī)模開發(fā)地下水所導(dǎo)致的地下水位下降、蒸發(fā)減弱等因素有關(guān);雖然深層高氟地下水濃度變化不大,但分布范圍擴(kuò)大,可能與淺層高氟地下水經(jīng)破損井管道滲漏進(jìn)入承壓含水層有關(guān).1995~2021年區(qū)內(nèi)中部與東北部淺層地下水砷濃度整體水平遠(yuǎn)高于世衛(wèi)組織飲用水的指導(dǎo)值(10μg/L),雖然深層高砷地下水濃度變化不大,但較淺層地下水高砷區(qū)分布更廣、富集區(qū)更多、污染更嚴(yán)重,這可能與承壓含水層的還原條件促進(jìn)砷的釋放有關(guān).因此,通榆縣深層地下水的開采,需要避開地下水高砷污染區(qū),且需要合理確定開采水平避免劣質(zhì)水入侵;同時(shí)也要關(guān)注開采井止水和開采井的破壞導(dǎo)致的串層污染問(wèn)題.

    高氟地下水;高砷地下水;富集機(jī)理;動(dòng)態(tài)變化

    地下水氟、砷污染是一個(gè)全球廣泛存在的環(huán)境問(wèn)題[1].例如,印度拉賈斯坦邦農(nóng)村地區(qū)超過(guò)60%的地下水水源氟超標(biāo),最高超標(biāo)程度達(dá)34倍[2];加納北部的氟中毒流行病區(qū),地下水氟濃度最高達(dá)67mg/L[3];中國(guó)華北平原的高氟水重災(zāi)區(qū),局部地下水氟濃度可達(dá)5.9mg/L[4].全球地下水砷污染問(wèn)題也很突出,例如:越南河內(nèi)市地下水砷污染嚴(yán)重的農(nóng)村地區(qū),作為飲用水的地下水砷平均濃度高達(dá)430μg/L[5];孟加拉國(guó)西南部恒河三角洲平原,地下水砷濃度范圍為43.1~1352μg/L,遠(yuǎn)高于世衛(wèi)組織飲用水的指導(dǎo)值(10μg/L)[6];在中國(guó)大同盆地中心,地下水砷濃度最高達(dá)366μg/L[7].飲用高氟、高砷地下水會(huì)對(duì)人類健康造成極大危害,是一系列地方性疾病的重要成因[27].識(shí)別地下水氟、砷濃度及釋放控制機(jī)制,對(duì)地下水氟、砷污染風(fēng)險(xiǎn)管控具有重要意義[8].

    在中國(guó)北方干旱半干旱地區(qū),慢性砷中毒與地氟病是常見的地方病,常見于呼和浩特盆地、大同盆地、運(yùn)城盆地、河套平原和松嫩平原.前人大量研究表明,干旱半干旱的封閉內(nèi)陸盆地、強(qiáng)還原條件下的含水層環(huán)境是地下水砷富集的有利場(chǎng)所[9-11];含氟礦物的溶解、離子交換作用和蒸發(fā)濃縮作用是高氟地下水形成的主要機(jī)制[12-15].吉林通榆縣位于松嫩平原西部,同時(shí)存在高氟和高砷地下水問(wèn)題,166個(gè)行政村一直飲用高氟水[16]、高砷水,是我國(guó)典型的氟、砷中毒地方病高發(fā)區(qū)之一[17].

    前人對(duì)通榆縣高氟地下水的分布特征、成因與富集機(jī)理做了大量的研究.早在1985年馬玉福等[18]認(rèn)為大興安嶺山地的沉積物提供了富含氟化物的礦物,而年降雨少、地下水埋深淺、徑流不暢和蒸發(fā)強(qiáng)等條件為高氟地下水形成提供了有利條件.李霄等[19]研究發(fā)現(xiàn)通榆縣1991~1995年間地下水氟含量有逐年增加的趨勢(shì).卞建民等[20]定量評(píng)價(jià)了通榆縣淺層地下水的脆弱性.張文強(qiáng)[17]認(rèn)為影響通榆縣氟中毒區(qū)域空間分布的諸地質(zhì)環(huán)境因素中,地下水動(dòng)力條件起著根本性作用.吳艷新[21]發(fā)現(xiàn)通榆縣農(nóng)村飲水安全問(wèn)題主要是高氟水、污染水等,其主要影響因素有氣候、降雨等.胡婧敏[22]發(fā)現(xiàn)通榆縣潛水水化學(xué)類型主要為HCO3-Na-Mg-Ca型,在影響氟中毒和病情輕重的各地質(zhì)環(huán)境因素中,地下水的動(dòng)力條件是至關(guān)重要的,其次是地貌因素.

    前人對(duì)通榆縣高砷地下水的研究,以高砷地下水空間分布特征、富集機(jī)理與健康風(fēng)險(xiǎn)評(píng)價(jià)等內(nèi)容為主.早期研究多為討論通榆縣高砷地下水空間分布特征或總結(jié)水文地球化學(xué)特征,例如卞建民等[23]發(fā)現(xiàn)通榆縣高砷地下水空間分布特征為:水平方向上在中部與北部的地勢(shì)低洼帶富集,垂向上在30~50m埋深的含水層水砷含量較高, Cl HCO3-Na型水的總砷含量最高.在此基礎(chǔ)上趙娟等[24]確定通榆縣地下水的砷安全閾值為0.05mg/L.隨著對(duì)地下水污染機(jī)制認(rèn)識(shí)的逐步深化,后續(xù)研究注重于高砷地下水富集機(jī)理或成因分析,例如張迪[10]發(fā)現(xiàn)鐵氧化物的還原性溶解是造成通榆縣Ⅰ、Ⅱ類地下水中高砷的主要原因,Ⅲ類地下水中硫化物的沉淀抑制了砷在深層地下水中的運(yùn)移和富集,人類農(nóng)業(yè)活動(dòng)是導(dǎo)致Ⅳ類地下水砷富集的主要原因.湯潔等[25]發(fā)現(xiàn)臭蔥石在地下水的溶解過(guò)程會(huì)形成鐵錳氧化物與砷化合物,隨著Eh降低,吸附于氧化物表面的含砷化合物會(huì)解吸進(jìn)入含水層.

    面對(duì)地氟病、慢性砷中毒等飲用水安全問(wèn)題,在國(guó)家“飲用水安全工程建設(shè)”項(xiàng)目的推動(dòng)下,通榆縣自2006年開始,相繼開展了“防病改水”重要舉措,開發(fā)埋藏較深的第三系大安組含水層地下水,有效改善了飲用水水質(zhì)狀態(tài).但在近年的調(diào)查中發(fā)現(xiàn),部分防病改水井中砷和氟的濃度出現(xiàn)增高趨勢(shì),使通榆縣供水安全再次出現(xiàn)隱患.本文旨在通過(guò)系統(tǒng)的地下水采樣和測(cè)試,分析通榆縣地下水中氟和砷的現(xiàn)狀污染水平,結(jié)合前人研究成果,對(duì)通榆縣地下水氟和砷濃度動(dòng)態(tài)變化特征及原因進(jìn)行分析,為進(jìn)一步改善通榆縣地下水環(huán)境質(zhì)量和提高飲用水安全保障程度提供科學(xué)依據(jù).

    1 材料與方法

    1.1 研究區(qū)概況

    通榆縣位于吉林省西北部,地處松嫩平原西部,地勢(shì)平坦,西北高東南低.研究區(qū)處于中溫帶半干旱大陸性氣候區(qū),多年平均氣溫5.5℃,區(qū)內(nèi)年平均降水量為287.4mm,多集中在5~8月,約占全年總降水量的90%,年平均蒸發(fā)量為1737.66mm,多集中在4~9月,約為年蒸發(fā)量的80%,年降水量小而集中,蒸發(fā)量為降水量的4~6倍.研究區(qū)內(nèi)主要河流為霍林河,發(fā)源自內(nèi)蒙古德福特勒罕的北麓,進(jìn)入向海水庫(kù)后分南北兩條支流進(jìn)入研究區(qū).本區(qū)主要出露的地層主要有第四系(Q)和新近系(N),其中新近系地層主要是泰康組和大安組.研究區(qū)為中生界以來(lái)的沉降區(qū),屬松遼斷陷盆地的西部隆起,主要構(gòu)造形跡為隆起和斷裂.

    圖1 研究區(qū)采樣點(diǎn)位置

    根據(jù)地下水賦存及巖性分布特征,本區(qū)包括2個(gè)主要含水層系統(tǒng):第四系松散巖類孔隙含水層系統(tǒng)和新近系碎屑巖類裂隙孔隙含水層系統(tǒng).研究區(qū)地下水主要補(bǔ)給方式為:大氣降水入滲補(bǔ)給、地表水滲漏補(bǔ)給、含水層間的越流補(bǔ)給、灌溉水回滲補(bǔ)給和側(cè)向地下徑流補(bǔ)給.地下水主要排泄方式為:蒸發(fā)排泄、地下徑流排泄和人工開采等.

    1.2 現(xiàn)狀調(diào)查

    2021年7月中旬共采集研究區(qū)地下水樣62份,采樣井深10~128m.其中,35份水樣來(lái)源于第四系承壓水井及新近系承壓水井(深層地下水),井深60~ 120m;27份水樣來(lái)源于第四系淺層潛水井(淺層地下水),井深在10~50m.采樣時(shí)根據(jù)《地下水環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》進(jìn)行操作,保存過(guò)程嚴(yán)格按照規(guī)范的要求進(jìn)行.水質(zhì)分析檢測(cè)項(xiàng)目包括有溶解性總固體(TDS),陽(yáng)離子:K+、Na+、Ca2+、Mg2+、Fe、Mn、As,陰離子:Cl-、HCO3-、CO3-、SO42-、F-.在采取水樣前,將水樣瓶用純水洗滌干凈,所有水樣采樣時(shí)用環(huán)境樣品潤(rùn)洗3次,水樣現(xiàn)場(chǎng)用Parafilm膜進(jìn)行密封,填寫標(biāo)簽,注明化驗(yàn)項(xiàng)目.用于金屬離子分析的樣品采用HNO3酸化至PH值小于2;用于砷分析的樣品加入濃HCL使其含量達(dá)到1%.樣品采集完畢后立即送往中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所實(shí)驗(yàn)室進(jìn)行水質(zhì)檢測(cè).

    1.3 歷史數(shù)據(jù)收集與分析

    主要通過(guò)知網(wǎng)文獻(xiàn)搜索的方式收集研究區(qū)淺深層地下水歷史砷和氟濃度數(shù)據(jù)、地下水位、地下水開采量、多年蒸發(fā)量等.收集到的水質(zhì)數(shù)據(jù)包括1995、2006、2013和2017年淺深層地下水氟濃度數(shù)據(jù),2006、2010、2013和2014年淺深層地下水砷濃度數(shù)據(jù).此外還收集到研究區(qū)2000~2010年地下水開采量、1998~2011年均蒸發(fā)量、1998~2011年第四系孔隙潛水含水層水位數(shù)據(jù).從收集到的水質(zhì)數(shù)據(jù)來(lái)看,研究區(qū)多數(shù)區(qū)域淺深層地下水砷和氟濃度普遍偏高.2000~2010地下水開采量整體呈上升趨勢(shì),1998~2011年均蒸發(fā)量整體呈下降趨勢(shì),1998~ 2011年淺層地下水位整體呈下降趨勢(shì).

    2 結(jié)果與討論

    2.1 水文地球化學(xué)特征

    (1)主要離子分布特征

    從地下水樣品的主要水化學(xué)指標(biāo)(表1)可知:

    ①淺層地下水TDS濃度變化范圍為288.42~ 1917.19mg/L,微咸水與淡水各占29.63%和70.37%;深層地下水TDS濃度變化范圍為314.35~ 1719.37mg/L,微咸水與淡水各占17.14%、82.86%.

    ②淺層地下水F-濃度變化范圍為0.04~6.24mg/ L,超標(biāo)率48.15%,變異系數(shù)為1.0,說(shuō)明樣本離散程度較高,說(shuō)明第四系淺層局部存在氟富集且污染程度高;深層地下水F-濃度變化范圍為0.25~1.79mg/ L,超標(biāo)率62.86%,變異系數(shù)為0.35,變異系數(shù)低但超標(biāo)率高,說(shuō)明高氟地下水在地下深層廣泛分布.

    ③淺層地下水As濃度變化范圍為0.86~ 114.58μg/L,超標(biāo)率18.52%,變異系數(shù)為2.23,變異系數(shù)高但超標(biāo)率低,表明第四系淺層局部存在氟富集且污染程度高;深層地下水As濃度變化范圍為0.79~109.73μg/L,超標(biāo)率51.43%,變異系數(shù)為1.18,變異系數(shù)與超標(biāo)率均較高,說(shuō)明高砷地下水在地下深層具局部富集且污染程度高.

    (2)典型污染物水化學(xué)特點(diǎn)

    根據(jù)水化學(xué)類型圖(圖2)可知:淺層高氟地下水化學(xué)類型主要為HCO3-Na和HCO3-Na·Ca型,淺層低氟地下水化學(xué)類型主要為HCO3-Ca·Na和HCO3-Ca·Mg型.深層高氟地下水化學(xué)類型主要為HCO3-Na型,深層低氟地下水化學(xué)類型主要為HCO3-Na·Ca型.

    淺層高砷地下水化學(xué)類型主要為HCO3-Na·Ca型,淺層低砷地下水化學(xué)類型主要為HCO3-Ca·Na和HCO3-Ca·Mg型.深層高砷地下水化學(xué)類型主要為HCO3-Na·Ca型,深層低砷地下水化學(xué)類型主要為HCO3-Ca·Na、HCO3-Na型.總體來(lái)說(shuō),淺層與深層地下水的水化學(xué)類型相似,但較深層地下水化學(xué)類型而言,淺層地下水樣本在piper三線圖的空間分布更趨向于HCO3-Ca·Na型水.

    表1 地下水主要化學(xué)成分統(tǒng)計(jì)

    注:As指標(biāo)單位為mg/L,其余指標(biāo)單位均為mg/L.

    2.2 地下水中氟和砷濃度分布特征

    從圖4可知,研究區(qū)整體上東南部地下水中氟含量較高,在瞻榆鎮(zhèn)與新華牛場(chǎng)一帶局部富集,最大值為6.22mg/L.從地貌地形、徑流條件等因素分析,研究區(qū)地處松遼平原西部,地勢(shì)平坦,西北高東南低,西北傾斜平原徑流條件較東南低平原好,氟含量低;而東南地形低洼,徑流條件較差,蒸發(fā)濃縮作用較強(qiáng)烈,氟含量高.淺層地下水砷含量從盆地邊緣至盆地中心逐漸升高,在烏蘭花鎮(zhèn)局部富集.總體來(lái)說(shuō)淺層地下水氟與砷含量分布特征在空間上具有相似性.從圖5可知,研究區(qū)深層地下水氟含量的分布特點(diǎn)為:西北部低,東北部高,在雙崗鹿場(chǎng)、雙崗鎮(zhèn)、鴻興鎮(zhèn)達(dá)到最高峰值.研究區(qū)深層地下水砷含量的分布特點(diǎn)為:西南部低,北部高,高砷點(diǎn)集中分布在雙崗鹿場(chǎng)、新興鄉(xiāng)一帶.

    圖4 淺層地下水氟、砷濃度空間分布

    圖5 深層地下水氟、砷濃度空間分布

    2.3 地下水中氟與砷的富集機(jī)理

    2.3.1 相關(guān)性分析 利用SOM將影響砷和氟釋放與轉(zhuǎn)化的化學(xué)變量(例如Na+、K+、Ca2+、Mg2+、HCO3?、SO42?、Cl?、TDS、Fe、Mn)在二維分量平面上表示出來(lái)(圖6、圖7),有助于量化砷、氟與多個(gè)變量之間的關(guān)系.從圖6可知,淺層地下水砷組分平面最大值位于圖的左下角,氟組分平面最大值位于圖的左下角與右下角;從圖7可知,深層地下水砷組分平面最大值位于圖的右下角,氟組分平面最大值位于圖的上方.多個(gè)變量之間的相關(guān)性可以通過(guò)在組成平面之間成對(duì)保持相同的相對(duì)拓?fù)浞植贾祦?lái)確定.為了根據(jù)SOM分析結(jié)果量化砷、氟與多個(gè)變量之間的相關(guān)性,基于SOM的輸出,采用皮爾遜相關(guān)系數(shù)來(lái)表示不同參數(shù)對(duì)砷、氟濃度的影響.從表2可知,淺層地下水F-與As沒有呈現(xiàn)明顯的顯著相關(guān)性,F-與Ca2+、Cl-均呈現(xiàn)一般相關(guān)性,說(shuō)明氟的富集和鹽效應(yīng)、螢石的沉淀溶解作用相關(guān).淺層地下水As與K+表現(xiàn)出較強(qiáng)相關(guān)性,這與地形低洼處的蒸發(fā)濃縮作用具有極大關(guān)系.從表3可知,深層地下水F-與Ca2+呈現(xiàn)中等相關(guān)性,As與Mn呈現(xiàn)一般相關(guān)性,與K+呈現(xiàn)較強(qiáng)相關(guān)性.總體來(lái)說(shuō)潛水-深層地下水典型污染物與其他水化學(xué)指標(biāo)的相關(guān)關(guān)系相似.

    圖6 淺層地下水各化學(xué)參數(shù)SOM平面展布(神經(jīng)元數(shù)量5×5)

    圖7 深層地下水各化學(xué)參數(shù)SOM平面展布(神經(jīng)元數(shù)量6×5)

    表2 淺層地下水化學(xué)成分相關(guān)系數(shù)

    2.3.2 氟的富集機(jī)理 (1)水巖作用 地下水含水介質(zhì)與水溶液接觸發(fā)生溶解沉淀作用,導(dǎo)致地下水水化學(xué)組份發(fā)生變化.研究區(qū)含水介質(zhì)主要含氟礦物為螢石(CaF2)、水氟(NaF)等.通過(guò)繪制主要離子關(guān)系圖,分析溶解沉淀作用對(duì)地下水氟富集的貢獻(xiàn).從圖8的Ca2+-SO42-關(guān)系圖可知,淺層低氟地下水較多樣本點(diǎn)位于=線上,說(shuō)明石膏(CaSO4)的溶解是該水樣點(diǎn)中離子的主要來(lái)源,而高氟地下水大多數(shù)樣本點(diǎn)位于=線以下,說(shuō)明除了石膏的溶解,含氟礦物螢石(CaF2)的溶解貢獻(xiàn)了大部分F-與Ca2+[26].從HCO3--Ca2+關(guān)系圖可知,淺層高氟、低氟地下水樣本點(diǎn)大多數(shù)位于=線上方,說(shuō)明淺層地下水化學(xué)反應(yīng)大致向方解石、石膏、螢石的溶解方向進(jìn)行.從Na+-SO42-關(guān)系圖可知,在相同SO42-濃度區(qū),低氟地下水樣本整體Na+濃度相較于高氟地下水高,說(shuō)明高氟地下水中氟濃度增大促進(jìn)水氟的析出,導(dǎo)致Na+濃度降低.從HCO3--Ca2++Mg2+關(guān)系圖可知,樣本點(diǎn)均位于=線附近,說(shuō)明白云石(CaMg(CO3)2)的溶解對(duì)研究區(qū)潛水化學(xué)演變具有一定貢獻(xiàn).

    表3 深層地下水化學(xué)成分相關(guān)系數(shù)

    從圖8的Ca2+-SO42-關(guān)系圖可知,深層低氟地下水全部樣本點(diǎn)均位于=線下方,說(shuō)明該水樣點(diǎn)中Ca2+離子的主要來(lái)源除了石膏的溶解,可能還存在白云石的溶解.從HCO3--Ca2+關(guān)系圖可知,深層高氟、低氟地下水樣本點(diǎn)大多數(shù)位于=線上方,說(shuō)明深層地下水化學(xué)反應(yīng)大致向方解石、石膏、螢石的溶解方向進(jìn)行.從Na+-SO42-關(guān)系圖可知,在相同SO42-濃度區(qū),高氟地下水樣本整體Na+濃度相較于低氟地下水高,說(shuō)明高氟地下水中水氟可能尚未析出.從HCO3--Ca2++Mg2+關(guān)系圖可知,樣本點(diǎn)均位于=線附近,說(shuō)明白云石的溶解對(duì)研究區(qū)承壓水化學(xué)演變具有一定貢獻(xiàn).

    圖9 地下水F-與Ca2+活度關(guān)系

    通過(guò)PHREEQC軟件進(jìn)行地球化學(xué)模擬,計(jì)算決定礦物溶解沉淀過(guò)程結(jié)果的飽和指數(shù)SI和Ca2+活度.由表4可知,淺、深層地下水石膏和螢石的SI平均值均小于-0.5,處于溶解階段,與離子比分析結(jié)果相似,說(shuō)明石膏和螢石的溶解對(duì)地下水氟與鈣離子的形成演化具有較大貢獻(xiàn).方解石與白云石SI平均值大于0,說(shuō)明方解石與白云石沉淀限制了Ca2+濃度升高.從圖9可知,山前區(qū)淺層地下水大部分樣品靠近趨勢(shì)線1,且F-活度升高時(shí)Ca2+活度變化不大,說(shuō)明山前淺層地下水氟富集主要受含氟礦物溶解控制.平原區(qū)淺層地下水大部分樣品位于趨勢(shì)線1上方,說(shuō)明平原區(qū)淺層地下水氟主要受螢石與方解石等礦物溶解控制.山前區(qū)深層地下水氟富集的主控因素與淺層相似,但平原區(qū)深層地下水的樣品多位于趨勢(shì)線1和2之間,說(shuō)明平原區(qū)深層地下水氟濃度受方解石、石膏與螢石等礦物溶解來(lái)源Ca2+的控制.

    (2)蒸發(fā)濃縮作用

    通過(guò)Gibbs圖分析巖石風(fēng)化作用與蒸發(fā)濃縮作用對(duì)淺層高氟地下水形成的影響,從圖10可知:淺層高氟地下水Na+/(Na++Ca2+)比值為0.1~0.8,TDS濃度約為300~1900mg/L,而部分水樣點(diǎn)氟值在一定程度上受巖石風(fēng)化作用影響,總體來(lái)講蒸發(fā)濃縮作用對(duì)大部分淺層的高、低氟地下水形成都有一定影響,但并不是高氟地下水形成的主導(dǎo)作用.

    導(dǎo)致以上現(xiàn)象發(fā)生的可能原因有兩個(gè):一是野外人為取樣、室內(nèi)水化學(xué)測(cè)試操作存在一定誤差;二是跟野外取樣時(shí)間有關(guān),取樣時(shí)正值夏季,研究區(qū)炎熱多雨,蒸發(fā)較強(qiáng)烈,有利于潛水中氟的富集,但潛水同時(shí)得到充足的大氣降雨補(bǔ)給而稀釋了水中的氟濃度,使氟濃度在1mg/L附近的地下水氟濃度降低而變?yōu)榈头叵滤?說(shuō)明蒸發(fā)濃縮作用對(duì)淺層大部分的高、低氟地下水形成都有一定影響.

    圖10 淺層地下水Gibbs圖

    2.3.3 砷的富集機(jī)理 地下水鐵、錳氧化物的還原溶解被廣泛認(rèn)為是促進(jìn)砷釋放的重要機(jī)制[5,27-28],通過(guò)As-Fe/Mn離子比圖與飽和指數(shù)SI等定性-定量方法探究砷的富集機(jī)理.

    從圖11可知,淺層及深層低砷地下水樣本點(diǎn)砷濃度與總鐵和錳基本無(wú)關(guān);而淺層和深層高砷地下水中砷濃度則與總鐵和錳濃度呈現(xiàn)一定的正相關(guān)關(guān)系,但并不明顯.通過(guò)PHREEQC計(jì)算了不同礦物的飽和指數(shù)SI(表4).由表4可知,鐵氧化物/氫氧化物礦物相的SI為正值,包括針鐵礦(淺層地下水SI平均值為6.78,深層地下水為6.80)、赤鐵礦(淺層地下水SI平均值為15.57,深層地下水為15.60)和Fe(OH)3(淺層地下水SI平均值為0.89,深層地下水為0.90),這表明水生地球化學(xué)環(huán)境有利于這些礦物相的沉淀[29].碳酸鹽礦物也可能是含水層中砷滯留/釋放的原因[30],從表4可知:鐵/錳碳酸鹽礦物的SI為負(fù)值,包括菱錳礦(淺層地下水SI平均值為-0.07,深層地下水為-0.28)和菱鐵礦(淺層地下水SI平均值為-0.62,深層地下水為-0.60),這說(shuō)明鐵/錳碳酸鹽礦物處于溶解階段,而淺層和深層高砷地下水中砷濃度則與總鐵和錳濃度呈現(xiàn)一定的正相關(guān)關(guān)系,說(shuō)明鐵/錳碳酸鹽礦物溶解可能對(duì)As的解吸有一定潛在作用.

    表4 不同礦物飽和指數(shù)(SI)描述性統(tǒng)計(jì)表

    結(jié)合上述分析可知,鐵氧化物/氫氧化物礦物沉淀與鐵/錳碳酸鹽礦物溶解可能是控制As釋放的主要因素,導(dǎo)致As從礦物中解吸到地下水中[31].

    2.4 地下水氟與砷富集主控因素

    利用主成分分析法(PCA)探討研究區(qū)地下水氟與砷富集的主控因素.基于各離子與F-的相關(guān)關(guān)系選取水質(zhì)參數(shù)F-、HCO3-、Cl-、SO42-、Na+、Mg2+、K+、Ca2+作為地下水氟的主成分分析變量.選取水質(zhì)參數(shù)As、Mn、Fe、HCO3-、Cl-、SO42-、TDS、Na+、Mg2+、K+、Ca2+作為地下水砷的主成分分析變量.KMO檢驗(yàn)值均大于0.5,<0.001,表明原始數(shù)據(jù)適合進(jìn)行主成分分析.基于特征值大于1的限定條件識(shí)別出3個(gè)主控因子,具體見圖12~圖13,方差累計(jì)貢獻(xiàn)度見表5和表6.

    表5 地下水氟主成分分析描述性統(tǒng)計(jì)表

    從圖12和表5可知,淺層地下水第一主成分(PC1)的貢獻(xiàn)率為32.91%,PC1與Na+、Mg2+、K+、Ca2+、Cl-呈強(qiáng)正相關(guān)關(guān)系,說(shuō)明PC1代表巖石風(fēng)化與蒸發(fā)濃縮作用.淺層地下水第二主成分(PC2)的貢獻(xiàn)率為28.85%,PC2與Ca2+呈強(qiáng)負(fù)相關(guān)關(guān)系、與F-呈中等正相關(guān)關(guān)系,說(shuō)明含鈣礦物沉淀在一定程度上促進(jìn)氟富集.淺層地下水第三主成分(PC3)貢獻(xiàn)率為18.08%,PC3與Na+、Ca2+呈弱正相關(guān)關(guān)系、與F-呈強(qiáng)正相關(guān)關(guān)系,說(shuō)明螢石與水氟的溶解對(duì)氟富集有積極作用.深層地下水第一主成分(PC1)與Na+、Mg2+、K+、HCO3-呈強(qiáng)正相關(guān)關(guān)系,代表碳酸鹽巖等的巖石風(fēng)化作用.PC2與Na+呈中等正相關(guān)、與F-呈強(qiáng)正相關(guān)關(guān)系,說(shuō)明水氟的溶解是深層地下水氟富集的主要機(jī)制.

    表6 地下水砷主成分分析描述性統(tǒng)計(jì)表

    注:VCR為方差貢獻(xiàn)率,CCR為累計(jì)貢獻(xiàn)率.

    圖12 含氟的地下水主要離子主成分分析

    從圖13和表6可知,淺層地下水第一主成分(PC1)的貢獻(xiàn)率為39.14%,PC1與Cl-、SO42-、Mg2+、Ca2+、TDS呈強(qiáng)正相關(guān)關(guān)系,說(shuō)明PC1代表巖石風(fēng)化與蒸發(fā)濃縮作用.第二主成分(PC2)的貢獻(xiàn)率為16.76%,PC2與Mn、總Fe呈強(qiáng)正相關(guān)關(guān)系、與As呈弱負(fù)相關(guān)關(guān)系,說(shuō)明鐵/錳碳酸鹽礦物的溶解對(duì)淺層As解吸釋放可能有一定限制作用.第三主成分(PC3)的貢獻(xiàn)率為16.46%,PC3與總Fe呈弱正相關(guān)關(guān)系、與Mn呈弱負(fù)相關(guān)關(guān)系、與As呈強(qiáng)正相關(guān)關(guān)系,說(shuō)明鐵碳酸鹽礦物溶解與錳氧化礦物沉淀促進(jìn)淺層As釋放.深層地下水第一主成分(PC1)與Na+、Mg2+、K+、Cl-、HCO3-呈強(qiáng)正相關(guān)關(guān)系,代表碳酸鹽巖等的巖石風(fēng)化作用.PC3與總Fe呈強(qiáng)負(fù)相關(guān)關(guān)系、與As呈中等正相關(guān)關(guān)系,說(shuō)明鐵氧化物/氫氧化物礦物的沉淀促進(jìn)As從礦物中解吸到地下水中.

    圖13 含砷的地下水主要離子主成分分析

    2.5 地下水中氟和砷濃度動(dòng)態(tài)變化

    宋偉[16](2018)在2016年對(duì)通榆縣防氟改水現(xiàn)況進(jìn)行了調(diào)查,707個(gè)地氟病病區(qū)均已實(shí)施改水,但部分改水工程因井內(nèi)管道滲漏破損等原因?qū)е鲁霈F(xiàn)水中氟超標(biāo)現(xiàn)象.通過(guò)對(duì)比前人研究成果,分析通榆縣進(jìn)行防氟改水后地下水中氟和砷濃度的變化特征,為下一步治理研究區(qū)地下水氟和砷污染提供科學(xué)依據(jù).

    (1)地下水氟濃度動(dòng)態(tài)變化

    淺層高氟地下水在通榆縣中部與東部低平原區(qū)富集,結(jié)合研究區(qū)地下水氟濃度文獻(xiàn)資料[18-19,32-33],在研究區(qū)中部與東部選取烏蘭花鎮(zhèn)-新華鎮(zhèn)-開通鎮(zhèn)-蘇公坨的四個(gè)代表不同水文地質(zhì)環(huán)境的行政分區(qū),繪制各分區(qū)1995年、2006年、2013年、2017年、2021年的地下水氟濃度浮動(dòng)柱狀圖(圖14),由圖14可知:

    ①除開通鎮(zhèn)以外,1995年以后各鄉(xiāng)鎮(zhèn)的淺層地下水氟濃度整體水平低于1995年的氟濃度值,可能與當(dāng)?shù)剡^(guò)度開發(fā)地下水資源導(dǎo)致地下水位下降、蒸發(fā)減弱等因素有關(guān)[34].通榆縣工業(yè)、農(nóng)業(yè)、生活用水均采自地下水,地表水資源較少,隨著城市化擴(kuò)張與耕地面積增加,地下水開采量逐年呈上升趨勢(shì)[35-36],導(dǎo)致地下水位下降.

    ②淺層地下水氟富集中心區(qū)有向東部轉(zhuǎn)移的趨勢(shì).1995至2021年五個(gè)年份,各鄉(xiāng)鎮(zhèn)淺層地下水氟濃度仍普遍高于世衛(wèi)組織飲用水的指導(dǎo)值(1mg/L).烏蘭花鎮(zhèn)、新華鎮(zhèn)和蘇公坨淺層地下水氟濃度整體水平(£2mg/L)呈現(xiàn)逐年下降的趨勢(shì),位于研究區(qū)東部的開通鎮(zhèn)淺層地下水氟濃度整體水平(32mg/L),高于其他三個(gè)鄉(xiāng)鎮(zhèn),這可能與開通鎮(zhèn)平緩低洼的地勢(shì)和活躍的人類活動(dòng)有關(guān),在今后的防氟治理工作中應(yīng)重點(diǎn)關(guān)注此類地區(qū)淺層地下水氟濃度的變化狀況.

    深層高氟地下水分布在研究區(qū)東北部,歷史濃度變化不大,為2~3mg/L;但空間分布范圍逐年擴(kuò)大,可能淺層地下水高氟地下水(中部最大值6.24mg/L)經(jīng)破損井管道滲漏進(jìn)入承壓含水層(中東部最大3.03mg/L),導(dǎo)致深層地下水氟濃度升高.

    (2)地下水砷濃度動(dòng)態(tài)變化

    淺層高砷地下水在研究區(qū)中部、北部的地勢(shì)低洼帶富集,結(jié)合研究區(qū)地下水砷濃度文獻(xiàn)資料[10,23-25,37],在研究區(qū)中部與北部選取烏蘭花鎮(zhèn)-新華鎮(zhèn)-開通鎮(zhèn)-四井子鄉(xiāng)的四個(gè)代表不同水文地質(zhì)環(huán)境的行政分區(qū),繪制各分區(qū)2006年、2010年、2013年、2014年、2021年的地下水砷濃度浮動(dòng)柱狀圖(圖15),從圖15可知:

    ①2006年至2021年各鄉(xiāng)鎮(zhèn)淺層地下水砷濃度整體水平遠(yuǎn)高于世衛(wèi)組織飲用水的指導(dǎo)值10μg/L.新華鎮(zhèn)淺層地下水砷濃度整體水平呈逐年下降趨勢(shì),其他三個(gè)鄉(xiāng)鎮(zhèn)地下水砷濃度整體水平未出現(xiàn)明顯變化趨勢(shì).

    ②2006年、2010年、2013年、2014年、2021年各年份淺層地下水砷濃度最高整體水平分別為新華鎮(zhèn)(50~100μg/L)、開通鎮(zhèn)(50~85μg/L)、開通鎮(zhèn)(90~140 μg/L)、四井子鄉(xiāng)(50~100μg/L)、烏蘭花鎮(zhèn)(53~114μg/L),而從地理方位來(lái)說(shuō)自南向北為新華鎮(zhèn)-開通鎮(zhèn)-烏蘭花鎮(zhèn)-四井子鄉(xiāng),說(shuō)明研究區(qū)淺層地下水砷富集中心區(qū)有自南向北部轉(zhuǎn)移的趨勢(shì).這可能是地下水從山前沖積平原向東北部低平原流動(dòng)時(shí),受研究區(qū)北部復(fù)興隆起與向陽(yáng)隆起帶的阻擋,徑流變緩的同時(shí)蒸發(fā)濃縮作用加強(qiáng),導(dǎo)致淺層高砷地下水在研究區(qū)東北部富集.

    深層高砷地下水分布在研究區(qū)中部與東北部,歷史濃度變化不大,為50~150μg/L,深層高砷地下水較淺層分布更廣、富集區(qū)更多、污染更嚴(yán)重,這可能與承壓含水層的還原條件促進(jìn)砷的釋放有關(guān).

    3 結(jié)論

    3.1 2021年7月中旬地下水樣本檢測(cè)結(jié)果顯示,新華牛場(chǎng)、瞻榆鎮(zhèn)第四系淺層氟富集且污染程度高,水中氟濃度最大值為6.2mg/L.深層地下水F-在東北雙崗鎮(zhèn)富集.淺層地下水As在烏蘭花鎮(zhèn)富集,深層地下水As在雙崗鹿場(chǎng)、新興鄉(xiāng)一帶富集.

    3.2 潛水氟含量受到含水介質(zhì)沉淀溶解作用的制約,蒸發(fā)濃縮作用對(duì)大部分淺層高氟與低氟地下水形成都有一定影響,但并不是高氟地下水形成的主導(dǎo)作用.水中鐵、錳氧化物的還原溶解促進(jìn)淺層、深層地下水As的釋放.

    3.3 除開通鎮(zhèn)以外,自2006年當(dāng)?shù)卣扇》婪乃胧┖?研究區(qū)中東部淺層地下水氟濃度整體水平有所降低,可能與當(dāng)?shù)亻_發(fā)地下水所導(dǎo)致的地下水位下降、蒸發(fā)減弱等因素有關(guān).雖然深層高氟地下水歷史濃度變化不大,但分布范圍擴(kuò)大,可能與淺層高氟地下水經(jīng)破損井管道滲漏進(jìn)入承壓含水層有關(guān).

    3.4 2006~2021年各鄉(xiāng)鎮(zhèn)淺層地下水砷濃度整體水平遠(yuǎn)高于世衛(wèi)組織飲用水的指導(dǎo)值,新華鎮(zhèn)淺層地下水砷濃度整體水平呈逐年下降趨勢(shì),其他三個(gè)鄉(xiāng)鎮(zhèn)地下水砷濃度整體水平均未出現(xiàn)明顯變化趨勢(shì).雖然深層高砷地下水歷史濃度變化不大,但較淺層分布更廣、富集區(qū)更多、污染更嚴(yán)重,這可能與承壓含水層的還原條件促進(jìn)砷的釋放有關(guān).

    [1] Sajil K P J. Groundwater fluoride contamination in Coimbatore district: a geochemical characterization, multivariate analysis, and human health risk perspective [J]. Environmental Earth Sciences, 2021,80(6).

    [2] Choubisa S L, Choubisa D, Choubisa A. Fluoride contamination of groundwater and its threat to health of villagers and their domestic animals and agriculture crops in rural Rajasthan, India [J]. Environmental Geochemistry and Health, 2023,45(3):607-628.

    [3] Sunkari E D, Adams S J, Okyere M B, et al. Groundwater fluoride contamination in Ghana and the associated human health risks: Any sustainable mitigation measures to curtail the long term hazards? [J]. Groundwater for Sustainable Development, 2022,16:100715.

    [4] 邢麗娜,郭華明,魏 亮,等.華北平原淺層含氟地下水演化特點(diǎn)及成因 [J]. 地球科學(xué)與環(huán)境學(xué)報(bào), 2012,34(4):57-67. Xing L N, Guo H M, Wei L, et al. Characteristics and causes of shallow fluoridated groundwater evolution in the North China Plain [J]. Journal of Earth Science and Environment, 2012,34(4):57-67.

    [5] Berg M, Tran H C, Nguyen T C, et al. Arsenic contamination of groundwater and drinking water in Vietnam: a human health threat. [J]. Environmental Science & Technology, 2001,35(13).

    [6] Halim M A, Majumder R K, Nessa S A, et al. Hydrogeochemistry and arsenic contamination of groundwater in the Ganges Delta Plain, Bangladesh [J]. Journal of Hazardous Materials, 2009,164(2):1335- 1345.

    [7] Guo Huaming, Wang Yanxin. Geochemical characteristics of shallow groundwater in Datong basin, northwestern China [J]. Journal of Geochemical Exploration, 2005,87(3):109-120.

    [8] 鄧安琪,董兆敏,高 群,等.中國(guó)地下水砷健康風(fēng)險(xiǎn)評(píng)價(jià) [J]. 中國(guó)環(huán)境科學(xué), 2017,37(9):3556-3565. Deng A Q, Dong Z M, Gao Q, et al. Health risk assessment of arsenic in groundwater in China [J]. China Environmental Science, 2017, 37(9):3556-3565.

    [9] 鄧婭敏.河套盆地西部高砷地下水系統(tǒng)中的地球化學(xué)過(guò)程研究 [D]. 武漢:中國(guó)地質(zhì)大學(xué), 2008. Deng Y M. Geochemical processes in the high arsenic groundwater system in the western Hetao Basin [D]. Wuhan: China University of Geosciences, 2008.

    [10] 張 迪,郭華明,倪 萍,等.氧化還原條件對(duì)地下水中砷釋放遷移的影響—以通榆縣高砷地下水為例 [J]. 第四紀(jì)研究, 2014,34(5): 1072-1081. Zhang D, Guo H M, Ni P, et al. Influence of redox conditions on the migration of arsenic release from groundwater-an example of high- arsenic groundwater in Tongyu County [J]. Quaternary Research, 2014,34(5):1072-1081.

    [11] 袁翰卿,李 巧,陶洪飛,等.新疆奎屯河流域地下水砷富集因素 [J]. 環(huán)境化學(xué), 2020,39(2):524-530. Yuan H Q, Li Q, Tao H F, et al. Arsenic enrichment factors in groundwater of Kuitun River Basin, Xinjiang [J]. Environmental Chemistry, 2020,39(2):524-530.

    [12] 劉圣鋒,高 柏,張海陽(yáng),等.海拉爾盆地地下水氟的分布特征及富集機(jī)理 [J]. 干旱區(qū)資源與環(huán)境, 2021,35(10):169-177. Liu S F, Gao B, Zhang H Y, et al. Distribution characteristics and enrichment mechanism of fluorine in groundwater of Hailar Basin [J]. Arid Zone Resources and Environment, 2021,35(10):169-177.

    [13] 劉春華,王 威,楊麗芝,等.山東省地下水氟富集規(guī)律及其驅(qū)動(dòng)機(jī)制 [J]. 地質(zhì)學(xué)報(bào), 2021,95(6):1962-1972. Liu C H, Wang W, Yang L Z, et al. Fluoride enrichment in groundwater and its driving mechanism in Shandong Province [J]. Journal of Geology, 2021,95(6):1962-1972.

    [14] 呂曉立,劉景濤,周 冰,等.塔城盆地地下水氟分布特征及富集機(jī)理 [J]. 地學(xué)前緣, 2021,28(2):426-436. Lv X L, Liu J T, Zhou B, et al. Distribution characteristics and enrichment mechanism of fluorine in groundwater of the Tacheng Basin [J]. Geological Frontiers, 2021,28(2):426-436.

    [15] 李予紅,宋曉光,胡 斌,等.張家口壩上地區(qū)高氟地下水分布與成因分析 [J]. 北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021,57(6):745-755. Li Y H, Song X G, Hu B, et al. Distribution and causes of high fluorine groundwater in Zhangjiakou dam area [J]. Journal of Beijing Normal University (Natural Science Edition), 2021,57(6):745-755.

    [16] 宋 偉.通榆縣防氟改水現(xiàn)狀調(diào)查分析 [J]. 中國(guó)地方病防治雜志, 2018,33(1):42. Song W. Investigation and analysis of the current situation of water fluoridation prevention and water reform in Tongyu County [J]. Chinese Journal of Endemic Disease Control, 2018,33(1):42.

    [17] 張文強(qiáng),董 冬,熊紹禮.吉林省通榆縣地方氟中毒與地質(zhì)環(huán)境關(guān)系的研究 [J]. 吉林地質(zhì), 2008,27(1):40-44. Zhang W Q, Dong D, Xiong S L. Study on the relationship between local fluorosis and geological environment in Tongyu County, Jilin Province [J]. Jilin Geology, 2008,27(1):40-44.

    [18] 馬玉福.吉林省通榆縣高氟地下水的形成條件與氟病防治措施 [J]. 吉林地質(zhì), 1985,(4):67-72. Ma Y F. Formation conditions of highly fluoridated groundwater in Tongyu County, Jilin Province and measures to prevent fluorosis [J]. Jilin Geology, 1985,(4):67-72.

    [19] 李 霄,邸志強(qiáng),柴 璐.吉林省通榆縣地下水化學(xué)動(dòng)態(tài)分析 [J]. 地質(zhì)與資源, 2008,17(4):292-296. Li X, Di Z Q, Chai L. Analysis of groundwater chemistry dynamics in Tongyu County, Jilin Province [J]. Geology and Resources, 2008,17(4): 292-296.

    [20] 卞建民,李立軍,楊 坡.吉林省通榆縣地下水脆弱性研究 [J]. 水資源保護(hù), 2008,(3):4-7. Bian J M, Li L J, Yang P. Groundwater vulnerability study in Tongyu County, Jilin Province [J]. Water Resources Protection, 2008,(3):4-7.

    [21] 吳艷新.通榆縣農(nóng)村飲水安全現(xiàn)狀與防治措施 [J]. 科技創(chuàng)新與應(yīng)用, 2014,(8):113. Wu Y X. Status quo of rural drinking water safety and preventive measures in Tongyu County [J]. Science and technology innovation and application, 2014,(8):113.

    [22] 胡婧敏.松嫩平原地方病嚴(yán)重區(qū)地下水氟的賦存特征及水質(zhì)安全評(píng)價(jià) [D]. 長(zhǎng)春:吉林大學(xué), 2016. Hu J M. Characteristics of fluorine in groundwater and evaluation of water quality and safety in endemic disease serious area of Songnen Plain [D]. Changchun: Jilin University, 2016.

    [23] 卞建民,湯 潔,封 靈,等.吉林西部砷中毒區(qū)水文地球化學(xué)特征 [J]. 水文地質(zhì)工程地質(zhì), 2009,36(5):80-83. Bian J M, Tang J, Feng L, et al. Hydrogeochemical characteristics of arsenic poisoning area in western Jilin [J]. Hydrogeology and Engineering Geology, 2009,36(5):80-83.

    [24] 趙 娟,李育松,卞建民,等.吉林西部地區(qū)高砷地下水砷的閾值分析及風(fēng)險(xiǎn)評(píng)價(jià) [J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2013,43(1):251-258. Zhao J, Li Y S, Bian J M, et al. Threshold analysis and risk evaluation of arsenic in high arsenic groundwater in western Jilin [J]. Journal of Jilin University (Earth Science Edition), 2013,43(1):251-258.

    [25] 湯 潔,卞建民,李昭陽(yáng),等.吉林省飲水型砷中毒區(qū)地下水砷的分布規(guī)律與成因研究 [J]. 地學(xué)前緣, 2014,(4):30-36. Tang J, Bian J M, Li Z Y, et al. Study on the distribution pattern and causes of arsenic in groundwater in the drinking water arsenic poisoning area of Jilin Province [J]. Geological Frontiers, 2014, (4):30-36.

    [26] 孫 龍,劉廷璽,段利民,等.平朔礦區(qū)不同水體水化學(xué)特征及氟分布成因 [J]. 環(huán)境科學(xué), 2022,43(12):5547-5559. Sun L, Liu T X, Duan L M, et al. Hydrochemical characteristics of different water bodies and the causes of fluorine distribution in the Pingshuo mining area [J]. Environmental Science, 2022,43(12):5547- 5559.

    [27] Berg M, Trang P T K, Stengel C, et al. Hydrological and sedimentary controls leading to arsenic contamination of groundwater in the hanoi area, vietnam; the impact of iron-arsenic ratios, peat, river bank deposits, and excessive groundwater abstraction [J]. Chemical Geology, 2008,249(1/2):91-112.

    [28] 王 翔,羅艷麗,鄧雯文,等.新疆奎屯地區(qū)高砷地下水DOM三維熒光特征 [J]. 中國(guó)環(huán)境科學(xué), 2020,40(11):4974-4981. Wang X, Luo Y L, Deng W W, et al. Three-dimensional fluorescence characteristics of DOM in high arsenic groundwater in Kuitun area, Xinjiang [J]. China Environmental Science, 2020,40(11):4974-4981.

    [29] Shakoor M B, Bibi I, Niazi N K, et al. The evaluation of arsenic contamination potential, speciation and hydrogeochemical behaviour in aquifers of Punjab, Pakistan [J]. Chemosphere, 2018,199.

    [30] Rasool A, Xiao Tangfu, Farooqi A, Muhammad Shafeeque, Sajid Masood, Salar Ali, Shah Fahad, Wajid Nasim. Arsenic and heavy metal contaminations in the tube well water of Punjab, Pakistan and risk assessment: A case study [J]. Ecological Engineering, 2017,39(4): 847-863.

    [31] Zhang Di, Guo Huaming, Xiu Wei, et al. In-situ mobilization and transformation of iron oxides-adsorbed arsenate in natural groundwater [J]. Journal of Hazardous Materials, 2017,321:228-237.

    [32] 孫述海,孫亞紅,遲寶明,等.吉林通榆氟中毒區(qū)地下水中氟來(lái)源及防改對(duì)策 [J]. 安徽農(nóng)業(yè)科學(xué), 2012,40(9):5428-5430. Sun S H, Sun Y H, Chi B M, et al. Sources of fluorine in groundwater and countermeasures for prevention and modification in Tongyu fluorosis area of Jilin [J]. Anhui Agricultural Science, 2012,40(9): 5428-5430.

    [33] 賈 璐,梁秀娟,肖長(zhǎng)來(lái),等.白城市淺層高氟地下水分布特征及富集機(jī)理 [J]. 水電能源科學(xué), 2021,39(1):51-54,123. Jia L, Liang X J, Xiao C L, et al. Distribution characteristics and enrichment mechanism of shallow high-fluoride groundwater in Baicheng City [J]. Hydropower Energy Science, 2021,39(1):51- 54,123.

    [34] 仲偉婉.通榆縣開通鎮(zhèn)擬建水源地地下水?dāng)?shù)值模擬及預(yù)測(cè) [D]. 長(zhǎng)春:吉林大學(xué), 2017. Zhong W W. Numerical simulation and prediction of groundwater in the proposed water source in Kaifaian Town, Tongyu County [D]. Changchun:Jilin University, 2017.

    [35] 周艷紅,尹 華,劉玉申.通榆縣水資源現(xiàn)狀及其節(jié)水潛力分析 [J]. 農(nóng)業(yè)與技術(shù), 2014,34(12):39-40. Zhou Y H, Yin H, Liu Y S. Analysis of the current situation of water resources and its water saving potential in Tongyu County [J]. Agriculture and Technology, 2014,34(12):39-40.

    [36] 李明彧.吉林省通榆縣“節(jié)水增糧行動(dòng)”水資源論證研究 [D]. 長(zhǎng)春:吉林大學(xué), 2017. Li M Y. Research on water resources demonstration of "water saving and grain increasing action" in Tongyu County, Jilin Province [D]. Changchun: Jilin University, 2017.

    [37] 邴智武.松嫩平原地下水氟、砷的富集規(guī)律及影響因素研究 [D]. 長(zhǎng)春:吉林大學(xué), 2009. Bing Z W. Study on the enrichment pattern of fluorine and arsenic in groundwater of Songnen Plain and the influencing factors [D]. Changchun: Jilin University, 2009.

    致謝:感謝吉林大學(xué)“凝聚青春力量,志愿服務(wù)賦能鄉(xiāng)村振興”項(xiàng)目對(duì)野外調(diào)研工作的支持!

    Enrichment mechanism and dynamic variation characteristics of fluorine and arsenic in groundwater of western Songnen Plain.

    DENG Yuan-dong1,2, YE Xue-yan1,2, WU Ya-min1,2, DU Xin-qiang1,2*

    (1.College of New Energy and Environment, Jilin University, Changchun 130021, China;2.Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China)., 2023,43(10):5277~5290

    In order to further reveal the the evolutionary trends of fluorine and arsenic in groundwater in the western part of the Songnen Plain, a comprehensive investigation was conducted on the current condition of fluorine and arsenic concentrations in groundwater, focusing on Tongyu County as a case study. Water quality data were analyzed using ion-ratio diagrams and self-organizing mapping networks. The results show that fluorine content of the phreatic water was governed by the precipitation and dissolution of the water-bearing medium, and evaporation and concentration had a certain influence on the formation of most of the shallow, high and low fluoride groundwater, but were not the dominant role in the formation of high fluoride groundwater. The dissolution of iron and manganese oxides from water promoted the release of As from shallow and deep groundwater. The overall level of the fluoride concentration in shallow groundwater in the mid-eastern part of the region has decreased since 1995, which may be attributed to factors such as the depletion of groundwater levels and decreased evaporation resulting from extensive local groundwater extraction. Although the deep, high fluoride groundwater did not show a significant change, its spatial coverage has expanded; This may be related to the leakage of shallow, high fluoride groundwater into the pressurized aquifers through the broken well pipes. The overall concentrations of arsenic in shallow groundwater in the central and northeastern regions from 1995to 2021 significantly exceed the WHO guideline values for drinking water. Although the concentrations of high-arsenic groundwater in deeper layers have remained relatively stable, they exhibit a wider distribution, more enriched zones, and greater contamination compared to high arsenic zone in shallow groundwater. This may be related to the reducing conditions in the pressurized aquifer that promoted the release of arsenic. Therefore, the exploitation of deep groundwater in Tongyu County should be conducted with cautions, taking into consideration of the areas affected by high levels of arsenic contamination. The exploitation level should be carefully determined to prevent the intrusion of polluted water. Additionally, attentions should be given to address the potentially cascading contaminations resulting from the leakage and destruction of wells.

    high-fluoride groundwater;high-arsenic groundwater;enrichment mechanism;change of dynamics

    X52

    A

    1000-6923(2023)10-5277-14

    2023-03-06

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41972247)

    * 責(zé)任作者, 教授, duxq@jlu.edu.cn

    鄧遠(yuǎn)東(1998-),男,湖南衡陽(yáng)人,吉林大學(xué)博士研究生,主要從事地下水資源與環(huán)境的研究工作.發(fā)表論文2篇.1875235909@qq.com.

    鄧遠(yuǎn)東,冶雪艷,吳亞敏,等.松嫩平原西部地下水氟和砷的富集機(jī)理與動(dòng)態(tài)變化特征 [J]. 中國(guó)環(huán)境科學(xué), 2023,43(10):5277-5290.

    Deng Y D, Ye X Y, Wu Y M, et al. Enrichment mechanism and dynamic variation characteristics of fluorine and arsenic in groundwater of western Songnen Plain [J]. China Environmental Science, 2023,43(10):5277-5290.

    猜你喜歡
    通榆縣淺層深層
    白城市通榆縣召開“十三五”教育信息化工作表彰暨教育2.0大數(shù)據(jù)平臺(tái)啟動(dòng)會(huì)
    淺層換填技術(shù)在深厚軟土路基中的應(yīng)用
    基于淺層曝氣原理的好氧顆粒污泥的快速培養(yǎng)
    考慮各向異性滲流的重力壩深層抗滑穩(wěn)定分析
    SAM系統(tǒng)對(duì)TDCS數(shù)據(jù)的優(yōu)化處理與深層應(yīng)用
    依靠農(nóng)業(yè)科技進(jìn)步推動(dòng)通榆縣農(nóng)業(yè)經(jīng)濟(jì)發(fā)展
    淺層地下水超采區(qū)劃分探究
    風(fēng)力發(fā)電發(fā)展SWOT分析
    包氣帶淺層地?zé)崛萘坑?jì)算方法商榷
    關(guān)于通榆縣風(fēng)力發(fā)電的調(diào)查分析報(bào)告
    天堂动漫精品| 激情在线观看视频在线高清| 国产精品亚洲一级av第二区| 麻豆久久精品国产亚洲av| 好男人电影高清在线观看| 2021天堂中文幕一二区在线观| 97人妻精品一区二区三区麻豆| 九九久久精品国产亚洲av麻豆| 高清日韩中文字幕在线| 亚洲内射少妇av| 免费大片18禁| 中亚洲国语对白在线视频| 日本一二三区视频观看| 亚洲成人久久爱视频| 99热6这里只有精品| 成人亚洲精品av一区二区| 国产综合懂色| 看黄色毛片网站| 一级毛片女人18水好多| 免费av观看视频| 欧美精品啪啪一区二区三区| 在线a可以看的网站| 国产99白浆流出| 99国产极品粉嫩在线观看| 美女免费视频网站| 国产精品免费一区二区三区在线| 欧美一区二区国产精品久久精品| 亚洲欧美一区二区三区黑人| 一边摸一边抽搐一进一小说| 在线观看日韩欧美| 日韩av在线大香蕉| 91在线精品国自产拍蜜月 | 18+在线观看网站| 天天一区二区日本电影三级| 九九在线视频观看精品| 亚洲中文字幕一区二区三区有码在线看| 午夜福利在线在线| 热99在线观看视频| 九色成人免费人妻av| 午夜亚洲福利在线播放| 最近最新免费中文字幕在线| 在线观看午夜福利视频| 一个人免费在线观看电影| 在线观看午夜福利视频| 成人性生交大片免费视频hd| 女人十人毛片免费观看3o分钟| 91在线精品国自产拍蜜月 | 在线播放国产精品三级| 精品一区二区三区av网在线观看| 天堂av国产一区二区熟女人妻| a级一级毛片免费在线观看| 欧美日韩国产亚洲二区| 国产精品98久久久久久宅男小说| 999久久久精品免费观看国产| 久久欧美精品欧美久久欧美| 国产精品香港三级国产av潘金莲| 国产精品一区二区三区四区久久| 听说在线观看完整版免费高清| 99久久无色码亚洲精品果冻| tocl精华| 夜夜夜夜夜久久久久| 午夜福利成人在线免费观看| 精品一区二区三区av网在线观看| 久久精品亚洲精品国产色婷小说| 亚洲成人免费电影在线观看| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av香蕉五月| 国产精品日韩av在线免费观看| 欧美zozozo另类| 在线观看日韩欧美| 亚洲男人的天堂狠狠| 尤物成人国产欧美一区二区三区| 午夜日韩欧美国产| 国产黄片美女视频| 国产精品影院久久| 丰满的人妻完整版| 国产又黄又爽又无遮挡在线| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区av网在线观看| 99国产精品一区二区三区| 精品一区二区三区人妻视频| av国产免费在线观看| 天堂√8在线中文| 亚洲国产精品成人综合色| 国产成人aa在线观看| svipshipincom国产片| 18+在线观看网站| 九九在线视频观看精品| 欧美日本亚洲视频在线播放| 日韩欧美三级三区| 深爱激情五月婷婷| 禁无遮挡网站| 久久久国产成人精品二区| 国产色爽女视频免费观看| 国产成+人综合+亚洲专区| 午夜两性在线视频| 午夜福利18| 十八禁网站免费在线| 欧美激情在线99| 国产成人av教育| 国产aⅴ精品一区二区三区波| 韩国av一区二区三区四区| 精品国产超薄肉色丝袜足j| 午夜福利18| 少妇的逼好多水| 色播亚洲综合网| 免费在线观看亚洲国产| 99热这里只有是精品50| 亚洲国产精品合色在线| 日韩精品青青久久久久久| 一进一出好大好爽视频| av在线天堂中文字幕| 亚洲av免费高清在线观看| 欧美日韩亚洲国产一区二区在线观看| 午夜视频国产福利| 国产精品亚洲一级av第二区| 美女大奶头视频| 久久精品综合一区二区三区| 十八禁人妻一区二区| 久久草成人影院| 可以在线观看毛片的网站| 一个人看的www免费观看视频| 亚洲精品一区av在线观看| 亚洲精品美女久久久久99蜜臀| 免费搜索国产男女视频| 97超级碰碰碰精品色视频在线观看| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区免费观看 | 在线国产一区二区在线| 日本黄色视频三级网站网址| 搡老岳熟女国产| 一个人免费在线观看的高清视频| 久久性视频一级片| 欧美极品一区二区三区四区| 丁香六月欧美| 99久久九九国产精品国产免费| 99国产极品粉嫩在线观看| 国产精品一区二区三区四区久久| 一个人观看的视频www高清免费观看| 99久久精品国产亚洲精品| 啦啦啦观看免费观看视频高清| 国语自产精品视频在线第100页| 最好的美女福利视频网| 麻豆国产97在线/欧美| 99riav亚洲国产免费| 人妻久久中文字幕网| 欧美日本视频| 亚洲成人久久性| 18禁国产床啪视频网站| 91久久精品国产一区二区成人 | 国产熟女xx| 国产极品精品免费视频能看的| 久久婷婷人人爽人人干人人爱| 我要搜黄色片| 非洲黑人性xxxx精品又粗又长| 午夜精品一区二区三区免费看| 一级黄色大片毛片| 国产不卡一卡二| 国产精品精品国产色婷婷| 日韩有码中文字幕| 久久国产精品人妻蜜桃| 国产精品99久久99久久久不卡| 精品久久久久久,| 中文字幕熟女人妻在线| 麻豆成人午夜福利视频| 真实男女啪啪啪动态图| 三级男女做爰猛烈吃奶摸视频| 国产一区在线观看成人免费| 香蕉丝袜av| 免费在线观看亚洲国产| 日日摸夜夜添夜夜添小说| 色哟哟哟哟哟哟| 97超视频在线观看视频| 免费看十八禁软件| svipshipincom国产片| 757午夜福利合集在线观看| 精品久久久久久,| 最近视频中文字幕2019在线8| 国产成人aa在线观看| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久电影 | 国内久久婷婷六月综合欲色啪| 岛国在线观看网站| 搡老岳熟女国产| 亚洲av一区综合| xxxwww97欧美| 日韩欧美在线乱码| 成人国产一区最新在线观看| 床上黄色一级片| 亚洲最大成人手机在线| 最近最新中文字幕大全电影3| 国产精品一区二区三区四区免费观看 | 久久精品人妻少妇| 婷婷六月久久综合丁香| 黄片小视频在线播放| 少妇的逼好多水| 嫩草影院精品99| 亚洲七黄色美女视频| 嫩草影院入口| 一本久久中文字幕| 99在线视频只有这里精品首页| 日韩有码中文字幕| xxxwww97欧美| 老司机福利观看| 丰满人妻熟妇乱又伦精品不卡| 九色国产91popny在线| 啦啦啦观看免费观看视频高清| 一区二区三区国产精品乱码| 欧美日本视频| 欧美日韩瑟瑟在线播放| h日本视频在线播放| 在线天堂最新版资源| 国产爱豆传媒在线观看| 亚洲精品日韩av片在线观看 | 国产久久久一区二区三区| 天天一区二区日本电影三级| 亚洲精品日韩av片在线观看 | 久久伊人香网站| 国产精品99久久99久久久不卡| 久久亚洲精品不卡| 国产精品 国内视频| 在线播放国产精品三级| 日本成人三级电影网站| 18禁国产床啪视频网站| 禁无遮挡网站| 中文字幕久久专区| 国产爱豆传媒在线观看| 丰满乱子伦码专区| aaaaa片日本免费| 亚洲一区二区三区不卡视频| 在线视频色国产色| 狠狠狠狠99中文字幕| 亚洲欧美精品综合久久99| 欧美激情久久久久久爽电影| av欧美777| 中国美女看黄片| 国产综合懂色| 国产一区二区三区视频了| 少妇的逼好多水| 狂野欧美白嫩少妇大欣赏| 亚洲午夜理论影院| 少妇人妻一区二区三区视频| 亚洲欧美日韩卡通动漫| 欧美国产日韩亚洲一区| 99视频精品全部免费 在线| 亚洲精品色激情综合| 欧美色视频一区免费| 国产伦精品一区二区三区四那| 亚洲精品亚洲一区二区| 亚洲无线在线观看| 日本黄色片子视频| 国产男靠女视频免费网站| 久久久精品欧美日韩精品| 成人鲁丝片一二三区免费| www日本在线高清视频| 国内揄拍国产精品人妻在线| 亚洲五月天丁香| 最近最新中文字幕大全电影3| 国内精品美女久久久久久| 久久久久久国产a免费观看| 老司机午夜十八禁免费视频| 内射极品少妇av片p| 搡老岳熟女国产| 欧美成人a在线观看| 最新美女视频免费是黄的| 国产免费一级a男人的天堂| 久久久国产精品麻豆| 亚洲人成电影免费在线| 特大巨黑吊av在线直播| 啦啦啦韩国在线观看视频| 9191精品国产免费久久| 国产精品99久久99久久久不卡| 国产高清视频在线播放一区| 啦啦啦免费观看视频1| 亚洲18禁久久av| 国产精品女同一区二区软件 | 亚洲在线自拍视频| 成人永久免费在线观看视频| 全区人妻精品视频| 欧美成人一区二区免费高清观看| 国产成人av教育| 法律面前人人平等表现在哪些方面| 免费人成在线观看视频色| 国产精品乱码一区二三区的特点| 色哟哟哟哟哟哟| 日韩 欧美 亚洲 中文字幕| 老汉色av国产亚洲站长工具| 欧美色视频一区免费| 人人妻人人看人人澡| 欧美av亚洲av综合av国产av| 欧美成人a在线观看| 亚洲av成人精品一区久久| 国产乱人视频| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 在线观看66精品国产| 一二三四社区在线视频社区8| 中文字幕精品亚洲无线码一区| 熟妇人妻久久中文字幕3abv| 国产精品av视频在线免费观看| 国产亚洲精品久久久com| 波多野结衣巨乳人妻| 国产成人系列免费观看| 欧美xxxx黑人xx丫x性爽| 看免费av毛片| 国产野战对白在线观看| 一个人观看的视频www高清免费观看| 国产 一区 欧美 日韩| 国内少妇人妻偷人精品xxx网站| 俄罗斯特黄特色一大片| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| a级毛片a级免费在线| 亚洲av不卡在线观看| 一级黄色大片毛片| 免费高清视频大片| 色视频www国产| 久久国产精品影院| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| av天堂中文字幕网| 麻豆国产97在线/欧美| 成年女人毛片免费观看观看9| 欧美日韩综合久久久久久 | 狂野欧美激情性xxxx| 精品不卡国产一区二区三区| 国产一区二区激情短视频| 蜜桃亚洲精品一区二区三区| 中文亚洲av片在线观看爽| 成人特级av手机在线观看| 欧美成人免费av一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲狠狠婷婷综合久久图片| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 亚洲欧美一区二区三区黑人| 九色国产91popny在线| 国产av在哪里看| 中国美女看黄片| 亚洲av中文字字幕乱码综合| 欧美在线黄色| 真人一进一出gif抽搐免费| 欧美乱色亚洲激情| 丰满人妻一区二区三区视频av | 又黄又粗又硬又大视频| 宅男免费午夜| 国产成人av教育| 精品福利观看| 亚洲精品国产精品久久久不卡| 老司机福利观看| 我的老师免费观看完整版| 观看美女的网站| 国产色爽女视频免费观看| 人人妻人人澡欧美一区二区| 久久久久国内视频| 中文字幕高清在线视频| 最新中文字幕久久久久| 亚洲成av人片免费观看| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 国产视频内射| 国产精品久久久久久久电影 | 韩国av一区二区三区四区| 香蕉av资源在线| 一级黄片播放器| 在线a可以看的网站| 综合色av麻豆| 国产在视频线在精品| 精品人妻1区二区| a在线观看视频网站| 搡老熟女国产l中国老女人| 国产伦精品一区二区三区视频9 | 看黄色毛片网站| 久久久成人免费电影| x7x7x7水蜜桃| 色老头精品视频在线观看| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 欧美三级亚洲精品| 超碰av人人做人人爽久久 | 国产免费一级a男人的天堂| 男女之事视频高清在线观看| 欧美绝顶高潮抽搐喷水| 97超视频在线观看视频| 首页视频小说图片口味搜索| 精品一区二区三区视频在线观看免费| 欧美一区二区国产精品久久精品| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| 国产精品 国内视频| 亚洲国产欧美网| 久久中文看片网| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 国产黄色小视频在线观看| 国产精品电影一区二区三区| 亚洲美女视频黄频| xxxwww97欧美| 精品一区二区三区人妻视频| 国产精品久久久久久人妻精品电影| 午夜免费男女啪啪视频观看 | 性色avwww在线观看| 特级一级黄色大片| 国产一区二区在线av高清观看| 高清在线国产一区| 无遮挡黄片免费观看| 男人和女人高潮做爰伦理| 免费看光身美女| 亚洲国产欧美人成| 尤物成人国产欧美一区二区三区| 国产男靠女视频免费网站| 亚洲成av人片免费观看| 身体一侧抽搐| 国产一区二区亚洲精品在线观看| 99久久成人亚洲精品观看| 亚洲国产精品sss在线观看| 在线观看日韩欧美| 欧美日韩综合久久久久久 | 国产激情欧美一区二区| 99久久精品热视频| 久久久国产成人精品二区| 久久久久久九九精品二区国产| 国产精品 国内视频| 国产成+人综合+亚洲专区| 哪里可以看免费的av片| 久久久久精品国产欧美久久久| 人人妻人人澡欧美一区二区| 国产亚洲av嫩草精品影院| 国产精品国产高清国产av| 欧美日韩中文字幕国产精品一区二区三区| 99热这里只有精品一区| 久久欧美精品欧美久久欧美| 欧美成狂野欧美在线观看| 手机成人av网站| 极品教师在线免费播放| 99热精品在线国产| 久久婷婷人人爽人人干人人爱| 又爽又黄无遮挡网站| 久久精品人妻少妇| 成年人黄色毛片网站| 波野结衣二区三区在线 | 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站| 十八禁人妻一区二区| 激情在线观看视频在线高清| 99国产综合亚洲精品| 国产av麻豆久久久久久久| 国产精品久久久久久精品电影| 中出人妻视频一区二区| av中文乱码字幕在线| 精品熟女少妇八av免费久了| 麻豆成人午夜福利视频| 小说图片视频综合网站| 国产高清有码在线观看视频| 欧美日本视频| 国产三级黄色录像| 亚洲在线观看片| 午夜免费观看网址| 国产成人啪精品午夜网站| 亚洲国产精品999在线| aaaaa片日本免费| 在线观看av片永久免费下载| 亚洲av免费高清在线观看| 在线a可以看的网站| 亚洲黑人精品在线| 伊人久久精品亚洲午夜| 欧美性猛交╳xxx乱大交人| 免费一级毛片在线播放高清视频| 亚洲精品色激情综合| 精品久久久久久久久久久久久| 日韩 欧美 亚洲 中文字幕| 欧美丝袜亚洲另类 | 内地一区二区视频在线| 美女大奶头视频| 日本黄色片子视频| 午夜福利在线在线| 日韩有码中文字幕| 人妻夜夜爽99麻豆av| 男人的好看免费观看在线视频| 婷婷精品国产亚洲av| av视频在线观看入口| 国产精品99久久久久久久久| 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 国产成年人精品一区二区| 午夜久久久久精精品| 国产蜜桃级精品一区二区三区| 少妇的丰满在线观看| 国产三级中文精品| 操出白浆在线播放| 老汉色∧v一级毛片| 成年人黄色毛片网站| 动漫黄色视频在线观看| 午夜免费激情av| 亚洲人成网站在线播放欧美日韩| 国产精品久久视频播放| 国产精品99久久99久久久不卡| 亚洲欧美一区二区三区黑人| 无限看片的www在线观看| 舔av片在线| 小蜜桃在线观看免费完整版高清| 九色成人免费人妻av| 免费在线观看日本一区| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 在线免费观看的www视频| 欧美大码av| 一进一出好大好爽视频| 色尼玛亚洲综合影院| 久久99热这里只有精品18| svipshipincom国产片| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 久久精品国产自在天天线| 人妻久久中文字幕网| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 男女下面进入的视频免费午夜| 一卡2卡三卡四卡精品乱码亚洲| 一a级毛片在线观看| 亚洲国产中文字幕在线视频| 亚洲片人在线观看| 亚洲不卡免费看| 一区福利在线观看| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 中亚洲国语对白在线视频| 亚洲无线观看免费| 男女床上黄色一级片免费看| or卡值多少钱| 18+在线观看网站| 黄色视频,在线免费观看| 久久久久免费精品人妻一区二区| 国产精品女同一区二区软件 | 久久九九热精品免费| 熟女电影av网| 久久久色成人| 日韩免费av在线播放| 三级毛片av免费| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| www日本黄色视频网| 法律面前人人平等表现在哪些方面| 黄色女人牲交| 国产又黄又爽又无遮挡在线| 五月玫瑰六月丁香| 国产激情欧美一区二区| 丰满乱子伦码专区| 好看av亚洲va欧美ⅴa在| 亚洲av第一区精品v没综合| 激情在线观看视频在线高清| 18美女黄网站色大片免费观看| 日韩中文字幕欧美一区二区| 亚洲欧美日韩高清专用| 大型黄色视频在线免费观看| 国产精品美女特级片免费视频播放器| 亚洲精品亚洲一区二区| 夜夜夜夜夜久久久久| 五月伊人婷婷丁香| 在线播放国产精品三级| 丁香欧美五月| 日本在线视频免费播放| www日本黄色视频网| 最近在线观看免费完整版| 好男人电影高清在线观看| 97超级碰碰碰精品色视频在线观看| 国产成人系列免费观看| 国产成人啪精品午夜网站| 久久中文看片网| 可以在线观看毛片的网站| 国产精品精品国产色婷婷| 757午夜福利合集在线观看| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区视频了| 午夜激情福利司机影院| 97人妻精品一区二区三区麻豆| 欧洲精品卡2卡3卡4卡5卡区| 欧美黑人欧美精品刺激| 小蜜桃在线观看免费完整版高清| 午夜亚洲福利在线播放| 变态另类成人亚洲欧美熟女| 欧美乱妇无乱码| 亚洲精品亚洲一区二区| 成年女人毛片免费观看观看9| 男女做爰动态图高潮gif福利片| 久久人妻av系列| 国产三级黄色录像| 欧美乱色亚洲激情| 欧美成人a在线观看| 成人特级黄色片久久久久久久| 精品熟女少妇八av免费久了| 午夜影院日韩av| 久久伊人香网站| 亚洲精品美女久久久久99蜜臀| 在线播放无遮挡| 亚洲黑人精品在线| 999久久久精品免费观看国产| 非洲黑人性xxxx精品又粗又长| 久久精品国产自在天天线| 真实男女啪啪啪动态图| 成年女人永久免费观看视频| 无限看片的www在线观看| 久久天躁狠狠躁夜夜2o2o| 12—13女人毛片做爰片一| 成人鲁丝片一二三区免费| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 国产伦人伦偷精品视频| 国产精品亚洲一级av第二区| bbb黄色大片| 色噜噜av男人的天堂激情| 天堂网av新在线| 免费看美女性在线毛片视频| 在线观看66精品国产| 亚洲人与动物交配视频|