CN109472398A - 基于相对鲁棒CVaR的电网投资项目组合优化方法 - Google Patents
基于相对鲁棒CVaR的电网投资项目组合优化方法 Download PDFInfo
- Publication number
- CN109472398A CN109472398A CN201811228672.2A CN201811228672A CN109472398A CN 109472398 A CN109472398 A CN 109472398A CN 201811228672 A CN201811228672 A CN 201811228672A CN 109472398 A CN109472398 A CN 109472398A
- Authority
- CN
- China
- Prior art keywords
- electric grid
- cvar
- grid investment
- investment
- projects combo
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000003064 k means clustering Methods 0.000 claims abstract description 6
- 230000009467 reduction Effects 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims abstract description 4
- 230000005611 electricity Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- 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/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种基于相对鲁棒CVaR的电网投资项目组合优化方法,包括以下几个步骤:(1)设置净现值率为不确定性随机变量,收集不确定型变量的预测结果,使用蒙特卡洛仿真法根据连续概率分布生成离散的随机样本,并使用K‑means聚类方法进行样本缩减。(2)建立电网投资组合相对鲁棒CVaR优化模型,根据不同项目关系、项目基本情况设置电网投资项目约束条件,保证模型的可拓展性。(3)电网投资项目组合决策,基于步骤(1)和步骤(2)的相关计算,得到最优项目组合和各电网投资项目相应CVaR风险度量,从而实现对电网投资项目组合决策的优化过程。本发明方法能有效减小投资风险,并减少决策的保守性,实现盈利与风险的平衡。
Description
技术领域
本发明属于优化技术领域,尤其涉及电网投资优化技术领域。
背景技术
传统的电网投资项目组合优化往往只是在给定期望投资回报率的前提下,求得满足投资能力、预期回报率等经济性约束和负荷需求、容载比等物理性约束,并使得电网投资利润或净现值等目标函数最优的电网投资项目组合,在新一轮电改的背景下,电网公司投资项目收益在外部需面对输配电量、不同监管周期输配电价和配售电业务竞争等的不确定性风险,在内部将面对管控趋严的资本性支出和成本性支出等的不确定性风险,落实精准投资是急需解决的突出难题。
现有研究中很少针对电网投资项目组合风险进行优化建模,在电网投资环境日益复杂的新形势下,该领域研究的重要性日益凸显。现有的基于风险度量的投资优化模型多应用于股票、债券等金融投资领域,并取得了较好的应用效果。条件风险价值(ConditionalValue at Risk,CVaR)在投资收益具有非对称分布性,尤其是厚尾的情况下,具有极大优越性,但传统模型的不稳定性,即参数的微小变化对最优解产生很大的影响,CVaR模型仍没有克服。鲁棒优化方法通过保证最坏情景下的结果最优化,能有效解决模型对参数敏感性过高、参数估计误差影响较大的问题,但往往使得优化结果过于保守。
在上述背景下,提出了一种基于相对鲁棒CVaR的电网投资项目组合优化方法,在无法得知何种随机变量概率分布更为精确的情况下,考虑最坏情景下相对基准组合的风险度量,提高电网投资的精准度。本发明提出的模型兼顾鲁棒性和随机性,满足次可加性、正齐次性和单调性,该模型因优良特性具有很好的实际意义。
发明内容
本发明提供了一种基于相对鲁棒CVaR的电网投资项目组合优化方法,包括以下几个步骤:
(1)设置净现值率为不确定性随机变量,收集不确定型变量的预测结果,使用蒙特卡洛仿真法根据连续概率分布生成离散的随机样本,并使用K-means聚类方法进行样本缩减。
(2)建立电网投资组合相对鲁棒CVaR优化模型,根据不同项目关系、项目基本情况设置电网投资项目约束条件,保证模型的可拓展性。
(3)电网投资项目组合决策,基于步骤(1)和步骤(2)的相关计算,得到最优项目组合和各电网投资项目相应CVaR风险度量,从而实现对电网投资项目组合决策的优化过程。
本发明方法为电网投资项目组合决策提供了优化模型,该模型不需要准确获知不确定性随机变量的分布特征,通过实现CVaR模型的相对鲁棒性,能有效减小投资风险,并减少决策的保守性,实现盈利与风险的平衡。
附图说明
图1是K-means方法的实现步骤图。
具体实施方式
为实现本发明的目的,下面结合具体实施方式,进一步阐述本发明。
针对电网投资收益存在不确定性和传统优化模型存在缺陷的问题,本发明提出了一种基于相对鲁棒CVaR的电网投资项目组合优化方法,通过蒙特卡洛仿真和K-means聚类方法进行随机样本的生成与削减,并构建了基于相对鲁棒CVaR的电网投资项目组合优化模型,本发明提出的内容与方法在电网投资收益不确定的情况下能够有效减小电网投资决策的风险,并相对减少投资决策的保守性,提高企业经济效益。
1、设置净现值率为不确定性随机变量,各电网投资项目的净现值率未知,对其预测时可分为离散性和连续型概率分布,未知何种分布最为精确。为减小计算的复杂性和保持模型的统一性,将连续型概率分布函数用离散样本来近似表示。使用蒙特卡洛仿真法根据连续概率分布生成离散的随机样本,并使用K-means聚类方法进行样本缩减。
样本生成与缩减方法如下:
本发明使用蒙特卡洛仿真法根据连续概率分布生成离散的随机样本,并使用K-means聚类方法进行样本缩减。K-means方法是一种比较典型的聚类划分方法,能有效地处理规模较大和高维的数据集合,能够高效分类大型数据集。该聚类方法对于给定G个向量和H个聚类数,按照误差平方和或准则函数的划分标准把数据分成N组,使同一类中的对象相似度尽可能高,而不同类之间的相似度尽可能低,按这样的方法所形成的每一组就称作一个聚类(或聚簇)。本文通过K-means方法寻找蒙特卡洛仿真生成的样本数据的H个聚类,使用每个聚类的最终聚类中心代表随机样本特征值,K-means方法的实现步骤如图1所示。
2、建立电网投资组合相对鲁棒CVaR优化模型,根据不同项目关系、项目基本情况设置电网投资项目约束条件。
电网投资组合相对鲁棒CVaR优化模型如下:
设f(X,Y)表示一个与决策向量X相关的损失函数,X∈Ω(Ω是Rn的一个子集,表示决策变量X的可行集),Y∈Rm是一个随机向量,表示不确定因素。假定随机向量Y的联合概率密度为p(Y),对于固定的X,f(X,Y)关于Y的不超过某一损失水平α的概率为:
对于任意X,ψ(X,α)为决策向量X下的损失累积分布函数,关于α非减和右连续。
αβ(X)=min{α∈R:ψ(X,α)≥β}
式中:αβ(X)和φβ(X)即为置信水平β∈(0,1)下的VaR和CVaR,在相同的置信水平下,CVaR≥VaR,这样最小化CVaR也就是最小化VaR,由于φβ(X)中含有αβ(X)难以直接求解,用Fβ(X,Y)来替代φβ(X),则:
[f(X,Y)-α]+=max[f(X,Y)-α,0]
式(4)中存在积分计算,计算较为复杂,为简化计算,可将积分项进行离散化处理求得期望值,简化后如式(6)所示:
式中,s=1,2,…,S代表某一概率分布下离散样本的个数。
定义相对鲁棒CVaR风险度量,公式如下:
其中且CVaRβ(z*(p),p)>0,相对鲁棒CVaR表示的实际意义为:设X为投资组合的可行集,p是随机变量y可能的概率分布,则CVaRβ(z*(p),p)表示在X上的投资组合关于分布p的极小CVaR。
以上分析可知,公式(8)表示投资组合x相对于基准投资组合z*(p)和基准组合相应的概率分布p的相对风险,而实际应用中,分布p不可能准确知道,我们只能确定p的置信域P。
因此,公式(9)表示在可能概率分布的置信域P中,相对风险最大时的随机变量概率分布。相对鲁棒CVaR相对于绝对鲁棒CVaR风险度量降低了模型的鲁棒性,并满足正齐次性、可加性和单调性。同时,相对于传统的CVaR模型,其不依赖于精确的随机变量概率分布,在对CVaR风险度量影响最坏的随机变量概率分布情况下,求得最优的投资组合,兼具鲁棒优化模型的优点,因此,其具有更好的实际应用意义。
对于置信域P的选择,考虑到电网投资的现实情况,假设随机变量项目投资回报率y的分布有有限种可能,对随机变量分布的估计可以用专家估计或经验估计等方法,但我们无法得知哪种估计值更为精确,即设这里pi表示第i种可能的概率分布,l代表可能存在的概率分布的数量。在PM中相对鲁棒CVaR可以表示为:
所以,在相对鲁棒CVaR风险度量下电网投资组合问题可以表示为:
为使上述模型更易于计算,引入人工变量θ,θ无实际含义,公式(11)可以等价为:
再次,引入人工变量根据公式(4),将上述简化后的模型变为:
这里需要指出,当γ*(pi)=1时,相对鲁棒CVaR优化模型变为绝对鲁棒CVaR优化模型,在此基础上,当l=1时,模型退化为常规CVaR优化模型,至此相对鲁棒CVaR模型已转化成线性规划,只需计算线性规划序列,当损失函数定义为线性函数时,该模型为凸规划问题,可使用内点法等多种有效算法进行求解。
净现值是反映电网投资项目经济可行性的重要财务指标,本文中损失函数定义为投资项目净现值的相反数,且为保证该优化问题为凸规划,设置损失函数如下:
式中,x=(x1,x2,…,xn)T,n∈N,N表示可投资的电网项目数,n表示项目编号;x1,x2,…xn表示投资电网项目的0-1变量,0代表不投资该项目,1代表投资该项目;表示电网项目的净现值向量,代表第n个电网投资项目在第i种分布情况下第s个离散样本中的不确定性净现值率,%;Un代表第n个项目的初始投资额,万元。
在新电改的形势下,电网公司逐步向公共服务提供者的角色转变,电网投资项目约束条件除了要考虑常规的经济性约束外,还需考虑电网投资的社会性和可靠性等,同时,不同项目之间的关系也需考虑在内,因此,本发明设置的目标函数包括以下几部分:
停电成本约束,电网企业因其公共服务提供者的身份,为用户供应安全、可靠的电能是其最为重要的任务,停电等事故将会带来巨大的经济损失和社会影响。因此,本文用停电成本来衡量项目的社会性。假设停电成本和容量成正比,则可以得到如下约束:
式中,cn代表第n个项目带来的新增容量,MW;c代表电网原有容量,MW;fmax代表可接受最大停电损失,万元;fmean代表平均停电成本,万元。
投资能力约束,在电网投资项目决策之前,电网公司需对规划期内的公司投资能力进行测算,投资能力和公司的经营活动、融资活动和投资活动等有关,投资项目的总投资额不能超过公司的投资能力,投资能力约束如下:
式中,Umax代表电网企业最大投资能力,万元。
电力需求约束,负荷预测是电网投资的基础性工作,需要考虑地区的经济发展状况、电量水平和历史数据等方面,电网建设水平需满足该地区电力需求,约束如下:
式中,D代表该地区的电力需求,MW。
可靠性约束,可靠性约束用容载比来表示。容载比是反映电网供电能力的重要技术经济指标之一,当容载比过高时,增加电网运行成本,造成投资的效率降低;反之,若过低,会降低电网调整的灵活性,抑制电量的消耗。因此,容载比应按政策规定控制在一定的范围内,约束如下:
式中,γmin为容载比下限;γmax为容载比上限;L为整个电网系统主线平均负荷。
投资预期收益约束,选用投资期望净现值代表投资收益,约束如下:
式中,代表第i种分布的期望净现值向量,万元;μ表示电网企业投资预期最小贴现值,万元。
投资项目关系约束,若项目之间相互独立,则有:
x1+x2+…xn≤N
若项目为互斥关系,则
x1+x2+…xs≤1
若项目之间相互依存,也就是说,只有选择了项目n,项目n'才有可能被选上;反之,若项目n未被选上,则项目n'也不可能被选上,则有:
xn-xn'≥0
若项目之间为关系极为紧密,互相依存,即必须同时选中或未选中这两个项目,则有:
xn-xn'=0
基于以上内容,将电网投资组合的损失函数和相关约束条件带入相对鲁棒CVaR模型,得到基于相对鲁棒CVaR的电网投资项目组合优化模型,模型定义如下:
(3)收集各电网投资项目净现值相关预测结果,通过步骤(1)中的样本生成与缩减方法获得随机概率样本,将随机概率样本输入步骤(2)得到的电网投资组合相对鲁棒CVaR优化模型中,得到最优项目组合和各电网投资项目相应CVaR风险度量,从而实现对电网投资项目组合决策的优化过程。
以上为本发明的具体计算方式。以上显示和描述了本发明的基本原理和具体细节,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。
Claims (2)
1.一种基于相对鲁棒CVaR的电网投资项目组合优化方法,其特征在于,包括以下几个步骤:
(1)设置净现值率为不确定性随机变量,收集不确定型变量的预测结果,使用蒙特卡洛仿真法根据连续概率分布生成离散的随机样本,并使用K-means聚类方法进行样本缩减;
(2)建立电网投资组合相对鲁棒CVaR优化模型,根据不同项目关系、项目基本情况设置电网投资项目约束条件,保证模型的可拓展性;
(3)电网投资项目组合决策,基于步骤(1)和步骤(2)的相关计算,得到最优项目组合和各电网投资项目相应CVaR风险度量,从而实现对电网投资项目组合决策的优化过程。
2.根据权利要求1所述的方法,其特征在于,所述的使用K-means聚类方法进行样本缩减包括:对于给定G个向量和H个聚类数,按照误差平方和或准则函数的划分标准把数据分成N组,使同一类中的对象相似度尽可能高,而不同类之间的相似度尽可能低,按这样的方法所形成的每一组就称作一个聚类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811228672.2A CN109472398A (zh) | 2018-10-22 | 2018-10-22 | 基于相对鲁棒CVaR的电网投资项目组合优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811228672.2A CN109472398A (zh) | 2018-10-22 | 2018-10-22 | 基于相对鲁棒CVaR的电网投资项目组合优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109472398A true CN109472398A (zh) | 2019-03-15 |
Family
ID=65666013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811228672.2A Pending CN109472398A (zh) | 2018-10-22 | 2018-10-22 | 基于相对鲁棒CVaR的电网投资项目组合优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109472398A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111460378A (zh) * | 2020-03-20 | 2020-07-28 | 四川大学 | 一种考虑风险测度的配电网精准投资项目优选方法 |
-
2018
- 2018-10-22 CN CN201811228672.2A patent/CN109472398A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111460378A (zh) * | 2020-03-20 | 2020-07-28 | 四川大学 | 一种考虑风险测度的配电网精准投资项目优选方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Quan et al. | A survey of computational intelligence techniques for wind power uncertainty quantification in smart grids | |
Zeng et al. | Hybrid probabilistic-possibilistic approach for capacity credit evaluation of demand response considering both exogenous and endogenous uncertainties | |
Alamsyah et al. | Artificial neural network for Indonesian tourism demand forecasting | |
CN104636834B (zh) | 一种改进的联合概率规划模型系统优化方法 | |
CN109919356A (zh) | 一种基于bp神经网络区间需水预测方法 | |
CN108960485A (zh) | 一种源-荷互动电力市场下的在线字典学习概率最优潮流方法 | |
CN108921425A (zh) | 一种资产项目投资分类的方法、系统及服务器 | |
Kasule et al. | Forecasting Uganda’s net electricity consumption using a hybrid PSO-ABC Algorithm | |
Saini et al. | Data driven net load uncertainty quantification for cloud energy storage management in residential microgrid | |
Rayati et al. | Stochastic optimization and Markov chain-based scenario generation for exploiting the underlying flexibilities of an active distribution network | |
Souhe et al. | Forecasting of electrical energy consumption of households in a smart grid | |
Song et al. | Risky multi-criteria group decision making on green capacity investment projects based on supply chain | |
Nammouchi et al. | Robust opportunistic optimal energy management of a mixed microgrid under asymmetrical uncertainties | |
CN117744501A (zh) | 考虑生态流量的水网系统调蓄节点优化调度与决策方法 | |
Wang et al. | A study of situation awareness-based resource management scheme in cloud environment | |
CN109472398A (zh) | 基于相对鲁棒CVaR的电网投资项目组合优化方法 | |
Xuesen et al. | Consideration of trends in evaluating inter-basin water transfer alternatives within a fuzzy decision making framework | |
Sang et al. | Privacy-preserving hybrid cloud framework for real-time TCL-based demand response | |
Dolanyi et al. | Risk-based constraints for the optimal operation of an energy community | |
Li et al. | Implicit stochastic optimization with data mining for reservoir system operation | |
Wen et al. | Short-term load forecasting with bidirectional LSTM-attention based on the sparrow search optimisation algorithm | |
Shen et al. | An interval analysis scheme based on empirical error and MCMC to quantify uncertainty of wind speed | |
Crespo | Developing a battery management system for self-consumption systems | |
Dai | A distributed and real-time machine learning framework for smart meter big data | |
CN117853154A (zh) | 一种基于多特征的电力交易规划方法、系统、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190315 |
|
WD01 | Invention patent application deemed withdrawn after publication |