СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДА В МНОГОПЛАСТОВОЙ СКВАЖИНЕ
METHOD FOR DETERMINING PROFILE OF FLUID INFLUX IN MULTI-PAY WELL

Изобретение относится к области геофизических исследований нефтяных и газовых скважин, а именно к определению профиля притока добываемого флюида в многопластовых скважинах с несколькими интервалами перфорации.
Определение профиля притока из многозонной скважины является важной задачей. Определение дебита отдельных интервалов перфорации необходимо, в частности, для принятия решения о необходимости проведения кислотной обработки, повторной перфорации и т.д..
Определение профиля притока обычно проводят во время промыслового каротажа добывающей скважины с помощью механических расходомеров (см., например, Hill, A.D.,. Production Logging - Theoretical and Interpretive Elements, SPE Monograph Series., 2002, стр. 61). Основными недостатками этого способа являются необходимость проведения специального каротажа скважины (в дополнение к операциям, проводимым в скважине во время перфорации и опробования скважины) и сложность определения дебитов малопродуктивных пластов.
Вклад различных интервалов перфорации может быть оценен также с помощью данных температурного каротажа добывающей скважины (см. Череменский Г.А., Прикладная геотермия, М. Недра, 224 стр., стр. №181) или из анализа нестационарных температурных данных, полученных при изменении дебита скважины (см. Чекалюк, Е.Б., Термодинамика нефтяного пласта, Москва, 1965, 234 стр., стр. №88, или Ramazanov, A., Valiullin, R.А., Shako, V., Pimenov, V., Sadretdinov, A., Fedorov, V., Belov, K., 2010. Thermal Modeling for Characterization of Near Wellbore Zone and Zonal Allocation, SPE 136256-MS). К недостаткам этих способов можно отнести необходимость анализа относительно небольших температурных сигналов и необходимость проведения специальных каротажей скважины или установки в скважине специального оборудования.
В соответствии с предлагаемым способом осуществляют измерения забойной температуры и забойного давления в скважине посредством датчиков, установленных на перфорационной колонне ниже всех интервалов перфорации, а также датчиков температуры, установленных на перфорационной колонне выше каждого интервала перфорации.
Измерения температуры и забойного давления осуществляют до проведения перфорации скважины и после перфорации до тех пор, пока температура добываемого флюида не вернется к первоначальной температуре пласта. Оценивают суммарный дебит скважины и рассчитывают избыточную тепловую энергию добываемого флюида для всех датчиков температуры, установленных на перфорационной колонне выше интервалов перфорации, после чего определяют дебит отдельных интервалов перфорации на основе рассчитанных избыточных тепловых энергий добываемого флюида и известного количества перфорационных зарядов в каждом интервале перфорации.
В соответствии с одним из вариантов осуществления изобретения суммарный дебит скважины определяют посредством измерения расхода на поверхности или в скважине.
В соответствии с другим вариантом осуществления изобретения суммарный дебит скважины определяют посредством расчета расхода по изменению забойного давления.
В соответствии с еще одним вариантом осуществления изобретения суммарный дебит скважины определяют посредством расчета расхода с использованием забойного давления и численного моделирования многопластовой добывающей скважины.
0
1
2
2
2
Изобретение поясняется чертежами, где на фиг. 1 показана схема скважины с двумя интервалами перфорации, на фиг. 2 приведен пример забойного давления P(t) и температур добываемого флюида T(t) и T(t) выше интервалов перфорации, на фиг. 3 показан дебит скважины, рассчитанный для давления, приведенного на Фиг. 2, на фиг. 4 приведена полная избыточная тепловая энергия добываемого флюида (рассчитанная по температуре Т, сплошная линия) и соответствующая энергия, рассчитанная по температуре T, на фиг. 5 показан алгоритм определения профиля притока с использованием численного моделирования многопластовой добывающей скважины.
Данное изобретение предлагает определять профиль притока в скважинах с несколькими интервалами перфорации с использованием результатов измерения скважинного давления и результатов измерения температуры с помощью датчиков, установленных на перфорационной колонне. Температуру надо измерять выше каждого интервала перфорации и на забое скважины, ниже всех интервалов перфорации.
0
0
i
Способ предусматривает измерение забойного давления P(t) и забойной температуры T(t), которая определяет среднюю температуру пород в рассматриваемом интервале глубин. Измерения осуществляют с помощью датчиков, установленных на перфорационной колонне в скважине ниже всех интервалов перфорации, а также измерения температуры T(t) добываемого флюида (i=1,2,..,m, m есть число интервалов перфорации) с помощью датчиков температуры, установленных на перфорационной колонне выше каждого интервала перфорации.
0
i
Измерения давления P(t) и температуры T(t) (i=0,1,..,m) начинают до перфорации (что позволяет определить пластовое давление и геотермальную температуру и продолжать измерение в течение нескольких часов после перфорации, до тех пор, пока температура добываемого флюида, нагретого благодаря энергии перфорационного взрыва, не вернется к первоначальной температуре пласта). При взрыве перфорационных зарядов часть энергии идет на испарение скважинного флюида и на энергию кумулятивной струи, но большая часть энергии идет на нагрев перфорационной колонны, обсадной трубы и породы вблизи скважины. Нагрев добываемого флюида происходит при его контакте с этими телами.
Затем оценивают суммарный дебит скважины Q(t), используя один из следующих способов:
- измерение расхода на поверхности или в скважине,
0
- расчет расхода по изменению забойного давления P(t) (если добываемый флюид не достигает поверхности),
0
- расчет расхода с использованием забойного давления P(t) и численного моделирования многопластовой добывающей скважины.
Параметры (проницаемости и скин факторы), определяющие продуктивность отдельных пластов, принимаются равными средним значениям, которые определяются в результате традиционного гидродинамического исследования скважины.
Рассчитывают избыточную тепловую энергию добываемого флюида для каждого температурного датчика
f
0
где T - средняя температура пород в рассматриваемом интервале глубин (определяемая T(t) и практически равная ей),  - объемная теплоемкость флюида.
i
Дебит отдельных интервалов перфорации рассчитывают по величинам E и известным количествам перфорационных зарядов в каждом перфорационном интервале.
Рассмотрим случай малодебитной скважины, когда в первые часы после перфорации нет излива добываемого флюида на поверхность.
2
1
0
0
Схема скважины с перфорационной колонной, пакером и двумя интервалами перфорации приведена на Фиг. 1, где показаны пакер - 1, клапан - 2, датчик температуры Т - 3, датчик температуры Т - 4, датчики забойной температуры и давления Т, Р - 5, вторая зона притока - 6, первая зона притока - 7, второй интервал перфорации - 8, первый интервал перфорации - 9.
0
1
2
На Фиг. 2 приведен синтетический пример забойного давления P(t) и температур добываемого флюида T(t) и Т(t) выше интервалов перфорации. Толстая кривая соответствует забойному давлению, которое равно ~50 бар перед перфорацией и увеличивается до пластового давления (около 85 бар) во время добычи в соответствии с тем, что поднимается уровень флюида в добывающей трубе. В данном случае принято, что интервалы перфорации имеют одинаковую протяженность и одинаковое количество перфорационных зарядов.
0
Если нет излива добываемого флюида на поверхность, суммарный дебит скважины Q(t) может быть рассчитан по забойному давлению P(t):
t
2
где r - внутренний радиус трубы, g=9.81 - м/с ускорение свободного падения,  - плотность флюида.
3
t
На Фиг. 3 показан дебит скважины, рассчитанный по этой формуле для давления, приведенного на Фиг. 2 (для =850 кг/м, r=0.038 м). Рассчитанный дебит далее используют для определения профиля притока.
В случае, если суммарный дебит скважины измерялся в скважине или на поверхности, этот дебит непосредственно используется для определения профиля притока.
1
2
f
Графики температуры T и Т (Фиг. 2) показывают, что сразу после перфорации температура потока добываемого флюида значительно больше (в данном случае на ~20 С), чем температура пород T (точки на Фиг. 2). Эта температура определяется нагревом скважинного флюида при взрыве и нагревом пластового флюида при его контакте с горячей породой, обсадной колонной и перфорационной колонной. Следует отметить, что температура породы может быть оценена по результатам измерения температуры в скважине перед перфорацией.
p
m
Поток пластового флюида охлаждает околоскважинную породу, обсадную и перфорационную колонну и через некоторое время (t=5÷10 час) после перфорации измеренные в скважине температуры приближаются к невозмущенной температуре пород (Фиг. 2). Это означает, что тепловая часть Е энергии взрыва перфорационных зарядов трансформировалась в избыточную тепловую энергию добытого флюида.
m
2
2
В данном случае m=2 и E≡E. Используя температуру Т, измеренную датчиком, который расположен выше всех интервалов перфорации, и дебит скважины Q(t), эту энергию можно рассчитать по формуле:
2
2
Сплошная линия на Фиг. 4 показывает избыточную тепловую энергию добываемого флюида для данных, приведенных на Фиг. 2. Видно, что через ~3 часа после перфорации Е достигает своего наибольшего значения Е≈16.5 МДж.
е
Полная энергия перфорационного взрыва, рассчитанная по удельной энергии взрыва и массе взрывчатого вещества, в рассматриваемом случае составляет Е≈28 МДж. Это означает, что приблизительно δ=60% от энергии взрыва было преобразовано в тепловую энергию породы, обсадной и перфорационной колонны:
m
e
E=δ⋅E.
Оставшаяся часть энергии взрыва (около 40%) была затрачена на разрушение породы, генерацию ударных волн в породе и в скважине или была быстро вынесена за пределы рассматриваемого интервала с газообразными продуктами взрыва.
Предлагаемая в данном изобретении процедура расчета дебита отдельных интервалов перфорации основана на следующих предположениях:
- величина δ одинакова для разных интервалов перфорации,
- флюиды, поступающие в скважину из разных интервалов перфорации, имеют одинаковые объемные теплоемкости,
- расстояние между интервалами перфорации невелико и можно пренебречь потерями тепловой энергии флюида в окружающие породы между интервалами перфорации,
- продолжительность добычи после перфорации и дебиты скважины достаточно велики, так что измеряемая датчиками температура флюида снижается до температуры невозмущенных пород.
Пусть m - число интервалов перфорации,
i
го
Q есть дебит из i интервала,  - суммарный дебит скважины,
m
есть дебит скважины из нижних i перфорационных интервалов, отнесенный к суммарному дебиту скважины (γ=1),
i
м
n есть число перфорационных зарядов в i интервале перфорации,
есть полное число перфорационных зарядов в скважине,
m
есть число зарядов в нижних i интервалах перфорации, отнесенное к полному число перфорационных зарядов в скважине (b=1),
i
го
Т(t) есть температура флюида, измеренная датчиком температуры, расположенным выше i интервала перфорации.
i
Дебит отдельных интервалов перфорации (на начальном этапе значения γ) рассчитывают с использованием закона сохранения энергии, который записывают для всех интервалов (i=1,2,..m):
или
где i=1,2,..m,
i
i
i
Искомые относительные продуктивности у (y=Q/Q, ) отдельных перфорационных интервалов рассчитывают по формулам:
1
1
В рассматриваемом случае двух интервалов перфорации (m=2) и одинакового числа перфорационных зарядов в интервалах (b=0.5) рассчитанная энергия Е(t) показана на Фиг. 4 пунктирной линией.
1
1
1
Расчетное значение безразмерного дебита γ(t) выходит на приблизительно постоянное значение через ~3 часа после перфорации: γ=y≈0.7.
0
i
i
В общем случае нестационарный дебит скважины Q(t) может быть рассчитан с использованием измеренного забойного давления P(t) и численной модели многопластовой добывающей скважины, в которую в качестве свободных параметров входят проницаемости {k} и скины {s} продуктивных пластов. Значения этих параметров могут быть найдены с использованием итерационной процедуры, алгоритм которой приведен на Фиг. 5.
i
i
0
ki
i
Первоначальный набор параметров, характеризующих продуктивные интервалы {k, s}, определяется с помощью традиционного гидродинамического исследования (ГДИ) скважины в предположении, что все пласты имеют одни и те же свойства. Для этих параметров с использованием измеренного забойного давления P(t) рассчитывают общий дебит скважины Q(t) и относительные дебиты отдельных пластов {y}. Затем, используя найденный дебит Q(t) и температуры {T(t)}, измеренными датчиками, расположенными выше продуктивных пластов, с использованием описанной выше процедуры находят относительные дебиты  и сравнивают два полученных набора чисел, характеризующих профиль притока, например, рассчитывают величину невязки S:
i
i
ki
Если S меньше заданной величины невязки ε: S<ε, то данный набор параметров принимается в качестве решения задачи. В противном случае значения параметров {k,s} изменяют, и вычисления продолжают до тех пор, пока векторы {y} и  с заданной точностью не совпадут.