一种用于分布式水文模型的复式河道断面概化及计算方法技术

本发明专利技术公开了一种用于分布式水文模型的复式河道断面概化及计算方法，包括：采用多个倒等腰梯形断面对实测复式河道断面进行拟合，得到各梯形断面的上底宽、下底宽和高的值；根据实测断面数据，拟合断面面积与流域汇流面积曲线，利用该曲线及各子流域汇流面积求各子流域河道断面面积；根据模拟断面和实测断面面积比例，等比缩放拟合梯形断面参数得各子流域河道断面的参数；求每个子流域断面的各梯形断面宽高比和面积及累积梯形断面深和面积；已知过水断面面积和累积梯形断面面积，可求过水断面的水深；已知过水断面水深和累积梯形断面深，可求过水断面面积、湿周、水力半径。本发明专利技术提高对复式断面的概化精度；进而提高模型汇流过程模拟精度。

【技术实现步骤摘要】
一种用于分布式水文模型的复式河道断面概化及计算方法
本专利技术涉及水文模型领域，具体涉及一种用于分布式水文模型的复式河道断面概化及计算方法。
技术介绍
分布式水文模型在水文模型领域应用越来越广泛；分布式水文模型一般需要对实际河道横断面进行概化，采用规则的断面形状(例如梯形或三角形)进行流量、过水断面面积与水深关系的计算。相比于实测断面而言，概化断面具有参数个数少、过水断面面积与水深计算关系简单，且可以根据不同子流域汇流面积进行类推等优点。一般情况下，现有分布式水文模型都采用一个单一形状(梯形或三角形)进行河道横断面概化；然而，实际河道往往是复杂，通常是呈阶梯状的复式断面，采用单一形状并不能完全地拟合实际河道断面形状，从而影响水文模型中汇流过程的计算，尤其是水深的计算准确性。
技术实现思路
本专利技术所要解决的技术问题是现有的分布式水文模型都采用单一形状进行河道横断面概化，不能完全地拟合实际河道断面形状，从而降低水文模型中汇流过程的计算精度的问题。为了解决上述技术问题，本专利技术所采用的技术方案是提供一种用于分布式水文模型的复式河道断面概化及计算方法，包括以下步骤：步骤S10、采集水文断面实测数据，并根据实测数据表现出的河道断面阶梯状特征，采用相应数量的倒等腰梯形断面对该复式河道断面进行拟合，并得到各梯形断面的上底宽T、下底宽B、高H的值；其中，上一个梯形断面的下底宽不小于下一个梯形断面的上底宽，且两者高程相同；步骤S20、根据采集的实测水文断面面积以及其对应的流域汇流面积，拟合得到实测水文断面面积与流域汇流面积的曲线，并利用该曲线及子流域的流域汇流面积，求得无实测数据的其他子流域的河道断面面积；步骤S30、采用河流上最近的上游或下游的实测水文断面的梯形断面概化形状作为各子流域梯形断面概化形状，并根据该子流域河道断面面积与最近的上游或下游的实测水文断面面积比，对该实测水文断面的各梯形断面参数进行等比例缩放，得到该子流域断面的各概化梯形断面的参数；步骤S40、以分布式水文模型程序读取的格式存储各子流域的梯形断面的概化参数；步骤S50、对每个子流域的概化梯形断面分别求解每个梯形断面宽高比、每个梯形断面面积、累积梯形断面深以及累积梯形断面面积；步骤S60、利用已知子流域河道断面的过水断面面积Ah和累积梯形断面面积，求解该子流域对应的过水断面的水深h′；步骤S70、利用已知子流域河道过水断面的水深h和累积梯形断面深SHk，求解该子流域对应的过水断面面积A′h、过水断面湿周χh、水力半径Rh。在上述方法中，对实际复式河道断面形状进行概化的倒等腰梯形断面数量为2个以上；最下面一个梯形断面的下底宽最小取值为0。在上述方法中，在步骤S40中，存储各子流域的梯形断面的概化参数的格式具体为：第一行标识整个流域中所有子流域所采用的梯形断面数量最大值，用于在模型应用中确定数组大小；第二行为标题行，其中，第1列为子流域编号，第2列表示当前子流域断面实际拟合梯形断面的个数，第3至5列依次表示该子流域断面对应的第一个梯形断面的上底宽、下底宽、高三个参数，第6至8列表示该子流域断面对应的第二个梯形断面的三个参数，依次类推，直到记录完该子流域断面对应的所有梯形断面的参数；如果没有对应的梯形断面，则所有参数为0；第三行以下各行按照子流域编号由大到小顺序依次记录各子流域概化断面对应的数据；每列对应第二行的一个标题，记录对应子流域概化断面的一个参数。在上述方法中，子流域第i个梯形断面宽高比θi计算公式为：子流域第i个梯形断面面积Ai计算公式为：Ai＝(θi·Hi+Bi)×Hi子流域的累积梯形断面深SHk的计算公式为：子流域的累积梯形断面深SAk的计算公式为：其中，i表示同一个子流域内具体梯形断面的编号；Ti表示该子流域第i个梯形断面的上底宽；Bi表示该子流域第i个梯形断面的下底宽；Hi表示该子流域第i个梯形断面的高；k为累计的梯形断面个数，0≤k≤N，N为该子流域河道断面概化的梯形断面个数。在上述方法中，步骤S60具体包括以下步骤：步骤S61、利用已知该子流域河道断面的过水断面面积Ah与该子流域的累积梯形断面面积逐个比较，确定使SAk＜Ah≤SAk+1成立时连续累积的梯形断面个数k；步骤S62、根据求得的k值和过水断面的水深h′的计算公式求解h′。在上述方法中，过水断面的水深h′的计算公式为：在上述方法中，步骤S70具体包括以下步骤：步骤S71、利用已知子流域河道过水断面的水深h与该子流域的累积梯形断面深SHk逐个比较，确定使SAk＜Ah≤SAk+1成立时连续累积的梯形断面个数k；步骤S72、根据求得的k值和过水断面面积A′h计算公式求解A′h；步骤S73、根据求得k值和过水断面湿周χh的计算公式求解χh；步骤S74、根据求得的过水断面湿周χh和水力半径Rh计算公式求解Rh。在上述方法中，过水断面面积A′h计算公式为：A′h＝SAk+(θk+1·(h-SHk)+Bk+1)×(h-SHk)。在上述方法中，过水断面湿周χh计算公式为：水力半径Rh计算公式为：在上述方法中，在步骤S30中，如果河流上没有实测水文断面，则采用临近河流的实测水文断面的梯形断面概化形状作为无实测水文断面河流上各子流域的梯形断面概化形状。与现有技术相比，本专利技术采用多个倒等腰梯形对复式断面进行拟合，并给出不同水位下断面面积、湿周、水力半径的计算公式，用于分布式水文模型汇流演算过程。相比于现有技术的缺点和不足，本专利技术具有以下有益效果：(1)采用多个梯形断面对实际河道横断面进行拟合，提高了对复式断面的概化精度；(2)给出了不同水深条件下，断面面积、湿周、水力半径等参数的计算公式，提高水文模型对断面水位-面积关系的模拟精度，进而提高模型汇流过程模拟精度。附图说明图1为本专利技术提供的一种用于分布式水文模型的复式河道断面概化及计算方法的流程图；图2为本专利技术中具体实施例的河道断面概化示意图；图3为本专利技术中一个子流域的所有梯形断面参数信息存储格式；图4为本专利技术中求解一个子流域的概化梯形断面的梯形断面宽高比、每个梯形断面面积、累积梯形断面深以及累积梯形断面面积的常量参数表。具体实施方式本专利技术提供了一种用于分布式水文模型的复式河道断面概化及计算方法，采用多个倒等腰梯形对复式断面进行拟合，并给出不同水位下断面面积、湿周、水力半径的计算公式，用于分布式水文模型汇流演算过程，提高水文模型对断面水位与面积关系的模拟精度，进而提高模型汇流过程模拟精度。下面结合说明书附图和具体实施方式对本专利技术做出详细说明。在进行断面概化的时候，需要确保多个梯形是通过上下堆积方式进行排列，而且需要确保上面梯形的下底宽要大于等于下面梯形的上底宽。多梯形断面概化使得水深～面积等关系为非连续的分段函数关系，因此在进行断面相关参数计算时，需要首先通过对比确定已知水深或者断面面积所位于的函数段，然后再根据计算公式进行计算如图1所示，本专利技术提供的一种用于分布式水文模型的复式河道断面概化及计算方法，包括以下步骤：步骤S10、采集水文断面实测数据(如图2中粗曲线所示)，并根据实测数据表现出的河道断面阶梯状特征，采用相应数量的倒等腰梯形断面对该复式河道断面进行拟合，得到各梯形断面的上底宽T、下底宽B、高H三个参数，并存储该梯形断面信息，存储格式本文档来自技高网...

【技术保护点】
1.一种用于分布式水文模型的复式河道断面概化及计算方法，其特征在于，包括以下步骤：步骤S10、采集水文断面实测数据，并根据实测数据表现出的河道断面阶梯状特征，采用相应数量的倒等腰梯形断面对该复式河道断面进行拟合，并得到各梯形断面的上底宽T、下底宽B、高H的值；其中，上一个梯形断面的下底宽不小于下一个梯形断面的上底宽，且两者高程相同；步骤S20、根据采集的实测水文断面面积以及其对应的流域汇流面积，拟合得到实测水文断面面积与流域汇流面积的曲线，并利用该曲线及子流域的流域汇流面积，求得无实测数据的其他子流域的河道断面面积；步骤S30、采用河流上最近的上游或下游的实测水文断面的梯形断面概化形状作为各子流域梯形断面概化形状，并根据该子流域河道断面与最近的上游或下游的实测水文断面面积比，对该实测水文断面的各梯形断面参数进行等比例缩放，得到该子流域断面的各概化梯形断面的参数；步骤S40、以分布式水文模型程序读取的格式存储各子流域的梯形断面的概化参数；步骤S50、对每个子流域的概化梯形断面分别求解每个梯形断面宽高比、每个梯形断面面积、累积梯形断面深以及累积梯形断面面积；步骤S60、利用已知子流域河道断面的过水断面面积Ah和累积梯形断面面积，求解该子流域对应的过水断面的水深h′；步骤S70、利用已知子流域河道过水断面的水深h和累积梯形断面深SHk，求解该子流域对应的过水断面面积A′h、过水断面湿周χh、水力半径Rh。...

【技术特征摘要】
1.一种用于分布式水文模型的复式河道断面概化及计算方法，其特征在于，包括以下步骤：步骤S10、采集水文断面实测数据，并根据实测数据表现出的河道断面阶梯状特征，采用相应数量的倒等腰梯形断面对该复式河道断面进行拟合，并得到各梯形断面的上底宽T、下底宽B、高H的值；其中，上一个梯形断面的下底宽不小于下一个梯形断面的上底宽，且两者高程相同；步骤S20、根据采集的实测水文断面面积以及其对应的流域汇流面积，拟合得到实测水文断面面积与流域汇流面积的曲线，并利用该曲线及子流域的流域汇流面积，求得无实测数据的其他子流域的河道断面面积；步骤S30、采用河流上最近的上游或下游的实测水文断面的梯形断面概化形状作为各子流域梯形断面概化形状，并根据该子流域河道断面与最近的上游或下游的实测水文断面面积比，对该实测水文断面的各梯形断面参数进行等比例缩放，得到该子流域断面的各概化梯形断面的参数；步骤S40、以分布式水文模型程序读取的格式存储各子流域的梯形断面的概化参数；步骤S50、对每个子流域的概化梯形断面分别求解每个梯形断面宽高比、每个梯形断面面积、累积梯形断面深以及累积梯形断面面积；步骤S60、利用已知子流域河道断面的过水断面面积Ah和累积梯形断面面积，求解该子流域对应的过水断面的水深h′；步骤S70、利用已知子流域河道过水断面的水深h和累积梯形断面深SHk，求解该子流域对应的过水断面面积A′h、过水断面湿周χh、水力半径Rh。2.根据权利要求1所述的方法，其特征在于，对实际复式河道断面形状进行概化的倒等腰梯形断面数量为2个以上；最下面一个梯形断面的下底宽最小取值为0。3.根据权利要求1所述的方法，其特征在于，在步骤S40中，存储各子流域的梯形断面的概化参数的格式具体为：第一行标识整个流域中所有子流域所采用的梯形断面数量最大值，用于在模型应用中确定数组大小；第二行为标题行，其中，第1列为子流域编号，第2列表示当前子流域断面实际拟合梯形断面的个数，第3至5列依次表示该子流域断面对应的第一个梯形断面的上底宽、下底宽、高三个参数，第6至8列表示该子流域断面对应的第二个梯形断面的三个参数，依次类推，直到记录完该子流域断面对应的所有梯形断面的参数；如果没有对...

【专利技术属性】
技术研发人员：刘佳嘉，严子奇，周祖昊，龚家国，桑学锋，牛存稳，贾仰文，王浩，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11