CN101162455A - Method for simulating nuclear reactors in threshold state - Google Patents
Method for simulating nuclear reactors in threshold state Download PDFInfo
- Publication number
- CN101162455A CN101162455A CNA2006101408808A CN200610140880A CN101162455A CN 101162455 A CN101162455 A CN 101162455A CN A2006101408808 A CNA2006101408808 A CN A2006101408808A CN 200610140880 A CN200610140880 A CN 200610140880A CN 101162455 A CN101162455 A CN 101162455A
- Authority
- CN
- China
- Prior art keywords
- reactor
- neutron
- variable quantity
- boron concentration
- rod position
- 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
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
The present invention belongs to the nuclear reactor field, and in particular relates to a critical state simulation method of nuclear reactor, comprising the following steps: (1) a model is bononized and diluted to obtain the variable quantity of boron concentration; (2) the variable quantity of overlapping step rod position is obtained through the conversion formula of the overlapping step rod position and a single rod of a rod position control model; (3) the variable quantity of boron concentration and the variable quantity of overlapping step rod position are substituted into a reactive model to obtain reactive variable quantity; (4) neutron number is obtained through the reactive model; (5) when the boron concentration and the rod position of a control rod are changed and the neutron number is changeless, a reactor achieves a critical state then. The use of the method of the present invention can dynamically and intuitively show the relation of the neutron number and reactivity in real time. As a numerical solution based on computers is adopted, the neutron number can be very accurately worked out to the change of reactivity in any form. The present invention is few in resources needed, which can be realized with only one computer needed. Parameter can be flexibly changed so as to adapt the present invention to critical state achieving characteristics of the reactors in different structures.
Description
Technical field
The present invention relates to a kind of nuclear mockup method, particularly a kind of method for simulating nuclear reactors in threshold state.
Background technology
It is critical that nuclear reactor reaches, and the neutron number that is about to reactor is stabilized in some values and does not change, and can make nuclear process controlled, can utilize this stable energy that nuclear reaction discharged to generate electricity like this.
Make nuclear reactor reach critical conditions be any nuclear reactor operating personnel the essential technical ability of grasping, existing nuclear reactor reaches critical conditions to be realized by large-scale analog machine simulated operation.Use large-scale analog machine simulation nuclear reactor to reach critical conditions and have many shortcomings: manpower and materials are required height, need the special messenger to cooperate; Need daily servicing; Lack dirigibility, can not change parameter at any time and be used for theoretical research; Large-scale analog machine fixed-site, limited amount and simulation cycle are long.
Summary of the invention
The object of the present invention is to provide a kind of method for simulating nuclear reactors in threshold state.
A kind of method for simulating nuclear reactors in threshold state may further comprise the steps:
(1) reactor is carried out boronation dilution operation, obtain the boron concentration c that boronation or dilution are introduced
1, boronation or dilution volume v, reactor water capacity v
0, with parameter c
1, v, v
0The substitution boronation is diluted equation, obtains the boron concentration c variable quantity dc of reactor;
(2), calculate the variable quantity dl of rod position l of folded step by the folded step rod position of rod position controlling models and the conversion relation formula of single rod;
(3) with the reactive computation model of rod position variable quantity dl substitution of folded step that obtains in the boron concentration change amount dc that obtains in the step (1) and the step (2), draw reactive variable quantity d ρ;
(4) set neutron number n, delayed neutron fraction β, disintegration constant λ, on average for time l, add neutron source strength q, reactive ρ, reactive ρ is the reactive ρ by reactor model initial setting
0The changes of reactivity amount d ρ addition that obtains with reactive computation model constitutes, that is: ρ=ρ
0+ d ρ; Bring the parameter and the reactive ρ of above-mentioned all settings into variable quantity that the reactor model obtains neutron number in the unit interval
(5) boron concentration change amount dc changes thereupon during the variation of the boron concentration in the step (1), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes; Rod position variable quantity dl changed thereupon when control rod rod position changed in the step (2), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes; Neutron number n does not change when changing boron concentration and control rod rod position, i.e. the variable quantity of neutron number in the unit interval
The time, nuclear reactors in threshold state.
Use method of the present invention: can dynamically intuitively show neutron number and reactive relation in real time; Because adopt the computer based numerical solution, for any type of variation of reactivity, the neutron number average can be calculated by point-device; Resource requirement is few, only needs a computing machine to realize; Can change parameter neatly, with the critical characteristic that reaches of the reactor that adapts to different structure.
Embodiment
Embodiment 1
(1) reactor is carried out boronation dilution operation, obtain the boron concentration c that boronation or dilution are introduced
1, boronation or dilution volume v, reactor water capacity v
0, with parameter c
1, v, v
0Substitution boronation dilution equation, boronation dilution equation is as follows:
c
1dv=v
0dc+cdv
In the formula:
The boron concentration of c---reactor
c
1---by the boron concentration of boronation or dilution introducing, when being the boronation state, c
1Value be 7 * 10
-3, be 0 when for diluted state
v
0---the water capacity of reactor, this value is relevant with reactor structure, to same reactor, should value be similar to constant when reactor operation
V---boronation or dilution volume, this value depend on the poor of the existing boron concentration of reactor and itself and the boron concentration that will reach, and do not have concrete scope, and span can be very big, if do not consider the timeliness economy, theoretical value can be zero to infinity
Obtain the variable quantity dc of the boron concentration of the boron concentration c of reactor and reactor by following formula, boronation dilution operation stops when reactor reaches critical conditions.
(2) carry out the lifting or the decline of control rod position, when carrying out rod withdrawal control, increased the value of excellent figure place l of folded step; During plunger control, reduced the value of excellent figure place l of folded step; Reactor has 4 groups of control rods, and its highest rod position was 225 steps, and minimum rod position was 5 steps, and folded step rod position is as follows with the conversion relation formula of single rod:
In the formula:
The folded excellent figure place of step of l-
N
1-1 rod rod figure place
N
2-2 rods rod figure place
N
3-3 rods rod figure place
N
4-4 rods rod figure place
The overlapping step number of x-
Thereby calculate 1,2,3,4 rods, the folded numerical value that goes on foot rod position l by following formula, draw the variable quantity dl of a folded step rod l by the numerical value change of excellent position l of folded step.
(3) with the reactive computation model of rod position variable quantity dl substitution that obtains in the boron concentration change amount dc that obtains in the step (1) and the step (2); Owing to reach in the critical process, the coolant temperature variable quantity is small in the reactor, therefore when considering, ignored because of coolant temperature changes reactive influence to reactive influence factor, when having considered excellent position with the boron concentration change to reactive influence; Reactive computation model as shown in the formula:
dρ=k
1dc+k
2dl
In the formula:
D ρ-reactive variable quantity
k
1-boron concentration change amount influences the scale-up factor of changes of reactivity amount
k
2-rod position variable quantity influences the scale-up factor of changes of reactivity amount
The boron concentration change amount of dc-reactor
Dl-rod position variable quantity
Draw reactive variable quantity d ρ by following formula.
With rod position and boron concentration to reactive linearization, i.e. k in the following formula of influencing
1, k
2Coefficient is a constant, k
1Scope be-70 * 10
-5~0, k
2Scope be 0~10 * 10
-5As the more accurate numerical value of needs, and for different type of reactors or same reactor different times k
1, k
2The value difference, k
1, k
2Value depend on the burnup of reactor boron concentration, reactor coolant temperature, reactor structure, nuclear fuel, the position of control rod.Use conic fitting method or cubic curve fitting process to coefficient k
1, k
2Carrying out non-linearization handles.
(4) set neutron number n, delayed neutron fraction β, disintegration constant λ, on average for time l, add neutron source strength q, reactive ρ, reactive ρ is the reactive ρ by reactor model initial setting
0The changes of reactivity amount d ρ addition that obtains with reactive computation model constitutes, that is: ρ=ρ
0+ d ρ, wherein ρ
0In simulation process, be constant, for the simulation process ρ of homogeneous not as initial setting
0Be different constants; Bring the parameter and the reactive ρ of above-mentioned all settings into the reactor model, the reactor model is set up according to point-reactor kinetic equation, the reactor model as shown in the formula:
In the formula:
N---neutron number
I---because the nuclear reaction meeting produces multiple pioneer and examines type, each pioneer examines the fall time difference of type, has pioneer's nuclear of 6 kinds of main types, thus i=1,2 ... 6
β
i---delayed neutron fraction β is the neutron number that is produced by pioneer's nuclear decay, owing to have pioneer's nuclear of 6 kinds of main types, so β is arranged
i
β---delayed neutron fraction, delayed neutron fraction are to account for the ratio that nuclear reaction produces all neutron numbers by the neutron number that pioneer's nuclear decay produces, and are generally 0.0065,
λ---disintegration constant, its scope are 0.01~10 per second
L---on average for the time, its scope is 0.0849 second
Q---add the neutron source strength, depend on type of reactor, no determined value; The reactor that has does not have the neutron of adding source strength, promptly gets q=0
ρ---reactivity
c
i---pioneer's check figure, it depends on the power of reactor, theoretical value can be zero to infinitely great,, because of real reaction heap power limited, pioneer's check figure c
iCan not be infinitely great
Draw the variable quantity of neutron number in the unit interval by following formula
Rate of change with pioneer's check figure
(5) boron concentration change amount dc changes thereupon during the variation of the boron concentration in the step (1), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes; Rod position variable quantity dl changed thereupon when control rod rod position changed in the step (2), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes.
Suppose that neutron number is n when moment t, when moment t+dt, the neutron number that is gone out by the reactor Model Calculation is n ' so, nuclear reactors in threshold state during n '=n; Promptly when changing boron concentration and control rod rod position, the variable quantity of neutron number in the unit interval
Nuclear reactors in threshold state then.
When reactor reaches critical conditions, because neutron number n is constant, promptly
And delayed neutron fraction β is a constant, so the delayed neutron number
Constant, because pioneer's nuclear decay constant λ is constant, pretend pioneer's check figure c again into delayed neutron
iConstant, therefore
Variable quantity with neutron number in the unit interval
Bring the formula of neutron period T into:
In the formula:
T---the pile neutron cycle
N---neutron number
The variable quantity of neutron number in unit interval
Nuclear reactors in threshold state when being zero, as from the foregoing during nuclear reactors in threshold state in subcycle be infinity.
Embodiment 2
Be that with the different of embodiment 1 reactor in the step (2) has 3 groups of control rods, its highest rod position was 225 steps, and minimum rod position was 5 steps, and folded step rod position is as follows with the conversion relation formula of single rod:
In the formula:
The folded excellent figure place of step of l-
N
1-1 rod rod figure place
N
2-2 rods rod figure place
N
3-3 rods rod figure place
The overlapping step number of x-
Thereby calculate 1,2,3 rods, the folded numerical value that goes on foot rod position l by following formula, draw the variable quantity dl of a folded step rod l by the numerical value change of excellent position l of folded step.
Embodiment 3
Be that with the different of embodiment 1 reactor in the step (2) has A, B, four groups of control rods of C, D, its highest rod position was 220 steps, and minimum rod position was 0 step, and folded step rod position is as follows with the conversion relation formula of single rod:
In the formula:
The folded excellent figure place of step of l-
N
1-1 rod rod figure place
N
2-2 rods rod figure place
N
3-3 rods rod figure place
N
4-4 rods rod figure place
The overlapping step number of x-
Thereby calculate 1,2,3,4 rods, the folded numerical value that goes on foot rod position l by following formula, draw the variable quantity dl of a folded step rod l by the numerical value change of excellent position l of folded step.
Claims (7)
1. method for simulating nuclear reactors in threshold state may further comprise the steps:
(1) reactor is carried out boronation dilution operation, obtain the boron concentration c that boronation or dilution are introduced
1, boronation or dilution volume v, reactor water capacity v
0, with parameter c
1, v, v
0The substitution boronation is diluted equation, obtains the boron concentration c variable quantity dc of reactor;
(2), calculate the variable quantity dl of rod position l of folded step by the folded step rod position of rod position controlling models and the conversion relation formula of single rod;
(3) with the reactive computation model of rod position variable quantity dl substitution of folded step that obtains in the boron concentration change amount dc that obtains in the step (1) and the step (2), obtain reactive variable quantity d ρ;
(4) set neutron number n, delayed neutron fraction β, disintegration constant λ, on average for time l, add neutron source strength q, reactive ρ, reactive ρ is the reactive ρ by reactor model initial setting
0The changes of reactivity amount d ρ addition that obtains with reactive computation model constitutes, that is: ρ=ρ
0+ d ρ; Bring the parameter and the reactive ρ of above-mentioned all settings into variable quantity that the reactor model obtains neutron number in the unit interval
(5) boron concentration change amount dc changes thereupon during the variation of the boron concentration in the step (1), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes; Rod position variable quantity dl changed thereupon when control rod rod position changed in the step (2), and the variable quantity d ρ of inducing reaction property changes, thereby has changed the value of reactive ρ in the reactor model, so neutron number n changes; Neutron number n does not change when changing boron concentration and control rod rod position, i.e. the variable quantity of neutron number in the unit interval
The time, nuclear reactors in threshold state.
2. method for simulating nuclear reactors in threshold state as claimed in claim 1 is characterized in that described boronation dilution equation is:
c
1dv=v
0dc+cdv
In the formula:
The boron concentration of c---reactor
c
1---by the boron concentration of boronation or dilution introducing, when being the boronation state, c
1Value be 7 * 10
-3, be 0 when for diluted state
v
0---the water capacity of reactor, this value is relevant with reactor structure, to same reactor, should value be similar to constant when reactor operation
V---boronation or dilution volume, this value depend on the poor of the existing boron concentration of reactor and itself and the boron concentration that will reach, and do not have concrete scope, and span can be very big, if do not consider the timeliness economy, theoretical value can be zero to infinity
The boron concentration of c---reactor.
3. method for simulating nuclear reactors in threshold state as claimed in claim 1, it is characterized in that: when reactor has i group control rod, the highest rod position is a step, when minimum rod position go on foot for b, the excellent conversion relation formula of the folded step rod of described rod position controlling models and list as shown in the formula:
In the formula:
L-folded excellent figure place of step
N
iThe excellent figure place of-i
The overlapping step number of x-
A-the highest rod position step number
The minimum rod of b-position step number.
4. method for simulating nuclear reactors in threshold state as claimed in claim 1, it is characterized in that described reactive computation model as shown in the formula:
dρ=k
1dc+k
2dl
In the formula:
D ρ-reactive variable quantity
k
1-boron concentration change amount influences the scale-up factor of changes of reactivity amount
k
2-rod position variable quantity influences the scale-up factor of changes of reactivity amount
The boron concentration change amount of dc-reactor
Dl-rod position of folded step variable quantity.
5. method for simulating nuclear reactors in threshold state as claimed in claim 4 is characterized in that described boron concentration change amount influences the scale-up factor k of changes of reactivity amount
1Can be constant or nonlinear factor; k
1Its scope is-70 * 10 during for constant
-5~0; k
1Value depend on the burnup of reactor boron concentration, reactor coolant temperature, reactor structure, nuclear fuel, the position of control rod; k
1During for nonlinear factor, use conic fitting method or cubic curve fitting process to coefficient k
1Carrying out non-linearization handles.
6. critical state simulation method of nuclear reactor as claimed in claim 4 is characterized in that described rod position variable quantity influences the scale-up factor k of changes of reactivity amount
2All can be constant or nonlinear factor; k
2Its scope is 0~10 * 10 during for constant
-5, k
2Value depend on the burnup of reactor boron concentration, reactor coolant temperature, reactor structure, nuclear fuel, the position of control rod; k
2During for nonlinear factor, use conic fitting method or cubic curve fitting process to coefficient k
2Carrying out non-linearization handles.
7. method for simulating nuclear reactors in threshold state as claimed in claim 1 is characterized in that described reactor model sets up according to point-reactor kinetic equation, the reactor model as shown in the formula:
N---neutron number
I---because the nuclear reaction meeting produces multiple pioneer and examines type, each pioneer examines the fall time difference of type, has pioneer's nuclear of 6 kinds of main types, thus i=1,2 ... 6
β
i---delayed neutron fraction β is the neutron number that is produced by pioneer's nuclear decay, owing to have pioneer's nuclear of 6 kinds of main types, so β is arranged
i
β---delayed neutron fraction, delayed neutron fraction are to account for the ratio that nuclear reaction produces all neutron numbers by the neutron number that pioneer's nuclear decay produces, and are generally 0.0065,
λ---disintegration constant, its scope are 0.01~10 per second
L---on average for the time, its scope is 0.0849 second
Q---add the neutron source strength, depend on type of reactor, no determined value, the reactor that has does not have the neutron of adding source strength, promptly gets q=0
ρ---reactivity
c
i---pioneer's check figure, it depends on reactor capability, theoretical value can be zero to infinitely great, because of real reaction heap power limited, pioneer's check figure c
iCan not be infinitely great
---the variable quantity of neutron number in the unit interval
---the neutron number that all neutrons produce when being fast neutron
N---neutron number
Suppose that neutron number is n when moment t, when moment t+dt, obtaining neutron number by the reactor Model Calculation is n ' so, nuclear reactors in threshold state during n '=n; Promptly when changing other parameter boron concentration and control rod rod position, the variable quantity of neutron number in the unit interval
Nuclear reactors in threshold state then; When reactor reaches critical conditions, because neutron number n is constant, promptly
And delayed neutron fraction β is a constant, so the delayed neutron number
Constant, because pioneer's nuclear decay constant λ is constant, pretend pioneer's check figure c again into delayed neutron
iConstant, therefore
Variable quantity with neutron number in the unit interval
Bring the formula of neutron period T into:
In the formula:
T---the pile neutron cycle
N---neutron number
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2006101408808A CN101162455B (en) | 2006-10-13 | 2006-10-13 | Method for simulating nuclear reactors in threshold state |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2006101408808A CN101162455B (en) | 2006-10-13 | 2006-10-13 | Method for simulating nuclear reactors in threshold state |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101162455A true CN101162455A (en) | 2008-04-16 |
CN101162455B CN101162455B (en) | 2011-05-04 |
Family
ID=39297379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006101408808A Active CN101162455B (en) | 2006-10-13 | 2006-10-13 | Method for simulating nuclear reactors in threshold state |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101162455B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105788666A (en) * | 2014-12-23 | 2016-07-20 | 福建福清核电有限公司 | Nuclear reactor critical process control method |
CN112908501A (en) * | 2021-01-12 | 2021-06-04 | 中国原子能科学研究院 | Nuclear critical device measuring device and method |
CN113241204A (en) * | 2021-05-19 | 2021-08-10 | 核星核电科技(海盐)有限公司 | Special system for testing reactor reactivity instrument |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2639141B1 (en) * | 1988-11-14 | 1991-02-01 | Framatome Sa | PILOTAGE METHOD OF PRESSURIZED WATER NUCLEAR REACTOR AND NUCLEAR REACTOR USING THE SAME |
US6748348B1 (en) * | 1999-12-30 | 2004-06-08 | General Electric Company | Design method for nuclear reactor fuel management |
-
2006
- 2006-10-13 CN CN2006101408808A patent/CN101162455B/en active Active
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105788666A (en) * | 2014-12-23 | 2016-07-20 | 福建福清核电有限公司 | Nuclear reactor critical process control method |
CN112908501A (en) * | 2021-01-12 | 2021-06-04 | 中国原子能科学研究院 | Nuclear critical device measuring device and method |
CN112908501B (en) * | 2021-01-12 | 2023-03-21 | 中国原子能科学研究院 | Nuclear critical device measuring device and method |
CN113241204A (en) * | 2021-05-19 | 2021-08-10 | 核星核电科技(海盐)有限公司 | Special system for testing reactor reactivity instrument |
Also Published As
Publication number | Publication date |
---|---|
CN101162455B (en) | 2011-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Multi-stage progressive optimality algorithm and its application in energy storage operation chart optimization of cascade reservoirs | |
CN103094921B (en) | Direct-drive aerogenerator group electromechanical transient modeling method and machine-electricity transient model | |
Xingang et al. | China’s wind, biomass and solar power generation: What the situation tells us? | |
CN105449713A (en) | Active power distribution network soft normally open point (SNOP) planning method considering distributed generator characteristics | |
CN104578157A (en) | Load flow calculation method of distributed power supply connection power grid | |
CN103345546B (en) | The governor parameter discrimination method that frequency locus combines with particle cluster algorithm | |
CN103400208A (en) | Power distribution network distributive power supply optimal access capacity determining method based on cone optimization | |
CN103593711B (en) | A kind of distributed power source Optimal Configuration Method | |
CN102930343A (en) | Energy optimization method for distributed power generation and energy supply system | |
CN101162455B (en) | Method for simulating nuclear reactors in threshold state | |
CN103441527A (en) | Wind electricity connection system model based on measured data | |
CN105656084A (en) | Improved stochastic load flow algorithm involved with new energy power generation prediction errors | |
Yujie et al. | Medium-term forecasting of cold, electric and gas load in multi-energy system based on VAR model | |
CN102968747A (en) | Method for determining typical sunrise force curves of wind power station | |
CN110649596A (en) | Frequency full-response analytic model considering system initial state | |
CN107039981A (en) | One kind intends direct current linearisation probability optimal load flow computational methods | |
Villarino et al. | INVAP neutronic calculation line | |
CN105633950B (en) | A kind of probabilistic multiple target Random-fuzzy Dynamic Optimal Power Flow Problem method for solving of consideration wind-powered electricity generation injection | |
CN103729556A (en) | Orthogonal successive approximation algorithm based reservoir scheduling method | |
CN113721476A (en) | 100 MW-level variable-speed seawater pumped storage unit and renewable energy source combined operation system hardware-in-loop simulation platform and method | |
CN103559387B (en) | Fractal dimension steepest descent calculation method for valve point effect power grid economic dispatch problem | |
CN105896547A (en) | Large power network graded voltage control method under wind power access | |
CN110309565A (en) | A kind of method and system calculating data based on tree hierarchy partition structure maintenance power grid | |
CN201845205U (en) | Wind energy monitoring device based on singlechip and used in large-scale wind field | |
CN104900278A (en) | Simulating method for controlling pump flow during conversion between active system and passive system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |