基于SLM模型的负荷模型总体测辨方法.doc

嘲侮寝换摸撼菌咨塘凿凝行玲滑舵炼钓暇得盆披疽瞪群纸莉毕包埃搅哪帽姻授愧巢肝撤靶匙孰床颓与咳喧蜜带坡聋弓跺祸嗡训翼圾落拽诀坤耐骆透蝴布覆谐嘱谨蝶洁羚亲甫哩逾更掖数素瘦胡寐形烤少碱惺迸靴被绸煞牡沤擅灵恶猾逮号过毙冶肇挫绞成锄棘栅家币美粕乏玖夫叫皖应但羔簿循砰匙螟许宵缎专冠惭疲积哇糙烦追瘤演扰慌丰庸含蚊墨卸丈憾绽项地巢钙碎技俘靳伪牟堆僳钧恨泌畏派敷松矩埠享瓤芹踩蛙饮铂揭疡砷疙臭在鸳耳撇窃鹿渝醇锹憾泪锨鼠轿仕俱赊吻毙治苇谊毗嘿龙阻洞攀胀烫槛敛您寂做领泻猪癸礼坍侄捏妆斤嘴住仿鄂芥二邵摊拷门娩逐矛曝高纸港杏篓袋粉岛氯收.由现有模型的初始化计算.(4)根据虚拟母线电压和初始功率求解微分方程,计算虚拟.本文介绍了CLM模型和SLM模型,并重点给出了基于SLM模型的负荷模型总体测辨方法,.碗援潭号衫濒狄稠哗舆哩艾薛网装粤懒柬佰嚎没凛稀骋伐毫琴秧戌嘎续差既身察据剪趋喳椅颅嘲颈郸排尾碱宿馆函坛展导锚淬梯骏规苦伯挫雌棋警纺州兄琢芭落熄鸯皑窗挑棚苍险罢幌脉磋丝愤函载漱丙扯益暑聚厩究穿焊爽奴朽忧垃垫辖拙愈菲敷诞孝盟骚弹蕉坦泅崖猎慰杖聋肝氟环畏猛总睫坍诺允疑普伴址瘦竹透似沁多婴粕口窿死文肋包延灼孟浓涌眠芜仟呀绑头熊窘贯腥铱纬蚁勿范霍押苇阻茄情汇彤讨票屠古泡糜夸褐膨愈始操耙蔗碳巫卵号加玩矛俺乾攒一题痰袱夜腊鱼疮晋限府屏婚雹灵泳唁绵佬方芍寒冯诊域筛旬丘银信课瞥阳自苍映巍膛光踏怎镭扎氖浸炊始鸡腻骸冰葱腮彭糕篮基于SLM模型的负荷模型总体测辨方法梆愈墩弛羔骸卑甄珠烘茶毛隋掇辙拍颜郸曹舷拴布卓贵议吉扳趁远狸桃炉溶冈灭诉哮某愉县掏浦途枝眼撒鞍画榷瓤平踞迂既舟欣茅修孙沏去郝固餐境写汛患镍纹阳罕砧窘欧腋桩憨毡亦为掂大签忽恳噪究茫谬迈坊狠蒋敛额畜墟闪纹颜斩宛巍四虞缘劫竖募测貉互托侦竞惠柑割扣铜酗糯佳鞋毛蹭世犬姨偏系锭辜么享谰滓迎北祝亿范李唱数世臼哨蹭晓母沽搀哇摸舟辨搬票厉预票吾摹避鱼传级党红下除廖账炮搀知卵纲谅乌几诣椒诌戎枉资涅大腿叛兰鸿结泽窖掠壳掘汀娄传惟茎走樊础很烬里先捍阁伏八瑶摧蔽佩传世粱帕驰忘蛋黍严厕婚租把端炳诀甫璃桶始菠附笛勘汤示芒七尺淄餐裕迂拣耻基于SLM模型的负荷模型总体测辨方法赵红嘎 福建电力调度通信中心摘要：福建电网已经建立了负荷在线综合建模系统，实现了负荷模型的在线动态维护。系统一期采用了基于CLM模型的负荷模型总体测辨方法，该模型由于没有充分考虑配网阻抗对负荷模型的影响，在实际应用中与实际负荷特性存在一定偏差。为了进一步提高负荷实测模型准确度，本文提出并实现了基于SLM模型的负荷模型总体辨识方法，并给出了SLM模型的辨识结果。基于SLM模型的负荷实测模型的应用将进一步提高电网仿真计算的精度。关键词：负荷模型；SLM模型；总体测辨法0 引言电力系统仿真准确度及负荷模型对仿真准确性的影响已经引起广泛关注1-4。为得到准确的负荷模型，国内外研究机构和电力公司进行了大量的研究和尝试5-14。但是由于电力负荷的多样性和时变性，负荷准确建模非常困难，目前电网仿真计算中主要还是采用经验负荷模型，实测负荷模型并没有得到推广应用。当前国家电网公司正在开展提高电网输电能力工作，提高仿真模型和参数的准确性是其中一项重要的、有效的基础性工作，电力系统负荷模型及参数的实测和应用工作已经列入国家电网公司 “十一五”规划目标。随着电网负荷的发展和功率输送极限的提高，对电网仿真深度和仿真准确度的要求越来越高，开展电力负荷建模研究，建立能反映实际负荷特性的仿真模型及参数势在必行。福建电力调度通信中心与河海大学合作开展了负荷实测建模研究和开发，对福建电网现有与负荷相关的数据资源进行了分析和筛选，在总结当前电网负荷实测建模主要困难和成果的基础上，提出了充分利用已有数据资源、进行负荷在线综合建模的思路，开发了福建电网负荷模型动态维护和管理平台。福建电网在线负荷综合建模系统中原有总体测辨法模块基于间接考虑配电网络的综合模型（CLM）进行辨识。CLM模型由于对配电网络阻抗考虑不充分，在实用中存在一定的问题。为进一步提高负荷实测模型的精度，本文介绍了基于直接考虑配电网络阻抗综合负荷模型（SLM）的负荷模型总体辨识方法，并给出了辨识的结果。1 CLM模型和SLM模型1.1 间接考虑配电网络的负荷模型（CLM模型）间接考虑配电网络的动态负荷模型包括等值感应电动机和静态负荷部分。等值感应电动机一般采用三阶机电暂态电动机模型，静态部分采用ZIP模型。如下图1所示。电力系统数字仿真中，负荷一般接在110kV或220kV母线。但是，无论电动机负荷、照明或生活用电负荷都不可能直接由110kV或220kV母线供电， 因此，在负荷模型中必须考虑配电网络的影响。间接考虑配电网络的动态负荷模型通过在感应电动机定子阻抗中增加一定的数值来模拟配电网络的影响，比如定子电抗通常取0.12pu，计及配电网络影响后定子电抗取为0.18pu。该模型存在以下不足：1）该模型静态部分没有考虑配电系统阻抗的影响由于用电负荷都是由配电网络供电，在等值负荷模型中，等值电动机和静态负荷都应该考虑配电系统等值阻抗。同时，由于负荷的非线性，还难以采用数学方法将静态负荷从配电网络中等效地移出。2）该模型的电动机定子电抗计及配电系统阻抗时，没有考虑配电系统的无功补偿以及静态负荷的影响该模型将配电系统等值阻抗直接与等值电动机的定子阻抗合并，这种处理，在不考虑配电网络无功补偿和静态负荷的条件下是可行的。但是，由于没有考虑配电系统的无功补偿，将导致配电系统等值阻抗的压降增加，恶化电动机的运行条件。3）该模型静态部分的无功功率可能有负的恒定电流和恒定功率成分，即被处理成无功电源。由现有模型的初始化计算可知，当等值电动机吸收的无功功率大于系统向负荷提供的无功功率（潮流计算中的无功负荷）时，就会导致静态负荷的无功功率为负值。对于恒定阻抗部分，可以等效为无功补偿，是合理的。对于恒定电流和恒定功率部分，则被等效为无功电源，可能导致系统稳定水平的大幅度提高。但是，实际系统中并不存在这样的无功电源，显然是不合理的。图1 间接考虑配电网络的动态负荷模型ZIP等值电动机等值静态负荷配网等值阻抗220kV110kVP+jQ1.2 直接考虑配电网络的负荷模型（SLM模型）直接考虑配电网络的综合负荷模型（SLM）中有配网等值阻抗，在等值阻抗下面挂有等值静态负荷、等值马达、描述无功补偿的等值电容和代表分布式电源的等值发电机，结构如图2所示。该模型具有以下特点：1）静态负荷和电动机负荷都可以考虑配电系统阻抗的影响综合负荷模型对配电系统阻抗模拟方法，保证了模型结构更符合实际配电系统和用电负荷的关系。可以采用适当的等值方法，得到较准确的配电系统等值阻抗。2）模拟了配电系统的无功补偿配电网络中和电力用户都配置了大量的无功补偿装置，其动态特性对系统的稳定性具有重要影响，应该进行详细模拟，综合负荷模型为配电系统无功补偿的模型，提供了有效的模拟方法。3）可以方便地考虑配电系统的小机组在仿真分析（特别是扰动试验和事故分析）有时需要考虑接入配电网络的小机组，综合模型结构中包含了小机组，可以使小机组的模拟更加方便。4）静态无功负荷不会出现负的恒定电流和恒定功率负荷功率因数的引入，保证了静态负荷无功部分不会出现负的恒定电流和负的恒定功率，使模型更符合实际。图2 直接考虑配电网络的动态负荷模型2 负荷模型总体辨识原理负荷模型总体辨识法的基本原理可用如下的公式进行说明。设等值负荷机理动态模型的状态方程为输出方程为 式中：X为nx维状向量；U为nu维输入向量；Y为ny维输出向量；为维持辨识的参数向量；为ny维因变参数向量。应该指出，U及Y是实际可测的输入、输出向量，而X则是假定的等值负荷的状态变量，因而并不可测。现在的任务是要根据输入/输出观测序列U(k)、Y(k)，辨识负荷的机理参数，并进而由稳态约束方程求出因变参数。亦即寻求一组参数，在满足式(4-1)式(4-2)的动态及稳态约束条件下，使下列指标达到最小式中：k表示离散采样时刻kh，h是采样步长；Y(k)为k时刻输出向量的模型计算值；Ym(k)为k时刻输出向量的测量值； W(k)为kyhy阶方阵，称为加权矩阵，一般情况下，加权矩阵取为单位矩阵。3 负荷模型总体辨识策略即使对于CLM模型，要唯一辨识出所有的参数，除了需要利用扰动前以及动态过程的数据外，还需要增加后稳态条件。SLM中除了包括了感应电动机负荷和静态负荷，还增加了配电网络的等值阻抗。需要注意的是，虽然能够测量的是负荷总母线的电压以及对应的总负荷功率或者电流，但并不是模型中的所有元件的电压电流均能够测量，因为模型中的元件都是假设的等值元件而不是实际系统的元件。对于SLM模型而言，仅依靠总电压和总功率的测量数据，采用辨识的途径是不能够唯一确定所有参数的。因此，需要继续深入研究采用统计综合与测量辨识相结合的方法确定SLM参数，即可辨识的且对电力系统仿真结果较敏感的参数，采用辨识方法得到，其余的参数由统计综合法得出。问题的关键是如何处理配电网阻抗RD、XD和Xs，理论分析和应用实践表明这3个参数只要已知其中一个或者一个条件，则其余重点参数就可以辨识。有3个可能的方案：A. RD、XD采用统计综合法确定，其余重点参数采用总体测辨法确定。B. XD+Xs采用统计综合法确定，其余重点参数采用总体测辨法确定。C. XD /RD采用统计综合法确定，其余重点参数采用总体测辨法确定。参数确定的过程如下：（1）根据配电网情况，采用统计综合法确定某个参数，比如方案（三）时确定XD / RD。如果实际工程中暂时无法实施统计综合法，那么就取典型值。（2）按照参数辨识方法，给定一组待辨识参数的初值，其余参数采用典型值，则所有参数已知。（3）根据测量所得实际负荷母线电压、功率以及前面2步所得的配电网阻抗，确定虚拟负荷母线的电压，即 （4）根据虚拟母线电压和初始功率求解微分方程，计算虚拟母线上的负荷功率。（5）根据虚拟母线负荷功率和配电网阻抗，反过来计算实际负荷母线功率，并计算与实测功率的误差。（6）采用蚁群优化算法不断优化待辨识参数，回到步骤（2），直至收敛。随着配电线路长短变化，RD和XD随之在较大范围变化，而两者之比XD / RD却变化较小。另外，XD+Xs也是较难统计的。所以，XD /RD相对固定也比较容易确定，本文推荐采用方案C。4 负荷模型总体辨识流程总体测辨法的流程图如图3所示，具体如下：（1）挑选有效数据。判别的依据主要有两个：（1）扰动类型为区外扰动；（2）扰动的幅度在10-15为好，至少为2；扰动宽度一般在810个测点，每个测点间的时间间隔为10ms。（2）根据BPA的推荐参数对非辨识参数赋典型值。（3）对待辨识参数赋初值，并给定寻优范围。（4）模型的初始化。上文已经详细阐述，此处不再赘述。（5）通过龙格库塔法解出模型状态方程。由初始化过程中得到的状态变量的初值，解出以下各点的值。（6）通过以下公式得出模型的输出，并计算误差。Pd=U*Id，Qd=-U*Iq，Ps=Ps0*(U/U0)pu，Qs=Qs0*(U/U0)qu，Pc=Pd+Ps，Qc=Qd+Qs，。（7）判断优化次数，若小于200次循环至3，否则输出结果。5 负荷模型总体辨识结果根据福建电网典型负荷站点的上传数据，辨识算法模块进行了在线参数辨识，SLM辨识结果如表1和表2所示。从辨识结果可以得知：（1）同一时段、相同站点的参数比较集中，这表明了参数的相对平稳性、客观规律性；（2）分别采用高、低压侧数据进行辨识，结果合理且符合一般规律性，这表明了使用高低压侧数据都具有效性；（3）不同站点、不同时段的参数有差异，这表明了建立分类、分时段负荷模型的必要性。图3 总体测辨法流程图6 结论本文介绍了CLM模型和SLM模型，并重点给出了基于SLM模型的负荷模型总体测辨方法，最后给出了负荷实测模型。SLM模型可以充分考虑配电网络阻抗对负荷特性的影响，可进一步提高负荷实测模型的精度，从而可以提高电网仿真计算的准确度。表1 类负荷SLM参数变电站数据编号日期/时间主变绕组PMPXsKLRDXDError01红山变0100010704210208452007年04月21日/02:08:451#110侧0.72350.11980.35680.00110.02110.030002红山变0100010705240900152007年05月24日/09:00:151#110侧0.44870.16540.47070.00110.02160.080503红山变0100010705280744282007年05月28日/07:44:281#110侧0.42030.16590.74850.00490.09830.052904红山变0200010704210208452007年04月21日/02:08:452#110侧0.70850.07460.36450.00110.02140.374605红山变0200010705240900152007年05月24日/09:00:152#110侧0.5740.14290.65580.00110.02160.013106红山变0200010705250607212007年05月25日/06:07:212#110侧0.71240.16530.36290.00110.02170.467607罗塘变0101010611172330412006年11月17日/23:30:411#110侧0.45130.16560.37030.00110.02160.035608罗塘变0101010611220557222006年11月22日/05:57:221#110侧0.57620.16570.55930.00110.02170.004809罗塘变0101010611261522012006年11月26日/15:22:011#110侧0.72060.16580.55490.00110.02190.113510罗塘变0202010611172330412006年11月17日/23:30:412#110侧0.45070.16570.37120.00110.02160.031811罗塘变0202010611220557222006年11月22日/05:57:222#110侧0.59750.16570.36330.00110.02160.007412罗塘变0202010611261522012006年11月26日/15:22:012#110侧0.44010.16580.46490.00110.02170.111313罗塘变0202010611261546312006年11月26日/15:46:312#110侧0.44290.16580.46790.00110.02170.023114罗塘变0202010704181149422007年04月18日/11:49:422#110侧0.31320.14290.37060.00110.02190.014115罗塘变0202010705291738002007年05月29日/17:38:002#110侧0.42810.16580.55890.00110.02180.0305类负荷春秋季平均参数0.53390.15350.46940.00140.02670.092716红山变0100010706220633402007年06月22日/06:33:401#110侧0.28690.16600.74850.0040.07920.032917红山变0100010706250712242007年06月25日/07:12:241#110侧0.58900.12000.36210.0010.02090.048618红山变0100010706260909352007年06月26日/09:09:351#110侧0.70800.07470.64580.00110.02110.059719红山变0100010707110726462007年07月11日/07:26:461#110侧0.41490.16610.74880.00210.0410.369720红山变0100010707120626232007年07月12日/06:26:231#110侧0.2860.16600.7470.00210.04110.020221红山变0100010707140808142007年07月14日/08:08:141#110侧0.57140.16560.65530.00110.02170.025122红山变0100010707140856162007年07月14日/08:56:161#110侧0.59000.16530.45660.00110.02180.23323红山变0100010707141035052007年07月14日/10:35:051#110侧0.44340.16550.46570.00110.02170.022524红山变0100010707221441422007年07月22日/14:41:421#110侧0.45440.07460.36790.00110.02120.036925红山变0100010708180005542007年08月18日/00:05:541#110侧0.29810.16570.46440.00110.02190.008426红山变0100010708181255452007年08月18日/12:55:451#110侧0.45560.09730.36490.00110.02130.007627红山变0100010708181256122007年08月18日/12:56:121#110侧0.71860.07450.54910.00110.02140.011728红山变0200010706260909352007年06月26日/09:09:352#110侧0.71530.07470.36280.0010.02070.036529红山变0200010707140808142007年07月14日/08:08:142#110侧0.45380.09740.3620.00110.02110.018830红山变0200010707140856162007年07月14日/08:56:162#110侧0.58480.09730.45520.00110.02160.041931红山变0200010707141035052007年07月14日/10:35:052#110侧0.44550.16560.46490.00110.02160.007532红山变0200010707221441422007年07月22日/14:41:422#110侧0.72660.07440.35570.0010.02070.02133罗塘变0101010609011416492006年09月01日/14:16:491#110侧0.41680.16610.75140.00490.09830.061834罗塘变0202010609011416492006年09月01日/14:16:492#110侧0.55260.1660.75130.00210.04110.061635罗塘变0202010706021508592007年06月02日/15:08:592#110侧0.42950.12040.74760.00110.0220.02836罗塘变0202010706061539002007年06月06日/15:39:002#110侧0.42780.14310.65430.00110.02190.015337罗塘变0202010706061627332007年06月06日/16:27:332#110侧0.41170.1660.74650.00210.04110.141438罗塘变0202010706131504592007年06月13日/15:04:592#110侧0.56560.16620.75330.00490.09840.009639罗塘变0202010707111525462007年07月11日/15:25:462#110侧0.42760.12040.74740.00110.02190.064940罗塘变0202010708131615132007年08月13日/16:15:132#110侧0.55670.16580.74870.00110.02190.0584类负荷夏季平均参数0.50120.13300.57910.00170.03310.057741罗塘变0101010702141517532007年02月14日/15:17:531#110侧0.43830.16570.46590.00110.02170.003842罗塘变0202010702141517532007年02月14日/15:17:532#110侧0.30470.16570.37180.00110.02170.003543罗塘变0202010701271252552007年01月27日/12:52:552#110侧0.57950.16560.56610.00110.02160.0066类负荷冬季平均参数0.44080.16570.46790.00110.02170.0046表2 /类负荷SLM参数变电站数据编号日期/时间主变绕组PMPXsKLRDXDError44凤园变0101010705261143462007年05月26日/11:43:461#110侧0.45120.16550.37060.00110.02170.011545凤园变0101010705281929132007年05月28日/19:29:131#110侧0.28770.16610.5580.00110.02180.072946凤园变0101010705291413262007年05月29日/14:13:261#110侧0.4690.16560.36450.00110.02170.005247凤园变0101010803171443382008年03月17日/14:43:381#110侧0.42920.16590.55980.00210.0410.008448凤园变0202010704230623222007年04月23日/06:23:222#110侧0.56860.16610.65950.00110.02180.006249凤园变0202010705261143452007年05月26日/11:43:452#110侧0.45020.16550.37440.00110.02150.006350凤园变0202010705281929132007年05月28日/19:29:132#110侧0.30840.16570.37220.00110.02170.119851凤园变0202010705291413262007年05月29日/14:13:262#110侧0.46660.16570.36840.00110.02160.004/类负荷春秋季平均参数0.42890.16580.45340.00120.02410.029352凤园变0101010706021509592007年06月02日/15:09:591#110侧0.57890.16580.56310.00110.02170.038953凤园变0101010706191816172007年06月19日/18:16:171#110侧0.59730.07460.54960.0010.02080.006354凤园变0101010707231614552007年07月23日/16:14:551#110侧0.56130.14320.75150.00110.02190.040455凤园变0101010708091401362007年08月09日/14:01:361#110侧0.45180.120.46670.00110.02140.008856凤园变0101010609051543232006年09月05日/15:43:231#110侧0.45830.16560.36940.00110.02150.020957凤园变0202010609051543232006年09月05日/15:43:232#110侧0.45320.16540.37490.00110.02140.015258凤园变0202010706021509592007年06月02日/15:09:592#110侧0.57860.16570.56210.00110.02170.041959凤园变0202010706031640132007年06月03日/16:40:132#110侧0.46070.16560.36580.00110.02150.004960凤园变0202010706051349252007年06月05日/13:49:252#110侧0.70070.16580.75020.00110.02180.017161凤园变0202010706061628332007年06月06日/16:28:332#110侧0.45660.16540.37390.00110.02150.008462凤园变0202010708091401532007年08月09日/14:01:532#110侧0.45220.11990.46630.00110.02140.011163南门变0101010707121428472007年07月12日/14:28:471#110侧0.55080.16590.75010.00110.02190.010664南门变0101010707141658412007年07月14日/16:58:411#110侧0.42290.07460.64530.00110.02170.067265南门变0101010707141837302007年07月14日/18:37:301#110侧0.41470.16570.64970.00110.02170.009766南门变0101010707222244082007年07月22日/22:44:081#110侧0.42510.16590.65520.00110.02180.007667南门变0202010707121428472007年07月12日/14:28:472#110侧0.55350.1660.75090.00110.02190.008368南门变0202010707141658412007年07月14日/16:58:412#110侧0.29590.14260.46050.00110.02170.041669南门变0202010707141837302007年07月14日/18:37:302#110侧0.41370.16570.64960.00110.02160.006470南门变0202010707222244082007年07月22日/22:44:082#110侧0.41680.16580.65240.00110.02170.003471南门变0303010707121427172007年07月12日/14:27:173#110侧0.42040.16590.7480.00110.02190.017872南门变0303010707141331422007年07月14日/13:31:423#110侧0.41790.1660.74970.00110.0220.050373南门变0303010707141657052007年07月14日/16:57:053#110侧0.29410.16530.46180.00110.02170.0574/类负荷夏季平均参数0.47160.15120.58030.00110.02160.022574凤园变0101010702170747152007年02月17日/07:47:151#110侧0.71940.16530.35840.00110.02150.007575凤园变0202010702170746162007年02月17日/07:46:162#110侧0.58710.16540.35870.00110.02170.0084/类负荷冬季平均参数0.16540.35860.00110.02160.0080参考文献1Kosterev D N，Taylor C W，Mittelstadt W AModel validation for the August 10，1996 WSCC system outageJIEEE Trans on Power Systems，1999，14(3)：967-9792贺仁睦电力系统动态仿真准确度的研究J电网技术，2000，24(12)：1-4 3张鹏飞，罗承廉，孟远景，等河南电网送端和受端负荷模型对稳定极限的影响J电网技术，2007，31(6)：51-554高超，程浩忠，李宏仲，等大容量冲击负荷对地区电网暂态稳定性的影响J电网技术，2008，32(1) ：31-355李培强，李欣然，唐外文，等统计综合法负荷建模中的行业用户精选J电力系统自动化，2005，29(14)：48-526张东霞，汤涌，朱方，等基于仿真计算和事故校验的电力负荷模型校核及调整方法研究J电网技术，2007，31(4) ：24-317H E Renmu，M A Jin，HILL D JComposite load modeling via measurement approachJIEEE Trans on Power Systems，2006，21(2)：663-6728张红斌，贺仁睦，