СПОСОБ ОПРЕДЕЛЕНИЯ ОБЛАСТИ ПРОТЕЧКИ РАДИОАКТИВНОГО АЗОТА В ПАРОГЕНЕРАТОРАХ ЯДЕРНЫХ РЕАКТОРОВ ТИПА КЛТ-40
METHOD FOR DETERMINING THE AREA OF LEAKAGE OF RADIOACTIVE NITROGEN IN STEAM GENERATORS OF NUCLEAR REACTORS OF THE KLT-40 TYPE

Изобретение относится к области радиационного контроля и обеспечения радиационной безопасности объектов использования атомной энергии и может применяться для обнаружения области протечки радионуклидов и оценки ее величины при эксплуатации водо-водяных ядерных реакторов [G01T 1/16, G01T 1/18, G01T 1/20].
Из уровня техники известен СПОСОБ УПРАВЛЕНИЯ ОСТАНОВОМ ВОДО-ВОДЯНОГО ЯДЕРНОГО РЕАКТОРА [RU 2706739 C2, опубл.: 20.11.2019], характеризующийся тем, что обнаружение утечки первого/второго контура осуществляют путем обнаружения одного или нескольких из следующих признаков:
- повышенная радиоактивность второго контура из-за загрязнения теплоносителем первого контура;
- повышенный запас воды во втором контуре;
- пониженный запас воды в первом контуре.
Недостатком аналога является невозможность определения области протечки в контурах реактора, а также низкая точность определения наличия протечки.
Также известно УСТРОЙСТВО И СПОСОБ ДЛЯ ОБНАРУЖЕНИЯ УТЕЧКИ [RU 2197718 C2, опубл.: 27.01.2003], в котором место утечки определяют, например, при известной скорости течения в коллекторном трубопроводе из промежутка времени, которое прошло между пиком давления в коллекторном трубопроводе и срабатыванием детектора для регистрации радиоактивности. Пик давления может, например, быть вызван открыванием, по меньшей мере, одного обратного клапана.
Недостатком аналога является невозможность определения области протечки в коллекторном трубопроводе.
16
7
Наиболее близкой по технической сущности является МОДЕЛЬ НАКОПЛЕНИЯ РАДИОНУКЛИДОВ В КОТЛОВОЙ ВОДЕ ПАРОГЕНЕРАТОРОВ АЭС С ВВЭР-440 И ВВЭР-1000 [Иванов Е.А., Пырков И.В., Хамьянов Л.П. / Атом. энергия. - 1994. - Т. 77, N 1. - С. 58-63, 85. - ISSN 0004-7163] и МЕТОДИКА ДИАГНОСТИКИ ПРОТЕЧКИ ТЕПЛОНОСИТЕЛЯ ПЕРВОГО КОНТУРА В КОТЛОВУЮ ВОДУ ПАРОГЕНЕРАТОРОВ АЭС С ВВЭР-440 И ВВЭР-1000 [Иванов Е.А., Пырков И.В., Хамьянов Л.П. / Атом. энергия. - 1994. - Т. 77, N 1. - С. 51-58, 85. - ISSN 0004-7163], физический смысл которых состоит в оценке проникновения указанных радионуклидов из первого контура в котловую воду парогенератора второго и степени их накопления. Модель предусматривает естественную линейную зависимость объемной активности указанных радионуклидов от мощности реактора. Определение накопленной активности указанных радионуклидов котловой воды парогенератора осуществляется путем измерения активности фильтров, но радиоактивный газ, например, радиоактивного азота N фильтрами не задерживается и выходит с паром, а изотопы йода остаются в воде (более чем на 99%).
16
7
Основная проблема прототипа состоит в том, что в его рамках удается лишь констатировать наличие протечки по активности изотопов, осевших на фильтры, и активности радиоактивного азота в помещении машинного зала при выходе пара на турбину. В лучшем случае можно оценить плотность потока гамма-квантов, бета-излучения и объемную активность радионуклида, используя показания сцинтилляционных детекторов. В указанном методе отсутствует возможность получения информации об области (на спиральном паропроводе) возникновения протечки, динамики ее развития, которая будет, очевидно, линейно зависеть от изменения мощности реактора, ее размерах и мощности «вброса» радиоактивного азота N из первого контура во второй. Отсутствие этой информации не дает возможности однозначно ответить на вопрос по какой причине происходит протечка, провести анализ металла водо-паропровода, его физико-механических свойств и их изменений при воздействии ионизирующего излучения и частой смене температурного режима в области протечки, изучить особенности способствующие возникновению микротрещин, через которые происходит протечка, с целью предотвращения подобных эффектов, а также прогнозировать величину внешней и ингаляционной дозы, которую может получить персонал, обслуживающий парогенератор и турбину, что приводит к нарушению условий радиационной безопасности в помещении машинного зала.
131
131
Исходя из этого в настоящем время запрещается эксплуатация парогенераторов при достижении суммарной удельной активности радионуклидов в продувочной воде 370 Бк/кг. Допустимая удельная активность I в продувочной воде каждого парогенератора не более 740 Бк/кг, допустимая средняя удельная активность I всех парогенераторов блока не должна превышать 185 Бк/кг. При этом протечка теплоносителя первого контура в котловую воду каждого парогенератора должна быть не более 5 кг/ч. В связи с этим важное значение приобретает разработка метода оценки протечек, в частности, разработка корректной модели, позволяющей найти причины протечки, ее область в парогенераторе и, таким образом, ответить на поставленные вопросы.
Задачей изобретения является устранение недостатков прототипа.
Технический результат изобретения заключается в обеспечении возможности обнаружения области протечки радиоактивного азота в парогенераторах ядерных реакторов КЛТ-40.
в
п
L
7
N
L
16
Указанный технический результат достигается за счет того, что способ определения области протечки радиоактивного азота в парогенераторах ядерных реакторов типа КЛТ-40, характеризующийся тем, что область протечки в трубопроводе парогенератора определяют путем определения плотности распределения воды ρ(Т), поступающей в водопровод парогенератора как функции температуры при её переносе вдоль оси X по водопроводу, плотности пара ρ(Т) как функции температуры или расстояния x при его переносе вдоль оси X по паропроводу при спрямлении его спиральной части, при этом эффективную длину водопровода, в пограничной области которой возникла протечка, определяют как расстояние x, при котором давление воды и пара в водо-паропроводе парогенератора уравновешено, активности радионуклида Nв области «вброса» из первого контура во второй, а также площади «вброса», отличающийся тем, что наличие радионуклида в паре определяют по измеренному значению мощности дозы ионизирующего излучения и плотности пара в точке выхода пара на турбину через приведенный размерный коэффициент α, характеризующий содержание радиоактивного азота в единице массы пара в точке выхода пара на турбину при его прохождении по паропроводу длиной L, расстояние x определяют путем решения уравнения
R
L
L
L
P
п
0
в
0
L
0
гр
L
где  – газовая постоянная, ρ(), ρ() плотности пара и воды в точке x≡ соответственно, – давление воды, а площадь «вброса» в области протечки находят по разности расстояний x, возникающих в результате флуктуации давления воды в области протечки, при которых также имеет место соответствующее равенство давлений воды и пара.
N
sf
γ
V
sf
γ
7
V
п
l
3
16
16
3
В частности, приведенный размерный коэффициент αопределяют по формуле , где D' = 2πKQ(L) – измеряемая величина мощности дозы в выделенном сферическом объеме радиоактивного пара V = π/6 диаметром l, K – γ-постоянная азота N, Q(L) – виртуальная объемная активность радионуклида N (Ки/м) на выходе из паропровода длиной L, возникающая за промежуток времени τ.
N
В частности, приведенный размерный коэффициент αопределяют, измеряя мощность дозы γ-β-излучения на выходе пара из паропровода, количество воды, образующейся при конденсации этого пара на выходе паропровода, температуру этого пара на выходе и определяя объем его выброса.
Краткое описание чертежей.
На фиг. 1 схематично изображена ядерно-энергетическая установка.
На фиг. 2 показана зависимость плотности воды от температуры.
На фиг. 3 показана зависимость скорости испарения водяного пара дистиллированной воды от температуры.
На фиг. 4 схематично изображен парогенератор с водо-паропроводом и протекающим теплоносителем.
На фиг. 5 показана зависимость удельной теплоты парообразования воды Λ(T) (МДж/кг) от температуры Т °C.
На фиг. 6 показана температурная зависимость скорости переноса пара в парогенераторе реактора КЛТ-40.
п
п
На фиг. 7 показана температурная зависимость отношения φ(Т) = u(Т)/v(Т): 1 – график отношения функций; 2 – аппроксимация отношения φ(Т) параболой в температурном диапазоне 170-290  °C.
п
п
На фиг. 8 показана зависимость скорости испарения водяного пара дистиллированной воды u(Т) от температуры: 1 – график функции u(Т); 2 – аппроксимация функции параболой в температурном диапазоне 170-290  °C.
0
гр
На фиг. 9 показана расчетная длина водного участка L в парогенераторе второго контура реактора КЛТ-40 при давлении Р: 1 – 6,56; 2 – 6,4; 3 – 6,24 МПа.
На фигурах обозначено: 1 – парогенератор, 2 – турбина, 3 – генератор, 4 – водопровод, 5 – паропровод, 6 – насос, 7 – конденсатор, 8 – защитная оболочка, 9 – стержни регулирования, 10 – трубопровод первого контура, 11 – реактор.
Осуществление изобретения.
Ядерная энергетическая установка – это устройство для получения тепловой, электрической или механической энергии в ходе управляемой ядерной реакции, осуществляемой в ядерном реакторе. Ядерная энергетическая установку, наряду с одним или несколькими ядерными (атомными) реакторами, включает парогенератор 1 (см. Фиг. 1), паровую турбину 2, приводимую ей в действие электрический генератор 3, а также трубопровод 4, паропровод 5, насос 6, конденсатор 7 и другое вспомогательное оборудование. В современных стационарных ядерных энергетических установках в качестве рабочего тела применяется в основном водяной пар, поэтому все теплообменники первого контура двухконтурных схем являются парогенераторами 1. Трубопровод первого контура 10 проходит через реактор 11, в защитной зоне которого, закрытой защитной оболочкой 8 смонтированы стержни регулирования 9 и сообщен с парогенератором 1 второго контура.
в
в
п
п
Общий принцип работы парогенератора 1 ядерной энергетической установки состоит в том, что во второй контур парогенератора 1 по трубопроводу 10 с помощью насоса 6 подают воду под давлением P, с температурой Т. Тепло от воды, циркулирующей в первом контуре, передают воде второго контура, тем самым нагревают и подают ее в парогенератор 1, где с образованием пара, выход которого через N спиральных паропроводов 5 того же внутреннего диаметра осуществляется с температурой Т под высоким давлением P. В процессе прохождения воды в парогенераторе 1 воду нагревают до температуры насыщения пара при соответствующем давлении, которая затем испаряется на внутренней поверхности трубок парогенератора 1, создавая эффективную пограничную область вода-пар, и, наконец, в виде перегретого пара по паропроводу 5 подают на турбину 2. Таким образом, на вход турбины 2 подают пар высокого давления, перегретый относительно температуры насыщения. Зависимость температуры кипения воды (парообразования) от её давления приведена в табл. 1.
Таблица 1
P
t<sub>k</sub>,°C
P
t<sub>k</sub>,°C
кПа
атм.
кПа
атм.
0,981
0,01
6,698
196,1
2,0
119,62
1,961
0,02
17,20
245,2
2,5
126,79
3,923
0,05
28,64
294,2
3,0
132,88
9,807
0,1
45,45
392,3
4,0
142,92
19,61
0,2
59,67
490,3
5,0
151,11
29,42
0,3
68,68
588,4
6,0
158,08
39,23
0,4
75,42
686,5
7,0
164,17
49,03
0,5
80,86
784,5
8,0
169,61
58,84
0,6
85,45
882,6
9,0
174,53
68,65
0,7
89,45
980,7
10,0
179,04
78,45
0,8
92,99
1961
20,0
211,38
88,26
0,9
96,18
2452
25,0
222,90
98,07
1,0
99,09
4903
50,0
262,70
101,3
1,033
100,00
9807
100,0
309,53
147,1
1,5
110,79
-
-
-
16
7
Автоматическое регулирование процесса подпитки водой парогенератора 1 приводит к повышению давления воды в водопроводе 4 и обратному смещению воды в первоначальное положение. В пограничной области на водопроводе 4 и паропроводе 5 температурный режим будет постоянно изменяться, что приводит к аналогичному изменению частоты механических напряжений водопровода 4 и паропровода 5 в этой области, последующей усталости металла и к вероятному появлению микротрещин, через которые из первого контура во второй 10 может проникать радионуклид N, создавая протечку. Поскольку плотность пара существенно меньше плотности воды, то радиоактивный азот будет распространяться в область паровой фазы, включая и выход пара на турбину 2.
16
16
16
7
7
7
в
п
0
0
0
0
L
L
0
Новым в изобретении является определение области протечки в водопроводе 4 и паропроводе 5, секундного расхода «вброса» радионуклида Nиз первого контура во второй, скорости «вброса», площади области «вброса», а также определение активности радионуклида Nв области «вброса» из первого контура во второй в заданной области паропровода 5, при этом наличие радионуклида в паре с учётом радиоактивного распада, в заданной точке (области) паропровода 5, определяют по измеренному значению плотности потока гамма, бета-излучения азота N или мощности дозы ионизирующего излучения, создаваемой указанным радионуклидом в точке выхода пара на турбину 2, для оценки которых определяют плотность распределения воды ρ(Т), поступающей в водопровод 4 парогенератора 1 , как функции температуры при её переносе вдоль оси x по водопроводу 5 (см. Фиг. 2), плотности пара ρ(Т) как функции температуры (см. Фиг. 3) или расстояния x при его переносе вдоль оси x по паропроводу 5 (см. Фиг.4) при спрямлении его спиральной части длиной L с внутренним радиусом R=1,5–2,0 см (0 ≤ r ≤ R) при длине водопровода 4 L < L (вертикальная область в левой части) и общей длине водопровода 4 и паропровода 5 L+L = Ls составляющей для парогенератора 1 от 10 до 20 м. Принимая, что рост температуры по водопроводу 4 и паропроводу 5 пропорционален расстоянию x, то при некотором x давление пара и воды может быть уравновешено. По значению x = L определяют эффективную длину водопровода 4, в пограничной области которой будет возникать протечка.
16
16
16
7
7
гр
7
½
1
1
2
2
Рассмотрим метод оценки протечки радиоактивного азота Nна примере парогенератора 1 реактора КЛТ-40. В теплоносителе (воде) первого контура находятся радионуклиды, включая и азот N, которые возникают в теплоносителе в результате физико-химического взаимодействия последнего с тепловыделяющими элементами (твелами) реактора 11. Этот теплоноситель служит источником тепла для спирального водопровода 4 и паропровода 5 парогенератора 1, в который подают воду под давлением P ≈ 6,1 ‒ 6,4 МПа и одновременно является источником радионуклидов, поступающих во второй контур. Вода, распространяющаяся по спиральному водопроводу 4 парогенератора 1, нагревается, что приводит к уменьшению её плотности, испарению и образованию насыщенного пара, давление которого увеличивается с ростом температуры, т.е. с ростом координаты x. Этот процесс происходит с определенной частотой. При этом в паропроводе 5 в силу разной плотности сред возникает область раздела пар-вода, которая в зависимости от флуктуаций давления воды или пара смещается от равновесного положения в ту или другую сторону. Изменение температурного режима металла паропровода 5 приводит к росту механических напряжений, в результате которых в этой области могут возникнуть микротрещины, через которые и могут проникать радионуклиды из теплоносителя первого контура во второй. Далее перемещаясь по паропроводу 5 пар парогенератора 1 в виде струи, поступающей на турбину 2, содержит радиоактивный азот N с периодом полураспада Т = 7,11 с, энергией гамма-излучения Emax = 6,134 MэВ и квантовым выходом ν, max = 69%; энергиями бета-излучения E = 4,288 МэВ, и выходом бета-частиц n = 68%; E = 10,419 МэВ, n = 26%, содержание которого в паре свидетельствует о признаке нарушения герметичности водопровода 4 или паропровода 5 второго контура парогенератора 1. Содержание радиоактивного азота при выходе пара на турбину 2 определяют с помощью сцинтилляционных детекторов с кристаллом NaJ(Tl).
в
в
0
0
2
Перенос воды в водопроводе 4 осуществляется с заданной скоростью, определяемой начальным давлением и плотностью, уменьшением плотности воды с ростом температуры и её испарением. В водопроводе 4 и паропроводе 5 небольшим радиусом масштаб турбулентности ограничен и не может быть больше внутреннего радиуса трубки, поэтому в уравнениях, описывающем перенос воды и пара турбулентность среды не учитывается. При этом, полагая, что поперечная скорость переноса воды на стенки трубопровода 4 равна нулю, стационарное уравнение, учитывающее эти физические процессы для частицы воды массой m = ρ(T)V, занимающей объём V = πRL, в цилиндрической геометрии будет выглядеть как:
,
=0
x
где  – средняя скорость переноса воды по водопроводу 4; G = const – генерация воды, нагнетаемой под давлением, [кг/с], определяемая в точке T|x = T выражением:
,
п
0
0
2
u(Т) – скорость испарения воды или скорость генерации пара, [кг/мс]; Sв = 2πRL – площадь испарения воды.
в
в
в
в
0
0
t
s
0
t
s
Площади испарения S воды и занимаемый ею объём V являются характеристиками одного и того же объёма среды, то их отношение S/V = 2/R. Спираль паропровода 5 представляют в виде горизонтального участка (Фиг. 4), а изменение температуры в среде второго контура парогенератора 1 аппроксимируют линейной функцией Т(x) = T + bx/L °С, в которой T = 170 °C, b = 120 °С, а x удовлетворяет неравенству 0 ≤ x ≤ L.
Скорость переноса воды массой m по водопроводу 4 определяют из условия равенства кинетической и потенциальной энергий давления воды на входе второго контура парогенератора 1:
где N – число водопроводов 4 в пакете.
в
В этом случае учитывают свойство воды как несжимаемой жидкости и считают ее скорость v в каждой трубке водопровода 4 постоянной величиной.
Окончательно среднюю скорость переноса воды в водопроводе 4 определяют через среднюю плотность воды из выражения:
в
в
0
t
s
2
3
где зависимость плотности воды от температуры ρ(Т) аппроксимируется параболой ρ(T) = а·T + b·T + c, [кг/м], где а, b, c – коэффициенты соответствующей размерности равные a = -0,002315; b = -0,39345; c = 1033,79, а изменение температуры с расстоянием x –Т(x) = T + bx/L.
в
0
в
0
t
s
3
По графику зависимости плотности воды ρ от температуры (см.Фиг.2) при температуре Т плотность воды ρ, которая изменяется как функция расстояния x в соответствии с принятой линейной зависимостью температуры Т от расстояния x, определяемой линейной функцией Т(x) = T + bx/L °С в водопроводе 4 и паропроводе 5 составит 730 и 900 кг/м при температуре 170 и 290 °С соответственно.
Выражение для среднего значения плотности будет выглядеть как:
п
п
Зависимости скорости испарения (скорости генерации пара) как функции температуры u(T) получают по формуле Клапейрона-Клазиуса, согласно которой резкий спад u(T) в области высоких температур обусловлен уменьшением удельной теплоты парообразования Λ(T), также зависящей от температуры (см. Фиг. 5).
,
0
где T – начальная температура.
п
0
п
0
3
-2
-1
Принимают для оценки u(T) в диапазоне температур 170-290°С T = 72 °С и u(T) = 1,375·10 кг·м·с (см. Фиг. 6).
Из фиг.2 и 3 следует, что в температурном диапазоне от 100 до 220 °C скорость испарения воды с ростом температуры резко увеличивается, а ее плотность, напротив, быстро уменьшается. Тогда, полагая постоянной скорость переноса воды по водопроводу 4, поскольку её перенос происходит без трения, которому способствует процесс парообразования на внутренней поверхности трубок парогенератора 1, давление воды с ростом x, т.е. с ростом температуры, будет также уменьшаться пропорционально плотности воды, а пара, напротив, возрастать.
Плотность воды определяют по формуле:
С
0
в
где постоянная  определяется через значение плотности воды из фиг. 2 при Т = Т, x = 0 или через параболическую зависимость плотности воды от температуры ρ(T):
п
п
п
п
п
0
2
При описании процесса скорости переноса пара также не учитывают скорость турбулентной диффузии пара, поскольку она значительно меньше его скорости, обусловленной адвективной составляющей и, обозначив массу пара через m = ρ·V, где V – внутренний объём паропровода 5 (V = πR·L), уравнение переноса, для плотности пара по паропроводу при его поперечной скорости равной нулю для стационарного процесса, определяют выражением:
п
п
п
0
0
п
п
где v – продольная скорость переноса пара по паропроводу 5; S – площадь парообразования (S ≈ 2πRL); τ = L/v – время «жизни» пара в паропроводе 5.
п
п
В последнем выражении скорости как переноса пара v, так и его генерации (при испарении воды) u зависят от его температуры и давления. В качестве оценки скорости переноса пара используют формулу Сен-Венана, определяющей истечение газа из резервуара в атмосферу как функции его температуры и давления:
,
п
ат
п
где k – постоянная адиабаты водяного пара; R – газовая постоянная (Дж/кг °K); T(L) – температура пара на выходе из паропровода (°K); P – атмосферное давление; P – давление пара на выходе из паропровода 5, МПа.
Общий график зависимости v(T) приведен на Фиг. 6.
п
0
0
Значение плотности пара в пограничной области ρ(L), т.е. при x = L, находят используя условия равенства давления воды и пара в этой области при заданной температуре, т.е.:
0
Рассматривая пар как идеальный газ и используя уравнение Клапейрона-Менделеева, граничное условие в области x = L получают в виде:
п
п
п
п
п
t
s
п
п
Поскольку плотность пара является сложной функцией температуры, обусловленной температурной зависимостью скорости испарения воды u(Т) и скоростью переноса пара в паропроводе 5 v(Т), которые, в свою очередь, зависят от координаты x, тоdρ/dx = (dρ/dT)(dT/dx) = (dρ/dT) (b/L), и учитывая, что отношение указанных функций φ(Т) = u(Т)/v(Т) также зависит от температуры, уравнение для плотности пара при условии стационарного процесса переписывают в виде:
где функция φ(Т) также аппроксимируется параболой:
f
f
f
-3
с относительной погрешностью температурной зависимости в интересующем диапазоне не хуже 1% (см. Фиг. 7), где а = -0,5524·10; b = 0,2717515; c, = -24,13 445.
При этих условиях решение уравнения для пара, определяющее его плотность, принимает вид:
s
t
s
0
где α = L/[b·(L - L)].
п
0
0
п
0
0
0
0
0
t
s
п
0
L
0
0
3
L
Для оценки начальной плотности пара ρ(Т) T(°K) используют уравнение Клапейрона – Менделеева и данными рис. 4, что дает ρ(Т) = P(T+273)/RT = 4,135 кг/м. Таким образом, для x = Lв соответствии с Т(x) = T + bx/L и решением уравнения для пара, получают выражение для значения плотности пара ρ(L) в пограничной области, что позволяет записать окончательное уравнение для граничного условия в области x≡  и, в конечном итоге, значение параметра L, определяющего координату x возможной протечки в паропроводе парогенератора выражают как:
п
п
п
п
Температурная зависимость функций u(Т), v(Т) и их отношения φ(Т) = u(Т)/v(Т) в диапазоне 170 ≤ Т ≤ 290 °С приведены на Фиг. 8, 6, 7 соответственно.
Из Фиг. 8 следует, что рост температуры определяет две области, в которых скорости генерации пара существенно различаются. В диапазоне температур 30 ≤ Т ≤ 240 °С скорость генерации с ростом температуры растет, достигая своего максимума, а затем уменьшается, что обусловлено резким спадом удельной теплоты парообразования воды Λ(T) в области температур 230 ≤ Т ≤ 380 °С (см. Фиг. 5). Поэтому наиболее значимой областью генерации пара из выделенной области температур является диапазон 440-530 °K (167-257 °С).
п
п
п
п
Далее вычисляют интеграл от скорости парообразования u(Т), значение которого находят, используя аппроксимацию параболой указанную кривую (см. Фиг. 8) в интересующем диапазоне температур 170 ≤ Т ≤ 290 °С. Чтобы упростить вычисления находят коэффициенты a = -0,7921; b=401,0046; c= -37083,2407, что позволяет провести расчет с точностью не хуже 1% и получить расчетную функцию:
.
0
l
0
r
0
Решение уравнения для граничного условия в области x = L находят графическим методом, строя графики зависимостей функций левой: U(L) и правой U(L) его частей, определяемых соответственно выражениями:
r
0
гр
0
Функция U(L) представляет собой, практически, постоянную, зависящую от параметра давления воды P. Поэтому абсцисса точки пересечения графиков и определяет величину параметра L – длину водного участка парогенератора 1.
s
l
0
гр
s
r
0
0
0
0
7
7
Результаты расчетов для парогенераторов 1, различающихся общей длиной L, приведены на Фиг. 9, на которой кривые 1, 2, 3 определяют зависимость функции U(L) при P= 6,4 МПа и общей длине парогенератора 1 L равной 20 м (кривая 1), 15 м (кривая 2) и 10 м (кривая 3), а функция U(L) (кривая 4) представляет собой величину, слабо изменяющуюся с ростом L от 2,098·10 (L = 0,1) до 2,101·10 (L = 0,9063).
гр
l
0
r
0
6
С ростом давления воды P во втором контуре значение функции U(L) также увеличивается и, напротив, уменьшается с падением давления (см. Табл. 2), значение постоянной U(L) также изменяется соответствующим образом, но характер решения остается прежним. Изменение давления определяют его флюктуацией и составляет ±2,5% от 6,4х10Па.
Таблица 2
Длина водного участка парогенератора <i>L</i><sub>0</sub>, м
<i>P</i>гр, МПа
6,24
6,4
6.56
<i>Ls</i>, м
10
0,7534
0.7636
0,7739
20
0,8944
0,9063
0,9124
Максимальная разность «хода» воды Δ<i>L</i><sub>0</sub>, м
<i>Ls</i>, м
10
0,0205
20
0,024
16
16
16
7
N
п
п
п
N
N
Поскольку период полураспада N весьма мал, то при формулировке уравнения переноса радиоактивного азота N по паропроводу 5 рассматривают его во времени, учитывая локальную генерацию в пограничной области вода-пар, уменьшение в результате радиоактивного распада и, кроме того, наличие радиоактивного азота в объёме пара, распространившегося в последнем в предыдущие моменты времени. При этом величину плотности пара определяют решением уравнения для плотности пара. Значение активности радиоактивного азота N, содержащегося в паровой фазе, и выходящего на турбину 2 Q(L), определяют путём её измерения на выходе за определённый промежуток времени τ, за который определяют и выход пара m (воды). Тогда отношение Q(L)/m(L) = α, определит размерный коэффициент α (Ки/кг), в области выхода пара на турбину 2, т.е. при x = L.
16
7
Генерация радиоактивного азота может быть представлена мощностью «вброса» Nв паровую фазу, которую определят следующим выражением:
в
V
N
P = Q G,
N
V
3
3
где G – секундный расход «вброса» [м/с]; Q – объёмная активность радиоактивного азота [Ки/м].
16
16
16
7
0
А
7
N
пг
0
vr
vr
0
г
Если объёмную активность N измерить на выходе из паропровода 5, то в его начальной точке (x = L), т.е. в области его генерации, значение начальной объёмной активности Qможет быть найдено с поправкой на радиоактивный распад N. Секундный расход G «вброса» радиоактивного азота N в паровую фазу в пограничной области вода-пар представляет собой произведение площади пограничной области паропровода 5S = 2πRl , где l = ΔL – ширина пограничной области вода-пар (см. Табл. 2), на скорость выхода радиоактивного азота из трещин U, являющейся искомой величиной:
N
0
vr
г
G = 2πRl U.
сп
сп
0
пг
пг
0
сп
сп
in
ex
in
ex
2
Если водо-паропровод изготовлен в виде спирали радиусом R с шагом h, то при длине пограничной области вода-пар равной ΔL, площадь этой области S определится выражением: S=(ΔL/h)·4π·R(R+R)/2, где R,R – внутренний и внешний радиусы трубки водо-паропровода соответственно.
vr
Величину l оценивают по показанию расходомера, т.е. по величине продвижения «хода» воды в трубе водопровода 4 на его начальном участке при флуктуации в нём давления воды, которое также может быть найдено путем измерения давления воды по показанию манометра и использованием приведенного расчетного алгоритма.
гр
vr
0
0
max
0
min
V0
16
3
16
При измерении давления воды с учетом ее флуктуации ±P величина l = ΔL = L, - L, будет зависеть от разности измеряемых давлений воды. При этом значение объёмной активности радиоактивного азота N Q (Ки/м) при его генерации через микротрещины, возникающие в пограничной области вода-пар паропровода 5 определяется объёмной активностью в объёме парогенератора 1. Исходя из этого, перенос радиоактивного азота N в паропроводе 5 описывают следующим уравнением:
0
где L – длина водопровода 4.
16
16
vr
0
В приведенном уравнении первый член в правой части описывает генерацию радиоактивного азота N в пограничной области вода-пар шириной l << L; η(x) – единичная функция; второй – убыль радиоактивного азота, содержащегося в паре, за счёт выноса последнего по паропроводу 5 на турбину 2, третий – убыль радионуклида за счёт радиоактивного распада азота N с постоянной распада λ. Вынос радиоактивного азота с паром из паропровода 5, представленного вторым членом в уравнении, зависит от адвективной скорости, определяемой формулой Сен-Венана, определяющей истечение газа из резервуара в атмосферу как функции его температуры и давления и для азота, должно быть несколько отличным от пара, поскольку постоянная адиабаты в упомянутой формуле для азота и пара отличаются (см. Табл. 3).
Таблица 3
Показатели адиабаты <i>k </i>для различных температур и газов
темп.,°C
газ
<i>k</i>
темп.,°C
газ
<i>k</i>
темп.,°C
газ
<i>k</i>
−181
H<sub>2</sub>
1,597
200
сухой воздух
1,398
20
NO
1,400
−76
1,453
400
1,393
20
N<sub>2</sub>O
1,310
20
1,410
1000
1,365
−181
N<sub>2</sub>
1,470
100
1,404
2000
1,088
15
1,404
400
1,387
0
CO<sub>2</sub>
1,310
20
Cl<sub>2</sub>
1,340
1000
1,358
20
1,300
−115
CH<sub>4</sub>
1,410
2000
1,318
100
1,281
−74
1,350
20
He
1,660
400
1,235
20
1,320
20
H<sub>2</sub>Oнасыщенный водяной пар
1,330
1000
1,195
15
NH<sub>3</sub>
1,310
100
1,324
20
CO
1,400
19
Ne
1,640
200
1,310
−181
O<sub>2</sub>
1,450
19
Xe
1,660
−180
Ar
1,760
−76
1,415
19
Kr
1,680
20
1,670
20
1,400
15
SO<sub>2</sub>
1,290
0
сухой воздух
1,403
100
1,399
360
Hg
1,670
20
1,400
200
1,397
15
C<sub>2</sub>H<sub>6</sub>
1,220
100
1,401
400
1,394
16
C<sub>3</sub>H<sub>8</sub>
1,130
k
ср
2
2
7
7
16
16
16
Значение постоянной адиабаты для заданного диапазона температур получают методом экстраполяции, используя полученную параболическую зависимость k(T), которой аппроксимировали данные по пару, приведенные в табл. 4, при этом среднее значение постоянной адиабаты составляло=  1,303. Поскольку в формуле Сен-Венана, определяющей истечение газа из резервуара в атмосферу как функции его температуры и давления, давление пара (см. фиг. 1) значительно больше атмосферного давления, то скорости выхода пара и азота N будут, в основном, определяться их температурой и значением постоянной адиабаты. Поэтому, используя данные табл. 3 для азота и насыщенного пара при T = 20 °C (HO) и T =15 °C (N), находят, что относительная погрешность скорости выброса азота N относительно пара составит не больше 7%, т.е. в пределах погрешности измерения скорости выброса пара. Последнее позволяет принять, это отличие незначительным и считать, что и вынос радиоактивного азота N из паропровода происходит со скоростью выноса пара. Подставляя правую часть уравнения:
м
в уравнение:
N
получают окончательное уравнение для активности азота Q(x, t):
,
0
решение которого при L ≤ x ≤ Ls принимает вид:
,
N
0
N
N
s
s
0
t
гр
s
0
L
L
в котором температура зависит от x, а постоянная Q(L) подлежит определению. При найденной функции активности пара Q(x, t) и плотности пара уточняют ранее введённый параметр α, который можно получить на выходе пара из паропровода в точке x = L (Т= Т + b °С), исключая протечку, поскольку, в момент времени, соответствующий выходу пара из паропровода 5 при P = 6,4 МПа,  = 20 м,  = 0,9063 м, и
При этом скорость пара при его выходе из паропровода 5 при указанных на фиг. 1 значениях температуры и давлении, согласно оценкам, составляет примерно 1512 м/с. При таких скоростях выхода пара на рабочие лопатки турбины 2 развивалась бы очень высокая радиальная скорость турбины 2, что повлекло бы значительный рост центробежной силы, которая бы привела к колоссальному росту напряжений в диске турбины 2 и особенно рабочих лопатках, к превышению допустимых напряжений указанных конструкций и, как следствие, к их разрушению. Эти возможные эффекты и приводят к необходимости снижению давления и температуры пара, путем использования соответствующих технических решений.
N
N
0
Выражение для параметра α[Ки/кг], которое позволяет определить искомое значение активности азота Q(L) в пограничной области вода-пар паропровода 5:
N
16
Величину параметра α определяют по оценке мощности дозы, создаваемой гамма-излучением азота N, из математического которой находят объемную активность.
16
-3
3
3
3
3
2
16
16
3
0
п
0
п
s
0
sf
γ
V
sf
sf
γ
7
V
п
Далее находят мощность дозы γ-излучения азота N при выходе пара в воздушную среду с образованием сферического облака диаметром l = v·τ, где v – скорость распространения пара в воздухе, которая по данным измерений составляет от 6 до 50 м/с; τ – время прохождения пара в паропроводе 5 длиной L=L-L, оценка которого составляет величину τп ~ 13·10 с. С этой целью используют формулу мощности дозы от виртуального шарового источника с равномерно распределенными по объему изотропными точечными источниками с удельной активностью Ки/м без учета многократного рассеяния при условии μr << 1, μ – коэффициент линейного ослабления γ-излучения; r = l/2– радиус сферы; D' = 2πKQ(L), где D' – измеряемая величина мощности дозы (мв/ч) в выделенном сферическом объеме радиоактивного пара V = πl/6; K = 14,652 [мв·м/ч·Ки] – γ-постоянная азота N, Q(L) – виртуальная объемная активность радионуклида N (Ки/м) на выходе из паропровода 5, возникающая за указанный промежуток времени τ, определяемая отношением:
V
sf
γ
Q (L) = D’ / 2πK l.
п
п
п
п
0
п
п
N
3
16
Значение плотности пара ρ(L) = 509,87 кг/м, а его массу на выходе из паропровода 5 в точке x = L определяют как m(L) = V·ρ(L) = πR2L·ρ(L), (m(L) = 10,905 кг. При этом параметр α (Ки/кг), оценку которого получили на основе измерений мощности дозы, создаваемой гамма-излучением азота N, находящимся в облаке пара, принимает вид:
0
N
пср
пср
0
п
п
N
N
N
0
п
-6
2
-6
-6
-6
Если принять, что в результате протечки, радиоактивный пар, выходящий из паропровода, создает мощность дозы сравнимую с радиационным фоном реактора 500 мкР/ч, то при заданных остальных параметрах и v = 50 м/с, величина α составит 1,102·10 Ки/кг и линейно будет возрастать или уменьшаться в соответствии с мощностью реактора. Оценку массы пара, выходящего из паропровода 5, определяют, как среднюю массу m пара со средней плотностью ρ и объемом πR·τ, выходящего из паропровода 5 за время τ со средней скоростью = 1415 м/с. Масса пара составит 6,976 кг, а величина α увеличится до 1,723·10 Ки/кг, т.е. эффективное значение параметра α будет находиться в диапазоне 1,102·10- 1,723·10 Ки/кг, а активность сосредоточенного изотопного источника Q(L), создаваемого радиоактивным азотом в области границы раздела вода-пар, за время τ находят из формулы для определения активности азота в пограничной области вода-пар паропровода 5:
V
0
3
Подставляют известные величины получают объемную активность Q(L) = 0,55 Ки/м.
Для измерения радиоактивности пара на радиационном фоне реактора используют гамма-детекторы, основанные на разностном методе ионизационных токов измерительной аппаратуры, состоящей, например, из проточной и непроточной ионизационных камер, одна из которых будет измерять радиационный фон реактора, а вторая – регистрировать как радиационный фон движущейся радиоактивной среды (радиоактивного пара), так и радиационный фон реактора.
16
7
В случае, рассматриваемом в настоящем изобретении содержание радионуклида радиоактивного азота N (Т½ = 7,11 с, с энергией бета-излучения 6,134 MэВ и квантовым выходом 69%; энергиями бета-излучения 4,288 МэВ, и выходом бета-частиц 68%, 10,419 МэВ и 26%) в струе пара парогенератора 1, поступающего на турбину 2, свидетельствующего о признаке нарушения герметичности водопровода 4 или паропровода 5 второго контура парогенератора 1 определяют путём использования сцинтилляционных детекторов с кристаллом NaJ(Tl), в силу высокой энергии гамма-, бета-излучения, по эффекту образования пара, при выходе пара на турбину 2.
Разность показаний измерительных приборов определит вклад радиационного фона движущейся радиоактивной среды в диапазоне скоростей ее переноса от 0,2 до 10 – 12 м/с.
v
3
Объемную активность азота Q [Ки/м], возникающего в паропроводе 5 в области границы раздела вода-пар, определяют по формуле:
v
N
0
0
vr
2
Q=Q(L)/(πRl),
0
vr
0
0,max
0,min
где R – внутренний радиус водопровода 4 и паропровода 5; l = ΔL = L - L [м] - разность хода измеряемых флуктуаций давления воды.
п
в
Используя найденные значения активности азота в области границы раздела вода-пар за время τ, находят величину мощности вброса радиоактивного азота Р из первого контура во второй:
в
V
0
0
vr
г
-4
P = Q(L)·2πRl·U ≈ 2,0·10 Ки/с.
N
в
v
3
Величину секундного вброса G [м/с] радиоактивного азота в область границы раздела вода-пар находят, используя выражение, определяемое отношением мощности вброса Р к объемной активности Q:
N
в
V
0
-4
3
G = P/Q(L) = 3,64·10 м/с.
Скорость вброса:
г
N
в
0
пср
U=G/S = Rv/2L=0,116 м/с.
Сравнивают полученные результаты измерения мощности дозы и предельной концентрацией радионуклида в помещении на выходе паропровода 5 на турбину 2 с нормативными данными, обеспечивающие безопасную работу персонала, принимают соответствующие решения: продолжить работу, поскольку выход радионуклида не приводит к превышению норм по радиационной безопасности; провести профилактические работы – провести дегазацию помещения; провести временный останов оборудования и провести ремонтные работы.