CN104598728A - Wind power generation-including power system state estimation method taking frequency change into consideration - Google Patents
Wind power generation-including power system state estimation method taking frequency change into consideration Download PDFInfo
- Publication number
- CN104598728A CN104598728A CN201510010384.XA CN201510010384A CN104598728A CN 104598728 A CN104598728 A CN 104598728A CN 201510010384 A CN201510010384 A CN 201510010384A CN 104598728 A CN104598728 A CN 104598728A
- Authority
- CN
- China
- Prior art keywords
- partiald
- state estimation
- power
- state
- model
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 6
- 206010008190 Cerebrovascular accident Diseases 0.000 claims 1
- 208000006011 Stroke Diseases 0.000 claims 1
- 239000003990 capacitor Substances 0.000 claims 1
- 230000005684 electric field Effects 0.000 claims 1
- 238000010248 power generation Methods 0.000 abstract description 9
- 230000001360 synchronised effect Effects 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000002347 injection Methods 0.000 abstract description 3
- 239000007924 injection Substances 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种计及频率变化的含风力发电的电力系统状态估计方法,随着大规模风电接入电网,风电出力的波动性会对电网频率产生显著影响,因此传统的状态估计模型已不再适用。本发明在同步发电机准稳态模型、负荷准稳态模型和异步风机简化RX模型的基础上,将频率偏差作为新的状态量引入到状态估计过程中,并且对发电机和负荷节点构造了新的零注入功率量测,最终建立了考虑频率偏差的状态估计模型。算例仿真结果表明该模型可有效计及系统的频率偏差,状态估计结果精度得到明显提高。
The invention discloses a method for state estimation of power systems including wind power generation that takes into account frequency changes. With large-scale wind power connected to the power grid, the fluctuation of wind power output will have a significant impact on the frequency of the power grid. Therefore, the traditional state estimation model has been No longer applicable. Based on the quasi-steady-state model of the synchronous generator, the quasi-steady-state model of the load and the simplified RX model of the asynchronous fan, the present invention introduces the frequency deviation as a new state quantity into the state estimation process, and constructs a model for the generator and load nodes The new zero-injection power measurement finally establishes a state estimation model considering frequency deviation. The simulation results of the example show that the model can effectively take into account the frequency deviation of the system, and the accuracy of the state estimation results is significantly improved.
Description
技术领域technical field
发明涉及一种计及频率变化的含风力发电的电力系统状态估计方法,属于电力系统运行和控制技术领域。The invention relates to a method for estimating the state of a power system including wind power generation in consideration of frequency changes, and belongs to the technical field of power system operation and control.
背景技术Background technique
作为能量管理系统(Energy Management System,EMS)的核心,电力系统状态估计通过对生数据的处理,获得状态量的最佳估计值。传统的加权最小二乘法(Weighted LeastSquares,WLS)状态估计算法估计质量和收敛性能很好,是状态估计的经典解法和理论基础,适应各种类型的量测系统。As the core of the Energy Management System (EMS), the state estimation of the power system obtains the best estimated value of the state quantity through the processing of raw data. The traditional Weighted Least Squares (WLS) state estimation algorithm has good estimation quality and convergence performance. It is a classic solution and theoretical basis for state estimation and is suitable for various types of measurement systems.
异步风机简化RX模型既充分考虑了异步发电机本身的特性,较详细地阐述了异步风力发电机的输出功率特性,又较传统RX模型的计算量小,精度满足计算要求。The simplified RX model of asynchronous wind turbine not only fully considers the characteristics of the asynchronous generator itself, but also elaborates the output power characteristics of the asynchronous wind turbine in more detail, and is less computationally intensive than the traditional RX model, and the accuracy meets the calculation requirements.
但现有技术中的状态估计模型仅仅考虑了风机的滑差,并未研究系统的频率变化。随着大规模风电接入电网,风电出力的波动性会对电网频率产生显著影响。However, the state estimation model in the prior art only considers the slip of the fan, and does not study the frequency change of the system. With the integration of large-scale wind power into the grid, the fluctuation of wind power output will have a significant impact on the frequency of the grid.
发明内容Contents of the invention
本发明提供一种计及频率变化的含风力发电的电力系统状态估计方法,将系统的频率偏差作为新的状态量建立新的状态估计模型,相关电力元件也采用考虑频率特性的模型,该方法能够有效计及系统的频率偏差,解决了现有技术中的问题,具有工程应用价值。The present invention provides a method for estimating the state of a power system including wind power generation that takes into account frequency changes. The frequency deviation of the system is used as a new state quantity to establish a new state estimation model, and the relevant power components also use a model that considers frequency characteristics. The method The frequency deviation of the system can be effectively taken into account, the problems in the prior art are solved, and the invention has engineering application value.
本发明为实现上述目的,采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:
一种计及频率变化的含风力发电的电力系统状态估计方法,依次包括以下步骤:A method for estimating the state of a power system including wind power generation in consideration of frequency changes, comprising the following steps in sequence:
(1)获得电力系统的网络参数和量测量;(1) Obtain network parameters and quantity measurements of the power system;
(2)利用步骤(1)获得的网络参数和量测量进行状态估计程序初始化;(2) utilize the network parameter and quantity measurement that step (1) obtains to carry out state estimation program initialization;
(3)建立计及频率变化的含风力发电的状态估计模型:(3) Establish a state estimation model including wind power generation that takes into account frequency changes:
min J(x)=[z-h(x)]TW[z-h(x)]min J(x)=[zh(x)] T W[zh(x)]
其中,J是目标函数;T表示矩阵的转置;W为对角权重矩阵;x为状态量,维数n=2N-1,N为节点数;z为量测量,维数m,m>n;h为m维非线性量测函数;Among them, J is the objective function; T represents the transposition of the matrix; W is the diagonal weight matrix; x is the state quantity, the dimension n=2N-1, N is the number of nodes; z is the quantity measurement, the dimension m, m> n; h is an m-dimensional nonlinear measurement function;
(4)考虑风机滑差s和系统频率偏差Δf建立新的雅克比矩阵:(4) Consider the fan slip s and the system frequency deviation Δf to establish a new Jacobian matrix:
其中,P和Q分别表示普通发电机对应的有功和无功;Pk和Qk分别为异步风机的接入节点k对应的有功和无功;其中和的维数与系统中风电场节点的数目相同;Among them, P and Q respectively represent the active power and reactive power corresponding to the ordinary generator; P k and Q k are respectively the active power and reactive power corresponding to the access node k of the asynchronous fan; where and The dimension of is the same as the number of wind farm nodes in the system;
(5)由初始各状态量V(0)、θ(0)、s(0)、Δf(0)计算量测量的计算值h(x(k))和雅克比矩阵H(x(k)),其中,θ表示电压相角,V表示电压幅值,s是滑差,Δf是频率偏差量,0表示初始状态量;(5) Calculated value h(x ( k) ) and Jacobian matrix H ( x ( k ) ), where θ represents the voltage phase angle, V represents the voltage amplitude, s is the slip, Δf is the frequency deviation, and 0 represents the initial state quantity;
(6)求解状态修正量Δx(k),判断是否满足收敛条件,若max{ΔV(k)|,|Δθ(k)|,|Δs(k)|,|Δf(k)|}>λ,迭代次数k=k+1,修正状态量V(k+1)=V(k)+ΔV(k),θ(k+1)=θ(k)+Δθ(k),s(k+1)=s(k)+Δs(k),Δf(k+1)=Δf(k)+ΔΔf(k),直到符合条件,输出结果,其中,k是迭代次数。(6) Solve the state correction amount Δx (k) to judge whether the convergence condition is satisfied, if max{ΔV (k) |,|Δθ (k) |,|Δs (k) |,|Δf (k) |}>λ , iteration times k=k+1, corrected state quantity V (k+1) = V (k) +ΔV (k) , θ (k+1) = θ (k) +Δθ (k) , s (k+ 1) =s (k) +Δs (k) , Δf (k+1) =Δf (k) +ΔΔf (k) , until the condition is met, output the result, where k is the number of iterations.
上述步骤(2)中,初始化的内容包括:设置迭代精度λ、最大迭代次数以及风力发电机滑差初值和频率偏差初值,形成节点导纳矩阵。In the above step (2), the initialization content includes: setting the iteration accuracy λ, the maximum number of iterations, and the initial value of the wind turbine slip and frequency deviation to form a node admittance matrix.
上述步骤(1)中,网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗,风电场的空气密度、风速,风机机型参数,系统初始频率。In the above step (1), the network parameters include: bus number, name, compensation capacitance, branch number of the transmission line, head node and end node numbers, series resistance, series reactance, parallel conductance, parallel susceptance, transformer transformation ratio and impedance, the air density and wind speed of the wind farm, the model parameters of the fan, and the initial frequency of the system.
上述步骤(1)中,量测量z包括:节点电压幅值、节点注入有功功率和无功功率,普通线路支路和变压器支路的有功功率和无功功率。In the above step (1), the quantity measurement z includes: node voltage amplitude, node injected active power and reactive power, active power and reactive power of ordinary line branch and transformer branch.
上述计及频率变化的含风力发电的电力系统状态估计是在同步发电机准稳态模型、负荷准稳态模型和异步风机简化RX模型的基础上,将频率偏差作为新的状态量引入到状态估计过程中,并且对发电机和负荷节点构造了新的零注入功率量测,最终建立了考虑频率偏差的状态估计模型。算例仿真结果表明该模型可有效计及系统的频率偏差,状态估计结果精度得到提高,具有工程应用前景。The above-mentioned state estimation of the power system including wind power generation considering the frequency change is based on the quasi-steady-state model of the synchronous generator, the quasi-steady-state model of the load and the simplified RX model of the asynchronous wind turbine, and introduces the frequency deviation as a new state quantity into the state In the estimation process, new zero-injected power measurements are constructed for generators and load nodes, and a state estimation model considering frequency deviation is finally established. The simulation results of the example show that the model can effectively take into account the frequency deviation of the system, and the accuracy of the state estimation results is improved, which has the prospect of engineering application.
附图说明Description of drawings
图1为本发明方法流程图。Fig. 1 is a flow chart of the method of the present invention.
图2为本发明提出的计及频率变化的含风力发电的电力系统状态估计所应用的算例系统,采用的是IEEE-14节点系统。Fig. 2 is a calculation example system applied to the state estimation of the power system including wind power generation in consideration of the frequency change proposed by the present invention, and the IEEE-14 node system is adopted.
图3为具体实施方式中异步电机的Γ形简化等效电路。Fig. 3 is a Γ-shaped simplified equivalent circuit of an asynchronous motor in a specific embodiment.
图4为两种状态估计模型电压平均估计误差。Figure 4 shows the average voltage estimation error of the two state estimation models.
具体实施方式Detailed ways
本发明以异步电机的准稳态模型为例分析风机的频率特性。为了简化计算,本发明采用RX模型。容量较大(大于40kW)的异步电机,由于其X1<<Xm,且R1和Rm可以忽略不计,可近似等效为图3所示的Γ形电路。The invention analyzes the frequency characteristics of the fan by taking the quasi-stable model of the asynchronous motor as an example. In order to simplify the calculation, the present invention adopts the RX model. An asynchronous motor with a large capacity (greater than 40kW), because its X 1 <<X m , and R 1 and R m can be ignored, can be approximately equivalent to the Γ-shaped circuit shown in Figure 3.
图3中R1、R2、Rm分别为定子电阻、转子电阻、励磁电阻,X1、X2和Xm分别表示定子电抗、转子电抗和励磁电抗,s是滑差,Δf是频率偏差量,PG和QG分别表示有功功率值和无功功率值。In Figure 3, R 1 , R 2 , and R m are stator resistance, rotor resistance, and excitation resistance, respectively; X 1 , X 2 , and X m represent stator reactance, rotor reactance, and excitation reactance; s is slip; Δf is frequency deviation PG and Q G represent active power value and reactive power value respectively.
为了考虑电网侧频率偏差对同步发电机的影响,同步发电机采用下面的准稳态模型:In order to consider the impact of grid-side frequency deviation on the synchronous generator, the following quasi-steady-state model is adopted for the synchronous generator:
式中:PG和QG分别表示同步发电机输出的有功功率值和无功功率值,PG_set和QG_set分别是同步发电机初始的有功功率值和无功功率值,PR是额定有功功率值,RR是对应同步发电机的调速率,aQ和bQ是同步发电机无功出力对应的调节系数,Δf表示频率偏差量,即系统稳态时频率与额定值的偏差。In the formula: PG and Q G respectively represent the active power value and reactive power value output by the synchronous generator, PG_set and Q G_set are the initial active power value and reactive power value of the synchronous generator respectively, PR is the rated active power Power value, R R is the regulation rate of the corresponding synchronous generator, a Q and b Q are the adjustment coefficients corresponding to the reactive output of the synchronous generator, Δf represents the frequency deviation, that is, the deviation between the frequency and the rated value in the steady state of the system.
负荷的准稳态数学模型采用考虑频率变化的静态模型,其多项式模型可表示如下:The quasi-steady-state mathematical model of the load adopts a static model considering the frequency change, and its polynomial model can be expressed as follows:
式中:PL和QL分别表示该负荷的有功和无功值,PL_set和QL_set分别表示该负荷的有功和无功的初始值,Kp和Kq分别表示负荷有功和无功对应的调节效应系数。pp、pc、pz和qp、qc、qz表示负荷模型静态电压特性系数,VL和VLB分别是该负荷的电压运行值和额定电压值,Δf是频率偏差量。In the formula: PL and Q L represent the active and reactive values of the load respectively, PL_set and Q L_set represent the initial values of active and reactive power of the load respectively, K p and K q represent the corresponding load active and reactive power adjustment effect coefficient. p p , p c , p z and q p , q c , q z represent the static voltage characteristic coefficients of the load model, V L and V LB are the voltage operating value and rated voltage value of the load respectively, and Δf is the frequency deviation.
电力系统状态估计的量测方程为:The measurement equation for power system state estimation is:
z=h(x)+εz=h(x)+ε
式中:x为状态量(维数n=2N-1,N为节点数);z为量测量(维数m,m>n);h为m维非线性量测函数;ε为m维量测误差。In the formula: x is the state quantity (dimension n=2N-1, N is the number of nodes); z is the quantity measurement (dimension m, m>n); h is the m-dimensional nonlinear measurement function; ε is the m-dimensional measurement error.
按最小二乘准则建立的目标函数如下:The objective function established according to the least squares criterion is as follows:
min J(x)=[z-h(x)]TW[z-h(x)]min J(x)=[zh(x)] T W[zh(x)]
其中J是目标函数,T表示矩阵的转置,W为对角权重矩阵,Wii=1/σi 2,σi为标准差。Where J is the objective function, T represents the transpose of the matrix, W is the diagonal weight matrix, W ii =1/σ i 2 , and σ i is the standard deviation.
一般情况下,h(x)为非线性函数,故采用迭代的方法求解。令x0是x的某一近似值,可以在x0附近对h(x)进行泰勒展开,保留一次项,并忽略二次以上的非线性项,得到:In general, h(x) is a nonlinear function, so an iterative method is used to solve it. Let x 0 be an approximate value of x, Taylor expansion can be performed on h(x) near x 0 , the first-order items are retained, and the non-linear items above the second order are ignored, to obtain:
h(x)≈h(x0)+H(x0)Δxh(x)≈h(x 0 )+H(x 0 )Δx
式中Δx=x-x0,H(x)为h(x)的雅克比矩阵。将此式代入目标函数中,可得到:In the formula, Δx=xx 0 , H(x) is the Jacobian matrix of h(x). Substituting this formula into the objective function, we can get:
J(x)=[Δz-H(x0)Δx]TW[Δz-H(x0)Δx]J(x)=[Δz-H(x 0 )Δx] T W[Δz-H(x 0 )Δx]
式中Δz=z-h(x0),将上式展开配方得到:In the formula, Δz=zh(x 0 ), expand the above formula to get:
J(x)=ΔzT[W-WH(x0)Σ(x0)HT(x0)W]ΔzJ(x)=Δz T [W-WH(x 0 )Σ(x 0 )H T (x 0 )W]Δz
+[Δx-Σ(x0)HT(x0)WΔz]TΣ-1(x0)[Δx-Σ(x0)HT +[Δx-Σ(x 0 )H T (x 0 )WΔz] T Σ -1 (x 0 )[Δx-Σ(x 0 )H T
×(x0)WΔz]×(x 0 )WΔz]
式中Σ(x0)=[HT(x0)WH(x0)]-1。In the formula, Σ(x 0 )=[H T (x 0 )WH(x 0 )] -1 .
上式中右边第一项与Δx无关。因此,欲使J(x)极小,第二项应为0,从而有:The first term on the right side of the above formula has nothing to do with Δx. Therefore, to make J(x) extremely small, the second term should be 0, thus:
Δx(l)=[HT(x(l))WH(x(l))]-1HT(x(l))W[z-h(x(l))]Δx (l) =[H T (x (l) )WH(x (l) )] -1 H T (x (l) )W[zh(x (l) )]
x(l+1)=x(l)+Δx(l) x (l+1) = x (l) +Δx (l)
其中l表示迭代次数,x按上式进行迭代修正,直到目标函数接近最小为止。Where l represents the number of iterations, and x is iteratively corrected according to the above formula until the objective function is close to the minimum.
因为风电的接入,本发明的状态估计模型在基本加权最小二乘法的基础上,不仅引入滑差s作为状态量引入修正方程,而且还考虑了频率偏移量Δf。所以状态估计的修正量扩展到Δx=[Δθ ΔV Δs ΔΔf]T,θ表示电压相角,V表示电压幅值,得到新的含有s和Δf的分块扩展雅可比矩阵为:Because of the access of wind power, the state estimation model of the present invention not only introduces the slip s as a state quantity into the correction equation on the basis of the basic weighted least square method, but also considers the frequency offset Δf. Therefore, the correction amount of the state estimation is extended to Δx=[Δθ ΔV Δs ΔΔf] T , θ represents the voltage phase angle, V represents the voltage amplitude, and the new block extended Jacobian matrix containing s and Δf is obtained as:
式中:P和Q分别表示普通发电机对应的有功和无功。Pk和Qk分别为异步风机的接入节点k对应的有功和无功。其中和的维数与系统中风电场节点的数目相同。In the formula: P and Q represent the active power and reactive power corresponding to the ordinary generator, respectively. P k and Q k are the active power and reactive power corresponding to the access node k of the asynchronous fan, respectively. in and The dimension of is the same as the number of wind farm nodes in the system.
系统中传统节点注入功率表示为:The injected power of traditional nodes in the system is expressed as:
式中:Pi和Qi分别表示节点i注入的有功和无功;Vi和Vj分别表示节点i和j的电压幅值;θij是节点i到节点j的电压相角差;Gij和Bij则表示节点导纳阵中对应节点i和j之间的电导和电纳;n是系统节点总数。In the formula: P i and Q i represent the active power and reactive power injected by node i respectively; V i and V j represent the voltage amplitudes of nodes i and j respectively; θ ij is the voltage phase angle difference from node i to node j; G ij and B ij represent the conductance and susceptance between the corresponding nodes i and j in the nodal admittance array; n is the total number of nodes in the system.
当计及发电机频率特性时,构建发电机节点零注入功率,此类发电机节点零注入功率可表达为:When the generator frequency characteristics are considered, the zero injection power of the generator node is constructed, and the zero injection power of this type of generator node can be expressed as:
式中:PGi和QGi分别表示发电机i注入的有功和无功。In the formula: P Gi and Q Gi represent the active power and reactive power injected by generator i, respectively.
当计及负荷频率特性时,构建负荷节点的零注入功率,此时负荷节点的零注入功率可表达为:When the load frequency characteristics are considered, the zero injected power of the load node is constructed. At this time, the zero injected power of the load node can be expressed as:
式中:PLi和QLi分别表示负荷i注入的有功和无功。In the formula: P Li and Q Li represent the active and reactive power injected by load i, respectively.
由前面各式可推导出雅可比矩阵H的各元素,其中部分分块矩阵元素如下:The elements of the Jacobian matrix H can be deduced from the previous formulas, and some of the elements of the block matrix are as follows:
式中:Pki和Qki分别表示风机节点i注入的有功和无功,si表示风机节点i对应的滑差。其中X12=X1+X2,下标i表示风机节点i对应的参数。In the formula: P ki and Q ki represent the active and reactive power injected by fan node i, respectively, and s i represents the slip corresponding to fan node i. Where X 12 =X 1 +X 2 , the subscript i represents the parameter corresponding to the fan node i.
根据上面的公式由初始各状态量V(0)、θ(0)、s(0)、Δf(0)计算量测量的计算值h(x(k))和雅克比矩阵H(x(k)),k是迭代次数,求解出状态修正量Δx(k),然后判断是否满足收敛条件,如果未达到收敛要求,修正状态量V(k+1)=V(k)+ΔV(k),θ(k+1)=θ(k)+Δθ(k),s(k+1)=s(k)+Δs(k),Δf(k+1)=Δf(k)+ΔΔf(k),重复上述操作,直到收敛精度达到要求。According to the above formula , the calculated value h(x ( k ) ) and the Jacobian matrix H ( x (k ) ) ), k is the number of iterations, solve the state correction quantity Δx (k) , and then judge whether the convergence condition is met, if the convergence requirement is not met, correct the state quantity V (k+1) = V (k) + ΔV (k) , θ (k+1) = θ (k) + Δθ (k) , s (k+1) = s (k) + Δs (k) , Δf (k+1) = Δf (k) + ΔΔf (k ) , repeat the above operations until the convergence accuracy meets the requirements.
上述方法具体步骤如下:The specific steps of the above method are as follows:
(1)获得电力系统的网络参数和量测量。网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗,风电场的空气密度、风速,风机机型参数,系统初始频率;量测量z包括:节点电压幅值、节点注入有功功率和无功功率,普通线路支路和变压器支路的有功功率和无功功率;(1) Obtain network parameters and quantity measurements of the power system. Network parameters include: bus bar number, name, compensation capacitance, branch number of transmission line, head node and end node numbers, series resistance, series reactance, parallel conductance, parallel susceptance, transformer ratio and impedance, wind farm air Density, wind speed, fan model parameters, system initial frequency; measurement z includes: node voltage amplitude, node injected active power and reactive power, active power and reactive power of ordinary line branches and transformer branches;
(2)利用上面获得的参数进行状态估计程序初始化。初始化的内容包括:设置迭代精度λ、最大迭代次数,风力发电机滑差初值和频率偏差初值,形成节点导纳矩阵;(2) Use the parameters obtained above to initialize the state estimation program. The initialization content includes: setting the iteration accuracy λ, the maximum number of iterations, the initial value of the slip of the wind turbine and the initial value of the frequency deviation, and forming the node admittance matrix;
(3)建立计及频率变化的含风力发电的状态估计模型:(3) Establish a state estimation model including wind power generation that takes into account frequency changes:
min J(x)=[z-h(x)]TW[z-h(x)]min J(x)=[zh(x)] T W[zh(x)]
(4)考虑风机滑差s和系统频率偏差Δf建立新的雅克比矩阵:(4) Consider the fan slip s and the system frequency deviation Δf to establish a new Jacobian matrix:
(5)由初始各状态量V(0)、θ(0)、s(0)、Δf(0)计算量测量的计算值h(x(k))和雅克比矩阵H(x(k));(5) Calculated value h(x ( k) ) and Jacobian matrix H(x ( k ) );
(6)求解状态修正量Δx(k),判断是否满足收敛条件,若max{|ΔV(k)|,|Δθ(k)|,|Δs(k)|,|Δf(k)|}>λ,迭代次数k=k+1,修正状态量V(k+1)=V(k)+ΔV(k),θ(k+1)=θ(k)+Δθ(k),s(k+1)=s(k)+Δs(k),Δf(k+1)=Δf(k)+ΔΔf(k),直到符合条件,输出结果。(6) Solve the state correction amount Δx (k) to judge whether the convergence condition is satisfied, if max{|ΔV (k) |,|Δθ (k) |,|Δs (k) |,|Δf (k) |}> λ, number of iterations k=k+1, corrected state quantity V (k+1) = V (k) +ΔV (k) , θ (k+1) = θ (k) + Δθ (k) , s (k +1) =s (k) +Δs (k) , Δf (k+1) =Δf (k) +ΔΔf (k) , until the condition is met, output the result.
本发明采用计及频率变化的含风力发电的电力系统状态估计,通过算例仿真,验证了本发明提出的模型效果显著,并且在电力系统估计结果的精度上,优于没有考虑频率偏差的模型。The present invention adopts the state estimation of the power system including wind power generation that takes into account the frequency change, and through the example simulation, it is verified that the model proposed by the present invention has a significant effect, and the accuracy of the power system estimation result is better than the model that does not consider the frequency deviation .
下面介绍本发明的实施例:Introduce the embodiment of the present invention below:
设风电场的平均空气密度为1.225kg/m3;风力发电机的扫掠面积为2642m2,初始滑差为-0.0044;风力发电机的切入风速、额定风速和切出风速分别为3m/s,16m/s和21m/s;风能利用系数Cp为0.1217。风力发电机的型号为V52-850,具体参数见表1:Suppose the average air density of the wind farm is 1.225kg/m 3 ; the swept area of the wind turbine is 2642m 2 , and the initial slip is -0.0044; the cut-in wind speed, rated wind speed and cut-out wind speed of the wind turbine are 3m/s respectively , 16m/s and 21m/s; wind energy utilization coefficient C p is 0.1217. The model of the wind turbine is V52-850, and the specific parameters are shown in Table 1:
表1V52-850机型参数Table 1 V52-850 model parameters
算例:Examples:
本发明采用图2所示的IEEE-14节点系统,为了对比两种模型的估计精度,当40台上述风力发电机构成的风电机组接入IEEE-14节点标准测试系统的3号节点,仿真结果如下表所示:The present invention adopts the IEEE-14 node system shown in Figure 2. In order to compare the estimation accuracy of the two models, when the wind turbines composed of 40 above-mentioned wind generators are connected to No. 3 node of the IEEE-14 node standard test system, the simulation results As shown in the table below:
表2估计结果及与真值的比较Table 2 Estimated results and comparison with the true value
由图4可知:未考虑频率偏差时,电压幅值和相角的平均估计误差分别为0.3570%和3.0257%;考虑频率偏差时,电压幅值和相角的平均估计误差分别为0.2662%和1.1151%。在估计结果上,未考虑频率的状态估计模型误差明显大于本文考虑频率的状态估计模型,充分说明了本文模型更能反映含风机系统的真实运行情况。It can be seen from Figure 4 that when the frequency deviation is not considered, the average estimation errors of the voltage amplitude and phase angle are 0.3570% and 3.0257% respectively; when the frequency deviation is considered, the average estimation errors of the voltage amplitude and phase angle are 0.2662% and 1.1151 respectively %. In the estimation results, the error of the state estimation model without considering the frequency is significantly larger than that of the state estimation model considering the frequency in this paper, which fully shows that the model in this paper can better reflect the real operation of the system with fans.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510010384.XA CN104598728B (en) | 2015-01-08 | 2015-01-08 | A kind of meter and the power system state estimation method containing wind-power electricity generation of frequency change |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510010384.XA CN104598728B (en) | 2015-01-08 | 2015-01-08 | A kind of meter and the power system state estimation method containing wind-power electricity generation of frequency change |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104598728A true CN104598728A (en) | 2015-05-06 |
CN104598728B CN104598728B (en) | 2017-11-03 |
Family
ID=53124507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510010384.XA Active CN104598728B (en) | 2015-01-08 | 2015-01-08 | A kind of meter and the power system state estimation method containing wind-power electricity generation of frequency change |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104598728B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095659A (en) * | 2015-07-27 | 2015-11-25 | 国电南瑞科技股份有限公司 | Method for province and region coordinative distributed state estimation based on cloud computing |
CN105633984A (en) * | 2016-03-08 | 2016-06-01 | 广州供电局有限公司 | Frequency state parameter detection method |
CN106199193A (en) * | 2016-06-30 | 2016-12-07 | 华北电力科学研究院有限责任公司 | Double-fed blower fan impedance hardware-in-the-loop test system and method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008109810A (en) * | 2006-10-26 | 2008-05-08 | Chugoku Electric Power Co Inc:The | Evaluation index analysis method for control of load frequency |
JP2009112148A (en) * | 2007-10-31 | 2009-05-21 | Toshiba Corp | Method and system of adjusting frequency of power system |
CN101707373A (en) * | 2009-11-20 | 2010-05-12 | 河海大学 | Automatic differentiation based power system state estimation method |
CN102522746A (en) * | 2011-12-13 | 2012-06-27 | 河海大学 | VSC-HVDC AC/DC optimal power flow method based on primal-dual interior point algorithm |
CN102545205A (en) * | 2011-12-22 | 2012-07-04 | 河海大学 | Distributed power system state estimating method based on automatic differentiation technology |
CN102623993A (en) * | 2012-04-12 | 2012-08-01 | 河海大学 | A State Estimation Method for Distributed Power System |
-
2015
- 2015-01-08 CN CN201510010384.XA patent/CN104598728B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008109810A (en) * | 2006-10-26 | 2008-05-08 | Chugoku Electric Power Co Inc:The | Evaluation index analysis method for control of load frequency |
JP2009112148A (en) * | 2007-10-31 | 2009-05-21 | Toshiba Corp | Method and system of adjusting frequency of power system |
CN101707373A (en) * | 2009-11-20 | 2010-05-12 | 河海大学 | Automatic differentiation based power system state estimation method |
CN102522746A (en) * | 2011-12-13 | 2012-06-27 | 河海大学 | VSC-HVDC AC/DC optimal power flow method based on primal-dual interior point algorithm |
CN102545205A (en) * | 2011-12-22 | 2012-07-04 | 河海大学 | Distributed power system state estimating method based on automatic differentiation technology |
CN102623993A (en) * | 2012-04-12 | 2012-08-01 | 河海大学 | A State Estimation Method for Distributed Power System |
Non-Patent Citations (4)
Title |
---|
唐西胜 等: "风力发电的调频技术研究综述", 《中国电机工程学报》 * |
崔明建 等: "考虑电网侧频率偏差的风电功率爬坡事件预测方法", 《电力系统自动化》 * |
庞博 等: "基于风力发电机简化RX模型的电力系统状态估计", 《电网技术》 * |
林今 等: "采用变速恒频机组的风电场有功功率波动对系统节点频率影响的动态评估模型", 《电力自动化设备》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095659A (en) * | 2015-07-27 | 2015-11-25 | 国电南瑞科技股份有限公司 | Method for province and region coordinative distributed state estimation based on cloud computing |
CN105095659B (en) * | 2015-07-27 | 2017-12-01 | 国电南瑞科技股份有限公司 | Coordinate distributed state estimation method to province based on cloud computing |
CN105633984A (en) * | 2016-03-08 | 2016-06-01 | 广州供电局有限公司 | Frequency state parameter detection method |
CN106199193A (en) * | 2016-06-30 | 2016-12-07 | 华北电力科学研究院有限责任公司 | Double-fed blower fan impedance hardware-in-the-loop test system and method |
CN106199193B (en) * | 2016-06-30 | 2020-12-08 | 华北电力科学研究院有限责任公司 | Doubly-fed fan impedance hardware-in-the-loop test system and method |
Also Published As
Publication number | Publication date |
---|---|
CN104598728B (en) | 2017-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512502B (en) | One kind is based on the normalized weight function the least square estimation method of residual error | |
CN107086606A (en) | A Parameter Identification Method for the Integrated Equivalent Asynchronous Motor Load Model of Distribution Network | |
Padrón et al. | Calculating steady-state operating conditions for doubly-fed induction generator wind turbines | |
CN103094900B (en) | Distributed generation power distribution network three-phase load flow calculation method taking phase sequence mixing method into consideration | |
CN104037763B (en) | A fast decomposition method power flow calculation method suitable for branch systems with small impedance | |
CN112214905B (en) | A wideband modeling analysis and simulation method and system for power system | |
CN104917174B (en) | Static voltage stability judgment method for active power distribution network | |
CN103414212B (en) | Consider the analogy method of the distribution network systems of double-fed wind generator machine access | |
CN111082427A (en) | A power flow calculation method for microgrid based on holomorphic function | |
CN114498678B (en) | An online evaluation method of power system inertia based on frequency-space correlation | |
Eminoglu | Modeling and application of wind turbine generating system (WTGS) to distribution systems | |
CN103532137A (en) | Method for estimating state of three-phase four-wire low-voltage distribution network | |
CN112234598B (en) | A method for initializing electromagnetic transient simulation | |
CN102545252B (en) | Voltage source commutation-high voltage direct current (VSC-HVDC) power flow computing method based on three-stage convergence Newton method | |
CN106446458A (en) | Weakly looped power distribution network load flow calculation method considering distributed power supplies | |
CN104794531A (en) | MDCOPF (modified direct current optimal power flow) method based on grid loss equivalent load model | |
CN104598728B (en) | A kind of meter and the power system state estimation method containing wind-power electricity generation of frequency change | |
CN110429648A (en) | Consider the small interference stability nargin probability evaluation method of failure of wind speed random fluctuation | |
CN108988391A (en) | The method for analyzing stability of double-fed fan rotor side converter based on revolving speed control | |
CN104934973B (en) | Power line loss calculation method considering temperature | |
CN106026082A (en) | Three-phase load flow calculation method for ungrounded power distribution network | |
CN105529736B (en) | Crow bar puts into double-fed fan motor group of planes fault current Equivalent Model computational methods | |
CN112800605B (en) | A Vector Wind Farm Modeling Method | |
CN109149583A (en) | Active power distribution network is succinctly pushed forward back substitution tidal current computing method | |
CN103199524B (en) | A kind of tidal current computing method adapting to the access of multiple distributed power source |
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 |