CN112861360B - 基于系统响应理论的马斯京根流量演算误差修正方法 - Google Patents
基于系统响应理论的马斯京根流量演算误差修正方法 Download PDFInfo
- Publication number
- CN112861360B CN112861360B CN202110188437.2A CN202110188437A CN112861360B CN 112861360 B CN112861360 B CN 112861360B CN 202110188437 A CN202110188437 A CN 202110188437A CN 112861360 B CN112861360 B CN 112861360B
- Authority
- CN
- China
- Prior art keywords
- river
- flow
- calculation
- section
- reach
- 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
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
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
- G01F25/10—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
Abstract
本发明公开了基于系统响应理论的马斯京根流量演算误差修正方法,主要包括:将研究流域的长河段划分多个子河段,构建所有子河段的马斯京根汇流演算方程;考虑各个子河段的出流变化是由河段区间的入流变化引起的系统响应,建立所有子河段的响应系统;收集历史预报和实测资料,计算各河段断面的流量预报误差;估计所有区间流量预报误差,反演修正后的区间流量;重新进行马斯京根汇流演算,得到未来时刻修正后的河道断面预报流量。本发明解决了长河段流域的河道洪水预报精度提升的难题,具有物理意义明确、计算效率高、修正结果客观合理等优点,适合在实际中推广和应用。
Description
技术领域
本发明属于水文预报及实时校正技术领域,特别涉及一种基于系统响应理论的马斯京根流量演算误差修正方法。
背景技术
河道洪水预报是一项重要的防洪减灾非工程措施,及时准确的河道洪水预报对河流系统的安全保障和水资源调度等意义重大。一般采用马斯京根法对河道汇流进行演算,对于长河段流域,需划分为多个子河段,再利用考虑区间入流的马斯京根矩阵方程进行连续演算。但是,演算过程中会存在各种来源的误差,包括断面流量预报误差、区间入流预报误差和汇流模型本身的误差,这些误差导致了流域出口断面的洪水预报值误差。当预报误差超过标准时,会给防洪决策带来风险,影响防洪效益。因此,需要对马斯京根流量演算过程中的诸多误差进行修正,以提高河道洪水预报精度,满足实际作业预报需求。
目前,河道洪水预报误差修正方法主要存在几方面不足:(1)由于人类认识的有限性,无法用水文模型或预报方案绝对准确地刻画流域的下垫面机制,难以从根本上降低预报误差,提升预报精度;(2)现有的误差修正方法,一般简化复杂汇流系统为单河段汇流系统进行预报误差修正,未考虑有区间入流的长河段连续汇流演算过程,忽略来源于上游断面和区间的误差,仅关注流域出口断面的误差。
发明内容
为了解决现有技术中存在的不足,本发明提供了基于系统响应理论的马斯京根流量演算误差修正方法,具有物理意义明确、计算效率高、修正结果客观合理等优点,可进一步提升长河段流域的河道洪水预报精度。
为解决上述问题,本发明具体采用以下技术方案:
基于系统响应理论的马斯京根流量演算误差修正方法,包括以下步骤:
步骤1,将一个长河段划分为n个子河段(图2),对于每个子河段都应用马斯京根法进行河道汇流演算,得到n个河段的马斯京根方程下断面预报出流的计算式;具体步骤为:
步骤1.1,对于每一个子河段:采用马斯京根汇流演算方程表达第i个子河段在t时刻的汇流过程:
其中,
式中,在t时刻,第i(=1,…,n)个河段的出流为n>1,第i-1个河段的出流为就是第i个河段的入流,第i个河段的支流汇入为参数组成为ai和bi,马斯京根法演算参数为和在t-Δt时刻,第i个河段的出流为入流为支流汇入为参数组成为ci和di,马斯京根法演算参数为和
整理后得到第i个河段的马斯京根方程预报出流的计算式为:
步骤1.2,同理步骤1.1,得到n个子河段的马斯京根方程下断面预报出流的计算方程组:
式中,在t时刻,第1,2,…,n个河段的出流为入流为支流汇入为参数组成为a1,a2,…,an,b1,b2,…,bn,在t-Δt时刻,第1,2,…,n个河段的出流为入流为支流汇入为参数组成为c1,c2,…,cn,d1,d2,…,dn。
步骤2,考虑河段出流的变化是由于区间入流的变化所引起的系统响应,可建立t时刻的n个河段和n个区间的响应系统;具体步骤为:
步骤2.1,对于每一级河段:认为t时刻第i(=1,…,n)个河段下断面出流的马斯京根法预报结果的误差全部来自于区间流量的预报误差,即考虑河段出流的变化是由于区间入流的变化所引起的系统响应:
上式的微分表达为:
步骤2.2,同理步骤2.1,得到t时刻的所有n级河段的响应系统:
上式写成矩阵形式为:
EQ=Η×Eq (2)
其中,
具体包括以下步骤:
本发明未提及的技术均参照现有技术。
本发明的有益效果:本发明提供的基于系统响应理论的马斯京根流量演算误差修正方法,能够反映多河段多区间入流的河段连续汇流演算过程,考虑了来源于上游断面和区间的误差,解决了长河段河道的洪水预报误差修正难题,进一步提升了河道洪水预报精度,且计算效率高、易于实现,具有较强的工程意义。
附图说明
图1是本发明基于系统响应理论的马斯京根流量演算误差修正方法的流程图;
图2是本发明长河段汇流系统的示意图;
图3是本发明大渡河丹巴站以上流域河段划分的概化图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,本发明提供的基于系统响应理论的马斯京根流量演算误差修正方法,包括以下步骤:
S1、如图2所示,对大渡河上游的一个流域,在其干流选取若干水文站,并将该长河段划分为多个子河段,对每个子河段都采用马斯京根法进行河道汇流演算,得到每个河段的马斯京根方程下断面预报出流的计算式,具体包括以下步骤:
1-1、在大渡河丹巴站以上流域的干流选取4个代表性水文站(日部,足木足,大金,丹巴),据此将该1个长河段划分为3个子河段(河段1:日部—足木足,河段2:足木足—大金,河段3:大金—丹巴),结果如图3所示;
1-2、对这3个子河段,都采用马斯京根汇流演算方程表达在t时刻的河道汇流过程:
其中,
式中,在t时刻,第1,2,3个河段的出流分别为入流分别为支流汇入分别为参数组成分别为a1,a2,a3,b1,b2,b3,马斯京根法演算参数分别为在t-Δt时刻,第1,2,3个河段的出流分别为入流分别为支流汇入分别为参数组成分别为c1,c2,c3,d1,d2,d3,马斯京根法演算参数为
1-3、整理后得到3个子河段的马斯京根方程下断面预报出流的计算式:
S2、考虑河段出流的变化是由于区间入流的变化所引起的系统响应,可建立t时刻的3个子河段和3个区间的响应系统:
上式写成各河段断面流量误差EQ与各区间入流误差Eq的响应系统的矩阵形式为:
EQ=Η×Eq
其中,
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (3)
1.基于系统响应理论的马斯京根流量演算误差修正方法,其特征在于,包括以下步骤:
步骤1,将一个长河段划分为n个子河段,n>1,对于每个子河段都应用马斯京根法进行河道汇流演算,得到n个河段的马斯京根方程下断面预报出流的计算式如式(1):
式中,为第t时刻第i个河段的出流,i=1,…,n,下同;为第t时刻第i个河段的入流;为第t时刻第i个河段的支流汇入;ai和bi为第t时刻的参数;为第t-Δt时刻第i个河段的出流;为第t-Δt时刻第i个河段的入流;为第t-Δt时刻第i个河段的支流汇入;ci和di为第t-Δt时刻的参数;
步骤2,认为河段出流的变化是由于区间入流的变化所引起的系统响应,基于步骤1的断面出流的计算式建立t时刻的n个河段和n个区间的响应系统如式(2)所示:
EQ=H×Eq (2)
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188437.2A CN112861360B (zh) | 2021-02-19 | 2021-02-19 | 基于系统响应理论的马斯京根流量演算误差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188437.2A CN112861360B (zh) | 2021-02-19 | 2021-02-19 | 基于系统响应理论的马斯京根流量演算误差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112861360A CN112861360A (zh) | 2021-05-28 |
CN112861360B true CN112861360B (zh) | 2021-10-26 |
Family
ID=75989684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110188437.2A Active CN112861360B (zh) | 2021-02-19 | 2021-02-19 | 基于系统响应理论的马斯京根流量演算误差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112861360B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115034160B (zh) * | 2022-06-28 | 2023-04-21 | 河海大学 | 一种基于马斯京根法参数转换为等效河道的计算方法 |
CN115017727B (zh) * | 2022-06-28 | 2023-04-21 | 河海大学 | 一种基于马斯京根法的汇污模拟方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
CN106803131A (zh) * | 2016-12-23 | 2017-06-06 | 华中科技大学 | 一种基于离散广义Nash汇流模型的河道洪水预报方法 |
US10078155B2 (en) * | 2016-06-24 | 2018-09-18 | Climacell Inc. | Real-time precipitation forecasting system |
CN111125969A (zh) * | 2019-12-25 | 2020-05-08 | 华中科技大学 | 一种跨水库流域河道径流演算方法和系统 |
CN112199901A (zh) * | 2020-09-14 | 2021-01-08 | 中国人民解放军火箭军研究院核技术研究所 | 一种无径流资料山区小流域山洪设计暴雨洪水计算方法 |
-
2021
- 2021-02-19 CN CN202110188437.2A patent/CN112861360B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
US10078155B2 (en) * | 2016-06-24 | 2018-09-18 | Climacell Inc. | Real-time precipitation forecasting system |
CN106803131A (zh) * | 2016-12-23 | 2017-06-06 | 华中科技大学 | 一种基于离散广义Nash汇流模型的河道洪水预报方法 |
CN111125969A (zh) * | 2019-12-25 | 2020-05-08 | 华中科技大学 | 一种跨水库流域河道径流演算方法和系统 |
CN112199901A (zh) * | 2020-09-14 | 2021-01-08 | 中国人民解放军火箭军研究院核技术研究所 | 一种无径流资料山区小流域山洪设计暴雨洪水计算方法 |
Non-Patent Citations (3)
Title |
---|
Multi-source error correction for flood forecasting based on dynamic system response curve method;Zhongmin Liang等;《Journal of Hydrology》;20210105;第1-12页 * |
复杂河道实时洪水预报的半自适应模型研究;李致家等;《水文计算与洪水预报》;20000731(第7期);第16-17、50页 * |
水库调度决策支持体系在北盘江梯级水库的应用;肖鹏等;《贵州电力技术》;20150630;第18卷(第6期);第18-22页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112861360A (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112861360B (zh) | 基于系统响应理论的马斯京根流量演算误差修正方法 | |
CN110991687B (zh) | 一种基于经验模型的水资源调度优化方法 | |
CN111914432B (zh) | 一种基于大数据的水文预报方法 | |
CN105243502B (zh) | 一种基于径流区间预测的水电站调度风险评估方法及系统 | |
CN108109076B (zh) | 一种考虑径流预报的梯级水电站群发电调度弃水风险分析方法 | |
CN108830419B (zh) | 一种基于ecc后处理的梯级水库群入库流量联合预报方法 | |
CN107045568B (zh) | 基于动态规划逐次逼近法的河道糙率反演方法 | |
CN108491974B (zh) | 一种基于集合卡尔曼滤波的洪水预报方法 | |
CN109783884B (zh) | 基于面雨量和模型参数同时校正的实时洪水预报误差修正方法 | |
CN114841417B (zh) | 一种高精度咸潮预报方法、系统和可读存储介质 | |
CN101899820A (zh) | 面向河流生态系统保护的流域地表水可利用量确定方法 | |
CN109214110B (zh) | 一种长距离提水工程优化调度方法 | |
CN117033888A (zh) | 一种基于分割基流的流域汇流单位线推求方法 | |
CN112861449B (zh) | 一种基于洪水预报误差反演的多河段联合校正模型 | |
CN113435631B (zh) | 一种洪水预报方法、系统、可读存储介质及计算设备 | |
CN116681180A (zh) | 一种基于两层相似性度量的湖泊水位相似预报方法 | |
CN109059875B (zh) | 一种驱动完美模型开展月尺度径流预报的方法 | |
CN113807745A (zh) | 一种电厂年度计划用水总量确定方法 | |
CN116628915A (zh) | 一种库区沿程流量预测方法 | |
CN112528462A (zh) | 一种基于动态响应曲线对主河道区间入流量修正的方法 | |
CN111932023A (zh) | 一种基于典型设计洪水过程线的小流域短期洪水预报方法 | |
CN112560230A (zh) | 一种提高河道演算稳定性的总变差递减方法 | |
CN116384124B (zh) | 一种基于过程分割的流域设计洪水计算方法 | |
CN113158394B (zh) | 基于蒸发误差补偿的流域水文模型误差修正方法和装置 | |
CN116484590A (zh) | 降雨中心在流域空间移动的情景模拟方法及系统 |
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 |