潮汐对地下水波动影响的数值模拟[J].pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 潮汐对地下水波动影响的数值模拟 陈娟 庄水英 李凌 河海大学水文水资源与水利工程科学国家重点实验室水生态环境模拟中心 江苏南京 摘要 本文从经典的 方程出发 在考虑浸润面影响的基础上 建立了基于有限差分法的滨海地区一维地 下水运动数值模型 通过将潮汐运动概化为正弦波 模拟了滨海地区地下水水位随潮汐波动的变化 结果显示 地下水水位的波动相对于潮汐波动来说具有不对称性 振幅衰减以及相位滞后等特征 通过比较受潮汐影响的地 下水平均水位与不考虑潮汐影响的地下水水位发现 实际的地下水平均水位要高于未考虑潮汐影响的地下水水 位 最后数值模拟了浸润面的形成和消失的过程 关键词 地下水水位 潮汐 数值模型 超高 浸润面 中图分类号 文献标识码 研究背景 在自然条件下 海平面的上下波动 潮汐 影响与其毗邻的滨海地下水系统的变化 这种变化反映在 潜水含水层内表现为潜水水位的变化 而反映在承压层内则表现为承压水头的变化 等 在 年 用边界元数值方法模拟了潮汐对滨海潜水水位波动的影响 以二维饱和地下水运动方程为控制方程 分 别考虑了自由边界 潜水面 和动边界 潮汐海平面边界 并且还考虑了沙滩倾角和浸润面 得出地下水 水位受潮汐影响的变化规律 振幅衰减 不对称性 相位滞后 在潮汐边界的影响下 滨海地下水水位 经周期平均后高于不受潮汐影响的地下水稳定水位 即滨海地下水水位在潮汐影响下存在超高 通过野外实验观测 总结了形成内陆水位超高的 个原因 浸润面的形成 非垂直沙滩导致的 边界条件不对称以及控制方程的非线性因素 等 进一步证实 由于非线性效应和沙滩倾角的影响 潮汐会产生内陆地下水水位超高 并且 当海岸沙滩角度较小时 这种由于潮汐产生的超高是非常明显 的 忽略潮汐造成的地下水平均水位超高将可能导致预测评估的地下水资源量出现误差 潮汐海平面的周期变动带动相邻的地下水水位波动 然而 地下水在沙滩斜坡上的 出点 与海平面 线并不总是重合的 见图 出点 与海平面线之间的区域即定义为浸润面 浸润面大多出现在落 潮时 对沙滩泥沙运移有重要的影响 浸润面的出现更是形成地下水平均水位超高的一个重要原 因 本文通过建立基于简单的一维地下水运动方程的有限差分数值模型 考虑浸润面条件 将潮汐运动 概化为正弦波 探索研究滨海潜水水位内陆波动随时间变化的关系 同时对比研究动态稳定后周期平 均的地下水水位与稳态解水位之间的关系 探寻两者之间存在的超高 并分析形成超高的主要原因 数值模型 控制方程 如图 所示的滨海潜水含水层 海平面的上下运动 潮汐 将引起地下水水位的变化 收稿日期 基金项目 国家杰出青年科学基金资助项目 教育部博士点基金资助项目 作者简介 陈娟 女 江苏张家港人 硕士生 研究方向为水生态与水环境模拟 为简化计算 作以下假定 含水层均质 各向同性 底部水平且不透水 忽略毛细现象和上表面 土壤水渗透 潜水流动符合 假设 潜水流动的连续性方程可表示为一维非线性 方程 式中 为潜水地下水水位 以底部为基准面 为有效孔隙度 为水力传导系数 为时 间 为横向坐标 边界条件 边界为向陆边界和向海边界 其 中 向陆边界条件模拟内陆地下水与河流或水库接 壤 取为常水头边界 需要注意的是 为研究地 下水水位受潮汐的影响 须保证内陆水头不对潮汐 波的传输产生影响 故尽量使横向距离 取值较 大 同时内陆水头 与平均海平面高度 相差 不会太大 向海边界受潮汐和浸润面的影响 在几 图 潮汐对滨海地下