“引江濟(jì)太”調(diào)水中望虞河水質(zhì)變化的規(guī)律

張 又，劉 凌，姚秀嵐，麻 林，趙 溦

(河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇南京 210098)

太湖流域具有獨(dú)特的地形地貌特征，在徑流和潮汐雙重作用下，人類活動(dòng)的強(qiáng)烈影響使之易形成重污染區(qū)。太湖流域水系與長(zhǎng)江、錢塘江相連，水利工程眾多，水體流動(dòng)性差，流域內(nèi)河網(wǎng)水質(zhì)普遍超標(biāo)，水污染問(wèn)題突出，湖泊富營(yíng)養(yǎng)化嚴(yán)重［1-2］，迫切需求改善太湖流域水系的連通性，以保障太湖流域的生態(tài)環(huán)境，實(shí)現(xiàn)太流流域水資源的可持續(xù)利用［3］。

望虞河位于太湖流域陽(yáng)澄淀泖區(qū)和武澄錫虞區(qū)的交界處，連接長(zhǎng)江與太湖，是太湖流域“引江濟(jì)太”的重要調(diào)水通道［4］。本研究旨在通過(guò)建立馬爾科夫模型來(lái)分析2008 年“引江濟(jì)太”調(diào)水對(duì)望虞河水質(zhì)的影響，探索江河湖連通性的改善對(duì)望虞河水質(zhì)影響的規(guī)律。

1 調(diào)水試驗(yàn)概況

2007 年5 月，太湖藍(lán)藻暴發(fā)導(dǎo)致無(wú)錫市供水危機(jī)，“引江濟(jì)太”應(yīng)急調(diào)度的實(shí)施，對(duì)改善無(wú)錫太湖水源地水質(zhì)起到了積極作用。為降低2008 年太湖藍(lán)藻大規(guī)模暴發(fā)的可能性，保障太湖流域供水安全，“引江濟(jì)太”調(diào)水于2008 年1 月提前啟動(dòng)，以進(jìn)一步探索冬春季引水對(duì)減少太湖藍(lán)藻暴發(fā)的作用。

2008 年1 月10 日開(kāi)啟常熟樞紐泵站，實(shí)施“引江濟(jì)太”;1 月22 日望亭水利樞紐開(kāi)閘引水入湖。由于入湖水質(zhì)受到影響，3 月22—25 日、3 月28 日至4 月2 日望亭水利樞紐關(guān)閘，待水質(zhì)改善后繼續(xù)開(kāi)閘引水。6 月9 日望亭水利樞紐關(guān)閘，引水結(jié)束。

此次調(diào)水試驗(yàn)，望虞河上共設(shè)置江邊閘下、虞義大橋、張橋、大橋角新橋和望亭閘下5 個(gè)監(jiān)測(cè)斷面，主要監(jiān)測(cè)指標(biāo)有DO、CODMn、NH3-N、TP 和TN。

2 馬爾科夫模型

馬爾科夫模型［5］是研究某一事件的狀態(tài)及狀態(tài)之間轉(zhuǎn)移規(guī)律的模型，它通過(guò)對(duì)t0時(shí)刻事件不同狀態(tài)的初始概率及狀態(tài)之間的轉(zhuǎn)移關(guān)系來(lái)研究t0+t 時(shí)刻狀態(tài)的變化趨勢(shì)。馬爾科夫模型具有無(wú)后效性，即狀態(tài)轉(zhuǎn)移概率僅與轉(zhuǎn)移出發(fā)狀態(tài)i、轉(zhuǎn)移步數(shù)k、轉(zhuǎn)移后狀態(tài)j 有關(guān)，而與轉(zhuǎn)移的初始時(shí)刻無(wú)關(guān)［6］。

馬爾科夫模型具有描述事物隨機(jī)變化的良好特性，已大量應(yīng)用于環(huán)境評(píng)價(jià)與預(yù)測(cè)方面。鐘政林等［7］將馬爾科夫模型應(yīng)用于河流綜合水質(zhì)預(yù)報(bào);馬占青等［8］應(yīng)用灰色馬爾科夫模型［9］預(yù)測(cè)城市污水排放量;何斌等［10］用馬爾科夫模型求得水質(zhì)不同級(jí)別之間的轉(zhuǎn)移概率，在此基礎(chǔ)上建立水質(zhì)動(dòng)態(tài)評(píng)價(jià)模型，評(píng)價(jià)水質(zhì)的動(dòng)態(tài)變化情況;馮為等［11］建立動(dòng)態(tài)的馬爾科夫水質(zhì)評(píng)價(jià)模型，對(duì)轉(zhuǎn)移概率賦權(quán)并計(jì)算絕對(duì)進(jìn)步度，在此基礎(chǔ)上引出相對(duì)進(jìn)步度的概念;馬惠群等［12］將基于Kendall τ 的指標(biāo)權(quán)重引入馬爾科夫模型，并應(yīng)用于大明湖的水質(zhì)評(píng)價(jià)。

筆者將模糊隸屬度與馬爾科夫模型相結(jié)合，應(yīng)用于江河湖連通條件下調(diào)水區(qū)水質(zhì)改善度的計(jì)算。通過(guò)計(jì)算望虞河調(diào)水試驗(yàn)中調(diào)水區(qū)水質(zhì)改善狀況的總進(jìn)步度，以及5 個(gè)斷面的水質(zhì)變化進(jìn)步度，以分析江河湖連通條件下調(diào)水區(qū)水體污染指標(biāo)質(zhì)量改善的規(guī)律。

2.1 水質(zhì)劃分

最常見(jiàn)的水質(zhì)劃分方法為:根據(jù)水環(huán)境質(zhì)量標(biāo)準(zhǔn)中的水質(zhì)類別，將水質(zhì)指標(biāo)值劃分為3 ～6 個(gè)狀態(tài)。這種方法無(wú)需將數(shù)據(jù)規(guī)格化，直接根據(jù)原始指標(biāo)值劃分等級(jí)狀態(tài)［7，10-11］。但是，當(dāng)某一指標(biāo)的檢測(cè)值大量集中在同一個(gè)水質(zhì)類別時(shí)，本方法達(dá)不到判斷水質(zhì)變化趨勢(shì)的目的。因此，需要重新規(guī)定隸屬范圍更小、更精細(xì)的狀態(tài)劃分標(biāo)準(zhǔn)。筆者綜合考慮各指標(biāo)監(jiān)測(cè)值的實(shí)際情況，對(duì)各指標(biāo)進(jìn)行規(guī)格化處理，將規(guī)格化后的指標(biāo)值劃分為11 個(gè)等級(jí)(即11 種狀態(tài))，狀態(tài)空間表示為E={1，2，…，11}，見(jiàn)表1。

表1 馬爾科夫水質(zhì)狀態(tài)劃分情況

2.2 指標(biāo)規(guī)格化

不同指標(biāo)的屬性值量綱是不同的，為消除量綱帶來(lái)的影響，需要對(duì)指標(biāo)值進(jìn)行規(guī)格化。規(guī)格化就是將一個(gè)屬性取值范圍投射到一個(gè)特定范圍之內(nèi)，以消除因數(shù)值大小不一而造成評(píng)價(jià)結(jié)果的偏差［13］。指標(biāo)的規(guī)格化一般選用Zadeh 公式:

越大越優(yōu)型

越小越優(yōu)型

式中:rij為第j 個(gè)樣本的第i 個(gè)指標(biāo)的標(biāo)準(zhǔn)化數(shù)值;xij為第j 個(gè)樣本的第i 個(gè)指標(biāo)的實(shí)測(cè)值;ximax為第i個(gè)指標(biāo)實(shí)測(cè)值的最大值;ximin為第i 個(gè)指標(biāo)實(shí)測(cè)值的最小值。

在評(píng)價(jià)指標(biāo)中，各指標(biāo)實(shí)測(cè)值變化范圍大小不一，傳統(tǒng)的Zadeh 公式使所有指標(biāo)實(shí)測(cè)值的最大、最小值均規(guī)格化為1 與0 或0 與1，夸大了實(shí)測(cè)值變化范圍較小的指標(biāo)在評(píng)價(jià)中的相對(duì)作用，同時(shí)削弱了實(shí)測(cè)值變化范圍較大的指標(biāo)在評(píng)價(jià)中的相對(duì)作用，使評(píng)價(jià)結(jié)果偏離實(shí)情［14］。

本文采用以下改進(jìn)的Zadeh 公式進(jìn)行計(jì)算:越大越優(yōu)型

越小越優(yōu)型

式中:x'imax為第i 個(gè)指標(biāo)的理想最大值;x'imin為第i 個(gè)指標(biāo)的理想最小值。

當(dāng)指標(biāo)特征值變化范圍較小時(shí)，指標(biāo)最大實(shí)測(cè)值與最小實(shí)測(cè)值相差不大，計(jì)算所得的指標(biāo)標(biāo)準(zhǔn)化數(shù)值也相差不大，反之亦然，因此能較好地反映指標(biāo)特征值變化的實(shí)際情況。同時(shí)，各指標(biāo)的最大實(shí)測(cè)值ximax與最小實(shí)測(cè)值ximin用其理想最大值x'imax和理想最小值x'imin代替，各指標(biāo)的理想最大值和理想最小值可根據(jù)樣本的實(shí)際情況選定。該方法可以有效避免異常的實(shí)測(cè)最大值或最小值對(duì)評(píng)價(jià)結(jié)果的影響。

根據(jù)表1 中劃分的質(zhì)量狀態(tài)，找到每個(gè)指標(biāo)標(biāo)準(zhǔn)化數(shù)值rij隸屬的狀態(tài)bij，從而得到指標(biāo)隸屬矩陣:

2.3 概率轉(zhuǎn)移矩陣的構(gòu)造

假設(shè)對(duì)m 個(gè)指標(biāo)進(jìn)行了n 次監(jiān)測(cè)，在這n 次監(jiān)測(cè)中相鄰的兩次監(jiān)測(cè)值由i 級(jí)轉(zhuǎn)變成j 級(jí)的總數(shù)為kij(i，j=1，2，…，11)，并且認(rèn)為第n 次監(jiān)測(cè)值轉(zhuǎn)變?yōu)楸旧?，i 級(jí)內(nèi)的原始樣本數(shù)為Ki，則i 級(jí)到j(luò) 級(jí)的轉(zhuǎn)移概率為

由此得到概率轉(zhuǎn)移矩陣:

2.4 進(jìn)步度計(jì)算

概率轉(zhuǎn)移矩陣P 已經(jīng)可以大致反映監(jiān)測(cè)期內(nèi)水質(zhì)的動(dòng)態(tài)發(fā)展情況［11］。為了避免馬爾科夫鏈發(fā)生多步轉(zhuǎn)移，以及滿足極限狀態(tài)下才能求得穩(wěn)定概率分布這一苛刻條件［15］，在一步轉(zhuǎn)移矩陣后，利用進(jìn)步度來(lái)動(dòng)態(tài)評(píng)估水質(zhì)變化情況。

進(jìn)步度sij有幾種不同表示方法，本研究選用以下方法計(jì)算進(jìn)步度sij［11］:

稱S=(sij)n×n為概率轉(zhuǎn)移矩陣P 的進(jìn)步矩陣，則絕對(duì)進(jìn)步度為

將PPD(s)進(jìn)行規(guī)格化，得到相對(duì)進(jìn)步度:

式中，N 即為所分狀態(tài)的個(gè)數(shù)，本文中N =11，即Ppd(s)=PPD(s)/100。

相對(duì)進(jìn)步度Ppd(s)是一個(gè)范圍在［-1，1］之間的數(shù)，可以直接反映水質(zhì)是變好還是變壞。若Ppd(s)＞0，說(shuō)明污染物指標(biāo)質(zhì)量是改善的，反之則是惡化的。水質(zhì)改善度最大時(shí)的相對(duì)進(jìn)步度為1，水質(zhì)惡化最嚴(yán)重時(shí)的相對(duì)進(jìn)步度為-1，相對(duì)進(jìn)步度的絕對(duì)值越大說(shuō)明改善或惡化的程度越大。由相對(duì)進(jìn)步度可以很直觀地看出單指標(biāo)或綜合指標(biāo)的變化情況，而且進(jìn)步度是一個(gè)概率統(tǒng)計(jì)參數(shù)，可以宏觀定量分析水質(zhì)與水量的關(guān)系，不受個(gè)別突變值的影響。

3 模型建立

根據(jù)太湖流域管理局提供的DO、CODMn、NH3-N、TP 和TN 這5 個(gè)指標(biāo)的檢測(cè)數(shù)據(jù)，所有斷面監(jiān)測(cè)到的DO 值均可達(dá)到Ⅱ類水質(zhì)及以上，評(píng)價(jià)時(shí)可不予考慮，故選擇CODMn(x1)、NH3-N (x2)、TP (x3)和TN (x4)這4 個(gè)指標(biāo)計(jì)算水質(zhì)改善情況的進(jìn)步度。

根據(jù)各指標(biāo)值的最大、最小值，結(jié)合水質(zhì)等級(jí)劃分標(biāo)準(zhǔn)［16］，選取各指標(biāo)的理想最大值和理想最小值，見(jiàn)表2。

表2 各指標(biāo)理想最大、最小值 mg/L

根據(jù)各監(jiān)測(cè)點(diǎn)代表河長(zhǎng)(江邊閘下7.8 km;虞義大橋13.4 km;張橋15.4 km;大橋角新橋16 km;望亭閘下7.4 km)，將各監(jiān)測(cè)點(diǎn)的CODMn、NH3-N、TP和TN 這4 個(gè)指標(biāo)的值加權(quán)平均，可計(jì)算得到望虞河全河1—6 月各指標(biāo)的綜合值。

根據(jù)式(3)、式(4)對(duì)望虞河各指標(biāo)的綜合數(shù)值以及各監(jiān)測(cè)點(diǎn)的指標(biāo)值進(jìn)行規(guī)格化處理，結(jié)果見(jiàn)表3 和表4。

表3 望虞河各指標(biāo)規(guī)格化結(jié)果

表4 各監(jiān)測(cè)點(diǎn)各指標(biāo)規(guī)格化結(jié)果

根據(jù)表1 所劃分的質(zhì)量狀態(tài)，將表3 和表4 中的數(shù)據(jù)按照式(5)計(jì)算，得到隸屬矩陣為

根據(jù)式(6)～式(8)，運(yùn)用MATLAB 軟件編程，將隸屬矩陣直接代入計(jì)算得到相對(duì)進(jìn)步度。

經(jīng)計(jì)算得到望虞河調(diào)水試驗(yàn)調(diào)水區(qū)1—6 月的水質(zhì)變化進(jìn)步度為:Ppd(s)=0.006 2;江邊閘下、虞義大橋、張橋、大橋角新橋及望亭閘下5 個(gè)監(jiān)測(cè)斷面1—6 月的水質(zhì)變化進(jìn)步度分別為:Ppd(s)1=0.057 5;Ppd(s)2= -0.014 2;Ppd(s)3=0.005 7;Ppd(s)4=0.023 6;Ppd(s)5=0.067 6。

4 計(jì)算結(jié)果分析

編程計(jì)算所得的水質(zhì)變化進(jìn)步度為水質(zhì)變化的相對(duì)進(jìn)步度Ppd(s)，其值可以直接反映水質(zhì)是變好還是變壞。

由計(jì)算結(jié)果可以看出，望虞河水質(zhì)在6 個(gè)月內(nèi)總體上是改善的，其進(jìn)步度為0.006 2，說(shuō)明其改善程度為0.006 2。通過(guò)分析計(jì)算結(jié)果可知，望虞河水質(zhì)從1—6 月總的水質(zhì)變化進(jìn)步度為0.037 2。從調(diào)水上游至下游，各監(jiān)測(cè)斷面1—6 月總的水質(zhì)變化進(jìn)步度分別為0.345 0、-0.085 2、0.034 2、0.141 6 和0.405 6(圖1)，即江邊閘下、張橋、大橋角新橋及望亭閘下這4 個(gè)監(jiān)測(cè)點(diǎn)的水質(zhì)改善程度分別為0.345 0、0.034 2、0.141 6 和0.405 6，而虞義大橋監(jiān)測(cè)點(diǎn)的惡化程度為0.085 2。

圖1 各監(jiān)測(cè)斷面水質(zhì)變化進(jìn)步度變化趨勢(shì)

每個(gè)斷面的水質(zhì)變化進(jìn)步度反映的是相應(yīng)斷面水質(zhì)隨時(shí)間的變化情況，即該點(diǎn)水質(zhì)較之前的變化程度，而不是沿程的變化情況。從圖1 可以看出，此次調(diào)水期間，望虞河河水在江邊閘下斷面處的水質(zhì)改善明顯，虞義大橋水質(zhì)效果最差，從張橋斷面以后水質(zhì)便開(kāi)始逐漸改善。

從時(shí)間跨度上分析水質(zhì)改善的原因。首先，1—6 月，由冬及夏，水體的自凈能力隨溫度的提升而有所增強(qiáng)，水質(zhì)因此在一定程度上得到改善;其次，這6 個(gè)月中長(zhǎng)江水量呈增長(zhǎng)趨勢(shì)，水質(zhì)隨水量增多而有所改善，作為望虞河調(diào)水的水源，其水質(zhì)的改善直接影響望虞河調(diào)水通道的第一個(gè)監(jiān)測(cè)點(diǎn)，因此江邊閘下的水質(zhì)改善明顯;第三，在調(diào)水前半程，江水流經(jīng)望虞河，一定程度上替換了原來(lái)的水體，在調(diào)水后半程，雖然調(diào)水對(duì)望虞河河水的替換效果減弱，但污染物沿程的沉積和擴(kuò)散也成為水質(zhì)改善的原因之一。由于這些現(xiàn)象多與調(diào)水水源和流域匯流的水量有關(guān)，故在夏季水量充沛時(shí)水質(zhì)改善效果更為顯著。

從沿程匯流的角度分析各監(jiān)測(cè)斷面水質(zhì)改善程度差別的原因。望虞河西岸支流污染較嚴(yán)重，以錫北運(yùn)河和九里河為甚，其水質(zhì)超標(biāo)達(dá)4 倍，其他西岸支流也不同程度的水質(zhì)超標(biāo)。錫北運(yùn)河和九里河的匯入，直接影響到虞義大橋和張橋兩個(gè)監(jiān)測(cè)斷面的水質(zhì)，因此，這兩個(gè)斷面的水質(zhì)改善情況較差，虞義大橋甚至呈水質(zhì)惡化的趨勢(shì)。由于河水自身的凈化作用等原因，流動(dòng)中的水質(zhì)逐漸好轉(zhuǎn)，至大橋角新橋斷面水質(zhì)已得到明顯改善。

5 結(jié) 論

a. 運(yùn)用馬爾科夫模型，得到調(diào)水后望虞河水質(zhì)改善情況的進(jìn)步度，結(jié)果顯示水質(zhì)從1 月到6 月總的改善程度為0.037 2。

b. 計(jì)算各斷面水質(zhì)變化進(jìn)步度，得到望虞河水質(zhì)的變化規(guī)律:望虞河江邊閘下斷面處的水質(zhì)改善明顯，虞義大橋水質(zhì)呈惡化趨勢(shì)，從張橋斷面以后改善程度逐漸增大。望虞河西岸的重度污染區(qū)嚴(yán)重阻礙望虞河水質(zhì)的改善，以虞義大橋和張橋兩個(gè)監(jiān)測(cè)點(diǎn)為甚。總體來(lái)講，調(diào)水對(duì)望虞河水質(zhì)的改善有明顯影響。

c. 采用改進(jìn)的Zadeh 公式對(duì)指標(biāo)實(shí)測(cè)值進(jìn)行規(guī)格化，并且根據(jù)樣本的實(shí)際情況選定各指標(biāo)的理想最大值和理想最小值，可有效避免各實(shí)測(cè)指標(biāo)最大值或最小值異常的問(wèn)題。

d. 調(diào)水雖然對(duì)望虞河水質(zhì)有影響，但是調(diào)水期間相鄰兩月的水質(zhì)指標(biāo)值變化并不大。本研究采用隸屬范圍更小、更精細(xì)的11 個(gè)等級(jí)的指標(biāo)質(zhì)量狀態(tài)劃分方法，比常規(guī)的劃分方法更精細(xì)，有效避免了因檢測(cè)值集中在同一水質(zhì)類別而難以判斷水質(zhì)變化趨勢(shì)的問(wèn)題。

e. 在計(jì)算水質(zhì)變化進(jìn)步度時(shí)，采用MATLAB 編程進(jìn)行計(jì)算，效率高，也方便將該方法應(yīng)用到其他樣本的估算中。

［1］任炳相，黃友璋，蔣澄宇. 江蘇省太湖污染防治概述［J］. 環(huán)境導(dǎo)報(bào)，1998(1):23-26. (REN Bingxiang，HUANG Youzhang，JIANG Chengyu. Introduction to the water pollution control of the Taihu Lake in Jiangsu Province［J］. Environment Herald，1998(1):23-26. (in Chinese))

［2］夏立忠，楊林章.太湖流域非點(diǎn)源污染研究與控制［J］.長(zhǎng)江流域資源與環(huán)境，2003(1):45-49. (XIA Lizhong，YANG Linzhang. Research on non-point source pollution in Tai Lake Region［J］.Resources and Environment in the Yangtze Basin，，2003(1):45-49.(in Chinese))

［3］劉克強(qiáng)，林澤新. 強(qiáng)化太湖流域江河湖連通促進(jìn)水資源優(yōu)化配置［J］. 中國(guó)水利，2011 (1):33. (LIU Keqiang，LIN Zexin.Strengthening connecting of river and lake in Taihu Region to promote optimal allocation of water resources［J］.China Water Resources，2011(1):33. (in Chinese))

［4］翟淑華，郭孟樸.望虞河引水對(duì)太湖影響前景分析［J］.水資源保護(hù)，1996 (4):12-15. (ZHAI Suhua，GUO Mengpu. Analysis of prospect of influence of water diversion of Wangyu River to Taihu Lake［J］. Water Resources Protection，1996(4):12-15.(in Chinese))

［5］汪榮鑫. 隨進(jìn)過(guò)程［M］. 西安:西安交通大學(xué)出版社，1987.

［6］陸傳賚. 排隊(duì)論［M］. 北京:北京郵電學(xué)院出版社，1994.

［7］鐘政林，曾光明，卓利，等.馬爾科夫過(guò)程在河流綜合水質(zhì)預(yù)報(bào)中的應(yīng)用［J］. 環(huán)境工程，1997(2):41-44.(ZHONG Zhenglin，ZENG Guangming，ZHUO Li，et al.Markov process applied to prediction for river comprehensive water quality ［J］. Environmental Engineering，1997(2):41-44.(in Chinese))

［8］馬占青，崔廣柏，楊宏杰，等.城市污水排放的灰色馬爾柯夫預(yù)測(cè)模型［J］. 河海大學(xué)學(xué)報(bào):自然科學(xué)版，2000，28(5):49-53. (MA Zhanqing，CUI Guangbo，YANG Hongjie，et al. Grey Markonian model for prediction of urban sewage drainage［J］. Journal of Hohai University:Natural Science，2000，28(5):49-53.(in Chinese))

［9］蔣承儀.灰色馬爾柯夫預(yù)測(cè)模型［J］.重慶建筑大學(xué)學(xué)報(bào)，1996(3):117-123. (JIANG Chengyi. Grey Markov forecasting model［J］. Journal of Chongqing Jianzhu University，1996(3):117-123.(in Chinese))

［10］何斌，陳燦，高登好. 水質(zhì)動(dòng)態(tài)評(píng)價(jià)的馬爾柯夫方法［J］.環(huán)境工程，2003，21(3):60-62. (HE Bin，CHEN Can，GAO Denghao. Markov method of dynamic assessment on water quality ［J］. Environmental Engineering，2003，21(3):60-62.(in Chinese))

［11］馮為，鄒志紅.河流水質(zhì)的動(dòng)態(tài)馬爾柯夫評(píng)價(jià)［J］.環(huán)境工程學(xué)報(bào)，2007，1 (8):132-135. (FENG Wei，ZOU Zhihong. Dynamic evaluation for water quality of rivers based on Markov process ［J］. Chinese Journal of Environmental Engineering，2007，1 (8):132-135. (in Chinese))

［12］ MA Huiqun，LIU Ling，CHEN Tao. Assessment model based on Markov chain［R］. Piscataway:IEEE Computer Society，2008.

［13］朱明.數(shù)據(jù)挖掘［M］.合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社，2002.

［14］崔振才，田文苓.模糊數(shù)學(xué)在水科學(xué)應(yīng)用中的幾個(gè)問(wèn)題［J］.河北工程技術(shù)高等專科學(xué)校學(xué)報(bào)，2002(2):8-12.(CUI Zhencai，TIAN Wenling.Analysis of some problems in hydroelectric energy with regard to fuzzy set［J］.Journal of Hebei Engineering and Technical College，2002(2):8-12.(in Chinese))

［15］黃寶宏，孫健，魏登云.馬爾可夫(Markov)鏈在體育教學(xué)評(píng)價(jià)中的應(yīng)用改進(jìn)［J］.天津體育學(xué)院學(xué)報(bào)，2005，20(2):33-35.(HUANG Baohong，SUN Jian，WEI Dengyun.Improve the application of the Markov chain method on teaching evaluation in physical education［J］. Journal of Tianjin Institute of Physical Education，2005，20(2):33-35.(in Chinese))

［16］GB3838—2002 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)［S］.