三峡库区库岸边坡地下水浸润线计算中的解析法.docx

三峡库区库岸边坡地下水浸润线计算中的解析法 1 库岸边坡地下水浸渍线解析法 库水波动是影响库岸变形破坏的主要动态因素。库水波动下库岸边坡地下水提取物的合理计算对水库边坡的稳定分析、预测和报告具有重要意义。 近几十年来,很多学者采用不同方法对边/滑坡地下水渗流问题进行了研究,如文献分别用有限元方法研究了降雨和库水波动下边坡地下水渗流场的变化趋势;文献用有限差分法研究了降雨和库水波动联合作用下的滑坡地下水浸润线。数值方法能较好解答复杂初始条件和边界条件下的边/滑坡地下水浸润线问题,但目前还未在工程实际中广泛使用,解析法是比数值法更容易推广应用的实用方法。文献基于含水层均质、侧向无限延伸并具水平不透水层,库岸垂直,初始地下水位水平和潜水流为一维流等假定,用解析解法研究了不考虑降雨和蒸发作用时库水位等速变化下的滑坡地下水浸润线。文献在上述假定的基础上研究了考虑降雨和蒸发影响时的情况,文献用解析解法研究了库水位呈正弦半波曲线变化时的滑坡地下水浸润线。上述解析解都是基于一定假定并对方程线性化处理后得到的,然而,用解析法计算的边/滑坡地下水浸润线经常与实际情况不符,这主要是因为没有系统地研究各种假定和线性化过程带来的误差并分析其适用条件。 结合前人研究成果,本文基于一定假定建立了不同的库岸边坡地下水渗流数学模型,并求出了库水位呈任意函数波动下的库岸边坡地下水浸润线解析解。通过对比分析不同数学模型的数值解和解析解,研究了库岸边坡垂直处理、潜水运动基本方程线性化处理、假定库水位等速变化、不考虑非饱和渗流和Dupuit假定带来的误差及其规律,从而分析计算库岸边坡地下水浸润线解析法的适用范围。 2 分析与解决 2.1 数学模型的基本假设 由于潜水运动基本方程(即布西涅斯克方程)是一个二阶非齐次偏微分方程,只有通过一定的假定并对方程线性化处理后才能给出浸润线的解析解。本文在建立用于解析解答的数学模型时,有以下基本假定: ①含水层均质、各向同性,侧向无限延伸并具有水平隔水层;②潜水流为一维流,不考虑垂直潜水面方向的渗流(dupuit假设);③ 库水波动范围的岸坡垂直;④不考虑非饱和渗流,即岩土介质的渗透性与其含水量无关;⑤由于降雨入渗和蒸发作用下岩土体中水的渗流过程十分复杂,目前的处理方法都有较大的误差,本文不考虑降雨入渗和蒸发情况。 2.2 含水层厚度变化模型 基于上述假设,可建立如图1所示的几何模型(将隔水底板高程定为0 m)。 不考虑降雨入渗和蒸发作用时,潜水一维非稳定运动的基本方程为 ??x(Κh?Η?x)=μ?Η?t。(1)??x(Kh?H?x)=μ?H?t。(1) 式中:潜水面上升时μ为重力给水度,下降时μ为饱和差;t为时间;K为渗透系数;x为位置坐标;h为含水层厚度;H为总水头(隔水底板水平时,h和H相等)。式(1)为一非线性方程,须对其线性化后才能解答,通常采用的线性化方法是用平均含水层厚度hm近似代替h。设水位传导系数a=Khm/μ,有 a?2h?x2=?h?t。(2)a?2h?x2=?h?t。(2) 设s(x,t)=h(x,t)-h(x,0),可得到以含水层厚度变化函数s(x,t)为未知函数的数学模型: {a?2s?x2=?s?t,0x∞,t0；s(x,0)=0,0x∞；s(0,t)=f(t),t0；s(∞,t)=0,t0。(3) 式中f(t)为库水位变动函数,可由库水位-时间曲线经过函数拟合得到。 2.3 平均含水层厚度的迭代计算 对于上述数学模型,通过对t进行拉普拉斯正变换,解答常微分方程,可得 S(x,p)=F(p)e-√pax,x0,p0。(4) 式中:S(x,p)为s(x,t)关于t的拉普拉斯变换,F(p)为f(t)关于t的拉普拉斯变换。运用拉普拉斯变换的微分性质和卷积定理对上式进行拉普拉斯逆变换,有 s(x,t)=L-1[F(p)e-√pax=x2a√π∫t0f(τ)ex24a2(t-τ)dτ?x0?t0。(5) 所以有: h(x,t)=h(x,0)+x2a√π∫t0f(τ)ex24a2(t-τ)dτ?x0?t0；h(x,t)=h(x,0)+f(t),x=0,t0。(6) 若假定库水位等速变化,即f(t)=vt,v为库水变化速度,上升为正,下降为负。则 S(x,p)=F(p)e-√pax=v0p2e-√pax。(7) 对上式(7)进行拉普拉斯逆变换,有 s(x,t)=L-1[v0p2e-√pax]=vt?4i2erfc(x2√at)=vtΜ(x2√at)。(8) 由于M(0)=1,h(x,t)可写成下列统一格式: h(x,t)=h(x,0)+s(x,t)=h(x,0)+vtΜ(x2√at)?x≥0,t0。(9) 其中, i2erfc(z)=2√π∫∞z∫∞z∫∞ze-u2dudzdz。(10) M(z)为一特殊函数,可