区域饱和-非饱和地下水流运动数值模拟

Saturated and unsaturated groundwater flow numerical simulation in large scale zone

根据多孔介质中区域饱和-非饱和水分运动的特征,建立基于迦辽金有限单元法的拟三维数学模型.通过对三维Richards方程沿饱和含水层垂向平均得到水平二维模型的基本方程,在各选定区域非饱和带沿水平方向平均得到多个一维非饱和水流运动方程.一维非饱和模型与二维饱和模型的耦合点为非饱和带计算所得的通过潜水面的流量项.最后,验证拟三维模型的模拟结果,说明模型合理可行,为研究大区域饱和-非饱和水分运动提供了有效的分析工具.

Based on the pattern of water flow in large scale saturated and unsaturated zones, a quasithree-dimensional finite element model was presented. By averaging 3-D flow equation over the aquifer thickness, the governing equation of the horizontal saturated 2-D model was obtained. The basic equation of the vertical unsaturated 1-D mode was received by averaging the 3-D saturated and unsaturated equation over the horizontal plane in all selected regions. The coupling point was a water flux term at the phreatic surface through which the unsaturated zone linked with the saturated one. The proposed model was testified; and it is shown that the model is accurate and computational effective.