職保平，王 宇，秦凈凈

(1.黃河水利職業(yè)技術學院, 河南 開封 475004；2.小流域水利河南省高校工程技術研究中心， 河南 開封 475004)

黑崗口水庫水沙特性的數(shù)值模擬

職保平1,2，王 宇1,2，秦凈凈1

(1.黃河水利職業(yè)技術學院, 河南 開封 475004；2.小流域水利河南省高校工程技術研究中心， 河南 開封 475004)

小浪底電站建成后，作為開封市水源地的黑崗口水庫不再進行排沙，為了研究黑崗口水庫的多年的水沙沉積及分布問題，通過對花園口水文站觀測數(shù)據(jù)及黑崗口引流量數(shù)據(jù)的分析，結合黑崗口地形數(shù)據(jù)，依據(jù)插值和擬合數(shù)學理論，利用Flow-3D軟件進行建模并進行數(shù)值計算得到水庫泥沙的空間分布，并推導預測黑崗口水庫未來的水沙分布及趨勢分析，對黑崗口水庫的管理及優(yōu)化引水工作方式等提供依據(jù)，具有一定的現(xiàn)實意義。

黑崗口水庫；水沙特性；水力特性

自小浪底建成后，黃河下游引水含沙量明顯下降，作為開封市水源地的黑崗口水庫總引沙量下降明顯，因此，基本不再進行排沙工作。雖然黃河下游含沙量下降，但經(jīng)過16年的長期引水，總引沙量依然十分可觀，目前引水口處明顯出現(xiàn)泥沙沉積效應。

水庫中水沙的長時間沉積和分布不均勻是泥沙研究中的重要工作內(nèi)容之一，而如何推算出水沙變化的過程和規(guī)律是治河泥沙工程中急需解決的問題。對于此類研究，已有很多學者利用物理模型實驗的方法進行了深入研究，給出了水庫中水沙分布變化規(guī)律的特性[1]。隨著研究工具的發(fā)展及數(shù)值模擬方法的完善，利用計算機數(shù)值模擬軟件可以更好的對該類問題進行研究，常用的水沙數(shù)學模型主要包括水流模塊和泥沙模塊。張紅藝等[2],楊程等[3]及馮淑萍等[4]通過對高含沙水庫和河道的數(shù)值模擬得出了泥沙分布相關規(guī)律，王勤香等[5]運用數(shù)值模擬方法研究了人工倒虹吸的水力特性，假冬冬等[6]通過三維數(shù)學模型系統(tǒng)的分析了水沙條件變化對河型河勢影響，程文等[7]和潘存鴻等[8]通過數(shù)值模擬分析了河道及水庫泥沙沖淤的分布。

為整體分析黑崗口水庫的水沙特性，以黑崗口水庫實測地形數(shù)據(jù)為基礎，建立完整的黑崗口水庫模型，以現(xiàn)有的原型觀測資料計算邊界條件，對黑崗口水庫進行了數(shù)值模擬，最終對未來水庫中水沙沉積和分布等問題進行研究，為庫區(qū)的運行和管理提供依據(jù)。

1 研究區(qū)域

本文研究區(qū)域以黑崗口水庫為核心，引用上游花園口水文站(1987年—2015年)數(shù)據(jù)平均值，見表1，該站是距黑崗渠首閘最近的上游水文站。以2015年徑流量和輸沙量為計算特征數(shù)據(jù)，見表2。

黑崗口水庫，又名黑池，位于小浪底樞紐下游180 km，花園口水文站下游75 km處，開封市西北部15 km處，緊臨黃河大堤，東窄西寬，長約5 km，寬500 m～1 000 m，水深3 m～8 m。黑崗口渠首閘位于黑崗口險工段30壩與31壩之間，大堤里程樁號為77+170處。該閘為鋼筋混凝土結構涵洞式水閘，共5孔，孔高2.0 m，寬1.8 m，閘室長6.85 m，閘底高程76.75 m，設計防洪水位82.64 m，設計引水流量50 m3/s，加大引水流量為64 m3/s，黑崗口1999年—2013年引水量見表3，其中2008年—2010年數(shù)據(jù)缺失。

表1 花園口水文站1987年—2015年水沙特征表

表2 花園口水文站2015年來水來沙年內(nèi)分配表

表3 黑崗口引水流量分布表

2 模擬計算模型

2.1 控制方程

本文涉及數(shù)值模型計算控制方程[7]包括連續(xù)性方程、動量方程、紊動能k方程，紊動能耗散率ε方程，泥沙輔助方程。根據(jù)計算水力學相關內(nèi)容[8-11]可知：

連續(xù)性方程：

(1)

動量方程：

(2)

紊動能k方程：

(3)

紊動能耗散率ε方程：

(4)

現(xiàn)有主要紊流模型中，Reynolds應力模型對脈動應力直接建立輸運方程，因此對于考慮各向異性及強曲率情況較為適合。但由于該類模型建模相對復雜且計算量大，在特殊情形下才考慮使用。渦黏模型中的RNG模型則通過對標準k-ε的湍流黏度進行修正，RNGk-ε模型是一種改進的k-ε模型，基于多尺度隨機過程重整化思想，方程形式同標準k-ε模型相似，模式常數(shù)由重整化理論算出，該模型可以更好的處理高應變率及流線彎曲程度較大的波動且計算量適中，基于上述原因，本研究中采用近年來逐漸成熟并被廣泛應用[12-14]的修正湍流黏度的RNG (Renormalization Group)k-ε模型。

泥沙模型的輔助方程[9]：

當來沙處于過飽和而床沙又較粗的條件下，床沙和懸沙就總體而言不發(fā)生交換，懸沙發(fā)生單向淤積。挾沙力級配是由來沙決定的，與床沙無關。

(5)

2.2 模型建立、計算范圍、網(wǎng)格剖分及邊界條件

本文地形建模資料采用2016年黑崗口水庫地形實測數(shù)據(jù)資料進行三維模型建立。水文資料采用實測2015年水位、流速、流量觀測資料，模型下邊界水位采用2016年實測流量與水位關系給定。以黑崗口水庫整體邊界為計算范圍。模型采用自適應邊界較好的正交網(wǎng)格進行模擬。網(wǎng)格均勻分布，模型計算時上游開邊界由實測流量控制，下游開邊界由實測水位控制。

圖1三維數(shù)值模型的建立

初始狀態(tài)按照黑崗口設計說明書進行設置，進出口條件按照現(xiàn)狀開封河務局提供的涵閘水文資料中黑崗口閘多年(1999年—2007年)平均引水量進行設定(Q=2.76 m3/s)模擬結果如下：

表5 模擬條件參數(shù)設置

3 模擬計算

3.1 計算工況

為了更好的了解目前以及未來黑崗口水庫的水沙變化趨勢，根據(jù)黑崗口水庫引水閘不同引水流量設定工況進行仿真計算，進口含沙量采用表2中多年平均值。模擬分析中根據(jù)不同引水量擬定工況，主要有在小引水量的情況下、在平均引水量的情況下以及大引水量的情況下，多年后黑崗口水庫水沙分布情況計算，各工況特征數(shù)據(jù)見表6。

表6 黑崗口水沙計算工況

3.2 計算結果分析

水沙仿真計算結果如下：

圖2 工況1情況下黑崗口水庫泥沙分布示意圖

圖3 工況2情況下黑崗口水庫泥沙分布示意圖

圖4工況3情況下黑崗口水庫泥沙分布示意圖

圖2～圖4為不同工況下黑崗口水庫泥沙分布模擬示意圖，其中，斷面A-A是黑崗口渠首閘處的典型分布情況，斷面B-B是湖中島西側(cè)的中心位置，斷面C-C是湖心島西側(cè)，斷面D-D是湖心島東側(cè)，斷面E-E是黑崗口至柳池的分水口斷面。分析可知：由于黃河中下游泥沙含量較大，因此黑崗口引水閘附近含沙量較大，含沙量達到2.6 kg/m3；同時由于湖心島的存在，使大量的泥沙沉積在進水端與湖心島之間(斷面A-A到C-C)，該區(qū)域沉沙量明顯高于湖心島東側(cè)的區(qū)域，這有利于向二級沉沙池柳池提供更為純凈的水源。黃河水經(jīng)渠首閘引水渠道口流入黑池，過流斷面迅速增加，流速降低顯著，因此形成較為明顯的沉沙廊道，見斷面A-A到斷面B-B之間的紅色區(qū)域，該廊道隨著時間的增長，沉沙淤積作用顯著。湖心島附近含沙量由于阻流的原因，島后形成一定的回旋作用，因此，島的東面(斷面D-D)含沙量高于其他水庫部位，在庫底形成淤積效應；柳池進口端(斷面E-E)整體懸沙、沉沙成分較少，未對輸水道形成危害。黑崗口引水閘附近常年有大量泥沙淤積，在不考慮挖沙等處理手段的情況下，仿真計算出的10年后的引水閘處(斷面標注)含沙量已經(jīng)達到8.5 kg/m3，在經(jīng)過計算后得出引水口處泥沙淤積深度達4 m左右。

4 結 論

研究得出黑崗口水庫在持續(xù)引水10年后，沖沙廊道不會至湖心島，形成鏈接，但由于水流作用，會在二者之間，形淤積并呈圓形分布，且與黑池北岸庫底聚集顯著。由于柳池進口末端(斷面E-E)水體流動較與整個黑池而言，非常緩慢，形成區(qū)域死水效應，懸移質(zhì)沙體逐步沉淀，該區(qū)域?qū)⑿纬蔀┩?，侵占水面。黑崗口引水閘附近10年后不采取排沙措施，造成引水效率下降甚至無法正常引水。該結論與郭維東等[15]研究結果不同，說明黑崗口水庫雖然引水渠流量較大，但目前并不同于河道型水庫的水力特性。

[1] 段志科,王延貴,戴 清.水庫泥沙數(shù)學模型與物理模型成果比較及其應用[J].泥沙研究,1994(4):76-85.

[2] 張紅藝,楊 明,張俊華,等.高含沙水庫泥沙運動數(shù)學模型的研究及應用[J].水利學報，2001，32(11):20-25.

[3] 楊 程,李春光,景何仿,等.寧夏沙坡頭河段二維水沙運動數(shù)值模擬研究[J].水力發(fā)電，2013，39(5):8-12.

[4] 馮淑萍,楊 程.沙坡頭庫區(qū)三維水沙運移數(shù)值模擬研究[J].水力發(fā)電，2014，40(7):22-25.

[5] 王勤香,潘繼順,王鵬濤.南水北調(diào)安陽河段倒虹吸水沙二維數(shù)值模擬[J].人民黃河,2010,32(9):120-122.

[6] 假冬冬,邵學軍,周建銀,等.水沙條件變化對河型河勢影響的三維數(shù)值模擬研究[J].水力發(fā)電學報,2014,33(5):108-113.

[7] 程 文,施練東,黃 晶,等.湯浦水庫泥沙沖淤分布數(shù)值模擬[J].水利水運工程學報，2015(6):107-111.

[8] 潘存鴻,曾 劍,唐子文,等.錢塘江河口泥沙特性及河床沖淤研究[J].水利水運工程學報,2013(1):1-7.

[9] 王福軍.計算流體動力學分析[M].北京:清華大學出版社,2004:58-76.

[10] 趙相航，解宏偉，郭 馨，等.基于VOF模型的臺階式溢洪道數(shù)值模擬[J].水利與建筑工程學報,2016,14(6):143-148.

[11] 關大瑋，程香菊．基于VOF模型的溢流壩三維水流數(shù)值模擬[J].中國水運,2015,15(6):51-55.

[12] 楊明遠,嚴以新,孔 俊,等.珠江口水流泥沙運動模擬研究[M].北京:海洋出版社,2008:50-56.

[13] 朱玲玲,牧振偉,龔厚建,等.阿拉溝水庫溢洪道水流特性數(shù)值模擬[J].水利與建筑工程學報,2012,10(6):22-25.

[14] 魯克恩,葛旭峰,李 琳.小石峽水電站溢洪道消力池結構優(yōu)化的三維數(shù)值模擬[J].水利與建筑工程學報,2011,9(5):15-17.

[15] 郭維東，李曉麗，胡 艷，等.河道型水庫泥沙沖淤數(shù)值模擬分析[J].沈陽農(nóng)業(yè)大學學報,2013,44(6):799-804.

NumericalSimulationofFlowandSedimentCharacteristicsinHeigangkouReservoir

ZHI Baoping1,2, WANG Yu1,2, QIN Jingjing1

(1.YellowRiverConservancyTechnicalInstitute,Kaifeng,He'nan475004,China；2.EngineeringTechnologyResearchCenterofSmallWatershedConservancyUniversityofHe'nanProvince,Kaifeng,He'nan475004,China)

After the completion of Xiaolangdi hydropower station, as the water source of Kaifeng city, Heigangkou reservoir no longer carried out desilting. In order to study the water and sediment deposition and distribution in the Heigangkou reservoir, based on the observation data of Huayuankou hydrological station and the data of Heigangkou drainage data, combined with the data of Heigangkou mouth terrain, based on the interpolation and fitting mathematical theory, this paper adopted Flow-3D software to model and calculate the reservoir spatial distribution of reservoir sediment. The distribution and trend of water and sediment in the future of Heigangkou reservoir was analyzed. This is the basis for the management of the Heigangkou reservoir and the way of optimizing the way of water diversion, and it has certain practical significance.

Heigangkoureservoir;waterandsedimentcharacteristics;hydrauliccharacteristics

10.3969/j.issn.1672-1144.2017.05.034

2017-05-11

2017-06-17

開封市科技局科技計劃項目“汴西湖水庫對開封水系的水沙變化影響研究”(1503037)

職保平(1983—)，男，河南新鄉(xiāng)人，博士，講師，主要從事水利水電工程方面的研究工作。E-mail: clownzhi@hotmail.com

TV651.1

A

1672—1144(2017)05—0188—04