考虑非线性环境量影响的混凝土重力坝坝基扬压力分析.pdf
《考虑非线性环境量影响的混凝土重力坝坝基扬压力分析.pdf》由会员分享,可在线阅读,更多相关《考虑非线性环境量影响的混凝土重力坝坝基扬压力分析.pdf(6页珍藏版)》请在咨信网上搜索。
1、第3 3卷 第8期 2 0 1 6年 8月 长江科学 院 院报 J o u r n a l o f Ya n g t z e R i v e r S c i e n t i fic R e s e a r c h I n s ti t u t e V0 1 33 No 8 Au g 2 0 1 6 d o i : 1 0 1 1 9 8 8 c k y y b 2 0 1 5 0 4 7 8 2 0 1 6 , 3 3 ( 0 8 ) : 4 2 - 4 6 , 5 8 考虑非线性环境量影响的混凝土重力坝 坝基扬压力分析 徐蔚 , 方春晖 , 王雪 , 朱凯 ( 1 浙江水利水电学院 建筑工程
2、学院, 杭州3 1 0 0 1 8 ; 2 浙江省水利河I = I 研究院, 杭州3 1 0 0 2 0 ; 3 浙江同济科技 职业学院 水利工程系, 杭州3 1 1 2 3 1 ; 4 河海大学 水文水资源与水利工程科学国家重点实验室, 南京2 1 0 0 9 8 ) 摘要: 基于混凝土重力坝上游库水位及降雨对其坝基扬压力影响的滞后效应及其非线性特征 , 建立 了考虑环境 量非线性影响过程的混凝土重力坝坝基扬压力监测模型, 将上游库水位及降雨对坝基扬压力的影响视为一个正态 分布过程 , 并采用量子遗传算法, 对相应滞后天数和影响天数进行智能寻优, 以提高拟合及预测精度。最后 , 以某 混凝土
3、重力坝为例, 采用上述方法对其 A 6 一 U P O 1 测点扬压力测值进行拟合和预测 , 通过与坝基扬压力实测值对 比分析, 验证了本文方法优于传统坝基扬压力统计模型, 具有一定科学和实践应用价值。 关键词: 混凝土重力坝; 扬压力 ; 滞后效应; 传统统计模型; 分组量子遗传算法 中图分类号: T V 6 9 8 1 文献标志码: A 文章编号 : 1 0 0 1 5 4 8 5 ( 2 0 1 6 ) 0 8 0 0 4 2 0 5 1 研究背 景 混凝土重力坝渗流安全性评价是大坝安全监控 的重要组成部分 , 其 中坝基扬压力作为作用在混凝 土重力坝上 的主要荷载 , 其 大小直接影
4、响坝体 的稳 定性。研究表 明, 坝高 1 0 0 m左右的重力坝 , 坝基 面上作用 的扬 压力大约是坝体 重量 的 2 0 _ 2 J 。而 大坝监测资料分析表 明, 坝身混凝土及基岩在上游 库水压力的作用下 , 其 防渗性能具有逐渐降低的趋 势 j 。因此 , 及时整理分析坝体、 坝基扬压力 的观 测资料 , 建立坝基扬压力统计模型并分析相应的变 化规律 , 对于验算大坝的稳定性和耐久性 , 掌握大坝 的安全性态 , 了解坝身混凝 土的抗渗性能 以及坝基 的帷幕 、 排水效应和坝基情况 的变化等 , 都具有十分 重要的意义 。 传统混凝土坝坝基扬压力统计模型中将前期库 水位或降雨对测压管
5、水位 的影 响视为一个 “ 平均 ” 的过程 , 即取观测 日前 i 天的库水位平均值对测 压管水位进行拟合 , 而实际上 , 上游水位对混凝土重 力坝坝基扬压力 的影响是一个先增 长后下 降的过 程 , 而不是一个平均 的过程 , 为此需要深人研究其因 子 的选择。 因此 , 本文结合混凝土重力坝渗流机理 , 综合考 虑上游库水位及降雨对坝基扬压力的非线性影响 , 将前期库水位及降雨对测压管水位的影响视为一个 正态分布过程 , 建立混凝土重力坝坝基 扬压力统计 模型 , 并以某混凝土重力坝为例 , 采用改进遗传算法 对上游库水位分量及降雨分量 的滞后天数与影响天 数进行智能寻优 , 确立最终
6、混凝 土重力 坝坝基扬压 力统计模型。最后 , 通过将本文方法拟合结果与混 凝土重力坝坝基扬压力传统统计模型拟合结果进行 对 比分析 , 验证了本文方法有效性 。 2 基本 理论 2 1 混凝土重力坝坝基扬压力统计模型 根据实测资料分析表 明, 混凝土重力坝坝基扬 压力主要受上游水位 、 降雨 、 基岩温度及 坝前 淤积 、 坝基帷幕 防渗 和排水效应 的时变 特性等影 响 。 。 为此 , 本文混凝土重力坝坝基扬压力统计模型采用 的形式为 P=P H+P R+P 。 +Pr 。 ( 1 ) 式 中l P为坝体测压管水位 ; P 是水位分量 ; P 是 降雨分量 ; P 是时效分量 ; P 为
7、温度分量。 2 1 1水位 分量 P 通过综合考虑混凝土重力坝坝体坝基扬压力与 上游库水位非线性关 系 , 不 同于传统统计模型将 上 游水位对坝基扬压力 的影 响视 为一个“ 平 均” 的过 程 , 如图 1 。根据上游水位对 坝基扬 压力 的影 响过 程 , 如 图 2 , 其 基本 服从 正 态分 布 过程 , 称 为影 响 曲线 。 收稿 日期 : 2 o 1 5 - o 6 0 6 ; 修 回日期: 2 o 1 5 - o 6 1 8 作者简介 : 徐蔚 ( 1 9 7 2 一 ) , 女, 浙江杭州人 , 副教授 , 硕士 , 主要 从事土木工程 防灾减灾研究 , ( 电话 ) 0
8、 5 7 1 8 6 9 2 9 0 5 5( 电子信箱 ) s h r g q 1 6 3 COI n 。 第8期 徐蔚 等考虑非线性环境量影响的混凝土重力坝坝基扬压力分析 4 3 , I X 3 尸 。 地 图 1 前 x a 天平均水位 图 2 正态分布影 响曲线 ( P一 影响权 ) ( P 一 影响权 ) F i g 1 Av e r a g e 、 v a t e r l e v e l Fi g 2 No r m a l d i s t r i b u t i o n o f fi r s t 3 d a y s ( P- I F) c u r v e ( P - I F) 基于
9、此 , 相应混凝土重力坝坝基扬压力上游水 位分量采用形式为 = 【 去e -(t-x 1)2 dt】 = m1 。 ( 2 ) 式中: 口 为水位分量的回归系数 ; 为水位分量的 滞后天数 ; 为水位分量的影响天数 ; ( t )为上游 水位 ; 为有效水位。 其 中, 。 , 需要通过试算求得 , 在求得实际滞 后天数和影响天数后 , 即可提高坝基扬压力统计模 型回归精度 , 为评价 防渗效果提供依据。此外 , 上游 水位一般每天只要一个测值 。因此 , 本文将连续性 积分改成离散型积分, 积分区间只需取 的2 3 倍 即可满足要求。 2 1 2 降雨分量 PR 在降雨过程 中, 有一部分人
10、渗产生地下水 , 地下 水主要通过节理裂隙渗流影响两岸坝段坝基的扬压 力 l 8 。在该过程 中呈 明显的非线性关 系, 并且具有 滞后效应, 采用前期降雨量的均值作为因子, 同样存 在难 以精确模拟降雨对扬压力 的影响。 单元裂隙岩体中渗流基本服从指数定律 , 即 =一 后 。 ( 3 ) 、 式中: 为渗流速度; k为渗透系数; 为孔隙水压 力 ; 为水 的重度 ;m为指 数形式渗 流定律 的指 数 , m O 5 , 3 2 5 。 根据降雨对地下水位的影响规律 和裂隙渗流的 指数定律 , 本文坝基扬压力降雨分量采用形式为 : e 】 = b 尺 。 r 】 。 ( 4 ) 式中: b
11、为降雨分量的回归系数; , 为降雨分量的滞 后天数; 为降雨影响权正态分 布标准差 ( 影响天 数 ) ; R( t )为 t 时刻的单位时段降雨量 ; 为有效降 雨 量 。 同样 , 和 经过试算求得 , 连续性积分可转化 为离散型积分 。 2 1 3温度分 量 P 渗流受地基裂 隙变化的影 响, 孔 隙变化受基岩 温度的作用 。而基 岩温度变化较小 , 基本呈年 周期变化规律 , 本文采用正弦波周 期函数作为温度 分量 , 即 P = c j i s i n z s ) 。( 5 ) 式 中 C C 为温度分量的回归系数。 2 1 4时效 分 量 时效分量是混凝土重力坝坝基扬压力的重要组
12、成部分 , 也是评价大坝渗流稳定性的重要依据 , 本文 将坝体测压管水位时效分量表示为 P 口=d 1 0+d 2 l n O 。 ( 6 ) 式 中: d , d 分别为时效分量线性项 、 对数项的回归 系数 ; 0为从蓄水初期或工程措施初期开始的天数 除以 l O 0, 即每 1 0 0 d增加1 0 。 综上 , 本文所取混凝土重力坝坝基扬压力统计 模型为 P + e f + C l i s i n + c s )+d l 0+d 2 l n O。 ( 7 ) c 咖 )+ + 。 ( 7 ) 然而 , 采用常规 方法计算库水位和降雨对坝基 扬压力影响的滞后天数及影响天数效率低 、 精度
13、差 , 为此本文采用改进遗传算法对上述参数进行智能寻 优 , 提高计算效率和精度 。 2 2 分组量子遗传算法 量子遗传算法( q u a n t u m g e n e t ic a l g o r i t h m, Q G A ) 结合量子计算与遗传算法 , 是一种新发展起来的概率 进化算法 引。该方法 以量子理论为基础, 采用量 子位概率编码表示染色体 , 通过不断更新量子旋转门 的作用来更新和优化种群, 达到搜索的目的。其具有 种群规模小 , 收敛迅速和全局寻优能力强等特点 , 并 在求解组合优化问题中取得显著成效 引 。 1 个量子位可以处于 0态 、 1态、 以及 1 0 ) 和
14、l 1 ) 之间的任意叠加态 。可以表示为 I )= I O )+ l 1 ) 。 ( 8 ) 式中: , 分别为量子位对应态的概率幅。I I 为 量子态达到 1 0 ) 态时相应的概率, l I 为量子态 达到 I 1 )态 的概率 , 并且满足归一化条件 , 即 l l + =1 。 ( 9 ) 因此 , l ) = c 。 s 导l 0 ) + e U s i n 导l 1 ) 。( 1 0 ) 长江科 学院 院报 2 0 1 6焦 式 中 : c 。 s 导 和 e sin 导 为 复 数 ; I c 。 s 0 I 和 I e sin 0 1 分 别 为 量 子 位 I 。 ) 和
15、I 1 ) 的 概 率 ,并 且 满足 J c 。s 詈 + e sin 导 I = 。 c , 传统量子遗传算法 ( Q G A) 通过量子位 进行编 码 ,通 过 概 率 幅 将 1 个 量 子 位 表 示 为 【 个量子位可以定义为 保 留本代最优 , 最后在每代计算 中 比较得 出全局最 优 。为了避免 目标 函数进行线性 、 非线性 等加 速适 应 函数出现的早熟现象 , 本文采用排序适应 函数 , 即 将 同一个种群的 m个染 色体 按 目标 函数值进行 从 小到 大 的 排 列。直 接 概 率 可 以 表 示 为 P ( i )= 2 , m 。 。 因此, n 3工程实例 。
16、其中 I i 2 + l =1 , i =1 , 2 , , n。任何量 子叠加态都可 以表示为这种形式。 在分组量子遗传算法 中, 量子旋转 门可以表示 为 吣 = , ( ) 。伽 I 。 ) 量子 门更新过程可以表示为 = -c 。 其 中,( i ) 和 ( 卢 ) 分别代 表染 色 体第 i 个量子比特旋转门更新前后的概率幅; 0 为 旋转角。( , 卢 ) 可以表示为 f : 。 。 i 一 in ;( 1 5 ) i O s i nO + c o s O 。 所 以, 可以得出: I i I +l卢 I =( c o s O 一 s i n O ) + ( s i n O + c
17、 o s O ) I I +I l =1 。 分组 量子遗传算法 流程如 图 3 。分组 量子 遗 传 算 法 通 过 运 行 代, 每一代相应随机产 生 r t 个 样 本 , 并 通 过每 一 层子种 群按类运行相 应 的量 子遗 传 算 法 , 相 应遗传算 法均采用不 同 的量 子 叠 加 态形 式 , 从 而保 持 种 群 的 多样 性 。 在每一代适应值计算 中 擘 生 掣 量 子 门 同 向 旋 转 l I 量 子门 异 向 旋 转 l l LI 产 生 新 的 叠 加 态 图 3 算法流程 Fi g 3 F l o w c h a r t o f al go r i t hm
18、3 1 工程简介 某水电站位 于云南省丽江市境 内, 电站设计总 装机容量为 2 4 0万 k W, 坝顶全长 6 4 0 I n , 最大坝高 1 6 0 m, 总库容9 1 3 亿 m 。枢纽工程 由碾压混 凝土 重力坝 、 右岸溢洪道 、 左右岸坝身泄洪冲沙孔和坝后 厂房等建筑物组成 , 右岸 台地布置 5孔开 敞式溢洪 道 , 其最大泄洪流量为 1 4 9 8 0 i n s , 坝后厂房安装 4 台单机容量 6 0 0 M w 水能发电机组 。水库正常蓄水 位 1 4 1 8 m, 校核水 位 1 4 2 1 0 7 n l , 相应库 容8 4 7 亿 r f l 。 ,具有周调节
19、能力 , 控制流域面积2 3 7 3 5 7 6 k m 。 电站 运 行 时 保 证 出力 4 7 3 7 MW, 年 平 均 发 电量 1 1 0 4 3 亿 k W h , 发电效益显著, 是“ 西电东送” 战 略 目标 的骨干电站之一 。 为了监测大坝渗流状况 , 根据监测施工单位统 计 , 大坝共 安装渗压计 1 0 6支 , 量水堰 9座 、 坝基 测 压管 4 6个 、 水位孔 1 1 个 、 水力学底座 1 3 1个 , 其 中 坝纵 0 + 1 5 7 5 0 i n 断面大坝渗流监测布置图如图 4 。 图 4 大 坝渗流监测布置 Fi g 4 La y o u t o f
20、mo n i t o r i n g p o i n t s f o r d a m s e e p a g e 第8期 徐 蔚 等 考虑非线性环境量影响的混凝土重力坝坝基扬压力分析 4 5 监测资料分析表明, 该坝段 A 6 - U P - 0 1 测点坝基扬 压力监测数据常年处 于较高水平 , 因此本文将对该 测点进行重点分析 。 3 2混凝 土重 力 坝坝基 扬压 力分 析 本文选取 A 6 - U P - 0 1测点 2 0 0 7年 2月 1 9日至 2 0 0 8年 1 2月 1 5日坝基扬压 力测值采用上述模 型 进行拟合。分组量子遗传算法控制参数设定如下: 取群体规模为3 0
21、0 0 , 量子位定为 1 5 , 初始搜索以等 概率叠加, 量子位概率 趔 幅均为 2 , 当迭代 次数 超过 1 0 0或者 复 蜊 相关系数超过 0 9 9时 即停止 迭代 。最 优个 体适应度值进化过程 如图 5 , 相应拟合结果 如图 6 。 进化代数 图5 最优个体适应度值 演化 曲线 Fi g 5 Ev o l u t i o n c u r v e o f o p t i ma l i n d i v i d u a l fit n e s s v a l u e 图 6 本 文方法拟合结果 F i g 6 Fi t t i n g r e s ul t o f t h e p
22、 r o p o s e d me t h o d 由图 5 、 图 6可知 , 采用分组量子遗传算法 最终 优化复相关系数为 0 9 8 1 , A 6 - U P一 0 1 测点扬压力实 测值与拟合值变化规律相同, 本文方法具有较高的 拟合精度。此外, 可以注意到坝基扬压力实测值在 个别 日期具有突变现象 , 初步分析为相应上游库水 位及降雨突变所致。 为进一步验证上述方法有效性 , 本文 同时采用 传统混凝土坝坝基扬压力统计模型 对 A 6 一 u P 一 0 1 测点测值进行拟合 , 该统计模型为 P 茎 + 6 + ( c sin + c o s )+d l 0+d 2 l n 。
23、。( 1 6 ) 式 中: 口 为上游水位分量的回归系数 ; 为观测值 为观测 日前 天平均库水位 ; 为观测 日前 i 天平均 降雨量。 取 m, =m =m =5 , 采用逐步 回归分析对 A 6 - U P - 0 1 测点坝基扬压力测值进行拟合 , 拟合结果如 图 7 。 日 期 图 7 传统方法拟合结果 Fi g 7 F i t t i n g r e s ul t o f t h e t r a d i t i o n a l m e th o d 由图 7可知 , 传统坝基扬压力统计模型拟合结 果与坝基扬压力实测值偏差较大 , 最终复相关系数 为 0 9 1 , 与本文方法相比,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 考虑 非线性 环境 影响 混凝土 重力坝 坝基 压力 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【z****6】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【z****6】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。