Обзор инженерных методов расчета дебитов горизонтальных и вертикальных скважин

 

1. Краткий обзор инженерных методов расчета дебитов

 горизонтальных скважин.

 

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

Рассмотрим задачу о квазистационарном  течении жидкости в пористой среде. Одиночная горизонтальная скважина длиной L дренирует область, ограниченную контуром питания с радиусом R_к. Толщина пласта - h, абсолютная проницаемость - K, динамическая вязкость жидкости - m, давление на контуре питания - p_к, давление на забое скважины - p_с, приведенный радиус скважины - r_с. Требуется определить дебит скважины.

Наиболее простое решение было предложено Ю.Т.Борисовым и В.П.Табаковым [1], [2].

Дебит горизонтальной скважины выражается формулой

    (1.1)

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

Внешне фильтрационное сопротивление  по форме совпадает с сопротивлением вертикальной скважины, отличаясь лишь тем, что вместо радиуса скважины r_с используется радиус r_экв = L / 4. Внутреннее сопротивление горизонтальной скважины принимается по методу эквивалентных фильтрационных сопротивлений Ю.П. Борисова [3] равным внутреннему сопротивлению батареи вертикальных скважин шириной L, расстояние между скважинами 2×d = h.

Формула (1.1) имеет тот недостаток, что вне зависимости от длины горизонтальной скважины контур питания предполагается радиальным. Точность данной формулы должна убывать с ростом отношения L/R_к.

 

 Для горизонтальной скважины контур нефтеносности должен иметь  эллипсообразный, а не круговой характер. С учетом этого Giger F. [4] представил формулу притока в горизонтальную скважину в виде

    (1.2)

где R_к - большая полуось эллипса, являющегося контуром питания.

Joshi S. [5] в развитие формулы (1.2) получил выражение

    (1.3)

где           (1.4)

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

Еще одна формула предложена в работе Renard G., Dupuy J. [6].

    (1.5)

где x = 2a / L и a вычисляются по формуле (1.4).

Формула (1.5) по внешнему фильтрационному  сопротивлению совпадает с формулой (1.3), а по внутреннему сопротивлению - с формулами (1.1) и (1.2).

Формулы (1.1) - (1.5) соответствуют случаю изотропного по проницаемости пласта.

С учетом анизотропии по проницаемости Joshi [5] предложена формула

 

    (1.6)

где   - коэффициент анизотропии;

K_г - проницаемость пласта в горизонтально направлении;

K_в - проницаемость по вертикали.

Формула Renard, Dupuy [6] для анизотропного пласта

    (1.7)

где 

Эта формула отличается от (1.6) расчетом внутреннего сопротивления скважины.

В работе Economides M. [7] дана поправка к формуле (1.6)

   (1.8)

Показано [8], что для больших значений коэффициента анизотропии формулы (1.7) и (1.8) являются более точными в сравнении с (1.6).

Для формул (1.2) - (1.8) имеют место некоторые ограничения.

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

L > h,        (1.9)

Для анизотропных пластов, помимо указанных, должно выполняться условие

L > b × h      (1.10)

Формулы (1.1) - (1.8) относятся к расчетам дебита, когда горизонтальная скважина находится в центре пласта относительно кровли и подошвы. Для случая асимметричного расположения скважины относительно кровли и подошвы пласта Joshi [5]

предложена формула

   (1.11)

 

где d - расстояние по вертикали от центра пласта до горизонтальной скважины. При d = 0 формула (1.11) переходит в (1.6).

Опираясь на точное решение П.Я. Полубариновой-Кочиной [9], Меркуловым В.П. [10] получена приближенная формула

  (1.12)

где    ;   ;  

;  ;  .

 

Среди полуэмпирических формул, использующих принцип добавочного фильтрационного сопротивления, отметим формулу Евченко В.С. [11]

     (1.13)

где   

при     ; 

 

Отличительная особенность формул (1.11) - (1.13) в том, что они учитывают асимметричность скважины относительно толщины пласта.

 

2. Сопоставительный  анализ формул притока

 в горизонтальные скважины

 

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

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

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

    (1.14)

Сопоставление расчетов по формулам для изотропного прласта можно сделать по таблице 1.

Поскольку при  взятии отношения дебитов многие параметры взаимно сокращаются, основными параметрами, определяющими величину отношения дебитов q_г/q_в остаются толщина пласта h, длина горизонтальной скважины L, радиус контура питания R_к. В таблице отражены результаты расчетов для h = 5м, h = 10м, h = 50м.  По  радиусу контура питания рассматривались варианты R_к = 100; 200; 500; 1000 метров.

Варианты по длине горизонтальной скважины L отвечают ограничениям (1.9).

Анализ таблицы  показывает, что наиболее близкие  результаты получаются по формулам Борисова-Табакова (1.1) и Renard, Dupuy (1.5). Расхождение результатов по этим формулам имеет место только для случаев L>R_k. Этого следовало ожидать, поскольку погрешность формулы (1.1), как было указано выше, возрастает с ростом отношения L/R_k.

Формула Joshi (1.3) дает несколько заниженные показатели, поскольку в ней внутреннее фильтрационное сопротивление горизонтальной скважины завышено, а формула Giger (1.2), наоборот, дает завышенные значения отношения q_г/q_в, поскольку в этой формуле занижается внешнее фильтрационное сопротивление для горизонтальной скважины.

 

 
 
 Сравнительно неплохо согласуются и результаты расчетов по формуле Меркулова (1.12). Для этой формулы отношение дебитов q_г/q_в постоянно выше, чем по формуле (1.5).

И совсем не согласующиеся  результаты получаются по формуле Евченко (1.13). Они сильно занижаются по сравнению  с другими формулами. Возможно, причиной является ошибка в формуле, допущенная в первоисточнике [11].

 В работе [7] для исходных данных h = 272 фута, L = 774 фута, R_к = 645фута, r_c=0,33 фута приводится значение отношения q_г/q_в, полученное в результате численного моделирования. Используя данное решение в качестве эталонного, была проверена точность инженерных формул (1.1) - (1.5), (1.12), (1.13).

Результаты представлены в таблице 2. По эталонному решению  значение отношения дебитов q_г/q_в при указанных исходных параметрах составляет 2,51. Наименьшее отклонение от эталонного значения достигается при вычислении дебита горизонтальной скважины по формуле Renard, Dupuy (1.5), наибольшее - по формуле Евченко (1.13). Промежуточное положение занимают формулы (1.1), (1.13), (1.2) и (1.3).

Таким образом, на основании данных таблиц 1 и 2 можно  сделать вывод, что наиболее точной инженерной формулой расчета дебита горизонтальной скважины является (1.5). Можно рекомендовать также формулу (1.1), дающую близкие результаты к (1.5).

 

 

3. Исследование  влияния на дебит горизонтальной  скважины параметров пластовой  системы и длины ствола.

 

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

 

 
 Влияние длины горизонтального  ствола однозначно - чем больше длина  ствола, тем больше дебит горизонтальной скважины. Однако зависимость дебита от длины L отнюдь не линейна.

К примеру, при  толщине пласта h =10 м и радиусе контура питания R_к = 500 м дебит горизонтальной скважины длиной L = 50 м превышает дебит вертикальной скважины в 2 раза, а при длине L = 500 м - в 5,9 раз. Таким образом, увеличение длины горизонтальной скважины в 10 раз приводит к увеличению дебита всего лишь в 2,95 раза. Это обстоятельство должно учитываться при технико-экономическом обосновании длины горизонтальной скважины в каждом конкретном случае проектирования.

Рассмотрим влияние  на эффективность горизонтальной скважины коэффициента анизотропии пласта по проницаемости .

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

Вместо формулы  Дюпюи для вертикальной скважины будем иметь

 

,    (1.15)

где           (1.16)

- эквивалентный радиус скважины.

Расчеты отношения дебитов q_г/q_в для анизотропных пластов представлены в таблице 3.

Толщина пласта принята равной h =10 м, радиус контура питания R_к = 500 м, длина скважины L = 300 м.

Расчеты выполнены для отношения  проницаемостей K_г/K_в = 10;5;2;1;0,5;0,2;0,1.

Для расчета дебита горизонтальной скважины использовались формулы (1.6), (1.7), (1.8) и (1.13).

 

 
 Поскольку формула Renard, Dupuy (1.7) является наиболее точной из них, достаточно проследить за изменением отношения q_г/q_в в зависимости от значения K_г/K_в на основе этой формулы.

Для изотропного пласта (K_г/K_в =1) отношение дебитов равно q_г/q_в = 4,28.

Для анизотропного пласта для случаев K_г/K_в > 1 отношение дебитов q_г/q_в возрастает, а при K_г/K_в < 1 - убывает.

Так, например, для отношения K_г/K_в = 10 отношение дебитов равно 11,84 , для случая K_г/K_в = 0,1 составляет 1,38.

Анализ показывает, что столь  резкое влияние анизотропии на отношение  дебитов горизонтальной и вертикальной скважин обусловлено в большей мере сильной зависимостью дебита вертикальной скважины от показателя анизотропии по формуле (1.15). Дебит же горизонтальной скважны по формуле (1.7) зависит от анизотропии не так существенно.

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

Из данных таблицы 3 следует основной вывод - эффективность горизонтальной скважины по отношению к вертикальной в анизотропном пласте выше, чем  в изотропном при K_г/K_в>1 и ниже при K_г/K_в < 1. Поэтому учет влияния анизотропии является обязательным при технико-экономическом обосновании варианта разработки.

Были выполнены расчеты и  по оценке влияния асимметричности  расположения горизонтальной скважины относительно кровли и подошвы пласта.

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

 

4. создание инженерной методики расчета

 нестационарного  притока в горизонтальную скважину

 

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