三峽水庫干流含沙量對不同形態(tài)磷濃度的影響分析

摘要：磷是水環(huán)境質(zhì)量評價和污染控制考核的重要指標，其在天然水體中以溶解態(tài)、顆粒態(tài)等形態(tài)存在。研究水中含沙量的改變對不同形態(tài)磷濃度的影響，對磷污染控制有著重要意義。選取三峽水庫干流為研究對象，采用Mann-Kendall非參數(shù)檢驗法和Pearson相關(guān)系數(shù)分析了該河段含沙量變化對總磷、溶解磷、顆粒磷等不同形態(tài)磷濃度的影響。結(jié)果表明：① 2003～2019年三峽水庫干流的平均流量無明顯變化，含沙量則明顯減少。② 總磷濃度在2003～2012年間起伏波動、緩慢上升，2013～2019年則持續(xù)下降；溶解磷濃度呈先上升后下降的趨勢，各斷面峰值出現(xiàn)在2013、2014年；顆粒態(tài)磷濃度呈波動下降趨勢。③ 與2003～2012年濃度均值相比，2013～2019年渾濁水樣磷濃度、沉降水樣磷濃度以及顆粒磷濃度分別下降了40%，25%和65%；溶解磷濃度則出現(xiàn)明顯上升，庫尾朱沱斷面、變動回水區(qū)寸灘斷面上升約60%，常年回水區(qū)清溪場和萬縣斷面分別上升約20%和12%；溶解磷/沉降水樣磷的比值快速上升，2013～2019年溶解磷/沉降水樣磷均達到68%以上。④ 在新水沙條件影響下，三峽水庫干流不同形態(tài)磷濃度已發(fā)生明顯改變，最易被浮游生物利用的溶解磷已成為三峽水庫干流磷的主要存在形式。研究成果可為三峽庫區(qū)營養(yǎng)鹽控制和治理提供科學(xué)依據(jù)。

關(guān) 鍵 詞：總磷；溶解磷；顆粒磷；含沙量；Mann-Kendall非參數(shù)檢驗；三峽水庫

中圖法分類號：TV697;X524

文獻標志碼：ADOI：10.16232/j.cnki.1001-4179.2024.10.010

0 引 言

磷是生物有機體生長和能量傳輸所必需的營養(yǎng)元素，也是自然水域中重要的營養(yǎng)或污染物質(zhì)^［1-3］。天然水體中磷的主要存在形態(tài)包括溶解態(tài)磷（dissolved phosphorus，DP）和顆粒態(tài)磷（particulate phosphorus，PP），且兩種形態(tài)之間會隨著地球化學(xué)過程的變遷而循環(huán)轉(zhuǎn)化^［4-6］。泥沙和水流是水體中磷遷移轉(zhuǎn)化的主要載體。河流水體磷輸移過程中，DP主要沿水流方向擴散遷移，而PP主要伴隨懸浮泥沙顆粒進行輸移，即部分PP沉降至水底、部分PP輸移更遠的距離至下游河道^［7-8］。

三峽水庫蓄水后，庫區(qū)水位抬升，干流水體流速減緩，輸沙能力降低。尤其自2013年開始，金沙江中下游梯級電站相繼運行并產(chǎn)生較強的攔沙作用，三峽水庫入庫泥沙大幅減少^［9-10］，梯級電站攔沙效應(yīng)對三峽水庫干流不同形態(tài)磷濃度的影響備受關(guān)注。丁帥等^［11］分析2009～2010年實測數(shù)據(jù)發(fā)現(xiàn)，三峽水庫沉積的總磷（total phosphorus，TP）、PP和DP分別占入庫的24.2%，30.8%和22.2%。秦延文等^［12］分析2015年1月和7月的實測數(shù)據(jù)發(fā)現(xiàn)，在豐水期和枯水期三峽水庫水體中磷均以DP為主要賦存形態(tài)，DP通量占TP通量的51%～96%，清溪場—萬州段是水、沙和磷的主要滯留區(qū)域。這些研究多為某一時段或某幾次監(jiān)測結(jié)果，覆蓋的時間范圍窄，無法完整展現(xiàn)梯級電站攔沙效應(yīng)對三峽水庫干流不同形態(tài)磷濃度的影響。也有學(xué)者利用長系列監(jiān)測數(shù)據(jù)開展相關(guān)研究，例如婁保鋒等^［13］研究了包括三峽水庫在內(nèi)的長江干流2001～2018年TP等水質(zhì)參數(shù)的變化趨勢。然而這些研究在監(jiān)測TP時往往按照GB 3838-2002《地表水環(huán)境質(zhì)量標準》要求對水樣靜置沉降30 min后再進行測定，該監(jiān)測方法無法完整地反映水體中的磷含量。

為彌補上述研究存在的不足，一些學(xué)者利用更有針對性的監(jiān)測數(shù)據(jù)開展研究。王殿常等^［14］發(fā)現(xiàn)，2003～2012年間，三峽水庫TP的攔截率為59.8%，DP的攔截率為14.6%，而2013～2017年間庫區(qū)TP和DP攔截率分別減小為14.4%和18.5%，但是該研究中大部分討論均圍繞寸灘（三峽水庫干流）、武?。踅⒛辖蜿P(guān)（三峽水庫壩下）3個斷面展開，研究范圍未集中在三峽水庫干流。周建軍等^［4］指出上游大量梯級水庫徹底改變了長江中下游泥沙狀況，攔截超過了75%總磷和大約6萬～10萬t潛在生物有效磷，但是該研究收集的不同形態(tài)磷監(jiān)測數(shù)據(jù)為2002年9月至2003年5月和2004年5～9月，覆蓋時間較早且序列較短。Zhang等^［15］收集了2004～2017年三峽庫區(qū)不同形態(tài)磷濃度數(shù)據(jù)，分析發(fā)現(xiàn)三峽庫區(qū)干流水體中TP和PP濃度與含沙量呈正相關(guān)，且該時期三峽庫區(qū)TP和PP濃度急劇下降，而DP濃度上升。值得注意的是，上述研究均以渾濁水樣磷濃度代表TP，對渾濁水樣TP和按照GB 3838-2002《地表水環(huán)境質(zhì)量標準》要求的沉降水樣TP間的差別，以及水體含沙量與該區(qū)別間的關(guān)系等方面并未涉及。

本文系統(tǒng)檢測了2003～2019年三峽水庫干流典型斷面渾濁水樣TP、沉降水樣TP、DP、PP等不同形態(tài)磷含量，利用M-K突變檢驗、季節(jié)性Kendall檢驗和Pearson相關(guān)系數(shù)等方法系統(tǒng)分析了不同形態(tài)磷濃度的變化趨勢及其與水沙條件的相關(guān)性，以期為全面了解三峽水庫干流水體磷遷移轉(zhuǎn)化提供基礎(chǔ)數(shù)據(jù)支撐。

1 研究區(qū)域與方法

1.1 研究區(qū)域

金沙江是三峽水庫入庫泥沙最主要的來源，其下游向家壩水電站從2012年10月初期蓄水，溪洛渡水電站從2013年5月初期蓄水。水電樞紐蓄水后，長江上游來水來沙均發(fā)生較大變化^［16］。本次研究的空間范圍為三峽水庫干流段，選取朱沱（干流入庫控制站）、寸灘（變動回水區(qū)，嘉陵江匯入口下游）、清溪場（常年回水區(qū)，烏江匯入口下游）、萬縣（常年回水區(qū)）4個典型斷面進行調(diào)查監(jiān)測。在監(jiān)測過程中，按照SL 219-2018《水環(huán)境監(jiān)測規(guī)范》要求，在每個斷面設(shè)置2～3條垂線、每條垂線設(shè)置2～3個采樣點，采用船載深水采樣器分層采集水樣。監(jiān)測斷面分布見圖1。

1.2 數(shù)據(jù)來源

水質(zhì)數(shù)據(jù)來源于長江水利委員會水文局每月1次的水質(zhì)監(jiān)測結(jié)果，含沙量數(shù)據(jù)來源于水文站同步監(jiān)測結(jié)果。為更好地說明三峽水庫干流水體中不同形態(tài)磷濃度變化及可沉降固體對其分布的影響，分別對渾濁水樣、沉降30 min水樣（根據(jù)GB 3838-2002《地表水環(huán)境質(zhì)量標準》，在開展總磷等參數(shù)監(jiān)測、評價時，水樣采集后需靜置沉降30 min，取上層非沉降部分測定^［17］）和0.45 μm濾膜過濾后水樣進行測定。測得的渾濁水樣磷濃度記為TCP（total cloudy phosphorus），沉降30 min水樣磷濃度記為TP（total phosphorus），0.45 μm濾膜過濾后水樣磷濃度記為TDP（total dissolved phosphorus）。TP與TDP之差即為顆粒態(tài)磷濃度^［12］，記為TPP（total particulate phosphorus）；TCP與TDP之差記為TNDP（total non-dissolvable phosphorus）。

研究時段為2003～2019年，每月同步測定水體中TCP、TP和TDP，其中TDP從2003年6月開始監(jiān)測。

1.3 研究方法

1.3.1 Mann-Kendall（M-K）非參數(shù)檢驗

該方法適用于徑流及氣溫、降水等氣象要素時間序列的趨勢變化分析和突變分析，其優(yōu)越性在于能夠檢驗線性或非線性的趨勢，計算簡便、不受少數(shù)缺失數(shù)值和異常值干擾、對樣本分布沒有要求^［18-19］。本次研究利用該方法分析水質(zhì)數(shù)據(jù)在時間序列上的變化趨勢及突變特征。

當M-K檢驗用于檢驗序列時，構(gòu)造一秩序列：

S=∑ki=1∑i-1jak=2，3，…，n（1）

其中：

a=1，x>x0，x≤x（2）

在其原序列隨機獨立等假設(shè)下，S的方差和期望分別為

E（S）=k（k-1）/4（3）

Var（S）=k（k-1）（2k+5）/72（4）

將S標準化得到：

UF=S-E（S）Var（S） k=1，2，…，n（5）

UF為統(tǒng)計量序列，構(gòu)成一條UF曲線，通過信度檢驗可得出所檢驗指標是否有明顯的變化趨勢。將時間序列x按逆序排列，通過此方法計算得到另一條曲線UB：

UB=-UF′k′=n+1-k（6）

兩條曲線在置信區(qū)間內(nèi)的交點確定為突變點。取顯著性水平α=5%時，其對應(yīng)的臨界值為1.96。UF>0，表示序列呈上升趨勢；UF<0，表明呈下降趨勢；UF>1.96或UF<-1.96，表示上升或下降趨勢明顯。

1.3.2 季節(jié)性Kendall（S-K）非參數(shù)檢驗

S-K檢驗是在M-K檢驗基礎(chǔ)上考慮了季節(jié)性變化，其主要優(yōu)點是不受可用資料系列存在漏測值、未檢出值的影響。該檢驗方法的原理是將歷年相同月（季）的水質(zhì)資料進行比較，如果后面的值（時間上）高于前面的值記為“+”號，否則記作“-”號。如果“+”號的個數(shù)比“-”號的多，則為上升趨勢；反之，則為下降趨勢，如果相等則為無趨勢。通常取顯著性水平α為0.10和0.01。α≤0.01指示檢驗具有高度顯著水平，0.01<α≤0.10指示檢驗具有顯著水平。

1.3.3 Pearson相關(guān)系數(shù)

首先利用Origin軟件計算各監(jiān)測斷面含沙量與不同形態(tài)磷濃度之間的Pearson相關(guān)系數(shù)，再依據(jù)相關(guān)系數(shù)的大小分析它們之間的相關(guān)關(guān)系^［20］。

2 結(jié)果與分析

2.1 三峽水庫干流水沙總體變化

三峽水庫范圍主要有嘉陵江、烏江、小江、大寧河、香溪河等幾十條河流匯入，其中以嘉陵江和烏江年徑流量較大，其余相對較小^［21］。嘉陵江入?yún)R口位于寸灘斷面上游，烏江入?yún)R口位于清溪場斷面上游。

三峽水庫干流各監(jiān)測斷面的水沙變化見圖2。三峽水庫蓄水以來，受支流入?yún)R和降雨面積增大的影響，庫區(qū)干流各監(jiān)測斷面沿程流量增大，位于常年回水區(qū)的清溪場和萬縣兩斷面平均流量差異較小。從年際變化來看，2003～2019年各斷面平均流量無明顯變化，而含沙量則明顯減少。各監(jiān)測斷面含沙量的減少主要發(fā)生在2013年以前，自2013～2019年，三峽水庫干流水體含沙量已降至較低水平，這一階段含沙量相對比較穩(wěn)定。2013～2019年與2003～2012年相比，朱沱、寸灘、清溪場、萬縣斷面的平均含沙量分別減小了76.1%，65.7%，63.6%，66.6%。

從年內(nèi)變化來看（圖3），各監(jiān)測斷面的含沙量時間分布基本一致，汛期（6～9月）要明顯大于蓄水期（10～12月）和消落期（1～5月），在消落期和蓄水期，各監(jiān)測斷面的含沙量基本都小于0.3 kg/m³。從空間分布來看，2003～2012年，各監(jiān)測斷面含沙量表現(xiàn)為沿程逐漸減小，即朱沱>寸灘>清溪場>萬縣；2013～2019年，各個水期的含沙量均明顯減小，在達到峰值的7月各監(jiān)測斷面的含沙量均小于0.5 kg/m³，僅為2003～2012年的30%～53%，在沿程分布上，萬縣斷面含沙量仍然最小，朱沱、寸灘、清溪場3個斷面的含沙量差異則明顯減小。

2.2 三峽水庫干流水體不同形態(tài)磷濃度總體變化

三峽水庫干流各監(jiān)測斷面2003～2019年不同形態(tài)磷濃度年均值變化趨勢見圖4。如圖4（a）、（b）所示，TCP和TP年際變化趨勢基本一致。2003～2012年，TCP和TP濃度起伏波動、整體呈緩慢上升趨勢，兩次峰值出現(xiàn)在2008年和2012年，各監(jiān)測斷面TCP濃度要明顯高于TP。該時間段內(nèi)，TCP濃度均值為0.203～0.285 mg/L，TP濃度均值為0.139～0.183 mg/L，TP濃度僅為TCP濃度的65%左右。自2013年起各監(jiān)測斷面TCP和TP濃度均開始下降，2014年開始各監(jiān)測斷面TP濃度均達到河流Ⅲ類水質(zhì)要求，TCP和TP濃度值逐漸接近，TP/TCP約為80%，萬縣斷面則達到90%。2013～2019年TCP濃度均值為0.121～0.166 mg/L，TP濃度均值為0.108～0.133 mg/L，與2003～2012年多年均值相比，TCP下降了約40%，TP下降了約25%。從沿程變化來看，除清溪場外，各監(jiān)測斷面TCP和TP濃度基本呈從上游到下游沿程降低趨勢。清溪場斷面TCP和TP濃度較高，可能是因為受上游入?yún)R支流烏江的影響，2009～2013年烏江TP顯著高于庫內(nèi)相關(guān)斷面，烏江匯入在這期間可使清溪場斷面TP年均值最大升高約0.057 mg/L^［22-23］。

TDP濃度年際變化趨勢與TP不同。如圖4（c）所示，2002～2007年，各監(jiān)測斷面TDP均呈窄幅波動，其中萬縣斷面2004年TDP較上年增大一倍。2008～2013年，各監(jiān)測斷面的TDP快速上升，2013年以后出現(xiàn)下降趨勢，但總體仍略高于2009年以前的濃度水平。與歷年最大值相比，從上游至下游各監(jiān)測斷面2019年TDP濃度分別下降62.4%，59.1%，65.2%，56.9%。2013～2019年TDP濃度均值較2003～2012年均值則有所上升，沿程分別上升56.8%，59.5%，20.5%，11.6%。

TPP濃度歷年最大值均出現(xiàn)在2008年（圖4（d）），說明在175.00 m試驗性蓄水以后，三峽水庫干流水體的TPP濃度已經(jīng)開始出現(xiàn)下降趨勢。自2013年以后，各監(jiān)測斷面的TPP急劇減小，2013～2019年TPP濃度均值降低至0.020～0.040 mg/L，較2003～2012年濃度均值下降了約65%。TNDP變化趨勢與TCP接近，如圖4（e）所示，2012年后各監(jiān)測斷面TNDP濃度開始顯著下降。相比而言，各監(jiān)測斷面TNDP濃度要明顯高于TPP，其中2003～2012年TNDP是TPP的1.9～2.2倍，2013～2019年TNDP則是TPP的1.6～1.8倍。從沿程變化來看，朱沱、寸灘、清溪場3個斷面的TPP和TNDP濃度變化不大，萬州斷面的TPP和TNDP濃度均顯著低于上游的各監(jiān)測斷面。

如圖5（a）～（d）所示，2003～2012年，各監(jiān)測斷面TCP和TP月均值呈倒“V”形，汛期（6～9月）>蓄水期（10～12月）>消落期（1～5月），其中各監(jiān)測斷面TCP汛期平均濃度約為蓄水期的2～3倍、消落期的2.7～3.3倍，TP汛期平均濃度約為蓄水期的1.8～2.0倍、消落期的1.7～2.4倍。TCP和TP濃度之差在汛期較為明顯，其中各監(jiān)測斷面TCP濃度約為TP濃度的1.7～1.8倍，而兩者間差異在1～5月和11～12月相對較小。2013～2019年，各監(jiān)測斷面TCP和TP月均值變化幅度降低，其中各監(jiān)測斷面TCP汛期平均濃度約為蓄水期的1.5～1.9倍、消落期的1.0～1.5倍，TP各月濃度近似呈“一”字形，不同水期間變化幅度較小。汛期懸浮泥沙含量高，單位體積水中懸浮泥沙所提供的對磷的有效吸附位數(shù)量越多^［24］，這可能是TCP汛期濃度顯著高于蓄水期和消落期的主要原因。

如圖5（e）～（f）所示，各監(jiān)測斷面TDP月均值變化趨勢與TCP、TP相反，基本呈汛期較低、蓄水期和消落期較高的趨勢。其中，朱沱、寸灘、清溪場3個斷面的TDP濃度均為消落期>蓄水期>汛期，萬縣斷面TDP則為消落期>汛期>蓄水期。2003～2012年各監(jiān)測斷面TDP多年月均值均顯著低于2013～2019年。

如圖5（g）～（j）所示，各監(jiān)測斷面TPP和TNDP濃度2003～2012年表現(xiàn)為汛期>蓄水期>消落期，2013～2019年蓄水期和消落期濃度差異逐漸縮小，表現(xiàn)為汛期>蓄水期≈消落期。

從沿程變化上看，除清溪場斷面外，各監(jiān)測斷面TCP、TP、TPP和TNDP均表現(xiàn)為從上游至下游沿程減小的趨勢，即朱沱＞寸灘＞萬縣，TDP則表現(xiàn)為從上游至下游沿程增加的趨勢，即朱沱≈寸灘<清溪場≈萬縣。

2.3 不同形態(tài)磷濃度和含沙量的趨勢分析

分別對各監(jiān)測斷面2003～2019年的不同形態(tài)磷濃度逐月監(jiān)測值以及含沙量變化情況進行M-K突變檢驗及趨勢分析，結(jié)果見表1。

各監(jiān)測斷面TCP和TP突變點均出現(xiàn)在2017年附近，表明TCP和TP濃度在2017年之后呈下降趨勢。這與圖2（a）、（b）的結(jié)果一致，即2017年開始，各監(jiān)測斷面的TP濃度基本能夠達到河流Ⅱ類水質(zhì)類別。張漫等^［18］研究發(fā)現(xiàn)2016年是三峽水庫總磷濃度變化的突變節(jié)點，在2016年之后有顯著下降趨勢，與本次研究基本一致。

各監(jiān)測斷面TDP突變點出現(xiàn)于2004～2008年，其中朱沱和寸灘斷面突變點均在2008年，清溪場斷面突變點在2007年，萬縣斷面突變點在2004年底，TDP上升趨勢出現(xiàn)節(jié)點表現(xiàn)為從上游到下游沿程提前。

各監(jiān)測斷面TPP和TNDP突變點出現(xiàn)于2011～2014年，其中朱沱斷面TPP和TNDP突變點均在2014年，寸灘斷面TPP和TNDP突變點均在2013年，清溪場斷面TPP突變點在2012年底、TNDP突變點在2013年，萬縣斷面TPP突變點在2011年、TNDP突變點在2012年初。TPP和TNDP下降趨勢出現(xiàn)節(jié)點表現(xiàn)為從上游到下游沿程提前。

朱沱和寸灘斷面含沙量突變點在2012年底，清溪場斷面突變點在2010年，萬縣斷面突變點則在2004年。這表明含沙量變化趨勢與TPP、TNDP類似，含沙量下降趨勢出現(xiàn)節(jié)點表現(xiàn)為從上游到下游沿程提前。尤其是2013年之后，各監(jiān)測斷面含沙量均快速下降。

采用季節(jié)性Kendall檢驗方法進行總磷變化趨勢分析，并采用水質(zhì)變化趨勢分析軟件PWQ Trend進行計算。總體來看，三峽水庫建成運行以來，各監(jiān)測斷面TCP、TP、TPP、TNDP和含沙量均呈顯著下降趨勢，而TDP呈顯著上升趨勢。其中，TCP、TP、TDP在2003～2012年呈上升趨勢，2013～2019年呈下降趨勢，TPP、TNDP和含沙量始終保持下降趨勢。這可能是因為2003～2010年三峽水庫試驗性蓄水期間，庫區(qū)新增淹沒區(qū)導(dǎo)致磷濃度有所上升，而步入高水位正常運行期后，面源新增污染趨于穩(wěn)定，隨著三峽大壩和金沙江上游梯級水庫攔沙效應(yīng)顯現(xiàn)，含沙量減小導(dǎo)致顆粒磷濃度持續(xù)降低，三峽水庫不同形態(tài)的磷濃度均呈下降趨勢。

2.4 不同形態(tài)磷濃度和含沙量的相關(guān)性分析

為進一步分析含沙量對不同形態(tài)磷濃度的影響，計算了2003～2019年含沙量與不同形態(tài)磷濃度的Pearson相關(guān)系數(shù)。從表2可以看出，2003～2012年，各監(jiān)測斷面TCP、TP、TPP和TNDP與含沙量的Pearson相關(guān)系數(shù)均大于0.6，具體表現(xiàn)為TNDP>TCP>TPP>TP，而TDP與含沙量相關(guān)性較差。2013～2019年，各監(jiān)測斷面不同形態(tài)的磷與含沙量相關(guān)性為TPP>TNDP>TCP>TP>TDP。劉尚武等^［25］根據(jù)2003～2016年實測水文資料研究發(fā)現(xiàn)三峽水庫入庫細顆粒泥沙的排沙比大于粗顆粒泥沙，呈含沙量明顯減少且泥沙細化的特點。TPP與含沙量的相關(guān)系數(shù)大于TNDP和TCP與含沙量的相關(guān)系數(shù)，可能是因為細顆粒泥沙與TPP的相關(guān)性更強，而粗顆粒的泥沙則與TNDP和TCP的相關(guān)性更強。與2003～2012年間相比，各監(jiān)測斷面不同形態(tài)的磷與含沙量相關(guān)性均不同程度減弱，僅清溪場斷面TPP和TNDP與含沙量的Pearson相關(guān)系數(shù)小幅上升。從監(jiān)測斷面來看，2013年以后，朱沱、寸灘、清溪場斷面的TPP、TNDP、TCP與含沙量的相關(guān)性呈現(xiàn)沿程增強的趨勢，而萬縣斷面相關(guān)系數(shù)則明顯減小。

3 討 論

3.1 含沙量對TDP/TP比例的影響

從各監(jiān)測斷面TDP與TP年均值之比來看（圖6），2003～2008年各監(jiān)測斷面TDP/TP均緩慢降低，2008年以后各監(jiān)測斷面TDP/TP快速升高，這與TDP濃度的M-K突變檢驗結(jié)果一致。其中，各監(jiān)測斷面TDP/TP在2012年均稍有回落，自2013年后基本維持在70%以上，2018～2019年為60%左右。2013～2019年，各監(jiān)測斷面的TDP/TP均值從上游到下游依次為68.5%，68.1%，69.1%，81.6%，這表明金沙江中下游梯級電站相繼運行后，三峽水庫萬縣斷面TDP/TP的比值較高，達到80%以上，朱沱、寸灘、清溪場3個斷面的TDP/TP則相差不大，均達到68%。

各監(jiān)測斷面TDP與TP月均值之比與含沙量月均值成反比。如圖7所示，當含沙量低于0.3 kg/m³時，TDP/TP比值隨含沙量的增加而快速下降，然后逐漸變緩；當含沙量超過1.0 kg/m³，TDP/TP比值趨于穩(wěn)定。這與報道的運用三峽水庫干流2004～2017年監(jiān)測數(shù)據(jù)的研究結(jié)果一致^［15］。當含沙量超過1.0 kg/m³，朱沱、寸灘兩斷面TDP/TP的占比穩(wěn)定在25.0%左右，清溪場斷面穩(wěn)定在30.0%左右，萬縣斷面穩(wěn)定在40.0%左右，各監(jiān)測斷面從上游到下游TDP/TP比值呈增大趨勢。

3.2 含沙量對顆粒態(tài)磷濃度的影響

泥沙是水中的重要磷庫，發(fā)揮著污染緩沖和營養(yǎng)調(diào)節(jié)等重要作用。泥沙顆粒的吸附作用可有效降低高磷污染期間水體溶解磷濃度，又可通過解吸補充低磷期間浮游生物對磷的需求^［26-27］。三峽水庫建成運行以來，朱沱—萬縣段的平均含沙量減小了63.6%～76.1%，TPP濃度下降了約63.6%～69.7%，TNDP濃度下降了63.2%～74.7%，TPP和TNDP與含沙量的下降幅度基本一致。泥沙減少導(dǎo)致三峽水庫非汛期水體透光深度比長江自然水體透光深度增加了8～10倍^［26］，同時水體中最易被浮游生物進行利用的溶解磷已成為磷的主要存在形式，兩者疊加之下，加大了水庫藻類暴發(fā)的風(fēng)險。

顆粒態(tài)磷濃度降低，也進一步減小了TCP與TP的濃度差。2003～2012年三峽水庫各監(jiān)測斷面TP/TCP約為65%左右，2013～2019年TP/TCP約為80%，其中萬縣斷面達到90%。各監(jiān)測斷面TP/TCP和含沙量的關(guān)系曲線如圖8所示。當含沙量小于0.15 kg/m³時，TP/TCP可以達到70%及以上。結(jié)合2.1節(jié)各監(jiān)測斷面水沙條件可以得出，對于三峽水庫干流河段，在蓄水期和消落期，按照現(xiàn)行的總磷監(jiān)測標準方法要求，水樣靜置30 min后測得的TP可以近似代表TCP；在汛期，若采用TP來代表TCP，則會有相當部分的磷被漏測，顆粒態(tài)磷被低估，導(dǎo)致泥沙對磷的影響機理和實質(zhì)不能全面反映出來。

4 結(jié) 論

（1）2003～2019年三峽水庫干流朱沱—萬縣段的平均流量無明顯變化，含沙量則明顯減少，含沙量的減少主要發(fā)生在2013年以前，2013～2019年含沙量則無明顯變化。

（2）三峽水庫建成運行以來，三峽水庫干流各監(jiān)測斷面TP濃度呈降低趨勢，其中，TPP濃度明顯降低，TDP濃度增大。2017年開始，各監(jiān)測斷面的TP濃度基本能夠達到河流Ⅱ類水質(zhì)類別。

（3）不同形態(tài)的磷濃度與含沙量的Pearson相關(guān)系數(shù)在2003～2012年為TNDP>TCP>TPP>TP，2013～2019年為TPP>TNDP>TCP>TP。三峽水庫細顆粒泥沙與TPP的相關(guān)性更強，而粗顆粒的泥沙則與TNDP和TCP的相關(guān)性更強。

（4）三峽水庫干流TDP/TP的比值與含沙量成反比。TDP/TP在2003～2012年為30%～56%，2013～2019年為68%～81%，其中，近壩斷面萬縣TDP/TP的比值最高，達到81%，朱沱、寸灘、清溪場3個斷面的TDP/TP則相差不大，均達到68%以上。

（5）2003～2012年三峽水庫干流TP/TCP約為65%左右，2013～2019年約為80%，其中萬縣斷面達到90%。在蓄水期和消落期，可以用TP近似代表TCP；在汛期，則必須測定TCP才能全面反映水體中的磷濃度。

參考文獻：

［1］ VITOUSEK P M，PORDER S，HOULTON B Z，et al.Terrestrial phosphorus limitation：mechanisms，implications，and nitrogen-phosphorus interactions［J］.Ecological Applications，2010，20（1）：5-15.

［2］ 胡佳，李艷華.天然水體中磷的存在形態(tài)及其對魚類影響研究進展［J］.水產(chǎn)學(xué)雜志，2020，33（4）：81-88.

［3］ 吳浩云，賈更華，徐彬等.1980年以來太湖總磷變化特征及其驅(qū)動因子分析［J］.湖泊科學(xué)，2021，33（4）：974-991.

［4］ 周建軍，張曼，李哲.長江上游水庫改變干流磷通量、效應(yīng)與修復(fù)對策［J］.湖泊科學(xué)，2018，30（4）：865-880.

［5］ 朱廣偉，秦伯強，張運林，等.近70年來太湖水體磷濃度變化特征及未來控制策略［J］.湖泊科學(xué)，2021，33（4）：957-973.

［6］ WORSFOLD P，MCKELVIE I，MONBET P.Determination of phosphorus in natural waters：a historical review［J］.Analytical Chimica Acta，2016，918：8-20.

［7］ YAN Q Z，DU J T，CHEN H T，et al.Particle-size distribution and phosphorus forms as a function of hydrological forcing in the Yellow River［J］.Environmental Science and Pollution Research，2016，23（4）：3385-3398.

［8］ 黃磊，方紅衛(wèi)，王靖宇，等.河流泥沙磷遷移過程的數(shù)學(xué)模型研究［J］.水利學(xué)報，2014，45（4）：394-402.

［9］ 王曉青，呂平毓，胡長霜.三峽庫區(qū)懸移質(zhì)泥沙對TP、TN等的吸附影響［J］.人民長江，2006，37（7）：15-17.

［10］李思璇，宋瑞，許全喜，等.長江上游總磷通量時空變化特征研究［J］.環(huán)境科學(xué)與技術(shù)，2021，44（5）：179-185.

［11］DING S，CHEN P P，LIU S M，et al.Nutrient dynamics in the Changjiang and retention effect in the Three Gorges Reservoir［J］.Journal of HydRS3BMjUbFuRXG6GNcTdUWw==rology，2019，574：96-109.

［12］秦延文，韓超南，鄭丙輝，等.三峽水庫水體溶解磷與顆粒磷的輸移轉(zhuǎn)化特征分析［J］.環(huán)境科學(xué)，2019，40（5）：2152-2159.

［13］婁保鋒，卓海華，周正，等.近18年長江干流水質(zhì)和污染物通量變化趨勢分析［J］.環(huán)境科學(xué)研究，2020，33（5）：1150-1162.

［14］王殿常，吳興華.三峽水庫磷輸移規(guī)律研究［J］.水利學(xué)報，2021，52（8）：885-895.

［15］ZHANG S Y，ZENG Y H，ZHA W，et al.Spatiotemporal variation of phosphorus in the Three Gorges Reservoir：impact of upstream cascade reservoirs［J］.Environmental Science and Pollution Research，2022，29：56739-56749.

［16］劉尚武，張小峰，許全喜，等.近50年來金沙江流域懸移質(zhì)輸沙特性研究［J］.泥沙研究，2020，45（3）：30-37.

［17］蘭靜，吳云麗，婁保峰，等.2004年以來長江中下游干流水體高錳酸鹽指數(shù)時空變化分析［J］.湖泊科學(xué)，2021，33（4）：1112-1122.

［18］張漫，張萬順，張瀟，等.三峽水庫175 m試驗性蓄水后庫區(qū)水質(zhì)時空變化規(guī)律［J］.人民長江，2022，53（3）：68-73，91.

［19］謝經(jīng)朝，趙秀蘭，何丙輝，等.漢豐湖流域農(nóng)業(yè)面源污染氮磷排放特征分析［J］.環(huán)境科學(xué)，2019，40（4）：1760-1769.

［20］李瓊芳，許樹洪，陳啟慧，等.環(huán)太湖各水資源分區(qū)入出湖河流總磷濃度與負荷變化分析［J］.湖泊科學(xué)，2022，34（1）：74-89.

［21］中國河湖大典編纂委員會.中國河湖大典［M］.北京：中國水利水電出版社，2010.

［22］翟婉盈，湛若云，卓海華，等.三峽水庫蓄水不同階段總磷的變化特征［J］.中國環(huán)境科學(xué)，2019，39（12）：5069-5078.

［23］卓海華，蘭秀薇，婁保鋒，等.烏江匯入對三峽水庫局部水體中磷質(zhì)量濃度的影響［J］.環(huán)境工程學(xué)報，2023，17（11）：3497-3504.

［24］婁保鋒，邱光勝.嘉陵江總磷通量變化及空間來源解析（2002～2019）［J］.中國環(huán)境科學(xué)，2024，44（3）：1448-1456.

［25］劉尚武，張小峰，許全喜，等.三峽水庫區(qū)間來沙量估算及水庫排沙效果分析［J］.湖泊科學(xué)，2019，31（1）：28-38.

［26］尹煒，王超，張洪.長江流域總磷問題思考［J］.人民長江，2022，53（4）：44-52.

［27］陳峰，李曉慧，王文靜，等.2018～2022年長江流域總磷時空變化特征研究［J］.水利水電快報，2023，44（8）：91-97，105.

（編輯：劉 媛）

Influence of sediment concentration on different forms of phosphorus concentration

in main stream of Three Gorges ReservoirWANG Wenjing，ZHAO Min，ZHAO Qichao

（Bureau of Hydrology，Changjiang Water Resources Commission，Wuhan 430010，China）

Abstract： Phosphorus，no matter dissolved and particulate forms，is an important index for water quality assessment and pollution control supervision.It is of great significance to study influence of sediment concentration variation on the concentration of different phosphorus forms.In this paper，the main stream of Three Gorges Reservoir （TGR） was selected as the research object，and the Mann-Kendall nonparametric test method and Pearson correlation coefficient were used to analyze the influence of sediment concentration on the concentration of different forms of phosphorus，such as total phosphorus，dissolved phosphorus and particulate phosphorus.The results showed that：① From 2003 to 2019，the average flow of monitoring sections in the main stream of the TGR hasn’t changed significantly，however the sediment concentration decYKMKNI9tm2OZTgu0uEKoNw==reased significantly.② The total phosphorus concentration in the mainstream of the TGR gently fluctuating upward from 2003 to 2012，and continued to decline from 2013 to 2019；the concentration of dissolved phosphorus increased first and then decreased，and the peak value of each section appeared in 2013 and 2014.The concentration of particulate phosphorus decreased in a fluctuated way.③ Compared with the average concentration from 2003 to 2012，the total phosphorus concentration in cloudy waters decreased by about 40%，the total phosphorus concentration in clarified waters decreased by about 25%，the particulate phosphorus decreased by about 65%.The dissolved phosphorus at Zhutuo Section at the end of the reservoir and Cuntan Section in the varied backwater area increased by about 60%，and the dissolved phosphorus in Qingxichang and Wanxian sections in the perennial backwater area increased by about 20% and 12% respectively.The ratio of dissolved phosphorus to clarified total phosphorus also increased rapidly，with dissolved phosphorus/clarified total phosphorus reaching over 68% from 2013 to 2019.④ Under the influence of new water and sediment regimes，phosphorus concentrations of different forms in the mainstream of the TGR have changed significantly.Dissolved phosphorus which is most easily used by plankton，has become the main form of phosphorus in the mainstream of the TGR.This study can provide scientific basis for nutrient control and management in the TGR area.

Key words： total phosphorus；dissolved phosphorus；particulate phosphorus；sediment content；Mann-Kendall nonparametric tests；Three Gorges Reservoir