(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111675848.0 (22)申请日 2021.12.31 (71)申请人 甘肃电通电力工程设计咨询有限公 司 地址 730050 甘肃省兰州市七里河区西津 东路628号 (72)发明人 宋军　詹伟　张海平　张海生　 王利平　王鑫　陈英普　蒋明华　 熊峰　徐建委　王仕俊　王文婷　 汪静　刘科　管虎　杨枢捷　 张江虹　张洪瑞　张亚明　俞骁　 郭化诚　邵生龙　孟佳佳　张小花　 王昕怡　王明浩　康倩倩　吕胤禛　 (74)专利代理 机构 北京众合诚成知识产权代理 有限公司 1 1246 代理人 陈波(51)Int.Cl. G06F 30/20(2020.01) G06Q 50/06(2012.01) G06F 17/18(2006.01) G06F 119/14(2020.01) (54)发明名称 一种系统侧谐波占主导情况下的谐波责任 定量计算方法 (57)摘要 本发明公开了一种系统侧谐波占主导情况 下的谐波责任定量计算方法， 具体是基于动态加 权阻抗模型进行谐波责任定量评估， 具体包括： 首先， 构建动态加权阻抗模型静态部分和动态部 分； 接着， 对动态加权阻抗模型参数进行优化； 然 后， 通过静态部分和动态部分阻抗值计算第h次 谐波频率下用户侧等效阻抗值； 最后， 求出谐波 电压责任。 本发 明基于动态加权阻抗模型计算等 效阻抗量化谐波源对关注母线的谐波责任， 解决 了系统侧谐波占主导情况下谐波责任计算的问 题， 扩展了谐波责任计算范围， 可 以高效地划分 供用电双方的谐波责任。 权利要求书4页 说明书8页 附图3页 CN 114330003 A 2022.04.12 CN 114330003 A 1.一种系统侧谐波占主导情况下的谐波责任定量计算方法， 基于动态加权阻抗模型， 其特征在于， 具体包括以下步骤： 步骤1： 构建动态加权阻抗模型静态阻抗部分和动态阻抗部分； 步骤2： 对动态加权阻抗模型动静态权 重进行优化； 步骤3： 通过静态部分和动态部分阻抗 值计算用户侧第h次谐波等效阻抗 值； 步骤4： 通过步骤3中求得的用户侧h次等效谐波阻抗 求出谐波电压责任。 2.根据权利要求1所述的一种系统侧谐波占主导情况下的谐波责任定量计算方法， 其 特征在于， 步骤1中， 对用户侧负荷构建动态加权阻抗模型， 该动态加权阻抗模型由静态阻 抗部分与动态阻抗部分加权串联而成， 其中所述静态阻抗部分是恒定阻抗模型， 所述动态 阻抗部分是感应电动机模型； 假定Xs是定子电抗， Rr是转子电阻， Xr是转子电抗， Xm是励磁电 抗， s是转子滑差； j表示虚数单位， C表示静态模型部 分， M表示动态模 型部分， 即感应电动机 模型； 则所述动态加权阻抗模型中静态阻抗的计算采用恒定阻抗模型， 表示为如式(1)所 示： 其中， U表示关注负荷基波电压测量值， P是关注负荷基波有功功率测量值， Q是关注负 荷基波无功 功率测量 值， 表示关注负荷基波相角值； 所述动态加权阻抗模型中动态阻抗被表示 为如式(2)所示： ZM＝jXs+(jXm)//(jXr+Rr/s)              (2)。 3.根据权利要求1所述的一种系统侧谐波占主导情况下的谐波责任定量计算方法， 其 特征在于， 步骤2中， 所述将动态加权阻抗模型中动静态权 重进行优化具体包括： 在所述动态加权阻抗模型中， 将 待识别的参数表示为{Xs,Xm,Xr,Td0}， 其中， Td0为电磁参 数， 通过最小二乘 ‑遗传算法优化动静态权 重； 将相邻N个样本为 一组进行优化， 具体优化模型如下： 对于第i个样本， 适应度函数f(i)表示该样本的功率误差平方和的相 反数， 如式(3)所 示： f(i)＝‑((Pp(i)‑P(i))2+(Qp(i)‑Q(i))2)          (3)， 其中P， Q分别为有功功率和无功功率的实测 值； Pp， Qp分别为有功功率和无功功率的预 测值， Pp的计算公式如式(4)所示： 式中， ZC、 X1、 Ud、 Uq、 Edp、 Eqp为待计算参量； 其中ZC是静态阻抗值， 通过恒定阻抗模型式 (1)计算； X1是电磁参数， 由待识别参数 Xm、 Xr、 Xs表示为如式(5)所示：权　利　要　求　书 1/4 页 2 CN 114330003 A 2式中， Ud,Uq分别表示样本的A相基波相量电压在p， q轴上的分量， 通过电压幅值U和相角 来计算， 表示 为如式(6)所示： 式中Edp,Eqp是电动势相量在d轴和q轴上的分量的预测值， 通过式(7)、 (8)进行计算： 公式(7)和(8)中的变量Ed、 Eq是电动势相量在d轴和q轴上的分量， 基于电磁动力学和显 示积分法求 解， 分别由式(9)、 (10)表示： 公式(7)和(8)中的X表示 为式(11)所示： X＝Xs+Xm                     (11)， 公式(7)和(8)中的Id,Iq分别表示样本的A相基波相量电流值在p， q轴上的分量， 通过电 流幅值I和相角 来计算， 表示 为式(12)所示: 公式(7)和(8)中的滑差s表示 为式(13)所示： 式中的T为辨识周期， 即相邻N个样本点的采样 时间； w0为同步转角速度， 由转子电阻Rr 和待识别参数 Xr,Xm,Td0计算， 表示 为式(14)所示：权　利　要　求　书 2/4 页 3 CN 114330003 A 3