抽水蓄能电站下库长系列年泥沙淤积特性研究

Journal of Yangtze River Scientific Research Institute

长江科学院院报

Vol.35 No.1Jan. 2 0 18

2018,35(1):5-10，15

doi ： 10.11988/ckyyb.20170074

抽水蓄能电站下库长系列年泥沙淤积特性研究

孙羽1，李向阳2，孙东坡1，张晓雷1，植凤娟3

(1.华北水利水电大学水利学院，郑州450046; 2.黄河水利委员会水文局，郑

3.凯里学院建工学院，贵州凯里556011)

摘要

州

450004;

:为了把握长系列年水沙过程对抽水蓄能电站下库区河床调整的影响，以天池抽水蓄能电站下库为例开展

了不同长系列年库区冲淤变形的数值模拟。根据下库水位周期性变化的运行特点，建立了基于MIKE的水沙数学 模型，并进行了模型参数调试与验证;模拟给出了长系列年的库区淤积总量及排沙比；同时也给出了库区淤积的沿

程分布与特征河段的淤积过程。成果分析表明:库区淤积发展状况与水沙系列相互响应，泥沙淤积主要呈现三角 洲形态，淤积重心主要集中于拦河坝上游0.5~ 1.8 km范围内;不同河段的淤积特性不同，淤积物组成也不同，主要 与其所处环境的水力与边界特性有关;频率洪水在长系列年水沙过程中虽然占据时间不长,但在库区河床冲淤量 中有重要影响;要提高排沙效率、减小淤积，需要进一步调整改善下水库的洪水调度与电站运行。研究成果与同步 物理模型试验结果基本一致，可供抽水蓄能电站水库调度运行参考。

关键词:抽水蓄能电站;下水库;长系列年;数值模拟;泥沙淤积;排沙比;淤积重心 中图分类号:TV147

文献标志码:A

文章编号= 1001-5485(2018)01-0005-06

1研究背景

近年来，我国各地相继建设了一批抽水蓄能电

针对库区水流泥沙数值模拟的研究有很多，如:

彭杨等[3] (2006)针对三峡库区建立了一维非恒定的 水沙数学模型，采用三级算法求解水流方程组，选用 汊点分沙模式与泥沙连续方程组进行非耦合求解，模 型验证结果良好;孙昭华等[4] (2007)建立了长江下游 近河口段(大通至徐六泾)具有河网特点的一维水沙 数学模型，实现了同时受径流、潮流双重作用下的非 恒定水沙运动的模拟;黄仁勇等[5] (2012)考虑水沙输 移过程中的非恒定性，建立了三峡水库一维非恒定树 状河网水沙数学模型，采用三峡水库蓄水运用后 2003—2009年实测水沙资料对模型进行了可靠性验 证，表明模型可在三峡水库优化调度研究中发挥预测 作用。但是针对抽水蓄能电站下库复杂水位及运用 条件，进行1 a乃至50 a的长系列水沙过程模拟还 不多见。鉴于泥沙问题的复杂性和抽水蓄能电站下 水库运行的特殊性，选择典型抽水蓄能电站，针对下 库水沙及边界特点，建立适宜的水沙数学模型，采用 数值模拟分析研究库区泥沙淤积规律是十分必要的。

站，以改善电网的运行状况。相对常规水电站，抽水 蓄能电站水库库容较小，且上、下水库循环抽放水运 行;在河流泥沙含量较高的地区兴建抽水蓄能电站， 泥沙淤积将直接影响到水库库容的损失。由于抽水 蓄能电站水库运行水位长期在正常蓄水位和死水位 之间变动，如果没有合理的处置措施，泥沙极易淤积 在正常蓄水位和死水位之间，占据水库调节库容。同 时由于抽水蓄能电站对过机泥沙的控制要求较高，如 果泥沙淤积近电站进出水口附近，过高的含沙水流会 磨蚀机电设备，甚至影响电站正常运行。因此，对下 边界条件复杂、运用工况多样的抽水蓄能电站下库区 开展泥沙淤积问题研究，是十分必要的[1]。

由于河流泥沙问题的复杂性，以及受工程边界条 件约束后河床再调整的不确定性，目前还没有成熟的 理论去预测复杂边界条件下的水库泥沙淤积量与淤 积分布。随着计算机迅速发展和相关理论的完善，逐 步发展起来的水沙数值模拟计算方法，在河道及水库 的水流泥沙运动模拟研究中取得了较好的效果[2]。

2工程概况

天池抽水蓄能电站位于河南省南召县境内的黄 鸭河上游，由上水库、下水库、输水系统及发电厂房4

收稿日期1017-01-17;修回日期1017-04-09

基金项目：国家自然科学基金重点项目（51579013);水利部公益性科研专项(201201074)作者简介:孙羽（1977-)，女，河南郑州人，助教，硕士，主要从事水利工程建设与管理研究。E-mail:Sunyu@ncu.eduxn 通信作者:孙东坡（195-)，男，河南叶县人，教授，硕士，博士生导师，主要从事河流动力学研究。E-mail: sundongpo@ncwu.eduxn

6

长江科学院院报

2018 年

部分组成（图1 ( a))。抽水蓄能电站下水库坝址布置 在黄鸭河干流一个牛轭型(•〇)河湾的下游峡谷(胳驼 头）出口，大坝为钢筋混凝土面板堆石坝，过水建筑物 有溢洪道、排沙洞、电站进出水口（图1(b))。拦河坝 坝顶高程540.6 m，最大坝高100.6 m;下水库正常蓄水 位537.5 m，死水位510.0 m，调节库容1 库容431万

店水文站、相邻白河干流白土岗水文站以及附近焦

园、马市坪等雨量站的观测资料。

3研究内容与研究方法

针对抽水蓄能电站下库区地形、来水来沙特点

197万m3,死

m3,总库容1 784.66万m3。电站总装机

容量120万kW，电站运行为周调节;发电时进出水口

排水，蓄能时进出水口抽水。建成后库区将型河

及工程运用条件，建立下库区及相关工程的水沙数 学模型;利用已有水沙资料，进行模型的验证与参数 调试。在此基础上，采用不同系列年水文泥沙条件 湾的狭颈处破开，形成右侧输水输沙通道，左侧河湾

下部修拦沙潜堰(基本不过流）的格局（图1 (b))。

(a)抽水蓄能电站上、下库平面位置

(b)下库鸟瞰效果图(建筑物与输水通道)

图1天池抽水蓄能电站布置 Fig.1 Layout of Tianchi pumped

storage power station

黄鸭河为白河支流，上游为石山林区，山高坡 陡，植被良好，流域内受人类活动影响较小，水土保 持较好。坝址以上黄鸭河河道长19.0 km、坡降 35.1%C、控制流域面积112.9 km2;坝址处河谷较开 阔，呈“ U”型，谷底高程468 ~ 472 m，10 m河床宽度70~

，下库区有马蹄河与天池岈沟2个支流。据水 沙资料分析，黄鸭河洪水峰高量小，陡涨陡落，坝址 处多年平均径流量为4 447万m3;平均含沙量为 0.99 kg/m3,但洪水期最大含沙量达到25 kg/m3;洪

水期推移质泥沙与悬移质均有，坝址处悬移质中值 粒径0.018 mm,推移质中值粒径为0.83 mm,卵砾石

为3.2 mm，推悬比达到45% ~ 50%。由于只有2006 年5月在坝址下游1.5 km临时设站，观测过当年水 位、流量、含沙量外，附近河段无实测水文资料，因此

该河段的水文泥沙分析还参考了下游30 km处李青进行下库区水沙数值模拟，分析统计库区泥沙淤积 量、纵向淤积形态;分析其成因，探寻淤积特性与电 站运用的关系;提出改善库区泥沙淤积的措施。研究长系列年的河道冲淤演变以一维水沙数值

模拟为宜，效率高、模拟成果相对可靠[6]。MIKE11 河流模拟系统是适用于一维河道、河网综合水沙数 值模拟软件，主要用于河口、河流的水力和泥沙传输 模拟。软件系统的核心是采用六点隐格式求解一维 河流非恒定流方程的水动力模块（HD)，同时还应

用泥沙输运模块（ST)和控制构筑物模块（SO)，用 以MI K解E11决下边界的复杂工程边界[7]。本文运用

进行库区河道一维水动力泥沙模型的搭 建，利用水动力模块与泥沙输移模块对下库区水沙 运动进行数值模拟；然后运用MIKE VIEW对数值 模拟结果进行对比分析。

4数学模型的建立

4.1模型的理论基础4.1.1水动力学模块（HD)

(1)水动力学控制方程。一维水动力学控制方

程为Saint-Venant方程组，艮口

一3,4 3

+ dQ

二 qdx ;)3lQ ^ 3 Q2 ^ 4dh(1 dt dx A dx ^ QIQI ）

C2AR

=0式中：为距离坐标;t为时间坐标;A为过水断面面 积;Q为流量;h为水位；为旁测入流量;C为谢才 系数;R为水力半径；为重力加速度。

(2)方程组的离散及开边界处理。方程组离散 格式采用Abbott六点隐式格式，在每一个网格节点按 顺序交替计算水位和流量，因为该格式为无条件稳 定，可以在较大的Courant数下保持计算稳定，故时间 步长可取较长。河道上开边界可给出水位时间变化 过程h =h(〇，或流量时间变化过程Q = Q(〇 ;河道下 开边界上给出流量-水位关系Q = Q(h)[89]。

第1期孙羽等抽水蓄能电站下库长系列年泥沙淤积特性研究7

4.1.2泥沙输运模块（

ST模块涉及输沙率、河床变形计算，涉及糙率

ST)

的动态调整变化，根据选用的泥沙运动公式，模拟计 算河床的动态冲淤过程。

推移质模拟计算选用Meyer-Peter & Mullet•公 式，即

少=贫、[(ys - Y槡y/y _ 1幺以]

=8(0-0J15; (2)0= (YsT-〇 Y)D= (ys调-y)D= (yspUl-y)D

。

⑶

式中:少为推移质无量纲单宽输沙率;y，ys为水和 泥沙的重度;P为水的密度；为推移质泥沙粒径;0 为无因次水流切应力，由式（3)确定；0。为无因次 起动切应力；〖为以水下重量计的推移质单宽输沙 率;U t为摩阻流速；为水力坡降;D是悬移质泥沙 特征粒径。

悬移质模拟计算采用无量纲的单宽输沙率形 式， 即

式中^T，T分别为悬移质无量纲和有量纲的单宽 输沙率。

选用Engelund &Fredsoe的悬移质模拟计算公

式， 即

0T = 0.3/A • @V@2 + 〇.15。

（5)

式中A是水流的沿程阻力系数，由Darcy-Wei/ach 公式确定。

泥沙输运计算与水动力学计算是相互连接、连 续进行的，泥沙输移计算时可将床面阻力和冲淤变 形反馈给水动力模拟，并输出输沙率、床面变形及阻 力变化值结果[10-4]。

4.2下库区水沙模型的建立及验证 4.2.1

模型范围、断面剖分与开边界处理

根据研究任务，考虑研究河段河道地形特点、洪 水泥沙特性等影响因素，选取模型范围:上边界距拟 建下库坝址约2.6 km，下边界为下库坝址断面。模拟 河段平面布置及断面布设见图2。为满足流场计算 精度要求，模拟河段共划分39个断面（CS30断面以 下，汊河段分输沙槽与原〇河湾2个断面，CS39为坝 前断面）；根据河道地形特点，断面间距为60~100 m。 模型上游开边界根据模拟工况采用黄鸭河相应的流 量过程与含沙量过程，计算区域下边界取在下水库坝 址断面。下游开边界即河道出口断面条件，按照下水 库坝前水位变化过程曲线作为控制下开边界条件。

图2

模拟河段计算断面的布设

Fig.2 Arrangement of sections for numerical simulation

4.2.2

模型调试与验证

4.2.2.1模型调试

模型采用2009年黄鸭河汛期实测水沙过程及

水位观测资料，河道地形采用未建工程前的实际河 道地形(！河湾），对数学模型进行4种水流条件的 河道洪水演进及河道阻力状况的调试。模型在调试 过程中，计算河段分别采用不同糙率模拟流场阻力； 另外在模拟时，考虑随流量及河床冲淤变化，糙率也 适当调整;在调试过程中，根据实测河道断面图对局 部河段地形进行了对位修正，更真实地反映河道断 面形态对过流的影响，更合理地校验水位。

4.2.2.2验证成果

数学模型经调试，曼宁糙率系数一般取0.025〜

0.060。通过2009年黄鸭河汛期实测水沙过程的验

证计算，将计算河段各测点实测水位与模拟洪水流 场对应点计算水位进行校验对比。对比结果表明， 洪水水位观测点处的计算水位与实测水位值基本吻 合，水位最大误差为

0.09 m，水位平均误差为

0.05 m，见图3;验证结果表明，数学模型对河道阻

力状况（糙率）的模拟是比较合适的。通过模拟河 段洪水水位、大河主流及河势与实际调查情况对比， 也基本正确反映了洪水的演进特点、河势变化及河 道阻力特征，河道洪水演进模拟也是可信的。同时 河床变形的验证也采用物理模型实测冲淤资料进行 了对比分析，淤积总量和淤积部位都比较吻合。综 上所述，数学模型在河道断面处理、糙率选择、其它 控制参数的选择上是49.9

合适的，所建一维数学一实测水位

模型能正确模拟计算

4985

计算水位

区域内的洪水演进与 ：5 7

流量/(m9 3-s

11

_1)

冲淤状况，可以进行下 图3

模型河道阻力验证

水库坝址上游计算河Fi3 Verifiecti°n°friver

段的洪水演进及河床

rStnc fr the modd

8

长江科学院院报2018 年

变形模拟，预测建坝后下库区水沙输移及河床冲淤 变化，作为设计参考。

图7中也可以清晰看出，淤积主要集中于坝前0.5〜 1.8 km之间，占库区淤积量的绝大部分;这与同步进 行的物理模型试验成果是基本一致的。下水库的淤

积主要包括2个区域:一是马蹄河出口附近上下游 河段，二是新开的右支排沙槽及以上河段，淤积重心 5 —维水沙数值模拟及成果分析5.1水沙长系列的模拟条件

模拟采用10〜

50 a 5个系列水沙过程，这里给

出计算采用的概化水沙50 a长系列过程，见图4。

为了便于和物理模型对比，长系列年采用概化处理 的方法为:基于水沙长系列中的总沙量不变，将长系 列年水沙过程中<0.3 kg/m3含沙量的小流量过程从 时间序列中排除，并将排除的总沙量依然分配到仍 保留的水沙过程中；计算下边界条件采用系列年中 下库水位变化按设计典型周调节过程，见图5。

P2 400-流量一含沙量

30Z 21 600 000252 _冬 1 200 800 150£芑镔400

nj 0

i^m. 1111^：[jlrfll终 5〇 _ 总 l 200 ln1 500 1 800 2 100 2 400 2 700

0初

时间/h

图4模拟河段上边界采用的50 a长系列流量

及含沙量过程

Fig.4 Process lines of 50-year series discharge

and sediment concentration for the upstream

of the numerical region

540a日535 运 530 525 #515 520 510

505〇300 600

900 1 200 1 500 /h

1 800 2 100 2 400 2 700

时间图5

模拟河段下边界采用的坝前运用水位调节曲线

Fig.5 Process lines of water stage for the downstream of

the numerical region

5.2长系列水沙过程模拟成果及分析5.2.1 淤积级向分布特性分析

水库库区泥沙淤积纵向形态的发展直接体现了 库区超饱和不平衡输沙的沿程发展过程，反映了库 容被淤积侵占的过程及淤积物的纵向分布。为了研 究该抽水蓄能电站下水库泥沙淤积的纵向发展特 点，从模拟计算数据成果中提取了下水库河床各断 面深泓点随水沙系列年运行时间的变化数据，绘制 出相应的库区河床纵剖面如图6所示。从图6中可 以直观地看出，下水库的淤积整体呈现典型的三角 洲淤积形态，随着长系列年的运行时间增长，淤积面 逐年a抬高，三角洲洲头缓慢向坝前推进。第1个 10 中有含沙量较大的几次洪水，故淤积量较其他 年相对较大;其余年份的淤积速率基本持平。从

和淤积总量基本与物理模型试验成果一致。

图6

不同系列年水沙条件下模拟河段 沿程纵剖面变化曲线

Fig.6 Longitudinal profiles of the numerical region

under different flow-sediment conditions

图7给出了 40 a末模拟河段各计算断面深泓 点的变化曲线，与原始床面对比可以看出纵向淤积 形态，库区泥沙淤积主要集中于下水库中下段，即拦 河坝以上2.0 km范围内；相对应于CS10 (距坝 1.84 CSSkm)断面以上的库尾段受回水变动影响;上游

1—CS08(距坝2.5 ~ I.% km)断面之间河床则

表现为有冲有淤（图6)，大洪水时还略有冲刷，基本

保持动态平衡。

0

0.5 1.0

1.5 2.0 2.5 3.0

距现里程/km

图7

40 a水沙长系列条件下下库区河段

沿程深泓点高程变化曲线

Fig.7 Variation of bed stage along the numerical region

under 40-year series flow-sediment condition

根据模拟计算成果分析，从50 a系列末下库区 淤积的纵向分布(见图8)和库区淤积量的沿时程分 布(见图％可看出：模拟河段的淤积发展规律与下 库区的河床边界特性有关，同时也与水沙系列相互 响应。淤积纵向分布特性表现为：CS25断面（距坝 0.95 km)以上壅水影响弱，沿程淤积很小;马蹄河出 口受回流影响有局部淤积区，CS25—CS30 (距坝 0.59 km)河段沿程水深迅速增大，泥沙淤积相应明 显递增。距坝约0.6 km以下岔河中右支输沙槽河段 沿程淤积较多，是库区淤积的主体，占近坝区淤积总

第1期孙羽等抽水蓄能电站下库长系列年泥沙淤积特性研究9

量的80%;而左支n河湾河段拦沙效果较好，淤积 相对很少，仅占淤积总量的1%;电站进水口附近 只有少量淤积，见图8。淤积沿时程分布特性可以 从

一遇洪水淤积量就占到82%。本文还给出了 20 a 长系列过程中，20 a —遇洪水前后的河床地形对比， 见图1。此次洪水下库区泥沙淤积量约为18.6万t，

50 a系列几个特征河段淤积的时间发展过程清

晰看出：一方面在长系列年过程中的丰沙时期，库区 淤积强度就比较高，少沙时期的淤积强度就比较低; 另一方面河床边界条件的影响也很重要，上游河段 CS21 — CS01 (距坝1.17~2.55 km)有冲有淤，故累积

t,水库次洪水排沙比约

为11.3%。从图10中可以看到距坝里程2.1 km以上 河河段的河床有冲有淤，而次洪水在距坝0.7 ~ 1.6 km范围内形成河床强烈淤积，坝前也有局部强

烈淤积，这主要与洪峰期的泄流能力有一定关系。 频率洪水过程的计算淤积量及分布和动床模型试验 排出水库的泥沙量为2.1万

淤积量少，且主要是粗颗粒淤积（粒径约为0.4 ~

0.8 mm);距坝600~ 1 500 m的库段泥沙分选淤积 严重，累积淤积量比较多;坝前段CS39—CS30淤积 量反而没有CS30 — CS21段的淤积量大，主要是细 颗粒沉积（粒径约为0.018 ~ 0.025 mm)。不同河段 的淤积特性不同，库区淤积重心在CS30 — CS21库

段，这是有支流入汇与回水影响开始明显的库段，今 后防淤排沙的重点也在CS30—CS21库段。数值模 拟下库区沿程淤积分布成果和淤积物颗粒组成成果 与物理模型试验成果基本一致。

140 000 □ Q河湾

120 000 「」输沙槽(CS30~CS39断面） 100 000

口 CS30断面以上河道

ft

8

80 0060 000 4000 20 000

0

1 3 5 7 9 11 13 15 17 19 21 3 5 7

断面号

2222931 33 35 3739

图8

50 a系列下库区沿程淤积分布

Fi!g.8 Deposition distribution along the numerical region under 50-year series flow-sediment condition of the downstream reservoir

(§/w置

图9

下库区特征河段淤积量沿时程变化曲线（50 a系列)

Fig. 9 Process lines of deposition amount incharacteristic sections of the downstream reservoir

(50-year series)

5.2.2 频率洪水在长系列年水沙过程中的重要

作用

频率洪水在长系列年水沙过程中虽然时间不 长，但对库区河床冲淤有重要影响。不同频率的典 型洪水均在库区形成强烈的河床变形，典型洪水出 现频率越低，次洪水淤积量越大。次洪水淤积量很 大，构成库区淤积的主体部分:1 a —遇洪水淤积量 占1 a长系列淤积总量（平均值）的30%，而50 a

实测值基本一致[12]。

560

■原始河床

540 . —20 a—遇洪水后河床

_/

日

/sw5^m

/km

0 10 20 a _遇洪水前后模拟河段沿程深泓点高程变化曲线

Fig. 10 Curves of bed stage (before and after a 20-year

flood ) along the numerical region

5.2.3库区淤积总量及排沙比分析

水库泥沙淤积量的大小不仅取决于来水来沙， 也受制于水库运用包括坝前水位的控制与电站运行 状况；同时，泥沙淤积量及其分布也直接影响到取水 发电设施的正常运行。表1为下水库在不同运行年 份的泥沙淤积总量与排沙比。通过表中数据可以看 出，随着水库运行年限的增长，下水库库区泥沙淤积 量逐年增加，而排出水库沙量很少，下水库的排沙比 未超过14%，还未达到水库泥沙淤积的平衡。

表1

下水库泥沙淤积量及排沙比

Table 1 Deposition amount and sediment delivery

ratio in downstream reservoir

运行年限/a

来沙量/(万m3)

淤积量/(万

mi3)

排沙比/%10124.43118.444.8220165.54155.266.2130213.76197.807.4740246.65222.989.5950318.54

271.76

14.69

6结论

(1)根据抽水蓄能电站的复杂运用条件及水沙 特性，建立了基于MIKE11的某抽水蓄能电站下库水 沙数学模型，并通过了模型验证。进行了 10~50 a长 系列水沙过程及典型洪水水沙过程的数值模拟，获得 了下库区河段水沙输移及河床冲淤变化数据。

10

长江科学院院报2018 年

移数值模拟[J].长江科学院院报，212, 2(1) :-12. 衣秀勇，关春曼，果有娜，等.DHI MIKE FLOOD洪水 模拟技术应用与研究[M].北京：中国水利水电出版 社,2014.

(2 )长系列年水沙过程在下库库区形成的泥沙 淤积主要呈三角洲形态，逐年增加;泥沙淤积主要集 中于拦河坝上游0.5 ~ 1.8 km范围内，淤积重心在支 流马蹄河以下至〇河湾分流处;粗沙主要淤积在距

[6]

[7] RAHMAN M M,ARYA D S,GOEL N K,et al. Design

Flow and Stage Computations in the Teesta River， Ban­坝1.5 km以上河段，细沙主要淤积在坝前段;淤积分

gladesh， Using Frequency Analysis and MIKE 11 Model­布及淤积物沿程分选特性与物理模型试验成果基本

ing [J ]. Journal of Hydrologic Engineering， 2011，

^'致。16(2), d〇i：10.1061/(ASCE)HE.1943-5584.0000299.

(3) 库区淤积发展与水沙系列是相互响应的，[ 8] PALIWALR，PATRA R R.Applicability of MIKE21 to

Assess Temporal md Spatial Variation in Water Quality of 不同河段的淤积特性不同，主要与其所处环境的水

an Estuary under the Impact of Effluent from an Industrial 力与边界特性有关；沿纵向淤积形态存在距坝

Estate [J]. Water Science and Technology，2011，63(9)，

0.59~ 1.4 km库段的淤积重心，这将是通过水库运

doi： 10.2166/wst.2011.193 • Source： PubMed.

用、调度来防淤排沙的关键区域。[9] PANDA R K，PRAMANIK N，BALA B. Simulation of

River Stage Using Artificial Neural Network and MIKE 11 (4) 在长系列水沙过程中，下库区持续淤积，下

Hydrodynamic Model [ J ]. Computers & Geosciences，

水库的排沙比没有超过15%。频率洪水在长系列

2010, 36(6) ： 735-745.

年水沙过程中虽然占据时间不长，但在库区河床冲 [10] ACKER S P，WHITE W R.Sediment Transport： New Ap­淤量中占据很大份额，有重要影响。要提高排沙效 proach and Analysis[J]. Journal of the Hydraulics Divi­

sion， ASCE，1973, 99(11) ：2041-2060.率、减小淤积，需要进一步调整改善下水库的洪水调

[11] ENGELUND F, FREDSOE J. Sediment Ripple and Dune 度与电站运行。上述研究成果可供抽水蓄能电站工

[J]. Annual Reviewof Fluid Mechanics，1982，14:13-

程设计及运行方案确定参考。

37.

[2] 华北水利水电大学.天池抽水蓄能电站泥沙防治技术

参考文献：

研究报告[R].郑州:华北水利水电大学，2016.

卢金友，范北林，等.三峡水库运用后荆江典 [1] 钱宁，万兆惠.泥沙运动力学[M].北京:科学出版 [3] 董耀华，

型河段冲淤变化计算分析[J].长江科学院院报， 社，1983.

2005,22(2) ：9-12.[2] 谢鉴衡.河流模拟[M].北京:水利水电出版社，1990.

[3] 彭杨,张红武.三峡库区非恒定流一维水沙数值模拟[4] 黄仁勇，黄悦.三峡水库干支流河道一维非恒定水

沙数学模型初步研究[J].长江科学院院报，2009, [J].水动力研究与进展A辑,2006,21(3)285-292.

26 (2) ： 9-13.[4] 孙昭华，陈飞，郭小虎.长江下游近河口段一维水沙

数值模拟[J].水利水运工程学报，2007, (3) :44-50.[]黄仁勇，李飞,张细兵.三峡水库运用初期库区水沙输(编辑:姜小兰)

Numerical Simulation on Long Series Sedimentation in Downstream

Reservoir of a Pumped Storage Power Station

SUN Yu1, LI Xiang-yang2, SUN Dong-po1, ZHANG Xiao-lei1 ,ZHI Feng-juan3

(1.School of Water Conservancy，North China University of Water Resources and Electric Power,Zhengzhou 450046 , China; 2.E>ureau of Hydrology，Yellow River Conservancy Commission , Zhengzhou 450004, China;

3.School of Architectural Engneering，Kaili University ，Kaili 556011，China)

Abstract ： In order to understand the impact of a long series water and sediment procesvoir bed，we simulated the scouring and deposition and deformation in the downstream reservoir of Tianchi pumpedstorage power station using MIKE based on water level fluctuation in the downstream resof parameters and model validation, we obtained the long series total deposition amount and sediment delivery ratio of the reservoir as well as the deposition profile along stream and deposition process of characteristic reaches. Results suggest that the sedimentation of reservoir corresponds with the water and sediment process. Sedimentation mainlydisplays a delta shape，focusing in the range of 0.5-1.8 km upstream the dam.A

(下转第15页）

第1期

李舜等斜交塔基局部冲刷规律研究15

性，下一步将结合数学模型等研究手段进行不同交 角下局部冲刷研究工作。参考文献：

[1] RICHARDSON J, PANCHANG V. Three-dimensional Sim­

ulation of Scour Inducing Flow at Bridge Piers[J]. Journal of Hydraulic Engineering, 1998，124(5) ：53〇-540.

[2] D’ ALESSANDRO C. Effect of Blockage on Cylindrical

Bridge Pier Local Scour [ D ]. Windsor： University of Windsor, 2013.

[3] PARK C W, PARK H I，CHO Y K. Evaluation of the Ap­

plicability of Pier Local Scour Formulae Using Laboratory and Field Data [ J ]. Marine Georesources & Geotechnolo- gy，2017, 35(1): 1-7.

[4] ISMAEL A, GUNAL M, HUSSEIN H. Effect of Bridge

Pier Position on Scour Reduction According to Flow^ Direc-

tion [ J ]. Arabian Journal for Science & Engineering, 2015, 40(6) : 1579-1590.

[]齐梅兰.采沙河床桥墩冲刷研究[J].水利学报，2005,

36(7) ： 835-839.[6] 高正荣，黄建维，卢中一.长江河口跨江大桥桥墩局部

冲刷及防护研究[M].北京：海洋出版社，2005.[7] 韩海骞，熊绍隆.潮流作用下桥墩局部冲刷规律研

究[J].浙江水利科技，214, 195(5) : 87-91.

[8] 于洋.斜交桥墩对洪水壅高影响的数值研究[D].青

岛：中国海洋大学，2013.[9]吴 新生，韩向东，黄卫东，等.新型ADV流速仪在桥墩

模型试验中的应用[J].长江科学院院报，2011， 28(11) : 104-109.

(编辑:陈

敏

)

Rules of Local Scour of Skewed Tower Footing

LI Shun^CHAI Zhao-hui2 3 ,LIU Tong-huan2 3, FENG Yuan2 3

(1.Elnergy China East China Electric Power Design Institute Co. , Ltd. , Shanghai 200331,

2.River Department, Yangtze River Scientific Research Institute, Wuhan 430010, China;3.ICey Laboratory of River Regulation and Flood Control of Ministry of Water Resources,

Yangtze River Scientific Research Institute, Wuhan 430010, China)

Abstract:The local scour of tower footing (pier) is a critical subject in planning and designject. Skew^ed tow^er footing has being gradually used in river-crossing engineering due to and economic conditions. But researchers paid more attention to the local scour of orthograther than that of skewed tower footing. In this article, the rules of local scour of a skew^ed tow^er footing are re­searched through generalized model test. Results reveal that the increment of flow velocity on deviational side of skewed tower footing is bigger than that on dorsal side; the maximum scour depth is largdeviational side; the scour pit displays an asymmetrical horseshoe shape, and the scour range on deviational side is larger than that on dorsal side. When bed surface around the tow^er footing is protected, this trend becomes weak. The results offer reference for researches on the local scour of skew^ed tow^er footing of river-crossing bridge or cablegallery construction.Key words: skewed tower footing; pile group; scour pit; scour depth; flowvelocity

(上接第1

页）

condition, the deposition features and deposit compositions of different reaches vary correspondingly. Floods of cer­tain frequency has a significant influence on deposition amount despite that flood does not occupy much of the flow^-sediment series. Flood dispatching and operation of powder station must be adjusted and optimizment delivery ratio and reduce sedimentation. The research results are in agreement with those from scale modeltest，hence could offer technical parameters for the reservoir dispatching operation.

Key words:pumped storage power station; downstream reservoir; long series; numerical simulation; deposition;

sediment delivery ratio; focus of deposition