High-voltage transmission line electromagnetic transient numerical calculation method based on dimension-increasing fine integration

Abstract

The invention discloses a high-voltage transmission line electromagnetic transient numerical calculation method based on dimension-increasing fine integration, and the method comprises the following steps: S100, obtaining high-voltage transmission line information data which comprises size information and feature information, and the feature information comprises insulator position information, ground wire position information and line material information; s200, performing calculation processing according to the size information and the line material information to obtain line resistance data, line capacitance data and line inductance data; s300, performing modeling processing on the high-voltage power transmission line according to the high-voltage power transmission line information data, the line resistance data, the line capacitance data and the line inductance data to obtain an electromagnetic transient mathematical model of the high-voltage power transmission line; and S400, performing analysis processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics.

Description

Technical Field

[0001] The invention relates to the technical field of electromagnetic transient calculation of power systems, in particular to a numerical calculation method of electromagnetic transient of high-voltage transmission lines based on dimension-increased precise integration. Background Art

[0002] High-voltage transmission lines are one of the important components of the power system. Using numerical calculation methods to analyze the electromagnetic transient process of high-voltage transmission lines can solve problems in many areas of the power system: it can simulate transmission line faults to obtain short-circuit currents and harmonic components; it can be used to study transmission line fault distance algorithms and calculate overvoltages when the transmission line is closed without load.

[0003] The existing numerical calculation method for electromagnetic transients of high-voltage transmission lines does not deal with singular integrals in electromagnetic transients. Since high-voltage transmission lines are slender structures, the accuracy of the numerical calculation results of electromagnetic transients of high-voltage transmission lines does not meet the requirements. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide a numerical calculation method for electromagnetic transient of high-voltage transmission line based on dimensional increase and precise integration, which solves the problem that the existing numerical calculation method for electromagnetic transient of high-voltage transmission line does not process the singular integral in electromagnetic transient, and the accuracy of the numerical calculation results of electromagnetic transient of high-voltage transmission line does not meet the requirements because the high-voltage transmission line is a slender structure.

[0005] In order to solve the above technical problems, the technical solution adopted by the present invention is: a numerical calculation method of electromagnetic transient state of high-voltage transmission line based on dimensional increase precise integration, comprising the following steps: S100, obtaining high-voltage transmission line information data, wherein the high-voltage transmission line information data includes size information and feature information, and the feature information includes insulator position information, ground wire position information, and line material information; S200, performing calculations and processing according to the size information and the line material information to obtain line resistance data, line capacitance data, and line inductance data; S300, modeling the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data to obtain an electromagnetic transient mathematical model of the high-voltage transmission line; S400, analyzing and processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics.

[0006] Preferably, the steps of performing calculation and processing according to the size information and the line material information in step S200 to obtain the line resistance data, the line capacitance data, and the line inductance data are as follows: S201, analyzing and processing the size information to obtain the transmission line length data, the transmission line cross-sectional area, and the number of turns of the transmission line; S202, analyzing and processing the line material information to obtain the resistivity, dielectric constant, and relative magnetic permeability of the transmission line; S203, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the resistivity of the transmission line to obtain the transmission line resistance data; S204, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the dielectric constant of the transmission line to obtain the transmission line capacitance data; S205, performing correction calculation according to the transmission line length data, the cross-sectional area of ​​the transmission line, the number of turns of the transmission line, and the relative magnetic permeability of the transmission line to obtain the transmission line inductance data.

[0007] Preferably, the calculation formula for transmission line resistance data is:

[0008] Where, R is the resistance data of the transmission line; is the resistivity of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; A is the resistance compensation parameter of the transmission line; The calculation formula for transmission line capacitance data is:

[0009] Where C is the capacitance data of the transmission circuit, is the dielectric constant of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; B is the capacitance compensation parameter of the transmission line; The calculation formula for transmission line inductance data is:

[0010] Where, L is the inductance data of the transmission line; is the magnetic permeability of the transmission line in vacuum; is the relative magnetic permeability of the transmission line; is the number of turns of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; and E is the inductance compensation parameter of the transmission line.

[0011] Preferably, in step S300, modeling the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data to obtain the electromagnetic transient mathematical model of the high-voltage transmission line specifically includes the following steps: S301, obtaining geographic information of high-voltage transmission lines; S302, performing data modeling operations on the high-voltage transmission line and the geographic information of the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data, to obtain an electromagnetic transient mathematical model of the high-voltage transmission line.

[0012] Preferably, the step S400 of analyzing and processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics specifically includes the following steps: S401, discretizing the electromagnetic transient mathematical model of the high-voltage transmission line to obtain a minimum unit, wherein the minimum unit includes an insulator minimum unit and a transmission line minimum unit; S402, analyzing the minimum unit through Maxwell equations to obtain an electromagnetic field equation model corresponding to the minimum unit; S403, calculating and processing the electromagnetic field equation model through dimensionality-increasing precise integration to obtain current data, charge distribution data and electromagnetic characteristics.

[0013] Preferably, in step S401, the electromagnetic transient mathematical model of the high-voltage transmission line is discretized to obtain a minimum unit, wherein the minimum unit includes an insulator minimum unit and a transmission line minimum unit, and specifically includes the following steps: S4011. Divide the transmission line in the electromagnetic transient mathematical model of the high-voltage transmission line by using the finite element method to obtain the minimum unit of the transmission line; S4012. Divide the insulators in the electromagnetic transient mathematical model of the high-voltage transmission line by the finite element method to obtain the minimum unit of the insulator.

[0014] Preferably, the step S403 of calculating and processing the electromagnetic field equation model by dimensional increase and precise integration to obtain current data, charge distribution data and electromagnetic characteristics specifically includes the following steps: S4031. Calculate and process the electromagnetic field equation model in the time domain to obtain the time domain integral of the transmission line and the insulator; S4032. Calculate and process the electromagnetic field equation model in the spatial domain to obtain the spatial integral of the transmission line and the insulator; S4033. Calculate and process the time domain integral and the space domain integral through dimensionality-increasing fine integration to obtain current data, charge distribution data, and electromagnetic characteristics.

[0015] Preferably, the step S4033 calculates and processes the time domain integral and the space domain integral by means of dimension-increased fine integration to obtain current data, charge distribution data, and electromagnetic characteristics, specifically including the following steps: S40331. Analyze and process the time domain integral and the space domain integral to determine the time domain singular integral and the space domain singular integral; S40332, transforming the time domain singular integral and the space domain singular integral by means of dimensionality-increasing fine integration to obtain the time domain regular integral and the space domain regular integral; S40333. Analyze and process the time domain regularized integral and the space domain regularized integral to obtain current data, charge distribution data, and electromagnetic characteristics.

[0016] The present invention provides a method for numerical calculation of electromagnetic transient of high-voltage transmission line based on dimensional increase and precise integration. The method analyzes the high-voltage transmission line to determine the capacitance data, resistance data and inductance data of the transmission line, and models the model in combination with the size information and characteristic information of the high-voltage transmission line to obtain the electromagnetic field equation model of the high-voltage transmission line. The electromagnetic field equation model is analyzed in the time domain and the space domain, and the singular integral in the electromagnetic field equation model is converted into a regular integral in combination with dimensional increase and precise integration, so that the numerical calculation result of electromagnetic transient is more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The present invention will be further described below in conjunction with the accompanying drawings and embodiments: Figure 1 is a flow chart of the method of the present invention; Figure 2 It is a flow chart of steps S201-S205 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration proposed by the present invention; Figure 3 It is a flow chart of steps S301-S302 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration proposed by the present invention; Figure 4 It is a flow chart of steps S401-S403 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration proposed by the present invention; Figure 5 It is a flow chart of steps S4011-S4012 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on dimension-increased precise integration proposed by the present invention; Figure 6 It is a flow chart of steps S4031-S4033 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration proposed by the present invention; Figure 7 It is a flow chart of steps S40331-S40333 in the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration proposed by the present invention; Figure 8 The electromagnetic transient calculation results (initial phase angle) of the high-voltage transmission line under the increased-dimensional precise integration ); Fig. 9The electromagnetic transient calculation results of the high-voltage transmission line under the CDA method (initial phase angle ). DETAILED DESCRIPTION

[0018] like Figure 1 As shown in the figure, the numerical calculation method of electromagnetic transient of high-voltage transmission line based on increased-dimensional precise integration includes: S100, obtaining high-voltage transmission line information data, wherein the high-voltage transmission line information data includes size information and feature information, and the feature information includes insulator position information, ground wire position information, and line material information; S200, performing calculations and processing according to the size information and the line material information to obtain line resistance data, line capacitance data, and line inductance data; S300, modeling the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data to obtain an electromagnetic transient mathematical model of the high-voltage transmission line; S400, analyzing and processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics; It will be understood by those skilled in the art that to perform numerical calculations of electromagnetic transients on high-voltage transmission lines, it is necessary to establish an electromagnetic field equation model of the high-voltage transmission lines, integrate the electromagnetic field equation model in the time domain and space domain, determine the time domain singular integrals and space domain singular integrals, convert the time domain singular integrals and space domain singular integrals through dimensionality-increasing fine integration, convert the singular integrals into regular integrals, and then process the regular integrals to obtain circuit data, charge distribution, and electromagnetic characteristics of the high-voltage transmission lines.

[0019] Reference Figure 2 As shown, the calculation and processing are performed according to the size information and the line material information to obtain the transmission line resistance data, the transmission line capacitance data, and the transmission line inductance data, which specifically includes the following steps: S201, analyzing and processing the size information to obtain the transmission line length data, the transmission line cross-sectional area, and the number of turns of the transmission line; S202, analyzing and processing the line material information to obtain the resistivity, dielectric constant, and relative magnetic permeability of the transmission line; S203, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the resistivity of the transmission line to obtain the transmission line resistance data; The calculation formula for transmission line resistance data is:

[0020] Where, R is the resistance data of the transmission line; is the resistivity of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; A is the resistance compensation parameter of the transmission line; S204, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the dielectric constant of the transmission line to obtain the transmission line capacitance data; The calculation formula for transmission line capacitance data is:

[0021] Where C is the capacitance data of the transmission line, is the dielectric constant of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; B is the capacitance compensation parameter of the transmission line; S205, performing correction calculation according to the transmission line length data, the cross-sectional area of ​​the transmission line, the number of turns of the transmission line, and the relative magnetic permeability of the transmission line to obtain the transmission line inductance data; The calculation formula for transmission line inductance data is:

[0022] Where, L is the inductance data of the transmission line; is the magnetic permeability of the transmission line in vacuum; is the relative magnetic permeability of the transmission line; is the number of turns of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; E is the inductance compensation parameter of the transmission line; It will be understood by those skilled in the art that the size and material of the high-voltage transmission line are analyzed to determine the resistance data, capacitance data and inductance data of the high-voltage transmission line. However, due to the different use environments of the high-voltage transmission lines and the impact of the external environment on the resistance data, inductance data and capacitance data of the high-voltage transmission line, the transmission line resistance compensation parameter A, the transmission line capacitance compensation parameter B and the transmission line inductance compensation parameter E are introduced to correct the resistance data, inductance data and capacitance data, so that the subsequent electromagnetic transient calculation results are more accurate.

[0023] Reference Figure 3 As shown, the high-voltage transmission line is modeled according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data, and the electromagnetic transient mathematical model of the high-voltage transmission line is obtained, which specifically includes the following steps: S301, obtaining geographic information of high-voltage transmission lines; S302, performing data modeling operations on the high-voltage transmission line and the geographic information of the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data, to obtain an electromagnetic transient mathematical model of the high-voltage transmission line.

[0024] Reference Figure 4 As shown, the analysis and processing based on the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics specifically includes the following steps: S401, discretizing the electromagnetic transient mathematical model of the high-voltage transmission line to obtain a minimum unit, wherein the minimum unit includes an insulator minimum unit and a transmission line minimum unit; The transmission lines and insulators are divided into units through the finite element method, the structure of the transmission lines, and the insulator structure. Each unit is analyzed through Maxwell's equations to determine the electromagnetic field equation model corresponding to each unit. The electromagnetic field equation model is then transformed based on the dimensionality-increasing fine integration to determine the current data, charge distribution data, and electromagnetic characteristics.

[0025] S402, analyzing the minimum unit through Maxwell equations to obtain an electromagnetic field equation model corresponding to the minimum unit; S403, calculating and processing the electromagnetic field equation model through dimensionality-increasing precise integration to obtain current data, charge distribution data and electromagnetic characteristics.

[0026] Reference Figure 5 As shown in FIG. 1 , the electromagnetic transient mathematical model of the high-voltage transmission line is discretized to obtain the minimum unit of the transmission line, which specifically includes the following steps: S4011. Divide the transmission line in the electromagnetic transient mathematical model of the high-voltage transmission line by using the finite element method to obtain the minimum unit of the transmission line; S4012. Divide the insulators in the electromagnetic transient mathematical model of the high-voltage transmission line by the finite element method to obtain the minimum unit of the insulator.

[0027] Reference Figure 6 As shown, the electromagnetic field equation model is calculated and processed by increasing the dimension and fine integration to obtain the current data, charge distribution data and electromagnetic characteristics, which specifically includes the following steps: S4031. Calculate and process the electromagnetic field equation model in the time domain to obtain the time domain integral of the transmission line and the insulator; S4032. Calculate and process the electromagnetic field equation model in the spatial domain to obtain the spatial integral of the transmission line and the insulator; S4033. Calculate and process the time domain integral and the space domain integral through dimensionality-increasing fine integration to obtain current data, charge distribution data, and electromagnetic characteristics.

[0028] Reference Figure 7 As shown, the time domain integral and the space domain integral are calculated and processed by the dimension-increased fine integration to obtain the current data, charge distribution data, and electromagnetic characteristics, which specifically includes the following steps: S40331. Analyze and process the time domain integral and the space domain integral to determine the time domain singular integral and the space domain singular integral; S40332, transforming the time domain singular integral and the space domain singular integral by means of dimensionality-increasing fine integration to obtain the time domain regular integral and the space domain regular integral; S40333. Analyze and process the time domain regularized integral and the space domain regularized integral to obtain current data, charge distribution data, and electromagnetic characteristics; Those skilled in the art will understand that the reason why the electromagnetic field equation model is integrated in the time domain and the space domain is that the current in the high-voltage transmission line is a conduction process, which changes with time and space. Therefore, the electromagnetic fields at different positions of different high-voltage transmission lines also change with time. Therefore, the electromagnetic field equation model is integrated in the time domain and the space domain, and then the time domain integral and the space domain integral are analyzed to determine the singular integrals in the time domain integral and the space domain integral, and then the singular integrals are converted through dimensional increase and fine integration to convert the singular integrals into regular integrals, and then the time domain integral and the space domain integral that do not contain the singular integrals are differentiated and restored to the electromagnetic field equation model, and the current data, charge distribution and electromagnetic characteristics are determined according to the electromagnetic field equation model.

[0029] Reference Figure 8 , 9 As shown in the figure, taking the sinusoidal signal as the input of the high-voltage transmission line as an example, the calculation formula is:

[0030] In the formula, ; a is a constant; is the phase angle, which is a known quantity; according to For sinusoidal information, new variables are introduced to increase the dimension of the state variables. The dimension increase formula is:

[0031] Using the dimension-increasing precise integration method, when When the voltage at the end of the high-voltage transmission line changes with time, the curve is as follows: Figure 8 ; Using the CDA method, when When the voltage at the end of the high-voltage transmission line changes with time, the curve is as follows: Fig. 9 ; In summary, compared with the CDA method, the increased-dimensional precise integration method is more effective in suppressing numerical oscillation problems.

[0032] The above embodiments are only preferred technical solutions of the present invention and should not be regarded as limiting the present invention. The protection scope of the present invention shall be the technical solutions recorded in the claims, including equivalent replacement solutions of the technical features in the technical solutions recorded in the claims. That is, equivalent replacement improvements within this scope are also within the protection scope of the present invention.

Claims

1. A numerical calculation method for electromagnetic transient of high-voltage transmission line based on dimensional-increased precise integration, characterized in that: The following steps are involved: S100, obtaining high-voltage transmission line information data, wherein the high-voltage transmission line information data includes size information and feature information, and the feature information includes insulator position information, ground wire position information, and line material information; S200, performing calculations and processing according to the size information and the line material information to obtain line resistance data, line capacitance data, and line inductance data; S300, modeling the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data to obtain an electromagnetic transient mathematical model of the high-voltage transmission line; S400, analyzing and processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics.

2. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 1 is characterized in that: In step S200, the steps of performing calculation and processing according to the size information and the line material information to obtain the line resistance data, the line capacitance data, and the line inductance data are as follows: S201, analyzing and processing the size information to obtain the transmission line length data, the transmission line cross-sectional area, and the number of turns of the transmission line; S202, analyzing and processing the line material information to obtain the resistivity, dielectric constant, and relative magnetic permeability of the transmission line; S203, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the resistivity of the transmission line to obtain the transmission line resistance data; S204, performing correction calculation according to the transmission line length data, the transmission line cross-sectional area, and the dielectric constant of the transmission line to obtain the transmission line capacitance data; S205, performing correction calculation according to the transmission line length data, the cross-sectional area of ​​the transmission line, the number of turns of the transmission line, and the relative magnetic permeability of the transmission line to obtain the transmission line inductance data.

3. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 2 is characterized in that: The calculation formula for transmission line resistance data is: In the formula, R is the resistance data of the transmission line; is the resistivity of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; A is the resistance compensation parameter of the transmission line; The calculation formula for transmission line capacitance data is: Where C is the capacitance data of the transmission circuit, is the dielectric constant of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; B is the capacitance compensation parameter of the transmission line; The calculation formula for transmission line inductance data is: Where, L is the inductance data of the transmission line; is the magnetic permeability of the transmission line in vacuum; is the relative magnetic permeability of the transmission line; is the number of turns of the transmission line; D is the length of the transmission line; S is the cross-sectional area of ​​the transmission line; and E is the inductance compensation parameter of the transmission line.

4. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 1 is characterized in that: In step S300, the high-voltage transmission line is modeled according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data to obtain the electromagnetic transient mathematical model of the high-voltage transmission line, which specifically includes the following steps: S301, obtaining geographic information of high-voltage transmission lines; S302, performing data modeling operations on the high-voltage transmission line and the geographic information of the high-voltage transmission line according to the high-voltage transmission line information data, line resistance data, line capacitance data, and line inductance data, to obtain an electromagnetic transient mathematical model of the high-voltage transmission line.

5. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 1 is characterized in that: The step S400 performs analysis and processing according to the electromagnetic transient mathematical model of the high-voltage transmission line to obtain current data, charge distribution data and electromagnetic characteristics, which specifically includes the following steps: S401, discretizing the electromagnetic transient mathematical model of the high-voltage transmission line to obtain a minimum unit, wherein the minimum unit includes an insulator minimum unit and a transmission line minimum unit; S402, analyzing the minimum unit through Maxwell equations to obtain an electromagnetic field equation model corresponding to the minimum unit; S403, calculating and processing the electromagnetic field equation model through dimensionality-increasing precise integration to obtain current data, charge distribution data and electromagnetic characteristics.

6. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 5 is characterized in that: In step S401, the electromagnetic transient mathematical model of the high-voltage transmission line is discretized to obtain the minimum unit, which includes the minimum insulator unit and the minimum transmission line unit. Specifically, the steps include: S4011. Divide the transmission line in the electromagnetic transient mathematical model of the high-voltage transmission line by using the finite element method to obtain the minimum unit of the transmission line; S4012. Divide the insulators in the electromagnetic transient mathematical model of the high-voltage transmission line by the finite element method to obtain the minimum unit of the insulator.

7. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 5 is characterized in that: The step S403 calculates and processes the electromagnetic field equation model by means of dimensional increase and precise integration to obtain current data, charge distribution data and electromagnetic characteristics, and specifically includes the following steps: S4031. Calculate and process the electromagnetic field equation model in the time domain to obtain the time domain integral of the transmission line and the insulator; S4032. Calculate and process the electromagnetic field equation model in the spatial domain to obtain the spatial integral of the transmission line and the insulator; S4033. Calculate and process the time domain integral and the space domain integral through dimensionality-increasing fine integration to obtain current data, charge distribution data, and electromagnetic characteristics.

8. The method for numerical calculation of electromagnetic transient state of high-voltage transmission line based on dimension-increased precise integration according to claim 7 is characterized in that: The step S4033 calculates and processes the time domain integral and the space domain integral by means of dimension-increasing fine integration to obtain current data, charge distribution data, and electromagnetic characteristics, specifically including the following steps: S40331. Analyze and process the time domain integral and the space domain integral to determine the time domain singular integral and the space domain singular integral; S40332, transforming the time domain singular integral and the space domain singular integral by dimensionality-increasing fine integration to obtain the time domain regular integral and the space domain regular integral; S40333. Analyze and process the time domain regularized integral and the space domain regularized integral to obtain current data, charge distribution data, and electromagnetic characteristics.

Images