北方碳酸岩类地区库区渗漏计算的数学模型.pdfVIP

北方碳酸岩类地区库区渗漏计算的数学模型 臧逸中1 刘满杰2 贾秀梅1 贾国臣2 1， (1中国地质科学院水文地质环境地质研究所，石家庄05005 2中水北方勘测设计研究有限责任公司勘察院，天津，300222) 摘要：根据中国北方碳酸盐岩地区的库区条件，在模型中考虑了多渗透层系统的渗流，考虑了河水和库水与地下 水的垂向交换，给出了符合北方碳酸岩类地区库区渗漏特点的数学模型，并以黄河龙口水利枢纽为实例，进行了不 同条件下库区渗漏和岩溶水运动的模拟计算。 关键词：地下水数学模型库区渗漏碳酸盐类地层 北方碳酸岩类含水层往往被河谷深切，河谷地带水库库盆多由强渗透含水岩层构成，库区渗漏 与岩溶水的关系十分密切，水库蓄水会对区域岩溶水运动产生影响，在进行这些库区的渗漏评价时， 建立恰当的岩溶裂隙水数学模型，用数值模拟方法研究库区渗漏和岩溶裂隙水运动，是一个理想的 选择。 1 数学模型 中国北方碳酸盐类地层主要含水空隙为构造裂隙，其次为溶隙、岩溶洞穴，赋存有丰富的岩溶 裂隙水。碳酸盐类地层受到岩性、区域地质和构造、具体构造形迹、挽近时期地壳运动、古岩溶发 育史、地形地貌、尤其是地下水流作用的控制，岩溶发育程度由补给区、迳流区到排泄区逐渐增强。 岩溶裂隙水有统一的区域性地下水流场，会形成由迳流区到排泄区的岩溶水强径流带、和排泄区的 岩溶水集中排泄带，水流运动符合达西定律。 用达西方法分析岩溶裂隙水时，意味着是在用一理想的孔隙介质的代表性连续体，来取代空隙 为构造裂隙、溶隙、岩溶洞穴的岩石的真实集合体。对于区域地下水系统而言，将裘布衣及福尔海 默流线为水平的而且等势线垂直底板这一假设作进一步的近似，认为流线顺系统含水层层面展布， 而且等势线垂直含水层系统的底板，层状介质水流的达西定律则可以由含水岩层单位宽度流量来定 义。潜水含水层单位宽度流量g由下式给出 PK(h-d)警倒 (1)、7 式中：K为含水层渗透系数L／T；h为含水层水力水头L；d为含水层底板标高L；dh／dl为水力梯度。 对于一个饱和厚度为常数脚(含水层厚度)的承压水含层来说，其单宽度流量q由下式确定 ’ g一胁等 (2) 在碳酸盐岩地层中，夹有一定厚度的泥灰岩、页岩，其透水性相对较弱，形成相对弱透水层， 使岩溶发育呈多层结构。水库蓄水后，随着库水位大幅度上升，原来处于地下水面以上的含水岩层 将相继淹没，尤其是一些强渗透含水岩层和强岩溶发育带随著水相继饱水后，渗透性明显转强。此 时，系统应由顶部为潜水，顶部以下为承压水的n个不同渗透性能的含水层来表述，多渗透层系统 第一作者简介：臧逸中，1949年生，研究员，地下水数值模拟研究，zangyizhong@sina．corn 943 的单位宽度流量q由下式给出 q=-(K。(向一d。)+∑K，％)等 (3) 式中：K为第i层渗透系数L／T；m，为第i层饱和厚度L；磊为第n层底板标高L。 通常情况下，岩溶裂隙水与河水的交换关系是确定的，补排量近似为常量，可以处理为模型的 输入或输出。水库蓄水后，岩溶裂隙水与河水和库水的水位差加大、通过河谷地段和库盆的补排量 增加或发生变化，且遂水库运行水位的改变而不断变化，已不再是常量，不宜再简单处理为模型的 输入和输出。此时，应该建立补排量与河水和库水水位差的函数关系，这一点尤为重要。河水(或 库水)与岩溶水垂向交换关系式可表示为 qz=A(h’一矗) (4) 式中：h为含水层水位L；h为河流或水库水位L；A为越流系数1／T。 根据上面对多渗透层系统的单位宽度流量、以及地表水与地下水垂向交换关系的定义，根据水 流连续性原理和水均衡原理，忽略岩溶水系统释水的延时作用，岩溶裂隙水系统的研究库区渗漏的 地下水数学模型可有如下