Аннотация

Аннотация

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

Полный текст

Введение В технике и различных технологиях современного производства самое широкое применение находят изделия, обладающие эффектом памяти формы [1]. В частности, электротехническая, кабельная и строительная индустрии являются наиболее крупными потребителями этой продукции [2-4]. Для создания молекулярно-сшитых полиэтиленов, наряду с чисто химическими технологиями и процедурой облучения направленными потоками ускоренных электронов, используется процедура гамма-облучения материала [3]. Технология создания дополнительных межмолекулярных связей с помощью ионизирующего гамма-облучения от источника Сo60 (кобальт 60) имеет не только ряд определенных преимуществ, связанных с качеством получаемого продукта, но и определенные конструктивные особенности. Радиационная доза, необходимая для возникновения эффекта сшивки, должна иметь вполне определенное значение, которое должно достигаться по возможности одновременно во всех точках облучаемого материала. В реальном производстве, в целях гарантированного исключения брака (недостаточность нерастворимой гель-фракции в материале), время экспозиции рассчитывается по минимальному значению мощности дозы. Она, в свою очередь, зависит, в том числе, от расположения источников ионизирующего излучения (ИИИ). Оптимальное размещение точечных ИИИ является предметом настоящих исследований. Следует сразу отметить, что гамма-излучение является изотропным. В существующих промышленных установках (ГУД 300, ГУД 300М) расположение источников ионизирующего излучения с незначительной погрешностью можно принять за прямолинейное. Облучаемые материалы находятся во вращающихся герметичных контейнерах, при этом ось вращения параллельна прямой линии, вдоль которой расположены источники гамма-квантов. Места расположения источников излучения фиксированы в пространстве и не могут быть изменены по конструктивным соображениям. В каждой точке локации источников излучения могут быть размещены несколько радиоактивных элементов различной мощности. Порядок расположения радиоактивных источников является параметром управления в задаче оптимизации. Постановка задачи и эвристический алгоритм решения Пусть имеется N источников гамма-излучения и М точек их размещений (N >> M). Обозначим через Аi активность i-го источника излучения. На момент t = 0 она известна для каждого i и обозначена как А0i. Имея в виду упомянутый изотропный характер гамма-излучения и его зависимость от времени (периода полураспада), мощность экспозиционной дозы Р(х,t) в точке с координатой х будет определятся как (1) где Гс - ионизационная гамма-постоянная; Т1/2 - период полураспада, Т1/2 = Тn, который для Co60 составляет 5,2 года; ri - расстояние от i-го источника до точки с координатой x, i - индекс суммирования по всем источникам от 1 до N. В каждый момент времени t реализуется некоторое распределение мощности облучения, для которого можно найти минимальное значение мощности: (2) Эта величина является базовой для расчета времени экспозиции облучаемого материала. Вследствие неравномерности распределения источников излучения некоторая часть материала получает излишнюю дозу, которую условно можно назвать потерянной и которая показана в виде заштрихованной области W на рис. 1. Рис. 1. Схема процесса облучения полимерного материала (1) с помощью источников излучения (2) Тогда очевидно, что оптимальное размещение N источников ионизирующего излучения должно обеспечивать минимум значению W - потерянной мощности: (3) где Т - период эксплуатации радиационной установки между двумя последовательными процедурами перезарядки; L - линейный размер зоны облучения. Отметим, что значение функционала W может быть пронормировано по L и T. При рассмотрении частного случая, когда все источники ИИ в формуле (1) распределены таким образом, что обеспечивается равномерность радиационной активности по длине L, для некоторого фиксированного момента времени t несложно получить аналитическое решение распределения мощности вида: p(x) = c/h (arctg ((L - x) / h - - (- arctg x / h)), где h - расстояние от линии размещения источников до линии размещения материала, а с - некоторая константа. Соответственно, pmin = c / h arctg (L / h) и, разумеется, p(x) - симметрична относительно L/2. Наличие подобных решений наводит на мысль о возможности сведения задачи оптимизации к поиску параметров некоторой непрерывной функции распределения активности А(х, t), заданной, например, в форме полинома и обеспечивающей выполнение условия (3) при существовании интегрального ограничения на А(х, t), суть которого сводится к обеспечению условия для суммарной активности источников ИИ. Следует отметить, что даже при наличии эффективного алгоритма по построению такой функции не удается избежать решения целочисленной задачи распределения активности дискретных точечных источников ИИ в соответствии с найденным решением. С учетом этого замечания минимизация функционала (3) проводилась путем решения задачи целочисленной оптимизации при наличии ряда технических ограничений: 1. Допускается размещение в одной точке локации с номером j нескольких источников гамма-излучения разной активности. Также допустимо и отсутствие источников в некоторых местах локации. 2. Полный набор источников излучения, предназначенных для зарядки в одну камеру для облучения, может содержать несколько меньших наборов источников с одинаковой активностью. Например, полный набор из N источников может содержать N1 источников с активностью А1, N2 источников с активностью А2 и т.д. При этом N1 + N2 + … + Nk = N. Как правило, число типов таких наборов k невелико, не более 3-4. 3. В силу технологических особенностей задача оптимального размещения источников имеет симметрию относительно точки х = L/2. Эта особенность существенно снижает вычислительные затраты на решение задачи. В общем случае, при наличии N источников ИИ и М точек их размещения (локаций) на отрезке длиной L, число возможных вариантов размещения чрезвычайно велико - M(N), и это обстоятельство следует учитывать при создании алгоритма оптимизации, который может иметь весьма серьезную вычислительную трудоемкость. Как правило, в одной точке локации физически может быть размещено не более 3 источников ИИ одновременно. Таким образом, количество возможных комбинаций в одной точке локации существенно сокращается. Кроме того, необходимо исключить из вычислительных процедур варианты расчетов W по формуле (3) с числом размещенных ИИИ, отличным от N. Оптимальное решение ищется в подпространстве N1 + N2 + … … + Nk - N = 0. Иные варианты, когда размещены не все ИИИ, очевидно, не имеет смысла рассматривать. Таким образом, задача дискретной оптимизации может быть сформулирована следующим образом: найти такое оптимальное распределение источников ионизирующего излучения в местах локации при которых целевая функция на заданном интервале времени от 0 до Т достигает минимума: (4) при следующих ограничениях: - в одной трубке может быть не более 3 ИИИ, количество трубок в камере (мест локации ИИИ) равно M; - общее количество источников различного типа равно N; - распределение должно быть симметричным относительно середины камеры Следует отметить, что задача (4) относится к задачам комбинаторной оптимизации [5] и является np-сложной. Для решения подобных задач обычно применяются эвристические алгоритмы [5], к которым можно отнести генетические [6], «муравьиные» [7] алгоритмы, метод Монте-Карло и т.п. Все эвристические алгоритмы направлены на уменьшение допустимого множества вариантов решения задачи комбинаторной оптимизации и ускорение сортировки оставшихся вариантов. В данном случае применим эвристику, основанную на аналитическом решении задачи при допущении о непрерывном распределении источников облучения, которое приведено выше. Как было показано, при равномерном распределении активности источников по длине A(x) функция мощности облучения P(x) имеет максимум по середине (графики 1 на рис. 2) и потери мощности достаточно велики. Чтобы добиться минимума потерь мощности, необходимо располагать источники таким образом, чтобы их активность на концах отрезка [0, L] возрастала, а функция мощности стремилась к константе (графики 2 на рис. 2). Рис. 2. Качественная картина распределения мощности излучения P(x) (пунктирная линия) при непрерывном распределении активности источников A(x) (сплошная линия) Теперь рассмотрим эвристический алгоритм приближенного решения поставленной дискретной задачи оптимизации (4), построенный на основе аналитического решения непрерывной задачи. Первым шагом в работе алгоритма является поиск комбинаторным способом множества допустимых вариантов распределения источников в местах локации с учетом ограничений задачи оптимизации. Как уже отмечалось, количество допустимых вариантов очень велико. Поэтому требуется разумно ограничить это множество таким образом, чтобы поиск оптимального распределения не занимал слишком много времени. Именно поэтому вторым шагом алгоритма является учет следующих эвристических условий: - распределение источников должно быть таким, чтобы их активность постепенно снижалась от края к середине отрезка [0, L], что соответствует функции 2 на рис. 2; - распределение источников должно быть таким, чтобы аппроксимирующая функция их активности была достаточно гладкой на заданном отрезке, из чего следует, что из допустимого множества решений можно исключить варианты, в которых суммарные активности источников в двух последовательных точках локации были бы максимальны или равнялись нулю. Последним шагом алгоритма является сортировка оставшихся вариантов распределения источников по значению целевой функции и поиск таких вариантов, в которых она минимальна. Демонстрационный пример В качестве примера решения рассмотрим вариант размещения в камере облучения 48 источников ИИ, 24 из которых имеют на момент зарядки в камеру активность А = 4,54·1013 Бк каждый (на схеме обозначены цифрой 1), а у 24 оставшихся активность каждого составляет А = 7,23·1013 Бк (на схеме обозначены цифрой 2). В любой точке локации число размещенных источников не должно превышать трех. Требуется разместить все источники в 34 местах локации. Схема типового размещения ИИИ представлена в виде таблицы (цифра 0 обозначает отсутствие ИИИ): 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Распределение потерь мощности при таком варианте размещения приведено на рис. 3. Соответствующее распределение активности источников приведено на рис. 5. Показанное на рис. 3 распределение потерянной мощности содержит характерное «западание мощности» по концам зоны облучения. Такой характер распределения весьма типичен для реального производства. Оптимальное, в смысле минимума потерь мощности, размещение 48 источников, полученное в результате расчетов, выглядит следующим образом: 2 2 2 1 2 2 2 1 1 1 1 0 2 2 1 1 2 2 1 1 2 2 0 1 1 1 1 2 2 2 1 2 2 2 1 2 2 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 2 1 2 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 2 Рис. 3. Мгновенные потери мощности (при t = 0 синий график, при t = 3 года красный график) Распределение потерь мощности по длине, соответствующее этому варианту, приведено на рис. 4. На рис. 6 приведено соответствующее этому варианту распределение активности источников излучения. Необходимо отметить, что значения параметров потерь мощности в двух сравниваемых вариантах фактически различаются в разы, кроме того, в оптимальном варианте, как и ожидалось, пропал упомянутый выше эффект «западания мощности» по концам зоны облучения. Рис. 4. Мгновенные потери мощности (при t = 0 синий график, при t = 3 года красный график) Рис. 5. Распределение активности при типовом размещении источников. Верхний график соответствует активности в момент t = 0, нижний t = 3 года Рис. 6. Распределение активности при оптимальном размещении источников. Верхний график соответствует активности в момент t = 0, нижний t = 3 года В целях оценки эффективности полученного решения следует отметить, что в момент времени t = 0 продолжительность экспозиции рассчитывается для типового варианта, исходя из мощности Pmin = = 0,22 Грей/с, а для оптимального варианта распределения Pmin = = 0,27 Грей/с. В момент времени t = 3 года показатели Pmin = 0,14 Грей/с и Pmin = 0,18 Грей/с соответственно. Таким образом, общая эффективность использования предложенного варианта размещения источников находится в диапазоне 22-28 %, что, несомненно, является экономически весьма привлекательным, поскольку не связано с какими-либо дополнительными затратами в оборудовании или технологиях. Из рис. 6 видно, что оптимальное распределение активности источников по длине камеры близко к тому, которое задано эвристикой, полученной из аналитического решения задачи оптимизации.