王晶娥

(遼寧省阜新水文局，遼寧 阜新 123000)

河流水沙演變規(guī)律分析是河流綜合治理規(guī)劃的重要依據(jù)。北方旱區(qū)河流存在水少沙多的特點(diǎn)，水沙演變規(guī)律分析對(duì)此類河流治理更為重要。當(dāng)前，對(duì)水沙演變規(guī)律分析主要是結(jié)合實(shí)測(cè)數(shù)據(jù)，利用趨勢(shì)檢驗(yàn)[1- 5]及突變檢驗(yàn)方法[6- 9]進(jìn)行多尺度的分析。柳河是遼河的一級(jí)支流，位于遼寧西北部干旱沙丘地區(qū)，屬于水土流失嚴(yán)重區(qū)域，區(qū)域水土流失量高達(dá)每年1560萬(wàn)t，近年來(lái)遼寧省開(kāi)展多沙河流的綜合治理，在綜合治理規(guī)劃時(shí)急需對(duì)多沙河流水沙演變規(guī)律進(jìn)行分析。為此本文結(jié)合柳河干流主要控制站彰武水文站1968～2016年實(shí)測(cè)水沙數(shù)據(jù)，對(duì)柳河近49年的水沙演變特征進(jìn)行定量分析。

1 區(qū)域概況及研究方法

1.1 區(qū)域概況

柳河從西北流向東南，流域內(nèi)水土流失嚴(yán)重，洪水過(guò)程陡漲陡落，洪水歷時(shí)短，河床沖淤變化特別劇烈。彰武水文站水文位于柳河中下游，在上游57km處建有鬧得海大型水利樞紐工程，鬧得海水庫(kù)至彰武站區(qū)間面積971km2，區(qū)間有兩條較小支流加入，還有大清溝中型水庫(kù)一座。具體見(jiàn)圖1。

圖1 柳河站點(diǎn)分布圖

1.2 研究方法

1.2.1趨勢(shì)非線性檢驗(yàn)方法

該方法假設(shè)時(shí)間數(shù)據(jù)序列(x1，x2，…，xn)為獨(dú)立、隨機(jī)同變量分布的樣本序列，檢驗(yàn)統(tǒng)計(jì)變量S計(jì)算方程為：

(1)

(2)

(3)

(4)

式中，S—正態(tài)分布的檢驗(yàn)統(tǒng)計(jì)變量；xi和xj—同一樣本中不同分布的兩個(gè)系列，其中E(ck)=k(k-1)/4；σs—標(biāo)準(zhǔn)方差；sign—運(yùn)算符號(hào)；n—樣本的總個(gè)數(shù)；Z—檢驗(yàn)值，若Z>0，檢驗(yàn)的時(shí)間序列為上升變化趨勢(shì)，Z<0，檢驗(yàn)的時(shí)間序列為上下降變化趨勢(shì)；Z的絕對(duì)值大于或等于2.32、1.64、1.28，表示通過(guò)置信度分別為99%、95%以及90%的顯著性檢驗(yàn)水平。

1.2.2突變分析方法

突變分定義一個(gè)統(tǒng)計(jì)量ck來(lái)分析n個(gè)獨(dú)立的、隨機(jī)變量同分布的樣本數(shù)據(jù)序列(x1，x2，…，xn)，計(jì)算方程分別為：

(5)

E(ck)=k(k-1)/4

(6)

σ(ck)=k(k-1)(2k+5)/72

(7)

(8)

式中，m—第i個(gè)樣本xi>xj(1≤j≤i)的累積數(shù)；E(ck)、σ(ck)—ck的均值和方差；UFk—ck的標(biāo)準(zhǔn)化；按逆序列數(shù)據(jù)(xn，xn-1，…，x1)重復(fù)上面過(guò)程，使|UBk=-UFk，k=n，n-1，…，UB1=0。如果UF和UB兩條曲線在置信區(qū)間內(nèi)出現(xiàn)交點(diǎn)，即為可能的突變點(diǎn)。

2 水沙演變特征分析

2.1 趨勢(shì)特征分析

(1)年尺度趨勢(shì)分析

結(jié)合彰武站1968～2016年實(shí)測(cè)水沙數(shù)據(jù)，分年尺度和汛期尺度對(duì)彰武站水沙趨勢(shì)演變特征進(jìn)行分析，如圖2、3，并結(jié)合趨勢(shì)統(tǒng)計(jì)方法對(duì)趨勢(shì)特征進(jìn)行檢驗(yàn)，檢驗(yàn)結(jié)果見(jiàn)表1。

圖2 水沙年尺度分析結(jié)果

時(shí)間尺度徑流量/108m3輸沙量/106t年尺度-5.49-6.65汛期-4.66-5.82非汛期-5.99-6.99

圖3 水沙月尺度分析結(jié)果

圖5 不同時(shí)間尺度的彰武站水沙突變特征結(jié)果

圖2為水沙年尺度趨勢(shì)特征分析結(jié)果，從圖2可看出，柳河彰武站在過(guò)去的49年中，徑流和輸沙量呈現(xiàn)明顯遞減趨勢(shì)，遞減的幅度為-4×106m3/a，輸沙量逐年遞減的幅度為-3×105t/a，年水沙趨勢(shì)檢驗(yàn)統(tǒng)計(jì)值為-5.49和-6.65，絕對(duì)值均高于趨勢(shì)檢驗(yàn)統(tǒng)計(jì)Z達(dá)到99%的值，經(jīng)檢驗(yàn)達(dá)到顯著遞減。從汛期和非汛期看，和年水沙變化相似，彰武站汛期和非汛期水沙下降趨勢(shì)也十分明顯，檢驗(yàn)統(tǒng)計(jì)值也均通過(guò)99%的遞減顯著性檢驗(yàn)，非汛期水沙遞減幅度分別為-3×105m3/a、3000t/a，而汛期水沙遞減幅度分別為-8×105m3/a、75000t/a?？梢?jiàn)，年水沙遞減幅度明顯高于汛期、非汛期，汛期水沙遞減幅度明顯高于非汛期。

2.2 水沙相關(guān)性特征

在水沙趨勢(shì)特征分析的基礎(chǔ)上，對(duì)柳河彰武水文站汛期和非汛期的水沙相關(guān)特征進(jìn)行分析，如圖4所示。

圖4 彰武站水沙汛期、非汛期相關(guān)性特征

柳河汛期和非汛期水沙關(guān)系均在0.5以上，汛期受水量影響，其相關(guān)系數(shù)低于非汛期，在汛期，彰武站受上游水庫(kù)控制影響，水沙影響較大，致使水沙相關(guān)性較差。而在非汛期，由于上游來(lái)水量的減少，上游水庫(kù)對(duì)下游彰武站水沙關(guān)系影響較弱，水沙相關(guān)性逐步增強(qiáng)。

3 水沙突變特征分析

結(jié)合突變分析方法對(duì)柳河干流彰武水文站1968～2016年水沙突變特征進(jìn)行分析，如圖5所示。

從彰武站各時(shí)間尺度突變特征分析可看出，彰武站年水沙突變年份較為一致，主要發(fā)生在2006年，經(jīng)過(guò)調(diào)查分析，突變?cè)蛟谟?006年擴(kuò)大城市引供水量后，水庫(kù)供水量增加，使得下游彰武站的徑流出現(xiàn)突變。徑流的變化也直接導(dǎo)致了彰武站輸沙量的突變。從汛期和非汛期水沙突變可看出，汛期徑流突變和年突變較為一致，突變?cè)蛑饕谟谌祟惢顒?dòng)的影響。而汛期水沙未出現(xiàn)突變，主要是汛期沙量較大，一部分輸沙量由于上游水庫(kù)供水影響，沉積在了水庫(kù)內(nèi)，而后在水庫(kù)放水后又重新進(jìn)入河道，使得汛期總體輸沙量為出現(xiàn)明顯突變。

4 水沙演變成因分析

為對(duì)柳河水沙演變的成因進(jìn)行分析，首先結(jié)合水沙雙累積曲線分析人類活動(dòng)對(duì)彰武站水沙演變影響，結(jié)果如圖6所示。

水沙雙累積曲線綜合反映水沙變化的特征，受降水影響，水沙雙累積曲線的斜率不發(fā)生變化，而當(dāng)受到人類活動(dòng)影響，水沙雙累

圖6 彰武水沙雙累積曲線分析結(jié)果

積曲線的斜率發(fā)生偏移。從圖6可看出，非汛期柳河主要控制站彰武站的水沙雙累積曲線斜率較為一致，為發(fā)生偏移。因此在非汛期水沙演變主要原因在于區(qū)域徑流量的遞減，而區(qū)域徑流量的遞減一部分在于氣候變化，另外由于近些年來(lái)，阜新地區(qū)加大了水保措施的治理，增加流域蓄水能力，使得區(qū)域徑流量出現(xiàn)明顯的遞減趨勢(shì)。而在汛期，其雙累積曲線各散點(diǎn)變化十分紊亂，無(wú)法定線，可見(jiàn)在汛期人類活動(dòng)對(duì)彰武站水沙關(guān)系影響十分顯著。彰武水文站上游有鬧德海水庫(kù)，該水庫(kù)除防洪功能外，重要的一項(xiàng)功能就是供水，在汛期水量較大的時(shí)候，供水量增加，使得下游河道水量和輸沙量呈現(xiàn)明顯變化，另一方面，由于水土保持措施治理影響，流域下墊面發(fā)生較為明顯的變化，使得彰武站的汛期水沙變化較大。

5 結(jié)論

(1)柳河水沙呈現(xiàn)顯著遞減趨勢(shì)，年水沙遞減幅度分別為-4×106m3/a、-3×105t/a，汛期遞減幅度小于年遞減幅度，但明顯高于非汛期的遞減幅度。

(2)汛期柳河受上游水庫(kù)控制影響，水沙相關(guān)性較弱，水沙非線性特點(diǎn)顯著，非汛期水沙相關(guān)性較好，呈現(xiàn)一定的線性變化關(guān)系。

(3)人類活動(dòng)是柳河汛期水沙演變的主因，氣候變化是非汛期水沙演變的主因。

(4)本文只定性的探討了北方典型多沙河流的水沙演變成因，在以后還需結(jié)合新方法定量探討多沙河流水沙演變的成因。

[1] 張洪義, 于海龍. 遼陽(yáng)市降雨徑流時(shí)空分布的差異性分析[J]. 吉林水利, 2012(03): 37- 39.

[2] 陸軍. 遼河中下游水沙關(guān)系變化研究[J]. 水利規(guī)劃與設(shè)計(jì), 2016(06): 38- 40+115.

[3] 王穎華, 張?chǎng)? 西寧市降水量特征及變化趨勢(shì)分析[J]. 水土保持研究, 2011, 18(05): 156- 160.

[4] 于國(guó)寶. 大凌河流域水沙特性及變化趨勢(shì)分析[J]. 水利規(guī)劃與設(shè)計(jì), 2017(02): 41- 43.

[5] 殷正雄. 水資源及其變化的區(qū)域環(huán)境影響機(jī)制[D]. 華中師范大學(xué), 2011.

[6] 李淑玲, 胡仲. 遼寧省水污染變化趨勢(shì)分析及治污舉措[J]. 現(xiàn)代農(nóng)業(yè)科技, 2010(23): 267+272.

[7] 陳水蓉. 趨勢(shì)分析在水質(zhì)管理中的應(yīng)用研究[D]. 天津師范大學(xué), 2012.

[8] 李文文, 傅旭東, 吳文強(qiáng), 等. 黃河下游水沙突變特征分析[J]. 水力發(fā)電學(xué)報(bào), 2014, 33(01): 108- 113.

[9] 王世鈞. 渭河上游水沙特性變化及其規(guī)律研究[J]. 水利規(guī)劃與設(shè)計(jì), 2013(09): 8- 10.