陳曉燕，何秉宇，劉江，茉莉得爾金斯汗，陳磊

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，新疆烏魯木齊830046；2.新疆大學(xué)智慧城市與環(huán)境建模自治區(qū)高校重點(diǎn)實(shí)驗(yàn)室，新疆烏魯木齊830046；3.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，新疆烏魯木齊830052；4.新疆環(huán)境保護(hù)科學(xué)研究院，新疆烏魯木齊830011，5.新疆環(huán)境污染監(jiān)控與風(fēng)險(xiǎn)預(yù)警重點(diǎn)實(shí)驗(yàn)室，新疆烏魯木齊830011)

水環(huán)境質(zhì)量下降的主要原因之一是有機(jī)污染物污染[1]．艾里克湖位于新疆克拉瑪依市東北部，對(duì)克拉瑪依市整體生態(tài)環(huán)境狀況有直接影響，對(duì)當(dāng)?shù)厣鐣?huì)經(jīng)濟(jì)的可持續(xù)發(fā)展有一定影響[2]．不同類型的天然水體對(duì)污染物自然降解有其特性[3?5]，20世紀(jì)80年代至今，雖然許多學(xué)者對(duì)國(guó)內(nèi)外多條河流的有機(jī)污染物降解速率進(jìn)行了大量研究[6]，但湖泊中污染物降解機(jī)理較河流中更為復(fù)雜，而針對(duì)湖泊污染物的降解僅有零星的研究[7,8]．有關(guān)新疆湖泊污染物降解系數(shù)的研究未見(jiàn)報(bào)道，考慮到新疆地理位置、氣候以及水體含鹽量高等特點(diǎn)，直接將《國(guó)家地表水水環(huán)境容量核定》中給定的污染物降解系數(shù)應(yīng)用于新疆湖泊中，其數(shù)值模擬結(jié)果或許存在一定的誤差，從而可能導(dǎo)致水環(huán)境容量測(cè)算誤差．

由于周圍融雪水的匯入，融凍期湖泊水量增加迅速，污染物入湖量增大，水溫、pH值、溶解氧等因素變化劇烈，導(dǎo)致湖泊污染物降解系數(shù)變化較大．本文選取化學(xué)需氧量（CODcr）、氨氮、總氮、總磷作為艾里克湖水環(huán)境中有機(jī)物污染的主要指標(biāo)[9]．研究有機(jī)污染物在融凍期（2月底至4月上旬）的綜合降解系數(shù)，探究降解系數(shù)的時(shí)空變化規(guī)律及其影響因素．

1 材料與方法

1.1 研究區(qū)概況

艾里克湖（85?450～85?520E，45?510～45?590N）位于新疆維吾爾自治區(qū)克拉瑪依市烏爾禾區(qū)境內(nèi)，地處準(zhǔn)噶爾盆地西部低山丘陵的東緣，是白楊河的尾閭湖．白楊河每年4-5月有積雪融水形成的洪水涌入艾里克湖，湖面也因此擴(kuò)大[10]．湖泊長(zhǎng)約12.4 km，寬約為3.5 km，面積約52.4 km2[11]．艾里克湖流域是典型的溫帶干旱區(qū)大陸性荒漠氣候，降水稀少，冬寒夏暖，氣溫年較差與日較差大，且春季、秋季多大風(fēng)，七八級(jí)大風(fēng)年均68天，平均風(fēng)速達(dá)3.5 m/s[12]．

1.2 研究方法

1.2.1 實(shí)驗(yàn)方法及時(shí)間

目前有關(guān)湖泊降解系數(shù)的確定方法主要有現(xiàn)場(chǎng)實(shí)際監(jiān)測(cè)法、室內(nèi)模擬實(shí)驗(yàn)法、歷史資料類比分析法、實(shí)測(cè)資料反推法、常規(guī)監(jiān)測(cè)資料估算法等[13]．室內(nèi)模擬實(shí)驗(yàn)法估算有機(jī)污染物降解系數(shù)屬于靜態(tài)實(shí)驗(yàn)，在模擬過(guò)程中無(wú)法體現(xiàn)水文、水力、溫度、河道特征、底泥性質(zhì)、季節(jié)變化等諸多綜合影響因素[14]．本文根據(jù)該湖泊的實(shí)際情況，采用現(xiàn)場(chǎng)實(shí)際監(jiān)測(cè)法來(lái)確定融凍期艾里克湖有機(jī)污染物的降解系數(shù)．

研究選擇艾里克湖融凍期（2月底至4月上旬）的有機(jī)污染物降解系數(shù)，其中采樣時(shí)間為2016年2月29日、3月7日、3月14日、3月21日、3月28日、4月4日．

1.2.2 采樣位置

根據(jù)以往的艾里克湖生態(tài)環(huán)境保護(hù)試點(diǎn)總體實(shí)施方案，以及對(duì)水質(zhì)監(jiān)測(cè)采樣時(shí)的難易程度，最終確定五個(gè)具有代表性的采樣點(diǎn)（即入湖口SH03-1、湖右岸SH03-2、湖左岸SH03-3、湖心部分SH03-4、湖尾處SH03-5），其中每個(gè)采樣點(diǎn)重復(fù)3次采樣，共計(jì)15個(gè)水質(zhì)采樣點(diǎn)．艾里克湖水質(zhì)采樣位置如圖1．

圖1 艾里克湖水質(zhì)采樣位置示意圖

1.2.3 現(xiàn)場(chǎng)監(jiān)測(cè)

現(xiàn)場(chǎng)采樣時(shí)，同時(shí)開(kāi)展GPS定位、水溫、pH值、溶解氧、流速等基礎(chǔ)資料的測(cè)量工作．采集水樣按《水和廢水監(jiān)測(cè)分析方法（第四版）》中的要求保存，帶至實(shí)驗(yàn)室待測(cè)．

1.2.4 降解系數(shù)的確定

國(guó)內(nèi)外研究表明，污染物的生物降解原理符合米氏方程，污染物濃度遠(yuǎn)大于米氏常數(shù)時(shí)，降解規(guī)律符合零級(jí)反應(yīng)，污染物濃度遠(yuǎn)小于米氏常數(shù)時(shí)，則符合一級(jí)反應(yīng)[15]．湖泊中污染物濃度一般小于米氏常數(shù)，所以其降解系數(shù)可用一級(jí)反應(yīng)動(dòng)力學(xué)模型進(jìn)行計(jì)算．依據(jù)研究[16,17]，本文對(duì)KCOD、K氨氮、K總氮、K總磷的確定，采用一級(jí)反應(yīng)動(dòng)力學(xué)模型進(jìn)行計(jì)算．一級(jí)反應(yīng)動(dòng)力學(xué)模型如下

式中：C為t時(shí)刻某污染物的濃度(mg/L)；K為污染物降解系數(shù)；t為反應(yīng)時(shí)間(d)．

將式（1）兩邊同時(shí)積分

式中：C0為某污染物的初始濃度(mg/L)．

由（3）得污染物降解系數(shù)K值[18]．

2 結(jié)果與討論

2.1 降解系數(shù)的時(shí)空變化規(guī)律

湖泊特性表明其水質(zhì)存在時(shí)空變化規(guī)律．本文從時(shí)間和空間兩個(gè)角度研究融凍期艾里克湖水環(huán)境中有機(jī)污染物的降解系數(shù)，探討其變化規(guī)律．

2.1.1 時(shí)間變化

融凍期的艾里克湖主要有機(jī)污染物KCOD的變化范圍為0.000 92 d?1～0.019 70 d?1，K氨氮的變化范圍為 0.000 58 d?1～0.037 19 d?1；K總氮變化范圍為0.001 84 d?1～0.038 15 d?1；K總磷變化范圍為0.000 59 d?1～0.073 85 d?1．

由圖2分析得出，KCOD、K氨氮、K總氮融凍初期呈遞增趨勢(shì)，但在3月21日或28日之后呈現(xiàn)下降趨勢(shì)，造成這種現(xiàn)象的原因可能是該時(shí)間段內(nèi)該地區(qū)氣溫回升過(guò)快．并且結(jié)合當(dāng)?shù)刈匀坏乩砬闆r，該湖區(qū)位于風(fēng)口，其風(fēng)力較大，湖水復(fù)氧能力較高，由此加快了有機(jī)物的降解．然而K總磷隨時(shí)間的變化規(guī)律并不明顯．

圖2 降解系數(shù)隨時(shí)間的變化

圖3 降解系數(shù)空間變化

2.1.2 空間變化

按照入湖口、湖右岸、湖左岸、湖心部分、湖尾處對(duì)該湖泊進(jìn)行空間分類，并進(jìn)行各個(gè)采樣點(diǎn)降解系數(shù)變化的討論分析，變化規(guī)律如圖3所示．

（1）KCOD：不同位置的KCOD變化呈現(xiàn)三種形態(tài)，即波動(dòng)型、穩(wěn)定型、增長(zhǎng)型．入湖口為波動(dòng)型，因入湖口位置水質(zhì)變化大，致使KCOD波動(dòng)較大．湖右岸為穩(wěn)定型，其原因是該位置地處艾里克湖西南，入湖的污染物與原水質(zhì)混合平緩，并未產(chǎn)生較大波動(dòng)．湖左岸、湖心部分、湖尾處為增長(zhǎng)型，即KCOD呈現(xiàn)有規(guī)律的增長(zhǎng)，然而3月28日-4月4日KCOD明顯下降，說(shuō)明入湖的污染物水量對(duì)該湖KCOD影響較弱，不足以影響總的降解系數(shù)．

（2）K氨氮：整體上3月7日-28日五個(gè)位置的K氨氮都趨于穩(wěn)定增長(zhǎng)，說(shuō)明湖水中污染物降解與水溫有直接關(guān)系，或由于入湖污染物量較少，不足以直接影響到氨氮的降解速率．4月4日五個(gè)位置的K氨氮都呈下降趨勢(shì)，造成這種現(xiàn)象的原因可能是該時(shí)間段內(nèi)匯入的污染物氨氮濃度偏高，對(duì)K氨氮的影響較大．湖尾處在21日的K氨氮就呈現(xiàn)下降趨勢(shì)，根據(jù)水力動(dòng)力學(xué)特性分析其位于入湖主流上，受入湖水流影響較大，該位置的降解系數(shù)轉(zhuǎn)折點(diǎn)出現(xiàn)較早．

（3）K總氮：入湖口K總氮變化較大，沒(méi)有穩(wěn)定的規(guī)律，說(shuō)明一直有污染源的匯入，對(duì)K總氮產(chǎn)生影響．3月7日-28日湖右岸、湖左岸、湖心部分、湖尾處符合常規(guī)降解趨勢(shì)，說(shuō)明入湖水量較少，K總氮呈現(xiàn)有規(guī)律的增長(zhǎng)，并在21日或28日出現(xiàn)K總氮的最大值．3月28日-4月4日五個(gè)位置的K總氮都出現(xiàn)下降趨勢(shì)，說(shuō)明入湖的污染物量超過(guò)了該湖的水體自凈能力．

（4）K總磷：可看出湖心部分的K總磷趨于穩(wěn)定，入湖口、湖右岸、湖左岸、湖尾處都存在較大的變化，規(guī)律不明顯，說(shuō)明一直都有污染物的匯入，入湖流量水質(zhì)與原水質(zhì)混合不均勻．

五個(gè)監(jiān)測(cè)點(diǎn)的KCOD、K氨氮存在一定規(guī)律性變化，可以歸納為波動(dòng)型、穩(wěn)定型、增長(zhǎng)型；K總氮、K總磷變化規(guī)律不明顯．說(shuō)明匯入艾里克湖的外界污染物水質(zhì)中化學(xué)需氧量、氨氮的含量相對(duì)穩(wěn)定；總氮、總磷含量不穩(wěn)定．

2.1.3 降解系數(shù)變化時(shí)空分析

利用SPSS軟件，以時(shí)間、空間為雙因子，各個(gè)降解系數(shù)為因變量，進(jìn)行雙因素方差分析，結(jié)果見(jiàn)表1．KCOD兩組間的時(shí)間變化非常顯著（P≤0.01），KCOD兩組間的空間變化為不顯著（P>0.05）；K氨氮兩組間的時(shí)間變化差異為非常顯著（P≤0.01），K氨氮兩組間的空間變化為不顯著（P>0.05）；K總氮兩組間的時(shí)間變化、空間變化差異為非常顯著（P≤0.01）；K總磷兩組間的時(shí)間變化為非常顯著（P≤0.01），K總磷兩組間的空間變化為不顯著（P>0.05）．

表1 降解系數(shù)與時(shí)空的雙因素方差分析

在融凍期內(nèi)，艾里克湖的KCOD、K氨氮、K總氮、K總磷都會(huì)隨著時(shí)間變化而遞增；對(duì)于不同空間位置上的KCOD、K氨氮、K總氮、K總磷并未呈現(xiàn)一定的規(guī)律性．

2.2 影響降解系數(shù)變化的因素

研究表明，水體中污染物質(zhì)的降解是一個(gè)很復(fù)雜的過(guò)程，受諸多因素的影響，主要包括水中污染物的特性、pH值、污染物初始濃度、水力特性、微生物群落、懸浮固體、溫度、溶解氧等幾個(gè)方面[19]．

利用R語(yǔ)言技術(shù)進(jìn)行降解系數(shù)與水溫、pH值、溶解氧、流速的相關(guān)性分析結(jié)果如表2所示．得出KCOD與流速間的差異為非常顯著；K氨氮與pH值、流速間的差異為非常顯著；K總氮與水溫、pH值、流速間的差異為非常顯著；然而K總磷分別與這四種因素間的差異為不顯著．

表2 降解系數(shù)與水溫、pH值、溶解氧、流速的相關(guān)性

利用R語(yǔ)言技術(shù)進(jìn)行多元線性回歸分析如表3所示（其中X1、X2、X3、X4分別為水溫、pH值、溶解氧、流速）．由于影響因子作用于湖泊水質(zhì)，因此在單因子分析的基礎(chǔ)上，進(jìn)行多變量相關(guān)分析，尋找影響該湖泊水質(zhì)降解系數(shù)變化的關(guān)鍵因子．

影響因子與KCOD、K氨氮、K總氮的多元線性回歸模型都達(dá)到顯著性水平(P<0.05)，其中KCOD、K氨氮、K總氮受四種因子影響較小(R2<0.5)，K總磷的多元線性回歸模型未達(dá)到顯著性水平(P>0.05)且受四種因子影響較小(R2<0.5)．KCOD、K氨氮、K總氮、K總磷都與水溫、溶解氧正相關(guān)；KCOD、K總磷與pH值正相關(guān)，K氨氮、K總氮與pH值負(fù)相關(guān)；KCOD、K氨氮、K總氮與流速正相關(guān)，K總磷與流速負(fù)相關(guān)．

圖4 Matlab的逐步回歸分析圖

從多元線性回歸分析中已經(jīng)得出降解系數(shù)與各影響因子間的相關(guān)關(guān)系，但需要進(jìn)一步確定影響該湖泊水環(huán)境中有機(jī)污染物降解系數(shù)的最主要因子，進(jìn)而利用Matlab的逐步回歸分析如圖4所示，最終建立逐步回歸模型：YCOD=0.004 271 7+0.011 431 5X4；Y氨氮=0.000 550 999+0.031 574 4X4；Y總氮=0.006 101 53+0.020 721 9X4；總磷的相關(guān)分析未建立成方程．綜上分析，流速是影響KCOD、K氨氮、K總氮的最主要因子，而K總磷的變化并不符合常規(guī)，未找出最主要的因子．

2.2.1 湖泊水溫的影響

水溫升高可以加速湖泊中有機(jī)污染物的分解、氧化反應(yīng)、微生物的活性，有利于污染物的降解．結(jié)果顯示：水溫對(duì)有機(jī)污染物的降解系數(shù)影響較大，一般來(lái)說(shuō)，北方水體夏季的降解系數(shù)要比冬季的降解系數(shù)高出1～2倍，這也是有些河流、湖泊枯水期水質(zhì)惡化的一個(gè)重要原因．本文研究時(shí)間段為融凍期，水溫變化從0～12?C遞增，每次采樣增幅僅有2～4?C．通過(guò)數(shù)據(jù)分析，KCOD、K氨氮、K總氮、K總磷都與水溫正相關(guān).

2.2.2 湖泊水體中pH值的影響

pH值是水質(zhì)中氫離子活度的一種標(biāo)度，也就是通常意義上水質(zhì)酸堿程度的衡量標(biāo)準(zhǔn)．根據(jù)監(jiān)測(cè)數(shù)據(jù)顯示，該時(shí)間段內(nèi)該水體中污染物pH值范圍為8.12～8.23，其變化并不明顯．KCOD、K總磷與pH值正相關(guān)，K氨氮、K總氮與pH值負(fù)相關(guān)．

2.2.3 湖泊水體中溶解氧的影響

溶解氧值是研究水自凈能力的一個(gè)指標(biāo)．研究時(shí)段水體中的溶解氧范圍8.01～9.66 mg/L，進(jìn)一步分析表明，KCOD、K氨氮、K總氮、K總磷與溶解氧濃度呈現(xiàn)正相關(guān)關(guān)系．

2.2.4 湖泊水文特征—流速的影響

湖泊水文特征包括流量、流速、湖寬、湖深、底泥等．由于研究的是融凍期的湖泊水質(zhì)，因而流量、湖寬、湖深、底泥等對(duì)該降解系數(shù)的影響不大．根據(jù)現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)，流速?gòu)?.1 m/s逐漸增長(zhǎng)至0.8 m/s，并且通過(guò)數(shù)據(jù)分析，KCOD、K氨氮、K總氮與流速正相關(guān)，K總磷與流速負(fù)相關(guān)，降解系數(shù)會(huì)隨著流速的增大而增加．

3 結(jié)論

（1）融凍期艾里克湖主要有機(jī)污染物降解系數(shù)變化較大．KCOD變化范圍為0.000 92 d?1～0.019 70 d?1,K氨氮變化范圍為0.000 58 d?1～0.037 19 d?1；K總氮變化范圍為0.001 84 d?1～0.038 15 d?1；K總磷變化范圍為0.00 059 d?1～0.073 85 d?1；

（2）融凍期艾里克湖主要有機(jī)污染物降解系數(shù)變化呈現(xiàn)兩種模式：規(guī)律性模式與非規(guī)律性模式．KCOD、K氨氮在融凍期表現(xiàn)出較好的規(guī)律性；K總氮、K總磷變化規(guī)律不明顯；

（3）融凍期艾里克湖主要有機(jī)污染物降解系數(shù)時(shí)間變化特征顯著，空間變化未呈現(xiàn)較好的規(guī)律性；

（4）融凍期艾里克湖主要有機(jī)污染物降解系數(shù)變化影響因素分析表明：水溫、溶解氧是主要影響因子，流速、pH值的影響有正負(fù)方向差異之分．

參考文獻(xiàn)：

[1]游雪靜,張玉珍,蘇玉萍,等.閩江流域水體氨氮降解系數(shù)實(shí)驗(yàn)?zāi)M研究[J].亞熱帶資源與環(huán)境學(xué)報(bào),2014(1):61-67.

[2]董莉莉.艾里克湖濕地生態(tài)保護(hù)指標(biāo)簡(jiǎn)析[J].河南水利與南水北調(diào),2015(5):25-26.

[3]倪浩清,李福田.水環(huán)境中生物化學(xué)時(shí)均反應(yīng)率的研究[J].水利學(xué)報(bào),2010,39(1):37-46.

[4]Leu H G,Ouyang C F,Pai T Y.Ef f ects of f l ow velocity and depth on the rates of reaeration and bod removal in a shallow open channel[J].Water Science&Technology,1997,35(8):57-67.

[5]Novotny V,Krenkel P A.A waste assimilative capacity model for a shallow,turbulent stream[J].Water Research,1975,9(2):233-241.

[6]周國(guó)勝.老府河水體中有機(jī)污染物降解規(guī)律研究[J].湖北理工學(xué)院學(xué)報(bào),2010,26(2):13-14.

[7]范成新,陳宇煒.太湖梅梁灣南部水體有機(jī)污染物降解表觀動(dòng)力學(xué)初步分析[J].湖泊科學(xué),1998,10(4):48-52.

[8]宋學(xué)宏,邴旭文,孫麗萍,等.陽(yáng)澄湖養(yǎng)殖水體COD降解動(dòng)力學(xué)研究[J].安徽農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,37(2):328-332.

[9]郭儒,李宇斌,富國(guó).河流中污染物衰減系數(shù)影響因素分析[J].氣象與環(huán)境學(xué)報(bào),2008,24(1):56-59.

[10]閆培鋒,周華榮,劉宏霞.白楊河—艾里克湖濕地土壤理化性質(zhì)的空間分布特征[J].干旱區(qū)研究,2008,25(3):406-412.

[11]胡云玲,焦鍵,劉江,等.分形理論在艾里克湖富營(yíng)養(yǎng)化評(píng)價(jià)中的應(yīng)用[J].新疆環(huán)境保護(hù),2015(1):8-12.

[12]姜亮亮,焦鍵,程艷,等.艾里克湖流域湖泊及植被生態(tài)需水研究[J].干旱區(qū)資源與環(huán)境,2016,30(4):108-114.

[13]楊亮秀,楊常亮,劉衛(wèi)紅,等.滇池流域盤龍江化學(xué)需氧量與總氮降解系數(shù)的測(cè)定研究[C]//2011中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第一卷).2011:298-401.

[14]劉洪燕.北方小基流河流氨氮綜合降解系數(shù)特征研究[J].能源與節(jié)能,2016(9):109-110.

[15]馮帥,李敘勇,鄧建才.平原河網(wǎng)典型污染物生物降解系數(shù)的研究[J].環(huán)境科學(xué),2016,37(5):1724-1733.

[16]馮帥,李敘勇,鄧建才.太湖流域上游河網(wǎng)污染物降解系數(shù)研究[J].環(huán)境科學(xué)學(xué)報(bào),2016,36(9):3127-3136.

[17]王雅鈺,劉成剛,黃天寅,等.不同處理技術(shù)對(duì)河流污染物降解系數(shù)影響的比較[J].環(huán)境工程學(xué)報(bào),2014,8(3):917-923.

[18]李國(guó)權(quán).三峽庫(kù)區(qū)排污口下游污染物降解規(guī)律的研究[D].重慶:重慶交通大學(xué),2013.

[19]中國(guó)環(huán)境規(guī)劃院.全國(guó)水環(huán)境容量核定技術(shù)指南[R].北京:中國(guó)環(huán)境規(guī)劃院,2003:55-56.