Determination of Penetration Parameters of a Bank or Dam Foundation with Sensitivity Analyzing Method
-
摘要: 堤坝工程中地基渗透系数的选取具较大随机性, 为了解决这个问题, 采用敏感性分析理论建立了求解渗透系数敏感性因子的数学模型, 该模型在河南省燕山水库水利工程中的应用表明: 覆盖层渗透系数的敏感性最大(> 10-2), 设计应取试验值的大值; 其他部分渗透系数的敏感性相对较小(10-6~10-3), 设计中可取试验值的均值; 对于防渗设施, 只要达到一定密实度, 再降低其渗透系数, 防渗效果不会再有明显改善Abstract: There appears to be random to design the penetration parameters for a bank or dam foundation. To solve this problem, this paper introduces a sensitivity analyzing method to calculate the penetration parameter's sensitivity factor of each part in the calculation model. The application of the method in Yanshan reservoir indicates that the covering layer's penetration parameter's sensitivity factor (> 10-2) is the biggest, whose design parameter should then be the biggest of the testing values. The sensitivity factors (10-6-10-3) of other parts are relatively small, whose design parameters can then be the mean values. As for the anti-penetration objects, if dense enough, the decrease of the penetration parameters will be insignificant.
-
堤防工程或水库坝基渗流模拟中, 要取得符合实际的模拟结果, 计算模型的建立与计算参数的选取同等重要.计算参数又可划分为几何尺寸参数和水力学参数两大类, 其中, 尺寸参数大多为人工控制参数, 易于量测; 而水力学参数, 如渗透系数等, 则因介质的不均一性和测量手段的不合理等原因, 往往具有离散性大、不易准确确定的特点.因此如何确定合理的渗透系数一直是一个困扰工程界的难题.目前比较成熟的做法是通过反演的手段确定合理的渗透系数组合.反演是以一定的实际水位观测资料为基础的, 通过调整计算模型各部分的渗透系数, 使模拟结果最终与实际观测结果吻合, 此时模型各部分的渗透系数即可用作实际模拟时的计算参数.然而, 对于那些尚处于规划设计阶段的水库工程, 显然不具备反演的条件, 此时, 即可采用敏感性分析的方法.
1. 敏感性分析的一般思路
敏感性分析(谭晓慧, 2001) 是系统分析方法的一种, 它以一定的数学模型为基础.对于一固定系统, 可设其影响变量的集合为(x1, x2, …, xn), 系统特征用变量Z表示; 影响变量的基准值集合为(x1′, x2′, …, xn′), 与之对应的系统特征变量取值为z′.分析变量xk对系统特征的影响时, 可令其余变量取基准值且固定不变, xk在其可能的范围内变动, 进而可得变量xk对系统特征的影响曲线
(1) 由此可得变量的敏感度函数
(2) 取xk=xk′, 即得变量xk的敏感因子Sk′
(3) 上述理论应用于水库坝基渗流分析中, 可令坝基最大水力比降J为系统特征值, 其影响变量则为计算模型各部分的渗透系数(k1, k2, …, kn).对于第i部分的渗透系数ki, 套用式(3) 可得其敏感性因子Si′
(4) 其中: ki′为第i部分渗透系数基准值(可取试验值的均值或大值均值); J′为各部分渗透系数取基准值时的坝基最大水力比降; φi′ (ki′) 为坝基最大水力比降表达式对ki取偏导后, 再令ki=ki′时的计算结果; 对于难以作微分处理的复杂系统而言, 可按以下差分格式作近似处理
(5) 由式(4)、式(5) 结合有限元分析方法, 即可确定计算模型各部分渗透系数对整个坝基渗流场的敏感性因子, 进而找出最不利的参数组合.
2. 计算模型及相关试验参数
燕山水库位于淮河流域支流沙颍河主要支流澧河上游甘江河上, 其不仅是历次淮河流域规划所选定的以防洪为主的大型水库, 而且还是南水北调中线引汉工程所选定的反调节水库之一.该坝基存在大规模顺河断层带, 宽100余m, 具渗透变形问题, 是人们关心的重要隐患.
于河槽位置选一垂直于坝轴线的剖面为计算剖面, 上游水位取一期工程正常蓄水位108.14 m, 下游水位取河床地表高程85.48 m.结合《燕山水库项目建议书》中的一期工程坝体结构和现场所作的勘查工作, 得计算模型, 如图 1所示.
模型中坝基各层的渗透系数由试验确定, 坝体各部分的渗透系数为设计控制值, 具体取值如表 1; 依此可确定敏感分析时的取值范围和基准值.
表 1 计算模型各部分渗透系数试验值Table Supplementary Table Calculation model's testing penetration parameters3. 渗透系数敏感性分析
编制有限元计算程序(毛昶熙, 1990; 陈崇希和唐仲华, 1993), 按前述敏感性分析思路对模型各部分渗透系数逐一作敏感性分析, 结果如下:
图 2对比列出了粘土斜墙取3种不同渗透系数时的渗流场模拟结果: 点划线对应渗透系数最小值6×10-7 cm/s, 虚线对应大值均值1×10-6 cm/s, 实线对应最大值3×10-6 cm/s.模拟中其他部分的渗透系数均取大值均值.不难看出, 当粘土斜墙渗透系数取3个不同值时, 其所对应的水头等值线模拟结果基本重合, 说明粘土斜墙只要达到了一定的密实度, 其渗透系数的改变对渗流场状态不会再产生明显的影响.
图 3对比列出了排水体取3种不同渗透系数时的渗流场模拟结果(对应关系如图 3); 模拟时其他部分的渗透系数亦取大值均值.可以看出, 随着排水体渗透系数的加大, 水头等值线在坝底的水平渗流段会向上游偏移, 但任意2组等值线间的间隔却变化甚小; 也就是说, 排水体渗透系数在其变化范围内的改变会影响水头值的分布, 但不会造成水力比降的显著提高或降低.
图 4对比列出的是任意料区渗透系数分别取5×10-5cm/s (填料方案1的大值均值)和3.5×10-1 cm/s (填料方案2) 时的水头等值线.虽然2个渗透系数间隔几个数量级, 但从图 4中可以看出, 渗流场的变化并不十分明显.而且与图 4的规律类似, 任意料区渗透系数的改变只会影响水头的分布, 而不会造成水力比降的提高或降低.
图 5为Q42卵石混合土渗透系数敏感性分析结果.从图 5中不难看出, 虽然其渗透系数只在同一数量级中变化(2.0×10-1~6.0×10-1 cm/s), 但整个坝基的渗流场却发生了显著的改变: 随着渗透系数的加大, 水头等值线间隔明显减小, 即对应位置的水力比降会有明显的提高.
图 6为Q3卵石混合土渗透系数敏感性分析结果.该层渗透系数的改变对整个坝基渗流场(尤其是坝头与坝脚部位) 亦会有较大的影响, 渗透系数的加大会导致坝基平均水力比降的提高.
图 7为断层破碎带渗透系数的敏感性分析结果.从图 7中可以看出3组等值线基本重合, 这说明破碎带渗透系数在其变化范围内的取值对整个坝基渗流场无显著影响.
为更好地指导设计, 还进行了垂直防渗体渗透系数的敏感性分析.模拟时, 上游水位取校核洪水位116.78 m, 下游水位仍取河床地表高程85.48 m; 防渗体厚4 m, 打入Q3卵石混合土层底部.当防渗体渗透系数取10-4, 10-5和10-8 cm/s时, 其水头等值线模拟结果如图 8所示, 当渗透系数小于10-5 cm/s后, 再减小其渗透系数, 防渗效果不再有明显改变, 图 8中10-5 cm/s所对应的等值线(虚线) 与10-8 cm/s所对应的等值线(点划线) 基本重合.表 2中, ki′为渗透系数基准值; J′为各部分渗透系数取基准值时坝基最大水力比降; φi′ (ki′) 为比降表达式对ki偏导后再令ki=ki′时的计算结果; Si′为敏感性因子.
表 2 计算模型各部分渗透系数敏感性分析结果Table Supplementary Table Sensitivity analyzing results of penetration parameters4. 结论
堤防工程及水库坝基渗流模拟时, 如不具备反演条件, 可在试验参数的基础上, 采用敏感性分析的方法确定计算模型各部分渗透系数的最不利组合, 然后以此作为实际渗流模拟时的计算参数.
计算模型不同部位渗透系数的改变对整个坝基渗流场的影响差异很大, 所以在模拟计算时其各自渗透系数的取值标准也应有所不同.对于燕山水库, 有如下分析结论: (1) 粘土斜墙及断层破碎带的渗透系数在其变化范围内发生改变时, 对整个坝基的渗流场特征不会产生显著影响.计算时其渗透系数可取大值均值. (2) 随排水体或任意料区的渗透系数的加大, 在坝底水平渗流段, 水头等值线会发生向上游的整体偏移, 但其偏移幅度很小.这说明排水体和任意料区渗透系数的改变亦不会对坝基渗流场特征产生显著影响, 在具体计算时其渗透系数亦可取大值均值. (3) 坝基覆盖层Q42及Q3卵石混合土渗透系数的改变对整个坝基渗流场会产生较大的影响.其中, Q42渗透系数的加大主要导致了坝基水平渗流段水力比降的大幅度增加, 而Q3渗透系数的加大则会导致坝前和坝后部分水力比降的提高.为了使评价结果偏安全, 在具体计算时两部分的渗透系数均应取最大值. (4) 就粘土斜墙和垂直防渗体而言, 只要它们的渗透系数小于一定值之后, 再降低其渗透系数, 其防渗效果不再会有明显提高.
-
表 1 计算模型各部分渗透系数试验值
Table 1. Calculation model's testing penetration parameters
表 2 计算模型各部分渗透系数敏感性分析结果
Table 2. Sensitivity analyzing results of penetration parameters
-
[1] Chen, C.X., Tang, Z.H., 1993. Numerical analysis of underground water seepage. China University of Geosciences Press, Wuhan(in Chinese). [2] Mao, C.X., 1990. Seepage calculation analyzing and controlling. Water Conservancy and Hydroelectric Power Press, Beijing(in Chinese). [3] Tan, X.H., 2001. Reliability calculation and sensitivity analysis of a plane sliding slope. Geology of Anhui, 11(1): 49-54 (in Chinese with English abstract). [4] 陈崇希, 唐仲华, 1993. 地下水流动问题数值分析. 武汉: 中国地质大学出版社. [5] 毛昶熙, 1990. 渗流计算分析与控制. 北京: 水利水电出版社. [6] 谭晓慧, 2001. 平面滑动边坡的可靠度计算及敏感性分析. 安徽地质, 11(1): 49-54. https://www.cnki.com.cn/Article/CJFDTOTAL-AHDZ200101010.htm -