Способ определения профиля притока флюидов многопластовых залежей в скважине

Изобретение относится к области геофизических исследований нефтяных и газовых скважин, а именно к определению профиля притока флюидов, поступающих в скважину из продуктивных пластов многопластовых коллекторов. Технический результат настоящего изобретения заключается в увеличении точности и достоверности определения профиля притока в многопластовой скважине на начальной стадии добычи, сразу после перфорации скважины. При этом нет необходимости в выстойке скважины между промывкой скважины и перфорацией, и не требуется измерять скорость изменения температуры в скважине до ее перфорации. Технический результат достигается тем, что осуществляют охлаждение призабойной зоны перед перфорацией скважины. Проводят перфорацию скважины и измеряют температуру потока в скважине над каждой зоной перфорации. Определяют дебит каждого продуктивного пласта, учитывая толщины зон перфорации и используя результаты измерений температуры, полученные в интервале между окончанием начальной стадии добычи, характеризующейся сильным влиянием объема скважины и быстрым изменением температуры потока в скважине, и временем, начиная с которого влияние охлаждения призабойной области скважины на измерения температуры становится незначительным. 7 з.п. ф-лы, 10 ил.

 

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

Обычно при оценке дебита отдельных продуктивных пластов по температурным данным производятся измерения температуры по всему стволу скважины при квазистационарных условиях добычи, а температура коллектора вблизи скважины считается близкой к температуре невозмущенного коллектора. При этом предполагается, что температура пород вблизи скважины равна невозмущенной температуре горного массива. Для определения дебита продуктивного пласта обычно используют результаты измерения температуры потока ниже и выше продуктивного пласта и температуру флюида, поступающего из пласта в скважину (см., например, Hill, A.D., 2002. Production Logging - Theoretical and Interpretive Elements. SPE Monograph Series).

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

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

В патенте РФ №2474687 описан способ определения профиля притока многопластовых залежей, использующий результаты измерения температур в скважине сразу после перфорации, когда температура пород вблизи скважины понижена предшествующими технологическими операциями (бурение, промывка скважины). В соответствии с этим способом дебиты отдельных пластов рассчитывают по измеренным в скважине скоростям изменения температур. Теоретической основой этого способа является установленная в изобретении линейная зависимость между скоростью изменения температуры добываемого флюида и удельным дебитом пласта, которая предполагает охлаждение призабойной зоны пласта и выстойку скважины между циркуляцией в скважине флюида и перфорацией.

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

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

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

Охлаждение призабойной зоны может быть осуществлено путем бурения или путем промывки скважины.

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

Начиная со второй зоны перфорации, последовательно снизу вверх рассчитывают относительные дебиты зон перфорации Yn,(n=2,3,…m) путем минимизации

S ( Y n ) = i F ( Y n , t i ) 2 min

где Y n = Q n Q 1 + Q 2 + .. + Q n , Qn(n=2,3,…m) - дебиты зон перфорации

F ( Y n , t ) = Y n a T n 1, a ( t ) T n , a ( t ) T n 1, a ( t ) T 1, a [ a n t ]

,

hn(n=1,2,…m) - толщина nй зоны перфорации, Tn,a(t) - температуры потока, измеренные в скважине выше nй зоны перфорации, ti - моменты времени, равнораспределенные в интервале t1<t<t2, где t1 - время окончания начальной стадии добычи, характеризующейся сильным вилянием объема скважины и быстрым изменением температуры потока в скважине, t2 - время, начиная с которого»влияние охлаждения призабойной области скважины на измерения температуры становится незначительным, причем количество моментов времени в этом интервале соответствует числу измерений температуры.

Дебиты Qn(n=2,3,…m) зон перфорации определяют по формулам:

Q 1 = Q 1 + Y 2 1 Y 2 + Y 3 ( 1 Y 3 ) ( 1 Y 3 ) + + Y m ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y m )

Q 2 = Q 1 Y 2 1 Y 2

Q 3 = Q 1 Y 3 ( 1 Y 2 ) ( 1 Y 3 )

Q n = Q 1 Y n ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y n )

где Q - суммарный объемный дебит скважины.

В соответствии с одним из вариантов реализации изобретения дополнительно измеряют температуру потока в скважине Tn,b(t) под каждой зоной перфорации, при этом F(Yn,t) определяют как

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

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

Изобретение поясняется чертежами, где

на фиг.1 показана схема скважины с двумя зонами перфорации и тремя датчиками температуры;

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

на фиг.2б приведены расчетные температуры флюидов, добываемых из зон с различными удельными дебитами (сплошные линии), и температура зоны с большим дебитом после масштабирования шкалы времени (маркеры);

на фиг.3 и фиг.4 приведены расчетные температуры для различных вариантов значений проницаемости пластов;

на фиг.5 и фиг.6 приведены зависимости невязки от относительного дебита пласта для различных вариантов значений проницаемости пластов;

на фиг.7а приведены расчетные температуры с нормальным шумом СКО=0.1 К. Вариант k1=100 мД, k2=30 мД.

на фиг.7б приведена зависимость невязки S от величины Y2. Решение Y2min=0.2.

на фиг.8а приведены расчетные температуры с нормальным шумом СКО=0.1 K. Вариант k1=30 мД, k2=100 мД.

на фиг.8б приведена зависимость невязки S от величины Y2. Решение Y2min=0.69.

на фиг.9а приведены значения подбираемых проницаемостей при решении обратной задачи. Вариант k1=100 мД, k2=30 мД.

на фиг.9б приведены заданные температуры и температуры, рассчитанные при решении обратной задачи. k1=96 mD, k2=29 mD, Y2≈0.23

на фиг.10а приведены значения подбираемых проницаемостей при решении обратной задачи. Вариант k1=30 мД, k2=100 мД.

на фиг.10б приведены заданные температуры и температуры, рассчитанные при решении обратной задачи. k1=29 mD, k2=96 mD, Y2≈0.75

Предлагаемый метод может быть использован в условиях типичной перфорации с использованием насосно-компрессорной трубы (НКТ).

Обычно охлаждение призабойной зоны происходит при бурении скважины. Если бурение было закончено недавно (за несколько несколько дней до перфорации), при этом не требуется проводить специальной промывки скважины Если скважина была пробурена давно, то перед перфорацией следует провести специальную промывку скважины, что обеспечит существенное (на 5-10 градусов Цельсия и более) снижение температуры околоскважинных по сравнению с невозмущенной температурой.

Предлагаемый способ использует то обстоятельство, что начальное радиальное распределение температуры приблизительно одинаково во всех рассматриваемых продуктивных пластах T≈T0(r).

Это предположение справедливо при выполнении следующих условий:

(a) продуктивные пласты расположены близко друг к другу (на расстоянии 30-50 м)

(b) все продуктивные пласты имеют приблизительно одинаковые тепловые свойства.

где

есть скорость фильтрации флюида, q[м3/м/с] - дебит, приходящийся на 1 м зоны перфорации, ρfcf - объемная теплоемкость флюида, ρrcr=ϕ·ρfcf+(1-ϕ)·ρmcm - объемная теплоемкость флюидонасыщенного резервуара, ρmcm - объемная теплоемкость матрицы горной породы, ϕ - пористость резервуара, r - расстояние до оси скважины..

В уравнении (1) не учитывается влияние кондуктивной теплопередачи, эффекта Джоуля-Томсона и адиабатического эффекта. Малое влияние кондуктивной теплопередачи в зоне притока при типичных удельных дебитах было показано путем детального численного моделирования. Относительное влияние эффекта Джоуля-Томсона и адиабатического эффекта определяется перепадом давления между пластом и скважиной и типичным понижением температуры пород в призабойной области. Поскольку на начальной стадии добычи перепады давления обычно невелики (10-30 Бар), а захолаживание пласта достигает 10 K и более, в первом приближении этими эффектами можно пренебречь. Более точное определение профиля притока по данным нестационарной термометрии может быть проведено с помощью численного моделирования (см. ниже).

Решение уравнения (1) имеет вид:

T ( r , t ) = T 0 ( r 2 + χ π q t ) ,               ( 3 )

где , T0(r) - радиальное распределение температуры в пласте перед началом добычи, q - удельный дебит пласта.

Из (3) следует выражение для зависимости от времени температуры поступающего в скважину флюида:

T i n ( t ) = T 0 ( r w 2 + χ π q t ) ,               ( 4 )

где rw - радиус скважины.

Рассмотрим скважину с двумя продуктивными зонами (Фиг.1). Температуры флюидов, поступающих из разных пластов с удельными дебитами q1 и q2, равны соответственно:

T i n ,1 ( t ) = T 0 ( r w 2 + χ π q 1 t )

и

T i n ,2 ( t ) = T 0 ( r w 2 + χ π q 2 t )

Очевидно, что температура Tin,2(t) может быть записана как T i n ,2 ( t ) = T 0 ( r w 2 + χ π q 1 t q 2 q 1 )

или

где

Согласно формуле (6), относительный дебит верхней продуктивной зоны Y2 может быть найден путем подбора такой величины масштабирующего коэффициента 'a 2', при котором будут совпадать зависящие от времени температуры флюидов, поступающих из разных продуктивных пластов.

Это утверждение иллюстрируют на Фиг.2 результаты численных расчетов, проведенных с использованием коммерческого пакета COMSOL Multiphysics 3.5. Расчетные радиальные распределения температуры после промывки скважины (24 часа промывки при начальной температуре породы 100°C и температуре циркулирующего флюида 50°C, пунктирная кривая и ее выстойки в течение 24 часов (сплошная кривая) показаны на Фиг.2а. Тонкая сплошная кривая на Фиг.2б показывает расчетную температуру поступающего в скважину флюида при удельном дебите q (q=10 м3/м/день), толстая сплошная кривая соответствует в два раза большему удельному дебиту 2·q. Маркеры, которые практически совпадают с тонкой кривой, показывают результат растяжения шкалы времени для кривой, соответствующей дебиту 2·q в два раза.

Формулы (5), (6) позволяют найти дебиты отдельных зон перфорации, если для всех зон измерены зависящие от времени температуры флюидов, поступающих из разных пластов в скважину. На практике измеряют только температуру потока в скважине. Температура, измеренная выше нижней зоны перфорации, приблизительно равна температуре флюида, поступающего из нижней зоны, но для всех остальных зон температура Tin обычно неизвестна.

В изобретении предлагается следующий способ решения этой проблемы.

Рассмотрим случай, при котором датчики температуры расположены на НКТ выше и ниже каждой зоны перфорации (Фиг.1).

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

где ρc - объемная теплопроводность флюида.

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

Из уравнения (7) следует формула для относительного дебита верхнего пласта Y2:

Учитывая, что

T1,a(t)≈Tin,1(t)

и используя формулу (5)

Tin,2(t)=Tin,1(a 2·t)≈T1,a(a 2·t),

находим уравнение относительно искомой величины Y2:

F ( Y 2 , t ) = 0                                     ( 9 )

где

Существенно, что в это уравнение не входит неизвестная температура добываемого флюида Tin2(t) и неизвестная температура невозмущенных горных пород Tf. Искомая величина Y2 определяется исключительно результатами измерения температуры в скважине: T1,a, T2,b, T2,a.

Если расстояние между пластами мало (порядка 10 м), то, как показывают расчеты, T2,b(t)≈T1,a(t) и соотношение (10) примет вид:

Относительный дебит верхнего пласта может быть найден из уравнения (9) для любого фиксированного момента времени 't', однако не все моменты времени эквивалентны.

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

При больших временах (t>t2), когда перестает сказываться влияние охлаждения призабойной зоны и флюид, поступающий в скважину, прогреет околоскважинные породы практически до температуры резервуара, становится более существенным влияние эффекта Джоуля-Томсона и влияние геотермы. Оба эти эффекта не учитываются описанной выше упрощенной моделью, что может привести к ошибкам при определении профиля притока. Поэтому при использовании описанной выше упрощенной модели определения профиля притока надо использовать данные, полученные при относительно небольших временах t<t2.

Кроме того, все величины T1a, T2b, и T2a измеряются с определенной погрешностью. Эта погрешность связана как с погрешностью датчиков, так и с погрешностью метода измерения. Действительно, в расчетную формулу (10) входят среднемассовые температуры потока, в то время как датчик измеряет температуру в определенной точке потока и имеется некоторое распределение температуры по поперечному сечению потока. Поэтому для уменьшения влияния различных факторов на результаты определения профиля притока целесообразно использовать все данные, полученные в интервале времени t1<t<t2.

В изобретении предлагается находить величину Y2 путем минимизации функции S(Y2) (11), которая рассчитывается как сумма квадратов невязок для всех моментов времени в интервале t1<t<t2, для которых имеются измеренные значения температуры:

S ( Y 2 ) = F ( Y 2 , t i ) 2 min                           ( 11 )

Предложенный способ определения профиля притока был протестирован на ряде типичных синтетических случаев. Было показано, что значения t1=1÷2 часа и t2=8÷10 часов обеспечивают удовлетворительную точность решения обратной задачи.

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

Начиная со второй зоны перфорации, последовательно снизу вверх рассчитывают относительные дебиты зон перфорации Yn,(n=2,3,…m) путем минимизации

S ( Y n ) = i F ( Y n , t i ) 2 min                      ( 12 )

где Y n = Q n Q 1 + Q 2 + .. + Q n , Qn(n=2,3,…m) - дебиты зон перфорации 0

hn(n=1,2,…m) - толщина nй зоны перфорации, Tn,a(t) - температура флюида, измеренная в скважине выше nй зоны перфорации, ti - моменты времени, равнораспределенные в интервале t1<t<t2, t1 - время окончания начальной стадии добычи, характеризующейся сильным вилянием объема скважины и быстрым изменением температуры скважинного флюида (по порядку величины это время определяется временем замены флюида в скважине в области продуктивных пластов), t2 - время, начиная с которого влияние охлаждения призабойной области скважины на измерения температуры становится незначительным (того же порядка, что и влияние геотермы и эффекта Джоуля-Томсона), причем количество моментов времени в этом интервале соответствует числу измерений температуры.

Дебит Qn(n=2,3,…m) зон перфорации определяют по формулам:

Q 1 Q 1 + Y 2 1 Y 2 + Y 3 ( 1 Y 3 ) ( 1 Y 3 ) + + Y m ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y m )                ( 15 )

Q 2 = Q 1 Y 2 1 Y 2

Q 3 = Q 1 Y 3 ( 1 Y 2 ) ( 1 Y 3 )

Q n = Q 1 Y n ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y n )

где Q - суммарный объемный дебит скважины.

Дебиты отдельных продуктивных пластов, определенные с помощью формул (12)-(15), можно рассматривать как приближенное решение поставленной задачи. В некоторых случаях, когда нарушаются обсуждавшиеся выше условия применимости аналитической модели, ошибка в определении индивидуальных дебитов может достигать 10-20% и более (см. ниже).

Возможность определения профиля притока с использованием предлагаемого способа была продемонстрирована на синтетических тестовых случаях, которые были генерированы с помощью нестационарной численной модели сопряженных процессов тепло- и массопереноса в пластах и в скважине. Функциональность модели позволяет воспроизводить произвольное чередование различных технологических процессов: промывку скважины, ее выстойку, добычу, нагнетание в скважину флюида. Численная модель скважины прошла многолетнее тестирование на аналитических решениях, с использованием коммерческих симуляторов (Comsol и Eclipse), успешно использовалась при моделировании и интерпретации сложных полевых случаев.

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

1. Промывка скважины в течение 24 часов. Температура циркулирующего в скважине флюида на глубине пласта принималась равной 100°C, невозмущенная температура пород - 123°C.

2. Выстойка скважины в течение 24 часов. Это этап не является обязательным для использования предлагаемого способа, но обычно существует некоторый интервал времени между промывкой и перфорацией скважины.

3. Перфорация скважины и добыча нефти в течение 24 часов с суммарным дебитом 100 м3/день.

Толщина зон перфорации h1=h2=10 м, расстояние между зонами 10 м.

Были рассмотрены два случая значений проницаемости пластов (скин-фактор принимался равным нулю):

( 1 )   k 1 = 100  мД , k 2 = 30  мД ( Y 0 .25 )

( 2 )   k 1 = 30  мД , k 2 = 100  мД ( Y 0 .75 ) .

В скобках даны относительные дебиты верхнего пласта в рассматриваемых случаях.

Фигуры 3, 4 показывают расчетные температуры T1,a, T2,b, T2,а и Tin,2 для этих вариантов. Температура T1,a равна температуре потока в точке, находящейся на 1 м выше верхней границы зоны перфорации, T2,b - на 1 м ниже нижней границы зоны перфорации. Видно, что в данном случае разница между температурами T1,a и T2,b очень мала, так как дебит скважины достаточно большой, пласты находятся близко друг к другу и потери энергии потока флюидов в окружающие породы малы. Это означает, что в данном случае нет необходимости в датчике температуры, находящемся ниже верхнего пласта (T2,b), и для функции F(Y2,t) может быть использована формула (11).

Расчеты показывают, что во всех случаях невязка S(Y2) (11) может быть рассчитана для интервала времени с 3 до 7 часов. Фигура 5 показывает, что в случае №1 невязка S(Y2) имеет единственный минимум при Y2≈0.233. Это значение очень близко к точному решению Y2=0.25 (погрешность ~7%). В случае №2 (Фиг.6) невязка S(Y2) также имеет единственный минимум при Y2min=0.68 (погрешность ~9%).

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

В некоторых случаях точность предлагаемого упрощенного способа определения профиля притока по температурным данным может быть недостаточна. Это может быть связано с нарушением условий, которые использовались при выводе расчетных формул. Например, при большом расстоянии между пластами будут различаться начальные радиальные распределения температуры в пластах, большое влияние будет оказывать геотермический градиент, при больших перепадах давления между пластом и скважиной нельзя будет пренебречь влиянием эффекта Джоуля-Томсона, при малых временах наблюдения нельзя будет пренебрегать влиянием объема скважины и т.д.

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

Фигуры 9, 10 иллюстрируют использование этого способа определения профиля притока. Для приведенных выше двух синтетических случаев автоматически подбирались проницаемости продуктивных пластов, обеспечивающие совпадение заданных и расчетных температур. Значения проницаемостей, полученные при использовании упрощенной модели, использовались как начальные приближения. Из Фиг.9, 10 видно, что практически полное совпадение температур было получено после приблизительно 150 решений прямой задачи. При этом найденные значения проницаемостей и дебитов практически совпадают с заданными:

( 1 )   k 1 = 96  мД , k 2 = 29  мД ( Y 2 0 .23 ) .

Заданные значения: k1=100 мД, k2=30 мД (Y2=0.25)

( 2 )   k 1 = 29  мД , k 2 = 96  мД ( Y 2 0 .75 ) .

Заданные значения: k1=30 мД, k2=100 мД (Y2=0.75).

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

1. Способ определения профиля притока флюидов многопластовых залежей в скважине, в соответствии с которым:
осуществляют охлаждение призабойной зоны,
проводят перфорацию скважины,
измеряют температуру потока в скважине над каждой зоной перфорации,
и определяют дебит каждого продуктивного пласта, учитывая толщины зон перфорации и используя результаты измерений температуры, полученные в интервале между окончанием начальной стадии добычи, характеризующейся сильным влиянием объема скважины и быстрым изменением температуры потока в скважине, и временем, начиная с которого влияние охлаждения призабойной области скважины на измерения температуры становится незначительным.

2. Способ по п.1, в соответствии с которым охлаждение призабойной зоны осуществляется за счет процесса бурения.

3. Способ по п.1, в соответствии с которым охлаждение призабойной зоны осуществляют путем промывки скважины.

4. Способ по п.1, в соответствии с которым температуру потока в скважине определяют с помощью датчиков, установленных на колонне НКТ, используемой для перфорации, над каждой зоной перфорации.

5. Способ по п.1, в соответствии с которым для определения дебита продуктивных пластов используют результаты измерений температуры, полученные предпочтительно в интервале от 1 до 10 часов с момента начала добычи.

6. Способ по п.1, в соответствии с которым, начиная со второй зоны перфорации, последовательно снизу вверх рассчитывают относительные дебиты зон перфорации Yn, (n=2,3,… m), путем минимизации функции S(Yn), которая рассчитывается как сумма квадратов невязок для всех моментов времени в интервале t1<t<t2, для которых имеются измеренные значения температуры: S ( Y n ) = i F ( Y n , t i ) 2 min
где Y n = Q n Q 1 + Q 2 + .. + Q n , Qn(n=2,3,…m) - дебиты зон перфорации
, -масштабирующий коэффициент, равный:
,
qn (n=l,2,…m) - удельный дебит nого пласта, hn(n=l,2,…m) - толщина nй зоны перфорации, Tn,a(t) (n=l,2,…m) - температуры потока, измеренные в скважине выше nй зоны перфорации, ti - моменты времени, равнораспределенные в интервале t1<t<t2, где t1 - время окончания начальной стадии добычи, характеризующейся сильным влиянием объема скважины и быстрым изменением температуры потока в скважине, t2 - время, начиная с которого влияние охлаждения призабойной области скважины на измерения температуры становится незначительным, причем количество моментов времени в этом интервале соответствует числу измерений температуры, и
определяют дебит Qn(n=2,3,…m) зон перфорации по формулам:
Q 1 Q 1 + Y 2 1 Y 2 + Y 3 ( 1 Y 2 ) ( 1 Y 3 ) + + Y m ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y m )
Q 2 = Q 1 Y 2 1 Y 2
Q 3 = Q 1 Y 3 ( 1 Y 2 ) ( 1 Y 3 )
Q n = Q 1 Y n ( 1 Y 2 ) ( 1 Y 3 ) ( 1 Y n )
где Q - суммарный объемный дебит скважины.

7. Способ по п.1, в соответствии с которым дополнительно измеряют температуру потока в скважине под каждой зоной перфорации, при этом F(Yn, t) определяют как

где Tn,b(t) и Tn,a(t) - температуры, измеренные в скважине ниже и выше nй зоны перфорации.

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



 

Похожие патенты:

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

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

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

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

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

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

Изобретение относится к шлангокабелям, предназначенным для работ в нефтяных и газовых скважинах и может быть использовано для перемещения предметов, в частности приборов в горизонтальных скважинах.

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

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

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

Изобретение относится к области газового машиностроения, в частности к устройствам исследования газовых и газоконденсатных месторождений на разных технологических режимах. Технический результат заключается в снижении массогабаритных характеристик устройства, выполнении его транспортабельным, для перевозки и размещения на автомобильном транспорте или прицепе без предварительной разборки и последующей сборки и исключение необходимости гидроиспытаний перед проведением исследований, что значительно сокращает время подготовительных работ перед измерением. Блок для исследования газовых и газоконденсатных скважин включает сепаратор газожидкостной смеси, устройства: сужающее, замера продукции сепарации, замера давления, температур и расхода газа, сбора жидкостей и механических примесей и запорно-регулирующую арматуру. Блок расположен в каркасе с размерами, вписывающимися в габариты транспортных средств, при этом сепараторы, и устройства, входящие в состав блока, в число которых дополнительно введен сверхзвуковой сепаратор, обвязаны трубопроводами. Сепаратор газожидкостной смеси по входу соединен с выходом(ами) сверхзвукового сепаратора, а выход очищенного газа сверхзвукового сепаратора соединен с сепаратором газожидкостной смеси или с трубопроводом выхода очищенного газа. 5 з.п. ф-лы, 2 ил.

Изобретение относится к нефтяной и газовой промышленности и используется для оптимизации процесса добычи нефти с помощью штанговых глубинных насосов. Техническим результатом является вывод скважины в автоматическом режиме на максимальный объем добычи нефти. Способ оптимизации процесса добычи нефти, в котором непрерывно контролируют параметры процесса добычи и регулируют частоту вращения электродвигателя насоса, причем в качестве основного контролируемого параметра процесса добычи выбирают коэффициент заполнения насоса, который рассчитывают по динамограмме, в качестве регулируемого параметра, подлежащего оптимизации, выбирают объем добываемой жидкости, а оптимизацию осуществляют за счет ступенчатого изменения частоты качаний насоса, при этом процесс управления осуществляют ступенчато в автоматическом режиме. Устанавливают частоту вращения электродвигателя, при которой частота качаний насоса много меньше номинальной, определяют площадь динамограммы работы качалки при максимальном наполнении насоса, значение которой принимают за базовое. После этого увеличивают частоту качаний насоса, непрерывно с заданной дискретностью фиксируют площадь динамограммы, а коэффициент заполнения насоса k3 рассчитывают как отношение измеренной площади динамограммы к базовой, при этом частоту качаний насоса не изменяют до тех пор, пока коэффициент заполнения насоса не перестанет изменяться, после этого определяют и фиксируют объем добываемой жидкости Si=nik3i, где n - число качаний насоса, k3 - коэффициент заполнения, после этого снова увеличивают частоту качаний насоса n и не изменяют ее до тех пор, пока коэффициент заполнения не перестанет изменяться, определяют и фиксируют объем добываемой жидкости, при новой частоте качаний, далее полученное значение Si=1 сравнивают с предыдущим Si и, если последующее значение окажется больше предыдущего, то следующий шаг изменения частоты осуществляют в ту же сторону, а если наоборот - в обратную сторону. 1 з.п. ф-лы, 1 ил.

Изобретение относится к горному делу, в частности к нефтегазовой промышленности, и может использоваться для замера профиля насосно-компрессорных и обсадных труб нефтегазовых скважин. Техническим результатом является сокращение трудоемкости замера профиля труб нефтегазовых скважин. Устройство содержит немагнитный корпус, состоящий из двух идентичных частей, выполненных из немагнитного материала и соединенных между собой. На каждой части корпуса установлены диаметрально противоположные кинематические подвижные пары шарнирно соединенных измерительных и вспомогательных рычагов, оснащенные опорными роликами для контакта с внутренней стенкой трубы и расположены в перпендикулярных друг другу плоскостях, пересекающихся вдоль центральной оси корпуса. Измерительные рычаги функционально соединены со своими преобразователями линейного перемещения в электрический сигнал, выполненными в виде цилиндрических катушек индуктивности, размещенных внутри каждой части корпуса, и ферромагнитных втулок, подвижно установленных снаружи частей корпуса в зоне действия магнитного поля цилиндрических катушек индуктивности. Ферромагнитные втулки шарнирно соединены с концами измерительных рычагов и подпружинены в направлении кинематических подвижных пар до упора в буртики, выполненные на каждой части корпуса, ограничивающие линейные перемещения ферромагнитных втулок на длине обмотки цилиндрических катушек индуктивности. Цилиндрические катушки индуктивности электрически связаны с блоком электроники, размещенным внутри герметичного корпуса, измеряющим суммарную напряженность магнитного поля цилиндрических катушек индуктивности и намагниченности ферромагнитных втулок с возможностью преобразования возбуждаемой эдс в электрические сигналы и передачи их на поверхность скважины по геофизическому кабелю, герметично закрепленному в корпусе. Вспомогательные рычаги шарнирно соединены с кольцами, подвижно установленными снаружи частей корпуса и подпружиненными в направлении кинематических подвижных пар до упора во вторые буртики, выполненные на каждой части корпуса, ограничивающие возвратно-поступательные перемещения колец при взаимодействии вспомогательных рычагов и вторых пружин. Последние выполнены с допустимым начальным усилием сжатия, превышающим допустимое конечное усилие сжатия первых пружин. Головная часть корпуса снабжена головкой, снабженной буртиком для упора второй пружины. 1 з.п. ф-лы, 1 ил.

Изобретение относится к способам исследования газовых и газоконденсатных скважин, определению их оптимальных технологических режимов, а именно к определению режимов максимального извлечения жидких продуктов при минимальных энергетических затратах, то есть минимальных потерях давления при различных режимах течениях газожидкостного потока. Технический результат предлагаемого технического решения заключается в расширении области исследования скважин при более широких режимах сепарации газожидкостных смесей. Способ включает сепарацию продукции скважины, замер дебита газа и объема, выносимых твердых и жидких фаз, при замере жидкой фазы под давлением сепарации и после дегазации в отдельной емкости. Замер устьевых давлений и температур на нескольких установившихся режимах, проводимый до стабилизации замеряемых параметров при утилизации отделяемых фаз. Замер объема жидкости и механических примесей осуществляется поочередно: на режиме сепарации без сброса давления газожидкостного потока перед сепарацией, на режиме сепарации со сбросом давления на сужающих устройствах (дросселях) перед сепарацией, на режиме сепарации со сбросом давления на сверхзвуковом сепараторе. Максимальный объем отсепарированной жидкости определятся по замеренным величинам при минимальных гидравлических потерях газожидкостной смеси. Отделенные фазы смешиваются с отсепарированным газовым потоком на всех режимах сепарации. Сепарация газожидкостной смеси осуществляется путем равномерного распределения сепарируемых фаз на вертикальной пористой структуре и последующим отводе с нее накопленной жидкости. Потоки газа и жидкости направляются на сепарацию продукции после режима сверхзвуковой сепарации. 5 з.п. ф-лы, 1 ил.

Изобретение относится к нефтяной и газовой промышленности, в частности к ликвидации скважин, выполнивших свое назначение. Способ ликвидации скважины с заколонными перетоками включает спуск колонны труб и установку цементных мостов в обсадной колонне скважины. Геофизическими исследованиями определяют наличие цементного кольца за эксплуатационной колонной и выявляют нарушения герметичности эксплуатационной колонны. Затем при отсутствии цементного кольца за эксплуатационной колонной и выявлении нарушений герметичности эксплуатационной колонны тампонированием производят наращивание цементного кольца за эксплуатационной колонной и герметизацию зон нарушения эксплуатационной колонны. Определяют качество тампонирования эксплуатационной колонны проведением повторных геофизических исследований. Затем в скважину от устья до забоя спускают заглушенную снизу колонну труб малого диаметра. Далее до забоя скважины спускают дополнительную колонну труб и производят установку цементного моста тампонированием под давлением от забоя до устья скважины с использованием термостойкого цемента с добавлением фиброволокна в количестве 0,2% от массы сухого цемента. Затем извлекают дополнительную колонну труб из скважины, доливают в ствол скважины термостойкий цемент до устья. После чего заполняют колонну труб малого диаметра незамерзающей жидкостью. Далее в колонну труб малого диаметра до забоя спускают оптоволоконный кабель. После ликвидации скважины периодически фиксируют температурное распределение в стволе скважины. Техническим результатом является повышение эффективности и надежности ликвидации скважины. 4 ил.

Группа изобретений относится к технике для исследования движения жидкостных потоков и сыпучих материалов, газожидкостных потоков, например процессов добычи газа в нефтегазовой отрасли, связанной с изучением процессов движения газожидкостных потоков в вертикальных, наклонных трубопроводах и отдельных устройствах. Техническим результатом является повышение точности проведения экспериментов и расширение функциональных возможностей за счет обеспечения возможности визуализации с улучшенным качеством происходящих в объеме и по длине лифтовой колонны труб и трубопровода процессов, а также за счет определения газосодержания потока как вдоль выбранного произвольно сечения, так и на произвольно выбранном ограниченном участке лифтовой колонны труб. Устройство содержит лифтовую колонну труб, блоки ввода в колонну и отвода из колонны, фото-видеорегистратор, импульсные вспышки с установленными на них софтбоксами, радиосинхронизатор фотовспышек, состоящий из источника сигнала и приемника сигнала, экран, сливной патрубок и разветвитель для синхрокорда. 2 н.п. ф-лы, 5 ил.

Изобретение относится к способу, устройству и машиночитаемому носителю данных, используемых при построении геологической модели нефтяного или иного месторождения. Технический результат - повышение точности вычислений параметров, используемых при построении геологической модели расположения нефтяных или иных месторождений. Изобретение позволяет для маркеров, выбранных в качестве начального решения, вычислить такие глубины маркера на каждой скважине, которые обеспечивают наилучшую суммарную корреляцию. Для каждого маркера, входящего в набор, определяется функционал, представляющий собой сумму коэффициентов корреляции комплекса методов ГИС для пар скважин, расположенных не далее заданного расстояния друг от друга. Для этого функционала вычисляются частные производные ,и полученный таким образом вектор сглаживается и используется для нахождения большего значения функционала на некотором отрезке вдоль этого вектора. Если большего значения не найдено, то последнее положение отметок маркера считается решением задачи, а если найдено, то производится сглаживание точки решения и процесс повторяется снова. На каждой итерации алгоритма производится сортировка глубин маркеров. 3 н.п. ф-лы, 9 ил.
Данное изобретение относится к способами оценки продуктивных пластов на нефтегазовых месторождениях, в частности к оценке их свойств. Технический результат заключается в более эффективной оценке свойств пористого пласта. Способ оценки свойств продуктивного пласта, пробуренного скважиной, включает закачку флюида с множеством индикаторных добавок субмикронного размера в ствол скважины и продуктивный пласт, ожидание обратного притока и определение свойств пласта. Данные свойства определяются посредством анализа изменений функции распределения индикаторов по размерам и типу в закачанном и добытом флюидах. 17 з.п. ф-лы.

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