CN109818346A - Method and device for transient calculation of closed loop current based on frequency domain analysis - Google Patents

Method and device for transient calculation of closed loop current based on frequency domain analysis Download PDF

Info

Publication number
CN109818346A
CN109818346A CN201711173997.0A CN201711173997A CN109818346A CN 109818346 A CN109818346 A CN 109818346A CN 201711173997 A CN201711173997 A CN 201711173997A CN 109818346 A CN109818346 A CN 109818346A
Authority
CN
China
Prior art keywords
pole
frequency
matrix
public pole
fdne
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
Application number
CN201711173997.0A
Other languages
Chinese (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.)
China Agricultural University
Original Assignee
China Agricultural University
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 China Agricultural University filed Critical China Agricultural University
Priority to CN201711173997.0A priority Critical patent/CN109818346A/en
Publication of CN109818346A publication Critical patent/CN109818346A/en
Pending legal-status Critical Current

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The invention discloses a kind of Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis, wherein method includes: to obtain N2One group of identical initial public pole in a element;Vector fitting is carried out one by one to each element of FDNE matrix according to initial public pole, pole and its number to guarantee each element are identical;If pole is not converged, one group of newly public pole is obtained after vector fitting;Vector fitting is re-started using new public pole as initial public pole, until the difference of initial public pole and new public pole in predetermined threshold level, obtains the FDNE matrix with public pole;Solve cyclization transient current.This method can realize Alloy White Iron Transient calculation by the rapid vector fitting process based on frequency-domain analysis, improve and calculate accuracy and speed, and further increase power supply reliability, ensure the normal production and living of user.

Description

基于频域分析的合环电流暂态计算方法及装置Method and device for transient calculation of closed loop current based on frequency domain analysis

技术领域technical field

本发明涉及电力系统中配电网合环电流暂态计算技术领域,特别涉及一种基于频域分析的合环电流暂态计算方法及装置。The invention relates to the technical field of closed loop current transient calculation in a power system, in particular to a closed loop current transient calculation method and device based on frequency domain analysis.

背景技术Background technique

合环暂态电流是合环能否成功的主要因素,因此,准确快速的计算合环暂态电流关系到能否合环成功,并且对于保障用户正常生产生活,提高供电可靠性具有重要意义。The loop closure transient current is the main factor for the success of the loop closure. Therefore, accurate and fast calculation of the loop closure transient current is related to the success of the loop closure, and is of great significance to ensure the normal production and life of users and improve the reliability of power supply.

相关技术中,有的模型采用简化的戴维南等值电路计算合环电流,但是精算精度不足;通过全网电磁暂态计算合环电流,虽然精度满足,但是因计算量太大无法实现;根据经验判断是否能够合环,只能针对特定的网络,不具有普遍性和代表性,亟待解决。In related technologies, some models use a simplified Thevenin equivalent circuit to calculate the closed loop current, but the actuarial accuracy is insufficient; although the closed loop current is calculated through the electromagnetic transient of the whole network, although the accuracy is satisfactory, it cannot be realized due to the large amount of calculation; according to experience To judge whether the loop can be closed, it can only be aimed at a specific network, which is not universal and representative, and needs to be solved urgently.

发明内容SUMMARY OF THE INVENTION

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。The present invention aims to solve one of the technical problems in the related art at least to a certain extent.

为此,本发明的一个目的在于提出一种基于频域分析的合环电流暂态计算方法,该方法可以提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。Therefore, an object of the present invention is to propose a closed loop current transient calculation method based on frequency domain analysis, which can improve the accuracy and speed of calculation, further improve the reliability of power supply, and ensure the normal production and life of users.

本发明的另一个目的在于提出一种基于频域分析的合环电流暂态计算装置。Another object of the present invention is to provide a closed loop current transient calculation device based on frequency domain analysis.

为达到上述目的,本发明一方面实施例提出了一种基于频域分析的合环电流暂态计算方法,包括以下步骤:获取N2个元素中一组相同的初始公共极点;根据所述初始公共极点对FDNE矩阵的每个元素逐个进行矢量拟合,以保证所述每个元素的极点及其个数相同;如果极点未收敛,则在矢量拟合结束后得到一组新公共极点;将所述新公共极点作为所述初始公共极点重新进行矢量拟合,直至所述初始公共极点与所述新公共极点的差值在预设门槛值内,得到具有公共极点的FDNE矩阵;求解合环暂态电流。In order to achieve the above object, an embodiment of the present invention proposes a closed loop current transient calculation method based on frequency domain analysis, including the following steps: obtaining a set of identical initial common poles in N 2 elements; The common poles perform vector fitting on each element of the FDNE matrix one by one to ensure that the poles and the number of each element are the same; if the poles do not converge, a new set of common poles will be obtained after the vector fitting; The new common pole is used as the initial common pole to perform vector fitting again until the difference between the initial common pole and the new common pole is within a preset threshold value, and an FDNE matrix with a common pole is obtained; transient current.

本发明实施例的基于频域分析的合环电流暂态计算方法,根据本发明实施例提出的基于频域分析的合环电流暂态计算方法,可以通过基于频域分析的快速矢量拟合法实现合环电流暂态计算,克服合环电流暂态计算中精度和速度无法完全满足的缺点,有效提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。The method for transient calculation of closed loop current based on frequency domain analysis according to the embodiment of the present invention and the method for transient calculation of closed loop current based on frequency domain analysis proposed in the embodiment of the present invention can be implemented by a fast vector fitting method based on frequency domain analysis. The closed-loop current transient calculation overcomes the shortcomings that the accuracy and speed cannot be fully satisfied in the closed-loop current transient calculation, effectively improves the calculation accuracy and speed, and further improves the reliability of power supply to ensure the normal production and life of users.

另外,根据本发明上述实施例的基于频域分析的合环电流暂态计算方法还可以具有以下附加的技术特征:In addition, the closed loop current transient calculation method based on frequency domain analysis according to the above-mentioned embodiment of the present invention may also have the following additional technical features:

进一步地,在本发明的一个实施例中,所述获取N2个元素中一组相同的初始公共极点,进一步包括:如果采样点数为N0,则求取N个N×N维的采样值矩阵;将所述N个N×N维的采样值矩阵中每个采样值矩阵中N2个元素值依次相加,以得到N0×1的采样值向量;对所述采样值向量进行矢量拟合,得到有理函数;获取所述有理函数的极点,得到所述初始公共极点。Further, in an embodiment of the present invention, the obtaining a set of identical initial common poles in the N 2 elements further includes: if the number of sampling points is N 0 , obtaining N N×N-dimensional sampling values matrix; add N 2 element values in each sample value matrix in the N N×N-dimensional sample value matrices in turn to obtain N 0 × 1 sample value vectors; perform vectorization on the sample value vectors Fitting to obtain a rational function; obtaining the pole of the rational function to obtain the initial common pole.

进一步地,在本发明的一个实施例中,所述具有公共极点的FDNE矩阵的每个元素都为一个频域函数:Further, in an embodiment of the present invention, each element of the FDNE matrix with a common pole is a frequency domain function:

其中,极点ai和留数c i为实数或为复数共轭对,d和h为实数,n为极点个数,s为复频率。Among them, the pole a i and the residue c i are real numbers or complex conjugate pairs, d and h are real numbers, n is the number of poles, and s is the complex frequency.

进一步地,在本发明的一个实施例中,所述求解合环暂态电流,进一步包括:Further, in an embodiment of the present invention, the solution for the closed loop transient current further includes:

建立FDNE(Frequency Dependent Network Equivalent,基于频域分析技术等值模型,以求解合环暂态电流,其中,输电网络电压和电流关系可以表示为:Establish FDNE (Frequency Dependent Network Equivalent, based on the equivalent model of frequency domain analysis technology to solve the closed loop transient current, where the relationship between the voltage and current of the transmission network can be expressed as:

其中,IB为边界母线电流,IE为除边界母线外其他母线电流,UE为除边界母线外其他母线电压,UB为边界母线电压,f为频率,下标B为边界母线,下标E为输电网络中其他母线;Among them, I B is the current of the boundary bus, I E is the current of the other buses except the boundary bus, U E is the voltage of the other buses except the boundary bus, U B is the voltage of the boundary bus, f is the frequency, the subscript B is the boundary bus, the lower Mark E is other busbars in the transmission network;

以得到收缩到边界后的FDNE矩阵和诺顿等值电流,并且所述诺顿等值电流为:To get the FDNE matrix and Norton equivalent current after shrinking to the boundary, and the Norton equivalent current is:

其中,f0为基频频率。Among them, f 0 is the fundamental frequency.

进一步地,在本发明的一个实施例中,还包括:将所述诺顿等值电流转化为三相电流源;通过隐式梯形法对所述三相电流源进行差分。Further, in an embodiment of the present invention, the method further includes: converting the Norton-equivalent current into a three-phase current source; and performing differential on the three-phase current source through an implicit trapezoidal method.

为达到上述目的,本发明另一方面实施例提出了一种基于频域分析的合环电流暂态计算装置,包括:获取模块,用于获取N2个元素中一组相同的初始公共极点;第一拟合模块,用于根据所述初始公共极点对FDNE矩阵的每个元素逐个进行矢量拟合,以保证所述每个元素的极点及其个数相同;判断模块,用于如果极点未收敛,则在矢量拟合结束后得到一组新公共极点;第二拟合模块,用于将所述新公共极点作为所述初始公共极点重新进行矢量拟合,直至所述初始公共极点与所述新公共极点的差值在预设门槛值内,得到具有公共极点的FDNE矩阵;计算模块,求解合环暂态电流。In order to achieve the above object, another embodiment of the present invention provides a closed loop current transient calculation device based on frequency domain analysis, including: an acquisition module for acquiring a set of identical initial common poles in N 2 elements; The first fitting module is used to perform vector fitting on each element of the FDNE matrix one by one according to the initial common pole, so as to ensure that the pole and the number of each element are the same; the judgment module is used if the pole is not Convergence, then a set of new common poles are obtained after the vector fitting is over; the second fitting module is used to re-fit the vector using the new common poles as the initial common poles until the initial common poles are the same as the all Assuming that the difference between the new common poles is within the preset threshold value, an FDNE matrix with common poles is obtained; the calculation module solves the closed loop transient current.

本发明实施例的基于频域分析的合环电流暂态计算装置,根据本发明实施例提出的基于频域分析的合环电流暂态计算方法,可以通过基于频域分析的快速矢量拟合法实现合环电流暂态计算,克服合环电流暂态计算中精度和速度无法完全满足的缺点,有效提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。The closed loop current transient calculation device based on frequency domain analysis according to the embodiment of the present invention, and the closed loop current transient calculation method based on frequency domain analysis proposed in the embodiment of the present invention, can be realized by the fast vector fitting method based on frequency domain analysis The closed-loop current transient calculation overcomes the shortcomings that the accuracy and speed cannot be fully satisfied in the closed-loop current transient calculation, effectively improves the calculation accuracy and speed, and further improves the reliability of power supply to ensure the normal production and life of users.

另外,根据本发明上述实施例的基于频域分析的合环电流暂态计算装置还可以具有以下附加的技术特征:In addition, the closed loop current transient calculation device based on frequency domain analysis according to the above-mentioned embodiment of the present invention may also have the following additional technical features:

进一步地,在本发明的一个实施例中,所述获取模块,进一步包括:第一获取单元,用于如果采样点数为N0,则求取N个N×N维的采样值矩阵;计算单元,用于将所述N个N×N维的采样值矩阵中每个采样值矩阵中N2个元素值依次相加,以得到N0×1的采样值向量;拟合单元,用于对所述采样值向量进行矢量拟合,得到有理函数;第二获取单元,用于获取所述有理函数的极点,得到所述初始公共极点。Further, in an embodiment of the present invention, the obtaining module further includes: a first obtaining unit, configured to obtain N N×N-dimensional sampling value matrices if the number of sampling points is N 0 ; a calculating unit , which is used to sequentially add N 2 element values in each sampling value matrix in the N N×N-dimensional sampling value matrices to obtain an N 0 × 1 sampling value vector; the fitting unit is used for pairing The sampling value vector is subjected to vector fitting to obtain a rational function; the second obtaining unit is used to obtain the pole of the rational function to obtain the initial common pole.

进一步地,在本发明的一个实施例中,所述具有公共极点的FDNE矩阵的每个元素都为一个频域函数:Further, in an embodiment of the present invention, each element of the FDNE matrix with a common pole is a frequency domain function:

其中,极点ai和留数ci为实数或为复数共轭对,d和h为实数,n为极点个数,s为复频率。Among them, the pole a i and the residue c i are real numbers or complex conjugate pairs, d and h are real numbers, n is the number of poles, and s is the complex frequency.

进一步地,在本发明的一个实施例中,所述计算模块,进一步包括:建立单元,用于建立FDNE等值模型,以求解合环暂态电流,其中,输电网络电压和电流关系可以表示为:Further, in an embodiment of the present invention, the calculation module further includes: a establishing unit for establishing an FDNE equivalent model to solve the closed loop transient current, wherein the relationship between the voltage and current of the transmission network can be expressed as :

其中,IB为边界母线电流,IE为除边界母线外其他母线电流,UE为除边界母线外其他母线电压,UB为边界母线电压,f为频率,下标B为边界母线,下标E为输电网络中其他母线;Among them, I B is the current of the boundary bus, I E is the current of the other buses except the boundary bus, U E is the voltage of the other buses except the boundary bus, U B is the voltage of the boundary bus, f is the frequency, the subscript B is the boundary bus, the lower Mark E is other busbars in the transmission network;

以得到收缩到边界后的FDNE矩阵和诺顿等值电流,并且所述诺顿等值电流为:To get the FDNE matrix and Norton equivalent current after shrinking to the boundary, and the Norton equivalent current is:

其中,f0为基频频率。Among them, f 0 is the fundamental frequency.

进一步地,在本发明的一个实施例中,还包括:将所述诺顿等值电流转化为三相电流源;通过隐式梯形法对所述三相电流源进行差分。Further, in an embodiment of the present invention, the method further includes: converting the Norton-equivalent current into a three-phase current source; and performing differential on the three-phase current source through an implicit trapezoidal method.

本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the present invention will be set forth, in part, from the following description, and in part will be apparent from the following description, or may be learned by practice of the invention.

附图说明Description of drawings

本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of embodiments taken in conjunction with the accompanying drawings, wherein:

图1为根据本发明实施例的基于频域分析的合环电流暂态计算方法的流程图;FIG. 1 is a flowchart of a method for transient calculation of closed loop current based on frequency domain analysis according to an embodiment of the present invention;

图2为根据本发明一个实施例的基于频域分析的合环电流暂态计算方法的流程图;2 is a flowchart of a method for transient calculation of closed loop current based on frequency domain analysis according to an embodiment of the present invention;

图3为根据本发明一个实施例的基于FDNE的等值模型示意图;3 is a schematic diagram of an equivalent model based on FDNE according to an embodiment of the present invention;

图4为根据本发明实施例的基于频域分析的合环电流暂态计算装置的结构示意图。FIG. 4 is a schematic structural diagram of a closed loop current transient calculation device based on frequency domain analysis according to an embodiment of the present invention.

具体实施方式Detailed ways

下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。The following describes in detail the embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein the same or similar reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary, and are intended to explain the present invention and should not be construed as limiting the present invention.

下面参照附图描述根据本发明实施例提出的基于频域分析的合环电流暂态计算方法及装置,首先将参照附图描述根据本发明实施例提出的基于频域分析的合环电流暂态计算方法。The method and device for calculating the closed loop current transient based on frequency domain analysis proposed according to the embodiments of the present invention will be described below with reference to the accompanying drawings. calculation method.

图1是本发明实施例的基于频域分析的合环电流暂态计算方法的流程图。FIG. 1 is a flowchart of a method for transient calculation of closed loop current based on frequency domain analysis according to an embodiment of the present invention.

如图1所示,该基于频域分析的合环电流暂态计算方法包括以下步骤:As shown in Figure 1, the closed loop current transient calculation method based on frequency domain analysis includes the following steps:

在步骤S101中,获取N2个元素中一组相同的初始公共极点。In step S101, a set of identical initial common poles in N 2 elements is obtained.

进一步地,在本发明的一个实施例中,获取N2个元素中一组相同的初始公共极点,进一步包括:如果采样点数为N0,则求取N个N×N维的采样值矩阵;将N个N×N维的采样值矩阵中每个采样值矩阵中N2个元素值依次相加,以得到N0×1的采样值向量;对采样值向量进行矢量拟合,得到有理函数;获取有理函数的极点,得到初始公共极点。Further, in an embodiment of the present invention, obtaining a set of identical initial common poles in N 2 elements further includes: if the number of sampling points is N 0 , obtaining N N×N-dimensional sampling value matrices; Add the N2 element values in each sampled value matrix of N N×N-dimensional sampled value matrices in turn to obtain an N0 ×1 sampled value vector; perform vector fitting on the sampled value vector to obtain a rational function ; Get the poles of the rational function and get the initial common pole.

可以理解的是,结合图1和图2,本发明实施例可以先为这N2个元素求取一组相同的初始公共极点。若采样点数为N0,则可求取N个N×N维的采样值矩阵,并将各个采样值矩阵中N2个元素值依次相加,对于N0个采样值矩阵则可得到一个N0×1的“采样值向量”,将这个“采样值向量”进行矢量拟合,可以得到一个有理函数,将这个有理函数的极点设为原始FDNE矩阵的初始公共极点。It can be understood that, with reference to FIG. 1 and FIG. 2 , in this embodiment of the present invention, a set of the same initial common poles may be obtained for the N 2 elements first. If the number of sampling points is N 0 , N N×N-dimensional sampling value matrices can be obtained, and the N 2 element values in each sampling value matrix can be added in turn. For N 0 sampling value matrices, an N A "sampled value vector" of 0 × 1, and by performing vector fitting on this "sampled value vector", a rational function can be obtained, and the pole of this rational function is set as the initial common pole of the original FDNE matrix.

在步骤S102中,根据初始公共极点对FDNE矩阵的每个元素逐个进行矢量拟合,以保证每个元素的极点及其个数相同。In step S102, vector fitting is performed on each element of the FDNE matrix one by one according to the initial common pole, so as to ensure that the poles and the number of each element are the same.

可以理解的是,结合图1和图2,基于步骤S101中的初始公共极点,对这个FDNE矩阵中的各个元素逐个进行矢量拟合,以保证FDNE矩阵各个元素的极点及其个数相同。It can be understood that, with reference to FIG. 1 and FIG. 2, based on the initial common pole in step S101, vector fitting is performed on each element in the FDNE matrix one by one to ensure that the poles and the number of each element of the FDNE matrix are the same.

在步骤S103中,如果极点未收敛,则在矢量拟合结束后得到一组新公共极点。In step S103, if the poles do not converge, a new set of common poles are obtained after the vector fitting is completed.

在步骤S104中,将新公共极点作为初始公共极点重新进行矢量拟合,直至初始公共极点与新公共极点的差值在预设门槛值内,得到具有公共极点的FDNE矩阵。In step S104, the new common pole is used as the initial common pole to perform vector fitting again until the difference between the initial common pole and the new common pole is within a preset threshold value, and an FDNE matrix with common poles is obtained.

可以理解的是,结合图1和图2,由矢量拟合法可知,当极点未收敛时,矢量拟合法结束后可以求出一组新极点,并且可以将这组新极点重新作为步骤S102中的初始公共极点,并重新进行矢量拟合,直到初始公共节点与算法结束后的新极点的差值在一定门槛值内,一般情况下迭代3-4次即可收敛并求得具有公共极点的FDNE矩阵。It can be understood that, with reference to Fig. 1 and Fig. 2, it can be seen from the vector fitting method that when the poles are not converged, a new set of poles can be obtained after the vector fitting method is completed, and the new set of poles can be re-used as the parameters in step S102. Initial common pole, and re-fit the vector until the difference between the initial common node and the new pole after the algorithm ends is within a certain threshold. Generally, iterate 3-4 times to converge and obtain the FDNE with the common pole. matrix.

可以理解的是,本发明实施例在矢量拟合的过程中,可以采取如下实用化处理方法:It can be understood that, in the process of vector fitting in this embodiment of the present invention, the following practical processing methods can be adopted:

首先在实际求取过程中,一般可以令h为0,对于实际网络而言,在无限大的频率下,实际网络的节点导纳不可能为无穷大。First, in the actual calculation process, h can generally be set to 0. For the actual network, at an infinite frequency, the node admittance of the actual network cannot be infinite.

因此,为了保证拟合的有理函数在基频和宽频下都呈现精确的频率特性,一般在拟合时对基频和宽频下的采样值设置相对较大的权重。Therefore, in order to ensure that the fitted rational function exhibits accurate frequency characteristics at both the fundamental frequency and the broadband frequency, relatively large weights are generally set to the sampled values at the fundamental frequency and the broadband frequency during fitting.

快速矢量拟合方法可以保证Y(s)各元素的极点是相同的,因此只有留数项和常数项不相同,则N×N维FDNE矩阵Y(s)可表示为:The fast vector fitting method can ensure that the poles of each element of Y(s) are the same, so only the residue term and the constant term are different, then the N×N-dimensional FDNE matrix Y(s) can be expressed as:

将上述矩阵写成传递函数的形式:Write the above matrix in the form of a transfer function:

Y(s)=C(sE-A)-1B+D,Y(s)=C(sE-A) -1 B+D,

式中为E单位矩阵,其维数与矩阵A相同。where E is the identity matrix, and its dimension is the same as that of matrix A.

A=diag(A1 … Ak …AN),A=diag(A 1 . . . A k . . A N ),

B=diag(B1 … Bk … BN),B=diag(B 1 ... B k ... B N ),

Bk=[(1 1 … 1)(1×n)]TB k = [(1 1 … 1) (1×n) ] T ,

k,m=1,2,…,N。k,m=1,2,...,N.

在步骤S105中,求解合环暂态电流。In step S105, the closed loop transient current is solved.

进一步地,在本发明的一个实施例中,求解合环暂态电流,进一步包括:建立FDNE等值模型,以求解合环暂态电流,其中,输电网络电压和电流关系可以表示为:Further, in an embodiment of the present invention, solving the closed loop transient current further includes: establishing an FDNE equivalent model to solve the closed loop transient current, wherein the relationship between the voltage and current of the transmission network can be expressed as:

其中,IB为边界母线电流,IE为除边界母线外其他母线电流,UE为除边界母线外其他母线电压,UB为边界母线电压,f为频率,下标B为边界母线,下标E为输电网络中其他母线;Among them, I B is the current of the boundary bus, I E is the current of the other buses except the boundary bus, U E is the voltage of the other buses except the boundary bus, U B is the voltage of the boundary bus, f is the frequency, the subscript B is the boundary bus, the lower Mark E is other busbars in the transmission network;

以得到收缩到边界后的FDNE矩阵和诺顿等值电流,并且诺顿等值电流为:To get the FDNE matrix and Norton equivalent current after shrinking to the boundary, and the Norton equivalent current is:

其中,f0为基频频率。Among them, f 0 is the fundamental frequency.

可以理解的是,FDNE等值模型如图3所示,本发明实施例可以运用快速矢量拟合法求解合环暂态电流。本发明实施例考虑到FDNE,因此输电网络电压和电流关系可以表示为:It can be understood that the FDNE equivalent model is shown in FIG. 3 , and the embodiment of the present invention can use the fast vector fitting method to solve the closed loop transient current. In the embodiment of the present invention, FDNE is considered, so the relationship between the voltage and current of the transmission network can be expressed as:

其中,下标B表示边界母线,而下标E表示输电网络中其他母线。对公式进行高斯消去,消去UE,从而得到收缩到边界后的FDNE矩阵和诺顿等值电流。由于220kV电压等级的输电网离合环点距离较远,因此诺顿等值电流可只考虑基频分量,故诺顿等值公式为:Among them, the subscript B refers to the boundary bus, and the subscript E refers to the other bus in the transmission network. Gaussian elimination is performed on the formula and UE is eliminated to obtain the FDNE matrix and Norton's equivalent current after shrinking to the boundary. Due to the long distance between the clutch and loop points of the transmission network with the voltage level of 220kV, the Norton equivalent current can only consider the fundamental frequency component, so the Norton equivalent formula is:

其中,f0是基频频率。where f 0 is the fundamental frequency.

进一步地,在本发明的一个实施例中,本发明实施例的方法还包括:将诺顿等值电流转化为三相电流源;通过隐式梯形法对三相电流源进行差分。Further, in an embodiment of the present invention, the method of the embodiment of the present invention further includes: converting the Norton equivalent current into a three-phase current source; and performing differential on the three-phase current source through an implicit trapezoidal method.

需要说明的是,为了使上述诺顿等值结果可以应用合环电流暂态计算中,需要做两步处理:首先将诺顿等值电流转化为三相电流源;其次使用隐式梯形法进行差分。由于本发明实施例的FDNE矩阵中各个元素的极点及其个数都相同,因此可以统一差分。并且根据仿真结果,计算合环暂态电流峰值与暂态特性。It should be noted that, in order to make the above Norton equivalent results applicable to the closed loop current transient calculation, two steps are required: first, the Norton equivalent current is converted into a three-phase current source; second, the implicit trapezoidal method is used to differentiate. Since the poles and the number of each element in the FDNE matrix in the embodiment of the present invention are the same, the difference can be unified. And according to the simulation results, the closed loop transient current peak value and transient characteristics are calculated.

进一步地,在本发明的一个实施例中,具有公共极点的FDNE矩阵的每个元素都为一个频域函数:Further, in an embodiment of the present invention, each element of the FDNE matrix with a common pole is a frequency domain function:

其中,极点ai和留数ci为实数或为复数共轭对,d和h为实数,n为极点个数,s为复频率。Among them, the pole a i and the residue c i are real numbers or complex conjugate pairs, d and h are real numbers, n is the number of poles, and s is the complex frequency.

可以理解的是,FDNE的实质是一个以频率为函数的节点导纳矩阵,它可以表示频率变化情况下的节点导纳矩阵。因此,N×N维FDNE可以表示为一个频域矩阵:It can be understood that the essence of FDNE is a nodal admittance matrix which is a function of frequency, which can represent the nodal admittance matrix in the case of frequency variation. Therefore, an N×N-dimensional FDNE can be expressed as a frequency-domain matrix:

矢量拟合法将FDNE矩阵的每个元素在一系列频率下的采样值拟合成一个连续的有理函数,并通过矢量拟合法求出的FDNE矩阵中的每个元素都可以表示为一个频域函数:The vector fitting method fits the sampled values of each element of the FDNE matrix at a series of frequencies into a continuous rational function, and each element in the FDNE matrix obtained by the vector fitting method can be expressed as a frequency domain function :

其中,极点ai和留数ci是实数或者分别以复数共轭对出现;d和h为实数;n为极点个数,FDNE矩阵中不同元素的ai、ci、d和h是不相同的。Among them, the pole a i and the residue c i are real numbers or appear as complex conjugate pairs respectively; d and h are real numbers; n is the number of poles, and a i , c i , d and h of different elements in the FDNE matrix are different identical.

可以理解的是,当N2个有理函数的极点不相同时,把FDNE矩阵中的N2个元素逐个应用到时域仿真,其计算量会急剧增加,同时时域实现也不方便,快速矢量拟合法可以有效地降低计算量。It can be understood that when the poles of the N 2 rational functions are different, applying the N 2 elements in the FDNE matrix to the time-domain simulation one by one, the calculation amount will increase sharply, and the time-domain implementation is inconvenient, and the fast vector The fitting method can effectively reduce the amount of calculation.

根据本发明实施例提出的基于频域分析的合环电流暂态计算方法,可以通过基于频域分析的快速矢量拟合法实现合环电流暂态计算,克服合环电流暂态计算中精度和速度无法完全满足的缺点,有效提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。According to the method for transient calculation of closed loop current based on frequency domain analysis proposed in the embodiment of the present invention, the transient calculation of closed loop current can be realized through the fast vector fitting method based on frequency domain analysis, which overcomes the accuracy and speed of transient calculation of closed loop current. The shortcomings that cannot be fully satisfied can effectively improve the accuracy and speed of calculation, and further improve the reliability of power supply to ensure the normal production and life of users.

其次参照附图描述根据本发明实施例提出的基于频域分析的合环电流暂态计算装置。Next, a closed loop current transient calculation device based on frequency domain analysis proposed according to an embodiment of the present invention will be described with reference to the accompanying drawings.

图4是本发明实施例的基于频域分析的合环电流暂态计算装置的结构示意图。FIG. 4 is a schematic structural diagram of a closed loop current transient calculation device based on frequency domain analysis according to an embodiment of the present invention.

如图4所示,该基于频域分析的合环电流暂态计算装置10包括:获取模块100、第一拟合模块200、判断模块300、第二拟合模块400和计算模块500。As shown in FIG. 4 , the closed loop current transient calculation device 10 based on frequency domain analysis includes: an acquisition module 100 , a first fitting module 200 , a judgment module 300 , a second fitting module 400 and a calculation module 500 .

其中,获取模100用于获取N2个元素中一组相同的初始公共极点。第一拟合模块200用于根据初始公共极点对FDNE矩阵的每个元素逐个进行矢量拟合,以保证每个元素的极点及其个数相同。判断模块300用于如果极点未收敛,则在矢量拟合结束后得到一组新公共极点。第二拟合模块400用于将新公共极点作为初始公共极点重新进行矢量拟合,直至初始公共极点与新公共极点的差值在预设门槛值内,得到具有公共极点的FDNE矩阵。计算模块500求解合环暂态电流。本发明实施例的装置10可以通过基于频域分析的快速矢量拟合法实现合环电流暂态计算,提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。Wherein, the acquisition modulo 100 is used to acquire a set of identical initial common poles in the N 2 elements. The first fitting module 200 is configured to perform vector fitting on each element of the FDNE matrix one by one according to the initial common pole, so as to ensure that the poles and the number of each element are the same. The judging module 300 is used to obtain a new set of common poles after the vector fitting ends if the poles are not converged. The second fitting module 400 is configured to perform vector fitting again using the new common pole as the initial common pole, until the difference between the initial common pole and the new common pole is within a preset threshold, and obtain an FDNE matrix with common poles. The calculation module 500 solves the closed loop transient current. The device 10 in the embodiment of the present invention can realize the closed loop current transient calculation through the fast vector fitting method based on frequency domain analysis, improve the calculation accuracy and speed, further improve the reliability of power supply, and ensure the normal production and life of users.

进一步地,在本发明的一个实施例中,获取模块100进一步包括:第一获取单元、计算单元、拟合单元和第二获取单元。Further, in an embodiment of the present invention, the obtaining module 100 further includes: a first obtaining unit, a calculating unit, a fitting unit and a second obtaining unit.

其中,第一获取单元用于如果采样点数为N0,则求取N个N×N维的采样值矩阵。计算单元用于将N个N×N维的采样值矩阵中每个采样值矩阵中N2个元素值依次相加,以得到N0×1的采样值向量。拟合单元用于对采样值向量进行矢量拟合,得到有理函数。第二获取单元用于获取有理函数的极点,得到初始公共极点。Wherein, the first obtaining unit is configured to obtain N N×N-dimensional sampling value matrices if the number of sampling points is N 0 . The computing unit is configured to sequentially add N 2 element values in each sample value matrix in the N N×N-dimensional sample value matrices, so as to obtain an N 0 ×1 sample value vector. The fitting unit is used to perform vector fitting on a vector of sampled values to obtain a rational function. The second obtaining unit is used to obtain the poles of the rational function to obtain the initial common pole.

进一步地,在本发明的一个实施例中,具有公共极点的FDNE矩阵的每个元素都为一个频域函数:Further, in an embodiment of the present invention, each element of the FDNE matrix with a common pole is a frequency domain function:

其中,极点ai和留数ci为实数或为复数共轭对,d和h为实数,n为极点个数,s为复频率。Among them, the pole a i and the residue c i are real numbers or complex conjugate pairs, d and h are real numbers, n is the number of poles, and s is the complex frequency.

进一步地,在本发明的一个实施例中,计算模块500进一步包括:建立单元。其中,建立单元用于建立FDNE等值模型,以求解合环暂态电流,其中,输电网络电压和电流关系可以表示为:Further, in an embodiment of the present invention, the computing module 500 further includes: a establishing unit. Among them, the establishment unit is used to establish the FDNE equivalent model to solve the closed loop transient current, where the relationship between the voltage and current of the transmission network can be expressed as:

其中,IB为边界母线电流,IE为除边界母线外其他母线电流,UE为除边界母线外其他母线电压,UB为边界母线电压,f为频率,下标B为边界母线,下标E为输电网络中其他母线;Among them, I B is the current of the boundary bus, I E is the current of the other buses except the boundary bus, U E is the voltage of the other buses except the boundary bus, U B is the voltage of the boundary bus, f is the frequency, the subscript B is the boundary bus, the lower Mark E is other busbars in the transmission network;

以得到收缩到边界后的FDNE矩阵和诺顿等值电流,并且诺顿等值电流为:To get the FDNE matrix and Norton equivalent current after shrinking to the boundary, and the Norton equivalent current is:

其中,f0为基频频率。Among them, f 0 is the fundamental frequency.

进一步地,在本发明的一个实施例中,本发明实施例的装置10还包括:将诺顿等值电流转化为三相电流源;通过隐式梯形法对三相电流源进行差分。Further, in an embodiment of the present invention, the apparatus 10 of the embodiment of the present invention further includes: converting the Norton equivalent current into a three-phase current source; and performing differential on the three-phase current source through an implicit trapezoidal method.

需要说明的是,前述对基于频域分析的合环电流暂态计算方法实施例的解释说明也适用于该实施例的基于频域分析的合环电流暂态计算装置,此处不再赘述。It should be noted that the foregoing explanation of the embodiment of the closed loop current transient calculation method based on frequency domain analysis is also applicable to the closed loop current transient calculation device based on frequency domain analysis in this embodiment, and will not be repeated here.

根据本发明实施例提出的基于频域分析的合环电流暂态计算装置,可以通过基于频域分析的快速矢量拟合法实现合环电流暂态计算,克服合环电流暂态计算中精度和速度无法完全满足的缺点,有效提高计算准确性和速度性,并且进一步提高供电可靠性,保障用户的正常生产生活。According to the closed loop current transient calculation device based on the frequency domain analysis proposed in the embodiment of the present invention, the closed loop current transient calculation can be realized through the fast vector fitting method based on the frequency domain analysis, which overcomes the accuracy and speed of the closed loop current transient calculation. The shortcomings that cannot be fully satisfied can effectively improve the accuracy and speed of calculation, and further improve the reliability of power supply to ensure the normal production and life of users.

在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", " Rear, Left, Right, Vertical, Horizontal, Top, Bottom, Inner, Outer, Clockwise, Counterclockwise, Axial, The orientations or positional relationships indicated by "radial direction", "circumferential direction", etc. are based on the orientations or positional relationships shown in the accompanying drawings, which are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying the indicated devices or elements. It must have a specific orientation, be constructed and operate in a specific orientation, and therefore should not be construed as a limitation of the present invention.

此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。In addition, the terms "first" and "second" are only used for descriptive purposes, and should not be construed as indicating or implying relative importance or implying the number of indicated technical features. Thus, a feature delimited with "first", "second" may expressly or implicitly include at least one of that feature. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise expressly and specifically defined.

在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。In the present invention, unless otherwise expressly specified and limited, the terms "installed", "connected", "connected", "fixed" and other terms should be understood in a broad sense, for example, it may be a fixed connection or a detachable connection , or integrated; it can be a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between the two elements, unless otherwise specified limit. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。In the present invention, unless otherwise expressly specified and limited, a first feature "on" or "under" a second feature may be in direct contact between the first and second features, or the first and second features indirectly through an intermediary touch. Also, the first feature being "above", "over" and "above" the second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is level higher than the second feature. The first feature being "below", "below" and "below" the second feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature has a lower level than the second feature.

在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。In the description of this specification, description with reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples", etc., mean specific features described in connection with the embodiment or example , structure, material or feature is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and combine the different embodiments or examples described in this specification, as well as the features of the different embodiments or examples, without conflicting each other.

尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。Although the embodiments of the present invention have been shown and described above, it should be understood that the above-mentioned embodiments are exemplary and should not be construed as limiting the present invention. Embodiments are subject to variations, modifications, substitutions and variations.

Claims (10)

1. a kind of Alloy White Iron Method for Transient Calculation based on frequency-domain analysis, which comprises the following steps:
Obtain N2One group of identical initial public pole in a element;
Vector fitting is carried out to each element of FDNE matrix according to the initial public pole, one by one to guarantee each member The pole and its number of element are identical;
If pole is not converged, one group of newly public pole is obtained after vector fitting;
Vector fitting is re-started using the new public pole as the initial public pole, until the initial public pole Difference with the new public pole obtains the FDNE matrix with public pole in predetermined threshold level;And
Solve cyclization transient current.
2. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that described to obtain Take N2One group of identical initial public pole in a element further comprises:
If sampling number is N0, then the sampling value matrix of N number of N × N-dimensional is sought;
By N in sampling value matrix each in the sampling value matrix of N number of N × N-dimensional2A element value is successively added, to obtain N0× 1 sampled value vector;
Vector fitting is carried out to the sampled value vector, obtains rational function;
The pole for obtaining the rational function obtains the initial public pole.
3. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that the tool The each element for having the FDNE matrix of public pole is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is complex frequency.
4. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that described to ask Cyclization transient current is solved, further comprises:
FDNE Equivalent Model is established, to solve cyclization transient current, wherein power transmission network voltage and current relationship can indicate Are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor other mothers in addition to the bus of boundary Line voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is other buses in power transmission network;
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
5. the Alloy White Iron Method for Transient Calculation according to claim 4 based on frequency-domain analysis, which is characterized in that also wrap It includes:
Three-phase current source is converted by the promise equal currents;
Difference is carried out to the three-phase current source by implicit trapezoid method.
6. a kind of Alloy White Iron Transient calculation device based on frequency-domain analysis characterized by comprising
Module is obtained, for obtaining N2One group of identical initial public pole in a element;
First fitting module, it is quasi- for carrying out vector one by one to each element of FDNE matrix according to the initial public pole It closes, pole and its number to guarantee each element are identical;
If judgment module obtains one group of newly public pole not converged for pole after vector fitting;
Second fitting module, for re-starting vector fitting for the new public pole as the initial public pole, directly To the difference of the initial public pole and the new public pole in predetermined threshold level, the FDNE with public pole is obtained Matrix;And
Computing module solves cyclization transient current.
7. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that described to obtain Modulus block further comprises:
First acquisition unit, if being N for sampling number0, then the sampling value matrix of N number of N × N-dimensional is sought;
Computing unit is used for N in sampling value matrix each in the sampling value matrix of N number of N × N-dimensional2A element value successively phase Add, to obtain N0× 1 sampled value vector;
Fitting unit obtains rational function for carrying out vector fitting to the sampled value vector;
Second acquisition unit obtains the initial public pole for obtaining the pole of the rational function.
8. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that the tool The each element for having the FDNE matrix of public pole is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is complex frequency.
9. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that the meter Module is calculated, further comprises:
Unit is established, for establishing FDNE Equivalent Model, to solve cyclization transient current, wherein power transmission network voltage and current Relationship can indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor other mothers in addition to the bus of boundary Line voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is other buses in power transmission network;
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
10. the Alloy White Iron Method for Transient Calculation according to claim 9 based on frequency-domain analysis, which is characterized in that also wrap It includes:
Three-phase current source is converted by the promise equal currents;
Difference is carried out to the three-phase current source by implicit trapezoid method.
CN201711173997.0A 2017-11-22 2017-11-22 Method and device for transient calculation of closed loop current based on frequency domain analysis Pending CN109818346A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711173997.0A CN109818346A (en) 2017-11-22 2017-11-22 Method and device for transient calculation of closed loop current based on frequency domain analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711173997.0A CN109818346A (en) 2017-11-22 2017-11-22 Method and device for transient calculation of closed loop current based on frequency domain analysis

Publications (1)

Publication Number Publication Date
CN109818346A true CN109818346A (en) 2019-05-28

Family

ID=66601412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711173997.0A Pending CN109818346A (en) 2017-11-22 2017-11-22 Method and device for transient calculation of closed loop current based on frequency domain analysis

Country Status (1)

Country Link
CN (1) CN109818346A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113051777A (en) * 2021-04-28 2021-06-29 北京华大九天科技股份有限公司 Method for correcting data by using vector fitting

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103106328A (en) * 2012-11-02 2013-05-15 南方电网科学研究院有限责任公司 Frequency correlation network equivalence generation method based on integral vector fitting method
CN106528934A (en) * 2016-10-10 2017-03-22 清华大学 Frequency-related network equivalence passivity correction method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103106328A (en) * 2012-11-02 2013-05-15 南方电网科学研究院有限责任公司 Frequency correlation network equivalence generation method based on integral vector fitting method
CN106528934A (en) * 2016-10-10 2017-03-22 清华大学 Frequency-related network equivalence passivity correction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XU ZHANG等: ""Implicitly Coupled Electromechanical and Electromagnetic Transient Analysis Using a Frequency-Dependent Network Equivalent"", 《IEEE TRANSACTIONS ON POWER DELIVERY》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113051777A (en) * 2021-04-28 2021-06-29 北京华大九天科技股份有限公司 Method for correcting data by using vector fitting
CN113051777B (en) * 2021-04-28 2021-08-31 北京华大九天科技股份有限公司 Method for correcting data by using vector fitting

Similar Documents

Publication Publication Date Title
CN108667048B (en) Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system
WO2021082040A1 (en) Constraint aggregation-based virtual power plant feasible output range identification method and apparatus
Zhou et al. Measurement-based harmonic modeling of an electric vehicle charging station using a three-phase uncontrolled rectifier
CN106777827B (en) Electromechanical-electromagnetic hybrid simulation method and system
CN109786897B (en) Lithium ion battery alternating excitation low-temperature heating method based on temperature change
CN108879664B (en) AC/DC system online voltage stability evaluation method based on wide area measurement
CN109900937A (en) A state-of-charge estimation method for lithium batteries with temperature compensation
CN103904643B (en) A kind of DC power flow computational methods considering network loss
CN109802406A (en) A method of analysis flexible DC transmission access system resonance stability
CN105701568B (en) A Heuristic Method for Rapid Optimization of Measuring Locations for Distribution Network State Estimation
CN107342586A (en) A kind of method and device for being used to determine power distribution network Static Equivalent
CN114935692A (en) A kind of converter impedance measurement method and device
CN115421044A (en) Lithium-ion battery state determination method, device, electronic equipment and storage medium
CN107947199B (en) Optimization model-based method for searching thermal stability security domain boundary of electric power system
CN109657276A (en) Electromagnetical transient emulation method and system
CN114626959B (en) A robust state estimation method for distribution network based on second-order cone programming considering multi-source measurements
CN109818346A (en) Method and device for transient calculation of closed loop current based on frequency domain analysis
CN105305392B (en) Short circuit calculation symmetrical component method suitable for the type IIDG power distribution network of control containing voltage
CN104917197B (en) A kind of method of parallel computation active distribution network three-phase unbalanced load flow
CN102709911B (en) Method for designing interface with harmonic characteristic hybrid simulation function
Gurusinghe et al. Efficient algorithms for real‐time monitoring of transmission line parameters and their performance with practical synchrophasors
CN107017631A (en) A kind of three-phase distribution net Observability analysis of power system based on linear circuit
CN104821577A (en) Three-phase four-wire system distribution network robust estimation method based on intelligent electric meter measurement
CN104062501B (en) Double-transformer substation harmonic wave state estimation method
Li et al. Modelling of the wide frequency equivalent circuit of the three-phase AC motor based on vector fitting method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190528

RJ01 Rejection of invention patent application after publication