CN105958483B - 考虑负荷特性的居民配电台区三相不平衡潮流计算方法 - Google Patents

考虑负荷特性的居民配电台区三相不平衡潮流计算方法 Download PDF

Info

Publication number
CN105958483B
CN105958483B CN201610383169.9A CN201610383169A CN105958483B CN 105958483 B CN105958483 B CN 105958483B CN 201610383169 A CN201610383169 A CN 201610383169A CN 105958483 B CN105958483 B CN 105958483B
Authority
CN
China
Prior art keywords
mrow
msub
mover
centerdot
load
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.)
Expired - Fee Related
Application number
CN201610383169.9A
Other languages
English (en)
Other versions
CN105958483A (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201610383169.9A priority Critical patent/CN105958483B/zh
Publication of CN105958483A publication Critical patent/CN105958483A/zh
Application granted granted Critical
Publication of CN105958483B publication Critical patent/CN105958483B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Peptides Or Proteins (AREA)

Abstract

本发明公开一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，在负荷模型方面，采用基于双笼型转子感应电动机模型的空调负荷计算模型与照明为主的静态负荷计算模型；在三相潮流计算模型方面，对负荷节点采用abc三相功率平衡方程和0相电流平衡方程，对配电变压器低压侧联络节点的abc三相和0相都采用电流平衡方程。本发明的方法能够降低三相潮流计算的编程求导计算量，易于实现；并准确地反映居民低压配电台区在各种运行工况下的三相不平衡程度并评估及其对配电网安全优质经济运行的影响。

Description

考虑负荷特性的居民配电台区三相不平衡潮流计算方法
技术领域
本发明涉及电力运行分析领域，具体涉及一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法。
背景技术
低压配电网是电力系统中直接连接用户的环节，其安全优质运行直接影响到用户的正常用电。由10kV降压到380V的配电变压器和其二次侧连接的380V出线及低压负荷所组成的低压配电台区，是低压配电网的基本形式。目前，我国的大多数低压配电台区，由于三相供电与单相供电并存，三相负荷接入不均，三相不平衡运行情况严重，馈线上负荷点的三相电压不平衡程度往往较大，严重影响到用户的供电电压质量。另外，三相不平衡程度的增大还会加大零线的电流和有功损耗，从而影响到低压配电台区运行的经济性。特别是对于众多的居民低压配电台区，由于台区中含有大量的单相空调负荷，而空调负荷在启动时的电流和功率比正常运行时要大几倍，启动瞬时往往会进一步加大低压配电台区的三相不平衡程度。
常规的配电网三相不平衡潮流计算中，负荷采用不随时间和运行状态变化的恒功率模型，没有反映实际负荷的变化特性，采用负荷节点恒功率模型的三相潮流计算评估配电台区运行的三相不平衡程度，与实际情况往往差别很大。特别是对于空调负荷启动的运行工况，由于其电流和功率达到正常运行时的几倍，采用正常运行功率值的恒功率模型进行三相潮流计算结果与实际电网状态差别很大，无法准确反映低压配电台区运行的三相不平衡程度。因此，对于含有大量单相空调负荷的众多居民低压配电台区，如何准确地反映台区在各种运行工况下的三相不平衡程度，亟需提出合理的三相潮流计算模型和方法。
发明内容
本发明的目的在于针对上述现有技术中存在的问题，提出一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，考虑居民配电台区中单相空调负荷的实际运行特性，能够模拟实际居民低压配电台区中空调负荷稳定运行和启动条件下台区的三相不平衡允许程度，并能够准确评估三相不平衡运行对于低压配电网供电安全性、电压质量和经济性的危害。
为达到上述发明的目的，本发明通过以下技术方案实现：
一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，包括如下步骤：
步骤S10，建立负荷节点的包括空调负荷模型和静态负荷模型的总负荷模型；所述空调负荷模型是空调负荷从配电网吸收的功率关于转差率的函数，空调负荷模型包括有空调稳定运行状态和空调启动状态；所述静态负荷模型是接入点静态负荷功率关于接入点电压幅值的幂函数模型；
步骤S20，先给定各个节点的电压初值，当前为第k次迭代，遍历N个负荷节点，所述k和N均为等于或者大于1的整数，依次判断各负荷节点的空调状态，并计算相应空调状态的接入点空调负荷功率；
步骤S30，依次计算各负荷节点的总负荷功率，并通过以配电变压器高压侧为平衡点的各负荷节点的三相潮流计算模型、以及通过配电变压器低压侧联络节点的三相潮流计算模型，计算得到负荷节点的节点电压和节点相角；
步骤S40，判断空调负荷模型和三相潮流计算模型是否均收敛，若是则当前节点电压和节点相角为最终解，否则将当前节点电压设定为迭代初值，k＝k+1后进入步骤S20。
进一步，所述空调负荷模型采用双笼型转子感应电动机等值电动模型。
进一步，所述空调稳定运行状态的空调负荷模型是接入点空调负荷功率关于空调稳定运行状态的转差率的函数。
进一步，所述空调启动状态的空调负荷模型是接入点空调负荷功率关于空调转差率为1的参数表达式。
进一步，所述步骤S20采用牛顿-拉夫逊算法计算相应空调稳定运行状态的转差率进而计算接入点空调负荷功率。
进一步，所述静态负荷模型表示如下：
其中，Ps、Qs分别为静态负荷接入网络的实际有功功率与无功功率，V为接入点的实际电压幅值，VN为额定电压，PN和QN为额定条件下负荷的有功功率和无功功率。
进一步，所述方法还包括有步骤S50，计算反映台区运行的三相不平衡程度的指标，所述指标包括台区网损、配电变压器负载率和各个节点的电压不平衡度。
进一步，步骤S40所述配电变压器高压侧为平衡点的各个负荷节点三相潮流计算模型如下：
其中，Pi,p和Qi,p分别表示第i个节点第p相从电网吸收的有功功率和无功功率，Vi,p表示第i个节点的第p相电压，Vj,p表示第j个节点的第p相电压，Zij表示第i个节点与第j个节点之间的abc相线路的阻抗，Zij,0表示第i个节点与第j个节点之间的0相线路的阻抗。
进一步，步骤S40所述配电变压器低压侧联络节点的三相潮流计算模型如下：
其中，表示配电变压器高压侧对低压侧的变比，为的共轭函数，θ为配电变压器高压侧对低压侧的电压移相角，Vn,p表示配电变压器高压侧的第p相电压，Vi,p表示配电变压器低压侧第p相电压，Vj,p表示非配电变压器侧第j个节点第p相电压，Zni表示配电变压器的等值阻抗，Zij表示配电变压器低压侧与第j个节点之间的线路阻抗，Zm表示配电变压器低压侧中性点对地阻抗。
进一步，步骤S30采用牛顿-拉夫逊算法求解以配电变压器高压侧为平衡点的各个负荷节点和配电变压器低压侧联络节点各自的三相潮流计算模型，求得配电网各节点电压和节点相角。
本发明的一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，具有如下有益效果：
在获取居民低压配电台区的结构和运行数据的基础上，通过求解本发明提出的考虑负荷特性的居民低压配电台区三相不平衡潮流模型，模型中考虑了居民低压配电台区以空调为主动态负荷和以照明为主静态负荷的运行特性，可以准确获得空调负荷正常运行及启动条件下台区的电压和功率分布情况，从而准确地反映居民低压配电台区在各种运行工况下的三相不平衡程度并评估及其对配电网安全优质经济运行的影响，为有效治理台区三相不平衡运行提供决策依据。
附图说明
图1为本发明的考虑负荷特性的居民配电台区三相不平衡潮流计算方法的步骤流程图。
图2为空调负荷的稳态等值电路示意图。
图3为本发明的考虑负荷特性的居民配电台区三相不平衡潮流计算方法的交替迭代算法流程图。
图4为低压配电网台区模型图。
图5为电压配电台区三相接线图。
具体实施方式
下面结合附图和实施例对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部实施例。
本发明提出考虑实际负荷特性的居民低压配电台区三相不平衡潮流计算模型，并评估居民低压配电台区三相不平衡运行的危害。在负荷模型方面，考虑基于双笼型转子感应电动机模型的空调负荷计算模型与照明为主的静态负荷计算模型；在三相潮流计算模型方面，对负荷节点采用abc三相功率平衡方程和0相电流平衡方程，对配电变压器低压侧联络节点的abc三相和0相都采用电流平衡方程。
一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，包括如下步骤：
步骤S10，建立负荷节点的包括空调负荷模型和静态负荷模型的总负荷模型；所述空调负荷模型是空调负荷从配电网吸收的功率关于转差率的函数，空调负荷模型包括有空调稳定运行状态和空调启动状态；所述静态负荷模型是接入点静态负荷功率关于接入点电压幅值的幂函数模型；
步骤S20，先给定各个节点的电压初值，采用交替迭代算法，当前为第k次迭代，遍历N个负荷节点，所述k和N均为等于或大于1的整数，依次判断各负荷节点的空调状态，并计算相应空调状态的接入点空调负荷功率；
步骤S30，依次计算各负荷节点的总负荷功率，并通过以配电变压器高压侧为平衡点的各负荷节点三相潮流计算模型和配电变压器低压侧联络节点的三相潮流计算模型，计算得到负荷节点的节点电压和节点相角；
步骤S40，判断空调负荷模型和三相潮流计算模型是否均收敛，若是则当前节点电压和节点相角为最终解，否则将当前节点电压设定为迭代初值，k＝k+1后进入步骤S20。
所述步骤S10的总负荷模型在本实施例中，居民低压配电台区的实际负荷主要以空调为主的动态负荷和以照明为主的静态负荷，因此总负荷模型采用两者并联的方式。
空调负荷模型包括有两种运行状态，分别为空调稳定运行状态和空调启动状态。
作为一个具体实施例，本发明的空调负荷模型采用双笼型转子感应电动机等值电路模型。
参看图1，空调负荷模型为稳定运行状态时，空调负荷模型如下描述：
其中为空调负荷接入配电网处的节点电压；rA、xA和xm分别为空调感应电动机的定子电阻、定子电抗和激磁电抗；r1、x1、r2和x2分别表示下笼转子电阻、互漏电抗、上笼转子电阻和上笼转子电抗；s为感应电动机的转差率；KH为容量折算比，即系统容量基值与感应电动机自身容量之比，表达式如式(2-1)所示：
KH＝SB/SBM (2-1)
空调负荷从配电网吸收的功率可由式(2-2)～(2-3)计算得到：
Rim(s)+jXim(s)＝KH(rA+jxA)+jKHxm//[jKHx1+(KHr1/s)//(KHr2/s+jKHx2)] (2-2)
其中：Rim+jXim是接入网络节点的空调感应电动机负荷的等值阻抗，是关于转差率s的函数；Pim+jQim是空调负荷从配电网吸收的功率值。
由上述分析可得，若要求得该负荷节点的注入功率，则需要得到转差率s的数值，可结合稳定运行时感应电动机的转矩平衡方程计算得到，转矩平衡方程如下所示[2]：
Pm-Pe＝0 (2-4)
其中，Pm和Pe分别为感应电动机在系统基准下的机械负载功率和电磁功率。
对于居民家用空调，压缩机带动的机械负载功率Pm常表示成转速的3次函数，如下：
Pm＝KL[(1-s)3]/KH (2-5)
其中，KL为负载率。
当感应电动机稳定运行时，双笼型转子的电磁功率可以等效为两个转子的电阻吸收功率。下文对电磁功率计算方式进行具体说明。
首先，通过对图1的虚线左边部分进行戴维南等值，则其等值电势为等值阻抗为Req+jXeq，表达式如下所示：
Req+jXeq＝KH(rA+jxA)//jKHxm (2-6)
因而图1中电压的表达式如下：
因此，下笼转子和上笼转子吸收的电磁功率表达式如下：
而转子吸收的总电磁功率表达式如下：
Pe＝Pe1+Pe2 (2-10)
因此，当空调负荷稳定运行时，只要空调负荷接入配电网处的节点电压已知，就可以采用牛顿-拉夫逊法，根据式(2-4)～(2-10)计算得到转子转差率s的大小，进而求得空调负荷节点从配电网吸收的功率大小。所述牛顿-拉夫逊法是常用的解非线性方程组的方法，将关于变量x的非线性方程通过泰勒级数展开，转变为关于变化量Δx的线性方程，再经过反复迭代，以求得足够精确的解。
空调负荷模型为空调启动状态时，空调负荷模型如下描述：
图1所示的为空调负荷稳定运行时的等效电路，完整的双笼型转子的等效电路中下笼中也存在漏电抗，但由于在稳定运行中，下笼漏电抗数值较小，故忽略不计。当电动机刚启动时，转子频率较高，下笼产生较大的漏电抗，故转子电流主要集中在漏电抗较小的上笼转子，因此，在对启动状态的计算中，忽略下笼转子，单独考虑上笼转子。由于电动机启动时，转差率s取值为1，从而得到启动状态下节点负荷的等值阻抗和功率表达式如下：
Rims+jXims＝KH(rA+jxA)+jKHxm//(jKHx1+KHr2+jKHx2) (2-11)
以上计算中空调负荷的参数取值如下表1所示：
表1空调负荷有关参数
在本实施例中的静态负荷模型如下：
照明为主的静态负荷采用幂函数模型表示，如下：
其中，Ps、Qs分别为静态负荷接入网络的实际有功功率与无功功率，V为接入点的实际电压幅值，VN为额定电压，PN和QN为额定条件下负荷的有功功率和无功功率。
在本实施例中的总负荷模型如下：
结合上述空调负荷模型和静态负荷模型，低压配电网中负荷节点的总负荷功率为空调负荷与静态负荷之和，表示如下：
其中，PL和QL分别为负荷节点吸收的总有功功率和总无功功率。
作为一个具体实施例，所述步骤S30的低压配电网三相潮流计算模型如下：
假定低压配电网有n个节点，配电变压器高压侧为第n个节点，选配电变压器高压侧为平衡节点，其abc三相电压幅值和相角给定。三相潮流计算模型中，负荷节点的abc三相采用功率平衡方程，0相采用基尔霍夫电流平衡方程，如下：
其中，Pi,p和Qi,p分别表示第i个节点第p相从电网吸收的有功功率和无功功率，Vi,p表示第i个节点的第p相电压，Vj,p表示第j个节点的第p相电压，Zij表示第i个节点与第j个节点之间的abc相线路的阻抗，Zij,0表示第i个节点与第j个节点之间的0相线路的阻抗。
三相潮流计算模型中，配电变压器低压侧联络节点的abc三相和0相都采用基尔霍夫电流平衡方程，如下：
其中，表示配电变压器高压侧对低压侧的变比，为的共轭函数，θ为配电变压器高压侧对低压侧的电压移相角，Vn,p表示配电变压器高压侧的第p相电压，Vi,p表示配电变压器低压侧第p相电压，Vj,p表示非配电变压器侧第j个节点第p相电压，Zni表示配电变压器的等值阻抗，Zij表示配电变压器低压侧与第j个节点之间的线路阻抗，Zm表示配电变压器低压侧中性点对地阻抗。
利用上述数学模型，分离实部与虚部后一共可以建立8(n-1)个平衡方程，若第n个节点为平衡节点，节点1～(n-1)中每个节点有abc三相及0相的电压与相角8个变量，共8(n-1)个变量，因此总变量个数与方程个数相同，可采用牛顿-拉夫逊算法进行居民低压配电网台区三相潮流计算。
所述步骤S20至步骤S30是的实现步骤，在本实施例中具体如下：
根据上述的三相潮流计算模型与空调负荷功率的计算模型，采用交替迭代算法反复迭代，先给定各个节点的电压初值，通过对式(2-4)～(2-14)采用牛顿-拉夫逊算法计算得到节点吸收功率，再通过式(2-15)和(2-16)采用牛顿-拉夫逊求得配网各节点电压，直至负荷转差率方程与潮流计算均收敛，图3为三相潮流交替迭代算法的流程图所示，具体如下所述：
步骤S21、获取居民低压配电台区的结构和运行数据；
步骤S22、给定各个节点的电压初值，根据空调负荷的运行状态，通过式(2-1)～(2-10)计算得到负荷节点中空调稳定运行状态的电动机转子转差率和abc三相功率，或者通过式(2-11)～(2-12)计算得到负荷节点中空调启动状态的abc三相功率，通过式(2-13)计算得到各个负荷节点中静态负荷的abc三相功率，最后通过式(2-14)计算得到节点abc三相功率；
步骤S31、通过式(2-15)和(2-16)求得配网各节点abc三相及0相电压及相角，判断节点电压是否满足收敛条件：
其中ξ为预先给定的小正数，和为本次和上一次迭代的节点电压；
步骤S41、返回到步骤S22和步骤S31交替迭代，直到居民低压配电台区三相潮流计算收敛，从而获得各节点电压和各支路功率的abc三相及0相的分布情况。
步骤S50、计算台区网损、配电变压器负载率和各个节点的电压不平衡度等指标，反映台区运行的三相不平衡程度。
下面举例说明本发明的考虑负荷特性的居民配电台区三相不平衡潮流计算方法的有益效果，为了验证考虑负荷特性的居民低压配电台区三相不平衡潮流计算模型正确性，选取某个实际居民低压配电台区简化后的7节点低压配网进行测试，配电台区接线如图4所示，低压配电台区的三相接线图如图5所示。
如图4中该低压配电台区共有7个节点，其中，配电变压器容量为SN＝200kVA，变比为10/0.4kV，并采用Dyn11接法，即高压侧对低压侧的电压移相角θ＝-π/6。变压器高压侧节点7为平衡节点，变压器低压侧为节点6，节点1、2、3位于同一条馈线上，节点4、5位于同一条馈线上。根据实际配电网数据，得到各线路的abc三相参数相等，如下表2所示，而0相线路参数为单相参数的1.5倍。另外，根据配电变压器型号S9-200，查找相关的铭牌参数，可计算得到电阻RT＝6.25Ω，电抗XT＝20Ω，计算中忽略了变压器励磁支路的电导和电纳。
表2各线路A、B、C三相参数表
在节点1～节点5对应的相序分别接入空调负荷与基于幂函数模型的静态负荷。同一节点的空调负荷与静态负荷的容量分别占该节点负荷总容量的60％与40％；其中，空调负荷的功率因数取0.85，额定功率与额定容量如下表3所示；静态负荷的功率因数取0.9，额定条件下的有功功率和无功功率如下表3所示。
表3各节点接入的空调负荷和静态负荷功率值
平衡节点7的abc三相电压幅值根据实际运行情况都设定为1.03p.u.，a相、b相和c相电压相角分布设定为0、-2π/3和-4π/3。为了分析考虑负荷特性的居民低压配电台区的三相不平衡潮流计算，分别对以下三种运行工况进行潮流计算：
1)所有节点空调负荷均稳定运行；
2)节点5空调负荷启动，其他节点空调负荷稳定运行；
3)节点4和节点5空调负荷同时启动。
计算得到各节点电压标幺值和空调负荷电动机转差率如表4所示。
表4不同运行状态下各节点电压与转差率
由上表4可以看出，由于各个节点接入的均为单相负荷，系统处于三相不平衡状态，各个节点的三相电压不相等，并且产生0相电压。相对于各节点空调稳定运行时的电压幅值，当节点5空调单独启动、以及节点4与节点5的空调负荷同时启动时，节点1～6的A相电压均明显降低，并且节点4与节点5同时启动状态节点1～6的A相电压幅值小于节点5单独启动状态的电压幅值。由于改变运行状态的空调负荷均在A相，因此对B、C相电压幅值和相角的影响较小。由于三相负荷不平衡，0线在各个节点处存在电位差，因而0线中有电流和功率损耗。节点4和5的空调运行状态变化对另一回馈线上的节点1～3的各相电压幅值与转差率的影响较小。
同时，由三相潮流计算结果得到3种运行工况下的配电台区线损与配电变压器负载率如表5所示，节点的电压不平衡度如表6所示。由表5可以看出，相对于各节点空调负荷稳定运行时，当节点5空调负荷单独启动、节点4与节点5的空调负荷同时启动，台区的线损、线损率、和配电变压器三相负载率增加明显。由于节点某一相空调负荷启动时的功率远大于稳定运行时的功率，因此配电变压器的该相负载率会明显增大，从而导致台区线损与线损率的大幅增加，而其他相的负载率则变化较小。由表6可以看出，当各节点空调负荷处于稳定运行状态时，各节点电压不平衡度均未超出国标规定的电网正常运行的允许值2％的上限。而当节点5单独启动、节点4与节点5的空调负荷同时启动，节点4、5的电压不平衡度均超出国标规定的电网正常运行的允许值上限，此时将会给配电台区的用户供电电压质量带来不利影响。
表5不同运行工况下的配电台区线损与配电变压器负载率
表6不同运行工况下各节点的电压不平衡度
通过考虑负荷特性的居民低压配电台区三相不平衡潮流计算，分析空调动态负荷不同运行状态对于低压配电台区运行的影响，得到如下结论：
1)采用具有空调动态负荷与幂函数静态负荷负荷特性的并联描述居民负荷特性，获得的居民低压配电台区三相不平衡潮流计算结果能够较为准确的反映台区的各种实际运行状态；
2)使用交替迭代法进行考虑负荷特性的居民低压配电台区三相不平衡潮流计算，具有较好的收敛性，易于编程实现；
3)相对于空调负荷稳定运行状态而言，当某个节点某一相空调负荷处于启动状态时，整个台区的线损、线损率均明显增加，并且配电变压器三相负载率与这一相负载率也明显增加；
4)相对于空调负荷稳定运行状态而言，当某个节点空调负荷处于启动状态时，会使该回馈线上各个节点产生较大的电压不平衡度，并且有可能导致节点电压不平衡度超出国标规定的允许上限值，影响用户的正常用电。
上述实施例仅用以说明本发明而并非限制本发明所描述的技术方案；因此，尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明，但是，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或者等同替换；而一切不脱离本发明的精神和范围的技术方案及其改进，其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，包括如下步骤：
步骤S10，建立负荷节点的包括空调负荷模型和静态负荷模型的总负荷模型；所述空调负荷模型是空调负荷从配电网吸收的功率关于转差率的函数，空调负荷模型包括有空调稳定运行状态和空调启动状态；所述静态负荷模型是接入点静态负荷功率关于接入点电压幅值的幂函数模型；所述空调稳定运行状态的空调负荷模型是接入点空调负荷功率关于空调稳定运行状态的转差率的函数；所述空调启动状态的空调负荷模型是接入点空调负荷功率关于空调转差率为1的参数表达式；
步骤S20，先给定各个节点的电压初值，当前为第k次迭代，遍历N个负荷节点，所述k和N均为等于或者大于1的整数，依次判断各负荷节点的空调状态，并计算相应空调状态的接入点空调负荷功率；
步骤S30，依次计算各负荷节点的总负荷功率，并通过以配电变压器高压侧为平衡点的各个负荷节点和配电变压器低压侧联络节点的三相潮流计算模型，计算得到负荷节点的节点电压和节点相角；
所述以配电变压器高压侧为平衡点的各个负荷节点三相潮流计算模型如下：
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>jQ</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>jQ</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>jQ</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <mfrac> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> </mrow> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> </mfrac> <mo>+</mo> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <mfrac> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> </mrow> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <mfrac> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> </mrow> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> </munder> <mfrac> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> </mrow> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mfrac> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中，Pi,p和Qi,p分别表示第i个节点第p相从电网吸收的有功功率和无功功率，Vi,p表示第i个节点的第p相电压，Vj,p表示第j个节点的第p相电压，Zij表示第i个节点与第j个节点之间的abc相线路的阻抗，Zij,0表示第i个节点与第j个节点之间的0相线路的阻抗；
所述配电变压器低压侧联络节点的三相潮流计算模型如下：
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> <mi>n</mi> </munderover> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> <mi>n</mi> </munderover> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> <mi>n</mi> </munderover> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>&amp;Element;</mo> <mi>i</mi> </mrow> </munder> <mrow> <mi>j</mi> <mo>&amp;NotEqual;</mo> <mi>i</mi> </mrow> <mi>n</mi> </munderover> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>j</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>+</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>/</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>+</mo> <mo>&amp;lsqb;</mo> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>a</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>b</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <msup> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>n</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>-</mo> <mover> <mi>k</mi> <mo>&amp;CenterDot;</mo> </mover> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mi>i</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中，表示配电变压器高压侧对低压侧的变比，为的共轭函数，θ为配电变压器高压侧对低压侧的电压移相角，Vn,p表示配电变压器高压侧的第p相电压，Vi,p表示配电变压器低压侧第p相电压，Vj,p表示非配电变压器侧第j个节点第p相电压，Zni表示配电变压器的等值阻抗，Zij表示配电变压器低压侧与第j个节点之间的线路阻抗，Zm表示配电变压器低压侧中性点对地阻抗；
步骤S40，采用交替迭代算法，判断空调负荷模型和三相潮流计算模型是否均收敛，若是则当前节点电压和节点相角为最终解，否则将当前节点电压设定为迭代初值，k＝k+1后进入步骤S20。
2.根据权利要求1所述的考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，所述空调负荷模型采用双笼型转子感应电动机等值电动模型。
3.根据权利要求1所述的考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，所述步骤S20采用牛顿-拉夫逊算法计算相应空调稳定运行状态的转差率进而计算接入点空调负荷功率。
4.根据权利要求1所述的考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，所述静态负荷模型表示如下：
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>P</mi> <mi>N</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>V</mi> <msub> <mi>V</mi> <mi>N</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mn>1.5</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>Q</mi> <mi>N</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>V</mi> <msub> <mi>V</mi> <mi>N</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2.5</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中，Ps、Qs分别为静态负荷接入网络的实际有功功率与无功功率，V为接入点的实际电压幅值，VN为额定电压，PN和QN为额定条件下负荷的有功功率和无功功率。
5.根据权利要求1所述的考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，所述方法还包括有步骤S50，计算反映台区运行的三相不平衡程度的指标，所述指标包括台区网损、配电变压器负载率和各个节点的电压不平衡度。
6.根据权利要求1所述的考虑负荷特性的居民配电台区三相不平衡潮流计算方法，其特征在于，步骤S30采用牛顿-拉夫逊算法求解以配电变压器高压侧为平衡点的各个负荷节点和配电变压器低压侧联络节点各自的三相潮流计算模型，求得配电网各节点电压和节点相角。
CN201610383169.9A 2016-05-31 2016-05-31 考虑负荷特性的居民配电台区三相不平衡潮流计算方法 Expired - Fee Related CN105958483B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610383169.9A CN105958483B (zh) 2016-05-31 2016-05-31 考虑负荷特性的居民配电台区三相不平衡潮流计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610383169.9A CN105958483B (zh) 2016-05-31 2016-05-31 考虑负荷特性的居民配电台区三相不平衡潮流计算方法

Publications (2)

Publication Number Publication Date
CN105958483A CN105958483A (zh) 2016-09-21
CN105958483B true CN105958483B (zh) 2018-06-01

Family

ID=56907476

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610383169.9A Expired - Fee Related CN105958483B (zh) 2016-05-31 2016-05-31 考虑负荷特性的居民配电台区三相不平衡潮流计算方法

Country Status (1)

Country Link
CN (1) CN105958483B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106844895B (zh) * 2016-12-30 2020-06-02 华南理工大学 一种冷热电联供微网能量流的解耦计算方法
CN106920012B (zh) * 2017-03-10 2021-02-19 国网山东省电力公司经济技术研究院 一种含潮流控制设备的参考电网评价方法
CN107968396B (zh) * 2017-11-23 2020-12-22 上海电力学院 一种基于本地相量测量信息的戴维南等值参数计算方法
CN109473983A (zh) * 2018-12-14 2019-03-15 中国航空工业集团公司西安飞机设计研究所 一种飞机交流系统稳态模型及计算方法
CN110601224B (zh) * 2019-09-03 2023-05-23 华南理工大学 基于家庭温控负荷的低压台区电压三相不平衡调节方法
CN114050578B (zh) * 2022-01-07 2022-04-15 国网北京市电力公司 一种配电网在线潮流计算方法、装置、设备及介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0923578A (ja) * 1995-07-03 1997-01-21 Jinbatsu Shiyuukichi 単相三線式配電需用家の構内配電改良方法

Also Published As

Publication number Publication date
CN105958483A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105958483B (zh) 考虑负荷特性的居民配电台区三相不平衡潮流计算方法
Kalcon et al. Small-signal stability analysis of multi-terminal VSC-based DC transmission systems
Liu et al. Transient model of air-conditioner compressor single phase induction motor
CN104362622B (zh) 一种带小水电机组地区电网的负荷模型和建模方法
WO2014089959A1 (zh) 一种基于工程应用的交直流混联大电网的动态等值方法
CN107086591A (zh) 一种基于无功补偿和智能换相的三相不平衡调节器
CN104201671B (zh) 一种含风电的三相不平衡配电网的静态电压稳定性评估方法
CN109921421B (zh) 双馈风电机组输出谐波电流模型建立方法
CN102403716A (zh) 一种多馈入交直流电网动态等值的方法
CN106549402B (zh) Tcsc抑制交直流并联运行系统次同步振荡的方法
CN115358623A (zh) 人工智能的区域电网电能质量配置方法和综合评价方法
CN104993711A (zh) 一种电压暂降过渡过程模拟装置及方法
Ghahremani et al. Synchronous generator third-order model parameter estimation using online experimental data
CN104701864A (zh) 一种采用svc作为无功补偿设备的无功规划方法
CN109586273A (zh) 基于风能hvdc外送系统的次同步振荡评估方法和系统
CN114301055A (zh) 基于宽频量测的电力系统间谐波潮流获取方法与系统
Flack et al. Controller hardware-in-the-loop (C-HIL) testing of decentralized EV-grid integration
Liu et al. Short Circuit Ratio analysis of multi-infeed HVDC system with a VSC-HVDC link
CN106849175A (zh) 一种双馈风电机组撬棒阻值整定方法
Jiang et al. Small-signal modeling and interaction analysis of LCC-HVDC systems based on harmonic state space theory
Shen et al. Modeling and Parameter Identification of the Photovoltaic Inverter based on VSG
CN105897092A (zh) 一种星接三相异步发电机实现单相运行的设计方法
CN111654026A (zh) 一种特高压并联电容器暂态分析方法、装置及存储介质
CN105281319B (zh) 基于加权最小二乘法的发电机参数动态聚合方法
Zhao et al. A General Integration Method for Small-Signal Stability Analysis of Grid-Forming Converter Connecting to Power System

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180601

Termination date: 20190531