沱河地下水库蓄水率定量评价.docxVIP

沱河地下水库蓄水率定量评价 地下水库是继山区水库和水库之后的一项重要节水工程。它具有安全可靠、方便、有效等特点。它可以改善该地区的地下水环境，缓解和抑制地下水开采引起的地下水位下降和地面沉降等环境地质问题，并在石家庄形成以下水位为中心的大规模地下水位下降。2006年，裂缝面积436.5km2，裂缝中心水位埋深52.4m，为沱河地下水库提供了巨大的储水空间。为了缓解中国北方平原日益严重的供水状况，国家实施了南北工程。石家庄是中国中南部最重要的河流之一。通过调蓄水源和当地雨水洪水资源，可实现约5.40.108m3（a）的补充水量，该系统为解决区域资源性缺水问题提供了水量保障。上述研究结果表明，石家庄北部的沱河干渠吞咽能很好地入渗，在动水条件下的入渗能力约为1.5m。d。这是一座大型地下水库，具有良好的建设条件。 由于具有占地少、安全性高、不易污染和多年调节等特点,地下水库在发挥外调水源的功能方面具备明显的优势[5~6].但滹沱河地下水库位于滹沱河冲洪积扇熵,超采形成的地下水位降落漏斗为一半封闭式的储水构造,人工补给会改变区域地下水动力场,填满降落漏斗后部分补给水会发生侧向渗流,势必降低了地下水库的有效蓄水率.在分析有效蓄水率概念的基础上,利用GIS技术计算不同丰枯遭遇条件下人工补给后地下水库的有效蓄水率,并与地表水库的有效蓄水率进行对比,分析二者对储水量变化的响应规律,讨论利用地下水库进行水资源调蓄的优势,为滹沱河地下水库的设计和运行管理提供科学指导. 1 入渗能力分析 为更加精确和定量化研究地下水库的储水效率问题,本次研究类似于入渗补给系数提出了有效蓄水率的概念.将地下水库人工补给水中能够有效补给地下水并储存在库区内的那部分水量称为有效需水量,其与入渗补给水量的比值定义为有效蓄水率.计算公式为 式中:γ为有效蓄水率(无量纲);ΔQ储为地下水的储变量,L3;Q补补为人工补给水量,L3. 有效蓄水率与入渗系数相似,都是反应地表水量能够有效补给地下水库的比例,但二者又有所不同:入渗系数反映的是入渗场或库区地表的入渗能力,即地下水库的补给功能;有效蓄水率反映的是地下水库的截留能力,即地下水库的储蓄功能. 一般地下水库的补给水源有:拦蓄雨水(Q雨水)、河流来水(Q河水)、地表水库弃水(Q弃水)、调引外来水(Q调水)、冬季河流闲水(Q闲水)、灌溉用水(Q灌水)、城市废水和污水(Q废污水)等,则地下水库有效蓄水率的计算公式可表达为 式中:β有效,地下水库的有效蓄水率;β无效,地下水库的无效蓄水率;Q蓄水,地下水库的有效需水量;Q补给,地下水库的可补给水量. Q蓄水作为地下水库储存效果的重要指标,以孔隙介质为储水空间的地下水库为例,其计算公式为 式中:V,地下水库的库容,m3;μ,重力给水度(无量纲);V*,蓄水体的体积,m3;S,地下储水空间面积,m2;H1,补给前地下水位,m;H2,补给后地下水位,m. 基于上述计算原理,在具体计算时,充分考虑到研究区含水层介质的非均质性,基于GIS技术,将研究区进行剖分,分别计算每个栅格补给前后的水位差,分别乘以该栅格的给水度和面积,最后加和得出地下水库的有效补给量. 2 不同丰枯条件下黄壁庄水库的有效蓄积效应 以滹沱河地下水库为例,建立地下水库的地下水流数值模型,预报不同丰枯条件下地下水位的恢复效果,对比分析地下水库的有效蓄水率和其上游黄壁庄水库不同入库水量条件下地表水库的有效蓄水率,总结有效蓄水率的变化规律,探讨各因素对有效蓄水率的影响机理. 2.1 数值模型的建立 本次研究模拟范围以滹沱河地下水库范围为依据,如图1所示,总面积913.5 km2.(1)虽然中更新统含水层局部具有微承压性,但与上层全新统-上更新统含水层具有密切水力联系,将二者概化为一层潜水含水层;(2)西部边界为山区-平原分界线,定为二类边界,北、东、南三侧由观测孔控制,为一类边界;(3)含水层介质特征概化为非均质各向同性,地下水流概化为平面二维非稳定流. 根据建立的水文地质概念模型,构建描述计算域水位分布的数学模型,地下水库库区地下水流的数学模型用如下微分方程的定解问题来描述. 其中:K,渗透系数,m/d;H1,水头(高程),m;B,底板标高,m;ε1,补给强度,m/d;ε2,排泄强度,m/d;μ,储水系数(给水度);H0,初始水位,m;D,计算区范围;Γ1,一类边界;Γ2,二类边界;φ,一类边界上的已知水头;q,单位宽度的侧向补给量. 上述偏微分方程表示潜水含水层数学模型,与初始条件和一类、二类边界条件共同组成滹沱河地下水库地下水流的定解问题.应用有限差分法将上述数学模型离散为有限差分方程组,利用GMS4.0(Groundwater Model System)软件包中的Modflow求解.由于2003年底黄壁庄水库副坝防渗加固工