CN114677031B - 一种基于数据包络分析的谐波治理需求评估方法 - Google Patents
一种基于数据包络分析的谐波治理需求评估方法 Download PDFInfo
- Publication number
- CN114677031B CN114677031B CN202210347350.XA CN202210347350A CN114677031B CN 114677031 B CN114677031 B CN 114677031B CN 202210347350 A CN202210347350 A CN 202210347350A CN 114677031 B CN114677031 B CN 114677031B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- loss
- efficiency
- transformer
- decision
- 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 42
- 238000004458 analytical method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000008901 benefit Effects 0.000 claims abstract description 5
- 239000003990 capacitor Substances 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 19
- 230000009467 reduction Effects 0.000 claims description 11
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 9
- 230000032683 aging Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 4
- 229910052802 copper Inorganic materials 0.000 claims description 4
- 239000010949 copper Substances 0.000 claims description 4
- 229910052742 iron Inorganic materials 0.000 claims description 4
- 238000004904 shortening Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 229910052731 fluorine Inorganic materials 0.000 claims description 3
- 125000001153 fluoro group Chemical group F* 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 238000004519 manufacturing process Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000001629 suppression Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241000220010 Rhode Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000003878 thermal aging Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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/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
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于数据包络分析的谐波治理需求评估方法,进行谐波损失分析和谐波治理需求分析,确定谐波治理需求综合评价指标集,包括输入指标和输出指标;采用DEA方法从输入指标与输出指标两个方面评估非线性用户和园区的谐波治理需求,对多个决策单元的相对优劣性进行评估;当出现大量或者全部决策单元都有效时,需考虑所有的可能权重来对决策单元进行评估,由此引入区间超效率模型;由于得到的是每个决策单元的效率区间,故采用基于极大极小后悔值法对决策单元进行比较和排序。本发明综合考虑多项影响指标,全面、系统、客观的评估非线性用户的用电环境以及谐波治理需求,进而提升用户的治理积极性,并对后续的治理起到指导作用。
Description
技术领域
本发明涉及电能质量技术领域,具体为一种基于数据包络分析的谐波治理需求评估方法。
背景技术
随着电力行业的逐步发展,电力系统的电力电子化率也逐渐提高。大量非线性负荷和冲击性负荷涌入电网,产生了严重的谐波污染问题。谐波对电力系统的影响主要包括谐波放大、谐振、使电力设备产生额外的附加损耗、缩短电力设备的使用寿命、干扰设备的正常运行及对测量带来误差和干扰等,更严重的是,谐波电流与系统侧阻抗相互作用进一步恶化公用电网的电压质量,影响与共用母线连接的其他电力负荷稳定运行,谐波污染已成为电网不可忽略的问题。目前,电力系统在谐波污染方面的两条解决途径是谐波管理和谐波治理,但是基于国家标准的谐波管理制度导致了非线性用户的治理积极性不足,电力谐波污染成为电力系统亟须解决的问题。
谐波问题给非线性用户和电网带来了许多困扰,尤其是大型工业园区用户。首先,谐波造成电能损耗,加大非线性用户的电费损失;其次,谐波会使得电力设备产生集肤效应,使得设备温度升高,缩短设备的使用寿命,严重时会造成设备损坏;再次,谐波电流会造成电容器组过载、引发谐振等风险,当发生谐振时,流入电网的谐波电流会高出数倍,对设备造成的电能损耗大大增加,甚至会造成非线性用户生产中断;最后,谐波电流不仅会导致电压谐波,而且其本身也会对电气设备产生影响,谐波可能会改变电网的电压过零点,一些高精密电子设备的运行会受到明显干扰,造成损失。
而多数高精尖制造业用户的生产线和生产设备,对供电可靠性和供电质量要求特别高,用电环境和生产设备运行参数稍有变化,就容易使其产品报废。多数工业用户的生产线都不同程度涵盖了类似电弧炉的冲击性和非线性负荷,对电网提供优质电能造成了消极影响。因此,引导用户主动参与治理电能质量问题,尽可能降低用户和电网的电能质量损失,为社会制造业、工业和地方经济的有力发展提供最大的保障。
目前,非线性用户谐波治理关键问题在于用户治理的需求分析不足。主要原因有两个:一是非线性用户是造成谐波污染问题的主导因素,但其自身对谐波问题的敏感度以及造成的经济损失不明确,谐波治理需求模糊,从而治理意愿薄弱;二是电能质量监测点大多分布在电网侧,导致用户获得的信息不对等,从而对自身的用电环境了解不足。
全面、系统地评估非线性用户的治理需求,指导谐波治理的实施,以实现让用户了解自身的用电环境并积极参与治理,解决日益严重的电力谐波污染问题。因此,亟须提出一种客观、全面、系统的谐波治理需求评估方法,以提升用户的治理积极性。
目前,对于电网用电环境的评估方法有层次分析法、模糊综合评判法、基于熵权的层次分析法等,但是传统方法确定指标的权重需依赖人为的经验和主观判断,缺乏系统性,因此需要更加客观、准确、稳定的评估方法。
发明内容
针对上述问题,本发明的目的在于为使非线性用户了解自身的用电环境,积极主动的参与治理,解决谐波污染问题,提出一种基于数据包络分析的谐波治理需求评估方法,综合考虑多项影响指标,全面、系统、客观的评估非线性用户的用电环境以及谐波治理需求,进而提升用户的治理积极性。技术方案如下:
一种基于数据包络分析的谐波治理需求评估方法,其特征在于,包括以下步骤:
步骤1:进行谐波损失分析和谐波治理需求分析,确定谐波治理需求综合评价指标集,包括输入指标和输出指标;
步骤2:采用DEA方法从输入指标与输出指标两个方面评估非线性用户和园区的谐波治理需求,对多个决策单元的相对优劣性进行评估;
步骤3:当出现大量或者全部决策单元都有效时,需考虑所有的可能权重来对决策单元进行评估,由此引入区间超效率模型;
步骤4:由于得到的是每个决策单元的效率区间,故采用基于极大极小后悔值法对决策单元进行比较和排序。
进一步的,所述进行谐波损失分析包括:
1)分析谐波引起的变压器和输电线路附加损耗:
输电线路上的附加损耗分为基波电流引起的损耗和谐波电流引起的损耗,而谐波损耗为基波损耗的kl倍;计算输电线路的谐波损耗:
其中:
式中,ΔPl为谐波引起的输电线路附加损耗;h为谐波电流次数;为最大谐波电流次数;R1为基波下的输电线路电阻;I1为基波电流;HRIh为h次谐波电流含有率,是谐波电流与基波电流的比率;Rh为h次谐波电阻;
变压器的损耗分为铜耗与铁耗,因谐波所产生的变压器附加损耗表示为:
式中,ΔPt为谐波引起的变压器附加损耗;h为谐波次数;为最大谐波次数;UAh、UBh、 UCh为A、B、C各相第h次谐波电压;IAh、IBh、ICh为A、B、C各相第h次谐波电流; 为A、B、C各相第h次谐波电压与谐波电流之间的相位差;
2)分析谐波引起的变压器降容
计算变压器的降容率:
其中:
式中,δt为谐波引起的变压器降容率;St、Vt为变压器的额定容量和额定电压;Sr、Vr为变压器在谐波条件下的运行容量和电压有效值;I*max为最大允许电流标幺值;P*LL为变压器负载损耗标幺值;P*EC为变压器绕组涡流损耗;FHL为谐波损耗因子;
3)分析谐波造成的电力设备寿命缩短
计算变压器的寿命损失率:
其中:
式中,LT为变压器在时间T内的寿命损失率;B为常数;θ(t)为t时刻的热点温度;FAA(t) 为t时刻的老化因子,需根据不同的变压器设计标准进行修正;
计算谐波条件下的电容器寿命:
式中,Lc为谐波条件下的电容器寿命;L0为正弦波形下标称电压时的电容器寿命;Kf为波形系数;Kp为峰值系数;Krms为均方根系数;a、b、c数值的大小表示三个因素对加速电容器老化的影响程度。
更进一步的,所述输入指标包括:
内部参数:变压器容量、输电线路长度、电容器容量;
负荷参数:负荷总容量、谐波源数量、谐波源容量和谐波源总谐波发射特性;
所述输出指标包括:
设备损耗:变压器附加损耗、输电线路附加损耗、变压器降容损耗、变压器寿命损耗、电容器寿命损耗;
国家标准:谐波电流含量裕度、谐波电压畸变率裕度;
潜在风险:谐振峰值数量、谐振程度。
更进一步的,所述谐振程度为:
式中,Ufh为h次谐波谐振情况下的谐波电压;U1为基波电压;
所述谐波电压畸变率裕度为:
式中,THDUGB为谐波电压畸变率限值;Uh为第h次谐波电压;
所述谐波电流含量裕度为:
式中,HRIGB为谐波电流含有率限值,Ih为第h次谐波电流。
更进一步的,所述超效率模型表示为:
式中,θ为输入相对于输出的有效利用率;ε为非阿基米德无穷小量;s-为输入指标的松弛变量;s+为输出变量的剩余变量;λj为决策单元的组合系数;n为评价对象的数量,即决策单元的数量;xj为第j个决策单元的输入指标,x0为第j0个决策单元的输入指标;yj为第j 个决策单元的输出指标,y0为第j0个决策单元的输出指标;e^T=(1,1,…,1)∈Em,eT=(1,1,…,1)∈Es;Em为m阶单位矩阵;Es为s阶单位矩阵。
若θ=1,s-=0,s+=0,则称为DEA有效;若θ=1,s-≠0或s+≠0,则称为DEA弱有效;若θ<1则称为非DEA有效;效率模型的效率指标没有小于等于1的约束,对于DEA有效的决策单元,其超效率值大于1,而对于非DEA有效的决策单元,其效率值不变;DEA有效表明谐波治理需求综合评估结果良好,非DEA有效为表明谐波治理需求综合评估结果差;
定义决策单元使用最优权重的交叉超效率值为效率区间上限;第d个决策单元DMUd的最大交叉超效率表示为:
其中,yrj为输出指标;xij为输入指标;s为输出指标数量;m为输入指标数量;μrd、ωid为相关联的输出和输入权重;θd为决策单元DMUd输入相对于输出的有效利用率;xid、yrd为被评估决策单元DMUd的输入指标和输出指标;为决策单元DMUd的效率区间上限;同理,其最小交叉超效率值可表示为:
其中,Edj为决策单元DMUd的效率区间下限;
由于得到的是每个决策单元的效率区间,故采用基于极大极小后悔值法,比较和排序决策单元的效率区间;
其中,rj为决策单元DMUj的最大损失效率;b为除自身外,其余决策单元的最大效率区间上限;
当E dj≥b时,决策者不会遭受损失,也不会后悔,此时rj=0;
结合上述两种情况,第j个决策单元的最大效率损失为:
因此,依据以下条件选择最佳效率区间:
进一步的,所述采用基于极大极小后悔值法对决策单元进行比较和排序,步骤如下:
Step1:计算每个DMU的最大效率损失,最小者记作Aj1,1≤j1≤n;
Step2:将Aj1排除,重新计算n-1个DMU的最大效率损失,最小者记作Aj2,1≤j2≤n,j2≠j1;
Step3:重复上述步骤,直到只剩下一个DMU;最终的排序为Aj1>Aj2>…>Ajn。
本发明的有益效果是:本发明提出一种基于数据包络分析的谐波治理需求评估方法,综合考虑多项影响指标,全面、系统、客观的评估非线性用户的用电环境以及谐波治理需求,进而提升用户的治理积极性;不需要预先估计参数、任何权重假设,避免了主观因素的影响;评估结果不仅可提高非线性用户的治理积极性,并对后续的治理起到指导作用。
附图说明
图1为谐波治理需求综合评估体系架构图。
图2为园区拓扑结构图。
图3为决策单元效率间隔图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。
1、谐波损失分析
(1)谐波引起的变压器和输电线路附加损耗
1)输电线路上的附加损耗可以分为基波电流引起的损耗和谐波电流引起的损耗,而谐波损耗为基波损耗的kl倍。输电线路的谐波损耗可表示为:
其中:
2)变压器的损耗分为铜耗与铁耗,其中铁耗是变压器铁芯内的励磁电流引起磁通变化产生的损耗,铜耗是负载电流通过配电变压器绕组产生的损耗。变压器谐波附加损耗的大小,一方面取决于运行中谐波电压和谐波电流的大小,另一方面还取决于配电变压器内部参数与频率的关系。因谐波所产生的变压器附加损耗可表示为:
式中,ΔPt为谐波引起的变压器附加损耗;h为谐波次数;为最大谐波次数;UAh、UBh、 UCh为A、B、C各相第h次谐波电压;IAh、IBh、ICh为A、B、C各相第h次谐波电流;为第h次谐波电压与谐波电流之间的相位差。变压器内部参数的频率特性与变压器的具体结构有关,目前尚无完全准确的测量方法。式(3)中相关参数可以在变压器低压侧通过测量获取,避免了在计算损耗时需要求出变压器的内部参数。
(2)谐波引起的变压器降容
谐波电流使得变压器油温或绕组产生过热现象,严重时,甚至会损坏变压器,因此在谐波环境下,为了确保变压器的安全运行,通常会降低变压器的实际使用容量,或在变压器定容时考虑一定的预留容量,从而造成了一定的损失。
变压器的降容率可表示为:
其中:
式中,δt为谐波引起的变压器降容率;St、Vt为变压器的额定容量和额定电压;Sr、Vr为变压器在谐波条件下的运行容量和电压有效值;I*max为最大允许电流标幺值;P*LL为变压器负载损耗标幺值;P*EC为变压器绕组涡流损耗;FHL为谐波损耗因子[IEEE C57.110标准]。
(3)谐波造成的电力设备寿命缩短
谐波不仅加剧了电力设备的损耗,而且使得其温度升高,绝缘受损,最终导致电力设备的使用寿命缩短。
1)谐波电流使得变压器热老化加速、绝缘寿命下降,根据IEEE 57.91标准中的定义,变压器的寿命损失率为:
其中:
式中,LT为变压器在时间T内的寿命损失率;B为常数;θ(t)为t时刻的热点温度;FAA(t) 为t时刻的老化因子,需根据不同的变压器设计标准进行修正。
2)谐波会使电容器电流、电压及容量增大,从而加速了电容器的热老化和电老化,最终导致了电容器的使用寿命缩短。谐波条件下的电容器寿命为:
其中:L0=130000h,a=-0.5,b=-6.1,c=-2.6。式中,Lc为谐波条件下的电容器寿命;L0为正弦波形下标称电压时的电容器寿命;Kf为波形系数;Kp为峰值系数;Krms为均方根系数; a、b、c数值的大小表示三个因素对加速电容器老化的影响程度。
(4)谐波产生的风险
谐波条件下很容易会引发谐波谐振或谐波放大现象,不仅使得电力设备的附加损耗大大增加、寿命缩短加剧,而且会导致电容器组爆炸、电压互感器损毁、开关跳闸、发电厂停机、停炉、异步电机反转等严重事故。
2、谐波治理需求分析
非线性用户参与治理的积极性不足,主要原因有两个:一是非线性用户是造成谐波污染问题的主导因素,但其自身对谐波问题的敏感度以及造成的经济损失不明确,谐波治理需求模糊,从而治理意愿薄弱;二是电能质量监测点大多分布在电网侧,导致用户获得的信息不对等,从而对自身的用电环境了解不足。
为全面、系统地评估非线性用户接受治理的必要性和治理需求,以实现让用户了解自身的用电环境并积极参与治理,需综合考虑多项谐波相关指标,如表1所示。
表1综合评价指标
其中,表1中指标的计算表达式为:
(1)谐振程度:
式中,Ufh为h次谐波谐振情况下的谐波电压;U1为基波电压。
(2)谐波电压畸变率裕度:
式中,THDUGB为谐波电压畸变率限值。
(3)谐波电流含量裕度:
式中,HRIGB为谐波电流含有率限值。
目前,对于电网用电环境的评估方法有层次分析法、模糊综合评判法、基于熵权的层次分析法等,但是传统方法确定指标的权重需依赖人为的经验和主观判断,缺乏系统性,因此需要更加客观、准确、稳定的评估方法。
3、数据包络分析(data envelopment analysis,DEA)
数据包络分析方法是由Charnes,Cooper和Rhodes在1978年提出,其主要采用数学规划模型评价具有多输入多输出的决策单元(decision making units,DMUs)之间的效益或相对有效性,是一种非参数的评估方法。DEA方法允许决策单元选择最有利的权重,以实现最佳的相对效率。DEA方法具有以下优点:1)不需要考虑输入与输出之间的函数关系;2)不需要预先估计参数、任何权重假设,以决策单元各输入输出的权重为变量,从最有利与决策单元的角度进行评价,避免了主观因素的影响;3)最终的评价结果与各输入输出数据的量纲无关,具有很强的客观性。
(1)DEA基本原理
假设有n个评价对象,即n个决策单元,每个决策单元都有m种类型的输入和s种类型的输出。每个DMUj(j=1,2,…,n)有m个输入指标xij(i=1,2,…,m)和s个输出指标yrj(r=1,2,…,s)。决策单元j的效率指数可表示为:
式中,V=(v1,v2,…,vm)T和U=(u1,u2,…,us)T为输入指标和输出指标的权重向量,均为实数向量。权重向量是待求参数,总可以求得合适的权系数使其满足:
Ej≤1,j=1,2,…,n (13)
以所有决策单元的效率指数为约束,第j0个决策单元的效率评价CCR模型可表示为:
式中,E0第j0个决策单元DUM0的相对效率。模型(14)是一个分式规划,通过Charnes-Cooper变换,可将其化为一个等价的线性规划问题。令
则有:
在线性规划(13)的最优解中,满足ω0>0,μ0>0,E0=1,则称DUM0为DEA有效,反之则成为非DEA有效。为有效检验决策单元的有效性,引入非阿基米德无穷小概念构建有效性判断数学模型,则(16)可表示为:
其对偶规划为:
其中
式中,ε为非阿基米德无穷小量;θ为输入相对于输出的有效利用率;s-为输入指标的松弛变量;s+为输出变量的剩余变量;λ为决策单元的组合系数;Em为m阶单位矩阵;Es为s阶单位矩阵。根据(14)可判断决策单元的有效性:若θ=1,s-=0,s+=0,则称为DEA有效;若θ=1,s-≠0或s+≠0,则称为DEA弱有效;若θ<1则称为非DEA有效。
(2)超效率模型
DEA的基本功能就是对多个决策单元的相对优劣性进行评估,但当出现大量或者全部决策单元都有效时,传统的DEA模型无法辨别有效决策单元的优劣。为解决这一问题,提出超效率模型,该模型的效率指标没有小于等于1的约束,对于DEA有效的决策单元,其超效率值可能大于1,而对于非DEA有效的决策单元,其效率值不变,这样就避免了效率值在1处的拥挤现象。超效率的数学模型可表示为:
式中,θ为输入相对于输出的有效利用率;ε为非阿基米德无穷小量;s-为输入指标的松弛变量;s+为输出变量的剩余变量;λj为决策单元的组合系数;n为评价对象的数量,即决策单元的数量;xj为第j个决策单元的输入指标,x0为第j0个决策单元的输入指标;yj为第j 个决策单元的输出指标,y0为第j0个决策单元的输出指标。
然而,传统的评估过程没有只考虑了对每个决策单元最有利的权重,没有将所有可能的权重考虑在内。在应用DEA模型进行评估时,其结果可能会随着权重的不同而发生变化,因此,应考虑所有的可能权重来对决策单元进行评估,由此引入区间超效率模型,以实现评估结果的客观性和准确性。
(3)超效率区间
定义决策单元使用最优权重的交叉超效率值为效率区间上限。DMUd的最大交叉超效率可表示为:
其中,yrj为输出指标;xij为输入指标;s为输出指标数量;m为输入指标数量;μrd、ωid为相关联的输出和输入权重;θd为决策单元DMUd输入相对于输出的有效利用率;xid、yrd为被评估决策单元DMUd的输入指标和输出指标;为决策单元DMUd的效率区间上限;同理,其最小交叉超效率值可表示为:
(4)极大极小后悔值排序方法
由于得到的是每个决策单元的效率区间,因此需要采用一种简单实用的方法对决策单元进行比较和排序。本发明采用一种基于极大极小后悔值法,该方法可以比较和排序决策单元的效率区间。
当E dj≥b时,决策者不会遭受损失,也不会后悔,此时rj=0。结合上述两种情况,第j 个决策单元的最大效率损失为:
因此,依据以下条件选择最佳效率区间:
基于上述分析,基于极大极小后悔值法对区间效率进行比较和排序,步骤如下:
Step1:计算每个DMU的最大效率损失,最小者记作Aj1(1≤j1≤n)。
Step2:将Aj1排除,重新计算n-1个DMU的最大效率损失,最小者记作Aj2(1≤j2≤n,j2≠j1)。
Step3:重复上述步骤,直到只剩下一个DMU。最终的排序为Aj1>Aj2>…>Ajn。
(5)算法应用
采用上述方法从输入与输出指标两个方面评估非线性用户的谐波治理需求,所考虑指标必须要全面反映相关的属性要素,其综合评估体系架构如图1所示。
DEA有效表明谐波治理需求综合评估结果良好;DEA弱有效为一般;非DEA有效为差。最终得到每个决策单元的效率区间和排序,可将决策单元的谐波治理需求分为三个等级,即非常需要、需要和良好,如表2所示。根据实际工程中需要关注的各项指标,该模型可加以扩展或简化,以满足不同的应用场合。
表2谐波治理需求等级
4、具体应用的最佳实施例
某大型工业园区,采用两路220kV供电,申请用电总容量为10000MVA,由5家典型的非线性用户构成,均为典型的6脉动整流谐波源电力用户,供电线路为电缆和铝包钢芯铝绞线,主要用电负荷包括电弧炉、不间断电源、变频器、照明空调等,其园区结构如图2所示。非线性设备固有的谐波特性在运行时会对电网造成谐波污染,根据监测数据和国家标准,该园区不能满足要求,需引导用户参与谐波治理,定义决策单元集为DMUj(j=1,2,3,4,5),该园区用户的指标数据如表3所示。
表3综合评价指标数据
采用所提方法对表3数据进行评估,输出指标求倒数转化为正指标,评估结果如表4所示,各决策单元的效率间隔如图3所示。
表4综合评价指标数据
从图表中可以看出,各非线性用户CCR模型的效率值为(1,0.6612,1,1,1),决策单元1、3、 4、5至少是弱有效的,DMU2是非DEA有效。为全面评估决策单元1、3、4的谐波治理需求,采用超效率模型进行综合评估,考虑其所有可能权重,各非线性用户超效率模型的最大效率值为(2.7857,0.6612,7.7879,10.106,1.5304),最小效率值为(1.2289,0.4156,2.0126,2.8212,0.6642),可以看出非DEA有效用户的计算结果与CCR模型相同,这样就避免了效率值在1处的拥挤情况,且每个决策单元的效率会随着权重的变化而变化,其中DMU4的影响最为明显。
采用所提方法对决策单元进行比较和排序,首先计算各个决策单元的最大损失效率:
R(DMU1)=max{max(0.6612,7.7879,10.106,1.5304)-1.2289,0}=8.877;
R(DMU2)=max{max(2.7857,7.7879,10.106,1.5304)-0.4156,0}=9.690;
R(DMU3)=max{max(2.7857,0.6612,10.106,1.5304)-2.0126,0}=8.093;
R(DMU4)=max{max(2.7857,0.6612,7.7879,1.5304)-2.8212,0}=4.967;
R(DMU5)=max{max(2.7857,0.6612,7.7879,10.106)-0.6642,0}=9.442.
显然,DMU4的最大损失效率最小,因此DMU4为最佳决策单元,将其排除,重新计算剩余决策单元的最大损失效率:
R(DMU1)=max{max(0.6612,7.7879,1.5304)-1.2289,0}=6.559;
R(DMU2)=max{max(2.7857,7.7879,1.5304)-0.4156,0}=7.372;
R(DMU3)=max{max(2.7857,0.6612,1.5304)-2.0126,0}=0.773;
R(DMU5)=max{max(2.7857,0.6612,7.7879)-0.6642,0}=7.124.
同理,DMU3的最大损失效率最小,因此DMU4为第二佳决策单元,将其排除,重新计算剩余决策单元的最大损失效率:
R(DMU1)=max{max(0.6612,1.5304)-1.2289,0}=0.302;
R(DMU2)=max{max(2.7857,1.5304)-0.4156,0}=2.370;
R(DMU5)=max{max(2.7857,0.6612)-0.6642,0}=2.122.
同理,DMU1的最大损失效率最小,因此DMU1为第三佳决策单元。重复上述过程,最终得到决策单元的排名顺序为:DMU4>DMU3>DMU1>DMU5>DMU2。
采用本方法可以更加直观的看出各个用户的谐波治理需求,其中,相较于其他用户,用户4的谐波用电环境最佳。该结果全面、系统、客观的评估非线性用户的用电环境以及谐波治理需求,可提升用户的治理积极性,并对后续谐波治理的实施起到指导作用。
Claims (5)
1.一种基于数据包络分析的谐波治理需求评估方法,其特征在于,包括以下步骤:
步骤1:进行谐波损失分析和谐波治理需求分析,确定谐波治理需求综合评价指标集,包括输入指标和输出指标;
步骤2:采用DEA方法从输入指标与输出指标两个方面评估非线性用户和园区的谐波治理需求,对多个决策单元的相对优劣性进行评估;
步骤3:当出现大量或者全部决策单元都有效时,需考虑所有的可能权重来对决策单元进行评估,由此引入区间超效率模型;
步骤4:由于得到的是每个决策单元的效率区间,故采用基于极大极小后悔值法对决策单元进行比较和排序;
所述超效率模型表示为:
Min[θ-ε(e^Ts-+eTs+)]
λj≥0,j=1,2,…,n
s-≥0,s+≥0 (12)
式中,θ为输入相对于输出的有效利用率;ε为非阿基米德无穷小量;s-为输入指标的松弛变量;s+为输出变量的剩余变量;λj为决策单元的组合系数;n为评价对象的数量,即决策单元的数量;xj为第j个决策单元的输入指标,x0为第j0个决策单元的输入指标;yj为第j个决策单元的输出指标,y0为第j0个决策单元的输出指标;e^T=(1,1,…,1)∈Em,eT=(1,1,…,1)∈Es;Em为m阶单位矩阵;Es为s阶单位矩阵;
若θ=1,s-=0,s+=0,则称为DEA有效;若θ=1,s-≠0或s+≠0,则称为DEA弱有效;若θ<1则称为非DEA有效;效率模型的效率指标没有小于等于1的约束,对于DEA有效的决策单元,其超效率值大于1,而对于非DEA有效的决策单元,其效率值不变;DEA有效表明谐波治理需求综合评估结果良好,非DEA有效为表明谐波治理需求综合评估结果差;
定义决策单元使用最优权重的交叉超效率值为效率区间上限;第d个决策单元DMUd的最大交叉超效率表示为:
其中,yrj为输出指标;xij为输入指标;s为输出指标数量;m为输入指标数量;μrd、ωid为相关联的输出和输入权重;θd为决策单元DMUd输入相对于输出的有效利用率;xid、yrd为被评估决策单元DMUd的输入指标和输出指标;为决策单元DMUd的效率区间上限;同理,其最小交叉超效率值表示为:
其中,E d j为决策单元DMUd的效率区间下限;
其中,rj为决策单元DMUj的最大损失效率;b为除自身外,其余决策单元的最大效率区间上限;
当E dj≥b时,决策者不会遭受损失,也不会后悔,此时rj=0;
结合上述两种情况,第j个决策单元的最大效率损失为:
因此,依据以下条件选择最佳效率区间:
2.根据权利要求1所述的基于数据包络分析的谐波治理需求评估方法,其特征在于,所述进行谐波损失分析包括:
1)分析谐波引起的变压器和输电线路附加损耗:
输电线路上的附加损耗分为基波电流引起的损耗和谐波电流引起的损耗,而谐波损耗为基波损耗的kl倍;计算输电线路的谐波损耗:
其中:
式中,ΔPl为谐波引起的输电线路附加损耗;h为谐波次数;l为最大谐波次数;R1为基波下的输电线路电阻;I1为基波电流;HRIh为h次谐波电流含有率,是谐波电流与基波电流的比率;Rh为h次谐波电阻;
变压器的损耗分为铜耗与铁耗,因谐波所产生的变压器附加损耗表示为:
式中,ΔPt为谐波引起的变压器附加损耗;UAh、UBh、UCh为A、B、C各相第h次谐波电压;IAh、IBh、ICh为A、B、C各相第h次谐波电流;为A、B、C各相第h次谐波电压与谐波电流之间的相位差;
2)分析谐波引起的变压器降容
计算变压器的降容率:
其中:
式中,δt为谐波引起的变压器降容率;St、Vt为变压器的额定容量和额定电压;Sr、Vr为变压器在谐波条件下的运行容量和电压有效值;I*max为最大允许电流标幺值;P*LL为变压器负载损耗标幺值;P*EC为变压器绕组涡流损耗;FHL为谐波损耗因子;
3)分析谐波造成的电力设备寿命缩短
计算变压器的寿命损失率:
其中:
式中,LT为变压器在时间T内的寿命损失率;B为常数;θ(t)为t时刻的热点温度;FAA(t)为t时刻的老化因子,需根据不同的变压器设计标准进行修正;
计算谐波条件下的电容器寿命:
式中,Lc为谐波条件下的电容器寿命;L0为正弦波形下标称电压时的电容器寿命;Kf为波形系数;Kp为峰值系数;Krms为均方根系数;a、b、c数值的大小表示三个因素对加速电容器老化的影响程度。
3.根据权利要求1所述的基于数据包络分析的谐波治理需求评估方法,其特征在于,所述输入指标包括:
内部参数:变压器容量、输电线路长度、电容器容量;
负荷参数:负荷总容量、谐波源数量、谐波源容量和谐波源总谐波发射特性;
所述输出指标包括:
设备损耗:变压器附加损耗、输电线路附加损耗、变压器降容损耗、变压器寿命损耗、电容器寿命损耗;
国家标准:谐波电流含量裕度、谐波电压畸变率裕度;
潜在风险:谐振峰值数量、谐振程度。
5.根据权利要求1所述的基于数据包络分析的谐波治理需求评估方法,其特征在于,所述采用基于极大极小后悔值法对决策单元进行比较和排序,步骤如下:
Step1:计算每个DMU的最大效率损失,最小者记作Aj1,1≤j1≤n;
Step2:将Aj1排除,重新计算n-1个DMU的最大效率损失,最小者记作Aj2,1≤j2≤n,j2≠j1;
Step3:重复上述步骤,直到只剩下一个DMU;最终的排序为Aj1>Aj2>…>Ajn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210347350.XA CN114677031B (zh) | 2022-04-01 | 2022-04-01 | 一种基于数据包络分析的谐波治理需求评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210347350.XA CN114677031B (zh) | 2022-04-01 | 2022-04-01 | 一种基于数据包络分析的谐波治理需求评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114677031A CN114677031A (zh) | 2022-06-28 |
CN114677031B true CN114677031B (zh) | 2022-11-25 |
Family
ID=82077184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210347350.XA Active CN114677031B (zh) | 2022-04-01 | 2022-04-01 | 一种基于数据包络分析的谐波治理需求评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114677031B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107203842A (zh) * | 2017-05-18 | 2017-09-26 | 西南交通大学 | 基于扩展云相似度与逼近理想解的谐波污染水平评估方法 |
CN112990695A (zh) * | 2021-03-11 | 2021-06-18 | 安徽大学 | 基于博弈dea-多目标决策的dg电能质量评估方法 |
CN112990710A (zh) * | 2021-03-18 | 2021-06-18 | 江苏方天电力技术有限公司 | 一种基于dea的风电场电能质量多层次综合评估方法 |
-
2022
- 2022-04-01 CN CN202210347350.XA patent/CN114677031B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107203842A (zh) * | 2017-05-18 | 2017-09-26 | 西南交通大学 | 基于扩展云相似度与逼近理想解的谐波污染水平评估方法 |
CN112990695A (zh) * | 2021-03-11 | 2021-06-18 | 安徽大学 | 基于博弈dea-多目标决策的dg电能质量评估方法 |
CN112990710A (zh) * | 2021-03-18 | 2021-06-18 | 江苏方天电力技术有限公司 | 一种基于dea的风电场电能质量多层次综合评估方法 |
Non-Patent Citations (1)
Title |
---|
电能质量治理设备运行状态识别及其治理效果评价;甄超等;《电测与仪表》;20200210(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114677031A (zh) | 2022-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109325694B (zh) | 基于承载能力的配电网优选方法 | |
Kansal et al. | Optimal placement of distributed generation in distribution networks | |
CN111861248B (zh) | 一种配电网的电能质量治理效果综合评估方法及装置 | |
CN114677031B (zh) | 一种基于数据包络分析的谐波治理需求评估方法 | |
Hamzeh et al. | Reliability evaluation of distribution transformers considering the negative and positive effects of rooftop photovoltaics | |
CN104578044A (zh) | 一种基于次同步振荡分析的确定机组作用系数的方法 | |
Angelim et al. | Probabilistic impact assessment of electric vehicles charging on low voltage distribution systems | |
Moreira et al. | Disturbing load classification based on the grey relational analysis method and load performance index | |
CN112016838B (zh) | 配电网能效指标体系的贡献率计算方法、系统及终端设备 | |
CN108229786B (zh) | 电网投资分配处理方法及装置 | |
Yin et al. | Global configuration strategy of SAPFs combined with comprehensive evaluation of power quality in distribution network | |
Menti et al. | Optimal sizing and limitations of passive filters in the presence of background harmonic distortion | |
Rastvorova | Assessment of the consumers’ contribution to the deterioration of the electrical power quality | |
Dalali et al. | Modified Thevenin‐based voltage instability indicator and load shedding approach for MCF connected network | |
CN111178634A (zh) | 一种配电网建设改造方案选定方法及系统 | |
Wang et al. | Evaluation Methods for the Development of New Power Systems Based on Cloud Model | |
CN113191675A (zh) | 多直流送端电网规划方案适应性评估方法及系统 | |
Urooj et al. | An Adaptive Neuro-Fuzzy based Methodology for Harmonic Analysis of a Power Transformer. | |
Ma et al. | Harmonie evaluation of grid with multiple harmonic sources based on DIgSILENT | |
Zou et al. | Discussion on Operation State Evaluation of Low Voltage Distribution Network | |
Sonwane et al. | Algorithms and methodology for development of OCP-PSO tool for optimal capacitor placement to enhance system reliability | |
Akbal et al. | PSO and CSA to Estimate of Parameter in Power Line | |
Zhang et al. | Comprehensive Evaluation and Improvement of Power Quality in Distribution Networks Based on DEA Model and Membership Function | |
Swapna et al. | Reactive power control in distribution network by optimal location and sizing of capacitor using Fuzzy and SFLA | |
Liu et al. | Evaluation of transformer state based on the life cycle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |