A mobile probe based reactivity measurement method, apparatus and medium

Abstract

The application discloses a kind of reactive measurement method, device and medium based on mobile detector, comprising: according to the control rod position to be measured, mobile detector follows the control rod height to be measured and moves based on optimal detector moving route, while measuring neutron signal, obtain original measurement neutron signal;During the process of reactivity measurement, the original measurement neutron signal is corrected using correction algorithm, and the corrected mobile detector neutron signal is obtained;According to the corrected mobile detector neutron signal, the final reactivity measurement result is obtained by reactivity measurement calculation.The device includes original signal acquisition unit, correction unit and reactivity measurement unit.The application moves the position of mobile detector according to the change of spatial effect of dynamic measurement process during the process of reactivity measurement, and compensates the spatial effect, so that the detector is in the position of minimum spatial effect during the measurement process, reduces the influence of neutron signal spatial effect, and improves the efficiency and accuracy of reactivity measurement.

Description

TECHNICAL FIELD

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

[0002] The reactivity measurement of the nuclear reactor core is mainly aimed at the measurement process of the control rod integral value, the subcriticality of the control rod, and the shutdown depth in the physical test process.

[0003] The commonly used reactivity measurement methods in nuclear power plants include the boron adjustment method, the rod exchange method, and the dynamic rod marking method. Among them, the boron adjustment method measures the control rod micro and integral value by adjusting the boron concentration, but the measurement time is long and the cost is high; the rod exchange method does not need to adjust the boron and takes slightly less time, and can measure the control rod integral value; the dynamic rod marking method can accurately measure the control rod integral value in a short time by inserting the control rod, combining the calculation analysis and the measured data, and is the method widely used in nuclear power plants at present. However, the above methods are generally used for reactivity measurement within about 2000 pcm, and have limitations for large reactivity measurement.

[0004] The commonly used large reactivity measurement method is the falling rod method. The basic principle of the falling rod method is that the control rod is instantly dropped into the reactor in a critical state, and the change of the neutron count rate with time is measured, and the measured reactivity is calculated by combining the neutron count rates before and after the falling rod. This method can realize the rapid measurement of the control rod integral value, the subcriticality of the control rod, and the shutdown depth by falling the rod, but this method is based on the point reactor model, and the measurement result is easily affected by the spatial effect of the neutron flux, and there is an error in the large reactivity measurement process.

[0005] In mature commercial reactors or power reactors, the detector is fixed in the measurement hole and cannot be moved. However, on the zero-power test device, the detector can be flexibly arranged and can be moved by a mechanical arm. In the reactivity measurement process, the neutron signal is a key parameter for reactivity measurement, and this quantity is affected by the spatial effect factor during the measurement process, and the measurement result has a deviation.

[0006] Currently, there is no report on the research of the reactivity measurement method based on the mobile detector or related patent technology by other units. Therefore, it is necessary to master the self-owned reactivity measurement method based on the mobile detector according to the differences between the self-owned core and the foreign core, and to provide the ability for efficient and accurate measurement of the core reactivity. SUMMARY

[0007] The technical problem solved by the present application is that in the existing reactivity measurement process, the neutron signal is a key parameter for reactivity measurement, which is affected by the spatial effect factor in the measurement process, and the measurement result has deviation and is not accurate enough.

[0008] The present application is implemented by the following technical solutions:

[0009] In the first aspect, the present application provides a reactivity measurement method based on a moving detector, which comprises:

[0010] When the reactor reaches a critical steady state for a period of time, reactivity is introduced by a control rod; according to the position of the control rod to be measured, the detector is moved to follow the height of the control rod to be measured and moved based on an optimal detector moving route, while measuring the neutron signal to obtain the original measured neutron signal;

[0011] In the reactivity measurement process, a correction algorithm is used to correct the original measured neutron signal to obtain a corrected moving detector neutron signal;

[0012] According to the corrected moving detector neutron signal, reactivity measurement calculation is performed to obtain the final reactivity measurement result.

[0013] Further, the acquisition of the optimal detector moving route comprises:

[0014] Based on three-dimensional space-time dynamics analysis, the actual measurement process is calculated and 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 moving detector at n measurement points in the movable area at several time points from 0-t; wherein i=1, 2, 3…n;

[0015] According to the neutron flux shape function, the signal correction factor of each measurement point is calculated;

[0016] Among the n measurement points, an optimal route is selected for the moving detector to move and measure; the optimal route is such that at t1, t2, t3, …, t m , the measurement points x1, x2, x3, …, x mSpatially adjacent, and selected on the route of all measurement point signal correction factor Cx1(t1), Cx2(t2), …, Cx m (t m ) relative to 1 the offset degree is minimum; wherein, m, n are positive integers.

[0017] Further, 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] Further, the expression of the offset degree is:

[0021]

[0022] Wherein, w i is the weight factor of the i-th time correction signal, which is determined according to the test uncertainty analysis and evaluation analysis; C xi is the signal correction factor of the x i th measurement point, t i is the i-th measurement time.

[0023] Further, the neutron flux shape function in the dynamic measurement process is obtained, comprising:

[0024] Using the Monte Carlo method analysis software based on high fidelity model to carry out three-dimensional space-time dynamic calculation analysis, the neutron flux shape function is obtained;

[0025] Or using the analysis software based on high fidelity model to calculate the response function of the contribution of each position in the core to the neutron flux at the outside moving detector; combined with the neutron flux distribution in the core obtained by the three-dimensional space-time dynamic analysis program, the neutron flux shape function is calculated.

[0026] Further, the expression of 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] wherein C(t i ) is a signal correction factor of each measuring point at time t i , corresponding to the correction factor of the selected route at the corresponding time and corresponding measuring point; I m (t) is the original measured neutron signal.

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

[0031]

[0032] wherein Λ is the neutron generation time; β is the effective share of delayed neutrons; β i is the i-th group delayed neutron share; λ 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 mobile detector neutron signal, and C is a further correction factor.

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

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

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

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

[0038] Further, the acquisition of the optimal detector moving route comprises:

[0039] a calculation simulation subunit, configured to perform calculation simulation on an actual measurement process based on three-dimensional space-time dynamics analysis to obtain a neutron flux shape function in a 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 a movable area at a plurality of time points from 0 to t; wherein i=1, 2, 3…n.

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

[0041] an optimal route selection subunit configured to select an optimal route from the n measurement points for the mobile detector to move and measure; the optimal route is such that the offset degree of the signal correction factors Cx1(t1), Cx2(t2), …, Cx m n of adjacent measurement points is the smallest; wherein m and n are positive integers. m m m

[0042] Further, 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] Further, the expression of the offset degree is:

[0046]

[0047] wherein w i is the weight factor of the i-th time correction signal, which is determined according to the test uncertainty analysis evaluation analysis; C xi is the signal correction factor of the x i -th measurement point, t i is the i-th measurement time.

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

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

[0050] or an analysis software based on a high-fidelity model is used to calculate the response function of each position in the core to the neutron flux at the mobile detector outside the core; combined with the neutron flux distribution in the core obtained by the three-dimensional space-time dynamic analysis program, the neutron flux shape function is calculated.

[0051] ​​​In a third aspect, the present application 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 mobile detector-based reactivity measurement method when executing the computer program.

[0052] In a fourth aspect, the present application further provides a computer-readable storage medium, which stores a computer program, wherein the computer program is executed by a processor to implement the above-mentioned mobile detector-based reactivity measurement method.

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

[0054] 1. The mobile detector-based reactivity measurement method, device and medium provided by the present application can reduce the influence of factors such as spatial effects in the reactivity measurement process.

[0055] 2. The mobile detector-based reactivity measurement method, device and medium provided by the present application use a three-dimensional analysis-based reactivity measurement detector movement algorithm to obtain an optimal detector movement route, and move the detector position according to the spatial effect change of the dynamic measurement process in the reactivity measurement process to compensate for the spatial effect, so that the detector is in a position with the minimum spatial effect during the measurement process, thereby reducing the influence of the neutron signal spatial effect. BRIEF DESCRIPTION OF DRAWINGS

[0056] The accompanying drawings, which are included to provide a further understanding of the embodiments of the present application and constitute a part of the application, illustrate embodiments of the present application and do not limit the present application. In the drawings:

[0057] Figure 1 The flowchart of the mobile detector-based reactivity measurement method of the present application;

[0058] Figure 2 The flowchart of the acquisition of the optimal detector movement route of the present application;

[0059] Figure 3 The structural block diagram of the mobile detector-based reactivity measurement device of the present application;

[0060] Figure 4 The structural block diagram of the acquisition of the optimal detector movement route of the present application. DETAILED DESCRIPTION

[0061] Hereinafter, the terms "include" or "may include" used in various embodiments of the present invention indicate the presence of the invented function, operation or element, and do 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, and should not be understood as excluding the presence of one or more other features, numbers, steps, operations, elements, components or combinations of the foregoing or the possibility of adding one or more features, numbers, steps, operations, elements, components or combinations of the foregoing.

[0062] The terms used in various embodiments of the present invention are only used to describe the purpose of specific embodiments and are not intended to limit the 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 those of ordinary skill in the art generally understood by the various embodiments of the present invention. The terms (such as those defined in generally used dictionaries) will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having idealized meaning or too 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 examples and drawings. The exemplary embodiments of the present invention and their descriptions are only used to explain the present invention and are not intended to limit the present invention.

[0064] In existing reactivity measurements, the neutron signal is a key parameter. However, this quantity is affected by spatial effects during the measurement process, resulting in biased and inaccurate results. Therefore, during reactivity measurements, it is necessary to compensate for spatial effects by moving the detector position according to the dynamic changes in spatial effects during the measurement process. This ensures that the detector is positioned to minimize spatial effects during the measurement process, optimizes the measurement data, reduces the impact of spatial effects, and thus improves measurement accuracy.

[0065] Therefore, the present invention designs a reactivity measurement method, device and medium based on a mobile detector. For dynamic reactivity measurement, the detector position is moved based on the optimal detector movement route during the reactivity measurement process according to the changes in the control rod position, thereby reducing the influence of the sub-signal spatial effect during the reactivity measurement process, making up for the shortcomings of the 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 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 I m (t);

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

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

[0071] Step A: Based on the 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 measurement points in the movable area at several time points from 0 to t; specifically, the neutron flux shape function ψ can be obtained by performing three-dimensional spatiotemporal dynamics calculation analysis based on the Monte Carlo method analysis software based on a high-fidelity model. n (t); or use analysis software based on the high straightness model to calculate the response function of each position in the core to the neutron flux 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, based on the neutron flux shape function, calculates the signal correction factor of each measurement point at n measurement points in the movable area of ​​the detector at a number of measurement moments within the time 0-t; the expression of the signal correction factor is:

[0073]

[0074] Where 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, selecting an optimal route from the n measurement points for the mobile detector to perform mobile measurement; the optimal route is such that t1, t2, t3, ..., t mThe measurement points x1, x2, x3, …, x of adjacent measurement time points in a measurement time point m The signal correction factors Cx1(t1), Cx2(t2), …, Cx m (t m ) of all measurement points on the selected route are adjacent in space.

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

[0077]

[0078] Wherein, w i is the weight factor of the i-th time correction signal, which is given according to the test uncertainty analysis evaluation or other analysis, or equal weight processing is adopted; C xi is the signal correction factor of the x i th measurement point, t i is the i-th measurement time.

[0079] The above technical scheme, the present application is based on the above optimal detector moving route to move the moving detector, while measuring the neutron signal, obtaining the original measurement neutron signal I m (t)。

[0080] Make the reactor reach the critical steady state for a period of time, introduce the reactivity through the control rod, at the same time, based on the optimal detector moving route obtained by the method for obtaining the optimal detector moving route, measure the neutron signal, obtain the original measurement neutron signal I m (t). During the measurement process, a small current amplifier is used to amplify the measurement signal of the moving detector, and the voltage signal is converted and transmitted to the data acquisition card, and the data is analyzed and processed through the processing terminal, and the actual measured neutron signal data is obtained.

[0081] Step 2, in the process of reactivity measurement, the original measurement neutron signal is corrected by using the correction algorithm, and the corrected moving detector neutron signal is obtained.

[0082] In this embodiment, the expression of 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] wherein, C(t i ) is a signal correction factor of each measuring point at t i time, corresponding to the correction factor of the selected route at the corresponding time and the corresponding measuring point; I m (t) is the original measured neutron signal.

[0086] The signal correction factor is analyzed according to the measurement result, 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, according to the corrected mobile probe neutron signal, the reactivity measurement calculation is carried out, and the final reactivity measurement result is obtained.

[0088] In this embodiment, the corrected mobile probe neutron signal I(t) obtained in step 2 is combined with the following or other reactivity measurement formula, and the final reactivity measurement result p is:

[0089]

[0090] In the formula, Λ is the neutron generation time; β is the effective share of delayed neutrons; β i is the i-th group delayed neutron share; λ 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; I(t) is the corrected mobile probe neutron signal, and C is a further correction factor.

[0091] The beneficial effects of the present application are that a nuclear reactor reactivity measurement method capable of reducing the influence of the above factors is established for the influence of the spatial effect in the reactivity measurement process. The present application adopts a reactivity measurement probe moving algorithm based on three-dimensional analysis, moves the probe position according to the spatial effect change of the dynamic measurement process during the reactivity measurement process, and compensates the spatial effect, so that the probe is in the position with the minimum spatial effect during the measurement process, reduces the influence of the neutron signal spatial effect, and improves the reactivity measurement efficiency and accuracy.

[0092] Embodiment 2

[0093] As shown in Figure 3 , the difference between this embodiment and embodiment 1 is that the present application further provides a reactivity measurement device based on a mobile probe, which is one-to-one corresponding to the function of the reactivity measurement method based on a mobile probe in embodiment 1; the device comprises:

[0094] ​The original signal acquisition unit is configured to introduce reactivity by controlling a control rod when a reactor reaches a critical steady state for a period of time; move a detector to follow a height of the control rod to be measured according to a position of the control rod to be measured and move the detector based on an optimal detector moving route while measuring a neutron signal to obtain an original measured neutron signal;

[0095] The correction unit is configured to correct the original measured neutron signal by using a correction algorithm during reactivity measurement to obtain a corrected moving detector neutron signal;

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

[0097] As a further implementation, as shown in Figure 4 The optimal detector moving route is obtained by:

[0098] The calculation simulation subunit is configured to perform calculation simulation on an actual measurement process based on three-dimensional space-time dynamics analysis to obtain a neutron flux shape function in a dynamic measurement process; the neutron flux shape function ψ i (0), ψ i (1), …, ψ i (t) is a neutron flux shape function of the moving detector at n measurement points in the movable area at a plurality of time points from 0-t; wherein i=1, 2, 3…n.

[0099] The correction factor calculation subunit is configured to calculate a signal correction factor of each measurement point according to the neutron flux shape function.

[0100] The optimal route selection subunit is configured to select an optimal route for the moving detector to move and measure among the n measurement points; the optimal route is such that the measurement points x1, x2, x3, …, xn at adjacent measurement times t1, t2, t3, …, tn in the 0-t time period are spatially adjacent, and the signal correction factors Cx1(t1), Cx2(t2), …, Cxn(tn) of all measurement points on the selected route have the smallest offset degree relative to 1; wherein m and n are positive integers. m m m m

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

[0102]

[0103] In the expression, C i (t) is the signal correction factor of the i-th measurement point, and ψ 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] wherein w i is the weight factor of the i-th moment correction signal, which is determined according to the test uncertainty analysis evaluation analysis; C xi is the signal correction factor of the x i th measurement point, t i is the i-th measurement moment.

[0107] As a further implementation, the neutron flux shape function in the dynamic measurement process is obtained, comprising:

[0108] The Monte Carlo method analysis software based on the high-fidelity model is used for three-dimensional space-time dynamic calculation analysis to obtain the neutron flux shape function.

[0109] Or the analysis software based on the high-fidelity model is used to calculate the response function of the contribution of each position in the core to the neutron flux at the out-of-core mobile detector; combined with the in-core neutron flux distribution obtained by the three-dimensional space-time dynamic analysis program, the neutron flux shape function is calculated.

[0110] Wherein, the execution process of each unit is executed according to the flow steps of the reaction measurement method based on the mobile detector described in embodiment 1, which will not be described one by one in this embodiment.

[0111] The present application establishes a nuclear reactor reactivity measurement method which can reduce the influence of the above factors such as spatial effect in the reactivity measurement process. The reactivity measurement method based on the mobile detector is used to move the detector position according to the control rod position change in the reactivity measurement process, reduce the influence of the neutron signal spatial effect, and improve the reactivity measurement efficiency and accuracy.

[0112] Meanwhile, the present application further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the above-mentioned reactivity measurement method based on the mobile detector when executing the computer program.

[0113] Meanwhile, the present application further provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the above-mentioned reactivity measurement method based on the mobile detector.

[0114] Those skilled in the art will appreciate that embodiments of the application can be readily used as software, hardware, or a combination of software and hardware. In one

[0115] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks.

[0116] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks.

[0117] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks in the flowcharts and / or a combination thereof. These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks.

[0118] The above detailed description has shown, described, and pointed out the aspects of the application in sufficient detail, that others skilled in the art can follow the present application. It is understood that various modifications can be made to the procedures, systems, and computer program products described herein, and such modifications are meant to be included within the scope of the application. Therefore, the scope of the application should be determined by the following claims, rather than by the detailed description and proce

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 movement route, while measuring the neutron signal to obtain the original measured neutron signal; During the reactivity measurement process, a correction algorithm is used to correct the original measured neutron signal to obtain a corrected moving detector neutron signal; performing a reactivity measurement calculation based on the corrected mobile detector neutron signal to obtain a final reactivity measurement result; The acquisition of the optimal detector movement 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 measurement 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 measurement point based on the neutron flux shape function; Among 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 moments, the measurement points at adjacent measurement moments 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 are m (t m ) has the smallest deviation from 1; where m and n are both positive integers.

2. The reactivity measurement method based on a mobile detector according to claim 1, characterized in that: The expression of the signal correction factor is: Where 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.

3. The reactivity measurement method based on a mobile detector according to claim 1, 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 based on 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.

4. The reactivity measurement method based on a mobile detector according to claim 1, characterized in that: Obtain neutron flux shape functions during dynamic measurements, including: The neutron flux shape function was obtained by performing three-dimensional spatiotemporal dynamics calculations and analysis using Monte Carlo analysis software based on a high-fidelity model. Alternatively, analysis software based on a high straightness model can be 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; combined with the neutron flux distribution in the core obtained by three-dimensional space-time dynamics analysis, the neutron flux shape function can be calculated.

5. The reactivity measurement method 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, corresponding 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; C xi For the xth i Signal correction factor for each measuring point, t i is the i-th measurement moment.

6. The reactivity measurement method based on a mobile detector according to claim 1, characterized in that: The final reactivity measurement result ρ is: Where, Λ is the progeny generation time; β is the effective fraction of delayed neutrons; β i is the delayed neutron share of group i; λ i is the decay constant of the delayed neutrons of group i, and e is the base of the natural logarithm. 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.

7. A reactivity measurement device based on a mobile detector, characterized in that: The device includes: an original signal acquisition unit, configured to move the detector to follow the height of the control rod to be tested and move based on an optimal detector movement route according to the position of the control rod to be tested, while measuring the neutron signal to obtain an original measured neutron signal; a correction unit, configured to correct the original measured neutron signal using a correction algorithm during the reactivity measurement process to obtain a corrected moving detector neutron signal; a reactivity measurement unit, configured to perform reactivity measurement calculations based on the corrected mobile detector neutron signal to obtain a final reactivity measurement result; The acquisition of the optimal detector movement route includes: A computational simulation subunit is used to computationally simulate the actual measurement process based on three-dimensional space-time dynamics to obtain a neutron flux shape function during the dynamic measurement process; the neutron flux shape function is the neutron flux shape function of the mobile detector at n measurement points in the movable area at a number of time points from time 0 to time t; a correction factor calculation subunit, configured to calculate a signal correction factor for each measurement point based on the neutron flux shape function; The optimal route selection subunit is used to select an optimal route for the mobile detector to perform mobile measurement among n measurement points; the optimal route is such that t1, t2, t3, ..., t m In the measurement moments, the measurement points at adjacent measurement moments 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 are m (t m ) has the smallest deviation from 1; where m and n are both positive integers.

8. The reactivity measurement device based on a mobile detector according to claim 7, characterized in that: The expression of the signal correction factor is: Where 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.

9. The reactivity measurement device based on a mobile detector according to claim 7, 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 based on 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.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the reactivity measurement method based on a mobile detector according to any one of claims 1 to 6 is implemented.

11. 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 6 is implemented.

Images