基于马斯京根参数K、X的等效河道计算方法研究.pdf

年 月水利学报 第 卷第 期文章编号：（）收稿日期：；网络首发日期：网络首发地址：?基金项目：国家重点研发计划项目（）；国家自然科学基金项目（）作者简介：王船海（），教授，主要从事水文物理规律模拟及水文预报研究。：通信作者：陈钢（），副教授，主要从事水文物理规律模拟及水文预报研究。：基于马斯京根参数 、的等效河道计算方法研究王船海，俞悦，吴金宁，陈钢，马腾飞，曾贤敏，（河海大学 水文水资源学院，江苏 南京 ；河海大学 水文水资源国家重点实验室，江苏 南京 ；江苏省水文水资源勘测局常州分局，江苏 常州 ）摘要：河道汇流演算本质上属于水力学问题，为了简便计算发展了一系列的水文学方法，马斯京根法便是其中最经典的方法之一。本文基于马斯京根法中的参数 、，与特征河长相结合，构建了不同类型的等效河道（矩形、抛物线型）断面下河宽、河道比降和最大水深的计算公式，并在实例中模拟出河道各断面不同时刻的水位、流速、流量等信息。通过在实际河道中验证，结果表明，该方法实现了河道大断面资料的扩充，且具有较高的精度，改进了传统的水文学河道汇流计算模式。该方法将传统水文方法与水动力方法结合起来，有利于实现少资料地区水位、流速过程模拟，为传统的马斯京根法转换为水动力模型计算提供新方法。关键词：马斯京根参数；河道汇流；等效河道；水文水动力中图分类号：文献标识码：?研究背景河道汇流演算属于水力学问题 ，但由于对河道实测断面资料的要求高、计算耗时长等问题，水文学家们研究和发展了一系列的水文学方法和简化的水动力学方法来进行演算 。马斯京根法便是最经典的水文学方法之一。水动力学方法 则主要是通过求解圣维南方程组，但由于圣维南方程组属于复杂的偏微分方程，求解过程较为复杂 。从首次面世以来，马斯京根法经历了从经验性方法上升为具有物理基础方法的发展历程 。在长期发展的过程中，研究者们在传统线性方法的基础上，发展出一系列的非线性方法，如马斯京根 康吉法 、可变参数马斯京根 康吉法 等。马斯京根法计算过程简便，所需资料少，且精度能满足一般科研、工程需求。目前，马斯京根法在我国已经取得了广泛应用，利用马斯京根法进行汇流演算的成果非常丰富 。同时，经过广大学者的不断探索，对圣维南方程组进行求解的水动力学模型发展已经日趋成熟，能够快速有效地模拟河道流量、水位、流速等要素。“等效河道”是一种将复杂的实际河流简化为具有相似水动力特性的“理想河道”的概念，能够将复杂的自然河道形态简化为等效的理想河道形状和参数，与原河道具有相似的水动力特性的虚拟河道。早在 年 等 就已提出了等效河道的概念，然而由于河道水力特性复杂，等效的河道需要考虑的要素过多。目前国内外对等效河道的研究较少，高学平等 通过在河道内分区设置植物，对含植物河道等效创面阻力进行了试验性研究；龚定等 根据原河道大断面形状对河道断面进行了梯形概化。可以看出，关于等效河道的研究大部分是采取的试验性研究，或是对河道断面直接进行简单概化处理，而通过简便的参数实现等效河道的公式推导的研究基本没有。随着国家数字孪生流域建设中“四预”的新要求，对洪水预报提出了新的挑战，包括预报河道水位、流速等信息。目前，对马斯京根法的相关研究大多聚集在对 和 参数值更精确的获取 ，以及马斯京根法演算结果精度的进一步提高中 。传统的马斯京根法输出结果类型单一，对河道内其他要素的获取需更多资料支撑，无法与水动力模型相结合进行运用。因此，本文提出了一种基于马斯京根参数 、的等效河道计算方法，通过构建等效河道，可求解出资料缺乏地区河道的河宽、河道比降、最大水深信息，通过水动力模型进行演算后实现了河道水位、流速、流量等要素的模拟。计算原理 传统马斯京根法水量平衡方程与槽蓄方程联立：()()()（）得：（）其中：；式中：为上断面流量，、分别为时段始末上断面流量，?；为下断面流量，、分别为时段始末下断面流量，?；为河段槽蓄量，、分别为时段始末河段槽蓄量，；为计算时段；和 为马斯京根法主要参数。为蓄量常数，表示稳定流时河段的传播时间，随流量的变化而变化；为楔蓄系数，由河道楔蓄因素和调蓄因素组成。的计算公式如下：（）式中：为特征河长，；为总河段长，。总河段总、总与分河段、关系如图 所示，将一个长河段分为 个子河段，建立马斯京根法总河段总、总与分河段、关系。总（）总 ()（）假设每个分河段的 ，则：总 总（）式中：总、总为马斯京根总河段参数；、为马斯京根分河段参数。图 马斯京根分河段示意 、与水动力模型河道要素转换关系结合特征河长，利用马斯京根法分河段参数、建立等效河道要素的转换关系。已知曼宁公式：槡（）波速 河长传播时间关系式：（）特征河长计算公式：（）令水力半径 （），联立式（）（），得到水力半径 、河道比降与马斯京根法分河段参数之间的关系：槡 （）断面面积 （）。将曼宁公式与流量公式（）联立，得洪峰流量 的表达式：（）槡（）对式（）中的 偏导得：槡 （）（）()（）通过联立式（）和式（），构建等效河道。将式（）代入式（），得：槡 （）（）()（）联立式（）与式（），得：槡 （）（）()()（）建立河道比降与马斯京根参数之间的关系。由式（）得到以下变形：（）槡（）将式（）代入式（）得：（）()（）（）()（）式中：为流速，?；为糙率；为稳定流时的波速，?；为平均流速转换波速系数；为河道比降；为分河段长，；为水深对洪峰流量的偏导；为特征河长，；为断面最大水深，；为洪峰流量，?；为水力半径，；为水力半径对水深的偏导。由式（）、式（）、式（）得到通用断面形状的马斯京根法参数、与等效河道要素之间的转换关系。具体河道断面计算方法 矩形断面如图 所示，假设河道断面形状为矩形时，经简化，水力半径 ，断面面积（），（），代入式（）得河道比降：（）式中：()（）将 带入式（），得断面最大水深：槡()（）由式（）得河宽：槡（）由式（）（）（）得到矩形断面的马斯京根法分河段参数、与水动力模型河宽、河道比降、最大水深 之间的转换关系。抛物线型断面如图 所示，当河道断面概化为抛物线型时，河宽 槡，则过水断面面积：（）槡 （）图 矩形断面示意图 抛物线型断面示意则 （）槡，。为表征抛物线曲线系数，需要另外求得。代入式（）得河道比降：（）式中：()（）将 带入式（），得断面最大水深：()槡（）由式（）得 的计算公式：()槡（）则河宽：()槡（）由式（）（）（）得到抛物线型断面的马斯京根法分河段参数、与水动力模型河宽、河道比降、最大水深 之间的转换关系。根据、与等效河道要素的转换关系，通过马斯京根参数、总河段长 、洪峰流量 、糙率 、平均流速转换波速系数 及分河段数 ，可以计算得到不同概化断面下等效河道的河宽 、河道比降、最大水深 等河道信息。本文只提出了抛物线型和矩形这两种概化断面下的转换关系，根据相同的原理，其他形状的概化断面（如三角形、梯形等）也可推导求出。基于 、等效河道的具体应用及合理性分析 基于 、等效河道的具体应用何惠等 提出了用最小二乘法对马斯京根法的 、参数进行最优估计，并 以 年 海河 流 域 南 运 河 称 沟 湾 至 临 清 段 的 一 次 洪 水 过 程 为 例，其 中 、。现以该方法率定的 、值运用于本文提出的等效河道方法，对该场次洪水进行演算。南运河位于河北省南部，南运河称沟湾至临清段全长 ，如图 （）所示。本文提出的等效河道方法对该河段的概化如图 （）（）所示。整个河道概化为一条长 的长河段，取 为 ，将河段等距离划分为 个子河段，取糙率 为 ，洪峰流量 为?。经计算，当河道断面概化为矩形时，河道全程比降为 ，河宽为 ；当河道断面概化为抛物线型时，河道全程比降为 ，河宽为 。、值对应的等效河道信息，转换后可用于水动力模型的边界条件。由临清断面形状信息，结合曼宁公式及流速流量关系可推出临清断面水位流量关系。取南运河称沟湾实测流量为上边界条件，临清断面的水位流量关系为下边界条件进行模拟。图 实际河道及下边界水位流量关系图 和表 为出口临清断面实测流量、参考文献 方法演算流量及本文方法演算流量对比。从表 统计结果来看，矩形断面下，均方根误差 为 ，比文献方法大 ；可决系数 为 ，与文献方法相近，均在 以上。抛物线型断面下，均方根误差 为 ，比文献方法大 ；可决系数 为 ，与文献方法相近，均在 以上。本算例下，两种概化方式演算结果相差不大，矩形的概化方式表现更好。从图 流量过程线来看，本文提出的两种概化方法演算的临清流量结果较好，与文献 方法演算结果相差不大，能反应临清出流过程。表 秤钩湾临清流量计算结果对比指标文献 演算临清流量本文方法演算临清流量（矩形）本文方法演算临清流量（抛物线型）注：可决系数 （）（）?（），式中 是出流的平均流量；均方根误差 （）（）槡本文方法演算出不同断面处的流速情况对比如图 ，可以看到，两种概化方式下，临清断面处流速情况与图 中的流量过程线一致，计算结果具有可靠性。由图 中不同断面的流速情况可以看出，矩形断面下，各断面流速于 月 日达到最大，其中临清断面流速为?。同时，抛物线型断面下，各断面在 月 日达到最大流速，其中临清处流速最大，为?。图 称沟湾临清流量过程线（年）图 不同等效河道断面处流速情况 合理性分析本文提出的方法使用的是洪峰流量，而非入流过程，不同流量下相应的河道概化断面也会出现偏差，选用洪峰流量进行演算，是考虑到洪峰流量是最值得关注的特殊情况。因此，不同入流及时间变化都是非线性分析的重要考虑因素。从表 和图 可知，当输入的洪峰流量分别为 ，?，概化断面为抛物线型时，种拟合效果均较好，均趋于 ，且 均小于 。但不同的洪峰流量在计算得到的洪水涨落过程也呈现出一定的差异性，体现了非线性的特征。从这三种计算结果来看，当洪峰流量为?时，误差最小。表 不同洪峰流量 下计算流量结果对比指标?表 临清断面水位计算结果对比指标矩形概化断面抛物线型概化断面 为了进一步验证本文方法计算水位、流速等信息的准确及合理性，将本文提出的两种概化河道计算出的结果与临清断面实测日平均水位进行对比。如表 所示，矩形断面下，均方根误差 为 ，可决系数 为 ；抛物线型断面下，均方根误差 为 ，可决系数 为 。如图 所示，本文提出的两种概化方式均能在一定程度上反应水位的变化情况，趋势一致。计算不同洪峰流量下，两种河道断面概化方式对水位的计算结果如表 和图 所示。可以看到，不同的洪峰流量对水位的计算结果呈现出一定的非线性特征，但从 和 这两个指标来看，模拟结果表现优秀，均小于 ，均趋于 。从水位过程线中可以清晰的看出，不同洪峰流量下，演算得出的结果差异较小。总的来看，相较于抛物线型断面，矩形断面计算结果更优。图 不同洪峰流量对应的模拟结果图 临清断面实测水位与计算水位对比表 不同洪峰流量下临清断面实测水位与计算水位结果对比指标矩形抛物线型?图 不同洪峰流量下临清断面实测水位与计算水位过程实测的河道大断面资料如图 所示，可以看到，称沟湾测站过水断面形状较不规则，主槽间距为 。在最大水位为 时，本文方法概化的矩形河道断面河宽 ，抛物线型河道断面河宽 。可以看到，这两种概化方式都十分对称，是一种理想化的河道概化方式，且大致能够描述河道断面。值得一提的是，本文提出的概化河道与实测大断面仍存在一定差距，这主要由于本文提出的概化方法使整个河道每个断面都均化为唯一的断面形状，称沟湾断面只是其中一个断面，因此存在一定差距是合理的现象。图 称沟湾站实测大断面形状与概化断面对比 结论本文提出了基于 、等效河道的计算方法，通过马斯京根法中参数 、构建等效河道，推演出不同的概化河道断面（矩形、抛物线型）下的河宽、河道比降和最大水深计算公式，并根据参考文献中率定好的参数在实例中进行验证。结果表明，该方法在传统马斯京根法的基础上，通过上断面流量、马法参数 、值、河道长度 、糙率 、平均流速转换波速系数 和洪峰流量 ，就可以推演出概化河道断面的河宽 、河道比降和最大水深 推导公式。利用这些信息，除了可以通过水动力学的方法对资料缺乏地区的河道进行流量预报外，在进一步的应用中，任意时刻、河道任意断面水深、流速信息也能在水动力模型中模拟出来。总体来说，本研究实现了河道大断面资料的扩充，改进了传统的水文学河道汇流计算模式，有利于实现河道水位、流速过程模拟，为传统的马斯京根法转换为水动力模型计算提供了新方法，具有物理意义明确、计算方法简便、计算精度高等优点。值得注意的是，传统的马斯京根法的分河段是为了解决长河段下出现的非线性的问题，而本文提出的分河段则是转换后用于水动力模型计算的，可以为任意值。本研究主要提出的是线性马斯京根法下的等效河道概化方法，但由于实际河道水力要素复杂，往往都会或多或少的受到非线性因素的影响，、参数的值会随着 的变化而变化，因此非线性马斯京根法转换为等效河道的方法也值得研究，以进一步减小计算误差，并缩小概化断面形状与实测断面形状之间的差距。参考文献：龚定，李致家，臧帅宏，等河道汇流演算方法研究与应用 河海大学学报（自然科学版），：（）：芮孝芳水文学原理 北京：中国水利水电出版社，：关志成，吴海龙，崔军，等马斯京根法洪水演进反演计算方法的探讨 水文，（）：芮孝芳，张超 法的发展及启示 水利水电科技进展，（）：杨开斌，许浩，李天庆，等马斯京根系列方法在金沙江“”白格堰塞湖溃决洪水演进中的应用 水电能源科学，（）：王浩骅，管光华，肖昌诚一维圣维南方程差分数值算法中稀疏矩阵求解方法比较及优选研究 灌溉排水学报，（）：胡庆云，王船海圣维南方程组 点线性隐格式的稳定性分析 河海大学学报（自然科学版），（）：芮孝芳，张超 法的发展及启示 水利水电科技进展，（）：（），（）：，（）：王浩正，邱依婷，韩冠宇，等基于马斯京根方法的概念模型在排水系统中的应用 中国给水排水，（）：阚光远，刘志雨，李致家，等新安江产流模型与改进的 汇流模型耦合应用 水科学进展，（）：冯艳，陈理想，王加虎改进的马斯京根法在西辽河的应用 东北水利水电，（）：，：，高学平，吕建璋，孙博闻，等含植物河道等效床面阻力试验研究 水利学报，（）：，邹强，何小聪，喻杉变指数非线性马斯京根模型及其参数率定方法 水电能源科学，（）：程银才，王军，李明华基于改进模式搜索法的马斯京根法模型参数率定方法 水电能源科学，（）：，郑金德，粟忠，范玉洁，等基于改进马斯京根法漓江上游河道流量演算模型探究 绿色科技，（）：，罗宇轩，陈华，林康聆，等变参数非线性马斯京根分段演算模型研究与应用 人民长江，（）：李致家，何蒙，闫凤翔，等河道洪水演算方法在淮河及滹沱河中游的应用 河海大学学报（自然科学版），（）：何惠，张建云马斯京根法参数的一种数学估计方法 水文，（）：，（，；，；，）：，（），：；（责任编辑：于福亮）