暴雨徑流推移過程中顆粒粒徑分布與污染物含量關系

關鍵詞：水源水庫；徑流演變；水質(zhì)；顆粒粒徑

中圖分類號：TV213 文獻標識碼：A 文章編號：1001-9235（2024）12-0059-15

水庫在社會發(fā)展以及人民幸福生活起到重要的作用，具有防洪、供水以及水力發(fā)電等多項功能［1-2］，水庫的建成改變了原有水體的水力條件，庫區(qū)大壩將會截留流域內(nèi)的氮磷等營養(yǎng)鹽以及其他污染物質(zhì)，增加了水體富營養(yǎng)化以及污染物內(nèi)源釋放的可能［3-4］，進一步造成水質(zhì)惡化。降雨徑流會推動流域內(nèi)的物質(zhì)轉移［5-6］，攜帶大量污染物的徑流會蓄積于水庫中，從而造成水庫在降雨期間水質(zhì)急劇惡化，威脅到城市供水的安全性以及水庫的生態(tài)環(huán)境。

目前研究者對氣候變化以及人為活動對徑流的影響［7-8］，對徑流后造成水庫水質(zhì)惡化［9-12］等方面研究較多，對徑流的沿程推移以及水質(zhì)演變研究較少。在汛期雨水沖刷山體，大量泥沙顆粒物隨著雨水匯流進入河流湖泊，泥沙顆粒通常攜帶大量污染物，對河流湖泊帶來沖擊影響，造成水源水庫水質(zhì)急劇惡化。徑流期進行野外全程監(jiān)測徑流變化情況，有助于及時掌握徑流情況，及時進行調(diào)度管理，有助于提高供水安全。

徑流雨水中的氮磷大部分以顆粒態(tài)存在［13-14］，徑流雨水中TP、TN與SS之間的相關系數(shù)在0. 85以上［15］，徑流雨水中氮磷顆粒物是城市水環(huán)境的重要影響因素［16］。陳吉平等［17］分析了黑河1945年后的75 a的徑流數(shù)據(jù)，發(fā)現(xiàn)徑流變化的主要因素是氣候變化，師衛(wèi)釗等［18］研究了嘉陵江上游徑流演變影響因素，發(fā)現(xiàn)人為活動是影響徑流減少主導因素，氣候是次要因素。

由于徑流入庫水體與水庫原水體密度不同，入流水體在進入水庫后，會形成一個明顯高濁度水層，該水層即為潛流層。本研究以西安金盆水庫為研究對象，對2022年10月、2023年4月以及2023年6月的3次徑流進行監(jiān)測分析，探究3場徑流的演變情況以及顆粒物攜帶污染物含量與粒徑大小的關系，為汛期水庫調(diào)度提供相應的方案，提高供水安全。

1研究對象與方法

1. 1黑河金盆水庫概況

金盆水庫位于陜西省周至縣內(nèi)，供水量占西安市供水總量的75%。金盆水庫樞紐工程位于黑河峪口以上約1. 5 km處，距西安市86 km，是一項以城市供水為主，兼有農(nóng)灌、發(fā)電、防洪等綜合利用的大（2）型水利工程。樞紐由攔河壩、泄洪洞、溢洪洞、引水洞、壩后電站及古河道防滲工程等建筑物組成。水庫按百年一遇洪水標準（Q=3600 m³/s）設計，2000年一遇洪水（Q=6 400 m³/s）校核。正常高水位594.0m，汛限水位593.0m，設計、校核洪水位分別為594.34、597.18m。泄洪洞位于左岸，洞身為圓拱直墻城門洞形，斷面10m×13m，進口高程545m，最大泄洪量2539 m³/s，最大流速41 m/s，采用挑流消能。溢洪洞位于右岸，圓拱直墻城門洞形，斷面10.0m×11.5m，進口高程578 m，最大泄洪量2200m³/s，最大流速41 m/s，采用挑流消能。引水洞位于右岸，一洞多用，為城市、農(nóng)業(yè)及發(fā)電供水，設計引水流量30.3m³/s，引水口高程514.3m，圓形斷面，直徑為3.5m，發(fā)電支洞以后（弧門閘室后）為無壓洞，圓拱直墻城門洞形，斷面6.0m×6.8m，進口放水塔頂高程600m，塔身高為91.7m，沿塔高分設上、中、下3個取水孔分層取水，高程分別為571.0、554.0、514.3m。洞身為直徑3.5m的壓力圓洞。

1.2采樣點布置與水樣采集

從上游退水位到主庫區(qū)在河道中泓線處選取取樣點及監(jiān)測斷面共10個：其中選取S1入流斷面為上游河道入流口，S10監(jiān)測斷面為庫區(qū)，S2—S9過水斷面為過渡斷面，各監(jiān)測點基本相距1km，系統(tǒng)分析徑流期沿程變化特點，及潛流推移過程。監(jiān)測點位置見圖1。

在主庫區(qū)S10附近，連續(xù)進行水樣的采集，每周一次，在徑流期間視徑流情況增加采樣頻率；對上游沿程S1—S9號點進行取樣和在線監(jiān)測，采樣頻率為每月2次，徑流期間根據(jù)徑流頻率、潛流層變化及水質(zhì)情況增加采樣頻率；采樣采用2.5L有機玻璃采水器對取樣點進行取水，其中對水深較淺的S1—S5號點分別取表層（水面下0.5m左右）、中層（取樣深度約水深一半）、底層（據(jù)庫底0.5m左右）3個不同深度進行取樣和在線檢測；對S6—S8號點，自表層向下間隔20 m取水，底層仍為距庫底0.5m進行取樣和在線檢測；對于S9—S10號點，自表層（水面下0.5m左右）向下每隔10 m進行水樣采集與在線監(jiān)測。所取水樣均裝入聚乙烯瓶及時送回現(xiàn)場實驗室處理，48 h內(nèi)完成總氮（Total Nitrogen，TN）、總磷（Total Phosphorus，TP）、氨氮（NH3-N）、硝態(tài)氮（NO₃-N）、高錳酸鹽指數(shù)和DOC（Dissolved OrganicCarbon）的測定。在徑流期間，對S1—S10所有采樣點的潛流層采樣，采用2.5 L有機玻璃采水器，取潛流層水樣于25 L 聚乙烯朔料桶內(nèi)用于粒徑分析。

在采集水樣同時，對各點位水溫（T）和溶解氧（Dissolved Oxygen，DO）等水質(zhì)理化指標掩參數(shù)選用HACH Hydrolab DS5型多參數(shù)水質(zhì)測定儀（美國哈希公司）垂向間隔為1 ～ 5 m進行原位監(jiān)測。

1. 3水質(zhì)分析方法

各個監(jiān)測點的分析主要包括現(xiàn)場實時水質(zhì)在線監(jiān)測和實驗室分析2部分?，F(xiàn)場實時在線監(jiān)測儀器選用HACH Hydrolab DS5型水質(zhì)分析儀（美國，哈希），用來實時監(jiān)測水深（H）、水溫（T）、濁度（NTU）、葉綠素a（Chl-a）、溶解氧（DO）、電導率（σ）和氧化還原電位（Oxidation-Reduction Potential，ORP）等水體理化指標；水樣采集后，儲存于4 ℃的冰箱中保存。

對于水樣中的營養(yǎng)鹽指標，如氮（N）、磷（P）以及有機污染物等的測定，依據(jù)國家頒布的標準方法執(zhí)行，并且所有測試工作確保在48h內(nèi)完成。TN 和DTN（Dissolved Total Nitrogen）采用堿性過硫酸鉀消解紫外分光光度法，TP和DTP（Dissolved Total Phosphorus）采用鉬酸銨分光光度法，NH3-N采用納氏試劑分光光度法，NO₃-N采用紫外分光光度法，總有機碳（Total Organic Carbon，TOC）將水樣經(jīng)0. 45 μm WhatmanGF/F濾膜過濾后，溶解性有機碳采用有機碳分析儀（日本，島津）進行測定，含沙量SS則將500mL水樣采用過濾—烘干—稱量方法測定。粒徑采用激光粒度儀進行測定，對于潛流層水樣分別采用30、20、10、5、3μm 等孔徑的濾膜進行連續(xù)分級過濾，測定過膜前和過膜后的水中的總磷、總鐵、總錳等指標。

2結果與分析

2.1徑流演變情況

2.1.1金盆水庫降雨情況

水庫實時庫容和主庫區(qū)水位的高度直接受到出入庫流量變化的影響，見圖2。同時，水庫的運行調(diào)度也會因上游來水量的變化和城市用水需求的變動而受到影響。3次降雨過程中，10、4、6月累計降雨強度分別為94.00、56.05、72.65mm，且10月的降雨強度高，造成水庫庫容迅速提升，4月的單日雨強最大，但由于是2023年初期降雨，土壤含水率較低，因此導致此次降雨產(chǎn)生的徑流規(guī)模有限。

2.1.2徑流演變情況（2022年10月3日）

2022年10月3—5日黑河流域平均面降雨量累計94mm，入庫洪峰流量1138m3/s，最大雨強在10月4日達到36mm/d，同期泄洪3706萬m3，出入庫水量基本持平，水庫水位穩(wěn)定在592m附近。金盆水庫調(diào)度運行情況見表1。

由圖3可見，2022-10-03徑流來水的溫度和濁度分別約為14℃、1178NTU，2號點底部含沙量最大為2.0g/L。10月3日S1全層水體濁度、水溫、DO分別處于729～934NTU、14℃、10.82～14.02mg/L范圍內(nèi)，全層水體相關理化指標基本一致，表明暴雨徑流潛流在S1處于全斷面完全混合徑流。隨著S2—S4水深增加，斷面各水層濁度基本沿水深增加而升高，暴雨徑流來水與原水層水體摻混，高濁水從底部潛流推進。相較于表層水體相關理化指標濃度基本穩(wěn)定，上游沿程監(jiān)測點底層水體濁度、DO濃度出現(xiàn)明顯提升，S5表層水體濁度為4.1NTU，溶解氧含量為7.61mg/L，而底層水體濁度達到1178NTU，溶解氧含量提高至9.88mg/L，此次徑流在S5開始以底部潛流的形式向壩前推進，沿程各監(jiān)測點底層DO含量得到大幅提高。

金盆水庫是一座峽谷型水庫，其特點是河道的調(diào)蓄能力相對較弱［19］。因此，在徑流推移過程中，洪峰決定其潛流特征。徑流洪峰流量大、含沙量高和密度大［20］，其潛流位置通常位于較深的水層，在洪峰成形前后，入庫流量、含沙量以及入流密度相對較小，導致潛流位置隨著這些參數(shù)的降低而上升。洪峰攜帶的顆粒物是外源污染物和營養(yǎng)鹽輸入的主要途徑［21］。當高濁洪峰出現(xiàn)時，它會導致相應潛流層的污染負荷急劇增加，這一現(xiàn)象表現(xiàn)為攜沙洪峰攜帶高濃度的污染物，沿著河床向壩前推移。

此次徑流的含沙量較高，導致此次徑流的來水密度大于水庫原水，在整個徑流沿河道推移的過程是一個完全底層潛流的過程，從底部一直向前推進，此次從洪峰到達時就開始進行監(jiān)測，可以發(fā)現(xiàn)10月3日的監(jiān)測洪峰到達5號點，以5號點的濁度最高，底部最大濁度達到1178NTU，此次徑流過程5號點0～10m水深范圍內(nèi)的濁度為4.1～10.9NTU，說明徑流經(jīng)過5號點時就以底部潛流的形式向前推進。底部潛流一直向前影響到入庫口的9號點，9號點的底部濁度最大達到389NTU，且以507～520m濁度升高尤為明顯，均在130NTU以上。

在這場徑流事件中，水流攜帶著顯著高于常規(guī)水平的溶解氧，其濃度超過10mg/L。當富含溶解氧的暴雨徑流潛入水下時，它與潛流層及其鄰近的水層迅速混合。這一過程促進了暴雨徑流水與庫區(qū)原有水體之間的物質(zhì)和熱量交換。因此，高濃度的溶解氧迅速被傳輸?shù)綆靺^(qū)原有水體中，有效地替換了之前存在于底部的厭氧區(qū)域，顯著提高了底層水域的溶解氧含量。在沿程直到9號點，各層水體均達到5mg/L以上，在10月3日主庫區(qū)水層基本未受到潛流影響，溶解氧濃度分布與9月29日基本一致，底部仍處于厭氧狀態(tài)。前面根據(jù)濁度分析已經(jīng)發(fā)現(xiàn)5號點之后徑流是以底部潛流的形式向前推進，在5號點之后，中上層水體中的溶解氧含量基本與9月21日（圖4）一致，從側面可以進一步驗證此次徑流為底部潛流。

2.1.3徑流演變情況（2023年4月4日）

2023年4月2—4日黑河流域平均面降雨量累計56.05mm，最大降雨強度在4月2日達到42mm/d，入庫洪峰流量為174.7m3/s。黑河金盆水庫在該次降雨產(chǎn)生的徑流未采取泄洪，對上游所有來水均滯蓄在水庫中，水庫水位從3月28日的560.10m升至4月4日的567.78m，隨著持續(xù)蓄水，黑河金盆水庫在4月6日水位上升至573.01m。金盆水庫調(diào)度運行情況見表2。2023-04-04的徑流來水水溫約7.65℃，來水濁度約130NTU，水體含沙量小于0.1mg/L，由于來水濁度和水體含沙量較低，暴雨徑流來水密度主要取決于來水水溫，形成等溫度層潛流。

圖5所示，4月4日2號點全層水體濁度、水溫、溶解氧含量分別為121.3～131.4NTU、7.65～7.66℃、11.99～14.42mg/L范圍內(nèi)，全層水溫基本一致，表層和底層水體濁度差為10NTU，表明暴雨徑流潛流在S2處于全斷面混合徑流。從圖5沿程濁度變化可以看出2—4號點表底水層濁度差僅為10NTU左右，表明水體以完全混合的形式向前推進，到5號點表底水層濁度差達到95.7NTU，表明暴雨徑流來水從深水層通過，潛流在S5開始下潛，S4—S5處于全斷面混合徑流向底部潛流過渡階段。

隨著河道延伸，水深逐漸增加，至S6點達到43m。在前期水質(zhì)監(jiān)測中，S6點以下至550m高程的水層溫度維持在7.39～7.80℃，其水體密度與徑流來水的密度大致相等。因此，暴雨徑流潛流主要在高程558m以下的水層中向大壩方向推進。后續(xù)監(jiān)測點水深繼續(xù)增加，延緩上下水體熱量傳遞，在到達8號點時，520m高程以下的水溫降至7.58℃，低于徑流來水，原水庫水體密度相對徑流來水較大，徑流潛流在S8難以繼續(xù)下潛，最終在高程527～560m水層繼續(xù)向壩前行進，暴雨徑流潛流在S8形成中部潛流，繼續(xù)推移到S9就在水庫的稀釋作用下與水庫原水完全混合，未能進入到主庫區(qū)。

在徑流前3月21日的沿程監(jiān)測數(shù)據(jù)看出，在3月21日，沿程各點水溫基本在6.89～12.00℃，高程550m左右以下的水層水溫低于8℃，此次徑流來水溫度基本與520m高程左右的水溫一致，到8號點隨著水深的增加，底部水溫有所下降，從而導致在8號點開始從中部向9號點推進。

在4月4日，水庫剛開始進入分層期，在3月21日的沿程監(jiān)測發(fā)現(xiàn)各點溶解氧均在3mg/L以上，而且3月是水庫水位最低時期，上游來水能在2—4號點完全與原水體混合使得高濃度溶解氧含量的水迅速占據(jù)全斷面，而在5號點之后隨著水深增加，底部溶解氧在徑流前開始向厭氧方向發(fā)展，此次的高濃度溶解氧快速補充了由于分層底部內(nèi)源消耗的溶解氧，由于徑流潛入到9號點中部時濁度已經(jīng)降至20NTU左右，說明潛入到9號點的徑流量少，對9號點底部溶解氧補充有限。

2.1.4徑流演變情況（2023年6月3日）

2023年6月1—5日黑河流域平均面降雨量累計80.41mm，6月3日降雨強度最大達到34.31mm/d，6月3日平均入庫流量達到198.84m3/s，6月4日日平均入庫流量342.5m3/s。金盆水庫在6月3—5日進行泄洪。金盆水庫調(diào)度運行情況見表3。

2023-06-03的徑流來水水溫約12.06℃，來水濁度約153.6～215.0NTU，水體含沙量小于0.2mg/L，6月5日來水的水溫為11.77℃，來水濁度37.9～46.1NTU，水體含沙量小于0.05mg/L。

圖6所示，6月3日S1、S2全層水體濁度、水溫、DO等相關理化指標濃度梯度較低，表明暴雨徑流潛流在S1、S2處于全斷面混合徑流，還未表現(xiàn)出異重流運動特征。

6月3日的徑流過程中，隨著S3和S4位置的水深逐漸增加，暴雨徑流來水與表層水體開始逐漸分離。此時，高濁度的水體開始下潛，但表層和底層水體的濃度差異并不顯著。相較于表層水體相關理化指標濃度升高有限，上游沿程監(jiān)測點底層水體濁度較表層稍高，此時表明S3、S4是底層潛流的過渡段，S5表層水體濁度最高僅為7.8NTU，而底層水體濁度分別升高至52.2NTU，由于S3—S7水深相對較淺，此時河床底部高程還處于潛流層范圍內(nèi)，暴雨徑流潛流在S3—S7以底部潛流向壩前行進，進入7號點以后，在5月18日的沿程監(jiān)測中可以看出，底部水溫降至10℃左右，此時徑流來水的密度小于原有底層水體的密度，故徑流潛入到7號點后開始以中部潛流的形式向壩前推進，到達8號點，潛流層為高程540～570m的范圍內(nèi)。潛流層推進到入庫口9號點時，由于此次來水量較少，加之與之前的低濁水（5月18日沿程各層水體濁度基本在5NTU以下）混合稀釋作用，入庫口的濁度稀釋到15NTU以下，此時主庫區(qū)未受到此次徑流的影響。

隨后在6月5日的檢測中（圖7），來水的水溫為11.77℃，來水濁度為37.9～46.1NTU，水體含沙量小于0.05mg/L。6月5日S1、S2全層水體濁度、水溫等相關理化指標濃度梯度較低，表明暴雨徑流潛流在S1、S2處于全斷面徑流，尚未表現(xiàn)出異重流運動特征。

3—7號點表層10以上的水層濁度均小于20NTU，而底部水層的濁度基本在50NTU以上，由于S3—S7水深相對較淺，此時河床底部高程還處于潛流層范圍內(nèi)，暴雨徑流潛流在S3—S7以底部潛流向壩前行進，而進入8號點開始，底部水溫開始低于此次徑流來水水溫，加之此次的含沙量不高，從而導致此次徑流來水密度小于原水底層水體密度，7號點后開始以中部潛流的形式向壩前推進，到達8號點，潛流層為高程510～560m的范圍內(nèi)，較6月3日潛流層有所下降。潛流層推進到入庫口9號點時，潛流層上升5m，潛流層高程為515～560m，由于在6月4日，水庫進行了泄洪處理，加之連續(xù)降雨導致徑流來水量大幅增加，此次監(jiān)測到徑流影響到主庫區(qū)，從高程530m處開始潛入主庫區(qū)，加之之前揚水曝氣系統(tǒng)的混合作用，使得對主庫區(qū)影響范圍較廣，530～580m高程范圍的濁度達到20～40NTU。

2.2徑流沿程水質(zhì)變化特征

2.2.1總磷變化

氮和磷是水生生態(tài)系統(tǒng)中重要的營養(yǎng)鹽，它們作為判斷水體富營養(yǎng)化程度的主要指標，并被視為藻類生長和繁殖的關鍵影響因素［22-23］。由圖8可以看出，整個徑流過程中，總磷基本呈現(xiàn)沿程遞減的趨勢，且隨著徑流不斷向前推進，徑流來水與黑河金盆水庫原水不斷混合稀釋，潛流層的總磷含量在不斷降低，10月3日的徑流為底部潛流的形式進入水庫，1—3號點由于水深較淺，表層水體中磷含量受來水的影響較大，導致表層水體中的磷含量較高，到4號點后，徑流對表層水體的影響較小，表層水體中磷含量在4—10號點基本變化不大。在4月4日以及6月3日的沿程監(jiān)測中可以看出，徑流來水對水深較淺的前2個監(jiān)測點影響較大，在這2次的徑流中，來水量較少，攜帶的顆粒物較少，對后面各層水體的影響較小，4號點之后的水體中磷含量變化浮動較?。ㄆ渲?月4日磷含量在0.01～0.06mg/L范圍內(nèi)波動，6月3日磷含量在0.04～0.08mg/L范圍內(nèi)波動）。對于這2次的沿程監(jiān)測，4月4日監(jiān)測時徑流入庫流量為174.7m³/s，6月3日的監(jiān)測時入庫流量198.84m³/s，2次基本相同的入庫流量，且基本都是在較淺的水域先以底部潛流后以中部潛流的形式向壩前推進。6月3日的磷含量在入庫口的完全混合斷面明顯高于4月4日的完全混合斷面，但在后面徑流推進過程中，6月3日的徑流磷含量在3號點就快速削減，且在后續(xù)水體中磷含量波動幅度較小，相比之下，4月4日的削減沒有6月3日快速，且后續(xù)水體中的磷含量波動浮動大，體現(xiàn)在潛流層總磷含量較高，表層水體波動較小，推測基本相同的入庫流量下，庫容對外來污染負荷產(chǎn)生正相關關系，6月3日的水位達到594m，而4月4日水位僅有567m，大庫容的情況下，可以快速稀釋徑流帶來的高污染水，使得水體中磷含量上升較慢，可以大幅減少徑流對主庫去供水的影響。對比3次徑流，10月3日、4月4日、6月3日入庫流量分別為534.85、174.70、198.84m3/s，10月3日的入庫流量是后2次的3.0、2.6倍，而入庫的磷含量是后者的5倍，說明入庫的磷含量與入庫流量呈正相關關系，入庫流量越大攜帶的磷含量越高。

2.2.2總氮變化

暴雨徑流攜帶大量氮，使主庫區(qū)水體尤其是潛流層水體水質(zhì)迅速惡化，由圖9可以看出：10月3日徑流來水的氮含量在3.2mg/L左右，4月4日徑流來水的氮含量在4mg/L左右，總氮在總體上在隨徑流來水的推進呈現(xiàn)下降的趨勢，10月3日的徑流推移到7號點后水體中的總氮削減到2mg/L左右，7號點之后的氮削減量較少，一方面徑流推移到7號點之后水深加大，對徑流來水有一定的稀釋作用，另一方面通過含沙量數(shù)據(jù)可以看出，徑流來水的顆粒物經(jīng)過沿程的沉降稀釋等作用，7號點的含沙量快速降低，進一步削減了此時來水中的氮含量；在4月4日的徑流來水流量遠小于10月3日徑流來水量的情況下，混合斷面的總氮含量甚至超過徑流來水氮含量，推測為在4月4日，混合斷面水深太淺，此次由于水庫在3月水位處于全年最低值，2號點水深僅僅只有30cm左右，在船能極限到達2號點附近的時候，混合斷面水深僅為0.46m，此次徑流來水將河底的沉積物沖刷起來導致2號點混合斷面的氮含量急劇升高，而到3號點后，水深為3.37m，加上來水濁度含沙量均低于10月3日的徑流來水，3號點總氮含量降至2.7mg/L，在受徑流影響較弱的8—10號點，總氮含量基本與徑流前一致；6月3日的徑流來水總氮含量低于前期5月29日監(jiān)測的總氮含量，導致此次徑流中總氮含量沿程基本無變化。

2.2.3硝氮變化

圖10所示，2022年10月3日、2023年4月4日、2023年6月3日3次徑流來水中硝氮含量分別為1.8、2.5、1.5mg/L左右，在3次徑流中混合斷面中的硝氮含量以4月4日為最高，4月4日作為2023年第一場較為大的入庫徑流，將山體上的污染物質(zhì)沖刷進入水庫造成硝氮含量高于其他2場徑流的硝氮，4月4日的徑流過程中，2—4號點處于完全混合狀態(tài)，其硝氮含量基本趨于一致，隨后徑流開始以潛流的形式向前推進，表層的硝氮含量基本不受影響，潛流層硝氮上升明顯，且隨著潛流層的推進，在水庫的稀釋作用下，硝氮含量沿程遞減，硝氮含量在4月4日明顯隨著潛流層的變化而隨之改變，對比4月4日總氮含量變化，可以看出此次徑流的氮污染物主要以硝氮的形式存在于其中；10月3日的硝氮濃度整體隨潛流層向前推移，由于此次徑流來水硝氮濃度為1.8mg/L，在徑流前的一次沿程分析中，硝氮濃度在1.5mg/L左右，因此來水硝氮濃度較水庫原水并不算高，對水聲較淺的1—4號點硝氮濃度提升較高，進入5號點后，隨著潛流層顆粒物的沉降以及水庫原水的稀釋，5號點之后的硝氮濃度基本沿程保持一致，浮動變化不大，對比10月3日總氮含量變化，可以看出此次徑流的氮污染物主要以顆粒態(tài)氮的形式存在于其中，該種形式的氮可以隨顆粒的沉降而快速去除；6月3日的硝氮含量基本維持在1.5mg/L左右波動，這是由于在6月3日的徑流中，上游來水的硝氮含量并不高（1.5mg/L），這與徑流前的5月23日沿程監(jiān)測的沿程各斷面水體中5號點及其之后硝氮含量基本一樣，5月23日的沿程監(jiān)測中，1號點硝氮較低為1.1mg/L左右，但1—3號點基本為完全混合斷面，帶來的硝氮迅速提高原有的水體，從而使6月3日的徑流中硝氮含量沿程基本不變。

2.2.4氨變化

圖11所示，10月3日這次徑流中水中氨含量平均值在沿程上呈現(xiàn)一個遞減的趨勢，這由徑流在推移的過程中受到水庫稀釋導致，而在4號點氨氮濃度急劇升高，從圖3可以看出，4號點左右的濁度最大，徑流攜帶的污染物導致該點的氨濃度異常升高。在4月4日以及6月3日的2場徑流中，氨濃度在國家地表水環(huán)境質(zhì)量標準中Ⅱ類水體（0.5mg/L）的要求范圍內(nèi)，因此這兩場徑流入庫不會導致黑河金盆水庫氨含量超標的問題（該水庫氨濃度執(zhí)行國家地表水環(huán)境質(zhì)量標準中Ⅱ類水體的要求），在4月4日6號點出現(xiàn)氨濃度異常升高的現(xiàn)象，這由于在6號點外來人員釣魚活動多在此點，留下了大量的釣魚垃圾，加上4月4日為2023年初次徑流，在2022年11月到2023年4月初5個月的時間內(nèi)基本無降雨產(chǎn)流，在此次降雨沖刷下釣魚垃圾中的氨大量進入水體，造成6號點氨濃度異常升高的現(xiàn)象。在6月3日這場徑流中，來水氨濃度在0.2mg/L，徑流前的5月23日沿程監(jiān)測的沿程各斷面水體中水中氨濃度在0.11mg/L左右，與來水氨濃度差異不大，氨濃度在3號點之后沿程變化不明顯。從6月3日的濁度分析可以看出，此場徑流在1—2號點完全混合，濁度在100NTU以上，而到3號點之后，3號點的濁度減到50NTU，在隨后的潛流過程中，從3號點的濁度為50NTU緩慢削減到8號點的20NTU，從而造成此次徑流氨濃度在3號點之后沿程變化浮動不大的現(xiàn)象。對比3場徑流可以發(fā)現(xiàn)，徑流入庫流量越大，攜帶的氨污染物濃度越高，在入庫流量基本相同的情況下，初次徑流攜帶的氨污染物濃度越高，庫容越大的情況下對污染物的削減越快，對水庫水質(zhì)的影響也較小。

2.2.5高錳酸鹽指數(shù)變化

圖12所示，對比3場徑流可以發(fā)現(xiàn)，徑流入庫流量越大，攜帶的有機污染物濃度越高。10月3日以及4月4日這兩場徑流中，高錳酸鹽指數(shù)基本隨著徑流的推進而呈現(xiàn)沿程遞減的趨勢，且以潛流層的高錳酸鹽指數(shù)較高，這與徑流來水攜帶的高濃度有機污染物有關；6月3日的高錳酸鹽指數(shù)在完全混合斷面（1、2號點）高錳酸鹽指數(shù)較高，而在3號點之后的高錳酸鹽指數(shù)沿程變化趨勢不大，這與2.2.4節(jié)氨濃度變化趨勢基本相同，其原因不再贅述。3號點濁度為50NTU左右，2號點濁度為110NTU，3號點是2號點的一半，高錳酸鹽指數(shù)下降43%，由2.1節(jié)濁度變化圖可以看出，6月3日4號點之后沿程的濁度變化較小，而這之后的高錳酸鹽指數(shù)呈小幅下降，趨勢一致；10月3日以及4月4日的濁度大致沿程遞減，高錳酸鹽指數(shù)在沿程也呈現(xiàn)遞減趨勢，這說明高錳酸鹽指數(shù)與濁度的相關性較大。

2.3顆粒粒徑與污染物關系

2.3.13次徑流潛流層顆粒粒徑分布百分比

水中無機顆粒物含量較高，表面積很大，既能吸附許多有機離子，還能吸附許多水中的有機物，因此水中無機顆粒物將影響水中各種物質(zhì)的存在形態(tài)、遷移和轉化［24］。隨著顆粒粒徑的減小，單位體積的表面積增加，表面電荷增強，以及金屬吸附能力提高，這些細小顆粒所攜帶的污染物含量相應增加［25］。3場徑流潛流層顆粒粒徑分布百分比見表4—6。由表可知，在不同徑流強度下，徑流入庫所攜帶的顆粒物的粒徑也各不相同，在2022年10月3日的徑流中，入庫流量達到1000m³/s以上，所攜帶的顆粒物中10μm以上顆粒物占比較大，能達到70%，而在2023年6月3日的徑流中，入庫流量在300m3/s左右，此時超過30μm的顆粒物明顯低于2022年10月3日的徑流。因此徑流量越大，由于沖刷作用攜帶進入庫區(qū)的大粒徑顆粒物占比越大。

2.3.2顆粒粒徑與污染物含量關系

對3次徑流顆粒物上攜帶的鐵錳磷進行分析（圖13—15）：2022年10月3日的徑流量明顯高于2023年6月3日的徑流量，導致10月3日水庫入庫水中攜帶的顆粒物明顯高于6月3日，從而使10月3日的入庫水流中的污染物含量也明顯高于6月3日；其次，在3次徑流中由于都攜帶大量顆粒物，且攜帶的顆粒物均以大顆粒為主，也導致入庫水中大顆粒為主要的污染物來源；粒徑攜帶的污染物含量與粒徑所占百分比成一定相關性，該粒徑范圍在水中占比越大，其攜帶的污染物含量也越高；分析單位百分比顆粒物攜帶的污染物含量可知，小于5μm和小于3μm顆粒上攜帶的污染物含量最高，小顆粒擁有更大的比表面積，可以吸附攜帶更多的污染物。

分析顆粒粒徑與污染物含量關系圖可以看出，入庫口代表徑流攜帶的顆粒情況，懸浮顆粒（顆粒粒徑大于10μm）占總粒徑的70%以上，其中10月3日懸浮顆粒所攜帶的鐵錳磷含量超過總量的78%，4月4日懸浮顆粒鐵錳含量超過總量的70%，10月3日懸浮顆粒所攜帶的鐵錳含量超過總量的62%、磷含量超過總量的57%，說明徑流攜帶的顆粒中10μm以上為主要污染物來源。

對徑流后的水質(zhì)進行分析，大顆粒沉降速度快，在徑流后濁度下降較快，小顆粒沉降速度慢，濁度下降到30NTU左右后水體中的濁度削減較慢，且單位質(zhì)量小顆粒攜帶污染物含量大，導致水庫污染物含量依舊較高，可通過揚水曝氣系統(tǒng)加速水庫水體流動，從而提高沉降速度，改善水質(zhì)。

在汛期來水污染物以大顆粒帶來的污染物含量最大，高濁水一般以潛流的形式進入水庫，對于高濁度水，可以通過泄洪排除一部分潛流層的水，減小來水對水庫水質(zhì)的影響，使用高水位的取水口，從而提高供水質(zhì)量。

3結論與建議

a）徑流沿程推移過程以及進入庫區(qū)位置主要由徑流量大小影響。徑流在上游會歷經(jīng)全斷面混合徑流到底部潛流的演變過程，小徑流主要受來水水溫影響，在入庫前即完成溫度的混合，徑流沿等溫度層向前推進；大徑流由于泥沙含量大，徑流來水密度大于水庫原水，徑流會完全以底部潛流方式向前推移。

b）對比3次徑流，10月3日、4月4日、6月3日入庫流量分別為534.85、174.70、198.84m3/s，10月3日的入庫流量是后2次的3.0、2.6倍，而入庫的TP含量是后者的5倍，CODMn是后2次的1.8、1.7倍，說明入庫TP、CODMn與入庫流量呈明顯正相關關系；4月4日與6月3日入庫流量基本相同，但4月4日為初期徑流，其TN、NO₃-N分別為6月3日的2.2、1.7倍，高于10月3日TN、NO₃-N濃度，說明TN、NO₃-N受初期徑流影響較大。

c）徑流攜帶的顆粒物中大顆粒沉降速度快，在徑流后濁度下降較快，小顆粒沉降速度慢，濁度下降到30NTU左右后水體中的濁度削減較慢，且單位質(zhì)量小顆粒攜帶污染物含量大，導致水庫污染物含量依舊較高，可通過揚水曝氣系統(tǒng)加速水庫水體流動從而提高沉降速度以及提高生物降解速，從而改善水質(zhì)。

d）在汛期水庫水質(zhì)急劇惡化的情況下，針對徑流來水可以采取上取水口進行取水，徑流進入庫區(qū)，及時采用泄洪洞快速排出徑流攜帶的高濁水，達到排濁蓄清，保證上層水體水質(zhì)較小受到影響。在后續(xù)研究中繼續(xù)深入分析顆粒粒徑與其他污染物如氮、COD之間的聯(lián)系，以及徑流后主庫區(qū)水質(zhì)變化情況。