CN108304668B - 一种结合水文过程数据和历史先验数据的洪水预测方法 - Google Patents
一种结合水文过程数据和历史先验数据的洪水预测方法 Download PDFInfo
- Publication number
- CN108304668B CN108304668B CN201810140277.2A CN201810140277A CN108304668B CN 108304668 B CN108304668 B CN 108304668B CN 201810140277 A CN201810140277 A CN 201810140277A CN 108304668 B CN108304668 B CN 108304668B
- Authority
- CN
- China
- Prior art keywords
- data
- hydrological
- historical
- knowledge
- variables
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 230000008569 process Effects 0.000 title claims abstract description 60
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 33
- 230000001364 causal effect Effects 0.000 claims abstract description 21
- 238000012549 training Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 238000004140 cleaning Methods 0.000 claims abstract description 4
- 230000001502 supplementing effect Effects 0.000 claims abstract description 4
- 239000002689 soil Substances 0.000 claims description 28
- 238000009826 distribution Methods 0.000 claims description 10
- 238000007619 statistical method Methods 0.000 claims description 8
- 230000008020 evaporation Effects 0.000 claims description 4
- 238000001704 evaporation Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 238000003064 k means clustering Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000007711 solidification Methods 0.000 description 5
- 230000008023 solidification Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000006424 Flood reaction Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Alarm Systems (AREA)
Abstract
本发明公开了一种结合水文过程数据和历史先验数据的洪水预测方法,利用加权聚类算法对影响降雨信息进行聚类分析,近似含水量信息;对数据进行清理,补充缺失的数据,修正错误数据;根据流域特征选取通用水文模型,获取水文过程中的变量之间的因果推理关系;对获取的一组组推理关系,采用有向图的形式进行表示,其中节点表示物理变量,边的方向表示因果推理关系,并以此拓扑结构构建贝叶斯网络;利用历史数据,结合通用水文模型,模拟产生中间结果数据,共同组成样本集;根据样本集获取历史先验知识,计算条件概率表,以实现对网络的训练,完成水文过程知识和历史先验知识的综合建模;输入测试数据,实现洪水预测。
Description
技术领域
本发明涉及一种结合水文过程数据和历史先验数据的洪水预测方法,属于洪水预报技术领域。
背景技术
山洪是山丘区中小河流由降水引发的突发性、暴涨暴落的洪水,而我国中小河流众多,各个流域具有地理地貌、气候的独特性,缺少统一便捷的建模手段。近年来,由于极端天气时间频多,地区持续性强降雨、突发性暴雨时常发生,而山丘地区下垫面呈现高陡坡、河流源短流急,在暴雨天气下极易发生山洪等山地灾害,对人民生命财产造成威胁。因此基于中小河流的洪水预报得到国家的高度重视。
在防洪领域中,常用的模型主要分基于物理(过程)驱动和基于数据驱动。物理驱动模型是基于水文概念模型的,该方法分析造成洪水的水文过程,建立适合某一流域的通用模型,模型参数需要依赖经验和人工交互不断迭代才能完成率定。模型预报的精度取决于建模者的知识和经验以及数据资料的完备情况。这类模型内包含的参数具备物理意义,具有良好的解释性,但中小河流复杂的环境因素加大了参数率定的难度。近年来,线性回归、神经网络、支持向量机等数据驱动模型被运用在洪水预报中,这类模型将水文过程视作黑箱子,不考虑系统内部的物理机制,通过建立输入输出样本间的映射关系以实现建模。由于无法获取精确的水文要素(降雨和流量),也无法获取流域内水循环要素的可靠时空分布,造成智能模型在预报性能上有一定的局限性。水文预报中,各个因素存在不确定性,因素之间相互关联,这使得物理驱动模型和数据驱动模型均有各自的不足。
发明内容
发明目的:针对现有技术中存在的问题,本发明提供一种结合水文过程数据和历史先验数据的洪水预测方法。该方法选定适合的通用水文模型,分解水文过程构建贝叶斯网络,通过通用水文模型产生中间变量结果以丰富样本集,并用统计学方法获得变量之间的先验知识,完成水文过程的知识固化。本方法选用通用水文模型,避免了精确参数的率定,但同时利用贝叶斯网络的特点,能对不确定知识有良好的表达。
技术方案:一种结合水文过程数据和历史先验数据的洪水预测方法,包括如下步骤:
第一步:收集水文要素(历史降雨量数据、历史流量数据、历史蒸发量数据等),根据降雨数据加权聚类近似推算土壤含水量,构成初始数据集。对数据进行清洗,对缺失数据进行补充;
第二步:根据流域特征选取适合该流域预报的水文模型,分析水文过程中的变量及变量之间的因果关系,收集水文过程知识,对水文过程知识进行建模,采用有向图的形式来表示水文过程,其中节点表示物理变量,边的方向表示因果推理关系,并以此拓扑结构构建贝叶斯网络;
第三步:利用历史数据,运行通用水文模型,模拟产生中间结果,历史水文要素数据和中间结果数据共同构成数据样本集,对样本集进行数理统计分析,分别计算每个变量的概率分布和变量之间的条件概率分布,以此获取历史先验知识,计算贝叶斯网络的每个节点的条件概率表,以实现对网络的训练,完成水文过程知识和历史先验知识的综合建模;
第四步:完成知识固化的贝叶斯网络可以替代水文模型,再输入新的水文要素,可以完成水文要素的预测。
第一步中加权聚类近似产生土壤含水量的步骤如下:
步骤1.1:收集整理历史降雨量;
步骤1.2:选取降雨特征,对土壤含水量计算点过去15天的降雨数据进行整理分析,抽取出多个时段的累计降雨量作为特征;
步骤1.4:对相对湿润度A进行排序,对聚类质心进行定级,评定土壤的含水百分比;
步骤1.5:在进行新一轮土壤含水量预测时,根据已有的降水特征,计算与各个聚类质心的距离,选取较近的质心评级近似该土壤含水量。
第二步中获取水文过程知识的步骤如下:
步骤2.1:综合分析目标流域的水文特性、气候、土壤特性等,选定适合的通用模型;
步骤2.2:在选定通用模型之后,分析水文过程中的变量及变量之间的因果关系,收集水文过程知识(洪水的物理过程中各个物理量之间的因果关系),具体的变量间因果关系如下所示:
(Et,Pt,WMt)→(Rt+1,WMt+1)
(St,Rt+1)→(St+1,RSt+1,RIt+1,RGt+1)
(TRSt+1,TRIt+1,TRGt+1)→(TRt+1)
(Qt-1,Qt,TRt+1)→(Qt+1,…,Qt+h)
其中下标t表示时刻,E表示蒸发量,P表示降雨量,WM表示土壤含水量(张力水),R表示产流量,S表示自由水含量,RS表示地表径流,RI表示壤中流,RG表示地下径流,TRS表示全流域的地表径流出流量,TRI表示全流域的壤中流出流量,TRG表示全流域的地下径流出流量,TR表示全流域出流,Q表示流量。
对水文过程知识进行建模,采用有向图的形式来表示水文过程,其中节点表示物理变量,边的方向表示因果推理关系;
步骤2.3:根据已有的网络拓扑图,构建贝叶斯网络,网络为有向图,有向边起点变量作为有向边终点变量的条件,节点保存变量的条件概率表。
第三步中利用贝叶斯网络以实现水文过程知识和历史先验知识的综合建模步骤如下:
步骤3.1:利用历史数据和通用的水文模型产生中间物理量的取值,需要产生的变量有R,RS,RI,RG,TRS,TRG,TRI,TRG,TR等,历史水文要素数据和中间结果数据共同构成数据集,对数据集进行离散化处理和规范化处理;R表示产流量,RS表示地表径流,RI表示壤中流,RG表示地下径流,TRS表示全流域的地表径流出流量,TRI表示全流域的壤中流出流量,TRG表示全流域的地下径流出流量,TR表示全流域出流。
步骤3.2:随机对数据集进行分组,其中2/3数据集构成训练集,1/3数据集构成测试集;
步骤3.3:针对网络中任意节点Y,其前置节点X1,X2,…,Xk,通过统计分析求解条件概率分布P(Y|X1,X2,…,Xk),完成对节点参数的训练;为增加网络的自适应性和泛化能力,可以根据一定的规则对数据集进行扩充,如假定前置节点Xi符合某种分布,根据分布特性产生数据对数据集进行扩充。
第四步中利用综合模型进行预报的步骤如下:
步骤4.1:输入t时刻的水文要素X1,X2,…Xn的具体取值X1=xt,1,X2=xt,2,…Xn=xt,n;
步骤4.2:计算argmaxyP(Y=y|X1=xt,1,X2=xt,2,…,Xn=xt,n),此时y作为网络的输出,即径流取值;
步骤4.3:步骤4.2中同时也可以获得中间变量的节点,作为t+1时刻的输入,继续进行推理预测,完成结果的连续预测。
有益效果:与现有技术相比,本发明提供的结合水文过程数据和历史先验数据的洪水预测方法,对影响降雨信息聚类获得难以实际测量的土壤含水量信息,再通过分析通用水文模型,提取物理量之间的因果关系,建立贝叶斯网络,通过对历史样本的统计分析,实现对网络的训练,以实现水文过程知识的固化,通过贝叶斯网络对洪水进行预报。
附图说明
图1为本发明的方法流程图;
图2为是实施例中综合模型和通用水文模型对比图,(a)、(b)、(c)、(d)分别是预见期1~4小时的结果展示;
图3为模型拓扑结构。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,结合水文过程数据和历史先验数据的洪水预测方法,主要包括以下几步:
一是利用加权聚类算法对影响降雨信息量进行聚类分析,根据以往15天的降雨情况来表征土地的干旱程度,以此来近似含水量信息;二是对数据进行清理,补充缺失的数据,修正错误数据,采用近邻数据的均值来平滑缺失值或修正错误值;三是根据流域特征选取适合的通用水文模型,分析水文过程中的变量及变量之间的因果关系,收集水文过程知识;四是对水文过程知识进行建模,采用有向图的形式来表示水文过程,其中节点表示物理变量,边的方向表示因果推理关系,并以此拓扑结构构建贝叶斯网络;五是利用历史数据,结合通用水文模型,模拟产生中间结果数据,加入到样本集;六是对样本集进行统计分析,获取历史先验知识,计算条件概率表,以实现对网络的训练,完成水文过程知识和历史先验知识的综合建模;七是输入测试数据,实现洪水预测。
由于中小河流所处地理环境复杂,土壤类型多样,以往测量土壤含水量方法效果不佳,再加上部分地区并没有土壤含水量的资料信息,因此本文提出一种用影响降雨信息来近似推算土壤含水情况的方法。以下是加权聚类计算土壤含水量的算法步骤:
步骤1.1:收集整理历史降雨量;
步骤1.2:选取降雨特征,对土壤含水量计算点过去15天的降雨数据进行整理分析,抽取出多个时段的累计降雨量作为特征;
步骤1.4:对相对湿润度A进行排序并定级,以此实现对聚类质心的评定(评定土壤的含水百分比);
步骤1.5:在进行土壤含水量推算时,根据已有的降水特征,计算与各个聚类质心的距离,选取较近的质心评级近似该土壤含水量。
物理驱动模型在传统水文预报中有很多发展,因此也衍生了很多种适合不同流域特性的水文模型,水文模型中蕴含了很多水文过程的知识,对水文模型分析可以提炼水文过程知识,并采用贝叶斯网络进行表达,其具体算法步骤如下:
步骤2.1:综合分析目标流域的水文特性、气候、土壤特性等,选定适合的通用水文模型;
步骤2.2:在选定通用模型之后,分析水文过程中的变量及变量之间的因果关系,收集水文过程知识,具体如下所示:
(Et,Pt,WMt)→(Rt+1,WMt+1)
(St,Rt+1)→(St+1,RSt+1,RIt+1,RGt+1)
(TRSt+1,TRIt+1,TRGt+1)→(TRt+1)
(Qt-1,Qt,TRt+1)→(Qt+1,…,Qt+h)
其中E表示蒸发量,P表示降雨量,WM表示土壤含水量(张力水),R表示产流量,S表示自由水含量,RS表示地表径流,RI表示壤中流,RG表示地下径流,TRS表示全流域的地表径流,TRI表示全流域的壤中流出流量,TRG表示全流域的地下径流出流量,TR表示全流域出流,Q表示流量。
对水文过程知识进行建模,采用有向图的形式来表示水文过程,其中节点表示物理变量,边的方向表示因果推理关系;
步骤2.3:根据已有的网络拓扑图,构建贝叶斯网络,网络为有向图,有向边起点变量作为有向边终点变量的条件,节点保存变量的条件概率表。
完成对水文过程知识的获取之后,需要采用历史先验知识对该网络进行训练。由于网络中部分中间节点并没有数据,为解决这个问题,利用通用水文模型,基于历史水文要素数据,产生中间物理量结果数据,加入数据集,再基于数理统计分析,获取历史先验知识,对网络进行训练,其具体的算法步骤如下:
步骤3.1:利用历史数据和通用的水文模型产生中间物理量的取值,需要产生的变量有R,RS,RI,RG,TRS,TRG,TRI,TRG,TR等,历史水文要素数据和中间结果数据共同构成数据集,对数据集进行离散化处理和规范化处理;
步骤3.2:随机对数据集进行分组,其中2/3数据集构成训练集,1/3数据集构成测试集;
步骤3.3:针对网络中任意节点Y,其前置节点X1,X2,…,Xk,通过统计分析求解条件概率分布P(Y|X1,X2,…,Xk),完成对节点参数的训练;为增加网络的自适应性和泛化能力,可以根据一定的规则对数据集进行扩充,如假定前置节点Xi符合某种分布,根据分布特性产生数据对数据集进行扩充。
在综合水文过程知识和历史先验知识建模后,可以用该模型进行水文预报,其具体的算法如下:
步骤4.1:输入t时刻的水文要素X1,X2,…Xn的具体取值X1=xt,1,X2=xt,2,…Xn=xt,n;
步骤4.2:计算argmaxyP(Y=y|X1=xt,1,X2=xt,2,…,Xn=xt,n),此时y作为网络的输出,即径流取值;
步骤4.3:步骤4.2中同时也可以获得中间变量的节点,作为t+1时刻的输入,继续进行推理预测,完成结果的连续预测。
为了验证本发明的效果,选取长湖汛期数据作为研究对象,昌化流域位于浙江省分水江流域上游,地势西北高东南高,属浙西山丘区,典型的中小河流水系。选取1998-2010年汛期场次洪水数据,数据时间间隔为1小时,其中1998年-2009年攻6790个数据作为训练样本,2010共671个数据作为测试样本。首先通过对雨量信息进行聚类,对土壤含水量进行近似估算,然后选定合适昌化流域的通用水文模型进行分析,建立贝叶斯网络模型,对样本数据进行数理统计分析,计算条件概率表,完成对网络的训练,最终完成水文过程的专家知识的固化,并投入预报。如图2所示,本文通过对比固化模型和通用水文模型,其确定性在不同预见期的对比结果如下:
预见期 | 固化模型 | 通用模型 |
1h | 0.939455 | 0.71 |
2h | 0.981491 | 0.67 |
3h | 0.801938 | 0.67 |
4h | 0.605704 | 0.66 |
实验结果表明:固化模型能对水文过程的专家知识进行良好的固化,并且具有良好的预测能力,这是因为固化模型不仅基于传统的通用水文模型,还借鉴了数据模型的特点,把严谨的数理推导转化成因果概率推导,弱化了计算性,降低了传统水文模型的应用难度。
Claims (8)
1.一种结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,包括如下步骤:
步骤1,收集水文要素,根据降雨数据加权聚类近似推算土壤含水量,构成初始数据集,对数据进行清洗,对缺失数据进行补充;
步骤2,根据流域特征选取通用水文模型,通过分析通用水文模型,获取水文过程中的变量之间的因果推理关系,对每组因果推理关系进行建模,采用有向图的形式进行表示,其中节点表示物理变量,边的方向表示因果推理关系,并以此拓扑结构构建贝叶斯网络;
步骤3,利用所述初始数据集中的历史数据,运行通用水文模型,模拟产生中间结果,历史水文要素数据和中间结果数据共同构成数据样本集,对样本集进行统计分析,获取历史先验知识,计算贝叶斯网络的每个节点的条件概率表,以实现对网络的训练,完成水文过程知识和历史先验知识的综合建模;
步骤4,对于综合建模的模型,输入新的水文要素,完成水文要素的预测。
2.如权利要求1所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,步骤1中加权聚类近似产生土壤含水量的步骤如下:
步骤1.1:收集整理历史降雨量;
步骤1.2:选取降雨特征,抽取出多个时段的累计降雨量作为特征;
步骤1.3:对多个特征初步赋权,再进行k-means聚类,获得多个聚类质心,计算相对湿润度A;
步骤1.4:对相对湿润度A进行排序,对聚类质心进行定级,评定土壤的含水百分比;
步骤1.5:在进行土壤含水量预测时,根据已有的降水特征,计算与各个聚类质心的距离,选取较近的质心评级近似该土壤含水量。
4.如权利要求1所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,步骤2中获取水文过程知识的步骤如下:
步骤2.1:选定通用模型;
步骤2.2:在选定通用模型之后,分析水文过程中的变量及变量之间的因果关系,收集水文过程知识,对水文过程知识进行建模,采用有向图的形式来表示水文过程,其中节点表示物理变量,边的方向表示因果推理关系;
步骤2.3:根据已有的网络拓扑图,构建贝叶斯网络,网络为有向图,有向边起点变量作为有向边终点变量的条件,节点保存变量的条件概率表。
5.如权利要求4所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,收集水文过程知识,具体如下所示:
(Et,Pt,WMt)→(Rt+1,WMt+1)
(St,Rt+1)→(St+1,RSt+1,RIt+1,RGt+1)
(TRSt+1,TRIt+1,TRGt+1)→(TRt+1)
(Qt-1,Qt,TRt+1)→(Qt+1,…,Qt+h)
其中,下标t表示时刻,h表示最大观测步长,E表示蒸发量,P表示降雨量,WM表示土壤含水量,R表示产流量,S表示自由水含量,RS表示地表径流,RI表示壤中流,RG表示地下径流,TRS表示全流域的地表径流出流量,TRI表示全流域的壤中流出流量,TRG表示全流域的地下径流出流量,TR表示全流域出流,Q表示流量。
6.如权利要求1所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,步骤3中利用贝叶斯网络以实现水文过程知识和历史先验知识的综合建模步骤如下:
步骤3.1:利用历史数据和通用的水文模型产生中间物理量的取值,需要产生的变量有R,RS,RI,RG,TRS,TRG,TRI,TR,历史水文要素数据和中间结果数据共同构成数据集,对数据集进行离散化处理和规范化处理;R表示产流量,RS表示地表径流,RI表示壤中流,RG表示地下径流,TRS表示全流域的地表径流出流量,TRI表示全流域的壤中流出流量,TRG表示全流域的地下径流出流量,TR表示全流域出流;
步骤3.2:随机对数据集进行分组,一组为训练集,一组为测试集;
步骤3.3:针对网络中任意节点Y,其前置节点X1,X2,…,Xk,通过统计分析求解条件概率分布P(Y|X1,X2,…,Xk),完成对节点参数的训练。
7.如权利要求6所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,步骤3.3中为增加网络的自适应性和泛化能力,根据前置节点Xi符合的规则对数据集进行扩充。
8.如权利要求1所述的结合水文过程数据和历史先验数据的洪水预测方法,其特征在于,步骤4中利用综合模型进行预报的步骤如下:
步骤4.1:输入t时刻的水文要素X1,t,X2,t,…Xn,t的具体取值X1,t=xt,1,X2,t=xt,2,…Xn,t=xt,n;其中,xt,i表示第i个观测点在时刻t具体的取值;
步骤4.2:计算argmaxyP(Y=y|X1,t=xt,1,X2,t=xt,2,…,Xn,t=xt,n),此时y作为网络的输出,即径流取值;
步骤4.3:步骤4.2中同时也可以获得中间变量的节点,作为t+1时刻的输入,继续进行推理预测,完成结果的连续预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810140277.2A CN108304668B (zh) | 2018-02-11 | 2018-02-11 | 一种结合水文过程数据和历史先验数据的洪水预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810140277.2A CN108304668B (zh) | 2018-02-11 | 2018-02-11 | 一种结合水文过程数据和历史先验数据的洪水预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108304668A CN108304668A (zh) | 2018-07-20 |
CN108304668B true CN108304668B (zh) | 2021-07-09 |
Family
ID=62865277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810140277.2A Active CN108304668B (zh) | 2018-02-11 | 2018-02-11 | 一种结合水文过程数据和历史先验数据的洪水预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108304668B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109491991B (zh) * | 2018-11-08 | 2022-03-01 | 四川长虹电器股份有限公司 | 一种无监督的数据自动清洗方法 |
CN109767264A (zh) * | 2018-12-20 | 2019-05-17 | 深圳壹账通智能科技有限公司 | 产品数据推送方法、装置、计算机设备和存储介质 |
CN110276477B (zh) * | 2019-05-29 | 2022-08-16 | 河海大学 | 一种基于分层贝叶斯网络和增量学习的洪水预报方法 |
CN110222823B (zh) * | 2019-05-31 | 2022-12-30 | 甘肃省祁连山水源涵养林研究院 | 一种水文流量波动情势识别方法及其系统 |
CN110428586B (zh) * | 2019-07-08 | 2021-01-05 | 大连理工大学 | 基于前期降雨和上下游拓扑关系的农村基层洪涝预警方法 |
US11386143B2 (en) | 2019-08-30 | 2022-07-12 | International Business Machines Corporation | Searching for analogue subsurface structures based on topological knowledge representation (TKR) |
CN111126699B (zh) * | 2019-12-25 | 2022-06-17 | 华中科技大学 | 一种考虑降雨量等级的水文预报方法及系统 |
CN111325400B (zh) * | 2020-02-20 | 2023-06-02 | 内蒙古自治区水利水电勘测设计院 | 一种高海拔长距离的输水定位方法及其定位系统 |
CN113159378B (zh) * | 2021-03-15 | 2022-04-05 | 中国科学院地理科学与资源研究所 | 无测站山区径流与遥感降水相结合的降水估算方法 |
CN113222368B (zh) * | 2021-04-28 | 2022-06-03 | 中国科学院城市环境研究所 | 一种基于雨水花园监测数据的雨洪预警方法 |
CN113377750B (zh) * | 2021-05-19 | 2022-03-29 | 山东大学 | 水文数据清洗方法及系统 |
CN113642794B (zh) * | 2021-08-16 | 2023-11-24 | 浙江大学 | 一种结合降雨和土壤水观测的山洪预报方法 |
CN113837475B (zh) * | 2021-09-27 | 2024-04-05 | 中水珠江规划勘测设计有限公司 | 有向图深度神经网络径流概率预报方法、系统、设备及终端 |
CN115081341B (zh) * | 2022-07-25 | 2022-11-11 | 江西武大扬帆科技有限公司 | 一种流域洪水仿真预警方法及系统 |
CN115115148B (zh) * | 2022-08-30 | 2022-12-30 | 武汉大学 | 基于过程-数据协同驱动的长期径流预报方法及系统 |
CN115146746B (zh) * | 2022-09-06 | 2022-11-22 | 武汉商学院 | 基于多源数据融合的区域滑坡灾害识别方法 |
CN116542124B (zh) * | 2023-03-13 | 2024-04-09 | 广东省科学院广州地理研究所 | 一种分布式水文模型的辅助建模方法 |
CN117828312B (zh) * | 2024-03-05 | 2024-05-24 | 中国科学院地理科学与资源研究所 | 一种流域水文环境的管理方法及相关设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298841A (zh) * | 2013-07-16 | 2015-01-21 | 杭州贵仁科技有限公司 | 一种基于历史数据的洪水预报方法和系统 |
KR101690181B1 (ko) * | 2015-07-16 | 2017-01-09 | 한국외국어대학교 연구산학협력단 | 돌발홍수지수 예측을 이용한 돌발홍수 모니터링 방법. |
CN106650767A (zh) * | 2016-09-20 | 2017-05-10 | 河海大学 | 基于聚类分析和实时校正的洪水预报方法 |
CN106845771A (zh) * | 2016-12-16 | 2017-06-13 | 中国水利水电科学研究院 | 一种基于前期雨量优选参数的洪水预报方法 |
CN106991278A (zh) * | 2017-03-21 | 2017-07-28 | 武汉大学 | 一种集合降水预报与实时洪水概率预报的耦合方法 |
-
2018
- 2018-02-11 CN CN201810140277.2A patent/CN108304668B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298841A (zh) * | 2013-07-16 | 2015-01-21 | 杭州贵仁科技有限公司 | 一种基于历史数据的洪水预报方法和系统 |
KR101690181B1 (ko) * | 2015-07-16 | 2017-01-09 | 한국외국어대학교 연구산학협력단 | 돌발홍수지수 예측을 이용한 돌발홍수 모니터링 방법. |
CN106650767A (zh) * | 2016-09-20 | 2017-05-10 | 河海大学 | 基于聚类分析和实时校正的洪水预报方法 |
CN106845771A (zh) * | 2016-12-16 | 2017-06-13 | 中国水利水电科学研究院 | 一种基于前期雨量优选参数的洪水预报方法 |
CN106991278A (zh) * | 2017-03-21 | 2017-07-28 | 武汉大学 | 一种集合降水预报与实时洪水概率预报的耦合方法 |
Non-Patent Citations (4)
Title |
---|
Bayesian flood forecasting methods: A review;Shasha Han 等;《Journal of Hydrology》;20170831;第340-351页 * |
Multimodel Combination Techniques for Analysis of Hydrological Simulations: Application to Distributed Model Intercomparison Project Results;Newsha K.Ajami 等;《Journal of Hydrometeorology》;20060831;第755-768页 * |
动态贝叶斯网络在水文预报中的应用;李维乾;《中国优秀硕士学位论文全文数据库工程科技II辑》;20111215;第C037-7卷(第S1期);第1-71页 * |
基于贝叶斯统计的水文模型不确定性研究;李明亮;《中国博士学位论文全文数据库基础科学辑》;20130616;第A012-4卷(第07期);第1-151页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108304668A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304668B (zh) | 一种结合水文过程数据和历史先验数据的洪水预测方法 | |
CN106650767B (zh) | 基于聚类分析和实时校正的洪水预报方法 | |
CN101480143B (zh) | 一种预测灌区作物单产量的方法 | |
CN113705877B (zh) | 基于深度学习模型的实时月径流预报方法 | |
CN117035201B (zh) | 平原河网水工程集群多目标调度规则制定方法及系统 | |
CN114254561A (zh) | 一种内涝预测方法、系统及存储介质 | |
CN111159149A (zh) | 一种基于三维卷积神经网络的河流流量预测方法 | |
CN113902580B (zh) | 一种基于随机森林模型的历史耕地分布重建方法 | |
CN109840587A (zh) | 基于深度学习的水库入库流量预测方法 | |
Sekertekin et al. | Short-term air temperature prediction by adaptive neuro-fuzzy inference system (ANFIS) and long short-term memory (LSTM) network | |
CN111027764B (zh) | 一种基于机器学习的适用于径流资料缺乏流域的洪水预报方法 | |
CN114529097B (zh) | 多尺度农作物物候期遥感降维预测方法 | |
CN113887635B (zh) | 一种流域相似性分类方法及分类装置 | |
CN110276477B (zh) | 一种基于分层贝叶斯网络和增量学习的洪水预报方法 | |
CN114372631A (zh) | 一种基于小样本学习和lstm的缺资料地区径流预测方法 | |
CN113836808A (zh) | 一种基于重污染特征约束的pm2.5深度学习预测方法 | |
CN116205136A (zh) | 基于径流滞后信息的大尺度流域深度学习洪水预报方法 | |
CN115759445A (zh) | 一种基于机器学习和云模型的分类洪水随机预报方法 | |
Nivesh et al. | Prediction of river discharge of Kesinga sub-catchment of Mahanadi basin using machine learning approaches | |
Zanella et al. | Internet of things for hydrology: Potential and challenges | |
El Ibrahimi et al. | Application of Neural Modeling and the SPI Index for the Prediction of Weather Drought in the Saïss Plain (Northern Morocco). | |
Nasseri et al. | A comparison between direct and indirect frameworks to evaluate impacts of climate change on streamflows: case study of Karkheh River basin in Iran | |
CN117010274A (zh) | 一种基于集成增量学习的地下水有害元素智能预警方法 | |
WO2023245399A1 (zh) | 基于土地系统和气候变化耦合的水稻生产潜力模拟方法 | |
Velasco et al. | Performance analysis of multilayer perceptron neural network models in week-ahead rainfall forecasting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |