观测井内水位方程幂级数解答.pdfVIP

观测井内水位方程的幂级数解答 葛国昌 河海大学岩土工程研究所，江苏南京（210098 ） E-mail:ggc2021@ 摘要：针对地下水位方程形式多样的特点，利用任意函数均可进行幂级数展开的性质，建立 了观测井内水位方程的幂级数解答通式。为说明算法思想，给出了该算法在地下水位的线性 和正弦函数变化中的应用算例。分析结果表明，本文算法用于观测井内水位方程的求解具有 适应性强、计算简单、精度高等优点。 关键词：观测井，水位方程，响应滞后，幂级数解答，近似多项式 1. 引言 随着土力学的发展及其在岩土工程勘察、设计与施工中的广泛应用，地下水位的控制问 题已越来越受到工程界的重视[ 1,2 ,3 ,4 ] 。鉴于地下水位的诸多影响因素及其不可视的特点，长 期以来，人们主要利用敞开型测压管或观测井对地下水位的动态变化进行长期的追踪监测， 以研究实际工程中所发生的渗流问题。但由于监测装置自身响应滞后的影响，通常使得观测 井内水位观测值并不能准确反映地下水位的真实变化特性。为此，Hvorslev[5 ]基于达西定律 和流量连续条件建立了考虑响应滞后影响的水位平衡微分方程，并针对地下水位随时间t 的 常水头、线性以及正弦函数变化，推导了考虑响应滞后影响的观测井内水位方程。显然， Hvorslev仅是分析了土层内比较理想状态下地下水位的变化形态，而在实际工程问题中，由 于地下水位的动态变化通常受到大气降雨、井点降水、时效、温度等外界因素[6 ,7 ] 的影响， 使得土层内地下水位的实际流态复杂多样。若是针对地下水位的每种变化形态，均采用求解 水位平衡微分方程的方式建立观测井内水位方程，显然过于繁琐，而对于复杂函数形式的地 下水位方程而言，还将面临无法求解的难题。为此，本文利用任意函数均可进行幂级数展开 的性质，建立了观测井内水位方程的幂级数解答通式，实用中，通过将地下水位原函数幂级 n 数展开后的 次替代多项式代入解答通式，即可得到观测井内水位方程。 2. 观测井水位平衡微分方程 考虑任意工况，设初始时刻，地下水位（高程）为g ，监测水位为m ；在任意时刻t， 0 0 地下水位为W ，监测水位为W ，如图 1 所示。为建立观测井内外水位平衡微分方程，假定： g m 土体均质、各向同性；渗流服从Darcy定律；土、水均不可压缩；忽略井壁的摩擦和井损。 - 1 - t = 0 g0 m0 sg sm W W m t 时刻 g 图 1 观测井内外水位状态 任意时刻 t，在内外水位差的作用下，通过观测井的渗流流量为 q kF (W =−W ) （1） m g 式中，k为土体的渗透系数；F 为观测井进水段的进水系数，又称形状系数，其物理意义为单 位渗透系数的土体当中在单位水头差作用下流入或流出观测井的水流量，它取决于进水段的 形状、尺寸、结构和位置[5,8,9 ] 。 设观测井进水段的内横截面积为A ，则流量连续条件可表示为 qdt =−AdWm