|
Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
Фильтрационные течения несжимаемой жидкости в неоднородных пластах
Лекция № 14
3.7.1 Общие замечания. В реальных условиях пористая среда редко бывает однородной. Неоднородной называется среда, у которой ее фильтрационные характеристики – пористость и проницаемость различны в различных точках. Однако часто даже в неоднородных пластах могут быть применены рассмотренные выше решения фильтрационных потоков, если эта неоднородность хаотичная (случайная). Тогда большие области пласта на макро уровне можно считать в среднем однородными. Но есть макро неоднородные пласты, в которых отдельные участки существенно различаются по фильтрационным характеристикам. В пластах-коллекторах нефти и газа выделяют следующие основные виды неоднородности. 1. Слоистая неоднородность. Фильтрационные характеристики в пределах слоев постоянны, а между собой различаются. При этом на границе пластов рассматривают два случая: слои гидравлически изолированы (границы непроницаемы) и между слоями существуют перетоки жидкости. Это случаи вертикальной неоднородности. 2. Зональная неоднородность, при которой пласт по площади распространения состоит из нескольких зон различной проницаемости. Это случай латеральной неоднородности. 3. Неоднородность, у которой проницаемость описывается непрерывной функцией от координат точек пространства k (x, y, z). Рассмотрим одномерные потоки несжимаемой жидкости, подчиняющиеся закону Дарси в таких неоднородных пластах. а) Слоистая неоднородность.
3.7.2 Прямолинейно-параллельный поток. Горизонтальный пласт постоянной толщиной h и шириной В состоит из n пропластков с толщинами h1, h2…hn, проницаемостью k1, k2…kn и пористостью m1, m2…mn. На контуре давление - Рк; в скважинах - Рr (рис.14.1).
Рис. 14.1. Разрез (а) и план пласта (б)
При отсутствии перетоков жидкости между пропластками распределение давления по координате х не будет зависеть от параметров среды и во всех пропластках будет одинаково. Оно будет аналогичным распределению давления в однородном пласте.
Скорость фильтрации в каждом i-м пропластке будет индивидуальной, т.к. зависит от проницаемости:
Дебит потока Q можно вычислить как сумму дебитов в отдельных пропластках Qi:
Движение частиц жидкости в каждом пропластке будет плоскопараллельным, но уравнения движения разные, из-за неодинаковой скорости фильрации
Для гидродинамических расчетов иногда бывает удобным заменить пласт со слоистой неоднородностью однородным пластом (h, B, Lk) со средневзвешенной проницаемостью, определенной на основе равенства дебитов.
б) Зональная неоднородность. Горизонтальный пласт (h, B, Lk, Рк, Рr) состоит из n зон: (k1, m1, l1), (k2, m2, l2), (ki, mi, li)…(kn, mn, ln); где: ki – проницаемость, mi – пористость, li – длина i-й зоны. Характеристики такого потока в пределах каждой однородной зоны будут рассчитываться по формулам
где: Р(i) (х) – распределение давления в i-й зоне, Рi-1 и Рi – давления на границах зон, играющие роль контурного и забойного давления в скважинах галереи соответственно. Градиент давления в каждой зоне постоянный, но различный в разных зонах
Дебит вследствие неразрывности потока несжимаемой жидкости будет постоянным в любом сечении потока (любой зоне).
Применяя к потоку в i-й зоне свойства пропорций, получим выражение дебита через обобщенные характеристики пласта и граничные значения давлений
Скорость потока также постоянна в любом сечении
При этом надо иметь в виду, что истинные скорости движения частиц будут меньше в зонах с большей пористостью и наоборот. Среднее значение проницаемости k ср такого неоднородного пласта определяется из равенства дебитов зонально - неоднородного и эквивалентного ему однородного пласта с проницаемостью k ср
Отсюда
Давление pi на границе раздела сред с различной проницаемостью, входящие в формулу р(х), определим из условия равенства скоростей фильтрации в этих зонах:
Например, если неоднородный пласт состоит из двух зон, как это часто бывает на практике, то
отсюда Теперь подставляя это решение в выражения распределения давления в зонах Р(1)(х) и Р(2)(х), найдем их выражения:
Если установившееся прямолинейное движение несжимаемой жидкости происходит в пласте, где проницаемость меняется непрерывно и задана функцией к = к (х), тогда дебит такого потока
Разделяя переменные и интегрируя, получим
3.7.3. Плоскорадиальный поток. а) Слоистая неоднородность.
Во всех пропластках распределение давления по цилиндрической координате r будет таким же как и для однородного пласта и подчиняться логарифмическому закону, поскольку граничные давления (Рк, Рс) в них одинаковы (рис. 14.3):
Градиент давления также одинаков:
Скорость фильтрации будет индивидуальной в каждом пропластке и пропорциональна его проницаемости - ki
Дебит будет равен сумме дебитов пропластков:
Средневзвешенное значение проницаемости можно определить из равенства дебитов слоисто-неоднородного и эквивалентного ему однородного пласта с проницаемостью кср:
Откуда б) Зональная неоднородность.
Рис. 14.4
В пласте имеет место установившийся плоскорадиальный поток несжимаемой жидкости по закону Дарси. Распределение давления в каждой i-й зоне будет подчиняться логарифмическому закону, где роль контурных давлений будут играть давления на внешних и внутренних границах зоны:
где: Р(i)(r) – распределение давления в i-й зоне, Рi и Рi-1 – давления на границах зон, играющие роль контурного и забойного давления в скважине, соответственно (рис. 14.5).
Градиент давления в каждой зоне также свой, но подчиняется гиперболическому закону
Дебит потока вследствие неразрывности потока жидкости будет постоянным через любую цилиндрическую поверхность: Скорость фильтрации в любой точке потока найдем как отношение дебита к соответствующему фильтрационному сечению пласта на координате r
Среднее значение проницаемости определим, как обычно, из сопоставления дебита зонально-неоднородного пласта с дебитом эквивалентного ему однородного пласта:
В практике разработки нефтяных месторождений значительный интерес представляет расчет параметров потока неоднородного пласта, состоящего из 2-х зон. Такая задача возникает при кислотной обработке призабойной зоны скважины, ее глинизации или парафинизации, установке гравийного фильтра и т.д. В этих случаях даже в однородном пласте вокруг скважины возникает кольцевая зона с отличными от внешней зоны параметрами (k, m). В этой связи очень важным бывает установить влияние проницаемости кольцевой призабойной зоны и остальной части пласта на ее продуктивность. Схема решения по определению распределения давления в зонах Р(1)(r) и Р(2)(r), а также дебита скважины аналогична рассмотренной ранее задачи для плоскопараллельного потока. Определяют давление Р1 на границе 2-х зон исходя из равенства скорости фильтрации на их границе
откуда
После этого находят Р(1)(r) и Р(2)(r) и дебит скважины. Интерференция скважин.
Лекция № 15 3.8.1 Общие положения
Введем некоторые новые понятия: - точечный сток – точка на плоскости пласта, поглощающая жидкость; сток можно рассматривать как гидравлически совершенную скважину бесконечно малого радиуса в пласте единичной мощности (рис. 15.2 (а)); - точечный источник – точка, выделяющая жидкость (модель нагнетательной скважины) (рис. 15.2 (б)).
Найдем потенциал Ф точечного стока на плоскости. Т.к. точечный сток является моделью добывающей скважины и течение вокруг него плоскорадиальное, то скорость фильтрации можно выразить через удельный дебит (дебит скважины на единицу мощности)
где: w - скорость фильтрации; Q- объемный дебит; q = Q/h - удельный дебит жидкости. Связь потенциала скорости фильтрации с вектором скорости
Если направление касательной
но для плоскорадиального течения Отсюда Потенциал в точке r = 0 и r = ¥ теряет смысл. Эквипотенциальные линии представляют собой семейство окружностей r = const. Для точечного источника выражение потенциала аналогичное, но q< 0. Найдем теперь потенциал точки стока не в плоскости, а в пространстве. Рассуждения аналогичные, что и для стока на плоскости, но движение вблизи такого рода стока радиально-сферическое, поэтому
Для потенциала точечного источника знак дебита меняется на противоположный. Модель точечного стока в пространстве будет в дальнейшем использована для решения различных задач притока жидкости к гидравлически совершенным и несовершенным скважинам. Отметим, что метод стоков и источников находит применение не только для решения задач фильтрации, но также задач теплопроводности, электричества и магнетизма. Вернемся к вопросам интерференции. Математический смысл метода суперпозиции заключается в том, что если имеется несколько источников фильтрационных потоков от скважин с потенциалами Ф1(x, y, z), Ф2 (x, y, z) … Фn (x, y, z), каждый из которых удовлетворяет уравнению Лапласа DФi = 0, то сумма их Гидродинамический смысл метода суперпозиции состоит в том, что давления (потенциалы) в любой точке пласта, вызванные работой каждой скважины (добывающей и нагнетательной) алгебраически суммируется, а вектор суммарной скорости фильтрации частицы жидкости в данной точке находится как геометрическая сумма векторов скоростей, вызванных работой каждой скважины. Пусть на неограниченной плоскости расположены n стоков, потенциал каждого из них в точке М равен Ф(i)M, где: i = 1, 2….n (рис. 15.3).
Каждая из функций потенциалов Ф(i)M удовлетворяет уравнению Лапласа, тогда и суммарный потенциал в точке М
также является его решением. Физически это означает, что фильтрационные потоки накладываются друг на друга. В этом и заключается принцип интерференции. Вектор суммарной скорости фильтрации
где: Метод суперпозиции можно использовать не только в пластах, имеющих круговой контур питания, или бесконечно больших пластах, но и имеющих контур питания иной формы или непроницаемую границу. В этом случае для выполнения граничных условий приходится вводить фиктивные точечные стоки и источники. При этом решение задач в таких пластах сводится к учету одновременной работы и реальных и фиктивных источников. Метод называется методом отображения стоков и источников. Рассмотрим, изложенные здесь принципы суперпозиции, при решении некоторых задачах, имеющих практическое применение в разработке нефтегазовых месторождений.
3.8.2. Приток жидкости к группе скважин в пласте с удаленным контуром питания. Пусть в горизонтальном пласте толщиной h расположена группа скважин А1, А2… Аn, c радиусами rci, работающая с различными забойными потенциалами Фci:
,
где: Так как контур питания достаточно удален от всех скважин, можно приближенно считать, что их расстояния до точек контура одинаковы и равны Rk (рис. 15.4). Потенциал в любой точке пласта, в том числе на забое любой скважины (Фсi), определяется как сумма потенциалов всех источников:
………………………………………………………….
Система состоит из n уравнений и содержит n+1 неизвестных (n дебитов и постоянную интегрирования с). Дополнительное уравнение получим, поместив точку М на контур питания:
Вычитая почленно уравнения системы из последнего уравнения (исключая тем самым с), получим новую систему из n неизвестных относительно qi:
в любой точке пласта М определяется как геометрическая сумма скоростей фильтрации, вызванных работой каждой скважины (рис.15.5).
Рис. 15.5
3.8.3. Приток жидкости к скважине с прямолинейным контуром питания. Пусть в полубесконечном пласте на расстоянии a от прямолинейного контура питания с контурным потенциалом Фk, работает в точке А одна добывающая скважина забойным потенциалом Фс. Требуется найти удельный дебит (q), скорость фильтрации (
Формула потенциала точечного стока Условие постоянства контурного потенциала Фk здесь не выполняется из-за конечного расстояния до контура питания. Для решения задачи используем рассмотренный метод отображения стоков и источников. Влияние прямолинейного контура приводит к появлению фиктивного зеркального источника – q* в точке А / на расстоянии aот прямолинейного контура питания. Потенциал в любой точке М пласта определяется как сумма потенциалов от двух источников:
где: q – действительный сток; -q – фиктивный источник. Потенциал на контуре получим, полагая r1 = r2: Фk = С = const. Постоянство потенциала свидетельствует о корректности применяемого метода. Для вычисления дебита скважины найдем ее забойный потенциал, переместив точку М на забой скважины, т.е. положив r1 = rс и r2 = 2a:
Формула совпала с формулой Дюпюи при условии Rк = 2а. В реальных условиях форма контура питания неизвестна, но вероятней всего она располагается между окружностью радиуса а и прямой, которой соответствует Rк =2а (рис. 15.7).
Поэтому удельный дебит q определяется из неравенства:
Найдем теперь гидродинамическое поле точечного источника возле прямолинейного контура как совокупность эквипотенциалей и линий тока. Уравнение линии равного потенциала можно поучить из выражения потенциала в любой точке М (х, у) пласта Положив этот потенциал постоянной величине и представив радиусы-векторы r1 и r2 в координатой форме, найдем уравнение линии равного потенциала, проходящей через точку М: Это уравнение можно преобразовать к уравнению семейства окружностей с центрами, лежащими на оси x:
Аналогично можно показать, что семейство линий тока также будет представлять окружности, но с центрами на оси у. Окружности будут перпендикулярными к эквипотенциалям и проходить через сток и фиктивный источник (рис.15.8).
3.8.4. Приток жидкости к скважине, расположенной вблизи непроницаемой прямолинейной границы.
где: 2а = r12, отсюда
Лекция №16 3.8.5. Приток жидкости к скважине, эксцентрично расположенной в круговом пласте. Пусть в плоском пласте мощностью h с круговым контуром питания радиуса Rк, на контуре которого поддерживается постоянный потенциал Фk, на расстоянии d от центра в т. A расположена скважина-сток, с забойным потенциалом Фс. Требуется определить дебит скважины и потенциал ФM (х, у) в любой точке пласта M (рис. 16.1).
Воспользуемся методом отображения стока в круге радиусом Rk. В этом случае отображением стока +q в т. A будет источник –q в т. A*, расположенной на продолжении ОА на расстоянии «а» от т. А. Найдем это расстояние из условия постоянства Фk на круге, в частности в 2-х его точках М1 и М2:
Для того, чтобы определить дебит скважины в т. А запишем выражение ее забойного потенциала:
Чтобы избавиться от константы вычтем полученное выражение забойного потенциала из выражения контурного потенциала в т. M1:
Подставляя сюда значение а, получим:
при d = 0 формула переходит в формулу Дюпюи. Выражение потенциала в любой точке М:
Вычитая из этого выражения уравнение ФМ1=Фk и учитывая выражение для «а», получим:
3.9. Метод электрогидравлических аналогий -метод эквивалентных фильтрационных сопротивлений. На примере притока жидкости к нескольким рядам (плоскопараллельный поток) или кольцеобразным батареям скважин (плоскорадиальный поток) ознакомимся с широко применяемым на практике при проектировании разработки месторождений методом эквивалентных фильтрационных сопротивлений. Метод предложен Ю.А. Борисовым и основан на аналогии движения жидкости в пористой среде и электрического тока в проводниках.
Рис. 16.2
Модель (L, Фк, Фс) – цепочка скважин на расстоянии 2s друг от друга и на расстоянии L от прямолинейного контура питания (рис16.2).
Рис. 16.3 При этом удельный дебит каждой скважины по методу отображения равен:
Рис. 16.3
П где:
где: Введение фильтрационных сопротивлений r и r/ позволяет записать удельный дебит в форме аналогичной закону Ома: где: q ® J; (Фk - Фc) ® Uk - Uc. Суммарный дебит всей прямолинейной цепочки из n скважин:
где: 2sn = В – длина цепочки скважин. Аналогично суммарный дебит круговой батареи из n скважин определяется выражением Введем аналогию между гидродинамическими характеристиками фильтрационного потока и характеристиками электрического тока: Q1/® I; (Pk-Pc)®DU; Тогда электрогидравлическая схема для одной цепочки (батареи) скважин будет иметь вид (рис. 16.4):
Рис. 16.4
Пусть теперь в полубесконечном пласте с прямолинейным контуром питания работают три параллельных цепочки добывающих скважин с числом n1, n2, n3 соответственно. Скважины в каждой цепочке имеют одинаковые радиусы Электрогидравлическая схема будет состоять из трех цепочек фильтрационных сопротивлений и будет выглядеть (рис. 16.5):
Рис. 16.5 Расчет схемы производится аналогично расчету разветвленных электрических цепей по законам Ома и Кирхгофа. Составляются алгебраические линейные уравнения по числу неизвестных: дебитов Q1/, Q2/, Q3/ (токов), если известны забойные давления (потенциалы), или наоборот. Внешние сопротивления рассчитываются по формуле
где: Внутренние сопротивления будут
Отметим, что приток жидкости к трем кольцевым батареям скважин с круговым контуром питания рассчитывается по такой же схеме электрических сопротивлений; при этом сохраняются и формулы расчета внутренних фильтрационных сопротивлений, а внешние сопротивления рассчитываются по формуле
При расчете фильтрационных сопротивлений следует учитывать, что номера прямолинейных или кольцевых батарей отсчитываются от контура питания. Контур питания первой батареи совпадает с истинным контуром, а каждой последующей, совпадает с положением линии предыдущей батареи (рис. 16.6). Расстояния Li для расчета внешних фильтрационных сопротивлений плоскорадиального потока показаны на рисунке 16.6-а, а для расчета ана-
Рис. 16.6
логичных сопротивлений круговых батарей соотношение 3.10. Приток жидкости к несовершенным скважинам. Лекция № 17
Скважина называется гидродинамически совершенной, если она вскрывает пласт на всю его мощность и на этом интервале скважина не обсажена (открыта), так что вскрытая поверхность забоя является фильтрующей. Скважина называется несовершенной по следующим условиям: - по степени вскрытия пласта (пласт вскрыт не на всю мощность), при этом ђ = b/h – называется относительным вскрытием пласта, где b – вскрытая мощность; - по характеру вскрытия; скважина обсажена (или в ней находится специальный фильтр) и она сообщается с пластом через перфорационные отверстия в трубе (цементе) или отверстия в фильтре. Встречаются скважины с двойным несовершенством. Гидродинамическое несовершенство скважин оценивается коэффициентом d = Q/Qсов. Приток жидкости к несовершенной скважине даже в бесконечном горизонтально - однородном пласте перестает быть плоскорадиальным. Строгие математические решения задач притока жидкости к несовершенным скважинам представляют большие трудности. Существует несколько подходов к решению таких задач: 1. Расчетный. Задача о притоке жидкости к несовершенным скважинам по степени вскрытия пласта математически исследовалась М. Маскетом. При этом получена следующая формула для дебита:
где:
- интеграл 2-го ряда Эйлера (Гамма-функция, для которой имеются таблицы в математическом справочнике). График j (ђ) имеет вид (рис 17.1).
Рис. 17.1
Иногда для расчета дебитов несовершенных по вскрытию пласта скважин используют формулу И. Козени:
Формула Козени удобна тем, что для рассчета не требует справочных данных Гамма-функции. 2. Электрическое моделирование, основанное на электрогидродинамической аналогии фильтрационных процессов. Ванна заполняется электролитом. В электролит помещается кольцевой электрод, имитирующий контур питания, а в центре его погружается штыревой электрод-скважина на глубину, соответствующую степени вскрытия пласта скважиной. К обоим электродам подается разность потенциалов тока. Возникающий ток является аналогом дебита скважины. Дебит гидравлически несовершенной скважины подсчитывается по формуле:
где: C = с1+с2 – дополнительное фильтрационное сопротивление, вызванное несовершенством скважины по степени вскрытия (с1) и характеру вскрытия (с2). Измеряя разность потенциалов и силу тока, можно подсчитать по закону Ома общее сопротивление, сделать пересчет и определить дополнительные фильтрационные сопротивления. Такие экспериментальные исследования проведены В.И. Щуровым. Им определены: Популярное:
|
Последнее изменение этой страницы: 2017-03-08; Просмотров: 850; Нарушение авторского права страницы