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

    珠江三角洲典型水產(chǎn)養(yǎng)殖區(qū)淺層地下水碘化物分布、來(lái)源及健康風(fēng)險(xiǎn)評(píng)估

    2024-01-01 00:00:00段磊曾經(jīng)趙顯林丘錦榮劉娜陶鈞實(shí)周建利
    關(guān)鍵詞:研究

    摘要:為了解珠江三角洲典型水產(chǎn)養(yǎng)殖區(qū)淺層地下水中碘化物的來(lái)源以及存在的飲用健康風(fēng)險(xiǎn),本研究綜合運(yùn)用數(shù)理統(tǒng)計(jì)分析、主成分分析、水質(zhì)評(píng)價(jià)和健康風(fēng)險(xiǎn)評(píng)估等方法對(duì)珠江三角洲典型淡水養(yǎng)殖區(qū)域內(nèi)魚塘、污水處理池、居民飲用井三種類型水源共21個(gè)代表性樣品進(jìn)行分析。結(jié)果表明:研究區(qū)淺層地下水水化學(xué)類型以Ca2+·HCO3-型為主,ρ(I-)介于2~343 μg/L之間,其中33.3%的監(jiān)測(cè)點(diǎn)位為高碘型水源,集中分布于研究區(qū)西部;淺層地下水中的碘易賦存于弱堿性還原環(huán)境中,而魚塘泥沙中的富碘有機(jī)質(zhì)的分解滲透作用可能會(huì)加重地下水的碘富集;城鎮(zhèn)化進(jìn)程中產(chǎn)生的還原性廢水和垃圾滲濾液伴隨著富含碳酸鹽巖體中有機(jī)質(zhì)的分解,都可能成為淺層高碘地下水的重要來(lái)源;研究區(qū)地表水(魚塘水、污水處理池水)中Ⅳ類水占比為16.7%,處于較差類別,其中化學(xué)需氧量(CODMn)與總氮(TN)為主要超標(biāo)因子;淺層地下水質(zhì)量類別最低為Ⅲ類水,整體水質(zhì)較好,部分點(diǎn)位存在pH和ρ(I-)超標(biāo)問題;成人的高碘飲用水健康風(fēng)險(xiǎn)較低,J4點(diǎn)位地下水源的兒童非致癌風(fēng)險(xiǎn)值大于1.0,建議加強(qiáng)對(duì)居民引用水井J4的ρ(I-)的常規(guī)監(jiān)控,以保障兒童飲用水健康。

    關(guān)鍵詞:珠江三角洲;水質(zhì)評(píng)價(jià);碘化物;健康風(fēng)險(xiǎn);淺層地下水;地表水

    doi:10.13278/j.cnki.jjuese.20230022

    中圖分類號(hào):X523;X820

    文獻(xiàn)標(biāo)志碼:A

    段磊,曾經(jīng)文,趙顯林,等.珠江三角洲典型水產(chǎn)養(yǎng)殖區(qū)淺層地下水碘化物分布、來(lái)源及健康風(fēng)險(xiǎn)評(píng)估.吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2024,54(5):16571674. doi:10.13278/j.cnki.jjuese.20230022.

    Duan Lei,Zeng Jingwen,Zhao Xianlin,et al. Distribution, Sources and Health Risk Assessment of Iodide in Shallow Groundwater in Typical Aquaculture Areas of the Pearl River Delta. Journal of Jilin University (Earth Science Edition),2024,54(5):16571674. doi:10.13278/j.cnki.jjuese.20230022.

    收稿日期:20230210

    作者簡(jiǎn)介:段磊(1996-),男,碩士研究生,主要從事農(nóng)業(yè)水土資源利用方面的研究,E-mail:936263632@qq.com

    通信作者:周建利(1971-),男,副教授,博士,主要從事水土環(huán)境方面的研究,E-mail:631334055@qq.com

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(41572217)

    Supported by the National Natural Science Foundation of China (41572217)

    Distribution, Sources and Health Risk Assessment of Iodide in Shallow Groundwater in Typical Aquaculture Areas of the Pearl River Delta

    Duan Lei 1, Zeng Jingwen2, Zhao Xianlin1, Qiu Jinrong2, Liu Na3, Tao Junshi2, Zhou Jianli1

    1." College of Agriculture, Yangtze University, Jingzhou 434000, Hubei, China

    2." South China Institute of Environmental Science, Ministry of Ecology and Environment, Guangzhou 510630, China

    3." College of Life Science and Technology, Jinan University, Guangzhou 510630, China

    Abstract:

    In order" to understand the sources of iodide in shallow groundwater in these typical aquaculture areas of the Pearl River delta as well as the health risks associated with drinking.

    This study analyzed a total of 21 representative samples from three different types of water sources: fish ponds, sewage treatment ponds, and groundwater in typical freshwater aquaculture areas of the Pearl River" delta.

    The results indicated that shallow groundwater in the study area had a water chemistry type primarily composed of Ca2+·HCO3- type with iodide concentrations ranging from 2343 μg/L of the monitored sites. 33.3% of the monitored sites were high iodine type water sources that were concentrated and dispersed in the western part of the study area. According to principal component analysis, iodine fugacity in groundwater was encouraged by the neutral to mildly alkaline reducing environment, and the degradation and percolation of iodine-rich organic materials from fish pond waste may have made the situation worse. The major causes of shallow, highly iodized groundwater may include degradation of organic matter in carbonate-rich rocks, urbanization with reduced effluent, and waste leachate leaks. The results of water quality evaluation showed that 16.7% of the surface water (fish ponds, sewage treatment ponds) in the study area was Class Ⅳ, and the chemical oxygen demand (CODMn) and total nitrogen (TN) were the main factors exceeding the standard. Although the pH and iodide levels in some areas of shallow groundwater exceed the recommended levels, overall water quality is satisfactory, with the worst groundwater quality being Grade Ⅲ. The examination of the health risks associated with groundwater revealed that while adults are at minimal risk from high iodine drinking water, children have a risk entropy of J4 groundwater sources that is larger than 1.0. It is suggested to strengthen the routine monitoring of residential drinking well J4 iodide concentration to ensure the health of children drinking water.

    Key words:

    Pearl River delta; water quality evaluation; iodide; health risk;shallow groundwater;surface water

    0" 引言

    近年來(lái)隨著經(jīng)濟(jì)高速發(fā)展,人們對(duì)地下水的開發(fā)利用愈加重視。我國(guó)地下水資源的空間分布特征是:北部水資源匱乏,南部水資源豐富,全國(guó)80%的水資源分布在南方;南方除山區(qū)較小的盆地有少量的地下水開采外,其余地下水開采活動(dòng)主要集中于人口稠密的城鎮(zhèn),在一些偏遠(yuǎn)的村莊,生活用水也都來(lái)自淺層地下水[13]。有研究[2,4]表明,碘傾向于在淺層地下水中富集,而淺層地下水容易受到工業(yè)點(diǎn)源污染、農(nóng)業(yè)面源污染、集約化養(yǎng)殖污染、居民生活污水和固廢(生活垃圾)污染等方面的影響。近年來(lái),隨著城鎮(zhèn)化進(jìn)程不斷加快,工業(yè)和農(nóng)業(yè)生產(chǎn)中大量含碘廢水和生活污水的排放對(duì)淺層地下水環(huán)境造成很大壓力,其中碘伏消毒劑在漁業(yè)養(yǎng)殖中通過溶劑或片劑等方式進(jìn)入水體環(huán)境中而對(duì)人體健康和水生態(tài)環(huán)境造成潛在威脅。因此,研究淺層地下水中碘的分布和來(lái)源已逐漸受到了國(guó)內(nèi)外學(xué)者的廣泛重視。除我國(guó)[57]外世界上很多國(guó)家存在高碘地下水(ρ(I-)>100 μg/L),如丹麥[8]、摩洛哥[9]、智利[10]、阿富汗[11]、俄羅斯[12]等。據(jù)文獻(xiàn)[13]報(bào)道,居民飲用水中ρ(I-)<10 μg/L為水源性缺碘地區(qū),ρ(I-)>100 μg/L為水源性高碘地區(qū)。碘作為最重要的微量元素之一,在能量代謝、體溫調(diào)節(jié)、身體和智力發(fā)育中都起著至關(guān)重要的作用,碘攝入不足或過量都將導(dǎo)致較高的患病機(jī)率[14],例如孕婦孕期碘攝入不足可增加死產(chǎn)、流產(chǎn)和嬰兒患癡呆癥的風(fēng)險(xiǎn),而當(dāng)?shù)鈹z入過量時(shí),可誘發(fā)甲狀腺功能亢進(jìn)、甲狀腺炎或智力喪失[1517]。自1980年以來(lái),沿海地區(qū)的許多甲狀腺腫流行都與飲用水中的碘過量有關(guān)[4]。世界衛(wèi)生組織(WHO)建議0~5歲兒童碘的攝入量為90 μg/d,6~12歲兒童碘的攝入量為120 μg/d,成年人碘的攝入量為150 μg /d,妊娠和哺乳期碘的攝入量為250 μg/d [18]。

    地下水中碘以碘化物(I-)、碘酸鹽(IO3-)和有機(jī)碘(OI)的形式存在[5]。地下水系統(tǒng)的有機(jī)組分、微生物活動(dòng)、pH值以及氧化還原環(huán)境等因素都會(huì)影響碘的遷移和轉(zhuǎn)化[19]。有研究[7, 2021]表明,pH值和氧化還原電位(Eh)作為碘地球化學(xué)行為的兩個(gè)重要參數(shù),共同揭示了碘在弱堿性還原條件下的地下水中更加易于富集的規(guī)律。另外,當(dāng)?shù)叵滤h(huán)境為氧化狀態(tài)時(shí),有機(jī)質(zhì)可以通過氧化作用將有機(jī)結(jié)合態(tài)的碘釋放到地下水中;而地下水在還原環(huán)境中又可誘發(fā)富碘有機(jī)質(zhì)的降解,富碘有機(jī)質(zhì)的氧化降解成為高碘地下水演變的重要水化學(xué)反應(yīng)之一[22]。同時(shí),氣候類型、沉積環(huán)境和地下徑流條件都影響著高碘地下水的賦存;此外,城鎮(zhèn)化進(jìn)程中所產(chǎn)生的污水和垃圾滲濾液的輸入也成為地下水中碘化物升高的重要影響因素[23]。Zhang等[24]在揭示華北滄州地區(qū)深層地下水中碘的富集機(jī)理時(shí)表明,地質(zhì)構(gòu)造差異產(chǎn)生的不同地下徑流條件是影響碘富集的主要因素;而呂曉立等[25]研究表明,珠江三角洲城鎮(zhèn)化過程中還原性垃圾滲濾液和污水輸入是高碘地下水的驅(qū)動(dòng)因素。

    目前,珠江三角洲地下水環(huán)境質(zhì)量評(píng)價(jià)顯示,該區(qū)地下水環(huán)境受到一定程度的有機(jī)污染[2627],氨氮、亞硝酸鹽、鐵、錳等重金屬存在普遍超標(biāo)現(xiàn)象[1516,20,28],區(qū)域地下水酸化(pH<6.5)較為普遍,“三氮”污染問題突出[17];但是,地下水中碘離子的質(zhì)量濃度尚未納入常規(guī)監(jiān)測(cè)指標(biāo),有關(guān)水源地碘離子的質(zhì)量濃度未引起公眾足夠重視。近年來(lái),有關(guān)我國(guó)地下水碘化物來(lái)源的研究多集中于我國(guó)黃淮海平原和內(nèi)陸盆地等區(qū)域[2931],但對(duì)珠江三角洲地區(qū),特別是在城鎮(zhèn)化和水產(chǎn)養(yǎng)殖活動(dòng)影響下淺層地下水中碘化物的分布、來(lái)源、水質(zhì)現(xiàn)狀及存在的健康風(fēng)險(xiǎn)關(guān)注較少。因此,本研究選擇珠江三角地區(qū)典型水產(chǎn)養(yǎng)殖區(qū)域內(nèi)魚塘水、污水處理池尾水和居民飲用井水為研究對(duì)象,對(duì)不同類型水源的水化學(xué)指標(biāo)進(jìn)行檢測(cè),分析其中碘化物的空間分布特征和來(lái)源,并對(duì)水質(zhì)進(jìn)行現(xiàn)狀分析以及對(duì)當(dāng)?shù)鼐用翊嬖诘娘嬘媒】碉L(fēng)險(xiǎn)進(jìn)行評(píng)估,以期為完善我國(guó)珠江三角洲典型水產(chǎn)養(yǎng)殖區(qū)地下水管理和農(nóng)村居民合理攝碘量提供科學(xué)依據(jù)。

    1" 材料與方法

    1.1" 研究區(qū)概況

    研究區(qū)為位于廣州市花都區(qū)赤坭鎮(zhèn)境內(nèi)的水產(chǎn)養(yǎng)殖區(qū)(圖1),屬南亞熱帶季風(fēng)氣候,降水充沛,多年平均降水量為1 818.7 mm,平均蒸發(fā)量為1 611.7 mm,水文地質(zhì)單元面積為4.82 km2。研究區(qū)的地勢(shì)大致為北部高、南部低,地貌類型以平原為主,丘陵和臺(tái)地次之。研究區(qū)平原地下水的主要來(lái)源為灌溉水、大氣降雨、地表水體的入滲和丘陵臺(tái)地等地下水的側(cè)向補(bǔ)給[20]。

    現(xiàn)場(chǎng)調(diào)查可知,研究區(qū)地下水主要有松散巖類孔隙水、層狀巖類裂隙水和覆蓋性碳酸鹽巖類巖溶水三大類。其中,松散巖類孔隙水呈條帶狀分布于養(yǎng)殖區(qū)西部及西南部,層狀巖類裂隙水零星分布于

    養(yǎng)殖區(qū)東北部及西北部,覆蓋性碳酸鹽巖類巖溶水則廣泛分布于養(yǎng)殖區(qū)。受地形地貌及地表水體影響,區(qū)域內(nèi)地下水的總體流向與地表水流動(dòng)方向基本一致,初步判斷為由北向南徑流、排泄。由圖2可知,區(qū)域內(nèi)土地利用類型主要為城鎮(zhèn)用地、水域(魚塘)和林地3類,其中城鎮(zhèn)用地面積占37.6%,水域面積占40.7%,林地面積占21.7%。

    1.2" 樣品采集與測(cè)試

    2022年7月,參考《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ 494—2009)[32]和《水質(zhì)樣品的保存和管理技術(shù)規(guī)定》(HJ 493—2009)[33],根據(jù)水產(chǎn)養(yǎng)殖區(qū)地表水與地下水分布情況,布設(shè)采樣點(diǎn)共21個(gè),包括12個(gè)地表水和9個(gè)淺層地下水。其中,地表水來(lái)源包括4個(gè)污水處理池和8個(gè)魚塘,地下水來(lái)源選取9個(gè)居民生活用井。具體點(diǎn)位分別命名為:W1、W2、W3、W4(4個(gè)污水處理池);Y1、Y2、Y3、Y4、Y5、Y6、Y7、Y8(8個(gè)魚塘);J1、J2、J3、J4、J5、J6、J7、J8、J9(9個(gè)地下水民井)。

    采集生活污水樣品時(shí),在污水處理池排水口處采集尾水;采集魚塘水樣品時(shí),在魚塘下游出水口附近取樣,將采樣所需貝勒管沒入水面下20 cm,盡量不要擾動(dòng)河底沉積物,以免污染水樣;在采集地下水樣品前,需開泵10~15 min去除水管與井中存留的積水,抽出水量達(dá)到井內(nèi)水體積的3~5倍后,每隔5 min再測(cè)一次現(xiàn)場(chǎng)理化指標(biāo),直至連續(xù)3次的相對(duì)誤差

    在±10%以內(nèi)再進(jìn)行取樣。所有樣品瓶在采集前先用超純水潤(rùn)洗,再用待采水樣潤(rùn)洗2遍后采集1 L水樣,所用采樣瓶為1 L高密度聚乙烯瓶。

    現(xiàn)場(chǎng)理化指標(biāo),包括溫度(T)、pH值、溶解氧(DO)、電導(dǎo)率(EC)和氧化還原電位(ORP),采用美國(guó)產(chǎn)便攜式哈希(portable HACH)多功能水質(zhì)分析儀測(cè)試。水樣用0.45 μm濾膜過濾后,添加相應(yīng)的保護(hù)劑(濃HNO3、H2SO4)?,F(xiàn)場(chǎng)記錄好取樣編號(hào)、取樣地點(diǎn)、日期和測(cè)定參數(shù),樣品采集完成后放入冷藏箱,盡快由專人護(hù)送交實(shí)驗(yàn)室完成測(cè)試。實(shí)驗(yàn)室監(jiān)測(cè)指標(biāo)包括總氮(TN)、硝態(tài)氮(NO3--N)、氨氮(NH4+-N)、碘化物(I-)和化學(xué)需氧量(CODMn)。采用DX120型離子色譜儀和ICPAES分別對(duì)陰離子(Cl-、SO42-、NO3-N、I-)和陽(yáng)離子(K+、Ca2+、Na+、Mg2+、NH4+N)質(zhì)量濃度進(jìn)行測(cè)定;采用酸堿指示劑滴定法測(cè)定CO32- 和HCO3-[34]質(zhì)量濃度;采用堿性過硫酸鉀消解紫外分光光度法(HJ 636—2012)對(duì)TN質(zhì)量濃度進(jìn)行分析;通過烘箱干燥分析溶解性總固體(TDS)質(zhì)量濃度;采用高錳酸鉀法測(cè)定CODMn質(zhì)量濃度。

    添加10%的平行樣品,且平行樣品之間的相對(duì)偏差需在±5%以內(nèi),通過電荷平衡誤差分析評(píng)估樣品整體測(cè)試的準(zhǔn)確性。

    ec=ca-cbca+cb×100%。(1)

    式中:ec為地下水的點(diǎn)位電荷平衡誤差;ca為陽(yáng)離子(K+、Ca2+、Na+、Mg2+)總濃度(mmol/L);cb為陰離子(HCO3-、Cl-、SO42-、CO32-)總濃度(mmol/L)。

    測(cè)試結(jié)果顯示,地下水的ec均在±5%之間,| ec |為0.006~0.046,平均值為0.025,說(shuō)明測(cè)試結(jié)果準(zhǔn)確。

    1.3" 評(píng)價(jià)方法

    1.3.1" 水質(zhì)分析與評(píng)價(jià)

    地表水評(píng)價(jià)標(biāo)準(zhǔn)以《淡水池塘養(yǎng)殖水排放要求》(SC/T 9101—2007)[35]中的Ⅱ類標(biāo)準(zhǔn)為主,未規(guī)定的指標(biāo)按《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002) [36]Ⅲ類標(biāo)準(zhǔn)(適用于水產(chǎn)養(yǎng)殖區(qū))指標(biāo)評(píng)價(jià),超過以上標(biāo)準(zhǔn)則為超標(biāo);地下水評(píng)價(jià)標(biāo)準(zhǔn)以《生活飲用水水源水質(zhì)標(biāo)準(zhǔn)》(CJ 3020—93)[37]中的Ⅱ類標(biāo)準(zhǔn)為主,未規(guī)定的指標(biāo)按《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T 14848—2017)[38] Ⅲ類標(biāo)準(zhǔn)作為超標(biāo)依據(jù);共同采用單因子污染指數(shù)法對(duì)pH值和DO、CODMn、NH4+N、NO3-N、I-質(zhì)量濃度進(jìn)行單項(xiàng)評(píng)價(jià),并用內(nèi)梅羅污染指數(shù)法對(duì)采樣點(diǎn)位水質(zhì)進(jìn)行整體評(píng)價(jià)[39]。該方法的運(yùn)用主要分為以下幾個(gè)步驟。

    1)選擇評(píng)價(jià)因子。根據(jù)研究區(qū)的監(jiān)測(cè)數(shù)據(jù),選取pH和DO、CODMn、NH4+N、NO3-N、I-質(zhì)量濃度指標(biāo)作為評(píng)價(jià)因子。

    2) 明確單項(xiàng)評(píng)價(jià)因子評(píng)分值[40]。各單項(xiàng)評(píng)價(jià)因子按地表水、地下水對(duì)應(yīng)的5類標(biāo)準(zhǔn)進(jìn)行評(píng)分。各單項(xiàng)評(píng)價(jià)因子的評(píng)分值Fi為0、1、3、6和10時(shí)分別對(duì)應(yīng)Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ類質(zhì)量類別[41]。

    3) 確定水質(zhì)綜合指數(shù)FN,計(jì)算公式為

    FN=" F2i,avg+F2i,max2 。(2)

    式中:FN為內(nèi)梅羅污染指數(shù);Fi,avg為Fi的平均值;Fi,max為Fi中最大值。

    4)進(jìn)行地下水質(zhì)量分級(jí)。地下水水質(zhì)按FN值分類,F(xiàn)N值的取值范圍為lt;0.80、0.80~2.50、2.50~4.25、4.25~7.20和gt;7.20時(shí),分別對(duì)應(yīng)地下水Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ類質(zhì)量類別。

    指標(biāo)超標(biāo)率按式(3)計(jì)算:

    S=ninj×100%。(3)

    式中:S為超標(biāo)率;ni為超標(biāo)樣本數(shù);nj為監(jiān)測(cè)樣本數(shù)。

    1.3.2" 健康風(fēng)險(xiǎn)評(píng)估

    根據(jù)世界衛(wèi)生組織(WHO)公布的致癌物清單,碘是一種非致癌物,因此本研究主要關(guān)注碘的非致癌風(fēng)險(xiǎn)潛力。飲用水源中的碘主要通過飲用過量攝入對(duì)人類健康構(gòu)成風(fēng)險(xiǎn)[42]。飲用攝入的非致癌風(fēng)險(xiǎn)評(píng)估采用美國(guó)環(huán)境保護(hù)署推薦的人體健康風(fēng)險(xiǎn)評(píng)估模型[43],人群根據(jù)生理和行為差異分為兒童、成人女性和成人男性,使用公式(4) [44]和(5) [45]計(jì)算各年齡組計(jì)算每日攝入量(estimated daily intake)和碘的非致癌風(fēng)險(xiǎn)(hazard quotient)。

    Di=ρiCdBw; (4)

    Qi=DiKi。 (5)

    式中:Di為各年齡組目標(biāo)污染物預(yù)計(jì)攝入量(mg·kg-1·d-1);ρi為飲用水中目標(biāo)污染物的質(zhì)量濃度(mg/L);Cd為飲用水的平均攝入量,兒童、成年女性、成年男性分別取0.7、1.5、1.5 L·d-1;Bw為平均體質(zhì)量,兒童、成年女性、成年男性分別取15、55、75 kg[45];Qi為非致癌風(fēng)險(xiǎn)值;Ki為通過飲用水?dāng)z入的污染物的參考劑量(mg·kg-1·d-1),碘的攝入?yún)⒖紕┝繛?.01 mg·kg-1·d-1[46]。當(dāng)Qi>1時(shí),非致癌風(fēng)險(xiǎn)超過可接受水平,說(shuō)明其可能對(duì)健康產(chǎn)生不利影響。

    2" 結(jié)果與討論

    2.1" 淺層地下水水化學(xué)特征

    圖3為研究區(qū)地下水水化學(xué)指標(biāo)分析結(jié)果。由圖3a—c可知:研究區(qū)淺層地下水pH值為4.85~7.67,平均值為6.47,其中J1、J2和J6水樣pH分別為6.09、4.85和4.85,均超出《生活飲用水水源水質(zhì)標(biāo)準(zhǔn)》(CJ 3020—93)[37]中的Ⅱ類標(biāo)準(zhǔn)的閾值(6.5~8.5);ρ(TDS)值為146.09~385.11 mg/L,平均值為280.56 mg/L(閾值1 000 mg/L);ρ(CODMn)為0.12~1.17 mg/L,平均值為0.55 mg/L(閾值6 mg/L)。ρ(TDS)和ρ(CODMn)均達(dá)標(biāo),研究區(qū)淺層地下水總體上呈弱酸性低溶解性總固體特征。

    由圖3d—k可知:主要的陽(yáng)離子平均質(zhì)量濃度由高到低為Ca2+(45.08 mg/L)、Na+(11.47 mg/L)、K+(6.48 mg/L)、Mg2+(2.39 mg/L),Ca2+是水中TDS的主要成分;主要的陰離子平均質(zhì)量濃度由高到低為HCO3-(121.06 mg/L)、SO42-(27.81 mg/L)、Cl-(20.50 mg/L)、NO3-N(7.20 mg/L)。其中:ρ(SO42-)為 8.12~59.39 mg/L,ρ(Cl-)為11.38~38.96 mg/L,ρ(SO42-)和ρ(Cl-)均達(dá)標(biāo)(閾值250 mg/L);ρ(NO3-N)為0.08~24.86 mg/L,超標(biāo)率為11%( 閾值20 mg/L)。

    2.2" 碘化物分布特征

    采用SPSS 27.0對(duì)研究區(qū)ρ(I-)進(jìn)行差異性分析,空間變異性大小取決于變異系數(shù)(CV)值,按照變異系數(shù)的劃分等級(jí),CVlt;10%為弱變異性,10%≤CV<100%為中等變異,CV≥100%則為強(qiáng)變異性[47]。研究區(qū)地表水ρ(I-)變異系數(shù)為84.4%,屬于中等變異;地下水ρ(I-)變異系數(shù)為158.2%,屬于強(qiáng)變異。采用地統(tǒng)計(jì)學(xué)方法中的反距離插值法分別對(duì)研究區(qū)地表水及地下水中ρ(I-)進(jìn)行空間插值分析(圖4),研究區(qū)地表水和地下水ρ(I-)在空間分布上具有一致性。地表水ρ(I-)為5~79 μg/L,未超過高碘水源上限閾值100 μg/L,地表水ρ(I-)分布呈現(xiàn)出整體分散、西部較高的空間分布特征,其中ρ(I-)最高點(diǎn)位為Y3(79 μg/L)和W2(74 μg/L)。淺層地下水源中碘化物的分布也呈現(xiàn)出整體分散、西部較高的空間分布特征,淺層

    地下水中ρ(I-)為2~343 μg/L,其中:J1、J2、J3、J6、J7為缺碘型水源(ρ(I-)<10 μg/L),占比55.6%,在研究區(qū)均勻分布;J8為適碘型水源,占比11.1%;J4、J5、J9為高碘型水源(ρ(I-)>100 μg/L),占比33.3%,主要集中在研究區(qū)西部,呈現(xiàn)出明顯的聚集性。雷萬(wàn)杉等[48]研究表明,區(qū)域尺度(如大洲、國(guó)家和?。┥?,地質(zhì)背景、氣候和地形因素對(duì)物質(zhì)的分布至關(guān)重要,在局部尺度(如縣、鄉(xiāng)和農(nóng)田)上,人類活動(dòng)特別是農(nóng)業(yè)活動(dòng)的影響更為顯著?,F(xiàn)場(chǎng)調(diào)查得知,J4、J5、J9民井周邊魚塘均養(yǎng)殖桂花魚,而聚維酮碘(PVPI)作為消毒劑被廣泛用于水產(chǎn)養(yǎng)殖中,推測(cè)淺層高碘地下水成因與聚維酮碘

    消毒劑的使用有關(guān)。已知聚維酮碘是單質(zhì)碘與聚乙烯基吡咯烷酮的無(wú)定型絡(luò)合物,碘在表面活性劑形成膠束中心而被運(yùn)輸,其中聚乙烯比咯烷酮作為有機(jī)載體,不僅能持續(xù)釋放游離碘,還可充當(dāng)?shù)獾膬?chǔ)存庫(kù)[49];此外,魚塘淤泥通常富含大量有機(jī)質(zhì),淤泥中碘的吸收與有機(jī)質(zhì)含量直接相關(guān)[50]。有研究表明,添加到富含有機(jī)物土壤中的碘離子和碘酸鹽會(huì)轉(zhuǎn)化為碘的有機(jī)形式,并被強(qiáng)烈保留在土壤中[51],而土壤環(huán)境中的pH、氧化還原電位(Eh)、鐵鋁(氫)氧化物、微生物等因素都對(duì)碘在土壤地下水系統(tǒng)中的遷移吸附產(chǎn)生重要影響[5255]。因此,推測(cè)研究區(qū)地表水Y3和淺層地下水J4、J5、J9中ρ(I-)較高的原因與頻繁使用聚維酮碘消毒劑的漁業(yè)活動(dòng)的有關(guān)。而魚塘未檢出超高質(zhì)量濃度碘化物的原因則可能是聚維酮碘穩(wěn)定性受溫度、光照、pH、有機(jī)質(zhì)等因素影

    響[5658],高溫強(qiáng)光照條件以及魚塘堿性環(huán)境(pHmean=8.01)使得聚維酮碘穩(wěn)定性變差,從而在短時(shí)間內(nèi)釋放出大量的游離碘,釋放后游離碘在水中發(fā)生歧化

    反應(yīng)最終產(chǎn)生大量的I-和IO3-。Fuge等[59]研究表明,在酸性還原條件下,游離碘以穩(wěn)定態(tài)的碘化物(I-)的形式存在,在堿性氧化條件下,游離碘轉(zhuǎn)化為更穩(wěn)定的碘酸鹽(IO3-)形式,而魚塘中長(zhǎng)期運(yùn)行的增氧設(shè)備使得魚塘氧化還原電位(Eh)介于91~197 mV之間;因此,魚塘在堿性氧化條件下均未檢出超高質(zhì)量濃度碘化物。

    由圖2可知,地下水點(diǎn)位J5、J9周邊土地利用類型主要為城鎮(zhèn)用地,J4中 ρ(I-)最高值(343 μg/L)出現(xiàn)在W2污水處理池附近;結(jié)合現(xiàn)場(chǎng)調(diào)查得知,W2污水處理設(shè)施已停運(yùn)但仍有污水排出,且周邊建有垃圾回收站,因此,推測(cè)J4超高ρ(I-)可能與周邊含碘污水排放和垃圾滲濾液的滲入有關(guān)。

    2.3" 地下水碘化物來(lái)源

    2.3.1" 碘化物與水化學(xué)指標(biāo)的相關(guān)性分析與主成分分析

    斯皮爾曼相關(guān)性分析(圖5)結(jié)果顯示,SO42-、Ca2+、Mg2+、HCO3-、Cl-、CODMn、Na+、I-為水中離子的主要成分,其質(zhì)量濃度都與ρ(TDS)呈正相關(guān)關(guān)系,其中:ρ(I-)與ρ(TDS)呈強(qiáng)相關(guān)關(guān)系(r2≥0.8, r為相關(guān)系數(shù)),表明環(huán)境中的蒸發(fā)、濃縮作用促進(jìn)了I-的富集;ρ(I-)與ρ(Ca2+)、ρ(Mg2+)、ρ(HCO3-)、ρ(CODMn)呈較強(qiáng)正相關(guān)關(guān)系(0.6≤r2<0.8),說(shuō)明I-可能與這些離子具有共同的來(lái)源。已有研究[60]表明,蒸發(fā)巖鹽溶解、海水入侵和城鎮(zhèn)化伴隨的污水泄漏可為地下水中的Ca2+、Mg2+ 和I-提供來(lái)源;此外,有報(bào)道[24, 61]表明礦物風(fēng)化和溶解有助于地下水中碘的富集,碳酸鹽巖性地區(qū)的碘含量相比非碳酸鹽巖性地區(qū)更豐富。由相關(guān)性分析(圖5)可知,ρ(HCO3-) 和pH值呈較強(qiáng)正相關(guān)關(guān)系,也與地下水還原因子ρ(NH4+N)呈弱正相關(guān)關(guān)系(0<r2≤0.3);說(shuō)明地下水中的HCO3-控制著pH值變化,二者共同反映了淺層地下水在弱堿性還原環(huán)境下有利于高碘地下水的賦存。此外,ρ(I-)與ρ(CODMn)呈較強(qiáng)正相關(guān)關(guān)系,說(shuō)明二者可能具有共同的來(lái)源;這與大同盆地地下水中碘的來(lái)源解析結(jié)果相似,即地下水環(huán)境沉積的有機(jī)物匯集了大量的碘,這些有機(jī)物中的碘又成為地下水中的來(lái)源[62]。

    根據(jù)水源類型將樣品分為地表水(魚塘、污水池)和地下水類型。通過運(yùn)用主成分分析法(PCA)進(jìn)一步對(duì)碘的來(lái)源進(jìn)行分析。通常條件下,碘的來(lái)源包含巖石的浸出、土壤雨水的淋溶以及受人類生產(chǎn)生活的影響[63]。如圖6a所示,地表水中ρ(NO3-N)、ρ(TN)與ρ(I-)均投影于主成分1(PC1)的正半軸,且距離I-較近,表明ρ(NO3-N)、ρ(TN)與ρ(I-)呈正相關(guān)關(guān)系,而ρ(NO3-N)、ρ(TN)又與頻繁的農(nóng)業(yè)活動(dòng)有關(guān)[6466];說(shuō)明農(nóng)業(yè)活動(dòng)是地表水I-的主要來(lái)源之一。

    地下水源的主成分分析(圖6b)中,PC1的方差貢獻(xiàn)率為47.327%(表1),pH值、ρ(TDS)與ρ(I-)在PC1上均為正載荷,pH與ρ(I-)呈正相關(guān)關(guān)系;說(shuō)明pH值影響著ρ(I-)的大小,在酸性較強(qiáng)的環(huán)境中,I-的氧化速度

    加快,導(dǎo)致I-損失加大,這與黃適爾等[42]研究結(jié)果一致。此外,也有研究[25]表明,地下水沉積物在弱堿性環(huán)境下會(huì)降低對(duì)I-的吸附能力,從而促進(jìn)I-釋放遷移至水相。主成分分析(表1)顯示,在PC1中,ρ(TDS)、pH、ρ(Ca2+)、ρ(HCO3-)、ρ(DO)、ρ(CODMn)、ρ(NH4+N)、EC、ρ(I-)荷載較大。有研究[67]表明,珠江三角洲沿海地區(qū)富含碳酸鹽的巖石通常富含碘有機(jī)質(zhì),研究區(qū)淺層地下水主要為重碳酸型水;表明弱堿性還原條件下,研究區(qū)孔隙含水層沉積地層中富碘有機(jī)質(zhì)的降解導(dǎo)致地下水中I-、NH4+N和HCO3-富集[61]。此外,也有研究[60, 68]指出,村鎮(zhèn)還原性污

    水通常含有較高質(zhì)量濃度的I-和NH4+N,而還原性污水中化糞池廢水以及生活污水的泄漏是地下水中的NH4+N和I-的來(lái)源之一。研究區(qū)地下水J4(ρ(I-)=343 μg/L)、J5(ρ(I-)=103 μg/L)和

    J9(ρ(I-)=107 μg/L)均為高碘地下水點(diǎn)位,污水處理池W2檢出ρ(I-)較高,且W2點(diǎn)距離點(diǎn)位J4較近,此外現(xiàn)場(chǎng)調(diào)查發(fā)現(xiàn),高碘地下水點(diǎn)位J4周邊還建有垃圾填埋場(chǎng);因此,J4點(diǎn)位檢出超高質(zhì)量濃度的碘化物則有可能是污水處理池尾水和垃圾滲濾液共同作用的結(jié)果,而J5、J9位于J4周邊,也可能受到J4高碘地下水遷移的影響。調(diào)查顯示,近些年來(lái)沿海地區(qū)水產(chǎn)養(yǎng)殖業(yè)日漸發(fā)達(dá),且河道人工挖沙現(xiàn)象越發(fā)普遍[69],挖沙過程破壞了孔隙含水層,使得歷史沉積層中富碘土壤淋溶解析,進(jìn)而促使孔隙含水層中碘富集[7071],同時(shí)聚維酮碘作為有機(jī)高分子消毒劑也在J5、J9周邊魚塘被廣泛使用;故J5、J9點(diǎn)位檢出超高質(zhì)量濃度的碘化物則有可能是J4高碘地下水的遷移與魚塘富碘底泥的淋濾共同作用的結(jié)果。因此,PC1代表弱堿性還原環(huán)境下,垃圾滲濾液和還原性富碘污水滲漏入滲以及魚塘富碘土壤的淋濾對(duì)研究區(qū)淺層高碘地下水的綜合影響。

    由研究區(qū)淺層地下水主要離子主成分分析(表1)知,主成分2(PC2)的方差貢獻(xiàn)率為16.934%,ρ(TN)、ρ(NO3 -N)、ρ(Na+)、ρ(Cl-)、ρ(SO42-)在PC2荷載較大。其中,ρ(TN)、ρ(NO3-N)與頻繁的農(nóng)業(yè)活動(dòng)有關(guān),ρ(SO42-)主要與珠江三角洲工業(yè)酸雨對(duì)淺層地下水的入滲有關(guān)[72]。主成分3(PC3)的方差貢獻(xiàn)率為15.412%,ρ(Mg2+)和ρ(K+)在PC3荷載較大。有調(diào)查[60]發(fā)現(xiàn),海水入侵可為地下水中的Na+、Mg2+、Cl-和I-提供來(lái)源證明,珠江三角洲孔隙含水層曾發(fā)生過三次海水入侵;因此,PC2和PC3代表農(nóng)業(yè)活動(dòng)、工業(yè)污染和海侵因素等綜合因素對(duì)該區(qū)地下水化學(xué)組分的影響。I-在PC2和PC3中荷載低,說(shuō)明農(nóng)業(yè)活動(dòng)、工業(yè)污染和海水入侵對(duì)淺層

    地下水中碘的富集影響較小。

    2.3.2" 水化學(xué)控制因素

    利用反映水化學(xué)組分特征的Gibbs圖分析ρ(Na+)和ρ(HCO3-)等水化學(xué)指標(biāo)對(duì)碘來(lái)源的影響,可進(jìn)一步研究碘來(lái)源與水化學(xué)控制因素的關(guān)系[7374]。水化學(xué)中陽(yáng)離子的質(zhì)量濃度比(圖7a)、陰離子的質(zhì)量濃度比(圖7b)與地下水中總?cè)芙饬康年P(guān)系能夠綜合反映地下水的水化學(xué)特征和組成,揭示地下水的主導(dǎo)因素。研究區(qū)高碘地下水ρ(I-)主要介于103~343" μg/L,陽(yáng)離子ρ(Na+)/(ρ(Na+)+ρ(Ca2+))為0.07~0.36,陰離子ρ(Cl-)/(ρ(Cl-)+

    ρ(HCO3-))為0.07~0.21,J4、J5、J9分布在巖石風(fēng)化區(qū)域,水化學(xué)驅(qū)動(dòng)因素受地質(zhì)因素影響較大,其主

    要受巖石風(fēng)化作用控制。從圖7a 可知,相較于其他點(diǎn)位(J1、J2、J3、J6、J7、J8),ρ(I-)質(zhì)量濃度較高的J4、J5和J9受到的蒸發(fā)、濃縮作用更強(qiáng)(ρ(TDS)值更大);說(shuō)明蒸發(fā)濃縮作用可能對(duì)水中I-的富集也起到了一定作用。

    2.3.3" 淺層高碘地下水化學(xué)特征

    Piper三線圖可以反映水體中主要離子組成的變化,從而反映水化學(xué)的組成特征[75]。如圖8所示,研究區(qū)淺層地下水水化學(xué)類型主要可分為Ca2+·HCO3-,Ca2+·Na+·HCO3-2種類型。研究區(qū)淺層地下水處于補(bǔ)給徑流區(qū),水化學(xué)類型以HCO3-占絕對(duì)優(yōu)勢(shì),說(shuō)明從補(bǔ)給區(qū)到排泄區(qū)地下水

    化學(xué)類型通常呈現(xiàn)從HCO3-→SO42-→Cl-變化的規(guī)律,這與朱海勇等[76]研究結(jié)果相符。研究區(qū)淺層地下水總體上呈弱酸性低溶解性總固體特征(圖8),陽(yáng)離子以Ca2+為主,陰離子以HCO3-為主,除J6為Ca2+·Na+·HCO3-水化學(xué)類型外,其余的點(diǎn)位都集中在第2區(qū)域,都屬于低溶解性總固體Ca2+·HCO3-為主導(dǎo)的水化學(xué)類型。民井J2、J6常年敞口開放,其淺層地下水酸化明顯可能受酸雨入滲補(bǔ)給影響。

    如圖9a所示,研究區(qū)地下水pH值范圍為4.85~7.67,J4、J5、J9高碘地下水pH范圍在7.19~7.67之間,為弱堿性。I-易被含水介質(zhì)中帶正電的物質(zhì)吸附[7778],pH值的增大會(huì)使膠體和黏土礦物帶更多負(fù)電荷,降低對(duì)I-的吸附[7980],所以當(dāng)pH>7.0時(shí)(圖9a),地下水中ρ(I-)明顯升高。相較于其他點(diǎn)位,J4、J5、J9高碘地下水中ρ(CODMn)和ρ(NH4+N)均較高(圖9b、c),而ρ(DO)、ρ(NO3-N)和Eh均較低(圖9d—f)。研究區(qū)地下水Eh集中于-43~286 mV,變化幅度較大;J4、J5、J9高碘地下水Eh均小于0 (圖9e),范圍在-43~-11 mV之間,研究區(qū)地下水J1、J2、J3、J6、J7和J8屬于偏氧化水(圖9e),J4、J5和J9屬于偏還原水,厭氧還原環(huán)境不利于氧化態(tài)的IO3- 形成[81],I-是研究區(qū)高碘地下水中穩(wěn)定存在的碘形態(tài)。綜上所述,氧化還原環(huán)境和pH值是影響研究區(qū)地下水中ρ(I-)高低的重要因素,弱堿性還原環(huán)境有利于地下水中碘的富集。

    2.4" 水質(zhì)污染評(píng)價(jià)

    目前,相關(guān)研究多采用單因子指數(shù)法、模糊綜合評(píng)價(jià)法和內(nèi)梅羅綜合指數(shù)法等進(jìn)行水質(zhì)評(píng)價(jià),或采用美國(guó)環(huán)保署(USEPA)健康風(fēng)險(xiǎn)評(píng)估(HHRA)模型進(jìn)行人類健康風(fēng)險(xiǎn)評(píng)價(jià),而將水質(zhì)評(píng)價(jià)及健康風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合來(lái)對(duì)水質(zhì)進(jìn)行綜合評(píng)價(jià)的研究較少。將水質(zhì)評(píng)價(jià)與健康風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合,可以更全面地了解水環(huán)境質(zhì)量以及地下水作為飲用水的安全性,也可為地下水污染防治和飲用水風(fēng)險(xiǎn)管理提供科學(xué)依據(jù)[82]。因此,本研究先采用單因子污染指數(shù)法和內(nèi)梅羅綜合指數(shù)法分別對(duì)研究區(qū)地表水和地下水的pH值、ρ(NH4+N)、ρ(NO3-N)、ρ(TN)、ρ(DO)、ρ(CODMn)、ρ(I-)和ρ(TDS)進(jìn)行單項(xiàng)評(píng)價(jià)和水質(zhì)整體評(píng)價(jià)(表2、3);然后將水質(zhì)評(píng)價(jià)與健康風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合,通過人類健康風(fēng)險(xiǎn)評(píng)估模型評(píng)估研究區(qū)居民飲用攝入碘的非致癌風(fēng)險(xiǎn)。

    2.4.1" 地表水水質(zhì)污染特征

    地表水水質(zhì)污染指數(shù)評(píng)價(jià)結(jié)果(表2)表明,碘化物因子均符合《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002) [36]Ⅲ類標(biāo)準(zhǔn),pH、ρ(DO)、ρ(TN)、ρ(NH4+N)、ρ(CODMn)和ρ(NO3-N)因子存在超標(biāo)。研究區(qū)超標(biāo)指標(biāo)的超標(biāo)率(表3)顯示,地表水中ρ(CODMn)、ρ(TN)、ρ(DO)、ρ(NH4+N)、pH和ρ(NO3-N)超標(biāo)率由大到小分別為75.0%、58.3%、50.0%、33.3%、25.0%、8.3%。由表2可知:

    W1和W2地表水質(zhì)量類別為Ⅳ類水,Ⅳ類水占比為16.7%,水質(zhì)較差(僅討論最低水質(zhì)類別);魚塘單因子污染指數(shù)均值>1的有ρ(DO) 、ρ(CODMn)和ρ(TN)。其中:魚塘ρ(DO)普遍較高的

    原因可能是,氧氣泵或其他氧氣供應(yīng)設(shè)備過于強(qiáng)勁,導(dǎo)致水中氧氣濃度升高;并且現(xiàn)場(chǎng)在采樣時(shí)發(fā)現(xiàn)魚塘存在大量藻類植物,藻類在白天進(jìn)行光合作用時(shí)會(huì)產(chǎn)生大量

    氧氣。ρ(CODMn)、ρ(TN)超標(biāo)是因?yàn)轸~塘中的代謝產(chǎn)物和殘餌是傳統(tǒng)水產(chǎn)養(yǎng)殖模式環(huán)境中有機(jī)質(zhì)污染負(fù)荷和ρ(TN)超標(biāo)的主要存在形式[83],益生菌的施用次數(shù)也是ρ(CODMn)、ρ(TN)較高的重要影響因素之一[40];此外,陳舊的魚塘富含有機(jī)質(zhì)的淤泥過厚,內(nèi)源沉積物釋放也很可能是造成ρ(CODMn)居高不下的原因。

    2.4.2" 地下水水質(zhì)污染特征

    地下水水質(zhì)污染指數(shù)評(píng)價(jià)結(jié)果(表4)表明,ρ(NH4+N)、ρ(CODMn)均符合《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T 14848—2017)[38]Ⅲ類標(biāo)準(zhǔn)。表4內(nèi)梅羅污染指數(shù)計(jì)算結(jié)果表明:J7、J8地下水質(zhì)量類別為Ⅰ類,共占比22.2%;J1、J3、J5、J9地下水質(zhì)量類別為Ⅱ類,共占比44.5%;J2、J4、J6地下水質(zhì)量類別為Ⅲ類,共占比33.3%,整體水質(zhì)較好。由表3可知,pH、ρ(I-)和ρ(NO3-N)因子存在超標(biāo),pH、ρ(I-)和ρ(NO3-N)超標(biāo)率由大到小分別為33.3%、33.3%、11.1%。表4單因子污染指數(shù)pH、ρ(I-)均值分別為1.83和0.89,均值較高,超標(biāo)情況較嚴(yán)重,超標(biāo)率均為33.3%。此外J2和J6的pH現(xiàn)場(chǎng)測(cè)得均為4.85,統(tǒng)計(jì)結(jié)果[84]表明,1982—2007年珠江三角洲地區(qū)降水pH值平均介于4.48~5.92,且呈逐年下降的態(tài)勢(shì);而J2和J6位于研究區(qū)水文地質(zhì)單元補(bǔ)給區(qū),降雨的影響可能使得補(bǔ)給區(qū)的pH整體呈偏酸性,而其他位于徑流區(qū)的點(diǎn)位pH未超標(biāo)原因可能是地下水從補(bǔ)給區(qū)到排泄區(qū)的流動(dòng)過程中跟地層間發(fā)生反應(yīng),酸性被中和[85]。

    2.4.3" 健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

    人類健康風(fēng)險(xiǎn)評(píng)估按照行為和生理差異將人群分為兒童、成年女性和成年男性。通過模型評(píng)估居民攝入碘的非致癌風(fēng)險(xiǎn),結(jié)果如圖10所示。由圖10可知:研究區(qū)非致癌風(fēng)險(xiǎn)值(Qi)最高人群出現(xiàn)在兒童中,兒童的Qi介于0.010~1.601之間,其中J4點(diǎn)

    位的Qi超過允許的限值(Qi=1.0),說(shuō)明該點(diǎn)位飲用水源對(duì)兒童存在非致癌風(fēng)險(xiǎn);成年女性的Qi值區(qū)間為0.050~0.935,成年男性的Qi值區(qū)間為0.010~0.686,均低于允許的限值;成年男性

    的Qi平均值

    (0.134)與成年女性的Qi平均值(0.182)均低于兒童的Qi平均值(0.312),說(shuō)明飲用水中的ρ(I-)不會(huì)通過飲用的途徑對(duì)成年女性和成年男性構(gòu)成非致癌風(fēng)險(xiǎn)。因此,在高碘地下水的治理修復(fù)方面,建議加強(qiáng)對(duì)民井J4 ρ(I-)的常規(guī)監(jiān)控,增加除碘設(shè)施,還可構(gòu)建新型材料以增加碘的吸附性能。例如,米生雍等[86]研究表明,異喹啉修飾金屬有機(jī)骨架配合物通過提高對(duì)CO2的吸附量,選擇性地分離甲烷中的CO2,進(jìn)而可對(duì)溶液中的碘產(chǎn)生良好的吸附性能。

    3" 結(jié)論

    1)研究區(qū)淺層地下水水化學(xué)類型主要為Ca2+·HCO3-和Ca2+·Na+·HCO3-。研究區(qū)ρ(I-)為2~343 μg/L,地表水和淺層地下水ρ(I-)在空間分布上存在一致性,呈現(xiàn)出整體分散、西部較高的空間分布特征。其中:地表水源(魚塘、污水處理池)ρ(I-)均未超標(biāo)(閾值100 μg/L),ρ(I-)較高點(diǎn)位Y2(79 μg/L)和W2(74 μg/L)均分布在研究區(qū)西部。9個(gè)淺層地下水點(diǎn)位中,有3個(gè)點(diǎn)位(J4、J5、J9)為高碘型水源(ρ(I-)>100 μg/L),占比為33.3%,主要集中在研究區(qū)西部,呈現(xiàn)出明顯的聚集性;有1個(gè)點(diǎn)位(J8)為適碘型水源,占比11.1%,分布于研究區(qū)東部;有5個(gè)點(diǎn)位(J1、J2、J3、J6、J7)為缺碘型水源(ρ(I-)<10 μg/L),占比為55.6%,在研究區(qū)均勻分布。

    2)相關(guān)性分析和主成分分析表明,弱堿性還原環(huán)境有利于碘的賦存,研究區(qū)西部地下水中I-富集可能源于孔隙含水層沉積地層中富碘有機(jī)質(zhì)的降解,村鎮(zhèn)居民還原性污水泄漏是地下水中的I-的重要來(lái)源之一;此外,魚塘富碘底泥土壤的淋濾可能是研究區(qū)淺層高碘地下水的另一重要來(lái)源。

    3)內(nèi)梅羅污染指數(shù)計(jì)算結(jié)果表明:研究區(qū)地表水中Ⅳ類水為較差類別,占比16.7%,其中,CODMn和TN污染最為嚴(yán)重;地下水整體水質(zhì)為Ⅲ類及以上類別,整體水質(zhì)較好,部分點(diǎn)位pH、碘化物污染較為嚴(yán)重。

    4)研究區(qū)J4點(diǎn)位的Qi超過允許的限值(Qi=1.0),該點(diǎn)位飲用水源對(duì)兒童存在非致癌風(fēng)險(xiǎn),若長(zhǎng)期飲用會(huì)對(duì)兒童健康帶來(lái)潛在威脅,應(yīng)引起重視。成年女性和成年男性不會(huì)通過飲用途徑造成攝入碘化物的非致癌風(fēng)險(xiǎn)。

    參考文獻(xiàn)(References):

    [1]" 王貴玲,陳德華,藺文靜,等.中國(guó)北方地區(qū)地下水資源的合理開發(fā)利用與保護(hù)[J].中國(guó)沙漠,2007,27(4):684689.

    Wang Guiling,Chen Dehua,Lin Wenjing,et al. Reasonable Exploitation and Utilization of Groundwater Resource in North China[J]. Journal of Desert Research,2007,27(4):684689.

    [2]" 賴健,張琳.淺談廣東省地下水監(jiān)測(cè)現(xiàn)狀[J].環(huán)境,2012,35(增刊1):111.

    Lai Jian,Zhang Lin. Discussion on the Status of Groundwater Monitoring in Guangdong Province[J]. Environmental,2012,35(Sup.1):111.

    [3]" Hang Z J,Luo G Z,Wang Z,et al. Study on Sustainable Utilization of Groundwater in North China Plain[J]. Resources Science,2009,31(3):335360.

    [4]" Andersen S,Guan H,Teng W,et al. Speciation of Iodine in High Iodine Groundwater in China Associated with Goitre and Hypothyroidism[J]. Biological Trace Element Research,2009,128(2):95103.

    [5]" Li J,Wang Y,Xie X,et al. Hydrogeochemistry of High Iodine Groundwater:A Case Study at the Datong Basin,Northern China[J]. Environmental Science:Rocesses amp; Impacts,2013,15(4):848859.

    [6]" Li J,Wang Y,Xue X,et al. Mechanistic Insights into Iodine Enrichment in Groundwater During the Transformation of Iron Minerals in Aquifer Sediments[J]. Science of the Total Environment, 2020,745:140922.

    [7]" Wei W,Nghiem A,Ma R,et al. Factors Controlling Iodine Enrichment in a Coastal Plain Aquifer in the North Jiangsu Yishusi Plain,China[J]. Journal of Contaminant Hydrology, 2021,243:103894.

    [8]" Voutchkova D D,Ernstsen V,Hansen B,et al. Assessment of Spatial Variation in Drinking Water Iodine and Its Implications for Dietary Intake:A New Conceptual Model for Denmark[J]. Science of the Total Environment,2014,493(15):432444.

    [9]" Azzakhnini I,Abdelouas A. Iodine Content in Groundwater of North Eastern Morocco and Its Relation with the Incidence of Goiter[J]. Materials Today: Proceedings,2019,13:11511160.

    [10]" lvarez" F,Reich M,Snyder G,et al. Iodine Budget in Surface Waters from Atacama:Natural and Anthropogenic Iodine Sources Revealed by Halogen Geochemistry and Iodine129 Isotopes[J]. Applied Geochemistry,2016,68:5363.

    [11]" Watts M J,Mitchell C J. A Pilot Study on Iodine in Soils of Greater Kabul and Nangarhar Provinces of Afghanistan[J]. Environmental Geochemistry and Health,2009,31(4):503509.

    [12]" Korobova E,Kolmykova L,Ryzhenko B,et al. Distribution and Speciation of Iodine in Drinking Waters from Geochemically Different Areas of Bryansk Region Contaminated After the Chernobyl Accident in Relation to Health and Remediation Aspects[J]. Journal of Geochemical Exploration,2018,184:311317.

    [13]" 賈清珍,張向東.對(duì)新修訂的《水源性高碘地區(qū)和高碘病區(qū)的劃定》標(biāo)準(zhǔn)的解讀[J].中華地方病學(xué)雜志,2017,36(3):226229.

    Jia Qingzhen, Zhang Xiangdong. Interpretation of the Newly Revised Standard of Demarcation of Water-Borne High Iodine Areas and High Iodine Disease Areas[J]. Chinese Journal of Endemiology, 2017,36(3):226229.

    [14]" Laurberg P,Pedersen I B,Knudsen N,et al. Environmental Iodine Intake Affects the Type of Nonmalignant Thyroid Disease[J]. Thyroid,2001,11(5):457469.

    [15]" Hetzel B S. Iodine Deficiency Disorders[J]. Lancet,1988, 331:13861387.

    [16]" Delange F. Risks and Benefits of Iodine Supplementation[J]. The Lancet,1998, 351:923924.

    [17]" Good World H O. Regional Office for the Eastern Mediterranean[J]. Eastern Mediterranean Health Journal,2013,27(2):111213.

    [18]" Zimmermann M B,Andersson M. Global Perspectives in Endocrinology:Coverage of Iodized Salt Programs and Iodine Status in 2020[J]. European Journal of Endocrinology,2021,185(1):R13R21.

    [19]" Wong G,Brewer P G. The Marine Chemistry of Iodine in Anoxic Basins[J]. Geochimica et Cosmochimica Acta,1977,41(1):151159.

    [20]" 龐園,張明珠,龐志研,等.廣花盆地地下水水質(zhì)評(píng)價(jià)及變化趨勢(shì)分析[J].水電能源科學(xué),2017,35(11):4043.

    Pang Yuan,Zhang Mingzhu,Pang Zhiyan,et al. Water Quality Assessment and Trend Analysis of Groundwater in Guanghua Basin[J]. Water Resources and Power,2017,35(11):4043.

    [21]" 劉君科,文俊,支兵發(fā).廣東中山地區(qū)地下水環(huán)境質(zhì)量現(xiàn)狀評(píng)價(jià)與分析[J].地下水,2009,31(4):7679.

    Liu Junke,Wen Jun,Zhi Bingfa. Evaluation & Analysis on Groundwater Environment Quality Status in Zhongshan Area in Guangdong Province[J]. Ground Water,2009,31(4):7679.

    [22]" Shimamoto Y S,Takahashi Y,Terada Y. Formation of Organic Iodine Supplied as Iodide in a Soil-Water System in Chiba,Japan[J]. Environmental Science amp; Technology,2011,45(6):20862092.

    [23]" 吳飛,王曾祺,童秀娟,等.我國(guó)典型地區(qū)淺層高碘地下水分布特征及其賦存環(huán)境[J].水資源與水工程學(xué)報(bào),2017,28(2):99104.

    Wu Fei, Wang Zengqi, Tong Xiujuan, et al. The Distribution Characteristics and Storage Environments of Rich Iodine in Shallow Groundwater of Typical Areas in China[J]. Journal of Water Resources and Water,2017,28(2):99104.

    [24]" Zhang Y,Chen L,Cao S,et al. Iodine Enrichment and the Underlying Mechanism in Deep Groundwater in the Cangzhou Region,North China[J]. Environmental Science and Pollution Research,2021,28(9):1055210563.

    [25]" 呂曉立,劉景濤,韓占濤,等.快速城鎮(zhèn)化三角洲地區(qū)高碘地下水賦存特征及驅(qū)動(dòng)因素:以珠江三角洲為例[J].環(huán)境科學(xué),2022,43(1):339348.

    Lü Xiaoli,Liu Jingtao,Han Zhantao,et al. Geochemical Characteristics and Driving Factors of High-Iodine Groundwater in Rapidly Urbanized Delta Areas:A Case Study of the Pearl River Delta[J]. Environmental Science,2022,43(1):339348.

    [26]" 詹志薇.珠江三角洲典型地區(qū)地下水有機(jī)污染物環(huán)境影響及其脆弱性評(píng)價(jià)[D].廣州:華南農(nóng)業(yè)大學(xué),2016.

    Zhan Zhiwei. Environmental Impacts and Vulnerability Assessment of Organic Pollutants in Ground Water in the Typical Area of Pearl River Delta[D]. Guangzhou:South China Agricultural University,2016.

    [27]" 丁紫榮.珠三角地區(qū)典型污染源及地下水中藥物與個(gè)人護(hù)理品的分布特征研究[D].長(zhǎng)沙:湖南師范大學(xué),2015.

    Ding Zirong. Occurrence and Distribution of Pharmaceutical and Personal Care Products in Typical Pollution Sources and Groundwater in the Pearl River Delta[D]. Changsha:Hunan Normal University,2015.

    [28]" 程新偉,孫繼朝.珠江三角洲地區(qū)酸性地下水分布特征及其影響因素研究[J].地下水,2017,39(5):2527.

    Cheng Xinwei, Sun Jichao. Study on Distribution Characteristics of Acid Groundwater and Its Influencing Factors in the Pearl River Delta[J]. Ground Water,2017,39(5):2527.

    [29]" Tang Q,Xu Q,Zhang F,et al. Geochemistry of Iiodine-Rich Groundwater in the Taiyuan Basin of Central Shanxi Province, North China[J]. Journal of Geochemical Exploration,2013,135(6):117123.

    [30]" Li J,Wang Y,Guo W,et al. Iodine Mobilization in Groundwater System at Datong Basin,China:Evidence from Hydrochemistry and Fluorescence Characteristics[J]. Science of the Total Environment,2014,468(15):738745.

    [31]" Huang G,Liu C,Li L,et al. Spatial Distribution and Origin of Shallow Groundwater Iodide in a Rapidly Urbanized Delta:A Case Study of the Pearl River Delta[J]. Journal of Hydrology, 2020,585: 124860.

    [32]" 水質(zhì)采樣技術(shù)指導(dǎo):HJ 494—2009[S].北京:中國(guó)環(huán)境科學(xué)出版社,2009.

    Water Quality Guidance on Sampling Techniques:HJ 494—2009[S]. Beijing:China Environmental Science Press,2009.

    [33]" 水質(zhì)樣品的保存和管理技術(shù)規(guī)定:HJ 493—2009[S]. 北京:中國(guó)環(huán)境科學(xué)出版社,2009.

    Water Quality Technical Regulation of the Preservation and Handling of Samples:HJ 493—2009[S]. Beijing:China Environmental Science Press,2009.

    [34]" 國(guó)家環(huán)保局本書編委會(huì).水和廢水監(jiān)測(cè)分析方法[M].北京:中國(guó)環(huán)境科學(xué)出版社,1989.

    China Environmental Protection Bureau Book Compilation Committee. Methods for Monitoring and Analysis of Water and Wastewater[M]. Beijing:China Environmental Science Press,1989.

    [35]" 淡水池塘養(yǎng)殖水排放要求:SC/T 9101—2007[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,2007.

    Requirement for Water Discharge from Freshwater Aquaculture Pond:SC/T 9101—2007[S]. Beijing:Standards Press of China,2007.

    [36]" 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn):GB 3838—2002[S].北京:中國(guó)環(huán)境科學(xué)出版社,2002.

    Environmental Quality Standards for Surface Water:GB 3838—2002[S]. Beijing:China Environmental Science Press,2002.

    [37]" 生活飲用水水源水質(zhì)標(biāo)準(zhǔn):CJ 3020—93[S].北京:華文出版社,1993.

    Water Quality Standard for Drinking Water Sources:CJ 3020—93[S]. Beijing:Sino-Culture Press,1993.

    [38]" 地下水質(zhì)量標(biāo)準(zhǔn):GB/T 14848—2017[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2017.

    Standard for Ground Water Quality:GB/T 14848—2017 [S]. Beijing:Standards Press of China,2017.

    [39]" 李曉航,張飛,王娟,等.秋季艾比湖流域水質(zhì)綜合分析與評(píng)價(jià)[J].環(huán)境污染與防治,2017,39(6):588593.

    Li Xiaohang,Zhang Fei,Wang Juan,et al. Analysis and Assessment of Water Quality of Ebinur Lake Basin in Autumn[J]. Environmental Pollution amp; Control,2017,39(6):588593.

    [40]" 楊坤,胡星明,盧文軒.巢湖流域水產(chǎn)養(yǎng)殖尾水水質(zhì)分析與評(píng)價(jià)[J].生物學(xué)雜志,2023,40(1):8590.

    Yang Kun,Hu Xingming,Lu Wenxuan. Analysis and Evaluation of Water Quality of Aquaculture Tail Water of Chaohu Lake Basin[J]. Journal of Biology,2023,40(1):8590.

    [41]" 陳綺,龐園.兩種水質(zhì)評(píng)價(jià)方法在廣州市從化區(qū)地下水水質(zhì)評(píng)價(jià)中的應(yīng)用[J].人民珠江,2016,37(7):7881.

    Chen Qi,Pang Yuan. Application of Two Kinds of Water Quality Evaluation Methods for Groundwater of Conghua District in Guangzhou City[J]. Pearl River,2016,37(7):7881.

    [42]" 黃適爾,邊睿,郭加汛,等.山東、安徽省村鎮(zhèn)飲用水源中碘化物的分布、來(lái)源及存在的健康風(fēng)險(xiǎn)[J].環(huán)境科學(xué)學(xué)報(bào),2022,42(12):104113.

    Huang Shier,Bian Rui,Guo Jiaxun,et al. Distribution,Sources and Health Risks of Iodide in Rural Drinking Water Sources in Shandong and Anhui Province[J]. Acta Scientiae Circumstantiae,2022,42(12):104113.

    [43]" Li P,He X,Li Y,et al. Occurrence and Health Implication of Fluoride in Groundwater of Loess Aquifer in the Chinese Loess Plateau:A Case Study of Tongchuan,Northwest China[J]. Exposure and Health,2019,11(2):95107.

    [44]" Li D,Gao X,Wang Y,et al. Diverse Mechanisms Drive Fluoride Enrichment in Groundwater in Two Neighboring Sites in Northern China[J]. Environmental Pollution,2018,237:430441.

    [45]" Xiao Y,Hao Q,Zhang Y,et al. Investigating Sources,Driving Forces and Potential Health Risks of Nitrate and Fluoride in Groundwater of a Typical Alluvial Fan Plain[J]. Science of the Total Environment,2022,802:149909.

    [46]" He X,Li P,Wu J,et al. Poor Groundwater Quality and High Potential Health Risks in the Datong Basin,Northern China:Research from Published Data[J]. Environmental Geochemistry and Health,2021,43(2):791812.

    [47]" 代述勇,雷加強(qiáng),趙景峰,等.塔里木盆地南緣策勒綠洲區(qū)地下水TDS空間變異及水化學(xué)特征分析[J].中國(guó)沙漠,2010,36 (1):818.

    Dai Shuyong,Lei Jiaqiang,Zhao Jingfeng,et al. Variation and Evolution of Groundwater Chemistry in the Middle and Lower Reaches of the Tarim River[J]. Journal of Desert Research,2010,36(1):818.

    [48]" 雷萬(wàn)杉,林鑫,段星星,等.土壤中硒元素來(lái)源和遷移作用研究現(xiàn)狀[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2024,54(1):264278.

    Lei Wanshan,Lin Xin,Duan Xingxing,et al. Source and Migration of Selenium in the Soil[J]. Journal of Jilin University (Earth Science Edition),2024,54(1):264278.

    [49]" 郎軼詠,王強(qiáng),劉蘭娣,等.聚維酮碘溶液穩(wěn)定性的研究[J].中國(guó)醫(yī)院藥學(xué)雜志,2008,28(17):15231525.

    Lang Yiyong,Wang Qiang,Liu Landi,et al. Study on the Stability of Povidone Iodine Solution[J]. Chinese Journal of Hospital Pharmacy,2008,28(17):15231525.

    [50]" Sheppard M I,Thibault D H. Chemical Behaviour of Iodine in Organic and Mineral Soils[J]. Applied Geochemistry,1992,7(3):265272.

    [51]" Shimamoto Y S,Takahashi Y,Terada Y. Formation of Organic Iodine Supplied as Iodide in a Soil-Water System in Chiba,Japan[J]. Environmental Science amp; Technology,2011,45(6):20862092.

    [52]" 曹寒,張?jiān)拢饾?,?土壤中碘的賦存形態(tài)及遷移轉(zhuǎn)化研究進(jìn)展[J].巖礦測(cè)試,2022,41(4):521530.

    Cao Han,Zhang Yue,Jin Jie,et al. Iodine Speciation,Transportation,and Transformation in Soils:A Critical Review[J]. Rock and Mineral Analysis,2022,41(4):521530.

    [53]" Duborská E,Urík M,Bujdo M,et al. Influence of Physicochemical Properties of Various Soil Types on Iodide and Iodate Sorption[J]. Chemosphere,2019,214:168175.

    [54]" Hong C,Weng H,Jilani G,et al. Evaluation of Iodide and Iodate for Adsorption-Desorption Characteristics and Bioavailability in Three Types of Soil[J]. Biological Trace Element Research,2012,146(2):262271.

    [55]" Qian K,Li J,Chi Z,et al. Natural Organic Matter-Enhanced Ttransportation of Iodine in Groundwater in the Datong Basin:Impact of Irrigation Activities[J]. Science of the Total Environment,2020,730:138460.

    [56]" 包亞明.不同pH調(diào)節(jié)方法對(duì)聚維酮碘溶液穩(wěn)定性的影響[J].黑龍江醫(yī)藥,2008,21(1):6263.

    Bao Yaming. Effect of Different pH Adjustment Methods on the Stability of Povidone Iodine Solution[J]. Heilongjiang Medicine Journal,2008,21(1):6263.

    [57]" 鄧永強(qiáng),邢啟德,孫陽(yáng)陽(yáng),等.聚維酮碘的特性與臨床應(yīng)用[J].中華醫(yī)院感染學(xué)雜志,2016,26(3):627629.

    Deng Yongqiang,Xing Qide,Sun Yangyang,et al. Properties and Clinical Applications of Povidone Iodine[J]. Chinese Journal of Nosocomiology,2016,26(3):627629.

    [58]" 巴劍波,劉玉明,殷明,等.聚維酮碘特性及其制劑研究進(jìn)展[J].中國(guó)消毒學(xué)雜志,2010,27(1):6769.

    Ba Jianbo,Liu Yuming,Yin Ming,et al. Research Progress of Povidone Iodine Properties and Preparations[J]. Chinese Journal of Disinfection,2010,27(1):6769.

    [59]" Fuge R,Johnson C C. Iodine and Human Health, the Role of Environmental Geochemistry and Diet: A Review[J]. Applied Geochemistry,2015,63:282302.

    [60]" 張昌延,何江濤,張小文,等.珠江三角洲高砷地下水賦存環(huán)境特征及成因分析[J].環(huán)境科學(xué),2018,39(8):36313639.

    Zhang Changyan,He Jiangtao, Zhang Xiaowen,et al. Geochemical Characteristics and Genesis Analyses of High-Arsenic Groundwater in the Pearl River Delta[J]. Environmental Science,2018,39(8):36313639.

    [61]" Korobova E. Soil and Landscape Geochemical Factors Which Contribute to Iodine Spatial Distribution in the Main Environmental Components and Food Chain in the Central Russian Plain[J]. Journal of Geochemical Exploration,2010,107(2):180192.

    [62]" 趙仕琳,劉文靜,孫丹陽(yáng),等.大同盆地地下水中有機(jī)質(zhì)對(duì)碘遷移轉(zhuǎn)化的影響[J].地球科學(xué),2023,48(12):46994710.

    Zhao Shilin,Liu Wenjing,Sun Danyang,et al. Effect of Organic Matter on Iodine Mobilization in Groundwater of Datong Basin[J]. Earth Science,2023,48(12):46994710.

    [63]" Zhang X Z, Yin L Y, Chen J, et al. A Study of the Multi-Tiered Risk Assessment Method of Site Groundwater Contamination Considering Transport-Induced Risk[J]. Hydrogeology amp; Engineering Geology, 2023, 50(2): 160170.

    [64]" 資惠宇,龐園,張明珠.基于主成分分析法的廣州市從化區(qū)地下水質(zhì)量評(píng)價(jià)[J].人民珠江,2017,38(10):7276.

    Zi Huiyu,Pang Yuan,Zhang Mingzhu. Application of Principal Component Analysis for Conghua District Groundwater in Guangzhou City[J]. Pearl River,2017,38(10):7276.

    [65]" Shi J T,Liu Q,Zou Y R,et al. Shallow Groundwater Surrounding the Likeng Landfill,Guangzhou,China:Major Ions and Elements Indicating the Contamination Sources[J]. Water Science and Technology,2016,74(9):22342244.

    [66]" Chen J,Zhu A,Tang C,et al. Nitrogen Aspects of Hydrological Processes:A Case Study in Likeng Landfill,Guangzhou,China[J]. Environmental Science:Processes amp; Impacts,2014,16(11):26042616.

    [67]" Fuge R,Johnson C C,Phillips W J. An Automated Method for the Determination of Iodine in Geochemical Samples[J]. Chemical Geology,1978,23(1/2/3/4):255265.

    [68]" Umezawa Y,Hosono T,Onodera S,et al. Sources of Nitrate and Ammonium Contamination in Groundwater Under Developing Asian Megacities[J]. Science of the Total Environment,2008,404(2/3):361376.

    [69]" 韓龍喜,計(jì)紅,陸永軍,等.河道采沙對(duì)珠江三角洲水情及水環(huán)境影響分析[J].水科學(xué)進(jìn)展,2005,5(16):685690.

    Han Longxi,Ji Hong,Lu Yongjun,et al. Impact of Sand Excavation in Pearl River Delta on Hydrology and Environment[J]. Advances in Water Science,2005,5(16):685690.

    [70]" Voutchkova D D,Ernstsen V,Kristiansen S M,et al. Iodine in Major Danish Aquifers[J]. Environmental Earth Sciences,2017,76:116.

    [71]" Xue X,Li J,Xie X,et al. Impacts of Sediment Compaction on Iodine Enrichment in Deep Aquifers of the North China Plain[J]. Water Research,2019,159:480489.

    [72]" 張濤,王明國(guó),張智印,等.然烏湖流域地表水水化學(xué)特征及控制因素[J].環(huán)境科學(xué),2020,41(9):40034010.

    Zhang Tao,Wang Mingguo,Zhang Zhiyin,et al. Hydrochemical Characteristics and Possible Controls of the Surface Water in Ranwu Lake Basin[J]. Environmental Science,2020,41(9):40034010.

    [73]" Marghade D,Malpe D B,Rao N S,et al. Geochemical Assessment of Fluoride Enriched Groundwater and Health Implications from a Part of Yavtmal District, India[J]. Human and Ecological Risk Assessment,2019,26(3):122.

    [74]" Miao Q,Li X,Xu Y,et al. Chemical Characteristics of Groundwater and Source Identification in a Coastal City[J]. Plos One,2021,16(8):e0256360.

    [75]" 胡春華,周文斌,夏思奇.鄱陽(yáng)湖流域水化學(xué)主離子特征及其來(lái)源分析[J].環(huán)境化學(xué),2011,30(9):16201626.

    Hu Chunhua,Zhou Wenbin,Xia Siqi. Characteristics of Major Ions and the Influence Factors in Poyang Lake Catchment[J]. Environmental Chemistry,2011,30(9):16201626.

    [76]" 朱海勇,陳永金,劉加珍,等.塔里木河中下游地下水化學(xué)及其演變特征分析[J].干旱區(qū)地理,2013,36(1):818.

    Zhu Haiyong,Chen Yongjin,Liu Jiazhen,et al. Variation and Evolution of Groundwater Chemistry in the Middle and Lower Reaches of the Tarim River[J]. Arid Land Geography,2013,36(1):818.

    [77]" Liu W,Yang H,Li B,et al. Determination of Bromine and Iodine Speciation in Drinking Water Using High Performance Liquid Chromatography-Inductively Coupled Plasma-Mass Spectrometry[J]. Geostandards and Geoanalytical Research,2011,35(1):6974.

    [78]" Yang H,Liu W,Li B,et al. Speciation Analysis for Iodine in Groundwater Using High Performance Liquid Chromatography-Inductively Coupled Plasma-Mass Spectrometry (HPLCICPMS)[J]. Geostandards and Geoanalytical Research,2007,31(4):345351.

    [79]" Hu Q,Zhao P,Moran J E,et al. Sorption and Transport of Iodine Species in Sediments from the Savannah River and Hanford Sites[J]. Journal of Contaminant Hydrology,2005,78(3):185205.

    [80]" Dai J L,Zhang M,Zhu Y G. Adsorption and Desorption of Iodine by Various Chinese Soils[J]. Environment International,2004,30(4):525530.

    [81]" 張偉娜.環(huán)境水體中碘的存在形態(tài)及其影響因素研究[D].長(zhǎng)春:吉林大學(xué),2012.

    Zhang Weina. Studies on Speciation and Influencing Factors of Iodine in Environmental Water System[D]. Changchun:Jilin University, 2012.

    [82]" 周巾枚,蔣忠誠(chéng),徐光黎,等.崇左響水地區(qū)地下水水質(zhì)分析及健康風(fēng)險(xiǎn)評(píng)價(jià)[J].環(huán)境科學(xué),2019,40(6):26752685.

    Zhou Jinmei,Jiang Zhongcheng,Xu Guangli,et al. Water Quality Analysis and Health Risk Assessment for Groundwater at Xiangshui,Chongzuo[J]. Environmental Science,2019,40(6):26752685.

    [83]" 孫振中,劉金金,張玉平.上海市淡水養(yǎng)殖水體中氮、磷的分布研究[J].水產(chǎn)科學(xué),2018,37(4):512521.

    Sun Zhenzhong,Liu Jinjin,Zhang Yuping. Distribution of Nitrogen and Phosphorus in Waters of Freshwater Aquaculture in Shanghai[J]. Fisheries Science,2018,37(4):512521.

    [84]" 張玉璽,孫繼朝,陳璽,等.珠江三角洲淺層地下水pH值的分布及成因淺析[J].水文地質(zhì)工程地質(zhì),2011,38(1):1621.

    Zhang Yuxi,Sun Jichao,Chen Xi,et al. Characteristics and Preliminary Analyses of the Formation of pH in Shallow Groundwater in the Pearl River Delta[J]. Hydrogeology amp; Engineering Geology,2011,38(1):1621.

    [85]" 張明珠,王卓微,黎坤,等.廣州市淺層地下水化學(xué)時(shí)空分布特征及其對(duì)土地利用類型的響應(yīng)[J].生態(tài)環(huán)境學(xué)報(bào),2017,26(9):15391546.

    Zhang Mingzhu,Wang Zhuowei,Li Kun,et al. Spatial and Seasonal Geochemical Characteristics of Shallow Groundwater in Response to Land Use Pattern in Guangzhou[J]. Ecology and Environmental Sciences,2017,26(9):15391546.

    [86]" 米生雍,郭華東.異喹啉修飾金屬有機(jī)骨架配合物的合成及其二氧化碳和碘的吸附性能[J].吉林大學(xué)學(xué)報(bào)(理學(xué)版),2023,61(5):12301236.

    Mi Shengyong,Guo Huadong. Synthesis of MOF Complex Modified by lsoquinoline and Its Adsorption Properties for Carbon Dioxide and Lodine[J]. Journal of Jilin University (Earth Science Edition),2023,61(5):12301236.

    猜你喜歡
    研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    關(guān)于遼朝“一國(guó)兩制”研究的回顧與思考
    EMA伺服控制系統(tǒng)研究
    基于聲、光、磁、觸摸多功能控制的研究
    電子制作(2018年11期)2018-08-04 03:26:04
    新版C-NCAP側(cè)面碰撞假人損傷研究
    關(guān)于反傾銷會(huì)計(jì)研究的思考
    焊接膜層脫落的攻關(guān)研究
    電子制作(2017年23期)2017-02-02 07:17:19
    av.在线天堂| 国产精品日韩av在线免费观看| 日本av手机在线免费观看| 亚洲四区av| 欧美极品一区二区三区四区| 欧美日韩国产亚洲二区| 亚洲第一区二区三区不卡| 久久欧美精品欧美久久欧美| 黑人高潮一二区| 91精品伊人久久大香线蕉| 国产亚洲一区二区精品| 色5月婷婷丁香| 你懂的网址亚洲精品在线观看 | 国产亚洲精品久久久com| 啦啦啦啦在线视频资源| 18+在线观看网站| 日本黄大片高清| 日本wwww免费看| 青春草亚洲视频在线观看| 男人和女人高潮做爰伦理| 高清毛片免费看| 日韩成人伦理影院| 91久久精品电影网| 亚洲国产欧美在线一区| 一个人观看的视频www高清免费观看| 久久99热6这里只有精品| 免费人成在线观看视频色| 熟女电影av网| 美女黄网站色视频| 69av精品久久久久久| 非洲黑人性xxxx精品又粗又长| 大又大粗又爽又黄少妇毛片口| 夫妻性生交免费视频一级片| 亚洲精品影视一区二区三区av| 99国产精品一区二区蜜桃av| 少妇熟女欧美另类| 国产亚洲精品av在线| 午夜福利网站1000一区二区三区| 亚洲欧洲日产国产| 日韩亚洲欧美综合| 成人美女网站在线观看视频| 成人毛片60女人毛片免费| 亚洲美女视频黄频| 日日摸夜夜添夜夜爱| a级毛色黄片| 看片在线看免费视频| 国产精品一区二区性色av| 国产精品一及| 精品免费久久久久久久清纯| 嘟嘟电影网在线观看| 久久久久久久久中文| 午夜福利高清视频| 久久人人爽人人爽人人片va| 午夜爱爱视频在线播放| 亚洲人成网站高清观看| 成年女人看的毛片在线观看| 精品久久久久久久久av| 中文字幕精品亚洲无线码一区| 亚洲av中文字字幕乱码综合| 亚洲精品日韩在线中文字幕| 视频中文字幕在线观看| 九九爱精品视频在线观看| 国产黄片美女视频| 亚洲美女视频黄频| 久久国内精品自在自线图片| 久久久a久久爽久久v久久| 最近中文字幕2019免费版| 亚洲av男天堂| 成人无遮挡网站| 99久久成人亚洲精品观看| 亚洲av二区三区四区| 国产伦精品一区二区三区视频9| 99久久无色码亚洲精品果冻| 国产精品麻豆人妻色哟哟久久 | 最近最新中文字幕免费大全7| 亚洲va在线va天堂va国产| 真实男女啪啪啪动态图| 在线天堂最新版资源| 国产高潮美女av| 亚洲精品成人久久久久久| 久久久久久久久大av| 久久久久九九精品影院| 国产伦精品一区二区三区四那| 日韩av不卡免费在线播放| 亚洲最大成人av| 麻豆乱淫一区二区| 午夜视频国产福利| 国产精品综合久久久久久久免费| 久久久久久久久久黄片| 日韩高清综合在线| 天堂av国产一区二区熟女人妻| 九九在线视频观看精品| 国产精品一及| 99在线视频只有这里精品首页| 欧美色视频一区免费| 免费黄网站久久成人精品| 色5月婷婷丁香| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 一本久久精品| 亚洲人成网站在线播| 国产午夜精品论理片| 男的添女的下面高潮视频| 精品午夜福利在线看| 六月丁香七月| 亚洲国产精品国产精品| 免费看光身美女| 亚洲成人精品中文字幕电影| 淫秽高清视频在线观看| 久久精品久久久久久久性| 小蜜桃在线观看免费完整版高清| 久久精品人妻少妇| 亚洲精品国产成人久久av| 毛片女人毛片| 亚洲丝袜综合中文字幕| 日韩精品青青久久久久久| 床上黄色一级片| 美女cb高潮喷水在线观看| 成年av动漫网址| 欧美激情国产日韩精品一区| 婷婷色综合大香蕉| 99国产精品一区二区蜜桃av| 成人av在线播放网站| 乱码一卡2卡4卡精品| 伦理电影大哥的女人| 亚洲欧美精品专区久久| 一夜夜www| 免费黄网站久久成人精品| 女人久久www免费人成看片 | 日本三级黄在线观看| 欧美人与善性xxx| 2021天堂中文幕一二区在线观| 色综合站精品国产| 亚洲国产成人一精品久久久| 亚洲在久久综合| 国产av在哪里看| 国产亚洲最大av| 亚洲最大成人手机在线| 欧美另类亚洲清纯唯美| 黄色配什么色好看| 国内揄拍国产精品人妻在线| 日韩亚洲欧美综合| 国产欧美日韩精品一区二区| 中文天堂在线官网| 综合色av麻豆| 如何舔出高潮| 少妇高潮的动态图| 亚洲激情五月婷婷啪啪| 乱人视频在线观看| 久久久久久久国产电影| 亚洲成av人片在线播放无| 黄片无遮挡物在线观看| 岛国在线免费视频观看| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 九色成人免费人妻av| 精品一区二区三区视频在线| 亚洲成av人片在线播放无| 在线天堂最新版资源| 色综合站精品国产| 韩国av在线不卡| 国产精品不卡视频一区二区| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 中文亚洲av片在线观看爽| 久久久成人免费电影| 一夜夜www| 久久人妻av系列| 国产成人精品久久久久久| 黑人高潮一二区| 免费观看的影片在线观看| 久久久a久久爽久久v久久| 国产精品久久久久久久电影| 丰满乱子伦码专区| 中文字幕人妻熟人妻熟丝袜美| 插逼视频在线观看| av免费在线看不卡| 一本一本综合久久| 欧美变态另类bdsm刘玥| 人体艺术视频欧美日本| 夜夜看夜夜爽夜夜摸| 亚洲不卡免费看| 国产探花在线观看一区二区| 内射极品少妇av片p| 久久久久国产网址| 久久精品91蜜桃| 亚洲激情五月婷婷啪啪| 婷婷色综合大香蕉| 国产一区二区亚洲精品在线观看| 成人午夜高清在线视频| 91精品伊人久久大香线蕉| 欧美色视频一区免费| 亚洲一区高清亚洲精品| 乱码一卡2卡4卡精品| 97超碰精品成人国产| 天天一区二区日本电影三级| 亚洲av日韩在线播放| 神马国产精品三级电影在线观看| 少妇丰满av| 精品久久国产蜜桃| 男女视频在线观看网站免费| 日韩一区二区三区影片| 内地一区二区视频在线| 免费一级毛片在线播放高清视频| 欧美一区二区精品小视频在线| 中文字幕熟女人妻在线| 亚洲欧美日韩无卡精品| 美女xxoo啪啪120秒动态图| 国产一级毛片在线| 三级经典国产精品| 最近视频中文字幕2019在线8| 久久99热6这里只有精品| 天天一区二区日本电影三级| 搡女人真爽免费视频火全软件| 插阴视频在线观看视频| 亚洲不卡免费看| 麻豆成人av视频| 99热这里只有是精品50| 大香蕉97超碰在线| 亚洲人与动物交配视频| 亚洲国产成人一精品久久久| 成人欧美大片| 欧美激情久久久久久爽电影| 国产伦一二天堂av在线观看| 九九在线视频观看精品| 成人三级黄色视频| 精品久久久久久久久av| 午夜福利在线观看吧| 精品无人区乱码1区二区| 亚洲精品乱码久久久v下载方式| 一级毛片aaaaaa免费看小| 三级男女做爰猛烈吃奶摸视频| 2021天堂中文幕一二区在线观| 亚州av有码| 免费观看性生交大片5| 国产精品久久电影中文字幕| 日韩欧美 国产精品| 精品久久久久久久久av| 中文在线观看免费www的网站| 午夜福利在线观看免费完整高清在| 麻豆国产97在线/欧美| 成年免费大片在线观看| 网址你懂的国产日韩在线| 美女黄网站色视频| 国产69精品久久久久777片| 欧美人与善性xxx| 久久这里只有精品中国| 午夜激情福利司机影院| 成人高潮视频无遮挡免费网站| 亚洲国产欧美在线一区| 国产v大片淫在线免费观看| 网址你懂的国产日韩在线| 国产成年人精品一区二区| 久99久视频精品免费| 国产午夜精品论理片| 国产一区二区在线观看日韩| 一级毛片我不卡| 十八禁国产超污无遮挡网站| 日韩欧美精品v在线| 国产精品久久视频播放| 亚洲人成网站在线观看播放| 国产一区二区在线观看日韩| 久久久久网色| 国产黄片美女视频| 少妇的逼好多水| 建设人人有责人人尽责人人享有的 | 欧美色视频一区免费| 青春草亚洲视频在线观看| 成人性生交大片免费视频hd| 日韩国内少妇激情av| 美女国产视频在线观看| 国产乱来视频区| 亚洲av日韩在线播放| 国产精品久久久久久av不卡| 天堂av国产一区二区熟女人妻| 三级经典国产精品| 日本三级黄在线观看| 好男人在线观看高清免费视频| 晚上一个人看的免费电影| 伦精品一区二区三区| 身体一侧抽搐| 热99在线观看视频| 91久久精品国产一区二区成人| 国产乱来视频区| 久久这里有精品视频免费| 纵有疾风起免费观看全集完整版 | a级一级毛片免费在线观看| 1000部很黄的大片| 国产精品久久久久久精品电影| 欧美精品国产亚洲| 丰满乱子伦码专区| 热99re8久久精品国产| 国产精品永久免费网站| 1000部很黄的大片| 久久亚洲精品不卡| 午夜免费激情av| 可以在线观看毛片的网站| 男人舔奶头视频| 日日摸夜夜添夜夜爱| 久久人妻av系列| 国产精品一区二区性色av| 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 午夜视频国产福利| 国产成人91sexporn| 欧美色视频一区免费| 91精品一卡2卡3卡4卡| 精华霜和精华液先用哪个| 久久久久久伊人网av| 三级毛片av免费| 看片在线看免费视频| 日本午夜av视频| 能在线免费观看的黄片| 久久久久久久久久成人| 欧美bdsm另类| 男人和女人高潮做爰伦理| 亚洲无线观看免费| 婷婷色麻豆天堂久久 | 在线免费观看不下载黄p国产| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 综合色丁香网| 午夜福利成人在线免费观看| 人人妻人人澡人人爽人人夜夜 | 久热久热在线精品观看| 在线播放无遮挡| 国产 一区精品| 日产精品乱码卡一卡2卡三| 观看美女的网站| 久热久热在线精品观看| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| 国产久久久一区二区三区| 身体一侧抽搐| 亚洲内射少妇av| 国产精华一区二区三区| 日韩视频在线欧美| 精品国产一区二区三区久久久樱花 | 一卡2卡三卡四卡精品乱码亚洲| 一个人看视频在线观看www免费| 亚洲欧洲国产日韩| 春色校园在线视频观看| 久久精品国产鲁丝片午夜精品| 久久鲁丝午夜福利片| 精品午夜福利在线看| 国产精品熟女久久久久浪| 久久亚洲国产成人精品v| av国产免费在线观看| 欧美激情在线99| 99热这里只有精品一区| 欧美日韩综合久久久久久| 日韩大片免费观看网站 | 免费黄色在线免费观看| 成人一区二区视频在线观看| 精品人妻视频免费看| 国产在视频线精品| 国产成人freesex在线| 久久久久久伊人网av| 久久精品夜色国产| 欧美+日韩+精品| 亚洲精品日韩av片在线观看| 国产美女午夜福利| 欧美+日韩+精品| 久久精品久久久久久噜噜老黄 | 国产麻豆成人av免费视频| 一级av片app| 久久久久久久久中文| 村上凉子中文字幕在线| 国产免费视频播放在线视频 | 啦啦啦啦在线视频资源| 欧美日韩国产亚洲二区| 国产一区二区三区av在线| 国产免费视频播放在线视频 | 国产一区二区在线av高清观看| 欧美色视频一区免费| 国产在线男女| 汤姆久久久久久久影院中文字幕 | 又爽又黄a免费视频| 国产女主播在线喷水免费视频网站 | 午夜福利在线观看免费完整高清在| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 久久99热这里只有精品18| 久热久热在线精品观看| 波多野结衣高清无吗| 一二三四中文在线观看免费高清| 中文资源天堂在线| 18禁在线播放成人免费| 国产午夜精品一二区理论片| 国产精品,欧美在线| 蜜桃亚洲精品一区二区三区| 蜜臀久久99精品久久宅男| 亚洲欧美日韩高清专用| av.在线天堂| 99久久无色码亚洲精品果冻| 日韩,欧美,国产一区二区三区 | 色5月婷婷丁香| 男女边吃奶边做爰视频| 国产综合懂色| 久久久成人免费电影| 男人狂女人下面高潮的视频| 欧美日韩在线观看h| 成人av在线播放网站| 熟女电影av网| 亚洲综合精品二区| 国产白丝娇喘喷水9色精品| 国产精品久久久久久av不卡| 九草在线视频观看| 午夜亚洲福利在线播放| 欧美激情国产日韩精品一区| 日本免费a在线| 国产日韩欧美在线精品| 久久精品影院6| 午夜精品一区二区三区免费看| 菩萨蛮人人尽说江南好唐韦庄 | 91狼人影院| 国产精品日韩av在线免费观看| 级片在线观看| 亚洲综合色惰| videos熟女内射| 在线免费观看的www视频| 亚洲精品一区蜜桃| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片| 亚洲乱码一区二区免费版| 国产麻豆成人av免费视频| 男的添女的下面高潮视频| 亚洲综合精品二区| 国产乱人视频| 亚洲成人久久爱视频| 国产成人精品久久久久久| 日韩欧美精品免费久久| 欧美最新免费一区二区三区| 99久久无色码亚洲精品果冻| 国产一级毛片七仙女欲春2| 18禁裸乳无遮挡免费网站照片| 国产黄片美女视频| 一个人观看的视频www高清免费观看| 黄色配什么色好看| 两个人视频免费观看高清| 一级黄色大片毛片| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 18禁在线播放成人免费| 欧美激情国产日韩精品一区| 亚洲av电影在线观看一区二区三区 | 日本与韩国留学比较| 国产精品久久久久久精品电影小说 | 国产精品国产高清国产av| 美女xxoo啪啪120秒动态图| 看黄色毛片网站| 亚洲欧洲日产国产| 日韩一区二区视频免费看| 偷拍熟女少妇极品色| 97超碰精品成人国产| 嫩草影院精品99| 午夜免费男女啪啪视频观看| 中文资源天堂在线| 少妇的逼好多水| 91狼人影院| 我要搜黄色片| 久久久久久久久久成人| 搡老妇女老女人老熟妇| 不卡视频在线观看欧美| 两个人视频免费观看高清| 国产精品野战在线观看| 亚洲18禁久久av| 精品国内亚洲2022精品成人| 97在线视频观看| 日韩一区二区视频免费看| 欧美一区二区国产精品久久精品| 97超碰精品成人国产| 麻豆国产97在线/欧美| 伦精品一区二区三区| 嫩草影院新地址| 夫妻性生交免费视频一级片| 春色校园在线视频观看| 久久欧美精品欧美久久欧美| 欧美3d第一页| 亚洲成人中文字幕在线播放| 国产高清三级在线| 国产男人的电影天堂91| 国产淫语在线视频| 麻豆国产97在线/欧美| 内射极品少妇av片p| 一区二区三区四区激情视频| 看十八女毛片水多多多| 亚洲高清免费不卡视频| 精品国产露脸久久av麻豆 | 18禁在线播放成人免费| 国产伦理片在线播放av一区| 久久精品久久精品一区二区三区| 国产精品.久久久| av国产久精品久网站免费入址| 尾随美女入室| 一边摸一边抽搐一进一小说| 国产av码专区亚洲av| 丰满少妇做爰视频| 国产免费福利视频在线观看| 水蜜桃什么品种好| 国产高潮美女av| 热99re8久久精品国产| 国产三级中文精品| 桃色一区二区三区在线观看| 亚洲精品成人久久久久久| 久久人妻av系列| 久久人人爽人人爽人人片va| 国产又色又爽无遮挡免| 国产黄色小视频在线观看| 成人午夜精彩视频在线观看| 亚洲内射少妇av| 亚洲国产日韩欧美精品在线观看| 亚洲18禁久久av| 高清毛片免费看| 爱豆传媒免费全集在线观看| 亚洲三级黄色毛片| 国产高清国产精品国产三级 | 国产三级在线视频| 亚洲人与动物交配视频| 2021天堂中文幕一二区在线观| 直男gayav资源| 精品国产三级普通话版| 免费人成在线观看视频色| 久久久久久大精品| 亚洲成人av在线免费| 亚洲五月天丁香| 日韩av在线免费看完整版不卡| 午夜激情福利司机影院| 日韩中字成人| 蜜桃久久精品国产亚洲av| 亚洲国产精品专区欧美| 国产精品人妻久久久久久| 久久精品久久精品一区二区三区| 日韩,欧美,国产一区二区三区 | 国产老妇女一区| 三级毛片av免费| 免费看a级黄色片| 在线免费十八禁| 国产不卡一卡二| 国产亚洲av片在线观看秒播厂 | 亚洲自偷自拍三级| 国产高清三级在线| 日本爱情动作片www.在线观看| av福利片在线观看| 简卡轻食公司| 日产精品乱码卡一卡2卡三| 大香蕉97超碰在线| 国产老妇女一区| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看 | 麻豆乱淫一区二区| 在线观看一区二区三区| 国产激情偷乱视频一区二区| 深爱激情五月婷婷| 国产精品一区二区三区四区免费观看| 水蜜桃什么品种好| 97人妻精品一区二区三区麻豆| 亚洲图色成人| 精品一区二区免费观看| 国产又色又爽无遮挡免| 欧美一区二区亚洲| 大香蕉久久网| 日韩av在线大香蕉| 亚洲精品乱码久久久久久按摩| 国产单亲对白刺激| 久久久成人免费电影| 亚洲三级黄色毛片| 久久韩国三级中文字幕| 麻豆久久精品国产亚洲av| 99热6这里只有精品| 日韩大片免费观看网站 | 欧美bdsm另类| a级一级毛片免费在线观看| 青春草亚洲视频在线观看| 色综合色国产| 国产麻豆成人av免费视频| 国产 一区精品| 中文在线观看免费www的网站| 亚洲国产欧美人成| 亚洲无线观看免费| 日韩欧美精品免费久久| 91久久精品电影网| 精品久久久久久久久久久久久| 国产亚洲91精品色在线| www.av在线官网国产| 2021少妇久久久久久久久久久| 久久精品国产自在天天线| 亚洲久久久久久中文字幕| 男女下面进入的视频免费午夜| 日本av手机在线免费观看| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 日韩 亚洲 欧美在线| 91aial.com中文字幕在线观看| 免费在线观看成人毛片| 国产av在哪里看| 国产精品国产三级国产专区5o | 日本av手机在线免费观看| 亚洲成色77777| 国产老妇伦熟女老妇高清| 人人妻人人澡欧美一区二区| 午夜精品国产一区二区电影 | 欧美高清性xxxxhd video| 亚洲自偷自拍三级| 亚洲国产日韩欧美精品在线观看| 日韩av在线大香蕉| 身体一侧抽搐| 国产69精品久久久久777片| 菩萨蛮人人尽说江南好唐韦庄 | 午夜日本视频在线| 国产淫片久久久久久久久| 国产精品福利在线免费观看| 成年免费大片在线观看| 日韩欧美精品v在线|