CN113033014B - 一种考虑多水源联合概率分布的区域可供水量估算方法 - Google Patents
一种考虑多水源联合概率分布的区域可供水量估算方法 Download PDFInfo
- Publication number
- CN113033014B CN113033014B CN202110381969.8A CN202110381969A CN113033014B CN 113033014 B CN113033014 B CN 113033014B CN 202110381969 A CN202110381969 A CN 202110381969A CN 113033014 B CN113033014 B CN 113033014B
- Authority
- CN
- China
- Prior art keywords
- water supply
- runoff
- distribution
- water
- correlation
- 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.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 295
- 238000009826 distribution Methods 0.000 title claims abstract description 112
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000005315 distribution function Methods 0.000 claims abstract description 49
- 241000039077 Copula Species 0.000 claims abstract description 39
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 20
- 238000010586 diagram Methods 0.000 claims description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 2
- 238000007726 management method Methods 0.000 abstract description 5
- 230000000875 corresponding effect Effects 0.000 description 11
- 239000004576 sand Substances 0.000 description 6
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 241001559589 Cullen Species 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013076 uncertainty analysis Methods 0.000 description 2
- 240000001949 Taraxacum officinale Species 0.000 description 1
- 235000005187 Taraxacum officinale ssp. officinale Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012546 transfer Methods 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- 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/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Mathematical Optimization (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Evolutionary Biology (AREA)
- Public Health (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Educational Administration (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Development Economics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
Abstract
本发明涉及水资源规划和管理领域,特别是涉及一种考虑多水源联合概率分布的区域可供水量估算方法。所述方法包括如下步骤:对设定区域的多个供水源进行分析,获取径流量数据;分别计算两两供水源之间的相关系数,判别供水源径流系列的相关性;针对相关性较强的供水源,利用Copula函数,得到联合概率分布函数,对多个供水源进行水文丰枯境遇分析,得到多个供水源不同丰枯境遇下的径流量;针对相关性较弱的供水源,基于P‑III型曲线绘制年径流量频率曲线,得到不同水文保证率下的径流量;结合水库蓄水要求及河流生态基流要求,确定上游水库在丰平枯年份的可供水量。本发明考虑了由于径流量变化引起供水变化的不确定性,为区域水资源的优化调配提供定量依据。
Description
技术领域
本发明涉及水资源规划和管理领域,特别是涉及一种考虑多水源联合概率分布的区域可供水量估算方法。
背景技术
区域可供水量的合理估算是区域水资源优化配置和规划管理的重要基础。目前大多数研究关于区域可供水量的估算一般采用典型年法,通过选取典型水文年流域的来自外调、内流径流量来估测可供水量,对于来水系列资料比较完整的区域,也有采用长系列调算法进行可供水量计算。对于区域内存在的多水源调水,一般根据河流的来水保证率按单项工程可供水量方法估算,扣除重复计算供水量后进行叠加,得到区域总的可供水量。这些方法计算简单,却存在着明显的缺点,即它的准确性低,只能进行粗略估计,未考虑水资源系统的不确定性,未能综合考虑和有效量化不同水源区之间的相关性,未能全面描述区域在多个供水源不同丰枯境遇联合概率下的可供水量,因而不能为流域内、外多水源径流的丰枯互补、相机补水提供定量化依据。
发明内容
有鉴于此,针对现有区域可供水量估算方法未能有效刻画水资源系统的不确定性,无法准确量化区域内多个供水源之间的相关性,以及未能综合考虑多个供水源不同丰枯境遇联合概率的不足,提出一种考虑多水源联合概率分布的区域可供水量估算方法,为进一步实现水资源优化配置提供可靠的基础。
为达到上述目的,本发明主要提供如下技术方案:
一种考虑多水源联合概率分布的区域可供水量估算方法,包括如下步骤:
步骤1:对设定区域的多个供水源进行分析,获取每个所述供水源的径流量数据;
步骤2:分别计算两两供水源之间的相关系数,判别多个供水源径流系列的相关性;
步骤3:针对径流系列相关性较强的供水源,利用单变量的参数估计方法构建其边缘分布函数,根据拟合度检验得到最优的边缘分布形式,再利用Copula函数,将不同的边缘分布连接起来,得到联合概率分布函数,对多个供水源进行水文丰枯境遇分析,得到多个供水源不同丰枯境遇下的径流量;针对和其它径流系列相关性较弱的供水源,基于P-III型曲线绘制年径流量频率曲线,得到不同水文保证率下的径流量;
步骤4:结合水库蓄水要求及河流生态基流要求,确定上游水库在丰平枯年份的可供水量。
上述的区域可供水量估算方法中,步骤3中,所述联合概率分布的表达式为:
C(F1(x1),F2(x2)),u1=F1(x1),u2=F2(x2)
式中:x为供水源径流量,x1为1号供水源径流量,x2为2号供水源径流量,u为单个供水源径流量的边缘分布满足的分布类型,u1为1号供水源径流量满足的分布类型,u2为2号供水源径流量满足的分布类型,C 为多个供水源径流量的联合分布满足的Copula函数分布类型,θ为Copula 函数的未知参数;
上述的区域可供水量估算方法中,所述分布类型为正态分布、对数正态分布、韦伯分布或伽马分布。
上述的区域可供水量估算方法中,获取所述联合概率分布函数的步骤包括:首先要利用极大似然法求解似然函数的最大值,确定单变量边缘分布函数的参数,通过AIC/BIC准则检验正态分布等备选分布函数的拟合优度,当AIC/BIC值越小,说明所选的边缘函数与经验分布拟合的效果越好,以此选定各个供水源径流量最优的边缘分布类型;再利用Copula函数,将不同的边缘分布函数连接起来,得到联合概率分布函数。
上述的区域可供水量估算方法中,步骤3中,选取不同的Copula函数来拟合两两供水源之间的联合概率分布,在平方欧式距离最小的原则下选定最优的Copula函数。
上述的区域可供水量估算方法中,所述单变量是指一个供水源的径流量数据;所述多变量是指多个供水源的径流量数据;
上述的区域可供水量估算方法中,所述Copula函数为BBV Copula、 ClaytonCopula或Frank Copula函数。
上述的区域可供水量估算方法中,步骤3中,得到多个供水源不同丰枯境遇下的径流量的步骤为:根据相关性较强的多个供水源分别处于特丰水年、偏丰水年、平水年、偏枯水年或特枯水年的不同组合形成的丰枯境遇下的径流量,分析所述多个供水源的两两之间径流系列的联合概率分布,从而得到一定概率下的径流量,
上述的区域可供水量估算方法中,所述P-III型曲线包括P=12.5%、 37.5%、62.5%、87.5%的P-III型曲线。
上述的区域可供水量估算方法中,步骤2中,通过Spearman系数、 Pearson系数、Kendall系数检验多个供水源径流系列的相关性,相关系数值的绝对值越大,变量间的相依性越好,其中:相关系数值小于0.2可视为极弱相关或无关,相关系数值大于0.6为强相关。
上述的区域可供水量估算方法中,步骤4中,所述水库蓄水要求为水库蓄水期间需保证河流生态水量和环境保护的基本要求,满足下游必需的供水、生态保护需要的最小流量;
上述的区域可供水量估算方法中,所述生态基流要求是指河道最小流量所应满足的最低要求。
上述的区域可供水量估算方法中,所述步骤3得到的相关性较强的多个供水源丰平枯境遇下的径流量与相关性较弱的单个供水源水文保证率下的径流量之和即为目标区域的预计可供水量。
上述的区域可供水量估算方法中,通过步骤4,根据联合概率分布图和预期目的,计算出特丰、偏丰、平、偏枯和特枯五个不同水文时期所对应概率的年径流量组合,得到不同联合分布概率下区域的预计年径流量区间范围,再结合水库蓄水要求和河流生态基流要求,进而估算得到区域实际可供水量。
借由上述技术方案,本发明区域可供水量估算方法至少具有下列优点:
1)利用Copula函数可求解多变量概率问题、适于表征变量间非线性关系的优势,本发明基于Copula函数研究径流相关性较好的不同供水源所在水文区的丰枯境遇进行联合分布概率分析,得到特丰、偏丰、平、偏枯和特枯五个不同水文时期所对应概率的年径流量组合,得到不同联合分布概率下区域的预计年径流量区间范围;结合基于P-III型曲线研究径流相关性较弱的供水源得到的不同水文保证率下的径流量;综合考虑上述径流量,结合区域规划政策等得到区域可供水量,有效量化了各供水源之间的相关性及多个供水源共存时不同丰枯遭遇下的联合分布概率,全面考虑了由于径流量变化引起供水变化的不确定性,准确性高,为区域水资源优化调配提供定量依据。
2)有效量化了不同水源区之间的相关性,准确刻画了区域在多个供水源不同丰枯境遇联合概率下的可供水量。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
图1是本发明的区域可供水量估算方法的流程图;
图2是本发明的实施例计算得到的黄河-长江径流系列联合概率密度;
图3是本发明的实施例计算得到的黄河-长江径流系列联合概率分布;
图4为Cullen&frey图((a)长江丹江口站径流系列、(b)黄河花园口站径流系列);
图5为大沙河年径流量频率曲线图;
图6为唐河年径流量频率曲线图。
具体实施方式
为更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明申请的具体实施方式、结构、特征及其功效,详细说明如后。
以雄安新区的区域供水为例,对本发明的区域可供水量估算方法进行进一步说明。
雄安新区位于大清河流域,水量来源分为流域内和流域外两部分,具体包括五个方面:上游清水产流、流域内水库联合调水、引黄入冀补淀工程、南水北调工程及再生水。因此,雄安新区的不同供水源分为三个水文区,即大清河流域水文区、长江流域水文区及黄河流域水文区。
如图1所示,本实施例所提供的一种考虑多水源联合概率分布的区域可供水量估算方法,包括如下步骤:步骤1:对设定区域的多个供水源进行分析,获取每个所述供水源的径流量数据;步骤2:分别计算两两供水源之间的相关系数,判别多个供水源径流系列的相关性;步骤3:针对径流系列相关性较强的供水源,利用单变量的参数估计方法构建其边缘分布函数,边缘分布函数各自服从某种分布类型(如对数正态分布、伽玛分布等),其横坐标为径流量数据,纵坐标为对应的概率,得到的是不同径流量所对应的概率,根据拟合度检验得到最优的边缘分布形式,再利用Copula函数,将不同的边缘分布连接起来,得到联合概率分布函数,对多个供水源进行水文丰枯境遇分析,得到多个供水源不同丰枯境遇下的径流量,通过联合概率分布函数,可以分析不同概率水平条件下的径流量,概率水平类比于不同丰枯境遇;针对和其它径流系列相关性较弱的供水源,基于P-III型曲线绘制年径流量频率曲线,得到不同水文保证率下的径流量;步骤4:根据联合概率分析图和预期目的,计算出大于等于某一概率时的年径流量组合,进而计算得到区域可供水量。
本实施例所提供的区域可供水量估算方法从分析区域不同供水源对应的径流系列的水文丰枯境遇入手,判别各径流系列的相关性,对相关性较强的分析边缘概率分布,进行联合概率分布分析,以此来充分考虑供水系统的不确定性;对相关性较弱的绘制径流P-III型曲线。综合各水源通过上述联合概率分布、P-III型曲线估算得到的来自外调、内流的可供水量,得到考虑气候变化、人为活动影响下区域基于不确定性的可供水量,为水资源优化配置和管理提供更加科学合理的可供水量基础数据。该方法有效量化了各供水源之间的相关性及多个供水源共存时不同丰枯遭遇下的联合分布概率,综合考虑了由于径流量变化引起供水变化的不确定性,准确性高,为区域水资源优化调配提供定量依据。
优选地,步骤3中,所述联合概率分布的表达式为:
C(F1(x1),F2(x2)),u1=F1(x1),u2=F2(x2)
式中:x为供水源径流量,x1为1号供水源径流量,x2为2号供水源径流量,u为单个供水源径流量的边缘分布满足的分布类型,u1为1号供水源径流量满足的分布类型,u2为2号供水源径流量满足的分布类型,C 为多个供水源径流量的联合分布满足的Copula函数分布类型,θ为Copula 函数的未知参数。所述分布类型为正态分布、对数正态分布、韦伯分布或伽马分布。
获取所述联合概率分布函数的步骤包括:首先要利用极大似然法求解似然函数的最大值,确定单变量边缘分布函数的参数,通过AIC/BIC准则检验正态分布等备选分布函数的拟合优度,当AIC/BIC值越小,说明所选的边缘函数与经验分布拟合的效果越好,以此选定各个供水源径流量最优的边缘分布类型;再利用Copula函数,将不同的边缘分布函数连接起来,得到联合概率分布函数。
极大似然估计法是通过一个随机样本满足一定的概率分布,但具体参数未知。参数估计是通过多次实验观察结果,并利用这些结果得出参数的近似值。最大似然估计是基于一个给定的参数可以最大化该样本发生概率的思想,放弃选择其他概率较小的样本,直接将该参数作为估计的真值。 AIC准则是一种统计方法,可以判断模型的优劣,包括函数拟合的偏差和参数变量引起的不可靠性。它适用于最大似然估计得到的Copula模型的检验。叶斯信息准则(BIC)虽然不及AIC方法,但它对高估的模型有更敏感的响应。两个变量AIC信息准则如下:
AIC=-2log(MLE)+2k
BIC=-2log(MLE)+klog n
式中,xi为供水源径流量,xi1表示1号供水源径流量,xi2表示2号供水源径流量;F为经验分布函数,其数值为经验频率值;ui为供水源径流量,ui1为1号供水源径流量,ui2为2号供水源径流量;C为选取的分布函数(理论分布函数),其数值为理论频率值;MSE为经验分布函数与理论分布函数之间的误差。
步骤3中,选取不同的Copula函数来拟合两两供水源之间的联合概率分布,在平方欧式距离最小的原则下选定最优的Copula函数。所述单变量是指一个供水源的径流量数据;所述多变量是指多个供水源的径流量数据;所述Copula函数为BBV Copula、ClaytonCopula或Frank Copula函数。
步骤3中,得到多个供水源不同丰枯境遇下的径流量的步骤为:根据相关性较强的多个供水源分别处于丰水年、平水年或枯水年的不同组合形成的丰枯境遇下的径流量,分析所述多个供水源的两两之间径流系列的联合概率分布,从而得到一定概率下的径流量。
步骤3中,得到多个供水源不同丰枯境遇下的径流量的步骤为:根据相关性较强的多个供水源分别处于特丰水年(概率小于等于12.5%)、偏丰水年(概率为12.5-37.5%)、平水年(概率为37.5-62.5%)、偏枯水年(概率为62.5-87.5%)或特枯水年(概率大于等于87.5%)的不同组合形成的丰枯境遇下的径流量,分析所述多个供水源的两两之间径流系列的联合概率分布,计算出特丰、偏丰、平、偏枯和特枯五个不同水文时期所对应概率的年径流量组合,得到不同联合分布概率下区域的预计年径流量区间范围,再结合水库蓄水要求和河流生态基流要求,进而估算得到区域可供水量实际值
步骤3中,所述P-III型曲线包括P=12.5%、37.5%、62.5%、87.5%的 P-III型曲线。
步骤2中,通过Spearman系数、Pearson系数、Kendall系数检验三个水文区的相关性,相关系数值的绝对值越大,变量间的相依性越好,其中:相关系数值小于0.2可视为极弱相关或无关,相关系数值大于0.6为强相关。
本实施例中,计算得到黄河-长江的径流系列的Spearman系数为0.73、 Pearson系数为0.82、Kendall系数为0.60,说明黄河-长江的径流系列强相关;大清河流域径流与这两者相关系数均小于0.07,说明长江-大清河、黄河-大清河的径流系列弱相关。本实施例中的相关系数通过SPSS或者 MATLAB软件,输入数据,相关分析即可得到相关系数的计算结果。
步骤4中,结合水库蓄水要求及河流生态基流要求,确定上游水库在丰平枯年份的可供水量。补水量就是可供水量,之所以称之为补水,是因为本地水不能自给自足,需要外调水。具体步骤为:首先根据水文学方法,计算出不同水文保证率(对应特丰、偏丰、平、偏枯和特枯年份)下的径流量,这个径流量不能全部放出,要综合考虑供水水库的蓄水量要求以及对于河流基本生态基流的要求,确定可以供给的水量范围。计算得到的丰枯境遇下的径流量和水文保证率下的径流量之和,其实就是可供水量。但是有时可能由于水库蓄水要求和河流生态基流要求,会对得到的径流量进行适当增减,进而确定可供水量。为了便于制定后面可操作的供水量分配方案,在这个范围中确定出分别对应于特丰、偏丰、平、偏枯和特枯年份的可供水量。
在实际应用中,得到多个供水源的Copula函数后,利用Copula函数的表达式和供水源的径流量频率分布特征,得到联合概率密度和联合分布图,并利用Matlab软件绘制它们的等高线图。观察联合概率分析图和等高线图,再结合供水源的实际情况,可以估算出大于等于某一概率时的年径流量组合,进而计算得到区域可供水量。
上述步骤详细说明如下:
(1)研究多个供水源的水文丰枯境遇
通过区域规划和水资源调查,确定区域多源供水的来源(即水文区),包括流域内、流域外的调水水源,运用水文水资源分析方法,分别研究多个供水源在不同年份的水文丰枯境遇。
(2)判别多个供水源径流系列的相关性
由于变量之间是否存在相关性,是判断能否采用Copula函数构造双变量联合分布模型的准则。通过Spearman系数、Pearson系数、Kendall系数等检验各水文区径流量之间的相关性。
(3)分析(2)中满足相关性要求的区域多个供水源之间的联合概率分布
针对径流系列相关性较强的供水源,利用单变量的参数估计方法构建其边缘分布函数,根据拟合度检验得到最优的边缘分布形式;根据平方欧式距离最小的原则选取适宜的Copula函数类型进行不同供水源径流系列的联合概率分布分析,以此来充分考虑供水系统的不确定性,得到多水源不同丰枯境遇下的径流量。
拟合度检验又称适应性检验,是对已建立的预测模型进行检验,比较其预测结果与实际发生的符合性。通常情况下,多个预测模型同时进行检验,并选择较好的拟合度进行试验,可以通过比较AIC和BIC值来完成参数的拟合优度评估。拟合优度检验是通过构建可以表征拟合程度的统计量来确定测试拟合模型与样本观察数据点的拟合程度;构建的测试统计量是样本观测值。相应的函数可以计算出所有测试对象的统计值,然后选择一个与之比较的标准,得到测试结论,然后判断拟合模型的拟合程度。常用评估拟合优度的方法有K-S检验、A-D检验等。K-S检验基于累积分布函数,该函数测试分布是否满足理论分布或比较两个经验分布的显著差异。单样本K-S检验用于测试观察到的数据的经验分布是否符合已知的理论分布。双样本K-S测试对两个样本的经验分布函数的位置和形状参数的差异敏感,使其成为比较两个样本的最有用和最常用的非参数方法之一。
(4)分析(2)中不满足相关性要求的单一供水源的P-III型曲线
针对和其它径流系列相关性较弱的供水源,选取20年及以上的径流数据系列,统计河流年径流量的样本均值、变差系数和偏态系数,基于P-III型曲线绘制年径流量频率曲线,得到不同水文保证率下的径流量。
(5)确定基于不确定性分析的区域可供水量
综合考虑多个供水源在不同水文保证率下的径流量,得到基于不确定性分析的区域可供水量。
基于各个供水源的径流量数据,构建它们的边缘分布函数(U和V),它们各自服从对数正态分布和伽玛分布,其横坐标为径流量数据,纵坐标为对应的概率,得到的是不同径流量所对应的概率;再构建它们的联合分布函数(Clayton Copula),0.9691是参数估计结果,联合分布函数得到的是多个供水源的组合概率,进而可分析不同的水文丰枯遭遇下的径流量。
通过上述区域可供水量估算方法计算的结果如下:
1)相关性判断
通过前期数据收集获得不同水文区径流数据有:黄河花园口站年径流系列(黄河水,1919-2013)、长江丹江口站年径流系列(中线水,1980-2005) 以及西大洋-王快站年径流系列(流域内,1956-2000)。取1980-2000年20 年的径流数据,计算流域内径流、黄河水、南水北调中线水两两间的秩相关系数(见下表1)。结果显示,黄河水和中线水径流的相关系数Spearman 系数为0.73、Pearson系数为0.82、Kendall为0.60,具有强相关性,可以用Copula函数进行联合概率分布;而西大洋-王快径流系列(流域内) 与黄河水、长江水都呈弱相关性,不适合用Copula函数进行联合概率分布计算,因而对其构建单变量分布函数,绘制P-III型概率密度曲线。
表1:不同水文区径流系列的秩相关系数
2)单变量边缘分布函数的构造
选取正态分布、对数正态分布、韦伯分布、伽马分布4 种分布函数,分别绘制长江径流系列、黄河径流系列的Cullen&frey图(见图4),统计 AIC/BIC的结果如表2所示。
对数正态分布的AIC/BIC最小且通过K-S检验,因而长江径流系列、黄河径流系列的边缘分布函数为对数正态分布。
表2:单变量分布函数AIC/BIC计算结果
3)长江-黄河径流系列联合概率分布
选择BBV Copula、Clayton Copula、Frank Copula、Gaussian Copula、 GumbelCopula这5种Copula函数进行拟合,计算平方欧式距离,确定最优Copula为ClaytonCopula。黄河-长江概率密度图如图2所示,联合概率累积分布如图3所示。分析黄河-长江径流的概率密度曲线可知,其存在年径流相同的点,即有点分布在y=x的曲线上,这一点比较特殊,但实际上长江的径流量远大于黄河的径流量,因而在y=x曲线的下部是更符合实际的概率分布。
进一步分析图3的联合概率累积分布图,可以得到大于等于某一概率时候的年径流量组合,进而可以得到特丰水年、偏丰水年、平水年、偏枯水年和特枯水年不同丰枯境遇下的径流量区间范围。当联合概率为12.5%时,黄河-长江可向雄安新区输入的年径流量为[6.7,8.9]亿m3;当联合概率为37.5%时,黄河-长江可向雄安新区输入的年径流量为[5.8,7.2]亿m3;当联合概率为62.5%时,黄河-长江可向雄安新区输入的年径流量为[4.3,6.1]亿m3;当联合概率为87.5%时,黄河-长江可向新区输入的年径流量为[3.2, 4.6]亿m3。因此,在不同丰平枯遭遇下,黄河-长江可向雄安新区输入的年径流量为[3.2,8.9]亿m3。实际规划决策结果表明,由于考虑水库蓄水要求和河流生态基流要求,现有国家引黄济淀工程分配2亿m3/年到新区,以及南水北调中线工程分配4亿m3/年到雄安新区,即黄河-长江向雄安新区输入的年径流量实际为6亿m3,符合本研究的研究结果。
4)唐河-大沙河径流频率统计
基于唐河(西大洋站)、大沙河(王快站)年径流数据(1956-2000年),统计2条河流年径流量的样本均值、变差系数Cv和偏态系数Cs。根据均值、Cv、Cs求P-III型密度函数中的参数α、β、a0,利用Matlab绘制年径流量频率曲线,结果如图5、图6所示。采用非线性最小二乘法优化适线,优化后两个径流系列的R2均达0.99。
大沙河径流P-III型曲线适线表
单位:亿m3
唐河径流P-III型曲线适线表
单位:亿m3
进一步对大沙河和唐河的径流平率曲线进行分析,大沙河的平均流量为6.24亿m3,Cv值为0.74,Cs/Cv为2,当保证率P=12.5%时,径流量为11.8亿m3;当保证率为37.5%时,径流量为5.6亿m3;当保证率为 62.5%时,径流量为3.9亿m3;当保证率为82.5%时,径流量为2.4亿m 3。而唐河的相关统计数据如下:样本值为45,均值为4.66亿m3,Cv为0.67,当保证率P取12.5%、37.5%、62.5%和82.5%时,径流量分别为8.7 亿m3、4.7亿m3、2.6亿m3、2.1亿m3。结合王快水库、西大洋水库的蓄水量要求以及对于河流基本生态基流的要求,确定上游水库在特丰、偏丰、平、偏枯和特枯年份可向雄安新区补水3.6亿m3、2.8亿m3、2亿m 3、1.5亿m3、0.4亿m3。
综上分析,通过对雄安新区不同供水源所在水文区的丰枯境遇进行分析,通过相关性检验结果显示,新区所在的大清河流域水文区与外调水所在的长江水文区、黄河水文区的径流系列相关性弱,而长江径流与黄河径流存在强相关性,优选出拟合度最好的Clayton copula对其联合概率分布进行描述。结果表明,在不同丰平枯遭遇下,黄河-长江可向雄安新区输入的年径流量为[3.2,8.9]亿m3现有国家引黄济淀工程以及南水北调中线工程分别分配2亿m3/年和4亿m3/年到新区,在合理的范围内,符合本研究的研究成果。而大清河流域径流选取大沙河和唐河进行P-III型曲线绘制,结合水库蓄水要求及河流生态基流要求,确定上游水库在特丰、偏丰、平、偏枯和特枯年份的可供水量分别为3.6亿m3、2.8亿m3、2亿m3、1.5 亿m3、0.4亿m3。
由于不同水源的径流量存在同丰、同枯、一丰一枯多种情况,运用传统方法仅考虑河流不同来水保证率下的径流量,扣除重复计算供水量后进行简单叠加,未考虑多水源共存时丰枯境遇不同而引起的联合概率分布,因而结果误差较大。本发明克服了这一局限,能客观反映多水源共存时联合概率分布的变化,为水资源优化配置和供水调度规划决策提供可靠的基础数据支持,从而提高了区域管理的科学性和有效性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.一种考虑多水源联合概率分布的区域可供水量估算方法,其特征在于,包括如下步骤:
步骤1:对设定区域的多个供水源进行分析,获取每个所述供水源的径流量数据;
步骤2:分别计算两两供水源之间的相关系数,判别多个供水源径流系列的相关性;
步骤3:针对径流系列相关性较强的供水源,利用单变量的参数估计方法构建其边缘分布函数,根据拟合度检验得到最优的边缘分布形式,再利用Copula函数,将不同的边缘分布连接起来,得到联合概率分布函数,对多个供水源进行水文丰枯境遇分析,得到多个供水源不同丰枯境遇下的径流量;针对和其它径流系列相关性较弱的供水源,基于P-III型曲线绘制年径流量频率曲线,得到不同水文保证率下的径流量;
步骤4:结合水库蓄水要求及河流生态基流要求,确定上游水库在丰平枯年份的可供水量;
其中,步骤3中,得到多个供水源不同丰枯境遇下的径流量的步骤为:根据相关性较强的多个供水源分别处于特丰水年、偏丰水年、平水年、偏枯水年或特枯水年的不同组合形成的丰枯境遇下的径流量,分析所述多个供水源的两两之间径流系列的联合概率分布,从而得到一定概率下的径流量,
所述P-III型曲线包括P=12.5%、37.5%、62.5%、87.5%的P-III型曲线。
2.根据权利要求1所述的区域可供水量估算方法,其特征在于,
步骤3中,所述联合概率分布的表达式为:
C(F1(x1),F2(x2)),u1=F1(x1),u2=F2(x2)
式中:x为供水源径流量,x1为1号供水源径流量,x2为2号供水源径流量,u为单个供水源径流量的边缘分布满足的分布类型,u1为1号供水源径流量满足的分布类型,u2为2号供水源径流量满足的分布类型,C为多个供水源径流量的联合分布满足的Copula函数分布类型,θ为Copula函数的未知参数;
所述分布类型为正态分布、对数正态分布、韦伯分布或伽马分布。
3.根据权利要求2所述的区域可供水量估算方法,其特征在于,
获取所述联合概率分布函数的步骤包括:首先要利用极大似然法求解似然函数的最大值,确定单变量边缘分布函数的参数,通过AIC/BIC准则检验正态分布等备选分布函数的拟合优度,当AIC/BIC值越小,说明所选的边缘函数与经验分布拟合的效果越好,以此选定各个供水源径流量最优的边缘分布类型;再利用Copula函数,将不同的边缘分布函数连接起来,得到联合概率分布函数。
4.根据权利要求1所述的区域可供水量估算方法,其特征在于,
步骤3中,选取不同的Copula函数来拟合两两供水源之间的联合概率分布,在平方欧式距离最小的原则下选定最优的Copula函数。
5.根据权利要求1所述的区域可供水量估算方法,其特征在于,
所述单变量是指一个供水源的径流量数据;多变量是指多个供水源的径流量数据;
所述Copula函数为BBV Copula、Clayton Copula或Frank Copula函数。
6.根据权利要求1所述的区域可供水量估算方法,其特征在于,
步骤2中,通过Spearman系数、Pearson系数、Kendall系数检验多个供水源径流系列的相关性,相关系数值的绝对值越大,变量间的相依性越好,其中:相关系数值小于0.2可视为极弱相关或无关,相关系数值大于0.6为强相关。
7.根据权利要求1所述的区域可供水量估算方法,其特征在于,
步骤4中,所述水库蓄水要求为水库蓄水期间需保证河流生态水量和环境保护的基本要求,满足下游必需的供水、生态保护需要的最小流量;
所述生态基流要求是指河道最小流量所应满足的最低要求。
8.根据权利要求1所述的区域可供水量估算方法,其特征在于,
所述步骤3得到的相关性较强的多个供水源丰平枯境遇下的径流量与相关性较弱的单个供水源水文保证率下的径流量之和即为目标区域的预计可供水量。
9.根据权利要求1所述的区域可供水量估算方法,其特征在于,通过步骤4,
根据联合概率分布图和预期目的,计算出多个供水源在特丰、偏丰、平、偏枯和特枯五个不同水文时期所对应概率的年径流量组合,得到不同联合分布概率下区域的预计年径流量区间范围,再结合水库蓄水要求和河流生态基流要求,进而估算得到区域实际可供水量。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110381969.8A CN113033014B (zh) | 2021-04-09 | 2021-04-09 | 一种考虑多水源联合概率分布的区域可供水量估算方法 |
LU501097A LU501097B1 (en) | 2021-04-09 | 2021-12-27 | An estimation method of regional available water supply amount with considering the joint probability distributions of multiple water sources |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110381969.8A CN113033014B (zh) | 2021-04-09 | 2021-04-09 | 一种考虑多水源联合概率分布的区域可供水量估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113033014A CN113033014A (zh) | 2021-06-25 |
CN113033014B true CN113033014B (zh) | 2022-05-03 |
Family
ID=76456395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110381969.8A Expired - Fee Related CN113033014B (zh) | 2021-04-09 | 2021-04-09 | 一种考虑多水源联合概率分布的区域可供水量估算方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113033014B (zh) |
LU (1) | LU501097B1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113688542A (zh) * | 2021-10-26 | 2021-11-23 | 长江水利委员会长江科学院 | 智能优化水资源配置方法、装置、计算机设备及存储介质 |
CN114611846B (zh) * | 2022-05-16 | 2022-07-26 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种城市新区多模式生态补水系统及方法 |
CN115860421B (zh) * | 2022-12-26 | 2023-08-29 | 长江水利委员会水文局 | 适配目标保证率的生态流量动态计算方法及系统 |
CN117575297B (zh) * | 2024-01-16 | 2024-04-23 | 水利部水利水电规划设计总院 | 区域水网系统调水受水区供水水源组成的识别方法及系统 |
CN117744501B (zh) * | 2024-01-19 | 2024-05-28 | 水利部水利水电规划设计总院 | 考虑生态流量的水网系统调蓄节点优化调度与决策方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335603A (zh) * | 2015-07-17 | 2016-02-17 | 清华大学 | 一种度量引水灌溉地区干旱程度的方法 |
CN106599562A (zh) * | 2016-12-06 | 2017-04-26 | 武汉大学 | 基于概率加权fdc法的河流生态需水量计算方法 |
CN109543147A (zh) * | 2018-10-26 | 2019-03-29 | 西安理工大学 | 一种流域降雨径流关系非线性快速诊断及归因分析的方法 |
CN110276150A (zh) * | 2019-06-27 | 2019-09-24 | 江西省水利科学研究院 | 一种基于Copula函数的山丘区河川基流量系列插补延长方法 |
CN111709564A (zh) * | 2020-06-08 | 2020-09-25 | 黄河水利委员会黄河水利科学研究院 | 一种气候变化条件下的流量历时曲线动态预测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867106A (zh) * | 2012-08-14 | 2013-01-09 | 贵州乌江水电开发有限责任公司 | 短期径流预报方法及系统 |
KR101705247B1 (ko) * | 2015-02-02 | 2017-02-13 | 한국건설기술연구원 | 사회기반시설물의 성능 평가 시스템 및 그 방법 |
CN111611692B (zh) * | 2020-04-26 | 2022-08-05 | 武汉大学 | 气候变化情景下基于等可靠度的设计洪水推求方法及系统 |
-
2021
- 2021-04-09 CN CN202110381969.8A patent/CN113033014B/zh not_active Expired - Fee Related
- 2021-12-27 LU LU501097A patent/LU501097B1/en active IP Right Grant
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335603A (zh) * | 2015-07-17 | 2016-02-17 | 清华大学 | 一种度量引水灌溉地区干旱程度的方法 |
CN106599562A (zh) * | 2016-12-06 | 2017-04-26 | 武汉大学 | 基于概率加权fdc法的河流生态需水量计算方法 |
CN109543147A (zh) * | 2018-10-26 | 2019-03-29 | 西安理工大学 | 一种流域降雨径流关系非线性快速诊断及归因分析的方法 |
CN110276150A (zh) * | 2019-06-27 | 2019-09-24 | 江西省水利科学研究院 | 一种基于Copula函数的山丘区河川基流量系列插补延长方法 |
CN111709564A (zh) * | 2020-06-08 | 2020-09-25 | 黄河水利委员会黄河水利科学研究院 | 一种气候变化条件下的流量历时曲线动态预测方法 |
Non-Patent Citations (3)
Title |
---|
"Monthly Runoff Prediction of Yunlong Reservoir Based on Copula Function";Chen JingGu等;《Water Resources and Power》;20160325;第34卷(第3期);24-27 * |
"考虑不确定性的流域水文过程及水库调度研究";郭爱军;《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》;20190815(第08期);A012-1,正文第33、36、40、118页 * |
"辽浑太流域大型水库供水联合调度模型及规则研究";王强;《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑》;20170315(第03期);C037-21,正文第22-31、68、76页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113033014A (zh) | 2021-06-25 |
LU501097B1 (en) | 2022-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113033014B (zh) | 一种考虑多水源联合概率分布的区域可供水量估算方法 | |
Yalcintas et al. | An energy benchmarking model based on artificial neural network method utilizing US Commercial Buildings Energy Consumption Survey (CBECS) database | |
Stephenson et al. | Bayesian inference for extremes: accounting for the three extremal types | |
CN108918815B (zh) | 一种土壤重金属风险预测方法 | |
Akay et al. | Flash flood potential prioritization of sub-basins in an ungauged basin in Turkey using traditional multi-criteria decision-making methods | |
CN107545103A (zh) | 煤矿区土壤重金属含量空间模型建立方法 | |
CN111898831B (zh) | 一种实时洪水概率预报实用化方法 | |
CN112506990A (zh) | 一种基于时空信息的水文数据异常检测方法 | |
CN110728409B (zh) | 一种洪水过程类型相似性挖掘与快速预测方法 | |
CN109190865B (zh) | 水质指标评价方法、预测方法、水域水质等级评价方法 | |
CN107886160B (zh) | 一种bp神经网络区间需水预测方法 | |
Zou et al. | Spatial interpolation of the extreme hourly precipitation at different return levels in the Haihe River basin | |
Hong et al. | Estimating extreme wind speed based on regional frequency analysis | |
Baasch et al. | Comparing gray box methods to derive building properties from smart thermostat data | |
CN116049340A (zh) | 考虑水位与水文连通性共存关系的湖泊最优水位确定方法 | |
CN117116382B (zh) | 引水工程影响下受水湖泊水质时空预测方法和系统 | |
CN115689051A (zh) | 一种基于GA算法耦合Morris和GLUE的SWMM模型参数自动率定方法 | |
CN111090831A (zh) | 一种湖泊面积变化关键驱动因子识别方法 | |
CN115907286A (zh) | 一种社会经济干旱评估方法 | |
CN106227965B (zh) | 一种顾及时空分布非平稳特征的土壤有机碳空间抽样网络设计方法 | |
CN112116200A (zh) | 基于动态比例代换和层次贝叶斯的缺灾情资料城市洪灾损失函数的构建方法 | |
CN117078077A (zh) | 一种高速公路路域的生态脆弱性评价方法 | |
Ye | Spatial variation and interpolation of wind speed statistics and its implication in design wind load | |
CN116777027A (zh) | 一种非正常日的负荷预测方法及系统 | |
CN115618720A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220503 |