秦歡歡，黃麗想，陳益平，高 柏，孫占學(xué),

(1.東華理工大學(xué)核資源與環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室，江西 南昌 330013；2.東華理工大學(xué)水資源與環(huán)境工程學(xué)院，江西 南昌 330013)

隨著我國(guó)經(jīng)濟(jì)發(fā)展和人口增長(zhǎng)，地表水污染已經(jīng)成為全社會(huì)不可忽視的重要生態(tài)環(huán)境問(wèn)題。受污染水體是許多疾病發(fā)生的重要誘發(fā)因素，飲水不良導(dǎo)致世界上約80%的疾病和約20%的癌癥病例[1-2]，我國(guó)約90%的癌癥病例是由于患者與致癌物質(zhì)直接接觸所致[3-4]。此外，全世界每年有超過(guò)80萬(wàn)人的死亡原因是飲用和接觸受污染水所致，亞、非、拉3大洲每年有近350萬(wàn)人的死亡原因是患上與水相關(guān)的疾病[5]。因此，水環(huán)境質(zhì)量與人體健康息息相關(guān)，對(duì)水環(huán)境的改善是人類社會(huì)的共同目標(biāo)。然而，國(guó)內(nèi)一般多關(guān)注水質(zhì)超標(biāo)情況，關(guān)于水質(zhì)指標(biāo)對(duì)人體健康的影響與評(píng)價(jià)則缺少關(guān)注[6-7]。與傳統(tǒng)水質(zhì)等級(jí)評(píng)價(jià)相比，水質(zhì)健康風(fēng)險(xiǎn)評(píng)價(jià)能更直觀地定量描述與評(píng)價(jià)水體各類污染物質(zhì)對(duì)人體健康可能產(chǎn)生的潛在風(fēng)險(xiǎn)與危害[8-10]，因而越來(lái)越受到全世界環(huán)境健康研究人員的重視。

一般來(lái)說(shuō)，水體中對(duì)人體健康有害的物質(zhì)包括化學(xué)致癌物、放射性物質(zhì)等基因毒物質(zhì)和非致癌物等軀體毒物質(zhì)[11-13]，人類通過(guò)直接攝入、口鼻吸入和皮膚接觸3種方式接觸這些有害物質(zhì)，而飲水是人體暴露于水體污染物的重要途徑之一[14]。相比于其他2種有毒物質(zhì)，化學(xué)致癌物對(duì)人體產(chǎn)生的健康風(fēng)險(xiǎn)和危害最大。As和Cd廣泛存在于水體、土壤和空氣中，其作為重金屬元素具有的致癌性、持久性、可富集和可放大等特性[15-17]使得它們可以通過(guò)許多不同途徑進(jìn)入人體，并對(duì)多個(gè)組織和器官產(chǎn)生不可逆的毒害作用[18-19]，容易引起多種心血管疾病和肺癌、肝癌等癌癥[20-21]，故而受到研究人員廣泛關(guān)注。目前，國(guó)內(nèi)有許多針對(duì)河流水體中砷、鎘和其他重金屬健康風(fēng)險(xiǎn)評(píng)價(jià)案例和研究成果[1,11-13,16-17,22]，這些研究采用美國(guó)環(huán)境保護(hù)署(US EPA)水環(huán)境健康風(fēng)險(xiǎn)評(píng)價(jià)模型，取得了較好的評(píng)價(jià)效果，展現(xiàn)了該方法的科學(xué)性和有效性。

拉薩河地處青藏高原，是西藏經(jīng)濟(jì)發(fā)展的“大動(dòng)脈”，其生態(tài)環(huán)境具有多變性和脆弱性。近年來(lái)受到沿線地?zé)崤c礦產(chǎn)開(kāi)發(fā)、城鎮(zhèn)化發(fā)展、人口急劇增長(zhǎng)等因素的劇烈影響，拉薩河遭受不可忽視的污染，水體中As和Cd等化學(xué)致癌物對(duì)流域居民健康造成潛在風(fēng)險(xiǎn)和危害[23-25]。然而，針對(duì)拉薩河水體中As和Cd健康風(fēng)險(xiǎn)的評(píng)價(jià)較少[26]。通過(guò)對(duì)拉薩河中下游和堆龍曲支流16個(gè)采樣點(diǎn)水樣進(jìn)行分析，查明拉薩河流域水體As和Cd濃度與空間分布特征，采用健康風(fēng)險(xiǎn)評(píng)價(jià)模型估算由飲水途徑所致的As和Cd健康風(fēng)險(xiǎn)，定量闡述As和Cd與人體健康之間關(guān)系，為拉薩河流域水質(zhì)安全保障、沿線居民健康風(fēng)險(xiǎn)防控和流域可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

拉薩河全長(zhǎng)為551 km，流域面積為32 588 km2，處于岡底斯-念青唐古拉地質(zhì)構(gòu)造帶中，呈現(xiàn)南北高中間低的地勢(shì)特點(diǎn)[27]。拉薩河谷基巖以中生代石灰?guī)r和喜瑪拉雅期花崗巖為主，以燕山期花崗巖和第三紀(jì)、第四紀(jì)火山噴出巖、碎屑巖為輔。拉薩河流域平均海拔為5 400 m，由高山和中高山河谷2種地貌類型組成，日照充足，氣溫較低，年均降水量和蒸發(fā)量分別為460和1 217 mm[23]。拉薩河流域分為上、中和下游，其中，上游是人口較少的牧區(qū)，中游耕地面積較小，但人類活動(dòng)相對(duì)較多，下游是農(nóng)業(yè)發(fā)達(dá)且人口密度最大的地區(qū)。拉薩河流域擁有十分豐富的礦產(chǎn)資源，包括有色金屬、地?zé)岷妥匀涣虻雀黝惖V產(chǎn)，大中型礦山有12個(gè)，其他礦山有57個(gè)。然而，研究區(qū)大部分選礦場(chǎng)建在拉薩河兩岸，從這些選礦場(chǎng)排出的廢水只經(jīng)過(guò)簡(jiǎn)單處理即排入拉薩河，對(duì)拉薩河生態(tài)環(huán)境產(chǎn)生嚴(yán)重影響。

1.2 樣品采集方法

于2017年7月23—29日沿著拉薩河流向選擇16個(gè)采樣點(diǎn)，其中，中游有6個(gè)(SP1～SP6)，下游有7個(gè)(SP7～SP13)，堆龍曲支流有3個(gè)(SP14～SP16)，利用瞬時(shí)采樣法對(duì)河水進(jìn)行采樣，水樣采集器采用預(yù)清洗過(guò)的高密度聚乙烯瓶，采樣點(diǎn)分布見(jiàn)圖1。由于上游人口稀少，人類活動(dòng)較弱，因此將采樣點(diǎn)SP1作為“源頭”，其他采樣點(diǎn)距SP1的距離為該采樣點(diǎn)沿拉薩河的距離。采樣時(shí)先潤(rùn)洗采集器2～3次，隨后逆水流方向進(jìn)行水樣采集，采用0.45 μm孔徑濾膜對(duì)水樣進(jìn)行過(guò)濾，然后加入6 mol·L-1超純硝酸將水樣酸化至pH<2后，包裹好水樣瓶，確保不泄露，同時(shí)做好標(biāo)記，樣品送至東華理工大學(xué)分析測(cè)試中心進(jìn)行測(cè)試。水樣經(jīng)預(yù)處理后，采用電感耦合等離子體質(zhì)譜儀(ICP-MS，Thermo Element 2型)，按照HJ 700—2014《水質(zhì) 65種元素的測(cè)定 電感耦合等離子體質(zhì)譜法》中標(biāo)準(zhǔn)方法對(duì)As和Cd濃度進(jìn)行測(cè)定。

圖1 研究區(qū)位置、水系和采樣點(diǎn)示意

1.3 健康風(fēng)險(xiǎn)評(píng)價(jià)模型及參數(shù)

采用US EPA確定的健康風(fēng)險(xiǎn)評(píng)價(jià)模型，以研究區(qū)不同年齡組人群為研究對(duì)象，考慮經(jīng)飲水途徑攝入的化學(xué)致癌物As和Cd對(duì)人體產(chǎn)生的致癌風(fēng)險(xiǎn)，具體計(jì)算公式為

(1)

式(1)中，Rc為As和Cd通過(guò)飲水途徑導(dǎo)致個(gè)人罹患癌癥的年總風(fēng)險(xiǎn)，a-1；Rci為化學(xué)致癌物i通過(guò)飲水途徑導(dǎo)致個(gè)人罹患癌癥的年風(fēng)險(xiǎn)， a-1；i為1和2分別表示As和Cd；Dci為經(jīng)飲水途徑攝入的化學(xué)致癌物i的單位體重日均暴露劑量，mg·kg-1·d-1；Qci為經(jīng)飲水途徑攝入的化學(xué)致癌物i的致癌強(qiáng)度系數(shù)，kg·d·mg-1；70為人均壽命，a；Cci為化學(xué)致癌物i質(zhì)量濃度，mg·L-1；UT為年齡組T人群日均飲水量，L·d-1；WT為年齡組T人群平均體重，kg。根據(jù)沈威等[28]和齊文[29]的研究以及原核工業(yè)部的推薦方法，將人群劃分為幼兒(≤7歲)、少年(>7～17歲)和成年(>17歲)，這3個(gè)年齡組人群日均飲水量分別為1.10、1.37和2.00 L·d-1[28-29]，平均體重分別為14.90、40.71和60.33 kg[30]。As和Cd的致癌強(qiáng)度系數(shù)分別為15和6.1 kg·d·mg-1[3]。

2 結(jié)果與討論

2.1 研究區(qū)As和Cd濃度及分布特征

拉薩河水體As和Cd濃度統(tǒng)計(jì)結(jié)果見(jiàn)表1。如表1所示，As質(zhì)量濃度范圍為0.650～4.270 μg·L-1，平均值為(2.280±0.969) μg·L-1，最大值和最小值相差5.57倍，范圍跨度較大；Cd質(zhì)量濃度范圍為0.067～0.134 μg·L-1，平均值為(0.092±0.017) μg·L-1，最大值和最小值相差1倍，范圍跨度較小。As和Cd濃度均未超過(guò)GB 5749—2022《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》規(guī)定限值，也均達(dá)到GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中I類和Ⅱ類水質(zhì)標(biāo)準(zhǔn)，說(shuō)明拉薩河水體水質(zhì)總體處于良好水平，可用于集中式生活飲用水水源。變異系數(shù)是衡量污染物濃度平均變異程度的指標(biāo)，變異系數(shù)大于50%表明濃度分布不均勻，可能存在因外源物質(zhì)流入導(dǎo)致的點(diǎn)源污染。研究區(qū)As和Cd變異系數(shù)分別為42.5%和18.9%，均小于50%，表明研究區(qū)As和Cd濃度分布相對(duì)均勻，拉薩河沿程出現(xiàn)As和Cd點(diǎn)源污染的可能性很小。

表1 拉薩河水體中As和Cd濃度統(tǒng)計(jì)

與國(guó)內(nèi)不同河湖中As和Cd平均濃度(表2[11,26,31-42])相比而言，拉薩河水體As和Cd濃度總體處于較低水平，但As濃度高于灤河、黃浦江(上游)、赤水河和贛江，Cd濃度高于漓江和九龍江。劉鳳等[26]在2010年4月測(cè)得拉薩河干流以及主要支流As和Cd平均質(zhì)量濃度分別為27.700和0.059 μg·L-1，As濃度顯著高于筆者研究結(jié)果，而Cd濃度則略低。

表2 國(guó)內(nèi)不同水體As和Cd平均濃度[11,26,31-42]

拉薩河水體中As和Cd濃度沿程變化情況見(jiàn)圖2。如圖2所示，拉薩河中下游和堆龍曲支流水體As濃度均沿程增長(zhǎng)，但呈不同趨勢(shì)。中游As濃度(y)與沿程距離(x)呈一定程度正相關(guān)關(guān)系，相關(guān)性方程為y=0.012 7x+1.765 3 (R2=0.323 2,P<0.05)，沿程距離每增加100 km，ρ(As)增加1.27 μg·L-1；下游As濃度與沿程距離呈較好正相關(guān)關(guān)系，相關(guān)性方程為y=0.014 9x+0.307 2 (R2=0.545 5,P<0.01)，沿程距離每增加100 km，ρ(As)增加1.49 μg·L-1；堆龍曲支流As濃度與沿程距離呈較弱正相關(guān)關(guān)系，相關(guān)性方程為y=0.011 6x-1.626 4 (R2=0.181 0,P<0.05)，沿程距離每增加100 km，ρ(As)增加1.16 μg·L-1。而Cd濃度則在中下游和堆龍曲支流呈不同變化趨勢(shì)，中游Cd濃度與沿程距離呈較強(qiáng)負(fù)相關(guān)關(guān)系，相關(guān)性方程為y=-0.000 3x+0.096 1 (R2=0.814 4,P<0.01)，沿程距離每增加100 km，ρ(Cd)下降0.03 μg·L-1；下游Cd濃度與沿程距離呈弱負(fù)相關(guān)關(guān)系，相關(guān)性方程為y=-0.000 2x+0.130 7 (R2=0.254 3,P<0.05)，沿程距離每增加100 km，ρ(Cd)下降0.02 μg·L-1；堆龍曲支流Cd濃度與沿程距離呈弱正相關(guān)關(guān)系，相關(guān)性方程為y=0.000 3x+0.019 2 (R2=0.2095,P<0.05)，沿程距離每增加100 km，ρ(Cd)增加0.03 μg·L-1。

圖2 拉薩河水體As和Cd濃度沿程分布曲線

拉薩河As和Cd濃度與采樣點(diǎn)經(jīng)緯度之間的關(guān)系見(jiàn)圖3。在地理分布上，拉薩河水體As濃度從南到北(緯度增加)整體呈下降趨勢(shì)，Cd濃度呈微弱上升趨勢(shì)，緯度每增加1°，ρ(As)減少1.778 μg·L-1，ρ(Cd)增加0.000 2 μg·L-1；拉薩河水體As濃度從西到東(經(jīng)度增加)整體呈上升趨勢(shì)，Cd濃度則呈下降趨勢(shì)，經(jīng)度每增加1°，ρ(As)增加0.456 μg·L-1，ρ(Cd)則減少0.010 7 μg·L-1。

圖3 拉薩河水體As和Cd濃度與經(jīng)緯度變化關(guān)系

2.2 As和Cd的健康風(fēng)險(xiǎn)評(píng)價(jià)

根據(jù)研究區(qū)水體As和Cd濃度數(shù)據(jù)，采用健康風(fēng)險(xiǎn)評(píng)價(jià)模型及相關(guān)參數(shù)，計(jì)算得到研究區(qū)各采樣點(diǎn)水體As和Cd經(jīng)由飲水途徑對(duì)不同年齡人群的健康風(fēng)險(xiǎn)，結(jié)果見(jiàn)表3。如表3所示，3個(gè)人群致癌總風(fēng)險(xiǎn)排序?yàn)橛變?3.647×10-5a-1)>少年(1.670×10-5a-1)>成年(1.645×10-5a-1)，幼兒經(jīng)飲水途徑的致癌總風(fēng)險(xiǎn)分別為少年和成年的2.18和2.22倍。這是因?yàn)橛變壕哂斜壬倌旰统赡旮用舾械纳眢w結(jié)構(gòu)，更容易受到化學(xué)致癌物的影響，因此其對(duì)幼兒的危害也更大，這與文獻(xiàn)[11]的研究結(jié)果一致。除幼兒組在采樣點(diǎn)SP4、SP12和SP13的總致癌風(fēng)險(xiǎn)超過(guò)國(guó)際輻射防護(hù)委員會(huì)(ICRP)推薦的最大可接受值(5×10-5a-1)外，其余各年齡組在不同采樣點(diǎn)的總致癌風(fēng)險(xiǎn)均低于US EPA推薦的最大可接受值(1×10-4a-1)和ICRP推薦值，說(shuō)明拉薩河水體由As和Cd導(dǎo)致的致癌風(fēng)險(xiǎn)總體上處于安全水平，但需要對(duì)幼兒在飲水安全方面給予更多關(guān)注。

表3 研究區(qū)各采樣點(diǎn)As和Cd的個(gè)人年風(fēng)險(xiǎn)水平

表3顯示，拉薩河水體化學(xué)致癌物經(jīng)飲水途徑對(duì)不同年齡人群的致癌風(fēng)險(xiǎn)在1×10-7～1×10-5范圍內(nèi)，As致癌風(fēng)險(xiǎn)值大于Cd。As和Cd致癌風(fēng)險(xiǎn)值范圍分別為0.462×10-5～6.714×10-5和0.019×10-5～0.086×10-5a-1，總致癌風(fēng)險(xiǎn)值范圍為0.485×10-5～6.767×10-5a-1。As經(jīng)飲水途徑的個(gè)人致癌健康風(fēng)險(xiǎn)最大值為6.714×10-5a-1，是ICRP推薦值的1.34倍。對(duì)拉薩河流域不同年齡人群，As致癌風(fēng)險(xiǎn)占As和Cd總致癌風(fēng)險(xiǎn)的比例大于98.3%，而Cd的致癌風(fēng)險(xiǎn)占比小于1.7%，表明拉薩河流域As致癌風(fēng)險(xiǎn)占據(jù)主導(dǎo)地位，貢獻(xiàn)了絕大多數(shù)致癌健康風(fēng)險(xiǎn)。由此可見(jiàn)，與Cd相比，As為拉薩河流域主要致癌污染物，對(duì)拉薩河流域進(jìn)行管理時(shí)需加強(qiáng)對(duì)As濃度的監(jiān)測(cè)，確保As引起的致癌風(fēng)險(xiǎn)處于可接受范圍。

從筆者研究結(jié)果可知，相比于國(guó)內(nèi)其他水體來(lái)說(shuō)，拉薩河流域水體As和Cd濃度水平較低(表2)，均不超過(guò)國(guó)家規(guī)定相關(guān)水質(zhì)標(biāo)準(zhǔn)限值，適合作為飲用水水源。相關(guān)致癌健康風(fēng)險(xiǎn)評(píng)價(jià)計(jì)算結(jié)果表明，As經(jīng)飲水途徑導(dǎo)致的致癌風(fēng)險(xiǎn)最大值略超過(guò)ICRP推薦的最大可接受風(fēng)險(xiǎn)值，需要給予一定關(guān)注和重視。因此，對(duì)于拉薩河流域水體水質(zhì)，在進(jìn)行常規(guī)監(jiān)測(cè)時(shí)需結(jié)合健康風(fēng)險(xiǎn)評(píng)價(jià)進(jìn)行管理，這有利于更全面地了解拉薩河作為飲用水水源地的環(huán)境質(zhì)量與健康安全狀況，以便有針對(duì)性地進(jìn)行風(fēng)險(xiǎn)管理，制定污染防治和治理策略。

2.3 與國(guó)內(nèi)其他水體比較與分析

采用健康風(fēng)險(xiǎn)評(píng)價(jià)模型計(jì)算得到國(guó)內(nèi)不同水體As和Cd對(duì)不同年齡人群經(jīng)飲水途徑的致癌風(fēng)險(xiǎn)見(jiàn)表4[11,26,31-42]。表4顯示，與國(guó)內(nèi)其他水體相比，拉薩河流域As和Cd總致癌風(fēng)險(xiǎn)(1.645×10-5～3.647×10-5a-1)處于較低水平。在幼兒、少年和成年3組人群中，拉薩河流域As和Cd總致癌風(fēng)險(xiǎn)在15條河流中均排第13位，僅高于灤河(1.470×10-5～3.260×10-5a-1)和赤水河(1.293×10-5～2.867×10-5a-1)。由此可見(jiàn)，拉薩河流域As和Cd經(jīng)飲水途徑的致癌風(fēng)險(xiǎn)總體處于可接受的安全水平。

表4 國(guó)內(nèi)不同水體As和Cd的平均個(gè)人年風(fēng)險(xiǎn)水平[11,26,31-42]

3 結(jié)論

(1)拉薩河流域水體As和Cd平均質(zhì)量濃度分別為(2.280±0.969)和(0.092± 0.017) μg·L-1，未超過(guò)GB 5749—2022規(guī)定限值，也達(dá)到GB 3838—2002中Ⅰ類和Ⅱ類水質(zhì)標(biāo)準(zhǔn)，變異系數(shù)分別為42.5%和18.9%，空間波動(dòng)性較小，存在As和Cd點(diǎn)源污染的可能性很小。與國(guó)內(nèi)其他水體相比，拉薩河水體As和Cd濃度處于較低水平。

(2)拉薩河中下游和堆龍曲支流As濃度均呈一定程度沿程增長(zhǎng)趨勢(shì)，Cd濃度在拉薩河中下游呈沿程下降趨勢(shì)，在堆龍曲支流呈沿程增長(zhǎng)趨勢(shì)。在地理分布上，隨著緯度增加，拉薩河水體As濃度呈下降趨勢(shì)，Cd濃度呈微弱上升趨勢(shì)；隨著經(jīng)度增加，拉薩河水體As濃度呈上升趨勢(shì)，Cd濃度呈下降趨勢(shì)。

(3)拉薩河流域水體As和Cd健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表明，與Cd相比，As為拉薩河流域水體起主導(dǎo)作用的化學(xué)致癌物，貢獻(xiàn)了超過(guò)98%的致癌風(fēng)險(xiǎn)，應(yīng)成為拉薩河流域水質(zhì)管理的優(yōu)先管控對(duì)象。拉薩河流域As和Cd經(jīng)飲水途徑的致癌風(fēng)險(xiǎn)總體處于可接受的安全水平。相比于少年和成年，幼兒為更敏感的健康風(fēng)險(xiǎn)受體，受As和Cd的危害更大，需要在幼兒飲水安全方面給予更多關(guān)注。