Подземная гидромеханика (ПГМ) Лектор: к.ф.-м.н., доцент Квеско Б.Б.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Уравнение неразрывности при установившаяся фильтрации : сжимаемой жидкости несжимаемой жидкости - в декартовой системе координат

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

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

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

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

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

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

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

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

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

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

Потенциальным течением будем называть течение, при котором проекции массовой скорости на оси ортогональной системы координат будут являться производными некоторой функции по направлениям данных осей. Уравнения потенциального движения ПОТЕНЦИАЛ ЗАКОН ДАРСИ через потенциал УРАВНЕНИЯ ЛАПЛАСА Установившееся течение ЗАКОН ДАРСИ

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

Замыкающие соотношения [, m, k, μ=f(p)] Зависимость плотности от давления или уравнения состояния а) Несжимаемая - =соnst.в) Упругая где c - коэффициент объёмного расширения,, V c - объём жидкости; c = (7-30) Па -1 - для нефти и (2,7-5) Па -1 для пластовой воды. с) Сжимаемая. р= R T - р пл < 9 Мпа; р < 1 Мпа р=z R T - р пл > 9 Мпа где R - газовая постоянная, Т - температура, z - коэффициент сверхсжимаемости. Изотермический процесс - или

Зависимость пористости от давления эф +р пл =р горн =const р горн = горн g H – горное давление Зависимость вязкости и проницаемости от давления При р < 10 Мпа показатель в выше приведенных экспоненциальных зависимостях меньше 1 и, следовательно, данные зависимости можно разложить в ряд Тейлора. Ограничиваясь первыми двумя членами, получаем где - общее обозначение выше приведённых параметров.

Начальные и граничные условия Начальные условия = о (x,y,z) при t=0 Если при t=0 пласт не возмущён, то = о =const. Граничные условия Внешняя граница : 1)постоянный потенциал (Г, t)= к =const - контур питания ; 2) постоянный расход G=F u=const или 3) переменный поток массы через границу 4) замкнутая внешняя граница 5) бесконечный пласт lim x (Г,t)= к =const y

Внутренняя граница 1) постоянный потенциал (r c, t)= c =const 2) постоянный массовый дебит 3) переменный потенциал на забое (r c,t)=f 2 (t) при r=r c ; 4) переменный массовый дебит 5) не работающая скважина

Установившаяся потенциальная одномерная фильтрация ВИДЫ ОДНОМЕРНЫХ ПОТОКОВ

ИССЛЕДОВАНИЕ ОДНОМЕРНЫХ ТЕЧЕНИЙ Задача исследования установившегося фильтрационного потока заключается в определении дебита (расхода), давления, градиента давления и скорости фильтрации в любой точке потока, а также в установлении закона движения частиц жидкости (или газа) вдоль их траекторий и в определении средневзвешенного по объёму порового пространства пластового давления.

Решение общего дифференциального уравнения Показатель формы потока Начало системы координат: галерея (для прямолинейно- параллельного потока); центр контура скважины в плоскости подошвы пласта (для плоско-радиального потока); центр полусферического забоя скважины (для сферически-радиального потока). Для укрупнённой трубки тока u= G /F( r ), где F=F( r ) - площадь эквипотенциальной поверхности прямолинейно-параллельный поток - F( r )=Bh; плоско-радиальный поток - F( r ) =2 h r; радиально-сферический поток - F( r ) = 2 r 2. G>0 - эксплуатационная скважина

Уравнение Дарси через расход прямолинейно-параллельный поток - A=Bh, j=0; плоско-радиальный поток - A =2 h, j=1; радиально-сферический поток - A = 2, j=2. j - показатель формы потока, т.к. характеризует вид одномерного течения. Уравнение для потенциала (j=0;2) Уравнение для потенциала (j=1)

Выражение для дебита при постоянных потенциалах на границах (j=0;2) (j=1 Уравнение для потенциала (j=0;2) Уравнение для потенциала (j=1)

Вид коллектораХарактери- стики Вид флюидаХарактеристикиПотенциал 1Недеформируемый (пористый) пласт k=const Несжимаемая жидкость =const; μ=const 2Трещиноватый (деформируемый) пласт Несжимаемая жидкость =const; μ =const 3Недеформируемый (пористый) пласт k=const Упругая жидкость μ =const 4Недеформируемый (пористый) пласт k=const Совершенный газ = cт р/ р ст - изотермическое течение; μ =const 5Недеформируемый (пористый) пласт k=const Реальный газ - изотермическое течение ПОТЕНЦИАЛЬНЫЕ ФУНКЦИИ р=z R T – общий случай; μ =const;

Вид коллектораВид флюидаПотенциал 1Недеформируемый (пористый) пласт Несжимаемая жидкость p 2Трещиноватый (деформируемый) пласт Несжимаемая жидкость p 4 3Недеформируемый (пористый) пласт Упругая жидкость 4Недеформируемый (пористый) пласт Газ p 2 5Трещиноватый (деформируемый) пласт Газ p 5

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

Функцииплоско-радиальное Распределение давления Градиент давления Уравнение притока Уравнение Дюпюи Уравнение движения Средневзве- шенное давление

1.Дебит не зависит от r, а только от депрессии d р. График зависимости Q от d р называется индикаторной диаграммой, а сама зависимость - индикаторной. Коэффициент продуктивности скважины 2. Градиент давления и, следовательно, скорость фильтрации обратно пропорциональны расстоянию и образуют гиперболу с резким возрастанием значений при приближении к забою. Анализ:

4. Изобары - концентрические, цилиндрические поверх- ности, ортогональные траекториям. 5. Дебит слабо зависит от величины радиуса контура r к для достаточно больших значений r к /r c, т.к. r к /r c входят в формулу под знаком логарифма. 3. Графиком зависимости р=р( r ) является логарифмическая кривая, вращением которой вокруг оси скважины образуется поверхность, называемая воронкой депрессии.

Течение совершенного газа через недеформируемый пласт по закону Дарси плоско-радиальное течение средневзвешенное давление уравнение движения интегрируем по времени от 0 до t и по расстоянию от R 0 до r, где R 0 - начальное положение частицы флюида Функция Лейбензона

Распределения плоско-радиальное Распределение давления Р=р 2 Градиент давления Уравнение притока Уравнение движения Средневзвешенн ое давление

Анализ – плоско-радиальное течение Распределение давления в недеформируемом пласте 1 - газ; 2 - несжимаемая жидкость Пьезометрическая кривая для газа имеет более пологий характер на большем своём протяжении, чем кривая несжимаемой жидкости; однако у неё более резкое изменение у стенки скважины, чем для несжимаемой жидкости. Распределение давления

Индикаторная зависимость при фильтрации газа по закону Дарси в переменных Q – Δp Уравнение притока Индикаторная зависимость для газа -- параболическая зависимость дебита Q ст от депрессии р к (с осью, параллельной оси дебитов) и линейная зависимость дебита от разницы квадратов пластового и забойного давлений. т.к. р к 2 - р с 2 = 2р к р с - ( р с ) 2 (где р с = р к - р с ) Индикаторная зависимость при фильтрации газа по закону Дарси в переменных Q – Δp 2

Распределение градиента давления Градиент давления вблизи забоя резко возрастает как за счёт уменьшения r, так и за счёт падения давления р, вызванного сжимаемостью газа. Изменение скорости фильтрации Скорость фильтрации слабо меняется вдали от скважины и резко возрастает в призабойной зоне несжимаемая жидкость газ

Течение несжимаемой жидкости в трещиноватом (деформируемом) пласте Потенциальная функция Трещиноватый (деформируемый) пласт Несжимаемая жидкость =const; μ =const Уравнение для потенциала

РАСПРЕДЕЛЕНИЕ ДАВЛЕНИЯ ГРАДИЕНТ ДАВЛЕНИЯ

Объёмный дебит Скорость фильтрации Кривые распределения давления 1- недеформируемый пласт 2 - трещиноватый пласт 1.Воронка депрессии для трещиноватого пласта более крутая, чем для пористого. В деформируемом трещиноватом пласте, за счет уменьшения раскрытости трещин, при снижении пластового давления возникают дополнительные фильтрационные сопротивления, вызывающие резкое понижение давления на сравнительно небольшом расстоянии от скважины, причем более резко снижается давление в пласте с большим *

2. Индикаторная кривая - парабола четвёртого порядка Вид индикаторной кривой при фильтрации несжимаемой жидкости в трещиноватом пласте

Анализ одномерных потоков при нелинейных законах фильтрации где

Течение несжимаемой жидкости в недеформируемом пласте Уравнение фильтрации при u=Q / (2 rh) Распределение давления в пласте Уравнение притока Дебит - положительный корень уравнения притока Индикаторная линия - парабола. Кривая распределения давления - гипербола и воронка депрессии - гипербола вращения. Крутизна воронки депрессии у стенки скважины больше, чем у чисто логарифмической кривой при течении по закону Дарси.

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

Уравнение притока Коэффициенты А и В определяют по данным исследования газовых скважин при установившихся режимах. или

Однородная несжимаемая жидкость в деформируемом (трещиноватом) пласте Закон фильтрации где Уравнение притока через давление и объемный дебит

Индикаторная кривая - результат сложения двух парабол: параболы четвёртого порядка, симметричной относительно оси, параллельной оси дебитов, и параболы второго порядка (относительно дебита Q) симметричной относительно оси, параллельной оси депрессий ( р с ) и отстоящей от последней на расстояние

Идеальный газ в деформируемом (трещиноватом) пласте Закон фильтрации в дифференциальной форме через потенциал Уравнение притока через давление и объемный дебит

Учет скин-эффекта при определении проницаемости в случае фильтрации по закону Дарси Q, cм3/сек ΔР,ат Проницаемость 0,250,19

Зависимость величины проницаемости от закона фильтрации Закон Дарси без скин-эффекта со скин-эффектом 0,25 дарси 0,19 дарси Нелинейный закон Δр=AQ +BQ дарси

ФИЛЬТРАЦИЯ В НЕОДНОРОДНЫХ СРЕДАХ Пласт называется макронеоднородным, если его фильтрационные характеристики (проницаемость, пористость) значительно, скачкообразно отличаются в разных областях. Виды макронеднородности Слоистая Зональная Общая

прямолинейно-параллельный поток плоско-радиальный поток закон Ома I =U / R для последовательных сопротивлений R = R i для параллельных -

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

ЗОНАЛЬНАЯ НЕОДНОРОДНОСТЬ Пласт по площади состоит из нескольких зон различных фильтрационных параметров, на границах которых данные параметры меняются скачкообразно. Массовый дебит постоянен и равен: при плоскорадиальном потоке Квазиоднородное приближение:

Парамет рr1 k1/k2 м 0,10,5210 0,2545, , , ,08472 Q/Q2*100%0,534, , , , , ,051586, , , , , ,2534 Ухудшение проницаемости призабойной зоны сильнее влияет на дебит, чем увеличение проницаемости в этой зоне

Увеличение (в %) дебита скважины при увеличении проницаемости призабойной зоны В случае фильтрации по закону Дарси увеличивать проницаемость призабойной зоны более, чем в 20 раз не имеет смысла, т.к. дальнейшее увеличение проницаемости практически не ведёт к росту дебита (при условии сохранения типа коллектора) Нарушение в пластовых условиях закона Дарси усиливает положительное влияние увеличенной проницаемости призабойной зоны на производительность скважины.

Схема притока к несовершенной скважине: а - по степени вскрытия; b - по характеру вскрытия Приток к несовершенным скважинам Виды несовершенств скважин. Приведённый радиус. Добавочное фильтрационное сопротивление Несовершенная скважина по степени вскрытия - это скважина с открытым забоем, вскрывшая пласт не на всю мощность, а частично (рис. а). Скважина, хотя и доведённая до подошвы пласта, но сообщающаяся с пластом только через отверстия в колонне труб, в цементном кольце или в специальном фильтре, называется несовершенной по характеру вскрытия пласта a) b)

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

Экспериментальные и теоретические исследования притока жидкости к гидродинамически несовершенной скважине Течение по закону Дарси Несовершенство по характеру вскрытия : В.И. Щуров С = С ( a, h) (a=h/D, h - мощность пласта, D - диаметр скважины; h=h вс /h, h вс - толщина вскрытия ). Несовершенство по степени вскрытия: И.М. Доуэлл, Маскет, Р.А. Ховард и М.С. Ватсон С = С (плотности перфорации, глубины прострела ) П лотность перфорации - число отверстий на 1 метр Дебит значительно зависит от плотности перфорации только до значений отверстий на 1 метр

Формула Маскета для дебита несовершенной по степени вскрытия скважины (основа - метод суперпозиции и отображения стоков) f - функция относительного вскрытия Формула Н.К.Гиринского - применяется если толщина пласта много больше радиуса скважины Коэффициент несовершенства

Если скважины несовершенны по характеру вскрытия, то коэффициент С увеличивается на величину сопротивления фильтра D - диаметр фильтрового отверстия в см; n - число отверстий на 1м перфорированной части.

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

1) R 1 (2-3) r c - из-за больших скоростей вблизи перфорации происходит нарушение закона Дарси и проявляется в основном несовершенство по характеру вскрытия; закон фильтрации - двухчленный ; С 3 - по графикам Щурова, а С 4 по формуле N - суммарное число отверстий; R 0 - глубина проникновения перфорационной пули в пласт.

2) R 2 h - линии тока искривляются из-за несовершенства по степени вскрытия;. фильтрация плоскорадиальна, но с переменной толщиной (от h вск до h); закон фильтрации - двухчленный. 3) R 2 < r< R к - действует закон Дарси и течение плоскорадиально

Общее уравнение притока к несовершенной скважине

Влияние радиуса скважины на её производительность Одиночная скважина r с - радиус 1 -ой скважины, r c / =xr c - радиус 2 -ой скважины; G - дебит 1 -ой скважины, G / =уG - дебит 2 -ой скважины;

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

Особенности упругого режима: Неустановившиеся процессы протекают тем быстрее, чем больше коэффициент проницаемости пласта k, и тем медленнее, чем больше вязкость жидкости и коэффициенты объёмной упругости жидкости и пласта. Параметры упругого режима Важнейшие параметры упругого режима: коэффициенты объёмной упругости жидкости и пласта. Коэффициент объёмной упругости жидкости ж характеризует податливость жидкости изменению её объёма и показывает, на какую часть первоначального объёма изменяется объём жидкости при изменении давления на единицу. ж - объём жидкости; знак минус указывает на то, что объём ж увеличивается с уменьшением давления; ж нефти - (7-30) м 2 /н; ж воды - (2,7-5) м 2 /н.

Коэффициент объёмной упругости пласта п - объём пласта; m - пористость; С слабо и сильно сцементированных горных пород находится в пределах (0,3-2) м 2 /н. Упругий запас з - это количество жидкости, высвобождающейся в процессе отбора из некоторой области пласта при снижении пластового давления до заданной величины, если высвобождение происходит за счет объёмного расширения жидкости и уменьшения порового пространства пласта. з = ж 0ж р + с 0 р= * 0 р., где 0ж - объём жидкости, насыщающей элемент объёма пласта 0 при начальном давлении р 0 ; р - изменение давления; * = m ж + с - коэффициент упругоёмкости пласта, показывающий долю объема жидкости от выделенного элемента объема пласта, высвобождающейся из элемента пласта при снижении давления на единицу

Коэффициент пьезопроводности пласта - характеризует скорость распространения изменения пластового давления В коллекторах – 1000см 2 /с æ 50000см 2 /c или 0.1м 2 /с æ 5м 2 /c. Параметр Фурье - определяет степень нестационарности процесса

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

Приток к скважине в пласте неограниченных размеров Вывод основного уравнения упругого режима Пласт - упругий, горизонтальный и большой протяженности и в нём имеется одна скважина, тогда движение жидкости в пласте можно считать плоскорадиальным. Уравнение пьезопроводности в цилиндрических координатах возмущение вызвано мгновенным стоком, существовавшим в момент t = t / Решение

Найдём значения постоянных. Будем считать, что в момент времени t = t / давление в пласте было р = р к = const. Тогда при r> 0 и при t = t / второй член правой части обращается в неопределённость типа / и определяется по правилу Лопиталя, что даёт С = р к. Для определения коэффициента А воспользуемся соотношением для упругого запаса з = * V п р для случая кольцевого элемента пласта с внутренним радиусом r, толщиной h и шириной dr, а также учтем падение давления р = p 0 - p по d з = * рd V п = После интегрирования в пределах от 0 до получим объём жидкости 3, выделившейся из всего пласта и, учитывая выражение для, определим коэффициент А:

Изменение давления во времени для скважины, введенной в неограниченный пласт в некоторый (начальный) момент времени и действующей мгновенно Изменение давления во времени для скважины, действовающей непрерывно с постоянным дебитом Q = Q 0 в течение времени dt / / через сток выделяется из пласта объём d 2 = Qdt /

Интегрально-показательная функция Свойства интегрально-показательной функции: -Ei(-u) изменяется от 0 до при изменении аргумента от 0 до ; функция -Ei(-u) представляется в виде сходящегося ряда Для малых u

Кривая КВД: погрешность не превышает 0,6% для бесконечного пласта при для конечного пласта погрешность расчета давления не превышает 1%, если r к > 1000r c и fo < 3, или Fo < 0,35. Пьезометрические кривые при пуске скважины в бесконечном пласте с постоянным дебитом Выводы: пьезометрические кривые представляют собой логарифмические линии. углы наклона касательных на забое скважины одинаковы для всех кривых. (1) Для точек вблизи забоя

Анализ основной формулы теории упругого режима 1. Основная формула строго справедлива лишь для точечного стока, т.е. при r с =0. Практические расчеты показывают, что ей можно пользоваться даже для укрупнённых скважин (r с 1км) и нельзя использовать только в первые доли секунды после пуска скважины. 2. Вскоре после пуска скважины вокруг неё начинает непрерывно увеличиваться область пласта, в которой для каждого момента времени давление распределяется так, как и при установившемся движении, т.е. давление оказывается квазиустановившимся и пьезометрические кривые будут кривыми логарифмического типа. 3. Стационарная скорость достигается очень быстро на небольших расстояниях от скважины.

Периодически работающая скважина Постановка задачи. В неограниченном пласте останавливается скважина, эксплуатирующаяся с постоянным дебитом Q в течении времени Т, сравнимого со временем проведения исследований. Понижение давления р / в момент времени Т: Повышение давления за счет работы источника с момента времени Т

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

Определение коллекторских свойств пласта по данным исследования скважин нестационарными методами Уравнение КВД гидропроводность пласта коэффициент проницаемости пласта по i = tg и радиусу r c скважины из коэффициента А можно определить коэффициент пьезопроводности пласта æ.

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

Приток к скважине в пласте конечных размеров в условиях упруго-водонапорного и замкнуто- упругого режима Приток к скважине в пласте конечных размеров с открытой внешней границей Постановка Пусть пласт имеет внешнюю границу радиусом r к, через которую может поступать вода при истощении упругого запаса. В центре пласта имеется скважина радиусом r с, которая мгновенно запускается в эксплуатацию с постоянным дебитом Q 0. Перед пуском скважины давление в пласте было р к. Исходные уравнения Уравнение упругого режима Формула Дюпюи Решаем совместно и получаем уравнение для определения давления

р у - установившееся давление в любой точке пласта или в реагирующей бездействующей скважине (при t = или Fo = ). Пьезометрические кривые при пуске скважины в конечном пласте с открытой внешней границей а - с постоянным дебитом; b - с постоянным забойным давлением р с Изменение дебита скважины с течением времени при постоянном забойном давлении р с

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

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

Пьезометрические кривые при пуске скважины в конечном пласте с закрытой внешней границей при постоянном забойном давлении Изменение дебита Q (кр.1) скважины и суммарной добычи Q cp (кр.2) с течением времени t Постоянное забойное давление

Неустановившееся фильтрация газа в пористой среде Уравнение Лейбензона Исходные соотношения Р=р 2, æ -- æ / =,

Пьезометрические кривые при неустановившемся притоке газа к скважине в разные моменты времени (а) и изменение давления с течением времени в фиксированных точках пласта (b) при малых r 2 /(4 æ t) Уравнение (1) используется для расчета коллекторских параметров газовых пластов методом обработки КВД. Принцип расчета такой же, что и в случае нефтяных скважин, но для получения линейной зависимости по оси ординат надо откладывать не депрессию, а разность квадратов пластового и забойного давлений изменение давления (1)

ПРИБЛИЖЕННЫЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧ ТЕОРИИ УПРУГОГО РЕЖИМА Аналитические решения большинства задач теории упругого режима представляются громоздкими формулами в виде бесконечного медленно сходящегося ряда или несобственного интеграла, содержащего специальные функции. В связи с этим были предприняты поиски приближенных эффективных решений задач неустановившейся фильтрации упругой жидкости в упругой пористой среде.

Метод последовательной смены стационарных состояний (ПССС) /развит И.А.Чарным/ Метод основан на предположении, что давление в пласте меняется во времени значительно медленнее, чем по координатам. Поэтому производную по времени можно в первом приближении отбросить, в результате чего для давления получается уравнение Лапласа, описывающее стационарный процесс. В каждый момент времени весь пласт условно разделяется на две области – возмущенную и невозмущенную. В возмущенной области пласта, начинающейся от стенки скважины, давление распределяется по закону установившегося движения жидкости и внешняя граница этой области служит в данный момент контуром питания. В невозмущенной области пласта давление всюду постоянно и равно начальному контурному давлению. Закон движения подвижной границы, разделяющей возмущенную и невозмущенную области, определяется при помощи уравнения материального баланса и граничных условий.

Прямолинейно-параллельный неустановившийся фильтрационный поток упругой жидкости 1.Приток к галерее, на которой поддерживается постоянный дебит Q. Пусть в момент времени t = 0 в горизонтальном пласте постоянной толщины h и ширины В пущена в эксплуатацию прямолинейная галерея, на которой поддерживается постоянный дебит Q. До пуска галереи давление во всем пласте было одинаковым и равным р к. К моменту времени t после пуска галереи граница возмущенной области распространится на длину l(t) (рис 1). Распределение давления в этой области считается установившимся, т.е. описывается линейной зависимостью:. (1) Рис. 1. Кривые распределения давления в прямолинейно-параллельном потоке по методу ПССС

Найдем закон перемещения во времени внешней границы возмущенной области l(t ). Для этого используем тот факт, что количество добытой продукции за время dt равно изменению упругого запаса жидкости в возмущенной зоне пласта за тот же промежуток времени. (2) где объем возмущенной зоны пласта, Согласно закону Дарси, отсюда Т.о. соотношение (2) можно переписать в виде или т.к Q=const Проинтегрируем полученное соотношение

Следовательно, формула для распределения давления в пласте (1) будет иметь вид Погрешность определения давления по сравнению с точной формулой составляет 25%

2. Приток к галерее, на которой поддерживается постоянное забойное давление p r = const В пласте в момент времени t = 0 пущена эксплуатационная галерея с постоянным забойным давлением p r = const. До пуска галереи давление во всем пласте было одинаковым и равным p к. Требуется найти распределение давления, закон перемещения границы возмущенной области l(t ) и изменение дебита галереи во времени Q(t ). Используем уравнение упругого режима (1) V(t), p, Q определим как и в предыдущем разделе:.,. Подставляя указанные величины в (1) подучим.

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

Плоскорадиальный неустановившийся фильтрационный поток упругой жидкости 1.Приток к скважине, на которой поддерживается постоянный дебит Q. Имеем: неограниченный горизонтальный пласт постоянной толщины h. В момент времени t = 0 пущена добывающая скважина радиусом r c с постоянным дебитом Q. До пуска скважины давление во всем пласте было одинаковым и равным p к. В соответствии с методом ПССС принимаем, что через время t после пуска скважины вокруг нее образуется возмущенная область радиусом R(t ), в которой давление будет распределено по стационарному закону. (3) В остальной части пласта сохраняется начальное пластовое давление p к. Требуется найти закон движения границы возмущенной области R(t).

Рис.2. Кривые распределения давления в плоскорадиальном потоке в разные моменты времени по методу ПССС (отбор осуществляется при ус- ловии Q = const) Кривые распределения давления в разные моменты времени приведены на рис. 2. Дебит скважины, очевидно, будет описываться формулой, аналогичной формуле Дюпюи,. Размеры возмущенной области найдем из уравнения материального баланса. (1) при Средневзвешенное пластовое давление в установившемся плоскорадиальном потоке определяется по формуле

Тогда Подставляя приведенные соотношения в уравнение материального балланса (1) получим и после интегрирования в пределах от 0 до t и от r c до R(t ) имеем Давление в любой точке пласта в любой момент времени t Депрессия в момент времени t: Относительная погрешность уменьшается с течением времени и составляет, по вычислениям, 10,6%, если 100; 7,5%, если fo = 103; 5,7%, если fo = 104.

2. Приток к скважине, на которой поддерживается постоянное давление p с = const Движения границы возмущенной области в этом случае можно определить по графику (рис. 3). Дебит скважины определяется по формуле Дюпюи при p с = const. Сравнение с результатами точных расчетов показывает, что погрешность определения дебита по методу ПССС составляет около 5%. Рис Зависимость безразмер- ного радиуса возмущенной облас- ти ( ) c r t R от безразмерного вре- мени fo при отборе жидкости с постоянным забойным давлени- ем c p = const

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

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

Прямолинейно-параллельный неустановившийся фильтрационный поток упругой жидкости 1.Приток к галерее, на которой поддерживается постоянный дебит Q. Пусть в момент времени t = 0 в горизонтальном пласте постоянной толщины h и ширины В пущена в эксплуатацию прямолинейная галерея, на которой поддерживается постоянный дебит Q. До пуска галереи давление во всем пласте было одинаковым и равным р к. К моменту времени t после пуска граница возмущенной области про- двинется на длину l(t ), при этом кривая распределения давления в этой области будет иметь вид параболы (рис.4). Рис. 4. Кривая распределения давления в прямолинейно-параллельном потоке по методу A.M.Пирвердяна

Уравнение параболы, задающей распределение давления в возмущенной области, определяется равенством. (2) Дебит галереи определяется по закону Дарси Продифференцируем выражение для давления и подставим х=0. В результате Т.о. выражение для дебита примет вид (3)

Закон движения границы возмущенной области определяется из уравнения материального баланса. (1) при Значение средневзвешенного пластового давления в возмущенной области к моменту времени t определим теперь, используя распределение (2) Тогда изменение давления и используя выражение (3) для дебита имеем Подставим полученные выражения в уравнение материального балланса (1)

Отсюда и после интегрирования в пределах от 0 до t и от 0 до l Формула для распределения давления (2) в возмущенной области пласта принимает вид Расчет депрессии дает погрешность по сравнению с точным решением примерно 9%, т. е. в 2,5 раза меньше, чем по методу ПССС.

2. Приток к галерее, на которой поддерживается постоянное забойное давление p r = const В пласте в момент времени t = 0 пущена эксплуатационная галерея с постоянным забойным давлением p r = const. До пуска галереи давление во всем пласте было одинаковым и равным p к. Требуется найти распределение давления, закон перемещения границы возмущенной области l(t ) и изменение дебита галереи во времени Q(t )

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

ОСНОВЫ ТЕОРИИ ФИЛЬТРАЦИИ МНОГОФАЗНЫХ СИСТЕМ Углеводородные системы ГомогенныеГетерогенные Составляющие (компоненты) размазаны по пространству и взаимодействуют на моле- кулярном уровне. Изменение физических и химических свойств непрерывно. Составляющие(фазы) - разделены отчетливыми геометрическими границами и взаимодействуют на поверхностях раздела. Изменение физических и химических свойств разрывно.

Характеристики многофазной среды НасыщенностьСкорость фазы Насыщенностью i порового пространства i –й фазы называется доля объема пор V i, занятая этой фазой в элементарном объеме: вектор скорости фильтрации u i фазы определяется как вектор, проекция которого на некоторое направление L равна отношению объемного расхода Q i данной фазы к площадке i, перпендику- лярной к указанному направлению:

Допущение: каждая фаза двигается под действием своего давления Закон фильтрации каждой из фаз: Зависимость относительных проницаемостей k i от насы- щенности Характерная несимметричная форма кривых относительной проницаемости объясняется тем, что при одной и той же насыщенности более смачивающая фаза занимает преимущественно мелкие поры и относительная проницаемость у неё меньше. Сумма относительных проницаемостей для каждого фиксированного значения Присутствие связанной смачивающей фазы мало влияет на течение не смачивающей жидкости, тогда как присутствие остаточной не смачивающей фазы значительно "стесняет" движение смачивающей фазы.

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

Капиллярное давление - р к =р 2 -р 1 Большее давление - на стороне жидкости, не смачивающей твердые зерна породы. Зависимость функции Леверетта от насыщенности : 1 - кривая вытеснения; 2 - кривая пропитки; А - остаточная насыщенность вытесняемой жидкости п - коэффициент межфазного поверхностного натяжения; - статический краевой угол смачивания между жидкостями и породой; m - пористость; J( ) безразмерная функция Леверетта.

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

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

Уравнения движения для многофазной фильтрации + Связь между давлениями Допущение: в любой точке каждая из фаз находится в состоянии термодинамического равновесия процессы однонаправленные Тогда: k i =k i ( ) и р к = р к ( ) + замыкающие отношения

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

Массовая скорость фильтрации капельно-жидкой фазы газированной жидкости Массовая скорость фильтрации свободного газа смеси Массовая скорость фильтрации газа, находящегося в растворе σ м (р) = G гр /G ж - массовая растворимость газа в жидкости, т. е. количество массы газа, раство- ренного в единице массы жидкости при давлении р

Массовая скорость всего газа (суммируем последние два соотношения) Объемный газовый фактор Г представляет собой отношение объемного газового дебита Q г, приведенного к давлению 1 ат, к объемному дебиту жидкого компонента Q ж, приведенному к тем же условиям. Газовый фактор Г в одномерном установившемся потоке сохраняется постоянным вдоль всего потока. г0 и ж0 - значения плотности газа и жидкого компонента - объемная растворимость газа в жидкости σ(р) = р Закон Генри растворимости газа в жидкости ( р

Формула газового фактора через физические параметры т.к. Объемный коэффициент нефти (р) характеризует изменение объема нефти вследствие изменений давления и количества растворенного газа. Величина (р) есть отношение удельных объемов нефти в пластовых и атмосферных условиях

Формула газового фактора через объёмный коэффициент где При постоянном газовом факторе Г данное уравнение, выражая зависимость между давлением р и насыщенностью s, служит уравнением состояния газированной жидкости Функции μ ж (р), μ г (p), (р) и σ(р) определяются по экспериментальным данным. Кривые зависимости коэффициента растворимости газа в нефти и объёмного коэффициента нефти от давления (1)

2. Определяем потенциальную функцию Определение массового дебита жидкой фазы газированной нефти 1. Находим насыщенность s из (1) 4. Подставляя граничные значения (р) в уравнение получаем формулу массового дебита жидкой фазы ε - показатель «несовершенства» жидкости 0

Вывод зависимости для объемного дебита жидкой фазы газированной смеси в плоско-радиальном потоке в виде формулы Дюпюи Допущения: k, ж и μ ж - постоянны 1. Получим выражения для потенциалов на забое и контуре 2. Вычитаем почленно полученные выражения и применяем теорему о среднем 3. Подставим Δ в и разделив на ж находим: k fж – некоторое среднее значение функции k f (р) в интервале изменения р от р с до р к.

При практических расчетах используют зависимости для дебита предложенные: Вогелем, где V c =0,8 Фетковичем Коэффициент С можно интерпретировать как индекс продуктивности пласта, а его значение зависит от подвижности фаз Значение коэффициента С увеличивается с ростом k и h и уменьшается с ростом скина

Резюме: дебит газированной жидкости равен дебиту воображаемой однородной несжимаемой жидкости, движущейся в пласте с коэффициентом проницаемости k' ж, меньшим k. Отличие от классической формулы Дюпюи Для газированной жидкости дебит зависит не только от депрессии р с, но и от величины давления р к или р с. Приближенные выражения для k' ж : 1.И. А. Чарный - k' ж = 0,65 k - для несцементированных песков 2.М. М. Глоговский и М. Д. Розенберг

ВЫВОДЫ по течению газированной жидкости 1.Дебит газированной жидкости при прочих равных условиях всегда меньше дебита однородной несжимаемой жидкости. 2. С повышением газового фактора при неизменяющейся депрессии р с дебит жидкой фазы уменьшается, а дебит газа увеличивается 3.При данной депрессии р с и газовом факторе Г более высокий дебит будет при более высоком пластовом давлении. 4. Для повышения добычи более эффективным средством является увеличение депрессии за счет повышения пластового (контурного) давления р к, но не путем снижения забойного давления р с.

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

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

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

Характерной особенностью данной системы является то, что её можно свести к одному уравнению для насыщенности. Полная система уравнений Схема одномерной двухфазной фильтрации с учетом силы тяжести u=u 1 +u 2 ; = ; функция Баклея Леверетта или функция распределения потоков фаз

Граничные условия: на входе в пласт (нагнетательная скважина или галерея): расход закачиваемой воды равенство нулю скорости фильтрации нефти из последнего условия вытекает, что k 2 = 0, следовательно, на этой поверхности = *. Начальные условия: задаются значения неизвестной функции в зависимости от пространственных координат при t = 0. Можно считать, что при t = 0 насыщенность всюду постоянна (например, = *).

Граничные условия: На выходе из пласта: пренебрегаем градиентом капиллярного давления по сравнению с градиентом давления в фазах, т. е. считаем, что при x = L, откуда следует, что при x = L. Экспериментально установлено, что вода не вытекает из гидрофильного пласта, а накапливается в выходном сечении, пока её насыщенность не достигнет значения *. В момент достижения значения * вода прорывается из пласта с сохранением на выходе этого значения насыщенности. Это явление получило название концевого эффекта. Математически оно приводится к сложному нелинейному граничному условию на выходе

Модель Рапопорта Лиса - для вытеснения без учета силы тяжести. Дифференциальное уравнение для насыщенности в данной модели – параболического типа. Модель Баклея Леверетта - для вытеснения без учета капиллярных сил. Уравнение насыщенности задач данного типа принадлежит к классу квазилинейных гиперболических уравнений первого порядка.

Функция Баклея Леверетта или функция распределения потоков фаз f( ) равна отношению скорости фильтрации вытесняющей фазы к суммарной скорости и равна объемной доле потока вытесняющей жидкости (воды) в суммарном потоке двух фаз определяет полноту вытеснения и характер распределения газоконденсатонефтенасы- щенности по пласту Вид функции Баклея-Леверетта и её производной Задачи повышения нефте- и газоконденсатоотдачи сводятся к применению таких воздействий на пласт, которые в конечном счете изменяют вид функции f( ) в направлении увеличения полноты вытеснения Задача Баклея Леверетта и ее обобщения

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

Задача Рапопорта – Лиса Распределение насыщенности в стабилизированной зоне l Cтабилизированная зона насыщенности перемещается, не изменяя своей формы, и распределение насыщенности в ней при постоянной скорости вытеснения – стационарно.

ОСНОВЫ ФИЛЬТРАЦИИ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ Рассматриваем нелинейные законы фильтрации, описывающие только безинерционные движения при условии, что фильтрующиеся жидкости обладают неньютоновскими свойствами.

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

Стационарно реологические жидкости при > 0, при 0. Вязкопластичные жидкости 0 - начальное (предельное) напряжение сдвига a) n < 1 Псевдопластичные жидкости b) n > 1 Дилатантные жидкости кажущаяся вязкость * убывает с возрастанием градиента скорости. * увеличивается с возрастанием градиента скорости. кр. 1 кр. 3 кр.4 Дилатантная - суспензии с большим содержанием твердой фазы. Псевдопластичная - растворы и расплавы полимеров

ЗАКОНЫ ФИЛЬТРАЦИИ Вязкопластичная жидкость в пористой среде - u>0;, u=0, где предельный (начальный) градиент Индикаторные линии: 1 - линейная аппроксимация неньютоновской жидкости; 2 - реальная неньютоновская жидкость; 3 – ньютоновская по закону Дарси Неньютоновские эффекты проявляются при малых скоростях фильтрации и в средах с малым размером пор, т. е. с малой проницаемостью Из-за неньютоновских свойств нефтей пропластки последовательно включаются в работу по мере превышения градиента давления предельного градиента сдвига.

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

Одномерные задачи фильтрации вязкопластичной жидкости Поток плоскорадиален (u>0); (u=0). Формула притока u=0, если dp/dr Установившееся течение вязкопластичной жидкости Интегрируем формулу притока при р(r c )=р c ; р(R к )=р к

Анализ Часть разности давлений в виде линейного слагаемого с угловым коэффициентом теряется на преодоление предельного градиента сдвига. При Q 0 давление не постоянно (как в случае фильтрации по закону Дарси), а изменяется по линейному закону. При тех же условиях наличие предельного градиента давления в пласте ведет к уменьшению дебита скважины по сравнению с фильтрацией по закону Дарси (формула Дюпюи). Индикаторная линия скважины Q( р с ) - прямолинейная, но не проходит через начало координат, а отсекает на оси депрессий отрезок, равный R к. Распределение давления в пласте Дебит скважины

Слоистый пласт Индикаторные линии при плоскорадиальном течении вязкопластичной жидкости через трёхслойный пласт (у каждого пропластка свои значениями толщин, проницаемости и начального градиента) Неустановившаяся фильтрация вязкопластичной жидкости При решении нестационарных задач на основе модели фильтрации с предельным градиентом в пласте образуется переменная область фильтрации, на границе которой (пока она не достигнет границы пласта) модуль градиента давления должен равняться предельному градиенту а давление - начальному пластовому. Уравнение пьезопроводности:

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

Образование застойных зон при вытеснении нефти водой - эффект фильтрации с предельным градиентом давления Схема образования застойных зон а - между двумя добывающими скважинами; b - при пятиточечной расстановке скважин (1 - нагнетательная скважина; 2 - добывающая скважина; 3 - зона застоя) Отношение незаштрихованных областей ко всей площади пятиточечной ячейки можно считать площадным коэффициентом охвата пласта заводнением. Величина застойной зоны и коэффициент охвата пласта зависят от параметра Коэффициент охвата пласта увеличивается с увеличением параметра

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

Рис Зависимость суммарного дебита от числа скважин Два вида задач: 1.Задаётся дебит скважин и требуется определить необходимое для этого дебита забойное давление и, кроме того, давление в любой точке пласта. 2.Задаётся забойное давление и требуется определить дебит. Последний вид условия встречается наиболее часто в практике разработки НГМ Методы решения Суперпозиции (потенциалов) Комплексного переменного Комформного отображения

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

Метод СУПЕРПОЗИЦИИ При совместном действии в пласте нескольких стоков (эксплуатационных скважин) или источников (нагнетательных скважин) потенциальная функция, определяемая каждым стоком (источником), вычисляется по формуле для единственного стока (источника).

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

Потенциал группы скважин по принципу суперпозиции Потенциал скважины при плоскорадиальном потоке Уравнение эквипотенциальных поверхностей Уравнение эквипотенциальных поверхностей при равенстве дебитов Линии тока образуют семейство кривых, ортогональ- ных изобарам

Приток к совершенной скважине Фильтрационный поток от нагнетательной скважины к эксплуатационной Схема расположения источника 0 1 и стока 0 2 знаки дебитов: источник G 1 = - G, сток G 2 = + G. Исходная формула Для данной постановки

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

Время обводнения Т (х=0; х 0 =2а) Площадь обводнения из равенства объёмов TQ и mh. Расстояние, пройденное частицей за время Т от нагнетательной скважины до эксплуатационной, вдвое больше расстояния пройденного другой частицей за это же время в положительном направлении оси х.

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

Приток к скважине в пласте с прямолинейным контуром питания Граничные условия: = к при r 1 =r 2,т.е. при r 1 /r 2 =1; = с при r 1 =r с, r 2 2а, т.е. при r 1 /r 2 r с /2а; Исходная формула МЕТОД ОТОБРАЖЕНИЯ ИСТОЧНИКОВ (СТОКОВ) - для выполнения тех или иных условий на границах вводятся фиктивные стоки или источники за пределами пласта

Приток к скважине, расположенной вблизи непроницаемой прямолинейной границы Данная задача может возникнуть при расположении добывающей скважины вблизи сброса или около границы выклинивания продуктивного пласта МЕТОД - отображения источника и стока Исходные формулы

Приток к скважине в пласте с произвольным контуром питания Схема видов контуров питания 1.При вычислении дебита скважины форма внешнего контура пласта не имеет сколько-нибудь существенного значения. 2. Чем дальше от внешнего контура пласта находится скважина, тем меньший дебит она имеет. Однако, так как величина расстояния входит под знаком логарифма, то даже значительное изменение этого расстояния мало влияет на величину дебита 3. В случае расположения скважины эксцентрично относительно контура поток можно считать плоско- радиальным и дебит рассчитывать по формуле Дюпюи если r к.>10 3 r c и эксцентриситет а 1 < r к /2.

Приток к бесконечным цепочкам и кольцевым батареям скважин Граничные условия: на контуре питания = к =const при r j =r к ; на контуре скважины = с =const при r 1 =r с ; r j (j 1)=2a sin[(n-1) /n]. Исходные формулы Приток к скважинам кольцевой батареи

При данных гр. условиях: Т.к., то Выражение для дебита одной скважины Область применения: размеры пласта во много раз больше площади внутри окружности батареи скважин (r к 10а ) - случай водонапорного режима. - rк 10а - случай режима растворенного газа

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

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

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

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

Формула Голосова П.П. для общего дебита скважин прямолинейной батареи: - для нечетного числа скважин 2n+1, где n - любое целое число - для четного числа скважин Здесь h - толщина пласта; - расстояние между скважинами; L – расстояние до контура. Ошибка в определении дебитов по данным формулам не превышает 3-4% при L=10км, r с =10см, при расстояниях между скважинами 100м 500м.

Фильтрационное поле бесконечной цепочки равностоящих скважин Формула дебита из формулы дебита скважин круговой батареи при r к = L + a; a = n /(2 ), где L = const - разность между радиусом контура питания и радиусом кольцевой батареи а ; = const - длина дуги окружности радиусом а между двумя соседними скважинами кольцевой батареи. Подставим значения r к, a

Где z= / (2 L), =е Массовый дебит скважин линейной батареи Здесь L - расстояние от контура питания до батареи; - расстояние между скважинами батареи; h - толщина пласта.

Массовый дебит батареи из n скважин Для несжимаемой жидкости

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

Метод эквивалентных фильтрационных сопротивлений (метод Борисова) Метод позволяет сложный фильтрационный поток в пласте при совместной работе нескольких батарей эксплуатационных и нагнетательных скважин разложить на простейшие потоки - к одиночно работающей скважине и к одиночно работающей батареи. закон Ома I =U / R Дебит прямолин. батареи сопротивления внешнеевнутреннее

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

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

«n» нагнетательных и эксплуатационных батарей a)b) Схема n-батарей с двумя контурами питания а) линейные батареи; b) кольцевые батареи

. Электрическая схема n-батарей с двумя контурами питания Сопротивления прямолинейная батарея круговая батарея

Законы Кирхгоффа для последовательных сопротивлений = i, а для параллельных - Приведенные формулы тем точнее, чем больше расстояние между батареями по сравнению с половиной расстояния между скважинами Схема замены соседних батарей скважин одной батареей

Интерференция несовершенных скважин 1) Определяется дебит совершенных скважин с радиусами r с по формулам теории интерференции для притока к стокам и источникам на плоскости. 2) Фильтрационное сопротивление каждой скважины увеличивается на величину коэффициентов несовершенства С i (i = 1,...,4). 3) Используется метод эквивалентных фильтрационных сопротивлений для исследования интерференции несовершенных скважин, в том числе при двухчленном законе фильтрации в виде - нелинейное сопротивление, добавляемое к внутреннему сопротивлению.

Взаимодействие скважин в неоднородно проницаемом и анизотропном пластах Исходные соотношения для дебитов: 1 -ая зона - 2-ая зона - = kФ+С, где Исключим 0 А) Кольцевая батарея во внутренней области

Анализ формулы: 1) При k 1 /k 2 = < 1 величина коэффициента суммарного взаимодействия (отношение суммарного дебита группы совместно действующих скважин к дебиту одиночной скважины) всегда выше, чем U батареи, действующей при тех же условиях в однородном пласте ( = 1). 2) Если же >1, то U будет меньше его значения в однородном пласте. Б) Кольцевая батарея во внешней области (а > R 0 ).

Анизотропный пласт Эффект взаимодействия будет значительно усиленным или ослабленным лишь при резком различии проницаемостей в двух определённых направлениях: в направлении линии расстановки скважин и в направлении перпендикулярном к этой линии. Ослабление взаимодействия наблюдается в случае более низкой проницаемости в направлении линии расстановки скважин по сравнению с проницаемостью в перпендикулярном направлении. Усиление эффекта взаимодействия происходит в обратном случае. Таким образом, для уменьшения эффекта взаимодействия при закладывании новых скважин следует выбирать направление, в котором пласт наименее проницаем.

Взаимодействие скважин В центре батареи действует нагнетательная скважина с дебитом равным дебиту батареи Сравнение дебитов скважин кольцевой батареи из n эксплуатационных скважин в двух случаях: 1)скважины имеют радиус r c и 2)скважины имеют радиус хr c.

Анализ 1) с увеличением числа эксплуатационных скважин кольцевой батареи влияние их радиуса на дебит уменьшается, если отсутствует нагнетание жидкости в пласт; 2) если в центре батареи находится нагнетательная скважина, то влияние радиуса скважины на дебит будет больше, чем при отсутствии центрального нагнетания жидкости в пласт. 3) радиус скважины влияет на производительность больше, чем при одиночной эксплуатационной скважине. Число скважин мало влияет на производительность.

Взаимодействие скважин при неустановившихся процессах По методу суперпозиции n - число скважин; Q j - объемный дебит стока (+) или источника(-) за номером j; р -понижение давления в какой либо точке пласта; r j - расстояние данной точки пласта от скважины за номером j Данная зависимость используется для расчета параметров пласта путем обработки кривой восстановления давления в случае скважины, эксплуатирующейся в течение длительного времени и остановленной для исследования.

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

Функцией комплексного переменного z= х + iy будет комплексное переменное F (z), если указан закон, позволяющий получить значение F (z) no заданному значению z. Рис. 1. Ортогональность изобар и линий тока Каждые две кривые, из которых одна принадлежит семейству кривых, определяемых уравнением (х, у) = С, а другая - семейству кривых (х, у) = С* (С и С* постоянные), пересекаются под прямым углом, т. е. два семейства кривых образуют ортогональную сетку в основной плоскости течения. 1. Общие положения теории функций комплексного переменного F (z) = F (х + iy) = (х, у) + i (х, у) (1) (х, у) и (х, у) - некоторые функции действительных переменных х и у; i мнимая единица. Задать функцию комплексного переменного - значит задать соответствие между парами чисел (х, у) и (, ).

Положения теории функций комплексного переменного Функции (х, у) и (х, у) удовлетворяют уравнению Лапласа, т, е. Положения 1 и 2 справедливы, если выполняются условия: Условия (3) называются уравнениями Коши Римана.

Характеристическая функция, потенциал и функция тока Кривые (х, у) = С - эквипотенциальные линии (для несжимаемой жидкости – изобары) Кривые (х, у) = С* взаимно ортогональны с эквипотенциальными линиями и, следовательно, направление векторов скорости фильтрации будет совпадать в любой данной точке М с направлением касательной к кривой семейства (х, у) = С*, т. е. кривые этого семейства можно считать линиями тока. (При установившемся движении линии тока и траектории частиц жидкости совпадают). Функция (х,у) называется функцией тока. Потенциальная функция течения и функция тока равны действительной и мнимой части некоторой функции F (z) комплексного переменного z Функция F (z) называется характеристической функцией течения (комплексным потенциалом).

Физический смысл функции тока (х, у Функцией тока можно назвать функцию, принимающую на линии тока (х, у) = С* значение (х, у) = С*, равное массе жидкости (газа), протекающей в единицу времени через поперечное сечение канала, построенного на линиях = 0 и (х, у) = С* 1 Функция тока определяется с точностью до произвольной постоянной, зависящей от выбора начальной линии тока = 0. Рис. 2. Распределение потока между двумя параллельными плоскостями 1 и 2 Для несжимаемой жидкости функция тока будет иметь значение объемного (а не массового) расхода жидкости через поперечное сечение канала, построенного на линиях тока = 0 и =С*.

1.Определяем характеристическую функцию, соответствующую данной задаче. 2.Отделяем в характеристической функции действительную часть от мнимой, т. е. определяем потенциальную функцию (х, у) и функцию тока (х, у). 3.Принимая различные значения функции, получаем уравнения семейства эквипотенциальных линий (х, у) = С, а придавая различные значения, находим уравнения семейства линий тока (х, у) = С*. 4. По эквипотенциальным линиям определяем распределение давлений в пласте, по линиям тока - направление движения и характер поля скоростей фильтрации. 5.Находим проекции вектора массовой скорости 6.Определяем массовую скорость Порядок исследования плоских течений с помощью комлексного переменного

Определение массовой скорости Массовая скорость фильтрации равна производной от характеристической функции F (z) по комплексному аргументу z. Вынося во второй скобке множитель i за знак скобки и воспользовавшись уравнениями Коши - Римана получим: Учитывая (4):

Модуль производной от характеристической функции течения равен модулю массовой скорости фильтрации Модуль производной от характеристической функции течения для несжимаемой жидкости будет равен скорости (а не массовой скорости) фильтрации жидкости u.

2. Характеристические функции некоторых основных типов плоского потока Исследуем течения, заданные характеристическими функциями вида F(z) = Az и F(z) = Alnz. I. Плоско-параллельный поток F(z) = Az, где z = x +iy, A - любое комплексное или действительное число (А = А1 + iA2 ) Отсюда Уравнение семейства эквипотенциальных линий - А 1 х – А 2 y = С. (7) Из (7): эквипотенциальные линии - прямые с угловым коэффициентом A1/А2. Уравнение семейства линий тока - А 1 у + А 2 х = С**. (8) Из (8): линии тока -- прямые с угловым коэффициентом (-A2/А1).

Рис. 4. Сетка, изображающая прямолинейно-параллельный поток в направлении, показанном стрелками Заданная характеристическая функция F(z) = Az соответствует прямолинейно- параллельному потоку. Фильтрационное поле представлено ортогональной прямолинейной сеткой, изображенной на рис. 4. Массовая скорость фильтрации Находим производную от F (z) no z Исходные формулы

II. Плоско-радиальный поток А) F(z) = A ln z, где А некоторое действительное число z– в полярных координатах: z = х +i y = r (cos θ + i sin θ) = re iθ, где г радиус - вектор точки; θ полярный угол. F(z) = A In (re iθ ) = A In r + iAθ. =Alnr; =Aθ. Уравнение эквипотенциальных линий - r=const(9) Эквипотенциальные линии - концентрическими окружностями с центром в начале координат Уравнение линии тока - θ = const. (10) Линии тока - прямые, проходящие через начало координат. Рис. 5. Карта эквипотенциальных линий и линий тока при плоско-радиальном (сходящийся или расходящийся) потоке. Центр скважины (сток или источник) находится в начале координат.

Из анализа одномерного плоско- радиального потока имеем Массовая скорость Находим производную от F (z) no z комплексное переменное модуль компл. числа = массовой скорости (11) (12) Приравниваем (11) и (12) находим Тогда (13) М>0 – сток (экспл. скважина) Функция (13) характеризует плоско-радиальное движение жидкости или газа в однородном горизонтальном пласте неограниченной протяженности. Скважина предполагается гидродинамически совершенной.

B) где а = а 1 + ia 2 Течение плоско-радиальное, но особая точка, в которой помещается точечный сток или точечный источник, сдвинута в направлении оси 0х на расстояние а 1, а в направлении оси-0y на расстояние a 2 и, следовательно, центр поперечного сечения скважины находится не в начале координат, а в точке а = а 1 + ia 2. z-а – в полярных координатах ( г расстояние любой точки плоскости потока не до начала координат, а до особой точки а = а 1 + ia 2, в которой помещается сток или источник; θ полярный угол с вершиной в особой точке ). Отсюда:

III. Несколько точечных стоков и источников Потенциальную функцию течения и функцию тока ψ, поддерживаемых всеми стоками и источниками, можно определить по методу суперпозиции, как алгебраическую сумму потенциальных функций течений или функций тока, поддерживаемых отдельными стоками и источниками, если бы каждый из них был единственным в пласте. Характеристическая функция сложного потока F j (z) - характеристическая функция, соответствующая стоку или источнику за номером j, находящемуся в точке а j :

Характеристическая функция течения при совместном действии источника и стока Характеристическая функция течения от нагнетательной скважины к эксплуатационной Рис. 6. Схема расположения источника 0 1 и стока 0 2 Отделяем действительную часть от мнимой

Тогда Уравнение семейства изобар Уравнение линий тока θ 1 - θ 2 =С* (15) (14) а 2 -a 1 =2a

1)Линии тока – окружности с координатами центров 2) Все окружности расположены на прямой, параллельной оси 0У и делящей расстояние между стоком и источником пополам. 3) Радиус окружности 4) Абсциссы точек пересечения Резуме: Линии тока представляют собой окружности, проходящие через центры обеих скважин, и ортогональны окружностям - изобарам. Центры всех этих окружностей расположены на прямой (эквипотенциальной линии), делящей расстояние между скважинами пополам

Характеристическая функция течения для кольцевой батареи скважин Характеристическая функция для n стоков Для плоско-радиального притока (а j комплексное число, определяю- щее положение стока за номером j). а j представим в тригонометрической форме, а –радиус батареи где

Подсчет времени движения частицы несжимаемой жидкости вдоль линии тока Жидкость несжимаема

Время движения частицы жидкости вдоль линии тока (20) (19) z=x-iy - сопряженное с z комплексное переменное. Разделим переменные и интегрируем вдоль линии тока

Стягивание контура нефтеносности (КН) к эксплуатационной кольцевой батарее Время движения частиц КН по линиям тока Характеристическая функция (23)

Время начала обводнения t о (по главной линии) Время движения по главной линии ( =0) ( =const ), (23) (22) Время движения по нейтральной линии ( = /n) Местоположение частицы контура нефтеносности на нейтральной линии тока r / в момент прорыва воды в скважины /из (25)/ /из(24)/

Анализ уравнения (27): 1)Величина г'/а возрастает с увеличением отношения r н /a; 2)Чем больше величина радиуса первоначального контура нефтеносности, тем больше отставание точек контура нефтеносности, движущихся по нейтральной линии тока, от точек контура, движущихся по главной линии тока. 3)При величине радиуса контура нефтеносности г н 2а можно пренебрегать членом уравнения (27) содержащим множитель а/г н и уравнение (27) можно записать в виде 4)Чем больше скважин в батарее, тем меньше отставание частиц контура нефтеносности от тех, которые движутся по главной линии тока т. е. тем равномернее стягивается контур. 5)Формы контура нефтеносности, которая первоначально была в виде окружности, искажается лишь в ближайшей окрестности скважин.

1) Общий объем добытой нефти за время безводной эксплуатации 2) Объем оставшейся в пласте нефти 3) Площадь ω н, занятая оставшейся в пласте нефтью в момент прорыва воды в скважины Вывод: при большом числе скважин в батарее нефтеносная площадь ω н не зависит от величины отношения а/г н.

Метод комформного отображения Функция F (z) комплексного переменного z = х +iy - поле некоторого плоского движения z = z (ς) Полученная из F функция F1 определяет некоторый плоский фильтрационный поток на плоскости ς. Изучив первый поток F, можно изучить поток F 1 (ς). Задаваясь той или иной преобразующей функцией z = z (ς), из одного потока F (z) плоскости z можно получить бесчисленное множество других потоков на плоскости ς. Таким образом, функция z=z(ς) реализует конформное отображение плоскости z на плоскость ς.

Вывод некоторых формул для притока к скважинам при помощи конформного отображения Отделим в функции z=z(ς)действительную часть от мнимой Отсюда x=x(ξ,η), y=y(ξ,η), ξ= ξ(x,y), η= η(x,y). (2) В зависимости от того, однозначна или многозначна преобразующая функция z=z(ς), каждой точке плоскости ς соответствует одна или несколько точек плоскости z. Точно так же каждой линии одной плоскости соответствует одна или несколько линий на другой плоскости. Таким образом, линиям тока и эквипотенциалям, т. е. сетке течения одной плоскости, будет соответствовать вполне определенная сетка течения на другой плоскости. При этом сами значения потенциала скорости Ф и функции тока Ψ будут одинаковыми на соответствующих друг другу линиях обеих плоскостей.

Производная dz/dς - некоторая функция комплексного переменного, определенная в соответствующих друг другу точках обеих плоскостей z и ς. Т.е. предел отношения не зависит от закона стремления к нулю отрезков Δξ и Δη. Отсюда следует, что в каждой точке плоскости ς и соответствующей (или соответствующих) ей точке плоскости z отношение соответствующих бесконечно малых отрезков dz и dς, постоянно Отсюда (4) Аргумент дроби равен разности аргументов числителя и знаменателя. argdz 1 это угол между направлениями элемента dz l и осью х. Таким образом, arg dz l arg dz 2 = arg dς 1 arg dς 2, т. е. углы между отрезками dz 1, dz 2 и отрезками dς 1, dς 2 равны. Преобразование z(ς) или ς (z) называется конформным, так как оно сохраняет подобие бесконечно малых элементов в соответствующих точках. (3)

Приток к скважине Плоскость z - скважина радиусом г с. Плоскость ς - скважина радиусом с При конформном отображении дебиты скважин стоков или источников сохраняются на обеих плоскостях l - произвольный замкнутый контур, окружающий скважину на плоскости z; λ - произвольный замкнутый контур, окружающий скважину на плоскости ς; dn и dl элементы нормали и касательной для контура l на плоскости z; dv и d λ элементы нормали и касательной для контура λ на плоскости ς. Доказательство Абсолютная величина дебита | Q | скважины на плоскости z составляющая скорости фильтрации по нормали к контуру.

По смыслу конформного преобразования, сохраняющего подобие бесконечно малых элементов в соответствующих точках обеих плоскостей, согласно формуле (3) имеем (7)(6) и сокращая на получим В правой части формулы (8) стоит абсолютная величина дебита скважины на плоскости ς, равная абсолютному значению дебита скважины на плоскости z.

Переход от плоско-параллельного к плоско-радиальному потоку методом комформных отображений - характеристическая функция плоско-параллельного потока (А положительная и действительная постоянная). Отделяем действительную и мнимую части Отсюда Эквипотенциали φ = Ах = const - семейство прямых, параллельных оси у Линии тока ψ = Ау = const прямые, параллельные оси х. (рис.) Проекции скорости фильтрации Вывод. Характеристическая функция течения F (z) = Az определяет плоско-параллельное течение в сторону отрицательной оси х с постоянной во всех точках скоростью u =- А.

Замена переменного Сравниваем действительные и мнимые части Прямым линиям х = const плоскости z соответствуют на плоскости ς кривые ln =const, = const, т. е. окружности с центром в начале координат, а прямым у=const лучи θ = const плоскости ς (рис. ). Сетке течения Ф = Ах = const, Ψ = Ay = const на плоскости z соответствует на плоскости ς сетка течения = const и θ = const, т. е. при А >0 приток к точечному стоку в начале координат с дебитом q = 2 А.

Приток к скважине в пласте с прямолинейным контуром питания Исходный поток - приток к точечному стоку на плоскости ς: Плоскость z: в точке х = 0, у = а находится скважина малого радиуса г с ; ось х является эквипотенциалью φ=φ к, а окружность малого радиуса г с эквипотенциалью φ=φ с ; пласт полубесконечен; контур питания прямолинейный Цель задачи: найти преобразование ς = ς (z) или обратное z =z (ς), которое реализует конформное отображение верхней полуплоскости z в круг = к плоскости ς, а точку z c = ia плоскости z, где расположен центр скважины радиусом г с, в начало координат ς = 0 плоскости ς Рис.7. Соответствие между плоскостями z и ς в задаче о притоке в скважину

Искомое преобразование: Доказательство: (11) 1)полагая z=ia, из формулы (11) получаем ς= 0, т. е. центру скважины на плоскости z соответствует начало координат ς= 0 на плоскости ς. 2)полагая в формуле (11) z = х любому вещественному числу, имеем Откуда следует действительная ось z = х перешла в окружность ρ н плоскости ς, а точка верхней полуплоскости z = ia в начало координат ς = 0.

Радиусы скважин обеих плоскостей связаны соотношением Из (11) Для комплексного потенциала на плоскости z получаем По формуле Дюпюи имеем Подставляем ρ с из формулы (12)

ОСНОВЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ МНОГОКОМПАНЕНТНОЙ ФИЛЬТРАЦИИ В основе всех современных методов прогнозирования показателей разработки месторождений природных углеводородов лежат численные методы интегрирования соответствующих дифференциальных уравнений в частных производных, описывающих процессы двухмерной или трехмерной многофазной фильтрации, и расчету картин двухмерных и трехмерных фильтрационных течений на их основе

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

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

Обратные задачи «Активные» обратные задачи - задачи управления, регулирования процесса разработки пласта или месторождения. «Пассивные» обратные задачи - распознавание объектов разработки и уточнение представления о состоянии и свойствах пластовой системы.

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

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

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

Актуальность математического моделирования грамотное использование результатов моделирования реально увеличивает коэффициент нефтеотдачи на %, дополнительная добыча нефти составляет %.

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

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