GIS electromagnetic field distribution simulation method, system and storage medium

Abstract

The invention discloses a GIS electromagnetic field distribution simulation method, system and storage medium. The method comprises the following steps: discretizing a three-dimensional simulation model of a GIS device into a spatial grid, and determining the corresponding relationship between the number and spatial coordinate data of each grid node; converting the Maxwell equations used to describe the GIS electromagnetic field into a fourth-order matrix, and using a four-step HIE-FDTD algorithm to solve the fourth-order matrix, and calculating the field strength data of each grid node; constructing a graph structure based on the spatial coordinate data and field strength data of the grid nodes; and generating new grid nodes and eigenvalues ​​of the new grid nodes between adjacent grid nodes using a node interpolation enhancement strategy based on the graph structure. The invention can quickly realize grid enhancement of a coarse grid simulation electromagnetic field distribution result, and also achieve a result of finer grid simulation.

Description

Technical Field

[0001] The present invention relates to the technical field of electric power equipment, and in particular to a GIS electromagnetic field distribution simulation method, system and storage medium. Background Art

[0002] Gas insulated switchgear (GIS) has the advantages of relatively small footprint, high reliability, and excellent breaking performance. It has been widely used in the world, especially in power systems. Years of operating experience of UHV GIS shows that despite its high operating reliability, various defects such as free metal particles, needle-like protrusions, fixed particles, insulating gas, etc. are still inevitable due to process problems in the manufacturing and assembly process. By detecting the distribution of electromagnetic fields in GIS, insulation faults can be discovered in time and the insulation state of GIS can be reliably evaluated so that appropriate measures can be taken to prevent further development and accidents. In this regard, there are currently two main methods: simulation experiments and simulation calculations. Due to the particularity of GIS equipment, experimental conditions are limited, and there is a lack of experimental data and characteristic research. Therefore, simulation has become an important means of studying GIS equipment failures.

[0003] At present, the main research methods of GIS simulation include Finite-Difference Time-Domain (FDTD) and Finite Element Method in Time Domain (FEMTD). Compared with FEMTD, FDTD is more widely used. The spatial sampling method of the electric field and magnetic field of the grid unit in FDTD directly satisfies Faraday's law and Ampere's law. Through time discretization and space discretization, the electric field at the current moment is obtained from the electric field at the previous moment and the magnetic field at the adjacent position, and the magnetic field at the current moment is obtained from the magnetic field at the previous moment and the electric field at the adjacent position. Therefore, only after the initial values ​​of the electric field and magnetic field of each grid unit are assigned, the calculation of the electromagnetic propagation process is completed by directly sampling the electric field and magnetic field alternately at half a time step.

[0004] However, during the simulation process, the coarseness of the grid division will affect the precision of the calculation results, but it will also make the time cost inversely proportional to the precision and increase exponentially. Therefore, the study of GIS time-domain finite-difference simulation and grid node data enhancement methods is important for GIS simulation, and can also provide a reference direction for the research and simulation of GIS equipment.

[0005] In the related technologies, the master's thesis "Simulation and Calculation of Transient Electric Field of GIS Equipment Based on Finite Difference Time Domain Method, Yang Rui, Huazhong University of Science and Technology" proposed a set of FDTD methods suitable for transient electric field calculation of GIS equipment. This paper selects appropriate grid generation algorithm and absorbing boundary conditions, divides GIS into two parts: coaxial bus pipeline and GIS complex equipment according to the structural characteristics during modeling, and studies key issues such as grid generation algorithm. Since GIS electromagnetic field simulation focuses on the accurate simulation of the interaction between equipment model and electromagnetic field distribution, it has high complexity and needs to consider both time cost and simulation accuracy. However, this paper mainly studies the improvement of simulation algorithm accuracy and ignores the time cost required in the actual simulation process. It pays more attention to theoretical research and lacks practical applicability. In addition, the main purpose of this paper is to analyze the electric field distribution simulation of GIS equipment under transient voltage. This simulation usually focuses on the response of equipment under voltage or current shock in a short period of time, and has weak generalization ability. In the patent application document with publication number CN107153721A, a symplectic finite-difference time-domain electromagnetic simulation method under a moving target is proposed. The Maxwell equations are constructed and the symplectic finite-difference time-domain method is used to solve the Maxwell equations. The algorithm solves the Maxwell equations and constructs the SFDTD method with 4th-order accuracy using the symplectic integral propagation factor. The essence is to approximate the time evolution matrix of the Maxwell equations. By adding the speed variable to the original Maxwell equations, more attention is paid to moving targets such as high-speed aircraft, missiles, and fast-moving cars. GIS electromagnetic field simulation focuses more on the accurate simulation of the interaction between the equipment model and the electromagnetic field distribution. GIS electromagnetic field simulation needs to consider the influence of the specific structure and material properties of the power equipment on the electromagnetic field, as well as how the electromagnetic field is distributed and changed inside these devices. In addition, GIS electromagnetic field simulation usually needs to consider the lossy and dispersion characteristics of the medium. In the actual electromagnetic field, the medium is often not an ideal perfect medium, but has a certain conductivity (causing loss) and frequency-related refractive index (causing dispersion). These characteristics have an important influence on the propagation and reflection of high-voltage and high-frequency electromagnetic fields. Summary of the invention

[0006] The technical problem to be solved by the present invention is how to efficiently obtain precise GIS electromagnetic field distribution simulation data.

[0007] The present invention solves the above technical problems by the following technical means:

[0008] The present invention proposes a GIS electromagnetic field distribution simulation method, which comprises:

[0009] Discretize the three-dimensional simulation model of GIS equipment into spatial grids, and determine the corresponding relationship between the number and spatial coordinate data of each grid node;

[0010] The Maxwell equations used to describe the electromagnetic field of GIS are converted into a fourth-order matrix, and the four-step HIE-FDTD algorithm is used to solve the fourth-order matrix and calculate the field strength data of each grid node;

[0011] Constructing graph structure based on spatial coordinate data and field strength data of grid nodes ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data;

[0012] The node interpolation enhancement strategy is used to generate new grid nodes and the eigenvalues ​​of the new grid nodes between adjacent grid nodes.

[0013] Furthermore, the three-dimensional simulation model of the GIS device is discretized into a spatial grid, and the corresponding relationship between the number and spatial coordinate data of each grid node is determined, including:

[0014] The three-dimensional simulation model of GIS equipment is discretized into rectangular space grids, and each grid node is numbered using the Yee grid concept;

[0015] Determine the corresponding relationship between the number of each grid node and its spatial coordinate data: ,in, It is along The discrete spatial step size in the axis direction, is the three-dimensional space coordinate of the grid node.

[0016] Furthermore, the Maxwell equations used to describe the electromagnetic field of the GIS are converted into a fourth-order matrix, and the fourth-order matrix is ​​solved using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node, including:

[0017] The Maxwell equations used to describe the electromagnetic field of GIS are converted into a sixth-order matrix form:

[0018]

[0019] In the formula, is the vector formed by the components of the electric field and magnetic field in the rectangular coordinate system, that is, T , is a sixth-order matrix;

[0020] Split the sixth-order matrix into four matrices:

[0021]

[0022] In the formula, for ,

[0023] for ;

[0024] The four-step HIE-FDTD algorithm is used to solve the fourth-order matrix form and calculate the field intensity data of each grid node.

[0025] Furthermore, the sixth-order matrix The formula is:

[0026]

[0027] In the formula, is the dielectric constant, is the magnetic permeability.

[0028] Furthermore, the four-step HIE-FDTD algorithm is used to solve the fourth-order matrix form and calculate the field strength data of each grid node, including:

[0029] The four-step HIE-FDTD algorithm is decomposed into , , , Four sub-steps, for each sub-step a semi-implicit difference scheme is used to calculate the three diagonal implicits of the field strength data about the grid nodes;

[0030] The pursuit method is used to solve the three diagonal implicit equations and obtain the grid node field strength data.

[0031] Furthermore, for The three diagonal implicits of the field strength data at the grid nodes are calculated step by step using a semi-implicit difference scheme, including:

[0032] Auxiliary vectors are introduced using a semi-implicit difference scheme , calculate the field strength data about the grid nodes The tridiagonal implicit form of :

[0033]

[0034] In the formula, is the point loss parameter of the medium, , , is the time step, Indicates that the electric field x , y , z The electric field component in the direction, Magnetic field x , y , z The magnetic field component in the direction, the superscript indicates the time step, for example Represents the time step n Ex, Represents the time step n and n The electric field E at the 1 / 4 position between +1 x .

[0035] Furthermore, the node interpolation enhancement strategy is used to generate new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes, and the formula is expressed as:

[0036]

[0037]

[0038]

[0039] In the formula, and are two adjacent grid nodes, is the new grid node, .

[0040] Furthermore, the field strength data of the new grid node is:

[0041]

[0042]

[0043] In the formula, and is the field strength data of two adjacent grid nodes, is the field strength data of the new grid node, is an undetermined coefficient, which is determined by λ in the above step and the field strength attribute relationship between the two nodes.

[0044] In addition, the present invention also proposes a GIS electromagnetic field distribution simulation system, the system comprising:

[0045] The grid discretization module is used to discretize the three-dimensional simulation model of the GIS equipment into a spatial grid and determine the corresponding relationship between the number and spatial coordinate data of each grid node;

[0046] The finite-difference time-domain simulation module is used to convert the Maxwell equations used to describe the electromagnetic field of GIS into a fourth-order matrix, and solve the fourth-order matrix using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node;

[0047] Graph structure representation module, used to build graph structure based on spatial coordinate data and field strength data of grid nodes ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data;

[0048] The node feature enhancement module is used to generate new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes by using the node interpolation enhancement strategy.

[0049] The advantages of the present invention are:

[0050] (1) The present invention adopts the finite difference time domain method in four steps to transform the Maxwell equations describing the GIS electromagnetic field into a fourth-order matrix to solve the problem, and then obtain the field strength data of each grid node. Then, a graph structure is constructed according to the spatial coordinate data and field strength data of the grid nodes. Based on the graph structure, a node difference enhancement strategy is adopted to interpolate between the samples of the minority class and the samples of their nearest neighbors to generate new samples, and finally obtain the simulation data of the fine grid, so as to quickly obtain the fine GIS electromagnetic field distribution simulation data.

[0051] 2) The present invention can significantly improve the accuracy and stability of the numerical solution by introducing a higher-order time integration scheme and a space integration scheme. The traditional FDTD method may be affected by numerical dissipation and numerical dispersion, and the present invention can reduce these problems through a more accurate integration method, which is particularly effective when dealing with long-term simulations or high-frequency electromagnetic wave propagation. For electromagnetic field simulations of complex structures, such as those with media, inhomogeneous media or inhomogeneous structures, this embodiment uses a four-step hybrid implicit explicit finite-difference time-domain (HIE-FDTD) method to better handle these complexities. Its accurate integration and numerical methods can effectively simulate electromagnetic field propagation and interactions between different media.

[0052] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] Figure 1 It is a flow chart of a GIS electromagnetic field distribution simulation method proposed in one embodiment of the present invention;

[0054] Figure 2 It is a principle block diagram of GIS electromagnetic field distribution simulation in one embodiment of the present invention;

[0055] Figure 3 It is a structural schematic diagram of a GIS electromagnetic field distribution simulation system proposed in one embodiment of the present invention. DETAILED DESCRIPTION

[0056] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in combination with the embodiments of the present invention. Obviously, the described embodiments are 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 creative work are within the scope of protection of the present invention.

[0057] like Figure 1 to Figure 2 As shown, the first embodiment of the present invention proposes a GIS electromagnetic field distribution simulation method, which includes the following steps:

[0058] S10, discretizing the three-dimensional simulation model of the GIS device into a spatial grid, and determining the corresponding relationship between the number and spatial coordinate data of each grid node;

[0059] It should be noted that, in this embodiment, a three-dimensional data model is established based on the GIS equipment and a simulation area of ​​the three-dimensional data model is discretized into a plurality of grid units.

[0060] S20, converting the Maxwell equations used to describe the electromagnetic field of GIS into a fourth-order matrix, and solving the fourth-order matrix using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node;

[0061] S30, constructing a graph structure based on the spatial coordinate data and field strength data of the grid nodes ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data;

[0062] S40, generating new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes using a node interpolation enhancement strategy.

[0063] This embodiment uses the time-domain finite difference method in four steps to transform the Maxwell equations describing the GIS electromagnetic field into a fourth-order matrix to solve the field strength data of each grid node, and then constructs a graph structure according to the spatial coordinate data and field strength data of the grid nodes. Based on the graph structure, a node difference enhancement strategy is used to interpolate between samples of the minority class and their nearest neighbor samples to generate new samples, and finally obtains fine grid simulation data, which can quickly obtain fine GIS electromagnetic field distribution simulation data.

[0064] As a further preferred technical solution, the step S10: discretizing the three-dimensional simulation model of the GIS equipment into a spatial grid and determining the corresponding relationship between the number and spatial coordinate data of each grid node includes the following steps:

[0065] S11, discretize the three-dimensional simulation model of GIS equipment into rectangular space grids, and use the Yee grid idea to number each grid node;

[0066] S12, determine the corresponding relationship between the number of each grid node and its spatial coordinate data: ,in, It is along The discrete spatial step size in the axis direction, is the three-dimensional space coordinate of the grid node.

[0067] As a further preferred technical solution, the step S20: converting the Maxwell equations used to describe the GIS electromagnetic field into a fourth-order matrix, and solving the fourth-order matrix using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node, includes the following steps:

[0068] S21. Convert the Maxwell equations used to describe the electromagnetic field of GIS into a sixth-order matrix form:

[0069]

[0070] In the formula, is the vector consisting of the components of the electric and magnetic fields in the rectangular coordinate system, ,in:

[0071]

[0072]

[0073]

[0074]

[0075]

[0076]

[0077] in, is a sixth-order matrix, expressed as:

[0078]

[0079] In the formula, is the dielectric constant, is the magnetic permeability.

[0080] It should be noted that this embodiment expresses the Maxwell equations in matrix form so that the electromagnetic field components are expressed in the form of vectors, and the differential operations in the equations are described by matrix operations. This form provides a more unified and compact description method, making the equations easier to understand and derive, and in the numerical solution process, numerical methods usually involve solving discretized equations, and the matrix form makes these solution methods more direct and effective.

[0081] S22. Convert the sixth-order matrix form into a fourth-order rectangular form:

[0082]

[0083] In the formula, and The representation forms are:

[0084]

[0085]

[0086] S23. The four-step HIE-FDTD algorithm is used to solve the fourth-order matrix form and calculate the field strength data of each grid node.

[0087] This embodiment can significantly improve the accuracy and stability of the numerical solution by introducing a higher-order time integration scheme and a space integration scheme. The traditional FDTD method may be affected by numerical dissipation and numerical dispersion, and this embodiment can reduce these problems through a more accurate integration method, which is particularly effective when dealing with long-term simulations or high-frequency electromagnetic wave propagation. For electromagnetic field simulations of complex structures, such as those with media, inhomogeneous media or inhomogeneous structures, this embodiment can better handle these complexities. Its accurate integration and numerical methods can effectively simulate electromagnetic field propagation and interaction between different media.

[0088] Specifically, the step S23: using a four-step HIE-FDTD algorithm to solve the fourth-order matrix form and calculate the field strength data of each grid node specifically includes:

[0089] The four-step HIE-FDTD algorithm is decomposed into , , , Four steps, including:

[0090] (1) :

[0091]

[0092] (2) :

[0093]

[0094] (3) :

[0095]

[0096] (4) Step by step:

[0097]

[0098] for Step by step: Using semi-implicit difference scheme to introduce auxiliary vectors , keeping the stability condition of the algorithm unchanged, we can get:

[0099]

[0100] in and is the mutual coupling equation, and the substitution method is used to eliminate The tridiagonal implicit equations are as follows:

[0101]

[0102] In the formula, is the point loss parameter of the medium, , , is the time step, represents the electric field components in the x, y, and z directions, The magnetic field components in the x, y, and z directions. The superscripts indicate the time step, e.g. Represents the time step n E x , represents the electric field Ex at the 1 / 4 position between time steps n and n+1, is the sixth-order identity matrix.

[0103] for , the auxiliary vector Substitute the expression into the expansion and get:

[0104]

[0105] for , the auxiliary vector Substitute the expression into the expansion and get:

[0106]

[0107] for , the auxiliary vector Substitute the expression into the expansion and get:

[0108]

[0109] According to the above equations, the components of the electric field and magnetic field in the x, y, and z directions at each step can be solved.

[0110] It should be noted that the present embodiment can obtain simulation calculation results through calculation, and the results are expressed as the spatial coordinates corresponding to each grid node and the distribution data of the electromagnetic field at that point. Since the coarseness of the simulation grid division is proportional to the calculation accuracy and inversely proportional to the calculation time, in order to solve the contradiction between accuracy and time, a real-time enhanced simulation method for grid data is designed, and grid enhancement is implemented for the coarse grid simulation electromagnetic field distribution results based on graph machine learning, and the result of finer grid simulation is also achieved. Furthermore, the step S40: using the node interpolation enhancement strategy to generate new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes is specifically used for:

[0111] Synthetic nodes are created by combining existing nodes through Node Interpolation, which is a generative enhancement strategy that generates new samples by interpolating between samples of the minority class and their nearest neighbors. The training distribution is as follows:

[0112]

[0113]

[0114]

[0115] In the formula, and are two adjacent grid nodes and both come from , is the new grid node, , that is, by the point and Point Generate a new node between , the new node position is given by The value of is determined.

[0116] As a further preferred technical solution, the field strength data of the new grid node is:

[0117]

[0118]

[0119] In the formula, and is the field strength data of two adjacent grid nodes, is the field strength data of the new grid node, is an undetermined coefficient. Δ is determined by λ and the field strength attribute relationship between the two nodes. The specific value needs to be confirmed in combination with the corresponding field strength, because the gradient of distance change is different from that of field strength change.

[0120] In this way, new nodes can be generated between adjacent nodes. The new node eigenvalues ​​also contain spatial data and field strength data, thereby achieving grid enhancement of the coarse grid simulation electromagnetic field distribution results and achieving the same results as the finer grid simulation results.

[0121] It should be noted that this embodiment involves spatial data and electromagnetic field data of existing nodes through an algorithm based on node proliferation between adjacent nodes. For grid data in GIS simulation, it has obvious topological structure characteristics. This embodiment can effectively capture the relationship between grids and retain the topological structure information of grid data. In addition, the algorithm of this embodiment generates the representation of new nodes through the characteristics of adjacent nodes, so that the algorithm can take into account the local characteristics of the nodes. In addition, this embodiment is an algorithm enhancement based on a graph structure for the coarse grid data obtained by simulation, which can adapt to grid data of different scales, including small-scale and large-scale grids, and is adaptable to GIS simulation data under different voltage levels of different devices, thereby realizing GIS simulation data enhancement and improving the quality of grid data.

[0122] In addition, if Figure 3 As shown, the second embodiment of the present invention proposes a GIS electromagnetic field distribution simulation system, characterized in that the system includes:

[0123] The grid discretization module 10 is used to discretize the three-dimensional simulation model of the GIS equipment into a spatial grid and determine the corresponding relationship between the number and spatial coordinate data of each grid node;

[0124] The finite difference time domain simulation module 20 is used to convert the Maxwell equations used to describe the electromagnetic field of the GIS into a fourth-order matrix, and solve the fourth-order matrix using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node;

[0125] A graph structure representation module 30 is used to construct a graph structure based on the spatial coordinate data and field strength data of the grid nodes. ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data;

[0126] The node feature enhancement module 40 is used to generate new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes by using a node interpolation enhancement strategy.

[0127] As a further preferred technical solution, the finite-difference time-domain simulation module 20 specifically includes:

[0128] The Maxwell equation conversion unit is used to convert the Maxwell equations used to describe the electromagnetic field of GIS into a sixth-order matrix form:

[0129]

[0130] In the formula, is the vector consisting of the components of the electric and magnetic fields in the rectangular coordinate system, is a sixth-order matrix;

[0131] The matrix conversion unit is used to convert the sixth-order matrix form into a fourth-order rectangular form:

[0132]

[0133] In the formula, for ,

[0134] for ;

[0135] The step-by-step solution unit is used to solve the fourth-order matrix form using the four-step HIE-FDTD algorithm and calculate the field strength data of each grid node.

[0136] As a further preferred technical solution, the node feature enhancement module 40 is specifically used to generate new grid nodes and feature values ​​of the new grid nodes between adjacent grid nodes using a node interpolation enhancement strategy, and the formula is expressed as:

[0137]

[0138]

[0139]

[0140] In the formula, and are two adjacent grid nodes, is the new grid node, .

[0141] It should be noted that other embodiments or implementation methods of the GIS electromagnetic field distribution simulation system described in the present invention can refer to the above-mentioned method embodiments and will not be described in detail here.

[0142] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0143] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0144] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A GIS electromagnetic field distribution simulation method, characterized in that: The method comprises: Discretize the three-dimensional simulation model of GIS equipment into spatial grids, and determine the corresponding relationship between the number and spatial coordinate data of each grid node; The Maxwell equations used to describe the electromagnetic field of GIS are converted into a fourth-order matrix, and the four-step HIE-FDTD algorithm is used to solve the fourth-order matrix to calculate the field strength data of each grid node, including converting the Maxwell equations used to describe the electromagnetic field of GIS into a sixth-order matrix form: In the formula, is the vector consisting of the components of the electric and magnetic fields in the rectangular coordinate system, is a sixth-order matrix; Convert the sixth-order matrix form to a fourth-order rectangular form: In the formula, and is a sixth-order matrix; The four-step HIE-FDTD algorithm is used to solve the fourth-order matrix form and calculate the field strength data of each grid node, including: the four-step HIE-FDTD algorithm is decomposed in the time domain into , , , Four steps, for each step, a semi-implicit difference format is used to calculate the three diagonal implicits about the grid node field strength data; the pursuit method is used to solve the three diagonal implicits to obtain the grid node field strength data, including the use of a semi-implicit difference format to introduce auxiliary vectors , calculate the field strength data about the grid nodes The tridiagonal implicit form of : In the formula, is the point loss parameter of the medium, , , is the time step, represents the electric field components in the x, y, and z directions, The magnetic field components in the x, y, and z directions. The superscripts represent the time steps. represents Ex at time step n, represents the electric field Ex at the 1 / 4 position between time steps n and n+1; Constructing graph structure based on spatial coordinate data and field strength data of grid nodes ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data; The node interpolation enhancement strategy is used to generate new grid nodes and the eigenvalues ​​of the new grid nodes between adjacent grid nodes. The formula is expressed as: In the formula, and are two adjacent grid nodes, is the new grid node, .

2. The GIS electromagnetic field distribution simulation method according to claim 1, characterized in that: The three-dimensional simulation model of the GIS device is discretized into a spatial grid, and the corresponding relationship between the number and the spatial coordinate data of each grid node is determined, including: The three-dimensional simulation model of GIS equipment is discretized into rectangular space grids, and each grid node is numbered using the Yee grid concept; Determine the corresponding relationship between the number of each grid node and its spatial coordinate data: ,in, It is along The discrete spatial step size in the axis direction, is the three-dimensional space coordinate of the grid node.

3. The GIS electromagnetic field distribution simulation method according to claim 1, characterized in that: The sixth-order matrix The formula is expressed as: In the formula, is the dielectric constant, is the magnetic permeability.

4. The GIS electromagnetic field distribution simulation method according to claim 1, characterized in that: The field strength data of the new grid node is: In the formula, and is the field strength data of two adjacent grid nodes, is the field strength data of the new grid node, is the coefficient to be determined.

5. A GIS electromagnetic field distribution simulation system, characterized in that: The system comprises: The grid discretization module is used to discretize the three-dimensional simulation model of the GIS equipment into a spatial grid and determine the corresponding relationship between the number and spatial coordinate data of each grid node; The finite-difference time-domain simulation module is used to convert the Maxwell equations used to describe the GIS electromagnetic field into a fourth-order matrix, and solve the fourth-order matrix using a four-step HIE-FDTD algorithm to calculate the field strength data of each grid node, including converting the Maxwell equations used to describe the GIS electromagnetic field into a sixth-order matrix form: In the formula, is the vector consisting of the components of the electric and magnetic fields in the rectangular coordinate system, is a sixth-order matrix; Convert the sixth-order matrix form to a fourth-order rectangular form: In the formula, and is a sixth-order matrix; The four-step HIE-FDTD algorithm is used to solve the fourth-order matrix form and calculate the field strength data of each grid node, including: the four-step HIE-FDTD algorithm is decomposed in the time domain into , , , Four steps, for each step, a semi-implicit difference format is used to calculate the three diagonal implicits about the grid node field strength data; the pursuit method is used to solve the three diagonal implicits to obtain the grid node field strength data, including the use of a semi-implicit difference format to introduce auxiliary vectors , calculate the field strength data about the grid nodes The tridiagonal implicit form of : In the formula, is the point loss parameter of the medium, , , is the time step, represents the electric field components in the x, y, and z directions, The magnetic field components in the x, y, and z directions. The superscripts represent the time steps. represents Ex at time step n, represents the electric field Ex at the 1 / 4 position between time steps n and n+1; Graph structure representation module, used to build graph structure based on spatial coordinate data and field strength data of grid nodes ,in, represents a collection of grid nodes, Represents an edge set, where the lines between grids are edges, and the characteristic values ​​of each grid node include spatial coordinate data and field strength data; The node feature enhancement module is used to generate new grid nodes and feature values ​​of new grid nodes between adjacent grid nodes using the node interpolation enhancement strategy. The formula is expressed as: In the formula, and are two adjacent grid nodes, is the new grid node, .

6. The GIS electromagnetic field distribution simulation system according to claim 5, characterized in that: The grid discretization module is specifically used for: The three-dimensional simulation model of GIS equipment is discretized into rectangular space grids, and each grid node is numbered using the Yee grid concept; Determine the corresponding relationship between the number of each grid node and its spatial coordinate data: ,in, It is along The discrete spatial step size in the axis direction, is the three-dimensional space coordinate of the grid node.

7. The GIS electromagnetic field distribution simulation system according to claim 5, characterized in that: The sixth-order matrix The formula is expressed as: In the formula, is the dielectric constant, is the magnetic permeability.

8. The GIS electromagnetic field distribution simulation system according to claim 5, characterized in that: The field strength data of the new grid node is: In the formula, and is the field strength data of two adjacent grid nodes, is the field strength data of the new grid node, is the coefficient to be determined.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 4 is implemented.

Images