CN105022931A - 一种直流电网计算方法及装置 - Google Patents
一种直流电网计算方法及装置 Download PDFInfo
- Publication number
- CN105022931A CN105022931A CN201510489135.3A CN201510489135A CN105022931A CN 105022931 A CN105022931 A CN 105022931A CN 201510489135 A CN201510489135 A CN 201510489135A CN 105022931 A CN105022931 A CN 105022931A
- Authority
- CN
- China
- Prior art keywords
- direct current
- bar
- current bus
- balance equation
- power balance
- 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.)
- Pending
Links
Abstract
本发明提供了一种直流电网计算方法及装置,该方法包括:确定每个直流汇流条线路总输出功率,其中,线路总输出功率是每个直流汇流条向与其相连的所有线路流出的总功率;根据线路总输出功率、负载功率确定每个直流汇流条的功率平衡方程,其中,负载是所述直流电网中的用电设备;根据每个直流汇流条的功率平衡方程得到功率平衡方程组,并通过联合求解功率平衡方程组得到每个直流汇流条的电压。通过本发明,解决了相关技术在计算中存在系统误差并且不便于编程计算的问题,进而达到了提高计算准确性和计算效率的效果。
Description
技术领域
本发明涉及飞机电源系统设计领域,具体而言,涉及一种直流电网计算方法及装置。
背景技术
相关技术中,飞机直流电网计算方法基本流程如下:统计各汇流条负载水平、估算负载电流、计算线路压降;重复上述步骤对每条线路进行计算即可完成整个电网计算。图1是根据相关技术中的直流电路示意图,如图1所示,其主要计算流程如下:
当已知汇流条I的电压幅值为VI、线路电阻为r。首先需要统计出汇流条II上总的负载功率P,此时电网计算目的是获取线路压降以及汇流条II上的电压幅值VII。
(1)忽略线路损耗,估算线路的电流I。一般情况下飞机直流电网额定电压为28V,因此线路中电流大小约为
(2)根据线路电流和电阻计算线路压降ΔV,由欧姆定律可知
ΔV=I×r
(3)根据始端电压和线路压降计算汇流条II的电压VII
VII=VI-ΔV
通过以上步骤就能得到给定负载水平时,单跟直流线路的压降以及线路末端电压,重复使用该过程即可完成整个飞机直流电网的计算工作。
可以看出,该方法在计算线路电流时首先忽略了线路损耗,并假设了线路电流为定值,因此会一定程度降低计算准确性;同时计算流程复杂,每计算一条线路都需要重复进行以上几个步骤;并且该方法是一种适用于人工计算的简化算法,不利于程序实现,当网络拓扑增大时计算量显著增加。
因此,常规飞机直流电网计算方法存在计算结果准确性低、计算流程复杂、人工计算量大等缺点。而随着计算机技术不断发展,对飞机直流电网进行精确计算称为可能。
发明内容
本发明提供了一种直流电网计算方法及装置,以至少解决存在一定系统误差,并且计算效率低的问题。
根据本发明的一个方面,提供了一种直流电网计算方法,该方法包括:确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压。
优选地,根据如下公式确定所述每个直流汇流条的所述线路总输出功率:其中,Pi是所述线路总输出功率,Vi是第i个直流汇流条的电压值,Vj是与所述第i个直流汇流条相邻的至少一个直流汇流条的电压值,Rij是所述第i个直流汇流条与所述第j个直流汇流条之间的线路电阻值。
优选地,根据所述线路总输出功率、所述负载功率确定所述每个直流汇流条的功率平衡方程为:其中,PLi是接在第i个直流汇流条上的负载功率,i=1,2,……,n,n是系统中直流汇流条的个数。
优选地,在包含n个汇流条的直流系统中,根据所述每个直流汇流条的功率平衡方程得到的所述功率平衡方程组为:
优选地,根据所述功率平衡方程组联合求解所述每个直流汇流条的电压包括:根据牛顿法求解所述功率平衡方程组,得到所述每个直流汇流条的电压。
根据本发明的另一方面,提供了一种直流电网计算装置,包括:线路总功率计算模块,用于确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;功率平衡方程确定模块,用于根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;功率平衡方程组确定模块,用于根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压。
优选地,所述线路总功率计算模块还用于根据如下公式确定所述每个直流汇流条的所述线路总输出功率:其中,Pi是所述线路总输出功率,Vi是第i个直流汇流条的电压值,Vj是与所述第i个直流汇流条相邻的至少一个直流汇流条的电压值,Rij是所述第i个直流汇流条与所述第j个直流汇流条之间的线路电阻值。
优选地,所述功率平衡方程确定模块确定的所述每个直流汇流条的功率平衡方程为:其中,PLi是与第i个直流汇流条相连的负载功率,i=1,2,……,n,其中n是直流汇流条的个数。
优选地,在包含n个汇流条的直流系统中,所述功率平衡方程组确定的所述功率平衡方程组为:
优选地,所述功率平衡方程组确定模块还包括电压确定单元,用于根据牛顿法求解所述
功率平衡方程组,得到所述每个直流汇流条的电压。
通过本发明,采用确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压,解决了相关技术中存在一定系统误差的问题,进而达到了提高计算准确性的效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据相关技术中的直流电路示意图一;
图2是根据本发明实施例中的直流电路示意图一;
图3是根据本发明实施例中的直流电路示意图二;
图4是根据本发明实施例中的直流电路示意图三;
图5是根据本发明实施例的民机直流电网局部示意图;
图6是根据本发明实施例的直流电网计算方法的流程图。
具体实施方式
下文中将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
本发明实施例提出了一种直流电网计算方法,该方法通过联立求解各汇流条功率平衡方程进行电网计算,避免了已有方法存在近似和简化的缺点,提高了计算结果准确性;同时便于程序编制,显著缩短了电网计算时间,避免了人工计算产生误差的可能性。该方法适用于各种电网构型以及不同负载模型下的直流电网计算问题,具有较高的实际应用价值。
图2是根据本发明实施例中的直流电路示意图一,如图2所示,在该直流线路中,已知该导线电阻为r,其两端汇流条编号分别为I和II,其电压幅值分别为VI和VII。
针对图2所示的直流电路,本发明实施例提供的直流电网计算方法可以通过步骤实现。
此时导线流过的电流大小为
则由汇流条I流向导线的功率为
PI=VIII-II (2)
将(1)式代入(2)式展开可得
由电路理论可知:当功率由汇流条流向导线时,(3)式计算结果为正;当功率由导线流向汇流条时,其计算结果为负。
图3是根据本发明实施例中的直流电路示意图二,如图3所示,该直流电路具有接有多根导线的汇流条i。图中与汇流条i相连的汇流条编号为k,l,……,m。并用j∈i表示这些与汇流条i连接的汇流条编号集合。
此时汇流条i向系统注入的功率应等于该汇流条向各条线路注入的功率之和,由(4)式可得该汇流条i向所有线路流入的总功率为
式中:Rij为汇流条i与j间的线路电阻。
将(4)式展开并整理如下
当连在汇流条i上的负载功率为PLi时,根据能量守恒原理,负载功率和流向线路的功率之和应为零,即应满足
此时,(6)式即为汇流条i的功率平衡方程。可以看出,在直流电网计算中每个汇流条都能得到(6)式所示的功率平衡方程,本方法将该方程作为电网计算的基本方程式,通过联立求解该方程组即可得完成电网计算。
图4是根据相关技术中的直流电路示意图三,图中节点0为调压点,电网计算目标是求解除调压点外其余n个汇流条电压值。
针对图4所示的直流电路示意图,本发明实例提供的直流电网计算方法可以通过以下步骤实现。
1)对于节点0,作为调压点该汇流条电压V在运行过程中基本不变,因此应作为已知量在电网计算中预先给定。
2)对于除调压点(节点0)以外的n个汇流条,此时电网计算待求量为这些流条的电压,因此共有n个待求量,将其按汇流条编号顺序排列为
x=[V1,V2,…,Vn]T (7)
由(6)式可知,每个汇流条均可推导得到一个功率平衡方程,在包含n个汇流条的直流系统中,将其按汇流条编号次序排列为
由(7)式可知,此时待求的未知量个数等于n,而(8)式中的方程数也等于n,因此方程数等于待求量数。通过求解(8)式就能得到x向量中的各汇流条电压。
可以看出,此时直流电网计算问题已转化为一个非线性方程组的求解问题,当线路参数和负载大小给定时采用合适方法即可求解。
此时采用了牛顿-拉夫逊法对该方程组进行了求解。该方法是非线性方程计算的经典方法,它将非线性方程组的求解问题转化为反复对相应的线性方程求解的过程,该过程通常称为逐次线性化。将(8)式对各待求量求偏导数得到的线性化方程为
式中:等式左项为功率不平衡量,n阶系数矩阵称为Jacobin矩阵,其计算公式为
根据牛顿-拉夫逊法计算流程可知,此时方程组计算步骤如下:
(1)给定各汇流条电压初值
(2)利用(12)式中的电压初值按(10)式计算功率不平衡量
(3)利用(12)式中的电压初值按(11)式计算Jacobin矩阵各元素,得到(9)式所示的线性化方程。
(4)求解线性修正方程(9)式,得到此时的电压修正量
(5)修正电压初值
V(1)=V(0)-ΔV(0) (15)
(6)将修正后的汇流条电压代入(10)式重新计算功率不平衡量,当功率不平衡量模值小于设定参数时认为计算结束并输出结果;否则令V(0)=V(1)返回第(2)步继续迭代计算。
图5是根据本发明实施例的民机直流电网局部示意图。此时设置节点0为调压点,直流电网额定电压为28V。
针对图5所示的示意图,本发明实施例提供的直流电网计算方法可以通过以下步骤实现。
(1)给定各汇流条电压初值,此时设定为V(0)=[28,28,28,28,28,28,28,28]T。
(2)利用此时电压初值按(10)式计算功率失配量
(3)由(11)式计算Jacobin矩阵,根据失配量ΔP(0)和Jacobin矩阵计算ΔV。
(4)修正电压为V(1)=V(0)-ΔV。
(5)将修正后的电压代入(10)式重新计算功率不平衡量,当不平衡功率小于1.0×10-6时认为迭代结束并输出结果,否则令V(0)=V(1)返回第(2)步重新计算。
此时,给定图5网络中各汇流条负载和各线路参数如表1和表2所示。
表1
表2
当系统参数如表1和表2所示时,采用本实施例的方法电网计算结果如表3所示。
表3
将本方法所得结果代入图5进行验证:以汇流条8为例,结果中该汇流条电压为26.357V,汇流条6电压为27.968V,这两汇流条间电阻为0.02846欧姆。由电路理论可知,此时由线路流向汇流条8的功率为
可以看出,此时由线路流向汇流条8的功率完全等于汇流条8上的负载功率,证明了本方法的准确性。由其它汇流条计算结果也能得到相同结论。
图6是根据本发明实施例的直流电网计算方法的流程图,如图6所示,该方法包括以下步骤:
步骤S602,确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;
步骤S604,根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;
步骤S606,根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压。
本发明实施例还提供了一种直流电网计算装置,包括:
线路总功率计算模块,用于确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;
功率平衡方程确定模块,用于根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;
功率平衡方程组确定模块,用于根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压。
其中,所述线路总功率计算模块还用于根据如下公式确定所述每个直流汇流条的所述线路总输出功率:其中,Pi是所述线路总输出功率,Vi是第i个直流汇流条的电压值,Vj是与所述第i个直流汇流条相邻的至少一个直流汇流条的电压值,Rij是所述第i个直流汇流条与所述第j个直流汇流条之间的线路电阻值。
其中,所述功率平衡方程确定模块确定的所述每个直流汇流条的功率平衡方程为:其中,PLi是与第i个直流汇流条相连的负载功率,i=1,2,……,n,其中n是直流汇流条的个数。
其中,所述功率平衡方程组确定的所述功率平衡方程组为:
其中,所述功率平衡方程组确定模块还包括电压确定单元,用于根据牛顿法求解所述功率平衡方程组,得到所述每个直流汇流条的电压。
本发明实施例具有以下技术效果:本发明提出了一种飞机直流电网计算的新方法,从功率平衡的原理出发得出汇流条功率平衡方程式,并通过联立求解这些方程完成电网计算。由于在建模和计算过程中不存在简化和近似,因此该方法计算结果准确性高。其次,便于程序编制,显著节省了计算时间,同时避免了人工计算存在错误的可能性。第三,该方法灵活性高,可适用于各种构型的电网计算问题,并且可以结合不同的负载模型或者与其它系统交联分析,具有较高的实用价值。
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种直流电网计算方法,其特征在于,包括:
确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;
根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;
根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解得到所述每个直流汇流条的电压。
2.根据权利要求1所述的方法,其特征在于,根据如下公式确定所述每个直流汇流条的所述线路总输出功率:其中,Pi是所述线路总输出功率,Vi是第i个直流汇流条的电压值,Vj是与所述第i个直流汇流条相连的至少一个直流汇流条的电压值,Rij是所述第i个直流汇流条与所述第j个直流汇流条之间的线路电阻。
3.根据权利要求1所述的方法,其特征在于,根据所述线路总输出功率、所述负载功率确定所述每个直流汇流条的功率平衡方程为:其中,PLi是接在第i个直流汇流条上的负载功率,i=1,2,……,n,n是系统中直流汇流条的个数。
4.根据权利要求3所述的方法,其特征在于,根据所述每个直流汇流条的功率平衡方程得到的所述功率平衡方程组:
5.根据权利要求1至4中任一项所述的方法,其特征在于,根据所述功率平衡方程组联合求解所述每个直流汇流条的电压,包括:根据牛顿-拉夫逊法求解所述功率平衡方程组,得到所述每个直流汇流条的电压。
6.一种直流电网计算装置,其特征在于,包括:
线路总功率计算模块,用于确定每个直流汇流条的线路总输出功率,其中,所述线路总输出功率是所述每个直流汇流条的向与其相连的所有线路流出的总功率;
功率平衡方程确定模块,用于根据所述线路总输出功率、负载功率确定所述每个直流汇流条的功率平衡方程,其中,所述负载是所述直流电网中的用电设备;
功率平衡方程组确定模块,用于根据所述每个直流汇流条的功率平衡方程得到功率平衡方程组,并根据所述功率平衡方程组联合求解所述每个直流汇流条的电压。
7.根据权利要求6所述的装置,其特征在于,所述线路总功率计算模块还用于根据如下公式确定所述每个直流汇流条向所述线路总输出功率:其中,Pi是所述线路总输出功率,Vi是第i个直流汇流条的电压值,Vj是与所述第i个直流汇流条相邻的至少一个直流汇流条的电压值,Rij是所述第i个直流汇流条与所述第j个直流汇流条之间的线路电阻值。
8.根据权利要求6所述的装置,其特征在于,所述功率平衡方程确定模块确定的所述每个直流汇流条的功率平衡方程为:其中,PLi是接在第i个直流汇流条上的负载功率,i=1,2,……,n,n是系统中直流汇流条的个数。
9.根据权利要求8所述的装置,其特征在于,在包含n个汇流条的直流系统中,所述功率平衡方程组确定的所述功率平衡方程组为:
10.根据权利要求6至9中任一项所述的装置,其特征在于,所述功率平衡方程组确定模块还包括电压确定单元,用于根据牛顿法求解所述功率平衡方程组,得到所述每个直流汇流条的电压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510489135.3A CN105022931A (zh) | 2015-08-11 | 2015-08-11 | 一种直流电网计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510489135.3A CN105022931A (zh) | 2015-08-11 | 2015-08-11 | 一种直流电网计算方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105022931A true CN105022931A (zh) | 2015-11-04 |
Family
ID=54412897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510489135.3A Pending CN105022931A (zh) | 2015-08-11 | 2015-08-11 | 一种直流电网计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105022931A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1929234A (zh) * | 2006-09-22 | 2007-03-14 | 天津大学 | 大规模电力系统网络分割潮流并行计算方法 |
CN103401234A (zh) * | 2013-07-23 | 2013-11-20 | 东北电力大学 | 一种基于发电机节点类型的潮流计算方法 |
-
2015
- 2015-08-11 CN CN201510489135.3A patent/CN105022931A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1929234A (zh) * | 2006-09-22 | 2007-03-14 | 天津大学 | 大规模电力系统网络分割潮流并行计算方法 |
CN103401234A (zh) * | 2013-07-23 | 2013-11-20 | 东北电力大学 | 一种基于发电机节点类型的潮流计算方法 |
Non-Patent Citations (3)
Title |
---|
姜舒婷等: "含潮流控制器的直流电网潮流计算方法", 《电网技术》 * |
由欣等: "电力系统潮流的有限元分析", 《南昌大学学报(工科版)》 * |
高荣等: "飞机电网仿真技术研究", 《沈阳航空工业学院学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5815382B2 (ja) | 平衡配電系統のための電力潮流解析 | |
CN104049175B (zh) | 具有多个多相母线的电网中的故障分析 | |
CN103825267B (zh) | 一种mmc-mtdc直流侧短路电流的计算方法 | |
CN101635457B (zh) | 一种基于状态估计残差参数灵敏度的电网参数估计方法 | |
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN103887790B (zh) | 输配电网协同潮流控制方法与系统 | |
CN106451456B (zh) | 一种基于环网电力系统的直算方法 | |
Samal et al. | A modified forward backward sweep load flow algorithm for unbalanced radial distribution systems | |
CN106602570B (zh) | 一种基于Matlab的快速分解法潮流计算方法 | |
CN103795057A (zh) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 | |
Semlyen et al. | Quasi-Newton power flow using partial Jacobian updates | |
CN107196306B (zh) | 基于Matlab稀疏矩阵的快速分解法潮流计算方法 | |
CN107069696A (zh) | 一种电力系统状态估计的并行计算方法 | |
CN110518590A (zh) | 考虑负荷电压静特性的配电网线性潮流计算方法 | |
CN105186492B (zh) | 一种基于潮流跟踪和电气剖分方法的网损分摊方法 | |
CN104899396A (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN103956741B (zh) | 一字链及支链式的三相对称多电源非环网电力系统直算法 | |
CN106532712A (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
Zhao et al. | Holomorphic embedding power flow for AC/DC hybrid power systems using Bauer's Eta algorithm | |
Jalili-Marandi et al. | Real-time electromagnetic and transient stability simulations for active distribution networks | |
Sameni et al. | A modified Newton-Raphson method for unbalanced distribution systems | |
Alsaadi et al. | An effective approach for distribution system power flow solution | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN105022931A (zh) | 一种直流电网计算方法及装置 | |
Samuel et al. | INVESTIGATING THE SELECTION OF A SUITABLE SLACK BUS: A CASE STUDY OF THE MULTI-GENERATING STATIONS OF THE NIGERIAN 330-KV POWER SYSTEM NETWORK. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151104 |