石 全 ,張建華,王會(huì)容,張玉珍,顏志俊,張 森 ,沙海明
(1.六安市水利局,安徽 六安 237061;2.江蘇省水利廳,江蘇 南京 210029;3.南京水利科學(xué)研究院 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210029;4.南京市水務(wù)局,江蘇 南京 210036)
河網(wǎng)區(qū)主要指呈網(wǎng)狀結(jié)構(gòu),同時(shí)具有一定空間范圍的網(wǎng)狀水系。由于我國地形多樣,地質(zhì)構(gòu)造較為復(fù)雜,河網(wǎng)區(qū)的類型也是多種多樣,河網(wǎng)區(qū)主要是由干流及其附近的支流組成,如有湖泊與河流相通,湖泊也是該河網(wǎng)區(qū)的一部分。按照河流位置與流向、流經(jīng)區(qū)域、支流分布、上下游情況及是否受潮汐影響等因素,我國河網(wǎng)區(qū)分為感潮河網(wǎng)區(qū)和非感潮河網(wǎng)區(qū)2大類。
(1)感潮河網(wǎng)區(qū)。河網(wǎng)內(nèi)部主要河流下游與海洋交匯、受潮汐流影響的河網(wǎng)區(qū)稱為感潮河網(wǎng)區(qū)。根據(jù)感潮河網(wǎng)區(qū)內(nèi)部主要河流規(guī)模、水文特征及其附近支流的密集程度等因素,感潮河網(wǎng)區(qū)又分為大型感潮河網(wǎng)區(qū)和一般感潮河網(wǎng)區(qū)2類。我國大型感潮河網(wǎng)區(qū)主要有太湖流域河網(wǎng)區(qū)和珠江三角洲河網(wǎng)區(qū);一般感潮河網(wǎng)區(qū)主要有閩江河網(wǎng)區(qū)、溫黃平原河網(wǎng)區(qū)、韓江三角洲河網(wǎng)區(qū)等。
(2)非感潮河網(wǎng)區(qū)。非感潮河網(wǎng)區(qū)指不受潮汐流影響,處于江河中上游地區(qū),一般有湖泊與河流相通的河網(wǎng)區(qū)。非感潮河網(wǎng)區(qū)由一個(gè)或若干個(gè)樹狀河網(wǎng)組成,對(duì)于一般樹狀分布河網(wǎng),上一級(jí)河道總是匯入至下一級(jí)河道,位于河道下游斷面的流量總是由其上游匯集而至。我國非感潮河網(wǎng)主要有鄱陽湖河網(wǎng)區(qū)、洞庭湖河網(wǎng)區(qū)、成都平原河網(wǎng)區(qū)等。
(1)河網(wǎng)自然特征。河網(wǎng)地區(qū)一般地勢(shì)平坦,河網(wǎng)密布,縱橫交錯(cuò),河漫灘發(fā)育,河床較淺,低坡平緩。網(wǎng)狀河流結(jié)構(gòu)是河網(wǎng)地區(qū)的典型地理特征,網(wǎng)狀河流具有比降較小、相互通連、河道穩(wěn)定及具有非常小的坡降和河道寬深比的基本特點(diǎn)。我國東部平原河網(wǎng)地區(qū),大部分地區(qū)受到潮汐的影響,受潮汐影響的區(qū)域,流向、流態(tài)呈現(xiàn)隨機(jī)性變化,易使污染物在河道往復(fù)傳播,充分混合,從而影響整片區(qū)域。
(2)水文水環(huán)境特征。河網(wǎng)區(qū)一般水量豐富、流速平緩、流向不定、水質(zhì)較差。河流流速平緩,水動(dòng)力條件不足,水體自凈能力較差,易發(fā)生水環(huán)境污染。河網(wǎng)地區(qū)人口密集,工農(nóng)業(yè)發(fā)達(dá),產(chǎn)生的污水也多。天然的地貌造就了河流縱橫交錯(cuò),水流相互通連,在潮汐作用下,水流往復(fù)運(yùn)動(dòng),污染物極易傳播,蜿蜒曲折,河道比降小,水流動(dòng)力嚴(yán)重不足,污染物難以迅速排出,容易累積水體遭受污染[1]。
(3)人工干預(yù)和調(diào)節(jié)影響大。河網(wǎng)區(qū)河流多受人工節(jié)制,水閘工程分隔、調(diào)控便利、流動(dòng)受限。高度的人工控制性,給環(huán)境流量的調(diào)節(jié)帶來了正負(fù)兩面的影響,正面影響是河流水位抬高、槽蓄水量增加等;負(fù)面影響是由于水閘的分隔,河流的連續(xù)性受到破壞,降低了水體的流動(dòng)性,自凈能力降低等。
(1)因果分析圖法。因果分析圖是逐層深入地分析問題產(chǎn)生原因的有效根據(jù),由質(zhì)量特征、要因、主干、支干等組成,主要通過對(duì)重點(diǎn)影響因素的篩選,結(jié)合各因素的特征加以分析,分類后進(jìn)行整理歸納,使之形成層次分明、條理清楚、一目了然的圖形。河網(wǎng)區(qū)水污染排放主要問題的原因,可能出現(xiàn)某種問題或呈現(xiàn)某種狀態(tài)的原因多種多樣,通過采用因果分析圖法,對(duì)河網(wǎng)區(qū)水污染排放可能影響因素或存在問題,采用圖示的方法,逐層深入排查可能原因,確定主要影響因素,進(jìn)而針對(duì)性地提出對(duì)策和措施,達(dá)到水污染物的有效控制,改善河網(wǎng)區(qū)水質(zhì)。
(2)河網(wǎng)區(qū)水污染控影響因果分析圖繪制。明確河網(wǎng)區(qū)水污染排放所要解決的問題和結(jié)果。首先由左至右畫出一條水平主干線,箭頭指向一個(gè)矩形框,框內(nèi)注明研究的問題,即結(jié)果;根據(jù)影響河網(wǎng)區(qū)水污染排放的各種因素,分析確定影響較大的方面原因;將每種大原因進(jìn)一步分解為中原因、小原因,直至分解的原因可以采取具體措施加以解決為止;檢查圖中的所列要素和原因是否齊全,對(duì)初步分析結(jié)果廣泛征求意見,進(jìn)行必要的補(bǔ)充和修改;選擇影響權(quán)重大的關(guān)鍵因素,以便重點(diǎn)采取措施加以控制。用因果分析圖來歸納影響河網(wǎng)區(qū)水污染排放要素和問題的因果關(guān)系。因果分析識(shí)別河網(wǎng)區(qū)水污染控制主要影響因素過程,見圖1。
河網(wǎng)區(qū)水污染控制的關(guān)鍵影響因素,是對(duì)區(qū)域水資源特征、開發(fā)利用、節(jié)約保護(hù)和污染排放等諸方面綜合分析的基礎(chǔ)上,經(jīng)過多要素分析而得出的,故需要對(duì)河網(wǎng)區(qū)水污染控制指標(biāo)體系中的影響因素作全面定量分析,計(jì)算因素的綜合影響貢獻(xiàn)程度,確定能綜合反映水資源、社會(huì)、經(jīng)濟(jì)、環(huán)境相互依存、相互協(xié)調(diào)關(guān)系的指標(biāo),使整體結(jié)構(gòu)更為簡(jiǎn)單。系統(tǒng)控制關(guān)鍵影響因素指標(biāo)定量診斷分析方法,主要有主成分分析法[2]等。
2.2.1 主成分分析法基本原理
主成分分析是一種模式識(shí)別中的降維映射方法,主要是將多維空間的信息在低維(二維或三維)空間表現(xiàn)出來,消除眾多信息相互重疊的部分。將原始變量進(jìn)行轉(zhuǎn)換,通過原始變量指標(biāo)的線性組合,優(yōu)化組合系數(shù),使新的變量指標(biāo)之間相互獨(dú)立且代表性好。
如在一個(gè)指標(biāo)體系中,有n個(gè)樣本m個(gè)變量,形成指標(biāo)矩陣X=[xij](i=1,2,...,n,j=1,2,...,m)。
圖1 因果分析法識(shí)別河網(wǎng)區(qū)水污染控制影響因素過程
它們的綜合指標(biāo)為z1,z2,...,zp(p≤ m),則:
式中:zi和zj互不相關(guān),z1為x1的一切線性組合的方差最大者,并且彼此互不相關(guān)。這些新變量z1,z2,...,zp就是原來變量的第一、第二、…、第p主成分。其中zi在總方差中占的比例最大,依次遞減。
2.2.2 主要分析步驟
選取指標(biāo)及指標(biāo)的同趨化處理。通過以上指標(biāo)選取,對(duì)與評(píng)價(jià)目的呈負(fù)相關(guān)的指標(biāo)進(jìn)行同趨化處理。用1/x或1-x代替原指標(biāo),保證所有指標(biāo)得分都對(duì)評(píng)價(jià)目標(biāo)作正貢獻(xiàn)。
(1)采集數(shù)據(jù)樣本。設(shè)樣本數(shù)據(jù)為n,對(duì)每個(gè)樣本收集到的指標(biāo)數(shù)據(jù)進(jìn)行趨同化處理后可得到矩陣為:
式中:Yij代表第i個(gè)樣本第j項(xiàng)評(píng)價(jià)指標(biāo)數(shù)值。
(2)指標(biāo)標(biāo)準(zhǔn)化處理。由于各指標(biāo)數(shù)據(jù)量綱都不一致,數(shù)量間的差異很大,需要將不同度量的指標(biāo)轉(zhuǎn)化為同度量的指標(biāo)。使各指標(biāo)具有可比性。標(biāo)準(zhǔn)化公式如下:
式中(i=1,2,...,n;j=1,2,...,m)。
(3)計(jì)算相關(guān)矩陣。將處理后的指標(biāo)進(jìn)行相關(guān)性分析,設(shè)rij為原來變量xi與xj的相關(guān)系數(shù),其計(jì)算公式為:
可得相關(guān)系數(shù)矩陣:
整理得到:R=XXT(/n-1),為實(shí)對(duì)稱矩陣(即rij= rji)。
(4)計(jì)算特征值和特征向量。利用上述相關(guān)系數(shù)矩陣R,求特征方程|λI-R|=0的m個(gè)特征值 λi(i=1,2,...,m),按其大小順序排列,即λ1≥λ2≥…≥λm≥0,然后分別求出特征向量e(ii=1,2,...,m)。
(5)計(jì)算主成分貢獻(xiàn)率及累計(jì)貢獻(xiàn)率。主成分貢獻(xiàn)率見式(6),累計(jì)貢獻(xiàn)率見式(7)。
(6)選擇主成分并計(jì)算主成分載荷。設(shè)定對(duì)主成分所包含的總體信息程度,即累計(jì)貢獻(xiàn)率。一般若前p個(gè)主成分提供了85%~95%的信息量,則取特征值 λ1,λ2,…,λp。所對(duì)應(yīng)的第一、第二、…、第p(p≤m)個(gè)主成分,變量數(shù)由m減至p個(gè),產(chǎn)生了降維的效果。
主成分載荷計(jì)算式為:
可進(jìn)一步計(jì)算主成分得分:
由此可對(duì)影響因素歸類,按照主成分代表的類別進(jìn)行系統(tǒng)診斷分析。
太湖流域?yàn)槲覈笮透谐焙泳W(wǎng)區(qū),涉及江蘇省、浙江省和上海市,總面積36895 km2,是我國經(jīng)濟(jì)社會(huì)發(fā)展最具活力的地區(qū)之一。2015年人口5997萬人,占全國總?cè)丝诘?.4%,GDP為6.69萬億元,占全國總GDP的9.9%,人均GDP11.2萬元,是全國人均GDP的2.3倍。
2015年太湖流域全年水質(zhì)達(dá)到或優(yōu)于Ⅲ類的比例僅為20.3%,主要為氨氮、總磷、高錳酸鹽指數(shù)、生化需氧量、化學(xué)需氧量和溶解氧等超標(biāo)[3]。
針對(duì)太湖流域河網(wǎng)區(qū)社會(huì)經(jīng)濟(jì)和水污染情況,利用因果分析圖法,河網(wǎng)區(qū)26個(gè)水污染控制的影響因素分析數(shù)據(jù),見表1[3]。
表1 太湖流域河網(wǎng)區(qū)水污染控制影響因素和分析數(shù)據(jù)
(續(xù)表1)
利用主成分分析方法,計(jì)算太湖流域河網(wǎng)區(qū)水污染控制影響因素特征值及主成分貢獻(xiàn)率,見表2。
由表2可知,太湖流域河網(wǎng)區(qū)水污染控制影響因素第一、第二、第三主成分的累計(jì)貢獻(xiàn)率已達(dá)到85.56%,達(dá)到主成分分析規(guī)定要求。根據(jù)主成分載荷計(jì)算公式得到主成分載荷矩陣,見表3。
表2 特征值及主成分貢獻(xiàn)率
表3 主成分載荷矩陣
(續(xù)表3)
根據(jù)表2的污染控制影響因素主成分貢獻(xiàn)率,結(jié)合表3的主成分載荷值,可以看出:
(1)第一主成分:對(duì)污染控制影響的貢獻(xiàn)率為65.85%,其中萬元GDP用水量(X1)、農(nóng)田畝均灌溉用水量(X4)、人均用水量(X6)、地表水源供水量(X11)、地下水源供水量(X12)、長(zhǎng)江引水總量(X13)、錢塘江引水總量(X14)、城鎮(zhèn)居民生活污水排放量(X15)、第二產(chǎn)業(yè)廢污水排放量(X16)、第三產(chǎn)業(yè)廢污水排放量(X17)、廢污水排放總量(X22)、城鎮(zhèn)廢污水處理率(X25)、農(nóng)田畝均化肥施用量(X26)等13個(gè)影響因素的載荷系數(shù)在0.9以上,與第一主成分都有較強(qiáng)的正相關(guān)或負(fù)相關(guān),表明這些因素對(duì)太湖流域河網(wǎng)區(qū)水污染控制具有很強(qiáng)的影響效力。
(2)第二主成分:對(duì)污染控制影響的貢獻(xiàn)率為12.57%,其中農(nóng)田灌溉水有效利用系數(shù)(X2)、工業(yè)用水重復(fù)利用率(X3)、GDP增長(zhǎng)率(X7)、人均GDP(X8)、人口(X9)、評(píng)價(jià)河長(zhǎng)水質(zhì)達(dá)到或優(yōu)于Ⅲ類(X19)、省界河流參評(píng)斷面水質(zhì)達(dá)到或優(yōu)于Ⅲ類(X20)、區(qū)水全年期水質(zhì)達(dá)標(biāo)率(X21)、農(nóng)業(yè)節(jié)水灌溉率(X23)、城鎮(zhèn)廢污水處理回用率(X25)等10個(gè)影響因素的載荷系數(shù)為0.8~0.9,與第二主成分都有較強(qiáng)的正相關(guān)或負(fù)相關(guān),表明這些因素對(duì)太湖流域河網(wǎng)區(qū)水污染控制具有一定的影響效力。
(3)第三主成分:指標(biāo)影響貢獻(xiàn)率僅為7.13%左右,其中人均水資源量(X5)、降雨量(X10)和耗水率(X2)等3個(gè)因素的的載荷系數(shù)為0.7~0.8,表明這些因素對(duì)太湖流域河網(wǎng)區(qū)水污染控制影響效力相對(duì)較小。
上述分析結(jié)果表明,供用水量、廢污水排放量等水資源開發(fā)利用和經(jīng)濟(jì)社會(huì)發(fā)展等因素,是太湖流域河網(wǎng)區(qū)水污染控制的關(guān)鍵影響因素。
河網(wǎng)地區(qū)是我國水資源比較豐富的地區(qū),密布的河網(wǎng)為社會(huì)經(jīng)濟(jì)的快速發(fā)展提供了寶貴資源。河網(wǎng)區(qū)人口密集、經(jīng)濟(jì)發(fā)達(dá),大量工農(nóng)業(yè)廢污水和生活污水直接或間接排放到河網(wǎng)中,對(duì)水體水質(zhì)和水生態(tài)環(huán)境構(gòu)成嚴(yán)重威脅[4]。河網(wǎng)區(qū)水污染控制的主要措施為:
(1)地表水污染控制。加強(qiáng)河網(wǎng)地區(qū)工業(yè)污染控制,做到增產(chǎn)減污或增產(chǎn)不增污;加大城市污水處理設(shè)施建設(shè)力度,提高污水處理水平;加強(qiáng)面源污染控制,減少氮、磷排放量;實(shí)施排污口整治工程,優(yōu)化入河排污口設(shè)置;加強(qiáng)監(jiān)測(cè)預(yù)警,提高應(yīng)急處置能力。
(2)地下水污染控制。調(diào)整河網(wǎng)地區(qū)供用水格局,將優(yōu)質(zhì)有限的深層地下水主要用于城鄉(xiāng)居民的生活供水;劃定地下水保護(hù)區(qū),設(shè)置衛(wèi)生防護(hù)帶;完善地下水資源管理制度,加強(qiáng)對(duì)地下水取水的監(jiān)督、檢查;完善地下水監(jiān)測(cè)網(wǎng)絡(luò)、動(dòng)態(tài)管理信息和預(yù)警系統(tǒng)。
(3)水生態(tài)保護(hù)和修復(fù)。合理配置河網(wǎng)地區(qū)水資源,優(yōu)化水利工程調(diào)度,建設(shè)生態(tài)用水保護(hù)工程,保障河流生態(tài)基流量和湖庫生態(tài)水位;建設(shè)河湖濕地生態(tài)工程,培育水生植被,恢復(fù)湖泊、河流水生態(tài)系統(tǒng)功能;健全生態(tài)監(jiān)測(cè)與預(yù)警網(wǎng)絡(luò)體系,維護(hù)河湖水生態(tài)系統(tǒng)良性循環(huán)。
參考文獻(xiàn):
[1]YU Miao,JIN Tong. Principal Component Analysis of Major Pollutants Discharge Amount in Major Cities,Agricultural Science&Technology[J] . 2017,18(7)∶1260-1262 .
[2]李玉珍,王宜懷.主成分分析及算法[J] . 蘇州大學(xué)學(xué)報(bào)(自然科學(xué)版),2005(01)∶32-36 .
[3]水利部太湖流域管理局 . 太湖流域及東南諸河水資源公報(bào)[R] . 2003-2016 .
[4]鄭石,王熒.基于主成分分析的農(nóng)業(yè)經(jīng)濟(jì)與環(huán)境污染問題研究 [J] . 江西農(nóng)業(yè)學(xué)報(bào),2017,29(6)∶125-130 .