(山西省水文水資源勘測(cè)分局，山西 太原 030000)

底泥是湖庫(kù)營(yíng)養(yǎng)物質(zhì)的重要蓄積庫(kù),污水的排入、地表徑流匯集以及水生生物殘骸,導(dǎo)致底泥中N、P元素逐步積累,易形成營(yíng)養(yǎng)物質(zhì)的內(nèi)負(fù)荷。當(dāng)外源排入湖庫(kù)的N、P營(yíng)養(yǎng)元素負(fù)荷量減少后,沉積物中的營(yíng)養(yǎng)元素會(huì)逐步釋放,成為湖泊富營(yíng)養(yǎng)化的主導(dǎo)因子[1-2]。磷被認(rèn)為是多數(shù)水體浮游藻類限制性營(yíng)養(yǎng)元素，湖泊沉積物中磷的特征、水土界面間的行為以及沉積物釋放規(guī)律研究，對(duì)探求湖泊富營(yíng)養(yǎng)化防治意義深遠(yuǎn)[3]。

漳澤水庫(kù)，又稱太行湖(屬于半封閉湖)，是人工湖泊，興建于1959年，位于山西省東南部的長(zhǎng)治市郊區(qū)，全湖區(qū)水面面積16.48km2(低水位時(shí))，最深處位于庫(kù)區(qū)泄水洞，約9m，庫(kù)尾深1m，庫(kù)中深4m，最高溫度出現(xiàn)在8月，約29～32℃，結(jié)冰期為12月上旬至次年3月，低溫為2～3℃，多年水體超過(guò)Ⅳ類水質(zhì)，處于富營(yíng)養(yǎng)化狀態(tài)。2018年初次對(duì)湖區(qū)下游至湖尾均勻布設(shè)監(jiān)測(cè)點(diǎn)1～16號(hào)(見(jiàn)圖1)，監(jiān)測(cè)底泥中總磷等污染物，并進(jìn)行了室內(nèi)靜態(tài)化學(xué)條件下總磷等營(yíng)養(yǎng)物質(zhì)釋放模擬試驗(yàn)，旨在湖區(qū)外源污染截?cái)嗪?，探究底泥?duì)水體的二次污染情況。由于室外模擬的難度較大，所以采集底泥進(jìn)行室內(nèi)分析及靜態(tài)模擬試驗(yàn)，通過(guò)近一年的底泥監(jiān)測(cè)和室內(nèi)靜態(tài)模擬試驗(yàn)，基本掌握了湖區(qū)底泥總磷的分布情況及總磷釋放初步特征，為解決水庫(kù)的富營(yíng)養(yǎng)化治理提供了參考。

圖1 太行湖(漳澤水庫(kù))底泥總磷分布 (單位：mg/L)

1 湖區(qū)底泥中總磷(TP)分布特征

1.1 時(shí)間變化特征

全湖監(jiān)測(cè)點(diǎn)總磷含量隨時(shí)間變化特征見(jiàn)圖2，最大值均出現(xiàn)在8月，其他月份變化不大，最小值均出現(xiàn)在10月，4月3～6號(hào)監(jiān)測(cè)點(diǎn)出現(xiàn)最小值，與10月含量基本一致?？傮w變化規(guī)律為8月驟增，9月開(kāi)始回落，到10月達(dá)到最低值，其后稍有回升。2號(hào)監(jiān)測(cè)點(diǎn)9月回落緩慢，16號(hào)監(jiān)測(cè)點(diǎn)10月回落值與7月起漲點(diǎn)達(dá)到水平。上述分布特征說(shuō)明，污染物的主要來(lái)源是上游河流洪水期(8月)攜帶的大量泥沙及陸源污染物進(jìn)入湖區(qū)，由于洪水期水動(dòng)力條件強(qiáng)、水流快，造成整個(gè)湖區(qū)底泥TP含量增高，2號(hào)區(qū)域9月比其他區(qū)域高也說(shuō)明大量泥沙等污染沉積后，一些細(xì)小的顆粒物繼續(xù)向下游低凹區(qū)域運(yùn)移沉積，16號(hào)監(jiān)測(cè)點(diǎn)位于主河槽入湖口，9月后仍會(huì)有部分泥沙輸入，但此時(shí)水動(dòng)力條件減弱，大部分?jǐn)y帶總磷的泥沙等停滯在此，造成TP比其他區(qū)域含量偏高。

圖2 底泥中總磷(TP)隨時(shí)間變化特征

1.2 空間變化特征

湖區(qū)總磷含量平均值為440mg/L,高值區(qū)位于壩址泄水洞區(qū)域和上游湖尾區(qū)，代表點(diǎn)為2號(hào)、14號(hào)、15號(hào)、16號(hào)，最大值為683mg/L,低值區(qū)在庫(kù)中心偏向庫(kù)區(qū)西側(cè)區(qū)域，代表點(diǎn)為1號(hào)、4號(hào)、6號(hào)、7號(hào)、9號(hào)、11號(hào)。底泥總磷的空間分布不均衡(見(jiàn)圖1)，湖區(qū)東側(cè)總磷平均含量整體比湖區(qū)西側(cè)高，有明顯的東西帶狀分布特征，此種特征可能是太行湖近湖尾處絳河支流水質(zhì)(Ⅱ類)好于南源(主源，劣Ⅴ類水質(zhì))水質(zhì)所致。

2 總磷釋放靜態(tài)模擬試驗(yàn)

底泥中的磷釋放與化學(xué)沉淀的形態(tài)有關(guān)，底泥中的磷主要以無(wú)機(jī)態(tài)的正磷酸鹽為主，一旦出現(xiàn)有利于鈣、鋁、鐵等不溶性磷酸鹽沉淀物溶解的條件，磷的釋放就可能發(fā)生[4]。

2.1 試驗(yàn)方法與試驗(yàn)計(jì)算公式

2.1.1 試驗(yàn)方法

將采集的湖區(qū)2號(hào)底泥原樣帶回實(shí)驗(yàn)室瀝干重力水后放入試驗(yàn)容器底部(平鋪5cm厚)，然后注入初始水樣，濃度為C0(采集湖區(qū)2號(hào)底泥上覆水，由于濃度較高，所以將其釋放至試驗(yàn)項(xiàng)目的Ⅲ類水質(zhì)以下，稱“稀釋水”)，見(jiàn)表1。每隔12h取200mL水樣檢測(cè)，并添加同體積稀釋水樣于容器內(nèi)，容器直徑10cm，水柱高25cm，采樣時(shí)于底部上方10cm處用虹吸管采取，以不擾動(dòng)底泥，加稀釋水時(shí)同樣不能擾動(dòng)底泥，如此反復(fù)，試驗(yàn)共進(jìn)行5天(120h)停止。本次共設(shè)計(jì)25組試驗(yàn)，測(cè)試不同T、pH值、DO化學(xué)條件下試驗(yàn)釋放水中的氨氮、總氮、總磷參數(shù)濃度變化，所有試驗(yàn)裝置均置于室內(nèi)自然光或日光燈照條件下。

表1 試驗(yàn)初始水樣濃度 單位：mg/L

2.1.2 試驗(yàn)計(jì)算公式

模擬試驗(yàn)考慮影響底泥釋磷主要因素：溫度、溶解氧、pH值作用。營(yíng)養(yǎng)鹽釋放試驗(yàn)不同文獻(xiàn)的公式大體一致，但是諸多文獻(xiàn)公式中各因素表述不一，下面列舉一些文獻(xiàn)中營(yíng)養(yǎng)鹽釋放試驗(yàn)公式及其概念表述，并推導(dǎo)營(yíng)養(yǎng)鹽模擬釋放試驗(yàn)公式。

2.1.3 文獻(xiàn)計(jì)算公式及概念表述列舉

a.計(jì)算沉積物的釋放強(qiáng)度γ和γ′[3]：

(1)

γ′=r/A

(2)

式中：γ為每次采水樣時(shí)的釋放量，mg；γ′為釋放強(qiáng)度，mg/m2；V為泥樣上層水體積，L，取V=1L；Cn為第n次采樣水中總磷濃度，mg/L；C0為泥樣上層水起始總磷濃度，mg/L；Ca為添水后泥樣上層水磷濃度，mg/L；n為采水樣的次數(shù)；Vn為每次采水樣體積，取Vn=35mL；A為與水接觸底泥面積，m2，A=0.0133m2。

b.沉積物的釋放強(qiáng)度[4-5]：

(3)

R=γ/(tA)

(4)

γ′=γ/A

(5)

式中:V為泥樣上方水的體積,L,取V=1L;Cn為第n次采樣時(shí)水中營(yíng)養(yǎng)鹽濃度,mg/L;C0為泥樣上方水的起始營(yíng)養(yǎng)鹽濃度,mg/L;Ca為添加水后泥樣上方營(yíng)養(yǎng)鹽濃度,mg/L;n為采水樣的次數(shù);Vn為每次采水樣的體積，取Vn=100mL;γ為每次采水樣時(shí)的釋放量,mg;γ′為釋放強(qiáng)度,mg/m2;R為釋放速度,mg/(m2·d);t為釋放時(shí)間,d;A為與水接觸的底泥面積,m2,A=0.0133m2。每次采樣后得出的R平均后即得到釋放速度。

c.底泥中污染物釋放率的計(jì)算公式[6-7]：

(6)

d.釋放速率計(jì)算[8]：

(7)

e.底泥氮磷的釋放通量采用下式計(jì)算[9]：

(8)

式中：Rn為到第n天為止底泥氮磷的釋放量,mg；V為反應(yīng)裝置中剩余水樣的體積,L；Ci為第i次采樣時(shí)氮磷的質(zhì)量濃度(i=1,2,…,n),mg/L；C0為初始氮磷的質(zhì)量濃度,mg/L；Vi為第i次采取的水樣體積,L。

f.柱狀模擬釋放營(yíng)養(yǎng)物的釋放速度按下式計(jì)算[10]：

(9)

2.1.4 營(yíng)養(yǎng)物的釋放量公式推導(dǎo)

設(shè)容器內(nèi)柱狀水體積為V(L)；各時(shí)刻釋放濃度值為C0、C1、C2、C3、C4、C5、C6、C7、C8、C9、C10、…、Cn-1、Cn(mg/L)(見(jiàn)表2)；每隔t時(shí)段取樣檢測(cè)，取樣體積Vi，同時(shí)補(bǔ)充同體積水樣至原刻度線，保證容器內(nèi)水柱體積V不變。因?yàn)槿∧骋粷舛人畼雍?，再將另一濃?設(shè)為Cai，i=1,2,3，…，n)的水樣(應(yīng)滿足低于取樣時(shí)的濃度)回加進(jìn)去，這時(shí)容器內(nèi)的濃度會(huì)發(fā)生變化，成為下一時(shí)段營(yíng)養(yǎng)鹽釋放的起始濃度。用后一時(shí)刻容器中某物質(zhì)的量(mg)減去前一時(shí)刻起始量(見(jiàn)表3)，即為某時(shí)段某營(yíng)養(yǎng)鹽的釋放量，將各時(shí)段的釋放量累加，即為試驗(yàn)期(nd)某營(yíng)養(yǎng)物質(zhì)釋放通量R(mg)，如果是氨氮可考慮容器的內(nèi)上層空間聚集的收集量。如果收集不到，計(jì)算通量可能會(huì)偏小。

表2 釋放試驗(yàn)期各時(shí)間點(diǎn)營(yíng)養(yǎng)物質(zhì)濃度

表3 某時(shí)刻取出再回填同體積(Vi)水樣后容器內(nèi)某物質(zhì)量

各時(shí)段釋放量計(jì)算如下：

0d釋放量：

0

0～1d釋放量：

VC1-VC0

1～2d釋放量：

VC2-[VC1-Vi(C1-Ca1)]

2～3d釋放量：

VC3-[VC2-Vi(C2-Ca2)]

3～4d釋放量：

VC4-[VC3-Vi(C3-Ca3)]

4～5d釋放量：

VC5-[VC4-Vi(C4-Ca4)]

……

(n-2)～(n-1)d釋放量：

VCn-1-[VCn-2-Vi(Cn-2-Can-2)]

(n-1)～nd釋放量：

VCn-[VCn-1-Vi(Cn-1-Can-1)]

將上述各時(shí)段釋放量相加得

R=VCn-[VC0+Vi(C1-Ca1)]+Vi(C2-Ca2)+Vi(C3-Ca3)+…+Vi(Cn-2-Can-2)+Vi(Cn-1-Can-1)

即：

i=1,2,3,…,n

(10)

上述公式即為底泥某營(yíng)養(yǎng)鹽釋放通量公式。則底泥釋放強(qiáng)度公式為

γ′=R/A

(11)

釋放速率公式為

γ=R/(tA)

(12)

式中：R為底泥某營(yíng)養(yǎng)鹽試驗(yàn)期的釋放通量，mg;V為泥樣上方水的體積，L;Cn為第n次(試驗(yàn)終止時(shí))采樣時(shí)容器內(nèi)水中營(yíng)養(yǎng)鹽濃度，mg/L;C0為泥樣上覆水初始某營(yíng)養(yǎng)鹽濃度，mg/L;Vi為每次采水樣體積，L;Ci-1為第i-1次采樣時(shí)容器內(nèi)上覆水中某營(yíng)養(yǎng)鹽濃度，mg/L;Cai-1為第i-1次所添加的水體中某營(yíng)養(yǎng)物質(zhì)濃度，mg/L;γ′為釋放強(qiáng)度，mg/m2;γ為釋放速率，mg/(m2·d)；t為釋放時(shí)間，d;A為與水接觸的沉積物表面積，m2。

經(jīng)過(guò)上面推導(dǎo)明確了各符號(hào)代表的含義，驗(yàn)證了底泥營(yíng)養(yǎng)鹽釋放量的經(jīng)典計(jì)算公式。

2.2 試驗(yàn)結(jié)果分析

2.2.1 溫度(T)對(duì)底泥總磷(TP)釋放影響試驗(yàn)

設(shè)計(jì)pH=7，DO=8.0mg/L，試驗(yàn)溫度15℃、20℃、25℃、30℃對(duì)底泥總磷(TP)釋放的影響。

根據(jù)上述推導(dǎo)的釋放計(jì)算公式計(jì)算出pH值和DO一定時(shí)不同溫度下總磷(TP)的釋放強(qiáng)度和釋放速率，見(jiàn)圖3、圖4?？偭揍尫艔?qiáng)度在不同溫度下均為上升趨勢(shì)，25℃、30℃釋放強(qiáng)度比15℃、20℃上升快，30℃時(shí)釋放強(qiáng)度上升最快，TP釋放速率初始下降非常快，48h后下降減緩，隨時(shí)間增加釋放速率趨于穩(wěn)定，不再下降，30℃時(shí)釋放速率最大。從圖3、圖4可以看出，溫度升高有利于總磷的釋放，并且溫度升高，微生物的活性增強(qiáng)，促進(jìn)生物擾動(dòng)、礦化作用和厭氧轉(zhuǎn)化，導(dǎo)致間隙水耗氧增加，使環(huán)境向還原態(tài)轉(zhuǎn)化，利于Fe3+向Fe2+轉(zhuǎn)化，促使底泥中磷釋放。此外，微生物的活動(dòng)還可使底泥中的有機(jī)態(tài)磷轉(zhuǎn)化成無(wú)機(jī)態(tài)的磷酸鹽而得以釋放。太行湖的最高溫出現(xiàn)在8月(30℃左右)，因此，底泥中總磷含量此月也最高，總磷釋放均高于其他月份。

圖3 不同溫度下總磷(TP)釋放強(qiáng)度(γ′)曲線

圖4 不同溫度下總磷(TP)速率(γ)曲線

2.2.2 pH值對(duì)底泥總磷(TP)釋放影響試驗(yàn)

設(shè)計(jì)溫度20℃，DO=8.0mg/L，試驗(yàn)不同pH值(5、7、9)時(shí)對(duì)底泥總磷(TP)的釋放影響(通過(guò)多年DO檢測(cè)，2號(hào)區(qū)DO較高，在5～13mg/L之間，DO隨溫度升高而降低，隨湖區(qū)水深而降低，通過(guò)現(xiàn)場(chǎng)監(jiān)測(cè)，檢測(cè)結(jié)果見(jiàn)表1)，通過(guò)太行湖區(qū)BOD5檢測(cè)，DO=8.0mg/L時(shí)，20℃條件下BOD5=0.5左右，DO降低可以接受。7月湖區(qū)2號(hào)位置DO檢測(cè)結(jié)果見(jiàn)表4。

表4 7月湖區(qū)2號(hào)位置DO檢測(cè)結(jié)果

據(jù)上述釋放計(jì)算式得20℃，DO=8.0mg/L時(shí)不同pH值條件下總磷的釋放強(qiáng)度和釋放速率(見(jiàn)圖5、圖6)?？偭揍尫艔?qiáng)度(γ′)在不同pH值條件下均呈上升趨勢(shì)；pH=10時(shí)平均釋放強(qiáng)度最大，明顯偏高于其他組，初始下降速度也最快，其次為pH=9時(shí)釋放速率較大,pH=7、pH=5平均釋放速率最小，趨勢(shì)基本一致。60h后釋放速率趨于穩(wěn)定。

圖5 不同pH值條件總磷(TP)釋放強(qiáng)度曲線

圖6 不同pH值條件總磷(TP)釋放速率曲線

2.2.3 溶解氧(DO)對(duì)底泥總磷(TP)釋放影響試驗(yàn)

設(shè)計(jì)溫度20℃，pH=7，試驗(yàn)不同DO條件下(DO指試驗(yàn)初始值，且常年DO監(jiān)測(cè)值較高，底泥耗氧試驗(yàn)，基本是15h后略有下降，文中不再列舉)對(duì)底泥總磷(TP)的釋放影響。由圖7、圖8可知，總磷釋放強(qiáng)度在不同DO條件下均呈上升趨勢(shì)，但處于低氧(DO=2.0mg/L)時(shí)釋放強(qiáng)度和釋放速率最大。好氧(DO=8.0mg/L)時(shí)釋放速率最小，60h后釋放速率均趨于穩(wěn)定。

圖7 不同DO條件下總磷釋放強(qiáng)度曲線

圖8 不同DO條件下總磷釋放速率曲線

但在好氧狀態(tài)下，不是所有形態(tài)的磷都不會(huì)釋放，底泥中的有機(jī)質(zhì)就會(huì)受到微生物的好氧分解，使不溶性的有機(jī)磷變成無(wú)機(jī)磷，易于釋放。只是由于底泥中的有機(jī)磷含量一般較低，或是產(chǎn)生的腐殖質(zhì)的吸附性釋放量不大，所以，總的影響結(jié)果是在好氧狀態(tài)下，底泥中磷的釋放量會(huì)降低。

2.3 總磷釋放強(qiáng)度數(shù)學(xué)模型推求

2.3.1 數(shù)學(xué)模型一

自變量為溫度(T)、pH值、溶解氧(DO)，因變量為總磷(TP)平均釋放速率(γ)，建立多元線性回歸模型，借助SPSS軟件實(shí)現(xiàn)，則數(shù)學(xué)模型如下：

γTP=-45.429+0.953(T)+4.2(pH)+1.295(DO)

(13)

模型復(fù)相關(guān)系數(shù)：

R=0.754a

式中：γTP為總磷(TP)的平均釋放速率，mg/(m2·d)。

方差檢驗(yàn)：假設(shè)顯著性水平a=0.05時(shí)(即1-a=95%)，自變量m=3，樣本n=25，則自由度(df)為(m,n-m-1)=(3，21)，由F檢驗(yàn)臨界值表(α=0.05(a))查得F0.05(3，21)=3.073，由表5得知方差檢驗(yàn)F=9.249>3.073，概率P(sig)=0.000a<0.05，說(shuō)明方程回歸顯著。

2.3.2 數(shù)學(xué)模型二

將釋放試驗(yàn)的初始濃度值(C0)引入模型，自變量為C0、溫度(T)、pH值、溶解氧(DO)，因變量為總磷(TP)的平均釋放速率(γ)，則多因子線性回歸數(shù)學(xué)模型見(jiàn)式(5)，借助軟件spss實(shí)現(xiàn)，其數(shù)學(xué)模型如下：

γTP=8.897+0.43(T)+2.087(pH)+0.064(DO)-1136.97(C0)

(14)

模型復(fù)相關(guān)系數(shù)：

R=0.844a

式中:γTP為總磷(TP)的平均釋放速率，mg/(m2·d)。

總磷(TP)多元線性回歸模型方差檢驗(yàn)：假設(shè)顯著性水平a=0.05時(shí)(即1-a=95%)，自變量m=4，樣本n=25，則自由度(df)為(m,n-m-1)=(4，20),由F檢驗(yàn)臨界值表,α=0.05(a),查得F0.05(4，20)=2.866，由表4得知方差檢驗(yàn)F=12.364>2.866，概率P(sig)=0.000a<0.05，說(shuō)明回歸方程顯著。數(shù)學(xué)模型一、模型二擬合結(jié)果見(jiàn)圖9、圖10。

圖9 模型一總磷(TP)釋放速率擬合曲線

2.4 模型外推預(yù)測(cè)

模型檢驗(yàn)合格后可用于外推預(yù)測(cè)。分別用模型一和模型二進(jìn)行外推預(yù)測(cè)，漳澤水庫(kù)底泥總磷(TP)年釋放量按下式計(jì)算：

圖10 模型二總磷(TP)釋放速率擬合曲線

(15)

式中：W為湖庫(kù)磷釋放量，t；γi為第i種條件下氮釋放速率，mg/(m2·d)；ΔTi為第i種條件所代表的時(shí)間段，d；A為湖庫(kù)水面面積，m2。

以2018年實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)，預(yù)測(cè)2018年漳澤水庫(kù)沉積物中總磷的月平均釋放速率γTN和釋放量W，見(jiàn)表5。

表5 2018年沉積物中TP釋放量預(yù)測(cè)

從表5可以看出，如采用模型一，漳澤水庫(kù)沉積物中總磷表現(xiàn)為全年釋放反應(yīng)，合計(jì)釋放量57.40t；如采用模型二，5月和8月沉積物中總磷表現(xiàn)為釋放反應(yīng)，合計(jì)釋放量5.20t，其余月份表現(xiàn)為負(fù)釋放，沉積量為182.00t，全年合計(jì)沉積量177.00t。

2.5 合理性驗(yàn)證

調(diào)查顯示，漳澤水庫(kù)2018年滲漏量和下泄流量總計(jì)1.25億m3，2018年底蓄水量為8487萬(wàn)m3，較2017年底減少3677萬(wàn)m3。庫(kù)區(qū)代表斷面總磷平均濃度由2017年的0.061mg/L降至2018年的0.054mg/L。外部點(diǎn)源污染物總磷為22.60t/a，面源污染物總磷為840.00t/a。

考慮到面源污染中氨氮以顆粒物攜帶為主，多在進(jìn)入水庫(kù)后先行沉淀，而釋放試驗(yàn)釋放速率則以監(jiān)測(cè)上覆水中TP濃度得出，為更為準(zhǔn)確地驗(yàn)證內(nèi)源釋放量，調(diào)查了入庫(kù)站點(diǎn)的污染量，其中濁漳河南源、絳河以及碧頭河2018年TP入庫(kù)量合計(jì)79.50t。通過(guò)質(zhì)量守恒驗(yàn)證，2018年漳澤水庫(kù)的外源總磷污染量為107.00t/a，通過(guò)下泄和滲漏僅排出TP污染物6.75t/a，由于水體總磷濃度和蓄水量減少，水體中溶解態(tài)TP量也較上年減少2.72t/a，其內(nèi)源釋放量約為-103.00t/a，即水體中TP向沉積物中沉淀103.00t/a，與模型二結(jié)論基本一致。因此采用模型二外推結(jié)論，即2018年漳澤水庫(kù)沉積物總磷釋放量小于沉積量，表現(xiàn)為水體中溶解態(tài)總磷向沉積物中的沉積吸附反應(yīng)，全年沉積總磷103.00t。

3 結(jié) 論

通過(guò)分析太行湖底泥中總磷含量的時(shí)空分布特征及在靜態(tài)試驗(yàn)條件下改變環(huán)境因子模擬太行湖底泥總磷(TP)釋放，得出如下結(jié)論：

a.太行湖底泥受洪水影響8月總磷含量最高，10月含量最低，庫(kù)區(qū)西側(cè)近絳河支流區(qū)域4月TP含量較低，比10月略低?？臻g分布為上游入湖口與下游壩址東側(cè)泄水洞區(qū)域TP含量為高值區(qū)，湖區(qū)TP含量存在明顯帶狀分布，總體為湖區(qū)東側(cè)高于西側(cè)區(qū)。

b.底泥釋放TP強(qiáng)度隨溫度升高而增加；隨pH值增加釋放強(qiáng)度增加；隨DO減小釋放強(qiáng)度增加。

c.太行湖底泥總磷(TP)平均釋放強(qiáng)度模型，試驗(yàn)25組樣本建模為γTP=-45.429+0.953(T)+4.2(pH值)+1.295(DO)，R=0.754a，試驗(yàn)組14組樣本建模為γTP=8.897+0.43(T)+2.087(pH值)+0.064(DO)-1136.97(CO)(DO)，R=0.844a。

d.經(jīng)模型二預(yù)測(cè)2018年底泥釋放磷量為5.20t,總體沉積177.00t，占外源TP輸入量的4.8%。

通過(guò)對(duì)湖區(qū)底泥總磷含量分析和基于化學(xué)條件下總磷靜態(tài)釋放模擬，外源(面源為主)仍是太行湖庫(kù)主要污染源，在外源磷輸入截?cái)嗪螅瑑?nèi)源底泥釋磷將成為其水體富磷主要因素，會(huì)引發(fā)二次水體磷污染，再度造成水體富營(yíng)養(yǎng)化。本次監(jiān)測(cè)湖區(qū)底泥總磷分布特征及釋磷模擬試驗(yàn)，為湖區(qū)污染治理及生態(tài)修復(fù)提供了技術(shù)參考。限于各方面條件，沒(méi)能進(jìn)行室外現(xiàn)場(chǎng)模擬，本次室內(nèi)試驗(yàn)只進(jìn)行了化學(xué)條件下靜態(tài)釋放模擬研究，是本文的缺陷所在，以后將進(jìn)一步對(duì)太行湖區(qū)底泥展開(kāi)動(dòng)態(tài)條件和微生物影響條件下的磷釋放進(jìn)行研究。