宮立剛

(瓦房店市水利勘測設(shè)計(jì)有限公司，遼寧 瓦房店 116300)

東風(fēng)水庫下游臨岸潛流帶水流場響應(yīng)規(guī)律研究

宮立剛

(瓦房店市水利勘測設(shè)計(jì)有限公司，遼寧 瓦房店 116300)

水庫下游臨岸潛流帶水流場響應(yīng)規(guī)律可用于預(yù)測不同水庫運(yùn)行情景下臨岸潛流帶的水頭與流量變化情況。本文以復(fù)州河?xùn)|風(fēng)水庫為例，基于測量數(shù)據(jù)得到其下游某斷面臨岸潛流帶的水頭及流量響應(yīng)過程，提出響應(yīng)強(qiáng)度指標(biāo)的概念，得到簡單直觀的響應(yīng)規(guī)律經(jīng)驗(yàn)公式，從而為該水庫下游臨岸系統(tǒng)的生態(tài)建設(shè)提供依據(jù)，并為水庫運(yùn)行方案的生態(tài)影響評價(jià)提供參考。

東風(fēng)水庫；臨岸潛流帶；水流場；響應(yīng)規(guī)律

1 概 述

潛流帶是水分飽和的沉積物層，它是地表水與地下水的互相作用區(qū)域，對河流生態(tài)系統(tǒng)起著重要的調(diào)控作用[1-2]。河流潛流帶的研究一般著重于河流地表水與河床地下水之間的垂向作用，而關(guān)于河流地表水與臨岸潛流帶之間的側(cè)向作用則研究得較少[3]。河道與河流臨岸潛流帶之間存在一定的物質(zhì)與能量交換過程，同時(shí)河流臨岸潛流帶又對河流生態(tài)系統(tǒng)內(nèi)的生物健康具有十分關(guān)鍵的影響[4-5]。近年來，水利生態(tài)研究已成為我國水利事業(yè)建設(shè)的一項(xiàng)重點(diǎn)內(nèi)容，而水庫的建設(shè)與運(yùn)行對流域生態(tài)系統(tǒng)的影響也成為當(dāng)前的研究熱點(diǎn)[6-7]。本文以東風(fēng)水庫為例，于2016年7月在水庫下游復(fù)州河設(shè)置的斷面進(jìn)行河道與臨岸潛流帶的測量，并基于測量數(shù)據(jù)分析該水庫下游臨岸潛流帶的水流場響應(yīng)規(guī)律，為該水庫的生態(tài)影響分析及該流域的生態(tài)承載力評價(jià)提供基礎(chǔ)。

2 測量與分析方法

2.1 工程概況

東風(fēng)水庫位于復(fù)州河中游，地處東經(jīng)121°54′、北緯39°48′。壩址以上控制面積663km2，設(shè)計(jì)標(biāo)準(zhǔn)100年一遇，校核標(biāo)準(zhǔn)2000年一遇，總庫容1.4億m3。該水庫是一座以灌溉和防洪為主，兼顧工業(yè)供水、發(fā)電和養(yǎng)殖等作用的大(2)型水庫。

水庫基本特征參數(shù)見表1，其中Ss為設(shè)計(jì)標(biāo)準(zhǔn)，Sj為校核標(biāo)準(zhǔn)，Vt為總庫容，Vn為正常庫容，Hn為正常蓄水位，Hs為設(shè)計(jì)水位，Hj為校核洪水位。

表1 東風(fēng)水庫基本特征參數(shù)

水庫所在流域?qū)俦睖貛Ъ撅L(fēng)區(qū)，多年平均降水量627.20mm，多年平均蒸發(fā)量1780.10mm，多年平均氣溫9.50℃，多年平均相對濕度66%，最大凍土深115cm，多年平均無霜日254天。

2.2 水流場響應(yīng)試驗(yàn)方法

2016年7月，在東風(fēng)水庫下游設(shè)置斷面進(jìn)行河道與臨岸潛流帶的測量，該斷面位于瓦房店太陽鄉(xiāng)，地處東經(jīng)121°54′11″、北緯39°46′27″，距離東風(fēng)水庫大壩約1.3km。沿河岸帶直線布置6個(gè)觀測井，分別編號為A、B、C、D、E、F，相對于河道與河岸交界處的距離分別為2m、4m、6m、8m、10m和12m。在各觀測點(diǎn)井及河道內(nèi)布置壓力傳感器以監(jiān)測水頭變化情況，該傳感器具備自動記錄與存儲功能，每0.5h捕捉1次數(shù)據(jù)。觀測時(shí)間為10天，共采集3360組水位數(shù)據(jù)?；谠摂?shù)據(jù)，可分析東風(fēng)水庫下游臨岸潛流帶水流場的響應(yīng)規(guī)律。

2.3 水流場響應(yīng)分析方法

水流場響應(yīng)可分為水頭響應(yīng)和流量響應(yīng)兩部分，這兩類響應(yīng)過程均源自河流與臨岸帶潛水層之間的側(cè)向作用，可用一維線性博欣內(nèi)斯克方程表示[8]：

(1)

式中h——水頭；

t——時(shí)間;

D——水力擴(kuò)散系數(shù);

x——距離。

水力擴(kuò)散系數(shù)D與河道和臨岸潛水帶的基本水力特性相關(guān)，表達(dá)式為

(2)

式中K——滲透系數(shù)，可基于測量數(shù)據(jù)采用數(shù)據(jù)擬合方法確定;

d——含水層飽水厚度;

ω——給水度。

在臨岸潛流帶水頭響應(yīng)問題中，式(1)具有兩個(gè)基本邊界條件:

(3)

(4)

式中h(x=0)——距離河岸0m處的水頭;

h0——河道水頭。

式(3)的含義為：河道河岸交界處含水層的水頭與河道處相同；式(4)的含義為：距離河岸無限遠(yuǎn)處的潛流帶水頭不受河道水頭的影響。

臨岸潛流帶流量響應(yīng)問題的分析公式為

(5)

式中Q(t)——某一時(shí)刻河岸單位長度上的體積流量。

3 結(jié)果與討論

3.1 水頭響應(yīng)過程

根據(jù)觀測傳感器數(shù)據(jù)可得到河道及各觀測井內(nèi)的實(shí)時(shí)水頭動態(tài)變化過程，對所收集的3360組數(shù)據(jù)進(jìn)行分類匯總，得到觀測時(shí)段內(nèi)日均水頭變化幅度值(見表2)。

表2 各位置處日均水頭變化幅度值 單位：mm

由表2數(shù)據(jù)可知，在觀測時(shí)間范圍內(nèi)東風(fēng)水庫下游觀測斷面處河道內(nèi)的日均水頭變化幅度值介于255～276mm之間，臨岸2 m處(觀測點(diǎn)A)至臨岸12m處(觀測點(diǎn)F)潛流帶日均水頭變化幅度值分別介于222～245mm之間、166～188mm之間、113～127mm之間、68～245mm之間、37～42mm之間和18～23mm之間。可見，東風(fēng)水庫下游臨岸潛流帶內(nèi)的水頭變化值低于河道本身，且距離越遠(yuǎn)該值越小。

河道內(nèi)日均水頭變化值的波動范圍為21mm，臨岸2 m處(觀測點(diǎn)A)至臨岸12 m處(觀測點(diǎn)F)潛流帶日均水頭變化值的波動范圍分別為23mm、22mm、14mm、9mm、5mm和5mm?？梢?，臨岸帶潛流帶的水頭變化值波動范圍可大于河道本身而不是呈線性減小的規(guī)律，說明水庫的運(yùn)行可引發(fā)較為復(fù)雜的臨岸潛流帶水頭響應(yīng)過程。各觀測點(diǎn)的波動值不盡相同，說明潛流帶的水頭響應(yīng)過程具有明顯的非均質(zhì)性；但總體而言，距離越大波動值越小。

綜上所述，東風(fēng)水庫下游臨岸潛流帶的水頭響應(yīng)具有明顯的非均質(zhì)性，總體而言，臨岸潛流帶的水頭響應(yīng)強(qiáng)度隨距離的增加而減弱。

3.2 流量響應(yīng)過程

表3 各觀測點(diǎn)日均側(cè)向潛流水體體積 單位：m2/d

基于監(jiān)測結(jié)果利用式(5)可得到各觀測井內(nèi)的實(shí)時(shí)流量動態(tài)變化過程，對所得的實(shí)時(shí)流量值Q(t)進(jìn)行積分，得到觀測時(shí)段內(nèi)各觀測井單位長度的日均側(cè)向潛流水體體積，見表3，其中，正值表示由河道流向臨岸潛流帶，負(fù)值反之。

由表3數(shù)據(jù)可知，在觀測時(shí)間范圍內(nèi)東風(fēng)水庫下游觀測斷面臨岸2 m處(觀測點(diǎn)A)至臨岸12 m處(觀測點(diǎn)F)單位長度日均側(cè)向潛流水體體積分別介于-6.58～3.41m2/d之間、-5.82～2.99m2/d之間、-4.4～2.3m2/d之間、-2.93～1.53m2/d之間、-1.64～0.85m2/d之間和-0.85～0.44m2/d之間?？梢?，在東風(fēng)水庫下游臨岸潛流帶內(nèi)，距離河道越遠(yuǎn)則流量越小。

臨岸2m處(觀測點(diǎn)A)至臨岸12 m處(觀測點(diǎn)F)潛流帶單位長度日均側(cè)向潛流水體體積的波動范圍分別為9.99m2/d、8.81m2/d、6.7m2/d、4.46m2/d、2.49m2/d和1.29m2/d?？梢姡▌臃秶瑯与S距離的增加而逐漸減小。

綜上所述，東風(fēng)水庫下游臨岸潛流帶的流量響應(yīng)也具有明顯的非均質(zhì)性，總體而言，在臨岸潛流帶內(nèi)，距離河道越遠(yuǎn)則流量響應(yīng)強(qiáng)度越弱。

4 響應(yīng)規(guī)律分析

提出響應(yīng)強(qiáng)度指標(biāo)的概念來分析水庫下游臨岸潛流帶水流場響應(yīng)規(guī)律，從而得到簡單直觀的經(jīng)驗(yàn)公式。該強(qiáng)度指標(biāo)可分為水頭響應(yīng)強(qiáng)度和流量響應(yīng)強(qiáng)度兩種。前者是各觀測點(diǎn)日均水頭變化幅度值與河道值的比值，后者是各觀測點(diǎn)單位長度日均側(cè)向潛流水體體積與觀測井A處值的比值。各位置處的響應(yīng)強(qiáng)度計(jì)算結(jié)果如下圖所示，其中，橫坐標(biāo)為河岸某點(diǎn)到河道與河岸交界處的距離，縱坐標(biāo)為響應(yīng)強(qiáng)度指標(biāo)值。響應(yīng)強(qiáng)度指標(biāo)值介于0～1之間，當(dāng)該值為0時(shí)表示無響應(yīng)，當(dāng)該值為1時(shí)表示臨岸潛流帶水流場變化與河道完全一致，因此，該值越大響應(yīng)強(qiáng)度越高。由圖可知，水頭響應(yīng)強(qiáng)度指標(biāo)和流量響應(yīng)強(qiáng)度指標(biāo)皆隨距離的增大而減小，其響應(yīng)規(guī)律皆可由對數(shù)曲線進(jìn)行擬合。水頭和流量響應(yīng)強(qiáng)度指標(biāo)的實(shí)際值序列與擬合曲線之間的確定系數(shù)R2值分別為0.9819和0.9855，該值較大，說明其擬合精確度滿足應(yīng)用要求。

東風(fēng)水庫下游臨岸潛流帶水流場響應(yīng)規(guī)律分析圖

在實(shí)際應(yīng)用中，可根據(jù)上圖所示擬合方程對東風(fēng)水庫下游臨岸潛流帶水流場響應(yīng)過程進(jìn)行推演，從而預(yù)測出不同水庫運(yùn)行情景下臨岸潛流帶不同位置處的水位與流量變化情況，為臨岸帶植物栽培種類與技術(shù)的選擇提供依據(jù)，并為不同水庫運(yùn)行方案的生態(tài)影響評價(jià)提供參考。

5 結(jié) 語

本文通過對復(fù)州河的東風(fēng)水庫下游設(shè)置的斷面河道及臨岸潛流帶水位進(jìn)行測量，得到不同觀測處的水頭及流量影響過程。觀測結(jié)果表明：該水庫下游臨岸潛流帶的水流場響應(yīng)具有明顯的非均質(zhì)性，且其強(qiáng)度隨距離的增加而逐漸減弱。提出響應(yīng)強(qiáng)度指標(biāo)的概念，得到簡單的經(jīng)驗(yàn)公式來量化地表示該案例中的水流場響應(yīng)規(guī)律。該經(jīng)驗(yàn)公式可用于初步預(yù)測不同水庫運(yùn)行情景下臨岸潛流帶不同位置處的水位與流量變化情況，但它只針對東風(fēng)水庫，在以后的研究中也可以對更多的水庫進(jìn)行分析，從而得到具有普遍適用性的臨岸潛流水流場響應(yīng)規(guī)律公式，以提高該領(lǐng)域的判斷與分析效率。

[1] 楊小全,金光球,李凌,等.河流潛流帶中膠體遷移的研究進(jìn)展[J].水利水電科技進(jìn)展,2010,30(6):78-83.

[2] 吳健,黃沈發(fā),唐浩,等.河流潛流帶的生態(tài)系統(tǒng)健康研究進(jìn)展[J].水資源保護(hù),2006,22(5):5-8.

[3] 金光球,李凌.河流中潛流交換研究進(jìn)展[J].水科學(xué)進(jìn)展,2008,19(2):285-293.

[4] 袁興中,羅固源.溪流生態(tài)系統(tǒng)潛流帶生態(tài)學(xué)研究概述[J].生態(tài)學(xué)報(bào),2003,23(5):956-964.

[5] 武玉峰.河流健康評價(jià)及保護(hù)對策分析——以牤牛河為例[J].水利規(guī)劃與設(shè)計(jì),2015(10):33-38.

[6] 賈瑞紅.水庫生態(tài)調(diào)度的措施及其發(fā)展前景[J].水利技術(shù)監(jiān)督,2011(5):42-44.

[7] 趙文耀.丹江口水庫流域生態(tài)環(huán)境保護(hù)現(xiàn)狀研究[J].水利技術(shù)監(jiān)督,2006,14(6),44-47.

[8] 朱蓓,趙堅(jiān),陳孝兵,等.水庫運(yùn)行對下游河岸潛流帶水位-溫度影響研究[J].水利學(xué)報(bào),2015,46(11):1337-1343.

Water Flow Field Response Regularity Research of the Near Shore Undercurrent Zone at Dongfeng Reservoir Downstream

GONG Ligang

(WafangdianCityHydraulicInvestigationandDesignCo.,Ltd.,Wafangdian116300,China)

Water flow field response regularity of near shore undercurrent zone at reservoir downstream can be used to forecast the changing condition of delivery head and flow rate in different reservoir operation circumstances. The article takes Fuzhou River Dongfeng Reservoir as an example, get the delivery head and flow rate response process of its downstream near shore undercurrent zone on the basis of measured data, put forward concept of response intensity index, get simple and visual response regularity empirical formula, thereby provide the basis for the ecological construction of the reservoir downstream near shore system and provide reference for ecological impact assessment of the reservoir operation scheme.

Dongfeng reservoir; near shore undercurrent; water flow field; response regularity

10.16617/j.cnki.11-5543/TK.2017.04.014

TV62

A

1673-8241(2017)04- 0053- 04