Термодинамика глубинных геофизических сред
В. Паньков, В. Ульман, Р. Хайнрих, Д. Краке

1. Введение

С термодинамической точки зрения, Земля является тепловой машиной, характеризуемой множеством параметров, которые определяются из уравнений состояния и моделей кондепарованных сред. В разработке этих моделей достигнут большой прогресс [напр., Jeanloz, 1983; Hemley et al., 1985, 1987; Wall et al., 1986; Catti, 1986; Cohen, 1987; Dovesi et al., 1987; Wolf and Bukowinsky, 1987, 1988; Wall and Price, 1988; Matsui et al., 1987; Matsui, 1988, 1989; Price et al., 1989; Catlow and Price, 1990; Isaak et al., 1990; Reynard and Price, 1990; Agnon and Bukowinsky, 1990a; Matsui and Price, 1991; D'Arco et al., 1991; Walzer, 1992; Silvi et al., 1993; Catti et al., 1993; Boison and Gibbs, 1993]. Однако практические исследования в геофизике в большой степени основаны на использовании полуэмпирических уравнений состояния [Birch, 1952, 1986; O. Anderson, 1966b, 1995; Pankov and Ullmann, 1979; D. Anderson, 1967, 1987, 1989; Stacey, 1981; Bina and Helffrich, 1992; Wall et al., 1993]. Прямое определение свойств геоматериалов на основе экспериментальных данных при высоких давлениях и температурах необходимо для решения многих задач геофизики и дает важные ограничения на структуру уравнений состояния.

Со времени фундаментальной работы Берча [Birch, 1952] была накоплена огромная информация по свойствам геоматериалов и их геофизическим приложениям [напр., Stacey, 1977a, 1977b; Jeanloz and Thompson, 1983; Brown and Shankland, 1981; Жарков и Калинин, 1968; Zharkov, 1986; Jeanloz and Knittle, 1989; O. Anderson et al., 1992a, 1992b, 1993; Kuskov and Panferov, 1991; D. Anderson, 1989; O. Anderson, 1988, 1995]. Настоящая работа посвящена обзору и анализу связей между основными термодинамическими характеристиками геоматериалов и их изменения с давлением и температурой. Прежде всего мы имеем дело с восемью параметрами 2-го порядка. Внимание обращается на их самосогласованное определение, связь с уравнениями состояния, сводку подходов к определению уравнений состояния. Приводится пример самосогласованной базы данных по термодинамическим свойствам минералов. Далее каждый из параметров 2-го порядка рассматривается отдельно: устанавливаются тождественные соотношения между их P-T производными (параметрами 3-го порядка) и анализируются практически полезные приближенные соотношения, включая примеры явных P-T зависимостей параметров 2-го порядка. Приводятся оценки некоторых параметров 4-го порядка. Соотношения между различными параметрами представляются в виде, удобном для практического использования экспериментальных данных и теоретического анализа. Наконец, даны некоторые приложения результатов к оценкам свойств нижней мантии. Наш анализ служит дополнением к обзорам, представленным в работах [O. Anderson, 1995; Stacey, 1995].


2. Основные термодинамические соотношения

В классической термодинамике простые системы, испытывающие обратимые изменения состояния, описываются множеством параметров, включая гидростатическое давление P, температуру T, объм V или плотность r и энтропию S. Отправным пунктом термодинамического анализа являются стандартные выражения для полных дифференциалов термодинамических потенциалов [напр., Callen, 1960; Morse, 1969; Kelly, 1973]

eqn001.gif(1)

eqn002.gif(2)

eqn003.gif(3)

eqn004.gif(4)

где E - внутренняя энергия, F - свободная энергия, G - энергия Гиббса и H - энтальпия.

Для описания термодинамических свойств в геофизике используются восемь параметров второго порядка: коэффициент объемного теплового расширения a, изобарическая CP и изохорическая CV теплоемкости, изотермический KT и адиабатический KS модули сжатия, термический коэффициент давления t и адиабатический градиент температуры tS. Эти параметры определяются формулами

eqn005.gif(5)

eqn006.gif(6)
eqn007.gif

eqn008.gif(7)

eqn009.gif(8)

eqn010.gif(9)

eqn011.gif(10)

Используя равенство вторых смешанных производных термодинамических потенциалов, из (1)-(4) получаем четыре соотношения Максвелла [см. также Stacеy, 1977a]

eqn012.gif(11)

eqn013.gif(12)

eqn014.gif(13)

eqn015.gif(14)

Легко показать, что вторые производные термодинамических потенциалов выражаются через введенные выше параметры или характеристики

eqn016.gif

eqn017.gif(15)

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

eqn018.gif(16)

eqn019.gif(17)

Таким образом, если из эксперимента известны, как обычно, a, CP, KS (или KT ), то пять упомянутых выше соотношений можно записать в виде

eqn020.gif

eqn021.gif(18)

Приведем два примера использования термодинамических соотношений в геофизике. Первый из них касается уравнения Вильямсона-Адамса для градиента плотности внутри Земли [Birch, 1952; D. Anderson, 1989], которое мы представим здесь в виде

eqn022.gif(19)

где r и P плотность и давление в Земле, F = KS/r - сейсмический параметр, t = dT/dP - tS - сверхадиабатический градиент температуры в Земле. Уравнение (19) легко получается из

eqn023.gif

с помощью (5), (16) и (14) или (18).

Другой пример - адиабатический градиент температуры по глубине l в Земле (см. (14))

eqn024.gif(20)

где g - ускорение силы тяжести и использовано уравнение равновесия dP/dl = rg. Согласно модели Земли PREM [D. Anderson and Dziewonski, 1981\link], F = 50, 80 и 117 км 2 /c 2 на глубинах, соответственно, 400, 1071 и 2740 км. Полагая T = 1700 K, 2200 К и 3000 К последовательно на тех же глубинах [D. Anderson, 1989; Паньков, 1989] и g = 1-1.5, получаем (dT/dl)S = 0.3-0.5 K/км - величина, обычно цитируемая в геофизической литературе.


This document was generated by TeXWeb (Win32, v.1.0) on August 10, 1998.