论文部分内容阅读
滨海地区海水入侵含水层将引起地区性水质恶化,导致淡水水资源短缺,同时引发土壤盐渍化及荒漠化等一系列生态环境问题,成为制约滨海地区经济社会可持续发展的重要因素。河流供应淡水量的减少、气候变化导致的海平面上升及人类活动影响使海水入侵形势愈发严峻。准确模拟及在不同情景下预测滨海地区海水入侵趋势对解决滨海地区水资源短缺与生态环境问题至关重要。在解决这一问题时,必须要考虑浓度等条件引起流体密度变化。因此,本文基于地下水动力学理论,依托分布式时变增益模型(DTVGM),构建二维变密度地下水渗流溶质运移模块,重点模拟海水入侵过程。以黄河三角洲南部广饶县为例,探讨不同气候变化情景及人类活动影响下研究区未来20年海水入侵趋势,分析降水丰枯与节水规划共同作用下海水入侵面积差异,并综合利用监测数据和模拟结果进行现状及未来海水入侵风险灾害评价,为区域水资源合理高效利用、规划与管理,提供科学依据和技术支撑。论文的主要结论如下: (1)以流场分布,速度分布、密度分布和浓度分布作为公共变量,将变密度地下水渗流模型与地下水溶质运移模型紧密耦合,利用强隐式(SIP)算法迭代求解,构建变密度地下水渗流溶质运移模块,并利用Henry问题对模块进行验证,结果合理。 (2)以分布式水文模型DTVGM模拟的土壤水入渗量作为变密度地下水渗流溶质运移模块的输入项,实现DTVGM与变密度地下水流溶质运移模块的松散式耦合。利用小清河流域水文站月流量资料及黄河三角洲南部广饶县海水入侵监测数据对耦合模型进行验证,耦合模型月流量模拟效率系数达到0.74以上,相对误差小于6.4%,相关系数大于0.86;地下水位模拟与实测资料相关系数达到0.959;cl-浓度模拟与观测数据相关系数为0.914,表明耦合模型对地表水流、地下水位及溶质运移过程模拟具有可靠性。 (3)根据气候变化情景、降水保证率以及节水规划,拟定12种不同预测情景,以2009年为基准年,利用耦合模型对研究区未来海水入侵变化趋势进行多情景预测分析。保持现状开采不变,在SERS A1B情景下,研究区2030年西北和东部地区海侵线南移,海水入侵面积分别增加了5.55km2和3.13km2。中部地区海侵线基本不变。12种预测方案中,到2030年最丰方案比最枯方案海水入侵面积少19.14 km2。 (4)最枯方案下,未来2015和2030年研究区北部除小清河河道沿线外大部分地区处于较高风险和高风险,其面积相对于基准年明显扩大,到2030年增加80%,而2030年高风险区比2015年高风险区面积增加20%。表明若该区持续干旱(90%降水保证率),高风险地区将随时间推移不断增加。由于淄河水库不断的补给地下水,研究区中部低风险区在2015年和2030年相比基准年扩大,面积增加10.4%。最丰方案下,研究区北部地区较高风险和高风险区域相比基准年面积减少70.3%,2030年比2015年减少43.7%,其原因是充足的雨水补给地下水,从而有效的防止海水入侵,使得研究区北部大部分地区由基准年的较高和高风险区变为了中等风险区。中部地区淄河水库周围低风险区的范围也随着时间推移不断扩大,表明河道型水库有效的防止了海水入侵。