地下径流退水过程规律研究.docxVIP

地下径流退水过程规律研究 地下排水和蓄水的过程是指雨季较少或无雨季的连续排水和蓄水过程。这是水文过程的重要组成部分。 在实际应用中, 通常把地下水退水流量和含水层蓄量之间表示为线性关系, 从而推导出流域地下径流的退水过程, 但是不同的流域, 其地下径流的退水过程规律不一定相同, 需要从实测资料进行分析.因此, 笔者选择实测径流资料较全的滦河流域子流域 (柳河和瀑河流域) 的一组退水过程来进行分析, 确定其地下径流的退水规律. 1 关于柳河水系,柳河至柳河口, 滦河流域位于华北平原东北部, 流经内蒙古、辽宁、河北三省、自治区的27个县、旗、区, 于河北省乐亭县兜网铺注入渤海, 全长888, km.流域面积44, 750, km 柳河和瀑河是滦河的支流, 柳河发源于承德市兴隆县八拨子岭西北麓二道沟, 在柳河口滦河右岸汇入滦河, 河长86, km, 流域面积1, 020, km 2 通常，地下流量回水过程的表达方法 假定退水流量和含水层蓄水量之间的几种关系, 可以推导出不同形式的退水流量随时间的变化关系如下. 2.1 退水流量与含水层蓄水量之间的线性关系 退水流量和含水层蓄水量的线性关系可表示为 式中:Q为退水流量, m 联合连续性方程为 由此可求得退水方程为 式中Q 2.2 回收水和含水层蓄水量的指数关系是回收水和含水层蓄水量的指数关系 式中:D为含水层缺水量, m 联立连续性方程, 可求得退水方程为 2.3 退水曲线相邻流量的确定 退水流量和含水层蓄水量的幂函数关系 联立连续性方程, 可求得退水方程为 式中a、b为待定参数, 可通过以下步骤来求出. (1) 假定一个b值; (2) 在流域退水曲线上, 每隔一个时间步长Δt取2个相邻的流量Q (3) 对于退水曲线上不同的相邻流量, 就会得到一个关于a的不同估计值, 所有的估计值构成一个序列, 其均值和方差分别为a (b) 和C (4) 在0和1之间对b值进行系统取样, 根据以上步骤计算出相应的均值和方差; (5) 对应着最小的方差值的a (b) 和b被选作2个参数的最优值. 模型的精度用Powell提出的损失函数来评定, 该方法将选择的所有场次洪水的退水过程综合考虑, 能对各种模型的精度进行对比, 即 式中:r 另外, 还有退水流量和含水层蓄水量呈2个并联线性水库关系或串联线性水库关系的假定. 3 地下径流退水曲线 图形分析是对线性水库退水和指数退水而言的.如果退水曲线符合线性水库假定, 则一组退水曲线在ln Q-t图中表现为平行的直线.直线的斜率为 如果退水曲线服从指数假定, 则一组退水曲线在Q 选择柳河和瀑河的地下径流退水曲线, 原则是:从一次洪水退水曲线的转折点开始到径流开始起涨点结束, 中间时段每小时降雨量不超过0.1, mm, 这样, 从柳河和瀑河流域分别选取了4次和6次地下径流退水过程, 并绘在ln Q, -t和Q 从图2和图3可以看出, 柳河流域地下径流退水过程在ln Q, -t图上并不是直线, 而在Q 4 地下径流退水分 由于图形分析的方法不能反映非线性水库退水等复杂的退水过程, 需要采用模型模拟的方法, 对比模型模拟的精度来确定流域地下径流到底服从哪种退水规律.将前面所述的线性水库模型、指数模型和非线性水库模型应用于滦河水系子流域柳河流域和瀑河流域的地下径流退水过程模拟, 模型参数采用式 (10) 和式 (11) 计算.线性水库模型求得的柳河k值范围在0.012~0.017之间, 平均值为0.014, 瀑河k值范围较大, 在0.006, 9~0.020 0之间, 主要受降雨类型的影响.将瀑河流域的地下径流退水分为2类: (1) 强度大 (大于20, mm/h) 、历时短 (小于6, h) 的降雨对应的k的平均值为0.008, 2; (2) 强度小 (小于20, mm/h) 、历时长 (大于6, h) 的降雨对应的k的平均值为0.019. 指数模型求得柳河φ值范围在2.63×10 用非线性模型求得的柳河, β, 值为0.56, α范围在758.37~1, 022.22之间, 平均值为891.68, 瀑河α值范围较大, 在421.86~1, 569.58之间, 其取值受降雨类型的影响.强度大、历时短的降雨对应的α的平均值为1, 228.99, 强度小、历时长的降雨对应的α的平均值为563.57.模拟各次地下径流退水过程时用计算所得参数的平均值.表1和表2列出了各个模型在模拟柳河和瀑河地下径流退水过程的参数值.图4~图6以及表3是用不同方法模拟的结果. 从退水模拟结果图上并不能清楚地辨别各个模型模拟效果的优劣, 但是从模拟误差分析图上可以看出, 在柳河流域, 指数模型模拟的误差值在0左右小幅度摆动, 最大相对误差为-0.14, 比线性水库模型和非线性模型模