Сравнительный анализ применения моделей Лихтенеккера и Лихтенеккера – Асаада для теоретической оценки эффективной теплопроводности горных пород

Вестник МГТУ. 2024. Т. 27, № 1. С. 103-112.

DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112

УДК 550.362

Сравнительный анализ применения моделей Лихтенеккера и Лихтенеккера - Асаада для теоретической оценки эффективной теплопроводности горных пород

С. А. Гусак

Горный институт КНЦ РАН, г. Апатиты, Россия; Филиал Мурманского арктического университета, г. Апатиты, Россия; e-mail: gusaksa@arcticsu.ru, ORCID: https://orcid.org/0000-0002-5028-326X

Реферат

Выполнен анализ условий и ограничений применения моделей Лихтенеккера и Лихтенеккера - Асаада для прогнозной оценки эффективной теплопроводности горных пород. Для сравнительной оценки рассмотрены экспериментальные данные по теплопроводности карбонатных пород (известняк, доломит) из бассейна Molasse на юго-западе Германии, импактной структуры Чиксулуб (полуостров Юкатан, Мексика) и баженовской свиты (Западная Сибирь, Россия). Сопоставительные оценки выполнялись на основе использования корреляционных соотношений, полученных авторами экспериментальных работ в результате регрессионного анализа петрофизических свойств карбонатных пород. Результаты расчетно-аналитической оценки показали, что наиболее высокая точность прогноза теплопроводности с использованием модели Лихтенеккера - Асаада может быть достигнута при наличии экспериментальных данных, обеспечивающих высокоточное определение структурного фактора. Изучены особенности влияния различных параметров модели Лихтенеккера - Асаада на погрешность оценки эффективной теплопроводности породы. Установлено, что на стадии предпроектных исследований для достаточно корректной теоретической оценки теплопроводности водонасыщенных слабопористых горных пород может быть использована модель Лихтенеккера.

Гусак С. А. Сравнительный анализ применения моделей Лихтенеккера и Лихтенеккера - Асаада для теоретической оценки эффективной теплопроводности горных пород. Вестник МГТУ. 2024. Т. 27, № 1. С. 103-112. DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112.

Comparative analysis of the application of the Lichtenecker and Lichtenecker - Asaad models for the theoretical assessment of the rocks effective thermal conductivity

Sergey A. Gusak

Mining Institute KSC RAS, Apatity, Russia; Apatity Branch of Murmansk Arctic University, Apatity, Russia; e-mail: gusaksa@arcticsu.ru, ORCID: https://orcid.org/0000-0002-5028-326X

Abstract

The analysis of the conditions and limitations of using the Lichtenecker and Lichtenecker -Asaad models for predictive assessment of the effective thermal conductivity of rocks have been carried out. For a comparative assessment, experimental data on the thermal conductivity of carbonate rocks (limestone, dolomite) from the South-West German Molasse basin, the Chicxulub impact structure (Yucatan Peninsula, Mexico) and the Bazhenov Formation (Western Siberia, Russia) have been considered. Comparative assessments were carried out based on the use of correlation relationships obtained by the authors of experimental works as a result of regression analysis of the petrophysical properties of carbonate rocks. The results of the computational and analytical assessment have shown that the highest accuracy of thermal conductivity prediction using the Lichtenecker - Asaad model can be achieved if experimental data are available that provide a highly accurate determination of the structure factor. The features of the influence of various parameters of the Lichtenecker - Asaad model on the error in estimating the effective thermal conductivity of the rock have been studied. It has been established that at the stage of pre-design research, the Lichtenecker model can be used for a fairly correct theoretical assessment of the thermal conductivity of water-saturated, slightly porous rocks.

Gusak, S. A. 2024. Comparative analysis of the application of the Lichtenecker and Lichtenecker -Asaad models for the theoretical assessment of the rocks effective thermal conductivity. Vestnik of MSTU, 27(1), pp. 103-112. (In Russ.) DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112..

Информация о статье Поступила в редакцию 29.01.2024;

принята к публикации 12.02.2024

Ключевые слова:

теплопроводность, горные породы, теоретическая оценка, модель, Лихтенеккер, Лихтенеккер - Асаад

Для цитирования

Article info

Received 29.01.2024;

accepted 12.02.2024

Key words:

thermal conductivity, rocks,

theoretical assessment, model, Lichtenecker, Lichtenecker - Asaad

For citation

Введение

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

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

Другой метод определения теплопроводности горных пород базируется на теоретической оценке эффективной теплопроводности пород с использованием различных моделей, связанных с определением теплопроводности на основе данных о содержании и составе породообразующих минералов и пористости, характеризующих многокомпонентную систему. Наиболее широкое применение получила модель среднего геометрического, которая часто называется моделью Лихтенеккера (Lichtenecker К.). В этой модели зависимость эффективной теплопроводности горной породы X от теплопроводности матрицы (минерального скелета) Xê, теплопроводности флюида в поровом пространстве Xf и пористости породы ф (отн. ед.) представлена в виде простого математического выражения

Ь = С  Ц.

Формула (1), которая применялась к широкому спектру смесей, включая биологические материалы, горные породы и минералы, а также композиционные материалы на основе пластмасс (Guo et al., 2017; Fuchs et al., 2015; Fuchs et al., 2013; Pechnig et al., 2007; Popov et al., 2003; Zhanga et al., 2018), долгое время считалась полуэмпирической по природе, поскольку, по мнению критиков, не имела какого-либо теоретического обоснования. В дальнейшем был выполнен ряд исследований, направленных на оценку и обоснование более прочной теоретической основы формулы Лихтенеккера (Goncharenko et al., 2000; Leäo et al., 2015; Simpkin, 2010; Zakri et al., 1998).

Однако, по мнению ряда российских и зарубежных специалистов, модель Лихтенеккера не может быть использована для достаточно корректной оценки теплопроводности горных пород, поскольку эта модель не учитывает текстурно-структурных особенностей неоднородных сред. Для учета этих особенностей была создана модифицированная теоретическая модель Лихтенеккера - Асаада, которая, по мнению различных специалистов, обеспечивает существенно более точные оценки и отличается от модели Лихтенеккера включением структурного фактора f следующим образом:

ь=^ • Ь ?.

В настоящей работе выполнена расчетно-аналитическая оценка, целью которой является анализ условий и ограничений применения моделей Лихтенеккера и Лихтенеккера - Асаада для прогноза теплопроводности горных пород. Для сравнительной оценки рассмотрены результаты экспериментальных исследований теплопроводности водонасыщенных карбонатных пород (известняк, доломит) из бассейна Molasse на юго-западе Германии, импактной структуры Чиксулуб (полуостров Юкатан, Мексика) и баженовской свиты (Западная Сибирь, Россия).

Экспериментальные исследования теплопроводности карбонатных пород

В работе немецких специалистов (Pechnig et al., 2007), посвященной оценке применимости теоретических моделей, было исследовано около 280 образцов горных пород из бассейна Molasse на юго-западе Германии. Керновые образцы изучались с помощью теплового и петрофизического сканирования с высоким разрешением по теплопроводности, плотности, пористости и скорости звука в породах в сухом и водонасыщенном состоянии. Результаты анализа минералогического состава пород позволили

Вестник МГТУ. 2024. Т. 27, № 1. С. 103-112.

DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112

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

Для учета эффектов геометрии пор и литификации осадочных пород (процесс превращения рыхлых осадков в твердые горные породы) в работе (Pechnig et al., 2007) использована модель Лихтенеккера - Асаада (формула (2)), на основе которой авторы выполнили численное моделирование для извлечения Xm и f из лабораторных данных. Результаты моделирования для теплопроводности матрицы и f-фактора приведены в таблице.

Таблица. Результаты моделирования методом Монте-Карло, определяющие параметры модели Лихтенеккера - Асаада (f-фактор (среднее ± стандартное отклонение) и теплопроводность матрицы

породы)

Тип породы Теплопроводность матрицы, Вт/(мК) f-фактор

Известняк 2,63 ± 0,13 0,77 ± 0,26

Доломит 3,71 ± 0,25 0,95 ± 0,13

В работе (Popov et al., 2011) приведены экспериментальные данные по теплопроводности доломитов и известняков импактной структуры Чиксулуб (полуостров Юкатан, Мексика). В этой работе описаны результаты измерений тепловых свойств кернов исследовательской скважины Yaxcopoil-1 (глубина 1 511 м), пробуренной в структуре Чиксулуб. На 451 сухом и водонасыщенном керне из интервала глубин 404-1 511 м были измерены теплопроводность, температуропроводность, объемная теплоемкость, коэффициент термической анизотропии, коэффициент термической неоднородности и, кроме того, пористость и плотность. Выполненные исследования показали, что осадочные породы, пересеченные скважиной Yaxcopoil-1, характеризуются значительными вариациями тепловых свойств, которые обусловлены изменчивостью пористости и литологии пород и демонстрируют разделение пласта на предударные и постударные участки по своим тепловым свойствам и пористости.

На основе корреляции между теплопроводностью и пористостью водонасыщенных пород авторами работы (Popov et al., 2011) получены следующие уравнения регрессии:

- для доломита:

Цф) = 4,87е-а018ф;

- для известняка:

Цф) = -0,028ф + 2,65.

На основе полученных уравнений регрессии авторами работы была оценена теплопроводность матрицы пород: 4,7 Вт/(мК) для доломитов и 2,61 Вт/(мК) для известняков.

В работе (Попов и др., 2017) приведены результаты теплофизических исследований образцов стандартного керна, целью которых, в частности, являлась оценка адекватности теоретических моделей теплопроводности реальным свойствам пород. На основе результатов изучения коллекции из 73 образцов карбонатных пород баженовской свиты (Западная Сибирь, Россия) получены уравнения регрессии, устанавливающие зависимость между теплопроводностью и пористостью для сухих, нефте-и водонасыщенных пород. При помощи этих уравнений выполнена оценка теплопроводности матрицы пород и подобраны значения структурного фактора f В частности, для водонасыщенных карбонатных пород баженовской свиты оцененное значение матрицы пород составило 3,1 Вт/(м-К), а подобранная величина фактора f - 0,51.

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

С целью оценки условий и ограничений применения модели Лихтенеккера (формула (1)) для прогнозных расчетов теплопроводности карбонатных пород в настоящей работе использованы результаты рассмотренных ранее экспериментальных исследований, по которым имелась достаточная информация для сравнительного анализа. Сопоставительные оценки выполнялись на основе корреляционных соотношений (в дальнейшем - экспериментальные данные), полученных авторами экспериментальных исследований в результате анализа петрофизических свойств карбонатных пород. При этом рассматривались водонасыщенные породы с теплопроводностью воды в поровом пространстве У = 0,6 Вт/(мК).

Так, с учетом данных таблицы, полученных авторами работы (Pechnig et al., 2007), формула Лихтенеккера - Асаада (2) для оценки зависимости теплопроводности пород У из бассейна Molasse от пористости ф и структурного фактора f может быть записана в следующем виде:

- для доломита:

У

- для известняка:

У

3,8

3,6

* 3,4

т

m

3,2 .

ь т с

1 3,0

о « 1

2,8 -

§ .

а 1

е

H 2,6

2,4 -

Результаты оценки теплопроводности пород по формулам (1), (5) и (6) приведены на рис. 1.

2,6 -

формула

доломит

0,04 0,08 0,12 0,16 Пористость, отн. ед.

«

2,5

т

m

ь,

Й 2,4

о н ч о в о

рп

по2,3 л п е

н

2,2

- формула

-

известняк

0,2

0,02 0,04 0,06 0,08 Пористость, отн. ед.

0,10 0,12

0

0

Рис. 1. Зависимость теплопроводности карбонатных пород бассейна Molasse от пористости, вычисленной

по модели Лихтенеккера

Аналогичные сравнительные оценки были выполнены на примере экспериментальных данных по теплопроводности доломитов и известняков импактной структуры Чиксулуб. Результаты сопоставительных расчетов по уравнениям регрессии (3 и 4) и с помощью модели Лихтенеккера при оцененной авторами работы (Popov et al., 2011) теплопроводности карбонатных пород приведены на рис. 2.

C учетом приведенных ранее значений теплопроводности матрицы пород и структурного фактора f полученных авторами работы (Попов и др., 2017) на основе регрессионного анализа образцов карбонатных пород баженовской свиты (Западная Сибирь, Россия), формула Лихтенеккера - Асаада для теоретической оценки теплопроводности таких пород в водонасыщенном состоянии может быть записана в следующем виде:

Цф) = 3,1

Результаты расчетов теплопроводности карбонатных пород баженовской свиты по формулам (1) и (7) приведены на рис. 3.

0 0,05 0,10 0,15 0,20 0,25 0 0,05 0,10 0,15 0,20 0,25 0,30 0,35 Пористость, отн. ед. Пористость, отн. ед.

Рис. 2. Зависимость теплопроводности карбонатных пород импактной структуры Чиксулуб от пористости, вычисленной по модели Лихтенеккера

Рис. 3. Зависимость теплопроводности карбонатных пород баженовской свиты от пористости, вычисленной по модели Лихтенеккера

Как видно из рис. 1-3, по мере увеличения пористости возрастают различия в теплопроводности пород, вычисленных по модели Лихтенеккера и экспериментальным данным. Применение формулы Лихтенеккера для сильнопористых пород может приводить к расхождению в расчетных и экспериментальных данных по теплопроводности пород более 14 %. При пористости менее 10 % это расхождение не превышает 8 %. При этом, как и в большинстве других исследований, модель Лихтенеккера прогнозирует более низкие значения эффективной теплопроводности горных пород. По результатам

исследований, выполненных на примере карбонатных пород баженовской свиты, в работе (Попов и др., 2017) делается вывод, что модель Лихтенеккера - Асаада дает существенно более точные оценки по сравнению с моделью Лихтенеккера.

В настоящей работе рассмотрено влияние различных параметров модели Лихтенеккера - Асаада на погрешность теоретической оценки эффективной теплопроводности породы.

Методика оценки

В работе (Попов и др., 2017) отмечается, что для успешного применения модели Лихтенеккера -Асаада необходимы предварительные экспериментальные исследования теплопроводности на представительной выборке изучаемого типа пород с применением прецизионного измерительного оборудования для определения фактора / соответствующего данному типу пород и конкретному поровому флюиду.

Последнее замечание имеет принципиальное значение для оценки условий и ограничений в применении модели Лихтенеккера. Из формул (1) и (2) со всей очевидностью следует, что различия в прогнозной оценке эффективной теплопроводности пород с помощью моделей Лихтенеккера и Лихтенеккера - Асаада определяются надежностью и точностью определения фактора /

Для анализа влияния различных параметров на погрешность определения фактора / в настоящей работе использован методический подход, предложенный в работах (Новиков, 2009; Столяров и др., 2007). Так, для оценки неопределенности 5/ в результатах определения фактора / в зависимости от относительной погрешности оценки теплопроводности минерального скелета (матрицы) породы 5Тт предложена следующая формула:

5/ = К  5Хт,

где

. т

Т7

К, = —г-Нт-,

1

Т - эффективная теплопроводность породы; Тт - теплопроводность матрицы породы; Т/ - теплопроводность порового флюида.

Согласно выражению (2), которое формализует модель Лихтенеккера - Асаада, неопределенность в оценке фактора / должна оказывать влияние на погрешность оценки эффективной теплопроводности породы 5Т. Для оценки взаимосвязи этих параметров использована зависимость, предложенная в работах

(Новиков, 2009; Столяров и др., 2007):

S— = К • Sf

где

v = Ф • in I —

Щ= ф •in

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

Из выражений (8) и (9) следует, что при заданной теплопроводности порового флюида (в частности, воды - для водонасыщенных пород), которая в несколько раз меньше теплопроводности матрицы и эффективной теплопроводности карбонатных пород, основным фактором, определяющим зависимость погрешности оценки фактора f от неопределенности в оценке —m, является соотношение — /—m.

На основе экспериментальных данных по эффективной теплопроводности рассматриваемых карбонатных пород, для которых были оценены теплопроводности матрицы пород (Попов и др., 2017; Pechnig et al., 2007; Popov et al., 2011), в настоящей работе с помощью формулы (9) выполнена оценка зависимости коэффициента К1 от соотношения —/—m. Общий характер такой зависимости иллюстрируют данные, приведенные на рис. 4, которые показывают резкое повышение погрешности в оценке фактора f с увеличением отношения —/—m, что для водонасыщенных горных пород является отражением обратной зависимости коэффициента К1 от пористости геоматериала. Так, например, при отношении — /—m для карбонатных пород баженовской свиты в диапазоне от 0,92 до 0,96 (уменьшение пористости от 10 до 5 %) значение коэффициента К увеличивается с 11,3 до 23,3. Это значит, что неопределенность в оценке теплопроводности матрицы на уровне 15-25 %, которая в большинстве реальных случаев является вполне удовлетворительным результатом (Столяров и др., 2007), может приводить к погрешности 170-280 % в оценке фактора f при пористости породы 10 % и 350-580 % при пористости 5 %.

Вестник МГТУ. 2024. Т. 27, № 1. С. 103-112.

DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112

25-

20-

15 —

10-

II ❖ о о доломит из бассейна Molasse ада доломит структуры Чиксулуб а а а известняк из бассейна Molasse + известняк структуры Чиксулуб • • • карбонатные породы баженовской свиты •

* А +

«

¿г а L Л.

/»if1 У

0,7

0,75

0,8

0,85

0,9

0,95

Рис. 4. Зависимость коэффициента К от отношения эффективной теплопроводности водонасыщенной породы к теплопроводности ее матрицы Fig. 4. Dependence of coefficient K on the ratio of the effective thermal conductivity of a water-saturated rock to the thermal conductivity of its matrix

Результаты расчетов коэффициента Kf в зависимости от соотношения X/Xm для карбонатных пород баженовской свиты и бассейна Molasse при оцененных значениях фактора f которые использованы в формулах (5-7), приведены на рис. 5.

' X/Jv

0,75 0,8

0,7

0,85

0,9

0,95

0-

г

•

о* о / карбонатные породы баженовской свиты . доломит из бассейна 9 9 9 Molasse известняк из бассейна Molasse

о о >

О

-0,05 -0,1 -0,15

£

-0,2 -0,25 -0,3-0,35

Рис. 5. Зависимость коэффициента Kf от отношения эффективной теплопроводности

водонасыщенной породы к теплопроводности ее матрицы Fig. 5. Dependence of the Kf coefficient on the ratio of the effective thermal conductivity of water-saturated rock to the thermal conductivity of its matrix

На рис. 5 видно, что в отличие от влияния неопределенности в оценке теплопроводности матрицы породы на погрешность оценки фактораf влияние структурного фактора f на неопределенность в оценке эффективной теплопроводности породы уменьшается при повышении соотношения X/Xm.

В дополнение к расчетным данным представляет интерес рассмотрение результатов теоретического анализа различных параметров модели Лихтенеккера - Асаада и, в частности, связи структурного фактора f с геометрическими характеристиками порового пространства (Новиков, 2009). В своем анализе автор этой работы базировался на теоретической модели эффективной теплопроводности многокомпонентной среды, основанной на теории эффективных сред, которая позволяет рассчитывать эффективное значение теплопроводности, исходя из данных о минеральном составе, пористости и геометрических характеристиках порового пространства (аспектное отношение и функция его распределения).

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

На основе использованной теоретической модели проведены расчеты эффективной теплопроводности пород с разным аспектным отношением пор/трещин, диапазон которого варьировался от 0,001 до 1 при вариации значений теплопроводности минерального скелета от 3 до 7 Вт/(м-К). Результаты расчетов показали, что влияние аспектного отношения пор/трещин на эффективную теплопроводность породы зависит от отношения теплопроводностей минерального скелета и флюида в поровом пространстве. При этом для водонасыщенных пород наблюдается минимальное влияние аспектного отношения пор/трещин. Для таких пород наибольшее влияние геометрических характеристик порового пространства на эффективную теплопроводность породы наблюдается при пористости от 10 до 40 %. При пористости менее 10 % расхождение в величине эффективной теплопроводности породы во всем диапазоне аспектного отношения пор/трещин не превышает 4 %. В этом случае теплопроводность матрицы выступает в качестве доминирующего фактора.

Заключение

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

Благодарности

Работа выполнена в рамках темы НИР № 0226-2019-0026 "Развитие методологии создания подземных комплексов для атомных станций малой мощности в Арктике".

Конфликт интересов

Автор заявляет об отсутствии конфликта интересов.

Библиографический список

Дучков А. Д., Соколова Л. С., Родякин С. В., Черныш П. С. Зависимость теплопроводности пород осадочного чехла Западно-Сибирской плиты от влажности и пористости // Геология и геофизика. 2014. Т. 55, № 5-6. С. 991-1000. DOI: http://dx.doi.org/10.15372/GiG20140520. EDN: SMXAZB. Зверев А. Б. Результаты натурных исследований устойчивости камерных сооружений подземной атомной станции // Использование подземного пространства страны для повышения безопасности ядерной энергетики : материалы междунар. конф., Апатиты, 20-22 октября 1992 г. В 3 ч. Ч. 2. Вопросы строительства и обеспечения безопасности подземных атомных станций. Апатиты : КНЦ РАН, 1995. С. 150-163.

Новиков С. В. Тепловые свойства терригенных коллекторов и насыщающих флюидов: дис. . канд. тех. наук. М., 2009. 155 с.

Попов Е. Ю., Ромушкевич Р. А., Попов Ю. А. Измерения тепловых свойств пород на стандартных образцах как необходимый этап теплофизических исследований месторождений углеводородов // Известия вузов. Геология и разведка. 2017. № 2. С. 56-70. EDN: YUINNB. Столяров М. М., Попов Ю. А., Тертычный В. В., Коробков Д. А. Особенности методики определения теплопроводности горных пород на основе теоретической модели Лихтенеккера - Асаада // Известия вузов. Геология и разведка. 2007. № 5. С. 69-72. EDN: LLZVFL. Fuchs S., Balling N., Förster A. Calculation of thermal conductivity, thermal diffusivity and specific heat capacity of sedimentary rocks using petrophysical well logs // Geophysical Journal International. 2015. Vol. 203, Iss. 3. P. 1977-2000. DOI: https://doi.org/10.1093/gji/ggv403.

BecTHHK MFry. 2024. T. 27, № 1. C. 103-112.

DOI: https://doi.org/10.21443/1560-9278-2024-27-1-103-112

Fuchs S., Schütz F., Förster H.-J., Förster A. Evaluation of common mixing models for calculating bulk thermal conductivity of sedimentary rocks: Correction charts and new conversion equations // Geothermics. 2013. Vol. 47. P. 40-52. DOI: https://doi.org/10.1016/j.geothermics.2013.02.002. Goncharenko A. V., Lozovski V. Z., Venger E. F. Lichtenecker's equation: Applicability and limitations // Optics Communications. 2000. Vol. 174, Iss. 1-4. P. 19-32. DOI: https://doi.org/10.1016/S0030-4018(99)00695-1.

Guo P. Y., Zhang N., He M. C., Bai B. H. Effect of water saturation and temperature in the range of 193 to 373 K on the thermal conductivity of sandstone // Tectonophysics. 2017. Vol. 699. P. 121-128. DOI: https://doi.org/10.1016/j.tecto.2017.01.024. Leäo T. P., Perfect E., Tyner J. S. Evaluation of Lichtenecker's mixing model for predicting effective permittivity of soils at 50 MHz // Transactions of the ASABE. 2015. Vol. 58, Iss. 1. P. 83-91. DOI: https://doi.org/10.13031/trans.58.10720. Pechnig R., Mottaghy D., Koch A., Jorand R. [et al.]. Prediction of thermal properties for mesozoic rocks of Southern Germany // Proceedings European Geothermal Congress 2007. Unterhaching, Germany, 30 May -1 June 2007. URL: https://pangea.stanford.edu/ERE/pdf/IGAstandard/EGC/2007/265.pdf. Popov Yu., Romushkevich R., Korobkov D., Mayr S. [et al.]. Thermal properties of rocks of the borehole Yaxcopoil-1 (Impact Crater Chicxulub, Mexico) // Geophysical Journal International. 2011. Vol. 184, Iss. 2. P. 729-745. DOI: https://doi.org/10.1111/j.1365-246X.2010.04839.x. Popov Yu., Tertychnyi V., Romushkevich R., Korobkov D. [et al.]. Interrelations between thermal conductivity and other physical properties of rocks: Experimental data // Pure and Applied Geophysics. 2003. Vol. 160. P. 1137-1161. DOI: https://doi.org/10.1007/pl00012565. Simpkin R. Derivation of Lichtenecker's logarithmic mixture formula from Maxwell's equations // IEEE Transactions on Microwave Theory and Techniques. 2010. Vol. 58, Iss. 3. P. 545-550. DOI: https://doi.org/10.1109/tmtt.2010.2040406. Zakri T., Laurent J.-P., Vauclin M. Theoretical evidence of Lichtenecker's mixture formulas based on the effective medium theory // Journal of Physics D: Applied Physics. 1998. Vol. 31, N 13. P. 1589-1594. DOI: https://doi.org/10.1088/0022-3727/31/13/013. Zhanga M., Bi J., Chen W., Zhang X. [et al.]. Evaluation of calculation models for the thermal conductivity of soils // International Communications in Heat and Mass Transfer. 2018. Vol. 94. P. 14-23. DOI: https://doi.org/10.1016/j.icheatmasstransfer.2018.02.005.

References

Duchkov, A. D., Sokolova, L. S., Rodyakin, S. V., Chernysh, P. S. 2014. Thermal conductivity of the sedimentary-cover rocks of the West Siberian Plate in relation to their humidity and porosity. Russian Geology and Geophysics, 55

Zverev, A. B. 1995. Results of in-situ studies of the stability of chamber structure in the underground nuclear power plant. Proceedings of Intern. conf. Utilization of countriy's underground space to increase nuclear power safety, October 20-22, 1992. Part 2. Apatity, pp. 150-163.

Popov, E. Yu., Romushkevich, R. A., Popov, Yu. A. 2017. Measurements of the rock thermal properties on the standard core plugs as a necessary stage of the thermal-physic investigations of the hydrocarbon fields.

Proceedings of Higher Educational Establishments. Geology and Exploration, 2, pp. 56-70. EDN: YUINNB.

Stolyarov, M. M., Popov, Yu. A., Tertuchny, V. V., Korobkov, D. A. 2007. Features of the method for determining the thermal conductivity of rocks based on the theoretical Lichtenecker - Asaad model.

Proceedings of higher educational establishments. Geology and Exploration, 5, pp. 69-72. EDN: LLZVFL.

Fuchs, S., Balling, N., Förster, A. 2015. Calculation of thermal conductivity, thermal diffusivity and specific heat capacity of sedimentary rocks using petrophysical well logs. Geophysical Journal International, 203(3), pp. 1977-2000. DOI: https://doi.org/10.1093/gji/ggv403. Fuchs, S., Schütz, F., Förster, H.-J., Förster, A. 2013. Evaluation of common mixing models for calculating bulk thermal conductivity of sedimentary rocks: Correction charts and new conversion equations. Geothermics, 47, pp. 40-52. DOI: https://doi.org/10.1016/j.geothermics.2013.02.002. Goncharenko, A. V., Lozovski, V. Z., Venger, E. F. 2000. Lichtenecker's equation: Applicability and limitations.

Optics Communications, 174(1-4), pp. 19-32. DOI: https://doi.org/10.1016/S0030-4018(99)00695-1. Guo, P. Y., Zhang, N., He, M. C., Bai, B. H. 2017. Effect of water saturation and temperature in the range of 193 to 373 K on the thermal conductivity of sandstone. Tectonophysics, 699, pp. 121-128. DOI: https://doi.org/ 10.1016/j.tecto.2017.01.024.

Leao, T. P., Perfect, E., Tyner, J. S. 2015. Evaluation of Lichtenecker's mixing model for predicting effective permittivity of soils at 50 MHz. Transactions of the ASABE, 58(1), pp. 83-91. DOI: https://doi.org/10.13031/ trans.58.10720.

Pechnig, R., Mottaghy, D., Koch, A., Jorand, R. et al. 2007. Prediction of thermal properties for mesozoic rocks of Southern Germany. Proceedings European Geothermal Congress 2007. Unterhaching, Germany, 30 May -1 June 2007. URL: https://pangea.stanford.edu/ERE/pdf/IGAstandard/EGC/2007/265.pdf.

Popov, Yu., Romushkevich, R., Korobkov, D., Mayr, S. et al. 2011. Thermal properties of rocks of the borehole Yaxcopoil-1 (Impact Crater Chicxulub, Mexico). Geophysical Journal International, 184(2), pp. 729-745. DOI: https://doi.org/10.1111/j.1365-246X.2010.04839.x.

Popov, Yu., Tertychnyi, V., Romushkevich, R., Korobkov, D. et al. 2003. Interrelations between thermal conductivity and other physical properties of rocks: Experimental data. Pure and Applied Geophysics, 160, pp. 1137-1161. DOI: https://doi.org/10.1007/pl00012565.

Simpkin, R. 2010. Derivation of Lichtenecker's logarithmic mixture formula from Maxwell's equations. IEEE Transactions on Microwave Theory and Techniques, 58(3), pp. 545-550. DOI: https://doi.org/10.1109/ tmtt.2010.2040406.

Zakri, T., Laurent, J.-P., Vauclin, M. 1998. Theoretical evidence of Lichtenecker's mixture formulas based on the effective medium theory. Journal of Physics D: Applied Physics, 31(13), pp. 1589-1594. DOI: https://doi.org/10.1088/0022-3727/31/13/013.

Zhanga, M., Bi, J., Chen, W., Zhang, X. et al. 2018. Evaluation of calculation models for the thermal conductivity of soils. International Communications in Heat and Mass Transfer, 94, pp. 14-23. DOI: https://doi.org/10.1016/j.icheatmasstransfer.2018.02.005.

Сведения об авторе

Гусак Сергей Андреевич - ул. Ферсмана, 24, г. Апатиты, Мурманская обл., Россия, 184209;

Горный институт КНЦ РАН, канд. техн. наук, доцент;

ул. Лесная, 29, г. Апатиты, Мурманская обл., Россия, 184209; филиал Мурманского арктического

университета в г. Апатиты;

e-mail: gusaksa@arcticsu.ru, ORCID: https://orcid.org/0000-0002-5028-326X

Sergey A. Gusak - 24 Fersmana Str., Apatity, Murmansk region, Russia, 184209;

Mining Institute KSC RAS, Cand. Sci. (Engineering), Associate Professor;

29 Lesnaya Str., Apatity, Murmansk region, Russia, 184209; Apatity Branch of Murmansk Arctic University;

e-mail: gusaksa@arcticsu.ru, ORCID: https://orcid.org/0000-0002-5028-326X