地下水建模方法和步骤剖析.ppt

? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 2 2 3 3 2 2 1 0 1 1 1 1 1 2 1 3 1 2 1 1 / / ... / / / ... 2 1 2 1 ... ... ... 2 1 2 1 2 1 n N n N e n N n N e n N n e n n e n n n e n n N n N n n n h t h t h t h t h h t h h h h h h ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4 / 1 ) 200 ( 004 . 0 25 . 0 20 8 ) ( 2 2 ? ? ? ? ? ? ? ? x t KM e ? ? 在这个例子中， ... 0,1,2,3,4, m, 10 , m 20 5 0 ? ? ? n h h n n m 10 , m 10 , m 10 , m 10 0 4 0 3 0 2 0 1 ? ? ? ? h h h h ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 5 4 3 2 1 0 1 1 4 1 3 1 2 1 1 25 . 0 25 . 0 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 n n n n n n n n n n h h h h h h h h h h 解：隐式格式一般方程为 于是有 根据初始条件得 根据边界条件得 ... 0,1,2,3,4, m, 10 , m 20 5 0 ? ? ? n h h n n m 10 , m 10 , m 10 , m 10 0 4 0 3 0 2 0 1 ? ? ? ? h h h h ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 5 0 4 0 3 0 2 1 0 0 1 1 4 1 3 1 2 1 1 25 . 0 25 . 0 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 h h h h h h h h h h 由初始条件和边界条件 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 5 . 12 10 10 15 10 25 . 0 10 10 10 20 25 . 0 10 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 25 . 0 0 0 25 . 0 2 1 4 1 3 1 2 1 1 h h h h 由此解得 t1 时刻的水头值为 m 10 , m 10 , m 10 , m 10 1 4 1 3 1 2 1 1 ? ? ? ? h h h h 在上述方程中取 n=0, 可以得到计算 t1 时刻水头值的方程 所以上述方程变成 同理，可计算 t2 时刻的水头值 2. 地下水数值模型建模步骤 模拟步骤 ? 建立概念模型 ? 建立数学模型 ? 数值方法及软件 （编程） ? 参数 ? 拟合模拟：模型校 正与检验 ? 参数敏感性分析 ? 预测模拟 软件 一、概念模型（模型概化） 根据详细的地形地貌、地质、水文地质、构造地质、水文地球化学、 岩石矿物、水文、气象、工农业利用情况等 ? 模拟的区域 : ? 含水层类型 : 潜水（无压）、承压、混合、多层 ? 维数 : 一维、二维、三维 ? 水流状态：稳定流 / 非稳定流、饱和流 / 非饱和流 ? 介质状况 : ? 均质和非均质 / 各向同性和各向异性 ? 孔隙 / 裂隙 / 双重介质 ? 流体的密度差 ? 边界条件和初始条件 必要时需进行一系列的室内试验与野外试验