李丹丹

(廣東省廣州生態(tài)環(huán)境監(jiān)測(cè)中心站，廣東 廣州 510000)

近年來(lái)，我國(guó)經(jīng)濟(jì)發(fā)展進(jìn)入了快車道。人民消費(fèi)水平得到了快速提升，生活質(zhì)量越來(lái)越高。水是生命之源，對(duì)日常生活、經(jīng)濟(jì)發(fā)展有著至關(guān)重要的作用，水質(zhì)的好壞更是影響著社會(huì)的方方面面。確保水質(zhì)達(dá)標(biāo)是政府部門的重要工作使命[1]，這就需要在水環(huán)境監(jiān)測(cè)中及時(shí)跟蹤水質(zhì)變化，及時(shí)發(fā)現(xiàn)異常監(jiān)測(cè)指標(biāo)，及早采取應(yīng)對(duì)措施，避免發(fā)生重大水質(zhì)污染問題[2]。

1 研究目的及意義

水環(huán)境監(jiān)測(cè)實(shí)踐工作中，經(jīng)常會(huì)遇到一些影響水體質(zhì)量的問題。如：地表水環(huán)境質(zhì)量達(dá)標(biāo)情況不理想；某個(gè)監(jiān)測(cè)指標(biāo)持續(xù)超標(biāo)但無(wú)法找到直接引起其污染變化的外部因素、某個(gè)時(shí)間段內(nèi)某些監(jiān)測(cè)指標(biāo)短時(shí)超標(biāo)后又能自行恢復(fù)正常等。水環(huán)境監(jiān)測(cè)指標(biāo)的變化是由多個(gè)因素綜合作用導(dǎo)致的，不僅僅可能由外部條件引起(如農(nóng)業(yè)面源污染、工業(yè)排污、生活排污等)[3]，還應(yīng)該考慮是否由其余某個(gè)監(jiān)測(cè)指標(biāo)的異常而引起。因此，了解某些監(jiān)測(cè)指標(biāo)之間的關(guān)聯(lián)性，摸清彼此間的相互影響程度，有助于在更多元化的角度和途徑去改善水質(zhì)狀況[4]。對(duì)于有效進(jìn)行水質(zhì)超標(biāo)預(yù)警，提前作出應(yīng)對(duì)措施，提高水環(huán)境治理效果具有積極意義。本文通過對(duì)流溪河流域上、中、下游斷面的數(shù)據(jù)分析發(fā)現(xiàn)，不同監(jiān)測(cè)斷面中的部分監(jiān)測(cè)指標(biāo)之間(主要為溫度、溶解氧、高錳酸鹽指數(shù)、氨氮、總磷)數(shù)值存在不同程度不同形式的關(guān)聯(lián)性。

2 監(jiān)測(cè)斷面選取原因

流溪河流域面積共2290 km2，其中2256.72 km2位于廣州市境內(nèi)(從化區(qū))。從化區(qū)位于廣州市北部，以綠色生態(tài)型、資源保護(hù)型發(fā)展方式為主[5]。流溪河流域源頭區(qū)域于從化區(qū)呂田鎮(zhèn)羅水樓村以上區(qū)域，含黃嶺水和羅水，源頭區(qū)域面積共 25.5 km2。現(xiàn)狀植被良好，人類活動(dòng)影響程度低。另外，流溪河流域由于涉及城市飲用水源保護(hù)區(qū)，故水環(huán)境治理及保護(hù)措施到位。

為了能夠更好地從監(jiān)測(cè)指標(biāo)之間的相互影響這個(gè)角度去分析數(shù)據(jù)，故本文選取了水質(zhì)較少受到人類活動(dòng)或其余外部因素影響的流溪河流域作為研究對(duì)象，采用運(yùn)行狀態(tài)穩(wěn)定且監(jiān)測(cè)頻次較高的廣州市水質(zhì)自動(dòng)監(jiān)測(cè)站的數(shù)據(jù)進(jìn)行分析，分別選取流溪河流域的上、中、下游各一個(gè)代表性斷面。

對(duì)其水質(zhì)自動(dòng)站2022年1-5月日均值監(jiān)測(cè)數(shù)據(jù)進(jìn)行相關(guān)性分析發(fā)現(xiàn)，在同一時(shí)間段內(nèi)，溶解氧濃度在同一流域不同空間上與溫度、高錳酸鹽指數(shù)、氨氮存在不同程度的關(guān)聯(lián)性。

3 溫度、溶解氧、高錳酸鹽指數(shù)、氨氮、總磷之間的關(guān)聯(lián)性

3.1 流溪河流域上游斷面

利用流溪河流域上游某斷面水質(zhì)自動(dòng)站2022年1-5月的監(jiān)測(cè)數(shù)據(jù)繪制折線圖發(fā)現(xiàn)，溫度-溶解氧曲線呈現(xiàn)高度對(duì)稱性，溫度先上升后下降，溶解氧先下降后上升，同一時(shí)間段內(nèi)兩者趨勢(shì)相反，如圖1所示，總磷-溶解氧曲線存在線性關(guān)系，在大部分之間段內(nèi)趨勢(shì)變化一致，總磷濃度上升溶解氧濃度也隨之上升，總磷濃度下降溶解氧濃度也隨之下降，如圖2所示。

圖1 流溪河流域上游某斷面水站溫度-溶解氧折線圖

圖2 流溪河流域上游某斷面水站總磷-溶解氧折線圖

利用Pearson相關(guān)系數(shù)法研究溶解氧分別和高錳酸鹽指數(shù)、氨氮、溫度、總磷之間的相關(guān)關(guān)系，結(jié)果表明溶解氧和溫度之間的相關(guān)系數(shù)值為-0.900，并且呈現(xiàn)出0.01水平的顯著性，因而說(shuō)明溶解氧和溫度之間有著顯著的負(fù)相關(guān)關(guān)系；溶解氧和總磷之間的相關(guān)系數(shù)值為0.669，并且呈現(xiàn)出0.01水平的顯著性，因而明溶解氧和總磷之間存在明顯的正相關(guān)關(guān)系，如表1所示，溶解氧和高錳酸鹽指數(shù)之間的p值大于0.05，兩者相關(guān)性不顯著；溶解氧與氨氮之間呈現(xiàn)出0.01水平的顯著性，但相關(guān)系數(shù)較小，相關(guān)性不明顯。

表1 流溪河流域上游某斷面水站

將溫度作為自變量，溶解氧作為因變量進(jìn)行線性回歸分析，一元線性回歸方程為溶解氧=13.292-0.231*溫度，R2值為0.834，意味著溫度可以解釋溶解氧的83.4%變化原因。對(duì)該一元線性回歸方程進(jìn)行F檢驗(yàn)時(shí)通過(F=661.362，p=0.000<0.05)，也即說(shuō)明溫度一定會(huì)對(duì)溶解氧產(chǎn)生影響關(guān)系。溫度的回歸系數(shù)值為-0.231(t=-25.717，p=0.000<0.01)，意味著溫度會(huì)對(duì)溶解氧產(chǎn)生顯著的負(fù)向影響關(guān)系，如表2所示。

表2 流溪河流域上游某斷面水站溫度-溶解氧線性回歸分析及F檢驗(yàn)結(jié)果

將總磷作為自變量，溶解氧作為因變量進(jìn)行線性回歸分析，一元線性回歸方程為：溶解氧=8.301+29.692×總磷，R2值為0.448，意味著總磷可以解釋溶解氧的44.8%變化原因。對(duì)該一元線性回歸方程進(jìn)行F檢驗(yàn)時(shí)通過(F=107.039， p=0.000<0.05)，也即說(shuō)明總磷一定會(huì)對(duì)溶解氧產(chǎn)生影響關(guān)系。總磷的回歸系數(shù)值為29.692(t=10.346，p=0.000<0.01)，意味著總磷會(huì)對(duì)溶解氧產(chǎn)生顯著的正向影響關(guān)系，如表3所示。

表3 流溪河流域上游某斷面水站

3.2 流溪河流域中游斷面

利用流溪河流域中游某斷面水質(zhì)自動(dòng)站2022年1-5月的監(jiān)測(cè)數(shù)據(jù)繪制折線圖發(fā)現(xiàn)，高錳酸鹽指數(shù)-總磷曲線趨勢(shì)高度吻合，總磷濃度上升，高錳酸鹽指數(shù)濃度也上升，總磷濃度下降，高錳酸鹽指數(shù)濃度也下降，如圖3所示。

圖3 流溪河流域中游某斷面水站總磷-高錳酸鹽指數(shù)折線圖

同理，利用Pearson相關(guān)系數(shù)法研究發(fā)現(xiàn)，高錳酸鹽指數(shù)和總磷之間的相關(guān)系數(shù)值為0.865，并且呈現(xiàn)出0.01水平的顯著性，因而說(shuō)明高錳酸鹽指數(shù)和總磷之間有著顯著的正相關(guān)關(guān)系，如表4所示，高錳酸鹽指數(shù)與溫度、氨氮、溶解氧之間的P 值均大于0.05，相關(guān)性不顯著。

表4 流溪河流域中游某斷面水站Pearson相關(guān)系數(shù)法分析結(jié)果

將總磷作為自變量，高錳酸鹽指數(shù)作為因變量進(jìn)行線性回歸分析，一元線性回歸方程為：高錳酸鹽指數(shù)=0.448+21.414×總磷，R2值為0.748，意味著總磷可以解釋高錳酸鹽指數(shù)的74.8%變化原因。對(duì)該一元線性回歸方程進(jìn)行F檢驗(yàn)時(shí)通過F(F=441.451，p=0.000<0.05)，也即說(shuō)明總磷一定會(huì)對(duì)高錳酸鹽指數(shù)產(chǎn)生影響關(guān)系?？偭椎幕貧w系數(shù)值為21.414(t=21.011，p=0.000<0.01)，意味著總磷會(huì)對(duì)高錳酸鹽指數(shù)產(chǎn)生顯著的正向影響關(guān)系，如表5所示。

表5 流溪河流域中游某斷面水站總磷-高錳酸鹽指數(shù)線性回歸分析及F檢驗(yàn)結(jié)果

3.3 流溪河流域下游斷面

利用李溪壩水質(zhì)自動(dòng)站2022年1-5月的監(jiān)測(cè)數(shù)據(jù)繪制折線圖發(fā)現(xiàn)，氨氮-溶解氧曲線存在線性關(guān)系，在大部分之間段內(nèi)趨勢(shì)變化相反，氨氮濃度上升，溶解氧濃度下降，氨氮濃度下降，溶解氧濃度上升，如圖4所示。

圖4 流溪河流域下游某斷面水站氨氮-溶解氧折線圖

同理，利用Pearson相關(guān)系數(shù)法研究發(fā)現(xiàn)，溶解氧和氨氮之間的相關(guān)系數(shù)值為-0.553，并且呈現(xiàn)出0.01水平的顯著性，因而說(shuō)明溶解氧和氨氮之間有著顯著的負(fù)相關(guān)關(guān)系，如表6所示，溶解氧和高錳酸鹽指數(shù)、溫度之間的P值均小于0.05，相關(guān)性不顯著；溶解氧和總磷之間的相關(guān)系數(shù)接近0，相關(guān)性較小。

將氨氮作為自變量，溶解氧作為因變量進(jìn)行線性回歸分析，一元線性回歸方程為為：溶解氧=9.247-7.216×氨氮，R2值為0.298，意味著氨氮可以解釋溶解氧的29.8%變化原因。對(duì)模型進(jìn)行F檢驗(yàn)時(shí)通過(F=63.105，p=0.000<0.05)，也即說(shuō)明氨氮一定會(huì)對(duì)溶解氧產(chǎn)生影響關(guān)系。氨氮的回歸系數(shù)值為-7.216(t=-7.944，p=0.000<0.01)，意味著氨氮會(huì)對(duì)溶解氧產(chǎn)生顯著的負(fù)向影響關(guān)系，如表7所示。

表6 流溪河流域下游某斷面水站Pearson相關(guān)系數(shù)法分析結(jié)果

表7 流溪河流域下游某斷面水站氨氮-溶解氧線性回歸分析及F檢驗(yàn)結(jié)果

4 研究結(jié)果分析與建議

4.1 分 析

通過分析流溪河流域上、中、下游監(jiān)測(cè)斷面的水質(zhì)自動(dòng)站監(jiān)測(cè)數(shù)據(jù)發(fā)現(xiàn)，其溫度、溶解氧、高錳酸鹽指數(shù)指數(shù)、氨氮、總磷之間存在關(guān)聯(lián)性。

(1)處于同一流域的不同斷面，各監(jiān)測(cè)指標(biāo)之間的關(guān)聯(lián)性大小存在差異，兩兩之間的線性關(guān)系也不完全相同，因此不能用某個(gè)斷面上某兩個(gè)監(jiān)測(cè)指標(biāo)之間的線性關(guān)系來(lái)代表整條流域的情況，更不能用來(lái)代表附近或者其他流域的情況。

(2)同一流域不同監(jiān)測(cè)斷面之間地質(zhì)條件、水文條件、生物本底條件等存在差異，污染情況也不同，因此各監(jiān)測(cè)指標(biāo)數(shù)值受到復(fù)雜因素的疊加影響，彼此之間的相互關(guān)聯(lián)性無(wú)法統(tǒng)一而論。必須在對(duì)同一流域相近斷面的大量有效數(shù)據(jù)進(jìn)行全面分析的基礎(chǔ)上，才能得出一個(gè)相對(duì)可靠的預(yù)測(cè)性信息。

4.2 建 議

(1)在進(jìn)行水質(zhì)預(yù)判預(yù)警時(shí)，可以采用合適的統(tǒng)計(jì)分析方法對(duì)同一斷面的大量監(jiān)測(cè)數(shù)據(jù)進(jìn)行演算，得出各主要監(jiān)測(cè)指標(biāo)(如溫度、溶解氧、高錳酸鹽指數(shù)、氨氮、總磷)之間的相關(guān)性，從而獲得其相對(duì)合理的線性規(guī)律，根據(jù)日常監(jiān)測(cè)的數(shù)據(jù)推算出其余相關(guān)指標(biāo)的數(shù)值，來(lái)推測(cè)出未來(lái)一段時(shí)間內(nèi)可能出現(xiàn)的變化情況，結(jié)合《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)的標(biāo)準(zhǔn)限值和各斷面水質(zhì)目標(biāo)，可預(yù)測(cè)該水體污染狀況；通過分析影響水質(zhì)的因子，找出關(guān)鍵控制因素，并制定相應(yīng)措施；對(duì)于不能達(dá)到預(yù)期要求或超出目標(biāo)值的監(jiān)測(cè)點(diǎn)，采取應(yīng)急處置預(yù)案；對(duì)可能存在超標(biāo)的監(jiān)測(cè)指標(biāo)提前關(guān)注。

(2)為提高監(jiān)測(cè)數(shù)據(jù)的合理性和可比性，在已驗(yàn)證某個(gè)斷面各監(jiān)測(cè)指標(biāo)之間的相關(guān)性后，將其按照一定的方法建立線性回歸方程，然后再用某個(gè)指標(biāo)數(shù)據(jù)，結(jié)合回歸方程計(jì)算與之相關(guān)性密切的監(jiān)測(cè)指標(biāo)的數(shù)值范圍，從而達(dá)到及時(shí)發(fā)現(xiàn)異常數(shù)據(jù)的目的，對(duì)異常數(shù)據(jù)進(jìn)行申訴，以便于使監(jiān)測(cè)數(shù)據(jù)更加科學(xué)、真實(shí)地反應(yīng)水質(zhì)狀況。

(3)在水環(huán)境監(jiān)測(cè)的實(shí)際工作中，對(duì)于某些需要重點(diǎn)關(guān)注的考核斷面，可以通過大量有效數(shù)據(jù)來(lái)分析各監(jiān)測(cè)指標(biāo)之間的關(guān)聯(lián)性，進(jìn)而得出針對(duì)該重點(diǎn)關(guān)注斷面的線性回歸模型和預(yù)測(cè)公式，從而對(duì)斷面水質(zhì)進(jìn)行預(yù)判，并結(jié)合當(dāng)前水環(huán)境監(jiān)測(cè)的具體情況提出一些改進(jìn)措施，提高水環(huán)境監(jiān)測(cè)結(jié)果的準(zhǔn)確性。

5 監(jiān)測(cè)指標(biāo)間關(guān)聯(lián)性分析帶來(lái)的水環(huán)境治理思路

(1)改變水體中氨氮、總磷濃度，間接使溶解氧濃度發(fā)生變化。在水質(zhì)監(jiān)測(cè)實(shí)際工作中可知，溶解氧較容易發(fā)生波動(dòng)，且容易受到外部不同因素的影響外，也是水質(zhì)治理中較難攻克的重要指標(biāo)。改善溶解氧超標(biāo)問題，最直接的方法是往水體中曝氣，但這一方法只能在小面積水域中操作，對(duì)于江河湖庫(kù)而言，不具有操作的可行性。經(jīng)過分析后發(fā)現(xiàn)部分?jǐn)嗝嬷邪钡?、總磷?duì)溶解氧有顯著負(fù)向影響，故我們還可以嘗試通過降低水體中氨氮、總磷的濃度，來(lái)達(dá)到一定程度上使水體溶解氧濃度上升的目的。具體操作為可以通過控制水域周邊的生活污水、工業(yè)廢水、農(nóng)業(yè)用水等排入水體，來(lái)減少水中有機(jī)物污染狀況，降低氮磷污染發(fā)生概率，從而間接提高溶解氧濃度。此外，由于溫度升高也可能引起水中溶解氧濃度降低，若在排除其他可能的因素后，溶解氧溶度依舊超標(biāo)，則可考慮是否由氣溫造成短時(shí)間內(nèi)的溶解氧超標(biāo)現(xiàn)象。

(2)遏制使總磷濃度升高的因素，一定程度上降低高錳酸鹽指數(shù)濃度。由于總磷與高錳酸鹽指數(shù)之間存在關(guān)聯(lián)性，面對(duì)部分?jǐn)嗝娓咤i酸鹽指數(shù)超標(biāo)現(xiàn)象，可考慮是否由河流兩岸及周邊農(nóng)田隨地表徑流形成的面源污染造成總磷超標(biāo)，造成水體受到大面積有機(jī)污染，造成水體中高錳酸鹽指數(shù)劇升高。因此可通過在水域中投放生物菌種加以曝氣增氧等措施，或者通過種植水生植物來(lái)吸收總磷，從而在一定程度上減少有機(jī)污染，降低高錳酸鹽指數(shù)，達(dá)到生態(tài)治理的目的。

(3)完善自動(dòng)監(jiān)測(cè)能力建設(shè)，為水環(huán)境監(jiān)測(cè)數(shù)據(jù)分析提供保障。環(huán)境監(jiān)測(cè)的質(zhì)量取決于數(shù)據(jù)測(cè)取的代表性、準(zhǔn)確性與精密性，監(jiān)測(cè)指標(biāo)之間的關(guān)聯(lián)性分析需要建立在大量的監(jiān)測(cè)數(shù)據(jù)上，考慮到人力物力的限制因素，利用自動(dòng)監(jiān)測(cè)設(shè)備獲取各監(jiān)測(cè)斷面數(shù)據(jù)是最高效的監(jiān)測(cè)方式。因此，需要提高自動(dòng)監(jiān)測(cè)能力建設(shè)，保障各類典型監(jiān)測(cè)斷面水質(zhì)自動(dòng)站數(shù)據(jù)的穩(wěn)定性、持續(xù)性和有效性，需要在自動(dòng)監(jiān)測(cè)方面加強(qiáng)能力培訓(xùn)和投入相關(guān)經(jīng)費(fèi)，建立并不斷完善自動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)體系，才能有效開展水環(huán)境質(zhì)量評(píng)價(jià)及管理，為打贏水污染治理這場(chǎng)硬仗奠定堅(jiān)實(shí)的基礎(chǔ)。

6 結(jié) 語(yǔ)

研究地表水各監(jiān)測(cè)指標(biāo)之間的相關(guān)性，對(duì)了解水質(zhì)變化趨勢(shì)具有重要意義。通過pearson相關(guān)系數(shù)法，并建立線性回歸模型分析了影響地表水主要污染因子間的相關(guān)關(guān)系及趨勢(shì)變化情況，得出相關(guān)數(shù)據(jù)的大致范圍，可作為指導(dǎo)決策的參考依據(jù)，檢查相關(guān)數(shù)據(jù)的合理性，及時(shí)發(fā)現(xiàn)異常數(shù)據(jù)，找出存在問題或原因所在，從而提出有針對(duì)性的治理對(duì)策和措施，為進(jìn)一步加強(qiáng)水環(huán)境管理提供科學(xué)依據(jù)。