Геохимическое моделирование - Geochemical modeling

Геохимическое моделирование это практика использования химическая термодинамика, химическая кинетика или оба, чтобы проанализировать химические реакции это влияет геологические системы, обычно с помощью компьютера. Используется в высокотемпературных геохимия для моделирования реакций, происходящих глубоко в недрах Земли, в магма, например, или для моделирования низкотемпературных реакций в водные растворы у поверхности Земли, предмет данной статьи.

Приложения для водных систем

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

Развитие геохимического моделирования

Гаррели и Томпсон (1962) впервые применили химическое моделирование к геохимии при 25 ° C и общем давлении в одну атмосферу. Их расчет, произведенный вручную, теперь известен как модель равновесия, который предсказывает распределение частиц, состояние минеральной насыщенности и летучесть газа на основе измерений объемного состава раствора. Удаляя небольшие аликвоты растворитель воды из уравновешенной родниковой воды и неоднократно пересчитывая распределение видов, Гаррелс и Маккензи (1967) моделировали реакции, происходящие при испарении родниковой воды.[3] Эта связь массопереноса с моделью равновесия, известная как модель пути реакции, позволил геохимикам моделировать процессы реакции.

Helgeson (1968) представил первую компьютерную программу для решения моделей равновесия и пути реакции,[4] которые он и его коллеги использовали для моделирования геологических процессов, таких как выветривание, осадок диагенез, испарение, гидротермальные изменения, и отложение руды.[5] Более поздние разработки в области геохимического моделирования включали переформулировку определяющих уравнений, сначала как обыкновенные дифференциальные уравнения, а затем как алгебраические уравнения. Кроме того, химические компоненты стали представлены в моделях водными веществами, минералами и газами, а не элементами и электронами, из которых состоят вещества, что упростило основные уравнения и их численное решение.[2]

Последние улучшения в мощности персональных компьютеров и программное обеспечение для моделирования сделали геохимические модели более доступными и более гибкими в их реализации.[6] Теперь геохимики могут построить на своих ноутбуках сложные пути реакции или реактивные транспортные модели для чего раньше требовался суперкомпьютер.[7]

Настройка геохимической модели

Водная система однозначно определяется своим химическим составом, температура, и давление.[8] Создание геохимических моделей таких систем начинается с выбора основы, набора водные виды, минералы, и газы которые используются для записи химических реакций и выражения состава. Требуемое количество базовых записей равно количеству составные части в системе, что фиксируется правило фазы термодинамики. Как правило, основу составляют вода, каждый минерал находится в равновесии с системой, каждый газ находится в известном состоянии. летучесть, и важные водные виды. Как только основа определена, разработчик модели может решить состояние равновесия, который описывается массовое действие и уравнения баланса массы для каждого компонента.[2]

При нахождении состояния равновесия разработчик геохимических моделей решает вопрос о распределении масс всех видов, минералов и газов, которые могут быть образованы из основы. Это включает Мероприятия, коэффициент активности, и концентрация водных видов, насыщенность состояние минералов и летучесть газов. Минералы с индексом насыщения (log Q / K), равным нулю, считаются находящимися в равновесии с жидкостью. Те, у кого положительный индекс насыщения, называются перенасыщенный, указывая на то, что они склонны к осаждению из раствора. Минерал является недонасыщенным, если его индекс насыщения отрицательный, что указывает на то, что он способствует растворению.[8]

Разработчики геохимических моделей обычно создают модели пути реакции, чтобы понять, как системы реагируют на изменения в составе, температуре или давлении. Настраивая способ задания массы и теплопередачи (т.е. открытые или закрытые системы), модели можно использовать для представления различных геохимических процессов. Пути реакции могут предполагать химическое равновесие, или они могут включать законы кинетической скорости для расчета времени реакций. Чтобы предсказать пространственно-временное распределение химических реакций, происходящих вдоль пути потока, геохимические модели все чаще объединяются с гидрологические модели переноса массы и тепла на формирование реактивные транспортные модели.[2] Специализированные программы геохимического моделирования, разработанные как перекрестно-связываемые реентерабельные программные объекты, позволяют создавать модели реактивного переноса любой конфигурации потока.[9]

Типы реакций

Геохимические модели могут моделировать множество различных типов реакции. Среди них:

Простой фазовые диаграммы или графики обычно используются для иллюстрации таких геохимических реакций. Eh-pH (Pourbaix) диаграммы, например, являются особым типом диаграмм активности, которые графически представляют кислотно-щелочной и окислительно-восстановительный химические процессы.

Неопределенности геохимического моделирования

Различные источники могут способствовать ряду результатов моделирования. Диапазон результатов моделирования определяется как неопределенность модели. Одним из наиболее важных источников, которые невозможно определить количественно, является концептуальная модель, которая разрабатывается и определяется моделистом. Дополнительными источниками являются параметризация модели в отношении гидравлических (только при моделировании транспорта) и минералогических свойств.[10] Параметры, используемые для геохимического моделирования, также могут вносить вклад в неопределенность модели. Это прикладная термодинамическая база данных и параметры кинетического растворения минералов.[11] Различия в термодинамических данных (т. Е. Константах равновесия, параметрах температурной коррекции, уравнениях активности и коэффициентах) могут привести к большим погрешностям. Кроме того, большой диапазон экспериментально полученных констант скорости для законов скорости растворения минералов может вызвать большие различия в результатах моделирования. Несмотря на то, что это хорошо известно, неопределенности не часто учитываются при проведении геохимического моделирования.[12]

Уменьшение неопределенностей может быть достигнуто путем сравнения результатов моделирования с экспериментальными данными, хотя экспериментальные данные существуют не для всех условий температуры и давления и для каждой химической системы.[12] Хотя такое сравнение или калибровка не может быть проведена, следовательно, геохимические коды и термодинамические базы данных являются современными и наиболее полезными инструментами для прогнозирования геохимических процессов.

Программы общего пользования

В Сайт USGS предоставляет бесплатный доступ ко многим из перечисленных выше программ. [34]

Смотрите также

дальнейшее чтение

  • Appelo, C.A.J. и Д. Постма, 2005 г., Геохимия, подземные воды и загрязнение. Тейлор и Фрэнсис, 683 стр. ISBN  978-0415364287
  • Бетке К.М., 2008, Моделирование геохимических и биогеохимических реакций. Издательство Кембриджского университета, 547 стр. ISBN  978-0521875547
  • Меркель, Б.Дж., Б. Планер-Фридрих, Д.К. Nordstrom, 2008 г., Геохимия подземных вод: Практическое руководство по моделированию природных и загрязненных водных систем. Спрингер, 242 стр. ISBN  978-3540746676
  • Элькерс, Э. и Дж. Шотт (ред.), 2009 г., Термодинамика и кинетика взаимодействия воды и породы. Обзоры по минералогии и геохимии 70, 569 с. ISBN  978-0-939950-84-3
  • Чжу, К. и Г. Андерсон, 2002 г., Экологические приложения геохимического моделирования. Издательство Кембриджского университета, 300 стр. ISBN  978-0521005777

Рекомендации

  1. ^ Чжу К. и Дж. Андерсон, 2002 г., Экологические приложения геохимического моделирования. Издательство Кембриджского университета, 300 стр.
  2. ^ а б c d Бетке, C.M., 2008, Моделирование геохимических и биогеохимических реакций. Издательство Кембриджского университета, 547 стр.
  3. ^ Гаррелс, Р. и F.T. Mackenzie, 1967, Происхождение химического состава некоторых источников и озер. Концепции равновесия в природных водах, успехи в химии. 67, Американское химическое общество, Вашингтон, округ Колумбия, стр. 222-242.
  4. ^ Helgeson, H.C., 1968, Оценка необратимых реакций в геохимических процессах с участием минералов и водных растворов, I. Термодинамические соотношения. Geochemica et Cosmochimica Acta 32, 853-877
  5. ^ Helgeson, H.C., R.M. Гаррелс и Ф. Mackenzie, 1969, Оценка необратимых реакций в геохимических процессах с участием минералов и водных растворов, II. Приложения. Geochemica et Cosmochimica Acta 33, 455-481
  6. ^ Чжу, К., 2009 г., Геохимическое моделирование реакционных путей и геохимических реакционных сетей. В E.H. Элкерс и Дж. Шотт (ред.), 2009 г., Термодинамика и кинетика взаимодействия воды и породы. Обзоры по минералогии и геохимии 70, 533-569
  7. ^ Брэди, П. и К. Бетке, 2000 г., Помимо подхода Kd. Грунтовые воды 38, 321-322
  8. ^ а б Андерсон, Г. 2009 г., Термодинамика природных систем. Издательство Кембриджского университета, 664 стр.
  9. ^ "ChemPlugin.ORG". ООО «Водные Растворы». Архивировано из оригинал 11 июля 2013 г.. Получено 3 мая 2013.
  10. ^ Детлефсен, Франк; Хаазе, Кристоф; Эберт, Маркус; Дамке, Андреас (01.01.2011). «Влияние вариаций входных параметров на взаимодействие воды и минералов при моделировании секвестрации CO2». Энергетические процедуры. 10-я Международная конференция по технологиям контроля парниковых газов. 4: 3770–3777. Дои:10.1016 / j.egypro.2011.02.311.
  11. ^ Хаазе, Кристоф; Детлефсен, Франк; Эберт, Маркус; Дамке, Андреас (01.06.2013). «Неопределенность в геохимическом моделировании растворения CO2 и кальцита в растворах NaCl из-за различных кодов моделирования и термодинамических баз данных». Прикладная геохимия. 33: 306–317. Дои:10.1016 / j.apgeochem.2013.03.001.
  12. ^ а б Хаазе, Кристоф; Эберт, Маркус; Детлефсен, Франк (2016-04-01). «Неопределенности геохимических кодов и термодинамических баз данных для прогнозирования воздействия двуокиси углерода на геологические образования». Прикладная геохимия. 67: 81–92. Дои:10.1016 / j.apgeochem.2016.01.008.
  13. ^ Muller, B., 2004, CHEMEQL V3.0, Программа для расчета равновесия химического состава, титрования, растворения, осаждения, адсорбции, кинетики, диаграмм pX-pY, диаграмм растворимости. Лимнологический исследовательский центр EAWAG / ETH, Кастаниенбаум, Швейцария
  14. ^ ван дер Ли, Дж. и Л. Де Виндт, 2000, ШАХМАТЫ, другой компьютерный код видообразования и комплексообразования. Технический отчет № LHM / RD / 93/39, Ecole des Mines de Paris, Фонтенбло
  15. ^ Рид М.Х., 1982, Расчет многокомпонентных химических равновесий и процессов реакции в системах, включающих минералы, газы и водную фазу. Geochimica et Cosmochemica Acta 46, 513-528.
  16. ^ Стифел, К. и A.C. Lasaga, 1994, Совместная модель переноса множества химических веществ и кинетических реакций осаждения / растворения в применении к реактивному потоку в однофазных гидротермальных системах. Американский журнал науки 294, 529-592
  17. ^ Стифел К.И., 2001, GIMRT, версия 1.2: Программное обеспечение для моделирования многокомпонентного, многомерного реактивного транспорта, Руководство пользователя. Отчет UCRL-MA-143182, Ливерморская национальная лаборатория Лоуренса, Ливермор, Калифорния.
  18. ^ Wolery, T.J., 1992a, EQ3 / EQ6, программный пакет для геохимического моделирования водных систем, обзор пакета и руководство по установке (версия 7.0). Отчет Национальной лаборатории Лоуренса Ливермора UCRL-MA-110662 (1).
  19. ^ Parker, D.R., W.A. Norvell и R.L. Chaney, 1995, GEOCHEM-PC, программа химического анализа для IBM и совместимых персональных компьютеров. В: R.H. Loeppert, A.P. Schwab и S. Goldberg (ред.), Химическое равновесие и модели реакций. Специальная публикация Американского общества почвоведов 42, 253-269
  20. ^ Бетке, К.М., и С. Йикель, 2014 г., Руководства пользователя Geochemist's Workbench, Версия 10.0. ООО «Водные растворы», Шампейн, Иллинойс
  21. ^ Кулик, Д.А., 2002, Подход к минимизации энергии Гиббса для моделирования сорбционных равновесий на границе раздела минеральная вода: термодинамические соотношения для многоузлового поверхностного комплексообразования. Американский журнал науки 302, 227-279
  22. ^ Ченг, Х. и Г. Yeh, 1998 г., Разработка трехмерной модели подземного потока, теплопередачи и химического переноса в реакциях: 3DHYDROGEOCHEM. Журнал гидрологии загрязнителей 34, 47-83
  23. ^ Вестолл, Дж. К., Дж. Л. Закари и Ф.Ф.М. Morel, 1976, MINEQL, компьютерная программа для расчета химического равновесного состава водных систем. Техническая записка 18, R.M. Лаборатория Парсонса, Департамент гражданской и экологической инженерии, Массачусетский технологический институт, Кембридж, Массачусетс.
  24. ^ Шерер, У.Д. и Д.К. МакЭвой, 1994, MINEQL +, Программа химического равновесия для персональных компьютеров, руководство пользователя, версия 3.0. Программное обеспечение для экологических исследований, Inc., Hallowell, ME.
  25. ^ Эллисон, Д.Д., Д.С. Браун, К.Дж. Ново-Градац, 1991, MINTEQA2 / PRODEFA2, модель геохимической оценки для экологических систем, руководство пользователя версии 3.0. Отчет Агентства по охране окружающей среды США EPA / 600 / 3-91 / 021.
  26. ^ Parkhurst, D.L., 1995, Руководство пользователя PHREEQC, компьютерной модели для видообразования, пути реакции, адвективного переноса и обратных геохимических расчетов. Отчет по исследованию водных ресурсов Геологической службы США 95-4227.
  27. ^ Parkhurst, D.L. и C.A.J. Appelo, 1999, Руководство пользователя PHREEQC (версия 2), компьютерной программы для видообразования, периодической реакции, одномерного переноса и обратных геохимических расчетов. Отчет о исследованиях водных ресурсов Геологической службы США 99-4259.
  28. ^ Леал А.М.М., Кулик Д.А., Смит В.Р. и Саар М.О., 2017, Обзор вычислительных методов для химического равновесия и кинетических расчетов для моделирования геохимического и реактивного переноса. Чистая и прикладная химия. 89 (5), 145–166.
  29. ^ Perkins, E.H., 1992, Интеграция интенсивных диаграмм переменных и фазовых равновесий жидкости с SOLMINEQ.88 pc / shell. В Ю.К. Харака, А.С. Маэст (ред.), Взаимодействие воды и камня, Балкема, Роттердам, стр. 1079-1081.
  30. ^ Сюй Т., Э. Sonnenthal, N. Spycher и K. Pruess, 2004, руководство пользователя TOUGHREACT: программа моделирования неизотермического многофазного реактивного геохимического переноса в геологических средах с переменной насыщенностью. Отчет LBNL-55460, Национальная лаборатория Лоуренса Беркли, Беркли, Калифорния.
  31. ^ hem.bredband.net/b108693/-VisualMINTEQ_references.pdf
  32. ^ Болл, Дж. и Д.К. Nordstrom, 1991, Руководство пользователя WATEQ4F, с пересмотренной базой термодинамических данных и тестовыми примерами для расчета видообразования основных, следовых и окислительно-восстановительных элементов в природных водах. Отчет открытого файла Геологической службы США 91-183.
  33. ^ Типпинг Э., 1994, WHAM - модель химического равновесия и компьютерный код для вод, отложений и почв, включающий дискретную локальную / электростатическую модель связывания ионов гуминовыми веществами. Компьютеры и науки о Земле 20, 973-1023.
  34. ^ «Программное обеспечение для геохимии водных ресурсов». water.usgs.gov. Получено 2020-09-25.