A method and device for estimating resistance parameters of a transmission line based on electrical quantity data

Abstract

The present invention relates to a method and device for estimating the resistance parameters of a transmission line based on electric quantity data, comprising: obtaining line loss electric quantity collected values ​​and calculated values ​​of the transmission line in each time period, determining a line loss electric quantity collected value matrix of the transmission line based on the line loss electric quantity collected values ​​of the transmission line in each time period, determining a line loss electric quantity calculated value matrix of the transmission line based on the line loss electric quantity calculated values ​​of the transmission line in each time period, and determining an estimated resistance value of the transmission line using the line loss electric quantity collected value matrix and the calculated value matrix of the transmission line. The present invention improves the accuracy and practicality of the estimated resistance value of the transmission line, thereby providing a guarantee for the safe operation of the power grid.

Description

Technical Field

[0001] The present invention relates to the field of power system calculation, and in particular to a method and device for estimating resistance parameters of a power transmission line based on electrical quantity data. Background Art

[0002] In recent years, with the interconnection and expansion of my country's power grid, transmission line parameters, as a crucial component of grid parameters, have become increasingly accurate, crucial for the safe and stable operation of the power system. Currently, transmission line parameters used in online applications such as state estimation and power flow calculations in dispatching automation systems are generally fixed and unadjusted year-round. However, in actual grid operation, transmission line parameters vary with the operating environment. In particular, transmission line resistance parameters increase or decrease with temperature fluctuations. Theoretical calculations show that if a transmission line's temperature fluctuates by 30°C, its actual operating resistance can vary by more than 10%.

[0003] Currently, lean grid operation generally aims to minimize active power losses. This lean grid operation, which aims to minimize losses, is based on transmission line resistance parameters. However, due to the low accuracy of traditional transmission line resistance estimation methods, grid loss analysis deviates significantly from actual results, significantly impacting dispatch plans, grid operation mode arrangements, and lean grid operation. Summary of the Invention

[0004] In view of the shortcomings of the existing technology, the purpose of the present invention is to provide a method and device for estimating the resistance parameters of a transmission line based on electricity data. According to the line loss electricity collection value matrix and the calculated value matrix of the transmission line, a transmission line line loss electricity error equation is established, and the error equation is solved to obtain the estimated value of the resistance of the transmission line, thereby improving the accuracy of the transmission line resistance estimation value and ensuring the safe operation of the power grid.

[0005] The purpose of the present invention is achieved by adopting the following technical solutions:

[0006] The present invention provides a method for estimating resistance parameters of a transmission line based on electrical quantity data, wherein the method comprises:

[0007] Obtain the collected and calculated values ​​of power loss of the transmission line in each period;

[0008] Determine a line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each period;

[0009] Determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each period;

[0010] The resistance estimation value of the transmission line is determined using the line loss electricity collection value matrix and the calculated value matrix of the transmission line.

[0011] Preferably, the acquisition of the collected and calculated values ​​of the line loss electricity of the transmission line in each time period includes:

[0012] Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h :

[0013] ΔW s,h =(W s,h,z,q -W s,h,f,m )+(W s,h,z,m -W s,h,f,q )

[0014] The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h :

[0015] ΔW j,h =A h ·R+ΔW c

[0016] In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,m is the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0017] Furthermore, the integral area A of the square curve of the impedance current in the hth period is determined by the following formula: h :

[0018]

[0019] In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

[0020] Furthermore, the transmission line impedance current I in the hth period is determined by the following formula: Z,h :

[0021]

[0022] In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, P m,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

[0023] Preferably, determining the line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each time period includes:

[0024] Determine the line loss power collection value matrix N of the transmission line according to the following formula:

[0025]

[0026] In the above formula, ΔW s,h is the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

[0027] Preferably, determining the line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each time period includes:

[0028] Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula:

[0029]

[0030] In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0031] Preferably, the determining of the estimated resistance value of the transmission line by using the line loss power collection value matrix and the calculated value matrix of the transmission line includes:

[0032] Establishing a line loss power error equation for transmission lines based on the line loss power collection value matrix and the calculated value matrix;

[0033] The least square method is used to solve the transmission line loss error equation to obtain the estimated resistance value of the transmission line.

[0034] Furthermore, the establishment of a transmission line loss power error equation based on the line loss power collection value matrix and the calculated value matrix of the transmission line includes:

[0035] Determine the transmission line line loss error equation as follows:

[0036] [E1,...,E h ,...,E H ] T =N-N'=N-MX

[0037] in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

[0038] Furthermore, the method of using the least squares method to solve the transmission line loss error equation to obtain the estimated resistance value of the transmission line includes:

[0039] Step 1. Solve using the least squares method Get an estimate of X in,

[0040] Step 2. Order but

[0041]

[0042] Step 3. Determine the estimated resistance R' of the transmission line as follows:

[0043] R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ;

[0044] In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,his the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0045] The present invention provides a device for estimating resistance parameters of a power transmission line based on electrical quantity data, wherein the device comprises:

[0046] The acquisition module is used to obtain the collected and calculated values ​​of the line loss electricity of the transmission line in each period;

[0047] A first determining module is used to determine a line loss power collection value matrix of a transmission line according to the line loss power collection value of the transmission line in each time period;

[0048] A second determining module is used to determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each time period;

[0049] The third determination module is used to determine the estimated resistance value of the transmission line by using the line loss power collection value matrix and the calculated value matrix of the transmission line.

[0050] Preferably, the acquisition module is used to:

[0051] Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h :

[0052] ΔW s,h =(W s,h,z,q -W s,h,f,m )+(W s,h,z,m -W s,h,f,q )

[0053] The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h :

[0054] ΔW j,h =A h ·R+ΔW c

[0055] In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,mis the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0056] Furthermore, the integral area A of the square curve of the impedance current in the hth period is determined by the following formula: h :

[0057]

[0058] In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

[0059] Furthermore, the transmission line impedance current I in the hth period is determined by the following formula: Z,h :

[0060]

[0061] In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, P m,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

[0062] Preferably, the first determining module is configured to:

[0063] Determine the line loss power collection value matrix N of the transmission line according to the following formula:

[0064]

[0065] In the above formula, ΔW s,his the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

[0066] Preferably, the second determining module is used to:

[0067] Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula:

[0068]

[0069] In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0070] Preferably, the third determining module includes:

[0071] An establishing unit is used to establish a line loss power error equation for a transmission line based on a line loss power collection value matrix and a calculated value matrix of the transmission line;

[0072] The acquisition unit is used to solve the transmission line line loss error equation using the least square method to obtain the estimated resistance value of the transmission line.

[0073] Furthermore, the establishing unit is used to:

[0074] Determine the transmission line line loss error equation as follows:

[0075] [E1,...,E h ,...,E H ] T =N-N'=N-MX

[0076] in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

[0077] Furthermore, the acquisition unit is used to:

[0078] Step 1. Solve using the least squares method Get an estimate of X in,

[0079] Step 2. Order but

[0080]

[0081] Step 3. Determine the estimated resistance R' of the transmission line as follows:

[0082] R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ;

[0083] In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0084] Compared with the closest prior art, the present invention has the following beneficial effects:

[0085] The present invention provides a method and device for estimating transmission line resistance parameters based on electricity data. This method obtains collected and calculated line loss values ​​for each time period, determines a line loss collection value matrix for the transmission line based on the collected line loss values ​​for each time period, and determines a line loss calculation value matrix for the transmission line based on the calculated line loss values ​​for each time period. The line loss collection value matrix and the calculated value matrix are then used to determine an estimated transmission line resistance. The technical solution provided by the present invention improves the accuracy and practicality of transmission line resistance estimation, ensuring the safe operation of the power grid. BRIEF DESCRIPTION OF THE DRAWINGS

[0086] Figure 1 This is a flow chart of a method for estimating resistance parameters of a transmission line based on electrical quantity data provided by the present invention;

[0087] Figure 2 This is a schematic diagram of an application scenario for obtaining line loss power collection values ​​and calculation values ​​in an embodiment provided by the present invention;

[0088] Figure 3 This is a schematic diagram of an application scenario for obtaining the impedance current of a transmission line in an embodiment provided by the present invention;

[0089] Figure 4 This is a structural diagram of a transmission line resistance parameter estimation device based on electrical quantity data provided by the present invention. DETAILED DESCRIPTION

[0090] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0091] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0092] The present invention provides a method for estimating the resistance parameters of a transmission line based on electrical quantity data. Figure 1 As shown, the method includes:

[0093] Obtain the collected and calculated values ​​of power loss of the transmission line in each period;

[0094] Determine a line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each period;

[0095] Determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each period;

[0096] The resistance estimation value of the transmission line is determined using the line loss electricity collection value matrix and the calculated value matrix of the transmission line.

[0097] In the preferred embodiment of the present invention, Figure 2 In the application scenario shown, the acquisition of the collected and calculated values ​​of the line loss electricity of the transmission line in each time period includes:

[0098] Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h :

[0099] ΔW s,h =(W s,h,z,q -W s,h,f,m )+(W s,h,z,m -Ws,h,f,q )

[0100] The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h :

[0101] ΔW j,h =A h ·R+ΔW c

[0102] In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,m is the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, and the above calculations are all performed using effective values.

[0103] Furthermore, the integral area A of the square curve of the impedance current in the hth period is determined by the following formula: h :

[0104]

[0105] In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

[0106] Further, in Figure 3 In the application scenario shown, the transmission line impedance current I in the hth period is determined by the following formula: Z,h :

[0107]

[0108] In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, Pm,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

[0109] In the most preferred embodiment of the present invention, the method of determining the line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each time period includes:

[0110] Determine the line loss power collection value matrix N of the transmission line according to the following formula:

[0111]

[0112] In the above formula, ΔW s,h is the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

[0113] In a preferred embodiment of the present invention, determining a transmission line loss electricity calculation value matrix according to the transmission line loss electricity calculation value in each time period includes:

[0114] Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula:

[0115]

[0116] In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0117] In a preferred embodiment of the present invention, the method of determining the estimated resistance value of the transmission line using the line loss power collection value matrix and the calculated value matrix of the transmission line includes:

[0118] Establishing a line loss power error equation for transmission lines based on the line loss power collection value matrix and the calculated value matrix;

[0119] The least square method is used to solve the transmission line loss error equation to obtain the estimated resistance value of the transmission line.

[0120] Furthermore, the establishment of a transmission line loss power error equation based on the line loss power collection value matrix and the calculated value matrix of the transmission line includes:

[0121] Determine the transmission line line loss error equation as follows:

[0122] [E1,...,E h ,...,E H ] T =N-N'=N-MX

[0123] in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

[0124] Furthermore, the method of using the least squares method to solve the transmission line loss error equation to obtain the estimated resistance value of the transmission line includes:

[0125] Step 1. Solve using the least squares method Get an estimate of X in,

[0126] Step 2. Order but

[0127]

[0128] Step 3. Determine the estimated resistance R' of the transmission line as follows:

[0129] R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ;

[0130] In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,his the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0131] The present invention provides a device for estimating resistance parameters of a transmission line based on electrical quantity data. Figure 4 As shown, the device includes:

[0132] The acquisition module is used to obtain the collected and calculated values ​​of the line loss electricity of the transmission line in each period;

[0133] A first determining module is used to determine a line loss power collection value matrix of a transmission line according to the line loss power collection value of the transmission line in each time period;

[0134] A second determining module is used to determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each time period;

[0135] The third determination module is used to determine the estimated resistance value of the transmission line by using the line loss power collection value matrix and the calculated value matrix of the transmission line.

[0136] In the most preferred embodiment of the present invention, the acquisition module is used to:

[0137] Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h :

[0138] ΔW s,h =(W s,h，z，q -W s,h,f,m )+(W s,h,z,m -W s,h,f,q )

[0139] The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h :

[0140] ΔW j,h =A h ·R+ΔW c

[0141] In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,m is the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW CIt is the corona loss of the transmission line.

[0142] The integral area A of the square curve of the impedance current in the hth period is determined by the following formula: h :

[0143]

[0144] In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

[0145] Furthermore, the transmission line impedance current I in the hth period is determined by the following formula: Z,h :

[0146]

[0147] In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, P m,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

[0148] In a preferred embodiment of the present invention, the first determining module is configured to:

[0149] Determine the line loss power collection value matrix N of the transmission line according to the following formula:

[0150]

[0151] In the above formula, ΔW s,h is the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

[0152] In a preferred embodiment of the present invention, the second determining module is configured to:

[0153] Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula:

[0154]

[0155] In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0156] In a preferred embodiment of the present invention, the third determining module includes:

[0157] An establishing unit is used to establish a line loss power error equation for a transmission line based on a line loss power collection value matrix and a calculated value matrix of the transmission line;

[0158] The acquisition unit is used to solve the transmission line line loss error equation using the least square method to obtain the estimated resistance value of the transmission line.

[0159] Furthermore, the establishing unit is used to:

[0160] Determine the transmission line line loss error equation as follows:

[0161] [E1,...,E h ,...,E H ] T =N-N'=N-MX

[0162] in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

[0163] Furthermore, the acquisition unit is used to:

[0164] Step 1. Solve using the least squares method Get an estimate of X in,

[0165] Step 2. Order but

[0166]

[0167] Step 3. Determine the estimated resistance R' of the transmission line as follows:

[0168] R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ;

[0169] In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

[0170] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0171] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0172] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0173] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0174] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered by the scope of protection of the claims of the present invention.

Claims

1. A method for estimating resistance parameters of a transmission line based on electrical quantity data, characterized in that: The method comprises: Obtain the collected and calculated values ​​of power loss of the transmission line in each period; Determine a line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each period; Determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each period; Determine the estimated resistance value of the transmission line using the line loss power collection value matrix and the calculated value matrix of the transmission line; The method of determining the estimated resistance value of the transmission line by using the line loss power collection value matrix and the calculated value matrix of the transmission line includes: Establishing a line loss power error equation for transmission lines based on the line loss power collection value matrix and the calculated value matrix; The least square method is used to solve the transmission line loss error equation to obtain the estimated resistance value of the transmission line; The method of using the least square method to solve the transmission line line loss error equation to obtain the estimated resistance value of the transmission line includes: Step 1. Solve using the least squares method Get an estimate of X in, Step 2. Order but Step 3. Determine the estimated resistance R' of the transmission line as follows: R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ; In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C A is the corona loss of the transmission line; h is the integrated area of ​​the square curve of the impedance current in the hth time period.

2. The method according to claim 1, wherein The acquisition of the collected and calculated values ​​of the line loss electricity of the transmission line in each time period includes: Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h : ΔW s,h =(W s,h,z,q -IN s,h,f,m )+(W s,h,z,m -IN s,h,f,q ) The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h : ΔW j,h =A h ·R+ΔW c In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,m is the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

3. The method according to claim 2, wherein Determine the integral area A of the square curve of the impedance current in the hth period by the following formula h : In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

4. The method according to claim 3, wherein the transmission line impedance current I in the hth period is determined by the following formula: Z,h : In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, P m,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

5. The method according to claim 1, wherein Determining the line loss power collection value matrix of the transmission line according to the line loss power collection value of the transmission line in each time period includes: Determine the line loss power collection value matrix N of the transmission line according to the following formula: In the above formula, ΔW s,h is the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

6. The method according to claim 1, wherein Determining the line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each time period includes: Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula: In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, h∈H, H is the total number of periods, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

7. The method according to claim 1, wherein The method of establishing a transmission line loss power error equation based on the line loss power collection value matrix and the calculated value matrix of the transmission line includes: Determine the transmission line line loss error equation as follows: [E1,...,E h ,...,HAVE BEEN H ] T =N-N'=N-MX in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

8. A device for estimating resistance parameters of a transmission line based on electrical quantity data, characterized in that: The device comprises: The acquisition module is used to obtain the collected and calculated values ​​of the line loss electricity of the transmission line in each period; A first determining module is used to determine a line loss power collection value matrix of a transmission line according to the line loss power collection value of the transmission line in each time period; A second determining module is used to determine a line loss electricity calculation value matrix of the transmission line according to the line loss electricity calculation value of the transmission line in each time period; A third determination module is used to determine the estimated resistance value of the transmission line using the line loss power collection value matrix and the calculated value matrix of the transmission line; The third determining module includes: An establishing unit is used to establish a line loss power error equation for a transmission line based on a line loss power collection value matrix and a calculated value matrix of the transmission line; An acquisition unit, configured to solve a transmission line line loss error equation using a least squares method to obtain an estimated resistance value of the transmission line; The acquisition unit is configured to: Step 1. Solve using the least squares method Get an estimate of X in, Step 2. Order but Step 3. Determine the estimated resistance R' of the transmission line as follows: R'=F 11 *ΔW s,1 +...+F 1h *ΔW s,h +...+F 1H *ΔW s,H ; In the above formula, E h is the line loss power error of the transmission line in the hth period, h∈H, H is the total number of periods, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the collected value of power loss of the transmission line in the hth period, F 1h is the hth element value of the first row of matrix F, X=[R,ΔW C ] T , N is the power loss matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

9. The device according to claim 8, wherein The acquisition module is used to: Determine the line loss power collection value ΔW of the transmission line in the hth period by the following formula: s,h : ΔW s,h =(W s,h,z,q -IN s,h,f,m )+(W s,h,z,m -IN s,h,f,q ) The calculated value of line loss ΔW of the transmission line in the hth period is determined by the following formula: j,h : ΔW j,h =A h ·R+ΔW c In the above formula, W s,h,z,q is the forward active power collected at the head end of the transmission line in the hth period, W s,h,f,q is the reverse active power collected at the head end of the transmission line in the hth period, W s,h,z,m is the forward active power collected at the end of the transmission line in the hth period, W s,h,f,m is the reverse active power collected at the end of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

10. The device according to claim 9, wherein Determine the integral area A of the square curve of the impedance current in the hth period by the following formula h : In the above formula, I Z,h is the transmission line impedance current in the hth period, T is the total number of moments in the hth period, I Z,h is the impedance current of the transmission line in the hth period, K is the total number of current measurement values ​​of the transmission line in the hth period, I Z,h,k+1 is the k+1th current measurement value of the transmission line in the hth period, I Z,h,k is the kth current measurement value of the transmission line in the hth time period, and Δt is the time interval between adjacent current measurements.

11. The device according to claim 10, wherein the transmission line impedance current I in the hth time period is determined by the following formula: Z,h : In the above formula, P q,h is the active power at the head end of the transmission line in the hth period, Q q,h is the active power at the head end of the transmission line in the hth period, U q,h is the bus voltage at the head end of the transmission line in the hth period, P m,h is the active power at the end of the transmission line in the hth period, Q m,h is the active power at the end of the transmission line in the hth period, U m,h is the bus voltage at the end of the transmission line in the hth period, Y c is the earth susceptance of the transmission line.

12. The device according to claim 8, wherein The first determining module is configured to: Determine the line loss power collection value matrix N of the transmission line according to the following formula: In the above formula, ΔW s,h is the collected value of the line loss electricity of the transmission line in the hth period, h∈H, H is the total number of periods.

13. The device according to claim 8, wherein The second determining module is configured to: Determine the line loss electricity calculation value matrix N' of the transmission line according to the following formula: In the above formula, ΔW j,h is the calculated value of the power loss of the transmission line in the hth period, A h is the integral area of ​​the square curve of the impedance current in the hth period, h∈H, H is the total number of periods, R is the rated resistance of the transmission line, ΔW C It is the corona loss of the transmission line.

14. The device according to claim 8, wherein The establishing unit is used to: Determine the transmission line line loss error equation as follows: [E1,...,E h ,...,HAVE BEEN H ] T =N-N'=N-MX in, X=[R,ΔW C ] T , N' is the line loss power calculation value matrix of the transmission line, N is the line loss power collection value matrix of the transmission line, ΔW s,h is the line loss power collection value of the transmission line in the hth period, h∈H, H is the total number of periods, A h is the integral area of ​​the square curve of the impedance current in the hth period, R is the rated resistance of the transmission line, ΔW C is the corona loss of the transmission line, E h is the transmission line loss error in the hth period.

Images