CN107316125A - 一种基于经济运行域的主动配电网经济运行评价方法 - Google Patents
一种基于经济运行域的主动配电网经济运行评价方法 Download PDFInfo
- Publication number
- CN107316125A CN107316125A CN201710336027.1A CN201710336027A CN107316125A CN 107316125 A CN107316125 A CN 107316125A CN 201710336027 A CN201710336027 A CN 201710336027A CN 107316125 A CN107316125 A CN 107316125A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mfrac
- power
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 212
- 238000011156 evaluation Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 90
- 238000005303 weighing Methods 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000013459 approach Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000005611 electricity Effects 0.000 claims description 47
- 238000004364 calculation method Methods 0.000 claims description 43
- 230000035699 permeability Effects 0.000 claims description 24
- 238000004146 energy storage Methods 0.000 claims description 21
- 230000001172 regenerating effect Effects 0.000 claims description 20
- 239000007789 gas Substances 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 17
- 238000007726 management method Methods 0.000 claims description 16
- 239000000446 fuel Substances 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 11
- 238000004891 communication Methods 0.000 claims description 11
- 238000012544 monitoring process Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 238000000205 computational method Methods 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 9
- 230000008901 benefit Effects 0.000 claims description 8
- 230000033228 biological regulation Effects 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 8
- 230000011218 segmentation Effects 0.000 claims description 8
- 230000006872 improvement Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 4
- 238000012423 maintenance Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000036642 wellbeing Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 239000003245 coal Substances 0.000 claims description 2
- 238000010835 comparative analysis Methods 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000010438 heat treatment Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims description 2
- 239000005364 simax Substances 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 230000005619 thermoelectricity Effects 0.000 claims description 2
- 239000002918 waste heat Substances 0.000 claims description 2
- 239000000567 combustion gas Substances 0.000 claims 2
- 230000005284 excitation Effects 0.000 claims 2
- 235000008429 bread Nutrition 0.000 claims 1
- 230000001276 controlling effect Effects 0.000 claims 1
- 230000002452 interceptive effect Effects 0.000 claims 1
- 230000009897 systematic effect Effects 0.000 claims 1
- 238000012797 qualification Methods 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 230000001131 transforming effect Effects 0.000 description 5
- 239000004020 conductor Substances 0.000 description 4
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 3
- 240000002853 Nelumbo nucifera Species 0.000 description 3
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 3
- 238000005265 energy consumption Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] 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
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明公开了一种基于经济运行域的主动配电网经济运行评价方法,其内容包括:在面向能源互联网的主动配电网的背景下定义14个二级评价指标,并与已有指标相结合,依据新的二级指标的特性分类处理确立4个新的一级指标,构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;利用熵权法修正,运用主客观赋权法相结合的方法改进区间层次分析法,并作为本发明的权重计算方法,利用配电网安全分析中“域”的方法论定义经济运行函数和经济运行域;基于经济运行域建立主动配电网各经济运行等级模型;根据经济运行函数是否越限判断待评样本是否属于对应等级。
Description
技术领域
本发明涉及主动配电网经济运行评价方法,尤其涉及一种基于经济运行域的主动配电网经济运行评价方法。
背景技术
配电网的经济运行是指通过调整不合理的网络结构,有目标的组织调整供用电负荷,使线路及设备在运行时间内,所输送的负荷尽量达到合理值,减少电能损耗,降低网损率,提高经济效益。
近年来,随着大规模的电动汽车(Electric Vehicle,EV)、储能系统(EnergyStorage System,ESS)、分布式可再生能源发电(Distributed Renewable EnergyGeneration,DREG)和冷热气电联产技术引入配电网,传统配电网在运行充裕性、安全性和经济性等方面面临较大挑战,因此主动配电网(Active Distribution Networks,ADN)成为解决这一问题的重要方案。主动配电网有更多的调控手段和能源途径促进可再生能源的消纳,对其经济运行进行评价不仅是节能降损及指导改造和规划的重要保证,也对协调管理可再生能源,提高能源利用率起到很好的指导作用,因此对其经济运行评价的研究具有非常重要的理论价值和实际意义。
目前已有多种评价方法被广泛应用于电网运行评价中,主要有层次分析法、灰色关联度法、BP神经网络法、模糊综合评价法、TOPSIS法和数据包络法等,近期的评价方法多是上述几种方法的组合或者改进,使得计算过程过于繁杂,且计算量较大,影响计算速度,使其在实际应用中受到限制。
目前国内外鲜有关于主动配电网经济运行评价的研究,主要都集中在电源规划、无功优化等方面。
综上所述,有必要发明一种主动配电网经济运行的评价方法,以解决现有主动配电网经济运行评价上存在的诸多问题。
发明内容
本发明目的在于提供一种精确性高、考虑因素全面、实用性强的基于经济运行域的主动配电网经济运行评价方法。
为实现上述目的,本发明是通过以下技术方案实现的:
一种基于经济运行域的主动配电网经济运行评价方法,该方法内容包括如下步骤:
步骤一在面向能源互联网的主动配电网的背景下定义14个二级评价指标;
步骤二将步骤一中定义的指标与已有指标相结合,并且依据定义的新的二级指标的特性分类处理,确立4个新的一级指标,构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;
步骤三利用熵权法修正,运用主客观赋权法相结合的方法改进区间层次分析法,并作为本发明的权重计算方法,利用配电网安全分析中“域”的方法论定义经济运行函数和经济运行域;
步骤四基于经济运行域建立主动配电网各经济运行等级模型;
步骤五根据经济运行函数是否越限判断待评样本是否属于对应等级。
由于采用上述技术方案,与现有技术相比本发明提供的一种基于经济运行域的主动配电网经济运行评价方法具有这样的有益效果:
1、本发明方法在面向能源互联网的大规模DREG和其他能源系统与主动配电网互联背景下定义新指标,并与原有指标相结合构建了主动配电网评价指标体系,并对部分已有指标进行了改进,使评价因素考虑更全面,应用背景更新,更接近实际情况;
2、本发明方法根据配电网安全性分析方法,利用“域”的方法论定义了经济运行函数和经济运行域,使得计算过程更加简单、快速,冗余性更小;
3、本发明方法充分考虑了大规模可再生能源和多种能源系统接入主动配电网的影响,并且在新背景下对部分已有指标进行了修正,将评价等级与经济运行域的超平面边界相结合,将抽象模糊的经济运行等级划分问题转化为判断经济运行函数是否越限,更加直观,更易于判断,实用性更强。对面向能源互联网的主动配电网经济运行进行评价不仅是节能降损及指导改造和规划的重要保证,也对协调管理多能源系统起到很好的指导作用。
附图说明
图1是本发明方法的面向能源互联网的主动配电网经济运行评价指标体系示意图;
图2是本发明方法中的经济运行域超平面二维断面示意图;
图3是本发明方法的总体流程图。
具体实施方式
下面结合附图对本发明做进一步说明:
本发明方法是在考虑大规模分布式电源接入主动配电网情况下提出的。
本发明的一种基于经济运行域的主动配电网经济运行评价方法,其总体流程图如图3所示,该方法的具体步骤如下:
(1)在面向能源互联网的主动配电网的背景下定义14个二级评价指标;参见图1所示的本发明方法的指标体系示意图,所述定义14个二级评价指标的具体内容如下:
(1-1)有功缓冲度
有功缓冲度定义为分布式可再生能源发电(Distributed Renewable EnergyGeneration,缩写为DREG。以下简称DREG)和多种能源接入系统前各个线路有功功率与接入后有功功率之差和接入前时有功功率之比的平均值,衡量DREG和多种能源系统接入后对线路输送有功功率的缓解程度,其计算公式为:
式中Pi、Pi DG分别为DREG和多种能源接入前、后时第i条线路有功功率,单位:MW;n为主动配电网线路的总数;
(1-2)潮流均衡度
潮流均衡度定义为线路的实际最大功率和极限容量之比的方差,大规模DREG的引入后系统潮流造成重要影响,潮流均衡度反映电网线路潮流分布是否均衡,其计算公式为:
式中, Simax分别为引入DREG后线路i的实际最大功率和极限容量,单位:MW;n为主动配电网线路的总数;
(1-3)可控负荷率
可控负荷率定义为实施直接负荷控制和可中断负荷管理控制的用电量占全部用电量的比例,对于风电、太阳能发电等可再生能源来说,通过负荷管理,可以有效改变负荷曲线形状,尤其是抬高低谷负荷,提高可再生能源消纳水平,其计算公式为:
式中,PCL为在负荷峰值时切断的可控负荷功率,单位:MW;Pav为日平均负荷,单位:MW;P(t)为t时刻的负荷,单位;MW;
(1-4)激励用电比例
激励用电比例定义为实施分时电价机制的用电量占全部用电量的比例,既包括普通负荷,也包括其他能源系统与主动配电网的能量交互,实施灵活电价机制的用电比例越高,用电侧与各种整体灵活性越好,有利于可再生能源的消纳,其计算公式为:
式中,PE为实施分时电价等灵活电价机制的用电量,单位:MW;
(1-5)智能监控通讯覆盖率
智能监控通讯覆盖率定义为拥有智能监测及通讯设备的配电变容量占全部配电变容量的比重,反映了主动配电网和通信系统的融合程度,也一定程度上反应了大数据采集技术的水平,主动配电网可以利用大数据海量处理技术进行负荷预测、需求侧管理等提高可靠性降低网损的措施,其计算公式为:
式中,Sk为第k台拥有智能监测及通讯设备的配电变容量,单位:MW;SDT为所有配电变总容量,单位:MW;
(1-6)用户参与度
用户参与度:用户参与度是指需求响应(Demand response,DR)项目中的响应负荷占总负荷的比例,用户不参与DR时,用电方式满意达最大值100%;用户参与DR后,调整用电方式、转移负荷,以用电满意度下降为代价,获得电费节约或经济补偿,其计算公式为:
式中,为参加DR的用户的负荷值,单位:MW;PL表示配网实际负荷值,单位:MW;
(1-7)DREG渗透率
DREG渗透率包含2种类型:分布式可再生能源静态渗透率λit和分布式能源有效渗透率λvx,其计算公式为:
式中,PDG为DREG的额定功率,MW;表示DREG的实际功率,单位:MW;表示配网最大负荷功率,单位:MW;基于分布式能源渗透率的定义,本发明方法DREG渗透率的定义如下:
DREG渗透率定义为静态渗透率λit和分布式能源有效渗透率λvx加权之和,其计算公式为:
式中,ωλ为权重,DREG渗透率反映了大规模DREG并网对主动配电网的影响程度;
(1-8)DREG消纳率
DREG消纳率定义为指主动配电网运行期间DREG的实际出力与最大允许出力之比,DREG消纳率反映了主动配电网对分布式能源的利用率,体现了主动配电网的兼容性及经济性,消纳率越高主动配电网越经济运行。
其计算公式为:
式中,表示实际环境下分布式能源的最大允许出力,单位:MW;
(1-9)调峰效应
调峰效应定义为负荷峰值时分布式发电、储能系统、冷热气电联产及可控负荷接入主动配电网的并网功率与负荷峰值之比,调峰效益能直观的反映出主动配电网协调管理分布式能源并网对调峰的影响,其计算公式为:
式中,为DREG在负荷峰值时的并网功率,单位:MW;为储能系统在负荷峰值时的放电功率,单位:MW;为冷热气电联产系统储热设备在负荷峰值时转带的热负荷功率,单位:MW;为电动汽车在负荷峰值时的放电功率,单位:MW;为冷热气电联供中燃气机组在负荷峰值时的发电功率,单位:MW;为在负荷峰值时切断的可控负荷功率;
(1-10)消纳贡献度
消纳贡献度定义为指电力系统以外其他能源系统消纳的DREG功率与DREG实际功率之比,多种能源系统接入主动配电网为其中大规模、高渗透率、随机性间歇性强的可再生能源提供更多的消纳选择空间。其计算公式为:
式中,为DREG工作时电动汽车充电对DREG的消纳量,单位:MW;为DREG工作时冷热气电热负荷对DREG的消纳量,单位:MW;冷热气电储热装置对DREG的消纳量,单位:MW;
(1-11)备用贡献度
备用贡献度定义为分布式可再生能源出力的备用容量与最大负荷之比,DREG的间歇性,用户侧消费负荷使用不确定性,对主动配电网备用容量提出了较高的要求,多种能源系统接入使得主动配电网的充裕性有了很大的改善。其计算公式为:
式中,为电动汽车放电的最大出力,单位:MW;为储能设备放电的最大功率,单位:MW;为冷热气电联供中燃气机组的最大出力,单位:MW;为DREG的最大出力,单位:MW;
(1-12)降损贡献度
多能源系统与主动配电网互联、DREG和ESS接入使得潮流流向由单向变为双向,在主动管理下系统潮流得到改善,从而引起电网损耗发生变化。本发明方法给出降损贡献度的定义如下:
降损贡献度定义为多能源系统与主动配电网互联引起的电网损耗变化量与并网前电网总损耗之比,其计算公式为:
δ=ΔPloss/ΔP0
式中,ΔPloss为DREG和多能源系统互联引起的网损变化量,单位:MW;ΔP0为线路总损耗,单位:MW;ΔPloss可表示为:
ΔPloss=ΔPDG+ΔPESS+ΔPEV+ΔPGAS
式中,ΔPDG、ΔPESS、ΔPEV和ΔPGAS分别为分布式可再生能源、储能系统、电动汽车和燃气机组引起的网损变化量;因此,有
式中,r为单位长度配电线路电阻值,单位:Ω/km;RT为变压器阻抗,单位:Ω;PL和QL分别为负荷的有功和无功功率;U为线路的电压,单位:V;PDG和QDG分别为并网DREG的有功和无功功率;PGAS和QGAS分别为并网冷热气电系统燃气机组的有功和无功功率;和分别为储能系统的充电和放电功率,单位:MW;和分别为电动汽车的充电和放电功率,单位:MW;M1、M2、M3和M4分别是DREG、储能设备、电动汽车充电站和燃气机组与变电站之间的配电线路长度,单位:m;
(1-13)消纳收益率
在多能源系统和DREG联合消纳过程中能带来很多方面的收益,主要包括延缓配电网络投资收益、节省网损收益、削峰填谷收益、提高可靠性收益及其他收益。本发明方法给出消纳收益率的定义如下:
消纳收益率定义为系统在进行分布式可再生能源消纳过程中带来的经济收益与所需的运行维护成本之比,其计算公式为:
式中,Ccap为延缓配电网络投资收益,单位:万元;Cgl为节省网损收益,万元;Cpls为削峰填谷收益,单位:万元;Crel为提高可靠性收益,单位:万元;Cres为其他收益,单位:万元;Cope为分布式能源运行维护费用,单位:万元;
Ccap为延缓配电网络投资时间所产生的收益,即年减少电网扩建容量的收益,大规模DREG接入主动配电网能起到很好的削峰作用,从而避免了因短暂的高峰负荷而追加的配电设备扩建容量的费用,其计算公式为:
Ccap=cdPNmax
式中,cd为配电设备的单位功率造价,单位:万元/MW;PNmax为分布式能源年消纳的最大功率,单位:MW;
Cgl为节省网损收益,其计算公式为:
Cgl=eΔPlossTmax
式中,e为平均电价,单位:万元/(MW·h);ΔPloss是DREG和多种能源系统并网贡献网损的变化量,单位:MW;Tmax为最大负荷年利用时间,单位:h;
Cpls为削峰填谷收益,此部分收益主要由储能系统和冷热气电联产系统产生,其公式为:
式中,n为年投运次数;分别表示第i小时段储能充、放电功率,单位:MW;分别表示第i小时段储能充、放电功率,单位:MW;PGAS(i)表示第i小时段冷热气电联供机组输出功率,单位:MW;表示第i小时段冷热气电联供储热装置的储、供热功率,单位:MW;表示第i小时段燃气机组余热回收装置效率,在同一时间内储能系统和电动汽车只会保持充电或者放电的一种状态;ei为第i小时段分时电价,单位:万元/(MW·h);eheat为第i小时段供热价格,单位:万元/(MW·h);
Crel为提高可靠性的收益即减少停电的损失,当发生停电事故时,分布式能源可作为备用电源投入配网,成为独立供电源,提高供电可靠性,从而减少停电损失;本发明方法采用电价倍数法计算停电损失,其计算公式为:
式中,k为分布式能源的平均容量系数;t为年平均停电时间,h;b为电价倍数,b=产电比/平均电价;
Cres为其他收益,主要包括节约燃料成本收益和政府补贴收益,分布式能源相比较传统能源而言,在发电过程中使用燃料较少,有些类型的分布式能源(太阳能、风能)甚至不消耗燃料,具有良好的经济效益,此外《电力需求侧管理城市综合试点工作中央财政资励资金管理暂行办法》中明确指出永久性节约和转移高峰电力负荷的补贴,东部地区补贴440元/kW,西部地区补贴550元/kW。
其计算公式为:
Cres=Cfuel+Cgon
式中,Cfuel为节约燃料成本收益,单位:万元;Cgon为政府补贴收益,单位:万元;Cgon单位功率的政府补贴资金,单位:万元/kW。T为各种清洁能源年利用小时数,h;α、β、γ分别表示燃料成本模型的系数;i表示各个清洁能源;
Cope为DREG和多种能源系统接入的运行维护费用,主要由其规模大小确定,其计算公式为:
Cope=cmPADN
式中,cm为单位功率的年运行维护成本,单位:万元/MW;PADN为DREG和多种能源系统的额定功率,单位:MW;
(1-14)环境收益率
环境收益率定义为主动配电网使用清洁能源生产的电能代替火力机组发电电能所带来的环境成本减少与所需的运行维护成本之比,其计算公式为:
式中,Ceco代替火力机组发电电能所减少环境成本,单位:万元。
(2)将步骤(1)中所定义的14个二级评价指标与已有二级评价指标相结合,并且依据定义的新的14个二级评价指标的特性分类处理,确立4个新的一级指标,构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;
所述构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;参见图1所示的本发明方法的指标体系示意图,其具体内容如下:
(2-1)构建面向能源互联网的主动配电网经济运行评价指标体系
在已有评价指标体系的基础上,针对面向能源互联网的主动配电网中多种能源形式参与可再生能源消纳的特点,加入步骤(1)中所定义的14个二级评价指标,并依据这14个二级评价指标的特性进行分类处理,确立4个新的一级评价指标,它们是系统灵活性、需求侧管理、多能源贡献和社会经济收益;
系统灵活度:可再生能源的接入为主动配电网甚至是整个电力系统的灵活度提升十分有利,但在目前高渗透率可再生能源接入的前提下,消纳率和调峰能力也成为配电网灵活度的有力补充,利用多种能源互联进行负荷转带能够显著提高配用电侧甚至整个能源互联网的灵活性,降低可再生能源利用率低造成的损失并增加经济收益;所述一级指标系统灵活度项下包括的二级指标为可控负荷率、激励用电比例、智能监控通讯覆盖率和用户参与度;
需求侧管理:主动配电网在管理模式上,基于智能通信平台,可实现对DSR的整合及对系统资产的分散式管理,实现电能品质优化控制,提高需求侧能源综合利用效率、多能源协同交互控制等目标。大规模DREG及DSR的接入后,用户侧管理,负荷控制,激励分时用电就成为了提高可再生能源消纳水平,增强配电系统运行经济性的重要手段;所述一级指标需求侧管理项下包括的二级指标为DREG渗透率、DREG消纳率和调峰效应;
多能源贡献:在多种能源系统接入面向能源互联网的主动配电网后,可再生能源的消纳需要考虑的不再仅仅是单一能源,即电能背景下的传统负荷就地消纳,而是要考虑多种能源系统形式参与到可再生能源的联合消纳过程中,对多能源系统在消纳过程、系统充裕性以及降低网损的量化考察成为体现运行经济性的重要手段;所述一级指标多能源贡献项下包括的二级指标为消纳贡献度、备用贡献度和降损贡献度;
社会经济收益:可再生能源接入配电网必然为减少火电煤耗,降低碳排放做出贡献,随着碳交易机制逐步引入电力市场,节能减排在决策过程中不再仅仅是一种单纯环境因素,而是可以带来明确的经济利益和社会效益的经济因素,其价值十分可观并将与日俱增;所述一级指标社会经济收益项下包括的二级指标为消纳收益率和环境收益率;
二级指标有功缓冲度和潮流均衡度依据其特性隶属于一级指标供电能力;
遵循系统性、科学性、客观性和实用性的原则,建立了包含一级评价指标和二级评价指标的递阶层次型主动配电网经济运行评价指标体系;
(2-2)完善传统配电网经济运行评价指标体系
本发明所述构建面向能源互联网的主动配电网经济运行评价指标体系是在已广泛应用的传统配电网经济运行评价指标体系基础上补充完善建立的;参见图1所示的本发明方法的面向能源互联网的主动配电网经济运行评价指标体系示意图,其中传统配电网经济运行评价指标包括如下内容:
供电能力,是衡量供电质量的主要一级指标之一,反映了供电部门向用户分配的电力是否合格。因此也是配电网建设改造评价的重要指标,以线路最大电压降落、母线电压合格率作为电压质量的二级指标:
线路最大电压降落ΔU,电压降落是指电力网中,首端相电压与末端相电压的差值。它反映该条线路末端用户的用电质量情况。线路最大电压降落直接反映了供电部门的供电稳定性,及用户用电负荷的安全性。其公式如下:
ΔU=U1-U2
式中,U1和U2分别为线路的首末端电压,单位:kV。
母线电压合格率ηB,是指母线电压合格时间与其运行时间的百分比。母线电压是用来衡量配电网经济运行情况的指标之一。待评价母线电压合格率公式如下:
ηB=tk/tΣ×100%
式中,tk为待评价母线电压合格时间;tΣ为待评价母线电压总运行时间。
电网损耗是衡量电网公司经济效益的一个重要一级指标。它在一定程度上能够反映电力网络运行耗费情况及其技术水平,因此也是电网经济技术评价的一个重要指标。以线路理论损耗率、配变理论损耗合格率作为电网损耗的二级指标:
线路理论损耗率ηL,线路理论损耗率是依据电网的结构参数(导线的规格型号、长度,设备的额定容量等)和运行参数(电压、电流、电量等)从理论计算中得出的线路损耗电量与线路供电量之比百分数。待评价线路理论损耗率公式如下:
式中,Ploss为待评价线路理论线损电量;Pg为待评价线路供电量。
配变理论损耗率合格率ηTL,配变理论损耗率ΔP是配电变压器损耗电量与输出电量之比的百分数。该指标反映了各条线路中变压器的损耗情况。ΔP的公式如下:
式中,ΔPg为待评价配电变压器的输入电量;ΔPs为待评价配电变压器的输出电量。待评价线路的配变理论损耗率合格率ηT公式如下:
ηTL=kP/kΣ×100%*
式中,kP为待评价线路中符合标准的配电变压器的个数;为待评价线路中配电变压器的总个数。
负荷是衡量线路安全经济运行的重要的一级指标之一。负荷不是恒定值,是随时间而变化的变动值,负荷状况能够反映电网的利用情况,从而反映电网运行的经济性。负荷情况会因为地区的不同而不同,伴随时间的变化负荷也存在高峰值与低谷值。因此考虑负荷情况应从多个角度去分析,在本发明中以线路功率因数、线路最佳负荷区域、线路负荷率作为负荷的二级指标:
线路负荷率cosφ,线路功率因数是指在交流配电线路中,线路电压与线路电流之间的相位差的余弦。线路功率因数反映整体线路的无功情况。待评价线路功率因数公式如下:
式中,Pl为待评价线路的有功功率;Sl为待评价线路的视在功率。
线路最佳负荷区域反映的是负荷的大小,根据负荷是否在线路的经济运行区域内来进行评价。线路总损耗ΔP0为固定损耗与可变损耗之和,公式如下:
式中,ΔPT为配电变压器的铁损之和;t为变压器运行时间;K为均方根负荷与平均负荷之比,取1.1;Req为配电网的等值电阻,包括线路与变压器两部分;P为计算期内平均负荷;U为配电网的额定电压;cosφ为功率因数。
配电网的网损率是指配网中损耗电量与供电电量的百分比,将上式两边同除以供电量P0=Pt可得网损率ΔP0%的公式如下:
将网损率对运行负荷求导,可知其一定存在极小值,经化简整理后可得网损率最低时所对应的运行负荷,称为配电网的经济运行负荷Pec,公式如下:
上式的意义是当配电网的运行负荷等于经济运行负荷Pec时,整个配电网的网损率出现极小值ΔP0min%,经济效果最佳。
由于实际运行负荷是一个变化的区域,而经济运行负荷只是区域中的一个点,很难做到恒定运行,因此,在将网损率控制在某一特定范围内的前提下,将配电网的运行负荷划定在一个区域内,此区域称为经济运行区域。
设配电网的网损率是最小网损率的α倍(α>1),即ΔP0%=αΔP0min%,则有:
化简整理可得网损率ΔP0%为最小网损率ΔP0min%的α倍时的运行负荷区域:
经济运行区域划分的合理与否,关键在于α的取值,通过对一些配电网的计算和分析,一般取1.2为宜。可将配电网的经济运行负荷区域[P]ec表示为[f1,f2],分别对应该区域的上下限:
线路负荷率Kf,线路负荷率是电网的平均负荷与最大负荷的比率。从经济运行角度考虑,负荷率越接近1,即峰谷差越小,电气设备的利用率越高,越有利于降损节能。因此,它是衡量用电均衡程度的指标,线路负荷是否合理直接影响了该线路的安全、合理、经济运行。待评价线路负荷率公式如下:
式中,Pl为待评价线路一定期间内的平均负荷Plmax为待评价线路一定期间内的最大负荷。
配网设备作为衡量配网经济运行技术的一级评价指标之一,配网设备的技术状况影响着整个配电网的技术水平。本发明以配变负载率合格率、线路导线截面积合格率、高耗能变压器占比作为配网设备的二级评价指标。
配变负荷率合格率ηT,配变负荷率是配电变压器在一定时间内的平均有功负荷与最高有功负荷的百分比,是反映发电、供电、用电设备是否充分利用的重要技术经济指标。配变负荷率公式如下:
式中,PT为待评价配电变压器一定期间内的平均负荷PTmax为待评价配电变压器一定期间内的最大负荷。给定KT>0.85为符合标准。待评价线路配变负荷率合格率公式如下:
式中,kh为待评价线路配电变压器符合标准的个数;k∑为待评价线路中所有配电变压器的个数和。
线路导线截面积合格率ηS,线路导线的截面积与线路的损耗情况有直接的关系。在输送负荷不变的情况下,使用截面积较粗的导线能够减少线路电阻,减少电能损耗和电压损失,达到降损节电效果,故以合格导线截面积所占比例来评价配电网线路的截面积选择是否经济。待评价线路导线截面积合格率公式如下:
式中,km为导线截面积合格的导线个数;kl为待评价线路中所有导线的个数和。
高耗能变压器占比ηH,在我国仍有SJ、SJL、SL等系列高损低效的配电变压器在电网中运行,不仅降低电网供电质量和运行的经济性,而且每年都造成巨额的电能损耗。为了减少电能的浪费,提高电网供电质量和运行经济性,应用低损高效的配电变压器替换高损耗变压器。待评价线路高耗能变压器占比公式如下:
式中,kg为待评价线路高耗能变压器个数;kΣ为待评价线路中所有配电变压器的个数和。
电网布局,随着电网负荷的不断变化,电网布局的合理性也成为影响配电网经济性的重要的一级指标。本发明以线路供电半径、运行方式、容载比作为电网布局的二级评价指标。
供电半径R是指以供电电源为中心,以供电距离为半径所画的圆为供电范围,这个半径即称为供电半径。低压供电半径指以配电变压器为中心的直线距离,它决定着配电变压器的合理布局,牵动配电网技术经济指标的优劣。供电半径的确定与负荷密度有很大关系,负荷密度越大,变电站容量就越大,供电半径就越小。《城市电力网规划设计导则》中规定,线路供电半径一般应满足:10kV线路的供电半径不大于15km。负荷密度小的地区,在保证电压质量和适度控制线损的前提下10kV线路供电半径可适当延长。
运行方式,改变运行方式就是通过改变网络结构,强行控制网络中的有功和无功流,实现负荷转移,达到降低网损的目的。本发明所涉及的运行方式都是在保证电网安全运行情况下的运行方式。由于运行方式的变化可以重新构造配电网络的运行结构,使负荷在各馈线之间相互转移而得到合理分配,进而降低损耗。因此,将运行方式的状态量化,以日平均负荷的变化ΔPav来表示运行方式的变化。
ΔPav=Pavmax-Pavmin
容载比Kc表示该线路、该站或该变压器的安装容量与最高实际运行容量的关系,反映容量备用情况。容载比过大,表示电网可供容量过多,反映电网建设早期投资增大;容载比过小,可供容量过小,又会使电网的适应性差,发生调度困难,甚至引起限电现象。要根据地区经济结构、发展速度、负荷增长情况以及企业的投资能力,来合理确定各规划期的容载比取值。线路容载比已成为确定该条线路容量的重要参数。待评价线路的容载比公式如下:
KC=∑S/ΣP
式中,ΣS为待评价线路所有配电变压器总容量;ΣP为待评价线路所有配电变压器的供电总负荷。
(2-3)对部分已有评价指标进行改进
对已有二级指标的计算方法增加在传统配电网中接入分布式能源后的变化量。现以二级指标线路最大电压降落为例介绍主动配电网接入分布式能源之后的计算方法,其它需要改进的二级指标计算方法与此类似。
传统配电网中没有接入分布式能源,须改进已有评价指标,使其反应主动配电网接入分布式能源后的运行特性,以下对其具体改进方法进行说明;
线路最大电压降落ΔU,改进计算方法中需要加入由于分布式能源引起的电压变化量,其计算公式如下:
ΔU=U1-U2+ΔUADN
式中,U1和U2分别为线路的首末端电压,单位:kV;ΔUADN为主动配电网加入分布式能源引起的电压变化量,单位:kV;
母线电压合格率ηB,改进计算方法中需要加入由于分布式能源引起的母线合格时间变化量,其计算公式如下:
式中,tk为待评价母线电压合格时间;tΣ为待评价母线电压总运行时间;ΔtADN为主动配电网加入分布式能源引起的电压变化量,单位:h;
线路理论损耗率ηL,改进计算方法中需要考虑由于接入分布式能源而引起的网络损失降低量,其计算公式如下:
式中,Ploss为待评价线路理论线损电量;Pg为待评价线路供电量;ΔPloss为DREG和多能源系统互联引起的网损变化量;
线路功率因数cosφ,改进计算方法中需要考虑由于接入分布式能源而引起的有功功率变化量,其计算公式如下:
式中,P为待评价线路的有功功率;S为待评价线路的视在功率,ΔPADN为分布式能源而引起的线路有功功率变化量。
线路负荷率Kf改进计算方法中需要考虑由于接入分布式能源而引起的线路负荷变化量。;其计算公式如下:
式中,Pload为待评价线路一定期间内的平均负荷;为DREG引起的负荷变化量。为待评价线路一定期间内的最高负荷。
(3)利用熵权法修正层次分析法,运用主客观赋权法相结合的方法改进区间层次分析法,并作为本发明的权重计算方法,利用配电网安全分析中“域”的方法论定义经济运行函数和经济运行域;
所述利用配电网安全分析中“域”的思想定义经济运行函数和经济运行域,具体定义内容如下:
(3-1)“域”的方法论是起源于输电网安全性研究的重要工具,相比于传统N-1逐点法,“域”能够描述运行点和安全边界的相对位置,提供系统整体的安全性裕度,有利于可视化监测和预防控制。近年来,一些学者尝试将“域”方法应用于配电网,提出了配电系统安全域(Distribution System Security Region,DSSR)的概念:DSSR是配电网中满足N-1安全约束的工作点集合,并建立了基于主变N-1的DSSR模型。
将“域”的方法论运用于主动配电网的经济运行评价中,在主动配电网中把所有满足相应程度的运行经济性的待评样本集合定义为配电网经济运行域;在此配电网安全分析方法的基础上结合超平面边界属性的特点划分经济等级;
(3-2)经济运行函数的定义
经济运行函数:在配电网经济运行评价中,把度量待评对象的经济运行程度函数称为经济运行函数,记为fE;
对于某一个一级指标,定义fE如下:
式中,xpq为第p个一级指标项下第q个二级指标,为第p个一级指标下第q个二级指标的第s个等级的越限值,r为第p个一级指标项下二级指标个数;经济运行函数的数值可以量化系统经济运行的程度:函数值越小,系统越经济;反之则越不经济;
第p个一级指标项下第q个二级指标的第s个等级的越限值的确定过程:根据不同二级指标的实际取值情况以及专家经验,为每个二级指标设置若干个准则值;为第p个一级指标项下第q个二级指标的第k个准则值;对第p个一级指标项下第q个二级指标的n个准则值求平均值,得:
大于均值的m个准则值再取均值;小于均值的l个准则值再取均值,得:
将区间按照一定的步长ξ均分为g(通常g=5)个子区间,每个区间对应一个等级,子区间边界值即为对应区间的越限值
(3-3)经济运行域的数学定义
经济运行域是在主动配电网中所有满足相应程度的运行经济性的待评样本集合,由经济运行函数的定义可得:
式中,为待评样本第s个等级的经济运行域;Θ为经济运行指标体系;yp为第p个一级指标,ωpq为第p个一级指标项下第q个二级指标的权重;
(3-4)权重计算方法
本方法中采用的权重计算方法为熵权法修正的改进区间层次分析法;区间层次分析法是以专家赋权法为计算基础的层次分析法;熵权法修正的区间层次分析法相对于数值层次分析法,充分考虑了专家对于指标比对主观因素判断的不确定性,能更好地改善层次分析法主观性;权重计算方法的具体步骤如下:
区间层次分析法是层次分析法在区间数上的扩展;
设为区间数矩阵,即记并记同样区间数向量即记并记对于给定的区间判断矩阵利用特征向量法分别求AL和AU的最大特征值所对应的具有正分量的归一化特征向量xL和xU;再由记α,β:
权重向量
本方法根据区间层次分析法得到各指标的权重区间,作为极大熵准则的约束条件,进行指标权重的非线性优化配置,其模型如下:
式中为区间层次分析法计算得到的第i个指标的权重区间。
熵权法确定客观权重方法如下:
指标归一化,将所有指标全部转化为相对应量数据,构建矩阵X=(aij)y×n,n为指标个数,y为待评样本个数。记xmax、xmin为同一指标最大和最小值,为平均值,得到B=(bij)y×n;
计算权重ωj
最终熵权法修正的改进区间层次分析法确定的权重为ω=Hjωi+(1-Hj)ωj,Hj为评价指标熵值。
(4)基于经济运行域建立主动配电网各经济运行等级模型,其具体建模过程如下:
(4-1)根据主动配电网中实测数据,计算各二级指标5个等级越限值;具体计算过程已在步骤(3)中已说明。
(4-2)由经济运行函数定义可知是r个形如的表达式的最大值,称每一个表达式为一个分段,则有r个分段:
fE在第n个分段具有初等函数形式,因此一定连续;对上式做恒等变形,可得
(4-3)fE所有分段的形式在指标体系的多维空间为一个超平面π(可由仿真工具拟合得出),一个超平面对应一个指标的等级准则;多个度量一级指标的二级指标形成的一组超平面包围的域即经济运行域,其超平面的二维断面示意图如图2所示。根据以上定义可得,基于经济运行域建立主动配电网各经济运行等级模型如下:
式中,为待评样本第s个等级的经济运行域;X为待评样本的指标向量,X=[x11x12 … x91 x92]T。
(5)根据经济运行函数是否越限判断待评样本是否属于对应等级;其具体判断过程如下:
(5-1)从数学上将经济运行等级划分问题转化为经济运行函数是否满足一个阈值的问题,从而判断待评样本或其中的某一一级指标是否处在某个等级被转化为经济运行函数与阈值比较;当fE=0时,待评样本(或其某一一级指标)处于该经济等级的边界上,故经济运行阈值取0:
fE≥0时,待评样本(或其某一级指标)处于该等级,
fE<0时,待评样本(或其某一级指标)不处于该等级;
(5-2)评价结果对比分析
运用模糊综合评价方法得到的评价等级与本发明方法所得结果一致,从而验证了本发明提出的评估方法的有效性和合理性。
以上所述的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (6)
1.一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:该方法内容包括如下步骤:
步骤一在面向能源互联网的主动配电网的背景下定义14个二级评价指标;
步骤二将步骤一中定义的指标与已有指标相结合,并且依据定义的新的二级指标的特性分类处理,确立4个新的一级指标,构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;
步骤三利用熵权法修正,运用主客观赋权法相结合的方法改进区间层次分析法,并作为本发明的权重计算方法,利用配电网安全分析中“域”的方法论定义经济运行函数和经济运行域;
步骤四基于经济运行域建立主动配电网各经济运行等级模型;
步骤五根据经济运行函数是否越限判断待评样本是否属于对应等级。
2.根据权利要求1所述的一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:在步骤一中,所述定义14个二级评价指标的具体内容如下:
(1-1)有功缓冲度
有功缓冲度定义为分布式可再生能源发电——DistributedRenewableEnergyGeneration,缩写为:DREG,和多种能源接入系统前各个线路有功功率与接入后有功功率之差和接入前时有功功率之比的平均值,其计算公式为:
<mrow>
<mi>d</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中Pi、Pi DG分别为DREG和多种能源接入前、后时第i条线路有功功率,单位:MW;n为主动配电网线路的总数;
(1-2)潮流均衡度
潮流均衡度定义为线路的实际最大功率和极限容量之比的方差,其计算公式为:
<mrow>
<mi>S</mi>
<mi>E</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
</mrow>
<mi>n</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
式中, Simax分别为引入DREG后线路i的实际最大功率和极限容量,单位:MW;n为主动配电网线路的总数;
(1-3)可控负荷率
可控负荷率定义为实施直接负荷控制和可中断负荷管理控制的用电量占全部用电量的比例,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>C</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>L</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mn>24</mn>
</msubsup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mn>24</mn>
</mfrac>
</mrow>
式中,PCL为在负荷峰值时切断的可控负荷功率,单位:MW;Pav为日平均负荷,单位:MW;P(t)为t时刻的负荷,单位;MW;
(1-4)激励用电比例
激励用电比例定义为实施分时电价机制的用电量占全部用电量的比例,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mi>E</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>P</mi>
<mi>E</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中,PE为实施分时电价等灵活电价机制的用电量,单位:MW;
(1-5)智能监控通讯覆盖率
智能监控通讯覆盖率定义为拥有智能监测及通讯设备的配电变容量占全部配电变容量的比重,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>S</mi>
<mi>k</mi>
</msub>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>D</mi>
<mi>T</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中,Sk为第k台拥有智能监测及通讯设备的配电变容量,单位:MW;SDT为所有配电变总容量,单位:MW;
(1-6)用户参与度
用户参与度:用户参与度是指需求响应(Demandresponse,DR)项目中的响应负荷占总负荷的比例,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中,为参加DR的用户的负荷值,单位:MW;PL表示配网实际负荷值,单位:MW;
(1-7)DREG渗透率
DREG渗透率包含2种类型:分布式可再生能源静态渗透率λit和分布式能源有效渗透率λvx,其计算公式为:
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>v</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
</mfrac>
</mrow>
式中,PDG为DREG的额定功率,MW;表示DREG的实际功率,单位:MW;表示配网最大负荷功率,单位:MW;基于分布式能源渗透率的定义,本发明方法DREG渗透率的定义如下:
DREG渗透率定义为静态渗透率λit和分布式能源有效渗透率λvx加权之和,其计算公式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mi>&lambda;</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>&lambda;</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>v</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mi>&lambda;</mi>
</msub>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mi>max</mi>
</msubsup>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>&lambda;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,ωλ为权重,DREG渗透率反映了大规模DREG并网对主动配电网的影响程度;
(1-8)DREG消纳率
DREG消纳率定义为指主动配电网运行期间DREG的实际出力与最大允许出力之比,其计算公式为:
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中,表示实际环境下分布式能源的最大允许出力,单位:MW;
(1-9)调峰效应
调峰效应定义为负荷峰值时分布式发电、储能系统、冷热气电联产及可控负荷接入主动配电网的并网功率与负荷峰值之比,其计算公式为:
<mrow>
<mi>&xi;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>^</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>^</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>^</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>^</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
<mo>^</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>L</mi>
</mrow>
<mo>^</mo>
</msubsup>
</mrow>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mfrac>
</mrow>
式中,为DREG在负荷峰值时的并网功率,单位:MW;为储能系统在负荷峰值时的放电功率,单位:MW;为冷热气电联产系统储热设备在负荷峰值时转带的热负荷功率,单位:MW;为电动汽车在负荷峰值时的放电功率,单位:MW;为冷热气电联供中燃气机组在负荷峰值时的发电功率,单位:MW;为在负荷峰值时切断的可控负荷功率;
(1-10)消纳贡献度
消纳贡献度定义为指电力系统以外其他能源系统消纳的DREG功率与DREG实际功率之比,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>E</mi>
<mi>A</mi>
<mi>T</mi>
</mrow>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msubsup>
</mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mfrac>
</mrow>
式中,为DREG工作时电动汽车充电对DREG的消纳量,单位:MW;为DREG工作时冷热气电热负荷对DREG的消纳量,单位:MW;冷热气电储热装置对DREG的消纳量,单位:MW;
(1-11)备用贡献度
备用贡献度定义为分布式可再生能源出力的备用容量与最大负荷之比,其计算公式为:
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mi>max</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mi>max</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
<mi>max</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mi>max</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mi>max</mi>
</msubsup>
</mrow>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中,为电动汽车放电的最大出力,单位:MW;为储能设备放电的最大功率,单位:MW;为冷热气电联供中燃气机组的最大出力,单位:MW;为DREG的最大出力,单位:MW;
(1-12)降损贡献度
降损贡献度定义为多能源系统与主动配电网互联引起的电网损耗变化量与并网前电网总损耗之比,其计算公式为:
δ=ΔPloss/ΔP0
式中,ΔPloss为DREG和多能源系统互联引起的网损变化量,单位:MW;ΔP0为线路总损耗,单位:MW;ΔPloss可表示为:
ΔPloss=ΔPDG+ΔPESS+ΔPEV+ΔPGAS
式中,ΔPDG、ΔPESS、ΔPEV和ΔPGAS分别为分布式可再生能源、储能系统、电动汽车和燃气机组引起的网损变化量;因此,有
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>rM</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>Q</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>rM</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>rM</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>rM</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>Q</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Q</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mi>L</mi>
<mo>+</mo>
<msub>
<mi>KR</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,r为单位长度配电线路电阻值,单位:Ω/km;RT为变压器阻抗,单位:Ω;PL和QL分别为负荷的有功和无功功率;U为线路的电压,单位:V;PDG和QDG分别为并网DREG的有功和无功功率;PGAS和QGAS分别为并网冷热气电系统燃气机组的有功和无功功率;和分别为储能系统的充电和放电功率,单位:MW;和分别为电动汽车的充电和放电功率,单位:MW;M1、M2、M3和M4分别是DREG、储能设备、电动汽车充电站和燃气机组与变电站之间的配电线路长度,单位:m;
(1-13)消纳收益率
消纳收益率定义为系统在进行分布式可再生能源消纳过程中带来的经济收益与所需的运行维护成本之比,其计算公式为:
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>l</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msub>
</mfrac>
</mrow>
式中,Ccap为延缓配电网络投资收益,单位:万元;Cgl为节省网损收益,万元;Cpls为削峰填谷收益,单位:万元;Crel为提高可靠性收益,单位:万元;Cres为其他收益,单位:万元;Cope为分布式能源运行维护费用,单位:万元;
Ccap为延缓配电网络投资时间所产生的收益,即年减少电网扩建容量的收益,其计算公式为:
Ccap=cdPNmax
式中,cd为配电设备的单位功率造价,单位:万元/MW;PNmax为分布式能源年消纳的最大功率,单位:MW;
Cgl为节省网损收益,其计算公式为:
Cgl=eΔPlossTmax
式中,e为平均电价,单位:万元/(MW·h);ΔPloss是DREG和多种能源系统并网贡献网损的变化量,单位:MW;Tmax为最大负荷年利用时间,单位:h;
Cpls为削峰填谷收益,此部分收益主要由储能系统和冷热气电联产系统产生,其公式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>l</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>V</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>P</mi>
</mrow>
</msubsup>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,n为年投运次数;分别表示第i小时段储能充、放电功率,单位:MW;分别表示第i小时段储能充、放电功率,单位:MW;PGAS(i)表示第i小时段冷热气电联供机组输出功率,单位:MW;表示第i小时段冷热气电联供储热装置的储、供热功率,单位:MW;表示第i小时段燃气机组余热回收装置效率,在同一时间内储能系统和电动汽车只会保持充电或者放电的一种状态;ei为第i小时段分时电价,单位:万元/(MW·h);eheat为第i小时段供热价格,单位:万元/(MW·h);
Crel为提高可靠性的收益即减少停电的损失,本发明方法采用电价倍数法计算停电损失,其计算公式为:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>kbetP</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mrow>
式中,k为分布式能源的平均容量系数;t为年平均停电时间,h;b为电价倍数,b=产电比/平均电价;
Cres为其他收益,主要包括节约燃料成本收益和政府补贴收益,其计算公式为:
Cres=Cfuel+Cgon
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
式中,Cfuel为节约燃料成本收益,单位:万元;Cgon为政府补贴收益,单位:万元;Cgon单位功率的政府补贴资金,单位:万元/kW。T为各种清洁能源年利用小时数,h;α、β、γ分别表示燃料成本模型的系数;i表示各个清洁能源;
Cope为DREG和多种能源系统接入的运行维护费用,主要由其规模大小确定,其计算公式为:
Cope=cmPADN
式中,cm为单位功率的年运行维护成本,单位:万元/MW;PADN为DREG和多种能源系统的额定功率,单位:MW;
(1-14)环境收益率
环境收益率定义为主动配电网使用清洁能源生产的电能代替火力机组发电电能所带来的环境成本减少与所需的运行维护成本之比,其计算公式为:
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>e</mi>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>e</mi>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msub>
</mfrac>
</mrow>
式中,Ceco代替火力机组发电电能所减少环境成本,单位:万元。
3.根据权利要求1所述的一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:在步骤二中,所述构建面向能源互联网的主动配电网经济运行评价指标体系,并对部分已有指标进行改进;其具体内容如下:
(2-1)构建面向能源互联网的主动配电网经济运行评价指标体系
在已有评价指标体系的基础上,针对面向能源互联网的主动配电网中多种能源形式参与可再生能源消纳的特点,加入步骤一中所定义的14个二级评价指标,并依据这14个二级评价指标的特性进行分类处理,确立4个新的一级评价指标,它们是系统灵活性、需求侧管理、多能源贡献和社会经济收益;
系统灵活度:可再生能源的接入为主动配电网甚至是整个电力系统的灵活度提升十分有利,但在目前高渗透率可再生能源接入的前提下,消纳率和调峰能力也成为配电网灵活度的有力补充,利用多种能源互联进行负荷转带能够显著提高配用电侧甚至整个能源互联网的灵活性,降低可再生能源利用率低造成的损失并增加经济收益;所述一级指标系统灵活度项下包括的二级指标为可控负荷率、激励用电比例、智能监控通讯覆盖率和用户参与度;
需求侧管理:主动配电网在管理模式上,基于智能通信平台,可实现对DSR的整合及对系统资产的分散式管理,实现电能品质优化控制,提高需求侧能源综合利用效率、多能源协同交互控制等目标。大规模DREG及DSR的接入后,用户侧管理,负荷控制,激励分时用电就成为了提高可再生能源消纳水平,增强配电系统运行经济性的重要手段;所述一级指标需求侧管理项下包括的二级指标为DREG渗透率、DREG消纳率和调峰效应;
多能源贡献:在多种能源系统接入面向能源互联网的主动配电网后,可再生能源的消纳需要考虑的不再仅仅是单一能源,即电能背景下的传统负荷就地消纳,而是要考虑多种能源系统形式参与到可再生能源的联合消纳过程中,对多能源系统在消纳过程、系统充裕性以及降低网损的量化考察成为体现运行经济性的重要手段;所述一级指标多能源贡献项下包括的二级指标为消纳贡献度、备用贡献度和降损贡献度;
社会经济收益:可再生能源接入配电网必然为减少火电煤耗,降低碳排放做出贡献,随着碳交易机制逐步引入电力市场,节能减排在决策过程中不再仅仅是一种单纯环境因素,而是可以带来明确的经济利益和社会效益的经济因素,其价值十分可观并将与日俱增;所述一级指标社会经济收益项下包括的二级指标为消纳收益率和环境收益率;
二级指标有功缓冲度和潮流均衡度依据其特性隶属于一级指标供电能力;
遵循系统性、科学性、客观性和实用性的原则,建立了包含一级评价指标和二级评价指标的递阶层次型主动配电网经济运行评价指标体系;
(2-2)完善传统配电网经济运行评价指标体系
本发明所述构建面向能源互联网的主动配电网经济运行评价指标体系是在已广泛应用的传统配电网经济运行评价指标体系基础上补充完善建立的;
(2-3)对部分已有评价指标进行改进
对已有二级指标的计算方法增加在传统配电网中接入分布式能源后的变化量。
4.根据权利要求1所述的一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:在步骤三中,所述利用配电网安全分析中“域”的思想定义经济运行函数和经济运行域,具体定义内容如下:
(3-1)将“域”的方法论运用于主动配电网的经济运行评价中,在主动配电网中把所有满足相应程度的运行经济性的待评样本集合定义为配电网经济运行域;在此配电网安全分析方法的基础上结合超平面边界属性的特点划分经济等级;
(3-2)经济运行函数的定义
经济运行函数:在配电网经济运行评价中,把度量待评对象的经济运行程度函数称为经济运行函数,记为fE;
对于某一个一级指标,定义fE如下:
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>E</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>r</mi>
</mrow>
</munder>
<mo>{</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>}</mo>
</mrow>
式中,xpq为第p个一级指标项下第q个二级指标,为第p个一级指标下第q个二级指标的第s个等级的越限值,r为第p个一级指标项下二级指标个数;经济运行函数的数值可以量化系统经济运行的程度:函数值越小,系统越经济;反之则越不经济;
(3-3)经济运行域的数学定义
经济运行域是在主动配电网中所有满足相应程度的运行经济性的待评样本集合,由经济运行函数的定义可得:
<mrow>
<msubsup>
<mi>&Omega;</mi>
<mi>E</mi>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mo>{</mo>
<mover>
<mi>y</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>&Element;</mo>
<mi>&Theta;</mi>
<mo>|</mo>
<mover>
<mi>y</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msubsup>
<mi>f</mi>
<mi>E</mi>
<mi>s</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mi>y</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</munderover>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
式中,为待评样本第s个等级的经济运行域;Θ为经济运行指标体系;yp为第p个一级指标,ωpq为第p个一级指标项下第q个二级指标的权重;
(3-4)权重计算方法
本方法中采用的权重计算方法为熵权法修正的改进区间层次分析法;区间层次分析法是以专家赋权法为计算基础的层次分析法;熵权法修正的区间层次分析法相对于数值层次分析法,充分考虑了专家对于指标比对主观因素判断的不确定性,能更好地改善层次分析法主观性;权重计算方法的具体步骤如下:
区间层次分析法是层次分析法在区间数上的扩展;
设为区间数矩阵,即记并记同样区间数向量即记并记对于给定的区间判断矩阵利用特征向量法分别求AL和AU的最大特征值所对应的具有正分量的归一化特征向量xL和xU;再由记α,β:
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>U</mi>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>,</mo>
<mi>&beta;</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>L</mi>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
权重向量
本方法将区间层次分析法得到各指标的权重区间,作为极大熵准则的约束条件,进行指标权重的非线性优化配置,其模型如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mi>F</mi>
<mo>=</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&omega;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>ln&omega;</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&omega;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<munder>
<mi>&omega;</mi>
<mo>&OverBar;</mo>
</munder>
<mi>i</mi>
</msub>
</mtd>
<mtd>
<msub>
<mover>
<mi>&omega;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中为区间层次分析法计算得到的第i个指标的权重区间。
5.根据权利要求1所述的一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:在步骤四中,所述基于经济运行域建立主动配电网各经济运行等级模型,其具体建模过程如下:
(4-1)根据主动配电网中实测数据,计算各二级指标5个等级越限值;
(4-2)由经济运行函数定义可知是r个形如的表达式的最大值,称每一个表达式为一个分段,则有r个分段:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>E</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
</mrow>
fE在第n个分段具有初等函数形式,因此一定连续;对上式做恒等变形,可得
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>&times;</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>E</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
(4-3)fE所有分段的形式在指标体系的多维空间为一个超平面π,可由仿真工具拟合得出,一个超平面对应一个指标的等级准则;多个度量一级指标的二级指标形成的一组超平面包围的域即经济运行域,根据以上定义可得,基于经济运行域建立主动配电网各经济运行等级模型如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>r</mi>
</mrow>
</munder>
<mo>{</mo>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>/</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msubsup>
<mo>}</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mi>X</mi>
<mo>&Element;</mo>
<msubsup>
<mi>&Omega;</mi>
<mi>E</mi>
<mi>s</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
8
式中,为待评样本第s个等级的经济运行域;X为待评样本的指标向量,X=[x11 x12... x91 x92]T。
6.根据权利要求1所述的一种基于经济运行域的主动配电网经济运行评价方法,其特征在于:在步骤五中,所述根据经济运行函数是否越限判断待评样本是否属于对应等级;其具体判断过程如下:
(5-1)从数学上将经济运行等级划分问题转化为经济运行函数是否满足一个阈值的问题,从而判断待评样本或其中的某一一级指标是否处在某个等级被转化为经济运行函数与阈值比较;当fE=0时,待评样本处于该经济等级的边界上,故经济运行阈值取0:
fE≥0时,待评样本处于该等级,
fE<0时,待评样本不处于该等级;
(5-2)评价结果对比分析
运用模糊综合评价方法得到的评价等级与本发明方法所得结果一致,从而验证了本发明提出的评估方法的有效性和合理性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710336027.1A CN107316125A (zh) | 2017-05-12 | 2017-05-12 | 一种基于经济运行域的主动配电网经济运行评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710336027.1A CN107316125A (zh) | 2017-05-12 | 2017-05-12 | 一种基于经济运行域的主动配电网经济运行评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107316125A true CN107316125A (zh) | 2017-11-03 |
Family
ID=60181568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710336027.1A Pending CN107316125A (zh) | 2017-05-12 | 2017-05-12 | 一种基于经济运行域的主动配电网经济运行评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107316125A (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108110790A (zh) * | 2017-12-15 | 2018-06-01 | 国网江苏省电力有限公司南京供电分公司 | 基于多源数据阶梯式动态模糊系统的主动配电网运行状态评估方法 |
CN108596476A (zh) * | 2018-04-23 | 2018-09-28 | 国网经济技术研究院有限公司 | 一种电源送出电网工程运行效益评价方法及系统 |
CN108649556A (zh) * | 2018-04-16 | 2018-10-12 | 东南大学 | 一种电网应急场景随机优化调度方法 |
CN108717597A (zh) * | 2018-04-23 | 2018-10-30 | 国网经济技术研究院有限公司 | 一种优化网架结构的电网工程运行效益评价方法及系统 |
CN108879659A (zh) * | 2018-06-22 | 2018-11-23 | 天津大学 | 一种弹性配电网灵活性的评价方法 |
CN109767080A (zh) * | 2018-12-23 | 2019-05-17 | 上海电力学院 | 一种社区综合能源系统的需求侧响应评估方法 |
CN109784573A (zh) * | 2019-01-24 | 2019-05-21 | 南方电网科学研究院有限责任公司 | 一种能源互联网多目标优化方法及装置 |
CN110046758A (zh) * | 2019-04-09 | 2019-07-23 | 湘潭大学 | 一种结合智能合约的微网用电调度方法 |
CN110417000A (zh) * | 2019-07-08 | 2019-11-05 | 国网天津市电力公司电力科学研究院 | 计及分布式电源和需求响应的配电安全域模型构建方法 |
CN110969347A (zh) * | 2019-11-22 | 2020-04-07 | 国网山东省电力公司临沂供电公司 | 一种输电网结构形态评估方法 |
CN110991715A (zh) * | 2019-11-22 | 2020-04-10 | 国网河北省电力有限公司电力科学研究院 | 一种基于改进灰色关联度的配电网能效评估体系和评估方法 |
CN111563691A (zh) * | 2020-05-18 | 2020-08-21 | 国网山西省电力公司经济技术研究院 | 一种接入新能源的交直流混合配电网性能评估方法 |
CN111724049A (zh) * | 2020-06-08 | 2020-09-29 | 国网河北省电力有限公司电力科学研究院 | 一种潜在电力能效服务客户的研判方法 |
CN113449241A (zh) * | 2020-03-25 | 2021-09-28 | 上海交通大学 | 含电力电子变压器的交直流混合配电网综合能效评估方法 |
CN113537708A (zh) * | 2021-06-09 | 2021-10-22 | 天津电力工程监理有限公司 | 基于组合赋权法的输变电工程造价管理成效综合评价方法 |
CN113595069A (zh) * | 2021-07-22 | 2021-11-02 | 江阴长仪集团有限公司 | 一种配电网光伏消纳容量提升方法 |
CN114943463A (zh) * | 2022-06-02 | 2022-08-26 | 西南石油大学 | 一种配电网节能技措评估方法 |
-
2017
- 2017-05-12 CN CN201710336027.1A patent/CN107316125A/zh active Pending
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108110790A (zh) * | 2017-12-15 | 2018-06-01 | 国网江苏省电力有限公司南京供电分公司 | 基于多源数据阶梯式动态模糊系统的主动配电网运行状态评估方法 |
CN108110790B (zh) * | 2017-12-15 | 2021-09-10 | 国网江苏省电力有限公司南京供电分公司 | 基于多源数据阶梯式的主动配电网运行状态评估方法 |
CN108649556B (zh) * | 2018-04-16 | 2021-06-01 | 东南大学 | 一种电网应急场景随机优化调度方法 |
CN108649556A (zh) * | 2018-04-16 | 2018-10-12 | 东南大学 | 一种电网应急场景随机优化调度方法 |
CN108596476A (zh) * | 2018-04-23 | 2018-09-28 | 国网经济技术研究院有限公司 | 一种电源送出电网工程运行效益评价方法及系统 |
CN108717597A (zh) * | 2018-04-23 | 2018-10-30 | 国网经济技术研究院有限公司 | 一种优化网架结构的电网工程运行效益评价方法及系统 |
CN108717597B (zh) * | 2018-04-23 | 2021-12-07 | 国网经济技术研究院有限公司 | 一种优化网架结构的电网工程运行效益评价方法及系统 |
CN108596476B (zh) * | 2018-04-23 | 2021-12-07 | 国网经济技术研究院有限公司 | 一种电源送出电网工程运行效益评价方法及系统 |
CN108879659A (zh) * | 2018-06-22 | 2018-11-23 | 天津大学 | 一种弹性配电网灵活性的评价方法 |
CN109767080A (zh) * | 2018-12-23 | 2019-05-17 | 上海电力学院 | 一种社区综合能源系统的需求侧响应评估方法 |
CN109784573A (zh) * | 2019-01-24 | 2019-05-21 | 南方电网科学研究院有限责任公司 | 一种能源互联网多目标优化方法及装置 |
CN109784573B (zh) * | 2019-01-24 | 2023-06-20 | 南方电网科学研究院有限责任公司 | 一种能源互联网多目标优化方法及装置 |
CN110046758A (zh) * | 2019-04-09 | 2019-07-23 | 湘潭大学 | 一种结合智能合约的微网用电调度方法 |
CN110417000A (zh) * | 2019-07-08 | 2019-11-05 | 国网天津市电力公司电力科学研究院 | 计及分布式电源和需求响应的配电安全域模型构建方法 |
CN110417000B (zh) * | 2019-07-08 | 2023-03-31 | 国网天津市电力公司电力科学研究院 | 计及分布式电源和需求响应的配电安全域模型构建方法 |
CN110969347A (zh) * | 2019-11-22 | 2020-04-07 | 国网山东省电力公司临沂供电公司 | 一种输电网结构形态评估方法 |
CN110969347B (zh) * | 2019-11-22 | 2023-08-18 | 国网山东省电力公司临沂供电公司 | 一种输电网结构形态评估方法 |
CN110991715A (zh) * | 2019-11-22 | 2020-04-10 | 国网河北省电力有限公司电力科学研究院 | 一种基于改进灰色关联度的配电网能效评估体系和评估方法 |
CN113449241A (zh) * | 2020-03-25 | 2021-09-28 | 上海交通大学 | 含电力电子变压器的交直流混合配电网综合能效评估方法 |
CN111563691A (zh) * | 2020-05-18 | 2020-08-21 | 国网山西省电力公司经济技术研究院 | 一种接入新能源的交直流混合配电网性能评估方法 |
CN111724049B (zh) * | 2020-06-08 | 2023-05-23 | 国网河北省电力有限公司电力科学研究院 | 一种潜在电力能效服务客户的研判方法 |
CN111724049A (zh) * | 2020-06-08 | 2020-09-29 | 国网河北省电力有限公司电力科学研究院 | 一种潜在电力能效服务客户的研判方法 |
CN113537708A (zh) * | 2021-06-09 | 2021-10-22 | 天津电力工程监理有限公司 | 基于组合赋权法的输变电工程造价管理成效综合评价方法 |
CN113595069A (zh) * | 2021-07-22 | 2021-11-02 | 江阴长仪集团有限公司 | 一种配电网光伏消纳容量提升方法 |
CN113595069B (zh) * | 2021-07-22 | 2023-10-20 | 江阴长仪集团有限公司 | 一种配电网光伏消纳容量提升方法 |
CN114943463A (zh) * | 2022-06-02 | 2022-08-26 | 西南石油大学 | 一种配电网节能技措评估方法 |
CN114943463B (zh) * | 2022-06-02 | 2024-04-26 | 西南石油大学 | 一种配电网节能技措评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107316125A (zh) | 一种基于经济运行域的主动配电网经济运行评价方法 | |
CN110866627B (zh) | 一种考虑阶梯气价的多区域电-气耦合综合能源系统优化调度方法 | |
CN105160451B (zh) | 一种含电动汽车的微电网多目标优化调度方法 | |
CN105591406B (zh) | 一种基于非合作博弈的微电网能量管理系统的优化算法 | |
CN106228258B (zh) | 一种计及需求侧管理的家庭能源局域网能量优化控制方法 | |
CN105071389B (zh) | 计及源网荷互动的交直流混合微电网优化运行方法及装置 | |
CN107301472A (zh) | 基于场景分析法和电压调节策略的分布式光伏规划方法 | |
CN102855592B (zh) | 多目标调峰优化方法和系统 | |
CN107437824A (zh) | 一种基于遗传算法的地区分布式电源准入容量的计算方法 | |
CN108471119B (zh) | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 | |
CN107769237A (zh) | 基于电动汽车接入的多能源系统协同调度方法及装置 | |
CN114155103A (zh) | 基于区块链合作博弈的能量共享联盟灵活性交易方法 | |
CN112365021A (zh) | 一种基于混合储能的区域综合能源系统规划优化方法 | |
CN112084629A (zh) | 一种基于两阶段鲁棒优化的多能虚拟电厂聚合方法 | |
Chen et al. | Information gap-based coordination scheme for active distribution network considering charging/discharging optimization for electric vehicles and demand response | |
CN115860413A (zh) | 基于负荷需求响应与双层可调鲁棒优化的并网冷热电联产微网系统经济调度方法 | |
CN115577909A (zh) | 考虑价格型需求响应和v2g的园区综合能源系统调度方法 | |
CN105279575B (zh) | 基于发电预测的多能源主体分布式博弈优化方法 | |
Tong et al. | An intelligent scheduling control method for smart grid based on deep learning | |
CN116402307A (zh) | 考虑可调度柔性资源运行特性的电网规划容量分析方法 | |
Wan et al. | Hierarchical multiobjective dispatching strategy for the microgrid system using modified MOEA/D | |
Qiu et al. | Local integrated energy system operational optimization considering multi‐type uncertainties: A reinforcement learning approach based on improved TD3 algorithm | |
CN114759616B (zh) | 一种考虑电力电子器件特性的微电网鲁棒优化调度方法 | |
CN105514986A (zh) | 一种基于虚拟电厂技术的der用户竞价并网方法 | |
CN115693797A (zh) | 考虑v2g与需求响应的配电网调度方法、介质及系统 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171103 |
|
WD01 | Invention patent application deemed withdrawn after publication |