A method for electromagnetic radiation suppression for data center flywheel energy storage systems

Abstract

The present application relates to flywheel energy storage electromagnetic suppression technology field, concretely for a kind of electromagnetic radiation suppression method for data center flywheel energy storage system, including the following steps, obtaining ground potential data generates response node set, comparing rotational speed and power event constructs synchronous topology, identifies coupling path and generates path atlas, adjusts magnetic shielding component, monitors magnetic flux change and generates electromagnetic compatibility state report.In the present application, by analyzing the mutation potential and recovery rate of ground node, screening the response node with disturbance characteristics, fusing the rotational speed and load power event, constructing the multi-source data response mapping in time domain, identifying the coupling path based on the physical connection relationship of node, collecting the magnetic flux density and magnetic hysteresis characteristics to identify the magnetic saturation trend, matching the excitation phase to implement dynamic adjustment, improving the response adaptability of shielding component, judging the stability by the change rate trend of magnetic flux density, controlling the reset operation of shielding structure, strengthening the shielding efficiency retention, realizing electromagnetic suppression under dynamic interference scene.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of flywheel energy storage electromagnetic suppression, and particularly relates to an electromagnetic radiation suppression method for a flywheel energy storage system of a data center. BACKGROUND

[0002] The technical field of flywheel energy storage electromagnetic suppression includes control and intervention means for the problem of electromagnetic radiation generated by the flywheel energy storage system during operation. The core content is to focus on the electromagnetic interference generated by the high-speed rotating components and their driving systems in the flywheel energy storage device during energy conversion and transmission to the external environment. The constraint of electromagnetic emission is realized through hardware layout, circuit optimization and shielding design, covering aspects such as energy storage system structure layout, electromagnetic compatibility design, electrical connection mode and grounding mode, focusing on solving the identification and control needs of electromagnetic radiation sources caused by the energy storage system under high-speed and high-frequency conditions, and ensuring its stable operation in complex power systems.

[0003] Among them, the electromagnetic radiation suppression method for the flywheel energy storage system of the data center refers to a method system for limiting or weakening the electromagnetic radiation interference generated by the flywheel energy storage device configured in the data center under the running state. It includes the identification of the radiation path in the flywheel motor driving circuit, the construction of the filter circuit of the power input and output port, the electromagnetic shielding structure design of the metal shell, and the adjustment and configuration of the system grounding mode. It mainly reduces the leakage of high-frequency interference signals by designing a low-pass filter network, constructs a closed magnetic field constraint area by deploying conductive shielding materials at key parts of the flywheel system, and forms a low-impedance leakage path by optimizing the electrical grounding form, thereby structurally suppressing electromagnetic radiation.

[0004] In the existing flywheel energy storage system electromagnetic radiation suppression process, the identification of electromagnetic interference sources mainly depends on the structure planning and preset parameters, and lacks real-time identification mechanism for dynamic disturbance behavior during operation, making it difficult to accurately track the actual propagation path of high-frequency disturbance. The synchronous comparison method based on event time sequence is not introduced in the processing process, which makes it impossible to establish an effective causal chain between different signal sources, and the identification result is prone to deviation or misjudgment. The topological connection relationship is not used as an auxiliary criterion to participate in path construction, and the actual formed shielding countermeasures lack structural logic constraints, resulting in the lack of precise spatial positioning ability in path construction. The dynamic adjustment mechanism based on magnetic state is not configured in the shielding structure design, which cannot identify the nonlinear response process of magnetic components caused by changes in excitation conditions during actual operation, and is prone to cause fluctuations in shielding effectiveness, especially in the unstable shielding area under high-speed and high-frequency power environment. The reset control mechanism cannot effectively evaluate the long-term trend of magnetic components, and the configuration adjustment stays in the static parameter setting stage, and the shielding ability presents a risk of degradation in long-term operation, resulting in weak recovery ability of the radiation path, which is prone to form new electromagnetic leakage paths under continuous or intermittent excitation. SUMMARY

[0005] The purpose of the present application is to solve the problems existing in the prior art, and a method for suppressing electromagnetic radiation of a flywheel energy storage system in a data center is proposed.

[0006] In order to achieve the above-mentioned purpose, the present application adopts the following technical scheme: a method for suppressing electromagnetic radiation of a flywheel energy storage system in a data center, comprising the following steps:

[0007] S1: Obtain the node grounding potential time series data in the flywheel energy storage system grounding network, perform potential difference division and mark the mutation points, screen the nodes whose potential change rate exceeds the surge threshold in a single sampling period, and restore to steady state in the next two periods, and generate a transient response node set combined with the time stamp;

[0008] S2: Statistically analyze the flywheel speed time series data and the load power time series data, detect the speed mutation event and the power transient event respectively, extract the event time, and compare the time window with the time stamp of the transient response node set, screen the nodes whose time difference is less than the time synchronization threshold, and obtain the synchronous response node topology;

[0009] S3: Based on the synchronous response node topology, retrieve the corresponding grounding module identifier, shielding connection point and rack grounding point, construct the electromagnetic coupling path, count the number of nodes in the path and trigger characteristics, and generate an electromagnetic radiation path map;

[0010] S4: Read the magnetic shielding component number in the electromagnetic radiation path map, collect the magnetic flux density data and the hysteresis loop characteristics, judge whether there is a magnetic saturation trend, screen the magnetic shielding components currently in the excitation phase, and perform magnetic circuit optimization adjustment, and generate a magnetic shielding optimization configuration table.

[0011] As a further scheme of the present application, the transient response node set includes mutation intensity label, steady state recovery identifier, time stamp feature, the synchronous response node topology includes event type corresponding relationship, time offset value between nodes, node physical layer structure, the electromagnetic radiation path map includes path node number, path excitation trigger feature, coupling module structure identifier, and the magnetic shielding optimization configuration table includes magnetic saturation prediction state, component excitation phase identifier, and optimization adjustment parameter.

[0012] As a further scheme of the present application, the specific acquisition step of the transient response node set is:

[0013] S111: Based on the grounding potential time series data of each node in the flywheel energy storage system grounding network, calculate the potential change rate of the node in the current sampling period, and perform mutation detection of the potential difference value, screen the node change rate value whose change rate exceeds the set surge threshold, and obtain the mutation potential change rate data;

[0014] S112: According to the mutation potential change rate data, according to the ground potential change rate sequence data of the corresponding node in the subsequent two periods, it is judged whether the potential is restored to the steady state interval, if the judgment standard of restoring the steady state is met, it is classified as a transient disturbance node, and the time stamp corresponding to the mutation moment is recorded, and a transient disturbance time tag set is obtained;

[0015] S113: According to the transient disturbance time tag set, the unique node identification information corresponding to each node at the mutation moment is extracted, and the disturbance response strength value of the node is calculated in combination with the node identification and time tag, and the time tag and node identification of the node whose disturbance response strength value is greater than the disturbance response threshold value are combined to generate a transient response node set.

[0016] As a further scheme of the application, the specific acquisition step of the synchronous response node topology is:

[0017] S211: Obtain the flywheel speed time series data and load power time series data in the flywheel energy storage system operation period, calculate the speed increment change value and increment change rate between adjacent sampling points, identify the speed mutation event and transient power mutation event, and record the corresponding time, generate a double-event time tag set;

[0018] S212: Based on the double-event time tag set, in combination with the time tag corresponding to each node in the transient response node set, determine the time difference between each event and node timestamp, set a time synchronization threshold, and record the nodes with a time difference absolute value less than or equal to the time synchronization threshold, to obtain a synchronous matching node index set;

[0019] S213: According to the synchronous matching node index set, identify and number each connection relationship between all marked nodes one by one, construct an adjacency matrix and mark whether there is a direct connection relationship between nodes, and establish a synchronous response node topology.

[0020] As a further scheme of the application, the specific acquisition step of the electromagnetic radiation path map is:

[0021] S311: Based on the physical connection relationship of each node in the synchronous response node topology, extract the ground module identification information, shielding connection point number and rack grounding point number corresponding to each node connection section, perform item-by-item filing and summarizing of the node corresponding grounding parameters in the path section, and establish a node grounding feature index set;

[0022] S312: According to the node grounding feature index set, perform structure recognition and trigger direction arrangement on the node sequence in each path, judge the electromagnetic coupling characteristics in the path, mark the shielding integrity in the path and record the structure information, and generate a coupling path feature set with the node sequence;

[0023] S313: According to the coupling path feature set, the node type of each node in the path is divided, the path propagation sequence identifier is established, the trigger starting order and the propagation chain direction of each node are identified, the number of nodes in each path and the proportion of each type of node are counted, the path structure difference is represented by combining the difference color line, and the electromagnetic radiation path atlas is generated.

[0024] As a further scheme of the present application, the specific acquisition step of the magnetic shielding optimization configuration table is:

[0025] S411: Read the magnetic shielding component number information marked in each path of the electromagnetic radiation path atlas, collect the magnetic flux density data, judge whether the saturation trend appears, select all magnetic shielding components that meet the saturation trend judgment and record the number, and acquire the magnetic saturation component number set;

[0026] S412: According to the magnetic saturation component number set, the current excitation state of each component is collected, whether each magnetic shielding component is in the rising edge time window is identified for the excitation waveform of each magnetic shielding component, the current magnetic flux density, loop reluctance, magnetic path length, coil turns and driving frequency are extracted, the magnetic circuit impact value of the magnetic shielding component is calculated, the magnetic shielding component whose magnetic circuit impact value is greater than the reference value is determined to establish the target list of magnetic permeance structure reconstruction, and the magnetic circuit adjustment structure set is established;

[0027] S413: Based on the magnetic circuit adjustment structure set, the corresponding magnetic permeance path of each magnetic shielding component to be adjusted is positioned in the structure diagram, the adjustment direction is judged by combining the magnetic circuit length and the magnetic permeability change trend, the component number, the adjustment direction, the target magnetic permeance value and the excitation frequency are recorded, and the magnetic shielding optimization configuration table is constructed.

[0028] As a further scheme of the present application, the method further comprises:

[0029] S5: The magnetic flux density change rate of each component in the magnetic shielding optimization configuration table is counted, the change amplitude in the moving average window is monitored, if the continuous detection period presents a convergence trend, a magnetic circuit reset instruction is triggered, the magnetic shielding component is adjusted to the initial configuration state, and an electromagnetic compatibility state report is generated;

[0030] The electromagnetic compatibility state report includes the magnetic flux density change trend, the shielding configuration stability index, and the reset adjustment trigger record.

[0031] As a further scheme of the present application, the specific acquisition step of the electromagnetic compatibility state report is:

[0032] S511: According to the magnetic shielding optimization configuration table recorded in all magnetic shielding component numbers, the magnetic flux density change rate between adjacent sampling periods is calculated, the absolute value of the continuous change rate of each magnetic shielding component is calculated, and it is judged whether they are all less than the preset change amplitude threshold value, then the corresponding component magnetic shielding is marked as tending to be stable state, and the magnetic flux convergence trend judgment result is obtained;

[0033] S512: Based on the magnetic shielding components identified in the magnetic flux convergence trend judgment result, a sliding average window is set and a magnetic flux density change value sequence is extracted in each window, the consistency and convergence characteristics of the change amplitude are evaluated, if the change trend is stable and the difference degree is within the range, the corresponding magnetic shielding component is judged to meet the magnetic circuit recovery condition, and is added to the reset execution candidate set according to the number, and the magnetic circuit reset target set is obtained;

[0034] S513: According to the magnetic circuit reset target set, the current magnetic permeance configuration of each magnetic shielding component is compared with the initial factory configuration parameters, the deviation evaluation of the current parameters and the initial values of each magnetic permeance path is performed, if the deviation range is within the specified threshold value, the magnetic circuit reset action can be triggered, the magnetic permeance path configuration is reset to the initial state, the state and path information of all restored components are summarized in order of component number, and an electromagnetic compatibility state report is generated.

[0035] Compared with the prior art, the advantages and positive effects of the present application are:

[0036] In the present application, by analyzing the mutation potential and recovery rate of the grounding node, the response node with disturbance characteristics is screened, the electromagnetic disturbance identification accuracy is improved, the speed and load power events are fused, the multi-source data response mapping in the time domain is constructed, the interference event correlation is enhanced, the coupling path is identified based on the node physical connection relationship, the spatial specificity of the propagation path construction is enhanced, the magnetic flux density and the magnetic hysteresis characteristics are collected to identify the magnetic saturation trend, the excitation phase is matched to implement dynamic adjustment, the shielding component response adaptability is improved, the stability is judged by the magnetic flux density change rate trend, the shielding structure reset operation is controlled, the shielding effectiveness retention is strengthened, the data-driven response regulation closed loop is constructed, and the electromagnetic suppression in the dynamic interference scene is realized. BRIEF DESCRIPTION OF DRAWINGS

[0037] Figure 1 The main step flowchart of the present application is shown in the figure;

[0038] Figure 2 The transient response node set acquisition flowchart of the present application is shown in the figure;

[0039] Figure 3 The synchronous response node topology acquisition flowchart of the present application is shown in the figure;

[0040] Figure 4 The electromagnetic radiation path map acquisition flowchart of the present application is shown in the figure;

[0041] Figure 5 A flowchart for obtaining an optimal configuration table of a magnetic shielding of the present application is shown in FIG. 1.

[0042] Figure 6 A flowchart for obtaining an electromagnetic compatibility state report of the present application is shown in FIG. 2. DETAILED DESCRIPTION

[0043] In order to make the purpose, technical solutions and advantages of the present application clearer, the present application is further described in detail below in combination with the drawings and examples. It should be understood that the specific examples described herein are only used to explain the present application and do not limit the present application.

[0044] In the description of the present application, it should be understood that the terms "length", "width", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer" and the like indicate the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as limiting the present application. In addition, in the description of the present application, the meaning of "a plurality of" is two or more, unless otherwise explicitly and specifically limited.

[0045] Please refer to Figure 1 An electromagnetic radiation suppression method for a data center flywheel energy storage system, comprising the following steps:

[0046] S1: Obtain the grounding potential time sequence data of each node in the grounding network of the flywheel energy storage system, perform potential difference on each node sequence and mark the mutation points, screen the nodes whose potential change rate exceeds the surge threshold (the transient overvoltage threshold set according to IEC61000-4-5 standard) in a single sampling period and returns to steady state in the subsequent two periods, and generate a transient response node set in combination with the time stamp;

[0047] S2: Statistically analyze the flywheel speed time sequence data and the load power time sequence data, detect the speed mutation event and the power transient event through the two groups of sequences respectively, extract the event time, and compare the time window with the time stamp of the transient response node set, screen the nodes with a time difference less than the time synchronization threshold (the standard time window (usually ≤10ms) for determining the synchronization event of the power system), and obtain the synchronization response node topology;

[0048] S3: Based on the physical connection relationship of the synchronization response node topology, retrieve the corresponding grounding module identifier, shielding connection point and rack grounding point, construct an electromagnetic coupling path (an electromagnetic radiation propagation model conforming to IEEE C95.3 standard), count the number of nodes and trigger characteristics in the path, and generate an electromagnetic radiation path map;

[0049] S4: Read the magnetic shielding component (magnetic conductive material component conforming to MIL-STD-461G standard) number in the electromagnetic radiation path map, collect magnetic flux density data and hysteresis loop characteristics (standard magnetization characteristic curve of ferromagnetic material), judge whether there is a magnetic saturation trend (magnetic flux density reaches more than 90% of the material saturation magnetization intensity), screen the magnetic shielding component currently in the excitation phase, and perform magnetic circuit optimization adjustment to generate a magnetic shielding optimization configuration table;

[0050] S5: Statistics the magnetic flux density change rate of each component in the magnetic shielding optimization configuration table, monitors the change amplitude in the moving average window (data smoothing window set according to EMI test standard CISPR16-1-1), if the continuous detection period presents a convergence trend (the change rate of three consecutive sampling periods decreases and the absolute value is <5%), the magnetic circuit reset instruction is triggered, the magnetic shielding component is adjusted to the initial configuration state, and the electromagnetic compatibility state report is generated.

[0051] The transient response node set includes mutation intensity label, steady-state recovery identification, timestamp feature, the synchronous response node topology includes event type corresponding relationship, node time offset value, node physical layer structure, the electromagnetic radiation path map includes path node number, path excitation trigger feature, coupling module structure identification, the magnetic shielding optimization configuration table includes magnetic saturation prediction state, component excitation phase identification, optimization adjustment parameter, the electromagnetic compatibility state report includes magnetic flux density change trend, shielding configuration stability index, reset adjustment trigger record.

[0052] Please refer to Figure 2 , the specific steps of S1 are:

[0053] S111: Based on the grounding potential time sequence data of each node in the flywheel energy storage system grounding network, calculate the potential change rate of the node in the current sampling period, and perform mutation detection of the potential difference value, screen the node change rate value that exceeds the set surge threshold, and obtain the mutation potential change rate data;

[0054] Based on the grounding potential time sequence data of each node in the flywheel energy storage system grounding network, the grounding potential value of each node at any two consecutive sampling time points is extracted, such as the first sampling point potential is 4.2V, the second sampling point potential is 5.5V, the potential change value is 1.3V, and the change rate is calculated as 130V / s by combining the sampling period Δt=0.01s. In the potential difference value mutation detection process, the absolute value of the difference between the current sampling time point potential and the previous sampling time point potential is judged. If the mutation change rate exceeds the threshold value, the node is recorded as a candidate mutation node. The threshold value is selected according to the transient overvoltage limit value set in IEC61000-4-5 standard, which is set to 850V / s. If the change rate exceeds this value, it is considered that the node has a mutation. For example: the potential of node A jumps from 2.1V to 10.7V in a certain period, and the sampling period is still 0.01s, so the change rate is 860V / s, which has exceeded the threshold value and is included in the candidate. At the same time, nodes with a change rate less than 600V / s are excluded, i.e. lower than the interference judgment standard. The change rate is in the 600-850V / s range, which is the interference judgment fuzzy area and is not recorded. According to this rule, after completing the node traversal, the change rate threshold node potential mutation rate data set is obtained, and the mutation change rate value is the most critical numerical result, representing the potential disturbance intensity experienced by the node, which is further used as the basis for subsequent identification.

[0055] S112: According to the mutation potential change rate data, according to the grounding potential change rate sequence data of the corresponding node in the subsequent two periods, it is judged whether the potential is restored to the steady state interval. If the judgment standard of restoring the steady state is met, it is classified as a transient disturbance node, and the time stamp corresponding to the mutation time is recorded to obtain the transient disturbance time label set.

[0056] Based on the data of the change rate of the abrupt potential, the trend is judged according to the sequence of the change rate of the ground potential in two consecutive cycles after the abrupt change of the node. The judgment criterion is whether the change rate of the potential decreases to the set steady-state range in two consecutive cycles. The steady-state recovery range is set as the change rate is less than 100V / s. If the change rate is 95V / s and 78V / s in the two cycles after the node abrupt change, it is considered to have recovered to the steady state, which meets the characteristics of transient disturbance, and the node is classified as a transient disturbance node. If the subsequent cycles still maintain a high change rate of potential, such as 340V / s and 400V / s, the node is removed. Whether the node meets the steady-state recovery is realized by comparing the change rate sequence. The timestamp of the node abrupt change is recorded with the sampling time t=0.25s. For example, node B, the abrupt change occurs at the 25th sampling point, and the corresponding timestamp is 25×0.01=0.25s. The node number and the abrupt change time are combined and recorded as the label "B-0.25". All nodes that meet the transient disturbance recovery criteria and their timestamps are sorted in turn to form a time-labeled set as the basic data source for response identification.

[0057] S113: Based on the transient disturbance time tag set, extract the unique node identifier information corresponding to each node at the time of the sudden change, and combine the node identifier with the time tag, using the formula:

[0058] ;

[0059] Calculate the first Disturbance response strength value of each node A transient response node set is generated by combining the time stamps and node identifiers of nodes whose disturbance response intensity values ​​are greater than the disturbance response threshold. Represents a node In the cycle The internal potential change value, Indicates the first The time sequence number of each cycle. Represents a node The standard deviation of the rate of change of potential between the two cycles before and after the mutation cycle. Represents a node In the The stable potential offset value of each non-mutation cycle. Represents a node The rate of change of the peak mutation, Represents a node The benchmark value for the steady-state rate of change of potential;

[0060] Based on the transient disturbance time label set, the node number and corresponding mutation time are extracted from each label item. A node-time mapping is established item by item to generate a transient disturbance node set. Then, the disturbance response intensity value is calculated for each node in the set. If the change values ​​of node C in the 1st to 5th periods are 1.2, 2.8, 0.5, 0.9, and 1.4V respectively, then the corresponding... Then the item is:

[0061] ;

[0062] Let the standard deviation of the potential change rate of this node during the period before and after the mutation be denoted. =12.5V / s, corresponding to the steady-state periodic offset potential value The voltages are 0.3V, 0.5V, and 0.7V respectively. Therefore, the sum of squares is 0.09 + 0.25 + 0.49 = 0.83. Taking the square root, we get:

[0063] ;

[0064] The peak change rate of node C is 870 V / s, the steady-state reference change rate is 90 V / s, and the absolute value of the difference is 780. The final calculation formula is:

[0065] ;

[0066] This value is compared with a disturbance response threshold of 0.035 (set based on statistical analysis of disturbance characteristics of multiple actual nodes in a flywheel energy storage system; 20 representative node samples were selected, and approximately 85% of the disturbance response intensity values ​​measured during their abrupt changes were within the range of 0.02 to 0.034; within this range, the node potential recovers to steady state after the abrupt change in no more than three cycles; when the intensity value exceeds 0.035, the disturbance exhibits characteristics of high instantaneous peak value, rapid potential decline, and significantly reduced subsequent fluctuations. Therefore, 0.035 is used as the distinguishing critical value. This value moderately increases with the increase of the peak value of the abrupt change potential rate and the sum of squares of the steady-state offset, possessing dynamic adaptability characteristics. It is suitable for actual working conditions with wide node distribution and complex electrical grounding grid structures, and can stably reflect the significant boundary between abrupt change response and steady-state recovery). If... If node C is included in the transient response node set, the result indicates that the mutation severity of node C meets the recording criteria and can be included in the analysis scope as a response node.

[0067] Table 1. Parameters of Nodes in Transient Perturbation Samples

[0068]

[0069] The key calculation parameters of node C are listed in Table 1, see Table 1, and the disturbance response intensity value is 0.0403, which meets the response node identification condition and is added to the transient response node set.

[0070] The disturbance response intensity value is a quantitative index for characterizing the degree of transient disturbance experienced by a node in the ground network of a flywheel energy storage system within a certain time period, which comprehensively reflects the cumulative degree of potential change rate of the node within multiple cycles before and after the mutation moment, the relative position of the mutation position in the time sequence, the degree of potential fluctuation before and after the mutation, and the degree of steady-state deviation, and is normalized by the difference between the disturbance peak value and the steady-state reference change rate. The specific embodiment is the aggregated performance of the node in the ground network, and the higher the value, the more dramatic and persistent the mutation of the node potential in a short period of time, the higher the transient response characteristics, and it is suitable for identifying and determining key disturbance sources, sensitive nodes and potential disturbance propagation paths. It is the core criterion for building a transient response node set.

[0071] Please refer to Figure 3 , the specific steps of S2 are:

[0072] S211: Obtain the flywheel speed time series data and load power time series data within the flywheel energy storage system operation cycle, calculate the speed increment change value and increment change rate between adjacent sampling points, identify the speed mutation event and transient power mutation event, and record the corresponding time, generate a double-event time tag set;

[0073] The flywheel speed time series data and the load power time series data in the flywheel energy storage system operation cycle are obtained. First, the sampling interval is uniformly standardized, and the sampling period is set to 10 ms, that is, 100 data points are recorded per second. For the flywheel speed sequence, the speed difference between two adjacent time points is extracted. For example, at t1=1.0 s and t2=1.01 s, the speed is 1500 rpm and 1810 rpm respectively, and the speed increment is 310 rpm, which exceeds the set threshold value of 300 rpm. Mark t2 as the speed mutation time. The threshold value refers to the typical characteristics of the dynamic response of the speed change in the high-inertia system, and is determined to be reasonable within the IEC60034 standard range. The load power sequence is processed in the same way. The power difference between each two adjacent time points is calculated. If the power increases from 5.2 kW to 7.0 kW at t3=2.04 s, the change is 1.8 kW, which exceeds the set power mutation threshold value of 1.5 kW, then t3 is recorded as the power mutation time. The threshold value is set according to the typical instantaneous rising power characteristics in the power load disturbance experiment to avoid misjudgment of power fluctuation. All time points that meet the conditions are marked with event type and stored in the event sequence. The mark “TS” represents speed mutation, and “PL” represents power disturbance. For example, “TS-1.01” and “PL-2.04”. Through this way, the time label extraction is completed, and a unified data set is arranged to form a double-event time label set.

[0074] S212: Based on the double-event time label set, the time difference between each event and the node timestamp is determined by combining the time label corresponding to each node in the transient response node set. A time synchronization threshold is set, and nodes with a time difference absolute value less than or equal to the time synchronization threshold are recorded to obtain a synchronization matching node index set.

[0075] According to the double-event time label set, each mutation event time point is traversed, and the time label corresponding to the node in the transient response node set obtained in the previous stage is compared. The difference between the node mutation time and the event time is extracted. If node R1 mutates at 1.013 seconds and event TS-1.01, the time difference between the two is 3 ms, and the absolute value is less than the set threshold value of 10 ms, which meets the power system synchronization event determination condition. The threshold value is derived from the maximum allowed drift standard of the synchronization clock error (IEEE1588). The judgment standard is whether it meets |t n -t e ≤0.01 s. After traversing all nodes and events, record the node numbers with a time difference within 10 ms, such as nodes R1, R3, and R6. Arrange to form an array [1, 3, 6], and mark it as a synchronization matching node index set. Exclude the synchronization failed nodes and mark the event association of each index, such as R1→TS and R3→PL, to ensure the unique mapping between nodes and events, and obtain the synchronization matching node index set.

[0076] S213: According to the synchronization matching node index set, the connection relationship between all marked nodes is identified and numbered one by one, the adjacency matrix is constructed, and whether there is a direct connection relationship between the nodes is marked. The nodes that have met the time synchronization condition are marked as synchronization state 1, and the nodes that are not involved are marked as state 0, and the synchronization response node topology is established;

[0077] According to the synchronization matching node index set, the system connection relationship parameters of each node are extracted, including bus number, feeder serial number, upstream and downstream node number and connection path identification, the network structure of the node serial number is found in combination with the preset power grid topology table, the adjacent node relationship matrix is arranged, for example, node 1 connects node 2, 5, node 3 connects node 4, 6, all connection relationships are uniformly marked, and the node state identification of the nodes in synchronization state is set to 1, and the others are set to 0. Finally, a topology table containing structure number, connection path and synchronization state is constructed, and a synchronization response node topology is generated.

[0078] Table 2 synchronization matching node structure and topology table

[0079]

[0080] Table 2 lists the structure topology parameter information of each node after time synchronization screening, referring to Table 2, the network propagation path in the synchronization response can be judged according to the node position and connection path, and the construction process of the synchronization response topology is supported.

[0081] Please refer to Figure 4 , the specific steps of S3 are:

[0082] S311: Based on the physical connection relationship of each node in the synchronization response node topology, the grounding module identification information, shielding connection point number and rack grounding point number corresponding to each node connection section are extracted, and the corresponding grounding parameters of the nodes in the path section are itemized and summarized, and a node grounding feature index set is established;

[0083] Based on the physical connection relationship of the synchronous response node topology, the connection path composed of nodes is deconstructed and the identification information is extracted. In the execution process, the topology connection matrix obtained in the previous step is first called, and the node pairs in the connection path are searched one by one. For example, the connection between nodes N1 to N4 is a path segment. For this path segment, the electrical structure port connected to the starting point and the ending point is extracted, and the connected ground module code is obtained through the association of the structure port, for example, N1 connects GND-01 module, which represents access to the bus bar, and N4 connects GND-02 module, which represents access to the shielded shell. The ground module identification needs to be confirmed according to the plant electrical installation drawing; Next, the interface parameters corresponding to the cables at both ends of the path are compared, such as the N1 connection point is a shielded joint with a model of SPX-M16, which is determined to be a shielded point with a number of SP-003 through the number rule, and the positioning coordinate value on the structure layout is extracted (such as N1 is installed on the third layer of rack K7, with a positioning point of X=1.2 m, Y=0.3 m) through the cabinet and rack structure diagram in the equipment installation drawing. The corresponding rack grounding point number J-K7-3 is obtained by matching the grounding reference point list. After all the parameters are extracted, the grounding parameter index entries including path number, starting and ending nodes, ground module number, shielded connection point number and rack grounding point number are merged to form. In order to improve the clarity of the structure, each path needs to be recorded separately and numbered to form a table, such as path segment number "P-001", associated nodes "N1-N4", corresponding ground modules "GND-01 to GND-02", shielded number "SP-003", grounding point "J-K7-3". Such processing is iterated through the entire topology network, and finally the node grounding feature index set is obtained.

[0084] S312: According to the node grounding feature index set, the node sequence in each path is structurally identified and the direction is arranged, the electromagnetic coupling characteristics in the path are judged, the shielding integrity in the path is marked and the structure information is recorded, and the coupling path feature set is generated according to the node sequence;

[0085] Based on the node grounding feature index set, the node pairs in each connection path are compared in structure and analyzed in relation to grounding. First, according to the known grounding module type, it is determined whether there is a cross-module connection in the path. If the starting grounding module of the path is GND-01, the ending is GND-03, and there is a relay node in the middle, it is marked as a multi-module path segment. The difference in the position of the grounding point in this structure may have a significant impact on the coupled path, and needs to be identified first. Then the node spacing is measured. Relying on the physical coordinates of the nodes in the process layout, the Euclidean distance between each two adjacent nodes is calculated. If the path "P-002" is composed of nodes N5→N6→N8, the node spacing is 0.6 meters and 0.9 meters respectively, and the total length of the path is 1.5 meters. Record it as the length parameter of this path segment. If the distance between the nodes exceeds 300 millimeters, mark it as a high interval path segment, and further record the shielding connection point number sequence and its corresponding cable number. For example, N5 uses shielding line SP-009, N6 uses SP-010, and the structural balance degree parameter of the path is established item by item. For each type of structure index, such as path length, grounding point spacing, shielding number difference, etc. Establish a basic feature parameter table to evaluate the structure level of the path in the overall topology. For example, path "P-002" is a connection within the same rack, and contains multiple grounding modules and multiple shielding points. Define its structure level as "complex structure" and record it in the feature field. Finally, a complete record table of structure type, length distribution, grounding point density, shielding sequence, etc. is formed as a multi-dimensional structure archive of the path, and a coupled path feature set is generated.

[0086] S313: According to the coupled path feature set, divide the nodes in each path into node types, establish path propagation sequence identification, identify the trigger starting order and propagation chain direction of each node, count the number of nodes in each path and the ratio of each type of node, represent the path structure difference with different color lines, and generate an electromagnetic radiation path map;

[0087] According to the coupling path feature set, the node number sequence involved in the path segment is traversed one by one, the path source node is determined according to the trigger timestamp of the starting node, and the propagation chain is constructed in ascending order of node number in the topology. If the node sequence in the path segment "P-003" is N10→N11→N13→N15, and the trigger order is N10=0.04s, N11=0.06s, N13=0.08s, N15=0.1s, then it is defined as the source, relay, relay, and terminal node in turn, and the propagation direction is indicated by an arrow on the path structure diagram, from left to right, indicating the starting to terminal path direction. The total number of nodes in the path structure is counted, the path contains 4 nodes, the structure sequence is complete and uninterrupted, and is marked as "continuous path". If there is a jump node (such as N10→N13) in the path, the timestamp of the intermediate node is missing, and it is recorded as "interrupted path". Then the distribution of the grounding module of each node is counted. If multiple nodes share a grounding module (such as N10 and N11 are both GND-01), the grounding sharing rate of the path is recorded as 50%. Combined with the node spacing and grounding point number, the propagation order table is arranged, and is summarized according to the path segment number, and is distinguished by color according to the path type, and the trigger direction is marked by a line. Finally, a complete path diagram is drawn. The path diagram should include node number, propagation direction, node trigger time, structure classification and other information, as an integrated presentation form of topology structure and electromagnetic propagation characteristics, and an electromagnetic radiation path atlas is generated.

[0088] Table 3 Grounding path structure parameter table

[0089]

[0090] Table 3 lists the grounding module combination, path length, node number and sharing rate data of typical path segments. Referring to Table 3, the differences in structure parameters of different paths can be clearly determined to support subsequent atlas generation.

[0091] Referring to Figure 5 , the specific steps of S4 are as follows:

[0092] S411: Read the magnetic shielding component number information marked in each path of the electromagnetic radiation path atlas, collect the magnetic flux density data, judge whether the saturation trend appears, select all magnetic shielding components that meet the saturation trend judgment and record the number, and obtain the magnetic saturation component number set;

[0093] The electromagnetic radiation path map is read to extract the magnetic shielding component number information. A list of components conforming to the magnetic permeable material requirements of the MIL-STD-461G standard is extracted. During the filtering operation, all component numbers with the prefix "M" are extracted based on the map identifier, such as M001, M002, etc. Material verification is performed on each component number, calling the material field from the component technical specifications table to determine if it is a standard magnetic permeable material such as iron-based alloy, permalloy, or pure nickel alloy. Numbers meeting the standard are retained and recorded in the list of valid components. Subsequently, magnetic flux density is acquired through an embedded magnetic flux sensing module. Measurement points are set up in the core and edge areas of the component for magnetic flux acquisition. The sampling time is set to 2ms. For example, the magnetic flux of component M001 is 0.79T at t=2.00s, and at t=2... At 01s, the magnetic flux is 0.81T, and the average value is 0.80T. This value is compared with its saturation magnetic flux density. If the component material is 1J85 alloy, its saturation magnetic flux density is 0.89T. Therefore, the current magnetic flux is 89.9% of the saturation magnetization. Combined with the magnetic saturation judgment benchmark value of 0.9×Bs=0.801T, since 0.80T is lower than the judgment threshold, the saturation mark is not triggered. Next, component M002 is measured, and the magnetic flux value is read as 0.84T, while its saturation magnetic flux density is 0.88T, corresponding to 90.9%, which is greater than the saturation judgment threshold of 0.792T. M002 is recorded as a magnetically saturated component. This type of data is archived and organized according to component number, magnetic flux value, and judgment result, and finally forms a set of component numbers containing all components that have reached the magnetic saturation trend condition, which is the magnetically saturated component number set.

[0094] S412: Based on the magnetic saturation component number set, collect the current excitation state of each component. For the excitation waveform of each magnetic shielding component, identify whether it is within the rising edge time window, and extract the current magnetic flux density, loop reluctance, magnetic path length, coil turns, and driving frequency using the formula:

[0095] ;

[0096] Calculate the first Magnetic circuit impact value of each magnetic shielding component To identify magnetically shielded components with magnetic circuit impact values ​​exceeding a reference value, a target list of components requiring magnetic permeability structure reconstruction is established, and a set of magnetic circuit adjustment structures is created. This indicates the current magnetic flux density of the component. The saturation magnetic flux density of its magnetically conductive material. This indicates the current driving frequency of the component, in Hz. This represents the resistance of the excitation circuit, in ohms. Indicates the number of turns in the excitation winding. Indicates the length of the magnetic permeation path, in meters;

[0097] According to the magnetic saturation component number set, the current state of each component is extracted whether it is in the excitation phase. The excitation state judgment is based on the driving signal frequency and the current waveform phase. The excitation frequency range is set to 1 kHz to 20 kHz, and the sampling frequency is 50 kHz. The power-on period of each component is analyzed. If the component M005 appears continuous current rise and synchronous detection voltage rise segment in the interval t=0.300s~0.301s, it is judged to be in the excitation rising edge, and the corresponding mark is the excitation state, and is recorded into the excitation component set. Extract its magnetic flux density value, winding turns, path length, resistance and frequency parameters. Take M005 as an example. Its magnetic flux density is 0.83T, saturated magnetic flux is 0.91T, frequency is 15kHz, loop resistance is 2.1Ω, coil turns is 20 turns, and magnetic path length is 0.9m. The above values are brought into the formula, and the operation process is as follows:

[0098] ;

[0099] The result Zg is 0.0367, which is lower than the reference value 0.05 (the setting basis is to compare the numerical distribution characteristics of the magnetic flux excess increase term and the magnetic circuit structure impedance term. When the magnetic flux density reaches more than 90% of the saturation strength, the frequency is more than 10kHz, the loop resistance is lower than 3Ω, and the winding turns are not higher than 25, the magnetic circuit impact value is concentrated between 0.04 and 0.06. Among them, those higher than 0.05 are generally accompanied by abnormal fluctuation of excitation end permeability or non-uniformity of magnetic field distribution. Therefore, the value 0.05 is taken as the critical reaction mark value of the magnetic circuit structure. This value increases with the increase of frequency and the decrease of winding density, and remains stable when the path length is less than 0.8m), so this component does not need to adjust the permeability structure. If the result of a component such as M007 is 0.0643, which exceeds the reference value, it is included in the adjustment list. After traversing all components, the adjustment number and target direction structure set are formed, which are the optimization execution objects of the magnetic circuit adjustment structure set.

[0100] The magnetic circuit impact value is used to measure the deviation between the electromagnetic response strength of the magnetic shielding component in the excitation state and the bearing limit of the magnetic circuit structure. Specifically, it reflects the transient magnetic response load of the component under the combined action of magnetic flux density, excitation frequency, resistance limit, coil structure and magnetic path length when it approaches the magnetic saturation boundary. The larger the value, the more likely the component is to produce nonlinear magnetization behavior or approach saturation state due to excitation coupling, resulting in a decrease in magnetic shielding effectiveness or abnormal accumulation of magnetic energy. Therefore, the magnetic circuit impact value not only comprehensively depicts the coordination matching degree of the current magnetic state and structure configuration of the component, but also provides a quantitative basis for the optimization adjustment of the magnetic path.

[0101] The formula is based on the square term of the magnetic flux density over-limit Reflect the difference between the current component and its saturation critical point, and multiply the difference with the excitation frequency The denominator is the loop resistance The numerator is the number of turns of the winding The electromagnetic impedance factor is formed, which represents the restriction effect of the loop on the flow of excitation current. The larger the impedance, the weaker the response. Therefore, the item is used as the divisor in the formula to attenuate the amplitude, forming a modulation structure inversely related to the magnetic saturation trend. Multiply the right side by The item combines the adjustment influence of the magnetic path length and the winding structure density on the entire magnetic field distribution path, where represents the length of the magnetic circuit body, and is the reciprocal of the coil density, which is used to weaken the adjustment effect of high-turn windings on length. The overall formula logic is as follows: the magnetic saturation critical difference is amplified by the excitation frequency, then weakened by the electromagnetic impedance, and then directionally corrected according to the structural path changes, so as to fully express the magnetic circuit impact strength of the magnetic shielding component under excitation in the case of approaching saturation.

[0102] S413: Based on the magnetic circuit adjustment structure set, locate the corresponding magnetic permeance path in the structure diagram for each magnetic shielding component to be adjusted, judge the adjustment direction based on the magnetic circuit length and the magnetic permeability change trend, record the component number, adjustment direction, target magnetic permeability value, and excitation frequency, and build a magnetic shielding optimization configuration table.

[0103] Based on the magnetic circuit adjustment structure set, the magnetic circuit optimization reconstruction is performed for the identified magnetic shielding component number to be adjusted. First, locate the magnetic permeance path section corresponding to the component in the structure diagram, record the start node and end node number, such as M007 connection path segment N14→N16→N19, total path length 1.2m, use equal interval node configuration, each segment interval 0.6m; judge its original magnetic permeability structure as single path configuration, magnetic permeability initial value μ is 4200H / m, combine the magnetic permeability change rule and the component current excitation state to judge the direction, if the magnetic permeability distribution along N14→N19 direction is 4200→3900→3600H / m, then set the adjustment direction as positive direction, the target value is set as increase amplitude 5%, record as target μ respectively 4410, 4095, 3780H / m; the adjusted magnetic circuit configuration needs to replace the magnetic permeability of each path segment and record the corresponding segment number, magnetic permeability before and after adjustment, target change amplitude and excitation frequency in the configuration table. Finally, build the table configuration structure with component number, path direction, target magnetic permeability value and corresponding frequency, and generate the magnetic shielding optimization configuration table.

[0104] Table 4 Magnetic shielding component optimization configuration table

[0105]

[0106] Table 4 lists the magnetic permeability settings and adjustment ratios before and after the optimization of the magnetic guide path structure of a typical excitation assembly. Referring to Table 4, this table serves as a configuration reference source for the execution of the magnetic guide path adjustment task.

[0107] Referring to Figure 6 The specific steps of S5 are as follows:

[0108] S511: According to the numbers of all magnetic shielding assemblies recorded in the magnetic shielding optimization configuration table, the rate of change of magnetic flux density between adjacent sampling periods is calculated, the absolute value of the continuous change rate of each magnetic shielding assembly is calculated, and it is judged whether they are all less than the preset change amplitude threshold. If so, the corresponding assembly magnetic shielding is marked as tending to a stable state, and the magnetic flux convergence trend judgment result is obtained.

[0109] The rate of change of magnetic flux density of the magnetic shielding assemblies listed in the magnetic shielding optimization configuration table is calculated. For each assembly, data is acquired for three consecutive sampling periods, the sampling period is 10 ms, time series data is constructed, and magnetic flux density difference value operation is performed point by point. The change rate is obtained by dividing the change between adjacent points by the magnetic flux value of the previous period. The magnetic flux densities of assembly M012 at t=0.100s, 0.110s, and 0.120s are 0.81T, 0.79T, and 0.77T, respectively, and the change rates are (0.79-0.81) / 0.81=-2.47%, (0.77-0.79) / 0.79=-2.53%, both are negative and the absolute values are less than 5%, the direction is consistent and the amplitude is in the convergence judgment interval. Then detect assembly M015, the sampling values are 0.84T, 0.82T, and 0.83T, and the corresponding change rates are -2.38% and -1.22%. The third segment change is positive and does not meet the continuous decrease condition, so it is not included. In this way, all assemblies in the magnetic shielding optimization configuration table are processed one by one. The numbers of assemblies that meet the three consecutive decrease conditions and have a change rate absolute value less than 5% are recorded, a magnetic flux convergence trend sequence is generated, and a data list is established according to the numbers. Finally, 8 numbers that meet the trend standard, such as M012, M018, and M021, are extracted from all assemblies and recorded, and a magnetic flux convergence trend judgment set is obtained.

[0110] S512: Based on the magnetic shielding assemblies identified in the magnetic flux convergence trend judgment result, a sliding average window is set and the magnetic flux density change value sequence is extracted in each window. The consistency and convergence characteristics of the change amplitude are evaluated. If the change trend is stable and the difference is within the range, the corresponding magnetic shielding assembly is determined to meet the magnetic circuit recovery condition, and is added to the reset execution candidate set according to the number. The magnetic circuit reset target set is obtained.

[0111] Based on the magnetic flux convergence trend, the component number in the set is determined, the change amplitude evaluation in the moving average window is performed, the window width is set to three sampling periods, the magnetic flux density sequence of M012 at t=0.100s to t=0.120s is [0.81T, 0.79T, 0.77T], the average value is 0.79T, the difference value of each segment in the window is normalized, and combined with the frequency, resistance, and permeability parameters to construct the convergence trend evaluation item. The frequency of M012 is 15kHz, the resistance is 2.1Ω, and the permeability is 4100H / m. Compared with the standard frequency of 10kHz and the standard permeability of 4000H / m, the change amplitude is within the preset deviation tolerance, which is marked as stable change. The window of M018 is [0.84T, 0.82T, 0.80T], the average value is 0.82T, the frequency is 16kHz, the resistance is 1.9Ω, and the permeability is 4300H / m. It is judged that the amplitude is also in the controllable interval, and is also included. If it is detected that M019 appears in the sequence [0.87T, 0.84T, 0.81T] with a change of more than 4%, it will be excluded. Finally, the component set that meets the reset judgment condition is obtained, including numbers M012, M018, M021, M022, M025, a total of 5, and the magnetic circuit reset target set is established.

[0112] S513: According to the magnetic circuit reset target set, compare the current magnetic shielding component permeability configuration with the initial factory configuration parameters, perform deviation evaluation on the current parameters and initial values of each magnetic permeability path, and if the deviation range is within the specified threshold, the magnetic circuit reset action can be triggered, the magnetic permeability path configuration is reset to the initial state, and the state and path information of all restored components are summarized in order of component number to generate an electromagnetic compatibility state report.

[0113] According to the magnetic circuit reset target set, components M012 and M018 are taken as examples, and the permeability of the current magnetic shielding configuration and the initial configuration at the factory are compared. The current permeability of component M012 path P23 is 4050H / m, the initial value is 4000H / m, the deviation is +1.25%, which is less than the maximum allowed offset of 5%, and it is determined that it can be directly reset. The corresponding magnetic permeability path length is 1.1m, and the number of turns is 25 turns. The current configuration of component M018 is 3850H / m, the initial value is 3800H / m, the deviation is +1.32%, the resistance is 2.3Ω, the frequency is 14.5kHz, and it is also within the allowed range, confirming the reset. After each reset, the component number, path number, original and current permeability, change amplitude, and state identification are recorded, and finally a table form structured data is generated, forming a unified record format, and an electromagnetic compatibility state report is generated.

[0114] Table 5 Magnetic circuit reset component state record table

[0115]

[0116] As shown in Table 5, the reset components M012 and M018 have a recovery of magnetic permeability in the adjustment range, and the table content is summarized according to the reset state data. See Table 5 for complete reset task archive information.

[0117] The above is only a preferred embodiment of the present application, and does not limit the present application in other forms. Any skilled person in the art can use the disclosed technical content to make changes or modifications as equivalent embodiments applied to other fields, but any simple modification, equivalent change and modification made to the above embodiments without departing from the technical solution content of the present application, according to the technical essence of the present application, still belongs to the protection scope of the technical solution of the present application.

Claims

1. A method of electromagnetic radiation suppression for a data center flywheel energy storage system, comprising: The method comprises the following steps: S1: acquiring time sequence data of node grounding potential in a flywheel energy storage system grounding network, performing potential difference division, marking mutation points, screening nodes with potential change rate exceeding a surge threshold in a single sampling period and returning to a steady state in the next two periods, and generating a transient response node set in combination with time stamps; S2: statistically analyzing flywheel speed time sequence data and load power time sequence data, detecting speed mutation events and power transient events respectively, extracting event time, and performing time window comparison with the time stamp of the transient response node set to screen nodes with time difference less than a time synchronization threshold to obtain a synchronous response node topology; S3: based on the synchronous response node topology, searching for corresponding grounding module identifiers, shielding connection points, and rack grounding points, constructing an electromagnetic coupling path, statistically analyzing the number of nodes in the path and triggering characteristics, and generating an electromagnetic radiation path map; S4: reading the magnetic shielding component number in the electromagnetic radiation path map, collecting magnetic flux density data and hysteresis loop characteristics, judging whether there is a magnetic saturation trend, screening the magnetic shielding components currently in the excitation phase, and performing magnetic circuit optimization adjustment to generate a magnetic shielding optimization configuration table.

2. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The transient response node set comprises mutation intensity labels, steady-state recovery identifiers, and time stamp features. The synchronous response node topology comprises event type correspondence, time offset values between nodes, and node physical layer structures. The electromagnetic radiation path map comprises path node numbers, path excitation triggering characteristics, and coupling module structure identifiers. The magnetic shielding optimization configuration table comprises magnetic saturation prediction states, component excitation phase identifiers, and optimization adjustment parameters.

3. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The specific acquisition steps of the transient response node set are as follows: S111: based on the grounding potential time sequence data of each node in the flywheel energy storage system grounding network, calculating the potential change rate of the node in the current sampling period, and performing mutation detection of the potential difference value, screening the node change rate value exceeding the set surge threshold, and acquiring the mutation potential change rate data; S112: according to the mutation potential change rate data, according to the grounding potential change rate sequence data of the corresponding node in the next two periods, judging whether the potential returns to the steady state interval, if the judgment standard of returning to the steady state is met, it is classified as a transient disturbance node, and the time stamp corresponding to the mutation time is recorded, and a transient disturbance time label set is obtained; S113: according to the transient disturbance time label set, extracting the unique node identifier information of each node at the mutation time, combining the node identifier and the time label, calculating the disturbance response intensity value of the node, and combining the time label and the node identifier of the node with the disturbance response intensity value greater than the disturbance response threshold value to generate a transient response node set.

4. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The specific acquisition steps of the synchronous response node topology are as follows: S211: acquiring flywheel speed time sequence data and load power time sequence data in a flywheel energy storage system operation period, calculating the speed increment change value and the increment change rate between adjacent sampling points, identifying speed mutation events and transient power mutation events, and recording the corresponding time to generate a double-event time label set; S212: Based on the double event time label set, the time difference between each event and the node timestamp is determined by combining the time label corresponding to each node in the transient response node set, a time synchronization threshold is set, and nodes with a time difference absolute value less than or equal to the time synchronization threshold are recorded to obtain a synchronization matching node index set; S213: According to the synchronization matching node index set, the connection relationship between all marked nodes is identified and numbered one by one to construct an adjacency matrix and mark whether there is a direct connection relationship between nodes, and a synchronous response node topology is established.

5. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The specific acquisition steps of the electromagnetic radiation path atlas are: S311: Based on the physical connection relationship of each node in the synchronous response node topology, the ground module identification information, shielding connection point number and rack grounding point number corresponding to each node connection section are extracted, and the node corresponding ground parameters in the path section are archived and summarized one by one to establish a node ground feature index set; S312: According to the node ground feature index set, the structure of each path node sequence is identified and the trigger direction is arranged, the electromagnetic coupling characteristics in the path are judged, the shielding integrity in the path is marked and the structure information is recorded, and the coupling path feature set is generated by the node sequence; S313: According to the coupling path feature set, the node type of each node in the path is divided, the path propagation sequence is established, the trigger starting order and propagation chain direction of each node are identified, the number of nodes in each path and the ratio of each type of node are counted, the path structure difference is represented by combining the difference color line, and the electromagnetic radiation path atlas is generated.

6. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The specific acquisition steps of the magnetic shielding optimization configuration table are: S411: Read the magnetic shielding component number information marked in each path of the electromagnetic radiation path atlas, collect the magnetic flux density data, judge whether there is a saturation trend, select all magnetic shielding components that satisfy the saturation trend judgment and record the number, and obtain a magnetic saturation component number set; S412: According to the magnetic saturation component number set, the current excitation state of each component is collected, whether each magnetic shielding component is in the rising edge time window is identified for the excitation waveform of each magnetic shielding component, the current magnetic flux density, loop inductance, magnetic path length, coil turns and driving frequency are extracted, the magnetic circuit impact value of the magnetic shielding component is calculated, the magnetic shielding component with a magnetic circuit impact value greater than a reference value is determined to establish a target list of magnetic guide structure reconstruction, and a magnetic circuit adjustment structure set is established; S413: Based on the magnetic circuit adjustment structure set, the corresponding magnetic guide path of each magnetic shielding component to be adjusted is positioned in the structure diagram, the adjustment direction is judged by combining the magnetic path length and the magnetic permeability change trend, the component number, adjustment direction, target magnetic guide value and excitation frequency are recorded, and the magnetic shielding optimization configuration table is constructed.

7. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 1, wherein, The method further comprises: S5: The magnetic flux density change rate of each component in the magnetic shielding optimization configuration table is counted, the change amplitude in the moving average window is monitored, if the continuous detection period presents a convergence trend, a magnetic circuit reset instruction is triggered, the magnetic shielding component is adjusted to the initial configuration state, and an electromagnetic compatibility state report is generated; The electromagnetic compatibility state report comprises a magnetic flux density change trend, a shielding configuration stability index, and a reset adjustment trigger record.

8. The electromagnetic radiation suppression method for a data center flywheel energy storage system of claim 7, wherein, The specific acquisition steps of the electromagnetic compatibility state report are as follows: S511: According to the numbers of all magnetic shielding components recorded in the magnetic shielding optimization configuration table, the magnetic flux density change rate between adjacent sampling periods is calculated, the absolute value of the continuous change rate of each magnetic shielding component is calculated, and it is judged whether they are all less than a preset change amplitude threshold. If yes, the corresponding component magnetic shielding is marked as tending to a stable state, and a magnetic flux convergence trend judgment result is obtained; S512: Based on the magnetic shielding components identified in the magnetic flux convergence trend judgment result, a sliding average window is set and a magnetic flux density change value sequence is extracted in each window, the consistency and convergence characteristics of the change amplitude are evaluated, if the change trend is stable and the difference degree is within the range, the corresponding magnetic shielding component is judged to meet the magnetic circuit recovery condition, is added to a reset execution candidate set according to the number, and a magnetic circuit reset target set is obtained; S513: According to the magnetic circuit reset target set, the current magnetic permeance configuration of each magnetic shielding component is compared with the initial factory configuration parameters, the deviation of the current parameters and the initial values of each magnetic permeance path is evaluated, if the deviation range is within the specified threshold, the magnetic circuit reset action can be triggered, the magnetic permeance path configuration is reset to the initial state, the state and path information of all recovered components are summarized in order of component number, and an electromagnetic compatibility state report is generated.

Images