Содержание к диссертации
Введение
Глава 1 Течение вязкой жидкости с примесью частиц по раскрывающейся трещине в пористой среде
1.1 Общие положения 19
1.2 Математическая постановка задачи 19
Допущения. Форма трещины. Кинематические соотношения. Уравнения движения.
1.3 Автомодельное движение 24
Осаждение частиц. Прилипание частиц.
1.4 Уравнения движения в безразмерной форме 28
Безразмерные переменные. Малый осадок частиц.
1.5 Остановка роста трещины 30
Постановка задачи. Анализ расчетов.
1.6 Формирование трещины гидроразрыва 34
Уравнения движения. Безразмерная форма. Начальные и краевые условия. Анализ расчетов.
1.7 Развитие трещины в ограниченном объеме 41
Система уравнений движения. Некоторые результаты расчетов.
1.8 Осаждение взвешенных частиц 46
Глава 2 Нестационарные течения при разрыве трубопровода
2.1 Истечение вязкой жидкости из трубы в затопленное пространство 51
2.1.1 Общие положения 51
2.1.2 Математическая постановка задачи 52
Допущения. Модель течения.
2.1.3 Безнапорное истечение 54
Автомодельное движение. Анализ решения.
2.1.4 Напорное истечение 59
Сила сопротивления. Безразмерная форма уравнений движения. Анализ расчетов.
2.1.5 Наклонная труба 63
Уравнения движения. Начало движения. Стабилизация течения.
2.1.6 Перетекание вязких жидкостей в трубе 66
Математическая постановка задачи. Анализ расчетов.
2.2 Растекание вязкой жидкости по плоской поверхности 69
2.2.1 Общие положения 69
2.2.2 Основные уравнения и допущения 69
Модель течения. Граничные режимы.
2.2.3 Предельные режимы протекания в пористую среду 72
Медленное протекание. Быстрое протекание. Свободное протекание.
2.2.4 Автомодельное движение 74
Постоянный расход. Постоянный объем.
2.2.5 Движение по наклонной плоскости 79
Математическая постановка задачи. Движение без протекания в грунт.
2.2.6 Растекание легкой жидкости по нижней поверхности ледяного покрова водоема 86
Постановка задачи. Постоянный расход. Постоянный объем.
Глава 3 Нестационарные непрерывные течения газа
3.1 Общие положения 89
3.2 Течение газа по системе труб при внешних воздействиях 91
Основные уравнения и допущения. Уравнения движения в безразмерных переменных. Сопряжение нестационарных потоков. Разностные уравнения. Политропический процесс сжатия.
3.3 Изотермическое течение газа 100
Уравнения движения. Согласование потоков. Разностные уравнения. Заданное поле температуры. Изотермическая модель компрессора.
3.4 Сравнение методов сопряжения нестационарных потоков 103
Постановка задачи. Анализ расчетов. Сравнение методов.
3.5 Движение газа Дюпре-Гирина 108
Основные уравнения и допущения. Уравнения движения в безразмерных переменных. Условия сопряжения потоков. Изотермическое течение газа. Неизотермическое течение газа.
3.6 Матричное представление ГТС 113
3.7 Начальные и краевые условия 115
3.8 Реализация разностной схемы 119
Выбор коэффициента искусственной вязкости. Определение значений Куранта. Коммерческий расход. Коэффициент сопротивления.
3.9 Движение реального газа 123
3.9.1 Основные уравнения 123
3.9.2 Упрощенные модели 126
3.9.3 Примеры неустановившихся течений 127
3.10 Программа «gas-stream» 131
Глава 4 Течение газа с разрывами
4.1 Общие положения 139
4.2 Изотермические ударные волны 139
4.2.1 Основные уравнения и допущения 139
Уравнения движения. Уравнения движения в безразмерных переменных. Общая постановка задачи.
4.2.2 Распространение граничного режима на полупрямой 144
Особенности безразмерной формы решения. Скачок давления на границе. Скачок расхода на границе.
4.2.3 Распространение граничного режима на отрезке 147
Скачок давления на границе. Скачок расхода на границе. Опорожнение трубки.
4.3 Ударные волны 158
4.3.1 Уравнение энергии 158
Уравнение для полной энергии. Теорема об изменении кинетической энергии. Уравнение для внутренней энергии. Дивергентная форма записи уравнений движения. Уравнения состояния. Уравнения движения в безразмерных переменных.
4.3.2 Распространение ударной волны по трубе 163
Постановка задачи. Сопротивление и теплообмен отсутствуют. Влияние теплообмена. Влияние сопротивления. Влияние сопротивления и теплообмена.
4.3.3 Гидравлический удар 168
Мгновенная остановка потока. Сравнение моделей течения. Влияние определяющих параметров.
4.3.4 Акустическое приближение 173
Линеаризация уравнений. Дисперсионное соотношение. Частные случаи.
4.4 Одномерные изэнтропические течения 176
4.4.1 Периодический граничный режим 176
4.4.2 Взаимодействие ударного импульса с контактным разрывом 181
4.5 Волновые двумерные течения газа 191
4.5.1 Торцевой разрыв трубы 191
Уравнения движения. Математическая постановка задачи. Результаты расчетов.
4.5.2 Разрушение участка трубы 198
Постановка задачи. Анализ расчетов.
Глава 5 Волновые движения газа с примесью частиц
5.1 Общие положения 204
5.2 Ударный выброс газа с примесью частиц 205
Уравнения движения фаз инертной газовзвеси. Постановка задачи. Некоторые результаты расчетов. Анализ результатов расчетов.
5.3 Отклонение потока газовзвеси в косой ударной волне 226
Соотношения на ударных скачках. Влияние объемного содержания частиц на параметры газовзвеси в равновесной области двухфазного течения за косой ударной волной. Влияние объемного содержания частиц на величину досжатия газа при нормальном отражении от стенки плоской ударной волны.
5.4 Модель решетчатого газа как система взаимодействующих частиц 240
Уравнения для корреляционных функций. Полубесконечная система частиц. Бесконечная система частиц во внешнем поле. Сглаживание осцилляции плотности.
5.5 Истечение в атмосферу реагирующей взвеси 251
Уравнения движения фаз реагирующей взвеси. Постановка задачи. Результаты расчетов.
Заключение 261
Литература 264
- Математическая постановка задачи
- Истечение вязкой жидкости из трубы в затопленное пространство
- Течение газа по системе труб при внешних воздействиях
- Изотермические ударные волны
- Ударный выброс газа с примесью частиц
Введение к работе
В работе исследуются нестационарные течения жидкостей, газов и смесей применительно к проблемам нефтегазового комплекса. Исследование выполнено на основе построения математических моделей с привлечением идей и подходов механики сплошных сред и вычислительной математики.
Задачи о неустановившихся течениях углеводородов привлекают внимание многих исследователей богатством затрагиваемых явлений и процессов, а также своей прикладной значимостью. В практике нефтедобычи одним из самых ярких гидродинамических процессов является гидроразрыв пласта. Применение гидроразрыва пласта началось в шестидесятых годах. К этому же времени относятся и первые теоретические исследования развития и формирования трещины гидроразрыва. С введением в разработку крупных месторождений Западной Сибири интерес к гидроразрыву пласта несколько снизился. В последние годы произошло существенное изменение структуры запасов углеводородов. Многие крупные месторождения в настоящее время находятся на вторичной стадии разработки, основные запасы нефти сосредоточены в низкопроницаемых пластах. Интерес к гидроразрыву вновь возрос.
На протяжении всей истории использования гидроразрыва нефтяного пласта ведутся экспериментальные и теоретические исследования. В результате совершенствуется технология проведения гидроразрыва и растет его эффективность. Ниже приведен краткий обзор теоретических работ зарубежных и отечественных авторов.
Одно из первых аналитических решений в существенно упрощенной постановке задачи о гидроразрыве пласта найдено R.Carter [162]. Основополагающими по теории гидроразрыва являются работы T.K.Perkins, L.R.Kern [176] и R.P.Nordgren [173]. Большинство предложенных моделей развития трещины обладают важным свойством существования автомодельных решений. Общий способ конструирования моделей гидроразрыва и поиска
7 определенных классов автомодельных решений дан M.Cleary et al [164, 165]. В работах S.Advani et al [155, 156, 157] найдены автомодельные решения для широкого класса моделей.
В отечественной литературе основополагающей в развитии теории гидроразрыва является работа Ю.П.Желтова, С.А.Христиановича [54]. Свой вклад внесли Г.И.Баренблатт [14, 15], Ю.В.Желтов (Ю.П.Желтов, Ю.В.Желтов [53]). Одна из последних работ по теоретическому моделированию гидроразрыва в русле идей С.А.Христиановича и Ю.П.Желтова принадлежит О.П.Алексеенко [3]. По теории и практике гидроразрыва имеется ряд монографий Ю.П.Желтова [52], А.С.Меликбекова [96], В.Г.Логинова, В.А.Блаженовича [89]. Распространение вертикальной трещины гидроразрыва в непроницаемой среде и с учетом свойств среды фильтровать жидкость рассмотрены Ю.Н.Гордеевым, А.Ф.Зазовским [40, 41]. О.Е.Ивашневым, Н.Н.Смирновым [59] описан процесс формирования трещины гидроразрыва в пористой среде. Рассматриваются вопросы и о притоке жидкости к скважине с трещиной (Р.Д.Раневская, Р.М.Кац [112]; M.J.Economides, K.G.Nolte [166] и
ДР-)-
Вторая часть работы посвящена изучению нестационарных процессов, связанных с транспортированием газа в разветвленных системах трубопроводов.
Современная технология транспорта газа превращает магистральный трубопровод в единую динамическую систему, требующую согласованной работы всех ее элементов. Увеличение протяженности и мощности трубопроводов вместе с усложнением его структуры предполагает развитую систему его математического обеспечения. Математическое моделирование неустановившихся течений в газотранспортных системах условно разделяется на два главных направления - аналитическое и численное.
Импульсом к развитию аналитического подхода послужило упрощение исходных уравнений предложенное И.А.Чарным [148]. Упрощенные
8 линеаризованные системы дифференциальных уравнений хорошо изучены. Привлекательность указанного подхода состоит в возможности построения локального решения, используя эффективные методы интегральных преобразований (Б.М.Мардонов, Г.А.Хромов [93]; Л.Х.Мирзаджанзаде, М.А.Гусейнзаде [98]; В.А.Юфин, А.И.Мамедов, В.А.Аллахвердиев [151, 152]). Проблема многомерности, порожденная структурой газотранспортной системы, в настоящее время нашла наиболее эффективное решение благодаря математической модели, разработанной С.А.Бобровским, С.Г.Щербаковым, М.А.Гусейнзаде, В.А.Юфиным и др. (С.А.Бобровский, С.Г.Щербаков, М.А.Гусейнзаде [19]; М.А.Гусейнзаде, В.А.Юфин [47]; М.А.Гусейнзаде, М.Ф.Степанова [46]; М.А.Гусейнзаде, Л.И.Другина, О.Н.Петрова, М.Ф.Степанова [45]). Предпринимаются попытки улучшить точность решения нестационарных задач магистрального транспорта газа путем более сложной линеаризации уравнений движения (В.В.Грачев и др., [43]). Вопросы линеаризации уравнений нестационарного неизотермического течения реального газа в трубопроводах рассматривались Б.Л.Кривошеином, В.П.Радченко [72]. Возможность получения более точных аналитических решений показана, например, в работах А.С.Трофимова, А.В.Судакова, А.В.Козлова [141]; А.С.Трофимова, Е.В.Кочаряна, В.А.Василенко [140].
Другое направление в решении задач магистрального транспорта газа -использование методов вычислительной математики. Достижения последних лет в решении этих задач связаны, как правило, с использованием ЭВМ. Результаты численных расчетов наряду с их самостоятельным значением могут быть использованы для обоснования упрощающих допущений математической модели процесса и оценки точности аналитических методов ее реализации. Как показывает анализ работ по изучению неустановившегося движения газа по трубам, математическое моделирование течения является основным способом исследования нестационарных процессов.
В связи с нестабильностью процесса транспортирования, неравномерностью загрузки и отбора газа необходим учет нестационарности движения. Ввиду значительной протяженности магистралей количество газа аккумулированного в них очень велико. Представляется важным умение определять текущее состояние газотранспортной системы при переменных граничных условиях и изменениях гидравлического режима.
Направление численного моделирования нестационарных течений в газотранспортной системе развито в работах Е.М.Минского, Ю.И.Максимова [97]. Предложены явные и неявные схемы расчета изотермического течения. Явные схемы приводят к простейшим алгоритмам, но обладают меньшей устойчивостью. В большинстве случаев отдают предпочтение неявным схемам. Несмотря на многократное увеличение объема вычислений на каждом временном слое общее число итераций уменьшается (И.С.Березин, Н.П.Жидков [18]).
Современные численные алгоритмы описаны в монографиях М.Г.Сухарева, Е.Р.Ставровского [125]; М.Г.Сухарева, А.М.Карасевича [124]. С целью увеличения быстродействия расчетов прибегают к линеаризации разностных уравнений. Каждый нелинейный член разбивается на произведение сомножителей, из которых линейный относится к новому временному слою, а все нелинейности к предыдущему слою. Полученная система алгебраических уравнений может иметь в итоге специфическую структуру, решение которой достигается методом прогонки, приводящем к устойчивому рекуррентному процессу.
Все более широкое распространение получают экономичные разностные схемы типа «предиктор-корректор» (М.Г.Сухарев, Е.Р.Ставровский [125]). Данный метод объединяет лучшие качества явных и неявных схем. Процесс счета является безитерационным, но обладает высокой устойчивостью. Процесс интегрирования разделен на два этапа, на первом из которых определяется
10 решение на промежуточном временном слое по упрощенным линеаризованным уравнениям, на втором - возмещается потеря точности.
При решении задач нестационарного неизотермического течения газа по магистралям существенно возрастает сложность конечно-разностной схемы ввиду добавления к системе уравнений энергии. Были предложены явные схемы расчета течения газа с учетом его теплообмена с окружающей средой. Неявные разностные схемы, решения полной системы уравнений движения газа по трубе, разработаны О.Ф.Васильевым, А.Ф.Воеводиным [33].
Одна из основных задач развития численного моделирования нестационарных течений в газотранспортной системе - максимальное сокращение времени счета. Использование современных численных методов при приемлемых значениях шагов интегрирования позволяет многократно сократить временные затраты. Возможности моделирования возрастают и в связи с совершенствованием самой вычислительной техники.
Идет поиск новых путей приближенного решения уравнений движения. Суть традиционных методов состоит в упрощениях, позволяющих получить аналитическое решение (И.А.Чарный [147]; И.Е.Ходанович [145]). М.Г.Сухарев, Е.Р.Ставровский [125] предлагают принципиально новый подход, основанный на использовании почти стабильных интегралов, связывающих граничные функции. Задачи о движении газов и газовых смесей в трубах рассматривались и многими другими авторами (О.Ф.Васильев, Э.А.Бондарев, А.Ф.Воеводин, М.А.Каниболотский [32]; Ю.В.Лапин, М.Х.Стрелец [88]; Г.В.Гогричиани, А.В.Шипилин [39]).
Быстрое стравливание газа из магистрали, а также нарушение ее герметизации приводит к возникновению интенсивных скоростных потоков, вовлекающих в движение пыль со стенок трубы. Струйное истечение газодисперсной смеси характеризуется сложной структурой. Картина течения заметно усложняется дополнительно возникающими разрывами, которые появляются в процессе дифракции ударной волны при ее выходе в открытое пространство.
Вопросами дифракции при малых углах отклонения потока занимался M.J.Lighthill [171]. Им показано, что в пренебрежении эффектами вязкости процесс дифракции автомоделей. При этом возникающая дифракционная картина существенно зависит от того, является ли течение сверх- или дозвуковым. К числу экспериментальных исследований дифракции на малых углах можно отнести работу I.I.Class [163]. В реальной картине течения даже при малых углах появляются дополнительные поверхности разрывов, не учитываемые в линейном приближении.
При больших углах отклонения потока картина дифракции меняется -происходит отрыв потока, поток поворачивается в веере Прандтля-Майера параллельно линии срыва. Дифракция сильных ударных волн при срывных течениях изучалась в работах B.W.Skews [180]; Г.М.Арутюняна [5]; Т.В.Баженовой, Л.Г.Гвоздевой, Б.С.Комарова, Б.Г.Сухова [13]; Т.В.Баженовой, Л.Г.Гвоздевой, Ю.Б.Жилина [12]. К одним из первых исследований численными методами фронта волны при дифракции, где получены срывные режимы, относятся работы Г.В.Тарнавского, В.И.Хоничева, В.И.Яковлева [126]; В.Н.Ляхова [91]. D.C.Pack [174] теоретически показал, что и при больших углах отклонения потока процесс дифракции должен проходить автомодельно. В работах T.V.Bazhenova, L.G.Gvozdeva, V.S.Komarov, B.G.Suchov [158]; S.Schuls [179] экспериментально обнаружено, что при увеличении интенсивности падающей ударной волны на дифрагирующей ударной волне возникает точка перегиба, которая при дальнейшем увеличении интенсивности превращается в маховскую конфигурацию.
В работе Е.Н.Васильева [31] проведено численное исследование нестационарного истечения сверхзвуковой струи, образующейся при дифракции ударной волны, выходящей из канала круглого сечения в затопленное пространство. Дано сравнение с экспериментом S.A.Sloan,
12 M.A.Nettleton [181]. Особенности ударного запуска сверхзвуковых струй, истекающих из плоских каналов постоянного сечения, анализируются в работе В.Б.Кислякова [67]. Структура стационарных сверхзвуковых струй с примесями взвешенных частиц исследуется в работах M.Sommerfeld [182-184]. Численное моделирование импульсных струй вязкого теплопроводного газа проводится в работе Н.М.Булгаковой [26].
Процесс ударного выброса огнетушащего порошка изучался в институте гидродинамики им. М.А.Лаврентьева. Разработан принципиально новый способ тушения пожаров на фонтанирующих скважинах. Возникающая при выбросе газопорошковая струя сворачивается в вихревое кольцо, которое сбивает пламя факела. Среди многочисленных исследований в этом направлении можно выделить работы Д.Г.Ахметова, О.П.Кисарова [9]; М.А.Лаврентьева, Б.В.Шабата [86]; Д.Р.Ахметова, Б.А.Луговцова, В.Ф.Тарасова [10].
В основе теоретических исследований двухфазной газовой динамики лежит феноменологический подход предложенный Х.А.Рахматулиным [113]. Анализ характеристических свойств уравнений двухжидкостных течений газа с диспергированными частицами и возникающих в таких течениях поверхностей разрывов выполнен А.Н.Крайко, Л.Е.Стерниным [71]. Дальнейшее развитие теории многоскоростного континуума с введением разрывов с поверхностными свойствами дано в работах А.Н.Крайко [69, 70].
Изучению нестационарных течений аэродисперсных смесей посвящены работы многих авторов (Ю.В.Казаков, А.В.Федоров, В.М.Фомин [63]; С.П.Киселев, Г.А.Руев, А.П.Трунев [66] и др.). Проведена оценка параметров потоков в рамках равновесного приближения с помощью построения инвариантов Римана в псевдогазе. Рассмотрена задача о распаде разрыва применительно к проблеме выброса частиц в свободную атмосферу. Изучена структура ударных волн и комбинированных разрывов в газовзвесях. Исследовано обтекание тел сверхзвуковым потоком запыленного газа (А.Н.Осипцов [107]; Л.А.Егорова, А.Н.Осипцов, В.И.Сахаров [50, 51]).
13 Построением математических моделей механики гетерогенных сред, описывающих течения смесей газа и твердых частиц с учетом вязкости газа и собственного давления дисперсной фазы занимались: А.В.Федоров [142]; Н.Н.Яненко, В.М.Фомин, Федоров и др. [154]; Е.В.Варламов, А.В.Федоров [30]; А.В.Федоров, Н.Н.Федорова [143].
Расчет разлета сжатого объема газовзвеси проведен Ю.В.Казаковым, А.В.Федоровым, В.М.Фоминым [62, 64]. В работе Б.Е.Гельфанда, Ю.В.Казакова, С.П.Медведева, А.Н.Поленова, А.В.Федорова, В.М.Фомина [37] рассмотрен процесс разлета сжатой стратифицированной газопылевой системы. Закономерности формирования и распространения ударных волн при разлете сжатого объема газовзвеси изучались Б.Е.Гельфандом, С.П.Медведевым, А.Н.Поленовым, А.М.Бартеневым [38]; Б.Е.Гельфандом, А.В.Губановым, С.П.Медведевым, Е.И.Тимофеевым, С.А.Цыгановым [36]. В исследовании В.М.Бойко, В.П.Киселева, С.П.Киселева [22] анализируется взаимодействие ударной волны с облаком частиц. Процесс разлета слоя порошкообразной среды рассматривался в работах А.Г.Кутушева, Д.А.Рудакова [80, 81]; А.Г.Кутушева [76]; Д.В.Садина [116, 117]. Систематическое изложение термогидродинамики гетерогенных смесей содержится в работах М.А.Лаврентьева, Б.В.Шабата [86]; Л.Е.Стернина, Б.Н.Маслова, А.А.Шрайбера, А.М.Подвысоцкого [123]; Г.М.Арутюняна [6]; Р.И.Нигматулина [101].
Развитие механики многофазных сред неразрывно связано с экспериментальными исследованиями. Ввиду сложности протекающих физических процессов и их скоротечности в практику исследований внедрены новые комплексы измерительной аппаратуры. В числе работ этого направления выделяются следующие: В.М.Бойко, В.Ф.Климкин [23]; В.М.Бойко, А.А.Карнаухов, В.Ф.Косарев, А.Н.Папырин [21]; В.М.Бойко, А.Н.Папырин [24]; V.M.Boiko, A.N.Papyrin [161]; V.M.Boiko et al [159]; V.M.Boiko, A.N.Malov [160].
Состояние проблемы
На данном этапе исследования процессов неустановившихся течений одно- и двухфазных сред в каналах, которыми занимались и продолжают заниматься многие исследователи: - недостаточно изучено влияние примеси частиц в жидкости гидроразрыва на динамику раскрытия трещины и ее конечную форму при изъятии нефти из низкороницаемых пластов; не сформировано общее мнение в подходе к решению нестационарных задач магистрального транспорта газа; не решено ряд задач прикладного характера.
Особое значение придается вопросам экологии, в связи с возможными случаями нарушения естественной структуры почвы, загрязнения окружающей среды при утечке перекачиваемой жидкости и газов, аварийных разрывов трубопроводов, возникновения стихийных очагов горения.
Цель работы: - решение ряда актуальных задач нефтегазового комплекса путем построения математических моделей нестационарных процессов движения жидкостей, газов и смесей; - исследование предложенных моделей численными и аналитическими методами; - получение количественных оценок и создание комплексов программ. Достижение цели обеспечивает решение следующих задач: изучение процесса формирования трещины гидроразрыва; рассмотрение вопросов связанных с утечкой перекачиваемой жидкости при разгерметизации или разрыве трубопровода; рассмотрение аспектов математического моделирования нестационарных течений газа в разветвленной системе труб с учетом внешних воздействий; - установление особенностей распространения ударных волн в шероховатых трубах, истечения газа при разрушении магистрали высокого давления; - исследование процесса ударного выброса аэродисперсной среды в атмосферу.
Научная новизна
1. Предложена математическая модель и рассмотрены особенности процесса течения вязкой жидкости с примесью частиц по раскрывающейся трещине в пористой среде. Изучено влияние примеси частиц в жидкости гидроразрыва на динамику раскрытия трещины и ее конечную форму.
2. Сделаны количественные оценки утечки перекачиваемой жидкости при гильотинном разрыве подводного трубопровода для случаев напорного и безнапорного истечения. Определена динамика растекания жидкости под ледяным покровом водоема и по поверхности грунта с учетом его влажности и наклона поверхности.
3. Сформулированы эффективные граничные условия сопряжения нестационарных газовых потоков в узлах разветвленной системы труб. Создан программный комплекс, позволяющий определять отклик системы на внешние воздействия.
Проведено исследование влияния сопротивления и теплообмена на характер распространения нестационарных ударных волн в трубе заполненной газом, ударного воздействия при разрушении участка магистрали высокого давления.
Проведено численное моделирование процесса ударного выброса аэродисперсной среды из цилиндрического канала в окружающее пространство. Исследовано влияние основных определяющих параметров на эффект образования вихревого дисперсного кольца и дальность разлета частиц.
Проведен теоретический анализ равновесного состояния смеси газа с твердыми частицами за косыми скачками уплотнения с учетом конечности объема занимаемого дисперсной фазой.
Исследовано влияние объемного содержания частиц в смеси на степень досжатия газа при нормальном отражении от стенки плоской ударной волны. Изучена одномерная модель решетчатого газа как система взаимодействующих частиц.
На защиту выносятся: - математическая модель течения вязкой жидкости с примесью частиц по раскрывающейся трещине в пористой среде. Результаты исследования влияния примеси частиц в жидкости гидроразрыва на динамику раскрытия трещины и ее конечную форму; - количественные оценки утечки перекачиваемой жидкости при гильотинном разрыве подводного трубопровода для случаев напорного и безнапорного истечения, динамики растекания жидкости под ледяным покровом водоема и по поверхности грунта с учетом его влажности и наклона поверхности; эффективные граничные условия сопряжения нестационарных газовых потоков в узлах разветвленной системы труб. Их реализация в программном комплексе, позволяющем определять текущее состояние системы при переменном граничном режиме; результаты исследования влияния сопротивления и теплообмена на характер распространения нестационарных ударных волн в трубе заполненной газом, ударного воздействия при разрушении участка магистрали высокого давления; результаты численного моделирования процесса ударного выброса аэродисперсной среды из цилиндрического канала в окружающее пространство,
17 анализа влияния основных определяющих параметров на эффект образования вихревого дисперсного кольца и дальность разлета частиц; - результаты теоретического анализа равновесного состояния смеси газа с твердыми частицами за косыми скачками уплотнения с учетом конечности объема занимаемого дисперсной фазой, влияния объемного содержания частиц в смеси на степень досжатия газа при нормальном отражении от стенки плоской ударной волны, одномерной модели решетчатого газа как системы взаимодействующих частиц.
Теоретическая значимость
Предложена теоретическая модель процесса течения вязкой жидкости с примесью частиц по раскрывающейся трещине в пористой среде. Предложена модель процесса двухслойного ламинарного течения несмешивающихся жидкостей в круглой трубе и плоском канале, обусловленного различием плотностей сред и внешними воздействиями. Разработаны специальные численные схемы расчета нестационарных газовых потоков в узлах разветвленной системы труб. Предложен метод решения бесконечной системы линейных алгебраических уравнений зацепляющегося типа.
Практическая значимость
Результаты исследования могут быть использованы в практике добычи и транспортировки углеводородного топлива. При проведении гидроразрыва нефтяного пласта, оценке последствий аварийной утечки перекачиваемой жидкости при разрывах подводного/наземного участков трубопроводов, определении зоны безопасности в условиях разрушения участка магистрали высокого давления. Представленный программный комплекс позволяет делать прогноз или вести технологический расчет нестационарных газовых потоков в разветвленной системе труб.
Публикации
По результатам исследования опубликовано 30 печатных работ, в том числе 16 - в рекомендованных ВАКом журналах, из которых 9 работ опубликовано без соавторов.
Личный вклад автора
Результаты, выносимые на защиту, получены лично автором. Автору принадлежит постановка задач по исследованию неустановившихся течений одно- и двухфазных сред в каналах, построение моделей процессов, анализ и интерпретация результатов исследования. В совместных публикациях, посвященных ударному выбросу газа с примесью частиц и утечке перекачиваемой жидкости, исследование выполнялось под руководством А.Г.Кутушева, которому выражаю глубокую признательность и благодарность.
Структура и объем работы
Диссертация состоит из введения, пяти глав, заключения и списка литературы. Общий объем работы составляет 280 страниц, включая 83 рисунка.
Математическая постановка задачи
Рассматривается процесс движения вязкой жидкости с примесью частиц по раскрывающейся трещине, с учетом потерь на просачивание в пористую среду. Предложены две модели развития трещины и исследованы их особенности.
Движения вязкой смеси по раскрывающейся трещине имеет место в технологии гидроразрыва пласта. С целью формирования трещины гидроразрыва, в скважину подается вязкая жидкость (жидкость гидроразрыва) с примесью твердых частиц (проппант). Закачиваемая вязкая смесь, надавливая на берега начальной трещины, расширяет ее. Ввиду просачивания жидкости гидроразрыва в пласт у берегов трещины скапливаются частицы и быстро осаждаются на «дно». После прекращения закачки смеси трещина вновь сужается. Окончательная ее форма определяется распределением по длине слоя осевших частиц. Асимптотика развития трещины указана R.P.Nordgren [173]. А.А.Поздняковым [111] представлена модель, единообразно описывающая весь ход процесса разрыва пласта, включая и зарождение щели. В работе О.Е.Ивашнева, Н.Н.Смирнова [59] проведено исследование роста трещины гидроразрыва в пористой среде. Цель данного исследования - определить влияние примеси частиц в жидкости гидроразрыва на динамику раскрытия трещины и ее конечную форму [133, 137].
Допущения. Движение жидкости с примесью частиц в трещине опишем в односкоростном квазиодномерном приближении. Трещина гидроразрыва предполагается вертикальной. Прикоснувшиеся к стенке частицы оседают на «дно», оставаясь в одном поперечном сечении. Осевший слой частиц не изменяет картину фильтрации. Объемным содержанием жидкости в осевшем слое пренебрегаем. Будем также считать, что закачиваемая смесь жидкости гидоразрыва и взвешенных частиц с объемными долями ах и а2 однородна. В силу принятых допущений объемное содержание взвешенных частиц остается неизменным по всей длине трещины. Размер частиц не учитывается.
Истечение вязкой жидкости из трубы в затопленное пространство
Рассматривается задача об утечке перекачиваемой жидкости при гильотинном разрыве подводного трубопровода. Предложена математическая модель процесса двухслойного ламинарного течения несмешивающихся жидкостей в круглой трубе и в плоском канале, обусловленного различием плотностей сред и внешними воздействиями.
Данная постановка известна в литературе как задача об утечке нефти при гильотинном разрыве подводного трубопровода (см., например, C.Kranenburg [167]; C.Kranenburg, E.Vegt [168]; В.И.Зоненко, Б.И.Ким, В.И.Минаев [55]). Полубесконечная горизонтальная труба заполнена одной жидкостью и полностью погружена в другую жидкость (рис. 2.1). Состояние покоя не может сохраняться, т.к. при этом невозможно согласовать давления по разные стороны от границы раздела сред. Вглубь трубы распространяется волна вытеснения, выдавливающая легкую фазу. Необходимо описать распространение волны вытеснения в трубе. В настоящей работе предлагается модель движения сред отличающаяся от (C.Kranenburg, E.Vegt [168]). Определяется форма автомодельного решения уравнений движения. Модель течения усложняется добавлением общего поступательного потока жидкостей. Основная цель исследования состояла в определении количественных оценок и качественных особенностей рассматриваемых процессов [85].
Течение газа по системе труб при внешних воздействиях
Рассмотрена численная схема расчета нестационарных течений газа в разветвленных системах трубопроводов с учетом внешних воздействий. Предложено моделирование движения газа в изотермическом приближении, а также с учетом изменения его температуры. Представлен программный комплекс расчета нестационарных течений в газотранспортной системе.
Интерес к численному моделированию нестационарных течений газа в трубах вызван возрастанием требований к точности расчета потоков в газотранспортной системе (ГТС). Проблема расчета состоит в нестабильности производственного процесса. Это обусловлено возможными изменениями граничного режима, компрессии, неравномерностью загрузки и отбора газа, отключением магистралей и т.д. Следующий простой пример иллюстрирует суть проблемы. Пусть ГТС состоит из одной магистрали, на ее выходе поддерживается постоянное давление, а на входе давление меняется по периодическому закону. Очевидно, что в зависимости от соотношения характерного времени стабилизации течения и периода колебаний реализуется тот или иной режим течения. Если период колебаний много больше указанного характерного времени, то движение газа квазистационарное; если много меньше, то на выходе имеем ровный поток; если же указанные времена соизмеримы, то на выходе нестабильный поток газа. В первом случае допустим гидравлический расчет. Состояние ГТС определяется граничными значениями, соответствующими данному моменту времени. В последних двух случаях расчета течения необходимо привлекать модель нестационарного движения газа.
Современные методы позволяют исследовать сложные математические модели движения газа в системе трубопровода. Так, например, в работе В.И.Зубова, В.Н.Котерова, В.М.Кривцова, А.В.Шипилина [56] представлены результаты моделирования нестационарных процессов, которые могут развиваться в трубопроводе, проложенном по морскому дну, в том числе и в нештатных ситуациях. Описанию алгоритмов расчета нестационарных движений газа посвящено много работ (см., например, F.Westerwelt [185]; М.Г.Сухарев, А.М.Карасевич [124]; Г.В.Гогричиани, А.В.Шипилин [39]; А.Ф.Воеводин, Р.И.Сафин [34]; В.А.Юфин, А.И.Мамедов, В.А.Аллахвердиев [151]). Особое внимание уделяется условиям устойчивости разностной схемы. Эти вопросы имеют важнейшее значение (А.А.Самарский, Ю.П.Попов [120]). Расчет движения газа в окрестности точки ветвления магистралей создает дополнительные сложности. В данном месте нарушается предположение о квазиодномерности течения, что приводит к нежелательному и существенному изменению шаблона разностной схемы [75, 129, 131]. Подобные проблемы сопряжения потоков возникают и в иных ситуациях (Т.Я.Грудницкая, В.И.Зубов, В.Н.Котеров, В.М.Кривцов, А.В.Шипилин, А.В.Щерпов [44]; С.В.Павлов, И.К.Якушев [108]; И.К.Якушев [153]).
В процессе создания программного комплекса, предназначенного для расчета нестационарных течений в ГТС, возникли две основные математические проблемы: а) произвольная геометрия (структура) ГТС, б) условия сопряжения нестационарных потоков. Составление программного комплекса для конкретной газотранспортной системы неоправданно спецификой ее эксплуатации. Проблема реализации в едином программном комплексе произвольной геометрии ГТС решается матричным способом задания схемы расположения всех ее элементов. В точках ветвления магистралей и локального внешнего воздействия (компрессии, сосредоточенного подвода/отвода газа) возникает разрыв гидродинамических параметров. В случае нестационарного движения, из-за специфики численного счета, связка параметров потока может «сковать» заметную часть расчетной области. Предложено заменить стандартные условия сопряжения потоков (баланса массы и энергии) на разностные уравнения для внутренних узлов ГТС. Указанные уравнения могут быть выведены, например, из интегральных соотношений для контрольного объема, содержащего внутренний узел.
Изотермические ударные волны
Рассматриваются особенности распространения ударных волн в шероховатых трубах. Исследуется влияние сопротивления и теплообмена на характер нестационарного течения газа.
Движение газа по трубе с разрывами встречается во многих случаях, в том числе может иметь место и в трубопроводной системе, рассмотренной в предыдущей главе. Отключение магистралей, изменение компрессии, отвода и подвода газа и т.д. - все это приводит к нестационарным режимам работы системы с возможным возникновением ударных волн. Основная цель исследования состояла в изучении общих закономерностей распространения таких волн по трубе [139].
Ударный выброс газа с примесью частиц
Изучение процессов движения дисперсных смесей газа с твердыми взвешенными частицами представляет научный и практический интерес. Такие смеси широкого используются в качестве двухфазных рабочих тел в различных аппаратах и установках современной техники, связанной с производством, переработкой и транспортировкой инертных и горючих пылевых и сыпучих материалов, с детонационным нанесением покрытий, с импульсным струйным направленным метанием огнетушащих порошков в очаги горения. Возрастающие требования к интенсификации и безопасности технологических процессов и к оптимизации используемого оборудования вызывают необходимость детального изучения поведения дисперсных сред в условиях переменных динамических нагрузок. При этом особое значение приобретают вопросы нестационарного волнового истечения газовзвесей во внешнюю среду из систем высокого давления в условиях их внезапной разгерметизации в аварийных ситуациях. Указанные вопросы актуальны и с точки зрения проектирования и разработки образцов новой техники, предназначенной для дальнего импульсного выброса газодисперсной среды из ствольных установок под действием сжатого газа или продуктов сгорания взрывчатых веществ.
Процесс разлета слоя дисперсных частиц в волнах разрежения и сжатия изучается, как правило, экспериментально на ударных трубах и численно на основе методов и уравнений механики многофазных систем. Аналитические исследования динамики газопылевых сред проводятся, как правило, в рамках упрощенной схематизации газовзвеси с использованием модели «эффективного» газа. До настоящего времени отсутствуют какие-либо полные аналитические решения задачи о нестационарном истечении сверхзвуковой струи, формирующейся при дифракции ударной волны, выходящей из канала трубы в окружающее газовое пространство. Некоторые результаты численного решения указанной задачи приведены Е.Н.Васильевым [31], В.Б.Кисляковым [67].
В работе Б.Е.Гельфанда, А.В.Губанова, С.П.Медведева и др. [36] представлены результаты экспериментального исследования динамики разлета облака инертной газовзвеси из камеры высокого давления (КВД) ударной трубы в газовое пространство камеры низкого давления (КНД). Ю.В.Казаковым, А.В.Федоровым, В.М.Фоминым[64] применительно к условиям опытов Б.Е.Гельфанда, А.В.Губанова, С.П.Медведева и др. [36] выполнено математическое моделирование плоского одномерного движения неравновесной дисперсной смеси газа с твердыми монодисперсными частицами мучной пыли. Б.Е.Гельфандом, С.П.Медведевым, А.Н.Поленовым, А.М.Бартеневым [38] на ударной трубе экспериментально изучен разлет горящей газовзвеси. А.Г.Кутушев, О.Н.Пичугин [79] численно исследовали волновое истечение реагирующей смеси газа и горящих частиц унитарного (содержащего окислитель) топлива из КВД в газ КНД испытательной трубы. В работе А.Г.Кутушева, Д.А.Рудакова [81] для условий опытов С.П.Медведева, А.Н.Поленова, Б.Е.Гельфанда, С.А.Цыганова [95] проведено математическое моделирование разлета слоя сжатой насыпной среды из КВД в КНД ударной трубы. В указанных выше работах, по существу, изучен лишь одномерный разгон дисперсной смеси в канале ударной трубы. Двумерное стационарное сверхзвуковое струйное истечение инертной газовзвеси из канала круглой трубы в окружающее газовое пространство теоретически (численно) и экспериментально исследовано M.M.Sommerfeld [182-184]. Нестационарное истечение сжатой смеси газа с частицами из цилиндрического канала конечных размеров в атмосферу анализируется А.С.Ивановым, В.В.Козловым, Д.В.Садиным [58], однако двухфазное течение за срезом канала не рассмотрено.
