Квеско Бронислав Брониславович, к.ф.-м-н., доцент, зав.кафедрой ПОДЗЕМНАЯ ГИДРОМЕХАНИКА. - презентация

1 Квеско Бронислав Брониславович, к.ф.-м-н., доцент, зав.кафедрой ПОДЗЕМНАЯ ГИДРОМЕХАНИКА

2 ОСНОВНЫЕ ПОЛОЖЕНИЯ Подземная нефтегазовая гидромеханика (ПГМ) наука о движении нефти, воды, газа и их смесей (флюидов) через коллектора. КОЛЛЕКТОРА горные породы, которые могут служить хранилищами флюидов и отдавать их при разработке

3 ФИЗИЧЕСКИЕ ОСНОВЫ ПОДЗЕМНОЙ ГИДРОМЕХАНИКИ Модели Абстрактные Физические Требования адекватности моделей реальным процессам: полнота - содержание достаточного числа признаков реального объекта; непротиворечивость - включенные признаки не должны противоречить друг другу; реализуемость - построенная математическая модель должна допускать аналитическое или численное решение, а физическая - реализацию в искусственных условиях; компактность и экономичность - процессы сбора информации, подготовка и реализация модели должны быть максимально просты, обозримы и экономически целесообразны. Теория осреднения Теория подобия

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

5 ПО СТЕПЕНИ СЖИМАЕМОСТИ ПО ЧИСЛУ ФАЗ РЕОЛОГИЧЕСКИЕ МОДЕЛИ ФЛЮИДОВ а) Несжимаемая - =соnst в) Упругая где c - коэффициент объёмного расширения, c = (7-30) Па -1 - для нефти и (2,7-5) Па -1 для пластовой воды. с) Сжимаемая. р=z R T - р пл > 9 Мпа R - газовая постоянная, Т - температура, z - коэффициент сверхсжимаемости. Гомогенные Гетерогенные Составляющие (компо- ненты) размазаны по пространству и взаимодействуют на молекулярном уровне. Изменение физических и химических свойств непрерывно. uxux Составляющие(фазы) – разделены отчетливыми геометрическими границами и взаимодействуют на поверхностях раздела. Изменение физических и химических свойств разрывно.

6 КОЛЛЕКТОРА ПО ОРИЕНТИРОВАННОСТИ ПАРАМЕТРОВ В ПРОСТРАНСТВЕ изотропныеанизотропные Изотропия - независимость изменения физических параметров от направления Анизотропия - различные изменения по отдельным направлениям. Упорядочные структуры - анизотропны по поверхностным параметрам.

7 МОДЕЛИ КОЛЛЕКТОРОВ ГЕОМЕТРИЧЕСКИЕ МЕХАНИЧЕСКИЕ

8 СМЕШАННЫЕ ТРЕЩИННЫЕ ВИДЫ КОЛЛЕКТОРОВ ПОРОВЫЕ (ГРАНУЛЯР- НЫЕ) трещиновато-пористые, трещиновато-каверновые и т.д. При этом первая часть в названии определяет вид пустот по которым происходит фильтрация. ГЕОМЕТРИЧЕСКИЕ

9 Фиктивный грунт Идеальный грунт ИДЕАЛИЗИРОВАННЫЕ МОДЕЛИ ПОРИСТЫХ КОЛЛЕКТОРОВ Слепок поровых каналов сцементированного песчаника

10 ГЕОМЕТРИЧЕСКИЕ ПАРАМЕТРЫ ПОРИСТЫХ КОЛЛЕКТОРОВ ГРАНУЛОМЕТРИЧЕСКИЙ СОСТАВ ПОРИСТОСТЬ ПРОНИЦАЕМОСТЬ УДЕЛЬНАЯ ПОВЕРХНОСТЬ

11 ГРАНУЛОМЕТРИЧЕСКИЙ СОСТАВ Гранулометрическим составом породы называют количественное (массовое) содержание в породе частиц различной крупности Эффективный диаметр – такой диаметр шаров, образующих эквивалентный фиктивный грунт, при котором гидравлическое сопротивление, оказыва- емое фильтрующейся жид- кости в реальном и эквивалентном грунте, одинаково.

12 ПОРИСТОСТЬ m о = V п /V ПОЛНАЯ ДИНАМИЧЕСКАЯОТКРЫТАЯ Просветность m s = F п /F Для газовых и нефтяных коллекторов в большинстве случаев m=15-22%, но может меняться в широких пределах: от нескольких долей процента до 52%.

13 ПРОНИЦАЕМОСТЬ - параметр породы, характе- ризующий её способность пропускать к забоям сква- жины флюиды. Проницаемость измеряется: в системе СИ - м 2 ; технической системе - дарси (д); 1д=1,02мкм 2 =1, м 2. Физический смысл проницаемости k заключается в том, что проницаемость характеризует площадь сечения каналов пористой среды, по которым происходит фильтрация. ВИДЫ ПРОНИЦАЕМОСТИ АБСОЛЮТНАЯ k ФАЗОВАЯ (ЭФФЕКТИВНАЯ) k i ОТНОСИТЕЛЬНАЯ k i

14 а) коэффициент насыщенности - отношение объёма Vf данного флюида, содержащегося в порах, к объёму пор Vп По виду флюида различают нефтенасыщенность, газонасыщенность, водонасыщенность. б) коэффициент связанности- отношение объёма, связанного с породой флюида V fс, к объёму пор ПАРАМЕТРЫ, СВЯЗАННЫЕ С НАЛИЧИЕМ ФЛЮИДОВ

15 ИДЕАЛИЗИРОВАННЫЕ МОДЕЛИ ТРЕЩИНОВАТО-ПОРИСТЫХ КОЛЛЕКТОРОВ Схема одномерной Схема пространственной модели трещинной среды - раскрытие; l - линейный размер блока породы

16 ПАРАМЕТРЫ ТРЕЩИННОЙ СРЕДЫ ТРЕЩИНОВАТОСТЬГУСТОТА РАСКРЫТОСТЬ т отношение объёма трещин V т ко всему объёму V трещинной среды. отношение полной длины l i всех трещин, находя- щихся в данном сечении трещинной породы к удвоенной площади сечения f Ширина трещины m т = т Г т,

17 т0 - ширина трещины при начальном давлении р 0 ; * т = п l / т0 - сжимаемость трещины; п - сжимаемость материалов блоков; l - среднее расстояние между трещинами. Для трещинных сред l/ т >100 и поэтому сжимаемость трещин высока.

18 МЕХАНИЧЕСКИЕ МОДЕЛИ ДЕФОРМАЦИЯ: 1. УПРУГАЯ ( S ); ПЛАСТИЧЕСКАЯ( S ); КРИП (ПОЛЗУЧЕСТЬ) - постепенное нарастание деформации при постоянном напряжении. 4. ХРУПКАЯ Реологические модели Кулона, Гука, Кельвина, Сен-Венана Абсолютно-твердое тело

19 ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ФИЛЬТРАЦИИ при отсутствии источников - стоков

20 Математическое описание гидродинамических процессов Смысл дифференциального уравнения Каждое из дифференциальных уравнений гидродинамики выражает определенный закон сохранения, в котором в качестве зависимой переменной используется некоторая физическая величина и отражен баланс между различными факторами, влияющими на эту переменную. Зависимыми переменными являются удельные свойства, т. е. свойства, отнесенные к единице массы, примерами являются: массовая концентрация, скорость (т. е. количество движения единицы массы), удельная энергия. Члены дифференциального уравнения такого типа выражают воздействия на единицу объема. Чистое истечение на единицу объема скорость изменения соответствующего свойства в единице объема. Дифференциальное уравнение состоит из членов, каждый из которых выражает воздействие на единицу объема, а сумма баланс этих воздействий.

21 1. Уравнение неразрывности Уравнение неразрывности при установившейся фильтрации : сжимаемой жидкости или G=ρUF=const несжимаемой жидкости или Q=UF=const Где в декартовой системе координат

22 2. Уравнение движения где р * =р+z g, u=dG/dt, G - расход массы жидкости в единицу времени через поверхность равного потенциала (массовый дебит); среда изотропна (k=const, μ =const) = 0 – течение медленное = 0 – изменение кинетической энергии мало - массовая сила сопротивления флюида о скелет горной породы Получаем уравнение движения в форме Дарси В декартовой системе координат

23 Первые экспериментальные наблюдения за движением воды в трубах, заполненных песком, провели французские инженеры А. Дарси (1856 г.) и Ж. Дюпюи ( гг.). Этими работами было положено начало теории фильтрации. Именем Дарси назван линейный закон фильтрации, который он установил, создавая первую систему водоснабжения в Европе. ЗАКОНЫ ФИЛЬТРАЦИИ пористой среды Закон Дарси

24 СКОРОСТЬ ФИЛЬТРАЦИИ: Q= w F п = w m F ЗАКОН ДАРСИ (ЛИНЕЙНЫЙ ЗАКОН ФИЛЬТРАЦИИ) Физический смысл скорости фильтрации - среднерасходная скорость фиктивного потока, в котором расход через любое сечение равен реальному расходу, поля давлений фиктивного и реального потоков идентичны, а сила сопротивления фиктивного потока равна реальной. Гидравлический уклон

25 Уравнение притока в форме Дюпюи Уравнение Дарси нефть газ z = (z c +z к ) / 2; μ = (μ c +μ к ) / 2; z с =z(p с ), μ с =μ (p с ), z к =z(p к ), μ к =μ (p к ). р* = р + pgz - приведенное давление

26 Границы применимости закона Дарси Верхняя граница инерционные силы Нижняя граница неньютоновские свойства Закон Дарси справедлив при соблюдении следующих условий: скорость фильтрации и градиент давления малы; изменение скорости фильтрации и градиента давления малы.

27 Число Рейнольдса Re=wa /μ ; w -характерная скорость течения: а - характерный геометри- ческий размер пористой среды; - плотность жидкости Зависимость Павловского Критическое число Рейнольдса Re кр =7,5-9. Зависимость Щелкачёва: Критическое число Рейнольдса Re кр =1-12. Скорость фильтрации u кр, при которой нарушается закон Дарси, называется критической скоростью фильтрации Верхняя граница

28 Нижняя граница модель с предельным градиентом начальные глинистые ячейки деформируемые ячейки

29 Законы фильтрации при Re > Re кр Одночленные законы: степенная зависимость C, n - постоянные, 1 n 2. Двухчленные зависимости ДарсиКраснопольского структурный коэффициент по Минскому (нефть) структурный коэффициент по Ширковскому (газ) (d – эквивалентный диаметр частиц)

30 Решая двухчленное уравнение фильтрации имеем уравнения притока: - для несжимаемой жидкости - для газа

31 ЗАКОНЫ ФИЛЬТРАЦИИ трещинной среды Линейный закон фильтрации Скорость фильтрации: u=m т w. Формула Буссинеска при представлении течения по трещинам, как течения между двумя плоскими параллельными пластинами Линейный закон фильтрации =k т –проницаемости трещиноватых сред Для трещиновато-пористой среды общая проницаемость определяется как сумма межзерновой и трещинной проницаемостей

32 Границы применимости линейного закона фильтрации трещинной среды Значения критических чисел Рейнольдса значительно зависят от шероховатости: для гладких трещин Re кр =500, для шероховатых - 0,4. Если величина относительной шероховатости меньше 0.065, то её ролью в процессе фильтрации можно пренебречь. Число Re для трещинной среды: