莊 楊， 楊 超， 張 興， 李一如

(1.常州市城市防洪工程管理處， 江蘇 常州 213165； 2.江蘇常州經濟開發(fā)區(qū)農業(yè)水利服務中心， 江蘇 常州 213165)

根據(jù)《水污染防治行動計劃》要求，地級及以上城市于2020年底前完成黑臭水體治理目標。為改善市區(qū)河道水質，常州市開展“暢流活水”系列工作，以加速水體更新為主要工作任務?！皶沉骰钏笨茖W引用長江優(yōu)質水源，“引清入城”的模式提高了城區(qū)河道水環(huán)境更新效率，水體有望在有序流動中恢復自凈能力。

河道淤積是導致黑臭水體發(fā)生的原因之一，淤泥會對河道水質構成持久的污染，且直接影響河道灌溉、排水功能的正常發(fā)揮以及河網生態(tài)效應，清淤是河道綜合治理的先行條件。目前，常見的河道淤積處理主要集中于工程措施，這一治理手段耗時長、投入高。“引清入城”的調水模式，加大了調水區(qū)域內河網水流量，加速了水流流速。調水期間，沿線河網水位有不同程度的抬高，主城區(qū)河道流速明顯提高，達到平均10～20 cm/s?；趥鹘y(tǒng)泥沙動力學理論，主城區(qū)在換水過程中，增大了水流的動能，水流挾沙能力也隨之加強。由此推斷，通過優(yōu)化水利工程調度，增加城區(qū)河道水動力條件，使其在提升地區(qū)水質過程中發(fā)揮前期效用，在改善水質的同時，降低地區(qū)淤積風險，為從源頭治理水體污染的總目標爭取時間。

為科學評估常州市主城區(qū)調水工作產生的減淤效果，目前分析沖淤效果的常見方法包括河道斷面勘測，檢測輸水斷面的懸移質含量，分析水流形態(tài)變化對沖淤能力的影響。這些方法雖然檢測結果較為準確，但是過程繁瑣，試驗要求和成本較高，檢測周期長。本文利用傳統(tǒng)泥沙動力學模型，分析調水前、后水體沖淤和挾沙能力的變化，對調水產生的沖淤效果進行初步評價，為探尋在調水中逐步達成地區(qū)沖淤平衡的狀態(tài)提供理論支撐。

1 工程設施

2017年7月，常州市審議通過《常州市運北片主城區(qū)“暢流活水”總體方案》，方案決定：活水工程引調長江水，采用雙源(德勝河、澡港河)供水，修建4座(盤龍苑、恐龍園、新市橋、洋橋)活動堰，增強了內部調控能力，提高小河道流動性，打造兩進三出的水流路徑(“兩進”指德勝河和澡港河，“三出”指肖龍港、老桃花港和新溝河)，配合水利工程的整體調度。工程建成以后，從澡港河引流(20～40 m3/s)的水資源全部流經主城區(qū)，按每天調運6 h，每年保證運行250 d估算，全程河網平均引水量約為10 m3/s。調水工作增大了市區(qū)河道流量，利用水流勢能與動能的轉換加速水流流動，加快了水循環(huán)，為通過非工程措施進行沖淤營造了契機。

2 水情變化

現(xiàn)有研究結果顯示，常州城區(qū)適宜引水規(guī)?？蛇_50 m3/s。在這一引水量下，“暢流活水”方案目前可以實現(xiàn)人為重構水位差20 cm以上。但由于南運河、串新河等南排泵站的啟用，通常換水區(qū)全河段較引水前抬高約7～11 cm，城南排水河道出水側水位抬高不明顯，以2020年12月的暢流活水調度運行水位數(shù)據(jù)為例，詳見表1。

表1 主城區(qū)平均水位變化

同時，引水期間，連通的河道基本形成由北往南、由西向東的有序流動(平均流速由1～7 cm/s提升至9～25 cm/s)。古城區(qū)內河道12～36 cm/s，其余暢通河道流速基本在10 cm/s以上，西北部區(qū)域總體流速7～20 cm/s。

3 水流沖淤和挾沙能力變化

挾沙水流的運動是一種相當復雜的現(xiàn)象，影響因素主要有：渾水性質、水流緩急程度、過流斷面寬深比、過流邊界糙率、水力半徑以及水沙綜合因素等。對于懸浮質泥沙，紊流理論認為懸浮在水體中的泥沙，受地心引力的影響有下沉趨勢，卻因受紊動擴散作用而上升，在運動過程中，水流和泥沙相互發(fā)生作用。一方面水流通過消耗其紊動能量挾帶和輸送泥沙，另一方面泥沙的存在又反過來改變水流的物理性質和水流結構，從而影響水流的能量、流速和含沙濃度?，F(xiàn)有實測資料顯示：長江水1982—1985年，含沙量范圍為0.20～3.76 kg/m3，2003—2004年，含沙量范圍為0.12～0.20 kg/m3。含沙量隨著沿岸綠化工程的推進還呈現(xiàn)逐年降低的態(tài)勢，按照現(xiàn)有的含沙量劃分標準[1]，可將長江水系劃分為低含沙量河流，本文假設調水區(qū)域內泥沙為均勻沙。

3.1 沖淤能力提升

水中泥沙常見的3種形態(tài)為懸移質、推移質、河床質，河流中泥沙以懸移質為主。對于水流挾沙能力的探討，多基于懸浮功的概念，判斷標準是水流對懸移質泥沙的運輸能力，在水流動能加大時會引發(fā)河床質泥沙轉化成懸移質、推移質，隨水流運動，3種形態(tài)沒有嚴格的界線，可以隨水流條件的變化而相互轉化。

懸移質泥沙隨水流移動主要依靠水流的能量轉換，要使得河床保持沖淤平衡狀態(tài)，單位體積挾沙水流總耗能需等于懸浮泥沙消耗的能量與其他途徑消耗能量之和。

河道斷面平均流速：

(1)

泥沙啟動流速：

(2)

泥沙揚動流速：

(3)

式中：ρ為水體密度，kg/m3；ρs為粒子密度，kg/m3；d為粒子直徑，m；h為水深,m；ω為泥沙顆粒沉速，m/s，可根據(jù)岡恰洛夫公式確定其取值；z為懸浮指標，泥沙粒徑對其影響較大,將泥沙視為均勻沙時，可以根據(jù)羅斯公式確定z值。

想通過調度實現(xiàn)沖淤效果，要加大U值。根據(jù)實測資料，長江水中懸移質泥沙最大粒徑常出現(xiàn)在夏汛期，懸移質粒徑范圍為0.0071～0.0276 mm，冬季則處于最小值狀態(tài),懸移質粒徑范圍為0.0064～0.0150 mm[2]，當d為8.5×10-7m時，z取0.7，ω取值0.0008 m/s。在12月份的調水過程中，水流平均流速由0.01～0.07 m/s提升至0.09～0.25 m/s，古城區(qū)內河道部分流段高達0.36 m/s，滿足U>Us>Uc。

3.2 水流挾沙能力變化

借鑒三峽工程庫區(qū)淤積測算和長江深水航道整治工程估算泥沙回淤等案例經驗，選取適用的挾沙能力計算公式[3]。張瑞瑾由能量守恒原理得出的流速、泥沙沉速、區(qū)域水深函數(shù)關系的一般表達式對低含沙量水流計算準確性較高，適用于研究常州地區(qū)調水期間的挾沙量變化規(guī)律，計算式為

(4)

式中：Sv為按體積比濃度計算的水流挾沙能力；k、m為系數(shù)；v為水流的垂線平均流速，m/s；g為重力加速度，m/s2。其中，k影響因素復雜，變化幅度較大，m的取值可以定位在0～3之間。

參考孔俊等[4]針對長江水流律定的挾沙能力公式系數(shù)的最佳結果，本文將張瑞瑾公式中k取值定為0.317，m值定為0.44。用張瑞瑾公式分析調水過程中水流挾沙能力變化，在k、m、ω取值均確定的情況下，影響水中懸移質數(shù)量的主要因素是流速和水深，就同一河段而言，調水過程中水流的挾沙能力遠超調水前。

分析可知，Sv值受流速的影響大于水深，對于水位差最大和最小的流段，調水前和調水時S值走向基本一致。調水開始后，隨著流量的增加，S值可以達到[0.0658，0.2577)，與調水前S值所處的[0.0037，0.0482)相比，水流挾沙能力增加效應明顯。

4 結 論

綜上所述，引流量較大的河道，水流加速明顯，挾沙能力也成倍速增加?！皶沉骰钏惫ぷ髟诟纳剖袇^(qū)水質的同時，能有效加大水流沖淤和挾沙能力，在水流沖刷下懸移質粒徑會隨之變大，過去因流速緩慢，從河水分離出來導致河道淤積的部分泥沙可以轉換為懸移質，并被過境水流帶走。

枯水期開展“暢流活水”在一定程度上增強了水流沖淤效果，適當利用可以改善淤積，從而降低工程措施清淤頻率，帶來一定的經濟價值。后期可以進一步通過長期的斷面輸沙量和檢測試驗，監(jiān)測水流挾沙量變化值和易淤流段的淤積改善情況，為實現(xiàn)非工程清淤提供數(shù)據(jù)支撐。也可以嘗試研究不同調水量、調水模式下，水流狀態(tài)與挾沙能力的變化，從而實現(xiàn)模型到實際應用的轉化，對“暢流活水”方案進一步提出優(yōu)化建議。