三峡库区非均质岸坡中浸润线的解析解.docx

三峡库区非均质岸坡中浸润线的解析解 由于地下水是影响水库边缘稳定性的主要因素，在库水位的作用下正确、适当地确定滑动线对边坡稳定性的评价非常重要。然而目前还没有合理的计算浸润线的公式, 设计人员往往依据经验确定浸润线。近年来, 国内外学者用数值模拟法对其进行了研究, 但该方法较为复杂, 使用极为不便。有少数学者研究了浸润线的解析解, 然而存在以下不足:①为简化边界条件而假设的初始潜水面水平, 与实际情况有较大出入;②以隔水底板水平的土坡或均质岩坡为研究对象, 所得公式的使用范围具有局限性, 难以推广。因此笔者拟根据三峡库区广泛存在的隔水底板倾斜的层状非均质岸坡建立基本模型, 并分别讨论在稳定渗流和非稳定渗流情况下, 库水位下降时岸坡中浸润线的确定方法。 1 基本方法 三峡库区广泛存在的缓倾角层状岸坡的简化地质剖面如图1所示。岸坡为粉砂岩、薄层泥岩互层结构。 1.1 在稳定增长过程中求解润湿线的方法 1.1.1 浸渍线方程的建立 吉林斯基势函数法求解浸润线的相关参数如图1所示, 同时假设库水位高于泥岩层。 将断面1的参数代入势函数基本公式得: Φg1=K1M1(H1-M1/2) +K2M2(H1- M1-M2/2) +0.5K1(H1-M1-M2)2(1) 同理可求出Φg2和Φgx。 式中:Φg1、Φg2和Φgx分别为断面1、断面2和其间任一断面处的势;K1、K2分别为粉砂岩、泥岩的渗透系数;M1、M2分别为粉砂岩层、泥岩层的厚度;H1为断面1处的水头值。 单宽流量q可表示为: q=Φg2-Φg1L=Φgx-Φg1x(2)q=Φg2?Φg1L=Φgx?Φg1x(2) 式中:L为断面1和断面2的间距;x为任取断面与断面1的间距。建立XOH坐标系, 解一元二次方程即得: H=1/K1[M2(K1-K2) +K2xtg θ+√(Κ2Μ22+2Κ1Μ2Μ1)(Κ2-Κ1)+2Κ1(Φg1+qx)](3)θ+(K2M22+2K1M2M1)(K2?K1)+2K1(Φg1+qx)???????????????????????????????????????√](3) 式中:θ为隔水底板与水平面的夹角;H为与断面1距离为x断面处的水头值。 当库水位低于泥岩层时浸润线方程的求法同上, 所得结果为一分段函数, 设x1、x2为浸润线与软弱夹层的两个交点的横坐标。 当0≤x≤x1时, 得到位于泥岩层之下的第一段浸润线的方程: Η=1/Κ1[Κ1xtgθ+√2Κ1(Φg1+qx)](4)H=1/K1[K1xtgθ+2K1(Φg1+qx)???????????√](4) 当x1≤x≤x2时, 得到位于泥岩层中的第二段浸润线的方程: Η=1/Κ2[Κ2xtgθ+Μ1(Κ2-Κ1)+√Κ1Μ21(Κ1-Κ2)+2Κ2(Φg1+qx)[JX-*2]][JX*2](5)H=1/K2[K2xtgθ+M1(K2?K1)+K1M21(K1?K2)+2K2(Φg1+qx)??????????????????????????√[JX?*2]][JX*2](5) 当x2≤x≤L时, 得到位于泥岩层之上的第三段浸润线的方程: H=1/K1[M2(K1-K2) +K1xtg θ+√(Κ2Μ22+2Κ1Μ2Μ1)(Κ2-Κ1)+2Κ1(Φg1+qx)[JX-*2]][JX*2](6)θ+(K2M22+2K1M2M1)(K2?K1)+2K1(Φg1+qx)???????????????????????????????????????√[JX?*2]][JX*2](6) 1.1.2 等效水头值的计算 裘布依假定结合等效厚度法的基本假设为库水位高于泥岩层。A、B端面的等效水头值Had、Hbd为: Had=H1-M2+K2M2/K1(7) Hbd=H2-M2+K2M2/K1(8) 单宽流量q可表示为: q=Κ1Ηbd-ΗadL?Ηbd+Ηad-Ltgθ2=Κ1Ηd-Ηadx?Ηd+Ηad-Ltgθ2q=K1Hbd?HadL?Hbd+Had?Ltgθ2=K1Hd?Hadx?Hd+Had?Ltgθ2 式中:Hadad、Hbdbd分别为断面1和断面2处水头值H1、H2的等效水头值;Hd为与断面1距离为x处水头值H的等效水头值。 解一元二次方程并还原等效厚度后得到的浸润线方程为: Η=1/2[a+√a2+4xLΗbd(Ηbd-Ltgθ)+4(Ηad)2(1-xL)]+H=1/2[a+a2+4xLHbd(Hbd?Ltgθ)+4(Had)2(1?xL)????????????????????????????????√]+ Μ2-Κ2Μ2Κ1(9) 式中:a=xtgθ。 1.2 内渗流基本方程 基本假设:①含水层均质、各向同性, 侧向无限延伸, 具有水平不透水层;②潜水流为一维流;③库水位降落前, 原始