向家坝导流底孔回填混凝土温度动态预测.pdf

第 3 2卷 第 2期 2 0 1 5年 2月 长江科学院院报 J o u r n a l o f Y a n g t z e R i v e r S c i e n t i fi c R e s e a r c h I n s t i t u t e Vn 1 3 2 No 2 F e b 2 0 1 5 DO I ： 1 0 3 9 6 9 j i s s n 1 0 0 1 5 4 8 5 2 0 1 5 0 2 0 2 4 向家坝导流底孔回填混凝土温度动态预测 周建兵 。 黄耀英 ， 何小鹏 ， 田开平 ( 1 中国长江三峡集团 向家坝工程建设部， 四川 宜宾6 4 4 6 0 0 ； 2 三峡大学 水利与环境学院， 湖北 宜昌4 2 3 0 0 2 ) 摘要： 导流底孔回填混凝土水泥含量高， 温控难度大， 一般需要埋设冷却水管进行通水冷却。针对导流底孔回填混 凝土的特点， 从动态预测角度出发， 将有热源水管冷却计算式结合混凝土浇筑仓实测温度， 根据浇筑仓当前实测温度 动态更新有热源水管冷却计算式中的 ， 然后进行未来 n天混凝土浇筑仓温度动态预测， 以指导和调控现场通水措 施。向家坝导流底孔回填混凝土温度动态预测的实例分析表明， 建议的温度动态预测式是可行的。 关键词 ： 导流底孔； 回填混凝土； 初期冷却 ； 动态预测； 向家坝水电站 中圈分类号 ： T V 5 4 4 ； T V 4 3 1 文献标志码： A 文章编号： 1 0 0 1 5 4 8 5 ( 2 0 1 5 ) 0 20 l 1 9一 o 4 1 研究背景 导流底孔混凝土 回填一般施工工期 紧， 施工难 度大， 其是影响电站蓄水和发电工期的关键节 点 。缩短底孔 混凝土 回填 工期具有重 要 的工程意 义 。缩短工期可从优化施工工艺和优化温控措施 2 个方面进 行。 目前 ， 对施 工工 艺 的优化 报道 比较 多 J ， 而对温控措施 的优化研究很少。对温控措 施进行优化 ， 必须要进行快速、 准确的温度预测。对 于混凝土水管冷却 的温度计 算 ， 朱伯芳 假设等效 冷却混凝土柱体的外表面为绝热边界 ， 分别给 出了 无热源水管冷却计算式和有热源水管冷却计算式 。 针对中后期冷却期间混凝土浇筑仓顶面已覆盖浇筑 若干层混凝土浇筑仓 ， 而且浇 筑仓侧面也 已经粘贴 保温苯板 ， 为此 ， 黄耀 英等 建议 采用无热 源水管 冷却 问题的混凝土平均温度计算式 ， 联合实测温度 ， 对无热源水管冷却问题的混凝土平均温度计算式 中 的 i 进行动态更新 ， 消除高掺粉煤灰缓慢放热以及 上下游表面不是绝热边界等引起 的误差 ， 动态预测 混凝土浇筑仓的降温过程线。对 于一期冷却期间 ， 混凝土存在水泥水化温升、 以及混凝 土新浇筑仓在 间歇期间受外界环境气温影响很 大等 问题 ， 导致直 接采用有热源水管冷却计算式预测 的精度较差 。 导流底孔回填混凝土可归属于大体积混凝土范 畴， 其水泥含量高 ， 温控难度大， 一般需要埋设冷却 水管进行通水冷却。对于混凝土新浇筑仓来说， 3 个面或 5个面暴露在环境温度 中， 新浇筑仓 内部的 温度受环境 温度 影响很大 。与混凝土新浇 筑仓 不 同 ， 对于导流底孔 回填的混凝土来说 ， 底孔周 围为温 度相对稳定的大体积混凝土， 且底孔内环境气温也 相对稳定 ， 考虑到混凝土是热不 良导体 ， 底孔 内环境 气温对回填混凝土内部温度影响较小 ， 即可以采用 有热源水管冷却计算式 ， 联合实测温度 ， 较好地进行 导流底孔回填混凝土 的温度预测 。为此 ， 本文采用 初冷有热源水管冷却温度预测式， 结合向家坝导流 底孔 回填混凝土实测温度 ， 进行导流底孔 回填混凝 土的温度动态预测。 2 初冷期间实时动态预测模型原理 2 1 初冷期间有热源水管冷却问题 对于导流底孔 回填的混凝土来说 ， 由于混凝土 是热不 良导体 ， 可以假设导流底孔 回填 的混凝土处 于近似绝热状态。此时 ， 可采用初冷有热源水管冷 却温度计算式。混凝土浇筑仓 内埋设冷却水管进行 通水冷却 ， 设等效冷 却水管直径为 D， 长度为 ， 有 热源 ， 混凝土初温为 ， 进 口水温为 ， 则混凝土平 均温度可表示为 ： T ( t )= +( 一 ) ( t )+0 。 ( t ) ；( 1 ) ( t )=e 叩 ； ( 2 ) P ：d k a D ； ( 3 ) 收稿 日期 ： 加1 3一o 92 4； 修 回日期 ： 2 0 1 3一l 1 2 4 基金项 目： 国家 自然科学基金项 目( 5 1 2 0 9 1 2 4 ) 作者简介： 周建兵( 1 9 8 5 一) ， 男， 湖北武汉人， 工程师， 硕士， 主要从事水利工程项目 管理工作， ( 电话) 1 8 6 0 8 3 1 9 3 5 0 ( 电子信箱) z h o u j i a n b i n g 3 0 9 1 6 3 c o m。 1 2 0 长江科学院院报 = 2 0 91 3 5 +0 3 2 ； ( 4 ) ： ； ( 5 ) C P q d=1 9 4 7 ( b ) ； ( 6 ) = 0 9 2 6 e x p - 0 0 3 1 4 ( 告_ 2 0 ) 】 ； ( 7 ) D =2 b =20 5 8 3 6 。 ( 8 ) 式中： t 为冷却时间； a为导温系数； D， b ， C 分别为等 效冷却柱体的直径 、 外半 径、 内半径； A为混凝土 的 导热系数 ； ， s ：分别 为 水管 水平 和垂 直 向间距 ； ( t ) 为与混凝土绝热温升有关的函数 。 混凝土绝热温升表达式有指数型、 双曲线型和 复合指数型等 。由于指数型绝热温升便于进行数学 运算 ， 为此 ， 朱伯芳 院士 给出了指数型绝热温 升 0 ( r )=O o ( 1一e ) 下的 ( t ) 为 ( ) ： ( e 一e - ra t ) 。 ( 9 ) ， n P 进一步分析表明， 指数型绝热温升表达式与试 验资料吻合得较差 ， 为此 ， 本文采用组合指数式绝热 温升进行函数 缈( t ) 的推导。 设 已经 获得 混 凝 土 双 曲线 型 式 的绝 热 温 升 ( 丁 )： 二 ，为分析 问题方便 ， 将双 曲线 型绝 热温 升转换为组合指数型绝热温升 ， 即 0 ( r ) =O o s ( 1一e 一 )+O o ( 1一s ) ( 1一e 一 ) 。 ( 1 0 ) 式中： s = 0 6 ； m1 =1 4 5 n ； m2 = 0 1 4 5 n 。 设仅 已知混凝土绝热温升试验值 ， 则采用优化 的方法来确定组合指数式中的 m ， m ： ， s 。 采用组合指数式绝热温升导 出的函数 哕( t ) 为 ( )： ( 。 叫 一。 )+ ， 1一P ( 。 一。 )。 ( 1 1 ) m2 P 如果在一期冷却期间， 采用多档水温进行通水 冷却 ， 则混凝土平均温度为 T ( t )=T w +( 一 ) ( t )+0 o q t ( t )；( 1 2 ) ( t )=e 。 ( 1 3 ) 绝热温升采用组合指数式 ， 即 ( f )： ( e - p it - m l t i e ( t + t i ) + m1一 p1 ( 。 “一e - m 2 ( t + t i ) )。( 1 4 ) m 2 一 p 函数中的 为第 i 档通水温度 ； T i 为第 i 一1档 水温通水结束且第 i 档水温开始通水时的混凝土温 度 ； 为第 i 档水温通水时的水冷 函数 ； t 为改变水 温或流量时刻， 当改变水温或流量时时间 t 需要从 0 开始。 2 2 基于实测温度的初冷期间浇筑仓温度动态预 测模式 由 2 1 节分析可知 ， 有热源水管冷却计算公式 隐含了等效冷却直径为 JD的混凝 土柱体 的外表面 为绝热边界 。虽然底孔周围为温度相对稳定的大体 积混凝土， 且底孑 L 内环境气温也相对稳定， 但对于底 孔 回填混凝土新浇筑仓来说 ， 仍然一定程度受底孔 内环境气温影响， 这导致采用有热源水管冷却计算 式 的预测值与实测温度存在一定 的差异。为此 ， 从 动态预测角度出发 ， 将有热源水管冷却计算式 和混 凝土浇筑仓实测温度相结合 ， 根据浇筑仓当前实测 温度动态更新有热源水管冷却计算式 中的 ， 然后 进行未来 n天混凝土浇筑仓温度动态预测 ， 以指导 和调控现场通水措施 。有热源水管冷却问题的一期 冷却期间混凝土浇筑仓温度预测如图 1 。 ! 未来 第 天时 I !混凝土浇 筑实I ： 测温度信息 I 一 期冷却时混凝土浇 筑仓当前温度信息 未来n 天拟通水水温、 通水流量 水管间距 、水管材质和 混凝土热学性能 有热源水管 冷却计算式 未来n 天混凝土浇筑 仓温度信息 图 1 一期冷却期间混凝土浇筑仓温度动态预测 Fi g 1 Dyna mi c t e m p e r a t ur e f or e c a s t i ng f or c o nc r e t e c as t i ng bl o c k i n fir s t s t ag e c oo l i ng 由图 1可见 ， 基于浇筑仓当前实测温度动态更 新有热源水管冷却计算式中的 ， 可以将浇筑仓层 面散热边界等引起的误差 ， 通过动态更新 来动态 实时修正 ， 从而克服有 热源水管冷却计算式温度预 测效果不理想的问题。 3 导流底孑 L 回填 混凝 土绝热温升优 化反演 由2 1节可知 ， 混凝土绝热温升是初冷期问有 热源水管冷却问题 的重要组成部分 ， 为此 ， 本文采用 初冷有热源水管冷却温度预测式并结合实测温度对 混凝土绝热温升参数进行 了优化反演 。由式 ( 1 0 ) 可见组合指数式与双曲线式实际上是等价 的， 二者 的参数可 以相互转换 ， 即可用上述方法反演最终绝 热温升 和双曲线型绝热温升常数 n 。 第2期 周建兵 等 向家坝导流底孔回填混凝土温度动态预测 1 2 1 对相同空间位置和相 同时间上的温度实测值与 计算值进行比较， 使两者误差最小 ， 可建立 目标函数 n m F ( X )=( 一 Mf ) ， ( 1 5 ) ：l i =1 使 ( ) 一m i n 。 式 中： 为 部位 时刻 的计算温度值 ， 该计算 温度 由初冷有热源水管冷却温度 预测式计算 获得 ； M ， 为 i 部位 时刻 的实测温度值 ； m为观测点总数 ； n为某 部位观测时间序列 向量的个数。 由于 目标 函数是设计变量的非线性函数 ， 所以 该优化模型属于无 约束 的非线性规划 ， 采用单纯形 法进行反演计算。单纯形法 在 维空 间的几何 图 形中有 +1个顶点 ， 顶点 之间用直线连接 。再 对 各顶点的 目标 函数逐一进行比较 ， 不断舍弃最坏点， 代之 以使 目标 函数值不断降低 的新点。不断循 环、 迭代 ， 直到 目标 函数满足 给定条件时终止计算 。本 文采用 V i s u a l F o r t r a n编制 了单 纯形法热学参 数反 分析程序。 4 实例分析 向家坝水电站是金沙江下游梯级开发 中最末的 一 个梯级， 坝址位于川滇两省交界的金沙江下游河 段上 ， 拦河大坝为混凝土重力坝 ， 坝顶 高程 3 8 4 m， 最大坝高 1 6 2 m， 水库正 常蓄水位 3 8 0 m， 电站总装 机容量为 6 4 0万 k W。左岸非溢流坝段和冲沙孔坝 段共设 6个宽 高为 1 0 m1 4 m 的导流底孔 。选 取典型导流底孑 L 回填混凝土浇筑仓进行分析， 首先 基于初冷有热源水管冷却计算式 ， 结合实测温度 ， 采 用单纯形法进行导流底孔回填混凝土绝热温升优化 反演 ， 然后基于优化反演参数进行导流底孔 回填混 凝土温度动态预测 。 4 1 导流底孔回填混凝土绝热温升参数优化反演 典型仓混凝土位于导流底孔顶拱处， 采用的是泵 送 自密实混凝土 ， 强度等级为 C 2 5 ， 混凝土配合 比为 水泥2 8 4 k g m ， 水 1 7 0 k g m 。 ， 粉煤灰9 4 k g m 。 ， 砂率 5 0 。其水泥含量比其他部位高， 采用 D N 2 5铁管进 行通水 冷却 ， 水 管外 直 径为3 3 7 m m， 水 管 间距 为 1 m 1 m( 水平向 垂直向) ， 通水工况如表 1 。 表 1 通水工 况 Tab l e 1 Ca s e s o f wat e r c oo l i ng l o ad 基于初冷有热源水管冷却计算式 ， 结合实测温 度 ， 采用单纯形法进行导流底孔 回填混凝土绝热温 升优化反演 ， 反演结果为 0 ( )=3 1 0 8 ( 1一e 。 。 。 叶)+ 2 O 7 2 ( 1一e 邮 )。 ( 1 6 ) 根据文献 7 中给出的绝热温升估算公式估算 的最终绝热温升为 5 4 ， 由式 ( 1 0 ) 和式 ( 1 6 ) 可知 ， 导流底孑 L 回填 混凝 土优化反演 的最终 绝热温升为 5 1 8 C， 与组合指数型绝热温升等价的双曲线型绝 热温升常数 n为1 8 6 d 。 4 2 导流底孔回填混凝土温度动态预测 基于优化反演的绝热温升 ， 采用初冷有热源水管 冷却温度预测式 ， 按 图 1的动态预测模式 ， 对典型仓 0 2 8 d 龄期混凝土温度进行了动态预测。具体参见 通水工况表 1 ， 在通水水温和通水流量变化时的第 3 ， 1 1 ， 1 4 ， 2 1 d ， 结合实测温度对式 ( 1 2 ) 中的 进行动态 修正。预测温度与实测温度对比如图2所示。 毳 兰 图 2基于反演参数的预测值与实测值对 比 F i g 2 Co m p a r i s o n b e t we e n p r e d i c t e d t e m p e r a t u r e b a s e d o n par a me t e r i nve r s i on a nd me a s ur e d t e m p e r a t ur e 由图 2可见 ， 基于优化反演绝热温升的计算温 度与实测温度的吻合效果令人满意 ， 这一方 面说 明 绝热温升的反演结 果是合理 的， 另一方面说 明采用 初冷有热源水管冷却温度动态预测方法可以很好进 行导流底孑 L 回填混凝土的温度预测。 5 结 论 本文采用初冷有热源水管冷却温度预测式结合 向家坝导流底孑 L 回填混凝土典型仓实测温度进行了 绝热温升优化反演 ， 并基 于反演参数对典型仓 0 2 8 d龄期 混凝 土 温度 进行 了动态 预 测 ， 得 到如 下 结论 ： ( 1 )针对导流底孔 回填混凝土 的特点 ， 从动态 预测角度出发， 建议将有热源水管冷却计算和混凝 土浇筑仓实测温度相结合 ， 基于浇筑仓 当前实测温 度动态更新有热源水管冷却计算式 中的 ， 然后进 1 2 2 长江 科学 院院报 2 0 1 5芷 行未来 r , 天混凝土浇筑仓温度动态预测 ， 以指导和 调控现场通水措施。 ( 2 )根据典型仓混凝土实测温度 ， 运用初冷有 热源水管冷却温度预测和优化反分析方法 ， 对 自密 实混凝土绝热温升进行 了反演 ， 优化反演结果与估 算结果接近且计算值与实测值吻合效果很好， 即绝 热温升的反演结果是合理的。 ( 3 )基于优化反演参数对典型仓混凝土 0 2 8 d 龄期进行了动态预测 ， 预测值与实测值吻合得较好 ， 即采用初冷有热源水管冷却温度动态预测方法 ， 可以 很好地进行导流底孔回填混凝土的温度预测。 参考文献 ： 1 2 3 4 何小 鹏 ， 许 光 明， 王 忠 民， 等 三 峡工程 导流 底孔 封堵 生产性试验研究 J 人民长江， 2 0 0 6 ， 3 7 ( 5 ) ： 1 51 8 ( H E X i a o p e n g ，X U G u a n g - m i n g ，WA N G Z h o n g m i n g ， e t a 1 Fu l l s c a l e Ex p e r i me n t a l S t u d y o f Pl ug g i n g Con s t r u e t i o n o f B o t t o m D i v e r s i o n O u t l e t o f T G P J Y a n g t z e R i v e r ， 2 0 0 6 ， 3 7 ( 5 ) ： 1 51 8 ( i n C h i n e s e ) ) 杨秀彬 小湾水 电站导流底孔 封堵混凝 土施 工 J 水 利水 电 施 工 ，2 0 1 2 ， ( 5 ) ： l 51 7 ( Y A N G X i u b i n P l u g g i n g Co n s t r u c t i o n o f Bo t t o m Di ve r s i o n Out l e t o f Xi a o w a n H y d r o p o w e r S t a t i o n J C o n s t ruc t i o n o f H y d r a u l i c E n g i n e e r i n g ， 2 0 1 2 ，( 5 ) ： 1 51 7 ( i n C h i n e s e ) ) 林毅辉 蒙州水电站拦河坝施工导流设计及导流底孔 封堵施 工 方法 的探 讨 J 小 水 电 ， 2 0 0 5 ，( 5 ) ： 6 6 6 7 f LI N Yi hu i Di s c u s s i o n o f Co ns t ru c t i o n Di v e r s i o n De s i g n a nd Pl ug g i n g Co ns t ruc t i o n Me t h o d o f Bo t t o m Di v e r s i o n O u t l e t o f Me n g z h o u H y d r o p o w e r S t a t i o n D a m J S ma l l Hy d r o P o w e r 2 0 0 5( 5 )： 6 66 7 ( i n C h i n e s e) ) 肖汉江， 崔建华， 李祖 民， 等 三峡大坝导流底孔封堵 温度应力仿 真计算 分析 J 长 江科学 院院报 ， 2 0 0 2 ， 1 9 ( 6 ) ： 1 7 2 0 ( X I A O H a n - j i a n g ， C U I J i a n h u a ， L I Z u 一 5 6 7 8 9 mi n g，e t a 1 Th e r ma l St r e s s Emu l a t i o n a l Ca l c ul a t i o n o n B l o c k a g e B o d y o f B o t t o m O u t l e t o f S p i l l wa y Da m S e c t i o n f o r T G P 1 J【 J o u r n a l o f Y a n g t z e R i v e r S c i e n t i fi c R e s e a r c h I n s t i t u t e ， 2 0 0 2 ，1 9 ( 6 ) ： 1 7 2 0 ( i n C h i n e s e ) ) 王 尚海， 陈剑 红 隔河岩大 坝导流底孑 L 封堵 的施 工 J 湖北水力发电，2 0 0 3 ，( 4 ) ： 1 92 1 ( WA N G S h a n g h a i ，C HE N J i a n h o n g P l u g g i n g C o n s t ruc t i o n o f B o t t o m Di v e r s i o n Ou t l e t o f Ge h e v a n Da m -Hu b e i Wa t e r P o we r ， 2 0 0 3，( 4)： 1 92 1 ( i n Ch i n e s e ) ) 李 然 小 湾 水 电工 程 坝 身 导 流底 孑 L 封 堵 施 工 技 术 J 水利水 电技术 ， 2 0 1 2 ， 4 3 ( 3 ) ： 4 2 4 6 ( L I R a n Co ns t ruc t i o n Te c hn o l o g y f o r P l u g g i ng o f Di v e r s i o n Bo t t o m Ou t l e t o f Da m f o r Xi a o wa n Hy d r o p o w e r P r o j e c t J1 Wa t e r Res o u r c e s a n d Hy d r o p o we r Eng i n e e r i ng， 201 2， 4 3( 3) ： 4 24 6 ( i n Ch i n e s e ) ) 朱伯芳 大体积混凝土温度应力与温度控制 M 北 京 ： 中 国 电 力 出 版 社 ， 1 9 9 9 ( Z H U B o f a n g T h e r m a l S t r e s s e s a n d T e mp e r a t u r e C o n t r o l o f Ma s s C o n c r e t e M B e i j i n g ：C h i n a E l e c t r i c P o w e r P r e s s 1 9 9 9 ( i n C h i n e s e ) ) 黄耀英， 周绍武， 付学奎， 等 中后期冷却期间混凝_ 浇筑仓温度动态预测模型 J 四川大学学报( 工程科 学 版 ) ，2 0 1 3 ， 4 5( 4 ) ： 3 43 8 ( H U A N G Y a o y i n g ， ZHOU S ha o wu， F U Xu e k u i ，e l a 1 Dy n a mi c Fo r e e a s t i n g Mo d e l o f Co n c r e t e Ca s t i n g Wa r e h o us e i n Mi d d l e k a t e C o o l i n g A g e J J o u r n a l o f S i c h u a n U n i v e r s i t y( E n g i n e e r i n g S c i e n c e E d i t i o n ) ，2 0 1 3 ，4 5( 4 ) ： 3 43 8 ( i n C h i n e s e ) ) 朱伯芳 混凝 土热 学力 学性能 随龄期变 化的组合 指数 公式 J 水利学报 ， 2 0 1 1 ， 4 2 ( 1 ) ： 1 7 ( Z HU B o f a n g Co mpo u nd Ex po n e n t i a l Fo r mu l a for Va r i a t i on o f Th e r mal a n d Me c h a n i c a l P r o p e r t i e s w i t h Ag e o f C o n c r e t e J 1 J o u r n a l o f Hy d r a u l i c E n g i n e e r i n g ，2 0 1 1 ，4 2( 1 )： 17 ( i n Ch i n e s e ) ) ( 编辑 ： 曾小汉) Dy n a mi c Pr e d i c t i o n o f Co n c r e t e Te mp e r a t ur e f o r t h e Pl u g g i n g o f B o t t o m Di v e r s i o n Ou t l e t o f Xi a n g j i a b a Hy d r o p o we r P r o j e c t Z HOU J i a n b i n g ，HU ANG Ya o y i n g ，HE Xi a o p e n g ，T I AN K a i p i n g ( 1 X i a n g j i a b a P r o j e c t C o n s t r u c t i o n D e p a r t m e n t ， T h r e e G o r g e s C o r p o r a t i o n ， Y i b i n 6 4 4 6 0 0 ， C h i n a ； 2 C o l l e g e o f H y d r a u l i c a n d E n v i r o n m e n t a l E n g i n e e r i n g ， T h r e e G o r g e s U n i v e r s i t y ， Y i c h a n g 4 4 3 0 0 2 ， C h i n a ) Ab s t r a c t ： Th e c o n c r e t e p l u g g i n g o f b o t t o m di v e r s i o n o u t l e t h a s hi g h p r o p o r t i o n o f c e me n t ，p o s i n g b i g c h a l l e n g e t o t h e t e mp e r a t u r e c o n t r o l ，t h e r e f o r e p i p e c o o l i n g i s ne e d e dI n t hi s r e s e a r c h，we r e c o mme n d t o c o mb i n e t h e h e a t p i p e c o o l i n g c a l c u l a t i o n f o r mu l a wi t h t h e me a s ur e d t e mp e r a t u r e o f c o n c r e t e b l o c k，whi c h me a n s t o d y n a mi c a l l y u p d a t e T i n t h e h e a t pi p e c o o l i n g c a l c u l a t i o n f o rm u l a a c c o r d i ng t o t he me a s u r e d t e mpe r a t u r e，a n d t h e n p r e d i c t t h e d y n a mi c t e mp e r a t u r e o f c o n c r e t e b l o c k o f n d a y s i n t he f u t ur e t o g u i d e a n d c o n t r o l t h e p i p e c o o l i n g Th e d y n a mi c p r e d i c t i o n f o r c o n c r e t e t e mp e r a t u r e i s p r o v e d t o b e f e a s i b l e b y t h e p r a c t i c a l e x a mp l e o f c o n c r e t e t e mpe r a t u r e o f p l u g g i n g o f b o t t o m d i v e r s i o n o u t l e t a t X i a n g j i a b a H y d r o p o w e r P r o j e c t K e y wo r d s ： b o t t o m d i v e r s i o n o u t l e t ； p l u g g i n g c o n c r e t e ；i n i t i a l p i p e c o o l i n g ； d y n a m i c p r e d i c t i o n ； X i a n g j i a b a h y d r o p o we r s t a t i o n