ФНГ / РЭНГМ / ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ГИДРОДИНАМИЧЕСКИХ МЕТОДОВ ИССЛЕДОВАНИЯ ПЛАСТОВ И СКВАЖИН.
(автор - student, добавлено - 6-05-2014, 21:56)
СКАЧАТЬ:
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ГИДРОДИНАМИЧЕСКИХ МЕТОДОВ ИССЛЕДОВАНИЯ ПЛАСТОВ И СКВАЖИН.
5.1. Задачи и методы изучения продуктивных пластов.
Для создания системы разработки того или иного эксплуатационного объекта требуется знать все свойства этого объекта, могущие оказать влияние на технологический процесс извлечения из него нефти или газа. Таких свойств много, изучение их является весьма сложной и трудоемкой задачей. Некоторые свойства пластов можно изучить путем прямого лабораторного исследования отбираемых из них породы, жидкости и газа (керны, пробы). В кернах изучаются пористость, проницаемость, содержание связанной воды и т.д. При исследовании проб определяются: вязкость, плотность, давление насыщения нефти газом, коэффициент растворимости, объемный коэффициент и др. Эти методы называются лабораторными. В нефтегазопромысловой практике применяется широкий комплекс косвенных методов, позволяющих давать количественную оценку интересующих нас свойств пласта на основании изучения других, связанных с ними свойств. Наиболее важными являются промыслово-геофизические и гидродинамические. По данным геофизических исследований (электрокаротаж, радиокаротаж акустические, механические каротажи и т.д.) определяются такие свойства пластов, как мощность, пористость, проницаемость, насыщенность и др. Геофизические методы являются одними из основных при изучении за- лежей. Но, так как зоны пласта в районе скважины обычно имеют поперечные размеры в несколько десятков сантиметров и их объем составляет лишь малую долю объема всего пласта, эти методы так же, как и лабораторные недостаточно надежны при изучении тех свойств пласта, которые резко могут изменяться на площади его распространения (например, проницаемость). Существенное преимущество имеют гидродинамические методы исследования, сущность которых заключается в определении некоторых важных гидродинамических характеристик пластов и скважин по данным измерений дебитов скважин и забойных давлений при установившихся и неустановившихся процессах фильтрации жидкостей и газов в пластах. Простейшим примером их является определение проницаемости образца породы, заключенной в трубе, на основании описанного выше опыта (2.2). Аналогичный смысл имеют и гидродинамические исследования реальных продуктивных пластов. Если, например, определить дебит нефтяной скважины в пластовых условиях и измерить величины пластового и забойного давления, то, воспользовавшись формулой Дюпюи, можно рассчитать величину средней проницаемости пласта в районе скважины: (5.1) Существует ряд формул гидродинамики, дающих связь между дебитами и давлениями при неустановившихся процессах фильтрации. Теоретические основы гидродинамических исследований будут рассмотрены ниже. Изучение залежей начинается сразу же после их открытия. Одной из главных задач исследований в начальный период является получение информации, необходимой для подсчета запасов нефти и газа: 1) выявление границ пласта по площади и разрезу, и распределения его общей и нефтеносной мощности; 2) сведения о пористости коллектора и его нефтегазонасыщенно-сти и об изменчивости этих характеристик; 3) оценка промысловых запасов залежи. Для этого кроме геологических и извлекаемых запасов важно знать товарные качества нефти и газа, а также свойства залежей, определяющих производительность скважин (мощность, проницаемость пласта, вязкость жидкости в пластовых условиях). Главную роль при решении первой задачи играют геофизические исследования, а при решении второй, третьей групп задач необходимо привлекать и лабораторные, и геофизические, и (в основном) гидродинамические исследования. После того, как установлены промышленные запасы и принято решение о вводе залежи в промышленную эксплуатацию, приступают к составлению технологической схемы или проекта разработки залежи. Для этой цели кроме той информации, которая уже имеется и используется в подсчете запасов, необходим комплекс данных об изменении проницаемости (или гидропроводности) пласта по площади залежи и в законтурной области, о продуктивности пласта в целом и отдельных его интервалов в различных частях залежи, об эффективности применяемых способов вскрытия пласта и перфорации скважин, об условиях работы скважин и др. Эти данные получаются на основании исследований. В процессе промышленной эксплуатации исследования пластов и скважин ведутся главным образом гидродинамическими методами. При этом решаются следующие важные задачи: а) уточнение гидродинамических характеристик пластов; б) контроль за ходом выработки пластов по площади и по разрезу; в) выявление действительной технологической эффективности отдельных элементов принятой системы разработки (система поддержания давления, схема расположения скважин; принятый способ вскрытия пластов, способ эксплуатации скважин и др.); г) выявление эффективности мероприятий по повышению или восстановлению производительности скважин.
На базе этих данных в систему разработки вносятся дополнения и усовершенствования (что продолжается в течении всего периода промышленной эксплуатации залежей).
5.2. Основные методы гидродинамических исследований.
Методы гидродинамических исследований скважин и пластов делятся на две основные группы: а) методы, основанные на промысловых измерениях дебитов и давлений при установившихся процессах фильтрации жидкостей и газов в пластах; б) методы, основанные на наблюдениях за изменениями дебитов и давлений во времени при неустановившихся процессах. К первой группе относятся метод установившихся отборов (метод пробных откачек) и метод карт изобар. Ко второй группе относятся метод давления и метод гидропрослушивания.
5.3. Метод установившихся отборов.
Для выявления зависимости между дебитом скважин и величинами забойных давлений при установившихся режимах эксплуатации проводятся специальные циклы исследований. В каждом цикле скважина эксплуатируется на нескольких установившихся режимах, отличающихся величиной дебита и забойного давления. После исследования скважины на одном режиме новый устанавливается не сразу (через несколько часов и суток). На каждом режиме измеряется дебит скважины и забойное давление. Если скважина дает многокомпонентную продукцию (нефть-вода, нефть-вода-газ, вода-газ), то измеряются дебиты отдельных компонентов. Нулевой режим: дебит-ноль а рзаб=рпл. По результатам строятся индикаторные диаграммы, (см. рис.32,33).
Эксплуатационные скважины могут иметь прямолинейные диаграммы, криволинейные - с выпуклостью, обращенные к оси дебитов, и диаграммы, одна часть которых прямолинейна, а другая - криволинейна, (см. рис.34).
Нагнетательные скважины могут иметь диаграммы прямолинейные, криволинейные с выпуклостью, обращенной к оси дебитов или к оси давлений, и диаграммы линейно-криволинейные (см.рис. 35).
Коэффициентом продуктивности эксплуатационной скважины называется отношение ее дебита к перепаду (депрессии) между пластовым и забойным давлением, соответствующим этому давлению:
(5.2) Коэффициент приемистости нагнетательной скважины: (5.2') Сопоставим эту формулу с формулой Дюпюи для случая притока жидкости:
откуда: (5.4) Введем обозначение kh/m=e. Параметр e называют гидропроводностью (д*м/сп или д*см/сп). Тогда (5.4’) Из этой формулы видно, что коэффициент продуктивности скважины зависит от гидропроводности пласта в ее районе, от величины приведенного радиуса и от расстояния между исследуемой скважиной и соседними ее окружающими. В формуле (5.4) величина k обозначает не физическую проницаемость коллектора, а фазовую проницаемость для жидкости, являющейся функцией насыщенности пласта этой жидкостью. При фильтрации однофазной жидкости (другая фаза находится в связанном состоянии) фазовая проницаемость для нее близка к физической проницаемости пласта. Нефтяные, нефтеводяные и водяные скважины (в пласте отсутствует газ) имеют, как правило, прямолинейные диаграммы. Угловой коэффициент прямой (тангенс угла j между индикаторной линией и осью давлений -
дисперсией) численно равен коэффициенту продуктивности скважины (5.5) Искривление линии (при отсутствии газа) может произойти по двум причинам: а) процесс фильтрации жидкостей в газе не подчиняется линейному закону (Дарси); б) изменяется физическая проницаемость пласта при переходе от одних режимов к другим. Существует нижний и верхний пределы скорости фильтрации, между которыми закон Дарси справедлив. Со скоростями меньше нижнего предела обычно не приходиться встречаться, а скорости выше верхнего предела встречаются в пластах вблизи забоев высокопроизводительных эксплуатационных, нагнетательных и особенно газовых скважин. Для оценки наличия этих скоростей обычно используют сами индикаторные линии (выпуклость в сторону дебитов). При увеличении величины давления могут быть расширения трещин в пласте. Искривление индикаторных линий в этом случае будет следующее: для эксплуатационных – выпуклость к оси дебитов; а для нагнетательных – выпуклость к оси давлений. Исходя из вышеизложенного, на основании индикаторных диаграмм устанавливаются оптимальные нормы отбора жидкости и газа а закачки рабочих агентов по скважинам. Сопоставление диаграмм, полученных по одной скважине в различное время, позволяет судить об изменении насыщенности пласта в ее районе, об изменении проницаемости пласта в призабойной зоне и т.п. Мероприятия, проводимые на скважине, также изменяют форму индикаторных диаграмм. Что касается гидродинамических характеристик пластов, то с помощью
индикаторных диаграмм можно непосредственно определить лишь величину средней гидропроводности пластов в районе скважины kh/m=e на момент исследования. Чтобы определить среднюю проницаемость пласта, необходимо по данным других (геофизических и лабораторных) исследований определить параметры: (5.6)
5.3.1. Фильтрация жидкости (нефти, воды, водонефтяной смеси).
Рассмотрим наиболее простые условия, а именно: а) фильтруется однофазная (однокомпонентная) жидкость (вода или нефть); б) индикаторная линия прямолинейная. Исходя из (4.18), приняв в ней p1=pпл, pс=рзаб, R1=dcp, можно найти (5.7) Здесь С - поправка на несовершенство скважины по степени и характеру вскрытия (см. 4.8); rс - действительный радиус скважины. Из уравнения (5.7) можно найти: (5.8) Если коэффициент продуктивности измерен в поверхностных условиях (hпов) и выражен в т/сут*(кгс/см2), то для перехода к размерностям (см3/с)*(кгс/см2) можно воспользоваться соотношением (5.9)
Тогда (5.10) Если коэффициент продуктивности h=tg j измерен в поверхностных условиях в единицах м3/сут*(кгс/см2), то (5.10’) Значения b (объемный коэффициент), gпов (удельный вес в новых условиях) определяется лабораторными исследованиями глубинных и поверхностных проб. dср - определяется из плана размещения скважин, за rс - принимается радиус долота. Для определения С пользуются способом В.И.Шурова. Поправку С можно рассматривать как сумму двух поправок: C=C1+C2, (5.11) где C1 - поправка на несовершенство скважины по степени вскрытия, C2 - то же по характеру вскрытия. Поправку C1 можно определить по графику рис. 36, где занесено несколько линий, каждая из которых характеризует зависимость С, от степени вскрытия d при определенном значении параметра a=h/D, где h - эффективная мощность пласта, D - диаметр скважины (долота). Для определения поправки C2 имеются аналогичные графики рис.37. Каждая кривая на графике отражает зависимость C2 от nD (n - среднее число отверстий на 1 м вскрытого интервала, D - диаметр скважины (в м)) при определенном значении параметра d=d'/D, где d' - средняя величина отношения глубины перфорационных каналов в породе пласта l' к диаметру скважины D равна 0,1 (l=l'/D=0,1). Для других значений d' будут другие диаграммы, где d' - диаметр перфорационных отверстий.
Гидропроводность пласта в районе скважин, по которым индикаторная
диаграмма криволинейна вследствии нарушения линейного закона фильтрации, можно определить по тем же формулам (5.10) или (5.10'). При этом, если начальный участок диаграммы прямолинейный, в формулы должна подставляться величина tg (p этого участка, если вся диаграмма криволинейна, то приближенное значение гидропроводности получают при подстановке в указанные формулы величины тангенса угла между осью давлений и касательной, проведенной к индикаторной линии в точке ее пересечения с осью давлений (см. рис.38). Если кривизна вследствии изменения степени раскрытости трещин в призабойной зоне в зависимости от давления, то гидропроводность пласта будет величиной переменной, зависящей от забойного давления. В этом случае строят несколько диаграмм для различных забойных давлений.
5.3.2.Фильтрация газонефтяной смеси.
В этом случае коэффициент продуктивности в принципе величина переменная и зависит от депрессии. Если обработку индикаторной диаграммы
проводить так же, как при фильтрации однофазной жидкости, то в итоге можно определить значения фазовой гидропроводности и фазовой проницаемости для нефти при разных забойных давлениях. Однако часто необходимо знать значения физической проницаемости коллектора и соответствующей ей величины гидропроводности пласта. Для этой цели применяется специальная методика обработки результатов исследований, основанная на использовании вспомогательной функции Н, имеющей размерность давления и учитывающей изменение фазовой проницаемости жидкости. Эта функция подставляется в функцию Дюпюи вместо давления , (5.12) где wн - дебит нефти в пластовых условиях в см3/с, k - физическая проницаемость коллектора в д, Нпл, Нзаб - значения функции Н при соответствующих давлениях. Если условия исследований не осложнены влияниями изменения размеров трещин в призабойной зоне и нарушениями линейного закона фильтрации, то комплекс 2pkh/mн(1n(dср/rс)+с) в (5.12) в процессе исследования (при переходах с режима на режим ) должен оставаться неизменным. Обозначим его h". Тогда (5.13) Как видно из (5.13), h" есть не что иное, как коэффициент продуктивности, который имела бы скважина в случае притока к ней однофазной жидкости (нефти). В связи с этим индикаторная диаграмма в координатах wн - DН" или q - DН" должна быть прямолинейной, а расчетная функция для определения гидропроводности будет иметь вид, аналогичный (5.10): (5.14) Здесь h" имеет размерность т/сут(кг/см2) (h"=q/DH). Функция Н определяется из выражения: , (5.15) где Н* - безразмерная вспомогательная функция, зависящая от давления, вязкости жидкости и газа, растворимости газа в нефти и от вида кривых относительных проницаемостей, характеризующих процесс фильтрации газожидкостной смеси в изучаемом коллекторе: , (5.16) где mг - вязкость газа в пластовых условиях, G - газовый фактор (G=Qг/Qн, м3/м3). Безразмерную функцию Н* можно найти по эмпирическим зависимостям, вид которых определяется формой диаграмм относительных проницаемостей и величиной параметра a= (mг/ mн*)l, где l - коэффициент растворимости газа в жидкости при пластовом давлении. В частности, для случая, когда коллектор представлен нецементированным песком и a »0.01, справедливы формулы: при 0£р*£15 Н*=0,39р*; при 15£р*£30 Н*=0,623 р*3,306; (5.17) при 30£р*£100 Н*=0,814 р*-10,03. Здесь р* - безразмерное давление. , (5.18) где р - измеренное давление, рат=1атм - атмосферное давление. В результате исследований скважины методом установившихся отборов должна быть получена информация о дебите жидкости, газовом факторе, забойном давлении в нескольких режимах эксплуатации, а также о
пластовом давлении. Кроме того, должны быть известны параметры для определения С, а также данные о вязкости нефти и газа в пластовых условиях и о коэффициенте растворимости l. Имеющаяся информация обрабатывается в следующей последовательности: 1) устанавливается тип коллектора и определяется параметр a. На основании и подбирается тип формулы (5.17) для определения Н*; 2) для каждого режима определяется: а) e ; б) р*заб и р*пл; в)Н*заб и Н*пл; г) Нзаб и Нпл; д) DН=Нпл-Нзаб. Данные заносятся в таблицу; 3) строится индикаторная диаграмма "q-DH"; 4) определяется h"=tgj=q/DH; 5) по (5.14) определяется параметр гидропроводности; 6) по (5.6) определяется проницаемость, если известна эффективная мощность пласта. 5.3.3.Фильтрация газа.
Теоретическими и экспериментальными исследованиями определено, что связь между дебитами и давлениями при установившейся фильтрации газа к скважине наилучшим образом описывается двухчленной формулой вида: , (5.19) где а и в - постоянные коэффициенты, зависящие от параметров пласта, wг - дебит газа в тыс.м3/сут при атмосферном давлении ратм (ратм=1,03 кгс/см2) и стандартной температуре Тст (293 К); рпл и рзаб - в кгс/см2. Произведение в w2г можно рассматривать как поправку на нелинейность закона фильтрации. Если положить b=0, то из (5.19) получим (5.19’) Сравнивая (5.19) с (4.23), предварительно введя в последнюю поправку на сжимаемость (z) и температуру (Тст/Тпл) и выразив в ней h в м, а wг в тыс.м3/сут, можно убедиться, что (5.20) Если внести поправку на несовершенство скважины через C, то получим , (5.20’) где z - средний коэффициент сжимаемости газа в диапазоне давлений рзаб¸рпл при температуре Тпл. Размерности: k - в д, h - в м, m - в сп, dср, rс, rпр - в м (или см), Тпл и Тст - в К. Поскольку коэффициент а пропорционален величине mг/kh=1/e, его определение в процессе исследования газовых скважин имеет такое же значение, как и определение коэффициента продуктивности при исследовании нефтяных скважин. Индикаторная диаграмма газовых скважин имеет вид рис. 39 . Однако для определения а и в более удобны графики зависимостей, изображенные на рис.40.
Из (5.19) следует, что (5.21) Полученный график должен иметь вид прямой линии. Получим линейную диаграмму, где АО=а; b=tgj; BC=a. 5.4. Метод восстановления давления.
5.4.1. Фильтрация жидкости, (нефти, воды, водо-нефтяной смеси).
Метод восстановления давления применяется в условиях проявления в чистом виде упругого режима эксплуатации. Такие условия возникают всегда при остановках, пусках скважин, а также при изменении режима их эксплуатации (см. рис. 3.3). Если по скважине, длительно эксплуатируемой при установившемся режиме, мгновенно изменить дебит на величину Dw, то давление в любой точке пласта (однородного), отстоящей от центра скважины на расстояние R, начнет изменяться в соответствии со следующим законом: (5.22) Здесь Dр - изменение давления в точке наблюдения в кгс/см2, Dw - изменение дебита скважины в пластовых условиях в см3/с, e - гидропроводность пласта в дсм/сп; t - время, отсчитываемое с момента изменения режима скважины в с, R - расстояние до точки наблюдения в см, c - так называемый коэффициент пьезопроводности пласта, имеющий размерность см2/с, Ei - символ специальной функции, аргументом которой является величина -R2/(4ct). Значения Ei определяются по таблицам. При отрицательных значениях аргумента функция Ei положительна, поэтому величины Dр и Dw всегда имеют разный знак. Условимся величину Dw считать положительной в тех случаях, когда дебит скважины (приток из пласта) увеличивается. Приемистость нагнетательных скважин и их увеличение принято считать величинами отрицательными и с увеличением количества жидкости, закачиваемой в скважину, давление в окружающих точках пласта возрастает. Коэффициент пьезопроводности пласта c, характеризует его способность к передаче возмущений (изменений давления), вызываемых изменениями режима эксплуатации скважин. Для однородного пласта , (5.23) где bс , bж - коэффициенты сжимаемости жидкости и пласта в (кгс/см2)-1; b* -коэф. упругоемкости пласта в кгс/см2, m -эффективная пористость в долях единицы. Величина c меняется в пределах от 102 до 105 см2/с. Если наблюдения проводить непосредственно на забое не самой скважины, то вместо R следует подставлять величину приведенного радиуса скважины rпр. При этом аргумент уже через несколько секунд становится очень малым по абсолютной величине. При с
высокой точностью соблюдается приближенное равенство (5.24) С учетом (5.24) формула (5.22) примет вид (5.25) На основе (5.25) строится методика обработки результатов исследований скважин методом восстановления давления. Учитывая, что lnx=2,31gx и lg(x*y)=lgx+lgy, запишем (5.25) в следующем виде: , (5.26) обозначаем: , (5.27) который также не зависит от времени и остается постоянной при восстановлении давления. Обозначим его буквой i. , (5.28) который также не зависит от времени. Переходя от Dw в см3/с в пластовых условиях к Dq в т/сут в поверхностных условиях и от размерности для e - д×см/сп к размерности д×м/сп, найдем: (5.28’) С учетом (5.28) вместо (5.27) можно записать: , (5.27’) Вводя (5.27¢) и (5.28¢) в (5.26), получим: (5.29) Формула (5.29) является уравнением прямой линии в
полулогарифмических координатах (см. рис.41). Приведенные расчетные данные и формулы - для однородного пласта. Но они успешно решают задачи реальных скважин. При этом e, k, c - средние значения этих параметров пласта в районе скважины. Порядок обработки данных исследований скважины методом восстановления давления.
2. Линейный участок графика экстраполируется (продолжается) до пересечения с осью ординат. Определяются А и i. Очевидно, (5.30)
3. Определяется параметр гидропроводности. Из формулы (5.28) следует, что (5.31) 4. Если известны величины h и m, то зная e, определяется проницаемость пласта k, по формуле (5.6). 5. Из формулы (5.27) определяется , (5.32) При этом c определяется из (5.23).
(5.33)
(5.34) Формула (5.34) вытекает из (5.4). Мы рассмотрим случай фильтрации однофазной нефти. Если в пласте имеется связанная вода, то, полученный из расчета, k характеризует фазовую проницаемость пласта для нефти, соответствующую начальной нефтенасы-щенности. Она может отличаться от величины физической проницаемости. А также получим e, характеризующую фазовую гидропроводность пласта при начальном значении нефтенасыщенности. Для водоносных скважин полученные k и e характеризуют соответственно физические проницаемость и гидропроводность. При исследовании нагнетательных скважин k характеризует фазовую проницаемость для воды, a e - фазовую гидропроводность пласта, соответствующие нефтенасыщенности. Эта же методика используется для обработки результатов исследования при фильтрации по пласту водонефтяной смеси, но определяемые при этом величины фазовой проницаемости kж для смеси и гидропроводности eж будут характеризовать некоторые осредненные значения этих параметров, зависящие от величины текущей нефтенасыщенности пласта в районе скважины по мере снижения нефтенасыщенности изменяются и величины kж и eж При наличии диаграммы относительных проницаемостей для данного пласта можно установить четкую зависимость между определяемыми при исследовании kж и eж и средней нефтенасыщенности пласта. Тангенс угла наклона линейного участка кривой восстановления давления (КВД) в полулогарифмических координатах определяется по формуле, аналогичной (5.28). , (5.35) где Dwж - приращение расхода жидкости в пластовых условиях в см/с, eж -осредненное значение фазовой гидропроводности пласта при фильтрации водонефтяной смеси в д×см/сп. Из (5.35) следует, что , (5.36) , (5.37) где wн и wв - изменение суммарного дебита соответственно для нефти и воды. Если известны Dqн пов и Dq в пов, то из (4.14) можно получить: , (5.38) , (5.38’) Таким образом:
Из (5.36) получим: (5.39) Величины Dqн и Dqв здесь имеют размерность т/сут; a gн и gв - т/м3 в поверхностных условиях.
|
|