CN103514967A - Intelligent reactor control method - Google Patents

Abstract

The invention relates to an intelligent reactor control method sequentially comprising the steps that: 1, parameters are collected; 2, neutron fluence rate of the reactor is sampled; 3, the neutron fluence rate is subjected to average processing with a calculation time interval delta t, such that average neutron fluence rate phi of each time point is obtained; 4, an intermediate variable of the t time is obtained; 5, a reactivity of the t time is obtained; 6, a neutron fluence rate expectation value of a t+1 time is set, such that intermediate variable of t+1 second is obtained; 7, the virtual reactivity of the t+1 time is obtained; 8, a control rod adjustment amount of the t+1 time is obtained; and 9, steps of a stepper motor driving the control rod is obtained. With the method provided by the invention, the control manner completely satisfies an objective law of reactor power growth or attenuation. During actual application, a stable value can be rapidly achieved according to a preset operation curve, such that overshoot and ringing caused by over-adjustment are avoided. During stable operation, neutron fluence rate can be stabilized within 0.3% under control.

Description

Intelligent response heap control method

Technical field

The present invention relates to a kind of intelligent response heap control method, particularly relate to a kind of intelligent response heap control method for meeting kinetics equation reactor.

Background technology

Traditional reactor control theory is according to current industry control PID(ratio, integration, differential) regulate, present most of reactor control systems are all according to this theory design, and do not take the inherent characteristic of reactor capability increase and decrease into account.

Regardless of different kinds of any reactor, all can move according to reactor kinetics equation, according to reactor kinetics equation, can accurately calculate current reactivity, if can set power or the neutron fluence rate of next second, adopt identical computing method just can calculate reactive variable quantity, and then the regulated quantity of controlled rod, calculate and make control rod run to accurate desired location in real time, make the real power of next second can approach very much predetermined value.This control model and traditional control model are very different, classic method only regulates the lifting of control rod according to the difference of performance number and setting value, and new control model has default intelligence, make the power of reactor reach fast and be stabilized in setting value, thereby avoided overshoot overshoot and vibration.Therefore need badly a kind of novel intelligent response heap control method is provided.

Summary of the invention

The technical problem to be solved in the present invention is to provide a kind of intelligent response heap control method that makes reactor reach fast stationary value.

For solving the problems of the technologies described above, intelligent response heap control method of the present invention, comprises the following steps successively:

The first step, gathers following parameter;

λ _ithe disintegration constant of-i group delayed-neutron precursor; Unit is s ^-1; Totally 15 groups;

β _effeffective share that-delayed neutron is total;

β _ieff-i organizes effective delayed neutron fraction, totally 15 groups;

Λ-neutron generation time; Unit is s;

The outer neutron source strength of Q-; Unit is neutron/s;

γ _sthe relative worth of-outer neutron source;

The efficiency of ε-outer neutron source;

K _hthe step-length of the stepper motor of-driving control rod, unit is mm/ step;

Intermediate variable $A_i=e^{-{\mathrm{\λ}}_i\mathrm{\Δt}};$

Intermediate variable $B_i=a_i(1-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Intermediate variable $C_i=a_i(e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Wherein: Δ t is interval computing time; a _i=β _ieff/ β _eff;

the differential efficiency of-control rod; Unit is mk/mm;

Second step, samples to the neutron fluence rate of reactor;

The 3rd step, in computing time interval of delta t, averages processing to neutron fluence rate, obtains each average neutron fluence rate Φ constantly;

The 4th step, obtains t intermediate variable I constantly _i(t);

I _i(t)＝I _i(t-1)·A _i+Φ(t)·B _i-Φ(t-1)·C _i；

i=1～15；I _i(0)＝0；

The 5th step, obtains t reactive ρ (t) constantly;

$\mathrm{\ρ}\left(t\right)=\mathrm{\Λ}\·\frac{\mathrm{\Φ}\left(t\right)-\mathrm{\Φ}(t-1)}{\mathrm{\Φ}\left(t\right)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{\mathrm{\Φ}\left(t\right)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_i\left(t\right)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{\mathrm{\Φ}\left(t\right)};$

The 6th step, sets t+1 neutron fluence rate expectation value Φ constantly _v(t+1), obtain the t+1 intermediate variable I of second _vi(t+1);

I _vi(t+1)＝I _i(t)·A _i+Φ _v(t+1)·B _i-Φ(t)·C _i；

i=1～15；

The 7th step, obtains t+1 virtual reactive ρ constantly _v(t+1);

${\mathrm{\ρ}}_v(t+1)=\mathrm{\Λ}\·\frac{{\mathrm{\Φ}}_v(t+1)-\mathrm{\Φ}\left(t\right)}{{\mathrm{\Φ}}_v(t+1)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{{\mathrm{\Φ}}_v(t+1)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_{\mathrm{vi}}(t+1)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{{\mathrm{\Φ}}_v(t+1)};$

The 8th step, obtains the regulated quantity Δ H of t+1 control rod constantly;

$\mathrm{\ΔH}=\frac{[{\mathrm{\ρ}}_v(t+1)-\mathrm{\ρ}\left(t\right)]}{\frac{\mathrm{d\ρ}}{\mathrm{dH}}};$

The 9th step, obtains driving the step number N of the stepper motor action of control rod;

$N=\frac{\mathrm{\ΔH}}{K_h}.$

Δt＝1s。

The invention enables control mode to meet the objective law of reactor-up or decay completely, this rule is difficult to common PID(ratio, integration, differential) regulate and realize, in actual use, can reach fast stationary value according to default operation curve (comprise and increase or reduce power), overshoot and the concussion of having avoided overshoot to cause, when stable operation, the present invention can compensate external interference rapidly and accurately, and neutron fluence rate is stabilized in 0.3%.See Fig. 1.

Accompanying drawing explanation

Fig. 1 is neutron fluence rate lifting curve figure.

Embodiment

Below in conjunction with drawings and Examples, the present invention is further detailed explanation.

The basis of this patent is to have adopted stepper motor to drive control rod lifting, makes the control accuracy of control rod within 0.2mm, and so just making accurately to control reactor according to reactivity becomes possibility.

The present invention comprises the following steps successively:

The first step, gathers following parameter;

λ _ithe disintegration constant of-i group delayed-neutron precursor; Unit is s ^-1; Totally 15 groups;

β _effeffective share that-delayed neutron is total;

β _ieff-i organizes effective delayed neutron fraction, totally 15 groups;

Λ-neutron generation time; Unit is s;

The outer neutron source strength of Q-; Unit is neutron/s;

γ _sthe relative worth of-outer neutron source;

The efficiency of ε-outer neutron source;

K _hthe step-length of the stepper motor of-driving control rod, unit is mm/ step;

Intermediate variable $A_i=e^{-{\mathrm{\λ}}_i\mathrm{\Δt}};$

Intermediate variable $B_i=a_i(1-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Intermediate variable $C_i=a_i(e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Wherein: Δ t is interval computing time, preferred Δ t=1s; a _i=β _ieff/ β _eff;

the differential efficiency of-control rod, is obtained in advance by Physical Experiment; Unit is mk/mm;

Second step is sampled to the neutron fluence rate of reactor on time, and the sample frequency of ADC is 5000 times/second, and every 50ms uploads to PC after neutron fluence rate data are averaged; This sampling is circulation execution incessantly on time;

The 3rd step, in computing time interval of delta t, averages processing to neutron fluence rate, obtains each average neutron fluence rate Φ constantly; Preferred Δ t=1s;

The 4th step, obtains t intermediate variable I constantly _i(t);

I _i(t)＝I _i(t-1)·A _i+Φ(t)·B _i-Φ(t-1)·C _i；

i=1～15；I _i(0)＝0；

The 5th step, obtains t reactive ρ (t) constantly;

$\mathrm{\ρ}\left(t\right)=\mathrm{\Λ}\·\frac{\mathrm{\Φ}\left(t\right)-\mathrm{\Φ}(t-1)}{\mathrm{\Φ}\left(t\right)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{\mathrm{\Φ}\left(t\right)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_i\left(t\right)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{\mathrm{\Φ}\left(t\right)};$

The 6th step, sets t+1 neutron fluence rate expectation value Φ constantly _v(t+1), obtain the t+1 intermediate variable I of second _vi(t+1);

I _vi(t+1)＝I _i(t)·A _i+Φ _v(t+1)·B _i-Φ(t)·C _i；

i=1～15；

The 7th step, obtains t+1 virtual reactive ρ constantly _v(t+1);

${\mathrm{\ρ}}_v(t+1)=\mathrm{\Λ}\·\frac{{\mathrm{\Φ}}_v(t+1)-\mathrm{\Φ}\left(t\right)}{{\mathrm{\Φ}}_v(t+1)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{{\mathrm{\Φ}}_v(t+1)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_{\mathrm{vi}}(t+1)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{{\mathrm{\Φ}}_v(t+1)};$

The 8th step, obtains the regulated quantity Δ H of t+1 control rod constantly;

$\mathrm{\ΔH}=\frac{[{\mathrm{\ρ}}_v(t+1)-\mathrm{\ρ}\left(t\right)]}{\frac{\mathrm{d\ρ}}{\mathrm{dH}}};$

The 9th step, obtains driving the step number of the stepper motor action of control rod;

$N=\frac{\mathrm{\ΔH}}{K_h}.$

Claims (2)

1. intelligent response is piled control method, comprises the following steps successively:

The first step, gathers following parameter;

λ _ithe disintegration constant of-i group delayed-neutron precursor; Unit is s ^-1; Totally 15 groups;

β _effeffective share that-delayed neutron is total;

β _ieff-i organizes effective delayed neutron fraction, totally 15 groups;

Λ-neutron generation time; Unit is s;

The outer neutron source strength of Q-; Unit is neutron/s;

γ _sthe relative worth of-outer neutron source;

The efficiency of ε-outer neutron source;

K _hthe step-length of the stepper motor of-driving control rod, unit is mm/ step;

Intermediate variable $A_i=e^{-{\mathrm{\λ}}_i\mathrm{\Δt}};$

Intermediate variable $B_i=a_i(1-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Intermediate variable $C_i=a_i(e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}-\frac{1-e^{-{\mathrm{\λ}}_i\mathrm{\Δt}}}{e^{{\mathrm{\λ}}_i\mathrm{\Δt}}});$

Wherein: Δ t is interval computing time; a _i=β _ieff/ β _eff;

the differential efficiency of-control rod; Unit is mk/mm;

Second step, samples to the neutron fluence rate of reactor;

The 3rd step, in computing time interval of delta t, averages processing to neutron fluence rate, obtains each average neutron fluence rate Φ constantly;

The 4th step, obtains t intermediate variable I constantly _i(t);

I _i(t)＝I _i(t-1)·A _i+Φ(t)·B _i-Φ(t-1)·C _i；

i=1～15；I _i(0)＝0；

The 5th step, obtains t reactive ρ (t) constantly;

$\mathrm{\ρ}\left(t\right)=\mathrm{\Λ}\·\frac{\mathrm{\Φ}\left(t\right)-\mathrm{\Φ}(t-1)}{\mathrm{\Φ}\left(t\right)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{\mathrm{\Φ}\left(t\right)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_i\left(t\right)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{\mathrm{\Φ}\left(t\right)};$

The 6th step, sets t+1 neutron fluence rate expectation value Φ constantly _v(t+1), obtain the t+1 intermediate variable I of second _vi(t+1);

I _vi(t+1)＝I _i(t)·A _i+Φ _v(t+1)·B _i-Φ(t)·C _i；

i=1～15；

The 7th step, obtains t+1 virtual reactive ρ constantly _v(t+1);

${\mathrm{\ρ}}_v(t+1)=\mathrm{\Λ}\·\frac{{\mathrm{\Φ}}_v(t+1)-\mathrm{\Φ}\left(t\right)}{{\mathrm{\Φ}}_v(t+1)}+{\mathrm{\β}}_{\mathrm{eff}}\·[1-\frac{1}{{\mathrm{\Φ}}_v(t+1)}\underset{i=1}{\overset{15}{\mathrm{\Σ}}}{I}_{\mathrm{vi}}(t+1)]-\frac{{\mathrm{\ϵ\γ}}_{s}Q}{{\mathrm{\Φ}}_v(t+1)};$

The 8th step, obtains the regulated quantity Δ H of t+1 control rod constantly;

$\mathrm{\ΔH}=\frac{[{\mathrm{\ρ}}_v(t+1)-\mathrm{\ρ}\left(t\right)]}{\frac{\mathrm{d\ρ}}{\mathrm{dH}}};$

The 9th step, obtains driving the step number N of the stepper motor action of control rod;

$N=\frac{\mathrm{\ΔH}}{K_h}.$

2. intelligent response heap control method according to claim 1, is characterized in that: described Δ t=1s.

Images