для игрока 2
V2 = d − |
1 |
. |
|
(X * , Lm ) |
|||
|
|
(Последние соотношения выводятся при доказательстве теоре-
мы 2.2.)
всистеме ядерных реакторов
2.2.1.Физическая постановка задачи
Вусловиях широкого ввода ядерно-энергетических мощностей и растущего разуплотнения графиков нагрузок энергосистем использование атомных электростанций (АЭС) в базовой части графика нагрузок будет представлять все большие трудности. Поэтому часть электростанций вынуждена будет работать в соответствии с суточными и сезонными колебаниями потребности в электроэнергии.
При этом решение проблемы маневренности энергоблоков АЭС возможно только после проведения комплекса научно-исследова- тельских и опытно-конструкторских работ. Основные из них следующие:
1)снятие ограничений на число циклов пуск-остановка для всего оборудования АЭС;
2)улучшение схем пусков и остановок АЭС и снижение потерь тепла при расхолаживании блоков за счет рационального использования остаточного энерговыделения ядерного топлива и теплоаккумулирующей способности графитовой кладки;
3)решение проблемы выбора оптимального типа топлива, защитных оболочек для него и конструкций тепловыделяющих элементов с учетом работы реакторов при переменных нагрузках;
4)приведение в соответствие с требованиями эксплуатации систем управления и защиты реакторов и обеспечение реакторов необходимым запасом реактивности.
Первые три проблемы носят в основном теплотехнический и прочностной характер, в то время как последняя связана с физикой ядерного реактора и обусловлена нестационарным отравлением реактора ксеноном.
131
Обеспечение реактора оперативным запасом реактивности, позволяющим компенсировать нестационарное отравление ксеноном, приводит к снижению энерговыработки реактора. Потерю энерговыработки реактора, работающего в переменном графике нагрузки по сравнению с энерговыработкой при работе на номинальной мощности, можно связать с резервируемым запасом реактивности соотношением
Q = aρ ,
где a – темп выгорания, 1
кВт сут; ρ – запас реактивности, отн. ед.; Q = Qm −Q – разность между энерговыработкой при работе в базисном режиме ( Qm ) и режиме переменных нагрузок ( Q ), кВт сут.
Величина резервируемого запаса реактивности, а следовательно, и потеря энерговыработки реактора, зависит от требуемой длительности работы реактора на пониженной мощности W, от степени
снижения мощности реактора ε = W , где Wн – номинальная мощ-
Wн
ность реактора, от величины плотности потока нейтронов при работе на номинальной мощности и других характеристик реактора.
В рамках точечной модели запас реактивности, обеспечивающий работу реактора на пониженной мощности εWн в течение произвольного времени, после снижения мощности равен
ρ = xм − xр ,
где xм и xр – максимальная и равновесная концентрации ксенона,
нормированные на |
ν f Σ f |
; |
Σ f |
– макроскопическое сечение деле- |
|
σX |
|||||
|
|
|
|
ния активной зоны реактора; ν f – среднее число вторичных ней-
тронов на акт деления; σX – сечение поглощения ксенона, см2.
На АЭС, как правило, устанавливается несколько энергоблоков. В общем случае энергоблоки могут отличаться электрической мощностью, темпом выгорания, другими характеристиками и рабо-
132
тать независимо друг от друга. Условием, связывающим реакторы, входящие в состав АЭС, является выработка заданного количества электроэнергии атомной электростанцией.
При ограниченном оперативном запасе реактивности, например, для реакторов с непрерывной перегрузкой или в конце кампании для реакторов с дискретной перегрузкой возникает задача оптимизации распределения запасов реактивности в системе реакторов с целью минимизации суммарной потери энерговыработки при удовлетворении заданному графику суточного снижения мощности.
2.2.2. Математическая постановка и решение оптимизационной задачи
Математически задача формулируется следующим образом: найти
|
|
|
N |
ρ |
i |
|
min |
∑ |
|
|
|||
ai |
|
|||||
ρ1 ... |
ρN i=1 |
|
|
|||
при ограничениях: |
|
|
|
|
|
|
N |
|
|
|
|
|
|
∑δi εi ( |
ρi ) = α; |
|
|
|
||
i=1 |
|
|
|
|
|
|
0 ≤ ρi |
≤ |
|
ρi m , i =1, ..., N, |
|||
где N – число реакторов на станции; Δρi – оперативный запас реактивности i-го реактора; εi – степень снижения мощности i-го реак-
N
∑Wi
тора; α – заданная степень снижения мощности АЭС, α = Ni=1 ;
∑Wi н
i=1
δi – доля электрической мощности i-го реактора, δi = NWi н .
∑Wi н
i=1
Заметим, что оптимизация возможна только при < 0 < α <1 , так как при α = 0 и α =1 значения оперативных запасов реактивности определены:
133
ρi = ρi м,
ρi = 0.
Поэтому качественно понятно, что эффект от оптимизации будет сильнее всего сказываться в средней части диапазона изменения α . Вид зависимости ε( ρ) возможной степени снижения
мощности реактора от запаса реактивности, вообще говоря, определяется режимом изменения мощности реактора.
Наиболее простым режимом для анализа является снижение мощности реактора с максимально возможной скоростью до определенного уровня и поддержание реактора на данном уровне до момента выхода АЭС на номинальную мощность. Характер зависимости ε( ρ) при различных плотностях потоков нейтронов в
этом случае показан на рис. 2.7.
Возможная степень снижения мощности,ε
|
ϕ=1. 1013H/см2с |
1,0 |
ϕ=3. 1013H/см2с |
0,8 |
ϕ=5. 1013H/см2с |
|
|
0,6 |
|
0,4 |
|
0,2 |
|
0,0
0 |
1 |
2 |
3 |
4 |
Величина резервируемого запаса пеактивности, Δρ%
Рис. 2.7. Зависимость возможной степени снижения мощности от величины резервируемого запаса реактивности
134
Как видно из рисунка, при малых плотностях потока нейтронов зависимость ε( ρ) представляет собой практически линейную
функцию. С увеличением плотности потока нейтронов нелинейность в характере зависимости ε( ρ) увеличивается.
В большинстве энергетических реакторов уровень плотности
потоков нейтронов больше, чем 1 1013 нейтр. при этом зависисм2 с
мость ε( ρ) , как видно из рисунка, носит нелинейный характер.
Для получения качественных результатов рассмотрим систему из двух реакторов. Зависимость ε( ρ) будем аппроксимировать по-
линомом второй степени.
|
|
|
|
|
|
|
|
N |
|
ρ |
i |
|
|
|
|
|
|||||
|
|
min |
|
∑ |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
ai |
|
|
|
|
|||||||||||||||
|
ρi ... |
ρN i=1 |
|
|
|
|
|
|
|||||||||||||
при ограничениях |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
ρ |
i |
|
|
|
|
ρ |
i |
|
|
2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
∑δi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= α; |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
1 − 2 |
|
ρim |
+ |
|
|
|
|
|
|
|
|
||||||||||
i=1 |
|
|
|
|
|
ρim |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
0 ≤ |
ρ1 ≤ |
ρ1m ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.88) |
||||
0 ≤ |
ρ2 ≤ |
|
ρ2m . |
|
|
|
|
|
|
|
|
|
|
||||||||
Введя новые переменные |
|
|
|
|
|
|
ρ1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
z |
= |
|
|
|
|
, |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
1 |
|
|
|
ρ1m |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
z2 = |
|
|
ρ2 |
, |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
ρ2m |
|
|
|
|
|
|
|
|||||
запишем задачу (2.88) в виде: найти |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
z1 |
|
|
|
|
|
|
|
z2 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
+ |
|
|
|
|
|
|
|||||||||||
|
min |
|
|
|
|
|
a2 ϕ |
|
|
|
|||||||||||
|
z1,z2 a1 ϕ1 |
2 |
|
||||||||||||||||||
135