Введение к работе
Актуальность темы.
Общепризнанным подходом к проведению теоретических исследований в прикладных проблемах естествознания является вычислительный эксперимент. Успех вычислительного эксперимента определяют
удачно выбранная математическая модель явления, учитывающая многообразие протекающих процессов, реальную геометрию, физические свойства среды,
использование эффективного численного метода, построенного на основе глубоко проработанной теории, адекватного в определенном смысле выбранной математической модели,
способ реализации алгоритма на ЭВМ.
При этом разработка эффективных численных методов наряду с увеличением мощности ЭВМ расширяет границы применимости вычислительного эксперименета и обусловливает широкое распространение ориентированных на конкретные прикладные области программных средств, реализованных в виде программных комплексов и пакетов прикладных проградім, которые расчитаны на широкий класс пользователей.
Одной из актуальных проблем, решение которой может быть получено средствами вычислительного экспетримента является проблема прогнозирования оптимальных режимов нефтедобычи и определения их техногенных последствий. В настоящее время в технологиях нефтедобычи и других геопользовательских работах математиче-
ское моделирование фильтрации в геологических средах со сложной структурой с учетом напряженно - деформированного состояния во многих случаях носит чисто качественный характер. Это обусловлено тем, что геомеханические модели являются неполными, а также ввиду неточности и субъективности геолого - геофизической информации о пласте. Более того, структурно многомасштабные среды вообще плохо поддаются моделированию. Рассмотрение более полных моделей, учитывающих реальное многообразие процессов и сложную геологическую структуру среды позволит усовершенствовать существующие технологии геопользования и в первую очередь технологии в нефтедобывающей промышленности с помощью численного моделирования процессов нагнетания жидкости в пласты через скважины. Это позволит оптимизировать параметры закачки жидкости с целью увеличения продуктивности пластов или для решения других технологических задач, а также предсказывать и предотвращать ряд негативных явлений. В частности, Нефтегорское землетрясение 27 мая 1995 года «было «подготовлено» технологическим нарушением свойств геологической среды вследствие отработки нефтяного месторождения»1, так как деформации горных пород, вызываемые перемещением флюидных масс, в тектонических нарушениях достигают уровней, которые близки к критическим, то есть к уровням неустойчивой деформации, превышение которых вызвает разрушение пород с высвобождением сейсмической энергии.
1Рудаков В.П. Геодинамические предпосылки Нефтегорского землетрясения 27 мая 1995 года // ДАН. 1995. Т. 345. N б. С. 81S-823
Состояние вопроса.
Задачи, целью решения которых является изучение процессов фильтрации в насыщенных пористых геологических средах, сопровождаемых изменением напряженно - деформированного состояния в этих средах, описываются уравнениями Био [Biot М.А., Rice J.R., Cleary М.Р., Николаевский В.Н. и др.]. Математическая модель Био наряду с уравнениями теории фильтрации — балансным уравнением для массы жидкости и выражением для потока согласно закону Дарси — включает уравнения линейной теории упругости для насыщенных пористых сред. Особенностью рассматриваемого класса задач является необходимость учитывать геологическую неоднородность среды, обусловленную историей развития: она состоит из ряда разновозрастных слоев сложной геометрии, резко отличающихся по своим свойствам, порой на несколько порядков. Таким образом, для изучения процессов фильтрации в насыщенных пористых геологических средах необходимо развитие численных методов на нерегулярных сетках. Под нерегулярными сетками здесь понимаются разностные сетки, которые не могут быть получены из прямоугольной гладким преобразованием координат.
Наряду с другими методами, в частности, методом конечных элементов [Сьярле Ф., Стренг Г., Фикс Дж., Heinrich В.], широко распространенным подходом к численному решению уравнений математической физики на нерегулярных сетках является метод опорных операторов, используемый в настоящей работе для построения дискретных алгоритмов. Метод опорных операторов позволяет построить согласованные разностные аналоги дифференциальных опе-
раторов тензорного анализа и на их основе — разностные схемы для уравнений механики сплошной среды, сохраняющие такие качества уравнений, как консервативность, симметричность (если она была) и положительность эллиптического оператора. Это дает возможность получать физически корректные результаты даже на грубых сетках, а также обеспечивает необходимый порядок аппроксимации в расчетной области в случае негладких коэффицентов и правых частей уравнений. Метод опорных операторов в том виде в котором он используется в настоящей работе, был предложен Самарским А.А., Тишкиным В.Ф., Фаворским А.П., Шашковым М.Ю. в 1981 году. В последующих работах этих и других авторов [Колдоба А.В., Короби-цин В.А., Коршия Т.К., Повешенко Ю.А.] были построены и исследованы разностные схемы для различных уравнений математической физики. Идея метода состоит в том, что один из рассматриваемых дифференциальных операторов первого порядка аппроксимируется непосредственно, а остальные таким образом, чтобы удовлетворить разностным аналогам некоторых интегральных тождеств. В частности, таким тождеством, содержащим тензорные поля, является
где Uj, рік — достаточно гладкие функции в области О. Его дискретные аналоги позволяют конструировать разностные схемы для уравнений линейной стационарной теории упругости. Так как метод опорных операторов не ограничен предположениями о структуре разностной сетки, он является адекватным подходом к построению разностных схем для уравнений линейной стационарной теории упру-
гости на нерегулярных сетках.
Непосредственным обобщением метода опорных операторов на случай учета различных физических нелинейных процессов (например, магнитогидродинамических, гравитационных, вязких и т.д.), разворачивающихся на фоне движения сплошной среды является метод диссипативных функций [Попов Ю.П. и др.]. Кратко остановимся на основных идеях метода. При моделировании механики движения такой «осложненной» сплошной среды обычно существует тензор напряжений, обусловленный тем или иным физическим полем, отвечающий за физический процесс (максвелловский тензор напряжений электромагнитного поля, тензор напряжений гравитационного поля, тензор вязких напряжений и т.д.) и тензор деформации (или скоростей деформации), отвечающий за кинематическое описание движения сплошной среды. Свертка симметризованного тензора деформации и «физического» тензора напряжений — диссипативная функция — определяет внутреннее (то есть не работающее через поверхность тела) преобразование энергии при соответствующем физическом процессе. Успешная пространственная дискретизация такой свертки в значительной мере и определяет эффективность применения метода опорных операторов и метода диссипативных функций, так как дальнейшее вычисление соответствующей силы (Лоренца, тяготения, вязкой и т.д.) в результате сопряженного интегрального преобразования является делом техники. Не такова модель Бно, так как в ней одновременно присутствуют два различных вида движения (жидкости и скелета) и соответственно две различных диссипативных функции единого энергетического преобразования. Поскольку
прямое дискретное моделирование такого дуализма представляется нецелесообразным, то имеет смысл уже на континуальном уровне воспользоваться обобщенным силовым тензором модели Био, отвечающим только за силовые характеристики процесса напряженного деформирования полной системы жидкость - скелет. Обобщенный силовой тензор и нестационарные объемные деформации рассмотренной вариационно - энергетической формулировки модели Био движения жидкости в деформируемом пористом теле легко дискретизируемы на основе метода опорных операторов.
Широкое распространение программных средств, опирающихся на эффективные алгоритмы и ориентированных на конкретные прикладные области возможно лишь при условии их реализации в виде пакетов прикладных программ. Использованные в работе над диссертацией инструментальные средства системы ТЕКОН позволили в значительной мере сократить трудозатраты на программирование алгоритмов. Система ТЕКОН представляет собой комплекс инструментальных средств, на основе которых возможно создание широкого класса пакетов прикладных программ для решения задач математической физики. Средствами ТЕКОНа создаются дискрипторы нерегулярных сеток и организовываются циклические и нециклические операции над ними. При этом используется динамическая организация памяти, которая в отличие от статической состоит в хранении информации в виде сегментов, упорядоченных посредством организации ссылок в список. Такая организация дает возможность возложить на ЭВМ часть громоздких символьных выкладок, позволяет производить декомпозицию матриц произвольной структуры. Инструмен-
тальные средства системы ТЕКОН позволили генерировать версию пакета программ для решения краевых задач теории фильтрации на основе разностных схем метода опорных операторов [Повещен-ко Ю.А., Дремов О.А. и др]. Эта, хорошо зарекомендовавшая себя в процессе решения практических задач, версия пакета послужила основой для создания программного комплекса предназначенного для решения более сложного класса задач — численного моделирования процессов фильтрации в насыщенных геологических средах, сопровождаемых изменением напряженно - деформированного состояния, в рамках вариационно - энергетической формулировки модели Био движения жидкости в пористом теле.
Цель работы. Настоящая диссертация посвящена разработке и реализации алгоритма численного решения задач пороупругости, описываемых уравнениями Био на основе разностных схем метода опорных операторов.
Научная новизна. Разработаны вычислительные алгоритмы для численного интегрирования уравнений механики Био в областях неоднородного геологического строения. Эти алгоритмы, основанные на идеях метода опорных операторов, сохраняют консервативность дифференциальных уравнений, что дает возможность получать физически корректные результаты даже на грубых сетках и в случае негладких коэффицентов и правых частей уравнений. Указанные алгоритмы реализованы в виде комплекса программ для персональных компьютеров. С помощью этого комплекса проведены расчеты в ряде
прикладных задач геофизики.
Апробация результатов диссертации. Основные результаты работы докладывались на научно - исследовательских семинарах ИПМ РАН им. М.В.Келдыша, ИММ РАН, НИИ механики МГУ им. М.В. Ломоносова, ВНИИГЕОСистем, НИИ океанологии РАН, на XIV Губкинских чтениях (ГАНГ им. И.М. Губкина, 15-17 октября 1996г., участие в стендовом докладе).
Публикации. Излагаемый в диссертации материал опубликован в 3 печатных работах, перечисленных в конце автореферата.
Структура И обьем диссертации. Диссертация состоит из введения, трех глав, заключения и списка литературы, включающего 51 наименований.
