一种分布式河北模型构建方法及其应用技术

本发明专利技术涉及一种分布式河北模型构建方法及其应用，包括下述步骤：步骤1、流域概化；步骤2、以格网为研究单元，利用河北模型，建立格网产流模块；步骤3、以滞后演算法以及马斯京根法建立逐格网分布式汇流模块。本发明专利技术利用成熟的河北模型进行格网的产流计算，再进行逐格网汇流计算，模型充分考虑了下垫面等因素的空间分布不均匀的情况，并可实现与大气模式进行耦合，从而在提高暴雨洪水的预报精度的同时，延长遇见周期，对我国北方半干旱、半湿润地区暴雨洪水预报、流域水循环机理以及水资源规划、水资源统一管理调度等具有重要意义。

【技术实现步骤摘要】

本专利技术涉及一种新型的分布式河北模型的构建方法及其应用，更具体地说是一种有效利用河北模型，以格网为研究单元进行产流分析，再进行逐格网汇流计算的分布式水文模型。
技术介绍
河北省是我国受气候变化和人类活动共同影响的区域之一。近几十年来，气候变暖和水资源的高度开发与利用使该地区的水循环过程严重偏离了其自然演化规律，而水循环过程又深刻影响着该地区的水资源系统、生态环境系统以及社会经济系统。气候变化的影响主要表现为大气环流改变，引起降水时空分布、强度和总量的变化等，人类活动的影响主要表现在对流域下垫面的影响和水资源的取、用、耗、排。河北省境内，人类活动对水循环的影响很大，不仅地表水资源的利用率超过100％，而且地下水资源已经超采超过30年，人类社会取用水资源不仅包括天然水源，而且包括大量的再生水。一方面，地下水位持续下降，土壤含水量较低，河道漏损量较大，常出现“降雨不产流”的情况；另一方面，气候变化等提高了历时短、雨量大的暴雨发生频率，使洪涝灾害风险加大。河北模型是结合河北省的流域特性和地理特点，根据已有的流域产汇流成果和预报经验，在对本地区的产汇流条件和特性进行专门分析研究的基础上，研制出的适合于河北省特点的洪水预报模型，在我国北方半干旱、半湿润地区尤其是海河流域得到广泛应用，预报精度较高、应用效果良好。但目前应用的河北模型属于集总式模型，无法考虑下垫面等因素的空间分布不均匀的情况，模型参数较少，实用性很强，但也使得河北模型在模拟降雨径流过程时存在较大局限性。当前数值大气模式与水文模型的耦合是水文预报的发展方向之一，具有提升暴雨洪水的预报精度和延长遇见期等优点，数值大气模式的降雨输出是以格网为单元的，客观上对构建具有良好模拟精度的分布式的水文模型提出要求，以便将水文模型与大气模式进行耦合。而目前河北模型无法直接与大气模式耦合，不能提供格网化的模拟预报结果。
技术实现思路
本专利技术设计了一种新型的分布式河北模型构建方法及其应用，其解决的技术问题是将集总式的河北模型改成分布式的水文模型，考虑了水文现象和要素的空间分布不均问题，实现与大气模式进行耦合，延长预报周期。为了解决上述存在的技术问题，本专利技术采用了以下方案：一种新型的分布式河北模型构建方法，它包括下述步骤：步骤1、流域概化；步骤2、以格网为研究单元，利用河北模型，建立产流模块；步骤3、以滞后演算法以及马斯京根法建立逐格网分布式汇流模块。进一步，所述步骤1中流域概化包括格网的水流方向概化、子流域概化、流域水系概化与流域边界概化，具体来说：a、格网的水流方向概化是对原始的DEM数据进行填洼处理，根据最陡坡度原则确定每一个格网的水流方向，并对格网进行编号，得到格网演算次序矩阵；b、子流域概化是按照一定的子流域出口断面格网阈值，确定属于子流域的格网，形成各子流域的边界；所述阈值是指格网的演算次序；c、流域水系概化是连接各子流域出口断面格网，形成概化的流域水系；d、流域边界概化是连接各子流域的边界，形成流域的边界。进一步，所述步骤2中以格网为研究单元，利用河北模型，建立产流模块，将天然径流分为地表径流和地下径流两种水源，先用下渗曲线和流域下渗能力分配曲线计算每一个格网的地表径流，再根据时段累积下渗量和流域蓄水容量分配曲线判断格网是否产生地下径流,并计算格网的地下径流量。进一步，所述步骤3中以滞后演算法以及马斯京根法建立逐格网分布式汇流模块，首先针对每一个子流域，通过格网演算次序矩阵，将每一个格网单元的地表径流、地下径流按照滞后演算法逐格网演算至子流域的出口断面格网；然后通过马斯京根法将各子流域的地表、地下径流演算至流域出口断面格网。在进行逐格网演算时，如果格网的土壤含水量未达到田间持水量，则上游网格的出流量首先补充格网的土壤含水量。上述构建方法的一种应用，其特征在于：新型的分布式河北模型对当前流域的径流量进行精细的模拟，并按照模拟结果进行水资源合理配置。上述构建方法的另一种应用，其特征在于：新型的分布式河北模型对未来径流量进行预测，从而为水资源管理、洪水预报提供参考。该新型的分布式河北模型构建方法及其应用具有以下有益效果：(1)本专利技术利用成熟的河北模型进行格网的产流计算，再进行逐格网汇流计算，实现对未来径流量的预测，模型对我国北方半干旱、半湿润地区暴雨洪水预报、流域水循环机理以及水资源规划、水资源统一管理调度等具有重要意义。(2)本专利技术充分利用GIS、DEM技术概化了流域，形成了数字化流域；(3)本专利技术利用河北模型，以单个格网为研究单元建立了格网产流模块，既利用了河北模型的结构和原理，又充分考虑了水文现象和要素的空间分布不均问题；(4)本专利技术大气模式模拟的降雨数据可作为分布式河北模型的输入数据，从而实现耦合，延长预报周期。附图说明图1：本专利技术新型的分布式河北模型计算流程图。具体实施方式下面结合图1，对本专利技术做进一步说明：一种新型的分布式河北模型，它包括下述步骤：步骤1、流域概化；步骤2、以格网为研究单元，利用河北模型，建立产流模块；步骤3、以滞后演算法以及马斯京根法建立逐格网分布式汇流模块。步骤1中流域概化步骤包含：格网的水流方向概化、子流域概化、流域水系概化与流域边界概化。a.格网的水流方向概化。①对原始的数字高程模型数据(DEM数据)进行填洼处理；DEM是“Digital Elevation Model”的缩写，是用一组有序数值阵列形式表示地面高程的一种实体地面模型。DEM数据是指各类DEM产品，例如SRTM数据是由美国太空总署(NASA)和国防部国家测绘局(NIMA)联合测量得到的数据。数字高程模型(DEM)表面往往会存在一些凹陷的区域，即洼地，洼地有可能是DEM对真实地形(如湖泊)的反映，也有可能是由于DEM数据测量误差所致。在进行水流方向计算时，洼地的存在会导致计算产生错误的结果，无法形成连续的汇流网络，因此首先需要对DEM数据进行“填洼处理”。“填洼处理”是指先利用水流方向数据计算出DEM中的洼地区域及其深度，由于与真实地形相比，DEM数据误差产生的洼地深度不会很大，我们可依据这些洼地深度数据设定深度阈值，将低于这个阈值的洼地格网高程增加一定的数值。②根据最陡坡度原则确定每一个格网的水流方向，并对格网进行编号，得到格网演算次序矩阵。这里“坡度”是指中心格网与邻域格网的高程差除以两格网间的距离。根据最陡坡度原则确定每一个格网的水流方向，即当进行坡度计算时，得到中心格网与某一邻域格网的坡度最大，那么中心格网的产流流向该邻域格网。b.子流域概化：按照一定的子流域出口断面格网阈值(阈值是指格网的演算次序)，确定属于子流域的格网，形成各子流域的边界。上述“演算次序”是指模型计算的顺序。本专利技术中演算次序为1的格网属于第一级格网，演算次序为2的格网为第二级格网，依此类推。阈值是指格网的演算次序，即格网所属的级别。若将子流域出口断面格网阈值定义为18，那么就将针对演算次序小于等于18的格网进行所属子流域的归类。演算次序等于18的格网代表一个子流域出口断面，有多少个18级格网就有多少个子流域。再根据水流流向判断其他格网属于哪一个特定的子流域，判定原则是：若水流流向某一特定的18级格网，那么该格网就属于这一特定子流域。c.流域水系概化：连接各子流域出口断面格网，形成本文档来自技高网...

【技术保护点】
一种新型的分布式河北模型构建方法，包括下述步骤：步骤1、流域概化；步骤2、以格网为研究单元，利用河北模型，建立格网产流模块；步骤3、以滞后演算法以及马斯京根法建立逐格网分布式汇流模块。

【技术特征摘要】
1.一种新型的分布式河北模型构建方法，包括下述步骤：步骤1、流域概化；步骤2、以格网为研究单元，利用河北模型，建立格网产流模块；步骤3、以滞后演算法以及马斯京根法建立逐格网分布式汇流模块。2.根据权利要求1所述新型的分布式河北模型构建方法，其特征在于：所述步骤1中流域概化依次进行格网的水流方向概化、子流域概化、流域水系概化以及流域边界概化，具体来说：a、格网的水流方向概化是对原始的DEM数据进行填洼处理，根据最陡坡度原则确定每一个格网的水流方向，并对格网进行编号，得到格网演算次序矩阵；b、子流域概化是按照一定的子流域出口断面格网阈值，确定属于子流域的格网，形成各子流域的边界；所述阈值是指格网的演算次序；c、流域水系概化是连接各子流域出口断面格网，形成概化的流域水系；d、流域边界概化是连接各子流域的边界，形成流域的边界。3.根据权利要求1或2所述新型的分布式河北模型构建方法，其特征在于：所述步骤2中以格网为研究单元，利用河北模型，建立产流模块，将天然径流分为地表径流和地下径流...

【专利技术属性】
技术研发人员：刘佳，李传哲，田济扬，聂汉江，于福亮，穆文彬，邱庆泰，范玲雪，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京;11