Научно-техническое обоснование рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки

Пашкевич, Роман Игнатьевич

Бесплатный автореферат и диссертация по наукам о земле на тему
Научно-техническое обоснование рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки
ВАК РФ 25.00.20, Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика

Автореферат диссертации по теме "Научно-техническое обоснование рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки"

На правах рукописи

Пашкевич Роман Игнатьевич

□03400258

НАУЧНО-ТЕХНИЧЕСКОЕ ОБОСНОВАНИЕ РАЦИОНАЛЬНЫХ ПАРАМЕТРОВ ТЕПЛОПЕРЕНОСА И ФИЛЬТРАЦИИ ДВУХФАЗНОГО ТЕПЛОНОСИТЕЛЯ ПРИ ОСВОЕНИИ ГЕОТЕРМАЛЬНЫХ МЕСТОРОЖДЕНИЙ КАМЧАТКИ

Специальность 25.00.20 - Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика

2 2 ОКТ км

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук

Хабаровск - 2009

003480258

Работа выполнена в Научно-исследовательском геотехнологическом центре Дальневосточного отделения Российской академии наук (НИГТЦ ДВО РАН)

Научный консультант:

доктор геолого-минералогических наук Трухин Юрий Петрович (НИГТЦ ДВО РАН, Петропавловск-Камчатский)

Официальные оппоненты:

доктор физико-математических наук, профессор

Попов Юрий Анатольевич (МГГРУ, Москва)

доктор технических наук

Галушкин Юрий Иванович (МГУ, Москва)

доктор технических наук

Курилко Александр Сардокович (ИГДС ДВО РАН, Якутск)

Ведущая организация: Институт физики Земли РАН, Москва

Защита состоится "05" нолЛ^Я 2009 г. в ч. 00 мин. на заседании диссертационного совета ДМ005.009.01 при Институте горного дела ДВО РАН по адресу: 680000, г. Хабаровск, ул. Тургенева, 51. тел./факс: (4212) 32-79-27.

С диссертацией можно ознакомиться в библиотеке Института горного дела ДВО РАН.

Автореферат разослан " СК-Т^Л 2009 г.

Ученый секретарь диссертационного совета, к.т.н., доцент

Корнеева С.И.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Многолетние теоретические и натурные исследования термогидродинамического и гидрогеохимического режима Мутновской магматогенной геотермальной системы выполнены Е.А. Вакиным, И.Т. Кирсановым, Б.Г. Поляком, Ю.П. Трухиным, P.A. Шуваловым, Ю.А. Тараном, В.П. Пилипенко, труды которых внесли значительный вклад в решение проблемы формирования магма-тогенных геотермальных систем. Численные модели теплопереноса в океанической и континентальной литосфере разработаны Ю.И. Галушкиным. Численные термогидродинамические модели гидротермальных месторождений Камчатки и Курильских островов разработаны A.B. Кирюхиным.

Методика расчета параметров пароводяных потоков в стволе эксплуатационных геотермальных скважин развивалась в работах В.А. Дрознина, В.Д. Бе-лодеда и А.Н. Шулюпина. Методы измерения теплофизических свойств горных пород при высоких температурах и давлениях разработаны Ю.А. Поповым. Вопросы технологии скважинных систем извлечения геотермальной энергии при докритических термодинамических условиях рассмотрены И.Т. Аладьевым, H.H. Ароновой, В.Л. Артемьевой, Э.И. Богуславским, В.А. Васильевым, Л.Ф. Вознюком, С.Г. Гендлером, Ю.Д. Дядькиным, А.Г. Егоровым, Г.Н. Забарным, Г.Н. Кононенко, O.A. Кремневым, Э.И. Мерзляковым, Ю.П. Морозовым, И.А. Павловым, Ю.М. Парийскнм, Т.Ю. Пискачевой, A.M. Пудовкиным, В.А. Романовым, И.А. Рыженко, А.Н. Саламатиным, H.H. Смирновой, В.Н. Трусовым, A.C. Цырульниковым, А.Н. Щербанем и A.B. Шурчковым.

В настоящее время для выработки энергии используется теплоноситель в виде пара, воды или их смеси, находящийся в естественных условиях продук-

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

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

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

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

Задачи исследований:

- обобщение и анализ современного состояния проблемы моделирования и управления термогидродинамическими процессами в эксплуатационных скважинах и породах геотермальных месторождений и магмагеотермальных систем;

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

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

- установление характера изменения и скорости волны разрежения в пористой среде в процессе ее быстрой декомпрессии;

- установление структуры профилей давления, скорости и плотности газа и их динамики в процессе быстрой декомпрессии пористой среды;

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

- разработка численных моделей Мутновского месторождения парогидро-тсрм, а также Мутновской, Авачинской и Кудрявой магмагеотермальных систем;

- разработка метода оценки глубины очага магмагеотермалыюй системы по геотермическим данным;

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

Объект исследовании: процессы теплопереноса и фильтрации в горных породах геотермальных месторождений Камчатки.

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

Методы исследований. При выполнении диссертационной работы использовался комплекс методов исследования:

- обобщение и анализ ранее выполненных исследований и практики эксплуатации геотермальных промыслов;

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

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

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

- математическое и экономико-математическое моделирование процессов теплопереноса в породных массивах;

- методы вычислительной математики при разработке численной модели фильтрации газа при быстрой декомпрессии пористой среды;

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

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

Основные научные положения, защищаемые автором.

1. Теоретические основы моделирования теплопереноса при фильтрации двухфазного теплоносителя в виде пароводяной смеси в продуктивных зонах геотермальных месторождений должны базироваться на методе объемного осреднения, дополненном сформулированными в работе новыми теоремами ос-

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

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

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

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

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

Научная новизна работы заключается в том, что:

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

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

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

- обоснованы рациональные способы управления слабогазирующими эксплуатационными скважинами геотермальных месторождений;

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

- разработан ряд новых численных термогидродинамических моделей Мутновского месторождения парогидротерм, Мутновской, Авачинской и Кудрявой магмагеотермальных систем, отличающихся от разработанных ранее учетом в явном виде и анализом фазовых переходов теплоносителя в полном природном диапазоне от надкритического флюида до перегретой жидкости;

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

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

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

Достоверность и обоснованность научных положений и выводов обеспечивается:

- удовлетворительной сходимостью со средним расхождением 9% результатов численного моделирования процесса быстрой декомпрессии газонасыщенной пористой среды с полученными в представляемой работе данными лабораторных экспериментов;

- представительным объемом исходной информации по 63-м скважинам Мутновского месторождения парогидротерм (термометрия и точечные замеры высокотемпературным глубинным термометром);

- достаточной точностью калибровки модели естественного состояния Мутновского месторождения парогидротерм: среднее отклонение от фактического распределения температур 11% при максимальном 17%;

- сходимостью результатов модели эксплуатации Мутновского месторождения парогидротерм с фактическими данными мониторинга: среднее отклонение 11% по температурам добычных и наблюдательных скважин и 16% по энтальпиям на забое добычных скважин;

- удовлетворительной сходимостью со средним расхождением 7,5% результатов моделирования теплопереноса в Мутновской магмагеотермальной

системе с фактическими данными замеров и оценок температуры и дебита фу-марольных газов;

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

Практическое значение работы состоит в том, что использование полученных результатов исследований обеспечило возможность:

- выполнить оценку естественных ресурсов Мутновского месторождения парогидротерм;

- увеличить объем добычи термоминеральных вод Мапкинского геотермального месторождения, снизив при этом эксплуатационные затраты и выполнив требование рационального отбора теплоносителя в соответствии с текущей потребностью;

- установить рациональные параметры перспективных геотермальных циркуляционных систем, функционирующих в начальных надкритических термодинамических условиях и создаваемых с целью использования тепловых ресурсов очагов магмагеотермальных систем, которые целесообразно применять при проектировании;

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

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

Реализация результатов работы.

Методические разработки в области численного моделирования процесса теплопереноса использованы ООО "Аква" в практике проектирования разработки месторождений термальных вод, а также при оценке запасов Мутновского месторождения парогидротерм, имеется акт внедрения.

Способы управления эксплуатационными геотермальными скважинами и способ автоматизации геотермального промысла внедрены в производственную деятельность ООО "Аквариус", имеется акт внедрения.

Результаты численного моделирования магмагеотермальных систем и режима эксплуатации геотермальных циркуляционных систем в надкритических условиях использованы в исследованиях теплопереноса в породах геотермальных месторождений Камчатки и Курильских островов.

Результаты экспериментальных и численных исследований процесса быстрой декомпрессии газонасыщенной пористой среды использованы в учебном процессе в КамчатГТУ (1997-2003) и в филиале ДВГТУ (ДВПИ им. В.В. Куйбышева) в г. Петропавловске-Камчатском (2001-2007) при преподавании дисциплин "Гидрогазодинамика", "Гидравлика", "Основы технологии и проектирования ТЭК".

Личный вклад автора состоит:

- в постановке цели, формулировании задач и разработке методики исследований;

- в систематизации и обработке исходных данных;

- в формулировании и доказательстве теорем осреднения по межфазной поверхности;

- в разработке и численной реализации математических моделей фильтрации и теплопереноса;

- в анализе результатов лабораторных и численных экспериментов и установлении динамических характеристик исследованных процессов фильтрации и теплопереноса;

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

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

Представляемая работа основана на результатах исследований, выполненных в 1994-2008 гг. при непосредственном участии и под руководством автора, в том числе, в 2004-2008 гг. - в рамках плановых тем НИР в НИГТЦ ДВО РАН: "Исследование высокотемпературных (Т>400°С) магматогенных геотермальных систем (современных формирующихся локальных рудЕю-магматических центров) и гидротермальных месторождений Камчатки и Курильских островов как нетрадиционных возобновляемых источников минеральных и энергетических ресурсов" (2002-2005 гг.) и "Энергосберегающие процессы и геотехнологии извлечения ценных компонентов из гидротермальных растворов геотермальных систем" (2006-2008 гг., № ГР 01.2.006 08893).

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

- Международной конференции "Тепловое поле Земли и методы его изучения" (Москва, 2008);

- расширенном семинаре лабораторий прикладной геофизики и вулканологии и сейсмотектоники Института физики Земли РАН (Москва, 2009);

- семинаре кафедры нефтегазовой и подземной гидромеханики Российского государственного университета нефти и газа имени М.И. Губкина (Москва, 2009);

- Международной научно-практической конференции "Эффективные энергетические системы и новые технологии" (Казань, 2001);

- Международном полевом Курило-Камчатском семинаре "Геотермальные и минеральные ресурсы областей современного вулканизма" (Петропавловск-Камчатский, 2005);

- Международном семинаре "Проект научного бурения на Мутновском вулкане. Исследование связи магматической и гидротермальной систем" (Петропавловск-Камчатский, 2006);

- региональной научно-практической конференции "Минерально-сырьевые ресурсы как фактор развития промышленной и социальной инфраструктуры Камчатского края" (Петропавловск-Камчатский, 2007);

- научно-технических советах ОАО "Геотерм" (Петропавловск-Камчатский, 2006-2007);

- Ученых советах НИГТЦ ДВО РАН (2006-2008);

- семинарах лаборатории геохимии и геотехнологии и лаборатории физико-химической гидродинамики НИГТЦ ДВО РАН (2004-2008);

- семинарах кафедры естественнонаучных и технических дисциплин филиала ДВГТУ (ДВПИ им. В.В. Куйбышева) в г. Петропавловске-Камчатском (20012007);

- конференциях профессорско-преподавательского состава КамчатГТУ (1994-2003).

Публикации. По теме диссертации опубликовано 40 печатных работ, в том числе 17 - в журналах и изданиях "Перечня ведущих рецензируемых научных журналов и изданий" Высшей аттестационной комиссии Минобрнауки России, из них 15 - в рекомендованных экспертным советом ВАК по наукам о Земле.

Объем и структура работы. Диссертация состоит из введения, семи глав, заключения, списка литературы из 231 наименования, иллюстрирована 138 рисунками, содержит 29 таблиц, общий объем - 303 страницы.

Автор выражает глубокую благодарность научному консультанту Ю.П. Трухину и Г.В. Секисову за научно-методическое консультирование; благодарит сотрудников НИГТЦ ДВО РАН: A.C. Латкина, В.А. Степанова и P.A. Шувалова за ценные советы и поддержку в работе; Т.П. Белову и В.В. Потапова -за организационную помощь. За сотрудничество при проведении проектных и опытных работ автор благодарен сотрудникам производственных организаций:

A.A. Демченко, С.М. Каратаеву, М.А. Петрову (ООО "Аквариус"), С.Ю. Киселеву, К.Н. Мальцевой (ООО "Аква"), И.И. Черневу (ОАО "Геотерм"). За сотрудничество в проведении лабораторных и численных экспериментов автор благодарит В.К. Панова, В.В. Таскина и A.B. Шадрина. Материалы диссертации обсуждались с М.А. Алидибировым, О.Л. Басмановым, Э.И. Богуславским,

B.Я. Вороновицким, С.Г. Гендлером, В.Г. Греем, Н.С. Жатнуевым, Г.Н. Забар-ным, В.В. Кадетом, A.B. Кирюхиным, М.Н. Кравченко, JI.K. Москалевым, A.A. Никитиным, Б.Е. Паршиным, Б.Г. Поляком, О.Б. Селянгиным, М.Д. Сидоровым, А.Л. Собисевичем, В.М. Сугробовым, И.С. Уткиным, С.А. Федотовым,

C.М. Хассанизаде, И.И. Черневым, P.A. Шуваловым, Ю.В. Шуваловым и А.Н.

Шулюлииым, которых автор благодарит за ценные советы и критические замечания. Автор благодарен за консультации по программным комплексам: НУ-Б110ТНЕ1Ш - сотрудникам национального центра Геологической службы США К. Киппу, П.А. Сиху, С.Е. Ингебритсену, Д.О. Хайбе, а также О. Абрамову из университета штата Аризона; Т01ТСН2 и \YinGridder - сотрудникам национальной лаборатории им. Лоуренса в Беркли (США) К. Пруссу и Л. Пану.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

В главе 4 обоснованы новые методы и средства управления процессами гидрогазодинамики и теплопереноса в стволе добычных геотермальных скважин, а также способа автоматизации геотермального промысла на примере Малкинского месторождения термоминеральных вод.

В главе 5 представлены результаты численного термогидродинамического моделирования Мутновского месторождения парогидротерм, разработаны методические рекомендации по выполнению численного моделирования геотермальных месторождений трещинно-жилыюго типа.

В главе 6 представлены результаты численного моделирования магмагеотермальных систем Камчатки и Курильских островов (вулканы Мутновский, Авачинский, Кудрявый), разработаны методические рекомендации по моделированию магмагеотермальных систем, предложен метод оценки глубины залегания очага Мутповской системы по геотермическим данным.

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

Основные результаты исследований отражены в следующих защищаемых положениях.

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

Используя подход Слеттери, на основе теоремы переноса Рейнольдса для движущейся материальной поверхности, для произвольных вектора а и скаляра <р (например, скорости и температуры), определенных на межфазной поверхности "теплоноситель - горная порода", сформулированы теоремы, выражающие градиент осредненного по межфазной поверхности потока вектора и градиент осредненного скаляра через изменения величин на микроуровне, в масштабе пор или трещин:

Щаща)ца = — £[траа] сИ + — £([тР(Та]у) щей1 - Ч1пАра(апр„)р„, (1)

Сре Сре

v<<p)^^ет = — | <Р*Ре М + — | (Р^реЧ) Прей1 - ШАра(<р)ро. (2)

Сре Сре

где средняя величина по межфазной поверхности Ар„ определена в виде:

(ЬЬ«=т- (3)

Ара

Сре- контур, ограничивающий поверхность /1^ ; /?=/ - жидкость; ¡¡=у - пар; а -твердая фаза (скелет породы); Я - радиус объема осреднения V, рис. 1.

Доказаны теоремы, связывающие средние по межфазной поверхности производные по времени от величин на микроуровне, с производными по времени от средних по межфазной поверхности: д да 1 С. , дЫАпа ■¿¡¡(апр^ро = (^Пра)ра + ((Уа)^рапра)ра + - Ф^0а\апс11--(апра)ра, (4)

Сре

д д<р 1 Г дЫАра

= - №Кмра^раПра)ра + Ф <Р^раУ0аМ--^—(<р)ра, (5)

с0е

где у/ра и КМра - скорость и средняя кривизна межфазной поверхности, п -единичный вектор, в зависимости от направления вектора скорости межфазной границы, сонаправленный или противоположный вектору N.

Аналогичные теоремы для пространственных и временных производных доказаны для случая переменного объема осреднения, когда наблюдаются существенные изменения свойств среды в макромасштабе.

Рис. 1. Объем осреднения Vс центром О' на произвольной кривой 5 в многофазной системе. О - центр макросистемы коорди-5 нат, г0 - радиус-вектор центра объема осреднения, г - радиус-вектор произвольной точки в макросистеме, р - радиус-вектор в системе, связанной с центром объема осреднения, А - поверхность объема осреднения, Лре - поверхность входа и выхода фазы /3 в V, /1/к, и >1/;„ - нормали к поверхностям А^ и Ар„, Т/к, и Т/!„ - единичные касательные векторы к контурам, ограничивающим поверхности Аре и Ар„ и ориентированные с Нрс и Пр„ по правилу правого винта, и \>ра - единичные векторы, нормальные к контурам, ограничивающим поверхности Аре и Ар„ и касательные к ним, N - нормаль к Л.

Из теорем (1) и (2) для случая двухфазного (жидкость-пар) потока следует: - | 2КтущуйА = гфУви (6)

где г=2/й, КШу - средняя кривизна межфазовой поверхности жидкость-пар, п(„ - единичный вектор нормали к поверхности жидкость-пар, ф - пористость, я/ - объемная насыщенность жидкостью.

Формула (6) является дополнением к теории объемного осреднения для двухфазного потока и позволяет выразить капиллярные силы в уравнениях сохранения импульса через макроскопический параметр - градиент насыщенности ЖИДКОСТЬЮ Я/.

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

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

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

тI ~ КкГ1

£, 1 Г

V, - —- - vv) - 7£,Р! - - I РуП^А + агфЧБ1 + = О,

АIV

где £(=ф5(, К- - тензор абсолютной проницаемости, уь vv - скорости жидкости и пара, р1 - давление жидкой фазы, р'у - микроскопическое давление пара, А1у -межфазная поверхность жидкость-пар, кг1 - относительная проницаемость, к1у -параметр сопротивления на межфазной границе, а - коэффициент поверхностного натяжения.

В медленном потоке силой трения на границе жидкость-пар можно пренебречь и уравнение (7) упрощается:

=--- [Уод + - | р>1исМ - агфУз1 - (8)

¿IV

В квазиравновесном случае уравнение для смеси имеет вид:

ЕуЦу £1{11

"V ~~ ~ 7£<Р< + + (^Ру + £1Рдё = о.

(9) (10)

Гидростатическое уравнение для пароводяной смеси:

Уф( 1 - + = <Я(1 - + ЗД^ + 5,.

При использовании стандартного определения капиллярного давления в виде рс = ру — р, уравнение гидростатики в случае постоянной пористости и постоянного давления газообразной фазы имеет вид:

ЩРс = -[(1 - 5,)ру + - (11)

Для умеренной плотности газа:

Щ =Р1ё + [аг + - р,)] —■ (12)

Для модели пористой среды в виде гексагональной укладки сфер диаметром В из уравнения (12) следует:

• в случае постоянного радиуса кривизны поверхности жидкость-газ

Рс/Рстак = ехр[-(51 - з^.^р&фО2/0,\Ъ$па\. (13)

• в случае, когда радиус кривизны линейно растет с высотой над уровнем нулевого капиллярного давления

рцфо2/0,155шт

Рс/Рс тах = [(1 - *,)/(! - *|т,Л

(14)

Расчеты по формулам (13) и (14) удовлетворительно согласуются с пересчитанными в осях рс/рс тех ~ Ъ экспериментальными данными, рис. 2.

Рс/Рс тах

1,00 0.75 0,50 0,25

•<

) 2

1/ 4/

\ 1

Рис. 2. Зависимость относительного капиллярного давления от насыщенности. 1,2- формула (13) соответственно при й/З и 1>/2; 3 - формула (14); 4 - экспериментальные данные [Реггагк! <?/ а/., 1990].

0,4

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

Лабораторные эксперименты по изучению процесса фильтрации пористой среды при ее быстрой декомпрессии проводились автором, совместно с В.К. Пановым, на вертикальной установке типа "ударная труба" квадратного сечения 40x40 мм, длиной 2060 мм. Камера высокого давления (КВД) заполнялась пористой средой и отделялась от камеры низкого давления (КНД) целлофановой диафрагмой, разрушаемой устройством прожига. Выбранные длины КВД (1025 мм) и КНД (1035 мм) обеспечивали необходимую длительность невозмущенного процесса истечения газа до отражения волн разрежения и давления от торцов трубы. В качестве модельной пористой среды использовался перли-то-пластобетон (коммерческое наименование "пластиприн") - материал на основе фенолформальдегидной смолы с добавлением вспененного перлита. Измеренная активная пористость материала составляла 0,8-0,95; средняя измеренная при градиентах давления до 32,5 кПа/м проницаемость - 6,67-10"" м2. Измерение давления производилось пьезодатчиками диаметром 3 мм на основе ЦТС-19, установленными заподлицо в боковой стенке КВД, а также в ее глухом торце. Сигналы датчиков регистрировались цифровыми запоминающими осциллографами С9-8. Начало отсчета времени связывалось с моментом разрыва диафрагмы, определяемом по сигналу датчика, установленного в КВД на расстоянии 23 мм от диафрагмы. Для корректировки искажения сигналов датчика вследствие стекания электрического заряда в ходе длительного, до 10 мс, процесса фильтрации использовалась составленная автором программа численного интегрирования второго закона Кирхгофа для цепи "датчик - повторитель сигнала". Постоянная времени измерительной цепи устанавливалась в специальных экспериментах с КВД, заполненной воздухом, путем сравнения сигналов датчиков с результатами численного интегрирования уравнений Эйлера по составленной автором программе на языке ФОРТРАН.

Установлено существование двух типов волн, сопровождающих процесс: волны упругой деформации, распространяющейся по твердому скелету пористой среды, и волны разрежения, распространяющейся по газу в порах. Скорость деформационной волны составляла 410±18 м/с. При этом темп падения давления в пористой среде в диапазоне избыточного начального давления газа в КВД от 1 до 2 бар составлял 84-21 бар/с. Инструментальная погрешность определения скорости волн не превышала 2%. Разброс, превышающий эту величину,

объяснялся неоднородностью материала пористой среды и ее неустранимыми продольными колебаниями в процессе фильтрации.

Для теоретического изучения процесса декомпрессии пористой среды на основе метода объемного осреднения разработана математическая модель фильтрации и теплопереноса в консервативной форме, включающая, по сравнению со стандартной, инерционные члены в уравнении импульса и энергии:

|р+|(риО = 0, (15)

д,.д цт

— + = (16)

е+"2

+ YW =°> С17)

где р - плотность газа, w - действительная скорость фильтрации, р - давление, р - динамическая вязкость, е - внутренняя энергия газа, ф и К- соответственно пористость и проницаемость среды, t - время, х - координата, направленная вдоль потока газа слева направо. Начальные условия имеют вид: н{х,0)=0, 0<х<1; p(x,0)=pi, р(х,0)=р;, ()<х<ха; р(х,0)=р2, р(х,0)=р2, x0<x<L, (18) где .v(l - расстояние от левого торца КВД до диафрагмы, L - общая длина трубы, индексы 1 и 2 относятся к значениям соответственно в КВД и КНД.

Система уравнений (15) - (17) решалась численно с помощью явного метода Мак-Кормака по составленной автором программе на языке ФОРТРАН.

Рассматривалось три случая расположения начального разрыва давления: а) на поверхности торца пористой среды, что соответствует условиям экспериментов; б) в газе, что соответствует случаю неполного заполнения КВД пористой средой; в) внутри пористой среды. На рис. 3-5 приведены расчетные профили давления фильтрующегося газа. В результате анализа данных численных экспериментов сделаны следующие выводы.

Распад разрыва на границе пористой среды. Ударной волны не образуется, поток газа запирается высоким фильтрационным сопротивлением пористой среды, на границе давление газа практически равно начальному. На границе "пористая среда - газ" формируется контактный разрыв, движущийся в сторону камеры низкого давления со скоростью, меньшей скорости волны разрежения в пористой среде. Скорость контактного разрыва растет с ростом проницаемости среды. Скорость газа на границе равна критическому значению, в расчетном диапазоне - порядка 60-150 м/с. Критическая скорость газа на границе растет с ростом начального давления в КВД и ростом проницаемости. От границы пористой среды в сторону КНД формируется область непрерывного течения. Волна разрежения имеет более крутой фронт, чем волна разрежения по чистому воздуху в случае отсутствия пористой среды. Скорость волны разрежения по газу в пористой среде ниже скорости звука, с которой распространяется волна разрежения по чистому газу, растет с ростом коэффициента проницаемости и при расчетном значении 6,67-10"11 м2 составляла 142,5 м/с. Скорость контактного разрыва растет с ростом проницаемости.

лл4

ЗА 4/

5/ 6/

0,36 0,43 0,50 0,57 0,64

Рис. 3. Эпюры безразмерного давления в случае расположения начального разрыва на границе пористой среды. 1 - 0,5; 2 - 1,0; 3 - 1,5; 4 - 2,0; 5 - 2,5;

6 - 3,0 мс с начала распада разрыва.

7 - положение начального разрыва.

Х'=Х/Т, р'=р/р2-

1,8 1,6

1,4 1,2 1

п Л 1 \ /

/ ■■■■■>■ /....... / / /Л,

0,29 0,36 0,43 0,50 0,57 0,64 .г'

Рис. 4. Эпюры безразмерного давления при неполном заполнении КВД. 1 -1; 2-2; 3-3; 4-4; 5 -5; 6-6 мс с начала распада разрыва.

7 - граница "пористая среда- газ".

8 - положение начального разрыва.

у *

У'^Ч

0.43

0,46 0,50

0.54

Рис. 5. Эпюры безразмерного давления при распаде разрыва внутри пористой среды. 1 - 1; 2 - 2; 3 - 3; 4 - 4; 5 - 5; 6 - 6 мс с начала распада разрыва.

7 - положение начального разрыва.

8 - граница "пористая среда- газ".

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

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

Распад разрыва внутри пористой среды. На границе первоначального разрыва критических давления и скорости не достигается, процесс идет плавно. Контактный разрыв формируется на границе начального разрыва, двигаясь в область КНД. Задний фронт разрыва пологий, в отличие от случая разрыва в чистом воздухе. Скорость максимальна в момент разрыва диафрагмы, ее эпюра симметрична относительно положения начального разрыва. С течением времени скорость падает, от границы пористой среды в область КНД формируется зона постоянной скорости газа.

Сравнение расчетного и экспериментального темпа падения давления газа. Теоретический фронт волны разрежения на временах меньше 0,8 мс с начала процесса практически совпадает с экспериментальным, рис. 6. На больших временах темп падения давления становится выше измеренного. С момента времени 3 мс темпы падения давления в численных и лабораторных экспериментах практически равны. Максимальное отклонение значений давления приходится на время около 3 мс и составляет 9%. При малых временах отклонение уменьшается: 0,8 мс - 2,3%, 1,6 мс - 4,4%. При больших временах отклонение также меньше максимального, для 11,2 мс составляет 6,7%. Более резкий фронт волны разрежения и быстрый спад давления, по сравнению с полученными в лабораторных экспериментах, может быть вызван нелинейностью действительного фильтрационного закона сопротивления. Нелинейность в математической модели не учитывалась ввиду отсутствия экспериментальных значений проницаемости на больших расходах газа.

1,9 1,8 1,7 1,6 1,5 1,4 1,3

1,6 3,2 4,8 6,4 8 9,6 *,мс

Рис. 6. Изменение безразмерного давления во времени в точке установки датчика. 1,2- лабораторные эксперименты №№ 474,475; 3 - численный эксперимент для условий №№ 474, 475; 4 - лабораторный эксперимент № 486; 5 - численный эксперимент для условий № 486. Начальное давление в КВД 0,2 МПа.

Сравнение расчетной и экспериментально определенной скорости волны разрежения. Экспериментальная скорость волны разрежения по газу определялась по времени прохода между: 1) датчиками, установленными на расстоянии 200 мм друг от друга в боковой стенке КВД; 2) боковым датчиком и датчиком, установленным в глухом торце КВД. Теоретическая скорость волны определялась по времени прохода между расчетными боковыми "датчиками". В табл.1 сведены результаты сравнения теоретически и экспериментально определенной скорости.

Таблица 1 - Скорость волны разрежении в пористом среде по данным численных и лабораторных экспериментов

№ лабораторного эксперимента 480 481 486 487* 488*

Измеренная скорость, м/с 126,7 126,1 138,4 179 185,2

Теоретическая скорость*, м/с 142,5

Отклонение от теории, % -12,5 -12,1 -3,2 20,4 23,1

* - значения определены по сигналам боковых датчиков

Экспериментальные значения скорости удовлетворительно согласуются с теоретическими. Среднее значение экспериментальной скорости составляет 151,1 м/с, теоретическое - 142,5 м/с. Результаты теоретических расчетов падения давления в точке установки датчика согласуются с экспериментальными значениями со средним расхождением 9%. Расхождения объясняются отклонением коэффициента проницаемости образца, участвовавшего в конкретном эксперименте, от принятого в расчетах значения вследствие неоднородности материала.

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

До 2003 г. на Малкинском геотермальном месторождении суммарный дебит добычных скважин был ниже требуемого по условию бесперебойного теплоснабжения, что было вызвано потерями давления на транспорт от устьев скважин до газосепаратора, а также потерями напора на преодоление высоты перегородки в газосепараторе. В целях увеличения добычи был реализован разработанный автором, совместно с М.А. Петровым, метод управления гидрогазодинамическими процессами в эксплуатационных колоннах добычных скважин. Вблизи устьев скважин №№12, 13 были установлены насосы, снимающие гидравлические сопротивления трубопроводов и газосепаратора. При работающих насосах скважины эксплуатируются в режиме свободного самоизлива с нулевым избыточным уровнем на устье, что обеспечивает увеличение их дебита. Разработана и реализована схема обвязки добычной слабогазирующей геотермальной скважины, в которой высота встраивания в эксплуатационную колонну трубопровода от скважины к насосу выбиралась по результатам расчетов и

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

На месторождении в качестве добычной используется также скважина №1, эксплуатирующаяся с помощью насоса, установленного на водоподъемной трубе выше оголовка. При работе скважины с требуемым дебитом уровень теплоносителя в ней опускается ниже устьевой отметки. Для установления условий устойчивой работы скважины в автоматическом режиме и предложенной схеме обустройства и автоматизации геотермального промысла были проведены эксперименты в производственных условиях. Устье скважины было герметизировано, кран на сброс газа из затрубья полностью закрыт. При дебите скважины 2,5 л/с расход газа составил 0,012 л/с. Статический уровень воды до начала откачки над колонным фланцем скважины составлял +2,30 м. вод. ст. Режим потока смеси вода-газ по номограмме Уоллиса был установлен снарядным, что подтверждалось прямыми наблюдениями: серия газовых снарядов чередовалась потоком воды.

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

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

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

Разработана и внедрена схема обустройства при автоматизированной эксплуатации геотермального месторождения, рис. 7.

^__Потребитель ЛЫ ~ Малкинский

лососевый рыбоводный завод

Рис. 7. Схема обустройства и автоматизации геотермального промысла. 1 - добычные скважины; 2 - запорная арматура; 3 - насос; 4 - щит управления насосом; 5 - щит управления скважинами; 6 - датчик давления; 7 - емкость-газосепаратор.

Предложен и реализован комбинированный технологический режим автоматизированной эксплуатации геотермальных месторождений. Добычные скважины №№12 и 13 работают основную часть времени в режиме самоизлива, с избыточным напором на устье, на трубопроводы до газосепаратора №1. В пиковый режим водоразбора по сигналу от датчика давления на магистральном трубопроводе от газосепаратора №1 включаются насосы, установленные на скважинах №№12 и 13, в результате чего последние переходят в режим свободного самоизлива. Дебит водоотбора увеличивается. Режим является комбинированным из двух: самоизлива и кранового. Скважина №1 является резервной, насос на ней автоматически включается и выключается по датчику давления, регистрирующего низкий или высокий уровень в газосепараторе №2. Основную часть времени работают скважины №№12 и 13, более глубокие, чем скважина №1, и перехватывающие глубинный поток теплоносителя с более высокой температурой. В результате выполняется одно из условий охраны недр, запрещающее разубоживание ресурсов геотермальных месторождений.

Представленные выше рекомендации позволили увеличить объем добычи на месторождении до величины защищенных запасов и обеспечить беспере-

^ s S s

& s

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

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

Численная трехмерная термогидродинамическая модель Мутновского месторождения иарогидротерм была реализована на базе программного комплекса TOUGH2, версия 2.0 [Pruess et ah, 1999], разработанного в Берклиевской лаборатории имени Лоурепса (LBNL, США). Вычислительная сетка разработана на базе программного комплекса WinGridder (LBNL), версия 2.0.1 [Pan, 2003]. Программный комплекс TOUGH2 численно реализует интегро-дифференциальным методом уравнения сохранения массы и энергии в проницаемой пористой или трещиноватой срсде, в приближении локального термодинамического равновесия, когда температуру флюида и скелета пород в представительном объеме можно считать равными. Использовался модуль реализации уравнения состояния EOS1 для чистой воды с максимальной температурой 350°С в соответствии с международной формуляцией свойств воды и водяного пара IF-67.

Для установления основных параметров геометрии численной модели автором, совместно с И.И. Черневым и A.B. Шадриным, на основе интерпретации и интерполяции данных по 12 геологическим разрезам месторождения, была разработана трехмерная геологическая модель месторождения с учетом фактического рельефа. Для выполнения численных экспериментов модель была упрощена и схематизирована до плоской трехмерной геометрической. Использовался подход одной пористости. В плане модель ограничена прямоугольником площадью 87 км2 (7,5x11,6 км) в координатной системе Гаусса-Кригера. Использовалась нерегулярная вычислительная сетка с размерами наиболее мелких блоков, включающих источники и скважины в горизонтальной плоскости, 25x25 м. Координаты центров блоков со скважинами и источниками соответствуют фактическим. Модель содержит 41472 блока.

По глубине модель распространяется от средней топографической поверхности участка (абс. отм. 700 м) до 2000 м (абс. отм. -1300 м) и включает пять слоев толщиной 75, 25, 800, 900 и 200 м. Первый, второй и третий слои соответствуют Среднему гидрогеологическому ярусу месторождения. Локальные

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

Гидрогеологические и теплофизические параметры пород были первоначально введены в модель равными данным, установленным по результатам откачек и опытных работ в процессе разведки и опытно-промышленной эксплуатации. Затем параметры уточнялись по результатам моделирования естественного состояния и эксплуатации. Теплоемкость и теплопроводность пород во всех блоках модели приняты постоянными: 1 кДж/кг-К и 2 Вт/м-К.

Границы модели в плане задавались неактивными. Объем граничных блоков был установлен значительным (1050 м3), поэтому их термодинамические параметры в численных экспериментах оставались постоянными. Границы модели являются границами постоянного давления, равного гидростатическому до начала эксплуатации, и постоянной температуры, в соответствии с геотермическим градиентом 70°С/км, полученным интерполяцией данных термометрии при бурении скважин №5 и №7, в точках, приходящихся на граничные блоки модели. Границы модели в плане определяют "стоки" теплоносителя за пределы участка, если давление теплоносителя внутри модели превышает давление в граничных блоках. В случае снижения давления внутри области ниже, чем на границах, из граничных блоков начинается поступление теплоносителя.

Глубинные "источники" геотермального теплоносителя задавались на нижнем слое модели. Они представлены 453 блоками с общим расходом 520 кг/с с постоянной энтальпией 1,6 МДж/кг, что соответствует температуре воды 341°С на линии насыщения при давлении 14,8 МПа, примерно равному гидростатическому давлению на глубине 2 км. Расход питания в модели задавался различным на участках: Дачный и Восточный блок - каждый по 140, Верхне-Мутновский и Вулканный - по 60 кг/с и распределялся между отдельными блоками равномерно, пропорционально площади горизонтальной поверхности блоков. Распределение блоков-"источников" представлено на рис. 8. На поверхности нижнего слоя, свободной от блоков-"источников", задавался тепловой поток с удельной плотностью 140 мВт/м", соответствующей начальному геотермическому градиенту 70°С/км и теплопроводности пород 2 Вт/м-К.

Достигнуто соответствие распределения в плане модельных и фактически измеренных температур, рис. 9-10. С увеличением глубины сходимость фактических и модельных геоизотерм повышается. Максимальное отклонение фактических и расчетных температур составило 16,8% при среднем отклонении 11%. На рис. 10 представлено сопоставление типичных фактической и расчет-нон термограмм. Сравнение фактических и расчетных термограмм производилось по 63 скважинам. Далее калибровка модели производилась по распределению температур в ходе эксплуатации месторождения. Окончательная настройка модели производилась с помощью изменения проницаемости отдельных блоков

Рис. 8. Распределение блоков питания и областей модели Мутновского месторождения паро-гидротерм. Числа величина питания на участках, кг/с. Цвет и римские цифры: I - высокопроницаемая область, II - среднепроницаемая, III - низкопроницаемая, IV - блоки с источниками основного питания, V- блоки периферийного питания. Линии: жирная сплошная граница площади водоупора, прямые жирные - разрезы, красные - установленные разломы.

Рис. 9. Геоизотермы Мутновского месторождения парогидротерм на абсолютной отметке 0 м: слева - построенные по результатам термометрии на стадии разведки; справа - полученные в численном моделировании естественного состояния.

Скв. 09

• \

1.....1 •• \ - д V

н. м

1500

О 50 100 150 200 250 Л "С 0 100 150 200 250 300 /."С

Рис. 10. Термограммы скважин Мутновского месторождения парогидротерм до начала эксплуатации: точки - термометрия в ходе бурения; линии - данные модели.

модели, приходящихся на эксплуатационные скважины, при этом использовался параметр РМХ во входном файле TOUGH2.

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

Выполнено сравнение фактических и расчетных энтальпий в блоках модели, соответствующих верхнему интервалу зон водопритоков в добычных скважинах. Давление, энтальпия и температура в блоках, содержащих скважины, были рассчитаны с помощью программного комплекса HOLA [Bjornsson et al., 1993]. Энтальпия пароводяной смеси на устье скважин рассчитывалась по таблицам, составленным по уравнениям свойств воды и водяного пара Международной системы уравнений 1997 г. для промышленных расчетов. Максимальное отклонение фактических и расчетных параметров не превосходит 16,3%. Отклонения объясняются отсутствием данных о величине скин-эффекта скважин.

По результатам моделирования были сделаны следующие выводы:

- естественные ресурсы месторождения можно оценить в 520 кг/с теплоносителя с энтальпией 1,6 МДж/кг, что соответствует температуре воды 341°С на линии насыщения при давлении 14,8 МПа. В тепловом эквиваленте указанная величина естественных ресурсов соответствует 832 МВт;

- возвратная закачка отработанного теплоносителя в полном объеме (сепа-рат и конденсат) в интервале глубин 600-1500 м оказывает существенное негативное влияние на термогидродинамический режим месторождения, особенно на Верхне-Мутновском участке.

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

Численное моделирование теплопсреноса в породах магмагеотермальных систем выполнялось на базе программного комплекса HYDROTHERM, Геологическая служба США, версия 2.2 и 3.0 [Hayba and Ingebritsen, 1994; Kipp et al., 2008], который конечно-разностным методом численно реализует уравнения сохранения массы и энергии в том же приближении, что и TOUGH2.

Дополнительно накладывается условие пренебрежения обратимой работой сжатия теплоносителя, что практически всегда выполняется в двухфазной и всегда выполняется в надкритической области состояния воды. HYDRO-THERM использует таблицы свойств чистой воды до максимальных значений температуры и давления 1200°С и 1 ГПа, что определило применение его для условий магмагеотермальных систем с надкритическими термодинамическими параметрами.

Концептуальная модель Мутновекой магмагеотермальной системы включает вертикальную проницаемую зону, проходящую через центр активного кратера Мутновского вулкана вдоль предполагаемого разлома в северо-северо-восточном направлении [Вакин, 1966]. Толщина зоны принята равной размеру активного кратера (150 м). Поверхность активного кратера гидравлически связана посредством вертикальной проницаемой флюидопроводящей зоны с областью дегазации магматического очага. Граница модели: по простиранию - 9 км от центра кратера, по вертикали - определяется варьируемой глубиной залегания очага. Буровые и геофизические работы в непосредственной близости к центру системы не проводились. Исходя из отсутствия надежных данных по геометрическим, теплофизическим, гидрогеологическим и геотермическим параметрам и свойствам трех областей модели (флюидопроводящая зона, магматический очаг и окружающие горные породы), автором, совместно с В.В.Таскиным, производились численные эксперименты на четырех различных типах моделей с разными начальными и граничными условиями, а также в широком диапазоне теплофизических и гидрогеологических параметров, табл. 2.

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

В численных экспериментах на модели А была установлена невозможность достижения расчетных температур газов в активном кратере, близких к измеряемым при глубинах залегания очага от 3 до 4,5 км, табл. 3. Адекватная модель системы должна включать в явном виде источники дегазирующегося флюида на входе во флюидопроводящую зону. В ходе численных экспериментов на модели В установлено превышение расчетных температур газов над измеренными. В экспериментах на модели С установлена общая динамика процесса в двух вариантах доминирующего теплопереноса в горных породах -кондуктивном и конвективном, а также изучено влияние на интенсивность процесса и его характеристики глубины залегания очага и давления флюида на входе во флюидопроводящую зону.

Таблица 2 - Типы моделей и параметры, варьируемые в численных экспериментах

Тип модели А В С D А, В, С, О

Форма очага (эллипс) горизонтальны й вертикальный, горизонтальный горизонтальный горизонтальный вертикальный, горизонтальный

Размеры большой и малой полуосей очага, км 2,2; 1,1 2,5; 1,5 •ч") 1 у1 2,2; 1,1 2,2...2,5; 1,1...1,5

Глубина очага, км 3,0...4,5 1,0 0,75...4,5 1,5 0,75...4,5

Температура апикальной части очага, "С 1200 900... 1200 1200 900 900... 1200

Проницаемость пород флюидопроводящей зоны, мД 10Л..102 10"3...1 1...103 1...102 юЛ.ло2

Проницаемость окружающих пород, мД 10"*...1 10° 10Л..1 10ц Ю-4...!

Проницаемость очага, мД 10'8...10"4 - - - юЛ.ло-4

Региональный тещговой поток, мВт/м* 80...120 120 120 120 80... 120

Тепловой поток на оси очага, мВт/м" 1200 - - - 1200

Расход дегазирующегося флюида (С), кг/с - 3...600 71 71 0...600

Давление дегазации - - Р*гст---Р*лст Ргст---2Ргст РгсТ' • -Рлст

Доминирующий тип теплопереноса в окружающих породах кондуктивный, конвективный кондуктивный, конвективный кондуктивный, конвективный кондуктивный кондуктивный, конвективный

Ргст. Рлст - гидростатическое и литостатическое давление на заданной глубине

Таблица 3 - Измеренная и полученная в численных экспериментах температура ("С) газов активного кратера Мутновской магмагеогермалыюй системы

Период наблюдений, годы Измеренная температура Температуры, полученные в численных экспериментах

А В С О

1963-2005 520-950 50-100 900-1000 500-700 525-840

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

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

Численные эксперименты с моделью £> позволили установить ожидаемые значения градиентов температуры (до 260°С/км) и давления (до 14 кПа/км) во флюидопроводящей зоне при варьировании ее проницаемости от 1 до 100 мД и изменении давления дегазации магмы от 15 до 30 МПа.

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

По результатам численных экспериментов с моделью С в 27-ми вариантах набора исходных параметров для случая доминирующего кондуктивного теп-

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

При малых глубинах очага (0,75 и 1,5 км) расчетная температура флюида вблизи поверхности кратера с ростом давления дегазации Рм сначала растет, а затем стабилизируется и начинает снижаться. Максимумы температуры ТАВ с ростом давления дегазации снижаются и сдвигаются в сторону его больших значений. При глубинах залегания магматического очага больше 2 км расчетная температура фумарол с увеличением давления дегазации в пределах от РГст до Рлст в принятом в расчетах диапазоне Рм только растет. Температура фумарол, соответствующая средним величинам давления дегазации между РГСт до Рлст для фиксированной глубины залегания магматического очага (линия Рср), почти линейно снижается с ростом Рм.

1000

Рис. 11. Зависимость температуры газов активного кратера Мутновской магмаге-отермалыгой системы от давления на входе во флюидопроводящую зону при разных

глубинах очага.

Задавая давления на рис. 11, можно по измеренной температуре фумарол оценить глубину залегания очага тремя следующими способами. • По среднему давлению дегазации в диапазоне от гидростатического до лито-статического на заданной глубине очага (Р|=32,7 МПа). Расчетная температура фумарол, определяемая по точке пересечения линий Рср и Рь составит Т1=810°С. Среднее значение фактически наблюдаемой температуры фумарол составляет 720°С за период 1963-2005 гг. Отклонение температуры фумарол,

определенной по графику от средней фактической составляет 13%. При этом расчетная глубина очага составляет Н1=2 км.

• По среднему значению давления в диапазоне от гидростатического на глубине очага 0,75 км до литостатического на глубине очага 3 км - Рг = 41,3 МПа. По графику Т2=775°С, отклонение от средней фактической - 8%, Нг=2,5 км.

• По средней температуре фумарол за период наблюдения - Тз=720°С и давлению Р2=41,3 МПа. Н3=2,8 км.

Расхождение между оценками глубин тремя различными способами не превышает 40% при среднем значении 2,5 км. Вероятной глубиной залегания магматического очага можно считать 2,5 км от поверхности активного кратера.

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

На базе программного комплекса НУОКОТНЕЯМ, автором, совместно с В.В.Таскиным, выполнено численное моделирование ГЦС типа "триплет" - одна нагнетательная и две добычных скважины в начальных надкритических термодинамических условиях, соответствующих ожидаемым в магмагеотермальных системах. В качестве начальных надкритических условий на забое добычных скважин приняты: температура 380°С, давление 22,5 МПа. Дебит нагнетательной скважины варьировали в пределах 5, 10 и 15 кг/с; глубина добычных скважин - 2,7 км; высоту забоя нагнетательной скважины над линией забоев добычных скважин (Н„) задавали 0, 250, 450, 550 м, расстояние между нагнетательной и добычной скважиной (Ь) - 200, 250, 450, 500 м; температуру теплоносителя на забое нагнетательной скважины принимали равной 100°С. Расчетная область размером 1,95x1x3 км включала продуктивную зону горных пород с проницаемостью 1 мД и две зоны малопроницаемых пород (10"5 мД): верхнюю мощностью 1 км и нижнюю - 0,2 км. Начальное распределение давления -гидростатическое; на верхней границе задавали атмосферное давление и температуру 10"С. Градиент температуры - 137°С/км, начальная температура на нижней границе - 420°С. Мощность открытых интервалов скважин - 50 м.

Изменение параметров теплоносителя в зависимости от Ь при Н„=0. Для установления рациональных параметров надкритической ГЦС важным является анализ изменения фазового состояния теплоносителя в ходе эксплуатации, а также исследование закономерностей этого изменения в зависимости от взаимного расположения добычных и нагнетательной скважин. В случае расположения скважин «в ряд» (Нн = 0) с Ь = 500 м, через 1 год с начала эксплуатации, в результате отбора теплоносителя и вызванного этим снижения давления, в продуктивной зоне вблизи забоев добычных скважин формируется зона влажного пара с водонасыщенностью, увеличивающейся в направлении от границ области влажного пара к забою. В остальной части продуктивной зоны состояние теплоносителя остается надкритическим. Через 2 года с начала эксплуатации фронт конденсации распространяется на большую часть сечения продуктивной

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

На забое нагнетательной скважины и вблизи нее давление теплоносителя остается выше критического. Зона перегретого пара формируется в области ниже и на периферии забоя добычных скважин, расширяясь в ходе эксплуатации к границам продуктивной зоны. При уменьшении расстояния между нагнетательной и добычными скважинами до Ь = 250 (200) м процесс фазовых переходов в продуктивной зоне развивается аналогично, однако его интенсивность и средняя водонасыщенность зоны увеличиваются. При фиксированных Н„ расстояние между скважинами определяет линейный размер области возмущенного состояния продуктивной зоны, темп охлаждения горных пород и величины градиентов давления, температуры и водонасыщенности в окрестности забоя скважин. Чем меньше X, тем интенсивнее охлаждаются породы и тем выше величина депрессии в призабойной части добычных скважин. Линейный размер области существенных градиентов водонасыщенности и ее средняя величина, наоборот, уменьшаются при сокращении Ь.

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

Установившиеся значения на забое добычных скважин при дебите 5, 10, 15 кг/с составляли: давления - 19, 15, 8,5 МПа; температуры - 360, 340, 300°С; энтальпии - 1900, 1700, 1400 кДж/кг; массовой водонасыщенности - 85, 92, 96%; требуемое давление нагнетания - 26, 34, 43 МПа. В целях установления рациональных геометрических параметров ГЦС, соответствующих эффективной работе системы, для разных наборов геометрических параметров взаимного расположения забоев добычной и нагнетательной скважин вычислялся отбор теплоты из продуктивной зоны. В результате сделаны следующие выводы.

Дебит скважин, проницаемость продуктивной зоны, относительная высота забоя нагнетательной скважины над линией добычных, расстояние между добычными и нагнетательной скважинами влияют на технологические параметры ГЦС. При Н„~0 рационально располагать забои добычных и нагнетательной скважин на расстоянии 250 м, при этом параметры добычных скважин на конец срока разработки выше, чем в случае с ¿=500 м, а требуемое давление нагнетания ниже, чем при ¿=200 м. Отбор тепла из продуктивной зоны ГЦС типа "три-

плет" зависит от относительной высоты забоя нагнетательной скважины над линией добычных и расстояния между скважинами.

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

В пределах Мутновской магмагеотермалыюй системы возможно расположение до 5-ти групп скважин типа "триплет", которые могут вскрыть на глубинах до 4-х км теплоноситель с температурой 400°С. На одной системе "триплет" возможно строительство ГеоЭС электрической мощностью 34 МВт при дебите теплоносителя 140 кг/с. В табл. 4 сведены показатели пилотной ГеоЭС, рассчитанные на базе программного комплекса ЕС5, разработанного в Массачусетсом технологическом институте (США), в реальных ценах 2007 г.

Таблица - 4 Расчетные технико-экономические показатели проекта ГеоЭС,

использующей тепловые ресурсы магматического очага (в ценах 2007 г.)

Капитальные затраты (тыс. долл. США)

Разведочные работы 6764,87

Бурение 14066,89

Стимуляция резервуара (гидроразрыв и др.) 7732,89

ГеоЭС 24555,36

Система транспорта теплоносителя 8215,55

Трансформаторная подстанция 110/35/10 кВ 2004,43

Линия В Л 110 кВ 1995,89

Суммарные капзатраты 65335,88

Эксплуч анионные расходы (тыс. долл. США/год)

Система скважин 545,33

ГеоЭС 2099,53

Теплоноситель 42,60

Подстанция и линия 80,01

Суммарные эксплутационные расходы 2767,47

Расчетный отпускной тариф (центы США/кВт ч) 9,70

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

ЗАКЛЮЧЕНИЕ

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

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

Основные теоретические положения и практические результаты работы заключаются в следующем.

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

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

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

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

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

6. Разработаны и внедрены на Малкинском месторождении термоминеральных вод новые способы управления работой добычных геотермальных скважин, выводящих теплоноситель с незначительным газовым фактором:

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

• скважины, работающей с понижением уровня в затрубном пространстве.

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

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

8. Разработаны численные термогидродинамические модели магма-геотермальных систем вулканов Мутновский, Авачинский, Кудрявый, а также Мутновского месторождения парогидротерм, откалиброванные по всей совокупности доступных данных геотермического и гидродинамического режимов. Произведена оценка величины естественных ресурсов Мутновского месторождения парогидротерм - 520 кг/с теплоносителя с энтальпией 1,6 МДж/кг, что соответствует температуре воды 341°С на линии насыщения при давлении 14,8 МПа и тепловому потоку 832 МВт.

9. Установлены расчетные распределения температуры, давления, фазового состояния и скоростей фаз флюида в массиве пород Мутновской магмагеотер-малыюй системы в четырех вариантах краевых условий, определяющих начальное термодинамическое состояние пород и условия на границах модели, двух вариантах доминирующего типа теплопереноса в породах - конвективного и кондуктивного, а также для двух сценариев теплового режима очага - остывающего и конвектирующего.

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

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

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

Основные научные результаты диссертации опубликованы в следующих работах:

Монография

1. Пашкевич Р.И., Таскин В.В. Термогидродинамическое моделирование теп-лопереноса в породах Мутновской магматогенной геотермальной системы. -Владивосток: Дальнаука, 2009.-250 с.

Статьи в журналах и изданиях "Перечня ведущих рецензируемых научных

журналов и изданий" ВАКМинобрнауки России

2. Пашкевич Р.И. Гидростатические парадоксы в геотермальных коллекторах // Вулканология и сейсмология.-1996.-№3.-С. 101-106.

3. Пашкевич Р.И., Таскин В.В. Численное исследование динамики паровых зон магмагеотермальных систем на примере вулкана Мутновский, Южная Камчатка // Геофизика.-2007.-№5.-С. 68-71.

4. Пашкевич Р.И., Таскин В.В. Термогидродинамическое моделирование магмагеотермальной системы вулкана Мутновский // Геоинформатика-2007.-№ З.-С. 57-61.

5. Пашкевич Р.И., Таскин В.В. Моделирование теплопереноса при фильтрации во флюидопроводящей зоне Мутновской магмагеотермальной системы: роль проницаемости пород и давления дегазации магмы // Вопросы современной науки и практики. Университет им. Вернадского. - 2007. №4(10)-С. 156-165.

6. Пашкевич Р.И., Таскин В.В. Термогидрогазодинамические процессы в массиве пород вулкана Мутновский при остывании его магматического очага: численный эксперимент // Естественные и технические науки.-2007-№5.-С. 169-174.

7. Пашкевич Р.И., Таскин В.В. Теплоперенос в породах Мутновской геотермальной системы: влияние формы, расхода дегазации очага и проницаемости пород //Горный информационно-аналитическийбюллетень-2008-№12.-С. 178-185.

8. Таскин В.В., Пашкевич Р.И. Теплоперенос в горном массиве при эксплуатации геотермальной циркуляционной системы // Горный журнал.-2008.-№11.-С. 40-42.

9. Пашкевич Р.И., Таскин В.В. Зависимость технологических параметров геотермальных циркуляционных систем в начальных надкритических условиях от относительной высоты забоя нагнетательной скважины // Горный информационно-аналитический бн>ллетень.-2008.-№12.-С. 186-190.

10. Таскип В.В., Пашкевич Р.И. Исследование термогидродинамического режима надкритической геотермальной циркуляционной системы //Известия высших учебных заведений. Горный журнал.-2008.-№8.-С. 162-164.

11. Пашкевич Р.И., Таскин В.В. Влияние расстояния между скважинами и глубины нагнетательной скважины надкритической геотермальной циркуляционной системы на отбор тепла // Горный информационно-аналитический бюллетень.-2008.-№12.-С. 191-194.

12. Пашкевич Р.И. Численное моделирование фильтрации при быстрой декомпрессии неподвижной газонасыщенной пористой среды // Газовая промышленность.-2009.-№7.-С. 40-42.

13. Пашкевич Р.И. Теплоперенос в горных породах магмагеотермалыюй системы Кудрявая (Курильские острова) // Горная промышленность-2009 — №3.-С. 60-61.

14. Пашкевич Р.И., Чернев И.И., Шадрин A.B. Термогидродинамическое моделирование Мутновского месторождения парогидротерм // Разведка и охрана недр.-2009.-№7.-С. 37-43.

15. Пашкевич Р.И. Осреднение по межфазной поверхности в фильтрационном потоке // Геофизика.-2009.-№4.-С. 63-65.

16. Пашкевич Р.И., Петров М.А. Способы управления и автоматизации добычи на геотермальных месторождениях // Маркшейдерия к недропользова-ние.-2009.-№4.-С. 36-37.

17. Паршин Б.Е., Муратов П.В., Пашкевич Р.И. Пути повышения эффективности работы воздушной конденсационной установки Верхне-Мутновской геотермальной электростанции на Камчатке // Электрические станции-2007.-№6.-С. 40-43.

18. Москалев Л.К., Таскин В.В., Пашкевич Р.И. Повышение эффективности работы сепараторов первой ступени Мутновской геотермальной электростанции путем предотвращения минералоотложений из теплоносителя И Гидротехническое строительство.-2006.-№12.-С. 20-24.

Статьи в других изданиях

19. Панов В.К., Пашкевич Р.И. Экспериментальное исследование мгновенной декомпрессии неподвижной газонасыщенной пористой среды // Тепломассообмен гетерогенных сред: Труды ПКВМУ. Вып. 1. -Петропавловск-Камчатский: РИО ПКВМУ, 1996.-С. 81-89.

20. Pashkevich R.I. About convective heat transfer in geothermal systems // Geo-thermal Resources Council TRANSACTIONS, Vol. 20, September/October 1996, Davis, California, USA, pp. 829-832.

21. Pashkevich R.I. On environmental aspects of geothermal development // Geothermal Resources Council TRANSACTIONS, Vol.20, September/October 1996, Davis, California, USA, pp. 241-243.

22. Pashkevich R.I. The contradictions in continua modelling of flow and heat transfer in fractured and/or porous reservoirs // Geothermal Resources Council TRANSACTIONS, Vol. 21, September/October 1997, Davis, California, USA, pp. 617-620.

23. Пашкевич Р.И. О методах включения капиллярных эффектов в модели двухфазного потока в геотермальных резервуарах // Инженерно-физические исследования на Камчатке: Труды КГАРФ. Вып. 9,-Петропавловск-Камчатский: Изд-во КГАРФ, 1999.-С.46-64.

Материалы конференций

24. Pashkevich R.I. Kamchatka geothermal resources development: problems and perspectives // Proc. Stanford Geothermal Workshop 18, Stanford University, California, USA, 1996, pp. 539-541.

25. Pashkevich R.I. Results and interpretations of hot and cold water injection experiments on Pauzhetka water-dominated geothermal field in Kamchatka // Proc. Stanford Geothermal Workshop 18, Stanford University, 1996, pp. 543-547.

26. Pashkevich R.I. On a capillary pressure in two-phase reservoirs //Proc. 18th New Zealand Geothermal Workshop, New Zealand, Auckland, 1996, pp. 291-294.

27. Pashkevich R.I. On possibility of HDR project in near-by region of Petropav-lovsk-Kamchatsky, Russia // Proc. 3rd International HDR Forum, May 13-16, 1996, Santa Fe, New Mexico, USA, pp. 144-145.

28. Pashkevich R.I. Method of capillary effects inclusion in two-phase flow modelling of a geothermal reservoir // Proc. Stanford Geothermal Workshop 20, Stanford University, California, USA, 1998, pp. 420-427.

29. Пашкевич Р.И. Перспективы ресурсо- и энергосбережения в геотермальном теплоснабжении на примере Камчатской области // Труды I Международной научно-практической конференции "Эффективные энергетические системы и новые технологии", Казань, 4-6 декабря 2001 г.-Казань: Изд-во Казан, гос. тех. ун-та, 2002. -С. 205-210.

30. Трухин Ю.П., Пашкевич Р.И., Таскин В.В. Термогидродипамическое моделирование магматогенно-гидротермалыюй системы вулкана Мутновский и его магматического очага // Минерально-сырьевые ресурсы как фактор развития промышленной и социальной инфраструктуры Камчатского края: материалы региональной научно-практической конференции, 1-2 марта 2007 г. Отв. ред. Ю.П. Трухин.-Петропавловск-Камчатский: НИГТЦ ДВО РАН, 2007.-С. 41-43.

31. Пашкевич Р.И., Таскин В.В. О возможности оценки глубины залегания очага Мутновской магмагеотермалыюй системы по геотермическим данным // Тепловое поле Земли и методы его изучения. Сб. науч. трудов X Международной конф., Москва, 18-19 июня 2008 г. /Отв. ред. Ю.А.Попов. -М.: РИО РГГРУ, 2008.-С. 176-181.

32. Pashkevich R.I., Taskin V.V. Heat transfer in geothermal system of Mutnovsky volcano: the influence of the form, discharge of magma chamber degassing and rocks permeability// Proc. Stanford Geothermal Workshop 34, Stanford University, Stanford, California, USA, 2009, pp. 88-93.

33. Pashkevich R.I., Taskin V.V. Numerical simulation of exploitation of supercritical enhanced geothermal system // Proc. Stanford Geothermal Workshop 34, Stanford University, Stanford, California, USA, 2009, pp. 215-223.

Подписано в печать 07.08.09. Заказ 145. Печать цифровая. Тираж 100 экз. Оперативная полиграфия ИП Романенко М.И. 683000, г.Петропавловск-Камчатский, ул. Ленинская, 46. тел. +7-(4152)-423-008.

Содержание диссертации, доктора технических наук, Пашкевич, Роман Игнатьевич

Основные условные обозначения.

Введение.

1. Анализ проблемы обоснования рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки.

1.1. Моделирование теплопереноса при фильтрации с учетом капиллярных сил.

1.2. Процесс быстрой декомпрессии проницаемой среды.

1.3. Проблема установления параметров теплопереноса в породах магмагеотермальных систем.

1.3.1. Численное моделирование магмагеотермальных систем.

1.3.2. Характеристика Мутновской магмагеотермальной системы.

1.3.3. Существующие концептуальные модели Мутновской магмагеотермальной системы.

1.4. Направления освоения тепловых ресурсов близповерхностных магматических очагов.

1.4.1. Извлечение тепловой энергии из магмы.

1.4.2. Добыча теплоты горных пород с помощью геотермальных циркуляционных систем.

1.4.3. Добыча геотермального теплоносителя высоких параметров из глубокозалегающих геотермальных коллекторов.

1.5. Постановка задач исследования.

2. Теоретические основы моделирования теплопереноса при фильтрации двухфазного потока в продуктивных зонах геотермальных месторождений.

2.1. Теоремы осреднения по межфазной поверхности в системе двухфазный (жидкость-пар) геотермальный теплоноситель-порода».

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

2.2.1. Новая система уравнений.

2.2.2. Сравнение существующих и предлагаемой модели.

2.2.3. Уравнения сохранения импульсов.

2.2.4. Градиент насыщенности как функция кривизны межфазной поверхности.

2.3. Преодоление гидростатических парадоксов в геотермальных коллекторах.

2.4. Выводы к главе 2.

3. Экспериментальное и численное исследование процесса быстрой декомпрессии неподвижной газонасыщенной пористой среды.

3.1. Эксперименты в лабораторных условиях.

3.2. Численные эксперименты.

3.3. Сопоставление результатов лабораторных и численных экспериментов.

3.3.1. Характер и темп снижения давления газа.

3.3.2. Скорость волны разрежения. 124'

3.4. Выводы к главе 3.

4. Методы и средства управления гидрогазотермодинамическими процессами в стволе добычных скважин геотермальных месторождений.

4.1. Анализ режима эксплуатации Малкинского геотермального месторождения.

4.2. Технологические решения по реконструкции промысла на Малкинском геотермальном месторождении.

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

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

4.5. Технологический режим и схема автоматизации геотермального промысла.

4.6. Выводы к главе 4.

5. Численное термогидродинамическое моделирование Мутновского месторождения парогидротерм.

5.1. Краткие сведения о месторождении.

5.2. Разработка концептуальной и численной модели месторождения.

5.3. Результаты моделирования.

5.3.1. Моделирование естественного состояния месторождения.

5.3.2. Моделирование эксплуатации.

5.3.3. Калибровка модели эксплуатации месторождения.

5.4. Выводы к главе 5.

6. Численное моделирование теплопереноса в породах магмагеотермальных систем.

6.1. Концептуальная и численная модель Мутновской магмагеотермальной системы.

6.2. Численные эксперименты по установлению параметров теплопереноса.

6.2.1. Остывающий магматический очаг в отсутствии дегазации.

6.2.2. Конвектирующий дегазирующий очаг.

6.2.3. Конвектирующий дегазирующий очаг, варьирование давления дегазации.

6.2.4. Конвектирующий дегазирующий очаг при постоянной начальной температуре в породах.

6.2.5. Метод оценки глубины залегания очага Мутновской магмагеотермальной системы по геотермическим данным.

6.2.6. Динамика паровых зон в породах Мутновской магмагеотермальной системы.

6.3. Кондуктивный теплоперенос в магмагеотермальной системе вулкана Авачинский.

6.4. Теплоперенос в породах Кальдеры Медвежья и вулкана Кудрявый.

6.5. Выводы к главе 6.

7. Технологические параметры и технико-экономическая оценка эффективности геотермальной циркуляционной системы при освоении ресурсов магматического очага.

7.1. Зависимость параметров циркуляционной системы от относительной глубины нагнетательной скважины.

7.2. Динамика параметров продуктивной зоны в ходе эксплуатации.

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

7.4. Влияние геометрических параметров системы на выработку тепловой энергии.

7.5. Технико-экономическая оценка эффективности освоения тепловых ресурсов близповерхностных магматических очагов.

7.6. Выводы к главе 7.

Введение Диссертация по наукам о земле, на тему "Научно-техническое обоснование рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки"

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

Многолетние теоретические и натурные исследования термогидродинамического и гидрогеохимического режима Мутновской магмагеотермальной системы выполнены Е.А. Вакиным, И.Т. Кирсановым, Б.Г. Поляком, Ю.П. Трухи-ным, P.A. Шуваловым, Ю.А. Тараном, В.П. Пилипенко, труды которых внесли значительный вклад в решение проблемы формирования магмагеотермальных систем. Численные модели теплопереноса в океанической и континентальной литосфере разработаны Ю.И. Галушкиным. Численные термогидродинамические модели гидротермальных месторождений Камчатки и Курильских островов разработаны A.B. Кирюхиным.

Методика расчета параметров пароводяных потоков в стволе эксплуатационных геотермальных скважин развивалась в работах В.А. Дрознина, В.Д. Бе-лодеда и А.Н. Шулюпина. Методы измерения теплофизических свойств горных пород при высоких температурах pi давлениях разработаны Ю.А. Поповым. Вопросы технологии скважинных систем извлечения геотермальной энергии при докритических термодинамических условиях рассмотрены И.Т. Аладьевым, H.H. Ароновой, B.JI. Артемьевой, Э.И. Богуславским, В.А. Васильевым, Л.Ф. Вознюком, С.Г. Гендлером, Ю.Д. Дядькиным, А.Г. Егоровым, Г.Н. Забарным, Г.Н. Кононенко, O.A. Кремневым, Э.И. Мерзляковым, Ю.П. Морозовым, И.А. Павловым, Ю.М. Парийским, Т.Ю. Пискачевой, A.M. Пудовкиным, В.А. Романовым, И.А. Рыженко, А.Н. Саламатиным, H.H. Смирновой, В.Н. Трусовым, A.C. Цырульниковым, А.Н. Щербанем и A.B. Шурчковым.

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

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

10 носителя. Проблема установления рациональных технологических параметров геотермальных циркуляционных систем в надкритических условиях в настоящее время не решена.

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

Основные научные положения, защищаемые автором:

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

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

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

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

Достоверность и обоснованность научных положений и выводов обеспечивается: удовлетворительной сходимостью со средним расхождением 9% результатов численного моделирования процесса быстрой декомпрессии газонасыщенной пористой среды с полученными в представляемой работе данными лабораторных экспериментов; представительным объемом исходной информации по 63-м скважинам Мутновского месторождения парогидротерм (термометрия и точечные замеры высокотемпературным глубинным термометром); достаточной точностью калибровки модели естественного состояния Мутновского месторождения парогидротерм: среднее отклонение от фактического распределения температур 11% при максимальном 17%; сходимостью результатов модели эксплуатации Мутновского месторождения парогидротерм с фактическими данными мониторинга: среднее отклонение 11% по температурам добычных и наблюдательных скважин и 16% по энтальпиям на забое добычных скважин; удовлетворительной сходимостью со средним расхождением 7,5% результатов моделирования теплопереноса в Мутновской магмагеотермальной системе с фактическими данными замеров и оценок температуры и дебита фумарольных газов; продолжительной проверкой, в течение шести лет, разработанных способов управления работой эксплуатационных скважин и технологической схемы автоматизации промысла в производственных условиях на Малкинском месторождении термоминеральных вод.

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

Реализация результатов работы. Методические разработки в области численного моделирования процесса теплопереноса использованы ООО "Аква" в практике проектирования разработки месторождений термальных вод, а также при оценке запасов Мутновского месторождения парогидротерм, имеется акт внедрения. Способы управления эксплуатационными геотермальными скважинами и способ автоматизации геотермального промысла внедрены в производственную деятельность ООО "Аквариус", имеется акт внедрения. Результаты численного моделирования магмагеотермальных систем и режима эксплуатации геотермальных циркуляционных систем в надкритических усло

15 виях использованы в исследованиях теплопереноса в породах геотермальных месторождений Камчатки и Курильских островов. Результаты экспериментальных и численных исследований процесса быстрой декомпрессии газонасыщенной пористой среды использованы в учебном процессе в КамчатГТУ (19972003) и в филиале ДВГТУ (ДВПИ им. В.В. Куйбышева) в г. Петропавловске-Камчатском (2001-2007) при преподавании дисциплин "Гидрогазодинамика", "Гидравлика", "Основы технологии и проектирования ТЭК".

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

Апробация работы. Результаты исследований на различных этапах докладывались на следующих конференциях, семинарах и советах: Международной конференции "Тепловое поле Земли и методы его изучения" (Москва, 2008); расширенном семинаре лабораторий прикладной геофизики и вулканологии и сейсмотектоники Института физики Земли РАН (Москва, 2009); семинаре кафедры нефтегазовой и подземной гидромеханики Российского государственного университета нефти и газа имени М.И. Губкина (Москва, 2009); Международной научно-практической конференции "Эффективные энергетические системы и новые технологии" (Казань, 2001); Международном полевом Курило-Камчатском семинаре "Геотермальные и минеральные ресурсы областей современного вулканизма" (Петропавловск-Камчатский, 2005); Международном семинаре "Проект научного бурения на Мутновском вулкане. Исследование связи магматической и гидротермальной систем" (Петропавловск-Камчатский,

2006); региональной научно-практической конференции "Минерально-сырьевые ресурсы как фактор развития промышленной и социальной инфраструктуры Камчатского края" (Петропавловск-Камчатский, 2007); научно-технических советах ОАО "Геотерм" (Петропавловск-Камчатский, 2006-2007); Ученых советах НИГТЦ ДВО РАН (2006-2008); семинарах лаборатории геохимии и геотехнологии и лаборатории физико-химической гидродинамики НИГТЦ ДВО РАН (2004-2008); семинарах кафедры естественнонаучных и технических дисциплин филиала ДВГТУ (ДВПИ им. В.В. Куйбышева) в г. Петро-павловске-Камчатском (2001-2007); конференциях профессорско-преподавательского состава КамчатГТУ (1994-2003).

Публикации. По теме диссертации опубликовано 40 печатных работ, в том числе 17 - в журналах «Перечня ведущих рецензируемых научных журналов и изданий» Высшей аттестационной комиссии Минобрнауки России, из них 15 - в рекомендованных экспертным советом ВАК по наукам о Земле.

Объем и структура работы. Диссертация состоит из введения, семи глав, заключения, списка литературы из 231 наименования, иллюстрирована 119 рисунками, содержит 18 таблиц, общий объем - 303 страницы.

Заключение Диссертация по теме "Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика", Пашкевич, Роман Игнатьевич

7.6. Выводы к главе 7

1. Выполнено численное моделирование ГЦС типа "триплет" — одна нагнетательная и две добычных скважины в начальных надкритических термодинамических условиях, соответствующих ожидаемым в магмагеотермальных системах. В качестве начальных надкритических условий на забое добычных скважин приняты: температура 380°С, давление 22,5 МПа.

2, Дебит скважин, проницаемость продуктивной зоны, относительная высота забоя нагнетательной скважины над линией добычных, расстояние между добычными и нагнетательной скважинами влияют на технологические параметры ГЦС.

3. При Нн=0 рационально располагать забои добычных и нагнетательной скважин на расстоянии 250 м, при этом параметры добычных скважин на конец срока разработки выше, чем в случае с ¿=500 м, а требуемое давление нагнетания ниже, чем при ¿=200 м. Отбор тепла из продуктивной зоны ГЦС типа "триплет" зависит от относительной высоты забоя нагнетательной скважины над линией добычных и расстояния между скважинами.

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

5. В районе Авачинского вулкана на доступной площади бурения возможно строительство 8-ми систем скважин типа "триплет", а в районе Мутновского вулкана, на расстоянии 2-х км от центра активного кратера - 5-ти систем скважин, которые могут вскрыть на глубинах до 4-х км теплоноситель с температурой 400°С.

6. На одной системе "триплет" возможно строительство ГеоЭС электрической мощностью 34 МВт при дебите теплоносителя 140 кг/с. Полная мощность ГеоЭС, использующих тепловые ресурсы указанных выше магмагеотермаль-ных систем составит 272 и 170 МВт.

ЗАКЛЮЧЕНИЕ

Основные теоретические положения и практические результаты работы заключаются в следующем.

1. Сформулированы новые теоремы осреднения по межфазной поверхности "двухфазный теплоноситель — порода", позволяющие выразить градиенты и временные производные от среднего по межфазной поверхности потока произвольного вектора (например, потока тепла), а также средней по поверхности скалярной величины (давления, температуры и др.) через производные и изменения величин в масштабе пор или трещин. Межфазная поверхность при этом рассматривается как геометрическая поверхность раздела фаз, не обладающая особыми свойствами. Указанные формулы осреднения являются очередным шагом в решении центральной проблемы механики многофазных сред - проблемы замыкания и выражения на макроскопическом уровне межфазных потоков.

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

• скважины, работающей с понижением уровня в затрубном пространстве.

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

7. Разработан и внедрен новый способ автоматизации геотермального промысла, заключающийся в использовании комбинированного режима работы

272 скважин - самоизлива, крановом и с понижением уровня в затрубном пространстве — обеспечивающего отбор теплоносителя по текущей потребности, при этом выбор режима устанавливается по уровню в газосепараторах.

8. Разработаны численные термогидродинамические модели магмагеотер-мальных систем вулканов Мутновский, Авачинский, Кудрявый, а также Мут-новского месторождения парогидротерм, откалиброванные по всей совокупности доступных данных геотермического и гидродинамического режимов. Произведена оценка величины естественных ресурсов Мутновского месторождения парогидротерм - 520 кг/с теплоносителя с энтальпией 1,6 МДж/кг, что соответствует температуре воды 341°С на линии насыщения при давлении 14,8 МПа и тепловому потоку 832 МВт.

9. Установлены расчетные распределения температуры, давления, фазового состояния и скоростей фаз флюида в массиве пород Мутновской магмагеотер-мальной системы в четырех вариантах краевых условий, определяющих начальное термодинамическое состояние пород и условия на границах модели, двух вариантах доминирующего типа теплопереноса в породах - конвективного и кондуктивного, а также для двух сценариев теплового режима очага - остывающего и конвектирующего.

Библиография Диссертация по наукам о земле, доктора технических наук, Пашкевич, Роман Игнатьевич, Хабаровск

1. Азиз, X. Математическое моделирование пластовых систем Текст. / Х.Азиз, Э.Сеттари. - М.: Недра,- 1982. - 407 с.

2. Александров, A.A. Таблицы теплофизических свойств воды и водяного пара: Справочник. Рек. Гос. службой стандартных справочных данных ГСССД Р-776-98 Текст. / А.А.Александров, Б.А.Григорьев. М.: Изд-во МЭИ, 2003. -168 с.

3. Алидибиров, М.А. Экспериментальное исследование процесса разрушения пористой среды при декомпрессии: отчет о НИР Текст. / М.А.Алидибиров. -ИВГиГ ДВО РАН: Петропавловск-Камчатский, 1994. 154 с.

4. Баренблатт, Г.И. Движение жидкостей и газов в природных пластах. -Текст. / Г.И.Баренблатт, В.М.Ентов, В.М.Рыжик. М.: Недра, 1984. - 208 с.

5. Берд, Р.Б. Явления переноса Текст. / Р.Б.Берд, В.Стьюарт, Е.Лайтфут. М.: Химия, 1974.-687 с.

6. Богуславский, Э.И. Добыча и использование тепла Земли. Физико-химическая геотехнология Текст. / Э.И.Богуславский, В.Ж.Аренс, Ю.Д.Дядькин. М.: МГГУ, 2001. - С.583-628.

7. Богуславский, Э.И. Технико-экономическая оценка освоения тепловых ресурсов недр Текст. / Э.И.Богуславский. Л.: ЛГИ, 1984. - 168 с.

8. Богуславский, Э.И. Экономико-математическое моделирование геотермальных циркуляционных систем Текст. / Э.И.Богуславский. Л.: ЛГИ, 1981. -104 с.

9. Вакин, Е.А. Термальные поля и горячие источники Мутновского вулканического района Текст. / Е.А.Вакин, И.Т.Кирсанов, Т.П.Кирсанова // Гидротермальные системы и термальные поля Камчатки. ДВНЦ АН СССР: Владивосток, 1976.-С. 85-114.

10. Вакин, Е.А. Активная воронка Мутновского вулкана Текст. / Е.А.Вакин, И.Т.Кирсанов, А.А.Пронин // Бюллетень вулканологических станций. 1966. -№40.-С. 25-35.

11. Вакин, Е.А. Общая характеристика Мутновского месторождения и прогнозная оценка ресурсов Текст. / Е.А.Вакин, Г.Ф.Пилипенко, В.М.Сугробов // Геотермические и геохимические исследования высокотемпературных гидротерм. М.: Наука, 1986. - 207 с.

12. Вакин, Е.А. Мутновский геотермальный район на Камчатке Текст. / Е.А.Вакин, Г.Ф.Пилипенко // Изучение и использование геотермальных ресурсов в вулканических областях. М.: Наука, 1979. - С. 36-46.

13. Вереина, О.Б. Тот^112-моделирование естественного состояния Мутновского геотермального резервуара Электронный ресурс. / О.Б.Вереина // Материалы Международного геотермального семинара. Россия, Камчатка, 9-15 августа 2004. 02 09.pdf, CD-ROM.

14. Гольдштик, М.А. Вязкие течения с парадоксальными свойствами Текст. / М.А.Гольдштик, В.Н.Штерн, Н.И.Яворский. Новосибирск: Наука, 1989. -336 с.

15. Гультен Э. Некоторые географические заметки к карте Ю. Камчатки Текст. / Э.Гультен // Изв. Русск. геогр. об-ва. т. 57, СПб. - 1925.

16. Гуревич, А.Е. Теоретические основы нефтяной гидрогеологии Текст. / А.Е.Гуревич, Л.И.Капченко, И.М.Кругликов. Л.: Наука, 1972. - 272 с.

17. Дитмар, фон К. Поездка и пребывание в Камчатке в 1851-1855 гг. т. 1 Текст. / К.Дитмар фон. СПб. -1901. — 156 с.

18. Дитмар, К. Поездки и пребывание в Камчатке в 1851-1854 гг. Текст. / К.Дитмар фон // Исторический очерк по путевым дневникам. — СПб, 1901.

19. Дядькин, Ю.Д. Тепловые процессы в горных выработках Текст. / Ю.Д.Дядькин, Ю.В.Шувалов, С.Г.Гендлер. Л.: ЛГИ, 1978. - 114 с.

20. Дядькин, Ю.Д. Процессы тепломассопереноса при извлечении геотермальной энергии Текст. / Ю.Д.Дядькин, С.Г.Гендлер. Л.: ЛГИ, 1985. - 92 с.276

21. Дядькин, Ю.Д. Геотермальная теплофизика Текст. / Ю.Д.Дядькин, С.Г.Гендлер, Н.Н.Смирнова. -Санкт-Петербург: Наука, 1993.-255 с.

22. Дядькин, Ю.Д. Разработка геотермальных месторождений Текст. / Ю.Д.Дядькин. М.: Недра, 1989. - 228 с.

23. Жатнуев, Н.С. Гидротермальные системы с паровыми резервуарами Текст. / Н.С.Жатнуев, А.Г.Миронов, С.Н.Рычагов, В.И.Гунин. — Новосибирск: СО РАН, 1996.-183 с.

24. Злобин, Т.К. Вулканы и вулканические извержения Текст. / Т.К.Злобин // Курильские острова (природа, геология, землетрясения, вулканы, история, экономика). Сахалинское книжн. изд-во: Южно-Сахалинск. - 2004.

25. Зеленский, М.Е. Извержение вулкана Мутновский (Камчатка) 17 марта 2002 г. Текст. / М.Е.Зеленский, А.А.Овсянников, Г.М.Гавриленко, С.Л.Сенюков // Вулканология и сейсмология. 2002. - №6. - С.25-28.

26. Келль, Н.Г. Карта вулканов Камчатки Текст. / Н.Г.Келль. СПб, 1928.

27. Кирсанов, И.Т. Состояние фумарол вулканов Мутновского и Горелого в период с октября 1959 по октябрь 1960 гг. Текст. / И.Т.Кирсанов // Бюллетень вулканологических станций. 1964. - №35. - С.34-43.

28. Кирюхин, A.B. Моделирование эксплуатации геотермальных месторождений Текст. / А.В.Кирюхин. Владивосток: Дальнаука, 2002. - 216 с.

29. Кирюхин, A.B. Модели теплопереноса в гидротермальных системах Камчатки Текст. /А.В.Кирюхин, В.М.Сугробов -М.: Наука, 1987. 152 с.

30. Конради, С.А. Геологический отдел Камчатской экспедиции 1908-1911 гг. Текст. / С.А.Конради, Н.Г.Келль // Изв. Русск. геогр. об-ва. 1925. - LVII. -вып. 1.

31. Корн, Г. Справочник по математике для научных работников и инженеров Текст. / Г.Корн, Т.Корн. М.: Наука, 1984. - 832 с.

32. Кочин, Н.Е. Векторное исчисление и начала тензорного исчисления Текст. / Н.Е.Кочин. М.: Наука, 1965. - 424 с.

33. Лыков, A.B. Тепломассообмен. Справочник Текст. / А.В.Лыков. М.: Энергия, 1972.-560 с.

34. Лыков, A.B. Тепломассообмен. Справочник Текст. / А.В.Лыков. М.: Энергия, 1978.-480 с.

35. Магма малоглубинных камер Текст. АН СССР СО ОТКЗ ИВ. - М.: Наука, 1970.-199 с.

36. Маделунг, Э. Математический аппарат физики Текст. / Э.Маделунг. М.: ГИФМЛ, 1960.-620 с.

37. Мальцева, К.И. Отчет по объекту: пересчет эксплуатационных запасов теплоносителя по Мутновскому месторождению парогидротерм (по состоянию на 31.12.2006). Текст. / К.И.Мальцева, Р.И.Пашкевич. ООО «Аква»: Елизово.-2007.

38. Маренина, Т.Ю. Геолого-петрографический очерк Мутновского вулкана Текст. / Т.Ю.Маренина // Тр. лаб. вулканологии АН СССР. 1956. - Вып.12. С. 1-52.

39. Можен, Ж. Механика электромагнитных сплошных сред Текст. / Ж.Можен. -М.: Мир.-1991.-560 с.

40. Муравьев, A.B. Повторная оценка тепловой мощности фумарольной деятельности на вулкане Мутновский (Камчатка) Текст. / А.В.Муравьев,278

41. Б.Г.Поляк, В.П.Турков, С.В.Козловцева // Вулканология и сейсмология. -1983.-№5. С.51-64.

42. Нигматулин, Р.И. Динамика многофазных сред : в 2-х т. Текст. / Р.И.Нигматулин. М.: Наука, 1987. - 464 с.

43. Никольский, A.JI. О волнах разрушения газированных пород Текст. / А.Л.Никольский // ДАН СССР, 1954. т. 91. - №5. - С. 1035-1038.

44. Никольский, A.A. О волнах внезапного выброса газированных пород Текст. / А.А.Никольский // ДАН СССР, 1953. т. LXXXVIII. - №4. - С.623-626.

45. Пашкевич, Р.И. Гидростатические парадоксы в геотермальный коллекторах Текст. / Р.И.Пашкевич // Вулканология и сейсмология. 1996. - №3. - С. 101105.

46. Пашкевич, Р.И. Теплоперенос в горных породах магмагеотермальной системы Кудрявая (Курильские острова) Текст. / Р.И.Пашкевич // Горная промышленность. 2009. - №3. - С.60-61.

47. Пашкевич, Р.И. Численное моделирование фильтрации при быстрой декомпрессии неподвижной газонасыщенной пористой среды Текст. / Р.И.Пашкевич II Газовая промышленность. 2009. - №7. - С.40-42.

48. Пашкевич, Р.И. Термогидродинамическое моделирование магма-геотермальной системы вулкана Мутновский Текст. / Р.И.Пашкевич, В.В.Таскин // Геоинформатика. 2007. - № 3. - С.56-61.

49. Пашкевич, Р.И. Численное исследование динамики паровых зон магма-геотермальных систем на примере вулкана Мутновский, Южная Камчатка Текст. / Р.И.Пашкевич, В.В.Таскин // Геофизика. 2007. - №5. - С.68-71.

50. Пашкевич, Р.И. Термогидрогазодинамические процессы в массиве пород вулкана Мутновский при остывании его магматического очага: численный эксперимент Текст. / Р.И.Пашкевич, В.В.Таскин // Естественные и технические науки. 2007. - № 5. - С. 169-173.

51. Поляк, Б.Г. Геотермические особенности области современного вулканизма Текст. / Б.Г. Поляк. М.: Наука, 1966. - 180 с.

52. Потапов, В.В. Физическая модель тепломассопереноса в магматогенной геотермальной системе под вулканом Мутновский Текст. /В.В. Потапов // Вулканология и сейсмология. 2002. - №2. - С. 21-30.

53. Селянгин, О.Б. Петрогенезис базальт-дацитовой серии в связи с эволюцией вулканоструктур Текст. / О.Б. Селянгин М.: Наука, 1987. - 148 с.

54. Селянгин, О.Б. К вулканам Мутновский и Горелый: вулканологический и туристический путеводитель Текст. / О.Б. Селянгин — Петропавловск-Камчатский: Новая книга. 2009. - 108 с.

55. Селянгин, О.Б. Новое о вулкане Мутновский: строение, развитие, прогноз Текст. / О.Б. Селянгин // Вулканология и сейсмология. 1993. - №1. - С.17-35.

56. Слеттери, Дж. Теория переноса импульса, энергии и массы в сплошных средах Текст. / Дж. Слеттери М.: Энергия, 1978. - 448 с.

57. Смирнов, В.И. Курс высшей математики для техников и физиков Текст. / В.И. Смирнов-Л.: ГНТИ, 1931.-519 с.

58. Таран, Ю.А. Геохимические исследования в кратере вулкана Мутновский (Камчатка) Текст. / Ю.А. Таран, В.П. Пилипенко, A.M. Рожков, Е.А. Вакин // Вулканология и сейсмология. 1991. - №5. - С. 37-55.

59. Таскин, В.В., Пашкевич, Р.И. Исследование термогидродинамического режима надкритической геотермальной циркуляционной системы Текст. / В.В. Таскин, Р.И. Пашкевич // Известия высших учебных заведений. Горный журнал. 2008. - №8. - С. 21-24.

60. Трухин, Ю.П. Некоторые закономерности современного гидротермального процесса Текст. / Ю.П. Трухин, В.В. Петрова М.: Наука, 1976. - 178 с.

61. Трухин, Ю.П. Ртуть в современном гидротермальном процессе Текст. /Ю.П. Трухин, И.И. Степанов, P.A. Шувалов М.: Наука, 1986. - 199 с.

62. Трухин, Ю.П. Геохимия современных геотермальных процессов и перспективные геотехнологии Текст. / Ю.П. Трухин М.: Наука, 2003. 375 с.

63. Трухин, Ю.П. Вероятный источник флюидов, формирующих современные гидротермы Камчатки Текст. / Ю.П.Трухин, P.A. Шувалов // Флюидный режим земной коры и верхней мантии. Иркутск, 1977. - С. 1—44.

64. Уткин, И.С. Динамика роста и развития проточных магматических очагов

65. Мутновско-Гореловской группы вулканов, их тепловые поля и накопленное281ими подземное тепло Текст. / И.С. Уткин, С.А. Федотов, И.Ф. Делемень, Л.И. Уткина // Вулканология и сейсмология. 2005. - №6. - С. 11-30.

66. Уткин, И.С. Об эволюции и размерах магматических очагов вулканов Текст. / И.С. Уткин, С.А. Федотов, И.Ф. Делемень, Л.И. Уткина // Вулканология и сейсмология. 1999. - №3. - С. 7-18.

67. Федотов, С.А. Магматические питающие системы и механизм извержений вулканов / С.А. Федотов Текст. М.: Наука, 2006. - 247 с.

68. Флетчер, К. Вычислительные методы в механике жидкостей Текст.: в 2 т. / К. Флетчер М.: Мир, 1991.-552 с.

69. Христианович, С.А. О волне дробления Текст. / С.А. Христианович // Известия АН СССР. ОТН. 1953. - №12. - С.1689-1699.

70. Христианович, С.А. Свободное течение грунтовой массы, вызванное расширением содержащегося в порах газа высокого давления Текст. / С.А. Христианович // Препринт. 1979. - №128. - М., ИПМ АН СССР. - 62 с.

71. Христианович, С.А. Неустановившееся течение грунтовой массы содержащей в порах газ высокого давления Текст. / С.А. Христианович // Физико-технические проблемы разработки полезных ископаемых. 1982. - № 3. — С.5-19.

72. Щербань, А.Н. Системы извлечения тепла земной коры и методы их расчета Текст. / А.Н. Щербань, А.С. Цырульников, Э.И. Мерзляков, И.А. Рыженко-Киев: Наукова думка, 1986. 236 с.

73. Albertsson, A. Fluid Handling and Evaluation Электронный ресурс./ A. Al-bertsson, J. Bjarnason, T. Gunnarsson, C. Ballzus, K. Ingason // IDDP feasibility report, Part III, 2003. 32 P., CD-ROM.

74. Ambusso, W. Steam-Water Relative Permeability Текст. / W. Ambusso, C. Satik, R.N. Home // Geoth. Res. Coun. Trans., 1996. 20. - pp 783-795.

75. Anderson, T.B. A fluid mechanical description of fluidized beds Текст. / T.B. Anderson, R. Jackson // Ind. Engng. Chem. Fundam. 1967. - 6(4). - pp. 527-539.

76. Aris, R. Vectors, tensors, and the basic equations of fluid mechanics Текст. / R. Aris -Prentis-Hall: Englewood Cliffs, N.J., 1962.

77. Asanuma, H. Microseismic Monitoring of Hydraulic Stimulation at the Australian HDR Project in Cooper Basin Электронный ресурс. / H. Asanuma, N. Soma, H. Kaieda, D. Wyborn // Proceedings of the World Geothermal Congress. Turkey, 2005. pdf.l615, CD-ROM.

78. Asaulov, S.G. A conceptual model and reservoir assessment for the Mutnovsky geothermal field, Kamchatka, Russia Текст./ S.G. Asaulov // Geothermal Training in Iceland 1994. UNU G.T.P., Iceland, 1994. Report 1. - pp. 1-30.

79. Bachmat, Y. Spatial macroscopization of processes in heterogeneous systems Текст. / Y. Bachmat//Israel J. Technol., 1972. 10(5). pp. 391—403.

80. Baria, R. Creation and mapping of 5000 m deep HDR/HFR reservoir to produceelectricity Электронный ресурс. / Baria R., Michelet S., Baumgartner J., Dyer В.,

81. Nicholls J., Hettkamp Т., Teza D., Soma N., Asanuma H., Garnish J., Megel T. //283

82. Proceedings of the World Geothermal Congress. Turkey, 2005. pdf.l627, CD-ROM.

83. Batini, F. Geophysical well logging a contribution to the fractures characterization Текст. / F. Batini, R. Bertani, B. Ciulli, A. Fiordelisi, P. Valenti // Proceedings, 27-th Workshop on Geothermal Reservoir Engineering. Stanford, 2002. - pp. 371-378.

84. Baumgartner, J. Electricity production from hot rocks Электронный ресурс. / J. Baumgartner, D. Teza, T. Hettkamp, G. Homeier, R. Baria, S. Michelet // Proceedings of the World Geothermal Congress. Turkey, 2005. -pdf.l624, CD-ROM.

85. Bear, J. Physical Principles of Water Percolation and Seepage Текст. / J. Bear, D. Zaslavsky, S.Irmay-UNESCO, Paris, 1968.

86. Bertani, R. High temperature adsorption measurements Текст. / R. Bertani, L. Pa-risi, R. Perini, B. Tarquini // Proc. Stanford Geoth. Workshop, 1996. 21. - pp 523-529.

87. Boehm, R.F. Modeling of a magma energy geothermal open cycle power plant / R.F. Boehm, J. Berg, A. Ortega Текст. // J. of Energy Res. Technology, 1989. V. Ш. - №. 1. pp. 239-245.

88. Bjornsson, G. A Multi-Feedzone Geothermal Wellbore Simulator Текст. / G. Bjornsson/ Lawrence Berkeley Laboratory - report LBL-23546, 1987.

89. Bjornsson, G. The wellbore simulator HOLA. Version 3.1. User's guide Текст. /G. Bjornsson, P. Arason, G.S. Bodvarsson National energy authority: Reykjavik, Iceland, 1993.

90. Brikowski, Т.Н. Modeling supercritical systems with TOUGH2: preliminary results using EOS ISC equation of state module Текст. / Т.Н. Brikowski // Proceedings, 26-th Workshop on Geothermal Reservoir Engineering. Stanford, 2001. pp. 208-215.

91. Brikowski, Т.Н. Modeling supercritical systems with Tough2: the EOSlsc equation of state module and a basin and range example Текст. / Т.Н. Brikowski //Geothermal Resources Council Transactions, 2001. 25. - pp. 285-289.

92. Brown, D. Fluid circulation and heat extraction from engineered geothermal reservoirs Текст. / D. Brown, R. DuTeaux, P. Kruger, D. Swenson, T. Yamaguchi // Geothermics. 1999. - 28. - pp. 553-572.

93. Cathles, L.M. An analysis of the cooling of intrusive by ground water convection which includes boiling Текст. / L.M. Cathles // Econ. Geol, 1977. -72. №5. -pp.804-826.

94. Cheng, P. Heat transfer in geothermal systems Текст. / P. Cheng //Adv. Heat Transfer, 1979.-14:1-105.

95. Collins, R.E. Flow of Fluids through Porous Materials Текст. / R.E. Collins -Reinhold: New York, 1964.

96. Coip, J.L. Utilization of magma energy-project summary. Energy resources of the pacific region Текст. / J.L. Colp, H.M. Stoller // AAPG Studies in Geology № 12, ed. Halbouty M.T. the American association of petroleum geologists, 1981. pp. 541-551.

97. Colp, J.L. Final report magma energy research project Текст. / J.L. Colp // Sandia National laboratories, Sand 2-2377, Albuquerque, NM, 1982.

98. Crapiste, G. H. A general closure scheme for the method of volume averaging Текст. / G. H. Crapiste, E.Rotstein, S. Whitaker // Chem. Engng. Sci., 1986. -41(2).-pp. 227-235.

99. Chu, T.Y. The magma energy program Текст. / T.Y. Chu, J.L. Dunn, J.T. Finger, J.D. Roudle, H.R. Westrich // Geothermal Res. Council Bulletin, 1990. 19. - №2. pp.42-52.

100. Cushman, J. H. Proofs of the volume averaging theorems for multiphase flow Текст. / J.H. Cushman // Adv. Water Res., 1982. 5(4). - pp. 248-253.

101. Cushman, J.H. Multiphase transport equations: 1. General equation for macroscopic statistical, local, space-time homogeneity Текст. / J.H. Cushman // Trans. Theory. Stat. Phys., 1983. - 12(1). - pp. 35-71.

102. Drew, D.A. Averaged field equations for two-phase media Текст. / D.A. Drew //Stud. appl. Math., 1971.- 50(2). pp. 133-166.

103. Drew, D.A. Averaged equations for two-phase flows Текст. // D.A. Drew, L. Segel // Stud. appl. Math., 1971. 50(3). - pp. 205-231.

104. Duchane, D. Hot dry rock geothermal energy in the USA moving toward practical use Текст. / D. Duchane // Proceedings of the World Geothermal Congress -Florence, 1995.-4.-pp. 2613-2617.

105. Dunn, J.C. Magma energy extraction-annual report for FY88 Текст. / J.C. Dunn // Sandia National laboratories, 1989. Sand 89-0567 - Albuquerque, NM.

106. Dunn, J.C. Magma energy extraction Текст. / J.C. Dunn, A. Ortega, C.E. Hicox, T.Y. Chu, R.P. Wemple, R.F. Boehm // Transaction of 12-th Workshop on Geothermal Reservoir Engineering. Stanford, 1987.— pp. 13-20.

107. Dunn, J.C. Energy extraction from crustal magma bodies Текст. / J.C. Dunn // ASME/JSME Thermal Engineering Joints Conference Proceeding. -II. pp. 93100.

108. Edlefsen, N.E. Thermodynamics of soil moisture Текст. / N.E. Edlefsen, A.B. Anderson // Hilgardia, 1943. 15. - pp. 31-298.

109. Ehara, S. Thermal structure beneath Kuju volcano and heat extraction from Kuju-Iwoyama solfatara field Текст. / S. Ehara // Ext. abstr. of workshop on deep-seated geoth. syst., 1994. March 8-10. - Tsukuba. Japan. - pp. 227-235.

110. Elders, W.A. Iceland deep drilling project update Текст. / W.A. Elders // Drilling observation and sampling of the Earth's continental crust, (DOSECC), 2006. — 1. -pp. 10-11.

111. Faust, C.R. Geothermal Reservoir Simulation: 1. Matematical models for liquid and vapour dominated hydrothermal system Текст. / C.R. Faust, J.W. Mercer // Water Resources Research, 1979. 15. -№ 1. -pp. 23-30.

112. Ferrand, L.A. A comparison of capillary pressure-saturation relations for drainage in two and three-fluid porous media Текст. / L.A. Ferrand, P.C.D. Milly, G.E. Pinder, R.P. Turrin //Adv. Water Resour, 1990. 13. - № 2. - pp. 54-63.

113. Fridleifsson, G.O. Progress of the Iceland deep drilling project: march 2003 Текст. / G.O. Fridleifsson, W.A. Elders, S. Saito // ICDP Newsletter. 2006. - 5. -pp. 8-12.

114. Fridleifsson, G.O. The Iceland deep drilling project: a search for deep unconventional geothermal resources Текст. / G.O. Fridleifsson, W.A. Elders // Geother-mics. 2005. - 34. - pp. 269-285.

115. Fujimitsu, Y. Hydrothermal system after the 1990-95 eruption near the lava dome of Unzen volcano, Japan Электронный ресурс. / Y. Fujimitsu, R. Kanou, J. Nishijima, S. Ehara // Proceedings of the World Geothermal Congress Turkey, 2005, CD-ROM

116. Gerlah, T.M. Fuels from magma-potential energy resources? Energy resources of the pacific region Текст. / T.M. Gerlah // AAPG Studies in Geology №12, ed. Halbouty M.T. The American association of petroleum geologists, 1981. pp. 553556.

117. Gianelli, G. Contact metamorphism in the Larderello geothermal system Текст. / G. Gianelli, G. Rugieri // Proceedings of the World Geothermal Congress 2000. Kyushu-Tohoku, 2000. pp. 1163-1168.

118. Gray, W.G. Paradoxes and realities in unsaturated flow theory Текст. / W.G. Gray, S.M. Hassanizadeh // Water Resour. Res., 1991. 27. - №8. - pp. 18471854.

119. Gray, W.G. Unsaturated flow theory including interfacial phenomena Текст. / W.G. Gray, S.M. Hassanizadeh // Water Resour. Res., 1991. 27. - № 8. - pp. 1855-1863.

120. Gray, W.G. A derivation of the equations for multi-phase transport Текст./ W.G. Gray // Chem. Engng. Sci. 1975. - 30. - pp. 229-233.

121. Gray, W.G. On the theorems for local volume averaging of multiphase systems Текст./ W.G. Gray, P.C.Y. Lee // Int. J. Multiphase Flow. 1977. - 3. - pp. 333340.

122. Gray, W.G. Local volume averaging of multiphase systems using a non-constant averaging volume Текст./ W.G. Gray / Int. J. Multiphase Flow. 1983. - 9(6). -pp. 755-761.

123. Gray, W.G. Averaging theorems and averaged equations for transport of interface properties in multiphase systems Текст. / W.G. Gray, S.M. Hassanizadeh // Int. J. Multiphase Flow. 1989. - 15(1). - pp. 81-95.

124. Gray, W.G. Boundary and interface conditions in porous media Текст. / W.G. Gray, S.M. Hassanizadeh // Water Reservoir Research, 1989.

125. Gray, W.G. Mathematical tools for changing spatial scales in the analysis of physical systems: / W.G. Gray, A. Leijnse, R.L. Kolar, C.A. Blain. Boca Raton, CRC Press. - 1993.

126. Gray, W.G. Thermodynamically constrained averaging theory approach for modeling flow and transport phenomena in porous medium systems: 3. Single-fluid-phase flow // W.G. Gray, S.M. Hassanizadeh // Water Resour. Res., 1991. 27(8). - pp. 1855-1863.

127. Gruszkiewicz, M.S. Measurements of water adsorption on the Geysers rocks Текст. / M.S. Gruszkiewicz, J. Horita, J.M. Simonson, R. Mesmer // Proc. Stanford Geoth. Workshop, 1996. 21. - 481-487.

128. Gulick, V.C. Some ground water considerations regarding the formation of small Martian gullies Текст. / V.C. Gulick // Abstracts of Papers Submitted to the Lunar and Planetary Science Conference, 2001. 32. - abstr. no. 2193.

129. Haar, L. NBS/NRC steam tables: thermodynamic and transport properties and computer programs for vapor and liquid states of water in SI units Текст. / L. Haar, J.S. Gallagher, G.S. Kell : New York, Hemisphere Publishing Corp., 1984. -320 P.

130. Hanao, M. Review of recent development of the Kakkonda deep reservoir, Japan Текст. / M. Hanao, M. Takanihashi // Proceedings 18-th Workshop on Geothermal Reservoir Engineering. Stanford, 1993. pp. 26-34.

131. Hassanizadeh, S.M. General conservation equations for multiphase systems: 1. Averaging procedure Текст. / S.M. Hassanizadeh, W.G.Gray // Adv. Wat. Re-sours, 1979.-2(3).- 131-144.

132. Hassanizadeh, S.M. General conservation equations for multiphase systems: 2. Mass, momentum, energy, and entropy equations Текст. / S.M. Hassanizadeh, W.G.Gray // Adv. Wat. Resours. 1979. - 2. - pp. 191-208.

133. Hassanizadeh, S.M. General Conservation Equations for Multiphase Systems: 1. Averaging Procedure Текст. /S.M. Hassanizadeh, W.G.Gray //Flow Through Porous Media. Recent Dev. Ed. G.Pinder, 1983. Southampton, - pp. 1-14.

134. Hassanizadeh, S.M. General conservation equations for multiphase systems: 3. Constitutive theory for porous media Текст. /S.M. Hassanizadeh, W.G.Gray //Adv. Wat. Resours., 1980. 3. - pp. 25-40.

135. Hayba, D.O. Multiphase groundwater flow near cooling plutons Текст. / D.O. Hayba, S.E. Ingebritsen // J. Geophys. Res. 1997. - 102. - pp. 12235-12252.

136. Hicox, C.E. Preliminary considerations for extraction of thermal energy from magma Текст./ C.E. Hicox, J.C. Dunn // Geothermal Resource Council Transactions, 1985. 9. - №2. - pp. 319-324.

137. Hiroyuki, A. Present status and remaining problems of HDR/HWR system design Текст./ A. Hiroyuki, D.V. Duchane, R.H. Parker, M. Kuriyagawa // Geothermics. -1999.-28.-pp. 573-590.

138. Hori, Y. Ogachi project with multi-layer fracturing method for HDR geothermal power outline and future plan Текст. / Y. Hori, H. Kaieda, K. Kitano // Proceedings of the World Geothermal Congress. Florence, 1995. - 4. - pp. 2691-2694.

139. Howes, F.A. The spatial averaging theorem revisited Текст. / F.A. Howes,

140. Whitaker // Chem. Engng. Sci. 40(8). - pp. 1387-1392.

141. Hurwitz, S. Groundwater flow, heat transport, and water-table position within volcanic edifices: Implications for volcanic processes in the Cascade Range Текст. /

142. Ishii, M. Thermo-fluid dynamic theory of two-phase flow Текст. / M. Ishii.- Ey-rolles: Paris. 1975.

143. Jupp, T. A thermodynamic explanation for black smoker temperatures Текст. / Т. Jupp, A. Schultz //Nature. 2000. - 403. - pp. 880-883.

144. Kaieda, H. Review of the Ogachi HDR project in Japan Электронный ресурс. / H. Kaieda, H. Ito, К. Kiho, К. Suzuki, H. Suenaga, К. Shin // Proceedings of the World Geothermal Congress. Turkey, 2005. 1601.pdf, CD-ROM.

145. Kim, Е.К. The economic value of sustainable development as applied to enhanced geothermal power generation : Master of Science in civil and environmental engineering Thesis Текст. / Е.К. Kim Massachusetts Institute of Technology, 2001. -81 P.

146. Kipp, K.L. Guide to the revised ground-water flow and heat transport simulator: HYDROTHERM Version 3 Текст. / K.L. Kipp, Jr., P.A. Hsieh, S.R. Charlton //U.S. Geological Survey Techniques and Methods 6-A25, 2008.

147. Kissling, W.M. Modeling of cooling plutons in the Taupo volcanic zone, New Zealand Текст. / W.M. Kissling // Proceedings, 24-th Workshop on Geothermal Reservoir Engineering. Stanford, 1999.

148. Kissling, W.M. Supercritical TOUGH2 code description and validation Текст. // W.M. Kissling, S.P. White - Industrial Research Limited Report 905, 1999.

149. Kitsou, O.I. Economic modeling of HDR enhanced geothermal systems Текст. / O.I. Kitsou, H.J. Herzog, J.W. Tester // Proceedings of the World Geothermal Congress 2000. Kyushu-Tohoku, 2000. pp. 3779-3784.

150. Lazzarotto, A. Technological developments in deep drilling in the Larderello area Текст. / A. Lazzarotto, F. Sabatelli [Электронный ресурс] // Proceedings of the World Geothermal Congress. Turkey, 2005. 1018.pdf, CD-ROM.

151. Ledru, P. Enhanced geothermal innovative network for Europe Текст./ P. Ledru, A. Genter //Proceedings European Geothermal Congress 2007. Unterhaching, Germany, 30 May-1 June 2007.

152. Manning, C.E. Permeability of the continental crust: implications of geothermal data and metamorphic systems Текст./ C.E. Manning, S.E. Ingebritsen // Reviews of Geophysics. 1999. - 37, l/February. - pp. 127-150.

153. Matsushima, N. Mathematical simulation of magma-hydrothermal activity associated with the 1977 eruption of Usu volcano Текст./ N. Matsushima // Earth Planets Space, 2003. 55. - pp. 559-568.

154. Matsunaga, I. Review of the HDR development at Hijiori site, Japan Электронный ресурс. / I. Matsunaga, H. Niitsuma, Y. Oikawa // Proceedings of the World Geothermal Congress. Turkey, 2005. 1635.pdf, CD-ROM.

155. Mercer, J.W. Review of simulation techniques for aquifer thermal energy storage (ATES) / J.W. Mercer, C.R. Faust, W. J. Miller, F.J. Pearson // Advances in hydroscience.- 1982.- 13.-pp. 1-129.

156. Moeckel, J.P. Thermodinamics of an interface / J.P. Moeckel //Archs Ration. Mech. Analysis. 1975. - 57(3). - pp. 255-280.

157. Moench, A.F. Simulation of steam transport in vapor-dominated geothermal reservoirs Текст. / A.F .Moench // Open File rep. 76-607, 43 pp., U.S. Geol. Surv., Menlo Park, Calif., 1976.

158. Nakada, S. Unzen Scientific Drilling Project (USDP): hot conduit ahead Текст. / S.Nakada // Newsletter. 2003. - 5. - pp.3-4.

159. Norton, D.L. Preliminary numerical analysis of the magma-hydrothermal history of the Geysers geothermal system, California, USA Текст. / D.Norton, J.B.Hulen // Geothermics. 2001. - V.30. - pp. 617-625.

160. Norton, D. Complex behavior of magma-hydrothermal processes: role of supercritical fluid Текст. / D.Norton, B.Dutrow // Geochimica et Cosmochimica Acta, 2001. V.65. - №. 21. - pp. 4009-4017.

161. Norton, D. Transport phenomena in hydrothermal systems: Cooling plutons Текст. / D. Norton, J.Knight // Amer. J. Sci., 1977. V.277. - pp. 937-981.

162. Norton, D. Sourcelines, sourceregions and pathlines for fluids in hydrothermal systems related to cooling plutons Текст. / D. Norton // Econ. Geol., 1978. 73. - pp. 21-28.

163. Official DOSECC (Drilling observation and sampling of the Earth's continental crust) website. http://www.dosecc.org/html/hawaii.html

164. Okubo, A. Numerical simulation of volcanomagnetic effects due to hydrothermal activity Текст. / A.Okubo, W.Kanda, K.Ishihara // Annuals of Disas. Prev. Res. Inst., Kyoto Univ., 2006. No. 49 C.

165. Pan, L. WinGridder. An interactive grid generator for TOUGH. Version 1.0. (User's Manual) Текст. / L.Pan, J.Hinds, C.Haukwa, Y.-Sh.Wu, G.Bodvarsson. -Lawrence Berkeley Laboratory Report LBNL-42957, 2001.

166. Pan, L. User's manual (UM) for WinGridder V2.0 Текст. / L.Pan Lawrence

167. Berkeley Laboratory Report LBNL-2001-131. 2002.293

168. Parker, R. The Rosemanowes HDR project 1983-1991 Текст. / R.Parker // Geo-thermics. 1999. V.28. -pp.603-615.

169. Persoff, P. Hydrologie Characterization of Four Cores from the Geysers Coring Project Текст. / P.Persoff, J.B.Hulen // Proc. Stanford Geoth. Workshop, 1996. -21. -pp.327-334.

170. Polyansky, O.P. Convection of two-phase fluid in a layered porous medium driven by the heat of magmatic dikes and sills Текст. / O.P.Polyansky, V.V.Reverdatto, A.V.Khomenko, V.G.Sverdlova // Geochemistry International, 2002. 40. Suppl. l.-pp. S69-S81.

171. Pritchett, J.W. STAR: Age of thermal reservoir simulation system Текст. / J.W.Pritchett // Proc. World Geothermal Congress 1995, Florence, Italy, May 1831, 1995. International Geothermal Association. 2959-2963.

172. Pruess, K. Grid Orientation and Capillary Pressure Effects in the Simulation of Water Injection into Depleted Vapor Zones Текст. / K.Pruess // Geothermics. -1991.-20 (5/6). 257-277.

173. Pruess, K. Numerical Simulation of Water Injection into Vapor-Dominated Reservoirs Текст. / K.Pruess // Proc. World Geoth. Cong., Florence, Italy, 1995. 3. -pp. 1673-1679.

174. Pruess, K. Injection Plume Behavior in Fractured, Vapor-Dominated Reservoirs Текст. / K.Pruess // Proc. Stanford Geoth. Workshop, 1996. 21. - pp. 413-^120.

175. Pruess, K. TOUGH2 User's Guide, Version 2.0 Текст. / K.Pruess, C.Oldenburg, G.Moridis // Lawrence Berkeley National Laboratory Report LBNL-43134, Berkeley, CA, November 1999.

176. Rathbun, J.A. Interaction of groundwater with impacts on Mars: Possible hydrothermal systems Текст. / J.A.Rathbun, S.W.Squyres //Abstracts of Papers Submitted to the Lunar and Planetary Science Conference, 2000. -31.- abstr. no. 1111.

177. Reid, M.E. Massive collapse of volcano edifices triggered by hydrothermal pressu-rization Текст. / M.E.Reid // Geology. 2004. - 32. - №5. - pp.373-376.

178. Rose, P. An enhanced geothermal system at Coso, California recent accomplishments Электронный ресурс. / P.Rose, J.Sheridan, J.McCulloch, J.N.Moore 1,

179. K.Kovacl, R.Weidler, S.Hickman // Proceedings of the World Geothermal Congress. Turkey, 2005. 1603.pdf.

180. Sanyal, S.K. An analysis of power generation prospects from enhanced geothermal systems Текст. / S.K.Sanyal, S.J.Butler // Geothermal Resources Council Transactions, 2005.-29.

181. Satik, C. Adsorption Characteristics of Rocks from Vapor-Dominated Geothermal Reservoir at the Geysers Текст. / С.Satik, M.Walters, R.N.Horne // Proc. Stanford Geoth. Workshop, 1996. 21. - 469-479.

182. Satik, C. An Experimental Study of Boiling in Porous Media Текст. / C.Satik, R.N.Horne // Geoth. Res. Coun. Trans., 1996. 20. - pp. 839-843.

183. Sato, Y. Status of Japanese HDR project at Hijiory Текст. / Y.Sato, K.Ishibashi, T.Takada, Yamaguchi // Proceedings of the World Geothermal Congress. Florence, 1995. IV. - pp. 2677-2678.

184. Schroeder, R. On deep seated geothermal resources Текст. / R.Schroeder, S.Takasugi, K.Osato // Proceedings, 23-th Workshop on Geothermal Reservoir Engineering. Stanford, 1998.-pp.451-^159.

185. Sengers, J.V. Representative equations for the viscosity of water substance Текст. / J.V.Sengers, B.Kamgar-Parsi // Journal of Physical and Chemical Reference Data. 1984.- 13. -pp.185-205.

186. Slattery, J.C. Flow of viscoelastic fluid through porous media Текст. /

187. J.C.Slattery // A.I.Ch.E. Jl, 1967. 13(3). - pp. 1067-1071.295

188. Smith, Т. An investigation of boiling processes in hydrothermal eruptions Текст. / T.Smith, R.McKibbin // Proceedings of the World Geothermal Congress 2000. -Kyushu-Tohoku, Japan. May 28-June 10. pp.699-703.

189. Shigeno, H. Evolution history of the Kakkonda magma-hydrothermal system, Japan, estimated though simplified-model numerical simulations Текст. / H.Shigeno // Proceedings of the World Geothermal Congress 2000. Kyushu-Tohoku, 2000. -pp.135-142.

190. Sta.Maria, R.B. Optimization of Water Injection into Vapor-Dominated Geothermal Reservoirs Текст. / R.B.Sta.Maria, R.N.Horne // Geoth. Res. Coun. Trans., 1996.-20.-859-869.

191. Sta.Maria, R.B. Simulating the Effects of Adsorption and Capillary Forces in Geothermal Reservoirs Текст. / R.B.Sta.Maria, A.S.Pingol // Proc. Stanford Geoth. Workshop, 1996.-21.-pp. 165-173.

192. Taran, Yu.A. A geochemical model for fumaroles of the Mutnovsky volcano, Kamchatka, USSR Текст. / Yu.A.Taran, V.P.Pilipenko, A.M.Rozhkov, E.A.Vakin // J. of Volcanol. and Geotherm. Res. 1992. - V.49. - pp. 269-283.

193. Tecplot. User's Manual. Version 9.2, Release 1. //Amtec Engineering, Inc. Belle-vue, Washington, September, 2002.

194. Tecplot. Reference Manual. Version 9.2, Release 1. //Amtec Engineering, Inc. Bel-levue, Washington, September, 2002.

195. Tesla, N. Our Future Motive Power Текст. / N.Tesla // Everyday Sciences and Mechanics. 1931.

196. The Future of Geothermal Energy Impact of Enhanced Geothermal Systems (EGS) on the United States in the 21st Century // MIT-led interdisciplinary panel, Massachusetts Institute of Technology, 2006.

197. Tomiya, A. Depth of magma chamber determined by experimental petrologic methods Текст. / A.Tomiya // Proceedings of the World Geothermal Congress 2000. Kyushu-Tohoku, 2000. pp. 1855-1859.

198. Wallroth, T. Hot dry rock research experiments at Fjallbacka, Sweden Текст. / T.Wallroth, T.Eliasson, U.Sundquist // Geothermics. 1999. - 28. - pp.617-625.

199. Watanabe, К. Numerical study on heat extraction from supercritical geothermal reservoir Текст. / К.Watanabe, Y.Niibori, T.Hashida // Proceedings World Geothermal Congress 2000, Kyushu Tohoku, Japan, May 28 - June 10, 2000. - pp. 3957-3961.

200. White, S.P. Modeling the Deep Geothermal System of the Uenotai Reservoir Электронный ресурс. / S.P.White, T.Okabe, T.Sato, M.Sato, T.Shiga, Y.Takahashi // Proceedings World Geothermal Congress 2005, Antalya, Turkey. -1183.pdf.

201. Whitaker, S. Diffusion and dispersion in porous media Текст. / S.Whitaker // A.I.Ch.E. Jl, 1967. 13(3). - pp. 420-427.

202. Yanagisawa, N. Scale variation of the production wells from deep reservoir in Kakkonda field Текст. / N. Yanagisawa, K.Fujimoto, Y.Hishi // Proceedings, 27-th Workshop on Geothermal Reservoir Engineering. Stanford, 2002. — pp.327-332.

203. Yano, Y. Numerical investigation of production behavior of deep geothermal reservoirs at super-critical conditions Текст. / Y.Yano, T.Ishido // Geothermics. -1998. 27. - pp.705-721.

204. СПИСОК ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

205. Пашкевич Р.И. Гидростатические парадоксы в геотермальных коллекторах Текст. / Р.И.Пашкевич // Вулканология и сейсмология. — 1996. №3. С. 101-106.

206. Pashkevich R.I. On environmental aspects of geothermal development Текст. / RJ.Pashkevich // Geothermal Resources Council TRANSACTIONS, 1996. Vol.20. September/October. — pp.241243.

207. Pashkevich R.I. About convective heat transfer in geothermal systems Текст. / R.I.Pashkevich // Geothermal Resources Council TRANSACTIONS, 1996. Vol. 20. September/October.- pp.829832.

208. Pashkevich R.I. Kamchatka geothermal resources development: problems and perspectives Текст. / RJ.Pashkevich // Proc. Stanford Geoth. Workshop, 1996. 18. —pp. 539-541.

209. Pashkevich R.I. Results and interpretations of hot and cold water injection experiments on Pauz-hetka water-dominated geothermal field in Kamchatka Текст. / RJ.Pashkevich // Proc. Stanford Geoth. Workshop, 1996. 18. - pp.543-547.

210. Pashkevich R.I. On a capillary pressure in two-phase reservoirs Текст. / R.I.Pashkevich // Proceedings 18th New Zealand Workshop, 1996. pp.291-294.

211. Pashkevich R.I. On possibility of HDR project in near-by region of Petropavlovsk-Kamchatsky, Russia Текст. / R.I.Pashkevich // Proceedings 3rd International HDR Forum, May 13-16. -Santa Fe, New Mexico, USA, 1996. pp.144-145.

212. Pashkevich R.I. The contradictions in continua modelling of flow and heat transfer in fractured and/or porous reservoirs Текст. / R.I.Pashkevich // Geothermal Resources Council TRANSACTIONS. Vol.21. September/October. - 1997. - pp.617-620.

213. Pashkevich R.I. Method of capillary effects inclusion in two-phase flow modeling of a geothermal reservoir Текст. / R.I.Pashkevich // Proc. Stanford Geoth. Workshop. 1998. - 20 - pp.420-427.

214. Пашкевич Р.И. Численное исследование динамики паровых зон магмагеотермальных систем на примере вулкана Мутновский, Южная Камчатка Текст. / Р.И.Пашкевич,

215. B.В.Таскин // Геофизика. — 2007. — №5. С.68-71 (соискателем выполнена постановка задачи, анализ результатов, сформулированы выводы).

216. Пашкевич Р.И. Термогидродинамическое моделирование магмагеотермальной системы вулкана Мутновский Текст. / Р.И.Пашкевич, В.В.Таскин // Геоинформатика. 2007. — № 3.

217. C.57-61 (соискателем выполнена постановка задачи, анализ результатов, сформулированы выводы).

218. Пашкевич Р.И. Численное моделирование фильтрации при быстрой декомпрессии неподвижной газонасыщенной пористой среды Текст. / Р.И.Пашкевич // Газовая промышленность. 2009. - №7. - С.40^2.

219. Пашкевич Р.И. Теплоперенос в горных породах магмагеотермальной системы Кудрявая (Курильские острова) Текст. / Р.И.Пашкевич // Горная промышленность. — 2009. №3. -С.60-61.

220. Пашкевич Р.И. Осреднение по межфазной поверхности в фильтрационном потоке Текст. / Р.И.Пашкевич // Геофизика. 2009.- №4. - С.63-65.

Информация о работе

Пашкевич, Роман Игнатьевич
доктора технических наук
Хабаровск, 2009
ВАК 25.00.20

Диссертация

Научно-техническое обоснование рациональных параметров теплопереноса и фильтрации двухфазного теплоносителя при освоении геотермальных месторождений Камчатки - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы