专利名称：：基于计算模型拼接的互联系统分布式潮流计算方法
技术领域：
：本发明属电力系统分析领域，更准确地说本发明涉及一种电力系统分布式潮流计算方法。技术背景随着科技的发展和调度集中化的需要，电网日益互联为一体，但是其和电力系统分层分区的调度和控制体制之间存在矛盾。目前各个电力调度中心的EMS多采用独自计算潮流的方式，外网对内网的影响不能准确反映到内网的计算中来，虽然采用了外网等值方法，但计算结果仍不够准确。因此，需要研究如何在潮流计算中将外网对内网的影响准确反映出来。其中一种思路是在多调度控制中心之间进行协调和组织，使用地理上广域分布的各调度中心的数据来计算全网的一体化潮流，从而准确反映外网的电气变化对内网的影响，这是分布式潮流计算的研究内容。具体如何在调度中心之间进行分布式潮流的协调和组织有多种方法，根据分析计算时所采用的模型，可概括为采用集中式建模和采用分布式建模两条思路。各电力调度中心EMS对各自所管辖或关心的电网均进行了建模采用集中式建模的协调和组织工作主要是在将各EMS的分散模型拼接起来建立一个全网统一的一体化模型；而采用分布式建模的协调和组织工作主要是在各EMS之间构造某种与协调算法相对应的外层迭代格式，并据其在各EMS之间进行计算数据的交互与迭代。文献一《电力系统多区域网络模型拼接方法》(中国专利申请号200610166302.1，公开日2007年7月25日)披露了一种电力系统多区域网络模型拼接方法。使用该方法，可以将各区域实际的网络模型拼接成集中式的全模型。虽然文献一并不涉及使用全模型来进行电力系统分析的内容，但得到全模型的目的是使用其进行电力系统的监控和一体化分析。若使用集中式全模型来计算各区域子网的潮流，即称之为采用集中式建模的分布式潮流计算。该方法存在的主要问题是在进行实际网络模型拼接时存在巨大的工作量。此外，由于电网模型变动频繁，在维护中也存在同样的问题。采用分布式建模思路目前所见有如下方法文献二《发输配全局电力系统分析》(电力系统自动化2000年第24巻第1期第17页)在发输配电网一体化分析背景下，披露了一种分布式潮流计算的"主从分裂算法"。在该算法中，各子网地位不平等，有主从(Master-Slave)之分。在外层迭代步中主网一般被看作是从网的广义电源端，而从网一般被看作是主网的广义负荷端。该算法视输电网侧为广义电源端，视配电网侧为广义负荷端，进行了分布式潮流计算。文献三《基于互联网的DTS联合培训仿真系统》(继电器2006年第34巻第2期第70页)披露了一种将如上"主从分裂算法"的思想运用到输电网互联分析中的分布式潮流算法。该算法在计算中将地调侧视为省调侧的广义负荷端来进行输电网背景下的分布式潮流计算。使用该算法求解输电网互联系统分布式潮流的主要问题有如下几点①难以处理地调侧110KV及以下系统合环操作的情况；②形成互联系统潮流计算的统一初始条件时，地调侧子DTS系统自身的不平衡功率难以处理；③该算法需对各子DTS系统核心潮流算法进行改动，从而降低了其对不同厂家平台的适应性。文献四《联合电力网潮流计算的分解协调法》(武汉水利电力大学(宜昌)学报1997年第19巻第1期第35页)披露了一种基于无功和电压的关系来构造外层迭代格式的分布式潮流分解协调算法。文献五《基于异步迭代的多区域互联系统动态潮流分解协调计算》(电力系统自动化2003年第27巻第24期第1页)及其补充文献六《分布式潮流计算异步迭代模式的补充和改进》(电力系统自动化2007年第31巻第2期第12页)披露了一种使用戴维南等值理论来构造外层迭代格式的分布式潮流分解协调算法。有别于主从式算法，这两种算法不从子网的特殊性质和结构出发去看问题，计算中各子网地位平等，其基本思路是各子网独立进行潮流计算后，在外层迭代中构造边界点状态变量的不动点迭代格式，迭代至子网间的边界点状态变量差值达到收敛精度，力图使计算结果和一体化计算结果相同。其共同存在的主要问题是算法组织模式较复杂，且其计算结果存在一定误差。
发明内容为解决现有技术的不足，本发明的目的在于提供一种既能避免集中式建模或者进行实际网络模型拼接时的巨大工作量，又能得到和采用集中式建模计算时完全一致的计算结果，同时在分布式潮流计算中避免相互外网等值，不在子网间进行协调和迭代，简化其组织模式，从而在工程上易于实施的基于计算模型拼接的互联系统分布式潮流计算方法。为了实现上述目的，本发明是采取以下的技术方案来实现的一种基于计算模型拼接的互联系统分布式潮流计算方法，包括下列步骤(1)、定义各区域电网之间的联络支路，所述的区域电网简称为子网，建立《联络支路定义表》，并定义各子网的边界拼接节点集合、内节点集合、外节点集合；(2)、各子网分别形成包括内节点集合和边界拼接节点集合的节点导纳矩阵、节点功率注入向量和节点类型信息向量；(3)、选择某一子网为计算主网，简称为主网，其余子网为计算从网，简称为从网；(4)、主网将自身的节点导纳矩阵、节点功率注入向量和节点类型信息向量和各从网的相对应的矩阵和向量拼接起来，分别得到全局节点导纳矩阵，全局节点功率注入向量和全局节点类型信息向量，由这三种全局信息量得到全局潮流计算方程组；所述的节点导纳矩阵、节点功率注入向量和节点类型信息向量称为计算模型；(5)、各子网在计算模型拼接过程中生成各自的节点编号对照表，或者在主网侧生成所有子网的节点编号对照表；(6)、统一选择解算全局潮流计算方程组的电压迭代初始值，或者伴随计算模型的拼接过程，通过拼接各子网电压迭代初始值向量的方法得到该值；(7)、在主网侧完成全局潮流计算方程组的求解，得到全局一体化潮流结果；(8)、主网侧按全传送模式或者部分传送模式将从网的潮流数据传回；若所有子网的节点编号对照表都生成在主网侧，本步骤中包括主网侧将各从网的节点编号对照表传回相对应的从网。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于两个子网间的全部联络支路恰能对两个子网进行直接的电气分割。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的联络支路选择为在EMS的网络模型中阻抗非零的线路或者变压器。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于当子网间的边界设备为三巻变时，所述的联络支路选择为其某一侧的等值支路。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的各子网分别形成的节点导纳矩阵、节点功率注入向量和节点类型信息向量中不包括外节点集合中的节点。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的全局导纳矩阵的拼接，总体上分两步进行先通过阶数扩充将各子网节点导纳矩阵拼接起来，再在该矩阵中加入联络支路，根据加入联络支路对该矩阵的所产生的影响对该矩阵进行修正。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于各子网节点编号对照表，其中记录了该子网除去外节点集合之外的所有节点在拼接过程中的不同编号之间的对应关系。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于在所述的全传送模式下，主网向从网传送回全部所需潮流数据，包括节点编号、节点电压幅值、相角、节点注入功率、支路潮流，子网自身不进行任何潮流计算。前述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于在所述的部分传送模式下，主网仅向从网传送回节点编号、节点电压幅值、相角，支路潮流和其他电气量计算由子网自己完成。本发明的有益效果和优点采用本发明的方法进行分布式潮流计算，解决了既有方法存在的问题。首先，采用本发明的算法无需全网一体化建模或进行实际的模型拼接，从而避免了由此带来的大量的模型拼接和此后的维护工作。并且，采用本发明的算法计算分布式潮流在子网间也无任何外层协调迭代计算，从而避免了复杂的算法组织模式。最后，采用本发明的算法能得到和采用集中式建模时完全一致的计算结果。总之，本发明的算法将集中建模计算和目前所提出的各种分布式协调迭代算法的优点结合起来，是一种准集中式的分布式潮流算法，既不进行实际的模型拼接，也不在子网间进行协调与迭代，而其计算的结果和全网一体化计算完全一致。图1是以主网S、从网&两个子网为例的本发明采用的基于计算模型拼接的互联系统分布式潮流计算流程框图。图2是用于说明计算模型拼接方案的两个子网模型拼接示意图。图3是以IEEE9节点系统为原型系统的计算模型拼接示意图。图4是忽略变压器励磁回路时的变压器的等值电路图。具体实施方式以下结合附图对本发明作具体的介绍如下在本发明中，披露了将节点导纳矩阵、节点注入功率向量和节点类型信息向量(称这三种信息量为计算模型)拼接得到全局潮流计算方程组的方法和组织模式。极坐标下的潮流方程组如式(l)所示<formula>formulaseeoriginaldocumentpage8</formula>式中"l，...,W，7'e/表示Z号后的标号为/的节点必须直接和节点/相联，并包括7=/的情况，f表示节点/的有功注入功率，0,表示节点/的无功注入功率，t/,表示节点/的电压幅值，f/,表示节点/的电压幅值，A和^表示节点Z和节点7在节点导纳矩阵中对应的互电导和互电纳，是节点/和节点/两节点电压的相角差。用常见的Newton法或FastDecoupled法来求解潮流方程组时，所需全网信息矩阵可分为以下3种节点导纳矩阵、节点初始注入功率列矩阵/向量、节点类型信息列矩阵/向量。已知这3个矩阵加上选定的迭代初值，即可求解常规潮流。在分布式建模的背景下，若能将各子网的这3个矩阵拼接起来，便得到可进行一体化解算的全局潮流计算方程组，其计算结果将和采用集中式建模时的结果完全一致。迭代初值可通过拼接得到，也可在全局潮流机算方程组拼接完成后再统一选择。在分布式潮流计算中，一般选择在具有最丰富计算资源的子网侧来完成拼接与计算，称为计算主网，简称为主网；其余子网称为计算从网，简称为从网。计算模型拼接的实质是将各子网潮流方程组拼接为全局潮流计算方程组，计算模型拼接方法包括节点导纳矩阵拼接方法、节点编号对照表生成方法、节点功率注入向量拼接方法和节点类型信息向量拼接方法四方面的内容。在一体化潮流计算完成后，由于拼接和计算在主网侧完成，主网侧潮流结果自然获得，但从网潮流结果需要由主网传送回去，设计有两种模式①全传送模式主网向从网传送回全部所需潮流数据，主要包括节点编号、节点电压幅值、相角、节点注入功率、支路潮流，子网自身不进行任何潮流计算；②部分传送模式主网仅向从网传送回节点编号、节点电压幅值、相角，支路潮流和其他电气量计算由子网自己完成。全传送模式优点在于子网无需计算，缺点在于数据传输量大；部分传送模式优缺点与之相反。传送回从网的节点编号一般已还原为拼接前的编号，从网获取后可直接使用，编号转化工作亦可由从网自身完成。此外，从网和主网有重复建模的部分，这部分的潮流数据也要传送回从网，其数据对应关系用主、从网统一命名设备ID的办法来解决。以子网S，为主网、以子网52从网为示例的算法组织模式见图1。在本发明中，披露了将各子网节点导纳矩阵拼接起来得到全局节点导纳矩阵的方法。虽然本发明中并不进行任何实际的网络模型拼接，但为了直观说明抽象的节点导纳矩阵的拼接方案，这里首先对计算模型拼接的整体方案在实际网络模型拼接这一层面上说明如下。图2是两个子网模型拼接示意图。国内各调度中心的电力系统应用软件普遍存在着建模部分重复的情况，图2表示了两个子网存在重复建模时的一般情况，图中S表示全网，S,、^为全网S的两个子网，图2(a)左、右图分别表示子网&、A的完整模型，S。表示两者重复建模的部分，n,、.、、为节点，L和是S。中的两条支路且恰能将S和S2进行电气分割，即/,,和是S和&的联络支路。本文以图2为例，提出模型拼接方案如下(1)定义子网间联络支路，其应恰能对子网进行电气分割。图2(a)中，假定/,,和被定义为联络支路。近些年来有关垂直通孔的研究R益引起人们重视的原因在于，随着特征尺寸的减少，金属线网的电阻在不断增大，而其中增长速度最快的就是垂直通孔。在特征尺寸为90纳米的芯片中，一个垂直通孔的电阻值可以相当于IOO微米的金属线。垂直通孔的电阻值在整个多层线网占据的比重越来越大，这对芯片设计的各个阶段都提出了新的挑战。特别地，这一现象对多层供电网络的设计具有尤其重要的影响。附图1中显示的是一种应用广泛的供电网络结构。在网络中，同一层的金属线都是相互平行，而相邻层的金属线相互垂直，一般情况下，高层的金属线具有较大的宽度(width)同时相邻金属线间的间距(pitch)也相对较大。在这种多层的供电网络中，自然需要通孔来实现不同层间金属线的互连。如附图1所示，为了不占用额外的布线资源，对信号线的布线造成影响，通孔一般都设置在相邻层间金属线的交叉点(intersection)上。垂直通孔电阻值的增加必然也会带来供电网络上电压降的增大。解决这一问题的基本方法就是增加垂直通孔的数目，如附图1所示，在有些金属线交叉点处就放置了多个通孔。当然，放置太多的垂直通孔也会造成制造成本的增加，并且受几何尺寸等因素影响，在各个交叉点处可以放置的通孔个数也有最大值限制。这样，随之而来的一个问题就是如何更好地在交叉点处配置通孔的数目。最简单的方法自然是在同一层的交叉点处分配相同数目的通孔，但是很显然，对大多数电路来讲这种方法并不是一个最佳的方法。与线宽优化类似，优化配置通孔的问题也可以抽象为一个有约束的数学优化问题。但是由于各个交叉点处只能配置整数个通孔，该问题是一个离散的数学优化问题，因此，直接用数学规划的方法求解将非常困难。到目前为止，尚没有提出一种成熟的方法解决这个问题。
发明内容针对目前尚没有一种能对多层供电网络中的垂直通孔进行优化配置以降低网络中电压波动的方法的现状，本发明提出了一种快速高效的优化配置垂直通孔的启发式方法。其特征在于该方法是一种以灵敏度信息作为指导，采用循环累加的方式实现多层供电网络中垂直通孔的优化配置的启发式方法，优化的目标是更加合理地配置各个金属线交叉点处垂直通孔的数目，以消除多层供电网络最底层节点中电压降超过规定最大允许值的违规节点，减少最底层节点的最大电压降。该方法是在计算机中依次按以下歩骤实现的-从网，拼接得到阶数扩充后的导纳阵(为<formula>formulaseeoriginaldocumentpage11</formula>(2)式中zeros(m，n)表示^x"阶的零矩阵，zeros(n，m)类推。易知，导纳阵(是m+"阶的方阵。需要说明，主网导纳阵和从网导纳阵在拼接时的位置不一定分别在左上角和右下角，可以任意放置，文中加以限制是为了能够清楚加以表述。(3)在导纳阵C中增加与联络支路对应的互导纳元素，并修正原各子网边界拼接节点对应的自导纳如下a)在导纳阵(中增加与联络支路对应的互导纳元素。图3(c)中，以联络线/,为例，假定其两端节点在导纳阵(中对应行号(亦即编号)为e、f，其导纳为;v，对地导纳为;v。，则在t中添加2个非对角元素如式(3)。联络线/2类同。<formula>formulaseeoriginaldocumentpage11</formula>(3)b)在导纳阵t中修正各边界拼接节点对应的自导纳。以拼接后节点、",为例，对应的自导纳&_和",对应的自导纳rm。#有如下增量<formula>formulaseeoriginaldocumentpage11</formula>(4)至此，.K。拼接完毕。c)对导纳阵C各行所对应节点进行编号优化，即将导纳阵《。进行相应行列对调初等变换得到全局导纳矩阵i;。。当联络支路为联络变而非联络线时，拼接思路并无不同，只是式(3)、(4)或有所变化。忽略变压器励磁回路时，变压器支路可用其漏抗串联一个理想变压器来模拟，如图4(a)所示。其n型等值电路如图4(b)所示。其中，^是变压器漏抗，JV=l/^，尺为变压器变比，节点e和节点/为变压器等值电路两侧端点，^和^v分别是节点e和节点/的端电压向量，h和)/分别是节点e和节点/的注入电流向量。据图4，易知通过联络变拼接导纳阵时，式(3)形式变化为<formula>formulaseeoriginaldocumentpage11</formula>而式(4)中修正量变化为<formula>formulaseeoriginaldocumentpage12</formula>(6)在本发明中，披露了一种建立节点编号对照表的方法。伴随着对导纳阵C进行优化编号得到全局导纳矩阵L的过程，应对每个子网都建立节点编号对照表,表中记录了各子网除去外节点之外的所有节点在拼接过程中的编号对应关系。以图1中主网S,为例，表l(a)即其节点编号对照表。拼接前主网S,的节点编号自l递增至w，形成与导纳阵&的行号对应的列向量A，放在第l歹!h拼接成导纳阵(后，得到与其行号对应的节点编号列向量"，，放在第2列；进行编号优化得到全局导纳矩阵乙后，得到与其行号对应的节点编号列向量。,，放在第3列，建表完成。从网&建表方法类同，但由于导纳阵&在导纳阵C的右下方，节点编号列向量。2和《2并不相同，其每个分量都加入了m(m为主网《的节点数，称为拼接偏移量，下文同)，如表l(b)所示。表1s,与&节点编号对照表(附图l)ff2"，111l+m乂22x222+m3333+m附附/+附(a)(b)在本发明中，还披露了一种将各子网节点注入功率向量拼接起来得到全局节点注入功率向量的方法。假设主网s的潮流方程组节点注入功率列向量为4，与导纳阵&,对应；从网52相应列向量为&，则拼接可分两步进行(1)将节点注入功率列向量^、A进行维数扩充拼接，得到与导纳阵《。对应的未经编号优化的全局节点注入功率向量&。如下&。=^(7)(2)查节点编号对照表，根据节点注入功率列向量A、A原节点在拼接前后的编号对应关系对&。进行相应的行调换，使各分量行号和原所对应节点优化后的编号一致，便得到与全局导纳矩阵乙对应的全局节点注入功率向量&。。或者亦可在对导纳阵(进行优化编号的同时，对&。进行相应调换得到全局节点注入功率向量^。。在本发明中，还披露了一种将各子网节点类型信息向量拼接起来得到全局节点类型信息向量的方法。在潮流计算中电力系统的节点分为P2节点、i^节点和平衡节点三种类型。若拼接前主网S中节点类型信息向量为惟一的平衡节点编号放在一元平衡节点号列向量t;,,中，P2节点号放在其尸2节点号列向量^.,中，^节点号放在尸F节点号列向量；，中，从网&中相应向量为<formula>formulaseeoriginaldocumentpage13</formula>节点类型信息向量的拼接按如下步骤进行。(i)将列向量7^2、rw2、r—各分量均加入拼接偏移量m。由于全网仅需保留一个平衡节点，从网&的平衡节点在拼接完成后将转化为全局潮流计算中的尸r节点，所以取<formula>formulaseeoriginaldocumentpage13</formula>。(2)取<formula>formulaseeoriginaldocumentpage13</formula>均和导纳阵t相对应。(3)由主网s、从网&节点编号对照表，査得列向量t;,'、t;:、；'各分量优化编号后的新编号，各自组成和全局导纳矩阵L相对应的全局平衡节点号列向量z;,、全局PS节点号列向量^和全局/^节点号列向量7;,这三个全局节点类型信息向量，拼接完成。或者，全局节点类型信息向量亦可用先查表后拼接的方式来得到。实施例下面是本发明的一个实施例，包括了采用本发明的算法以IEEE9节点系统为原型系统进行的分布式潮流计算的过程和结果。本发明的其它的特征、目的和优点也可以从实施例的说明和附图中看出。附图3说明了以IEEE9节点系统为原型系统使用本发明的算法的过程及结果比对，其网络和计算参数均采用IEEE9节点系统标准数据。图3(b)表示了两个具有重复建模联络支路/,、/2的两个子网5,、S2，各自具有从l开始的节点编号。实施例步骤如下。(1)在《联络支路定义表》中定义/,、/2为子网间的联络支路；根据联络支路的定义和本发明的技术方案，定义子网S、^的边界拼接节点集合^、内节点集合乂、外节点集合A^。S边界节点集合A^="2},M.,("3，"4,"5)，iV￡1={",，}，&边界节点集合^,"2}，A^={"3,"4}，"}。(2)子网S,和&各自拓扑到边界拼接节点，各自形成节点导纳矩阵、节点注入功率向量和节点类型信息向量。子网5;和&各形成导纳阵^和&.2。导纳阵^阶数为5，导纳阵^阶数为4。子网S的潮流方程组节点注入功率列向量为4，与导纳阵&,对应，子网&相应列向量为&。子网S中节点类型信息向量为惟一的平衡节点编号放在一元列向量7;,1中，尸0节点号放在列向量r一中，尸r节点号放在列向量r—中。子网&中相应向量为r<2、7;9.2、t;v2。(3)选择子网&为主网，子网&为从网。(4)主网s，将自身的和接收从网&传送过来的节点导纳矩阵、节点功率注入向量和节点类型信息向量(称这三种信息量为计算模型)按照本发明的技术方案拼接起来，最终得到全局导纳矩阵乙，全局节点注入功率向量&。，全局节点类型信息向量z;,、rP9、rpv，完成计算模型的拼接工作。其过程在直观层面上如图3(c)、(d)所示。其中，图3(d)中节点编号由半动态优化法得到。(5)伴随着计算模型拼接过程，子网s,和&各自生成节点编号对照表，分别见表2(a)、(b)。(6)统一选择全局潮流计算方程组电压迭代初始值。(7)由以上拼接过程得到了全局潮流计算方程组，在主网S,侧完成该全局潮流方程组的求解，得到全局一体化潮流计算结果。关键计算结果量全网电压相量记为f)，见表3。(8)主网S,侧按全传送模式将从网S2的潮流数据传回到从网S2。结果验证为了说明所得到的计算结果的正确性，我们再在如图3(a)所示IEEE9节点系统标准数据下进行常规潮流计算，得到全网电压相量记为W，将这个计算结果看作是该实施例分布式潮流计算结果的真值，见表3。注意到图3(a)、(d)中的节点编号对应关系，可以看到，O和W完全一致，即采用本发明算法得到的节点电压和真值完全一致。表2S,与&节点编号对照表<table>tableseeoriginaldocumentpage15</column></row><table>(a)(b)表3拼接算法潮流结果与真值(标幺值)比较<table>tableseeoriginaldocumentpage15</column></row><table>注A是拼接后优化后编号，见附图3(d)。2是IEEE9节点系统编号，见附图3(a)。本发明已以较佳实施例公布如上，但上述实施例不以任何形式限定本发明，凡采用等同替换或等效变换的形式所获得的技术方案，均落在本发明的保护范围之内。权利要求1、基于计算模型拼接的互联系统分布式潮流计算方法，包括下列步骤(1)、定义各区域电网之间的联络支路，所述的区域电网简称为子网，建立《联络支路定义表》，并定义各子网的边界拼接节点集合、内节点集合、外节点集合；(2)、各子网分别形成包括内节点集合和边界拼接节点集合的节点导纳矩阵、节点功率注入向量和节点类型信息向量；(3)、选择某一子网为计算主网，简称为主网，其余子网为计算从网，简称为从网；(4)、主网将自身的节点导纳矩阵、节点功率注入向量和节点类型信息向量和各从网的相对应的矩阵和向量拼接起来，分别得到全局节点导纳矩阵，全局节点功率注入向量和全局节点类型信息向量，由这三种全局信息量得到全局潮流计算方程组；所述的节点导纳矩阵、节点功率注入向量和节点类型信息向量称为计算模型；(5)、各子网在计算模型拼接过程中生成各自的节点编号对照表，或者在主网侧生成所有子网的节点编号对照表；(6)、统一选择解算全局潮流计算方程组的电压迭代初始值，或者伴随计算模型的拼接过程，通过拼接各子网电压迭代初始值向量的方法得到该值；(7)、在主网侧完成全局潮流计算方程组的求解，得到全局一体化潮流结果；(8)、主网侧按全传送模式或者部分传送模式将从网的潮流数据传回；若所有子网的节点编号对照表都生成在主网侧，本步骤中包括主网侧将各从网的节点编号对照表传回相对应的从网。2、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于两个子网间的全部联络支路恰能对两个子网进行直接的电气分割。3、根据权利要求1或2所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的联络支路选择为在EMS的网络模型中阻抗非零的线路或者变压器。4、根据权利要求3所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于当子网间的边界设备为三巻变时，所述的联络支路选择为其某一侧的等值支路。5、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的各子网分别形成的节点导纳矩阵、节点功率注入向量和节点类型信息向量中不包括外节点集合中的节点。6、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于所述的全局导纳矩阵的拼接，总体上分两步进行先通过阶数扩充将各子网节点导纳矩阵拼接起来，再在该矩阵中加入联络支路，根据加入联络支路对该矩阵的所产生的影响对该矩阵进行修正。7、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于各子网节点编号对照表，其中记录了该子网除去外节点集合之外的所有节点在拼接过程中的不同编号之间的对应关系。8、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于在所述的全传送模式下，主网向从网传送回全部所需潮流数据，包括节点编号、节点电压幅值、相角、节点注入功率、支路潮流，子网自身不进行任何潮流计算。9、根据权利要求1所述的基于计算模型拼接的互联系统分布式潮流计算方法，其特征在于在所述的部分传送模式下，主网仅向从网传送回节点编号、节点电压幅值、相角，支路潮流和其他电气量计算由子网自己完成。全文摘要本发明涉及一种基于计算模型拼接的互联系统分布式潮流计算方法，包括定义各子网之间的联络支路和边界节点集合；各子网分别形成信息量；选择主网和从网；主网和从网的信息量拼接得到全局信息量，从而得到全局潮流计算方程组；在主网侧生成或各子网生成各自的节点编号对照表；统一选择电压迭代初始值；在主网侧求解，得到全局一体化潮流结果；主网侧按全传送模式或者部分传送模式将从网的潮流数据传回，完成分布式潮流计算。本发明避免了全网一体化建模或进行实际的模型拼接带来的大量的模型拼接和此后的维护工作；避免了复杂的算法组织模式，同时还可以得到和采用集中式建模时完全一致的计算结果。文档编号G06F19/00GK101221599SQ20071019154公开日2008年7月16日申请日期2007年12月12日优先权日2007年12月12日发明者卢建刚,姚建国,锋李,杨胜春,明高申请人:国电南瑞科技股份有限公司;广东电网公司