CN112085156A - 一种分压线损率预测及降损潜力评估方法及*** - Google Patents

一种分压线损率预测及降损潜力评估方法及*** Download PDF

Info

Publication number
CN112085156A
CN112085156A CN202010682242.9A CN202010682242A CN112085156A CN 112085156 A CN112085156 A CN 112085156A CN 202010682242 A CN202010682242 A CN 202010682242A CN 112085156 A CN112085156 A CN 112085156A
Authority
CN
China
Prior art keywords
loss rate
line loss
line
sample
input
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010682242.9A
Other languages
English (en)
Other versions
CN112085156B (zh
Inventor
常征
胡鹏
汤恒
马瑞
吴龙雨
赵爽
苏同厚
曹冬梅
孙智慧
李奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaiyuan Power Supply Co Of State Grid Anhui Electric Power Co ltd
Bengbu Power Supply Co of State Grid Anhui Electric Power Co Ltd
Original Assignee
Huaiyuan Power Supply Co Of State Grid Anhui Electric Power Co ltd
Bengbu Power Supply Co of State Grid Anhui Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaiyuan Power Supply Co Of State Grid Anhui Electric Power Co ltd, Bengbu Power Supply Co of State Grid Anhui Electric Power Co Ltd filed Critical Huaiyuan Power Supply Co Of State Grid Anhui Electric Power Co ltd
Priority to CN202010682242.9A priority Critical patent/CN112085156B/zh
Publication of CN112085156A publication Critical patent/CN112085156A/zh
Application granted granted Critical
Publication of CN112085156B publication Critical patent/CN112085156B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06395Quality analysis or management
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Molecular Biology (AREA)
  • Evolutionary Computation (AREA)
  • Educational Administration (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种分压线损率预测及降损潜力评估方法及***,包括以下步骤:获取电网的运行数据,包括发售电量统计、输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指标体系;采用箱线图排除原始数据中的异常值,并通过多月累加的方法得到新样本;把处理过得到的新样本数据,输入事先训练好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降损潜力进行评估。本发明提出的分压线损率预测及降损潜力评估方法及***,根据电压等级的不同,量化线路长度、变压器容量等因素对线损率的影响程度,充分挖掘电网的降损潜力,为电网未来的降损工作提出了合理建议和科学理论指导。

Description

一种分压线损率预测及降损潜力评估方法及***
技术领域
本发明涉及电力***技术领域,具体涉及一种分压线损率预测及 降损潜力评估方法及***。
背景技术
节能降损工作是响应国家建设节约型社会的具体行动,也是电网 公司提高企业经营效益的内在要求。节能降损有助于合理利用能源, 提高能源利用效率和经济效率,降低能源损耗,保护环境,是电网公 司社会责任和企业责任的集中体现。电网降损规划是在建设网架坚强、 电网协调发展、结构优化、电力供应充足的同时,合理的采取各种切 实可行的有效的节能降损措施,尽可能的减小线损,降低损耗,提升 电网供电可靠性及电力企业的经验效益。
然而,理论线损计算以大负荷日实测的数据为基础(SCADA和用 采***等数据),基于潮流计算和等值电阻法等来计算真实的线路损耗 和变压器损耗,工作量巨大且计算值与统计线损存在偏差,同时线损 率与各类影响因素间存在复杂的非线性关系,不同因素对线损率的影 响程度各不相同,定性的分析不足以指导实际的线损管理工作,研究 如何确定影响线损率的关键因素并量化其降损潜力,对于提高电网节 能降损工作的效率及经济性具有重要意义。考虑到线损率与及影响因 素间的关系难以通过解析的方式进行描述,有研究将回归分析法应用 于低压台区线损率的计算,但是回归方程难以确定,计算准确度低。
近年来,机器学习理论的发展与应用为线损率计算提供了新的途 径。采用机器学习的方法无须建立反映实际物理机理的数学模型,可 以利用机器学习模型强大的自学能力、推广能力以及非线性处理能力 来挖掘线损率与特征参数之间的关系。
发明内容
本发明提出的一种分压线损率预测及降损潜力评估方法及***, 可解决上述背景技术的不足。
为实现上述目的,本发明采用了以下技术方案:
一种分压线损率预测及降损潜力评估方法,包括以下步骤,
S100、获取电网的运行数据,包括发售电量统计、输电线路参数 和变压器参数,建立不同电压等级的线损率影响因素指标体系;
S200、采用箱线图排除原始数据中的异常值,并通过多月累加的 方法得到新样本;
S300、把处理过得到的新样本数据,输入事先训练好的线损率预 测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降 损潜力进行评估。
进一步的,还包括S400、改变输入模型数据的取值范围,进行灵 敏度分析,评估各降损措施的降损潜力。
进一步的,所述S100、获取电网的运行数据,包括发售电量统计、 输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指 标体系;
包括:
建立不同电压等级的线损率影响因素指标体系,以国家电网的“规 划计划信息管理平台”样本数据的月度线损率报表为基础,选用的原始 表格及关键特征包括:
“公司分压线损统计表”:电厂上网电量、转入电量、转出外省电 量、售电量;
“35kV及以上线路情况”:线路总条数、线路总长度、20年及以上 线路条数及长度;
“变压器情况”:变压器总台数、变压器总容量、20年及以上变压 器台数及容量;
“新能源电厂并网运行情况表”:风、光以及其他新能源的上网电 量;
通过“公司分压线损统计表”计算线损率y表示为:
Figure BDA0002586270710000031
式中,Q1为电厂上网电量;Q2为转入电量;Q3为转出省外电量;Q4为售电量。
进一步的,所述S200、采用箱线图排除原始数据中的异常值,并 通过多月累加的方法得到新样本;
包括:
采用箱线图排除原始数据中的异常值,箱线图定义了异常值的数 值范围是[下四位数-1.5*四分位数间距,上四分位数+1.5*四分位数间 距],对检测出的异常数据进行剔除;
采用多月累加的方式对原始样本进行数据处理,从而得到新样本。
进一步的,S300所述线损率预测模型的训练步骤如下:
根据输入***的样本,以实际线损率与预测线损率的误差最小为 目标训练BP神经网络的连接权值、偏置相关网络参数,具体如下:
建立训练集D={(x1,y1),(x2,y2),...,(xm,ym)},
Figure BDA0002586270710000032
测试集 共有m个样本,每个样本的输入有d个特征,输出为1维实值向量; 输出层神经元的阈值为θ={θ1,θ2,...,θl},隐藏层的阈值为 δ={δ1,δ2,...,δq},输入层第i个神经元与隐藏层第h个神经元的连接权 值为vih,隐藏层第h个神经元与输出层第j个神经元之间的连接权值 为whj,隐藏层和输出层神经元都选用sigmoid函数作为传递函数,表 示为:
Figure BDA0002586270710000033
对于第k个输入样本(xk,yk):
隐藏层第h个神经元接收到的净输入ak h和输出Ak h分别为:
Figure BDA0002586270710000041
Figure BDA0002586270710000042
输出层第j个神经元从隐藏层接收到的净输入和输出分别为:
Figure BDA0002586270710000043
Figure BDA0002586270710000044
该神经网络在输入样本(xk,yk)上的均方误差为:
Figure BDA0002586270710000045
通过误差反向传播算法采用梯度下降的策略通过不断迭代对网络 进行训练,使得网络输出与期望的输出间的误差Ek最小,以隐藏层到 输出层的连接权值whj为例进行具体推导:
以目标函数的负梯度方向对参数进行调整,目标函数是网络输出 与期望输出的误差,即:
wh,j←wh,j-ηΔwh,j
Figure BDA0002586270710000046
式中,η表示学习率,取值范围为0到1,用于控制算法每一轮迭 代中的更新步长,若太大则会发生振荡,太小则影响训练的收敛速度;
根据链式求导法则,误差Ek对whj的梯度表示为:
Figure BDA0002586270710000047
由于
Figure BDA0002586270710000048
Sigmoid函数作为传递函数满足:
f′(x)=f(x)(1-f(x)) (11)
得到
Figure BDA0002586270710000051
最后得到参数whj的更新公式为:
Figure BDA0002586270710000052
同理得到输入层到隐藏层权值以及隐藏层、输出层偏置的更新公 式:
Figure BDA0002586270710000053
Figure BDA0002586270710000054
Figure BDA0002586270710000055
对于容量为m的样本集D={(x1,y1),(x2,y2),...,(xm,ym)},对应的误 差为整个样本集中各样本误差的和:
Figure BDA0002586270710000056
各参数更新时的负梯度为所有样本负梯度的和,即:
Figure BDA0002586270710000057
Figure BDA0002586270710000058
Figure BDA0002586270710000059
Figure BDA00025862707100000510
选取实际线损率与BP神经网络输出预测线损率的平均绝对百分 误差eMAPE作为评价预测效果的依据:
Figure BDA0002586270710000061
式中,
Figure BDA0002586270710000062
和yi分别表示测试集中第i个样本的预测线损率和实际线 损率,M表示测试集中样本的总数;
一次训练集和测试集划分得到的性能指标可能具有特异性,为此 这里重复20次实验,每次随机选取训练集和测试集,将20次实验的 eMAPE的平均值作为预测模型的性能评价指标。
进一步的,所述S300、把处理过得到的新样本数据,输入事先训 练好的线损率预测模型中,得到电网分压线损率的预测值,即对分压 线损率预测及降损潜力进行评估;
包括:
假设各电压等级与外省的电量输送情况基本不变,即转入电量和 转出外省电量不变,电厂上网电量的年增长率为η,则该电压等级的负 荷功率Pload为:
Figure BDA0002586270710000063
式中,Q1为电厂上网电量;Q2为转入电量;Ymax为年最大负荷利用 小时数;根据电网投资期望容载比和负荷功率确定新建变电站的总容 量Pnew,并假设变电站的平均容量基本不变,进而确定新建变电站的 座数N1
Figure BDA0002586270710000064
Pnew=N1·Pmean (24)
式中,K为期望容载比;Plast为上一年变电站容量;Pmean为变电站 平均容量;[.]为向上取整函数;
假设各年每回线路输送的平均功率及该电压等级线路的平均长度 保持不变,进而确定新建线路的条数N2及长度lnew
Figure BDA0002586270710000071
lnew=N2·lmean (26)
式中,Pline为单回线路平均输送功率;lmean为线路平均长度;
将预测的上网电量、线路条数、变电站容量等参数输入分压线损 率预测模型中,得到电网未来的分压线损率。
进一步的,所述S400、改变输入模型数据的取值范围,进行灵敏 度分析,评估各降损措施的降损潜力;包括:
基于灵敏度分析方法给出各电压等级线损率的关键影响因素,针 对模型y=f(x1,x2,...,xi),xi为模型的第i个输入属性值,令每个输入属 性值在可能的取值范围内变动,确定这些输入属性的变动对模型输出 值的影响程度。
另一方面,本发明还公开一种分压线损率预测及降损潜力评估系 统,包括以下单元:
指标体系单元,用于获取电网的运行数据,包括发售电量统计、 输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指 标体系;
新样本处理单元,用于采用箱线图排除原始数据中的异常值,并 通过多月累加的方法得到新样本;
预测评估单元,用于把处理过得到的新样本数据,输入事先训练 好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线 损率预测及降损潜力进行评估。
进一步的,还包括:
灵敏度分析单元,用于改变输入模型数据的取值范围,进行灵敏 度分析,评估各降损措施的降损潜力。
进一步的,还包括以下子单元:
线损率预测模型的训练单元,用于根据输入***的样本,以实际 线损率与预测线损率的误差最小为目标训练BP神经网络的连接权值、 偏置相关网络参数。
由上述技术方案可知,本发明的分压线损率预测及降损潜力评估 方法及***,针对当前电网输入电量与输出电量的结算不同步导致线 损计算不准确,线损率与及影响因素间的关系难以建模,电网的降损 潜力难以定量评估等一系列问题,本申请的一种分压线损率预测及降 损潜力评估方法,通过多月累加的数据处理方法,解决了输入电量与 输出电量的结算不同步的问题;通过建立基于BP神经网络的线损率预 测模型,无须建立反映实际物理机理的数学模型,解决了线损率与及 影响因素间的关系难以通过解析的方式进行描述的技术难题;通过预 测模型对影响线损率的各项因素进行灵敏度分析,解决了电网的降损 潜力难以定量评估的技术难题。
本发明提出的分压线损率预测及降损潜力评估方法及***,根据 电压等级的不同,量化线路长度、变压器容量等因素对线损率的影响 程度,充分挖掘电网的降损潜力,为电网未来的降损工作提出了合理 建议和科学理论指导。
附图说明
图1为本发明方法的总流程图;
图2为分压线损率影响指标体系图;
图3为同期***月度分压电量传导关系图;
图4为500kV原始数据箱线图;
图5为多月累加数据处理原理图;
图6为BP神经网络结构图;
图7为基于BP神经网络的线损率预测模型参数训练图;
图8为本实施例500kV线损率预测结果;
图9为分压线损率预测结果;
图10为500kV线损率随各影响因素变化情况。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结 合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、 完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是 全部的实施例。
如图1所示,本实施例所述的分压线损率预测及降损潜力评估方 法,包括:
针对当前电网输入电量与输出电量的结算不同步导致线损计算不 准确,线损率与及影响因素间的关系难以建模,电网的降损潜力难以 定量评估等一系列问题,本发明提出了一种分压线损率预测及降损潜 力评估方法,该方法如图1所示,主要包括如下4个步骤:
S100、获取电网的运行数据,包括发售电量统计、输电线路参数 和变压器参数,建立不同电压等级的线损率影响因素指标体系;
S200、采用箱线图排除原始数据中的异常值,并通过多月累加的 方法得到新样本;
S300、把处理过得到的新样本数据,输入事先训练好的线损率预 测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降 损潜力进行评估。
还包括S400、改变输入模型数据的取值范围,进行灵敏度分析, 评估各降损措施的降损潜力。
以下针对上述步骤,具体说明:
1)建立不同电压等级的线损率影响因素指标体系,采用箱线图排 除原始数据中的异常值,采用多月累加的方法得到新样本;
建立不同电压等级的线损率影响因素指标体系,以国家电网的“规 划计划信息管理平台”样本数据的月度线损率报表为基础,选用的原始 表格及关键特征包括:
“公司分压线损统计表”:电厂上网电量、转入电量、转出外省电 量、售电量;
“35kV及以上线路情况”:线路总条数、线路总长度、20年及以上 线路条数及长度;
“变压器情况”:变压器总台数、变压器总容量、20年及以上变压 器台数及容量;
“新能源电厂并网运行情况表”:风、光以及其他新能源的上网电 量。
通过“公司分压线损统计表”计算线损率y表示为:
Figure BDA0002586270710000101
式中,Q1为电厂上网电量;Q2为转入电量;Q3为转出省外电量;Q4为售电量。
采用箱线图排除原始数据中的异常值,箱线图定义了异常值的数 值范围是[下四位数-1.5*四分位数间距,上四分位数+1.5*四分位数间 距],对检测出的异常数据进行剔除。
为了减小由于抄表不同期带来的误差,采用多月累加的方式对原 始样本进行数据处理,从而得到新样本,除首、尾两月外,其他各月供、 售电抄表不同步对线损率虚增或虚降的影响相互抵消。
2)根据输入***的样本,以实际线损率与预测线损率的误差最小 为目标训练BP神经网络的连接权值、偏置等网络参数;
建立训练集D={(x1,y1),(x2,y2),...,(xm,ym)},
Figure BDA0002586270710000102
测试集 共有m个样本,每个样本的输入有d个特征,输出为1维实值向量。 输出层神经元的阈值为θ={θ1,θ2,...,θl},隐藏层的阈值为 δ={δ1,δ2,...,δq},输入层第i个神经元与隐藏层第h个神经元的连接权 值为vih,隐藏层第h个神经元与输出层第j个神经元之间的连接权值 为whj,隐藏层和输出层神经元都选用sigmoid函数作为传递函数,表 示为
Figure BDA0002586270710000111
对于第k个输入样本(xk,yk):
隐藏层第h个神经元接收到的净输入
Figure BDA0002586270710000112
和输出
Figure BDA0002586270710000113
分别为:
Figure BDA0002586270710000114
Figure BDA0002586270710000115
输出层第j个神经元从隐藏层接收到的净输入和输出分别为:
Figure BDA0002586270710000116
Figure BDA0002586270710000117
该神经网络在输入样本(xk,yk)上的均方误差为:
Figure BDA0002586270710000118
通过误差反向传播算法采用梯度下降的策略通过不断迭代对网络 进行训练,使得网络输出与期望的输出间的误差Ek最小,以隐藏层到 输出层的连接权值whj为例进行具体推导:
以目标函数的负梯度方向对参数进行调整,目标函数是网络输出 与期望输出的误差,即:
wh,j←wh,j-ηΔwh,j
Figure BDA0002586270710000119
式中,η表示学习率,取值范围为0到1,用于控制算法每一轮迭 代中的更新步长,若太大则会发生振荡,太小则影响训练的收敛速度。
根据链式求导法则,误差Ek对whj的梯度可以表示为:
Figure BDA0002586270710000121
由于
Figure BDA0002586270710000122
Sigmoid函数作为传递函数满足:
f′(x)=f(x)(1-f(x)) (11)
得到
Figure BDA0002586270710000123
最后得到参数whj的更新公式为:
Figure BDA0002586270710000124
同理得到输入层到隐藏层权值以及隐藏层、输出层偏置的更新公 式:
Figure BDA0002586270710000125
Figure BDA0002586270710000126
Figure BDA0002586270710000127
对于容量为m的样本集D={(x1,y1),(x2,y2),...,(xm,ym)},对应的误 差为整个样本集中各样本误差的和:
Figure BDA0002586270710000128
各参数更新时的负梯度为所有样本负梯度的和,即:
Figure BDA0002586270710000129
Figure BDA0002586270710000131
Figure BDA0002586270710000132
Figure BDA0002586270710000133
选取实际线损率与BP神经网络输出预测线损率的平均绝对百分误 差eMAPE作为评价预测效果的依据:
Figure BDA0002586270710000134
式中,
Figure BDA0002586270710000135
和yi分别表示测试集中第i个样本的预测线损率和实际线损 率,M表示测试集中样本的总数。
一次训练集和测试集划分得到的性能指标可能具有特异性,为此这 里重复20次实验,每次随机选取训练集和测试集,将20次实验的eMAPE的平均值作为预测模型的性能评价指标。
3)对于输入线损率预测模型的影响因素进行灵敏度分析,预测未 来上网电量、线路数量、变电站容量等参数,输入模型得到未来线损 率的预测值;
基于灵敏度分析方法给出各电压等级线损率的关键影响因素,针 对模型y=f(x1,x2,...,xi)(xi为模型的第i个输入属性值),令每个输入 属性值在可能的取值范围内变动,研究这些输入属性的变动对模型输 出值的影响程度。
假设各电压等级与外省的电量输送情况基本不变,即转入电量和 转出外省电量不变,电厂上网电量的年增长率为η,则该电压等级的负 荷功率Pload为:
Figure BDA0002586270710000136
式中,Q1为电厂上网电量;Q2为转入电量;Tmax为年最大负荷利用 小时数。根据电网投资期望容载比和负荷功率可以确定新建变电站的 总容量Pnew,并假设变电站的平均容量基本不变,进而可以确定新建 变电站的座数N1
Figure BDA0002586270710000141
Pnew=N1·Pmean (24)
式中,K为期望容载比;Plast为上一年变电站容量;Pmean为变电站 平均容量;[.]为向上取整函数。
假设各年每回线路输送的平均功率及该电压等级线路的平均长度 保持不变,进而可以确定新建线路的条数N2及长度lnew
Figure BDA0002586270710000142
lnew=N2·lmean (26)
式中,Pline为单回线路平均输送功率;lmean为线路平均长度。
将预测的上网电量、线路条数、变电站容量等参数输入分压线损 率预测模型中,得到电网未来的分压线损率。
以下举例说明本发明实施例的具体应用:
以某地区对未来五年的分压线损率进行预测;
对于500kV主网,主要考虑电厂上网电量、转入电量、转出外省 电量以及电网设备等因素对该电压等级线损率的影响。以2016年为基 准年,社会总用电量的年增长率为5%,年最大负荷利用小时数取4000 小时,实际线损率、BP神经网络的输入变量和输出线损率由表1给出。
表1 2016年线损分析参数表
Figure BDA0002586270710000143
Figure BDA0002586270710000151
(1)模型参数优化
确定最优BP神经网络结构,以2016年的数据进行模型的训练, 表2给出了不同网络结构下的线损率预测性能。
表2不同网络结构下500kV线损率预测模型性能指标
Figure BDA0002586270710000152
根据表2,500kV线损预测模型的最优结构为四层,其中第1和第 2隐藏层的神经元数分别为15和5。在该结构下,选取一组训练集和 测试集,得到如图8所示的计算结果。
如图8所示,测试集共有15个样本,预测模型的平均百分误差为 0.8098%,说明该模型能够很好的根据上网电量、输入电量、转出外省 电量以及电网结构等数据对线损率做出预测。
采用与上述500kV电压等级相同的分析计算方法可以确定220kV、 110kV、35kV、10kV及以下各电压等级线损率预测模型的最佳BP神 经网络结构并计算其性能指标,如表3所示。
表3各电压等级最优预测模型
Figure BDA0002586270710000161
(2)分压线损率预测
根据上述方法可以获取未来10年各电压等级的BP神经网络预测模 型的输入数据,输入预测模型后,得到如图9所示的未来10年分压线 损率预测结果。
根据图9,与2016年相比,2022年的线损率将下降0.047%;220kV的线损率 呈现缓慢上升趋势,将达到0.94%且趋于平稳;110kV及35kV的线损率则呈现 先上升后下降的趋势,这主要受新能源渗透率变化的影响,前期新能源受负荷 需求和网架结构的影响难以全部就地消纳,存在倒送现象,线损将有所升高。 后期随着电网公司的投资建设,网架结构更趋合理,运行方式更加灵活,另一 方面,用电负荷的增长幅度大于新能源渗透率的增长幅度,110kV和35kV的新 能源就地消纳能力增强,有助于降损。
(3)灵敏度分析和降损潜力评估
基于前文建立的BP神经网络模型,保持其他参数不变,分别将表 中的电厂上网电量、转入电量、转出外省电量、线路总条数、20年及 以上线路长度、变压器容量以及20年及以上变压器容量调高 0.5%~10%,再次输入模型得到相应的线损率变化情况如图10所示。
以500kV为例可以看出各输入属性对线损率降损潜力的影响程度 从高到底依次为:转入电量、上网电量、转出外省电量、20年及以上 变压器容量、20年及以上线路长度、平均线路长度、变压器总容量。
另一方面,本发明还公开一种分压线损率预测及降损潜力评估系 统,包括以下单元:
指标体系单元,用于获取电网的运行数据,包括发售电量统计、 输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指 标体系;
新样本处理单元,用于采用箱线图排除原始数据中的异常值,并 通过多月累加的方法得到新样本;
预测评估单元,用于把处理过得到的新样本数据,输入事先训练 好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线 损率预测及降损潜力进行评估。
进一步的,还包括:
灵敏度分析单元,用于改变输入模型数据的取值范围,进行灵敏 度分析,评估各降损措施的降损潜力。
进一步的,还包括以下子单元:
线损率预测模型的训练单元,用于根据输入***的样本,以实际 线损率与预测线损率的误差最小为目标训练BP神经网络的连接权值、 偏置相关网络参数。
可理解的是,本发明实施例提供的***与本发明实施例提供的方 法相对应,相关内容的解释、举例和有益效果可以参考上述方法中的 相应部分。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系 统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全 软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请 可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存 储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实 施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机 程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指 令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和 /或方框图中的流程和/或方框的结合。可提供这些计算机程序指令 到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设 备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理 设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程 和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据 处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算 机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现 在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指 定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设 备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生 计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提 供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多 个方框中指定的功能的步骤。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管 参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员 应当理解:其依然可以对前述各实施例所记载的技术方案进行修改, 或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不 使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种分压线损率预测及降损潜力评估方法,其特征在于:包括以下步骤,
S100、获取电网的运行数据,包括发售电量统计、输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指标体系;
S200、采用箱线图排除原始数据中的异常值,并通过多月累加的方法得到新样本;
S300、把处理过得到的新样本数据,输入事先训练好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降损潜力进行评估。
2.根据权利要求1所述的分压线损率预测及降损潜力评估方法,其特征在于:还包括S400、改变输入模型数据的取值范围,进行灵敏度分析,评估各降损措施的降损潜力。
3.根据权利要求1所述的分压线损率预测及降损潜力评估方法,其特征在于:所述S100、获取电网的运行数据,包括发售电量统计、输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指标体系;
包括:
建立不同电压等级的线损率影响因素指标体系,以国家电网的“规划计划信息管理平台”样本数据的月度线损率报表为基础,选用的原始表格及关键特征包括:
“公司分压线损统计表”:电厂上网电量、转入电量、转出外省电量、售电量;
“35kV及以上线路情况”:线路总条数、线路总长度、20年及以上线路条数及长度;
“变压器情况”:变压器总台数、变压器总容量、20年及以上变压器台数及容量;
“新能源电厂并网运行情况表”:风、光以及其他新能源的上网电量;
通过“公司分压线损统计表”计算线损率y表示为:
Figure FDA0002586270700000011
式中,Q1为电厂上网电量;Q2为转入电量;Q3为转出省外电量;Q4为售电量。
4.根据权利要求3所述的分压线损率预测及降损潜力评估方法,其特征在于:所述S200、采用箱线图排除原始数据中的异常值,并通过多月累加的方法得到新样本;
包括:
采用箱线图排除原始数据中的异常值,箱线图定义了异常值的数值范围是[下四位数-1.5*四分位数间距,上四分位数+1.5*四分位数间距],对检测出的异常数据进行剔除;
采用多月累加的方式对原始样本进行数据处理,从而得到新样本。
5.根据权利要求1所述的分压线损率预测及降损潜力评估方法,其特征在于:S300所述线损率预测模型的训练步骤如下:
根据输入***的样本,以实际线损率与预测线损率的误差最小为目标训练BP神经网络的连接权值、偏置相关网络参数,具体如下:
建立训练集D={(x1,y1),(x2,y2),...,(xm,ym)},
Figure FDA0002586270700000021
测试集共有m个样本,每个样本的输入有d个特征,输出为l维实值向量;输出层神经元的阈值为θ={θ1,θ2,...,θl},隐藏层的阈值为δ={δ1,δ2,...,δq},输入层第i个神经元与隐藏层第h个神经元的连接权值为vih,隐藏层第h个神经元与输出层第j个神经元之间的连接权值为whj,隐藏层和输出层神经元都选用sigmoid函数作为传递函数,表示为:
Figure FDA0002586270700000022
对于第k个输入样本(xk,yk):
隐藏层第h个神经元接收到的净输入ak h和输出Ak h分别为:
Figure FDA0002586270700000023
Figure FDA0002586270700000024
输出层第j个神经元从隐藏层接收到的净输入和输出分别为:
Figure FDA0002586270700000031
Figure FDA0002586270700000032
该神经网络在输入样本(xk,yk)上的均方误差为:
Figure FDA0002586270700000033
通过误差反向传播算法采用梯度下降的策略通过不断迭代对网络进行训练,使得网络输出与期望的输出间的误差Ek最小,以隐藏层到输出层的连接权值whj为例进行具体推导:
以目标函数的负梯度方向对参数进行调整,目标函数是网络输出与期望输出的误差,即:
wh,j←wh,j-ηΔwh,j
Figure FDA0002586270700000034
式中,η表示学习率,取值范围为0到1,用于控制算法每一轮迭代中的更新步长,若太大则会发生振荡,太小则影响训练的收敛速度;
根据链式求导法则,误差Ek对whj的梯度表示为:
Figure FDA0002586270700000035
由于
Figure FDA0002586270700000036
Sigmoid函数作为传递函数满足:
f′(x)=f(x)(1-f(x)) (11)
得到
Figure FDA0002586270700000037
最后得到参数whj的更新公式为:
Figure FDA0002586270700000041
同理得到输入层到隐藏层权值以及隐藏层、输出层偏置的更新公式:
Figure FDA0002586270700000042
Figure FDA0002586270700000043
Figure FDA0002586270700000044
对于容量为m的样本集D={(x1,y1),(x2,y2),...,(xm,ym)},对应的误差为整个样本集中各样本误差的和:
Figure FDA0002586270700000045
各参数更新时的负梯度为所有样本负梯度的和,即:
Figure FDA0002586270700000046
Figure FDA0002586270700000047
Figure FDA0002586270700000048
Figure FDA0002586270700000049
选取实际线损率与BP神经网络输出预测线损率的平均绝对百分误差eMAPE作为评价预测效果的依据:
Figure FDA00025862707000000410
式中,
Figure FDA00025862707000000411
和yi分别表示测试集中第i个样本的预测线损率和实际线损率,M表示测试集中样本的总数;
重复20次实验,每次随机选取训练集和测试集,将20次实验的eMAPE的平均值作为预测模型的性能评价指标。
6.根据权利要求1所述的分压线损率预测及降损潜力评估方法,其特征在于:所述S300、把处理过得到的新样本数据,输入事先训练好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降损潜力进行评估;
包括:
假设各电压等级与外省的电量输送情况基本不变,即转入电量和转出外省电量不变,电厂上网电量的年增长率为η,则该电压等级的负荷功率Pload为:
Figure FDA0002586270700000051
式中,Q1为电厂上网电量;Q2为转入电量;Tmax为年最大负荷利用小时数;根据电网投资期望容载比和负荷功率确定新建变电站的总容量Pnew,并假设变电站的平均容量基本不变,进而确定新建变电站的座数N1
Figure FDA0002586270700000052
Pnew=N1·Pmean (24)
式中,K为期望容载比;Plast为上一年变电站容量;Pmean为变电站平均容量;[.]为向上取整函数;
假设各年每回线路输送的平均功率及该电压等级线路的平均长度保持不变,进而确定新建线路的条数N2及长度lnew
Figure FDA0002586270700000053
lnew=N2·lmean (26)
式中,Pline为单回线路平均输送功率;lmean为线路平均长度;
将预测的上网电量、线路条数、变电站容量等参数输入分压线损率预测模型中,得到电网未来的分压线损率。
7.根据权利要求2所述的分压线损率预测及降损潜力评估方法,其特征在于:
所述S400、改变输入模型数据的取值范围,进行灵敏度分析,评估各降损措施的降损潜力;包括:
基于灵敏度分析方法给出各电压等级线损率的关键影响因素,针对模型y=f(x1,x2,…,xi),xi为模型的第i个输入属性值,令每个输入属性值在可能的取值范围内变动,确定这些输入属性的变动对模型输出值的影响程度。
8.一种分压线损率预测及降损潜力评估***,其特征在于,包括以下单元:
指标体系单元,用于获取电网的运行数据,包括发售电量统计、输电线路参数和变压器参数,建立不同电压等级的线损率影响因素指标体系;
新样本处理单元,用于采用箱线图排除原始数据中的异常值,并通过多月累加的方法得到新样本;
预测评估单元,用于把处理过得到的新样本数据,输入事先训练好的线损率预测模型中,得到电网分压线损率的预测值,即对分压线损率预测及降损潜力进行评估。
9.根据权利要求8的一种分压线损率预测及降损潜力评估***,其特征在于,还包括:
灵敏度分析单元,用于改变输入模型数据的取值范围,进行灵敏度分析,评估各降损措施的降损潜力。
10.根据权利要求8的一种分压线损率预测及降损潜力评估***,其特征在于,还包括以下子单元:
线损率预测模型的训练单元,用于根据输入***的样本,以实际线损率与预测线损率的误差最小为目标训练BP神经网络的连接权值、偏置相关网络参数。
CN202010682242.9A 2020-07-15 2020-07-15 一种分压线损率预测及降损潜力评估方法及*** Active CN112085156B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010682242.9A CN112085156B (zh) 2020-07-15 2020-07-15 一种分压线损率预测及降损潜力评估方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010682242.9A CN112085156B (zh) 2020-07-15 2020-07-15 一种分压线损率预测及降损潜力评估方法及***

Publications (2)

Publication Number Publication Date
CN112085156A true CN112085156A (zh) 2020-12-15
CN112085156B CN112085156B (zh) 2024-05-03

Family

ID=73734685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010682242.9A Active CN112085156B (zh) 2020-07-15 2020-07-15 一种分压线损率预测及降损潜力评估方法及***

Country Status (1)

Country Link
CN (1) CN112085156B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114362134A (zh) * 2021-11-06 2022-04-15 国网河南省电力公司新乡供电公司 一种基于线损合格率的中压线路降损方法
CN114609480A (zh) * 2022-05-16 2022-06-10 国网四川省电力公司电力科学研究院 一种电网损耗异常数据检测方法、***、终端及介质
CN115081747A (zh) * 2022-07-27 2022-09-20 国网浙江省电力有限公司 一种基于知识图谱技术的数据处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103942727A (zh) * 2014-04-10 2014-07-23 国家电网公司 一种基于电网特征差异的线损水平评价方法
WO2015172560A1 (zh) * 2014-05-16 2015-11-19 华南理工大学 基于bp神经网络的中央空调冷负荷的预测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103942727A (zh) * 2014-04-10 2014-07-23 国家电网公司 一种基于电网特征差异的线损水平评价方法
WO2015172560A1 (zh) * 2014-05-16 2015-11-19 华南理工大学 基于bp神经网络的中央空调冷负荷的预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
田宏杰;: "线损分析预测在供电管理中的应用", 电力***保护与控制, no. 07 *
苏志雄;曾民星;孔灿;林毓;张祥华;曾江;: "基于线损影响因素的配电网线损计算方法研究", 电气应用, no. 2 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114362134A (zh) * 2021-11-06 2022-04-15 国网河南省电力公司新乡供电公司 一种基于线损合格率的中压线路降损方法
CN114362134B (zh) * 2021-11-06 2023-11-21 国网河南省电力公司新乡供电公司 一种基于线损合格率的中压线路降损方法
CN114609480A (zh) * 2022-05-16 2022-06-10 国网四川省电力公司电力科学研究院 一种电网损耗异常数据检测方法、***、终端及介质
CN115081747A (zh) * 2022-07-27 2022-09-20 国网浙江省电力有限公司 一种基于知识图谱技术的数据处理方法
CN115081747B (zh) * 2022-07-27 2022-11-11 国网浙江省电力有限公司 一种基于知识图谱技术的数据处理方法

Also Published As

Publication number Publication date
CN112085156B (zh) 2024-05-03

Similar Documents

Publication Publication Date Title
CN112488395B (zh) 一种配电网线损预测方法及***
CN112085156B (zh) 一种分压线损率预测及降损潜力评估方法及***
CN112990500B (zh) 基于改进加权灰色关联分析的台区线损分析方法及***
CN110910026B (zh) 一种跨省输电线路线损智能管理与决策方法及***
CN104573857A (zh) 一种基于智能算法优化组合的电网负荷率预测方法
CN114936742A (zh) 一种供水***调度代理决策方法
Dang et al. A novel multivariate grey model for forecasting periodic oscillation time series
CN110796293A (zh) 一种电力负荷预测方法
CN111311026A (zh) 一种顾及数据特征、模型和校正的径流非线性预测方法
CN110533247A (zh) 一种采用气温数据异常点补偿的月度用电量预测方法
CN101685964A (zh) 基于大规模并行处理的在线调度辅助决策方法和***
CN105205572A (zh) 基于EMD分解和Elman算法的光伏出力预测方法
Xu et al. Investigating the determinants of the growth of the new energy industry: using quantile regression approach
CN107679651A (zh) 一种基于灰色预测模型与bp人工神经网络模型的月用电量预测方法
Shiyun et al. Data analysis and cloud computing of power grid infrastructure projects under the background of informatization
Sun et al. Reliability prediction of distribution network based on PCA-GA-BP neural network
CN109871998B (zh) 一种基于专家样本库的配电网线损率预测方法及装置
Ban Energy decision making of steel company based on energy management system
Ma et al. Research on substation project prediction method in power transmission and transformation by improved neural network intelligent model
Li et al. The Potential of Electrifying China's Energy Use: an Analysis Based on STIRPAT Model.
Gao et al. Investment forecast of power network infrastructure project based on BP neural network
CN111369083A (zh) 一种面向电网项目投资效益的评价方法与***
Duan et al. Heavy Overload Prediction Method of Distribution Transformer Based on GBDT
Sha et al. Cost Risk Assessment of Power Transmission and Transformation Project Based on MCS-GA-ELM
CN115062871B (zh) 基于多智能体强化学习的智能电表状态评估方法

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