2010年_基于统计模型的白石水库坝顶位移分析(韩卫许小华)

基于统计模型的白石水库坝顶位移分析韩卫，许小华辽宁省白石水库管理局，朝阳122000摘要：利用统计模型原理，对白石水库大坝坝顶水平位移和垂直位移进行解析计算，定量分析了影响坝顶位移的主要因素及影响程度，并根据分析成果，对今后大坝安全监测的工作重点给出了建议。关键词：白石水库；坝顶位移；统计模型；分析。1工程概况白石水库位于辽宁省北票市上园镇附近的大凌河干流上，总库容16.45亿m3，是干流上唯一的大（I）型控制性骨干工程。大坝为混凝土重力坝，部分采用RCD碾压混凝土技术。最大坝高49.3m，坝顶长513m，分为32个坝段。水库1996年9月正式开工，1999年9月下闸蓄水。用于坝顶变形监测的真空激光准直系统于2001年1月调试完成，由激光发射端设备、接收端设备、测点设备、真空管道和抽真空设备等组成，在1#～30#坝段坝顶上游侧中部，各设置1个测点，共计30个。2统计模型的基本原理根据实测资料可定性地分析本段各监测量的变化规律，为进一步分析影响监测量的主要因素及影响程度大小，需根据原因量（如库水位、温度等）与效应量（如变形、渗流等）的之间的关系建立数学模型，对影响监测量变形的各因素进行解析计算，解释其变形规律，预测监测量今后的变化值及趋势。根据历史实测资料分析，坝内任一监测量变化的统计模型表达式可表示为：SBYYYYTH2ˆ0式中：HY为水压分量；TY为温度分量；Y为时效分量；0B为常数项；S为剩余标准差。（1）水压分量上游水压分量HY主要是由水荷载引起的，水荷载包括：水重、水平静水压力等，它们均是水位H的函数。上游水压分量与H呈非线性关系，可写成通式如下：iiliHHaaY10式中：0a、ia分别为回归系数，H为上游水位或上游水深，l为水位幂次因子，本次计算取l=4。（2）温度分量由于混凝土材料变形及基岩裂隙的展开与闭合基本上是呈周期性变化，采用正弦波周期函数作为温度因子，温度分量如下：)(10sfbbYkinjT式中：0b、jb分别为温度因子回归系数；)(sfk为温度函数及周期函数，计算中温度因子S为观测日期前不同时段的气温平均值，本次气温S的选取分别为当天值、前15天、前30天、前60天、前90天值。（3）时效分量时效分量是监测量形变的一个重要因素，是分析大坝变形关心的首要问题之一。时效变形的一般规律是蓄水初期或某一工程措施初期变化较快，然后随时间而逐渐趋于平稳，但也会出现随时间变化而增大，应根据坝体实际情况认真分析。时效分量表达式为：)(10tfccYkkmk式中：C0、Ck分别为时效分量的回归系数，t为从起测日开始的天数。引起时效形变的原因较复杂，它的产生可能与大坝混凝土徐变、塑性变形、裂缝变形、基岩徐变，以及外界条件的影响等诸多因素有关，本次计算时效函数)(tfk选用：线性函数)(tL、对数函数)(tLn。3坝顶监测资料分析本文资料分析时段为2001年1月至2007年7月。3.1库水位白石水库为多年调节水库，库水位变化较为平缓，多年在107.0～118.52m左右之间变化，最高水位118.52m出现在2006年4月10日，远未达到正常高水位127.0m；多年库水位变幅为11.52m，各年度库水位变幅较小，在3.18m～6.72m之间变化，多年平均库水位为114.72m。3.2日常平均气温白石水库地处寒冷地区，坝址多年平均气温6.8℃，每年7、8月温度较高，月平均最高气温分别为22.3℃，21.6℃；每年12月、次年1、2月温度最低，月平均最低气温分别为-16.7℃、-20.9℃。3.3坝顶水平位移监测成果本文选取6#、16#、28#坝段为典型坝段，对坝顶水平位移观测成果进行统计计算，在显著水平α=0.05下，统计模型见表1，解析成果见表2。表1坝顶激光准直测点水平位移统计模型成果表测点数学模型相关系数R标准差S(mm)样本数量JD06-XY=-4.292171+0.06612t+4.37333E-05H3-0.05053T1-0.08688T15-0.1133764T60.96880.8283222JD16-XY=-0.89904+0.16585t-1.2608Ln(t+1)-0.08259T1-0.1885414T150.89001.849216JD28-XY=-2.872177-0.06244t+2.46988Ln(t)+2.696677E-05H3-0.0295195T1+0.033579T15-0.046973T600.94960.4620196式中：1、“H”为相对坝基的水深，坝基高程取为85.0m；2、t为时间。表2坝顶激光准直测点水平位移统计解析成果表测点实测变幅(mm)解析变幅(mm)时效分量(%)水压分量(%)温度分量(%)最大拟合残差(mm)相应时间JD06-X15.7711.717.227.1875.607.11401.01.04JD16-X30.813.4120.940.0079.0618.79101.01.04JD28-X6.615.7867.078.5124.42-3.565001.12.133.4坝顶垂直位移监测成果本文选取4#、16#、28#坝段为典型坝段，对坝顶垂直位移观测成果进行统计计算，在显著水平α=0.05下，统计模型见表1，解析成果见表2。表3坝顶激光准直垂直位移统计模型成果表测点数学模型相关系数R标准差S(mm)样本数量JD04-ZY=-2.895297-0.01789t+0.707524Ln(t)-0.011265T30-0.1295T900.96030.4487220JD16-ZY=-0.89904+0.16585t-1.2608Ln(t)-0.08259T1-0.1885414T150.97770.5181217JD28-ZY=-1.484533+0.0331472t-0.83681Ln(t)-0.137007T15+0.2443T60-0.2853624T900.92490.7940202表4坝顶激光准直垂直位移统计解析成果表测点实测变幅(mm)解析变幅(mm)时效分量(%)水压分量(%)温度分量(%)最大拟合残差(mm)相应时间JD04-Z7.424.8924.32075.68-4.66101.01.04JD10-Z12.357.1229.96070.04-8.16301.01.04JD16-Z10.747.9216.30083.70-2.29101.11.06JD22-Z12.5912.238.55091.453.05701.12.19JD28-Z12.898.1618.33081.67-2.490401.12.134大坝运行性状综合评述4.1大坝坝顶水平位移选取6#坝段，绘制坝顶水平位移成果统计解析曲线如图1所示。JD06-X水平位移模型分量图-8-404812位移值(mm)实测值模型值-8-404812位移值(mm)水压分量-8-404812位移值(mm)温度分量-12-8-404801-101-702-102-703-103-704-104-705-105-706-106-707-107-7日期位移值(mm)时效分量（1）由解析计算成果知，坝顶水平位移主要受温度的影响，温度分量占总量的24.4～88.04%；时效影响次之，占总量的11.96～67.07%；受水压的影响较小，占总量的0～13.19%。坝顶各测点位移过程曲线平滑，连续性好。坝顶水平位移受温度影响明显，坝顶水平位移随季节呈周期性变化，温度升高坝体倾向上游位移；温度降低坝体倾向下游位移，坝体呈弹性变化，约有1月左右的滞后，符合一般重力坝变形规律。（2）坝顶各测点上、下游向的相对位移值在-8.6～+24.2mm之间变化，最大值+24.2mm发生在16#坝段，最小值-8.6mm发生在11#坝段；各测点平均位移值在-2.4～+5.6mm之间变化。（3）坝顶各测点均有不同程度的向下游方向的不可逆位移，向下游不可逆位移为1.2～5.5mm，从1#坝段至28#坝段不可拟变形呈逐渐增大趋势，最大值5.5mm出现在28#坝段。4.2大坝坝顶垂直位移选取10#坝段，绘制坝顶水平位移成果统计解析曲线如图2所示。JD10-Z垂直位移模型分量图-6-4-2024位移值(mm)实测值模型值-4-20246位移值(mm)水压分量-4-20246位移值(mm)温度分量-9-7-5-3-1101-101-702-102-703-103-704-104-705-105-706-106-707-107-7日期位移值(mm)时效分量（1）坝顶垂直位移主要受温度的影响，温度分量占总量的70.0～91.5%；受时效影响小，占总量的8.5～30.0%；不受水位影响。温升时坝顶上抬，温降时坝顶下沉，坝体呈弹性变化，约有2月左右的滞后，呈周期性变化，符合重力坝坝体垂直位移变化的一般规律。（2）坝顶相对垂直位移值在-21.64～+7.11mm之间变化，最大值+7.11mm出现在29#坝段，最小值-21.64mm出现在15#坝段。多年平均相对垂直位移值在-0.27～+5.9mm之间变化，变化趋势较为平缓。（3）坝顶各测点均有不同程度呈下沉的不可逆位移，时效分量曲线平缓，其下沉不可逆位移为0.2-3.5mm。4.3建议白石水库大坝坝顶位移过程曲线平滑，呈周期性变化，连续性好，符合一般重力坝变形规律。在目前低水位运行期间，坝顶已产生向下游不可逆位移，今后应重点监测其位移的发展趋势，以及坝踵应力是否有受拉的情况，并加强巡视检查工作。作者简介：韩卫（1971—），男，河北保定人，高级工程师，工程硕士，主要从事大坝安全监测运行管理及资料分析工作，E-mail：[e-mail]。