Method and device for fitting xenon oscillation test data of nuclear power plant

Abstract

The present disclosure belongs to the technical field of nuclear power and particularly relates to a method and device for fitting xenon oscillation test data of a nuclear power plant. The present disclosure finds a peak value or a valley value of an AO(t) curve changing with time, determines an oscillation period T through a time difference of a peak value in one period or a peak-valley value in half a period, selects an arbitrary time t, and obtains a corresponding time t1=-T / 4 or t2=T / 4 when AO(t)-AO(t+T / 2)=0, thereby determining AO(t) corresponding to t1 and t2 as an oscillation center C, a midpoint time between t1 and t2 being a point at t0=0 time, and further determining parameters of a free oscillation function. Thus, a reliable xenon oscillation test data processing method is provided, so that correct oscillation period, amplitude, stability index, oscillation center and other parameters are obtained, and it is verified whether a reactor core meets design requirements.

Description

Technical Field

[0001] This invention belongs to the field of nuclear power technology, specifically relating to a method and apparatus for fitting xenon oscillation test data in nuclear power plants. Background Technology

[0002] In related technologies, neutron flux is relatively high (e.g., neutron flux greater than 10¹³ neutrons / cm²). 2 In a large reactor core (0.sec), if the power in a certain local area decreases while the total power remains constant, the power in other areas will increase. However, the decrease in power in that local area leads to increased xenon toxicity, further reducing its power and iodine production, while the opposite occurs in other areas. Due to the initial decrease in iodine production and the delayed decay of iodine, the xenon toxicity in that area weakens over time, causing the power to increase, and iodine production also increases, while the opposite occurs in other areas. This process results in power oscillations within the reactor core. Because these oscillations are related to xenon toxicity, they are called xenon oscillations. Xenon oscillations can cause excessively high local power, leading to cladding burn-out and compromising the safe operation of the reactor. Furthermore, long-term local power fluctuations can cause material fatigue, threatening the safe operation of the reactor. Therefore, the design requires that the axial xenon oscillations of the reactor be convergent or suppressable by control rods. Xenon oscillation tests artificially induce axial xenon oscillations using control rods, and the changes in axial power deviation are measured to verify whether the reactor core meets the design requirements. This experiment is usually limited by time constraints, and it is generally impossible to complete the measurement of a full oscillation cycle; only half a cycle's measurement data can be obtained. Traditional data processing methods assume that the time from the peak to the trough is half a cycle T / 2, thus determining the point t=0. However, the entire xenon oscillation is decaying, and the determined t=0 time point and period have significant errors, resulting in a large discrepancy between the obtained tritium oscillation parameters and the actual values. Therefore, how to accurately obtain tritium oscillation parameters during xenon oscillation experiments has become an urgent problem to be solved. Summary of the Invention

[0003] To overcome the problems existing in related technologies, a method and apparatus for fitting xenon oscillation test data in nuclear power plants are provided.

[0004] According to one aspect of the present disclosure, a method for fitting xenon oscillation test data in a nuclear power plant is provided, the method comprising:

[0005] The in-pile power distribution measurement data were processed by power distribution processing software to obtain the in-pile axial power offset at multiple target times. ;

[0006] For each target time, based on the power Pr and axial power deviation ΔI of the off-core power range channel measured at that target time, the axial power offset corresponding to that target time is determined using Equation 1. ;

[0007] Formula 1

[0008] According to each target time corresponding The time corresponding to the target moment Perform correction to obtain the corrected result. ;

[0009] Determine each The curve changes over time, and the fitted curve yields the free oscillation function after the control rod is raised to the initial position, as shown in Equation 2.

[0010] Formula 2

[0011] In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; and C is the oscillation center.

[0012] Find the peak or trough value of the curve of AO(t) changing with time, and determine the oscillation period T by the time difference between the peak value of one cycle or the peak and trough value of half a cycle;

[0013] Choose any time t. When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. The AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint between t1 and t2 is the point at time t0=0.

[0014] When t0=0, AO(t0)=A+C, then we get A= AO(t0)-C;

[0015] When t3 = T / 2, AO(t3) = -A Bt +C, then we get .

[0016] In one possible implementation, the method further includes: based on the target time corresponding to and The difference between them is obtained by fitting. and The calculation formula between them, and the application of this calculation formula to all Perform correction to obtain the corrected result. .

[0017] According to another aspect of the present disclosure, a data fitting device for a nuclear power plant xenon oscillation test is provided, the device comprising:

[0018] The first analysis module is used to process the in-pile power distribution measurement data using power distribution processing software to obtain the in-pile axial power offset at multiple target times. ;

[0019] The second analysis module is used to determine the axial power offset corresponding to each target time based on the power Pr and axial power deviation ΔI of the off-pile power range channel measured at that target time, using Equation 1. ;

[0020] Formula 1

[0021] The correction module is used to adjust the values ​​corresponding to each target time. The time corresponding to the target moment Perform correction to obtain the corrected result. ;

[0022] The fitting module is used to determine each The curve changes over time, and the fitted curve yields the free oscillation function after the control rod is raised to the initial position, as shown in Equation 2.

[0023] Formula 2

[0024] In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; and C is the oscillation center.

[0025] The first determining module is used to find the peak or valley value of the AO(t) curve over time, and to determine the oscillation period T by the time difference between the peak value of one cycle or the peak and valley value of half a cycle.

[0026] The second determining module is used to select any time t. When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. The AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint time between t1 and t2 is the point at time t0=0.

[0027] The third determining module is used to determine that when t0=0, AO(t0)=A+C, and obtain A= AO(t0)-C;

[0028] The fourth determining module is used to determine AO(t3) = -A when t3 = T / 2. Bt +C, to get .

[0029] In one possible implementation, the correction module further includes: a correction submodule, configured to, based on the target time corresponding to... and The difference between them is obtained by fitting. and The calculation formula between them, and the application of this calculation formula to all Perform correction to obtain the corrected result. .

[0030] According to another aspect of the present disclosure, a data fitting device for a nuclear power plant xenon oscillation test is provided, the device comprising:

[0031] processor;

[0032] Memory used to store processor-executable instructions;

[0033] The processor is configured to execute the above-described method.

[0034] According to another aspect of the present disclosure, a non-volatile computer-readable storage medium is provided, on which computer program instructions are stored, which, when executed by a processor, implement the above-described method.

[0035] The beneficial effects of this disclosure are as follows: This disclosure finds the peak or trough value of the AO(t) change curve over time, and determines the oscillation period T by the time difference between the peak value of one cycle or the peak and trough value of half a cycle; by selecting any time t, when AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4 is obtained, thereby determining that the AO(t) corresponding to t1 and t2 is the oscillation center C, and the midpoint time between t1 and t2 is the point at time t0=0, thus determining the parameters of the free oscillation function. This provides a reliable method for processing xenon oscillation test data, thereby obtaining the correct oscillation period, amplitude, stability index and oscillation center parameters, thereby correctly verifying whether the reactor core meets the design requirements. Detailed Implementation

[0036] The present invention will be further described in detail below with reference to specific embodiments.

[0037] In one possible implementation, a method for fitting xenon oscillation test data from a nuclear power plant is provided. This method can be executed by a terminal device, which can be, for example, a server, desktop computer, laptop computer, etc. This disclosure does not limit the type of terminal device. The method may include:

[0038] Step 10: The in-pile power distribution measurement data is processed by power distribution processing software to obtain the in-pile axial power offset at multiple target times. ;

[0039] Step 20: For each target time, based on the power Pr and axial power deviation ΔI of the off-core power range channel measured at that target time, determine the corresponding axial power offset using Equation 1. ;

[0040] Formula 1

[0041] Step 30, based on the corresponding target time... The time corresponding to the target moment Perform correction to obtain the corrected result. ;

[0042] For example, it can be based on the corresponding target time. and The difference between them is obtained and The calculation formula between them, and the application of this calculation formula to all Perform correction to obtain the corrected result. .

[0043] For example, for each target moment, it is also possible to... and When the difference between them is greater than a preset threshold, and The median or average between them is used as the corrected value. .

[0044] Step 40, determine each The curve changing over time is fitted to obtain the free oscillation function AO(t) after the control rod is raised to the initial position, as shown in Equation 2:

[0045] Formula 2

[0046] In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; and C is the oscillation center.

[0047] Find the peak or trough value of the curve of AO(t) changing with time. The oscillation period T can be determined by the time difference between the peak value of one cycle or the peak and trough value of half a cycle.

[0048] Subtract AO(t) from the AO value after T / 2, i.e., AO(t)-AO(t+T / 2). When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint time between t1 and t2 is the point at t0=0.

[0049] When t0=0, AO(t0)=A+C, then we get A= AO(t0)-C;

[0050] When t3 = T / 2, AO(t3) = -A Bt +C, then we get .

[0051] In related technologies, nuclear power technicians rely on subjective judgment to fit the AO curve over time t, lacking clear data processing guidelines and exhibiting significant subjectivity and objectivity. This disclosure, through rigorous procedures and mathematical methods, accurately determines the time point t0=0, thereby correctly determining parameters such as the oscillation center C, amplitude A, and stability index B. Table 1 compares the parameters of the fitted curves calculated based on this disclosure, related technologies, and theoretical expectations. Referring to Table 1, the stability index B of this disclosure is closer to the theoretical expectation, resulting in a better fitting effect.

[0052] Table 1

[0053]

[0054] In one possible implementation, a data fitting device for a nuclear power plant xenon oscillation test is provided, the device comprising:

[0055] The first analysis module is used to process the in-pile power distribution measurement data using power distribution processing software to obtain the in-pile axial power offset at multiple target times. ;

[0056] The second analysis module is used to determine the axial power offset corresponding to each target time based on the power Pr and axial power deviation ΔI of the off-pile power range channel measured at that target time, using Equation 1. ;

[0057] Formula 1

[0058] The correction module is used to adjust the values ​​corresponding to each target time. The time corresponding to the target moment Perform correction to obtain the corrected result. ;

[0059] The fitting module is used to determine each The curve changes over time, and the fitted curve yields the free oscillation function after the control rod is raised to the initial position, as shown in Equation 2.

[0060] Formula 2

[0061] In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; and C is the oscillation center.

[0062] The first determining module is used to find the peak or valley value of the AO(t) curve over time, and to determine the oscillation period T by the time difference between the peak value of one cycle or the peak and valley value of half a cycle.

[0063] The second determining module is used to select any time t. When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. The AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint time between t1 and t2 is the point at time t0=0.

[0064] The third determining module is used to determine that when t0=0, AO(t0)=A+C, and obtain A= AO(t0)-C;

[0065] The fourth determining module is used to determine AO(t3) = -A when t3 = T / 2. Bt +C, to get .

[0066] In one possible implementation, the correction module further includes: a correction submodule, configured to, based on the target time corresponding to... and The difference between them is obtained by fitting. and The calculation formula between them, and the application of this calculation formula to all Perform correction to obtain the corrected result. .

[0067] The description of the above-mentioned device can be found in the description of the above-mentioned method, and will not be repeated here.

[0068] This disclosure can be a system, method, and / or computer program product. A computer program product may include a computer-readable storage medium having computer-readable program instructions loaded thereon for causing a processor to implement various aspects of this disclosure.

[0069] Computer-readable storage media can be tangible devices capable of holding and storing instructions for use by an instruction execution device. Computer-readable storage media can be, for example—but not limited to—electrical storage devices, magnetic storage devices, optical storage devices, electromagnetic storage devices, semiconductor storage devices, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of computer-readable storage media include: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), static random access memory (SRAM), portable compact disc read-only memory (CD-ROM), digital multifunction disc (DVD), memory sticks, floppy disks, mechanical encoding devices, such as punch cards or recessed protrusions storing instructions thereon, and any suitable combination of the foregoing. The computer-readable storage media used herein are not to be construed as transient signals themselves, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through waveguides or other transmission media (e.g., light pulses through fiber optic cables), or electrical signals transmitted through wires.

[0070] The computer-readable program instructions described herein can be downloaded from computer-readable storage media to various computing / processing devices, or downloaded via a network, such as the Internet, local area network, wide area network, and / or wireless network, to an external computer or external storage device. The network may include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers, and / or edge servers. A network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards them to the computer-readable storage media in the respective computing / processing device.

[0071] Computer program instructions used to perform the operations of this disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-dependent instructions, microcode, firmware instructions, status setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Smalltalk, C++, etc., and conventional procedural programming languages ​​such as the "C" language or similar programming languages. The computer-readable program instructions may execute entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or may be connected to an external computer (e.g., via the Internet using an Internet service provider). In some embodiments, electronic circuitry, such as programmable logic circuitry, field-programmable gate arrays (FPGAs), or programmable logic arrays (PLAs), is personalized by utilizing the status information of the computer-readable program instructions to implement various aspects of this disclosure.

[0072] The various embodiments of this disclosure have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen to best explain the principles, practical application, or improvement of the technology in the market, or to enable others skilled in the art to understand the embodiments disclosed herein.

Claims

1. A method for fitting xenon oscillation test data in a nuclear power plant, characterized in that, The method includes: The in-pile power distribution measurement data were processed by power distribution processing software to obtain the in-pile axial power offset at multiple target times. ; For each target time, based on the power Pr and axial power deviation ΔI of the off-core power range channel measured at that target time, the axial power offset corresponding to that target time is determined using Equation 1. ; Formula One; According to each target time corresponding Axial power offset corresponding to the target at that moment The correction is performed to obtain the axial power offset corresponding to the target time after correction. The corrected axial power offset corresponds to the target time. The correction process is as follows: According to each target time corresponding and The difference between them is obtained by fitting. and The calculation formula between them, and the application of this calculation formula to all The correction is performed to obtain the axial power offset corresponding to the target time after correction. ; Determine the axial power offset corresponding to each corrected target time. The curve changes over time, and the fitted curve yields the free oscillation function after the control rod is raised to the initial position, as shown in Equation 2. Formula 2; In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; C is the oscillation center; and AO(t) is the free oscillation function. Find the peak or trough value of the curve of AO(t) changing with time, and determine the oscillation period T by the time difference between the peak value of one cycle or the peak and trough value of half a cycle; Choose any time t. When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. The AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint between t1 and t2 is the point at time t0=0. When t0=0, AO(t0)=A+C, then we get A= AO(t0)-C; When t3 = T / 2, AO(t3) = -A Bt +C, then we get .

2. A data fitting device for xenon oscillation tests in nuclear power plants, characterized in that, The device includes: The first analysis module is used to process the in-pile power distribution measurement data using power distribution processing software to obtain the in-pile axial power offset at multiple target times. ; The second analysis module is used to determine the axial power offset corresponding to each target time based on the power Pr and axial power deviation ΔI of the off-pile power range channel measured at that target time, using Equation 1. ; Formula 1; The correction module is used to adjust the values ​​corresponding to each target time. Axial power offset corresponding to the target at that moment The correction is performed to obtain the axial power offset corresponding to the target time after correction. The corrected axial power offset at the target time. The correction module further includes: a correction submodule, used to correct the error according to each target time. and The difference between them is obtained by fitting. and The calculation formula between them, and the application of this calculation formula to all The correction is performed to obtain the axial power offset corresponding to the target time after correction. ; The fitting module is used to determine the axial power offset corresponding to each corrected target time. The curve changes over time, and the fitted curve yields the free oscillation function after the control rod is raised to the initial position, as shown in Equation 2. Formula 2; In the formula, t is time; T is the oscillation period; A is the amplitude; B is the stability index; C is the oscillation center; and AO(t) is the free oscillation function. The first determining module is used to find the peak or valley value of the AO(t) curve over time, and to determine the oscillation period T by the time difference between the peak value of one cycle or the peak and valley value of half a cycle. The second determining module is used to select any time t. When AO(t)-AO(t+T / 2)=0, the corresponding time t1=-T / 4 or t2=T / 4. The AO(t) corresponding to t1 and t2 is the oscillation center C. The midpoint time between t1 and t2 is the point at time t0=0. The third determining module is used to determine that when t0=0, AO(t0)=A+C, and obtain A= AO(t0)-C; The fourth determining module is used to determine AO(t3) = -A when t3 = T / 2. Bt +C, to get .

3. A data fitting device for xenon oscillation tests in nuclear power plants, characterized in that, The device includes: processor; Memory used to store instructions executed by the processor; The processor is configured to execute the nuclear power plant xenon oscillation test data fitting method according to claim 1.

4. A non-volatile computer read-only storage medium storing computer program instructions thereon, characterized in that, When the computer program instructions are executed by the processor, they implement the data fitting method for a xenon oscillation test in a nuclear power plant as described in claim 1.

Images