《工学论文-土石坝渗流观测资料分析模型及方法》由会员分享，可在线阅读，更多相关《工学论文-土石坝渗流观测资料分析模型及方法（9页珍藏版）》请在金锄头文库上搜索。
1、土石坝渗流观测资料分析模型及方法摘要：为了有效地根据土石坝的原型观测资料来坝体和坝基中存在的渗流，本文以渗流为基础，结 合工程实际，考虑了水位、降雨量、坝基排水以及时效等因素对渗流的，建立了土石坝坝体和坝基 测压管以及通过大坝渗流量等观测资料的统计分析模型和。经对青山和对河口水库等工程的实际， 说明其效果是令人满意的。 关键词：土石坝 渗流分析 渗流观测 分析模型 在土石坝坝体和坝基适当部位，有计划地设置一些测压管或渗压计，以及在其下游适当部位设置观 测渗流量的量水堰，并进行观测，可及时了解水库在运行过程中坝体的浸润线位置和渗流区各点渗 透压力的大小，以及通过坝体和坝基渗流量的变化情况，这对大
2、坝的渗流和稳定分析都具有很大的 实际意义。对土石坝各部位的测压管水位和渗流量，选用合理的分析模型进行及时的分析是监测土 石坝运行安全的重要。本文从渗流的支配方程入手，建立了土石坝中有压、无压渗流及其渗流量观 测资料的分析模型。经过实际应用表明，它可较好地解决实际工程问题。1 土石土石坝坝渗流的支配方程渗流的支配方程 忽略地下水流动方程中的惯性项，土石坝渗流的支配方程1,2为 (1)渗流场为均质各向同性时，式（1）变为 或(2) 式中：kx、ky、kz分别为x、y、z方向上的渗透系数，h为水头，=-kh为渗流速度势。 对稳定渗流而言，它的解实际上可归结为在满足某特定边界条件下，求解上述方程式。对
3、无压渗流 问题，由于浸润面事先为未知边界，故在求解过程中，先假定浸润面边界，然后需通过反复试算， 才可以对问题进行求解。根据 .和 .的推导结果1，对具有自由面的缓变渗流，当坐标轴位于不透水层面时，其不稳定渗流的方程形式为：。在稳定渗流时，则渗流方程的 形式为：。以上式中：H为水深函数；ne为有效孔隙率；t为时间。在这种情况下，浸 润线位置即是方程中的一个变量，故它无需作为边界条件来考虑。由于这时地下水流水深函数H的平方项亦满足拉普拉斯方程，故只需以H2为基本变量，就可求解有压渗流一样的方法解决无压渗流 问题。因此人们常将上述方程应用于无压渗流问题中。2 坝坝基有基有压压渗流渗流观测资观测资料
4、分析料分析 根据上述渗流支配方程的基本特性，当渗流场固定时，各点的位势应不随时间而变。位势可用下式表示：。式中：hi为测压管水位，H1、H2分别为上下游水位。对坝体基本不透水 的坝基有压渗流而言，由于渗流场是不变的，故渗流场中某固定点的位势与上下游水位变化无关 ，基本上为一常数。若渗流场中某点的位势随时间增大，则意味着上游铺盖遭受局部破坏或在坝基 内部产生渗透变形；相反，若位势减小，则意味着水库淤积或坝基内部淤填等等。因此对坝基有压 渗流问题，通常可采用位势法来分析其坝基渗流情况。 对于渗流场中有中间位势的有压渗流，则场中各点位势不仅是空间位置的函数，还与中间位势有关 ，中间位势固定，渗流场中
5、各定点的位势亦固定。所以，在对坝基中设有减压井情况下的渗流分析中，应结合减压井的位势进行分析，其合理分析模型为：i=a+bj。式中：i为坝基测压管位势；j为坝基排水减压井位势；a、b为回归系数。 青山水库是浙江省的一座大(2)型水利工程，坝基属双层地基，并采取了铺盖、减压井和反滤盖重 等渗流控制措施。其坝基的渗流是属于有中间位势的有压渗流。作者采用了测压管位势与减压井位 势的相关分析，然后对各个特定减压井位势下各测压管位势进行分析，进一步搞清了坝基所存在的渗透变形问题。现以典型的0+400观测断面（观测设施的断面布置见图1）为例3进行分析。该断 面历年测压管位势与减压井位势的相关系数见表1。减
6、压井特定位势下的测压管历年位势过程线见 图2。 表1 历年测压管位势与减压井位势的相关系数管号19891990199119921993199419951996199719981999B4d4E4F40.8360.9830.9340.9900.9080.9700.9620.9890.8880.9750.9870.9960.9500.9740.8960.9890.8490.9350.9660.9650.4410.9070.9820.9770.0320.9780.9800.9770.5610.9550.9800.9910.3660.9650.9810.9910.5190.9430.9560.9720
7、.0370.9480.8980.993图1 青山水库0+400断面测压管分布示意 （单位：m）由表1和图2所示，在不同的井位势下，B4管在1993年以前其位势与井位势的相关性较好，管位势变 化基本正常。但1993年以后，管位势与井位势关系曲线的相关系数由原来的0.8360.950(19891 993年)降至0.0320.561，不同井位势下的测压管位势变化曲线出现明显异常，说明该管位势基本 上已不受减压井位势所控制，坝基渗流场出现明显变化。从历年库水位变化过程线来看，B4管位势 变化出现上述异常现象的原因是由于大坝经过1993年、1995年、1996年、1999年库水位屡创新高的 影响，导致坝
8、基土体产生一些局部冲刷和淤积现象所形成的。D4、E4管的位势变化则出现历年下降 的趋势，很显然这两管位势下降的原因是由于测压管与减压井之间的坝基土体遭受局部冲刷所致。由此可见，坝基内部已出现较为明显的渗透变形。由坝基渗透坡降的表明3，在不同水位条件下 ，B4、D4管之间的渗流比降为0.1630.277，它们不仅远大于允许比降值0.12，并且还超过了 该地基土的破坏比降0.24。这与上述分析结果相吻合。 图2 各特定减压井位势下的测压管位势历年变化曲线3 坝坝体无体无压压渗流渗流观测资观测资料料 对于坝体无压渗流而言，测压管水位与其所处的位置、上下游水位变化幅度、降雨量强度和筑坝材 料的渗透性等
9、因素有关。由于库水位变化传递到测压管位置存在一定的滞后时间，所以测压管水位 变化还与前期库水位有关。此外，象土体固结，坝前淤积，下游冲刷等因素对渗流状态亦会产生。 归纳起来，影响测压管水位的因素主要包括库水位、降雨量和时效3个方面，根据具有自由面的无 压渗流支配方程的基本特性，库水位与测压管水位并不成线性关系，但它们的平方项则基本满足线 性关系。故库水位与测压管水位(以相对不透水层面为基面)之间合理的回归方程型式可表示为： (3)式中：h为测压管水位；H为当日库水位；Hi-j为前期库水位的均值；R0为当日降雨量；Ri-j为前期平均降雨量；为时效。 针对青山水库0+400断面的4支坝体测压管，参
10、照上述模型，其回归方程3为 (4) 式中：H0为观测日当天的库水位；Hi-j为观测日前第i天到第j天的平均库水位；R0为观测日当天的降雨量；Ri-j为观测日前第i天到第j天的平均降雨量；为观测日至正常观测起始日天数的百分之一；ak(k=0 21)为待定系数； 根据选定的统计模型，取各测压管所有的观测数据进行逐步回归，并求出其最佳的回归方程，简要 回归成果见表2，测压管水位实测值和回归值的过程线见图3。从图表中看出回归曲线与实测过程 线拟合得较好。 图3 测压管实测与模型拟合曲线表2 各坝体测压管回归分析简要成果测压管号样本数复相关系数标准差F检验值F0.01临界值A4C4G41G4235385
11、2843490.9500.8810.8940.9450.3580.3610.2390.166187451223442.133.242.252.85 4 渗流量渗流量观测资观测资料的分析料的分析 渗流量观测不仅能了解水库的渗漏损失，更重要的是它能有效地监视土石坝的安全情况，如果渗流 量随时间发生了较大的变化，则意味着渗流场内可能己产生局部的渗透破坏现象，就应及时找出其 原因并采取必要的补救措施。大坝总渗流量通常由3部分组成：(1)通过坝体的渗流量；(2)通过 坝基的渗流量；(3)通过两岸绕渗或两岸地下水补给的渗流量。为了监测各部分的渗透稳定性，应 尽量做到分区观测。 渗流量的大小除直接受上下游水
12、位变化的影响，还与大坝断面轮廓、坝体和坝基土渗透系数大小以 及其防渗排水设施等因素有关。此外，由于土体固结、坝前淤积，坝体或坝基产生渗透变形而导致 渗流出口土体被冲刷。由于温度变化而会改变水的粘滞性和影响基岩裂隙的张开度等，以上这些因 素均会对渗流量产生影响。综上所述，影响渗流量的因素主要包括；上下游水位、降雨量、时效及 温度等。但在土石坝中，温度通常影响较小，可不加以考虑。由前述渗流支配方程可见，对于有压与无压渗流共存情况，其渗流量应为大坝作用水头H和H2的函数。由库水位和渗流量过程线可见， 渗流量随库水位的变化同样亦存在滞后时间，故它亦与前期库水位有关。其合理的回归方程型式为 ： (5)
13、式中：Q为渗流量；H0为当日库水位；Hi-j为前期库水位的均值；R0为当日降雨量；Ri-j为前期平均降雨量；为时效。对河口水库是浙江省的大(2)型水库，主坝和副坝坝区均设置了量 水堰。主坝1号堰测量主坝坝体渗流量，2号堰测量绕坝渗流量。副坝1号堰测量副坝坝体渗流量，2号堰测量绕坝渗流量。针对对河口水库的具体情况，可参照上述回归方程，其具体形式4取为 (6)表3 各渗流量回归分析简要成果量水堰号样本数复相关系数标准差F检验值F0.01临界值主坝1主坝2副坝1副坝234893050360731610.8750.8200.8560.8580.2600.2250.3750.10070748065654
14、92.000.142.042.00 式中：Q为渗流量；H0为观测日当天库水位；Hi-j为观测日前第i天到第j天的平均库水位；R0为观测日当天的降雨量；Ri-j为观测日前第i天到第j天的平均降雨量；为观测日至正常观测起始日天数的百分之一；ak(k=0 18)为待定系数。 根据选定的统计模型，选入和剔除因子的F检验值取2.50，通过逐步回归的求出最佳回归方程，回 归简要成果见表3。从表中可以看出，4个量水堰所测得渗流量的复相关系数都在0.82以上，剩余 标准差为0.10.375，说明回归方程具有一定的精度。 根据回归方程选入因子，对各分量进行偏相关分析，求得渗流量与各分量的偏相关系数，结果见表 4
15、。从表4可见，影响渗流量的主要因素是降雨量，其次是库水位与时效分量、说明量水堰测出的 渗流量受降雨影响较大。故在坝址设置量水堰时应尽量减少或消除地面径流和降雨对渗流量的影响 ，通常可在坝面及两岸设置排水沟，将非渗流水另行排除，此外设置量水堰的位置亦应考虑到尽量 不受建筑物泄水的影响。 表4 偏相关分析成果量水堰号复相关系数偏相关系数分量的影响程度库水位降雨量时效分量第一第二第三主坝1主坝2副坝1副坝20.8750.8200.8560.8580.5060.4380.5540.5950.6420.5990.5630.5630.5410.4770.2280.419降雨量降雨量降雨量库水位时效分量时效
16、分量库水位降雨量库水位库水位时效分量时效分量 5 结结 语语 (1)在具有排水减压井的有压渗流中，由于渗流场中存在中间位势。故在位势分析时，应采用在不 同中间位势（井位势）条件下，来分析各测压管位势的变化。(2)对具有自由表面的无压渗流而言 ，由于库水位与测压管水位并不成线性关系，但它们的平方项则基本满足线性关系，故可用h的平方项为变量来建立回归方程。(3)对有压与无压并存的渗流情况，渗流量为大坝作用水头H和H2的 函数，并且和前期库水位有关，故可根据上述影响因素来建立其回归方程。(4)根据本文所提供的 分析模型和方法，分别对青山水库和对河口水库等工程的原型观测资料进行了分析，分析结果表明 ，其效果良好。参参 考考 文文 献献： 1 , . M. ,1953.
