学位论文 > 优秀研究生学位论文题录展示

海域地表水与地下水相互作用下水量与水质运移规律数值模拟研究

作　者: 王全荣
导　师: 唐仲华
学　校: 中国地质大学
专　业: 地下水科学与工程
关键词: 地表水地下水耦合数值模拟溶质运移有限差分法
分类号: P641
类　型: 硕士论文
年　份: 2010年
下　载: 484次
引　用: 0次
阅　读: 论文下载

内容摘要

地下水与地表水的交互作用既是水资源研究的核心和热点,也是陆水相互作用中重要的物理过程。国际水文科学协会(IAHS)和国际水文地质学家协会(IAH)分别于1986年和1994年将地表水和地下水相互作用的问题正式提交会议讨论,于是地表水和地下水相互作用的研究便成了水文学及水文地质学方面研究的热点问题,然而,由于各种条件的限制,现有的模拟或者以单独研究地表水为主,或者以单独研究地下水为主,而且在模拟地表水与地下水耦合作用时,常常将地表水刻画成一维或者二维流,实际上,海域地表水体的流动往往是三维的。随着经济的快速增长和计算机速度的不断扩大,对数值模拟结果精度的要求越来越高,因此,本文旨在将地表水与地下水作为一个整体,进一步研究海域地表水与地下水相互作用下水量与数值分布规律的数值模拟问题。论文的第一部分在忽略海域地表水流动性的情况下研究地表水与地下水耦合数值模拟。首先,在忽略咸淡水密度差异的条件下,采用有限差分法建立控制方程的差分格式与半解析有限差分格式,讨论有限差分格式对时间网格剖分、空间网格剖分以及初始条件选择的敏感性,同时研究有限差分格式、半解析有限差分格式与解析解的联系与区别；其次,在考虑咸淡水密度差异的情况下,研究密度差异对地下水流场的影响；最后,通过算例对比分析传统有限差分法、半解析有限差分法和解析解,揭示这些模型的应用条件以及范围。结果表明(1)半解析数值解是采用解析法消去时间项,然后采用数值法计算空间项。数值实验结果表明,空间步长对半解析数值解的影响不敏感,对于传统的数值解,时间步长不同计算结果差异很大,因此,对于海潮波动的数值模拟,时间的剖分对数值结果的影响很大。(2)在潜水含水层中,当海潮波动相对含水层厚度比较小时,模型的线性化给计算结果带来的误差较小,可以忽略不计,反之,当波动比较大时,误差随着也增大。(3)求解海潮波动引起地下水流动的问题时,由于海潮的周期性,解析解和半解析数值解都不需要初始条件,然而,采用传统有限差分法求解非稳定流问题,初始条件必不可少,数值模拟结果表明采用平均海平面作为初始条件,效果最佳。(4)储水系数对海水入侵引起地下水流场变化比较敏感,储水系数越小,忽略咸淡水密度差异造成的误差越大。对于储水系数较大的潜水含水层,海水与淡水的密度差异只影响临海地区的流场,但是在远离海域的地方,流场几乎不受影响。(5)海水的浓度随着时间的推移慢慢向内陆延伸,但是由于水头海域水头的波动很小,因此浓度的传播很慢。论文的第二部分对应于本文的第三章,是在考虑地表水流动的情况下,研究海域地表水与地下水耦合数值模拟问题,其中,地表水采用三维的Navier-Stokes方程组来描述。本文在Katerina差分格式的基础上,分析其收敛性条件以及差分方程的特点,针对该差分方程的缺陷,推导了上游加权半隐式半显式差分格式,然后,提出了改进的方法——加权因子预测法,并推导了加权因子的计算公式,最后,通过两个数值算例,验证了改进法的合理性以及有效性,并研究了含水层性质对耦合流场下水质和水量空间分布的影响。结果表明(1)现有的半显式半隐式有限差分格式虽然计算量小,然而由于其收敛性条件难以把握,很容易造成部分区域质量不守恒或者波动现象,于是,本文提出了加权因子预测法,数值实验表明该法效果很好。(2)当含水层的渗透性较差和储水系数较大时,地表水与地下水的交互作用比较弱,即地表水水位的强烈变化在含水层中的衰减太快,反之,当储水系数比较小时,地表水与地下水的耦合作用强烈,地表水水位变化在含水层中的衰减变慢。论文的第三部分对应于本文的第四章,首先建立控制方程的有限差分格式,然后通过数值算例研究污染物在含水层与地表水中的迁移规律,结果表明：与地下水流动速度相比,地表水的流动速度很快,因此,溶质在含水层的传播速度慢,而当运移到地表水体时,传播速度很快。若某地区含水层受到了污染,则尽量将污染物圈定在某个区域内,比如加隔水墙,或者在被污染的含水层附近采用低渗透的介质使污染物的扩散减慢,同时在含水层中注入某些生物化学溶剂,使含水层慢慢净化。第四部分是应用实例,采用数值模拟手段研究雷州半岛地表水与地下水耦合作用下水质和水量的空间分布规律。首先通过分析总结现有的勘察资料,弄清楚该地区的水文地质背景,并概化出该地区的水文地质概念模型,然后,针对该模型,建立对应的数学模型与数值模型,最后,采用自制模拟程序,模拟该地区的水量与水质的分布规律。数值模拟结果表明,随着海水水位的波动,含水层的水头也在波动,然而,波动的范围很小,即只在距离海域100m处以内的含水层中的水流受影响,远离海岸的地方几乎不受海水潮汐动态的影响。对于海水入侵情况,在刚开始,海域的咸水入侵含水层,然而,随着时间的推移,咸水的入侵速度变慢,到14小时后,入侵基本结束,即咸淡水界面保持稳定。因此,在雷州半岛,如果没有高强度的开采,是不会导致海水入侵问题。

全文目录

摘要  6-8
ABSTRACT  8-13
第一章绪论  13-26
  1.1 研究目的和意义  13-14
    1.1.1 研究目的  13
    1.1.2 研究意义  13-14
  1.2 国内外研究现状及存在的主要问题  14-24
    1.2.1 国内外研究现状  15-23
    1.2.2 主要存在的问题  23-24
  1.3 主要研究内容及研究思路  24-26
    1.3.1 研究内容  24
    1.3.2 研究思路  24-25
    1.3.3 研究的创新点  25-26
第二章忽略流动性的地表水与地下水相互作用的数值模拟研究  26-46
  2.1 恒定密度地表水与地下水耦合数值模拟研究  26-38
    2.1.1 数学模型  26-27
    2.1.2 传统的有限差分数值模型  27-28
    2.1.3 半解析有限差分数值模型  28-30
    2.1.4 网格剖分对数值模拟结果的影响  30
    2.1.5 初始条件对数值模拟结果的影响  30
    2.1.6 数值解的检验  30-38
  2.2 变密度地表水与地下水耦合数值模拟研究  38-45
    2.2.1 数学模型  38-39
    2.2.2 数值方法  39-40
    2.2.3 数值实验  40-45
  2.3 小结  45-46
第三章三维地表水与地下水相互作用下水流运动规律的数值模拟研究  46-84
  3.1 三维地表水与地下水耦合的数学模型  46-51
    3.1.1 地表水的控制方程  46-48
    3.1.2 地下水流的控制方程  48-50
    3.1.3 地表水与地下水耦合模型  50-51
  3.2 三维地表水与地下水耦合的数值方法  51-62
    3.2.1 地表水流数学模型的有限差分格式  51-58
    3.2.2 地下水流数学模型的有限差分格式  58-59
    3.2.3 地表水和地下水耦合数学模型的数值模型  59-61
    3.2.4 算法以及流程  61-62
  3.3 差分格式的讨论与改进  62-66
    3.3.1 对流项的有限差分格式  62-65
    3.3.2 有限差分格式的讨论  65
    3.3.3 有限差分格式的改进  65-66
  3.4 数值模拟中特殊问题的处理  66-69
    3.4.1 单元格的疏干(drying)与再充水(rewetting)的处理  66-67
    3.4.2 潜水含水层自由表面的处理  67-68
    3.4.3 水跃处理  68-69
  3.5 数值实验  69-83
    3.5.1 算例1:地表水半隐式半显式有限差分数值算法的研究  69-74
    3.5.2 算例2:地表水与地下水耦合作用对含水层参数敏感性的研究  74-83
  3.6 小结  83-84
第四章三维地表水与地下水相互作用下溶质迁移规律的数值模拟研究  84-91
  4.1 溶质运移耦合数学模型  84-86
    4.1.1 溶质运移数学模型  84-85
    4.1.2 溶质运移与水流方程耦合数学模型  85-86
  4.2 数值方法  86
  4.3 数值实验  86-89
  4.4 小结  89-91
第五章雷州半岛地表水与地下水耦合数值模拟  91-100
  5.1 区域地质概况  91-96
    5.1.1 自然地理概况  91-92
    5.1.2 地质构造  92-95
    5.1.3 水文地质特征  95-96
  5.2 概念模型以及数值模型  96-97
    5.2.1 概念模型  96
    5.2.2 数学模型以及数值方法  96-97
  5.3 结果与讨论  97-98
  5.4 小结  98-100
第六章结论及建议  100-103
  6.1 结论  100-101
  6.2 建议  101-102
  6.3 不足与展望  102-103
致谢  103-104
参考文献  104-112

海域地表水与地下水相互作用下水量与水质运移规律数值模拟研究

内容摘要

全文目录

相似论文