李春培，王有寧，劉 婕，田 嶶，李丙官，朱天云

(湖北工程學(xué)院 生命科學(xué)技術(shù)學(xué)院，湖北 孝感 432000)

隨著社會(huì)經(jīng)濟(jì)的發(fā)展，我國農(nóng)資產(chǎn)品的消耗量日益增加，從而產(chǎn)生較多農(nóng)業(yè)面源污染問題[1]，如水體富營養(yǎng)化[2-5]、土壤板結(jié)[6]等。掌握氮磷排放狀況是解決農(nóng)業(yè)面源污染問題的前提之一，但目前關(guān)于農(nóng)場(chǎng)排水溝渠氮磷排放規(guī)律的研究報(bào)道還較為少見。

本文以湖北朱湖農(nóng)場(chǎng)排水溝渠為研究對(duì)象，通過遙感影像確定研究溝渠的地理位置，設(shè)定采樣點(diǎn)，判斷水流的流速和流向，然后對(duì)排水溝渠中N、P濃度進(jìn)行監(jiān)測(cè)采樣分析并運(yùn)用灰色關(guān)聯(lián)水質(zhì)評(píng)價(jià)方法，對(duì)排水溝渠進(jìn)行綜合評(píng)價(jià)，進(jìn)一步掌握溝渠排水的氮磷污染狀況，以期為農(nóng)業(yè)面源污染中氮磷的吸附攔截提供支持，為長江中下游水體富營養(yǎng)化治理提供參考。

1 材料與方法

1.1 研究區(qū)域概況

孝感市孝南區(qū)朱湖農(nóng)場(chǎng)地處江漢平原北部，地理坐標(biāo)為東經(jīng)113°51′～113°57′，北緯30°47′～30°51′，總面積43 km2，耕地面積2 333 hm2，養(yǎng)殖水面800 hm2，常住人口2.4萬。具有獨(dú)特的區(qū)位優(yōu)勢(shì)和貨暢其流的交通優(yōu)勢(shì)(見圖1)。

圖1 朱湖農(nóng)場(chǎng)溝渠斷面布設(shè)與位置

1.2 采樣點(diǎn)的布設(shè)

排水溝渠內(nèi)存在嚴(yán)重水華現(xiàn)象，排水溝渠植物生長較多，溝渠兩側(cè)建有綠化帶，部分溝渠兩側(cè)進(jìn)行水泥硬化。溝渠兩側(cè)在相同的管理制度下進(jìn)行管理，排水溝渠中部分水直排入府河，部分通過排入菱角湖，再匯入府河，最終匯入長江。

根據(jù)朱湖農(nóng)場(chǎng)排水體系分布特征，結(jié)合采樣布設(shè)的代表性和準(zhǔn)確性原則，設(shè)置控制斷面1、8，分別位于朱湖排水系統(tǒng)的入口以及出口，其次設(shè)置監(jiān)測(cè)斷面 2、3、4、5、6、7 于朱湖農(nóng)場(chǎng)主要排水溝渠(見圖1)，其目的是從時(shí)間和空間上監(jiān)測(cè)排水體系中的水質(zhì)。

1.3 樣品采集

定時(shí)定點(diǎn)進(jìn)行水樣采集。本研究采樣時(shí)間：第1次采樣為2019年5月16日，第6次采樣為2019年6月28日，其中第4次與第6次均是雨后采樣。每隔7天進(jìn)行一次采樣監(jiān)測(cè)，每次取樣時(shí)間為上午9時(shí)至11時(shí)。取樣時(shí)，在水體未被擾動(dòng)的情況下，取水面下10 cm水體，將其裝入預(yù)先清洗的500 mL玻璃瓶中，加入硫酸酸化保存樣品，并滴加兩滴氯仿以抑制微生物的活動(dòng)。

1.4 樣品處理及測(cè)定

將野外采集的水樣帶回實(shí)驗(yàn)室進(jìn)行指標(biāo)分析。主要監(jiān)測(cè)指標(biāo)為化學(xué)需氧量(COD)、總磷 (TP)、銨態(tài)氮(NH4+-N)、硝態(tài)氮(NO3--N)，各項(xiàng)指標(biāo)監(jiān)測(cè)時(shí)，都應(yīng)經(jīng)過0.45 μm濾膜進(jìn)行過濾后，在進(jìn)行下一步監(jiān)測(cè)。監(jiān)測(cè)分析方法參見《水和廢水監(jiān)測(cè)分析方法(第四版)》[7]。

1.5 數(shù)據(jù)處理及方法

數(shù)據(jù)處理采用Excel 2010和 IBM SPSS Statistics 24等數(shù)據(jù)處理軟件。

1.6 模型建立

灰色關(guān)聯(lián)法可以在數(shù)據(jù)少、樣本小的情況下進(jìn)行數(shù)據(jù)分析，通過對(duì)實(shí)測(cè)數(shù)據(jù)形成的數(shù)列與構(gòu)建的理想標(biāo)準(zhǔn)數(shù)列的關(guān)聯(lián)程度，來研究系統(tǒng)中子系統(tǒng)的關(guān)聯(lián)程度[8]。在水質(zhì)評(píng)價(jià)過程中，存在許多不可控制的自然因素影響實(shí)驗(yàn)監(jiān)測(cè)結(jié)果。本文選取傳統(tǒng)的灰色關(guān)聯(lián)法對(duì)湖北朱湖農(nóng)場(chǎng)排水溝渠進(jìn)行水質(zhì)評(píng)價(jià)。

首先，構(gòu)建標(biāo)準(zhǔn)數(shù)列和實(shí)測(cè)數(shù)列，標(biāo)準(zhǔn)數(shù)列以GB3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》[9]為基礎(chǔ)建立，實(shí)測(cè)數(shù)列以斷面監(jiān)測(cè)的數(shù)據(jù)為基礎(chǔ)建立。其次，進(jìn)行數(shù)據(jù)無量綱化處理，本文采用均值化[10]進(jìn)行無量綱化處理。再次，采用傳統(tǒng)的灰色關(guān)聯(lián)表達(dá)式進(jìn)行關(guān)聯(lián)系數(shù)與關(guān)聯(lián)度的計(jì)算。最后，進(jìn)行關(guān)聯(lián)度的排序，定義實(shí)測(cè)數(shù)列與標(biāo)準(zhǔn)數(shù)列的關(guān)聯(lián)度越高，其相關(guān)性越密切。

2 結(jié)果與分析

2.1 COD時(shí)間維度、空間維度差異分析

通過室內(nèi)實(shí)驗(yàn)指標(biāo)監(jiān)測(cè)，所測(cè)定的COD最大值為844.78 mg/L，最小值為51.20 mg/L，其平均值為249.82 mg/L。從時(shí)間維度上看，朱湖農(nóng)場(chǎng)排水溝渠的COD呈現(xiàn)逐漸增加的趨勢(shì)。第六次采樣前發(fā)生了較強(qiáng)烈的降雨事件，促使第1次采樣、第5次采樣與第6次采樣所測(cè)定的COD值存在顯著性差異。在第1次采樣至第2次采樣均為晴天，第4次采樣為降雨事件發(fā)生之后，由于農(nóng)田施肥和降雨事件雙重作用，使得第2次采樣與第4次采樣所測(cè)定值存在顯著性差異(表1)。由于第4次采樣與第6次采樣均為雨后采樣，所測(cè)定的值均有大幅度的上升，這表明降雨事件會(huì)促進(jìn)地表污染物質(zhì)通過地表徑流和地下滲透的方式進(jìn)入排水溝渠。

表1 COD時(shí)間維度差異分析

從空間維度上看，朱湖農(nóng)場(chǎng)排水溝渠各斷面的COD呈現(xiàn)波動(dòng)式變化趨勢(shì)。根據(jù)野外調(diào)查，其水流方向?yàn)橛蓽S河進(jìn)入控制斷面8，再到控制斷面1，通過控制斷面1匯入府河。其中控制斷面1、監(jiān)測(cè)斷面3、5、7較前一斷面所測(cè)定COD平均值較高，在這三個(gè)斷面存在生活污水排入口，這是導(dǎo)致COD值較高的原因之一。由表2可知，各斷面之間所測(cè)定的COD值無顯著差異。

表2 COD空間維度差異分析

2.2 NH4+-N時(shí)間維度、空間維度差異分析

通過室內(nèi)實(shí)驗(yàn)指標(biāo)監(jiān)測(cè)，所測(cè)定的NH4+-N最大值為0.87 mg/L，最小值為 0.04 mg/L，其平均值為0.22 mg/L。從時(shí)間上看，由施肥活動(dòng)、微生物硝化與反硝化以及降雨事件綜合影響，使得朱湖農(nóng)場(chǎng)排水溝渠的NH4+-N呈現(xiàn)先減小后增加的趨勢(shì)。其中第1次采樣至第4次采樣期內(nèi)，NH4+-N呈現(xiàn)減小趨勢(shì)，這是排水溝渠自凈能力的表現(xiàn)。第4次采樣至第6次采樣NH4+-N呈現(xiàn)逐漸增加的趨勢(shì)，這是在降雨事件發(fā)生之后，施肥活動(dòng)也隨之進(jìn)行，未被土壤膠體固定的NH4+-N隨地表徑流匯入排水溝渠，促使水體中NH4+-N增加。

由表3可知，第1次采樣與其他采樣之間存在顯著差異，第2次采樣與其他采樣之間存在顯著差異，第3次采樣與第4次采樣存在顯著差異，第4次采樣與第6次采樣存在顯著差異。

表3 NH4+-N時(shí)間維度差異分析

從空間維度上看，朱湖農(nóng)場(chǎng)排水溝渠各斷面的NH4+-N呈現(xiàn)波動(dòng)式變化趨勢(shì)。其中監(jiān)測(cè)斷面6與控制斷面1所監(jiān)測(cè)的NH4+-N增加峰值，這與施肥量有關(guān)。由于排水溝渠是厭氧環(huán)境，存在底泥微生物的硝化與反硝化作用，促使監(jiān)測(cè)斷面6至監(jiān)測(cè)斷面3呈現(xiàn)消減作用。由表4可知，監(jiān)測(cè)斷面3、4與其他斷面呈現(xiàn)顯著性差異，這兩個(gè)監(jiān)測(cè)斷面均值較小，這是由消減作用導(dǎo)致。

表4 NH4+-N空間維度差異分析

2.3 NO3--N時(shí)間維度、空間維度差異分析

通過室內(nèi)實(shí)驗(yàn)指標(biāo)監(jiān)測(cè)，所測(cè)定的NO3--N最大值為6.12 mg/L，最小值為0.49 mg/L，其平均值為2.40 mg/L。從時(shí)間上看，由施肥活動(dòng)、微生物硝化與反硝化以及降雨事件綜合影響，使得朱湖農(nóng)場(chǎng)排水溝渠的NO3--N呈現(xiàn)逐漸增加的趨勢(shì)。在降雨事件與施肥活動(dòng)發(fā)生之后，未被土壤膠體固定的NO3--N隨地表徑流進(jìn)入排水溝渠，使得排水溝渠中NO3--N大幅增加。由表5可知，第4、5、6采樣所測(cè)定的NO3--N無顯著差異，而第1次采樣與其他采樣所測(cè)定的NO3--N存在顯著差異。

表5 NO3--N時(shí)間維度差異分析

從空間維度上看，朱湖農(nóng)場(chǎng)排水溝渠各斷面的NO3--N呈現(xiàn)波動(dòng)式變化趨勢(shì),其中監(jiān)測(cè)斷面3所測(cè)的NO3--N均值最大，這與生活廢水排入量有關(guān)。由于排水溝渠是厭氧環(huán)境，存在底泥微生物的硝化與反硝化作用，促使各斷面呈現(xiàn)波動(dòng)式變化。由表6可知，監(jiān)測(cè)斷面2、4、5、6、7與控制斷面1之間無顯著差異，監(jiān)測(cè)斷面2、4、6、7與控制斷面1以及監(jiān)測(cè)斷面2之間無顯著差異，監(jiān)測(cè)斷面3與控制斷面8存在顯著差異。

表6 NO3--N空間維度差異分析

2.4 TP時(shí)間維度、空間維度差異分析

通過室內(nèi)實(shí)驗(yàn)指標(biāo)監(jiān)測(cè)，所測(cè)定的TP最大值為0.36 mg/L，最小值為0.01 mg/L，其平均值為0.16 mg/L。從時(shí)間上看，由施肥活降雨事件綜合影響，使得朱湖農(nóng)場(chǎng)排水溝渠的TP 呈現(xiàn)逐漸增加的趨勢(shì)。在降雨事件與施肥活動(dòng)發(fā)生之后，未被土壤膠體固定的TP隨地表徑流進(jìn)入排水溝渠，使得排水溝渠中TP大幅增加(表7)。由于施肥活動(dòng)與降雨事件的發(fā)生，使得第4次采樣、第6次采樣與其他采樣存在顯著差異。這表明降雨事件與排水溝渠中TP有較高的相關(guān)性。

表7 TP時(shí)間維度差異分析

從空間維度上看，朱湖農(nóng)場(chǎng)排水溝渠各斷面的TP呈現(xiàn)波動(dòng)式變化趨勢(shì)。其中監(jiān)測(cè)斷面2、6較前一監(jiān)測(cè)斷面呈現(xiàn)驟增現(xiàn)象，這與研究區(qū)內(nèi)施肥量有關(guān)。 由表8可知，監(jiān)測(cè)斷面2、3、4、5、6、7以及控制斷面8之間無顯著差異，監(jiān)測(cè)斷面3、7以及控制斷面1、8與監(jiān)測(cè)斷面2、4、5、6之間存在顯著差異。

表8 TP空間維度差異分析

2.5 灰色關(guān)聯(lián)法水質(zhì)評(píng)價(jià)

2.5.1 評(píng)價(jià)結(jié)果的確定

評(píng)價(jià)結(jié)果的確定是實(shí)測(cè)數(shù)列與標(biāo)準(zhǔn)數(shù)列相關(guān)性的排序，按照大到小的順序進(jìn)行排序，其排序在第一位的評(píng)價(jià)級(jí)別為水質(zhì)評(píng)價(jià)等級(jí)(表9)。

表9 斷面水質(zhì)評(píng)價(jià)結(jié)果表(節(jié)選)

2.5.2 基于灰色關(guān)聯(lián)分析法的農(nóng)場(chǎng)排水溝渠水質(zhì)綜合評(píng)價(jià)結(jié)果

采樣灰色關(guān)聯(lián)法和內(nèi)梅羅法對(duì)排水溝渠進(jìn)行水質(zhì)評(píng)價(jià)，評(píng)價(jià)結(jié)果見表10。

表10 朱湖農(nóng)場(chǎng)水質(zhì)評(píng)價(jià)結(jié)果

在無雨條件下，運(yùn)用灰色關(guān)聯(lián)法所得結(jié)果是控制斷面8、監(jiān)測(cè)斷面7、2以及控制斷面1均為V類水質(zhì)，監(jiān)測(cè)斷面5、6是最優(yōu)水質(zhì)斷面，均為I類水質(zhì)，運(yùn)用內(nèi)梅羅法所得結(jié)果，除了監(jiān)測(cè)斷面5、6為IV，其他斷面為V類水質(zhì)。在有雨條件下，運(yùn)用灰色關(guān)聯(lián)法所得結(jié)果是控制斷面8和控制斷面1均為V類水質(zhì)，而監(jiān)測(cè)斷面2是最優(yōu)水質(zhì)斷面，為I類水質(zhì)，運(yùn)用內(nèi)梅羅法所得結(jié)果是各斷面均為V類水質(zhì)。結(jié)果表明，在水質(zhì)綜合評(píng)價(jià)中，灰色關(guān)聯(lián)法與內(nèi)梅羅法所得結(jié)果有所差異，由于水質(zhì)綜合評(píng)價(jià)選擇指標(biāo)較少，這是形成差異的主要原因。總體上來看，灰色關(guān)聯(lián)法適應(yīng)于排水溝渠的水質(zhì)綜合評(píng)價(jià)。

3 小結(jié)與討論

本文以灰色關(guān)聯(lián)法為基礎(chǔ)，建立水質(zhì)綜合評(píng)價(jià)模型，這是一種可繼續(xù)探索的方向。傳統(tǒng)的內(nèi)梅羅法[11]在水質(zhì)綜合評(píng)價(jià)的精確度較低，在建立灰色關(guān)聯(lián)法的水質(zhì)綜合評(píng)價(jià)模型時(shí)，應(yīng)采用更加精確的水質(zhì)綜合評(píng)價(jià)法，例如模糊數(shù)學(xué)法[12]、改進(jìn)內(nèi)梅羅法等[13]。在水質(zhì)綜合評(píng)價(jià)過程中，也應(yīng)增加水質(zhì)綜合評(píng)價(jià)指標(biāo)數(shù)，指標(biāo)較少也是影響灰色關(guān)聯(lián)法精度的主要原因。

根據(jù)野外調(diào)查、實(shí)地走訪與室內(nèi)指標(biāo)監(jiān)測(cè)，得出以下結(jié)論：

1)朱湖農(nóng)場(chǎng)排水溝渠的水流流向?yàn)闁|南向西北，在其上游設(shè)有水閘，其目的是調(diào)節(jié)研究區(qū)域內(nèi)的灌溉水量。排水溝渠中水的主要用途為農(nóng)田灌溉、養(yǎng)殖魚蝦等。在研究區(qū)域內(nèi)菱角湖建有濕地公園，濕地公園水質(zhì)與排水溝渠隔開，濕地公園的水經(jīng)由處理之后，再引入濕地公園內(nèi)。排水溝渠水污染來源主要是人們排放的生產(chǎn)生活廢水，在排水溝渠下游存在水華現(xiàn)象，但其中仍有水生生物存活，部分水華區(qū)域有人為清理的痕跡。

2)從時(shí)間上來看，排水溝渠中所測(cè)定COD、NO3--N、TP等指標(biāo)的含量逐漸增加。這是降雨事件、施肥活動(dòng)、生活污水排入量等因素綜合作用的結(jié)果。而排水溝渠中NH4+-N呈現(xiàn)波動(dòng)變化趨勢(shì)，這與排水溝渠中微生物的硝化與反硝化作用有關(guān)。

3)從空間上看， 在施肥活動(dòng)、生活污水排入量以及微生物硝化與反硝化作用的共同作用下，使得排水溝渠中COD、NO3--N、TP以及NH4+-N各指標(biāo)呈現(xiàn)波動(dòng)式變化。

4)在水質(zhì)綜合評(píng)價(jià)中，灰色關(guān)聯(lián)法與內(nèi)梅羅法所得結(jié)果有所差異，由于水質(zhì)綜合評(píng)價(jià)選擇指標(biāo)較少，是形成差異的主要原因。總體上來看，灰色關(guān)聯(lián)法適應(yīng)于排水溝渠的水質(zhì)綜合評(píng)價(jià)。