非稳态井模型的理论推导及应用.docxVIP

非稳态井模型的理论推导及应用

非稳态井模型是描述地下水或油气在多孔介质中流动时，井的流量、水头（或压力）随时间动态变化的理论模型，是水文地质学、油气开发工程等领域分析流体运移规律的核心工具。与稳态模型（不考虑时间因素）不同，非稳态模型更贴近实际工程中流体流动的瞬时性，其理论推导和应用对资源开发、环境评价等具有重要意义。

一、基本假设与核心定律

非稳态井模型的推导建立在一系列简化假设和经典渗流定律的基础上，这些假设为理论分析提供了可行性，同时也明确了模型的适用范围：

均质各向同性假设：含水层或储层的渗透性能（渗透系数）在空间上均匀分布，且在各个方向上相同。

径向流假设：流体在井周围仅沿径向流动（从远处向井中心汇聚），忽略垂向流速，简化为平面径向流问题。

完整井假设：井贯穿整个含水层厚度，且井壁全部进水，保证水流在含水层全厚度上均匀汇入井内。

无限含水层假设：含水层在水平方向无限延伸，边界对流动的影响可忽略（适用于抽水初期或远离边界的情况）。

弹性储水假设：含水层中的水和介质本身具有弹性，水头变化会引起介质压缩和水的膨胀，从而释放或储存水（储水效应）。

核心定律包括达西定律和连续性方程：

达西定律：描述渗流速度与水头梯度的关系，表达式为v=-K\nablah，其中v为渗流速度，K为渗透系数，h为水头，负号表示水流从高水头向低水头流动。对于径向流，可简化为v_r=-K\frac{\partialh}{\partialr}（r为径向距离）。

连续性方程：基于质量守恒，流体流入与流出微元体的质量差等于微元体内质量的变化率，表达式为\nabla\cdot(\rhov)+\frac{\partial(\rhon)}{\partialt}=0（\rho为流体密度，n为孔隙度，t为时间）。在不可压缩流体和弱压缩介质假设下，可简化为与水头相关的形式。

二、非稳态渗流基本方程的推导

1.变量定义

导水系数T=KM（M为含水层厚度），表征含水层导水能力；

储水率S_s：单位体积含水层水头降低1单位时释放的水量；

储水系数S=S_sM：单位面积含水层水头降低1单位时释放的水量（无量纲）。

2.基本方程推导

结合达西定律和连续性方程，对平面径向流场景（仅考虑r和t的变化）进行推导：

渗流速度与流量关系：通过半径r的圆柱面的流量Q满足Q=2\pirMv_r，代入达西定律得Q=-2\pirT\frac{\partialh}{\partialr}（定流量抽水时Q为常数）。

连续性方程简化：对于非稳态流动，单位时间内圆柱体内释放的水量等于流量随半径的变化率，即2\pirMS_s\frac{\partialh}{\partialt}=-\frac{\partialQ}{\partialr}。

联立方程：将流量表达式代入连续性方程，整理得：

\frac{\partial^2h}{\partialr^2}+\frac{1}{r}\frac{\partialh}{\partialr}=\frac{S}{T}\frac{\partialh}{\partialt}

该方程为非稳态径向渗流的基本微分方程，形式与热传导方程一致，属于抛物线型偏微分方程，描述了水头h随径向距离r和时间t的变化规律。

三、泰斯模型：无限含水层中非稳态井的经典解

泰斯模型（TheisModel）是1935年由泰斯（C.V.Theis）提出的无限含水层中定流量抽水的非稳态解，是最基础的非稳态井模型。

1.定解条件

初始条件：抽水前含水层水头均匀分布，即h(r,0)=h_0（h_0为初始水头）。

边界条件：

井壁处：当r\to0时，流量为常数Q（定流量抽水）；

无限远处：当r\to\infty时，水头不受抽水影响，即h(\infty,t)=h_0。

2.方程求解

引入降深s(r,t)=h_0-h(r,t)（水头降低值），基本方程转化为：

\frac{\partial^2s}{\partialr^2}+\frac{1}{r}\frac{\partials}{\partialr}=\frac{S}{T}\frac{\partials}{\partialt}

采用拉普拉斯变换求解：对时间t进行拉普拉斯变换，消去时间导数，将偏微分方程转化为常微分方程，结合边界条件得到变换域内的解，再通过反演得到原函数。最终解得：

s(r,t)=\frac{Q}{4\piT}W(u)

其中，u