Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.geol.msu.ru/deps/geochems/soft/ru/OptimC.doc Дата изменения: Mon Dec 28 17:51:57 2015 Дата индексирования: Sun Apr 10 01:18:10 2016 Кодировка: koi8-r

Программа OptimC


Руководство пользователя


О программе 1
Установка программы 2
Файлы и папки 2
О создании файлов данных 2
Начало работы с программой OptimC 2
Исходные данные 3
Первая часть данных 3
Вторая часть данных 4
Результаты расчётов 4
Оптимизация RB-параметров 5
Задание набора параметров для оптимизации 5
Задание дополнительных ограничений 5
Дополнительные возможности 6
Таблица значений pK диссоциации комплекса 6
Экспорт для Unitherm 6
Литература 7
Приложение. Сообщения об ошибках 8
Ошибки, возникающие в начале работы 8
Ошибки, возникающие при чтении файла данных 8
Ошибки, возникающие при расчётах или после них 8



О программе

Программа OptimC является Windows версией программы UT-RYZ для DOS.
Назначение программы - нахождение параметров модифицированной модели
Рыженко-Брызгалина (RB) [1, с. 74], описывающей термодинамические свойства
комплексных частиц водного раствора. Параметры этой модели pK298, A(zz/a),
B(zz/a) рассчитываются в программе методом наименьших квадратов по
экспериментальным значениям константы диссоциации комплекса, измеренным при
различных значениях температуры и давления (или температуры и плотности
воды). Оптимизироваться могут как все параметры одновременно, так и
произвольное их подмножество.
Помимо безусловной оптимизации указанных параметров возможна и их
условная оптимизация, при которой учитываются дополнительные ограничения, а
именно: в некоторых точках вместо константы диссоциации может быть
определено лишь её максимально возможное значение.
Программа OptimC работает под управлением операционной системы Windows
(от 98 и выше) и требует наличия установленной программы Microsoft Excel
(версии 8.0 и выше). Используемые программой алгоритмы расчёта
термодинамических свойств воды [2] и других термодинамических функций
полностью соответствуют используемым в программном пакете HCh (в частности,
в программе Unitherm). По этой причине оптимизированные программой
параметры могут непосредственно вноситься в базу термодинамических данных
Unitherm. Тем не менее, для работы программы OptimC наличие установленного
на компьютере пользователя пакета HCh необязательно.

Установка программы

Для установки программы на компьютер запишите в отдельную («домашнюю»)
папку на жестком диске два файла - исполняемый файл OptimC.exe и шаблон для
создания файлов данных OptimC.xls. Если на компьютере установлен пакет HCh,
рекомендуется записать эти файлы в папку C:\HCh_Win\Main (это облегчит
включение оптимизированных параметров частиц в базу термодинамических
данных Unitherm).
Создайте ярлык исполняемого файла OptimC и поместите его на рабочий
стол или в любое другое удобное для вас место. Если на компьютере
установлен пакет HCh, можно также поместить ярлык OptimC в группу HCh for
Windows меню Пуск. В дальнейшем с помощью этого ярлыка вы сможете запускать
программу OptimC обычным образом.

Файлы и папки

В качестве файла исходных данных программа OptimC использует рабочую
книгу Excel (точнее, один из её рабочих листов) и туда же записывает
результаты работы. Вы можете самостоятельно подготовить исходные данные с
помощью программы Excel, а затем запустить программу OptimC, передав ей
подготовленный файл. Однако, возможно, удобнее поступить по-другому:
запустить программу OptimC, нажать в ней кнопку "Create new data file" и
заполнить требуемые поля в открывшемся шаблоне. В обоих случаях записать
файл на диск можно только с помощью программы Excel (меню Файл, команда
«Сохранить как...»). Никогда не записывайте файл данных в домашнюю папку
программы OptimC. Используйте для таких файлов другую, «рабочую» папку.
Программа OptimC запоминает рабочую папку и при следующем запуске откроет
её автоматически.

О создании файлов данных

Для создания файла данных удобнее всего воспользоваться шаблоном.
Шаблон - это рабочая книга Excel (OptimC.xls), состоящая из одного листа,
содержащего все нужные ключевые слова и размеченного под данные для
программы OptimC. Ключевые слова - это слова, которые описывают данные и
которые вы не должны изменять (если вы измените какое-либо ключевое слово,
программа OptimC не будет работать с таким файлом данных). Каждая ячейка
шаблона, содержащая ключевое слово, снабжена комментарием, разъясняющим
смысл данных, описываемых этим ключевым словом. Шаблон имеет атрибут
«только для чтения», поэтому после заполнения всех необходимых полей файл
данных следует записать под другим именем (и расширением .xls).
Вы можете создать файл данных для программы OptimC в виде рабочей книги
Excel c несколькими рабочими листами, каждый из которых соответствует своей
частице водного раствора. Создайте дополнительный рабочий лист в книге
обычным образом командой Excel «Добавить рабочий лист», после чего
скопируйте шаблон (или уже заполненный рабочий лист книги), вставьте его в
новый рабочий лист и отредактируйте.

Начало работы с программой OptimC

При запуске программа OptimC первым делом проверяет наличие в памяти
компьютера работающей программы MS Excel и, если её нет, запускает её, не
открывая никакой рабочей книги. Если программа Excel уже работает во время
старта OptimC, новой копии Excel не создаётся. Если в системе не
установлена программа Excel или версия этой программы слишком старая,
программа OptimC не стартует. Если же с программой Excel всё в порядке, на
экране появляется сокращённое окно программы OptimC содержащее всего две
кнопки, которые позволяют открыть нужный файл данных или создать новый.
Кнопка "Open existing data file" позволяет открыть уже существующую на
диске рабочую книгу Excel с исходными данными задачи. Эта кнопка доступна
всегда. При нажатии этой кнопки на экране появляется стандартный диалог
открытия файла, в котором вы можете выбрать папку и необходимую рабочую
книгу Excel (расширение .xls). После нажатия кнопки «Открыть» выбранная
рабочая книга открывается, и программа OptimC начинает читать информацию с
текущего рабочего листа. При чтении данных программа проверяет их на
наличие ошибок. Если ошибка обнаружена, выдается диагностирующее сообщение
(с указание строки и столбца, где обнаружена ошибка) и чтение данных
прекращается - в этом случае необходимо исправить указанную ошибку на
текущем рабочем листе и вернуться к программе OptimC для повторного
считывания данных. Если при открытии файла данных ошибок не обнаружено,
программа OptimC сообщает, какое количество экспериментальных точек
содержится в файле, и открывает своё окно полностью.
Кнопка "Create new data file" позволяет создать новый файл данных,
вводя исходную информацию в шаблон рабочего листа Excel. Эта кнопка
доступна только в том случае, если в домашней папке программы OptimC
присутствует файл OptimC.xls (см. раздел «Установка программы»). При
нажатии этой кнопки шаблон загружается в Excel и начинает считываться
программой как обычный файл данных. Так как файл шаблона не содержит данных
(в том числе и обязательных), то при его чтении программа обнаруживает
ошибку - отсутствие обязательных данных - и выдает соответствующее
сообщение. Вам нужно исправить эту «ошибку», то есть ввести необходимые
данные и вернуться к программе OptimC. Далее вы работаете с шаблоном так
же, как и с обычными файлами данных. Единственное, что нужно помнить:
записывать полученный файл на диск следует из Excel (меню Файл, команда
«Сохранить как...»).
Признаком того, что рабочий файл открыт, и в нём нет ошибок, является
полное окно программы OptimC; при этом заголовок окна содержит имя текущей
рабочей книги и текущего рабочего листа. Полное окно программы содержит
дополнительно несколько элементов управления: набор переключателей,
позволяющих выбрать RB-параметры, подлежащие оптимизации, и командную
кнопку для её запуска. Если текущий рабочий лист содержит дополнительные
ограничения на величины констант диссоциации, то в окне программы OptimC
имеется также переключатель режима условной оптимизации.

Исходные данные

Исходные данные для программы OptimC можно разделить на две части: RB-
параметры и другие свойства комплекса (строки 1-12) и таблица
экспериментальных значений константы (pK) его диссоциации (остальные
строки). Ниже описаны требования, которым должны удовлетворять эти данные.

Первая часть данных

Первая строка рабочего листа предназначена для заголовка - общего
описания данных на рабочем листе. Программа заголовок не использует.
Вторая строка: ячейка C2 должна содержать название комплекса, свойства
которого представлены на данном рабочем листе. Заполнение этой ячейки
обязательно. Если вы планируете включить оптимизированные параметры этого
комплекса в базу данных Unitherm, его название должно удовлетворять
соглашениям, принятым в этой базе для частиц водного раствора, то есть
представлять собой химическую формулу комплекса и не содержать более 14
символов.
Строки 4-7: введите в ячейки B4-B7 и C4-C7 соответственно
стехиометрические коэффициенты и названия составляющих комплекс простых
частиц. Тем самым вы зададите реакцию диссоциации комплекса, параметры
которой будут представлены вами ниже в этом файле. Максимальное число
частиц, на которые может диссоциировать комплекс, равно четырем (без учета
стехиометрических коэффициентов). Если комплекс диссоциирует на меньшее
число частиц, неиспользованные ячейки должны оставаться пустыми. Если вы
планируете в дальнейшем включить свойства данного комплекса в базу данных
Unitherm, помните, что все образующие комплекс частицы должны содержаться в
этой базе данных, и каждая из них должна принадлежать к одному их следующих
типов: 1) встроенные - частицы H2O, H+ и OH-; 2) базисные - частицы,
описываемые моделью HKF (Helgeson-Kirkham-Flowers); или 3) водорастворённые
неэлектролиты - частицы, описываемые моделью AD (Akinfiev-Diamond).
Информация о стехиометрии реакции диссоциации комплекса не используется
программой OptimC при оптимизации её параметров, однако она полезна как
конкретизация реакции диссоциации и необходима для правильного
функционирования базы данных Unitherm.
Строки 9-11: введите в ячейки B9-B11 известные вам значения RB-
параметров комплекса (pK298, A, B), а в следующие ячейки (столбец C) -
ссылки на источники этих данных. Если значение какого-либо из параметров
вам неизвестно, можно оставить соответствующую ячейку пустой, но лучше всё-
таки ввести какое-нибудь число (напр., 0).

Вторая часть данных

Вторая часть исходных данных содержит экспериментальные значения
константы диссоциации (точнее, её отрицательного логарифма) при различных
условиях. Строки 13 и 14 содержат заголовок таблицы; здесь вводить ничего
не надо.
Собственно данные начинаются со строки 15; одна строка описывает одну
экспериментальную точку. Каждая точка определяется следующими параметрами:

- температура в градусах Цельсия; должна задаваться обязательно;

- давление в барах (значение "0" задает давление насыщенного пара воды);
если давление неизвестно (не измерялось), ячейку нужно оставить пустой (см.
ниже);

- плотность воды в граммах на кубический сантиметр; плотность воды задаётся
как альтернатива заданию давления (см. выше) - для каждой экспериментальной
точки может (и должна) быть задана лишь одна из этих величин;

- вес измерения - число, обратное стандартной ошибке измерения; допускается
нулевое (или, эквивалентно, пустое) значение - такие точки не учитываются
при оптимизации; вместо веса измерения в ячейки этой графы можно вводить
символ "<" (без кавычек) - этот знак говорит программе OptimC, что
следующая ячейка содержит не значение pK диссоциации, а лишь верхнюю
границу для него;

- экспериментальное значение pK диссоциации (или его верхний предел - см.
выше); ввод этого значения обязателен;

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

Результаты расчётов

Результаты работы выводятся программой OptimC на тот же рабочий лист в
ячейки, расположенные справа от исходных данных. Ниже приводится их краткое
описание.
Ячейки E9-E11 ("init.") содержат значения RB-параметров комплекса перед
оптимизацией. Если при открытии файла данных эти ячейки пусты, они
заполняются программой OptimC путем копирования соответствующих данных из
ячеек B9-B11. Если же в какой-либо из этих ячеек уже имеется значение, оно
сохраняется неизменным. При оптимизации RB-параметров именно эти значения
являются исходными. При желании вы можете сами ввести сюда любые значения.
Ячейки F9-F11 ("opt.") содержат значения RB-параметров комплекса после
оптимизации. Они заполняются программой результатами расчётов, выполненных
после нажатия кнопки "Optimize.". Вводить сюда значения вручную не имеет
смысла.
Ячейки E12 и F12 содержат взвешенные среднеквадратические ошибки
аппроксимации экспериментальных значений pK диссоциации комплекса с помощью
исходного и оптимизированного набора RB-параметров соответственно. Первая
из них вычисляется программой OptimC до, а вторая - после оптимизации
параметров. Эти ячейки не содержат значений, если в файле отсутствует
таблица экспериментальных значений pK.
Ячейки G9-G11 содержат доверительные интервалы (на 95% уровне) для
значений оптимизированных параметров, если соответствующие параметры
оптимизировались, и не содержат ничего (пусты), если они не
оптимизировались (сравните в этом случае соответственные числа в колонках E
и F).
Столбцы F и G справа от таблицы экспериментальных данных содержат
величины, полученные программой OptimC при последней оптимизации параметров
комплекса:

- F ("opt.") - значения pK, рассчитанные по набору оптимизированных RB-
параметров;

- G ("error") - разности между оптимизированными (F) и экспериментальными
(E) значениями pK. Именно эти значения используются при вычислении
взвешенной среднеквадратической ошибки аппроксимации для оптимальных данных
(F12).

Оптимизация RB-параметров

Когда кнопка "Optimize RB parameters" в окне программы OptimC доступна,
её нажатие запускает процесс поиска таких значений выбранных вами RB-
параметров комплекса, которые в наибольшей степени соответствуют
экспериментальным данным (при условии выполнения дополнительных
ограничений, если это требуется). Таким образом, перед нажатием этой кнопки
необходимо сделать две вещи: задать набор параметров, которые требуется
оптимизировать и указать, требуется ли учёт дополнительных ограничений
(если они есть).

Задание набора параметров для оптимизации

Параметры, которые можно оптимизировать, представлены в главном окне
программы OptimC. Для того, чтобы включить параметр в список
оптимизируемых, достаточно его отметить (поставить галочку слева от него).
Соответственно, снятие галочки исключает параметр из списка.
Замечание. Справа от параметра в главном окне программы OptimC указано
его текущее значение. Если значение отсутствует (т. е. стоит знак "?"), то
исключить такой параметр из списка оптимизируемых не удастся.
Рекомендации. Обращайте внимание на величину доверительного интервала
для параметров, полученных в результате оптимизации. Например, если
доверительный интервал для параметра B оказывается больше его значения (по
модулю), то это говорит о том, что параметр B не значим, и, следовательно,
имеет смысл сделать его равным нулю и исключить из процесса оптимизации.
Выбор оптимизируемых параметров можно также значительно облегчить, если на
том же рабочем листе средствами Excel построить график типа XY Scatter (без
соединительных линий) по данным столбцов E и F таблицы экспериментальных
данных (на этом графике вдоль оси X будут отложены экспериментальные, а
вдоль оси Y - расчётные значения pK). По тому, как множество точек на этом
графике отклоняется от диагонали Y=X, хорошо видно влияние каждого из
параметров на качество аппроксимации. Экспериментируйте, благо параметров
немного!

Задание дополнительных ограничений

Если среди ваших экспериментальных данных имеются точки, в которых
значение pK диссоциации комплекса не определено из-за слишком низкого
значения (т. е. его значение оказалось ниже предела чувствительности
метода), то даже такие точки можно использовать с пользой для дела. Вы
можете включить такие данные в общую таблицу, введя вместо веса измерения
символ "<", а вместо значения pK - его верхний предел (эти «неполноценные»
экспериментальные данные будут выделены в таблице цветом). Если не
указывать особо, то при оптимизации программа OptimC не будет использовать
такие точки (так же, как и точки с нулевым весом), однако расчётные
значения pK и их отклонения от заданных значений будут выдаваться для всех
точек без исключения. Поэтому вы всегда можете проверить, удовлетворяются
ли для оптимизированных параметров условия pKрасч < pKmax или нет (если
удовлетворяются, то значения в графе "error" будут отрицательными для всех
таких точек). Если в некоторых точках ограничения "<" не удовлетворяются
(значения в графе "error" положительны), и вас это не устраивает, вы можете
повторить оптимизацию, предварительно включив требование жёсткого
выполнения всех таких условий. Для этого достаточно поставить галочку у
"Activate the "<" restrictions" и щёлкнуть на кнопке "Optimize.". В
результате этой операции все ограничения на значения pK будут
удовлетворены, однако, конечно, суммарная погрешность аппроксимации
возрастёт. Какой из вариантов расчёта оставить в качестве окончательного
результата, конечно, решать вам.

Дополнительные возможности

В данном разделе описаны возможности, недоступные через графический
интерфейс программы OptimC.

Таблица значений pK диссоциации комплекса

Оптимизировав RB-параметры, вы, возможно, захотите узнать, как эти
параметры позволяют предсказать pK диссоциации вашего комплекса в более
широком диапазоне условий, чем в исходных экспериментальных точках.
Программа OptimC предоставляет возможность получить стандартную (для
программы Unitherm) таблицу таких значений (для температур до 500 њС и
давлений до 5 кбар).
Для того, чтобы увидеть эту таблицу, достаточно нажать клавишу F9
(программа OptimC должна быть в этот момент активна). Прежде, чем закрыть
окно с таблицей, вы можете скопировать всю таблицу в буфер обмена либо с
помощью локального меню (правая кнопка мыши), либо с помощью «горячих»
клавиш. При нажатии комбинации клавиш Ctrl+Insert таблица помещается в
буфер обмена в формате, удобном для последующей вставки в MS Excel; при
нажатии клавиш Ctrl+C таблица помещается в буфер обмена в формате, удобном
для последующей вставки в MS Word. В последнем случае, после вставки
таблицы в Word для лучшего представления рекомендуется задать для неё
непропорциональный шрифт, например, Courier.

Экспорт для Unitherm

Оптимизированные параметры комплекса можно записать в файл data.txt в
формате, пригодном для непосредственного импорта в базу данных Unitherm.
Для этого достаточно нажать клавишу F2 - данные будут добавлены к этому
файлу, если он присутствует в «домашней» папке программы OptimC; если же
файла data.txt в этой папке нет, он будет создан. Если «домашней» папкой
программы OptimC является папка C:\HCh_Win\Main (см. раздел «Установка
программы»), то последующий импорт данных в Unitherm не потребует никаких
дополнительных манипуляций с этим файлом (разумеется, если имя частицы
удовлетворяет соглашениям, принятым в HCh).
Замечание. Программа OptimC учитывает региональную настройку вашего
компьютера на использование символа-разделителя целой и дробной части
чисел. Так, если ваш компьютер настроен на использование десятичной
запятой, все числа в программе представляются именно с ней, так как Excel
использует то же представление чисел. Однако, в файл data.txt программа
всегда выводит числа с десятичной точкой, так как пакет HCh понимает только
это представление.



Если необходимо сослаться на эту программу, можно воспользоваться ссылкой
[3].




Литература


1. М. В. Борисов, Ю. В. Шваров. Термодинамика геохимических процессов:
Учебное пособие // М., Изд-во МГУ, 1992, 256 с.
2. J. Kestin, J. V. Sengers, B. Kamgar-Parsi, and J. M. H. Levelt Sengers
(1984). Thermophysical properties of fluid H2O. // J. Phys. Chem. Ref.
Data, 13, #1, pp. 175-183.
3. Shvarov Yu.V. A suite of programs, OptimA, OptimB, OptimC, and OptimS
compatible with the Unitherm database, for deriving the thermodynamic
properties of aqueous species from solubility, potentiometry and
spectroscopy measurements // Applied Geochemistry, 2015, v. 55, p. 17-27.


Приложение. Сообщения об ошибках

|Сообщение об ошибке |Комментарий |
|Ошибки, возникающие в начале работы |
|Could not run Excel |На компьютере не установлен MS Excel или его |
| |версия ниже 8.0 (Office 97). Установите |
| |современный Excel |
|Could not open the file|Ошибка открытия файла. Проверьте, не занят ли |
| |файл другой программой |
|Do not work in the home|Вы пытаетесь записать рабочий файл в домашнюю |
|folder! |папку, что не разрешается. Запишите файл в |
| |рабочую папку |
|Ошибки, возникающие при чтении файла данных |
|Format error or |Программа обнаружила, что в указанной ячейке |
|incorrect keyword |содержится недопустимое значение. Исправьте |
| |ошибку |
|Missing mandatory data |Указанная ячейка должна содержать данные, однако |
| |она пуста. Введите необходимое значение |
|T/P/D out of range |Температура, давление или плотность воды имеет |
| |значение, недопустимые в модели HKF. Исправьте |
| |параметры или запустите OptimC с опцией /ie (на |
| |свой страх и риск) |
|P(sat) does not exist |При указанной температуре давления насыщенного |
| |пара воды не существует. Исправьте температуру |
| |или давление |
|Weight is negative |Обнаружен отрицательный вес экспериментальной |
| |точки. Исправьте ошибку |
|Invalid data format |Указанная ячейка содержит значение, которое не |
| |является числом. Исправьте ошибку |
|Cell must be empty |Указанная ячейка не должна содержать никаких |
| |значений. Очистите её и продолжите работу |
|Wrong complex |В описании реакции диссоциации комплекса задана |
|composition |либо частица без стехиометрического коэффициента,|
| |либо стехиометрический коэффициент без частицы |
|No complex composition |Реакция диссоциации комплекса в файле не задана. |
| |Задайте реакцию диссоциации |
|Define pressure or |В указанной строке не задано ни давления, ни |
|density |плотности воды. Задайте нужный параметр |
|Could not calculate |Не удалось вычислить плотность воды при заданных |
|density |температуре и давлении. Проверьте правильность |
| |данных |
|OptimC closed: could |После возврата управления программа OptimC не |
|not find an open |нашла ни одной открытой рабочей книги Excel. Для |
|worksheet |нормального завершения работы сначала закрывайте |
| |OptimC, и лишь потом Excel |
|Unexpected error |Непредвиденная ошибка. Сообщите об этом |
| |разработчику |
|Ошибки, возникающие при расчётах или после них |
|Could not find solution|Попробуйте уменьшить число оптимизируемых |
| |параметров |
|Could not copy the |Системная ошибка копирования данных. Код ошибки |
|table. Error code ... |указан в сообщении |
|Could not open file |Программа не смогла открыть файл data.txt. |
|'data.txt' (is it |Проверьте, не используется ли этот файл другой |
|busy?) |программой |
|Could not write to file|Программа не может записать данные в файл |
|'data.txt' (is it |data.txt. Проверьте, не защищен ли файл от записи|
|protected?) | |