专利名称：基于虚拟恒电流负荷的环状配电网潮流获得方法
技术领域：
本发明涉及一种基于虚拟恒电流负荷的环状配电网潮流获得方法，应用于电力监控领域。
背景技术：
在配电网监控管理系统中，配电网潮流是高级应用软件提供给调度员和计划人员的最重要的功能，用于分析配电网络电气状态和检查越界范围等，同时，它也是配电网络分析的基础软件，用于电压/无功优化、操作模拟和结线变化分析等。配电网络具有不同于输电网络的显著特点1)配电网基本是放射形结构，或者是接近辐射状网络结构(弱环结构)；2)支路参数!■和X相差不大，有时甚至会出现r > X 的情况；3)三相支路参数不平衡和三相负荷不对称问题比较突出。配电网络的这些特点造成对传统的潮流算法来说实际上属于病态条件，例如由于配电网支路参数r/x比值较大， 使原来在高压输电网中行之有效的算法，如快速解耦法等，在配电网中不再有效；其次，由于配电网中不对称元件，如未换位的输电线路、三相三铁心柱变压器、电力机车负荷及其它三相不对称负荷的存在，使配电网的三相电压、电流不再对称，因而不能像对称系统那样， 只计算单相的情况，而必须进行三相潮流计算。因此，获得配电网潮流必须采用不同于输电网且适合配电网特点的、有效的、鲁棒性好的潮流计算方法。目前常用的三相潮流计算可分为回路阻抗法、前推回代法。S.K Goswami和S. K Basu在1991年提出的回路阻抗法以回路电流为变量，从馈线根节点到每一个负荷点形成一条回路，根据基尔霍夫电压定律，列回路电流方程，该方法处理网孔能力较强，但节点和支路编号处理复杂。前推回代法最初是从手算算法发展而来的，当用来进行辐射状配电网的潮流计算时，该算法的效率是所有算法中最高的，占用内存也很少，前推回代法由于自身的优点在国外已被不断改进且应用于与生产实际中，它的缺点是，当应用于环状网络时需要进行特殊的处理。
发明内容
本发明的目的在于克服上述不足，提供一种能够对环状配电网络进行分析的基于虚拟恒电流负荷的环状配电网潮流获得方法。本发明的目的是这样实现的一种基于虚拟恒电流负荷的环状配电网潮流获得方法，所述方法的步骤为
步骤I、计算机从配电网管理系统中获取配电网的网络参数；
步骤2、形成配电网的树枝集合和连枝集合，其具体形成步骤如下
步骤2. I、通过步骤I获取的网络参数得到节点集合和边集合；
步骤2. 2、置节点集合中所有节点的是否被处理标志位为False (即表示还没有处理过)；初始化堆栈，将节点集合中的根节点压入堆栈，并置根节点的是否被处理标志位为 True，初始化树枝集合和连枝集合为空；步骤2. 3、若堆栈为空，跳转至步骤2. 6，否则执行步骤2. 4 ；
步骤2. 4、从堆栈中弹出节点，该节点为节点一，遍历与该节点相连的边，与该边相连的另一节点为节点二，若节点二的是否被处理标志位为False，则将该边加入到树枝集合，将节点二压入堆栈，并置节点二的是否被处理标志位为True ；
步骤2. 5、跳转至步骤2.3 ;
步骤2. 6、将边集合中除树枝集合外的边放入连枝集合中；
步骤2. 7、遍历每条连枝，利用步骤I得到的网络参数，获得连枝的支路导纳矩阵；
步骤3、重新编号，编号步骤为
步骤3. I、将连枝从边集合中移除；
步骤3. 2、初始化数据
置队列为空，置当前整数为1，置根节点的新编号为1，将节点集合中的根节点放入队列中，置根节点的是否被处理标志位为True ；
步骤3. 3、如果队列为空，跳转至步骤4，否则，从队列的对首取出节点，并将该节点设为当前节点；
步骤3. 4、从边集合中取出与当前节点相连的边，通过改变获得与当前节点相连的节点，将该节点设置为子节点，若子节点的是否被处理标志位为False，则将当前整数累加1， 并将累加后的当前整数置为该子节点的新编号，同时将该子节点存入队列的队尾；
步骤3. 5、跳转至步骤3. 3 ；
步骤4、初始化遍历所有连枝，在有连枝相连的节点上设置一个虚拟恒电流负荷，电流值置为0，使用数组initialV存储节点电压，置数组initialV中元素为对应节点的额定电压，预置ε、ε i分别为足够小的常数，
步骤5、置各节点初始电压为initialV ；
步骤6、对由树枝组成的辐射状网络使用前推回代法计算各节点电压和支路电流； 步骤6. I、按节点编号从大到小的顺序遍历所有节点，设当前节点为currentN,如果 currentN没有父节点，即currentN是根节点，贝U转到步骤6. 2 ;否则,设currentN与父节点之间的支路为br,利用公式I=(S/V)*计算currentNode流入负荷的电流,其中S是负荷功率，V为节点电压，*号表示共轭，br的电流为流入负荷的电流加上以currentN为父节点的树枝上的电流，遍历以currentN为父节点的所有树枝,计算得到br上的电流；
步骤6. 2、置实数deltaV为无穷小，按节点编号从小到大的顺序遍历所有节点， 设当前节点为currentNo ,从节点集中获得currentNo的子节点，遍历currentNo的子节点，设当前子节点为sonNo,设currentNo与sonNo之间支路为bra, Ibra为6. I计算所得bra的电流，设sonNo的当前节点电压为VsonNol,更新sonNo的节点电压为 VsonNo=VcurrentNo-Zbr^Ibra,其中 VcurrentNo 为 currentNo 的电压，Zbra 为支路 bra 的阻抗矩阵，设VsonNo与VsonNol对应相位差值的模最大值为a,如果deltaV〈a,则置 de I 3,ν_&;
步骤6. 3、如果deltaV〈 ε χ，转到步骤7，否则转到6. I ；
步骤7、置最大偏差delta为无穷小,遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2，nl上虚拟恒电流负荷loadl，n2上虚拟恒电流负荷为load2，步骤6得到的计算结果中取出branch两端的电压vl和v2,设branch导纳矩阵为Y,利用公式I=Y* (vl_v2)，计算出支路电流I，设Ioadl电流值为Iloadl, load2电流值为Iload2，设Iloadl与 Iload2对应相位差值的模的最大值为b,如果delta〈b,置delta=b ；
步骤8、如果delta〈e，输出节点电压和支路电流结果，退出；否则，进入步骤9;
步骤9、遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2,更新nl上虚拟恒电流负荷的电流值为I，n2上虚拟恒电流负荷的电流值为-I ；
步骤10、将步骤6计算所得各节点电压存入initialV中；
步骤11、转到步骤5。与现有技术相比，本发明的有益效果是
配电网正常情况下呈现辐射状网络结构，但在某些特殊的运行方式下(如故障处理和网络重构)，可能会出现短时间的环网运行，最常见的是单环网形式。为了计算环状配电网的潮流，本发明提出一种基于补偿电流的环状配电网潮流获得方法。该方法首先分析出环状配电网络的树枝和连枝，接着通过两层迭代来解决潮流问题，内层迭代用前推回代法计算有树枝组成的纯辐射状网络潮流，外层迭代不断修正连枝两端节点的注入电流达到整体收敛；从而实现对环状配电网络的分析计算，有助于提高配电网络的整体监控能力。
图I为本发明基于虚拟恒电流负荷的环状配电网潮流获得方法的流程图。图2为本发明环状测试配电网络的结构图。图3为本发明图2所示环状测试配电网络的树枝部分的结构4为辐射状配电网络局部。
具体实施例方式参见图1，本发明涉及的一种基于虚拟恒电流负荷的环状配电网潮流获得方法，该方法的步骤为
步骤I、计算机从配电网管理系统中获取配电网的网络参数，
步骤2、形成配电网的树枝集合和连枝集合，其具体形成步骤如下
步骤2. I、通过步骤I获取的网络参数得到节点集合和边集合；
步骤2. 2、置节点集合中所有节点的是否被处理标志位为False (即表示还没有处理过)；初始化堆栈，将节点集合中的根节点压入堆栈(节点数据从配电网管理系统中获得时计算机已经知晓哪个节点为根节点)，并置根节点的是否被处理标志位为True，初始化树枝集合和连枝集合为空；
步骤2. 3、若堆栈为空，跳转至步骤2. 6，否则执行步骤2. 4 ；
步骤2. 4、从堆栈中弹出节点，该节点为节点一，遍历与该节点相连的边，与该边相连的另一节点为节点二，若节点二的是否被处理标志位为False，则将该边加入到树枝集合，将节点二压入堆栈，并置节点二的是否被处理标志位为True ；
步骤2. 5、跳转至步骤2. 3 ；
步骤2. 6、将边集合中除树枝集合外的边放入连枝集合中；
步骤2. 7、遍历每条连枝，利用步骤I得到的网络参数，获得连枝的支路导纳矩阵；
步骤3、重新编号，编号步骤为步骤3. I、将连枝从边集合中移除；
步骤3. 2、初始化数据
置队列为空，置当前整数为1，置根节点的新编号为1，将节点集合中的根节点放入队列中，置根节点的是否被处理标志位为True ；
步骤3. 3、如果队列为空，跳转至步骤4，否则，从队列的对首取出节点，并将该节点设为当前节点；
步骤3. 4、从边集合中取出与当前节点相连的边，通过改变获得与当前节点相连的节点，将该节点设置为子节点，若子节点的是否被处理标志位为False，则将当前整数累加1， 并将累加后的当前整数置为该子节点的新编号，同时将该子节点存入队列的队尾；
步骤3. 5、跳转至步骤3.3 ;
步骤4、初始化遍历所有连枝，在有连枝相连的节点上设置一个虚拟恒电流负荷，电流值置为0，使用数组initialV存储节点电压，置数组initialV中元素为对应节点的额定电压，预置ε、ε i分别为足够小的常数，
步骤5、置各节点初始电压为initialV ；
步骤6、对由树枝组成的辐射状网络使用前推回代法计算各节点电压和支路电流； 步骤6. I、按节点编号从大到小的顺序遍历所有节点，设当前节点为currentN,如果 currentN没有父节点，即currentN是根节点，贝U转到步骤6. 2 ;否则,设currentN与父节点之间的支路为br,利用公式I=(S/V)*计算currentNode流入负荷的电流,其中S是负荷功率，V为节点电压，*号表示共轭，br的电流为流入负荷的电流加上以currentN为父节点的树枝上的电流，遍历以currentN为父节点的所有树枝,计算得到br上的电流；
步骤6. 2、置实数deltaV为无穷小，按节点编号从小到大的顺序遍历所有节点， 设当前节点为currentNo ,从节点集中获得currentNo的子节点，遍历currentNo的子节点，设当前子节点为sonNo,设currentNo与sonNo之间支路为bra, Ibra为6. I计算所得bra的电流，设sonNo的当前节点电压为VsonNol,更新sonNo的节点电压为 VsonNo=Vcu;rrentNo-Zb;r*Ibra,其中 VcurrentNo 为 currentNo 的电压，Zbra 为支路 bra 的阻抗矩阵，设VsonNo与VsonNol对应相位差值的模最大值为a,如果deltaV〈a,则置 de I 3,ν_&;
步骤6. 3、如果deltaV〈 ε χ，转到步骤7，否则转到6. I ；
步骤7、置最大偏差delta为无穷小,遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2，nl上虚拟恒电流负荷loadl，n2上虚拟恒电流负荷为load2，步骤6得到的计算结果中取出branch两端的电压vl和v2,设branch导纳矩阵为Y,利用公式I=Y* (vl_ v2)，计算出支路电流I，设Ioadl电流值为Iloadl, load2电流值为Iload2，设Iloadl与 Iload2对应相位差值的模的最大值为b,如果delta〈b,置delta=b ；
步骤8、如果delta〈 ε，输出节点电压和支路电流结果，退出；否则，进入步骤9 ;
步骤9、遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2,更新nl上虚拟恒电流负荷的电流值为I，n2上虚拟恒电流负荷的电流值为-I ；
步骤10、将步骤6计算所得各节点电压存入initialV中；
步骤11、转到步骤5。以图2所示的配电网络为例，具体步骤说明如下步骤I、初始化
从配电网管理系统中获得配电网的网络参数。该配电网络只有馈线和负荷，馈线的支路号、首节点编号、末节点编号和阻抗如表I所示，负荷所在节点以及负荷功率功率如表2 所示。表I馈线参数表表I馈线参数表
权利要求
1.一种基于虚拟恒电流负荷的环状配电网潮流获得方法，其特征在于所述方法的步骤为步骤I、计算机从配电网管理系统中获取配电网的网络参数，步骤2、形成配电网的树枝集合和连枝集合，其具体形成步骤如下步骤2. I、通过步骤I获取的网络参数得到节点集合和边集合；步骤2. 2、置节点集合中所有节点的是否被处理标志位为False (即表示还没有处理过)；初始化堆栈，将节点集合中的根节点压入堆栈，并置根节点的是否被处理标志位为 True，初始化树枝集合和连枝集合为空；步骤2. 3、若堆栈为空，跳转至步骤2. 6，否则执行步骤2. 4 ；步骤2. 4、从堆栈中弹出节点，该节点为节点一，遍历与该节点相连的边，与该边相连的另一节点为节点二，若节点二的是否被处理标志位为False，则将该边加入到树枝集合，将节点二压入堆栈，并置节点二的是否被处理标志位为True ；步骤2. 5、跳转至步骤2.3 ;步骤2. 6、将边集合中除树枝集合外的边放入连枝集合中；步骤2. 7、遍历每条连枝，利用步骤I得到的网络参数，获得连枝的支路导纳矩阵；步骤3、重新编号，编号步骤为步骤3. I、将连枝从边集合中移除；步骤3. 2、初始化数据置队列为空，置当前整数为1，置根节点的新编号为1，将节点集合中的根节点放入队列中，置根节点的是否被处理标志位为True ；步骤3. 3、如果队列为空，跳转至步骤4，否则，从队列的对首取出节点，并将该节点设为当前节点；步骤3. 4、从边集合中取出与当前节点相连的边，通过改变获得与当前节点相连的节点，将该节点设置为子节点，若子节点的是否被处理标志位为False，则将当前整数累加1， 并将累加后的当前整数置为该子节点的新编号，同时将该子节点存入队列的队尾；步骤3. 5、跳转至步骤3.3 ;步骤4、初始化遍历所有连枝，在有连枝相连的节点上设置一个虚拟恒电流负荷，电流值置为0，使用数组initialV存储节点电压，置数组initialV中元素为对应节点的额定电压，预置ε、ε i分别为足够小的常数，步骤5、置各节点初始电压为initialV ；步骤6、对由树枝组成的辐射状网络使用前推回代法计算各节点电压和支路电流； 步骤6. I、按节点编号从大到小的顺序遍历所有节点，设当前节点为currentN,如果 currentN没有父节点，即currentN是根节点，贝U转到步骤6. 2 ;否则,设currentN与父节点之间的支路为br,利用公式I=(S/V)*计算currentNode流入负荷的电流,其中S是负荷功率，V为节点电压，*号表示共轭，br的电流为流入负荷的电流加上以currentN为父节点的树枝上的电流，遍历以currentN为父节点的所有树枝,计算得到br上的电流；步骤6. 2、置实数deltaV为无穷小，按节点编号从小到大的顺序遍历所有节点， 设当前节点为currentNo ,从节点集中获得currentNo的子节点，遍历currentNo的子节点，设当前子节点为sonNo,设currentNo与sonNo之间支路为bra, Ibra为6. I计算所得bra的电流，设sonNo的当前节点电压为VsonNol，更新sonNo的节点电压为 VsonNo=VcurrentNo-Zbr*Ibra，其中 VcurrentNo 为 currentNo 的电压，Zbra 为支路 bra 的阻抗矩阵，设VsonNo与VsonNol对应相位差值的模最大值为a，如果deltaV〈a，则置 de I 3,ν_&;步骤6. 3、如果deltaV〈 ε χ，转到步骤7，否则转到6. I ；步骤7、置最大偏差delta为无穷小,遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2，nl上虚拟恒电流负荷loadl，n2上虚拟恒电流负荷为load2，步骤6得到的计算结果中取出branch两端的电压vl和v2,设branch导纳矩阵为Y,利用公式I=Y* (vl_ v2)，计算出支路电流I，设Ioadl电流值为Iloadl, load2电流值为Iload2，设Iloadl与 Iload2对应相位差值的模的最大值为b,如果delta〈b,置delta=b ；步骤8、如果delta〈 ε，输出节点电压和支路电流结果，退出；否则，进入步骤9 ;步骤9、遍历所有连枝,设当前连枝为branch, branch两端节点为nl和n2,更新nl上虚拟恒电流负荷的电流值为I，n2上虚拟恒电流负荷的电流值为-I ；步骤10、将步骤6计算所得各节点电压存入initialV中；步骤11、转到步骤5。
全文摘要
本发明涉及一种基于虚拟恒电流负荷的环状配电网潮流获得方法，通过两层迭代来解决潮流问题，内层迭代用前推回代法计算有树枝组成的纯辐射状网络潮流，外层迭代不断修正连枝两端节点的注入电流达到整体收敛。本发明基于虚拟恒电流负荷的环状配电网潮流获得方法，能够对环状配电网络进行分析。
文档编号H02J3/00GK102593825SQ20121003804
公开日2012年7月18日 申请日期2012年2月20日 优先权日2012年2月20日
发明者张建兴 申请人:江苏方程电力科技有限公司