Подготовка математических программ Deca и Radrisk к проведению экспертизы безопасности проекта Белорусской АЭС

Размещено на http://www.allbest.ru/

Подготовка математических программ Deca и Radrisk к проведению экспертизы безопасности проекта Белорусской АЭС

Н.В. Горбачева,

Н.В. Кулич,

Н.В. Береснева

ГНУ «ОИЭЯИ-СОСНЫ» НАН БЕЛАРУСИ, Минск, Беларусь

После получения эксплуатирующей организацией лицензии на размещение атомной станции на площадке в Гродненской области актуальной задачей для регулирующего органа Республики Беларусь является экспертиза безопасности проекта энергоблока АЭС-2006 с реактором ВВЭР-1200. Предполагается, что защитная оболочка реактора в проекте белорусской АЭС будет иметь инновационное инженерное решение, что позволит смягчить последствия аварийного выброса при авариях с повреждением активной зоны.

ГНУ «ОИЭЯИ-Сосны» НАН Беларуси оказывает научно-техническую поддержку деятельности регулирующего органа при проведении экспертизы безопасности белорусской АЭС. Ведется подгототовка специалистов, освоение тяжелоаварийных кодов анализа безопасности блока АЭС. Планируется получить результаты к моменту ввода в эксплуатацию первого энергоблока. Для решения текущих задач экспертизы материалов, обосновывающих безопасность проекта энергоблока, в институте ведется работа по развитию собственных программных средств.

Требования нормативного регулирования высшего уровня [1, 2] предписывают обязательность выполнения вероятностного анализа безопасности при проектировании блока белорусской АЭС и устанавливает целевые ориентиры для вероятностных показателей безопасности в отношении возможного ущерба здоровью населения. Так, в целях исключения необходимости эвакуации населения за пределами зоны планирования защитных мероприятий, следует стремиться к тому, чтобы оцененное значение вероятности предельного аварийного выброса не превышало 10-7 на реактор в год. Величина выброса основных дозообразующих нуклидов в окружающую среду при запроектных авариях должна быть ограничена в такой степени, чтобы дозы облучения населения на границе зоны планирования защитных мероприятий и за ее пределами не превышали значений, требующих принятия решений о мерах защиты населения. безопасность аварийный выброс

Введенные с 01.01.2013 г. Санитарные нормы и правила «Требования к радиационной безопасности» приведены в соответствие с рекомендациями МАГАТЭ, в которых требуется учитывать дисперсный состав выброса при оценке дозы облучения персонала, а оценку воздействия на население необходимо проводить для различных возрастных групп [3, 4].

В требованиях к содержанию ООБ в части анализа аварийных процессов [5] указывается, что при выполнении оценки выхода радионуклидов при авариях необходимо учитывать глубину выгорания топлива, изменение во времени мощности энергоблока, активность по изотопному составу топлива, относительное содержание изотопов йода с учетом органической и элементарной фракций.

Обоснование безопасности новых проектов АЭС объективно сопряжено с рядом трудностей. МАГАТЭ и регулирующий орган США обращают внимание, что для нового блока АЭС характеристики выброса радионуклидов в атмосферу в ходе развития тяжелой аварии могут существенно отличаться от базового проекта. Поэтому, при обосновании безопасности необходимо более тщательно подойти к вероятностному анализу безопасности уровня 3, хотя в настоящее время и для действующих АЭС не всегда результаты ВАБ-2 заканчиваются выполнением ВАБ-3 [6, 7].

В данном отчете предложен инженерный, практический подход к прогнозированию характеристик источника выброса при анализе безопасности проекта блока в поддержку деятельности регулирующего органа. Так, для того, чтобы иметь возможность уточнять последствия выбросов (например, образование органических соединений йода), необходимо разработать параметрическую модель аварийного источника. При разработке модели все радионуклиды, которые могут оказать негативное воздействие в результате выброса, группируются по сходству физико-химических свойств в стандартные группы, как это принято при выполнении ВАБ-2. В отношении наименее изученных процессов, влияющих на формирование аварийного источника, неопределенность которого в результате может оказаться весьма существенной, следует привлекать экспертные оценки. Так, в данном отчете экспертным путем получены оценки доли от активности продуктов деления, накопленной в активной зоне, по группам с различными летучими свойствами.

Для более эффективного обращения с неопределенностями, имеющими место при обосновании безопасности при лицензировании блока АЭС, хорошо зарекомендовал себя подход к разработке и использованию результатов кодов, основанный на улучшенной оценке (“best-estimated”) [8]. В своем развитии данный подход имеет несколько приближений к модели, входным данными и граничным условиям, а также предположения, касающегося выполнения функции критическими системами безопасности. Для целей лицензирования предпочтение отдается кодам улучшенной оценки, основные требования к которым следующие. Для кодов, выполняющих улучшенную оценку аварийных процессов, необходимо, чтобы модель достаточно подробно отражала ведущую закономерность, а при его разработке использовались реалистичные исходные данные и граничные условия, где адекватным для определения «реалистичные данные» является «статистические данные». При использовании результатов расчетов кодами улучшенной оценки в целях лицензирования ключевым элементом является анализ неопределенности.

В настоящее время для анализа сложных аварийных процессов в ядерно-энергетических установках зарекомендовала себя практика разработки и применения расчетных комплексов, выполняющих сопряженный расчет с привлечением кодов улучшенной оценки. В данной работе представлен опыт разработки расчетного комплекса программ на базе собственных, разработанных в ОИЭЯИ-Сосны, кодов DECA и RADRISK для обоснования безопасности проекта блока АЭС-2006 в части оценки радиологических последствий запроектных аварий. Информация, представленная ниже, показывает, что данные коды соответствуют понятию “коды улучшенной оценки”, в которых в качестве неопределенных входных использованы реалистичные приближения.

Базовая версия кода DECA выполняет расчет изотопного состава топлива, концентраций и удельных активностей актиноидов и продуктов деления в активной зоне ВВЭР-1000 в зависимости от глубины выгорания. [9]. Изотопный состав продуктов деления насчитывает 650 ядер, образующихся в изобарных цепочках с атомными весами от 72 до 166. Изотопный состав актиноидов насчитывает 58 ядер от изотопов Th до Fm. Математическая модель описывает выход продуктов деления в процессе деления тепловыми нейтронами первичного топлива - 235U и вторичного топлива - 239Pu и 241Pu, ядерные и нейтронно-физические превращения актиноидов. При разработке кода DECA решена проблема сложности моделируемого процесса, связанная со структурой модели. Для описания структуры генетических связей трансурановых элементов и продуктов деления в модели использован математической аппарат теории графов, что позволяет учесть реакции, приводящие к обратным превращениям от дочернего ядра к материнскому в схеме трансмутации актиноидов, например, реакцию 238Pu(n,)239Pu и другие из пар (n,2n), (n,).

Сопоставление результатов расчетов кодом DECA с расчетами кодом ORIGEN-S показало, что при различии в глубине выгорания 2,7%, отличие данных следующее: по активностям продуктов деления Kr-85 - 3%, Zr-95 - 3%, Сs-137 - 5%, Cs-134 - 15%, Sr-90 - 10%, Ce-144 - 9%; по актиноидам: Pu-238 - 2%, Pu -239 - 1%, Np- 237 - 17%, Am-241 - 20%, Am-243 - 4.4%, Cm-242 - 1,5%, Cm-244 - 43% [9]. Проведенный анализ качества модели выгорания показал, что модель адекватно отражает закономерности формирования радионуклидного состава в процессе выгорания топлива, а численный метод эффективно и с высокой точностью выполняет расчет концентраций и активностей радионуклидов. Неопределенность результатов расчетов радиационных характеристик топлива определяется погрешностями используемых в модели выгорания кода DECA ядерно-физических констант, взятых из библиотеки JNDC-V2 (JAERI-1320, 1990), неопределенностями нейтронно-физических расчетов, неоднородностью энерговыделения по объему активной зоны в точечном приближении. Особо следует отметить, что разработанная модель обеспечивает надежную взаимосвязь рассчитываемых характеристик с глубиной выгорания топлива, а также соотношение активностей радионуклидов внутри каждой фракции химических элементов.

В процессе разработки расчетного комплекса выполнена модернизация кода DECA по следующим направлениям. Разработана версия кода DECA, выполняющего вероятностный анализ неопределенности результатов расчета концентраций и активностей при неопределенных входных параметрах модели выгорания. Предварительно проанализированы и получены оцененные интервалы изменения 40 стохастических параметров модели выгорания кода DECA. Использована рекомендуемая МАГАТЭ вычислительная технология, выполняющая вероятностно-статистический анализ неопределенности с применением процедур Монте-Карло [10]. Вычислительная процедура, оформленная в виде программного модуля CUB, подключена к коду DECA. В качестве меры неопределенности результатов моделирования кодом DECA концентраций Ci и активностей Ai радионуклидов, образующихся в процессе выгорания топлива в реакторе, принят 90%-й доверительный интервал [Сi,q=5%, Сiq=95%], где Сi,q=5% , Ci,q=95%, - значения 5% и 95%-ых квантилей концентрации i- го радионуклида.

Как говорилось выше, для обоснования безопасности на этапе лицензирования блока АЭС целесообразно выполнять ВАБ-3. Для этих целей в ГНУ «ОИЭЯИ-Сосны» разработан программный продукт - код RADRISK. Код RADRISK выполняет расчет радиологических последствий аварийных выбросов радионуклидов для населения и окружающей среды и дает оценку радиационного риска. Для расчета вероятности воздействия на население аварийных выбросов в код RADRISK включен программный модуль CUB, который создает выборки случайных значений входящих в вероятностную модель параметров, таких как статистические характеристики метеорологических параметров, активность радионуклидов, высоты подъема облака и других неопределенных параметров. Результаты расчета показателей ущерба здоровью и уровня загрязнения окружающей среды представлены в виде 5-, 50-, 95%-ых квантилей плотности поверхностного загрязнения, приземной концентрации радионуклидов, индивидуальной дозы облучения населения и радиационного риска, что соответствует методики выполнения ВАБ-3 [11, 12].

Известно, что в ходе развития тяжелой аварии образуется огромное количество аэрозолей, в том числе горячих частиц. Радиоактивные аэрозоли, образующиеся при эксплуатации ядерных реакторов, отличаются в зависимости от типа реактора сложностью дисперсного и многокомпонентностью радиоизотопного состава [13]. За исключением РБГ и газообразных форм йода (I2, CH3I), все радионуклиды представлены в выбросе в форме частиц аэрозоля. Поэтому полное описание аварийного источника должно включать физические характеристики аэрозольного выброса, такие как дисперсный состав, плотность аэрозоля, концентрацию частиц. Однако, эта информация редко бывает доступна из результатов ВАБ-2. Поэтому при выполнении расчетов радиологических последствий в большинстве кодов, выполняющих ВАБ-3, принято считать, что все частицы аэрозольного выброса имеют аэродинамический диаметр (АМАД) 1 мкм. Однако необходимость учета реальных аэродисперсных характеристик радиоактивного аэрозоля заключается в том, что при одной и той же концентрации и времени экспозиции его аккумуляция в различных отделах органов дыхания значительно зависит от аэродинамического диаметра.

В данной работе для того, чтобы получить более реалистичные оценки радиологических последствий, разработана и введена в код RADRISK дискретная вероятностная модель источника выброса.

Для описания дисперсного состава наиболее оправдано теоретически установленное А.Н. Колмогоровым логарифмически-нормальное распределение для функции плотности распределения F(X) числа частиц по размерам диаметров, X:

, (1)

где м - медиана и - стандартное отклонение величины ln X.

Предпочтение, отдаваемое логарифмически-нормальному распределению связано с удобством решения практических задач, а также с тем фактом, что оно связано с физической сущностью образования аэрозольных систем в результате специфических процессов диспергирования [13]. В работе [14] авторами показано, что статистические характеристики размеров топливного зерна могут служить замещающим критерием (признаком) для выражения вероятностной меры связи между результатами ВАБ-1 и ВАБ-3.

Значения µ и могут быть получены из результатов статистической обработки выборки измерений размеров частиц по следующим формулам:

, (2)

где Х10 и Х90 - 10- и 90%- ные квантили, соответственно.

Параметры логарифмически-нормального распределения ХG и ?G могут быть использованы для характеристики средней поверхности и среднего объема частиц с тем же стандартным отклонением ?G на основании соотношения

(3)

при v=2 и v=3.

В работе [13] показано, что не только распределение частиц по размерам, но и связанные с ними физические характеристики подчиняются логнормальному закону. Так, при н=2 активность аккумулируется на поверхности частицы, а при н=3 активность содержится в объеме частицы:

(4)

Для реализации дискретной модели аварийного источника в виде совокупности радиоактивных частиц, статистические размеры которых распределены по логнормальному закону (1), использованы процедуры Монте-Карло. Пусть Хi - случайный размер частицы из интервала [Xmin, Xmax]. Исходя из предположения о виде функций распределения, с помощью метода Монте-Карло построим выборку размера n:

X_min Х₁, Х₂, Х₃, …, Х_n X_max (5)

Из выборки (5) можно найти все характеристики, включая среднее значение, моду, медиану, стандартное отклонение на интервале [Xmin, Xmax]. Построив на основании (5) интегральную функцию распределения, можем найти численные оценки для квантилей распределения любого порядка.

Например, если найдено, что 95%-ая квантиль не превосходит Хi, то факт, что лишь 5% частиц могут иметь размер, превосходящий Хi, может быть полезно использован при анализе неопределенности внесенной в результат предположением о виде распределения.

При наличии в выбросе нескольких типов аэрозолей дискретная вероятностная модель аварийного источника представляется суммой по всем n-частицам каждой из фракций (Qk ):

(6)

Размеры частиц, плотность аэрозоля и концентрация частиц влияют на поведение радиоактивных материалов в процессе переноса и осаждения при атмосферной диффузии. Процессы, описывающие поведение аэродисперсных частиц в процессе атмосферной диффузии и оседания слишком сложны, чтобы использовать детальные модели в практических целях. В данной работе использованы эмпирические данные для скоростей оседания частиц в зависимости от диаметра, полученные в работе [15] (таблица 1).

Таблица 1 Скорости оседания частиц в зависимости от диаметра

Прогноз воздействия выброса на население при запроектной тяжелой аварии на АЭС-92

Для проведения расчетных исследований воздействия выброса радионуклидов на население при запроектной аварии, результаты расчета наработки радионуклидов в активной зоне кодом DECA сопряжены с входными данными кода RADRISK посредством файла обмена данными.

Код DECA выполняет расчет активностей и масс продуктов деления, накапливающихся в процессе выгорания топлива, с последующей сортировкой по номенклатуре радионуклидов с различными физико-химическими свойствами. Результаты расчетов наработки масс и активностей в летучих фракциях радионуклидов в топливе реактора ВВЭР-1000 при глубине выгорания 11 МВт*сут/кг на момент аварии приведены в таблице 2.

Таблица 2 Активности и массы радионуклидов (50% квантиль), содержащиеся в различных фракциях радионуклидов, полученные кодом DECA

Ранний отказ защитной оболочки при тяжелой аварии

В качестве параметра, характеризующего величину выброса активности Qi i-го радионуклида из разрушенной активной зоны при различных авариях, используется коэффициент фракционирования fi. Прогноз величины выброса радионуклидов при j-классе аварии определяется через параметры источника выброса по следующей формуле:

; (7)

, (8)

где Ai - активность i-го радионуклида в активной зоне к моменту аварийного сброса мощности реактора, I 1,…I; qji - доля активности i-го радионуклида, высвобождающаяся из топливной матрицы в зависимости от температуры в активной зоне при j-аварии; ej - степень разрушения активной зоны при j-аварии; Lj - протечка, доля активности выходящая в окружающую среду при разгерметизации системы защитной оболочки реактора, fi. - коэффициент фракционирования радионуклидов при тяжелых авариях. Коэффициенты фракционирования радионуклидов при запроектной аварии на реакторах типа ВВЭР-1000, относящихся к группам химических элементов с различными физико-химическими свойствами, получены по материалам работы [6].

Оцененный выброс изотопов РБГ (85mKr, 85Kr, 87-88Kr, 133Xe, 135Xe), цезия (134Cs, 137Cs) и йода (131-135I) при отказе защитной оболочки в первые часы после начала развития аварийного процесса составил величину Q = 1.685•1016Бк.

Дисперсный состав и вклад в активность различных фракций выброса радионуклидов, полученные с использованием дискретной вероятностной модели аварийного источника, приведены в таблице 3.

Таблица 3 Дисперсный состав и распределение активности аварийного выброса

На рисунках 1, 2 приведены результаты расчета модернизированным кодом RADRISK, соответственно, плотности поверхностного загрязнения, С(x,y) и эффективной индивидуальной дозы облучения персонала, E(x,y), для следующих характеристик выброса: а) линия ____ ? частицы со средним диаметром 1 мкм, б) линия_ _ _? частицы со средним диаметром 5 мкм; с) линия ……. ? без учета дисперсных свойств выброса (скорости осаждения изотопов йода - 2 см/с, цезия - 0,8 см/с). Результаты получены при предельно неустойчивой категории погоды (А), а результаты оценки дозы персонала в соответствии с требованиями Гигиенического норматива «Критерии оценки радиационного воздействия» Республики Беларусь.

Рисунок 1 - 95% квантиль плотности поверхностного загрязнения (в логарифмическом масштабе)

Рисунок 2 - 95 % квантиль дозы облучения персонала.

Результаты расчетов с помощью модернизированного кода RADRISK показывают, что в случае выброса аэрозолей со средним диаметром частиц 5 мкм оценка доз облучения персонала примерно на порядок выше, чем при выбросе аэрозолей со средним диаметром частиц 1 мкм, но ниже чем консервативный расчет без учета дисперсных свойств выброса.

Более реалистичный прогноз воздействия выброса со сложным аэрозольным составом может быть получен и для населения, применяя рекомендованные МКРЗ дозовые коэффициенты.

Для выполнения независимой проверки результатов обоснования безопасности блока Белорусской АЭС с учетом отечественных нормативных требований выполнена модернизация собственных программных средств DECA и RADRISK и последующее сопряжение их в программный комплекс. Для модернизации кодов использован перспективный для анализа безопасности при лицензировании проектов АЭС методический подход, обеспечивающий «улучшенную оценку» по сравнению с консервативным подходом результатов моделирования и анализ неопределенности прогноза воздействия на население аварийных выбросов. Разработана и внедрена в код RADRISK дискретная вероятностная модель, учитывающая сложный дисперсный состав выброса, применение которой обеспечивает реалистичные оценки при наличии полидисперсного состава выброса.

Сопряжение кода RADRISK с кодом DECA позволяет оперативно учитывать реальное выгорание топлива и возможные изменения изотопного состава, вызванного переходными и аварийными процессами, и другие требования регулирующего органа при выполнении независимой. Работоспособность комплекса программ DECA и RADRISK продемонстрирована на примере запроектной аварии ВВЭР-1000.

Список литературы

1 Общие положения обеспечения безопасности атомных станций (ОПБ АС). ТКП 170-2009 (02300)

2 Правила ядерной безопасности реакторных установок атомных станций (ПБЯ РУ АС) ТКП 171-2009 (02300)

3 Санитарные нормы и правила «Требования к радиационной безопасности». Гигиенический норматив «Критерии оценки радиационного воздействия» Пост. Министерства здравоохранения Республики Беларусь № 213 от 28.12.12.

4 Международные основные нормы для защиты от ионизирующих иэлучений и опасного обращения с источниками излучения. Серия изданий по безопасности № 115. МАГАТЭ, Вена, 1997 г.

5 Требования к содержанию отчета по обоснованию безопасности атомной электростанции с реактором типа ВВЭР. ТКП 294-2010 (02300)

6 Procedures for conducting probabilistic safety assessment of nuclear power plants (level 2) // Safety series № 50-P-8. - Vienna: IAEA, 1995. - 132 p.

7 L. Soffer, S. B. Burson/ Accident Source Terms for Light-Water Nuclear Power Plants. Final Report U.S. Nuclear Regulatory Commission. Office of Nuclear Regulatory Research NUREG1465.

8 Р. Уголева Экспертиза и аттестация ПС, применяемых для расчетов и обоснования безопасности в области использования атомной энергии. Труды НТЦ ЯРБ. - М.: НТЦ ЯРБ. 2011. - 65.с.

9 Н.В. Горбачева, Н.В. Береснева, Н.В. Кулич, В.В. Скурат. Верификация и использование математических программ для оценки запасов радиоактивности в обеспечение радиационной безопасности комплекса ядерных установок на площадке АЭС/6-я МНТК «Обеспечение безопасности АЭС с ВВЭР» , 26 - 29 мая 2009 г. Подольск, ФГУП ОКБ «ГИДРОПРЕСС».

10 Evaluating the Reliability of prediction made using environmental transfer models // Safety Series №. 100. - Vienna: IAEA, 1989. - 106 р.

11 N. Harbachova, N. Kulich. Model Advancement In Level 3 PSA Application with the Use of Code RADRISK /Output of the IAEA Technical Meeting on Level 3 Probabilistic Safety Assessment. Working Materials. IAEA Headquarters, Vienna, Austria July 2-6, 2012.

12 Н. В. Горбачева, Н.В. Кулич. Разработка программных средств для проведения работ по независимой проверке результатов анализа безопасности белорусской АЭC в части ВАБ-3/ 7-я МНТК «Обеспечение безопасности АЭС с ВВЭР», 17-20 мая 2011 г. Подольск, ФГУП ОКБ «ГИДРОПРЕСС».

13 Быховский, А.В. Горячие аэрозольные частицы при техническом использовании атомной энергии. М.: Атомиздат, 1974. - 256 с.

14 N. Harbachova, N. Kulich. The code RADRISK Enhancement for NPP safety justification: fuel design changes/ Proc. of Technical Meeting on Probabilistic Safety Assessment for New Nuclear Power Plants' Design. IAEA, Vienna, Austria, October 1-5, 2012.

15 Пути миграции искусственных радионуклидов в окружаеющей среде. Радиоэкология после Чернобыля. /Под. ред. Ф. Уоррена и Р. Харрисона. М.: Мир,1999. - 512 с.

Размещено на Allbest.ru