CN103473604B - 一种日前发电计划综合性能评价方法 - Google Patents
一种日前发电计划综合性能评价方法 Download PDFInfo
- Publication number
- CN103473604B CN103473604B CN201310370259.0A CN201310370259A CN103473604B CN 103473604 B CN103473604 B CN 103473604B CN 201310370259 A CN201310370259 A CN 201310370259A CN 103473604 B CN103473604 B CN 103473604B
- Authority
- CN
- China
- Prior art keywords
- index
- unit
- membership function
- value
- few days
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000005611 electricity Effects 0.000 claims description 49
- 230000009194 climbing Effects 0.000 claims description 42
- 239000003245 coal Substances 0.000 claims description 32
- 230000002950 deficient Effects 0.000 claims description 27
- 238000010977 unit operation Methods 0.000 claims description 19
- 230000005283 ground state Effects 0.000 claims description 18
- 230000014509 gene expression Effects 0.000 claims description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 230000006978 adaptation Effects 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 6
- GOLXNESZZPUPJE-UHFFFAOYSA-N spiromesifen Chemical compound CC1=CC(C)=CC(C)=C1C(C(O1)=O)=C(OC(=O)CC(C)(C)C)C11CCCC1 GOLXNESZZPUPJE-UHFFFAOYSA-N 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 230000005619 thermoelectricity Effects 0.000 claims description 5
- 230000007812 deficiency Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000000205 computational method Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000007726 management method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及电力工业领域的评价方法,具体涉及一种日前发电计划综合性能评价方法。该方法包括下述步骤:①建立评价指标体系;②确定各评价指标的隶属度函数;③确定日前发电计划的总体评分。本发明提供的方法制定出一套适合日前发电计划的评估指标体系,并根据各指标的特点提出每类指标的隶属度函数计算方法,结合各指标最终权重计算获得评估结果,本发明所提供的评估方法将客观评价每种计划方案的优劣性,使计划的编制具有可供对比分析的参考体系。
Description
技术领域
本发明涉及电力工业领域的评价方法,具体涉及一种日前发电计划综合性能评价方法。
背景技术
电力工业中节能经济调度是提高电网运行可靠率,实现节能减排,提高调度计划精细化管理水平的有效技术手段。根据实际需求,目前节能经济调度已在我国得到全面推广,各级调度中心根据本地需求构建了适应本地特色的节能经济调度模型。然而截至目前为止尚未出现对日前发电计划进行优劣评估的方法,因此日前发电计划缺少进行对比分析的参照体系。
发明内容
针对现有技术的不足,本发明的目的是提供一种日前发电计划综合性能评价方法,该方法制定出一套适合日前发电计划的评估指标体系,并根据各指标的特点提出每类指标的隶属度函数计算方法,结合各指标最终权重计算获得评估结果,本发明所提的评估方法将客观评价每种计划方案的优劣性,使计划的编制具有可供对比分析的参考体系。
本发明的目的是采用下述技术方案实现的:
本发明提供一种日前发电计划综合性能评价方法,其改进之处在于,所述方法包括下述步骤:
①建立评价指标体系;
②确定各评价指标的隶属度函数;
③确定日前发电计划的总体评分。
优选的,所述步骤①中,评价指标包括三个一级指标:安全性能指标、经济性能指标和适应性能指标。
较优选的,所述安全性指标包括二级指标:机组运行安全性指标和电网运行安全性指标;
所述机组运行安全性指标包括三级指标:机组限额指标、机组禁止运行区指标和机组爬坡指标;所述电网运行安全性指标包括三级指标:系统备用指标、系统爬坡指标、系统基态潮流及N-1潮流指标;其中基态潮流及N-1潮流包括四级指标:支路潮流指标和断面潮流指标;
所述经济性能指标包括二级指标:机组运行经济性指标和电网运行经济性指标;
所述机组运行经济性能指标包括三级指标:机组平均发电煤耗指标和机组平均负载率指标;所述电网运行经济性指标包括三级指标:平均网损指标、弃风电量指标和弃水电量指标;
所述适应性能指标包括二级指标:模式适应性指标和运行适应性指标;
所述模式适应性指标包括三级指标:“三公”电量偏差比指标;所述运行适应性指标包括三级指标:机组出力平滑性指标和机组计划执行率指标。
优选的,所述步骤②中,二级指标、三级指标和四级指标统称为中间指标,所述中间指标的由下层指标累积求取,底层指标按照隶属度函数求取;
分别确定机组限额指标、机组禁止运行区指标、机组爬坡指标、系统备用指标、系统爬坡指标、基态潮流及N-1潮流、机组平均发电煤耗指标、机组平均负载率指标、平均网损指标、弃风电量指标、弃水电量指标、“三公”电量偏差比指标、机组出力平滑性指标和机组计划执行率指标的隶属度函数;
为随着评价值的增大隶属度趋于1的类型,称为递增隶属度;为随着评价值减少隶属度趋于1的类型,称为递减隶属度;
对递增隶属度采用正弦函数表示,对递减隶属度采用余弦函数表示,表达式分别如下:
其中:n为自变量,n0和n1为常数;当n趋近于n1时,μ趋近于1;当n趋近于n0时,μ趋近于1。
较优选的,机组限额指标是考核日前计划中有多少机组是超过机组限额安排计划的;确定机组限额指标隶属度函数包括:机组限额指标隶属度函数自变量n用机组限额的不合格数cu表示:
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,max为机组i的出力上限,pi,min为机组i的出力下限;cu越少越好,根据cu计算隶属度函数时选择余弦函数, 其中n0=0,n1根据系统计算规模决定,取值为机组总台数的1/10。
较优选的,机组禁止运行区指标是考核日前计划中有多少机组被安排在禁止运行区;确定机组运行区指标隶属度函数包括:机组运行区指标隶属度函数自变量n用运行在禁止运行区的机组的不合格数cu表示:
cu=cu+1 if pi,t∈([pi,1,pi,2],[pi,3,pi,4]...) (4);
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,1,pi,2,pi,3,pi,4为禁止运行区间出力值;cu越少越好,根据cu计算隶属度函数时选择余弦函数,其中n0=0,n1=1。
较优选的,机组爬坡约束指标是考核日前计划中有多少机组计划临近两点出力变化超出机组爬坡速率;确定机组爬坡约束指标隶属度函数包括:机组爬坡约束指标隶属度函数的自变量n用机组爬坡不合格数cu表示:
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,up为上爬坡限值,pi,down为下爬坡限值;cu越少越好,根据cu计算隶属度函数时选择余弦函数;其中n0=0,n1为机组总台数的1/20。
较优选的,系统备用指标是考核日前计划能否满足系统正负备用,当不满足系统备用要求时,将使系统抗风险能力差,降低系统安全系数,确定系统备用指标隶属度函数包括:系统备用指标隶属度函数的自变量n用备用考核不合格数cs表示:
式中,cs为系统备用不合格累计值,初始值为0;pi,t,max为机组i在时刻t的上限发电能力;pi,t,min为机组i在时刻t的下限发电能力;pi,t为机组i在时刻t的出力;Rt,up为系统在时刻t的上备用要求,Rt,down为系统在时刻t的下备用要求;cs越少越好,在隶属度函数计算时采用余 弦函数;其中n0=0,n1的取值根据评估对象电网构成特性取值。
较优选的,系统爬坡指标是考核日前计划下系统总出力的爬坡能力,当系统爬坡能力不足时,影响系统抗波动能力;确定系统爬坡指标隶属度函数包括:系统爬坡指标隶属度函数的自变量n用不合格数cs表示:
式中,cs为系统爬坡不合格累计值,初始值为0;pi,t为机组i在时刻t的出力;PSt,up,PSt,down分为系统是时刻t的上下爬坡要求;cs越少越好,在隶属度函数计算时采用余弦函数;将爬坡能力少于系统总负荷的5%的点作为爬坡不合格点;其中n0=0,n1的取值根据评估对象电网构成特性取值。
较优选的,定义计划日基态或N-1潮流在任意时刻出现超过支路限额5%的情况,则该支路为计划潮流不合格支路,确定基态潮流及N-1潮流隶属度函数包括:基态潮流及N-1潮流隶属度函数的自变量n用计划潮流不合格支路表示:
式中,pt为各个节点在时刻t注入有功功率;X为网络的节点阻抗矩阵;lt,i,j为节点i,j之间支路在时刻t的有功潮流;xi,j为节点i,j之间支路的电抗值;li,j,max,li,j,min分别为支路i,j的上下限额;cs为系统潮流不合格支路累计值,不合格支路越少越好,采用选择余弦函数作为隶属度函数;其中n0=0,n1根据系统计算规模决定,取值为支路总数的1/20,基态潮流及N-1潮流隶属度函数的计算同样适用于断面潮流评估。
较优选的,机组平均发电煤耗指标针对火电机组;评估火电机组的节能减排效果是计算机组的平均发电煤耗,即系统每千瓦时所消耗的标准煤耗,确定机组平均发电煤耗指标的隶属度函数包括:机组平均发电煤耗指标的隶属度函数的自变量n用系统平均发电煤耗cm表示:
式中,M表示系统总标准煤耗;Qu表示系统火电总发电量;cm表示系统平均发电煤耗;cm越低越好,采用选择余弦函数作为隶属度函数;取n0=320,n1=360。
较优选的,机组负载率指标针对火电机组,火电机组负载率反映机组运行经济性的辅助指标;确定机组负载率指标隶属度函数包括:机组负载率指标隶属度函数的自变量n用机组负载率cm表示:
式中,Qu表示火电机组总发电量;Qmax表示火电机组最大可发电量,cm越大越好,采用选择正弦函数作为隶属度函数,取n0=0.7,n1=0.9。
较优选的,确定平均网损指标的隶属度函数包括:平均网损隶属度函数的自变量n用平均网损表示:
式中,Snet为全网计划周期内损耗;Qnet为全网计划周期内机组上网总电量;网络损耗越低则表明电网运行越经济,采用选择余弦函数作为隶属度函数;对500kV电网取n0=0.01,n1=0.03,对220kV电网取n0=0.03,n1=0.05。
较优选的,确定弃风电量和弃水电量指标隶属度函数包括:弃风电量和弃水电量指标隶属度函数自变量n用下述表达式表示:
式中,pf,r,t,pp,r,t分别为新能源在时刻t的预测电力和计划电力;R代表新能源总数;风电、水电弃用越少越好,采用余弦隶属度函数;其中n0=0,n1的取值根据风电和水电的预测发电量确定,n1为风电和水电的预测发电量的30%。
较优选的,“三公”电量偏差指标衡量“三公”电量偏差一致性,确定“三公”电量偏差指标隶属度函数包括:“三公”电量偏差指标隶属度函数的自变量n用下述表达式表示:
式中,Qp,QH分别表示计划电量和期望电量;εi表示偏差最小电厂的比例;当出现一个电厂的偏差大于最小偏差的20%时,则认为该电厂不合格,不合格电厂越少越好,采用余弦隶属度函数;其中n0=0,n1的取值根据“三公”电厂总数确定,取值为“三公”电厂总数的20%。
较优选的,确定机组出力平滑性指标隶属度函数包括:机组出力平滑性指标隶属度函数的自变量n用下述表达式表示:
式中,ci , p表示机组的毛刺点;pi,t代表机组i在时刻t的出力;cs代表机组出力平滑不合格机组数目;cs越少越好,采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,取值为火电机组总数的20%。
较优选的,机组计划执行率指标是指机组实际执行发电与计划电力的偏差;确定机组计划执行率指标隶属度函数包括:机组计划执行率指标隶属度函数的自变量n用下述表达式表示:
式中,pi,t,plan,pi,t,real分别表示机组在时刻t的计划出力和实际出力;slt,plan,slt,real分别表示在时刻t的系统负荷预测和实际负荷;ci,p表示机组的偏差不合格点;cs表示计划执行整体不合格点,越少越好,采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,取值为火电机组总数的20%。
优选的,所述步骤③中,确定日前发电计划的总体评分包括:
根据步骤②中各指标的隶属度函数,并结合各指标在总体中所占的评分比重得出日前发电计划的总体评分,其表达式如下:
式中,N表示评分指标总数,μn表示指标n的隶属度值,λn表示指标n所占的计分比值各指标的λ之和为1;κ为总体评分;对于λn的取值根据各地的实际情况赋值。
与现有技术比,本发明达到的有益效果是:
1、本发明首先提供了针对日前发电计划综合性能的一套指标评估体系,该指标体系从计划安全性、经济性及实用性三个大方面全方位的描述了日前计划各类性能,能够较好的对日前发电计划的结果性能进行全面概括。
2、本发明提出的针对日前计划结果性能的评估方法可以根据各地的电网实际情况,调度计划区域业务特色调整隶属度计算参数及指标权重,从而计算获得各具参考及指导意义的性能评分结果。树状的多层级指标管理体系可以对指标进行分类管理灵活删减,为指标体系的进一步完善提供了良好的扩充性。
3、本发明针对日前发电计划结果性能所提出的评估方法,既为计划制定部门提供了提高计划质量的参考体系,也为第三方监督机构提供了有效的监控手段,为进一步推进电网节能减排提供了技术辅助手段。
附图说明
图1是本发明提供的日前发电计划综合性能评价方法的流程图;
图2是本发明提供的安全性能指标示意图;
图3是本发明提供的经济性能指标示意图;
图4是本发明提供的适应性能指标示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本发明的目的在于对现有的节能经济调度日前发电计划综合性能做出合理性评价。电力工业中节能经济调度是提高电网运行可靠率,实现节能减排,提高调度计划精细化管理水平的有效技术手段。根据实际需求,目前节能经济调度已在我国得到全面推广,各级调度中心根据本地需求构建了适应本地特色的节能经济调度模型。然而截至目前为止尚未出现对日前发电计划进行优劣评估的方法,因此日前发电计划缺少进行对比分析的参照体系。本发明的目的还在于通过分析发电计划自身模型特性与工程实用需求,在构建评估分析指标基础上建立指标评分细则,从而获取整体评价体系,最终给出综合评估结论。
本发明提供的日前发电计划综合性能评价方法的流程图如图1所示,包括下述步骤:
①建立评价指标体系:
日前发电计划是根据电网系统负荷预测,在综合考虑火电、水电、风电等各种能源运行特性,以电网安全运行为基础,以满足系统整体运行最优性为目标的优化问题。日前发电计划是电力系统调度计划类的核心业务,其涉及面较广,为确保指标体系的全面性,本发明将指标进行分层分级管理。
首先划分三类一级指标:安全性能指标,经济性能指标和适应性能指标。
安全性能下含二级指标:机组运行安全指标和电网运行安全指标;
机组运行安全下含三级指标:机组限额,机组禁止运行区,机组爬坡;电网运行安全下含三级指标:系统备用,系统爬坡,系统基态潮流、系统N-1潮流;基态潮流和N-1潮流包含四级指标:支路潮流、断面潮流,本发明提供的安全性能指标示意图如图2所示。
经济性能下含二级指标:机组运行经济性和电网运行经济性;机组运行经济性包含三级指标:机组平均发电煤耗,机组平均负载率;电网运行经济性包含三级指标:平均网损,弃风电量,弃水电量。本发明提供的经济性能指标示意图如图3所示。
适应性性能下含二级指标:模式适应性和运行适应性;模式适应性下含三级指标:“三公”电量偏差比;运行适应性下含三级指标:机组出力平滑性,机组计划执行率。本发明提供的适应性能指标示意图如图4所示。
②确定各评价指标的隶属度函数:
由于指标分层分级管理,因此中间指标将由下层指标累积求取,而底层指标则按照隶属度函数求取。
底层指标隶属度函数。隶属度函数是用来表征不同自变量取值对因变量的影响程度,本文通过隶属度函数来评价各项指标在不同计划结果下得分。在本发明中所有指标都可以分为两种类型,一种为随着评价值的增大隶属度趋于1的类型,称为递增隶属度;一种为随着评价值减少隶属度趋于1的类型,称为递减隶属度。根据实际情况,本发明对递增隶属度采用正弦函数表示,如式(1)所示;而递减隶属度采用余弦函数表示,如式(2)所示。
计算公式中n为自变量,n0和n1为常数。从计算公式可见,对于公式(1),当n趋近于n1时,μ趋近于1;对于公式(2),当n趋近于n0时,μ趋近于1。对于不同的评价指标只需要选择不同的隶属度函数,并确定合适的n0和n1即可计算出隶属度。
计算各个指标的隶属度主要是确定各个指标的自变量n,具体计算如下:
一)机组限额评估。该指标主要是考核日前计划中有多少机组是超过机组限额安排计划的,不合格数计算公式如式(3)所示。
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,max为机组i的出力上限,pi,min为机组i的出力下限。cu越少越好,根据cu计算隶属度函数时选择余弦函数,其中n0=0,n1可以根据系统计算规模决定,在此取值为机组总台数的1/10。
二)机组禁止运行区评估。该指标主要是考核日前计划中有多少机组被安排在禁止运行区。不合格数计算公式如式(4)所示。
cu=cu+1 if pi,t∈([pi,1,pi,2],[pi,3,pi,4]...) (4);
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,1,pi,2,pi,3,pi,4为禁止运行区间出力值。cu越少越好,根据cu计算隶属度函数时选择余弦函数,其中n0=0,n1=1。
三)机组爬坡约束评估。该指标主要是考核日前计划中有多少机组计划临近两点出力变化超出机组爬坡速率。爬坡不合格数计算公式如式(5)所示。
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,up为上爬坡限值,pi,down为下爬坡限值。cu越少越好,根据cu计算隶属度函数时选择余弦函数。实际运行时,超过爬坡约束的部分将由其他机组分担或由联络线分担,因此会影响计划执行。其中n0=0,n1为机组总台数的1/20。
四)系统备用评估。该指标主要是考核日前计划能否满足系统正负备用,当不满足系统 备用要求时,将使系统抗风险能力较差,降低系统安全系数。备用考核不合格数目计算公式如式(6)所示。
式中,cs为系统备用不合格累计值,初始值为0;pi,t,max为机组i在时刻t的上限发电能力;pi,t,min为机组i在时刻t的下限发电能力;pi , t为机组i在时刻t的出力;Rt,up为系统在时刻t的上备用要求,Rt,down为系统在时刻t的下备用要求。cs越少越好,在隶属度函数计算时采用余弦函数。其中n0=0,n1的取值需要根据评估对象电网构成特性取值,当评估对象与外界电网联系较为紧密时,该值可以取得稍大一点。当评估对象为孤岛电网时,该值可以取得较小。
五)系统爬坡评估。该指标主要是考核日前计划下系统总出力的爬坡能力,当系统爬坡能力不足时,将影响系统抗波动能力。爬坡不合格数计算公式如式(7)所示。
式中,cs为系统爬坡不合格累计值,初始值为0;pi,t为机组i在时刻t的出力;PSt,up,PSt,down分为系统是时刻t的上下爬坡要求。cs越少越好,在隶属度函数计算时采用余弦函数。本发明将爬坡能力少于系统总负荷的5%的点作为爬坡不合格点。其中n0=0,n1的取值需要根据评估对象电网构成特性取值,当电网包含较多随机负荷或电源时,该值需要取得较小;当电网主要以火电、水电构成并与外网联系紧密时,该值可以取得较大。
六)基态潮流及N-1潮流评估。电力系统计划编制结果首先需要满足电网基态运行下的潮流约束,其次需要满足电网在N-1方式下的潮流约束。本文定义计划日基态或N-1潮流在任意时刻出现超过支路限额5%的情况,则定义该支路为计划潮流不合格支路。不合格支路数计算公式如式(8)所示。
式中,式中,pt为各个节点在时刻t注入有功功率;X为网络的节点阻抗矩阵;lt,i,j为节点i,j之间支路在时刻t的有功潮流;xi,j为节点i,j之间支路的电抗值;li,j,max,li,j,min分别为支路i,j的上下限额;cs为系统潮流不合格累计值,因此不合格支路越少越好,所以采用选择余弦函数作为隶属度函数。其中n0=0,n1可以根据系统计算规模决定,在此取值为支路总数的1/20。上述评价方法同样适用于断面潮流评估。
七)机组平均发电煤耗评估。该指标只针对火电机组。评估机组的节能减排效果主要是计算机组的平均发电煤耗,即系统每千瓦时所消耗的标准煤耗。系统平均发电煤耗计算公式如式(9)所示。
式中M表示系统总标准煤耗;Qu表示系统火电总发电量;cm表示系统平均发电煤耗;cm越低越好,所以采用选择余弦函数作为隶属度函数。根据煤耗记录,世界最先进的百万千瓦级超临界机组的煤耗在300克/千瓦时左右,本发明取n0=320,n1=360。
八)机组负载率评估。该指标只针对火电机组。机组负载率同样是反映机组运行经济性的辅助指标。由于机组的发电煤耗更多决定于机组自身的物理属性,对于大容量高性能机组比例少的电网该指标很难达到优质,但是机组负载率则只跟计划安排相关,因此也是评估一个电网计划编制经济性的辅助指标。机组负载率的计算公式如式(10)所示。
式中,Qu表示火电机组总发电量;Qmax表示火电机组最大可发电量。通常情况下,cm越大越好,所以采用选择正弦函数作为隶属度函数。考虑到系统正备用与爬坡要求,系统平均负载率很难达到1,本发明取n0=0.7,n1=0.9。
九)平均网损评估。不同的发电计划对应不同的潮流分布,不同的潮流分布将致使不同的网络损耗,平均网损计算公式如式(11)所示。
式中,Snet为全网计划周期内损耗;Qnet为全网计划周期内机组上网总电量。网络损耗越 低则表明电网运行越经济,所以采用选择余弦函数作为隶属度函数。不同的电压等级其网损是不一样的,考虑到目前输电技术条件下网损是难以避免的,本发明对500KV电网取n0=0.01,n1=0.03,对220KV电网取n0=0.03,n1=0.05。
十)弃风、弃水电量评估。风电和水电作为清洁可再生能源,电网在保证安全运行条件下应该给予全额接纳,因此对风电和水电的弃用将作为电网运行经济性的考核指标。弃风弃水电量计算方法如式(12)所示。
式中,pf,r,t,pp,r,t分别为新能源在时刻t的预测电力和计划电力;R代表新能源总数;风电、水电弃用越少越好,所以采用余弦隶属度函数。其中n0=0,n1的取值根据风电和水电的预测发电量确定,n1为风电和水电的预测发电量的30%。
十一)“三公”电量偏差比评估。“三公”调度是我国主要的调度模式,因此计划对“三公”电量的把握是对计划结果适应性的重要参考指标。本发明采用的是“三公”电量偏差比指标,它是指计划电量和从中长期得到的期望电量的偏差。由于期望电量没有考虑计划日的实际电网运行情况,因此出现计划电量和期望电量不一致是正常现象,但是这种偏差对于每一个“三公”电厂而言应该是一致的,所以该指标主要是衡量“三公”电量偏差一致性。计算公式如式(13)所示。
式中,Qp,QH分别表示计划电量和期望电量;εi表示偏差最小电厂的比例;当出现一个电厂的偏差大于最小偏差的20%时,则认为该电厂不合格,不合格电厂越少越好,因此采用余弦隶属度函数。其中n0=0,n1的取值根据“三公”电厂总数确定,可以取值为“三公”电厂总数的20%。
十二)机组出力平滑性评估。由于发电机组有着自己的运行要求,因此机组的计划出力曲线应该尽可能的平滑,而不应该出现大量毛刺从而影响计划的实际执行。机组出力平滑性作为对机组计划适应性水平指标之一,计算公式如式(14)所示。
式中,ci,p表示机组的毛刺点;pi,t代表机组i在时刻t的出力;cs代表机组出力平滑不合格机组数目。cs越少越好,所以采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,可以取值为火电机组总数的20%。
十三)机组计划执行率评估。所谓机组计划执行率是指机组实际执行发电与计划电力的偏差。这个偏差综合反映了计划编制的适应性与合理性。当然,偏差的出现也可能是由于负荷预测的偏差以及其他人为原因所造成。该偏差指标的计算公式如式(15)所示。
式中,pi,t,plan,pi,t,real分别表示机组在时刻t的计划出力和实际出力;slt,plan,slt,real分别表示在时刻t的系统负荷预测和实际负荷;ci,p表示机组的偏差不合格点;cs表示计划执行整体不合格点。越少越好,因此采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,可以取值为火电机组总数的20%。
③确定日前发电计划的总体评分:
步骤②描述了各类指标的隶属度计算方法,此时只要再给出各指标在总体中所占的评分比重即可计算获得总体评分,总体评分计算公式如式(16)所示。
式中,N表示评分指标总数,μn表示指标n的隶属度值,λn表示指标n所占的计分比值,各指标的λ之和为1;κ为最终得分。对于λn的取值可以根据各地的实际情况进行给值,主要在于侧重的评估内容,当更重视安全性能时可以对安全指标给予较大的值;当更重视经济性能时可以对经济指标给予较大的值。例如我们选取系统备用,系统支路基态潮流作为评估的安全指标,选择系统平均煤耗,系统弃水量、弃风量作为评估的经济指标。假设现已经计算得到系统备用的隶属度为0.995,支路基态潮流隶属度为0.99,平均煤耗隶属度为0.95,弃水量隶属度为1,弃风量隶属度为0.98。当评估人员认为安全性与经济性同等重要时各指标的λn均取为0.2,此时最终评分:
κ=0.995×0.2+0.99×0.2+0.95×0.2+1×0.2+0.98×0.2=0.983;
当评估人员认为系统安全性更为重要,给系统安全指标赋予更大权重时,此时最终评分 为κ=0.995×0.4+0.99×0.3+0.95×0.1+1×0.1+0.98×0.1=0.988。由此可见对于更为注重安全性的计划,该计划结果具有更高的评估得分。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (15)
1.一种日前发电计划综合性能评价方法,其特征在于,所述方法包括下述步骤:
①建立评价指标体系;
②确定各评价指标的隶属度函数;
③确定日前发电计划的总体评分;
所述步骤①中,评价指标包括三个一级指标:安全性能指标、经济性能指标和适应性能指标;
所述安全性能指标包括二级指标:机组运行安全性指标和电网运行安全性指标;
所述机组运行安全性指标包括三级指标:机组限额指标、机组禁止运行区指标和机组爬坡指标;所述电网运行安全性指标包括三级指标:系统备用指标、系统爬坡指标、系统基态潮流及N-1潮流指标;其中基态潮流及N-1潮流包括四级指标:支路潮流指标和断面潮流指标;
所述经济性能指标包括二级指标:机组运行经济性指标和电网运行经济性指标;
所述机组运行经济性指标包括三级指标:机组平均发电煤耗指标和机组平均负载率指标;所述电网运行经济性指标包括三级指标:平均网损指标、弃风电量指标和弃水电量指标;
所述适应性能指标包括二级指标:模式适应性指标和运行适应性指标;
所述模式适应性指标包括三级指标:“三公”电量偏差比指标;所述运行适应性指标包括三级指标:机组出力平滑性指标和机组计划执行率指标;
所述步骤②中,二级指标、三级指标和四级指标统称为中间指标,所述中间指标的由下层指标累积求取,底层指标按照隶属度函数求取;
分别确定机组限额指标、机组禁止运行区指标、机组爬坡指标、系统备用指标、系统爬坡指标、基态潮流及N-1潮流、机组平均发电煤耗指标、机组平均负载率指标、平均网损指标、弃风电量指标、弃水电量指标、“三公”电量偏差比指标、机组出力平滑性指标和机组计划执行率指标的隶属度函数;
为随着评价值的增大隶属度趋于1的类型,称为递增隶属度;为随着评价值减少隶属度趋于1的类型,称为递减隶属度;
对递增隶属度采用正弦函数表示,对递减隶属度采用余弦函数表示,表达式分别如下:
其中:n为自变量,n0和n1为常数;当n趋近于n1时,μ趋近于1;当n趋近于n0时,μ趋近于1。
2.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组限额指标是考核日前计划中有多少机组是超过机组限额安排计划的;确定机组限额指标隶属度函数包括:机组限额指标隶属度函数自变量n用机组限额的不合格数cu表示:
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,max为机组i的出力上限,pi,min为机组i的出力下限;cu越少越好,根据cu计算隶属度函数时选择余弦函数,其中n0=0,n1根据系统计算规模决定,取值为机组总台数的1/10。
3.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组禁止运行区指标是考核日前计划中有多少机组被安排在禁止运行区;确定机组运行区指标隶属度函数包括:机组运行区指标隶属度函数自变量n用运行在禁止运行区的机组的不合格数cu表示:
cu=cu+1 if pi,t∈([pi,1,pi,2],[pi,3,pi,4]...)(4);
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,1,pi,2,pi,3,pi,4为禁止运行区间出力值;cu越少越好,根据cu计算隶属度函数时选择余弦函数,其中n0=0,n1=1。
4.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组爬坡约束指标是考核日前计划中有多少机组计划临近两点出力变化超出机组爬坡速率;确定机组爬坡约束指标隶属度函数包括:机组爬坡约束指标隶属度函数的自变量n用机组爬坡不合格数cu表示:
式中,cu为机组台数累计值,初始值为0;pi,t为机组i在时刻t的出力;pi,up为上爬坡限值,pi,down为下爬坡限值;cu越少越好,根据cu计算隶属度函数时选择余弦函数;其中n0=0,n1为机组总台数的1/20。
5.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,系统备用指标是考核日前计划能否满足系统正负备用,当不满足系统备用要求时,将使系统抗风险能力差,降低系统安全系数,确定系统备用指标隶属度函数包括:系统备用指标隶属度函数的自变量n用备用考核不合格数cs表示:
式中,cs为系统备用不合格累计值,初始值为0;pi,t,max为机组i在时刻t的上限发电能力;pi,t,min为机组i在时刻t的下限发电能力;pi,t为机组i在时刻t的出力;Rt,up为系统在时刻t的上备用要求,Rt,down为系统在时刻t的下备用要求;cs越少越好,在隶属度函数计算时采用余弦函数;其中n0=0,n1的取值根据评估对象电网构成特性取值。
6.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,系统爬坡指标是考核日前计划下系统总出力的爬坡能力,当系统爬坡能力不足时,影响系统抗波动能力;确定系统爬坡指标隶属度函数包括:系统爬坡指标隶属度函数的自变量n用不合格数cs表示:
式中,cs为系统爬坡不合格累计值,初始值为0;pi,t为机组i在时刻t的出力;PSt,up,PSt,down分为系统是时刻t的上下爬坡要求;cs越少越好,在隶属度函数计算时采用余弦函数;将爬坡能力少于系统总负荷的5%的点作为爬坡不合格点;其中n0=0,n1的取值根据评估对象电网构成特性取值。
7.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,定义计划日基态或N-1潮流在任意时刻出现超过支路限额5%的情况,则该支路为计划潮流不合格支路,确定基态潮流及N-1潮流隶属度函数包括:基态潮流及N-1潮流隶属度函数的自变量n用计划潮流不合格支路表示:
式中,pt为各个节点在时刻t注入有功功率;X为网络的节点阻抗矩阵;lt,i,j为节点i,j之间支路在时刻t的有功潮流;xi,j为节点i,j之间支路的电抗值;li,j,max,li,j,min分别为支路i,j的上下限额;cs为系统潮流不合格支路累计值,不合格支路越少越好,采用选择余弦函数作为隶属度函数;其中n0=0,n1根据系统计算规模决定,取值为支路总数的1/20,基态潮流及N-1潮流隶属度函数的计算同样适用于断面潮流评估。
8.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组平均发电煤耗指标针对火电机组;评估火电机组的节能减排效果是计算机组的平均发电煤耗,即系统每千瓦时所消耗的标准煤耗,确定机组平均发电煤耗指标的隶属度函数包括:机组平均发电煤耗指标的隶属度函数的自变量n用系统平均发电煤耗cm表示:
式中,M表示系统总标准煤耗;Qu表示系统火电总发电量;cm表示系统平均发电煤耗;cm越低越好,采用选择余弦函数作为隶属度函数;取n0=320,n1=360。
9.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组负载率指标针对火电机组,火电机组负载率反映机组运行经济性的辅助指标;确定机组负载率指标隶属度函数包括:机组负载率指标隶属度函数的自变量n用机组负载率cm表示:
式中,Qu表示火电机组总发电量;Qmax表示火电机组最大可发电量,cm越大越好,采用选择正弦函数作为隶属度函数,取n0=0.7,n1=0.9。
10.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,确定平均网损指标的隶属度函数包括:平均网损隶属度函数的自变量n用平均网损表示:
式中,Snet为全网计划周期内损耗;Qnet为全网计划周期内机组上网总电量;网络损耗越低则表明电网运行越经济,采用选择余弦函数作为隶属度函数;对500kV电网取n0=0.01,n1=0.03,对220kV电网取n0=0.03,n1=0.05。
11.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,确定弃风电量和弃水电量指标隶属度函数包括:弃风电量和弃水电量指标隶属度函数自变量n用下述表达式表示:
式中,pf,r,t,pp,r,t分别为新能源在时刻t的预测电力和计划电力;R代表新能源总数;风电、水电弃用越少越好,采用余弦隶属度函数;其中n0=0,n1的取值根据风电和水电的预测发电量确定,n1为风电和水电的预测发电量的30%。
12.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,“三公”电量偏差指标衡量“三公”电量偏差一致性,确定“三公”电量偏差指标隶属度函数包括:“三公”电量偏差指标隶属度函数的自变量n用下述表达式表示:
式中,Qp,QH分别表示计划电量和期望电量;εi表示偏差最小电厂的比例;当出现一个电厂的偏差大于最小偏差的20%时,则认为该电厂不合格,不合格电厂越少越好,采用余弦隶属度函数;其中n0=0,n1的取值根据“三公”电厂总数确定,取值为“三公”电厂总数的20%。
13.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,确定机组出力平滑性指标隶属度函数包括:机组出力平滑性指标隶属度函数的自变量n用下述表达式表示:
式中,ci,p表示机组的毛刺点;pi,t代表机组i在时刻t的出力;cs代表机组出力平滑不合格机组数目;cs越少越好,采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,取值为火电机组总数的20%。
14.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,机组计划执行率指标是指机组实际执行发电与计划电力的偏差;确定机组计划执行率指标隶属度函数包括:机组计划执行率指标隶属度函数的自变量n用下述表达式表示:
式中,pi,t,plan,pi,t,real分别表示机组在时刻t的计划出力和实际出力;slt,plan,slt,real分别表示在时刻t的系统负荷预测和实际负荷;ci,p表示机组的偏差不合格点;cs表示计划执行整体不合格点,越少越好,采用余弦隶属度函数,其中n0=0,n1的取值根据火电机组总数确定,取值为火电机组总数的20%。
15.如权利要求1所述的日前发电计划综合性能评价方法,其特征在于,所述步骤③中,确定日前发电计划的总体评分包括:
根据步骤②中各指标的隶属度函数,并结合各指标在总体中所占的评分比重得出日前发电计划的总体评分,其表达式如下:
式中,N表示评分指标总数,μn表示指标n的隶属度值,λn表示指标n所占的计分比值,各指标的λ之和为1;κ为总体评分;对于λn的取值根据各地的实际情况赋值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310370259.0A CN103473604B (zh) | 2013-08-22 | 2013-08-22 | 一种日前发电计划综合性能评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310370259.0A CN103473604B (zh) | 2013-08-22 | 2013-08-22 | 一种日前发电计划综合性能评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103473604A CN103473604A (zh) | 2013-12-25 |
CN103473604B true CN103473604B (zh) | 2016-08-10 |
Family
ID=49798447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310370259.0A Active CN103473604B (zh) | 2013-08-22 | 2013-08-22 | 一种日前发电计划综合性能评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103473604B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104090560B (zh) * | 2014-05-06 | 2017-02-08 | 内蒙古云谷电力科技股份有限公司 | 一种监测供电综合环境评价指标的设备 |
CN104239967A (zh) * | 2014-08-29 | 2014-12-24 | 华北电力大学 | 一种含风电场的电力系统多目标经济调度方法 |
CN104657910B (zh) * | 2014-12-29 | 2017-11-24 | 国网湖南省电力公司 | 一种电网每日运行评价方法 |
CN105005833A (zh) * | 2015-08-20 | 2015-10-28 | 长江水资源保护科学研究所 | 水电梯级开发对流域生态安全影响的评价方法 |
CN105184672A (zh) * | 2015-08-31 | 2015-12-23 | 同济大学 | 一种针对三公调度发电计划的评估方法 |
CN105205737A (zh) * | 2015-10-21 | 2015-12-30 | 东南大学 | 一种基于电站机组的经济性评价系统 |
CN107506895A (zh) * | 2017-07-26 | 2017-12-22 | 国网能源研究院 | 一种含新能源发电的电力系统运行经济性综合评价系统 |
CN108575497B (zh) * | 2018-04-02 | 2020-07-28 | 河海大学 | 一种基于模糊物元的刚性防浪林的树种选样方法 |
CN109754225B (zh) * | 2018-12-29 | 2023-08-29 | 昆明电力交易中心有限责任公司 | 一种日前发电计划制定流程的监控系统的设计方法及装置 |
CN115496417A (zh) * | 2022-10-31 | 2022-12-20 | 国网山东省电力公司青岛供电公司 | 一种海岛综合供能系统评价方法及系统 |
CN117154852B (zh) * | 2023-10-27 | 2024-01-02 | 华南理工大学 | 考虑n-1潮流约束下基于交替迭代的电网经济调度方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254091A (zh) * | 2011-06-10 | 2011-11-23 | 华北电力大学 | 20千伏中压配电方案评估方法 |
CN102411735A (zh) * | 2011-09-09 | 2012-04-11 | 河海大学常州校区 | 可重构装配系统重构规划方案评价方法 |
CN102982388A (zh) * | 2012-10-31 | 2013-03-20 | 中国电力科学研究院 | 一种日前电力系统经济调度的方法 |
CN103051001A (zh) * | 2013-01-07 | 2013-04-17 | 中国电力科学研究院 | 一种小周期实时发电计划 |
CN103077435A (zh) * | 2012-12-24 | 2013-05-01 | 江西省电力科学研究院 | 一种基于组合加权综合评价模型的sec指标评估方法 |
-
2013
- 2013-08-22 CN CN201310370259.0A patent/CN103473604B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254091A (zh) * | 2011-06-10 | 2011-11-23 | 华北电力大学 | 20千伏中压配电方案评估方法 |
CN102411735A (zh) * | 2011-09-09 | 2012-04-11 | 河海大学常州校区 | 可重构装配系统重构规划方案评价方法 |
CN102982388A (zh) * | 2012-10-31 | 2013-03-20 | 中国电力科学研究院 | 一种日前电力系统经济调度的方法 |
CN103077435A (zh) * | 2012-12-24 | 2013-05-01 | 江西省电力科学研究院 | 一种基于组合加权综合评价模型的sec指标评估方法 |
CN103051001A (zh) * | 2013-01-07 | 2013-04-17 | 中国电力科学研究院 | 一种小周期实时发电计划 |
Also Published As
Publication number | Publication date |
---|---|
CN103473604A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103473604B (zh) | 一种日前发电计划综合性能评价方法 | |
CN104331844B (zh) | 一种电网基建项目投资决策方法 | |
CN108428045A (zh) | 一种配电网运行健康状态评估方法 | |
CN109409702A (zh) | 一种电网运行指标分析系统 | |
CN107611966A (zh) | 一种考虑差异可靠性的有源配电网供电能力评估方法 | |
CN106127377A (zh) | 一种智能电网多能源综合协调水平评价方法 | |
CN103530823A (zh) | 基于otfn-ahp的微电网规划设计综合评价体系 | |
CN107634532A (zh) | 基于改进v2g和优先级调度的充放储一体站控制方法 | |
CN104751246A (zh) | 一种基于随机机会约束的有源配电网规划方法 | |
CN103577891B (zh) | 一种含分布式电源的多孤岛微网优化合作运行方法 | |
CN110334914A (zh) | 一种基于风险思想的源网荷储协调水平评价方法 | |
CN107194574A (zh) | 一种基于负荷损失的电网安全风险评估方法 | |
CN110119888A (zh) | 一种基于分布式电源接入的有源网格化规划方法 | |
CN106600136A (zh) | 一种电力断面越限控制效率评估方法 | |
CN111799794B (zh) | 一种考虑输电阻塞盈余的输电网扩展规划方法 | |
CN106096285B (zh) | 一种储能系统应对高风电渗透率系统调频需求效用评估方法 | |
CN110061496A (zh) | 源网荷多目标规划模型的可信度优化方法 | |
CN110245794A (zh) | 考虑灵活性的多能源汇集中心火储容量双层优化方法 | |
CN106549420A (zh) | 考虑风险和风力发电的电力系统运行备用优化方法 | |
CN114781159A (zh) | 分层分区的配电网光伏接入现状及承载能力评估方法 | |
CN106600130A (zh) | 基于模糊神经网络模型的独立微网安全性分析方法 | |
CN105678415A (zh) | 一种分布式电源配电网的净负荷预测方法 | |
CN105844359A (zh) | 一种基于可靠性成本效益分析的城市电网结构优化方法 | |
CN107563641A (zh) | 一种计及灾害偏好的配电网抗灾多场景差异化规划方法 | |
CN116014797A (zh) | 一种配网新能源接纳能力提升的评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |