地下水模型介绍.ppt

Build Computer Model 设置流动模型参数: 边界条件 定水头边界 抽水井边界 第三十页，共七十七页，2022年，8月28日 GW-recharge (Source 10-4 m/d) Wells 1000 m3/d Fixed head 32.0 m (1st BC) Fixed head 46.0 m (1st BC) 初始水头: demo_head_ini.trp 渗流系数 : conduc2d.trp 含水层底板标高 : bot_san1.trp 含水层顶板标高 : demo_relief.trp 含水层顶板 Aquifer top 含水层底板 Aquifer bottom 渗流系数 Conductivity m/s 定水头边界 (1st BC) 模型数据 Model data 地下水补给 定水头边界 (1st BC) 地下水抽水井 2D 地下水流动模型 2D Model scenario flow model 第三十一页，共七十七页，2022年，8月28日 Boundary conditions flow: Sign Boundary condition Geometry Unit Head (1st kind) – Dirichlet boundary condition -> fixed hydraulic head node m Flux (2nd kind) – Neumann boundary condition -> in/outflux line (2D) area (3D) m2/d m/d Transfer (3rd kind) – Cauchy boundary condition -> groundwater/surfacewater interaction line (2D) area (3D) m Well (4th kind) -> nodal extraction/injection node m3/d Note: By flux BC + values mean outflow fluxes, whereas – values correspond influxes By well BC + values for pumping wells and – values for injection wells 水流模型边界条件类型: 第三十二页，共七十七页，2022年，8月28日 3rd Kind Boundary Condition (Cauchy Type) q = ? ? (hOW – hGW) q (Infiltrations-/ Exfiltrationsrate) [m/d] ? Transferrate (Leakage-Factor) [10-4 1/d] hOW Water table in River hGW Water table in Groundwater For Transferrate ? : ? = kf / d kf conductivity of colmation bed layer d thickness of colmation bed layer. 第三十三页，共七十七页，2022年，8月28日 3rd Kind Boundary Condition (Cauchy Type) 第三十四页，共七十七页，2022年，8月28日 Constraint Conditions Example: Well with minimum hydraulic head constraint Qbc [m3/d] hmin [m] Q2 [m3/d] < Qbc h [m] = hmin Boundary condition (4th kind) Constraint condition Resulting flux Boundary condition (1st kind) Note: the algebraic signs for constraint conditions are contrary to the boundary conditions, i.e., influxes are positive, outfluxes negative! t0 = 0 t1 > t0 第三十五页，共七十七页，2022年，8月28日 输入、内插流场模型参数: 渗流系数k值 : conduc2d.trp 含水层底板标高: bot_san1.trp 含水层顶板标高: demo_relief.trp Poros