Reactivity measurement method and device based on mobile detector and medium

Abstract

The invention discloses a reactivity measurement method and device based on a mobile detector and a medium, and the method comprises the steps: enabling the mobile detector to follow the height of a to-be-measured control rod according to the position of the to-be-measured control rod and move based on an optimal detector movement path, measuring a neutron signal, and obtaining an original measurement neutron signal; in the reactivity measurement process, correcting an original measurement neutron signal by adopting a correction algorithm to obtain a corrected neutron signal of the mobile detector; and according to the corrected neutron signal of the mobile detector, carrying out reactivity measurement calculation to obtain a final reactivity measurement result. The device comprises an original signal acquisition unit, a correction unit and a reactivity measurement unit. In the reactivity measurement process, according to the spatial effect change condition of the dynamic measurement process, the position of the detector is moved for spatial effect compensation, so that the detector is located at the position with the minimum spatial effect in the measurement process, the influence of the neutron signal spatial effect is reduced, and the reactivity measurement efficiency and accuracy are improved.

Description

Technical Field

[0001] The invention relates to the field of nuclear reactor reactivity measurement, and in particular to a reactivity measurement method, device and medium based on a mobile detector. Background Art

[0002] The reactivity measurement of nuclear reactor core is mainly aimed at the measurement processes such as control rod integral value, rod stuck subcriticality, shutdown depth, etc. during physical testing.

[0003] Common reactivity measurement methods for nuclear power plants include: boron adjustment method, rod replacement method and dynamic rod engraving method. Among them, the boron adjustment method measures the micro-integral value of the control rod by adjusting the boron concentration, but the measurement time is long and the cost is high; the rod replacement method does not require boron adjustment, takes a slightly shorter time, and can measure the integral value of the control rod; the dynamic rod engraving method inserts the control rod, combines calculation analysis with measured data, and can measure the integral value of the control rod more accurately in a shorter time. It is a method widely used in nuclear power plants now. However, the above methods are generally used for reactivity measurements within about 2000pcm, and have limitations for large reactivity measurements.

[0004] The commonly used method for measuring large reactivity is the rod drop method. The basic principle of the rod drop method is to instantly drop a control rod into a reactor in a critical state, measure the change of the neutron count rate over time, and calculate the reactivity to be measured by combining the neutron count rate before and after the rod drop. This method can achieve rapid measurement of large reactivity such as the control rod integral value, rod stuck subcriticality, and shutdown depth by dropping the rod. However, this method is based on a point reactor model, and the measurement results are easily affected by the spatial effect of the neutron flux, and there are errors in the large reactivity measurement process.

[0005] In mature commercial reactors or power reactors, the detectors are fixed in the measurement channel and cannot be moved. However, in the zero-power test device, the detectors can be flexibly arranged and can be moved by grabbing with a robotic arm. In the reactivity measurement process, the neutron signal is the key parameter for reactivity measurement. This quantity is affected by the spatial effect factor during the measurement process, and the measurement results are biased.

[0006] So far, no other unit has conducted research on mobile detector-based reactivity measurement methods or related patented technologies. Therefore, it is necessary to master the autonomous mobile detector-based reactivity measurement method based on the differences between domestic cores and foreign cores, so as to provide the ability to efficiently and accurately measure core reactivity. Summary of the invention

[0007] The technical problem to be solved by the present invention is that in the existing reactivity measurement process, the neutron signal is a key parameter for reactivity measurement. This quantity is affected by the spatial effect factor during the measurement process, and the measurement result has deviations and is not accurate. The purpose of the present invention is to provide a reactivity measurement method, device and medium based on a mobile detector. For dynamic reactivity measurement, according to the change of the control rod position, the detector position is moved based on the optimal detector moving route during the reactivity measurement process, the influence of the sub-signal spatial effect during the reactivity measurement process is reduced, the shortcomings of the existing measurement method are compensated, and the accuracy and efficiency of the reactor reactivity measurement are improved.

[0008] The present invention is achieved through the following technical solutions:

[0009] In a first aspect, the present invention provides a reactivity measurement method based on a mobile detector, the method comprising:

[0010] When the reactor reaches a critical steady state for a period of time, reactivity is introduced through the control rods; according to the position of the control rod to be tested, the mobile detector follows the height of the control rod to be tested and moves based on the optimal detector movement route, while measuring the neutron signal to obtain the original measured neutron signal;

[0011] In the reactivity measurement process, the original measured neutron signal is corrected by using a correction algorithm to obtain a corrected mobile detector neutron signal;

[0012] A reactivity measurement calculation is performed according to the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

[0013] Furthermore, the acquisition of the optimal detector movement route includes:

[0014] Based on three-dimensional spatiotemporal dynamics analysis, the actual measurement process is simulated to obtain the neutron flux shape function in the dynamic measurement process; the neutron flux shape function ψ i (0),ψ i (1)…ψ i (t) is the neutron flux shape function of the mobile detector at n measuring points in the movable area at a number of time points from 0 to t; where i = 1, 2, 3…n;

[0015] Calculating a signal correction factor for each measuring point according to the neutron flux shape function;

[0016] Among the n measurement points, an optimal route is selected for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t m In the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x mThe signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the selected route. m (t m ) has the smallest deviation relative to 1; where m and n are both positive integers.

[0017] Furthermore, the expression of the signal correction factor is:

[0018]

[0019] In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i (0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

[0020] Furthermore, the expression of the offset degree is:

[0021]

[0022] Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined and given according to the experimental uncertainty analysis and evaluation; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

[0023] Furthermore, a neutron flux shape function in a dynamic measurement process is obtained, including:

[0024] The Monte Carlo method analysis software based on a high-fidelity model is used to perform three-dimensional space-time dynamics calculations and analysis to obtain the neutron flux shape function;

[0025] Alternatively, the response function of the contribution of each position in the core to the neutron flux at the mobile detector outside the core can be calculated using analysis software based on the high straightness model; the neutron flux shape function can be calculated by combining the neutron flux distribution in the core obtained by the three-dimensional space-time dynamics analysis program.

[0026] Furthermore, the corrected moving detector neutron signal I(t i ) is:

[0027] I(t i )=I m (t i )C(t i )

[0028] C(t i )=C xi (t i )

[0029] Among them, C(t i ) is t i The signal correction factor of each measuring point at the time corresponds to the correction factor of the corresponding measuring point at the corresponding time on the selected route; I m (t) is the original measured neutron signal.

[0030] Further, the final reactivity measurement result ρ is:

[0031]

[0032] Where Λ is the neutron generation time; β is the effective fraction of delayed neutrons; β i is the delayed neutron fraction of the i-th group; i is the i-th group delayed neutron decay constant, e is the base of natural logarithm, and the above parameters are calculated by physical calculation software or other calculation methods;

[0033] I(t) is the corrected neutron signal of the moving detector, and C is a further correction factor.

[0034] In a second aspect, the present invention further provides a reactivity measurement device based on a mobile detector, the device comprising:

[0035] The original signal acquisition unit is used to introduce reactivity through the control rod when the reactor reaches a critical steady state for a period of time; according to the position of the control rod to be tested, the mobile detector follows the height of the control rod to be tested and moves based on the optimal detector movement route, while measuring the neutron signal to obtain the original measured neutron signal;

[0036] A correction unit, used to correct the original measured neutron signal by using a correction algorithm during the reactivity measurement process to obtain a corrected mobile detector neutron signal;

[0037] The reactivity measurement unit is used to perform reactivity measurement calculation according to the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

[0038] Furthermore, the acquisition of the optimal detector movement route includes:

[0039] The computational simulation subunit is used to perform computational simulation on the actual measurement process based on three-dimensional spatiotemporal dynamics analysis to obtain the neutron flux shape function in the dynamic measurement process; the neutron flux shape function ψ i (0),ψ i (1)…ψ i (t) is the neutron flux shape function of the mobile detector at n measuring points in the movable area at a number of time points from 0 to t; where i = 1, 2, 3…n;

[0040] A correction factor calculation subunit, used to calculate the signal correction factor of each measurement point according to the neutron flux shape function;

[0041] The optimal route selection subunit is used to select an optimal route from n measurement points for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t m In the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x m The signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the selected route. m (t m ) has the smallest deviation relative to 1; where m and n are both positive integers.

[0042] Furthermore, the expression of the signal correction factor is:

[0043]

[0044] In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i (0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

[0045] Furthermore, the expression of the offset degree is:

[0046]

[0047] Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined and given according to the experimental uncertainty analysis and evaluation; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

[0048] Furthermore, a neutron flux shape function in a dynamic measurement process is obtained, including:

[0049] The Monte Carlo method analysis software based on a high-fidelity model is used to perform three-dimensional space-time dynamics calculations and analysis to obtain the neutron flux shape function;

[0050] Alternatively, the response function of the contribution of each position in the core to the neutron flux at the mobile detector outside the core can be calculated using analysis software based on the high straightness model; the neutron flux shape function can be calculated by combining the neutron flux distribution in the core obtained by the three-dimensional space-time dynamics analysis program.

[0051] In a third aspect, the present invention further provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned reactivity measurement method based on a mobile detector when executing the computer program.

[0052] In a fourth aspect, the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned reactivity measurement method based on a mobile detector is implemented.

[0053] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0054] 1. The present invention provides a reactivity measurement method, device and medium based on a mobile detector, and establishes a nuclear reactor reactivity measurement method that can reduce the influence of factors such as spatial effects during the reactivity measurement process. The present invention moves the detector position based on the optimal detector movement route during the reactivity measurement process according to the change of the control rod position, reduces the influence of the sub-signal spatial effect during the reactivity measurement process, and improves the accuracy and efficiency of the reactor reactivity measurement.

[0055] 2. The present invention provides a reactivity measurement method, device and medium based on a mobile detector, which adopts a reactivity measurement detector movement algorithm based on three-dimensional analysis to obtain the optimal detector movement route. During the reactivity measurement process, the position of the mobile detector is compensated for the spatial effect according to the change of the spatial effect of the dynamic measurement process, so that the detector is in the position with the minimum spatial effect during the measurement process, thereby reducing the influence of the spatial effect of the neutron signal. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings:

[0057] Figure 1 This is a flow chart of a reactivity measurement method based on a mobile detector of the present invention;

[0058] Figure 2 The flowchart of obtaining the optimal detector moving route of the present invention;

[0059] Figure 3 This is a structural block diagram of a reactivity measurement device based on a mobile detector of the present invention;

[0060] Figure 4 This is a structural block diagram of obtaining the optimal detector moving route of the present invention. DETAILED DESCRIPTION

[0061] Hereinafter, the term "include" or "may include" used in various embodiments of the present invention indicates the presence of the invented function, operation or element, and does not limit the addition of one or more functions, operations or elements. In addition, as used in various embodiments of the present invention, the terms "include", "have" and their cognates are intended only to indicate specific features, numbers, steps, operations, elements, components or combinations of the foregoing items, and should not be understood as first excluding the presence of one or more other features, numbers, steps, operations, elements, components or combinations of the foregoing items or the possibility of adding one or more features, numbers, steps, operations, elements, components or combinations of the foregoing items.

[0062] The terms used in various embodiments of the present invention are only used for the purpose of describing specific embodiments and are not intended to limit various embodiments of the present invention. As used herein, the singular form is intended to also include the plural form, unless the context clearly indicates otherwise. Unless otherwise limited, all terms used here (including technical terms and scientific terms) have the same meaning as the meaning generally understood by those of ordinary skill in the art to which the various embodiments of the present invention belong. The terms (such as the terms defined in the dictionary generally used) will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having an idealized meaning or an overly formal meaning, unless clearly defined in various embodiments of the present invention.

[0063] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments and drawings. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention.

[0064] In the existing reactivity measurement process, the neutron signal is the key parameter for reactivity measurement. This quantity is affected by the spatial effect factor during the measurement process, and the measurement results are biased and inaccurate. Therefore, in the reactivity measurement process, it is necessary to move the detector position to compensate for the spatial effect according to the change of the spatial effect in the dynamic measurement process, so that the detector is at the position with the minimum spatial effect during the measurement process, optimize the measurement data, reduce the influence of the spatial effect, and thus improve the measurement accuracy.

[0065] Therefore, the present invention designs a reactivity measurement method, device and medium based on a mobile detector. For dynamic reactivity measurement, according to the change of the control rod position, the detector position is moved based on the optimal detector moving route during the reactivity measurement process, thereby reducing the influence of sub-signal spatial effects during the reactivity measurement process, making up for the shortcomings of existing measurement methods, and improving the accuracy and efficiency of reactor reactivity measurement.

[0066] Example 1

[0067] like Figure 1 As shown, the present invention provides a reactivity measurement method based on a mobile detector, the method comprising:

[0068] Step 1: When the reactor reaches a critical steady state for a period of time, reactivity is introduced through the control rod; according to the position of the control rod to be tested, the mobile detector follows the height of the control rod to be tested and moves based on the optimal detector movement route, while measuring the neutron signal to obtain the original measured neutron signal I m (t);

[0069] The present invention adopts a neutron detector to collect measurement signals inside or outside the pile.

[0070] In this embodiment, Figure 2 As shown, the acquisition of the optimal detector moving route includes:

[0071] Step A, based on three-dimensional space-time dynamics analysis, the actual measurement process is simulated to obtain the neutron flux shape function ψ during the dynamic measurement process. n (t); the neutron flux shape function ψ n (t) is the neutron flux shape function of the mobile detector at n measuring points in the movable area at a number of time points from 0 to t; specifically, the Monte Carlo method analysis software based on the high-fidelity model can be used to perform three-dimensional space-time dynamics calculation analysis to obtain the neutron flux shape function ψ n (t); or use the analysis software based on the high straightness model to calculate the response function of each position in the core to the neutron flux contribution at the mobile detector outside the core; combine the neutron flux distribution in the core obtained by the three-dimensional space-time dynamics analysis program to calculate the neutron flux shape function ψ n (t).

[0072] Step B, according to the neutron flux shape function, based on n measuring points in the movable area of ​​the moving detector at a number of measuring moments within the time 0-t, calculate the signal correction factor of each measuring point; the expression of the signal correction factor is:

[0073]

[0074] In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i (0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

[0075] Step C, among the n measuring points, select an optimal route for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t mIn the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x m The signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the selected route. m (t m ) relative to "1" has the smallest deviation σ (the deviation can be measured by the sum of deviations or the sum of variances, etc.); where m and n are both positive integers.

[0076] Specifically, the expression of the offset degree is:

[0077]

[0078] Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined based on the experimental uncertainty analysis and evaluation or other analysis, or is processed by equal weight; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

[0079] The above technical solution, the present invention moves the mobile detector based on the above optimal detector moving route, and measures the neutron signal at the same time to obtain the original measured neutron signal I m (t).

[0080] The reactor is allowed to reach a critical steady state for a period of time, and reactivity is introduced through control rods. At the same time, the neutron signal is measured based on the optimal detector movement route obtained by the method for obtaining the optimal detector movement route, and the original measured neutron signal I is obtained. m During the measurement process, a small current amplifier is used to amplify the measurement signal of the mobile detector, which is converted into a voltage signal and then transmitted to the data acquisition card. The data is analyzed and processed by the processing terminal to obtain the actual measured neutron signal data.

[0081] Step 2, in the reactivity measurement process, using a correction algorithm to correct the original measured neutron signal to obtain a corrected mobile detector neutron signal;

[0082] In this embodiment, the corrected moving detector neutron signal I(t i ) is:

[0083] I(t i )=I m (t i )C(t i )

[0084] C(t i )=C xi (ti )

[0085] Among them, C(t i ) is t i The signal correction factor of each measuring point at the time corresponds to the correction factor of the corresponding measuring point at the corresponding time on the selected route; I m (t) is the original measured neutron signal.

[0086] The signal correction factor is analyzed according to the measurement results, and an appropriate signal correction factor is selected so that the tail of the corresponding reactivity measurement curve after all the reactivity is introduced is a horizontal straight line.

[0087] Step 3: Perform reactivity measurement calculation based on the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

[0088] In this embodiment, the final reactivity measurement result ρ is obtained by combining the corrected mobile detector neutron signal I(t) obtained in step 2 with the following or other reactivity measurement formulas:

[0089]

[0090] Where Λ is the neutron generation time; β is the effective fraction of delayed neutrons; β i is the delayed neutron fraction of the i-th group; i is the decay constant of the i-th group of delayed neutrons, e is the base of the natural logarithm, and the above parameters are calculated by physical calculation software or other calculation methods; I(t) is the corrected mobile detector neutron signal, and C is a further correction factor.

[0091] The beneficial effect of the present invention is that: in view of the influence of factors such as spatial effect in the reactivity measurement process, a nuclear reactor reactivity measurement method capable of reducing the above influence is established. The present invention adopts a reactivity measurement detector movement algorithm based on three-dimensional analysis, and in the reactivity measurement process, according to the change of spatial effect in the dynamic measurement process, the position of the moving detector is compensated for the spatial effect, so that the detector is at the position with the minimum spatial effect during the measurement process, reducing the influence of the spatial effect of the neutron signal, and improving the efficiency and accuracy of the reactivity measurement.

[0092] Example 2

[0093] like Figure 3 As shown, the difference between this embodiment and embodiment 1 is that the present invention further provides a reactivity measurement device based on a mobile detector, which has a one-to-one functional correspondence with a reactivity measurement method based on a mobile detector in embodiment 1; the device comprises:

[0094] The original signal acquisition unit is used to introduce reactivity through the control rod when the reactor reaches a critical steady state for a period of time; according to the position of the control rod to be tested, the mobile detector follows the height of the control rod to be tested and moves based on the optimal detector movement route, while measuring the neutron signal to obtain the original measured neutron signal;

[0095] A correction unit, used to correct the original measured neutron signal by using a correction algorithm during the reactivity measurement process to obtain a corrected mobile detector neutron signal;

[0096] The reactivity measurement unit is used to perform reactivity measurement calculation according to the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

[0097] As a further implementation, Figure 4 As shown, the acquisition of the optimal detector moving route includes:

[0098] The computational simulation subunit is used to perform computational simulation on the actual measurement process based on three-dimensional spatiotemporal dynamics analysis to obtain the neutron flux shape function in the dynamic measurement process; the neutron flux shape function ψ i (0),ψ i (1)…ψ i (t) is the neutron flux shape function of the mobile detector at n measuring points in the movable area at a number of time points from 0 to t; where i = 1, 2, 3…n;

[0099] A correction factor calculation subunit, used to calculate the signal correction factor of each measurement point according to the neutron flux shape function;

[0100] The optimal route selection subunit is used to select an optimal route from n measurement points for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t m In the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x m The signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the selected route. m (t m ) has the smallest deviation relative to 1; where m and n are both positive integers.

[0101] As a further implementation, the expression of the signal correction factor is:

[0102]

[0103] In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i(0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

[0104] As a further implementation, the expression of the offset degree is:

[0105]

[0106] Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined and given according to the experimental uncertainty analysis and evaluation; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

[0107] As a further implementation, a neutron flux shape function during a dynamic measurement is obtained, including:

[0108] The Monte Carlo method analysis software based on a high-fidelity model is used to perform three-dimensional space-time dynamics calculations and analysis to obtain the neutron flux shape function;

[0109] Alternatively, the response function of the contribution of each position in the core to the neutron flux at the mobile detector outside the core can be calculated using analysis software based on the high straightness model; the neutron flux shape function can be calculated by combining the neutron flux distribution in the core obtained by the three-dimensional space-time dynamics analysis program.

[0110] The execution process of each unit can be performed according to the process steps of the reactivity measurement method based on a mobile detector described in Example 1, and will not be described in detail in this embodiment.

[0111] The present invention aims at the influence of factors such as spatial effect in the reactivity measurement process, and establishes a nuclear reactor reactivity measurement method that can reduce the above influence. The reactivity measurement method based on mobile detector is adopted, and the position of the detector is moved during the reactivity measurement process according to the change of the control rod position, so as to reduce the influence of the spatial effect of the neutron signal and improve the reactivity measurement efficiency and accuracy.

[0112] At the same time, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned reactivity measurement method based on a mobile detector when executing the computer program.

[0113] At the same time, the present invention also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned reactivity measurement method based on a mobile detector is implemented.

[0114] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0115] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0116] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0117] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0118] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A reactivity measurement method based on a mobile detector, characterized in that: The method includes: According to the position of the control rod to be tested, the mobile detector follows the height of the control rod to be tested and moves based on the optimal detector moving route, while measuring the neutron signal to obtain the original measured neutron signal; In the reactivity measurement process, the original measured neutron signal is corrected by using a correction algorithm to obtain a corrected mobile detector neutron signal; A reactivity measurement calculation is performed according to the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

2. The method for measuring reactivity based on a mobile detector according to claim 1, characterized in that: The acquisition of the optimal detector moving route includes: Based on three-dimensional space-time dynamics, the actual measurement process is simulated to obtain the neutron flux shape function in the dynamic measurement process; the neutron flux shape function ψ i (0),ψ i (1)…ψ i (t) is the neutron flux shape function of the mobile detector at n measuring points in the movable area at a number of time points from 0 to t; where i = 1, 2, 3…n; Calculating a signal correction factor for each measuring point according to the neutron flux shape function; Among the n measurement points, an optimal route is selected for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t m In the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x m The signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the route. m (t m ) has the smallest deviation relative to 1; where m and n are both positive integers.

3. The method for measuring reactivity based on a mobile detector according to claim 2, characterized in that: The expression of the signal correction factor is: In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i (0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

4. The method for measuring reactivity based on a mobile detector according to claim 2, characterized in that: The expression of the offset degree is: Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined and given according to the experimental uncertainty analysis and evaluation; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

5. The method for measuring reactivity based on a mobile detector according to claim 2, characterized in that: Obtain neutron flux shape functions during dynamic measurements, including: The Monte Carlo method analysis software based on a high-fidelity model is used to perform three-dimensional space-time dynamics calculations and analysis to obtain the neutron flux shape function; Alternatively, analysis software based on a high straightness model is used to calculate the response function of the contribution of each position in the core to the neutron flux at the mobile detector outside the core; the neutron flux shape function is calculated by combining the neutron flux distribution in the core obtained by three-dimensional space-time dynamics analysis.

6. The method for measuring reactivity based on a mobile detector according to claim 1, characterized in that: The corrected moving detector neutron signal I(t i ) is: I(t i )=I m (t i )C(t i ) C(t i )=C xi (t i ) Among them, C(t i ) is t i The signal correction factor of each measuring point at the time corresponds to the correction factor of the corresponding measuring point at the corresponding time on the selected route; I m (t) is the original measured neutron signal.

7. The method for measuring reactivity based on a mobile detector according to claim 1, characterized in that: The final reactivity measurement result ρ is: Where, Λ is the progeny time; β is the effective fraction of delayed neutrons; β i is the delayed neutron share of the i-th group; λ i is the decay constant of the i-th group of delayed neutrons, e is the base of the natural logarithm, and the above parameters are calculated by physical calculation software; I(t) is the corrected neutron signal of the moving detector, and C is a further correction factor.

8. A reactivity measurement device based on a mobile detector, characterized in that: The device includes: The original signal acquisition unit is used to move the detector to follow the height of the control rod to be tested and move based on the optimal detector moving route according to the position of the control rod to be tested, and measure the neutron signal at the same time to obtain the original measured neutron signal; A correction unit, used to correct the original measured neutron signal by using a correction algorithm during the reactivity measurement process to obtain a corrected mobile detector neutron signal; The reactivity measurement unit is used to perform reactivity measurement calculation according to the corrected mobile detector neutron signal to obtain a final reactivity measurement result.

9. The reactivity measurement device based on a mobile detector according to claim 8, characterized in that: The acquisition of the optimal detector moving route includes: A computational simulation subunit is used to perform computational simulation on the actual measurement process based on three-dimensional space-time dynamics to obtain a neutron flux shape function in the dynamic measurement process; the neutron flux shape function is a neutron flux shape function of n measurement points of the mobile detector in the movable area at a number of time points from time 0 to time t; A correction factor calculation subunit, used to calculate the signal correction factor of each measurement point according to the neutron flux shape function; The optimal route selection subunit is used to select an optimal route from n measurement points for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t m In the measurement time, the measurement points at adjacent measurement time are x1, x2, x3, ..., x m The signal correction factors Cx1(t1), Cx2(t2),…, Cx are spatially adjacent and all measurement points on the route. m (t m ) has the smallest deviation relative to 1; where m and n are both positive integers.

10. The reactivity measurement device based on a mobile detector according to claim 9, characterized in that: The expression of the signal correction factor is: In the formula, C i (t) is the signal correction factor of the i-th measurement point, ψ i (0) is the neutron flux shape function at time 0, ψ i (t) is the neutron flux shape function at time t.

11. The reactivity measurement device based on a mobile detector according to claim 9, characterized in that: The expression of the offset degree is: Among them, w i is the weight factor of the correction signal at the i-th moment, which is determined and given according to the experimental uncertainty analysis and evaluation; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

12. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the reactivity measurement method based on a mobile detector as described in any one of claims 1 to 7 is implemented.

13. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method for measuring reactivity based on a mobile detector according to any one of claims 1 to 7 is implemented.

Images