CN112687411A - Reactivity measurement method based on multi-detector relay signals - Google Patents
Reactivity measurement method based on multi-detector relay signals Download PDFInfo
- Publication number
- CN112687411A CN112687411A CN202011551550.4A CN202011551550A CN112687411A CN 112687411 A CN112687411 A CN 112687411A CN 202011551550 A CN202011551550 A CN 202011551550A CN 112687411 A CN112687411 A CN 112687411A
- Authority
- CN
- China
- Prior art keywords
- neutron
- reactivity
- sequence
- groups
- sequences
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a reactivity measurement method based on multi-detector relay signals, which comprises the following steps of S1: respectively setting proper and fixed small current amplification gears for the n neutron detectors, enabling the signal amplitudes of the detectors passing through the small current amplifiers to be mutually partially covered, and performing S2 reactivity control to obtain n groups of amplified neutron measurement signal sequences; s3, acquiring a sliding window: finding out an effective amplitude value crossing area and recording as a sliding window; s4, physical modeling: n groups of neutron flux value sequences; s4, calculating an evaluation factor: obtaining n groups of evaluation factor sequences; s5, obtaining a simulated gear shifting splicing point through optimal evaluation: the time point corresponding to the optimal evaluation value is regarded as a simulated gear shifting splicing point; s6, constructing a sequence required by reactivity calculation: constructing n groups of reactivity signal sequences to be detected; cutting off n groups of reactive signal sequences to be tested by taking the simulated gear shifting splicing points as nodes, recombining and sequencing to obtain sequences required by reactivity calculation; and S7, calculating the reactivity.
Description
Technical Field
The invention relates to the field of nuclear reactor reactivity measurement, in particular to a reactivity measurement method based on multi-detector relay signals.
Background
The reactivity measurement of the reactor core of the nuclear reactor mainly aims at the measurement processes of control rod micro, integral value, rod clamping subcritical degree, shutdown depth and the like in the physical test process, and has very important significance for calibrating the core nuclear design and providing reference basis for the reactor operation by measuring the reactivity.
The existing reactivity measurement methods are: a single neutron detector is adopted to shift and control a small current amplifier; reactivity calculations were performed on the obtained measurement signals.
Specifically, certain reactivity is introduced after the reactor reaches a critical state in the process, the neutron signal changes along with time in the measurement process, and the small current amplifier selects an amplification gear according to the amplitude of the neutron signal, so that the amplified neutron signal has enough linearity. However, when the gear is shifted, due to the electronic characteristics of circuit hardware, a part of amplified neutron signals generate certain fluctuation, and the measurement accuracy is reduced. During some negative reactivity measurements, the neutron signal may span 3-4 steps, and signal fluctuations from gear shifts will affect the reactivity measurement.
Disclosure of Invention
The invention aims to provide a reactivity measurement method based on multi-detector relay signals.
The invention is realized by the following technical scheme:
a reactivity measurement method based on multi-detector relay signals comprises the following steps:
s1, layout of neutron detectors: arranging n neutron detectors at a reactivity measuring station of a target reactor, and setting proper and fixed low-current amplification gears for the n neutron detectors respectively, so that signal amplitudes of the detectors passing through respective low-current amplifiers can be partially covered with each other and can be connected with each other, and a wider measuring range is realized;
s2, reactivity control: performing reactivity control on the target reactor from the control start time t0To the control end time tendContinuously monitoring by using n neutron detectors to obtain n groups of amplified neutron measurement signal sequences corresponding to the n neutron detectors one by one;
s3, acquiring a sliding window: observing the signal amplitude of each time point of each group of amplified neutron measurement signal sequences, finding an effective amplitude crossing region of one group of amplified neutron measurement signal sequences and the other group of amplified neutron measurement signal sequences, and recording a time segment corresponding to the effective amplitude crossing region as a sliding window;
s4, physical modeling: according to the physical parameters of the target reactor and the neutron detector in the step S1, three-dimensional physical modeling is carried out in a computer, and the three-dimensional physical modeling is calculated by using the built model to obtain: from the control start time t0To the control end time tendN groups of neutron flux value sequences corresponding to the n neutron detectors one by one;
s4, calculating an evaluation factor: calculating to obtain an evaluation factor E aiming at the neutron flux value corresponding to each moment in the n groups of neutron flux value sequencesi(t), obtaining n groups of evaluation factor sequences;
s5, obtaining a simulated gear shifting splicing point through optimal evaluation:
taking 1 time point from the sliding window as an estimated simulated gear shifting splicing point, truncating n groups of evaluation factor sequences by taking the estimated simulated gear shifting splicing point as a node, recombining and sequencing the truncated evaluation factor sequences of different segments according to time sequence to obtain an estimated evaluation factor sequence, and carrying out average value calculation on the estimated evaluation factor sequence to obtain an estimated value; repeating the process until all the time points in the sliding window are subjected to sliding calculation;
comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point;
s6, constructing a sequence required by reactivity calculation:
firstly, constructing n groups of reactivity signal sequences to be detected based on n groups of amplified neutron measurement signal sequences;
cutting off n groups of reactive signal sequences to be detected by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the cut-off n groups of reactive signal sequences to be detected in different sections from small to large according to the absolute amplitude of output signals to obtain a sequence required by reactivity calculation;
s7, reactivity calculation: reactivity calculations were performed based on the sequence required for reactivity calculations.
The invention aims to solve the problems that when the change range span of the neutron signal is large and gear shifting exists, the measurement accuracy is improved, the measurement efficiency can be improved while the wide-range measurement is realized, and the measurement accuracy can also be improved. The design principle of the invention is as follows:
adopt a plurality of neutron detectors, the gear that the undercurrent of every neutron detector of fixer launched, these gears are when the earlier stage is experimental, ensure that it can cover each other to set up, thereby realize the wide range, then, the signal synchronization with these neutron detectors is recorded simultaneously, adopt backstage computer to splice these signals, keep every effective data of surveying in its gear measuring range, like this to the system, shift gears the action in its measuring process, consequently, there is not signal fluctuation, thereby guarantee the measurement accuracy. All that is needed in the present invention is to find a suitable splice point. For each neutron detector, the measured signal is a continuous signal, and thus, it can be regarded as a signal sequence. The invention provides a processing method for searching an optimal splicing point through theoretical data, which has the following principle: performing physical modeling on physical parameters based on the target reactor and the neutron detector, for example, adopting Monte Carlo method physical calculation analysis software based on a high-fidelity model or other physical calculation analysis software; thus, theoretical neutron flux value data can be obtained, and the neutron flux value data and an actual neutron measurement signal are combined for evaluation; the best splice point is found. In order to simplify the calculation, the gear covering area between the neutron detectors is used as a sliding window for calculation, and the calculated evaluation values are analyzed and compared to find the optimal splicing point. And then splicing a plurality of actually measured neutron signal sequences, wherein the spliced sequence is substantially equal to a signal obtained by continuously shifting 1 neutron detector. Compared with the signals obtained by 1 neutron detector, the sequence adopted by the invention has no signal fluctuation caused by physical gear shifting. The inventor names the method as a multi-detector relay signal splicing technology. At present, no other unit develops the research of a reactivity measurement method based on multi-detector relay signals or reports of related patent technologies in China. Therefore, it is necessary to master an autonomy reactivity measurement method based on multi-detector relay signals according to the difference between the autonomous reactor core and the foreign reactor core, and to provide capability for efficient and accurate measurement of reactor core reactivity.
Preferably, factor E is evaluatedi(t) is obtained by calculation according to the formula (1),
wherein E isi(t) represents the evaluation factor at time t of the ith neutron detector,
represents the neutron flux value at the time t of the ith neutron detector,
indicating the ith neutron probeNeutron flux value at time 0 of the detector.
Preferably, the specific process of acquiring the sliding window is as follows:
when n is 2, the amplified neutron measurement signal sequences of the 2 neutron detectors are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
The level of a small current amplification gear of the 1 st neutron detector is I, the level of a small current amplification gear of the 2 nd neutron detector is II,
after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 1 st neutron detector and the corresponding small current amplification gear level I, it can be determined that: vm,1(t0) To Vm,1(tk) As a valid segment sequence, Vm,1(tk+1) To Vm,1(tend) Is an invalid segment sequence;
after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 2 nd neutron detector and the corresponding small current amplification gear level II, it can be determined that: vm,2(t0) To Vm,2(tL-1) For invalid segment sequences, Vm,2(tL) To Vm,2(tend) Is an invalid segment sequence;
the neutron detector sets up suitable and unchangeable undercurrent amplification gear separately and satisfies: the signal amplitudes of the detectors passing through the respective low-current amplifiers can be mutually partially covered, so that Vm,1(tk) At a time later than Vm,2(tL) Of time, therefore, Vm,2(tL) To Vm,1(tk) As the crossing region of effective amplitude, corresponding to tLTo tkThe time period of (a) is denoted as a sliding window.
Preferably, the specific process of obtaining the simulated gear shifting splicing point by the optimal evaluation is as follows:
when n is 2, the sequences of 2 groups of evaluation factors are respectively:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(tend)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(tend)；
From the sliding window: t is tLTo tkGet thThe time point is used as an estimated gear shifting imitating splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated gear shifting imitating splicing point as a node, and the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence to obtain the estimated evaluation factor sequences as follows:
E1(t0)、E1(t1)、...E1(tL)、...、E1(th)、E2(th)、...、E2(tk)、...、E2(tend)；
wherein E is1(t0)、E1(t1)、...E1(tL)、...、E1(th) Evaluation factor sequence assigned to the 1 st neutron detector, E2(th)、...、E2(tk)、...、E2(tend) The evaluation factor sequence attributed to the 2 nd neutron detector,
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will thSequentially setting the next time point, and repeating the process until all time points in the sliding window are subjected to sliding calculation;
and comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point.
Preferably, the specific process of obtaining the simulated gear shifting splicing point by the optimal evaluation is as follows:
when n is 2, the sequences of 2 groups of evaluation factors are respectively:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(tend)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(tend)；
From the sliding window: t is tLTo tkGet thThe time point is used as an estimated gear shifting imitating splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated gear shifting imitating splicing point as a node, and the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence to obtain the estimated evaluation factor sequences as follows:
E1(tL)、...、E1(th)、E2(th)、...、E2(tk)；
wherein E is1(tL)、...、E1(th) Evaluation factor sequence assigned to the 1 st neutron detector, E2(th)、...、E2(tk) The evaluation factor sequence attributed to the 2 nd neutron detector,
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will thSequentially setting the next time point, and repeating the process until all time points in the sliding window are subjected to sliding calculation;
and comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point.
In the method, two specific processes of obtaining the simulated gear shifting splicing point through optimal evaluation are provided, namely two evaluated evaluation factor sequences are provided, data in a time period of a full measurement process adopted in the first method participate in optimal evaluation calculation, and data in a sliding window section participate in optimal evaluation calculation adopted in the second method. Both can be achieved theoretically, wherein the first one has the advantages that: the evaluation result is more fidelity, but the calculation amount is large, so that the variance is small, and obvious gear shifting points are difficult to find; the first of these has the advantage that: the variance is large, obvious gear shifting points are easy to find, the calculation amount is small, and the uncertainty of distortion exists in the evaluation result.
Preferably, the specific process for constructing the sequence required for reactivity calculation of S6 is as follows:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, directly taking the amplified neutron measurement signal sequences as 2 groups of reactivity signal sequences to be detected without processing;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
Start time: t is t0And the control ending time: t is tendThe sliding window is as follows: t is tLTo tkSetting the imitated gear shifting splicing point as tXWherein, tLTo tkAt t0To tendWithin the range of tXAt tLTo tkWithin the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
The sequence required for reactivity calculation was:
Vm,1(t0)，Vm,1(t1)，，…，Vm,1(tX)，Vm,2(tx)，…，Vm,2(tend)。
Vm,1(tX) Is a simulated gear shifting splicing point t in an amplified neutron measurement signal sequence of a 1 st neutron detectorXThe amplified neutron measurement signal of (1); vm,2(tx) Imitated gear shifting splicing point t in amplified neutron measurement signal sequence of No. 2 neutron detectorXThe amplified neutron measurement signal of (1).
When reactivity is calculated by V in the desired sequencem,1(tX) And Vm,2(tx) When the temperature of the water is higher than the set temperature,
Vm,1(tX) Is retained in the sequence required for reactivity calculation, Vm,2(tx) The device is used for signal synchronization processing;
preferably, the specific process for constructing the sequence required for reactivity calculation of S6 is as follows:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, optimizing the amplified neutron measurement signal sequences to obtain 2 groups of reactivity signal sequences to be detected;
the formula of the optimization process is as follows:
wherein, Vi(t) is the optimized processing result of the amplified neutron measurement signal of the ith neutron detector at the time t, phiDet,i(t) is the neutron flux value of the ith neutron detector at time t, Vm,i(t) is an amplified neutron measurement signal of the ith neutron detector at the time t;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
Start time: t is t0And the control ending time: t is tendThe sliding window is as follows: t is tLTo tkSetting the imitated gear shifting splicing point as tXWherein, tLTo tkAt t0To tendWithin the range of tXAt tLTo tkWithin the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
After optimization, 2 groups of reactive signal sequences to be detected are respectively as follows:
The sequence required for reactivity calculation was:
V1(t0)，V1(t1)，…，V1(tx)，V2(tx)，…，V2(tend)；
V1(tx) For the simulated gear shift splicing point t in the 1 st group of reactive signal sequences to be testedXTo be measured for a reactivity signal; v2(tx) Imitated gear shifting splicing point t in group 2 reactive signal sequence to be testedXTo be measured for a reactivity signal.
When reactivity is calculated by V in the desired sequence1(tx)，V2(tx) When the temperature of the water is higher than the set temperature,
V1(tx) Is retained in the sequence required for reactivity calculation, V2(tx) For signal synchronization processing.
In the method, the invention provides 2 to-be-detected reactive signal sequences, and the to-be-detected reactive signal sequences can be obtained by directly splicing original amplified neutron measurement signal sequences; or the amplified neutron measurement signal sequence can be optimized and then spliced to obtain a reactivity signal sequence to be detected.
In the above method, the setting of n to 2 only expresses the processing method performed in the case of 1 shift position, and is not a limitation of the present invention. When n is another number, if there are a plurality of shift positions, the processing may be performed according to the above-described equivalent principle.
The invention has the beneficial effects that: aiming at the influence of factors such as shift fluctuation of a neutron measurement signal in the reactivity measurement process, a reactor reactivity measurement method capable of reducing the influence is established. By adopting a method based on multi-detector relay measurement, signal shift fluctuation is reduced, errors introduced by a space effect of a detector are reduced, and the efficiency and the accuracy of the reactivity measurement are improved.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a schematic flow chart of the present invention.
FIG. 2 is a schematic diagram of an enlarged neutron measurement signal sequence.
FIG. 3 is a schematic diagram of the sequence of the evaluation factors.
FIG. 4 is a diagram showing the sequence of the evaluated evaluation factor after recombination.
FIG. 5 is another schematic diagram of the sequences of the evaluated evaluation factors after recombination.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
Example 1
As shown in fig. 1-5:
a reactivity measurement method based on multi-detector relay signals comprises the following steps:
s1, layout of neutron detectors: arranging n neutron detectors at a reactivity measuring station of a target reactor, and setting proper and fixed low-current amplification gears for the n neutron detectors respectively, so that signal amplitudes of the detectors passing through respective low-current amplifiers can be partially covered with each other and can be connected with each other, and a wider measuring range is realized;
s2, reactivity control: performing reactivity control on the target reactor from the control start time t0To the control end time tendContinuously monitoring by using n neutron detectors to obtain n groups of amplified neutron measurement signal sequences corresponding to the n neutron detectors one by one;
s3, acquiring a sliding window: observing the signal amplitude of each time point of each group of amplified neutron measurement signal sequences, finding an effective amplitude crossing region of one group of amplified neutron measurement signal sequences and the other group of amplified neutron measurement signal sequences, and recording a time segment corresponding to the effective amplitude crossing region as a sliding window;
s4, physical modeling: according to the physical parameters of the target reactor and the neutron detector in the step S1, three-dimensional physical modeling is carried out in a computer, and the three-dimensional physical modeling is calculated by using the built model to obtain: from the control start time t0To the control end time tendN groups of neutron flux value sequences corresponding to the n neutron detectors one by one;
s4, calculating an evaluation factor: calculating to obtain an evaluation factor E aiming at the neutron flux value corresponding to each moment in the n groups of neutron flux value sequencesi(t), obtaining n groups of evaluation factor sequences;
s5, obtaining a simulated gear shifting splicing point through optimal evaluation:
taking 1 time point from the sliding window as an estimated simulated gear shifting splicing point, truncating n groups of evaluation factor sequences by taking the estimated simulated gear shifting splicing point as a node, recombining and sequencing the truncated evaluation factor sequences of different segments according to time sequence to obtain an estimated evaluation factor sequence, and carrying out average value calculation on the estimated evaluation factor sequence to obtain an estimated value; repeating the process until all the time points in the sliding window are subjected to sliding calculation;
comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point;
s6, constructing a sequence required by reactivity calculation:
firstly, constructing n groups of reactivity signal sequences to be detected based on n groups of amplified neutron measurement signal sequences;
cutting off n groups of reactive signal sequences to be detected by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the cut-off n groups of reactive signal sequences to be detected in different sections from small to large according to the absolute amplitude of output signals to obtain a sequence required by reactivity calculation;
s7, reactivity calculation: reactivity calculations were performed based on the sequence required for reactivity calculations.
The invention aims to solve the problems that when the change range span of the neutron signal is large and gear shifting exists, the measurement accuracy is improved, the measurement efficiency can be improved while the wide-range measurement is realized, and the measurement accuracy can also be improved. The design principle of the invention is as follows:
adopt a plurality of neutron detectors, the gear that the undercurrent of every neutron detector of fixer launched, these gears are when the earlier stage is experimental, ensure that it can cover each other to set up, thereby realize the wide range, then, the signal synchronization with these neutron detectors is recorded simultaneously, adopt backstage computer to splice these signals, keep every effective data of surveying in its gear measuring range, like this to the system, shift gears the action in its measuring process, consequently, there is not signal fluctuation, thereby guarantee the measurement accuracy. All that is needed in the present invention is to find a suitable splice point. For each neutron detector, the measured signal is a continuous signal, and thus, it can be regarded as a signal sequence. The invention provides a processing method for searching an optimal splicing point through theoretical data, which has the following principle: performing physical modeling on physical parameters based on the target reactor and the neutron detector, for example, adopting Monte Carlo method physical calculation analysis software based on a high-fidelity model or other physical calculation analysis software; thus, theoretical neutron flux value data can be obtained, and the neutron flux value data and an actual neutron measurement signal are combined for evaluation; the best splice point is found. In order to simplify the calculation, the gear covering area between the neutron detectors is used as a sliding window for calculation, and the calculated evaluation values are analyzed and compared to find the optimal splicing point. And then splicing a plurality of actually measured neutron signal sequences, wherein the spliced sequence is substantially equal to a signal obtained by continuously shifting 1 neutron detector. Compared with the signals obtained by 1 neutron detector, the sequence adopted by the invention has no signal fluctuation caused by physical gear shifting. The inventor names the method as a multi-detector relay signal splicing technology. At present, no other unit develops the research of a reactivity measurement method based on multi-detector relay signals or reports of related patent technologies in China. Therefore, it is necessary to master an autonomy reactivity measurement method based on multi-detector relay signals according to the difference between the autonomous reactor core and the foreign reactor core, and to provide capability for efficient and accurate measurement of reactor core reactivity.
Preferably, factor E is evaluatedi(t) is obtained by calculation according to the formula (1),
wherein E isi(t) represents the evaluation factor at time t of the ith neutron detector,
represents the neutron flux value at the time t of the ith neutron detector,
represents the neutron flux value at time 0 of the ith neutron detector.
Preferably, the specific process of acquiring the sliding window is as follows:
when n is 2, the amplified neutron measurement signal sequences of the 2 neutron detectors are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
Specifically, as shown in FIG. 3, let t end100, wherein the 1 st neutron detector small current amplification step level is set as I, therefore, the step is at t0To t60For an effective measurement range, the 1 st neutron detector small current amplification gear level is set to II, so that it is at t40To t100Is the effective measurement range; their 2 neutron detectors are at t40To t60The effective measurement ranges of (1) are mutually covered;
therefore, after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 1 st neutron detector and the corresponding small current amplification gear level I, it can be determined that: vm,1(t0) To Vm,1(tk) As a valid segment sequence, Vm,1(tk+1) To Vm,1(tend) Is an invalid segment sequence;
after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 2 nd neutron detector and the corresponding small current amplification gear level II, it can be determined that: vm,2(t0) To Vm,2(tL-1) For invalid segment sequences, Vm,2(tL) To Vm,2(tend) Is an invalid segment sequence;
the neutron detector sets up suitable and unchangeable undercurrent amplification gear separately and satisfies: the signal amplitudes of the detectors passing through the respective low-current amplifiers can be mutually partially covered, so that Vm,1(tk) At a time later than Vm,2(tL) At the time of the day,thus, Vm,2(tL) To Vm,1(tk) As the crossing region of effective amplitude, corresponding to tLTo tkIs recorded as a sliding window, t40To t60Is a sliding window.
Example 2
Preferably, on the basis of the above-mentioned embodiments,
let t end100, wherein the 1 st neutron detector small current amplification step level is set as I, therefore, the step is at t0To t60For an effective measurement range, the 1 st neutron detector small current amplification gear level is set to II, so that it is at t40To t100Is the effective measurement range; their 2 neutron detectors are at t40To t60The effective measurement ranges of (1) are mutually covered; t is t40To t60Is a sliding window.
The specific process for obtaining the simulated gear shifting splicing point through optimal evaluation comprises the following steps:
when n is 2, as shown in FIG. 3, the sequences of 2 groups of evaluation factors are:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(tend)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(tend)；
As shown in FIG. 4, ERecombination (t45)To get t at50The time point is used as a recombined estimated evaluation factor sequence obtained by the estimated simulated gear shifting splicing point; eRecombination (t50)To get t at50The time point is used as a recombined estimated evaluation factor sequence obtained by the estimated simulated gear shifting splicing point;
specifically, from the sliding window: t is t40To t60Get t45The time point is used as an estimated simulated gear shifting splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated simulated gear shifting splicing point as a node, the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence,as shown in fig. 4, the evaluated evaluation factor sequence was obtained as:
E1(t0)、E1(t1)、...E1(t40)、...、E1(t45)、E2(t45)、...、E2(t60)、...、E2(t100)；
wherein E is1(t0)、E1(t1)、...E1(t40)、...、E1(t45) Evaluation factor sequence assigned to the 1 st neutron detector, E2(t45)、...、E2(t60)、...、E2(t100) The evaluation factor sequence attributed to the 2 nd neutron detector,
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will t40The next time point in sequence, as shown in FIG. 4, will be t45Change to t50And repeating the process to obtain the sequence of the evaluated evaluation factors as follows:
E1(t0)、E1(t1)、...E1(t40)、...、E1(t50)、E2(t50)、...、E2(t60)、...、E2(t100)；
wherein E is1(t0)、E1(t1)、...E1(t40)、...、E1(t50) Evaluation factor sequence assigned to the 1 st neutron detector, E2(t50)、...、E2(t60)、...、E2(t100) A sequence of evaluation factors attributed to the 2 nd neutron detector;
repeating the process until all the time points in the sliding window are subjected to sliding calculation;
comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point, and setting txImitating a gear shifting splicing point.
Preferably, on the basis of the above-mentioned embodiments,
let t end100, wherein the 1 st neutron detector small current amplification step level is set as I, therefore, the step is at t0To t60For an effective measurement range, the 1 st neutron detector small current amplification gear level is set to II, so that it is at t40To t100Is the effective measurement range; their 2 neutron detectors are at t40To t60The effective measurement ranges of (1) are mutually covered; t is t40To t60Is a sliding window.
The specific process for obtaining the simulated gear shifting splicing point through optimal evaluation comprises the following steps:
when n is 2, as shown in FIG. 3, the sequences of 2 groups of evaluation factors are:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(t100)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(t100)；
As shown in FIG. 5, ERecombination (t45)To get t at50The time point is used as a recombined estimated evaluation factor sequence obtained by the estimated simulated gear shifting splicing point; eRecombination (t50)To get t at50The time point is used as a recombined estimated evaluation factor sequence obtained by the estimated simulated gear shifting splicing point;
specifically, from the sliding window: t is t40To t60Get t45The time point is used as an estimated simulated gear shifting splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated simulated gear shifting splicing point as a node, the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence,
E1(t40)、...、E1(t45)、E2(t45)、...、E2(t60)；
wherein E is1(t40)、...、E1(t45) Ascribed to the 1 st neutron detectorEvaluation factor sequence of (E)2(t45)、...、E2(t60) A sequence of evaluation factors attributed to the 2 nd neutron detector;
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will thThe next time point is the sequence, t45Change to t50And repeating the process to obtain the sequence of the evaluated evaluation factors as follows:
E1(t40)、...、E1(t50)、E2(t50)、...、E2(t60)；
wherein E is1(t40)、...、E1(t50) Evaluation factor sequence assigned to the 1 st neutron detector, E2(t50)、...、E2(t60) A sequence of evaluation factors attributed to the 2 nd neutron detector;
until all time points in the sliding window are subjected to sliding calculation;
and comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point.
In the method, two specific processes of obtaining the simulated gear shifting splicing point through optimal evaluation are provided, namely two evaluated evaluation factor sequences are provided, data in a time period of a full measurement process adopted in the first method participate in optimal evaluation calculation, and data in a sliding window section participate in optimal evaluation calculation adopted in the second method. Both can be achieved theoretically, wherein the first one has the advantages that: the evaluation result is more fidelity, but the calculation amount is large, so that the variance is small, and obvious gear shifting points are difficult to find; the first of these has the advantage that: the variance is large, obvious gear shifting points are easy to find, the calculation amount is small, and the uncertainty of distortion exists in the evaluation result.
Example 3
As shown in fig. 1-5:
based on the above example, S6, the specific process for constructing the sequence required for reactivity calculation is:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, directly taking the amplified neutron measurement signal sequences as 2 groups of reactivity signal sequences to be detected without processing;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
Start time: t is t0And the control ending time: t is t100The sliding window is as follows: t is t40To t60Setting the imitated gear shifting splicing point as tXWherein, t40To t60At t0To t100Within the range of tXAt t40To t60Within the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(t100)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(t100)；
The sequence required for reactivity calculation was:
Vm,1(t0)，Vm,1(t1)，，…，Vm,1(tX)，Vm,2(tx)，…，Vm,2(t100)。
Vm,1(tX) Is a simulated gear shifting splicing point t in an amplified neutron measurement signal sequence of a 1 st neutron detectorXThe amplified neutron measurement signal of (1); vm,2(tx) Imitated gear shifting splicing point t in amplified neutron measurement signal sequence of No. 2 neutron detectorXThe amplified neutron measurement signal of (1).
When reactivity is calculated by V in the desired sequencem,1(tX) And Vm,2(tx) When the temperature of the water is higher than the set temperature,
Vm,1(tX) Is retained in the sequence required for reactivity calculation, Vm,2(tx) The device is used for signal synchronization processing;
preferably, the specific process for constructing the sequence required for reactivity calculation of S6 is as follows:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, optimizing the amplified neutron measurement signal sequences to obtain 2 groups of reactivity signal sequences to be detected;
the formula of the optimization process is as follows:
wherein, Vi(t) is the optimized processing result of the amplified neutron measurement signal of the ith neutron detector at the time t, phiDet,i(t) is the neutron flux value of the ith neutron detector at time t, Vm,i(t) is an amplified neutron measurement signal of the ith neutron detector at the time t;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
Start time: t is t0And the control ending time: t is t100The sliding window is as follows: t is t40To t60Setting the imitated gear shifting splicing point as tXWherein, t40To t60At t0To t100Within the range of tXAt t40To t60Within the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
in 1 stAmplified neutron measurement signal sequence of the sub-detectors: vm,1(t0)，Vm,1(t1)，…，Vm,1(t100)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(t100)；
After optimization, 2 groups of reactive signal sequences to be detected are respectively as follows:
The sequence required for reactivity calculation was:
V1(t0)，V1(t1)，…，V1(tx)，V2(tx)，…，V2(t100)；
V1(tx) For the simulated gear shift splicing point t in the 1 st group of reactive signal sequences to be testedXTo be measured for a reactivity signal; v2(tx) Imitated gear shifting splicing point t in group 2 reactive signal sequence to be testedXTo be measured for a reactivity signal.
When reactivity is calculated by V in the desired sequence1(tx)，V2(tx) When the temperature of the water is higher than the set temperature,
V1(tx) Is retained in the sequence required for reactivity calculation, V2(tx) For signal synchronization processing.
In the method, the invention provides 2 to-be-detected reactive signal sequences, and the to-be-detected reactive signal sequences can be obtained by directly splicing original amplified neutron measurement signal sequences; or the amplified neutron measurement signal sequence can be optimized and then spliced to obtain a reactivity signal sequence to be detected.
In the above embodiment, setting n to 2 only expresses the processing method performed in the case of 1 shift position, and is not a limitation of the present invention. When n is another number, if there are a plurality of shift positions, the processing may be performed according to the above-described equivalent principle.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A reactivity measurement method based on multi-detector relay signals is characterized by comprising the following steps:
s1, layout of neutron detectors: arranging n neutron detectors at a reactivity measuring station of a target reactor, and setting proper and fixed low-current amplification gears for the n neutron detectors respectively, so that signal amplitudes of the detectors passing through respective low-current amplifiers can be partially covered with each other and can be connected with each other, and a wider measuring range is realized;
s2, reactivity control: performing reactivity control on the target reactor from the control start time t0To the control end time tendContinuously monitoring by using n neutron detectors to obtain n groups of amplified neutron measurement signal sequences corresponding to the n neutron detectors one by one;
s3, acquiring a sliding window: observing the signal amplitude of each time point of each group of amplified neutron measurement signal sequences, finding an effective amplitude crossing region of one group of amplified neutron measurement signal sequences and the other group of amplified neutron measurement signal sequences, and recording a time segment corresponding to the effective amplitude crossing region as a sliding window;
s4, physical modeling: according to the physical parameters of the target reactor and the neutron detector in the step S1, three-dimensional physical modeling is carried out in a computer, and the three-dimensional physical modeling is calculated by using the built model to obtain: slave controlSystem start time t0To the control end time tendN groups of neutron flux value sequences corresponding to the n neutron detectors one by one;
s4, calculating an evaluation factor: calculating to obtain an evaluation factor E aiming at the neutron flux value corresponding to each moment in the n groups of neutron flux value sequencesi(t), obtaining n groups of evaluation factor sequences;
s5, obtaining a simulated gear shifting splicing point through optimal evaluation:
taking 1 time point from the sliding window as an estimated simulated gear shifting splicing point, truncating n groups of evaluation factor sequences by taking the estimated simulated gear shifting splicing point as a node, recombining and sequencing the truncated evaluation factor sequences of different segments according to time sequence to obtain an estimated evaluation factor sequence, and carrying out average value calculation on the estimated evaluation factor sequence to obtain an estimated value; repeating the process until all the time points in the sliding window are subjected to sliding calculation;
comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point;
s6, constructing a sequence required by reactivity calculation:
firstly, constructing n groups of reactivity signal sequences to be detected based on n groups of amplified neutron measurement signal sequences;
cutting off n groups of reactive signal sequences to be detected by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the cut-off n groups of reactive signal sequences to be detected in different sections from small to large according to the absolute amplitude of output signals to obtain a sequence required by reactivity calculation;
s7, reactivity calculation: reactivity calculations were performed based on the sequence required for reactivity calculations.
2. The method for reactivity measurement based on multi-probe relay signals according to claim 1, wherein the evaluation factor E isi(t) is obtained by calculation according to the formula (1),
wherein E isi(t) represents the evaluation factor at time t of the ith neutron detector,
represents the neutron flux value at the time t of the ith neutron detector,
represents the neutron flux value at time 0 of the ith neutron detector.
3. The method for measuring the responsiveness based on the multi-detector relay signal according to claim 1, wherein the specific process for acquiring the sliding window is as follows:
when n is 2, the amplified neutron measurement signal sequences of the 2 neutron detectors are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
The level of a small current amplification gear of the 1 st neutron detector is I, the level of a small current amplification gear of the 2 nd neutron detector is II,
after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 1 st neutron detector and the corresponding small current amplification gear level I, it can be determined that: vm,1(t0) To Vm,1(tk) As a valid segment sequence, Vm,1(tk+1) To Vm,1(tend) Is an invalid segment sequence;
after observing the relationship between the signal amplitude of the amplified neutron measurement signal sequence of the 2 nd neutron detector and the corresponding small current amplification gear level II, it can be determined that: vm,2(t0) To Vm,2(tL-1) For invalid segment sequences, Vm,2(tL) To Vm,2(tend) Is an invalid segment sequence;
the neutron detector sets up suitable and unchangeable undercurrent amplification gear separately and satisfies: the signal amplitudes of the detectors passing through the respective low-current amplifiers can be mutually partially covered, so that Vm,1(tk) At a time later than Vm,2(tL) Of time, therefore, Vm,2(tL) To Vm,1(tk) As the crossing region of effective amplitude, corresponding to tLTo tkThe time period of (a) is denoted as a sliding window.
4. The method for reactivity measurement based on multi-probe relay signals according to claim 1,
the specific process for obtaining the simulated gear shifting splicing point through optimal evaluation comprises the following steps:
when n is 2, the sequences of 2 groups of evaluation factors are respectively:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(tend)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(tend)；
From the sliding window: t is tLTo tkGet thThe time point is used as an estimated gear shifting imitating splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated gear shifting imitating splicing point as a node, and the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence to obtain the estimated evaluation factor sequences as follows:
E1(t0)、E1(t1)、...E1(tL)、...、E1(th)、E2(th)、...、E2(tk)、...、E2(tend)；
wherein E is1(t0)、E1(t1)、...E1(tL)、...、E1(th) Evaluation factor sequence assigned to the 1 st neutron detector, E2(th)、...、E2(tk)、...、E2(tend) The evaluation factor sequence attributed to the 2 nd neutron detector,
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will thSequentially setting the next time point, and repeating the process until all time points in the sliding window are subjected to sliding calculation;
and comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point.
5. The method for reactivity measurement based on multi-probe relay signals according to claim 1,
the specific process for obtaining the simulated gear shifting splicing point through optimal evaluation comprises the following steps:
when n is 2, the sequences of 2 groups of evaluation factors are respectively:
evaluation factor sequence of the 1 st neutron detector: e1(t0)、E1(t1)、...E1(tend)；
Evaluation factor sequence of the 2 nd neutron detector: e2(t0)、E2(t1)、...E2(tend)；
From the sliding window: t is tLTo tkGet thThe time point is used as an estimated gear shifting imitating splicing point, 2 groups of evaluation factor sequences are cut off by taking the estimated gear shifting imitating splicing point as a node, and the cut-off evaluation factor sequences of different sections are recombined and sequenced according to the time sequence to obtain the estimated evaluation factor sequences as follows:
E1(tL)、...、E1(th)、E2(th)、...、E2(tk)；
wherein E is1(tL)、...、E1(th) Evaluation factors attributed to the 1 st neutron detectorSubsequence, E2(th)、...、E2(tk) The evaluation factor sequence attributed to the 2 nd neutron detector,
carrying out average value calculation on the evaluated evaluation factor sequence to obtain an evaluation value; will thSequentially setting the next time point, and repeating the process until all time points in the sliding window are subjected to sliding calculation;
and comparing the evaluation value with '1', recording the evaluation value closest to '1' as an optimal evaluation value, and regarding a time point corresponding to the optimal evaluation value as a simulated gear shifting splicing point.
6. The method for reactivity measurement based on multi-probe relay signals according to claim 1,
s6, the specific process for constructing the sequence required by reactivity calculation is as follows:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, directly taking the amplified neutron measurement signal sequences as 2 groups of reactivity signal sequences to be detected without processing;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
7. The method for reactivity measurement based on multi-probe relay signals according to claim 6,
start time: t is t0And the control ending time: t is tendThe sliding window is as follows: t is tLTo tkSetting the imitated gear shifting splicing point as tXWherein, tLTo tkAt t0To tendWithin the range of tXAt tLTo tkWithin the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
The sequence required for reactivity calculation was:
Vm,1(t0)，Vm,1(t1)，，…，Vm,1(tX)，Vm,2(tx)，…，Vm,2(tend)。
Vm,1(tX) Is a simulated gear shifting splicing point t in an amplified neutron measurement signal sequence of a 1 st neutron detectorXThe amplified neutron measurement signal of (1); vm,2(tx) Imitated gear shifting splicing point t in amplified neutron measurement signal sequence of No. 2 neutron detectorXThe amplified neutron measurement signal of (1).
8. The method for reactivity measurement based on multi-probe relay signals according to claim 1,
s6, the specific process for constructing the sequence required by reactivity calculation is as follows:
when the number n is 2, the compound is,
firstly, constructing 2 groups of reactivity signal sequences to be detected based on 2 groups of amplified neutron measurement signal sequences; wherein, the process of constructing 2 groups of reactivity signal sequences to be detected is as follows: 2, optimizing the amplified neutron measurement signal sequences to obtain 2 groups of reactivity signal sequences to be detected;
the formula of the optimization process is as follows:
wherein, Vi(t) is an optimization processing junction of the amplified neutron measurement signal of the ith neutron detector at the time tFruit, phiDet,i(t) is the neutron flux value of the ith neutron detector at time t, Vm,i(t) is an amplified neutron measurement signal of the ith neutron detector at the time t;
and cutting off 2 groups of amplified neutron measurement signal sequences by taking the simulated gear shifting splicing points as nodes, and recombining and sequencing the 2 groups of amplified neutron measurement signal sequences of different sections after cutting off according to the sequence from small to large of the absolute amplitude of the output signals to obtain a sequence required by reactivity calculation.
9. The method for reactivity measurement based on multi-probe relay signals according to claim 8,
start time: t is t0And the control ending time: t is tendThe sliding window is as follows: t is tLTo tkSetting the imitated gear shifting splicing point as tXWherein, tLTo tkAt t0To tendWithin the range of tXAt tLTo tkWithin the range;
the 2 groups of amplified neutron measurement signal sequences are respectively as follows:
1, amplified neutron measurement signal sequence of a neutron detector: vm,1(t0)，Vm,1(t1)，…，Vm,1(tend)；
2, amplified neutron measurement signal sequence of the neutron detector: vm,2(t0)，Vm,2(t1)，…，Vm,2(tend)；
After optimization, 2 groups of reactive signal sequences to be detected are respectively as follows:
group 1 reactive signal sequences to be tested: v1(t0)，V1(t1)，…，V1(tend)；
Group 2 reactive signal sequences to be tested: v2(t0)，V2(t1)，…，V2(tend)；
The sequence required for reactivity calculation was:
V1(t0)，V1(t1)，…，V1(tx)，V2(tx)，…，V2(tend)；
V1(tx) For the simulated gear shift splicing point t in the 1 st group of reactive signal sequences to be testedXTo be measured for a reactivity signal; v2(tx) Imitated gear shifting splicing point t in group 2 reactive signal sequence to be testedXTo be measured for a reactivity signal.
10. The method for reactivity measurement based on multi-probe relay signals according to claim 7 or 9,
when reactivity is calculated by V in the desired sequencem,1(tX) And Vm,2(tx) When the temperature of the water is higher than the set temperature,
Vm,1(tX) Is retained in the sequence required for reactivity calculation, Vm,2(tx) The device is used for signal synchronization processing;
when reactivity is calculated by V in the desired sequence1(tx)，V2(tx) When the temperature of the water is higher than the set temperature,
V1(tx) Is retained in the sequence required for reactivity calculation, V2(tx) For signal synchronization processing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011551550.4A CN112687411B (en) | 2020-12-24 | 2020-12-24 | Reactivity measurement method based on multi-detector relay signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011551550.4A CN112687411B (en) | 2020-12-24 | 2020-12-24 | Reactivity measurement method based on multi-detector relay signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112687411A true CN112687411A (en) | 2021-04-20 |
| CN112687411B CN112687411B (en) | 2022-06-24 |
Family
ID=75452727
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011551550.4A Active CN112687411B (en) | 2020-12-24 | 2020-12-24 | Reactivity measurement method based on multi-detector relay signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112687411B (en) |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6181759B1 (en) * | 1999-07-23 | 2001-01-30 | Westinghouse Electric Company Llc | Method and apparatus for determining nearness to criticality of a nuclear fueled electric power generating unit |
| JP2001255391A (en) * | 2000-01-06 | 2001-09-21 | Toshiba Corp | Reactivity measuring method for reactor, measuring method for moderator temperature coefficient and measuring device of moderator temperature coefficient |
| US6430247B1 (en) * | 1999-07-05 | 2002-08-06 | Framatome | Method and system for monitoring at least one operating parameter of the core of a nuclear reactor |
| CN101105986A (en) * | 2007-08-03 | 2008-01-16 | 蔡光明 | Reactor reactivity measuring method |
| KR101216608B1 (en) * | 2012-02-13 | 2012-12-31 | 유저스(주) | Method and system for calculating online digital reactivity using field programmable gate array |
| ES2415031A1 (en) * | 2012-10-22 | 2013-07-23 | Centro De Investigaciones Energéticas, Medioambientales Y Tecnológicas (Ciemat) | System for monitoring reactivity in a subcritical nuclear reactor |
| CN103345950A (en) * | 2013-07-12 | 2013-10-09 | 中科华核电技术研究院有限公司 | Ex-core nuclear detection system and method for pressurized water reactor |
| CN103794256A (en) * | 2013-11-19 | 2014-05-14 | 国核（北京）科学技术研究院有限公司 | Method and system for monitoring abnormality of reactor |
| CN104215809A (en) * | 2014-09-24 | 2014-12-17 | 中国核动力研究设计院 | Method for switching power measuring span of digital reactivity instrument |
| CN104376884A (en) * | 2013-08-15 | 2015-02-25 | 中国原子能科学研究院 | Portable reactivity meter for carrying out full measuring range automatic monitoring on reactor signal, and reactivity correction method |
| CN104409120A (en) * | 2014-12-11 | 2015-03-11 | 中国核动力研究设计院 | Switching method for current linear amplification automatic range |
| JP2015094746A (en) * | 2013-11-14 | 2015-05-18 | 株式会社東芝 | Effective delayed neutron fraction measurement method and effective delayed neutron fraction measurement apparatus |
| CN107170498A (en) * | 2017-05-25 | 2017-09-15 | 中国核动力研究设计院 | Reactivity meter and control method with online dynamic scale control rod worth function |
| CN109117591A (en) * | 2018-09-13 | 2019-01-01 | 中国核动力研究设计院 | A kind of kinetic reactivity measurement method based on multi-detector measuring signal |
| CN109273119A (en) * | 2018-09-13 | 2019-01-25 | 中国核动力研究设计院 | Optimize the method for neutron detector position when measuring big reactivity on critical assembly |
| CN109961861A (en) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | Reactivity measurement instrument under a kind of nuclear power plant reactor depth subcritical state |
| CN111799005A (en) * | 2020-08-13 | 2020-10-20 | 福州奇正谷科技有限公司 | Different-range normalization method, storage medium and real-time online reactivity meter |
-
2020
- 2020-12-24 CN CN202011551550.4A patent/CN112687411B/en active Active
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6430247B1 (en) * | 1999-07-05 | 2002-08-06 | Framatome | Method and system for monitoring at least one operating parameter of the core of a nuclear reactor |
| US6181759B1 (en) * | 1999-07-23 | 2001-01-30 | Westinghouse Electric Company Llc | Method and apparatus for determining nearness to criticality of a nuclear fueled electric power generating unit |
| JP2001255391A (en) * | 2000-01-06 | 2001-09-21 | Toshiba Corp | Reactivity measuring method for reactor, measuring method for moderator temperature coefficient and measuring device of moderator temperature coefficient |
| CN101105986A (en) * | 2007-08-03 | 2008-01-16 | 蔡光明 | Reactor reactivity measuring method |
| KR101216608B1 (en) * | 2012-02-13 | 2012-12-31 | 유저스(주) | Method and system for calculating online digital reactivity using field programmable gate array |
| ES2415031A1 (en) * | 2012-10-22 | 2013-07-23 | Centro De Investigaciones Energéticas, Medioambientales Y Tecnológicas (Ciemat) | System for monitoring reactivity in a subcritical nuclear reactor |
| CN103345950A (en) * | 2013-07-12 | 2013-10-09 | 中科华核电技术研究院有限公司 | Ex-core nuclear detection system and method for pressurized water reactor |
| CN104376884A (en) * | 2013-08-15 | 2015-02-25 | 中国原子能科学研究院 | Portable reactivity meter for carrying out full measuring range automatic monitoring on reactor signal, and reactivity correction method |
| JP2015094746A (en) * | 2013-11-14 | 2015-05-18 | 株式会社東芝 | Effective delayed neutron fraction measurement method and effective delayed neutron fraction measurement apparatus |
| CN103794256A (en) * | 2013-11-19 | 2014-05-14 | 国核（北京）科学技术研究院有限公司 | Method and system for monitoring abnormality of reactor |
| CN104215809A (en) * | 2014-09-24 | 2014-12-17 | 中国核动力研究设计院 | Method for switching power measuring span of digital reactivity instrument |
| CN104409120A (en) * | 2014-12-11 | 2015-03-11 | 中国核动力研究设计院 | Switching method for current linear amplification automatic range |
| CN107170498A (en) * | 2017-05-25 | 2017-09-15 | 中国核动力研究设计院 | Reactivity meter and control method with online dynamic scale control rod worth function |
| CN109961861A (en) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | Reactivity measurement instrument under a kind of nuclear power plant reactor depth subcritical state |
| CN109117591A (en) * | 2018-09-13 | 2019-01-01 | 中国核动力研究设计院 | A kind of kinetic reactivity measurement method based on multi-detector measuring signal |
| CN109273119A (en) * | 2018-09-13 | 2019-01-25 | 中国核动力研究设计院 | Optimize the method for neutron detector position when measuring big reactivity on critical assembly |
| CN111799005A (en) * | 2020-08-13 | 2020-10-20 | 福州奇正谷科技有限公司 | Different-range normalization method, storage medium and real-time online reactivity meter |
Non-Patent Citations (5)
| Title |
|---|
| WANG WC DNEG: "The first application of modified neutron source multiplication method in subcriticality monitoring based on Monte Carlo", 《NUCLEAR ENGINEERING AND TECHNOLOGY》 * |
| 吴军 等: "数字化多通道宽量程ITER中子通量密度测量系统研制", 《原子能科学技术》 * |
| 王丽华等: "先进反应性测量系统的研发与应用", 《核动力工程》 * |
| 韩钰等: "核反应堆反应性测量技术研究", 《科技视界》 * |
| 黄礼渊 等: "反应堆物理试验用便携式数字反应性仪的研制", 《核电子学与探测技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112687411B (en) | 2022-06-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108333617B (en) | The quick Peak Search Method of radioactive material quality detection in a kind of seawater | |
| CN110929769B (en) | Reactor mechanical fault joint detection model, method and device based on vibration and sound | |
| CN110730435B (en) | Data drift blind calibration method for distributed wireless sensor network | |
| CN112687411B (en) | Reactivity measurement method based on multi-detector relay signals | |
| CN114492629A (en) | Abnormality detection method, abnormality detection device, electronic apparatus, and storage medium | |
| CN109839627B (en) | Sonar detection distance measurement and calibration system and method under pool condition | |
| CN114265001B (en) | Smart electric meter metering error evaluation method | |
| CN105407496B (en) | A kind of method of erroneous measurements in identification wireless sensor network | |
| CN112285650B (en) | Method, system and storage medium for positioning unknown wave velocity sound emission source in presence of abnormal TDOA | |
| CN115826070B (en) | Method for determining mining value position in white-sentry-type uranium ore to be detected by using core logger | |
| CN116522085A (en) | Full-automatic inhaul cable frequency extraction, fixed-order and cable force identification method and application | |
| Neveux et al. | Robust filtering for linear time-invariant continuous systems | |
| CN110702344A (en) | Closed-loop bridge structure damage diagnosis method based on distributed strain measurement technology | |
| Khosraviani et al. | Damage detection in continuous deck bridges using statistical cross-correlation function method | |
| CN116298670A (en) | Intelligent fault positioning method and system suitable for multi-branch distribution line | |
| CN111695501B (en) | Equipment soft fault detection method based on operating system kernel calling data | |
| CN117553938B (en) | Error correction method for temperature secondary instrument in high-low temperature environment | |
| CN103472435A (en) | Method for two-stage location of sound emission source of ocean platform complex structure system | |
| CN108152527B (en) | Digital speed measurement method based on median average filtering | |
| CN112967825A (en) | Reactivity measurement method based on correction signal uncertainty analysis | |
| CN113009107A (en) | Method for detecting unfrozen water content of frozen soil | |
| CN114737949B (en) | Sensor for measuring gas-liquid two-phase flow gas phase distribution and terminal system thereof | |
| CN116593830B (en) | Cable fault detection system and detection method | |
| CN117572853B (en) | Magnetic field controller performance test analysis management system | |
| CN113192659B (en) | Nuclear reactor reactivity measurement method based on primary neutron signal average sampling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |