An online measurement method and system for reactivity parameters based on a pulse source

By acquiring and processing pulse source measurement data in real time, calculating slow-range neutron background and dynamic parameters, and calculating reactivity values in real time, the data lag problem of traditional pulse source methods is solved, and the efficiency and safety of reactor tests are improved.

CN116189941BActive Publication Date: 2025-07-22NUCLEAR POWER INSTITUTE OF CHINA
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202211489818.5
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2022-11-25
Publication Date
2025-07-22
Estimated Expiration
2042-11-25

AI Technical Summary

Technical Problem

The offline measurement working mode of the traditional pulse source method leads to lag in the analysis and processing of reactive measurement data, affecting the progress and safety of the reactor test.

Method used

A reactive parameter online measurement method based on pulse source is provided. By obtaining measurement data in real time, processing a neutron counting array, calculating the background level of slow-out neutrons, dynamic parameters and instantaneous neutron amplitude, and calculating the reactive value in real time.

Benefits of technology

Real-time measurement of reactive parameters is realized, the problem of data analysis and processing and evaluation lag is solved, and the efficiency and safety of reactor tests are improved.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN116189941B_ABST
    Figure CN116189941B_ABST
Patent Text Reader

Abstract

An on-line measurement method and system for reactivity parameters based on a pulse source according to an embodiment of the present invention include: obtaining in real time the fundamental wave prompt neutron decay parameters at different rod positions of a control rod; generating a reference point for reactivity measurement of the control rod by fitting according to the obtained fundamental wave prompt neutron decay parameters at different rod positions of the control rod; calculating the subcriticality reactivity values of the control rod in different rod position states according to the obtained fundamental wave prompt neutron decay parameters at different rod positions of the control rod and the reference point for reactivity measurement; and the embodiment of the present invention solves the technical problem that the traditional pulse source method uses an off-line measurement working mode, resulting in a lag in the analysis, processing and evaluation of data.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to an online measurement method and system for reactivity parameters based on a pulse source. Background Art

[0002] With the development of computer technology and the improvement of data acquisition design performance, online measurement technology has become an important development field in modern metrology and testing, and online measurement has been introduced into various data acquisition measurement models. Existing online measurement modes for other reactivity measurement methods have been developed, while the online measurement method for reactivity parameters based on pulse source technology has not been reported. Due to the large differences in the principles and equipment of different reactivity measurement methods, the existing online measurement modes are not applicable to the pulse source method. The pulse source method is a reactor physics test method that uses a neutron tube to periodically inject neutrons into the reactor and obtains reactivity parameters by measuring the leakage neutrons. The neutron tube belongs to Class II radiation device, and the maximum neutron yield injected into the reactor can reach 2*10 8 c / s.

[0003] The traditional offline measurement working mode of the pulse source method is inefficient. Its working mode is: according to the test procedure, raise the control rod position, input the measurement parameters, start the neutron tube to emit pulsed neutrons for measurement, observe the measurement curve and save the data, then readjust the control rod position for the next measurement, repeat the above steps and save the measurement data until the entire test is completed and then perform data analysis. It can be seen that in the traditional offline measurement working mode, data acquisition and analysis are carried out separately, and the analysis result lags behind, resulting in the lag of the analysis, processing and evaluation of reactivity measurement data. Especially in the critical physics test carried out on the reactor, due to the limitations of reactor operation safety and test operation procedures, if repeated measurements occur due to incorrect data, it will not only be time-consuming and laborious, but also seriously affect the test progress. Summary of the Invention

[0004] To solve the technical problem that the traditional pulse source method using the offline measurement working mode leads to the lag of data analysis, processing and evaluation, the embodiments of the present invention provide an online measurement method and system for reactivity parameters based on a pulse source.

[0005] The embodiments of the present invention are implemented through the following technical solutions:

[0006] In a first aspect, the embodiments of the present invention provide an online measurement method for reactivity parameters based on a pulse source, including:

[0007] Real-time obtain the measurement data at the control rod position of the pulse source in the current measurement cycle;

[0008] Process the measurement data to obtain the delayed neutron background level in the current measurement cycle;

[0009] Process the measurement data to obtain the dynamic parameters of the current measurement cycle;

[0010] Calculate the prompt neutron amplitude of the pulse source within the current measurement cycle;

[0011] Calculate the real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude;

[0012] Further, calculating the real-time reactivity value based on the delayed neutron background level, dynamic parameters of the current measurement cycle, and prompt neutron amplitude includes:

[0013] Calculate the real-time reactivity value according to formula (7):

[0014]

[0015] where ρ is the real-time reactivity value of the pulse source in the current measurement cycle, Φ d is the delayed neutron background level of the pulse source in the current measurement cycle, α p is the dynamic parameter of the pulse source in the current measurement cycle, A p is the prompt neutron amplitude of the pulse source in the current measurement cycle.

[0016] Further, the measurement data includes a neutron count array.

[0017] Further, processing the measurement data to obtain the delayed neutron background level of the current measurement cycle includes:

[0018] Compare each channel count in the neutron count array with the background count discrimination threshold one by one according to the channel number, and average the channel counts in the neutron count array that are less than the background count discrimination threshold to obtain the delayed neutron background level of the pulse source in the current measurement cycle.

[0019] Further, processing the measurement data to obtain the dynamic parameters of the current measurement cycle includes:

[0020] After deducting all channel counts in the neutron count array that are less than the background count discrimination threshold according to the delayed neutron background level of the pulse source in the current measurement cycle, a new neutron count array is obtained;

[0021] Compare the maximum count value in the channel X where the maximum count value in the new neutron count array is located with the count values in the adjacent channels X + 1, X + 2,..., X + N1 of the channel where the maximum count value is located in turn. If the count values in the channels X, X + 1, X + 2,..., X + N1 decrease in turn, then the channel where the count value in the channel X + N1 is located is the high-order harmonic termination channel; X and N1 are both positive integers;

[0022] Compare each channel of the new neutron count array with the delayed neutron background level of the current measurement period of the pulse source one by one starting from the initial channel. If the count values of channels 1, 2, …, N2 and the delayed neutron background level increase sequentially, the channel where the count value of channel 1 is located is the last channel of the delayed neutron background; N2 is a positive integer.

[0023] Fit the dynamic parameters by the method of successive rejection according to the data between the last channel of the delayed neutron background and the termination channel of the high-order harmonics.

[0024] Further, process the measurement data to obtain the delayed neutron background level of the current measurement period, including:

[0025] Successively compare the sub-count array n i with the background count discrimination threshold n b , and take the average count of each channel in the tail range of the neutron count array n i < n b as the delayed neutron background level Φ of the current measurement period d ; where i is the channel number, i = 1, 2, 3…N3, and N3 is a positive integer greater than zero.

[0026] Further, process the measurement data to obtain the dynamic parameters of the current measurement period, including:

[0027] According to the delayed neutron background level Φ of the current measurement period d , subtract the count of each channel in the neutron count array n i that is less than the background count discrimination threshold n b to obtain a new neutron count array nc i ;

[0028] Search for the maximum count value nc i of the new neutron count array nc max , successively compare nc max , nc max+1 …nc max+i channel counts, and take the nc max in the adjacent N4 channels that satisfy nc max+1 > nc max+i …nc max+i as the termination channel of the high-order harmonics;

[0029] Search the array nc i sequentially from the initial channel, query the first channel count value that satisfies nc min < Φ d , successively compare nc min-i … < nc min-1 < nc min , and take the nc min-i … < nc min-1<nc min The adjacent M channels of min , confirm the last channel nc of the delayed neutron background min-i ;

[0030] Take nc max+i To nc min-i Channel data, and the dynamic parameters are obtained by fitting according to the channel-by-channel rejection method. Among them, min > i; max, min, and M are all positive integers greater than zero; i is the channel number, i = 0, 2, 3…N4, and N4 is a positive integer greater than zero.

[0031] Furthermore, calculate the prompt neutron amplitude of the pulse source within the current measurement period; including:

[0032] For nc max+i To nc min-i Fit the data between the channel data to obtain a fitting equation;

[0033] Back-calculate the fitting equation to obtain the prompt neutron amplitude A p .

[0034] Furthermore, it also includes: calculating the real-time reactivity value according to the delayed neutron background level, the dynamic parameters of the current measurement period, and the prompt neutron amplitude, including:

[0035] Calculate the current delayed neutron background level of the pulse source, the dynamic parameters of the current measurement period, and the prompt neutron amplitude period using formula (7), and perform error analysis on the calculation result. If the calculation result converges and is less than or equal to the critical error, the real-time reactivity value is obtained.

[0036] In a second aspect, an embodiment of the present invention provides an on-line measurement system for reactivity parameters, including:

[0037] An acquisition unit for real-time acquisition of measurement data at the control rod position of the pulse source in the current measurement period;

[0038] A first processing unit for processing the measurement data to obtain the delayed neutron background level of the current measurement period;

[0039] A second processing unit for processing the measurement data to obtain the dynamic parameters of the current measurement period;

[0040] A first calculation unit for calculating the prompt neutron amplitude of the pulse source within the current measurement period; and

[0041] A second calculation unit for calculating the real-time reactivity value according to the delayed neutron background level, the dynamic parameters, and the prompt neutron amplitude.

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

[0043] An online measurement method and system for reactivity parameters based on a pulse source according to an embodiment of the present invention obtain measurement data at the control rod position of the pulse source in the current measurement cycle in real time; process the measurement data to obtain the delayed neutron background level in the current measurement cycle; process the measurement data to obtain the dynamic parameters in the current measurement cycle; calculate the prompt neutron amplitude of the pulse source in the current measurement cycle; and calculate the real-time reactivity value according to the delayed neutron background level, dynamic parameters, and prompt neutron amplitude, solving the technical problem that the traditional pulse source method uses an offline measurement working mode, resulting in a lag in the analysis, processing, and evaluation of data. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts.

[0045] Figure 1 It is a schematic flow chart of an online measurement method for reactivity parameters based on a pulse source.

[0046] Figure 2 It is a schematic structural diagram of an online measurement system for reactivity parameters based on a pulse source. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0047] To make the objectives, technical solutions, and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the embodiments and the drawings. The illustrative embodiments and descriptions thereof of the present invention are only used to explain the present invention and are not intended to limit the present invention.

[0048] In the following description, a large number of specific details are set forth in order to provide a thorough understanding of the present invention. However, it is obvious to those of ordinary skill in the art that the present invention does not have to employ these specific details. In other embodiments, well-known structures, circuits, materials, or methods have not been specifically described to avoid obscuring the present invention.

[0049] Throughout the specification, references to "one embodiment", "an embodiment", "an example" or "an illustration" mean that a particular feature, structure, or characteristic described in connection with the embodiment or example is included in at least one embodiment of the present invention. Thus, the phrases "one embodiment", "an embodiment", "an example" or "an illustration" appearing throughout the specification do not necessarily all refer to the same embodiment or example. Additionally, the particular features, structures, or characteristics may be combined in any suitable combination and / or sub-combination in one or more embodiments or examples. Further, those of ordinary skill in the art will understand that the diagrams provided herein are for illustrative purposes only and are not necessarily drawn to scale. The term "and / or" used herein includes any and all combinations of one or more of the associated listed items.

[0050] In the description of the present invention, the orientation or positional relationship indicated by terms such as "front", "rear", "left", "right", "upper", "lower", "vertical", "horizontal", "high", "low", "inner", "outer", etc. is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the protection scope of the present invention.

[0051] Embodiment

[0052] To solve the technical problem that the traditional pulse source method uses an off-line measurement working mode, resulting in a lag in the analysis, processing, and evaluation of data, in a first aspect, an embodiment of the present invention provides an on-line measurement method for reactivity parameters based on a pulse source, as shown in Figure 1 and includes:

[0053] S1. Obtain in real time the measurement data at the control rod position in the current measurement cycle of the pulse source;

[0054] S2. Process the measurement data to obtain the delayed neutron background level in the current measurement cycle;

[0055] S3. Process the measurement data to obtain the dynamic parameters in the current measurement cycle;

[0056] S4. Calculate the prompt neutron amplitude of the pulse source in the current measurement cycle;

[0057] S5. Calculate the real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude.

[0058] Thus, in the embodiment of the present invention, by acquiring in real time the measurement data at the position of the control rod in the current measurement cycle of the pulse source; processing the measurement data to obtain the delayed neutron background level in the current measurement cycle; processing the measurement data to obtain the dynamic parameters in the current measurement cycle; calculating the prompt neutron amplitude of the pulse source within the current measurement cycle; and calculating the real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude, the technical problem that the traditional pulse source method uses an off-line measurement working mode, resulting in lag in the analysis, processing, and evaluation of data, is solved.

[0059] Further, calculating the real-time reactivity value based on the delayed neutron background level, the dynamic parameters of the current measurement cycle, and the prompt neutron amplitude includes:

[0060] Calculating the real-time reactivity value according to formula (7):

[0061]

[0062] where ρ is the real-time reactivity value of the pulse source in the current measurement cycle, Φ d is the delayed neutron background level of the pulse source in the current measurement cycle, α p is the dynamic parameter of the pulse source in the current measurement cycle, and A p is the prompt neutron amplitude of the pulse source in the current measurement cycle.

[0063] The derivation process of formula (7) is as follows:

[0064] Assume that the prompt neutrons and delayed neutrons in the reactor are represented by subscripts p and d respectively. Then the balance equations of prompt neutrons and delayed neutrons in the reactor can be expressed as:

[0065]

[0066]

[0067]

[0068]

[0069] In the formula, L is the neutron disappearance operator, F is the neutron generation operator, Xp and Xd are the prompt and delayed neutron fission terms respectively, Q is the external neutron source term, and Km and Cm are the delayed neutron fraction and precursor concentration respectively.

[0070] Combined with (3), perform time integration on equation (2)

[0071]

[0072] Substitute (4) into (5) and multiply both sides by the conjugate function Integrating gives the variations of prompt neutrons and delayed neutrons:

[0073]

[0074] Changing the upper limit of the above integral to the pulse interval time T gives the expression for reactivity in the subcritical state:

[0075]

[0076] Furthermore, the measurement data includes a neutron count array.

[0077] Furthermore, processing the measurement data to obtain the delayed neutron background level for the current measurement cycle; includes:

[0078] Comparing each channel count in the neutron count array with the background count discrimination threshold one by one according to the channel number, and averaging the channel counts in the neutron count array that are less than the background count discrimination threshold to obtain the delayed neutron background level for the current measurement cycle of the pulse source.

[0079] Furthermore, processing the measurement data to obtain the dynamic parameters for the current measurement cycle; includes:

[0080] After deducting all the channel counts in the neutron count array that are less than the background count discrimination threshold according to the delayed neutron background level for the current measurement cycle of the pulse source, a new neutron count array is obtained;

[0081] Comparing the maximum count value in the channel X where the maximum count value in the new neutron count array is located with the count values in the adjacent channels X+1, X+2,..., X+N1 of the channel where the maximum count value is located in sequence. If the count values in the channels X, X+1, X+2,..., X+N1 decrease in sequence, then the channel where the count value in the channel X+N1 is located is the high-order harmonic termination channel; X and N1 are both positive integers;

[0082] Comparing each channel in the new neutron count array with the delayed neutron background level for the current measurement cycle of the pulse source one by one according to the channel number. If the count values in the channels 1, 2,..., N2 and the delayed neutron background level increase in sequence, then the channel where the count value in the channel 1 is located is the last channel of the delayed neutron background; N2 is a positive integer;

[0083] The dynamic parameters are obtained by fitting using the data between the last channel of the delayed neutron background and the high-order harmonic termination channel by the channel-by-channel rejection method.

[0084] Furthermore, processing the measurement data to obtain the delayed neutron background level for the current measurement cycle; includes:

[0085] Successively comparing the neutron count array n i with the background count discrimination threshold n b , and the neutron count array ni <n b The average count value of each channel in the tail range is used as the delayed neutron background level Φ of the current measurement period d ; where i is the channel number, i = 1, 2, 3…N3, and N3 is a positive integer greater than zero.

[0086] Furthermore, process the measurement data to obtain the dynamic parameters of the current measurement period; including:

[0087] According to the delayed neutron background level Φ of the current measurement period d , subtract all the channel counts in the neutron count array n i that are less than the background count discrimination threshold n b to obtain a new neutron count array nc i ;

[0088] Search for the maximum count value nc i of the new neutron count array nc max , and successively compare nc max , nc max+1 …nc max+i channel counts, and take the nc max in the adjacent N4 channels that satisfy nc max+1 >nc max+i ...>nc max+i as the high - order harmonic termination channel;

[0089] Search the array nc i sequentially from the initial channel, query the first channel count value that satisfies nc min <Φ d , and successively compare nc min-i …<nc min-1 <nc min , and take the adjacent M channels that satisfy nc min-i …<nc min-1 <nc min to confirm the last channel nc min-i of the delayed neutron background;

[0090] Take the data from channel nc min+i to nc min-i and fit to obtain the dynamic parameters according to the channel - by - channel rejection method. Where min>i; max, min, and M are all positive integers greater than zero; i is the channel number, i = 0, 2, 3…N4, and N4 is a positive integer greater than zero.

[0091] Furthermore, calculate the prompt neutron amplitude of the pulse source within the current measurement period; including:

[0092] For nc max+i to nc min-iFitting the data between the channels of data to obtain a fitting equation;

[0093] Back-calculating the prompt neutron amplitude A from the fitting equation p .

[0094] Further, it also includes: calculating a real-time reactivity value based on the delayed neutron background level, the dynamic parameters of the current measurement period, and the prompt neutron amplitude, including:

[0095] Calculating the current delayed neutron background level of the pulse source, the dynamic parameters of the current measurement period, and the prompt neutron amplitude period using formula (7), and performing error analysis on the calculation results. If the calculation results converge and are less than or equal to the critical error, the real-time reactivity value is obtained.

[0096] The general steps of the online measurement method for the reactivity parameters of a pulse source include: First, analyzing the measurement data to obtain the delayed neutron background level Φ d , performing calculation and analysis on the measurement data to obtain the dynamic parameter α p , calculating the A p value in real time for each measurement period of the pulse source, and calculating the online measurement reactivity value according to equation (7).

[0097] Specifically, the online measurement method for the reactivity parameters of a pulse source includes:

[0098] Step 1. Delayed neutron background analysis and confirmation: Set the background count discrimination threshold nb, and search the neutron count array n i acquired by data collection by the algorithm, and compare the array n i with the threshold nb one by one. Take the average value of the counts of each channel in a certain range at the tail of the neutron count array n i <n b as the delayed neutron background value Φ d ;

[0099] Step 2. According to the delayed neutron background level Φ d calculated in Step 1, subtract the background count of the neutron count array n i to obtain a new neutron count array nc i ;

[0100] Step 3. Dynamic parameter α p Analysis and confirmation: Search for the maximum count value nc i of the array nc max , compare nc max , nc max+1 ...nc max+i channel counts one by one, and take the one that satisfies nc max >nc max+1 ...>nc max+iThe adjacent N channels (N is set according to experimental conditions), confirm nc max+i as the high-order harmonic termination channel;

[0101] Step Four: Delayed neutron background analysis and confirmation: Search the array nc sequentially from the initial channel i , query the first-channel count value that satisfies nc min < Φ d , compare nc min-i ... < nc min-1 < nc min in sequence, and take the adjacent M channels (M is set according to experimental conditions) that satisfy nc min-i ... < nc min-1 < nc min to confirm the last channel nc of the delayed neutron background; min-i ;

[0102] Step Five: Take the data from channel nc max+i to nc min-i , and according to the channel-by-channel rejection method, combine the algorithm criterion to fit and obtain the dynamic parameter α p ;

[0103] Step Six: Calculate the A p value in real time. For each measurement period of the pulse source, according to nc max+i determined in Step Three and nc min-i determined in Step Four, fit the data between nc max+i and nc min-i , and inversely deduce the A p value according to the fitting equation.

[0104] Step Seven: Analyze according to the calculation formula (7), set the error analysis criterion in the software parameters, and the algorithm automatically gives the error analysis result. When the error analysis result converges to a certain value, prompt whether to terminate the measurement, and calculate the subcritical reactivity value ρ i of this rod state.

[0105] Exemplarily, an online measurement method for reactivity parameters based on a pulse source includes

[0106] Step One: Delayed neutron background analysis and confirmation, set the background count discrimination threshold nb, search the neutron count array n i collected by the algorithm from the data acquisition, compare the array n i with the threshold nb successively, and take the average value of the counts of each channel in a certain range at the tail of the neutron count array n i < n b as the delayed neutron background level Φ d ;

[0107] Step 2: Based on the delayed neutron background level Φ calculated in Step 1 d , subtract the background count of the neutron count array n i to obtain a new neutron count array nc i ;

[0108] Step 3: Dynamic parameter α p Analysis and confirmation: Search for the maximum count value nc i of the array nc max , and successively compare nc max , nc max+1 ... nc max+i channel counts. Select 5 adjacent channels that satisfy nc max > nc max+1 … > nc max+i to confirm nc max+5 as the termination channel of the high-order harmonic;

[0109] Step 4: Delayed neutron background analysis and confirmation: Search the array nc i sequentially from the initial channel, query the first channel count value that satisfies nc min < Φ d , and successively compare nc min-i … < nc min-1 < nc min . Select 8 adjacent channels that satisfy nc min-1 … < nc min-1 < nc min to confirm the last channel nc min-8 of the delayed neutron background;

[0110] Step 5: Take the data from channel nc max+i to nc min-i , and based on the channel-by-channel rejection method, combine the algorithm criterion to fit and obtain the dynamic parameter α p ;

[0111] Step 6: Calculate the A p value in real time. For each measurement period of the pulse source, based on nc max+i determined in Step 3 and nc min-i determined in Step 4, fit the data between nc max+i and nc min-i , and inversely deduce the A p value according to the fitting equation.

[0112] Step 7: Analyze according to the calculation formula (7), set the error analysis criterion in the software parameters, and let the algorithm automatically give the error analysis result. When the error analysis result converges to meet a certain value, prompt whether to terminate the measurement, and calculate the subcritical reactivity value ρ i of this rod state.

[0113] In a second aspect, an embodiment of the present invention provides an online measurement system for reactivity parameters. Referring to Figure 2 as shown, it includes:

[0114] An acquisition unit, configured to acquire in real time the measurement data at the control rod position of the pulse source in the current measurement cycle;

[0115] A first processing unit, configured to process the measurement data to obtain the delayed neutron background level in the current measurement cycle;

[0116] A second processing unit, configured to process the measurement data to obtain the dynamic parameters in the current measurement cycle;

[0117] A first calculation unit, configured to calculate the prompt neutron amplitude of the pulse source within the current measurement cycle; and

[0118] A second calculation unit, configured to calculate a real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude.

[0119] The above specific embodiments further elaborate on the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. An on-line measurement method for reactivity parameters based on a pulse source, characterized in that, Including: Obtain the measurement data at the control rod position in the current measurement cycle of the pulse source in real time; Process the measurement data to obtain the delayed neutron background level in the current measurement cycle; Process the measurement data to obtain the dynamic parameters in the current measurement cycle; Calculate the prompt neutron amplitude of the pulse source within the current measurement cycle; Calculate the real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude; Process the measurement data to obtain the delayed neutron background level in the current measurement cycle; including: Sequential comparison of neutron count arrays With the background count discrimination threshold , the average count of each channel in the tail range of the neutron count array is used as the delayed neutron background level for the current measurement period ; where i is the channel number, i = 1, 2, 3... N3, and N3 is a positive integer greater than zero; Process the measurement data to obtain the dynamic parameters in the current measurement cycle; including: According to the delayed neutron background level of the current measurement cycle , subtract the counts of all channels in the neutron count array that are less than the background count discrimination threshold to obtain a new neutron count array ; Search for a new neutron counting array The maximum count value , and compare successively channel counts, and take the in the adjacent N4 channels that meet as the high-order harmonic termination channel; Search the array sequentially from the initial channel , query the first channel count value that satisfies , compare sequentially , take the adjacent M channels that satisfy , and confirm the last channel of the delayed neutron background ; Take to channel data, and obtain dynamic parameters by fitting according to the channel-by-channel rejection method; where min > i; max, min, and M are all positive integers greater than zero; i is the channel number, i = 0, 2, 3…N4, and N4 is a positive integer greater than zero; Calculate the real-time reactivity value based on the delayed neutron background level, dynamic parameters in the current measurement cycle, and prompt neutron amplitude, including: Calculate the real-time reactivity value according to formula (7): (7) Among them, is the real-time reactivity value of the current measurement cycle of the pulse source, is the delayed neutron background level of the current measurement cycle of the pulse source, is the dynamic parameter of the current measurement cycle of the pulse source, is the prompt neutron amplitude of the current measurement cycle of the pulse source.

2. The on-line measurement method of the reactivity parameter based on a pulse source according to claim 1, characterized in that The measurement data includes a neutron count array.

3. The online measurement method of reactivity parameters based on a pulse source according to claim 2, characterized in that, Calculate the prompt neutron amplitude of the pulse source within the current measurement cycle; Including: Pair to Fit the data between the data to obtain a fitting equation; Back-calculate the prompt neutron amplitude from the fitting equation .

4. The on-line measurement method of the reactivity parameter based on the pulse source according to claim 3, characterized in that, Further including: Calculate the real-time reactivity value based on the delayed neutron background level, dynamic parameters in the current measurement cycle, and prompt neutron amplitude, including: Calculate the current delayed neutron background level of the pulse source, the dynamic parameters in the current measurement cycle, and the prompt neutron amplitude using formula (7), and perform error analysis on the calculation result. If the calculation result converges and is less than or equal to the critical error, then obtain the real-time reactivity value.

5. An on-line measurement system for reactive parameters, characterized in that, Used to implement an online measurement method for reactivity parameters based on a pulse source according to any one of claims 1-4, including: An acquisition unit for obtaining the measurement data at the control rod position in the current measurement cycle of the pulse source in real time; A first processing unit for processing the measurement data to obtain the delayed neutron background level in the current measurement cycle; A second processing unit for processing the measurement data to obtain the dynamic parameters in the current measurement cycle; A first calculation unit for calculating the prompt neutron amplitude of the pulse source within the current measurement cycle; and A second calculation unit for calculating the real-time reactivity value based on the delayed neutron background level, dynamic parameters, and prompt neutron amplitude.

Citation Information

Patent Citations

  • Automated system for controlling neutron physical parameters of research nuclear plant

    RU2593389C1

  • System and method to determine reactivity

    US20220230769A1