A method for predicting electromagnetic response of rotating equipment in power plant under nuclear burst pulse excitation

Abstract

The application discloses a kind of electromagnetic response prediction methods of nuclear explosion pulse excitation under power plant rotating equipment, it is related to the technical field of computational electromagnetics.Nuclear explosion pulse signal is decomposed, and spectrum signal under multiple frequencies is obtained;According to the spectrum signal, generate the electric field excitation source of power plant rotating equipment at each frequency;Establish the multi-conductor transmission line model of power plant rotating equipment;At each frequency, the electric field excitation source of frequency is regarded as input, and the multi-conductor transmission line model is solved, and the electromagnetic response of power plant rotating equipment at frequency is obtained;According to the electromagnetic response at each frequency, determine the electromagnetic response of power plant rotating equipment, and the electromagnetic response of power plant rotating equipment under nuclear explosion pulse excitation can be accurately predicted by the method.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of computational electromagnetics, in particular to a method for predicting electromagnetic response of rotating equipment of a power plant under nuclear explosion pulse excitation. BACKGROUND

[0002] With the current international situation becoming increasingly complex, the global nuclear situation is not optimistic, and the nuclear confrontation between major powers is escalating and the nuclear race is intensifying, which poses a huge threat to key basic livelihood facilities. Among them, power plants as important power infrastructure are concerned about the impact of strong electromagnetic pulses caused by high-altitude nuclear explosions on their safety.

[0003] Power plants contain a large number of rotating equipment, involving multiple systems of various professions, including generators, coal mills, primary / secondary air fans, and main pump motors, etc. In the traditional design process of rotating machines in power plants, although the anti-electromagnetic interference capability in the conventional industrial environment and the protection capability of various operating overvoltages and lightning overvoltages have been considered, the impact of the strong electromagnetic pulse environment generated by high-altitude nuclear explosions on the equipment has not been taken into account in the design, and relevant research work has not been fully carried out. In fact, rotating machines are long-term in the strong electric field of working voltage, and under the coupling effect of strong electromagnetic pulses of high-altitude nuclear explosions, the rotating machines may be subjected to transient overvoltage impact, causing serious electrical damage to the inter-turn insulation and main insulation of the winding, thereby damaging the overall insulation level of the equipment, causing partial discharge and insulation breakdown, etc. accidents, ultimately leading to damage to the rotating machine, the power plant system cannot operate normally, and seriously affecting the rescue and recovery work after the high-altitude nuclear explosion.

[0004] Therefore, in order to improve the ability of rotating equipment of power plants to cope with electromagnetic pulse threats, it is urgent to provide a method for accurately predicting the electromagnetic response of rotating equipment of a power plant under nuclear explosion pulse excitation. SUMMARY

[0005] Therefore, in order to improve the ability of rotating equipment of power plants to cope with electromagnetic pulse threats, it is urgent to provide a method for accurately predicting the electromagnetic response of rotating equipment of a power plant under nuclear explosion pulse excitation.

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

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

[0008] The nuclear explosion pulse signal is decomposed to obtain a plurality of frequency spectrum signals under different frequencies;

[0009] According to the frequency spectrum signal, an electric field excitation source of the rotating equipment of the power plant under each frequency is generated;

[0010] A multi-conductor transmission line model of a rotating equipment of a power plant is established;

[0011] At each frequency, the electromagnetic response of the rotating equipment of the power plant at the frequency is obtained by solving the multi-conductor transmission line model with the electric field excitation source of the frequency as input;

[0012] The electromagnetic response of the rotating equipment of the power plant is determined according to the electromagnetic responses at the frequencies.

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

[0014] The nuclear explosion pulse signal is subjected to time-frequency domain conversion to obtain the spectral signals;

[0015] The spectral signals with amplitudes greater than or equal to a preset amplitude threshold value are obtained from the spectral signals, and the spectral signals with amplitudes greater than or equal to the preset amplitude threshold value are determined as the spectral signals at the 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] Wherein, E i (t) is the electric field excitation source of the i-th frequency in the spectral signal at time t, A i is the amplitude corresponding to the i-th frequency in the spectral signal, and f i is the i-th frequency in the spectral signal.

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

[0020] According to the electrostatic field theory, the parameters of the multi-conductor transmission line model are determined according to the state information of the rotating equipment of the power plant; the parameters include distributed capacitance parameters and distributed conductance parameters;

[0021] The calculation formula of the capacitance parameter is:

[0022]

[0023] Wherein, C ij is the capacitance parameter between the i-th conductor and the j-th conductor of 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 conductance parameter is:

[0025]

[0026] wherein G ij is the conductance between the i-th conductor and the j-th conductor, and sigma is the conductivity between the i-th conductor and the j-th conductor.

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

[0028] establishing a finite element model of the power plant rotating equipment, inputting state information of the power plant rotating equipment into the finite element model, and obtaining parameters of the multi-conductor transmission line model.

[0029] Preferably, the multi-conductor transmission line model is solved by taking the electric field excitation source of the frequency as input, and electromagnetic response of the power plant rotating equipment at the frequency is obtained, comprising:

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

[0031] The spatial partial differential operators in the voltage partial differential equation and the current partial differential equation are discretized by using a central difference format, and the discretized voltage partial differential equation and the current partial differential equation are combined into a time-domain ordinary differential equation;

[0032] According to the ordinary differential equation theory, a form function of the solution of the time-domain ordinary differential equation is determined;

[0033] The integral term in the form function is approximated by constructing a Gauss integral formula through two-point interpolation to obtain a discrete iterative recurrence formula;

[0034] The electric field excitation source of the frequency is substituted into the discrete iterative recurrence formula to obtain the electromagnetic response of the power plant rotating equipment at the frequency.

[0035] The present application provides a kind of electromagnetic response prediction device of power plant rotating equipment under nuclear explosion pulse excitation, comprising:

[0036] The decomposition module is used to decompose the nuclear explosion pulse signal to obtain a plurality of frequency spectrum signals at different frequencies.

[0037] The generation module is used to generate an electric field excitation source for the power plant rotating equipment at each frequency according to the frequency spectrum signals.

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

[0039] a solving module, configured to solve the multi-conductor transmission line model by taking the electric field excitation source of the frequency as input to obtain electromagnetic response of the rotating equipment of the power plant at the frequency;

[0040] a determining module, configured to determine the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each frequency.

[0041] The present application provides a computer readable storage medium, the storage medium stores a computer program, the computer program is executed by a processor to realize the above-mentioned electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation.

[0042] The present application provides a computer device, comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executes the program to realize the above-mentioned electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation.

[0043] The above-mentioned at least one technical scheme adopted by the present application can achieve the following beneficial effects:

[0044] In the present application, because the interaction of signals of different frequencies with the rotating equipment is different during transmission, the nuclear explosion pulse signal is decomposed into spectral signals at multiple frequencies, which can more clearly determine the contribution of each frequency component to the electromagnetic response of the rotating equipment, avoid the analysis difficulty caused by the mutual interference of different frequency components in the time domain signal, and improve the accuracy of predicting the electromagnetic response of the rotating equipment of the power plant. BRIEF DESCRIPTION OF DRAWINGS

[0045] The accompanying drawings, which are included to provide a further understanding of the present application, constitute a part of the present application and illustrate embodiments of the present application and its description, which serve to explain the present application and do not constitute improper limitations on the present application. In the drawings:

[0046] Figure 1 A nuclear explosion pulse excitation electromagnetic response prediction method flowchart of the rotating equipment of the power plant is provided in the present application;

[0047] Figure 2 A structure diagram of a unit length series resistance matrix, a series inductance matrix, a parallel conductance matrix and a parallel capacitance matrix is provided in the present application;

[0048] Figure 3 Another nuclear explosion pulse excitation electromagnetic response prediction method flowchart of the rotating equipment of the power plant is provided in the present application;

[0049] Figure 4 A nuclear explosion pulse excitation electromagnetic response prediction device schematic diagram of the rotating equipment of the power plant is provided in the present application;

[0050] Figure 5 A computer device schematic diagram for implementing a nuclear explosion pulse excitation electromagnetic response prediction method of a power plant rotating equipment is provided. DETAILED DESCRIPTION

[0051] In order to make the objects, technical solutions and advantages of the present application clearer, the technical solutions of the present application will be described below in connection with the specific embodiments of the present application and corresponding drawings. Obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative labor fall within the scope of protection of the present application.

[0052] In order to improve the ability of the rotating equipment of the power plant to cope with electromagnetic pulse threats, it is necessary to carry out research on the characteristics of nuclear explosion strong electromagnetic pulses and the damage mechanism to electrical and electronic equipment, and to develop targeted protection methods and protection equipment, and to propose relevant response measures for prevention, protection and recovery within a reasonable range.

[0053] Therefore, the present application proposes a nuclear explosion pulse excitation electromagnetic response prediction method of a power plant rotating equipment in view of the above research needs of the rotating equipment of the power plant. First, the nuclear explosion electromagnetic pulse is decomposed into signals of multiple frequencies, which are respectively used as excitation sources of different frequency points. Second, a multi-conductor transmission line model is established. Finally, the above model is solved based on the time domain precise integration method, and the electromagnetic response of the nuclear explosion strong electromagnetic pulse excitation in the rotating equipment of the power plant is accurately and efficiently predicted.

[0054] The technical solutions provided by the embodiments of the present application will be described in detail below in connection with the drawings.

[0055] Figure 1 A flowchart of a nuclear explosion pulse excitation electromagnetic response prediction method of a power plant rotating equipment in the present application is provided, which specifically includes the following steps:

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

[0057] In a specific embodiment, a nuclear explosion occurs at an altitude of more than 30 kilometers above the ground. The instantaneously emitted gamma rays interact with air molecules to produce Compton current. Under the influence of factors such as atmospheric density, X-rays and geomagnetic field, strong electromagnetic pulses are oscillated and radiated. The electromagnetic pulses generated by high-altitude nuclear explosions have the characteristics of very short rise time, slow decay rate and long duration. A double exponential pulse waveform is usually used to describe it, and the time domain expression is:

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

[0059] wherein E(t) is the electric field strength value at time t; E0 is the peak value of field strength, generally 50 kV / m; β is the pulse front parameter; α is the pulse trailing parameter; k is the correction parameter; t r is the rising time, the time for the pulse field strength to rise from 10% to 90%; t f is the falling time, the time for the pulse field strength to fall from 90% to 10%; τ FWHM is the pulse half-width, the time interval for the pulse field strength to reach 50% twice. The typical high-altitude nuclear explosion pulse model parameters are shown in Table 1.

[0060] Table 1

[0061]

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

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

[0064] Optionally, the time-frequency domain conversion of the explosion pulse excitation signal includes Fast Fourier Transform (FFT), wavelet analysis method, etc. The time-frequency domain conversion of the explosion pulse excitation signal is performed to obtain a set of discrete sequences X(k) (frequency spectrum signals) in the frequency domain, corresponding to frequencies f1, f2,..., f n .

[0065] The preset amplitude threshold A min of the frequency component amplitude in the frequency spectrum signal is set, the frequency components less than the preset amplitude threshold are discarded, the remaining components are extracted, and a new set of discrete sequences X(i) containing I components is obtained, and the discrete sequence X(i) is the final plurality of frequency spectrum signals at the plurality of frequencies.

[0066] It should be noted that according to the characteristics of the nuclear explosion pulse signal, the signal in a certain frequency range is considered as the characteristics of the nuclear explosion pulse signal, and the signals beyond this range are not considered in the nuclear explosion pulse. Therefore, the frequency spectrum signal with an amplitude greater than or equal to a preset amplitude threshold is determined as the frequency spectrum signal at a plurality of frequencies, and the electromagnetic response of the rotating equipment of the power plant is calculated based on the frequency spectrum signal at the plurality of frequencies.

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

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

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

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

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

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

[0073]

[0074] Wherein, U is an n×1 voltage matrix containing all voltage components; I is an n×1 current matrix containing all current components; U s is a voltage excitation source term; I s is a current excitation source term; R is an n×n unit length series resistance matrix; L is an n×n unit length series inductance matrix; G is an n×n unit length parallel conductance matrix; C is an n×n unit length parallel capacitance matrix. The structure diagram of the unit length series resistance matrix, the unit length series inductance matrix, the unit length parallel conductance matrix and the unit length parallel capacitance matrix is as shown in Figure 2 The specific expression form is as follows:

[0075]

[0076] In the formula, R in the R matrix ii is the conductor resistance, r is the return resistance of the grounded core; 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 ijσ 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 constructed according to actual conditions, and the magnetic circuit model of the rotating machine winding is solved by using the magnetic resistance network method.

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

[0089] The finite element model of the rotating machine winding is established by using the Q3D module of the ANSYS software, and the distributed parameter results under different frequencies can be obtained. Generally, in the case of low frequency, the finite element analysis method can be used to calculate the distributed parameters, and in other cases, the analytical method can be used to calculate the distributed parameters. By using the characteristics of the two methods, the respective advantages are maximized, and the distributed parameters of the multi-conductor model of the rotating machine are obtained efficiently and accurately.

[0090] In S104, the electromagnetic response of the rotating equipment of the power plant at the frequency is obtained by solving the multi-conductor transmission line model with the frequency electric field excitation source as input.

[0091] In an exemplary embodiment, the electromagnetic response of the rotating equipment of the power plant at the frequency is obtained by solving the multi-conductor transmission line model with the frequency electric field excitation source as input by using the time domain precise integration method.

[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 by using the central difference format, and the discretized voltage partial differential equation and the current partial differential equation are combined into a time domain ordinary differential equation; according to the theory of ordinary differential equation, the form function of the solution of the time domain ordinary differential equation is determined; the integral term in the form function is approximated by constructing the Gauss integral formula by two-point interpolation to obtain a discrete iterative recurrence formula; the electromagnetic response of the rotating equipment of the power plant at the frequency is obtained by substituting the frequency electric field excitation source into the discrete iterative recurrence formula.

[0093] The two control equations in formula (3) can be a voltage control equation and a current control equation, respectively. 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 by using the central difference format, and the two equations in formula (10) after spatial discretization are combined into one ordinary differential equation, and the expression is as follows:

[0096]

[0097] wherein, Y is a one-dimensional column vector containing all voltage components U and current components I in the calculation region, H is a coefficient matrix determined by the spatial step and the model distribution parameters and not changing with time, f(t) is a one-dimensional column vector introduced by the excitation source. Wherein, the calculation region can be considered as the region of the nuclear explosion point and its influence range, which should contain the analyzed rotating equipment of the power plant.

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

[0099]

[0100] The right integral term of formula (12) is approximated by using the Gauss integral formula, and the Gauss integral formula is constructed by using two-point interpolation, and a discrete iterative recursive formula can be obtained, as shown in formula (13).

[0101]

[0102] wherein, Y k+1 is the electromagnetic response at the frequency, Y k is the value of the voltage and current components Y(kΔt) at t k time, T is the exponential matrix e HΔt of the coefficient matrix H, and the excitation vector f refers to the electric field intensity at the position (such as the surface, a circuit port) of the device part due to the incidence of the nuclear explosion electromagnetic pulse, and the position and intensity of the electric field excitation source E are represented by the excitation vector f.

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

[0104] The electromagnetic response at each frequency obtained above can be superimposed to obtain the electromagnetic response of the rotating equipment of the power plant.

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

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

[0107] S302, establish a multi-conductor transmission line model of the rotating machine winding.

[0108] S303, determine the distributed parameters in the multi-conductor transmission line model by using an analytical method or a finite element method.

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

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

[0111] In the traditional design process of the rotating machine of the power plant, although the anti-electromagnetic interference capability in the conventional industrial environment and the protection capability against various types of operating overvoltage and lightning overvoltage have been considered, the influence of the strong electromagnetic pulse environment generated by the high-altitude nuclear explosion on the equipment has not been taken into account in the design. In fact, the rotating machine is long-term in a strong electric field of working voltage, and under the coupling action of the strong electromagnetic pulse of the high-altitude nuclear explosion, the rotating machine may be impacted by a transient overvoltage, which may cause serious electrical damage to the turn-to-turn insulation and the main insulation of the winding, and finally lead to damage of the rotating machine, failure of the power plant system to operate normally, and serious influence on the rescue and recovery work after the high-altitude nuclear explosion. The present application aims at the above problems, and establishes an electromagnetic response characteristic prediction algorithm system for the rotating equipment of the power plant under the nuclear explosion pulse excitation based on the multi-conductor transmission line model and the time domain precise integration method, which has the following advantages: (1) the electromagnetic response characteristics of the rotating equipment of the power plant under the nuclear explosion pulse excitation can be predicted, especially the overvoltage characteristics, which can guide the design personnel to carry out research and development of protection methods and protection equipment, and propose relevant countermeasures for prevention, protection and recovery within a reasonable range; (2) the prediction algorithm has the advantages of less iteration steps and high calculation efficiency due to the wide stability condition of the time domain precise integration method; (3) the calculation accuracy of the time domain precise integration method does not change with the change of the time step, so the prediction algorithm can use the largest time step possible while maintaining the calculation accuracy, and has the advantages of high calculation accuracy; and (4) the algorithm has good universality.

[0112] In the application of the electromagnetic response prediction method for the rotating equipment of the power plant under the nuclear explosion pulse excitation provided by the present application, the excitation source can not be selected according to Figure 1The sequence of each step shown is executed, and the execution sequence of each specific step can be determined as required, and the present application does not limit this.

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

[0114] Figure 4 The electromagnetic response prediction device of the rotating equipment of the power plant under the nuclear explosion pulse excitation provided by the present application is shown in the schematic diagram, and the device 400 comprises:

[0115] The decomposition module 401 is used for decomposing the nuclear explosion pulse signal to obtain a plurality of frequency spectrum signals under different frequencies.

[0116] The generation module 402 is used for generating the electric field excitation source of the rotating equipment of the power plant under each frequency according to the frequency spectrum signal.

[0117] The construction module 403 is used for establishing the multi-conductor transmission line model of the rotating equipment of the power plant.

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

[0119] The determination module 405 is used for determining the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response under each frequency.

[0120] The specific limitation of the electromagnetic response prediction device of the rotating equipment of the power plant under the nuclear explosion pulse excitation can be referred to the limitation of the electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation in the above, and will not be repeated here. Each module in the above electromagnetic response prediction device of the rotating equipment of the power plant under the nuclear explosion pulse excitation can be realized by software, hardware and combination thereof in whole or in part. The above each module can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory in the computer device in the form of software, so as to call and execute the operation corresponding to each module by the processor.

[0121] The present application also provides a computer readable storage medium, which stores a computer program, and the computer program can be used to execute the electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation provided above. Figure 1

[0122] The present application also provides a computer readable storage medium, which stores a computer program, and the computer program can be used to execute the electromagnetic response prediction method of the rotating equipment of the power plant under the nuclear explosion pulse excitation provided above. Figure 5 The structure schematic diagram of the computer device is shown as Figure 5 ​​As shown, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory, and of course can also include other hardware required by the business. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs to implement the above Figure 1 Provided is a method for predicting the electromagnetic response of rotating equipment in a power plant under a nuclear explosion pulse excitation.

[0123] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer readable storage medium, and when executed, can include the processes of the above-mentioned embodiment methods. In the embodiments provided by the present application, any reference to a memory, storage, database, or other medium can include at least one of a non-volatile and volatile memory. The non-volatile memory can include a read-only memory (ROM), a magnetic tape, a floppy disk, a flash memory, or an optical memory. The volatile memory can include a random access memory (RAM) or an external cache memory. As an illustration but not as a limitation, the RAM can be in various forms, such as a static random access memory (SRAM) or a dynamic random access memory (DRAM).

[0124] The technical features of the above embodiments can be combined in any manner. To make the description concise, not all possible combinations of the technical features in the above embodiments are described, but as long as the combinations of the technical features do not contradict, they should be considered within the scope of the present application.

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 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 to obtain the electromagnetic response of the rotating equipment of the power plant at the frequency; determining the electromagnetic response of the rotating equipment of the power plant according to the electromagnetic response at each of the frequencies; 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 based on electrostatic field theory and state information of the power plant rotating equipment; the parameters include distributed capacitance parameters and distributed conductance parameters; The calculation formula of the capacitance parameter is: ； in, The first rotating equipment in the power plant i The conductor and j The capacitance parameter between the conductors is For the said i a conductor and the j The dielectric constant of the insulating medium between the conductors, For the said i The surface area per unit length of a conductor, For the said i a conductor and the j When there is a potential difference U between the two conductors, i The electrostatic field strength on the surface of a conductor; The distributed conductivity parameter is: ； in, For the said i a conductor and the j The conductance between the conductors, For the said i a conductor and the j The conductivity between conductors.

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; Spectrum signals having frequency amplitudes greater than or equal to a preset amplitude threshold are obtained from the spectrum signals, and spectrum signals having frequency amplitudes greater than or equal to the preset amplitude threshold are determined as spectrum signals at the multiple frequencies.

3. The method according to claim 1, characterized in that The calculation formula of the electric field excitation source is: ； in, for t The spectrum signal at time The electric field excitation source of frequency, The first i The amplitude corresponding to the frequency, The first i frequencies.

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

5. The method according to claim 1, wherein The method of using the electric field excitation source of the frequency 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 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; discretizing the spatial partial differential operators in the voltage partial differential equation and the current partial differential equation using a central difference format, and merging the discretized voltage partial differential equation and current partial differential equation 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.

6. 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, configured to generate an electric field excitation source for rotating equipment in a power plant at each frequency according to data of the spectrum signal; A construction module is configured to establish a multi-conductor transmission line model of the power plant rotating equipment; a parameter acquisition process of the multi-conductor transmission line model of the power plant rotating equipment includes: determining parameters of the multi-conductor transmission line model based on electrostatic field theory and state information of the power plant rotating equipment; the parameters include distributed capacitance parameters and distributed conductance parameters; The calculation formula of the capacitance parameter is: ； in, The first rotating equipment in the power plant i The conductor and j The capacitance parameter between the conductors is For the said i a conductor and the j The dielectric constant of the insulating medium between the conductors, For the said i The surface area per unit length of a conductor, For the said i a conductor and the j When there is a potential difference U between the two conductors, i The electrostatic field strength on the surface of a conductor; The distributed conductivity parameter is: ； in, For the said i a conductor and the j The conductance between the conductors, For the said i a conductor and the j The conductivity between conductors; a solving module, configured to solve the multi-conductor transmission line model at each frequency by taking the electric field excitation source of the frequency as input, and obtain 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.

7. A computer-readable storage medium, characterized in that The computer-readable 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 5 is implemented.

8. 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 the processor implements the method according to any one of claims 1 to 5 when executing the computer program.

Images