Electromagnetic response prediction method of power plant rotating equipment under nuclear explosion pulse excitation

Abstract

The invention discloses an electromagnetic response prediction method for power plant rotating equipment under nuclear explosion pulse excitation, and relates to the technical field of computational electromagnetism. Decomposing the nuclear explosion pulse signal to obtain spectrum signals under multiple frequencies; generating an electric field excitation source for the rotating equipment of the power plant at each frequency according to the frequency spectrum signal; establishing a multi-conductor transmission line model of the power plant rotating equipment; under each frequency, an electric field excitation source of the frequency is used as input, the multi-conductor transmission line model is solved, and the electromagnetic response of the power plant rotating equipment under the frequency is obtained; according to the electromagnetic response under each frequency, the electromagnetic response of the power plant rotating equipment is determined, and the method can accurately predict the electromagnetic response of the power plant rotating equipment under nuclear explosion pulse excitation.

Description

Technical Field

[0001] The present invention relates to the technical field of computational electromagnetics, and in particular to a method for predicting electromagnetic responses of rotating equipment in a power plant under nuclear explosion pulse excitation. Background Art

[0002] As the current international situation becomes increasingly complex, the global nuclear situation is not optimistic. Escalating nuclear confrontation and intensified nuclear competition between major powers pose a significant threat to critical infrastructure and livelihoods. Power plants, as crucial electrical infrastructure, are a key concern for the safety of power plants worldwide, especially those subjected to the strong electromagnetic pulses (EMPs) caused by high-altitude nuclear explosions.

[0003] Power plants contain a large amount of rotating equipment, encompassing multiple systems across various disciplines. These include generators, coal mills, primary and secondary fans, main pump motors, and other rotating electrical machines. While traditional power plant rotating electrical machine design considers electromagnetic interference immunity in conventional industrial environments and protection against various operational overvoltages and lightning overvoltages, the impact of the strong electromagnetic pulse (EMP) environment generated by a high-altitude nuclear explosion on these machines has not been fully considered, and relevant research has yet to be fully conducted. Rotating electrical machines are exposed to strong electric fields at their operating voltages for extended periods. Coupled with the strong EMP from a high-altitude nuclear explosion, these machines can be subjected to transient overvoltage shocks, causing severe electrical damage to the inter-turn insulation and main insulation of the windings. This can compromise the overall insulation level of the equipment, leading to partial discharge and insulation breakdown. Ultimately, these rotating electrical machines can be damaged, disrupting the normal operation of power plant systems and severely impacting rescue and recovery efforts following a high-altitude nuclear explosion.

[0004] Therefore, in order to improve the ability of rotating equipment in power plants to cope with electromagnetic pulse threats, there is an urgent need for a method that can accurately predict the electromagnetic response of rotating equipment in power plants under nuclear explosion pulse excitation. Summary of the Invention

[0005] Based on this, it is necessary to provide a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation to address the above technical problems. This method can accurately predict the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation.

[0006] The present invention adopts the following technical solutions:

[0007] The present invention provides a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, comprising:

[0008] Decompose the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies;

[0009] Generate an electric field excitation source at each frequency for rotating equipment in the power plant according to the spectrum signal;

[0010] Build multi-conductor transmission line models for rotating equipment in power plants;

[0011] At each frequency, the electric field excitation source of the frequency is used as input to solve the multi-conductor transmission line model and obtain the electromagnetic response of the rotating equipment of the power plant at the frequency;

[0012] Based on the electromagnetic responses at various frequencies, the electromagnetic responses of rotating equipment in power plants are determined.

[0013] Preferably, the nuclear explosion pulse signal is decomposed to obtain spectrum signals at multiple frequencies, including:

[0014] Perform time-frequency domain conversion on the nuclear explosion pulse signal to obtain a spectrum signal;

[0015] Spectrum signals whose frequency amplitudes are greater than or equal to a preset amplitude threshold are obtained from the spectrum signals, and spectrum signals whose frequency amplitudes are greater than or equal to the preset amplitude threshold are determined as spectrum signals at multiple frequencies.

[0016] Preferably, the calculation formula of the electric field excitation source is:

[0017] E i (t) = A i sin2πf i t;

[0018] Among them, E i (t) is the electric field excitation source of the ith frequency in the spectrum signal at time t, A i is the amplitude corresponding to the i-th frequency in the spectrum signal, f i is the i-th frequency in the spectrum signal.

[0019] Preferably, the parameter acquisition process of the multi-conductor transmission line model of the rotating equipment of the power plant includes:

[0020] Based on the electrostatic field theory, the parameters of the multi-conductor transmission line model are determined based on the status information of the rotating equipment in the power plant; the parameters include distributed capacitance parameters and distributed conductance parameters;

[0021] The calculation formula of capacitance parameters is:

[0022]

[0023] Among them, C ij is the capacitance parameter between the i-th conductor and the j-th conductor in the rotating equipment of the power plant, ε is the dielectric constant of the insulating medium between the i-th conductor and the j-th conductor, S is the surface area per unit length of the i-th conductor, and E is the electrostatic field intensity along the surface of the i-th conductor when there is a potential difference U between the i-th conductor and the j-th conductor;

[0024] The distributed conductivity parameters are:

[0025]

[0026] Among them, G ij is the conductivity between the i-th conductor and the j-th conductor, and σ is the conductivity between the i-th conductor and the j-th conductor.

[0027] Preferably, the parameter acquisition process of the multi-conductor transmission line model of the rotating equipment of the power plant includes:

[0028] A finite element model of the rotating equipment in the power plant is established, and the state information of the rotating equipment in the power plant is input into the finite element model to obtain the parameters of the multi-conductor transmission line model.

[0029] Preferably, the electric field excitation source of the frequency is used as input to solve the multi-conductor transmission line model to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency, including:

[0030] The voltage control equation and the current control equation in the multi-conductor transmission line model are transformed to obtain the voltage partial differential equation and the current partial differential equation;

[0031] The spatial partial differential operators in the voltage partial differential equation and the current partial differential equation are discretized using the central difference scheme, and the discretized voltage partial differential equation and current partial differential equation are merged into time-domain ordinary differential equations.

[0032] According to the theory of ordinary differential equations, determine the formal function of the solution of time domain ordinary differential equations;

[0033] The integral term in the formal function is approximated by constructing the Gaussian integral formula through two-point interpolation, and a discrete iterative recursive formula is obtained;

[0034] Substituting the electric field excitation source of the frequency into the discrete iterative recursive formula, the electromagnetic response of the rotating equipment in the power plant at the frequency is obtained.

[0035] The present invention provides a device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, comprising:

[0036] A decomposition module is used to decompose the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies;

[0037] A generation module, used for generating an electric field excitation source for rotating equipment in a power plant at each frequency according to the spectrum signal;

[0038] Building blocks for modeling multi-conductor transmission lines for rotating equipment in power plants;

[0039] A solution module is used to solve the multi-conductor transmission line model at each frequency by taking the electric field excitation source of the frequency as input to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency;

[0040] The determination module is used to determine the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each frequency.

[0041] The present invention provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation is implemented.

[0042] The present invention provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation is implemented.

[0043] At least one of the above technical solutions adopted by the present invention can achieve the following beneficial effects:

[0044] In the present invention, since signals of different frequencies interact differently with rotating equipment during transmission, decomposing the nuclear explosion pulse signal into spectrum signals at multiple frequencies can more clearly determine the contribution of each frequency component to the electromagnetic response of the rotating equipment, avoid the analysis difficulties caused by mutual interference between different frequency components in the time domain signal, and improve the accuracy of the predicted electromagnetic response of the rotating equipment of the power plant. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0046] Figure 1 A schematic flow chart of a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by the present invention;

[0047] Figure 2 A schematic structural diagram of a series resistance matrix, a series inductance matrix, a parallel conductance matrix, and a parallel capacitance matrix per unit length provided by the present invention;

[0048] Figure 3 A schematic flow chart of another method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by the present invention;

[0049] Figure 4 A schematic diagram of a device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by the present invention;

[0050] Figure 5 A schematic diagram of a computer device for realizing a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by the present invention. DETAILED DESCRIPTION

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

[0052] In order to improve the ability of rotating equipment in power plants to cope with the threat of electromagnetic pulses, carry out targeted research and development of protection methods and equipment, and propose relevant prevention, protection and recovery measures within a reasonable range, it is necessary to conduct research on the characteristics of strong electromagnetic pulses from nuclear explosions and the destructive mechanism of electrical and electronic equipment.

[0053] To address these research needs for rotating equipment in power plants, this paper proposes a method for predicting the electromagnetic response of rotating equipment in power plants under nuclear explosion pulse excitation. First, the nuclear explosion electromagnetic pulse is decomposed into multiple frequency signals, each of which serves as the excitation source at different frequencies. Second, a multi-conductor transmission line model is established. Finally, this model is solved using a time-domain precise integration method to accurately and efficiently predict the electromagnetic response of rotating equipment in power plants under the strong nuclear explosion electromagnetic pulse excitation.

[0054] The technical solutions provided by various embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0055] Figure 1 The flowchart of the method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation in the present invention specifically includes the following steps:

[0056] S101, decomposing the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies.

[0057] In a specific example, a nuclear explosion occurs at an altitude of more than 30 kilometers. The instantaneously emitted gamma rays interact with air molecules to generate Compton currents. Under the influence of factors such as atmospheric density, X-rays, and the Earth's magnetic field, these currents oscillate and radiate strong electromagnetic pulses. The electromagnetic pulses generated by a high-altitude nuclear explosion have characteristics such as a very short rise time, a slow decay rate, and a long duration. They are typically described using a double-exponential pulse waveform, expressed in the time domain as:

[0058] E(t)=kE0(e -αt -e -βt) (1)

[0059] Where, E(t) is the electric field strength at time t; E0 is the peak value of the field strength, which is generally 50 kV / m; β is the pulse leading edge parameter; α is the pulse trailing edge parameter; k is the correction parameter; t r is the rising edge time, the time when the pulse field strength rises from 10% to 90%; t f is the falling edge time, the time it takes for the pulse field strength to drop from 90% to 10%; τ FWHM The pulse width at half maximum is the time interval between two pulse field strengths reaching 50%. The typical high-altitude nuclear explosion pulse model parameters are shown in Table 1.

[0060] Table 1

[0061]

[0062] By selecting a set of parameters as the high-altitude nuclear explosion pulse model parameters, the time-domain expression of the pulse excitation can be obtained according to formula (1). By selecting an appropriate time step, a set of discrete data x(n) of the nuclear explosion electromagnetic pulse is obtained, and this discrete data is used as the nuclear explosion pulse signal.

[0063] In an exemplary embodiment, a nuclear explosion pulse signal is decomposed to obtain spectrum signals at multiple frequencies, including: performing time-frequency domain conversion on the nuclear explosion pulse signal to obtain a spectrum signal; obtaining a spectrum signal whose frequency amplitude is greater than or equal to a preset amplitude threshold from the spectrum signal, and determining the spectrum signal whose frequency amplitude is greater than or equal to the preset amplitude threshold as the spectrum signal at multiple frequencies.

[0064] Optionally, the burst pulse excitation signal is subjected to time-frequency domain conversion including fast Fourier transform (FFT), wavelet analysis, etc., to obtain a set of discrete sequences X(k) (spectrum signals) in the frequency domain, corresponding to frequencies f1, f2..., f n .

[0065] Set the preset amplitude threshold A of the frequency component amplitude in the spectrum signal min , discard the frequency components smaller than the preset amplitude threshold, extract the remaining components, and obtain a new discrete sequence X(i), which contains I components. The discrete sequence X(i) is the spectrum signal at multiple frequencies required in the end.

[0066] It should be noted that, based on the characteristics of nuclear explosion pulse signals, signals within a certain frequency range are considered nuclear explosion pulse signals, while signals outside this range are not considered nuclear explosion pulse signals. Therefore, spectral signals with frequency amplitudes greater than or equal to a preset amplitude threshold are identified as spectral signals at multiple frequencies, and these spectral signals at multiple frequencies are used to calculate the electromagnetic response of the power plant's rotating equipment.

[0067] S102 : generating an electric field excitation source at each frequency for rotating equipment in the power plant according to the spectrum signal.

[0068] The corresponding frequency f of the discrete sequence X(i) i and amplitude A i Extract and obtain a set of components as the electric field excitation source. The calculation formula of the electric field excitation source is:

[0069] E i (t) = A i sin 2πf i t (2)

[0070] Among them, E i (t) is the electric field excitation source of the ith frequency in the spectrum signal at time t, A i is the amplitude corresponding to the i-th frequency in the spectrum signal, f i is the i-th frequency in the spectrum signal.

[0071] S103, establishing a multi-conductor transmission line model of rotating equipment in a power plant.

[0072] Taking the rotating equipment of a power plant as an example, a multi-conductor transmission line model of the rotating motor winding of the power plant is established. The control equation of the multi-conductor transmission line model of the rotating motor winding is as follows:

[0073]

[0074] Among them, U is the n×1 voltage matrix, which contains all voltage components; I is the n×1 current matrix, which contains all current components; U s is the voltage excitation source term; I s is the current excitation source term; R is the n×n unit length series resistance matrix; L is the n×n unit length series inductance matrix; G is the n×n unit length parallel conductance matrix; C is the n×n unit length parallel capacitance matrix. The structural diagram of the unit length series resistance matrix, series inductance matrix, parallel conductance matrix and parallel capacitance matrix is ​​shown in the figure below. Figure 2 The specific expression is as follows:

[0075]

[0076] In the formula, R in the R matrix ii is the conductor resistance, r is the grounded core return resistance; the main diagonal element L in the L matrix ii is the conductor's self-inductance, and the off-diagonal element L ij is the mutual inductance between the two corresponding conductors; the diagonal element of G is the conductivity of the conductor to the ground G ii The sum of the conductance between the conductor and other conductors, other elements G ij is the conductivity between the two corresponding conductors; the diagonal element C is the capacitance of the conductor to the ground ii The sum of the capacitances with other conductors, other elements C ij is the capacitance between the two corresponding conductors.

[0077] It should be noted that the excitation source term U s and I s It can be derived from formula (2) according to the form of the specific problem.

[0078] Based on the above multi-conductor transmission line model, the parameters of the multi-conductor transmission line model can be obtained in two ways: analytical method and finite element analysis method.

[0079] In one embodiment, a process for obtaining parameters of a multi-conductor transmission line model of rotating equipment in a power plant includes: determining parameters of the multi-conductor transmission line model based on state information of the rotating equipment in the power plant according to electrostatic field theory; the parameters include distributed capacitance parameters and distributed conductance parameters;

[0080] According to electrostatic field theory, the calculation formula of distributed capacitance parameters is:

[0081]

[0082] Among them, C ij is the capacitance parameter between the i-th conductor and the j-th conductor in the rotating equipment of the power plant, ε is the dielectric constant of the insulating medium between the i-th conductor and the j-th conductor, S is the surface area per unit length of the i-th conductor, and E is the electrostatic field intensity along the surface of the i-th conductor when there is a potential difference U between the i-th conductor and the j-th conductor.

[0083] Rotating equipment in power plants, such as motors, consists of multiple mutually insulated independent conductors, including stator cores, windings, and rotors. Let the total number of these conductors be N, and all conductors are numbered 1, 2, ..., N in sequence; then "i" and "j" are integers between 1 and N, referring to the corresponding numbered conductors.

[0084] According to electrostatic field theory, the calculation formula of distributed conductivity parameters is:

[0085]

[0086] Among them, G ijis the conductivity between the i-th conductor and the j-th conductor, and σ is the conductivity between the i-th conductor and the j-th conductor.

[0087] It should be noted that the distributed resistance and distributed inductance parameters can be solved by constructing a magnetic circuit model of the rotating motor winding according to actual conditions and using the reluctance network method.

[0088] In another embodiment, a parameter acquisition process of a multi-conductor transmission line model of rotating equipment in a power plant includes: establishing a finite element model of the rotating equipment in the power plant, inputting state information of the rotating equipment in the power plant into the finite element model, and obtaining parameters of the multi-conductor transmission line model.

[0089] Using the Q3D module in ANSYS software, a finite element model of a rotating motor winding is established, which can yield distributed parameter results at different frequencies. Typically, at low frequencies, finite element analysis is used to calculate distributed parameters, while analytical methods can be used for other calculations. By leveraging the characteristics of both methods and maximizing their respective strengths, the distributed parameters of a rotating motor multi-conductor model can be efficiently and accurately determined.

[0090] S104 , at each frequency, using the electric field excitation source of the frequency as input, solving the multi-conductor transmission line model to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency.

[0091] In an exemplary embodiment, a frequency electric field excitation source is used as input, and a multi-conductor transmission line model is solved by a time domain precise integration method to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency.

[0092] Specifically, the voltage control equation and the current control equation in the multi-conductor transmission line model are transformed to obtain the voltage partial differential equation and the current partial differential equation; the spatial partial differential operators in the voltage partial differential equation and the current partial differential equation are discretized using the central difference format, and the discretized voltage partial differential equation and the current partial differential equation are merged into the time-domain ordinary differential equation; according to the theory of ordinary differential equations, the formal function of the solution of the time-domain ordinary differential equation is determined; the integral term in the formal function is approximated by constructing a Gaussian integral formula through two-point interpolation to obtain a discrete iterative recursive formula; the electric field excitation source of the frequency is substituted into the discrete iterative recursive formula to obtain the electromagnetic response of the rotating equipment of the power plant under the frequency.

[0093] The two control equations in formula (3) can be respectively the voltage control equation and the current control equation. The two control equations in formula (3) are transformed to obtain the voltage partial differential equation and the current partial differential equation, as shown in formula (10).

[0094]

[0095] Based on formula (10), the time partial differential operator is kept unchanged, the spatial partial differential operator is discretized using the central difference format, and the two equations in formula (10) after spatial discretization are combined into an ordinary differential equation, which is expressed as follows:

[0096]

[0097] in, Y is a one-dimensional column vector containing all voltage components U and current components I in the calculation area, H is a coefficient matrix determined by the spatial step size and model distribution parameters and does not change with time, f(t) is a one-dimensional column vector introduced by the excitation source. The calculation area can be considered as the nuclear explosion point and the area within its impact range, which should include the rotating equipment of the power plant under analysis.

[0098] According to the theory of ordinary differential equations, the analytical solution of formula (11) can be obtained as formula (12).

[0099]

[0100] By using the Gaussian integral formula to approximate the integral term on the right side of formula (12), and applying two-point interpolation to construct the Gaussian integral formula, we can obtain the discrete iterative recursive formula, as shown in formula (13).

[0101]

[0102] Among them, Y k+1 is the electromagnetic response at frequency, Y k t k The values ​​of the voltage and current components at the moment Y(kΔt), T is the exponential matrix e of the coefficient matrix H HΔt The excitation vector f refers to the electric field intensity at some locations of the device (such as the surface or a circuit port) due to the incident electromagnetic pulse of the nuclear explosion. The position and intensity of the electric field excitation source E are represented by the excitation vector f.

[0103] S105 , determining the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each frequency.

[0104] The electromagnetic responses at various frequencies obtained above can be superimposed to obtain the electromagnetic responses of the rotating equipment in the power plant.

[0105] In an exemplary embodiment, the present invention also provides a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, such as Figure 3 As shown, this embodiment includes the following steps:

[0106] S301, determine the nuclear explosion electromagnetic pulse model, perform fast Fourier decomposition on the pulse excitation, and obtain the excitation source expression at different frequencies.

[0107] S302: Establish a multi-conductor transmission line model of the rotating electrical machine winding.

[0108] S303: Determine the distribution parameters in the multi-conductor transmission line model using an analytical method or a finite element method.

[0109] S304 , selecting an excitation source of a certain frequency as input, solving the established multi-conductor transmission line model using a time domain precise integration method, and obtaining the electromagnetic response characteristics of the rotating equipment of the power plant under the excitation of the frequency.

[0110] S305, repeating S304 until all the excitation source conditions at all frequencies in S301 are simulated, and finally obtaining the electromagnetic response characteristics of the rotating equipment of the power plant under the nuclear explosion pulse excitation.

[0111] While traditional power plant rotating motor design takes into account electromagnetic interference resistance in conventional industrial environments and protection against various operational overvoltages and lightning surges, the impact of the strong electromagnetic pulse environment generated by a high-altitude nuclear explosion on the equipment is not considered. In reality, rotating motors are exposed to strong electric fields at their operating voltages for extended periods. Coupled with the strong electromagnetic pulses of a high-altitude nuclear explosion, these motors can be subjected to transient overvoltage shocks, causing severe electrical damage to the inter-turn insulation and main insulation of the windings. This can ultimately damage the rotating motor, disrupting normal power plant operations and severely impacting rescue and recovery efforts following a high-altitude nuclear explosion. In response to the above problems, this application establishes an algorithm system for predicting the electromagnetic response characteristics of rotating equipment in power plants under nuclear explosion pulse excitation based on a multi-conductor transmission line model and a time-domain fine integration method. The algorithm system has the following advantages: (1) It can predict the electromagnetic response characteristics of rotating equipment in power plants under nuclear explosion pulse excitation, especially the overvoltage characteristics, guide designers to carry out research and development of protection methods and protective equipment, and propose relevant prevention, protection and recovery measures within a reasonable range; (2) Since the time-domain fine integration method has a loose stability condition, the prediction algorithm of this application has the advantages of a small number of iterations and high computational efficiency; (3) Since the computational accuracy of the time-domain fine integration method does not change with the change of the time step, the prediction algorithm of this application can use the largest possible time step while maintaining the computational accuracy, and has the advantages of high computational accuracy; (4) The algorithm proposed in this application has good versatility.

[0112] When applying the electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation provided by the present invention, it is not necessary to Figure 1The steps are executed in the order shown. The specific execution order of the steps can be determined according to needs, and the present invention does not limit this.

[0113] The above is a method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by one or more embodiments of the present invention. Based on the same idea, the present invention also provides a corresponding device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, such as Figure 4 shown.

[0114] Figure 4 This is a schematic diagram of a device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation provided by the present invention. The device 400 includes:

[0115] The decomposition module 401 is used to decompose the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies.

[0116] The generating module 402 is configured to generate an electric field excitation source for the rotating equipment of the power plant at each frequency according to the spectrum signal.

[0117] The construction module 403 is used to establish a multi-conductor transmission line model of rotating equipment in a power plant.

[0118] The solving module 404 is used to solve the multi-conductor transmission line model at each frequency by taking the electric field excitation source of the frequency as input to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency.

[0119] The determination module 405 is configured to determine the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each frequency.

[0120] The specific limitations of the device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation can be found in the aforementioned limitations of the method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, and will not be repeated here. Each module in the aforementioned device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation can be implemented in whole or in part through software, hardware, or a combination thereof. Each of the aforementioned modules can be embedded in or independent of a processor in a computer device in hardware form, or can be stored in a memory in the computer device in software form, so that the processor can call and execute the corresponding operations of each of the aforementioned modules.

[0121] The present invention also provides a computer-readable storage medium, which stores a computer program, which can be used to execute the above Figure 1 A method for predicting the electromagnetic response of rotating equipment in power plants under nuclear explosion pulse excitation is provided.

[0122] The present invention also provides Figure 5 The structural diagram of the computer equipment shown in FIG. Figure 5As shown in the figure, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory. Of course, it may also include other hardware required for the business. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to achieve the above Figure 1 A method for predicting the electromagnetic response of rotating equipment in power plants under nuclear explosion pulse excitation is provided.

[0123] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0124] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of the present invention.

Claims

1. A method for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, characterized in that: include: Decompose the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies; generating an electric field excitation source for rotating equipment in the power plant at each frequency according to the frequency spectrum signal; Establishing a multi-conductor transmission line model of rotating equipment in the power plant; At each frequency, taking the electric field excitation source of the frequency as input, solving the multi-conductor transmission line model, and obtaining the electromagnetic response of the rotating equipment of the power plant at the frequency; The electromagnetic response of the rotating equipment of the power plant is determined according to the electromagnetic response at each of the frequencies.

2. The method according to claim 1, characterized in that The nuclear explosion pulse signal is decomposed to obtain spectrum signals at multiple frequencies, including: Performing time-frequency domain conversion on the nuclear explosion pulse signal to obtain a spectrum signal; A spectrum signal whose frequency amplitude is greater than or equal to a preset amplitude threshold is obtained from the spectrum signal, and a spectrum signal whose frequency amplitude is greater than or equal to the preset amplitude threshold is determined as the spectrum signal at the multiple frequencies.

3. The method according to claim 1, characterized in that The calculation formula of the electric field excitation source is: E i (t)=A i sin2πf i t; Among them, E i (t) is the electric field excitation source of the ith frequency in the spectrum signal at time t, A i is the amplitude corresponding to the i-th frequency in the spectrum signal, f i is the i-th frequency in the spectrum signal.

4. The method according to claim 1, characterized in that The parameter acquisition process of the multi-conductor transmission line model of the power plant rotating equipment includes: Determining the parameters of the multi-conductor transmission line model according to the electrostatic field theory and the state information of the rotating equipment of the power plant; the parameters include distributed capacitance parameters and distributed conductance parameters; The calculation formula of the capacitance parameter is: Among them, C ij is the capacitance parameter between the ith conductor and the jth conductor in the rotating equipment of the power plant, ε is the dielectric constant of the insulating medium between the ith conductor and the jth conductor, S is the surface area per unit length of the ith conductor, and E is the electrostatic field intensity along the surface of the ith conductor when there is a potential difference U between the ith conductor and the jth conductor; The distributed conductivity parameter is: Among them, G ij is the conductivity between the i-th conductor and the j-th conductor, and σ is the conductivity between the i-th conductor and the j-th conductor.

5. The method according to claim 1, characterized in that The parameter acquisition process of the multi-conductor transmission line model of the power plant rotating equipment includes: A finite element model of the rotating equipment of the power plant is established, and state information of the rotating equipment of the power plant is input into the finite element model to obtain parameters of the multi-conductor transmission line model.

6. The method according to claim 1, characterized in that The method of using the electric field excitation source of the frequency as input, solving the multi-conductor transmission line model, and obtaining the electromagnetic response of the rotating equipment of the power plant at the frequency includes: Transforming the voltage control equation and the current control equation in the multi-conductor transmission line model to obtain a voltage partial differential equation and a current partial differential equation; The spatial partial differential operators in the voltage partial differential equation and the current partial differential equation are discretized using a central difference format, and the discretized voltage partial differential equation and current partial differential equation are combined into a time-domain ordinary differential equation; Determining the formal function of the solution of the time domain ordinary differential equation according to the theory of ordinary differential equations; The integral term in the formal function is approximated by constructing the Gaussian integral formula through two-point interpolation, and a discrete iterative recursive formula is obtained; Substituting the electric field excitation source of the frequency into the discrete iterative recursive formula, the electromagnetic response of the rotating equipment of the power plant at the frequency is obtained.

7. A device for predicting the electromagnetic response of rotating equipment in a power plant under nuclear explosion pulse excitation, characterized in that: include: A decomposition module is used to decompose the nuclear explosion pulse signal to obtain spectrum signals at multiple frequencies; A generating module, used for generating an electric field excitation source for rotating equipment of a power plant at each frequency according to data of the frequency spectrum signal; A construction module for establishing a multi-conductor transmission line model of rotating equipment of the power plant; A solution module, for solving the multi-conductor transmission line model at each frequency by taking the electric field excitation source of the frequency as input, and obtaining the electromagnetic response of the rotating equipment of the power plant at the frequency; The determination module is used to determine the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each of the frequencies.

8. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

9. A computer device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method according to any one of claims 1 to 6 is implemented.

Images