龙华口碾压混凝土重力坝表面永久保温对施工期温度应力的影响.pdf

第 8 卷第 5 期 2 0 1 0 年 1 O月 南水 北调与水利 科技 S o u t h - t o - N o r t h Wa t e r T r a n s f e r s a n d Wa t e r S c i e n c e T e c h n o l o g y Vo 1 8 No 5 0c t 2 0 1 0 d o i ： 1 0 3 7 2 4 S P d 1 2 0 1 2 0 1 0 0 5 0 5 2 龙 华 口碾压 混凝 土重力坝 表面永久保温对施工期温度应力的影响 金鑫， 李鹏辉 ， 薛 城 ， 孙国强 ( 清华大学 水沙科学与水利水电工程国家重点实验室， 北京 1 0 0 0 8 4 ) 摘要： 龙华口碾压混凝土重力坝位于寒冷地区， 冬季坝体内外温差大， 易产生较大的拉应力而导致表面开裂。现结合 大坝工程实际施工条件及进度安排， 通过三维有限元仿真分析， 主要研究了浇筑进度和表面永久保温措施对施工期大 坝混凝土温度应力的影响。仿真分析结果表明， 施工进度安排对大坝温度应力有较大影响， 采取大坝表面永久保温措 施可显著减少坝体内外温差， 有效控制大坝温度应力， 为施工期温控设计提供了参考依据。 关键词： R C C重力坝； 表面保温； 温度应力； 仿真分析 中图分类号： T V6 4 2 2 ； TV5 4 4 9 2 1 文献标识码： A 文章编号： 1 6 7 2 1 6 8 3 ( 2 0 1 0 ) 0 5 0 0 5 2 0 6 I n flu e n c e o f Pe r ma ne nt Su p e r fic i a l He a t Pr e s e r v a t i o n o n t he Te mpe r a t u r e S t r e s s d u r i n g t h e C o n s t r u c t i o n Pe r i o d o f L o n g h u a k o u RCC Gr a v i t y Da m J I N Xi n ， L I P e n g - h u i ，XUE C h e n g， S UN Gu o - q i a n g ( S t a t e Ke y L a b o r a t o r y o f Hy d r o s c i e n c e a n d En g i n e e r i n g，Ts i n g h u a Un i v e r s i t y， Be i j i n g，1 0 0 0 8 4 ，C h i n a ) Ab s t r a c t ：Lo n g hu a k o u RCC g r a v i t y d a m i s l o c a t e d a t a c o l d r e g i o n I n wi n t e r t he r e i s a l l o bs e r v a b l e t e mp e r a t u r e d i f f e r e n c e b e t we e n i n s i d e a n d o u t s i d e o f t he d a m，wh i c h ma y i n du c e l a r g e r t e ns i l e s t r e s s t r ig g e r i n g t he s u r f a c e c r a c k i n g Th e i n f l u e n c e o f v a r i o u s s c h e d ul e s a n d p e r ma n e n t s u p e r f i c i a l he a t p r e s e r v a t i o n o n t he t e mp e r a t u r e s t r e s s d u r i n g t h e c o ns t r u c t i o n p e r i o d i s s t u d i e d us i n g t h e 3 D f i n i t e e l e me n t s i mu l a t i o n a n a l y s i s ba s e d o n t h e a c t u a l c o nd i t i o n s Th e r e s u l t s s h o w t ha t v a r i o u s s c h e d ul e s h a v e g r e a t i n f l u e nc e O il t h e t e mp e r a t u r e s t r e s s o f t h e d a mTh e p e r ma n e n t s u p e r f i c i a l he a t p r e s e r v a t i o n me a s u r e s c a n s i g n i f i c a n t l y d e c r e a s e t h e t e mpe r a t u r e dif f e r e n c e a n d e f f e c t i v e l y c o n t r o l t h e t e mp e r a t u r e s t r e s s o f t h e d a m wh i c h c a n p r o v i d e t h e i mp o r t a n t r e f e r e nc e f o r t h e t e mp e r a t u r e c o n t r o l d e s ig n o f t h e d a m d u r i n g t h e c o n s t r u c t i o n p e r i o d Ke y wo r d s ：RCC gr a v i t y d a m； s u p e r f i c i a l h e a t p r e s e r v a t i o n； t e mp e r a t u r e s t r e s s ； s i mu l a t i o n a n a l y s i s 龙华 口水 电站 工程 位于 山西 省 盂县 ， 主坝 为 碾压 混 凝 土( RC C) 重力坝 ， 最大坝高 6 3 2 m， 坝顶长 3 5 3 I T I 。坝体 基本断面为三角形 ， 坝顶高程 5 7 5 I n ， 宽 5 m。上游面 自高 程 5 2 5 In处以 1： 0 2坡度放坡至建基面； 下游面坝坡 1： 0 7 5 ， 最大底宽 5 0 1 T I 。大坝共分 1 4个坝段， 包括 2个表孔 溢流坝段和 2个底孔泄流坝段。坝体不设纵缝， 横缝宽度为 1 5 30 In 。 为了清楚地了解龙华 口碾压混凝土重力坝的温度特性及 温度应力 ， 本文选取 6号典型坝段为研究对象， 使用清华大学 水利系自主开发的三维有限元仿真计算程序 S a p t s V1 0模 拟了施工期全过程的温度和徐变应力变化( 含累积 自重) ， 分 析了坝体上、 下游表面采取永久保温措施的效果， 为碾压混凝 土重力坝的温度控制及防裂设计提供依据【 1 。 l 计算模型及参数 1 1 计算模型 选取 6号坝段作为典型坝段进行建模分析， 地基深度方 向取 1 0 0 r n ( 约 1 5 倍坝高) ， 上、 下游方向分别取 1 0 0 r r l 。坝 基边界约束条件为： 基岩上下游计算边界给定顺河向约束 ， 其 他两向自由； 基岩左右边计算边界给定沿坝轴线方向的约束， 其他两向自由； 基岩底部计算边界面给予固定约束， 即z、 Y 、 三向约束。计算坐标系符合右手法则， 轴为横河向( 指向右 岸) ， y轴指向上游( 顺水流方向) ， z 轴向上( 表示高程) 。 6 号典型坝段( 含地基) 共离散为 1 2 9 8 0 个单元 ， 其中 坝体 1 1 1 6 0 个单元 ， 地基 1 8 2 0 个单元， 均采用 8节点 6 面体等参单元。坝体和地基节点一共 1 5 3 3 4个 。单元的 收稿 日期 ： 2 0 1 0 0 7 3 0 修回 日期： 2 0 1 0 0 8 2 6 作者简介 ： 金鑫( 1 9 7 6 一 ) ， 男 ， 河南信阳人 ， 助理研究员 ， 硕士 ， 主要从 事水利工程坝工结构设计及科研工作 。 强 * 强 a” |一 豫 5 2 试 验研究 金 鑫等 龙华 口碾压混凝土重力坝表面永久保温对施工期温度应力的影响 大小依据施工进度和层厚来确定。 分 区见 图 2 。 计算网格见图 1 ， 材料 1 3 计算边界条件 图 1 典型坝段有限元计算网格 Fi g 1 Fi n i t e e l e me nt me s h o f t y p i c a l d a m s e c t ion C 1 5常态 混凝土垫层 图 2 典型坝段材料分 区 Fi g 2 Ma t e r i a l p a r t i t io n o f t y pic a l d a m s e c t i o n 1 2 计算参数 根据资料统计， 坝址多年平均气温为 9 6， 逐月多年平 均气温见表 1 。 表 1 坝址逐月多年平均气温 Ta b l e 1 M o nt h l y a nn u a l a v e r a g e t e mp e r a t ur e o f t h e d a m s i t e 材料的热力学参数见表 2 。 表 2 材料的热力学参数 Ta b l e 2 Th e r mo d y n a mi c p a r a me t e r s O f ma t e r i a l s 在计算中， 材料的热学、 力学参数因分区的不同而采用相 应的参数。施工期各坝段的上、 下游面及施工浇筑面的气温 值参见 表 1 ， 坝段 的左 、 右 面取 为绝热面 。 坝体上、 下游表面按喷涂 5 c m厚聚氨酯闭孔保温材 料_ 3 进行永久保温， 折算后的等效表面散热系数选取为 0 5 W ( m K) 。 1 4坝 体应 力控 制 指标 根据混凝土重力 坝设计 规范 l_ 5 ， 混凝 土 的抗 裂 应力 控制 标准按照极限拉伸值法确定。结合相关工程经验， 应力控制 标准取值见表 3 。 表 3 混凝土容许拉 应力 Ta bl e 3 Al l o wa b l e t e n s i l e s t r e s s o f c on c r e t e s 2 计算分析工况 2 1 不 同施 工控 制进 度 安排 本工程坝体混凝土总方量约为 4 1 万 r n 3 ， 其中碾压混凝土 约 3 9 万 m3 。2 0 0 9 年 8月开始基 础垫层混凝 土浇筑 ， 1 2 月 一2 月气温在零度以下时停浇， 不采取大的制冷和加热系统， 夏季 高温季节采取简易温控措施控制混凝土浇筑温度不超过 2 3 ， 其余时间浇筑温度为月平均温度+2辐射热。 进度安排： 根据工程计划， 工程施工 总工期 2 5年， 2 0 1 0年 6 月浇筑至坝顶。平均每个月浇筑 1 0 r n 。 进度安排 ： 结 合实 际施 工能 力与 施工 进度 ，2 0 1 0年 8 月浇筑至坝顶 。平均 每个月浇筑 8 m。 进度安排： 2 0 1 0年 1 O月浇筑到坝顶。平均每个月浇 筑 6 m。 进度安排 ： 2 0 1 0 年 l 1 月浇筑 到坝 顶 ， 夏季停 工 。平均 每个月浇筑 8 m。 2 2不同表面永久保温方案 河床( 5 3 0 m高程) 以下强约束区的大坝混凝土上下游面 采取表面贴苯板然后覆土回填等工程措施 ， 保证冬季大坝上 下游面 的温度不 小于 1 O。而对 于河床 以上混 凝土 表面是 否采取保温措施 ， 本文分两种方案进行对比研究。 方案 1 ： 河床以上不采用保温措施 。 方案 2 ： 河床以上采用永久保温材料保温 6 。 。 结合上述不同施工进度安排 ， 本文共分析了 8种不同的 工况 ， 见表 4 。 表 4 仿 真分析工况 Ta b l e 4 S i mu l a t io n a n a l y s i s c a s e s 试 验 撅 尧 第 8 卷 总第5 O期 南水北调与水利科技 2 0 1 0年第 5 期 3 施工期温度应力仿 真结果与分析 施工期坝体温度应力随时间变化， 而工程中比较关心在 整个施工期大坝可能出现的最大拉应力。由于冬季外界气温 较低， 混凝土大坝内外温差较大， 此时在混凝土表面附近会产 生较大的拉应力 。 3 1 进度下工况 1 、 2的成果分析 工况 1 为河床以上不采用保温措施， 重力坝横断面中心 点温度沿高程过程线见图 3 ； 工况 2为河床以上采用永久表 面保温措施 ， 重力坝横断面中心点温度沿高程过程线见图 4 。 其中温度场单位为。 一 一 赠 5 0 1 0 一一 一 。 一 一 。 一 一 一 ， ： 0 9 7 0 9 1 0 l 0 一l 1 0 5 l ( 1 - 8 l O 一1 l 1 l 一 2 l l 一6 年 月 图 3 工况 1中心点沿高程温度过程线 Fig 3 Te mp e r a t ur e h y d r o g r a g h o f c e n t e r a l o n g t he e l e v a t i o n of 1 s t c a s e 年一 月 图 4 工况 2中心点沿高程温度过程线 Fi g 4 Te mp e r a t u r e hy dr o gr a g h o f c e n t e r a l o n g t h e e l e v a t ion o f 2 n d c a s e 对比可以发现： 采用永久保温后， 河床以上混凝土内部最 高温度比无保温措施时高约 3， 但是其内部温度变幅将明 显减小， 表面温度能保持在 1 0左右， 对降低混凝土表面拉 应力有利 。 工况 1的大坝轴向最大应力包络图和横断面冬季 1月坝 轴向应力见图 5和图 6 ； 工况 2的大坝轴向最大应力包络图 和横断面冬季 1 月坝轴向应力见图 7和图 8 。最大应力包络 图给出了坝段 中线横 剖 面上 距 上游 面 0 5 m( u p ) 、 正 中心 ( c e n t r e ) 、 距下游面0 5 m( d o w n ) 处在施工期的最大拉应力沿 坝高分布情况。其中应力场单位为 MP a 。 在进度安排下 ， 大 坝最大拉应力 出现在 5 4 5 1T I 高程 附 近的坝上下游表面处， 出现 的时间在 2 0 1 1年 1月。工况 1 没有保温措施时最大拉应力值达到 1 4 MP a 左右， 而工况 2 采用永久保温措施后最大拉应力值在 1 0 MP a以下。另 外 ， 在基础强约束， 大坝处于高应力状态， 尽管在大坝浇筑 过程中， 垫层出现微裂缝 ， 释放了部分拉应力 ， 但该部位仍值 得重点关注 。 5 哉 泰 麦 惶 - 0 2 Fi g 5 0 0 2 0 4 0 6 0 8 l 1 2 1 4 l 6 l 8 应 力， MP a 图 5 工况 1 坝轴向最大应力包络图 Ma x imu m a xia l s t r e s s e n v e l o p d i a g r a m o f 1 s t c a s e 图 6 工况 1 横断面冬季 1月坝轴向应力图 F i g 6 Ax i a l s t r e s s o f c r o s s - s e c t i o n i n J a n u a r y o f 1 s t c a s e 成 力 MP a 图 7 工况 2 坝轴 向最 大应力包络 图 Fi g 7 Ma x i mu m a x i a 1 s t r e s s e nv e l o p d ia g r a m o f 2 n d e a s e 图 8 工况 2横断面冬季 1月坝轴 向应 力图 Fi g 8 Ax i a l s t r e s s o f c r o s s - s e c t i o n i n J a n u a r y o f 2 n d c a s e 渤 渤 渤 册 湖 啪 渤 u v 龌 啪 跏 渤 第8卷 总第 5 O期 南水北调与水利科技 2 0 1 0年第5期 4 5 5 0 ；3 0 5 0 0 9 7 O 9 1 0 l 0 l l 0 5 1 0 8 1 0 一l 1 l 1 2 l 1 6 1 1 9 年一 月 图 1 5 工况 6中心点沿高程温度过程线 Fig 1 5 T e n p e r a t u r e h y g r a g h o f c e n t e r a l o n gt h ede v a fi o n o f 6 t h c a s e 幄 应 J MP a 图 1 6 工况 6坝轴向最大应力包络图 Fi g 1 6 Ma x i mu m a x i a l s t r e s s e n v e l o p d i a g r a m o f 6 t h c a s e 5 O 4 5 4 ( ) 一 3 5 p 3 0 2 5 籍20 1 5 l 0 5 O 0 9 7 0 9 一J 0 1 0 一l i O 一 5 1 0 8 l 0 1 】 l l 一 2 l 1 6 1 l 一 9 年一 月 图 1 7 工况 7中心点沿高程温度过程线 Fi g 1 7 Te mp e r a t u r e hy d r o g r a g h o f c e n t e r a l o n g t h e e l e v a t i o n o f 7 t h c a s e 图 1 8 工况 7 坝轴向最大应力包络图 F i g 1 8 Ma x i mu m a x i a l s t r e s s e n v e l o p d i a g r a m o f 7 t h c a s e 年一 月 图 1 9 工况 8中心点沿高程温度过程线 Fi g t 9 Te mp e r a t u r e hy d r o g r a g h o f c e n t e r a l o n g t he e l e v a t i o n o f 8 t h c a s e 哪 豫鬻 弹辨 | 警t 豫 5 6 试 验研究 鑫 姬 应力 M P a 图 2 O 工况 8坝轴向最大应力包络图 Fi g 2 0 Ma x i mu m a x i a l s t r e s s e n v e l o p d i a g r a m o f 8 t h c a s e 3 5 不 同工况温度 应 力成果 总结 不同进度不同工况下河床以上大坝施工期最大轴向拉应 力值见表 5 。 表 5 河床以上坝体施工期最大轴向拉应力值表 Ta bl e 5 Ma x i mu m a x i a 1 s t r e s s d u r i n g t h e c o n s t r u c t i o n p e r i o d o f t h e d a m a b o v e t h e r i v e r b e d 可见施工进度安排对大坝温度应力有较大影响， 在条件 允许的前提下 ， 尽量保证大坝在夏季来临前浇筑完成； 另外也 可利用低温季节适当加快施工进度， 而在夏季应尽量避免高 温浇筑。在河床高程以上的坝体上、 下游表面采用聚氨酯闭 孑 L 保温材料进行永久保温后， 坝面受外界温度影响减小， 可以 大大削减施工期和运行期的温差应力场。 4 结语 本文通过对龙华口碾压混凝土重力坝典型坝段在不同进 度安排下不同工况的温度场和应力场进行仿真分析， 研究了 不 同浇筑进度对大坝温度应力的影响 以及大坝上下游混凝土 表面采取永久保温措施的必要性。 施工进度安排对大坝温度应力有较大影响， 在条件允 许的前提下， 应尽量利用低温季节适当加快施工进度， 而避免 在夏季高温浇筑。 河床以上混凝土表面采用永久保温后 ， 其内部最高温 度将增加约 3， 但 由于其表面温度能保持在 1 O左右以 上， 混凝土内部及表面温度变幅将明显减小， 对降低混凝土表 面拉应力有利。 河床以上若无表面保温措施， 各工况下大坝在施工完 后的第一个冬季( 2 0 1 1 年 1 月) 均会产生较大的表面拉应力， 可能会引起坝体裂缝。 河床以上采取永久表面保温措施后， 坝体最大拉应力 可减小 0 4 O 8 MP a ， 均在坝体混凝土允许拉应力范围内， 大大改善了坝体应力条件。 参考文献 ： E 1 3 刘光廷 ， 胡昱 ， 王恩志 ， 等 石门子碾压混凝土拱坝温度场实测与 仿真计算 J 清华大学学报( 自然科学版) ， 2 0 0 2 ， 4 2 ( 4 ) ， 5 3 9 5 4 2 ( I I U Gu a n g - t i n g， HU Yu， WANG En - z h i ， e t a 1 An a l y s i s a n d Me a s u r e me n t o f t h e Te mp e r a t u r e Fi e l d i n t h e S h i me n z i 跏 m 姗 金 鑫等 龙华口碾压混凝土重力坝表面永久保温对施工期温度应力的影响 R C C Ar c h Da m J Ts i n g h u a Un i v e r s i t y ， 2 0 0 2 ， 4 2 ( 4 ) ， 5 3 9 5 4 2 ( i n Ch i n e s e ) ) 2 朱伯芳 大体积混凝土温度应力 与温度 控制 M 北 京 ： 中国 电 力出版社 ， 1 9 9 9 ( Z HU B o - f a n g ， T e mp e r a t u r e S t r e s s a n d C o n t r o i s o f Ma s s C o n c r e t e M 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 Ch i n e s e ) ) 1- 3 杜彬 聚氨酯硬质泡沫在大坝工程中的应用研究r -J 水利水电科 技进展 ， 2 0 0 2 ， 2 2 ( 4 ) ， 1 4 1 6 ( D u B i r L Ap p l i c a t io n o f P o l y u r e t h a n e R d f o a m t o d a m P mj e c t s - J A d v a n c e s i n S c i e n c e a n d T e c h n o l o g y o f Wa t e r Re s o u r c e s ， 2 0 0 2， 2 2 ( 4 ) ， 1 4 1 6 ( in Ch ine s e ) ) 4 郭继师 ， 寇福立 ， 温贵 明， 等 聚氨酯保温技术在北疆某 水电站的 研究与应用 J 水利水电工程设计， 2 0 0 8 ， 2 7 ( 3 ) ， 4 6 4 8 ( G U O J i s h i ， KO U F u l i ， W E N Gu i mi n g ， e t a 1 Th e Re s e a r c h a n d Ap p l i c a t i o n f o r Po l y me t h a ne He a t Pr e s e r v a t i o n Te c h n o l o g y U s e d i n O n e H y d r o p o w e r S t a t io n i n N o r t h P a r t o f U r u mc h i J D e s i g n o f W a t e r Re s o u r c e s 8 L Hy d r o e l e c t r i c En gine e r i n g， 2 0 0 8， 2 7 ( 3) 4 6 4 8 ( i n Ch i n e s e ) ) 1- 5 3 S L 3 1 9 - 2 0 0 5 ， 混凝土重力坝设计规范 S ( S L 3 1 9 2 0 0 5 ， D e s i g n s p e c i f i c a t io n f o r c o n c r e t e g r a v i t y d a m s s ( i n C h i n e s e ) ) 6 3 曹海， 夏世法， 林锋， 等 浅谈某严寒地区碾压混凝土坝施工质量 控制 J 南水北调与水利科技， 2 0 1 0 ， 8 ( 3 ) ， 1 4 9 1 5 2 ( C A O Ha i ， XI A S h i f a ，LI N F e n g，e t a 1 Th e Pr e l i mi n a r y S t u d y o f Qu a l i t y Co n t r o l f o r t h e Co n s t r u c t io n o f R C C Da ms i n Co l d a n d D r y R e g i o n D S o u t h - t o - N o r t h Wa t e r T r a n s f e r s a n d Wa t e r SCi e n c e& Te c hn o l o g y， 2 0 1 0， 8 ( 3 )， 1 4 9 1 5 2 ( i n Ch i n e s e ) ) 7 阮新民， 夏世法， 鲁一晖， 等 高寒地区某 R C C重力坝保温措施 ( 上接 第 4 8页) 参考文献： 1 徐宗学， 孟翠玲， 巩同梁， 等 西藏 自 治区气温变化趋势分析口 自然资 源学 报 ， 2 0 0 9 ， 2 4( 1 ) ， 1 6 2 1 7 0 ( XU Z o n g - x u e ， ME NG Cu i l i n g， GONG To n g - l i a n g， e t a 1 Lo n g - t e r m Tr e nd An a l y s i s f o r A i r T e m p e r a t u r e i n t h e T i b e t a n R e g i o n J J o u r n a l o f N a t u r a l Re s o ur c e s ， 2 0 09， 2 4 ( 1 )， 1 6 2 - 1 7 0 ( i n Chi n e s e ) ) 2 Z h e n g ， H ， Z h a n g ， L ， L i u ， C ， S h a o ， Q ， F u k u s h i m a ， ， C h a n g e s i n s t r e a m fl ow r e g i me i n h e a dwa t e r c a t c h me n t s o f t h e Ye l l o w Ri v e r b a s i n s i n c e t h e 1 9 5 0 s J Hr d r o l o g i c a l P r o c e s s e s ， 2 0 0 7 2 1， 8 8 6 - 89 3 3 赵少华 ， 杨永辉， 邱国玉， 等 河北平原 3 4 年来气候变化趋势分 析 J 资源科学， 2 0 0 7 ， 2 9 ( 4 ) ， 1 0 9 一 I 1 3 ( Z H A O S h a o - h u a ， YAN G Yo n g - h u i ， QI U Gu o - y u ， e t a 1 An a l y z i n g t h e C l i ma t i c Ch a n g e s i n He b e i P l a i n o v e r t h e L a s t 3 4 Ye a r s J Re s o u r c e s S c i e n c e， 2 0 0 7， 2 9 ( 4 )， 1 0 9 1 1 3 ( i n Ch i n e s e ) ) 及效果分析 J 南水北调 与水 利科技 ， 2 0 0 9 ， 7 ( 3 ) ， 2 7 3 0 ( R U AN Xi n - mi n，XI A S h i f a ，LU Yi - h u i ，e t a 1 Ap p l i c a t i o n a n d Ef f e c t Ev a l u a t i o n f o r t h e He a t Pr e s e r v a t i o n M e a s u r e s o f a RCC G r a v i t y D a m i n Col d R e g i o n J S o u t h - t o - N o r t h Wa t e r T r a n s f e r s a n d Wa t e r S c i e n c e& Te c h n o l o g y ， 2 0 0 9 ， 7( 3 ) ， 2 7 3 0 ( i n Chi n e s e ) ) 8 牛万吉， 王建平， 罗清萍， 等 严寒干旱地区 R C C重力坝的保温 防裂措施 J 水 利 水 电科 技 进 展 ， 2 0 0 9 ， 2 9 ( 1 ) ， 4 4 4 6 ( NI U Wa n - j i ， WANG J i a n - p in g ， L UO Qi n g p ing ， e t a 1 Re s e a r c h o n t h e he a t p r e s e r v a t io n a n d c r a c k p r e v e nt i o n me a s u r e s f o r RCC g r a v i t y d a m i n s e v e r l y c o l d a n d d r y a r e a J Ad v a n c e s i n S c i e n c e a n d Te c h n o l o gy o f W a t e r Re s o u r c e s ， 2 0 09 ，2 9( 1)， 4 4 4 6 ( i n Ch i n e s e ) ) 9 劳道邦， 王炳恒， 王乃超， 等 温泉堡水库碾压混凝土拱坝低温表 面保护研究 J 南水北调与水利科技， 2 0 0 7 ， 5 ( 6 ) ， 1 0 2 1 0 6 ( LAO Da o - b a n g，W ANG Bi 11 g h e n g，W ANG Na i e h a o，e t a 1 Lo w Te mp e r a t u r e S u r f a c e Pr o t e c t i o n Re s e a r c h o n t h e We n q u a n p u Re s e r v o i r RC C Ar c h Da m J S o u t h t o - No r t h Wa t e r Tr a n s f e r s a n d Wa t e r S c i e n c e Te c h n o l o g y ， 2 0 0 7 ， 5 ( 6 ) ， 1 0 2 1 0 6 ( i n Ch i n e s e ) ) 1 0 马涛， 吕琦 ， 孙胜利 ， 等 严寒地区碾压混凝土坝温度场和温度 应力仿真分析 J 水资源与水工程学报， 2 0 0 9 ， 2 0 ( 3 ) ， 1 1 9 1 2 3 ( MA Ta o ， L U Qi ， S UN S h e n g - l i ， e t a 1 S i mu l a t i o n An a l y s i s o f Te mp e r a t u r e Fi e l d a n d Th e r ma l St r e s s o f RCC Da m i n E x t r e me Co l d Re g i o n J J o u r n a l o f Wa t e r Re s o u r c e s& Wa t e r En gine e r i n g， 2 0 0 9， 2 0( 3) ， 1 1 9 1 2 3 ( i n Ch i ne s e ) 4 B o u c h e t R ， E v a p o t r a n s p i r a t i o n r e e l l e a t p o t e n t ie l l e ， s i g n i f i c a t io n c l i m a t i q u e J I AH S P u b l i c a t i o n ， 1 9 6 3 ， 6 2 ， 1 3 4 1 4 2 5 3 B r u t s a e r t W， P a r l a n g e ， L B ， H y d r o l o g ic c y c l e e x p l a i n s t h e e v a p o r a t i o n p a r a d o x J N a t u r e ， 1 9 9 8 ， 3 9 6 ， 2 9 3 0 6 孙福宝， 杨大文， 刘志雨， 等 基于B u d y k 0 假设的黄河流域水热 耦合平衡规律研究 J 水利学报， 2 0 0 7 ， 3 8 ( 4 ) ， 4 0 9 4 1 6 ( S U N Fu b a o， YANG Da - we n， LI U Zh i y u， e t a 1S t u d y o n c o u p l e d wa t e r - e n e r g y ba l a n c e i n Ye l l o w Ri v e r ba s i n b a s e d on Bu d y k o Hy p o t h e s i s J S h u i L i X u b a o ， 2 0 0 7 ， 3 8 ( 4 ) ， 4 0 9 4 1 6 ( i n C h i ne s e ) ) 7 陈福军 ， 沈彦俊 ， 胡乔利 ， 等 海河流域 N DVI 对气候变化的响应 研究 J 遥感学报( 待发) ， 2 0 1 0 ( C H E N F u - j u n ， S HE N Y a n - j a n ， HU Qi a o - l i ， e t a 1 Re s p o n s e s o f NDVI t o c l i ma t e c h a n g e i n H a i h e R i v e r B a s i n J J o u r n a l o f R e m o t e S e n s i n g( a c c e p t e d ) ， 2 O1 O ( i n Chi n e s e ) ) 试 硕 。 究