К оглавлению

УДК 550.822.3

Определение коэффициента сжимаемости пустотного пространства коллектора по данным исследований на неустановившихся режимах фильтрации

В.А. КОРОТЕНКО, М.Е. СТАСЮК (ЗапСибНИГНИ)

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

Методы определения g по данным исследований на установившихся режимах фильтрации были предложены ранее [2, 3, 5, 6]. Однако до настоящего времени им не уделялось должного внимания. В работе [3] для определения влияния деформаций на фильтрационные свойства коллектора предлагается использовать кривые падения давления (КПД). При этом допускается, что характер изменения его проницаемости при снижении давления описывается экспоненциальной зависимостью и обусловлен процессами деформации пустотного пространства. Методика [7], основанная на обработке КПД интегральным методом Умрихина - Чарного (упрощение Б.А. Богачева), позволяет определить начальную гидропроводность и коэффициент , характеризующий зависимость гидропроводности пласта от параметра . На основании этой методики был разработан способ определения g [4], однако его применение связано с графическими построениями, что вызывает некоторые затруднения. Во-первых, достоверность результатов зависит не только от экспериментальных данных, но и от точности обработки и аппроксимации кривой; во-вторых, определение искомых параметров возможно только в том случае, когда КПД имеют достаточно четкую криволинейную форму.

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

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

Если , получим общеизвестные формулы [5].

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

где  - гидропроводность начальная и текущая.

В работе [3] приводится зависимость гидропроводности от депрессии для КПД:

где - импульс депрессии,- приток жидкости в ствол скважины, - коэффициент зависимости гидропроводности от .

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

Для двух различных значений  и  после логарифмирования можно записать

Из (3)

где  - среднее значение в интервале, откуда

Из формулы (6) следует, что если КПД аппроксимируется прямой линией, то коэффициенты  и  постоянны. Другой способ определения основывается на методике [4, 7]. Из (4) и (5) следует:

Из (3)

Для пользования формулой (7) нужно знать . Для этого необходимо построить КПД в координатах , найти значение  в различных точках, построить зависимость  от  и по соотношению (8) рассчитать :

Значения  можно вычислить, пользуясь аналитическим методом определения  без учета притока. Согласно [1]

где  

Из (6) и (9) имеем

В случае притока жидкости к стволу скважины из (8) получим

Из приближенной формулы интегрального метода Умрихина - Чарного можно записать

После преобразования (11) получим

Пренебрегая членами, содержащими  и учитывая (7):

Если вторым слагаемым в (12) пренебречь, то

Коэффициент сжимаемости пустотного пространства  по известному значению  или  найдем по формуле [6]

где - коэффициент сжимаемости нефти, 2*10-3 МПа-1

В качестве примера определения  рассмотрим данные обработки КПД по скв. 25, 42, 125 Салымского месторождения. В первых двух они обрабатывались с учетом притока жидкости в ствол скважины после ее остановки, а в скв. 125 - без учета притока.

Обработка КПД по скв. 25 графоаналитическим методом с использованием формулы (8) затрудняется из-за большого радиуса кривизны КПД, поэтому определялся по формулам (12) и (14).

Запись КПД продолжалась 8,5 ч, депрессия возросла от нуля до 19,5 МПа. Если не рассматривать начальный участок КПД ввиду наибольшей неточности экспериментальных данных, то коэффициент  меняется в пределах 0,0455-0,0537 МПа-1,  - в интервале 0,0145-0,0172 МПа-1. При определении графоаналитическим методом указанных коэффициентов для скв. 42 по формулам (8), (14) использовали результаты обработки, приведенные в [7]. Значения этих коэффициентов находили также по формулам (8), (12) - (14):

Значения, полученные графоаналитическим методом по (8), (14), отличаются от полученных аналитически по (12), (14) в среднем на 16 %; результаты более грубого приближения по (13), (14) - в среднем на 8 %. При этом следует учитывать, что графоаналитическая обработка связана с точностью графических построений, что также увеличивает вероятность ошибки при определении . При обработке результатов исследования скв. 125 без учета притока в ствол скважины  рассчитывали по формулам (6), (10), (12), (14). Если начальный участок КПД не учитывать, то наиболее удовлетворительные результаты дает применение формулы (12):

Выводы

1. Разработан аналитический метод определения g по данным исследования скважин на неустановившихся режимах фильтрации.

2. В отличие от ранее предложенного [4,7], данный метод учитывает изменение сжимаемости в зависимости от депрессии на пласт и характеризуется более высокой точностью и разрешающей способностью.

СПИСОК ЛИТЕРАТУРЫ

1.      Бузинов С.П., Умрихин И.Д. Гидродинамические методы исследования скважин и пластов. М., Недра, 1973.

2.      Горбунов А.Т. Разработка аномальных нефтяных месторождений. М., Недра, 1981.

3.      Горбунов А.Т., Николаевский В.Н. О нелинейной теории упругого режима фильтрации.- НТС, Сер. Добыча нефти, М., Недра, 1964, с. 73-95.

4.      Горностаев С.Г., Нестеров В.Н., Пешков В.Е. Метод определения упругости продуктивных отложений баженовской свиты по результатам гидродинамических исследований (Тез. докл. 6-й научно-технической конференции молодых ученых и специалистов). Тюмень, ЗапСибНИГНИ, 1983, с. 85-86.

5.      Котяхов Ф.И. Физика нефтяных и газовых коллекторов. М., Недра, 1972.

6.      Медведский Р.И., Стасюк M.Е. Метод обработки результатов исследования на неустановившихся режимах фильтрации.- В кн.: Комплексирование геолого-геофизических исследований в скважинах. Тюмень, 1983, с. 25-30.

7.      Пешков В.Е., Нестеров В.Н., Асаева Е.Е. Геолого-экономическая эффективность применения рациональной технологии испытания скважин Салымского месторождения.- В кн.: Пути повышения эффективности бурения и испытания разведочных скважин. Тюмень, 1975, с. 87-92.