Научная статья на тему 'Полнопрофильный рентгенофазовый анализ состава электролита алюминиевого производства на основе мультипопуляционного эволюционного алгоритма и данных элементного анализа'

Полнопрофильный рентгенофазовый анализ состава электролита алюминиевого производства на основе мультипопуляционного эволюционного алгоритма и данных элементного анализа Текст научной статьи по специальности «Математика»

CC BY
183
22
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГЕНЕТИЧЕСКИЕ АЛГОРИТМЫ / GENETIC ALGORITHMS / САМОКОНФИГУРАЦИЯ / SELF-CONFIGURATION / РЕНТГЕНОВСКАЯ ПОРОШКОВАЯ ДИФРАКЦИЯ / X-RAY POWDER DIFFRACTION / FULL-PROFILE ANALYSIS / МЕТОД РИТВЕЛЬДА / RIETVELD METHOD / КОЛИЧЕСТВЕННЫЙ РЕНТГЕНОФАЗОВЫЙ АНАЛИЗ / QUANTITATIVE X-RAY PHASE ANALYSIS / ЭЛЕКТРОЛИТ АЛЮМИНИЕВОГО ПРОИЗВОДСТВА / ALUMINUM ELECTROLYTE

Аннотация научной статьи по математике, автор научной работы — Дубинин П.С., Залога А.Н., Бураков С.В., Гусев К.А., Андрющенко Е.С.

Разработан метод количественного рентгенофазового анализа (КРФА) на основе самоконфигурируемого мультипопуляционного параллельного генетического алгоритма (СГА), который успешно применен для бесстандартного автоматического полнопрофильного анализа образцов электролита алюминиевых электролизеров. Метод позволяет избежать работы по эмпирическому выбору настроек генетического алгоритма и показывает достаточную точность автоматического КРФА по методу Ритвельда под управлением СГА (полученные значения криолитового отношения соответствуют аттестованным значениям с точностью, характеризуемой стандартным отклонением всего 0,035 ед. КО), а также возможность его использования для оперативного контроля состава электролита на алюминиевом производстве.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Дубинин П.С., Залога А.Н., Бураков С.В., Гусев К.А., Андрющенко Е.С.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Full-Profile X-Ray Phase Analysis of the Composition of Aluminum Electrolyte on the Basis of a Multipopulation Evolutionary Algorithm and Elemental Analysis Data

We developed a method of self-configurable multipopulational parallel genetic algorithm (SGA) and successfully applied it to a standardless automatic full-profile X-ray phase analysis of aluminum bath samples. The method spares a user-analyst from working on the empirical choice of the genetic algorithm settings and shows the efficiency of the automatic QPA by the Rietveld method under the control of the SGA, as well as the possibility of using it for express monitoring of the bath composition in aluminum production.

Текст научной работы на тему «Полнопрофильный рентгенофазовый анализ состава электролита алюминиевого производства на основе мультипопуляционного эволюционного алгоритма и данных элементного анализа»

Journal of Siberian Federal University. Chemistry 2 (2018 11) 262-272

УДК 539.26; 519.65; 51-72; 543.44

Full-Profile X-Ray Phase Analysis of the Composition of Aluminum Electrolyte on the Basis of a Multipopulation Evolutionary Algorithm and Elemental Analysis Data

Petr S. Dubinin*3, Alexandr N. Zalogaa, Sergey V. Burakovb, Konstantin A. Gusevb, Evgeniy S. Andryushchenkoa, Igor S. Yakimova, Oksana E. Bezrukovaa and Alexandr S. Samoiloa

aSiberian Federal University 79 Svobodny, Krasnoyarsk, 660041 Russia bReshetnev Siberian State University of Science and Technology 31 Krasnoyarsky Rabochy, Krasnoyarsk, 660037, Russia

Received 28.03.2018, received in revised form 03.04.2018, accepted 12.05.2018

We developed a method of self-configurable multipopulational parallel genetic algorithm (SGA) and successfully applied it to a standardless automatic full-profile X-ray phase analysis of aluminum bath samples. The method spares a user-analyst from working on the empirical choice of the genetic algorithm settings and shows the efficiency of the automatic QPA by the Rietveld method under the control of the SGA, as well as the possibility of using it for express monitoring of the bath composition in aluminum production.

Keywords: genetic algorithms, self-configuration, X-ray powder diffraction, full-profile analysis, Rietveld method, quantitative X-ray phase analysis, aluminum electrolyte.

Citation: Dubinin P.S., Zaloga A.N., Burakov S.V., Gusev K.A., Andryushchenko E.S., Yakimov I.S., Bezrukovа O.E., Samoilo A.S. Full-profile x-ray phase analysis of the composition of aluminum electrolyte on the basis of a multipopulation evolutionary algorithm and elemental analysis data, J. Sib. Fed. Univ. Chem., 2018, 11(2), 262-272. DOI: 10.17516/1998-28360073.

© Siberian Federal University. All rights reserved

* Corresponding author E-mail address: dubinin-2005@yandex.ru

Полнопрофильный рентгенофазовый анализ состава электролита алюминиевого производства на основе мультипопуляционного эволюционного алгоритма и данных элементного анализа

П.С. Дубинин3, А.Н. Залога3, С.В. Бураков6, К.А. Гусев6, Е.С. Андрющенкоа, И.С. Якимова, О.Е. Безрукова8, А.С. Самойлоа

аСибирский федеральный университет Россия, 660041, Красноярск, пр. Свободный, 79 бСибирский государственный университет науки и технологий имени М.Ф. Решетнева Россия, 660037, Красноярск, пр. Красноярский рабочий, 31

Разработан метод количественного рентгенофазового анализа (КРФА) на основе самоконфигурируемого мультипопуляционного параллельного генетического алгоритма (СГА), который успешно применен для бесстандартного автоматического полнопрофильного анализа образцов электролита алюминиевых электролизеров. Метод позволяет избежать работы по эмпирическому выбору настроек генетического алгоритма и показывает достаточную точность автоматического КРФА по методу Ритвельда под управлением СГА (полученные значения криолитового отношения соответствуют аттестованным значениям с точностью, характеризуемой стандартным отклонением всего 0,035 ед. КО), а также возможность его использования для оперативного контроля состава электролита на алюминиевом производстве.

Ключевые слова: генетические алгоритмы, самоконфигурация, рентгеновская порошковая дифракция, метод Ритвельда, количественный рентгенофазовый анализ, электролит алюминиевого производства.

Введение

В настоящее время алюминий получают в электролизных ваннах разложением глинозема (Al2O3), растворенного в расплаве фторидов при температуре около 950 °С [1]. Основным компонентом расплава является натриевый криолит (Na3AlF6), к которому для улучшения технологических характеристик добавляют фторид алюминия (AlF3) и некоторые фториды щелочных и щелочноземельных металлов, чаще всего это фторид кальция (CaF2) и фторид магния (MgF2). В результате введения добавок изменяются важные для технологии электролиза свойства расплавов: электропроводность, плотность, температура ликвидуса, растворимость глинозема, давление паров над расплавом, а также коррозионное воздействие на материалы электродов, футеровку и другие элементы конструкции электролизера. В процессе работы электролизная ванна постоянно испытывает колебания и отклонения от оптимального химического состава.

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

Такимосновным технологическим параметром, характеризующим химический состав электролизной вынны, выступает криолитовое отношение (КО) - интегральнаи величина, выражающая величину мальнаго отношения фторида натрия к фториду алюминия или отношении ихмасс:

KO = C(NaF, маа.%) = 2 C(NaF, масс.%) C(AlF3, маа.%) C(AlF3, масс.% /

где С - содержания фторидов натрия и алюминия в молярных или массовых процентах.

Дополнительно, наряду с КО, контролируются концентрации добавок фторидов. Технологический экспрасс--анализ состава электролита производится количественным рвнтгенофа-зовым методсм (КРФА). Анализ проводится в лабораторных условияа на отобранных из ванн и охьыжденных пробах по измеряемым на аифрактометрах интенсивностям анаиитических линий криаталлическиа фаз. Необходимая частота контроля ванны - раз в 2-3 дня, точность контроля КО ~ 0,04, точность контроля фторидов ~ 10 % (отн.), время на анализ одной пробы -несколько минут.

Охлажденные пробы электролита российских алюминиевых заводов имеют достаточно сложный фазоаый состав, приведенный втабл. 1.

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

Таблица 1. Фазы проб охлажденного электролита российских алюминиевых заводов

Table 1. The electolyte at Russian aluminum smelters

№ Фаза Хим. формула Содержание (% масс.) Область КО

1 Криолит N a3AlF6 0-90 > 1.67

2 Хиолит Na5Al3F14 0-85 < 3.0

3 Фторид натрия NaF 0-5 > 3.0

4 Ca-криолит 1 NaCaAlF6 0-15 < 3.0

5 Ca-криолит 2 Na2CajAl2F14 0-20 < 2.95

6 Флюорит CaF2 0-9 > 2.45

7 Веберит Na2MgAlF7 0-15 < 2.85

8 Нейборит NaMgF3 0-6 > 2.5

9 а-, Р-, у- глинозем Al2O3 2-5

характеристикам, пробам промышленного электролита стандартных образцов; эта процедура достаточно сложная, трудоемкая и на практике не всегда возможная.

Нами для автоматического анализа фазового состава образцов электролита алюминиевого производства предлагается бесстандартный , структурно-чувствительный КРФА по методу Ритвельда на основе разработанного многопопуляционного эволюционного алгоритма [2]. Так как большинство аналитических лабораторий обеспечены соответствующим оборудованием и возможностью проведения химического анализа, например, рентгеноспектральным флуоресцентным методом, то для повышения точности КРФА привлекается информация о количественном элементном составе образцов. Разработанный метод является гибким и может использоваться для контроля состава любых кристаллических систем на любом современном промышленном производстве.

Экспериментальная часть

Метод самоконфигурирования генетических алгоритмов

Полнопрофильный анализ по методу Ритвельда широко используется для лабораторного количественного рентгенофазового анализа. Однако его приложения для производственного технологического контроля недостаточно развиты. Причина в том, что метод Ритвельда основан на нелинейном методе наименьших квадратов (МНК), который для сходимости требует установки корректных исходных приближений уточняемых параметров для каждого образца. В лабораторном варианте метода требования к точности исходных приближений ниже, так как имеется возможность их интерактивного поэтапного уточнения. Но для технологического контроля требуется высокая степень автоматизма и способность унифицированного анализа большого количества образцов. При этом автоматическое уточнение большого числа профильных и структурных параметров из их унифицированных исходных значений, не достаточно хорошо подходящих одновременно ко всем образцам, часто приводит к расходимости МНК. Одним из вариантов решения этой задачи является применение генетических алгоритмов (ГА) для выбора начальных значений параметров образца, эволюционного отбора перспективных параметров и их дальнейшего автоматического уточнения по МНК метода Ритвельда.

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

этапы самоконфигурируемого генетического алгоритма (СГА) безусловной оптимизации могут быть описаны следующим образом.

1. При старте алгоритма выбор того или иного варианта для каждого вида оператора (селекция, скрещивание, мутация) является равновероятным (вероятность применения каждого варианта, например, селекции равна p=1/z, где z - число операторов данного вида). То есть пока не накоплена статистика по эффективности применения, равновероятно применяются несколько вариантов селекции, несколько вариантов мутации, несколько вариантов скрещивания.

2. На каждом поколении оценка эффективности применения того или иного варианта операторов ГА производится по средней пригодности индивидов эволюции, получаемых с его помощью: averagefitnessi= f/n, i=1, 2, ..., z, где averagefitnessi - средняя пригодность порожденных потомков; f - сумма пригодностей индивидов, полученных i-м оператором; ni - количество индивидов, полученных i-м оператором (данного вида, например, селекции); z - число операторов.

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

4. Для порождения каждого потомка следующего поколения стохастически (согласно полученным распределениям вероятностей) выбираются операторы, при помощи которых он и будет сгенерирован.

Самоконфигурирование освобождает пользователя от подбора настроек генетического алгоритма, при этом эффективность решения задачи остается приемлемой (при наилучшем выборе параметров генетического алгоритма эффективность решения задач несколько выше, но подбор параметров ГА требует времени и высокой квалификации пользователя).

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

пространству поиска. На каждом из вычислительных узлов с помощью операторов рекомбинации и селекции эволюционно проирходит формирование потомков с меньшими значениями целевой функции. Операторы мутации случайным образом «разбрасывают» их по пространству поиска, иногда с уоеличением значе нил целевой функции. "Часть моделей с моньшим значением целевой функции уточняютая с помощью МНК метода Рртвельда. Затем, is результате общей селеьощи формируется новая популяция тесьовыр моделей-аотамков, в среднем, с лучшей пригодностью. Определанное количества лучших тестоаых моделей ив популяций на рабяяих узлах отправляатся нууправляющьй вычислитееьный уилл СУА. ЗВссс; эти решения, на-1сопле1И1иы^ на упрлоеяющем ьлле, сл<ср)тиуеиются л порядке умеьушения знярвния целевой функ-1сеи. Периодически члеть уичших решений!, накопленных к танном- по1толению эволюции на управляющем узле, пе°чайня отбирается и ряссышается обратно в популяции на рабочих узлах случайным образом. Такая умеренная миграцал обеспеиивает распространваие удачиых реша-ний ico по-1л-ля1^и.)-]\н и улуишает о(5п-ую с;нсод^]\-ост->.

СЛ^мо1-О]Е-(^^1Н0л^и0)(-]Зс^1си^ реализовяно для индиви-еуалене-х проуиесов у-еАЛ на рабочих узлах аогенсно вышяяпиааяояй лхеме СГА . ДД,.гия каждого ирэоинаиссЕ! оудяется ставдартный нибор генетических яперяторов - по несколько япереторов каждего типв 13 честности, одноточечное, двухточечное и равноморное скрещиванит; р-шговся и турнирная селекция с разными размерами турнира; низкая, средняя и вигсокаи веряятность мутации Распределение вероятностей для выбора нтнлеетнягo ьяпа оператора пприии :-<знЕаееряот.цп.ии:и следующего пяоллеиия данной популяции постоянно суею-иется в яользу более яффеитивных вьриaнroв cаII:<з;]элlт)0]рo:в зг счет минее эффективных. ВданнеН реаливтиии эффяктивность конкретных гепeтичтсних операторов оценивается только локально, для индивидуальных ГА на отдельных рабочих узлах, без учета их эффективности в процессах ГА на дру гих узлах. Последнее может являться способом повышения оящей эффективности эволюционнаго КРФА, но требкит изучения.

Метод) полнопрофильного КРФА в параллель-ям СГА

КРФА ли методу Ритвельда основан ни итерационней минимизации с помощью нелинейного МНК неоязки между расчетным и экспериментальным профилями дифрактограммы:

2

Ф(0р) = £( ) - -сир -е Арр0г )) -е 0, (2)

i

едя Yo, Yc - экспериментальные и выяисленные инoерoивнoгти в точих профиля 26р we - весовые коэффиииенеы; Рт - вектои зна-ееий пр-оффиjikjh:^.:?-, миироструктуеных и кристоллострук-турныгх параметров, пОлру31^ннь>г;х на &-а итррсции; АОт - пргиращения к ним, вычисляемые по МНК; при k=0 задаются исходные значения всех параметров.

Уточняамая пот методи Ритвельда ча-тт пареметров oекяoрa P определяат состая параметрических строк - индивидов эволюции по иУ. Полный н-бор параметров веяи^а P, включея фиксированные и уточняемые, представляет тео^ую модрло харии-ерш^к многофазного образца. При выпа-нении эворюцианносо КРФАнеобхидимо задавать область значений уточняемых параматров. Лучшие из еначиний, ааходимьнх в этой облетти в процессе эколюции по ГА, используются в качестве исходных значений для уточнения по метод- Ритвельда. Для аа-дания ипласти титря на -триметры накладываются огранттчения - для 1т^кдо^о параметра

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

Для автаматизации полнопрофильного КРФА важно то, что всегда известно целевое значение Я-фактора, предварительно определяемое в результате декомпозиции профиля дифрак-тограммы без данных о кристаллической структуре фаз (обычно по методу ЛеБэйла [5]). После достижения одной из тестовых моделей значения Я-фактора, близкого к целевому, выполняется ее финальное уаочнение по метода Ритвельда и производится расчет концентроции фаз Са в образце. При отсу тствии аморфной фазы для расчета используется нормированное выражение:

N

Са= %а1аы уа / £ sjzjмJvj, (3)

где 5а - масштабный коэффициент расчетной дифрантограммы Уса фазы а, который выбирается для ее оптимального соответствия экспериментальной дифрактограмме Го; Уа, 2а, Ма - объем ячейки, число формульных единиц в ячейке и молекулярная масса фазы а; N -число фаз многофазного образца.

В присутствии аморфной фазы используется вариант КРФА с внутренним стандартом [6].

Про(5оемой эволюционного КРФА может являться изость области поиска, в которой м-фактор изменяотоя недостаточно значительно, чтобы1 служить эффективным критерием пригодности для эволюционного поиска по ГА. Поэтому для нсиления чувствительности целевой функции пригодности к Я-фактору добавлена невязка между выгаисленнылми ио данных КРФА и измеренными значениями элементного состава:

Rp = 100-

I w -Е-сх]2

I w,-[Го, - Ye, (P)]2 ^ SaZ0MaVa '

(1 - WXA ) -

Iw, - (YO,)2 " Iw - (Cf )2

,

(4)

где C,Xji - измеренные значения элемента t в образце; Pta - массовая доля химического элемента t в фазе а; wXA - весовой вклад в Rwp данных о химическом составе образца (обычно 0,5).

Интеграция данной реализации полнопрофильного анализа по методу Ритвельда с параллельным СГА обеспечивает автоматическое выполнение КРФА. При этом целевой функцией СГА является профильный R-фактор метода Ритаельда.

Разработанный СГА был применен для анализа фазового состава и КО комплекта из 24-х отраслевых стандартных образцов (ОСО) электролита, к настоящему времени внедренного на пяти Российских алюминиевых заводах. Данные ОСО были выбраны по нескольким причи-

и

нам: во-первых, они приготовлены из материала самого электролита нескольких алюминиевых заводов, т.е. отвечают составу и микроструктуре фаз реального промышленного электролита; во-вторых, ОСО аттестованы по химическому и фазовому составу с высоким качеством, стандартное отклонение данных аттестации составляет, в среднем 0.008 ед. КО; в-третьих, фазовый состав комплекта варьируется в широком интервале, например КО меняется от ~1,9 до 3 [7].

Порошковые дифракционные данные получены с использованием CuKa - излучения на диф -рактометре Shimadzu-7000S. Образец растирался в агатовой ступке и подготавливался методом прямой набивки кюветы. Условия съемки: диапазон от 10 до 90° по шкале 29 с шагом 0,01°, At -2,4 с, скоростью - 0,25°/мин. Структурные модели фаз брались из базы данных ICSD [8].

Метод СГА полностью автоматический, поэтому для анализа всех образцов использовался идентичный избыточный список фаз (см. табл. 1) и одинаковые настройки СГА. Для заданных фаз определялись значения параметров, указанные в табл. 2 (всего 90 параметров). Если фаза отсутствовала в образце, ее масштабный коэффициент обнулялся и концентрация, соответственно, тоже. Для каждого образца проводили по три запуска СГА и далее брали усредненные значения их фазовых концентраций. Для КРФА одного образца по программе СГА затрачивалось около 4-5 мин.

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

Для работы СГА были выбраны следующие генетические операторы, вероятность использования которых адаптивно варьировалась в ходе эволюционного процесса при вычислении каждого образца:

Таблица 2. Уточняемые параметры по СГА для образцов электролита Table 2. Refinable SGA parameters for electrolyte samples

Фаза Перечень уточняемых параметров

Параметр коррекции сдвига нуля дифрактограммы

NasAlF6 S, a,b,c,p,U,W, Eta0, Asyml, координаты атомов (12 позиций), текстура по Марч-Далласу по направлению hkl [112]

NasAlsF14 S, a,c,U,W, Eta0, Asyml, координаты атомов (9 позиций)

NaCaAlF6 S, a,b,c,p,U,W, Eta0, Asyml

Na2Ca3Al2F14 S, a, U,W, Eta0, Asyml

CaF2 S, a, U,W, Eta0, Asyml

Na2MgAlF7 S, a,b,c,U,W, Eta0, Asyml

NaMgFs S, a,b,c,U,W, Eta0, Asyml

a-Al2O3 S, a,b,U,W, Eta0, Asyml

y-AlOs S, a,U,W, Eta0, Asyml

5 — масштабный коэффициент фазы; а,Ь,с,в — параметры кристаллической ячейки; и,Ж — параметры полуширины пика (по модели псевдо-Войгта); Eta0 и А.чут1— параметры формы пика и асимметрии соответственно.

- операторы селекции: турнирная с размером турнира 3, 5, 7, 9, ранговая;

- операторы скрещивания: двухточечное, трехточечное, равномерное;

- операторы мутации: низкий уровень, средний уровень, высокий уровень, с тремя значениями среднеквадратического отклонения каждый.

Кроме этого, была применена гибридизация генетического алгоритма с локальным спуском по Ламарку каждые 10 поколений [9].

Результаты и обсуждения

В качестве примера итоговых результатов, получаемых при полнопрофильном КРФА по параллельному СТА, на рис. 1 представлен типичный результат моделирования дифрактограм-мы одного из образцов электролита. Профильный Р-фактор составил 8,6 % отн.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

На рис. 3 сравниваются аттестованные значения КО с КО, рассчитанными из фазового состава, полученного по СТА. График в идеале должен проходить через нуль под углом 450, а отклонения от него будут свидетельствовать о наличии систематической межметодической погрешности. Для количественной оценки ее существования на рис. 3 представлено также регрессионное линейное уравнение (y=a+bx). Соответственно, стандартное отклонение SD показывает величину случайной межметодической погрешности. Эти расчеты и график выполнены с помощью программы статистической обработки аналитических данных ORIGIN.

Рис. 1. Результаты моделирования дифрактограммы одного из образцов электролита по параллельному СТА с разностбойкривой

Fig. 1. Model powder pattern calculated by the parallel SGA for an electrolyte sample and the difference curve

л Я 1П1Ь2лиззв<с«5в«ЬЕ:]«>га75 1ва1к»

Номер тсжоления

Рис. 2. Типичный график сходимости анализа образца электролита по СГА

Fig. 2. A typical SGA convergence graph obtained during the analysis of an electrolyte sample

KO OCO

Рис. 3. График соответствия расчетных и аттестованных КО отраслевых стандартных образцов электролита (КО OCO - значения аттестованных КО, КО SGA - значения КО, рассчитанных из фазового состава по СГА; SD - стандартное о тклонение (% масс.); R - коэф фициент корреляции ; регрессионное уравнение y=a+bx)

Fig. 3. Correspondence between the calculated and certified values of CR for branch reference samples of the electrolyte (KO OCO - certified values of CR, KO SGA - values of CR calculated from the; phase composition in SGA, SD - standard deviation (% mass), R - correlation coefficient, regression equation y = a + bx)

Вывод] ы

1. Определенные по СГА значения КО соответствуют аттестованным значениям с точностью, характеризуемой стандартным отклонением 0,035 ед. КО и укладываются в 95%-й интервал надежности, изображенный пунктиром на рис. 3.

2. Регрессионное уравнение хорошо соответствует уравнению y = a+bx , при этом коэффициент а статистически незначимо отличен от 0, но b завышен на 3,5 % отн. Последнее показывает, что имеет место небольшое систематическое завышение результатов

СГА. Это вызвано завышением концентрации фазы криолита Na3AlF6 и свидетельствует о неполноте автоматического уточнения ее кристаллической структуры по СГА, искажаемой вследствие неполного перехода высокотемпературной модификации криолита в низкотемпературную при неравновесной кристаллизации проб электролита в процессе пробоотбора на заводах.

3. В целом, полученные характеристики соответствуют технологическим требованиям, предъявляемым на алюминиевом производстве к точности анализа КО, и метод СГА может быть рекомендован к внедрению для оперативного контроля КО электролита после устранения причин отмеченной систематической погрешности.

Исследование выполнено при поддержке Российского фонда фундаментальных исследований, Правительства Красноярского края, Красноярского краевого фонда поддержки научной и научно-технической деятельности в рамках проекта № 16-43-243049.

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

1. Apisarov A., Frolov A., Gusev A. et al. Modiefid alumina-cryolite bath with high electrical conductivity and dissolution rate of alumina. Light metals 2007. P. 571-576.

2. Zaloga A., Akhmedova S., Yakimov I., Burakov S., Semenkin E., Dubinin P., Piksina O., Andryushchenko E. Genetic Algorithm for Automated X-Ray Diffraction Full-Profile Analysis of Electrolyte Composition on Aluminium Smelters. Lecture Notes in Electrical Engineering 2016. Vol. 383, P. 79-93.

3. Semenkin E.S., Semenkina M.E. The Choice of Spacecrafts' Control Systems Effective Variants with Self-Configuring Genetic Algorithm. In: Ferrier, J.L. et al (Eds.): Informatics in Control, Automation and Robotics: Proceedings of the 9th International Conference ICINCO '2012. Rome: Italy, 2012. Vol. 1, P. 84-93.

4. Semenkin E.S., Semenkina M.E. Self-configuring Genetic Algorithm with Modified Uniform Crossover Operator. Advances in Swarm Intelligence. Lecture Notes in Computer Science 7331. Springer-Verlag, Berlin Heidelberg 2012. P. 414-421.

5. Le Bail A. The Profile of a Bragg Reflection for Extracting Intensities. Chapter 5. Powder Diffraction Theory and Practice, ed. R.E. Dinnebier and S.J.L. Billinge. Cambridge, Royal Society of Chemistry, 2008. 571 p.

6. Bish D.L., Howard S.A. Quantitative phase analysis using the Rietveld method. Journal of Applied Crystallography 1998, Vol. 21, P. 86-91.

7. Якимов И.С., Дубинин П.С., Залога А.Н., Пиксина О.Е., Кирик С.Д. Разработка отраслевых стандартных образцов электролита алюминиевых электролизеров. Стандартные образцы, 2008 (4), С. 30-42 [Yakimov I.S., Dubinin P.S., Zaloga A.N., Piksina O.E., Kirik S.D. Development of industry standard electrolyte samples of aluminum electrolyzers. Certified Reference Materials 2008, Vol. 4, P. 30-42. (In Russ.)]

8. Inorganic Crystal Structure Database (ICSD) http://www2.fiz-karlsruhe.de/icsd_publications.

html.

9. M. Tim Jones. AI Application Programming (Programming Series). Charles River Media, 2003. 363 p.

i Надоели баннеры? Вы всегда можете отключить рекламу.