CN111814314A - 基于稳定性评价的水电站补偿分摊方案评价方法 - Google Patents

基于稳定性评价的水电站补偿分摊方案评价方法 Download PDF

Info

Publication number
CN111814314A
CN111814314A CN202010550888.1A CN202010550888A CN111814314A CN 111814314 A CN111814314 A CN 111814314A CN 202010550888 A CN202010550888 A CN 202010550888A CN 111814314 A CN111814314 A CN 111814314A
Authority
CN
China
Prior art keywords
hydropower station
value
hydropower
compensation allocation
evaluation
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.)
Granted
Application number
CN202010550888.1A
Other languages
English (en)
Other versions
CN111814314B (zh
Inventor
李英海
汪利
夏青青
郭家力
涂玉律
聂盼盼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN202010550888.1A priority Critical patent/CN111814314B/zh
Priority to CN202310964511.4A priority patent/CN117081080A/zh
Publication of CN111814314A publication Critical patent/CN111814314A/zh
Application granted granted Critical
Publication of CN111814314B publication Critical patent/CN111814314B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Educational Administration (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了基于稳定性评价的水电站补偿分摊方案评价方法,包括:建立水电站发电优化调度模型,比较不同联盟方式下各水电站发电量;利用不同的补偿分摊方法,进行梯级水电站补偿分摊;对不同补偿分摊方法的分摊结果进行稳定性评价;建立指标体系,构造评价指标值矩阵,制定梯级水电站补偿分摊方案;基于后悔理论,计算得出最优的补偿分摊方案。本发明对梯级水电站制定出多个补偿分摊方案,将稳定性作为指标之一,建立了指标体系,对补偿分摊方案进行综合评价,选出的最优补偿分摊方案既具有稳定性,又实现了梯级水电站的整体最优,可促进梯级水电站的联合调度与相互协作。

Description

基于稳定性评价的水电站补偿分摊方案评价方法
技术领域
本发明属于水电站调度领域,具体涉及一种基于稳定性评价的水电站补偿分摊方案评价方法。
背景技术
在我国水能资源大规模开发格局下,各大流域梯级水电站群已初具规模,水电站群不再是片面追求单个水电站的发电量最大,而是从梯级整体出发,合理利用有限的水能资源以产生尽可能大的综合效益。然而,在流域梯级开发中,梯级上下游各水电站可能分属不同业主,在联合调度运行时,上游龙头水电站及其他具有较大径流调节能力的水电站往往需要牺牲自身的部分利益来提高下游电站的发电效益,从而使整个梯级达到最优,为促进梯级水电站的联合调度与相互协作,这就涉及多业主梯级水电站间补偿效益分摊的问题。
目前,梯级水电站补偿效益分摊主要的方法包括:条约、协议或规则、单指标法、综合指标法、比例分配法、Shapley值法等。通过较多的实例分析可知,多种分摊方法均可应用于梯级水电站补偿效益分摊中,而如何从中选出最优方案任然是亟待解决的现实问题。
鉴于此,研究水电站补偿分摊方案评价方法,建立汉江中游南河流域上的三里坪、寺坪、过渡湾梯级水电站补偿效益分摊方案的指标体系,基于主客观综合赋权法和后悔理论对多种梯级水电站补偿效益分摊方案进行综合评价。
发明内容
本发明的技术问题是现有的梯级水电站补偿分摊方法有多种,根据这些方法制定的梯级水电站补偿分摊方案各有优缺点,目前缺乏对梯级水电站补偿分摊方案的评价方法,如何制定公平合理的补偿分摊方案并从制定的多个补偿分摊方案中选出最优方案仍然是亟待解决的问题。
本发明的目的是解决上述问题,提供一种基于稳定性评价的水电站补偿分摊方案评价方法,利用不同补偿分摊方法进行梯级水电站补偿分摊,并对补偿分摊结果进行稳定性评价,通过建立指标体系制定水电站的补偿分摊方案,基于后悔理论选出最优的补偿分摊方案,促进梯级水电站的联合调度与相互协作,使梯级水电站整体达到最优。
本发明的技术方案是基于稳定性评价的水电站补偿分摊方案评价方法,包括以下步骤,
步骤1:建立水电站发电优化调度模型,比较不同联盟方式下各水电站发电量;
步骤2:利用不同的补偿分摊方法,进行梯级水电站补偿分摊;
步骤3:对步骤2中的补偿分摊结果进行稳定性评价;
步骤4:建立指标体系,构造评价指标值矩阵,制定梯级水电站补偿分摊方案;
步骤5:基于后悔理论,根据步骤4的各补偿分摊方案的评价指标值,计算得出最优的补偿分摊方案。
优选地,步骤2中,所述不同的补偿分摊方法,包括贡献度法、变异系数法、Shapley值法、变异系数-Shapley值法。
优选地,步骤3中,所述对步骤2中的补偿分摊结果进行稳定性评价,采用分裂倾向法。
进一步地,步骤4还包括运用专家评分法和熵权法计算指标体系的各个指标的权重。
优选地,步骤4的指标体系包括分摊比例、梯级增加发电量、最小出力、分裂倾向值。
进一步地,步骤1包括以下子步骤,
步骤1.1:以单一水电站发电效益最大或者考虑不同联盟方式下发电效益最大为目标函数,分别建立发电优化调度模型:
1)单一水电站年发电效益最大:
Figure BDA0002542474670000021
2)不同联盟方式下年发电效益最大:
Figure BDA0002542474670000022
其中B为调度期内联盟水电站总发电效益;K为联盟水电站总数;i为水电站编号,表示梯级中第i个水电站;t、T分别为调度期各时段编号及总时段数;Ai为水电站i的综合出力系数;Qi,t为水电站i在时段t的发电引用流量;Hi,t为水电站i在时段t的发电水头;Δt为时段长;pi为水电站i的上网电价;
步骤1.2:确定约束方程,主要包括水量平衡方程、水电站蓄水量约束、水库水位约束、水电站出力约束、水库下泄流量约束以及上下游水电站间流量联系;
步骤1.3:采用动态规划方法进行发电调度优化计算,得到不同联盟方式下各水电站的发电量及发电效益。
进一步地,步骤5包括以下子步骤,
步骤5.1:构造补偿分摊方案的评价指标值矩阵并对其进行标准化处理;
步骤5.2:针对每项指标,取标准化后矩阵中指标评价值的最优值,建立理想方案;
步骤5.3:构建后悔值矩阵;
步骤5.4:构建感知效用矩阵;
步骤5.5:结合各指标的权重,计算各补偿分摊方案的综合感知效用值;
步骤5.6:比较各补偿分摊方案的综合感知效用值,选出最优的补偿分摊方案。
相比现有技术,本发明的有益效果:
1)本发明对梯级水电站制定出多个补偿分摊方案,将补偿分摊结果的稳定性作为指标之一,建立了补偿分摊方案的指标体系,对补偿分摊方案进行综合评价,便于选出梯级水电站最优补偿分摊方案,选出的最优补偿分摊方案既具有稳定性,又实现了梯级水电站的整体最优,促进了梯级水电站的联合调度与相互协作;
2)本发明提出的变异系数-Shapley值法既考虑了电站的个体特征又考虑了电站的边际贡献,在充分肯定了龙头电站在联盟中的作用与价值的同时,又促进了下游各级电站积极参与梯级联合调度;
3)本发明将后悔理论与效用相结合,真实地反应了水电站业主的决策行为,增加了最优补偿分摊方案的准确性、科学性。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1为基于稳定性评价的水电站补偿分摊方案评价方法的流程示意图。
图2为实施例的梯级水电站的示意图。
图3为实施例的不同补偿分摊方法的对比图。
具体实施方式
如图1所示,基于稳定性评价的水电站补偿分摊方案评价方法,包括以下步骤,
步骤1:建立水电站发电优化调度模型,计算出不同联盟方式下各水电站发电效益;
步骤1.1:以单一水电站发电效益最大或者考虑不同联盟方式下发电效益最大为目标函数,分别建立发电优化调度模型:
①单一水电站年发电效益最大:
Figure BDA0002542474670000031
②不同联盟方式下年发电效益最大:
Figure BDA0002542474670000032
其中B为为调度期内联盟水电站总发电效益;K为联盟水电站总数;i为水电站编号,表示梯级中第i个水电站;t、T分别为调度期各时段编号及总时段数;Ai为水电站i的综合出力系数;Qi,t为水电站i在时段t的发电引用流量;Hi,t为水电站i在时段t的发电水头;Δt为时段长;pi为水电站i的上网电价。
步骤1.2:确定约束方程,主要包括水量平衡方程、水电站蓄水量约束、水库水位约束、水电站出力约束、水库下泄流量约束以及上下游水电站间流量联系等。
步骤1.3:采用动态规划方法进行发电调度优化计算,得到不同联盟方式下各水电站的发电量及发电效益。
步骤2:分别采用贡献度法、变异系数法、Shapley值法、变异系数-Shapley值法进行补偿效益分摊;
贡献度法是指在梯级水电站联合优化调度的发电收益中,电站单位水能转化为电能进而转化为发电效益的效率越高,对联合调度发电收益的贡献度越大,计算式如下:
Figure BDA0002542474670000041
其中xi为第i个水电站的分摊比例;n为梯级水电站总数;piki为第i个水电站势能转化为发电收益的转化系数;ki发电效率系数,代表势能转化为电能的系数;pi为上网电价,代表电能转化为发电效益的系数。
变异系数法是一种客观赋权的方法,该方法根据水电站各个体特征参数在所有被评价的个体特征参数上观测值的变异程度大小来对其赋权。假设有n个水电站,每个水电站有m个个体特征参数,从而形成了一个参数值矩阵X=(Xij)n×m,计算过程如下:
①先对各项参数值进行归一化,再对参数值矩阵进行标准化处理,为使每个参数值均大于0,且在[0.6,1]范围内,对于正向参数有:
Figure BDA0002542474670000042
其中,Yij为第i个水电站的特征参数j的标准值,Yij∈[0.6,1];Xmax j、Xmin j为特征参数j的最大、最小特征值。
②计算第j个特征参数值的平均值和均方差:
Figure BDA0002542474670000051
Figure BDA0002542474670000052
③计算第j个特征参数值的变异系数:
Figure BDA0002542474670000053
④计算第j个特征参数权重值:
Figure BDA0002542474670000054
⑤计算第i个水电站的第j个特征参数值的比重:
Figure BDA0002542474670000055
⑥计算第i个水电站综合权重值:
Figure BDA0002542474670000056
其中,υi为第i个水电站的个体特征权重系数,即分摊比例。
Shapley值法是一种用于解决多方合作对策利益分配问题的数学方法,该方法按照参与人员边际费用的期望值进行效益分摊,以保证分摊的有效性与公平性。
Shapley值使用的基础公理有对称公理、有效公理、可加公理。当满足上述基础公理后,各水电站的Shapley值具有唯一解:
Figure BDA0002542474670000057
其中,s为联盟中水电站个数;N为所有可能的联盟方式;S为N的一个子联盟;(s-1)!(n-s)!/n!为联盟S出现的联盟概率;ν(S)-ν(S-i)为电站i对联盟S的边际贡献值。
结合变异系数法和Shapley值法各自优势,提出变异系数-Shapley值法,根据公式(4)~(10)得到各水电站个体特征权重系数υi,根据公式(11)得到各水电站的Shapley值σi,则变异系数-Shapley值法公式如下:
θi=υiσi (12)
Figure BDA0002542474670000061
其中xi为第i个水电站的分摊比例;υi为第i个水电站的个体特征权重系数;σi为第i个水电站的的Shapley值。
步骤3:用分裂倾向(Propensity to disrupt,PTD)法对步骤2的4种补偿分摊方法的分摊结果进行稳定性评价;
分裂倾向法是一种定量评价合作博弈解决方案稳定性的常用方法,是在某种分配结果下大联盟失去成员i所遭受的损失与成员i离开大联盟自身遭受的损失之比:
Figure BDA0002542474670000062
其中PTDi为某种分配结果下成员i的分裂倾向;ci为大联盟中成员i分配得到的收益;v(i)为成员i实行单库调度时取得的收益;V(N-{i})为失去成员i的大联盟总收益。
步骤4:以分摊比例、梯级增加发电量、最小出力、PTD值为指标建立指标体系,制定4种补偿分摊方案,构造指标评价值矩阵;
步骤5:基于后悔理论对步骤4制定的4种分摊方案进行评价,得出最优的分摊方案。
后悔理论是一种在放弃独立性公理的前提下,在进行决策过程时,将后悔和欣喜这两种不同的心理行为同时考虑到该项决策中的理论,具体表现为一种感知效用函数。此函数主要分为两部分:一是当前选择方案的效用函数,二是选择此方案的后悔-欣喜函数。令a和b分别为选择A方案和B方案得到的结果,那么决策者选择方案A的感知效用函数为:
q(a,b)=v(a)+R(v(a)-v(b)) (15)
函数v(a)和v(b)分别为选择A方案和B方案得到的效用函数,而函数R(v(a)-v(b))为后悔-欣喜函数,当后悔-欣喜函数值大于0时,表示为选择A方案和不选B方案感到喜悦;当后悔-欣喜函数值小于0时,表示为选择A方案和不选B方案感到后悔。
对于多个不同的方案而言,决策者的感知效用函数如下
qi=v(ai)+R(v(ai)-v(a*)) (16)
其中a*=max{ai,i=1,2,…,m};R(v(ai)-v(a*))<0,即后悔-欣喜函数值为负数,表示决策者的后悔程度。
分摊方案的评价过程包括:
假设有F个补偿效益分摊方案组成决策方案集M=(M1,M2,…,MF),L个决策指标组成指标集I=(I1,I2,…,IL),其中分裂倾向(PTD)为其中一个指标,指标权重为W=(w1,w2,…,wL)。①构造评价指标值矩阵及其标准化处理
方案集M对指标集I的评价值的评价指标值矩阵为:
Figure BDA0002542474670000071
由于各个决策指标之间存在性质或量级上的不同,所以需要先对样本矩阵进行标准化处理:为使每个参数值均大于0,且在[0.6,1]范围内,对于效益型指标有:
Figure BDA0002542474670000072
对于成本型指标有:
Figure BDA0002542474670000073
其中,xij表示决策者在Ij指标下对方案Mi的评价值;xmaxj表示决策者在Ij指标下对方案Mi的评价值中最大的数值;xminj表示决策者在Ij指标下对方案Mi的评价值中最小的数值;yij为Ij指标下方案Mi标准值,yij∈[0.6,1],(i=1,2,…,F;j=1,2,…,L)。
②建立理想方案
设xj +=max(xij),(i=1,2,…,F;j=1,2,…,L),称
xj +=(x1 +,x2 +,…,xL +) (20)
为决策的理想方案;
③构建后悔值矩阵
计算不同指标值xij的函数效应值vij,需要先构造效用函数v(x)。鉴于现实生活中进行风险决策时,决策者多尽量规避风险,故效用函数v(x)是单调递增的凹函数,满足v'(x)>0和v”(x)<0。运用幂函数代表效用函数,计算式如下
v(x)=xα (21)
其中1>α>0,α越小,表明决策者越规避风险。
④构建感知效用矩阵
分别计算各方案的后悔值,需要构建后悔-欣喜函数R(·)。鉴于决策者对风险性决策是规避风险的,故后悔-欣喜函数R(·)是单调递增的凹函数,即后悔-欣喜函数满足:R'(·)>0和R”(·)<0。决策者的后悔-欣喜函数值的计算式如下
Figure BDA0002542474670000081
其中ε是决策者的后悔规避系数,ε>0,ε越大,后悔规避程度就越大;vab为2个不同方案的效用值的差值。
再由式(8)和(9)构建决策者的效用感知函数矩阵为
qij=vij+Rij (23)
⑤通过专家打分法和熵权法确定主客观综合权重:
专家评分法是一种定性描述定量化方法,它首先根据评价对象的具体要求选定若干个评价指标,再根据评价指标制订出评价标准,聘请若干代表性专家凭借自己的经验按此评价标准给出各指标的评价分值,然后对其进行结集。专家评分法计算的指标j的权重表示为wzj
熵权法是综合指标法中的一种。信息是系统内有序程度的一种度量,熵则是系统内无序程度的一种度量,两者符号相反,绝对值相等。运用熵可以度量和评价指标数据所含的信息量,同时可以此确定各项指标的权重值。
熵权法计算权重的步骤如下:
a.求各指标的信息熵
Figure BDA0002542474670000082
其中ej表示评价指标j的信息熵,Uij表示第i个方案的第j个评价指标值的比例,
Figure BDA0002542474670000083
yij为第i个方案的评价指标j的标准值;
b.确定各指标权重
Figure BDA0002542474670000091
其中wkj表示熵权法计算的评价指标j的权重,gj表示第j个特征指标差异系数,gj=1-ej
综合主客观权重:
Figure BDA0002542474670000092
wj为评价指标j的综合主客观权重,参数
Figure BDA0002542474670000093
根据决策者的偏好信息确定。
⑥计算方案Mi在不同指标之下的综合感知效用值,计算公式如下:
Figure BDA0002542474670000094
实施例选择南河中下游干流上的三里坪水电站A、寺坪水电站B、过渡湾水电站C组成的梯级水电站作为实施对象,以三里坪上游50%平水年逐月径流为年径流过程,进行补偿分摊计算。3座电站的位置关系如图2所示。
实施结果如下:
(1)采用动态规划法对步骤1的发电优化调度模型进行求解,得到不同联盟方式下三里坪、寺坪、过渡湾3座水电站的发电量及发电效益,如表1所示。
表1不同联盟方式下各水电站的发电量及发电效益表
Figure BDA0002542474670000095
三里坪、寺坪、过渡湾3座水电站联合运行的梯级总发电量相比3座水电站单独运行的总发电量增加0.738亿kwh,总发电效益增加了0.290亿元,按照步骤2的4个补偿分摊方法来进行补偿分摊比例计算,结果如图3和表2所示。
表2各水电站补偿分摊比例表
Figure BDA0002542474670000096
从表2可以看出,三里坪水电站作为龙头电站,因其具有补偿调节能力,从而大幅提高了梯级整体发电效益,在进行补偿效益分摊时也获得了最高的补偿效益,这充分肯定了龙头电站的补偿调节作用;寺坪、过渡湾水电站在获得比各电站单独运行时高的补偿效益的同时,二者的分摊差异也很大,由此可见次级电站间的个体差异,以及选择有效的分摊方法有利于促进下游各级电站积极参与联合调度。
(2)计算分裂倾向值,结果如表3所示。水电站i的PTD值越大,则它分得的收益相较于它对联盟所做的贡献越小,它越倾向于离开联盟;相反,PTD值越小,则水电站i越倾向于留在大联盟中。若各水电站的PTD值均较接近且较小,表明该分摊方案下各个电站的谈判能力接近,且都倾向于留在大联盟中,即该方案为稳定性较好、较为接受的分摊方案。
表3不同分摊方法下各水电站分裂倾向值表
Figure BDA0002542474670000101
(3)构造评价指标值矩阵,选择以分摊比例、不同联盟方式下相对于单独调度时梯级增加的发电量、最小出力、PTD值为评价指标,4种方案的评价指标值如表4所示。
表4评价指标值表
Figure BDA0002542474670000102
(4)通过专家打分法和熵权法确定主客观综合权重,结果如表5所示。
表5各评价指标值主客观综合权重表
Figure BDA0002542474670000103
(5)通过后悔理论进行方案评价,其中后悔值矩阵为:
Figure BDA0002542474670000104
后悔-欣喜函数矩阵为:
Figure BDA0002542474670000111
感知效用矩阵为:
Figure BDA0002542474670000112
各方案的综合感知效用值以及方案评价结果如表6所示。
表6各方案综合感知效用值表
Figure BDA0002542474670000113
通过上述评价过程,得出如下结论:
①三里坪水电站作为龙头电站获得了最高的补偿,通过变异系数-Shapley值法获得补偿相对于其他方法的分摊结果均有不同程度的增加,寺坪、过渡湾水电站也获得比单独运行时更高的补偿效益且二者的分摊差异较大。这充分肯定了龙头电站在联盟中的作用与价值,更具公平性,在促进下游各级电站积极参与梯级联合调度的同时也凸显了次级电站的个体差异。
②通过PTD值比较分析发现,变异系数-Shapley值法是最为稳定较为大家接受的分摊方法,其中三里坪、寺坪水电站的PTD值均为最小,充分肯定了其作为上游电站的补偿调节能力,过渡湾水电站PTD值较大,主要原因是分摊中充分考虑了各水电站的贡献率,对贡献率较大牺牲较多的上游电站进行了较大的补偿。
③4种评价方案中,方案M1的综合感知效用值为0.822,方案M2的综合感知效用值为0.832,方案M3的综合感知效用值为0.877,方案M4的综合感知效用值为0.948,方案M4>方案M3>方案M2>方案M1,>表示优于,方案M4最优,方案M1最劣。

Claims (8)

1.基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,包括以下依次执行的步骤,
步骤1:建立水电站发电优化调度模型,比较不同联盟方式下各水电站发电量;
步骤2:利用不同的补偿分摊方法,进行梯级水电站补偿分摊;
步骤3:对步骤2中的补偿分摊结果进行稳定性评价;
步骤4:建立指标体系,构造评价指标值矩阵,制定梯级水电站补偿分摊方案;
步骤5:基于后悔理论,根据步骤4的各补偿分摊方案的评价指标值,计算得出最优的补偿分摊方案。
2.根据权利要求1所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤2中,所述不同的补偿分摊方法,包括贡献度法、变异系数法、Shapley值法、变异系数-Shapley值法。
3.根据权利要求1所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤3中,所述对步骤2中的补偿分摊结果进行稳定性评价,采用分裂倾向法。
4.根据权利要求1所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤4还包括运用专家评分法和熵权法计算指标体系的各个指标的权重。
5.根据权利要求1所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤4中,所述指标体系包括分摊比例、梯级增加发电量、最小出力、分裂倾向值。
6.根据权利要求2所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤1包括以下子步骤,
步骤1.1:以单一水电站发电效益最大或者考虑不同联盟方式下发电效益最大为目标函数,分别建立发电优化调度模型:
1)单一水电站年发电效益最大:
Figure FDA0002542474660000011
2)不同联盟方式下年发电效益最大:
Figure FDA0002542474660000012
其中B为调度期内联盟水电站总发电效益;K为联盟水电站总数;i为水电站编号,表示梯级中第i个水电站;t、T分别为调度期各时段编号及总时段数;Ai为水电站i的综合出力系数;Qi,t为水电站i在时段t的发电引用流量;Hi,t为水电站i在时段t的发电水头;Δt为时段长;pi为水电站i的上网电价;
步骤1.2:确定约束方程,主要包括水量平衡方程、水电站蓄水量约束、水库水位约束、水电站出力约束、水库下泄流量约束以及上下游水电站间流量联系;
步骤1.3:采用动态规划方法进行发电调度优化计算,得到不同联盟方式下各水电站的发电量及发电效益。
7.根据权利要求4所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,步骤5包括以下子步骤,
步骤5.1:构造补偿分摊方案的评价指标值矩阵并对其进行标准化处理;
步骤5.2:针对每项指标,取标准化后矩阵中指标评价值的最优值,建立理想方案;
步骤5.3:构建后悔值矩阵;
步骤5.4:构建感知效用矩阵;
步骤5.5:结合各指标的权重,计算各补偿分摊方案的综合感知效用值;
步骤5.6:比较各补偿分摊方案的综合感知效用值,选出最优的补偿分摊方案。
8.根据权利要求6所述的基于稳定性评价的水电站补偿分摊方案评价方法,其特征在于,所述贡献度法是指在梯级水电站联合优化调度的发电收益中,电站单位水能转化为电能进而转化为发电效益的效率越高,对联合调度发电收益的贡献度越大;
贡献度法的计算式如下:
Figure FDA0002542474660000021
其中xi为第i个水电站的分摊比例;n为梯级水电站总数;piki为第i个水电站势能转化为发电收益的转化系数;ki发电效率系数,代表势能转化为电能的系数;pi为上网电价,代表电能转化为发电效益的系数;
所述变异系数法根据水电站各个体特征参数在所有被评价的个体特征参数上观测值的变异程度大小来对其赋权,假设有n个水电站,每个水电站有m个个体特征参数,从而形成了一个参数值矩阵X=(Xij)n×m,计算过程如下:
①先对各项参数值进行归一化,再对参数值矩阵进行标准化处理,为使每个参数值均大于0,且在[0.6,1]范围内,对于正向参数有:
Figure FDA0002542474660000022
其中,Yij为第i个水电站的特征参数j的标准值,Yij∈[0.6,1];Xmax j、Xmin j分别为特征参数j的最大、最小特征值;
②计算第j个特征参数值的平均值和均方差:
Figure FDA0002542474660000031
Figure FDA0002542474660000032
③计算第j个特征参数值的变异系数:
Figure FDA0002542474660000033
④计算第j个特征参数权重值:
Figure FDA0002542474660000034
⑤计算第i个水电站的第j个特征参数值的比重:
Figure FDA0002542474660000035
⑥计算第i个水电站综合权重值:
Figure FDA0002542474660000036
其中,υi为第i个水电站的个体特征权重系数,即分摊比例;
所述Shapley值法按照参与人员边际费用的期望值进行效益分摊,以保证分摊的有效性与公平性,当满足有对称公理、有效公理、可加公理后,各水电站的Shapley值具有唯一解:
Figure FDA0002542474660000037
其中s为联盟中水电站个数;N为所有可能的联盟方式;S为N的一个子联盟;(s-1)!(n-s)!/n!为联盟S出现的联盟概率;ν(S)-ν(S-i)为电站i对联盟S的边际贡献值;
结合变异系数法和Shapley值法各自优势,提出变异系数-Shapley值法,根据式(4)~(10)得到各水电站个体特征权重系数υi,根据公式(11)得到各水电站的Shapley值σi,变异系数-Shapley值法的计算式如下:
θi=υiσi (12)
Figure FDA0002542474660000041
其中xi为第i个水电站的分摊比例;υi为第i个水电站的个体特征权重系数;σi为第i个水电站的的Shapley值。
CN202010550888.1A 2020-06-16 2020-06-16 基于稳定性评价的水电站补偿分摊方案评价方法 Active CN111814314B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010550888.1A CN111814314B (zh) 2020-06-16 2020-06-16 基于稳定性评价的水电站补偿分摊方案评价方法
CN202310964511.4A CN117081080A (zh) 2020-06-16 2020-06-16 以发电效益最大为目标的水电站补偿分摊方案评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010550888.1A CN111814314B (zh) 2020-06-16 2020-06-16 基于稳定性评价的水电站补偿分摊方案评价方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202310964511.4A Division CN117081080A (zh) 2020-06-16 2020-06-16 以发电效益最大为目标的水电站补偿分摊方案评价方法

Publications (2)

Publication Number Publication Date
CN111814314A true CN111814314A (zh) 2020-10-23
CN111814314B CN111814314B (zh) 2023-08-29

Family

ID=72845191

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202010550888.1A Active CN111814314B (zh) 2020-06-16 2020-06-16 基于稳定性评价的水电站补偿分摊方案评价方法
CN202310964511.4A Pending CN117081080A (zh) 2020-06-16 2020-06-16 以发电效益最大为目标的水电站补偿分摊方案评价方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202310964511.4A Pending CN117081080A (zh) 2020-06-16 2020-06-16 以发电效益最大为目标的水电站补偿分摊方案评价方法

Country Status (1)

Country Link
CN (2) CN111814314B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114677064A (zh) * 2022-05-27 2022-06-28 长江水利委员会水文局 一种耦合最优性与稳定性的梯级水库调度决策支持方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105787269A (zh) * 2016-02-25 2016-07-20 三明学院 一种基于后悔理论的异质多属性变权决策方法
CN105931131A (zh) * 2016-04-18 2016-09-07 华电电力科学研究院 梯级水电站群联合调度效益分配的Shapley值法
CN107153975A (zh) * 2016-10-28 2017-09-12 大连理工大学 一种基于联盟博弈的多控制性水库梯级电站补偿效益分摊方法
CN108988398A (zh) * 2018-07-25 2018-12-11 武汉大学 一种级联水电站群多利益主体调峰效益分配方法
CN110032770A (zh) * 2019-03-18 2019-07-19 华中科技大学 抽水蓄能机组双机相继开机规律的多目标优选方法及系统
CN110110965A (zh) * 2019-04-04 2019-08-09 浙江大学宁波理工学院 基于合格值自适应后悔理论的方案选择决策方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105787269A (zh) * 2016-02-25 2016-07-20 三明学院 一种基于后悔理论的异质多属性变权决策方法
CN105931131A (zh) * 2016-04-18 2016-09-07 华电电力科学研究院 梯级水电站群联合调度效益分配的Shapley值法
CN107153975A (zh) * 2016-10-28 2017-09-12 大连理工大学 一种基于联盟博弈的多控制性水库梯级电站补偿效益分摊方法
CN108988398A (zh) * 2018-07-25 2018-12-11 武汉大学 一种级联水电站群多利益主体调峰效益分配方法
CN110032770A (zh) * 2019-03-18 2019-07-19 华中科技大学 抽水蓄能机组双机相继开机规律的多目标优选方法及系统
CN110110965A (zh) * 2019-04-04 2019-08-09 浙江大学宁波理工学院 基于合格值自适应后悔理论的方案选择决策方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周彬彬;申建建;程春田;: "基于联盟博弈的多控制性水库梯级电站补偿效益分摊方法", 水电能源科学, no. 11, pages 58 - 62 *
曹云慧;王丽萍;王春超;朱艳霞;李继伟;周婷;: "基于熵权Shapley值法的梯级水电站补偿效益分摊", 水电能源科学, no. 02, pages 91 - 94 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114677064A (zh) * 2022-05-27 2022-06-28 长江水利委员会水文局 一种耦合最优性与稳定性的梯级水库调度决策支持方法
CN114677064B (zh) * 2022-05-27 2022-09-30 长江水利委员会水文局 一种耦合最优性与稳定性的梯级水库调度决策支持方法

Also Published As

Publication number Publication date
CN111814314B (zh) 2023-08-29
CN117081080A (zh) 2023-11-17

Similar Documents

Publication Publication Date Title
CN108599373B (zh) 高比例可再生能源电力系统输配协调调度目标级联分析法
CN110728410A (zh) 一种计及需求响应灵活性和不确定性的负荷聚合商经济调度方法
CN105160149B (zh) 一种模拟调峰机组的需求响应调度评估体系构建方法
CN106803171A (zh) 考虑日前市场价和风电不确定性的月度合约市场购电方法
CN107507078A (zh) 一种基于讨价还价博弈的分布式能源收益分配策略
CN106056307A (zh) 一种基于综合赋权法的电能替代技术综合效益评估方法
Ali Development and Improvement of Renewable Energy Integrated with Energy Trading Schemes based on Advanced Optimization Approaches
CN115358519A (zh) 一种虚拟电厂优化调度方法及装置
CN115313361A (zh) 一种规模化可调资源的联合优化控制方法和装置
CN111008463A (zh) 一种考虑发电侧储能的容量配置优化方法、装置及设备
CN114612001A (zh) 集群电动汽车参与电网削峰的调控指令分解方法和系统
CN111814314A (zh) 基于稳定性评价的水电站补偿分摊方案评价方法
Cheng et al. An innovative profit allocation to distributed energy resources integrated into virtual power plant
CN107154625B (zh) 基于模糊贝叶斯学习的电动汽车放电电价谈判方法
CN115952897A (zh) 一种城市能源系统的优化配置方法
CN112651544B (zh) 一种增量配电多主体协调优化方法
CN115833255A (zh) 一种灵活性供需平衡的源荷协同优化调度方法及终端
CN115514014A (zh) 一种含高比例风电的新型电力系统灵活性资源供需博弈优化调度方法
CN115811094A (zh) 一种考虑碳成本传导的新型电力系统调度方法与系统
CN114493227A (zh) 一种基于勾股模糊-topsis的新能源调度后评估方法及介质
Guo et al. Virtual Power Plant Operation Based on Entropy-Weight-Improved Cloud Risk Assessment Model.
CN113255957A (zh) 综合服务站不确定因素的定量优化分析方法及系统
CN105335799A (zh) 一种基于混合多属性评价的错峰用电预案编制方法
CN114649820B (zh) 一种风光储场站多环节容量分配方法及系统
Zhang et al. Review of Revenue Allocation Method for Virtual Power Plants

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