Способ определения характеристик скважины, призабоспособ определения характеристик скважины, призабойной зоны и пласта и устройство для его осуществлейной зоны и пласта и устройство для его осуществления ния - RU2179637C1

Чертежи

Описание

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

Известен способ разработки продуктивного пласта, включающий импульсную нестационарную закачку пластовой жидкости, замеры на устье давления закачки и расхода жидкости, определение накопленного расхода и производной функции репрессии, характеризующей работу, затрачиваемую на нестационарное течение в пласте единицы расхода пластовой жидкости, построение графика зависимости производной функции репрессии от накопленного расхода для диапазона значений гидропроводности, заведомо включающих искомую гидропроводность и выбор среди множества полученных кривых производной линии, ближе остальных отвечающей условию постоянства производной, соответствующей искомой гидропроводности пласта (Патент РФ 2151859, кл. Е 21 В 43/20, опублик. 2000 г.).

Известен способ эксплуатации скважины с одновременным определением параметров загрязненной призабойной зоны скважины, включающий импульсную нестационарную закачку пластовой жидкости со ступенчатыми изменениями расхода от минимальных до максимальных величин с задаваемым периодом, через каждые 5÷60 с замер и регистрацию давления, плотности и расхода закачиваемой пластовой жидкости, пересчет данных на забойные условия, для каждого замера в условиях импульсной нестационарной закачки пластовой жидкости во время каждого режима закачки определение функции репрессии, характеризующей нестационарное течение в призабойной зоне скважины на протяжении данного режима закачки жидкости, для каждого режима построение графика зависимости функции репрессии от логарифма времени закачки на данном режиме, на каждом из полученных графиков выделение начального наклонного прямолинейного участка, методом наименьших квадратов нахождение параметров выделенных прямолинейных участков, по которым определение гидропроводности и пьезопроводности загрязненной призабойной зоны, а также ее радиуса и коэффициента скин-эффекта (Патент РФ 2151858, кл. Е 21 В 43/20, опублик. 2000 г.).

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

Наиболее близким к изобретению по технической сущности является способ эксплуатации скважины, в результате осуществления которого происходит определение характеристик скважины, призабойной зоны и пласта, включающий импульсную нестационарную закачку реагента, на устье скважины замер давления закачки и расхода реагента, пересчет данных замеров на забойные условия, определение накопленного расхода и работы, затрачиваемой на нестационарное течение в призабойной зоне скважины единицы расхода реагента, по этим показателям с учетом текущей проводимости пласта рассчет коэффициента скин-эффекта, в качестве текущей проводимости использование величины, определенной по результатам кратковременного импульсного нестационарного испытания данной скважины на приемистость пластовой жидкостью, изменение режима закачки реагента при достижении необходимых фильтрационных свойств призабойной зоны скважины, определенных по скин-эффекту, рассчитанному по накопленному расходу и работе течения единицы расхода реагента в призабойной зоне скважины с учетом текущей проводимости пласта (Патент РФ 2151855, кл. Е 21 В 43/20, опублик. 2000 г. - прототип).

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

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

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

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

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

Известно контролирующее устройство для газовой скважины. Устройство, устанавливаемое на устье газовой скважины для определения величины давления на забое скважины, содержит фиксированный измерительный комплекс. Последний имеет датчики для измерения давления и температуры газа, которые введены в газовый поток, идущий через устье скважины. Для систематического контроля за правильностью измерений, проведенных с помощью датчиков, измерительный блок содержит автоматическое устройство, обеспечивающее периодическое взятие проб из газового потока, идущего через устье скважины. С этим устройством технологически связан процессор, обеспечивающий вычисление величины давления газа на дне скважины, которые он производит на основании данных, полученных с помощью датчиков, установленных в газовом потоке, идущем через устье скважины. К процессору подсоединен блок памяти, обеспечивающий прием, запоминание и хранение данных о давлении и температуре газа, которые поступают в блок памяти процессора через заданные интервалы времени. К блоку памяти подключен дисплей, на экране которого демонстрируется цифровая информация о давлении и температуре в газовом потоке, идущем через устье скважины, а также информация о давлении газа на забое скважины (Патент США 4414846, кл. 37-151, опублик. 1983).

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

Наиболее близким к изобретению по технической сущности является для определения характеристик скважины, призабойной зоны и пласта, включающее датчик давления и датчики расхода, подключенные к приспособлению для регистрации параметров среды. (А.С. СССР 1469110, кл. Е 21 В 47/00, опублик. 1989 г.)
Известное устройство позволяет замерять расход рабочего агента при закачке в скважину и определять направление его движения в скважине. Однако устройство не позволяет контролировать такие параметры, как давление и его изменение. Кроме того, устройство способно определять только параметры непосредственно в точке определения и не способно определять параметры на удалении, например, на забое скважины.

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

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

Кроме того, на измерительном участке дополнительно могут быть размещены датчики для определения плотности и температуры.

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

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

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

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

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

Устройство включает измерительный участок 1 с датчиками расхода 2, датчиком давления 3, дифференциальным манометром 4 с импульсными трубками 5, 6, соединенными соответственно с началом и концом участка. Устройство соединено через нагнетательную линию 7 со скважиной 8. Измерительный участок 1 выполнен длиной, позволяющей фиксировать перепады давлений при течении жидких сред с минимальными гидравлическими сопротивлениями.

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

Выходы датчика давления 3 и дифференциального манометра 4 соединены электрическими кабелями 9 и 10 с блоками искрозащиты 11 и 12, размещенными в выносном блоке 13, и далее с блоком сбора информации 14. Выходы датчиков расхода 2 соединены электрическими кабелями 15 и 16 с вторичными блоками расходомеров 17 и 18 и далее с блоком сбора информации 14. Блок сбора информации 14 соединен с компьютером 19.

Устройство работает следующим образом.

При нагнетании рабочего агента через измерительный участок 1 в скважину 8 аналоговые сигналы с датчика давления 3 и дифференциального манометра 4 посредством электрических кабелей 9 и 10 через блоки искрозащиты 11 и 12 поступают в выносной блок 13 и далее в блок сбора информации 14. В блоках искрозащиты 11 и 12 производится гальваническая развязка электрических цепей.

Частотные сигналы с датчиков расхода 2 посредством электрических кабелей 15 попадают на вторичные блоки расходомеров 17 и 18, откуда с помощью соединительных кабелей 16 поступают на вход соответствующих каналов блока сбора информации 14.

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

При проведении технологической операции по воздействию на нефтяной пласт для интенсификации отбора или изоляции водопритока, выравнивания профиля притока или поглощения расход нагнетаемой рабочей жидкости остается относительно неизменным лишь в течение отдельных весьма коротких промежутков времени и изменяется в широких пределах в течение всей операции. В предложенном способе изначально закладывают режим импульсной нестационарной закачки реагента как наиболее общий и в наибольшей мере отвечающий условиям производства. Реализуемый в особых условиях на практике стационарный режим закачки реагента является частным случаем общего импульсного нестационарного режима. При этом справедливы все расчеты и выводы предлагаемого способа. Процесс импульсной нестационарной закачки реагента характеризуется значительными колебаниями расхода и давления с хаотичными изменениями величин по амплитуде и частоте. Расход может изменяться по амплитуде от 0,084 до 7,6 л/с, по частоте - от 0,002 до 0,02 Гц, при этом максимальный расход обеспечивает недопущение развития искусственной трещиноватости в призабойной зоне (максимально допустимое забойное давление в процессе закачки жидкости должно быть меньше давления раскрытия трещин в призабойной зоне скважины). Устьевое давление нагнетания может изменяться по амплитуде от 1 до 10÷15 МПа при той же частоте.

При проведении технологических операций на скважине посредством информационно-измерительного комплекса измеряют и регистрируют устьевое давление, плотность, потери давления на измерительном участке и объемный расход нагнетаемого реагента с интервалом в 5÷60 с (т.е. с периодом опроса 5÷60 с). Для каждого замера по этим данным в реальном времени процесса рассчитывают забойное давление с учетом искривления ствола скважины, реологии жидкости и ее нагрева, изменениями в этой связи гидростатического давления и потерь на жидкостное трение в насосно-компрессорных трубах при закачке в насосно-компрессорные трубы обычных ньютоновских жидкостей, а также полимерных, глинистых, цементных растворов и других неньютоновских жидкостей. Забойное давление определяют в такой последовательности.

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

На α^-M этапе будет закачана α^-Я жидкость (при α^=1;2 и т.д. в зависимости от числа жидкостей для закачки). Для каждого замера расхода нагнетаемой α^-Й жидкости Q_α(t) и перепада давления на измерительном участке ΔP_ИЗМ(t) в реальном времени процесса рассчитывают вспомогательные параметры G_α, U_α:

где Q_α(t)- расход α^-Й нагнетаемой жидкости в момент времени t после начала закачки, м³ /сут;
d_ИЗМ - внутренний диаметр измерительного участка, м;
L_ИЗМ - длина измерительного участка (расстояние между осями импульсных трубок на измерительном участке), м;
ΔP_ИЗМ (t)- перепад давлений на измерительном участке (между осями импульсных трубок на измерительном участке) в момент времени t после начала закачки α^-Й жидкости, МПа;
α^- порядковый номер нагнетаемой жидкости.

Размерности вспомогательных параметров G_α, U_α таковы:
|G_α | = 1/сут; |U_α| = МПa.
Величины вспомогательных параметров G_α и U_α, рассчитанные по формулам (1) для текущего момента t, наносят на график зависимости
U_α= U_α(G_α). (2)
На фиг.3 представлен график зависимости вспомогательного параметра U_α от вспомогательного параметра G_α на примере закачке в скважину гелеобразующего состава: водного раствора сополимера "Комета" и смолы "ДЭГ", где по оси абсцисс отложены величины lgG, по оси ординат - величины lgU.

После получения первых 30÷40 значений U_α, G_α производят аппроксимацию полученного массива точек подбором функциональной зависимости U_α= U_α(G_α). В дальнейшем по мере поступления новых данных (значений U_α, G_α)) зависимость U_α= U_α(G_α) уточняют.

После установления функциональной зависимости (2) для каждого замера Q_α(t) расхода α^-Й жидкости в реальном времени процесса вычисляют вспомогательный параметр

где d_НКТ - внутренний диаметр насосно-компрессорных труб, м.

Размерность вспомогательного параметра:

Принимая

из функциональной зависимости U_α= U_α(G_α) определяют величину соответствующую величине

Размерность вспомогательного параметра:
Для каждого замера расхода Q_α(t) в реальном времени процесса рассчитывают коэффициент λ[(fluid_α),Δt] гидравлических сопротивлений потоку в насосно-компрессорных трубах α^-Й нагнетаемой жидкости:

где вспомогательный параметр, определяемый выражением (5), МПа;
d_НКТ - внутренний диаметр насосно-компрессорных труб, м;
ρ_УСT(fluid_α)- плотность α^-Й нагнетаемой жидкости в устьевых условиях, кг/м³;
Q_α(t)- расход α^-Й жидкости в момент времени t закачки, м³/сут;
λ [(fluid_α),Δ t]- коэффициент гидравлических сопротивлений потоку в насосно-компрессорных трубах α^-Й нагнетаемой жидкости, величина безразмерная.

Определенные из уравнения (6) значения коэффициента λ[(fluid_α),Δt] гидравлических сопротивлений потоку в насосно-компрессорных трубах α^-Й жидкости наносят на график зависимости λ[(fluid_α),Δt] от расхода α^-Й жидкости:
λ[(fluid_α),Δt] = Ф(Q_α (t)). (7)
На фиг.4 представлен график зависимости коэффициента гидравлических сопротивлений в насосно-компрессорных трубах λ[(fluid_α),Δt] от расхода Q_α(t) гелеобразующего состава: водного раствора сополимера "Комета" и смолы "ДЭГ", где по оси абсцисс отложены величины lgQ_α(t), по оси ординат - величины lgλ [(fluid_α),Δt].
После получения первых 30÷40 точек [значений λ и Q(t)] производят аппроксимацию полученного массива подбором корреляционной зависимости λ[(fluid_α),Δt] = Ф(Q_α(t)). В дальнейшем по мере поступления новых данных [значений λ и Q(t)] зависимость (7) уточняют.

Для каждого замера расхода Q_α(t) по корреляционной зависимости λ[(fluid_α),Δt] = Ф(Q_α(t)) в реальном времени процесса рассчитывают коэффициент λ [(fluid_α),Δt] гидравлических сопротивлений потоку в насосно-компрессорных трубах α^-Й нагнетаемой жидкости.

По этим данным для каждого замера расхода Q_α(t) в реальном времени процесса рассчитывают потери давления Р_ТР(t) на жидкостное трение в насосно-компрессорных трубах потока α^-Й нагнетаемой жидкости:

где Р_ТР(t) - потери давления потока α^-Й нагнетаемой жидкости в насосно-
компрессорных трубах вследствие жидкостного трения в момент времени t, МПа;
L - длина насосно-компрессорных труб (расстояние по стволу скважины от устья до башмака насосно-компрессорных труб), м;
λ[(fluid_α),Δt]- коэффициент гидравлических сопротивлений потоку в насосно-компрессорных трубах α^-Й нагнетаемой жидкости, определенный для каждого замера расхода Q_α(t) по корреляционной зависимости (7)
λ [(fluid_α),Δt] = Ф(Q_α (t));
d_НКТ - внутренний диаметр насосно-компрессорных труб, м;
ρ_УСT(fluid_α)- плотность α^-Й жидкости в устьевых условиях, кг/м³;
Q_α(t)- расход α^-Й жидкости в момент времени t закачки, м³/сут.

Динамическое забойное давление Р_С(t) на глубине башмака насосно-компрессорных труб L в текущий момент времени t равно:
Р_C(t)=P_УСТ(t)+Р_Г(t)-P_ТР(t) (9)
где Р_С(t) - забойное давление на глубине L башмака насосно-компрессорных труб в момент t нагнетания α^-Й жидкости, МПа;
Р_УСТ(t) - устьевое давление в момент t закачки α^-Й жидкости, МПа;
Р_Г(t) - гидростатическое давление, создаваемое столбом α^-Й жидкости в момент времени t, МПа;
Р_ТР (t) - потери давления потока α^-Й нагнетаемой жидкости в насосно-компрессорных трубах вследствие жидкостного трения в момент времени t, МПа, определены по формуле (8).

Репрессия на пласт ΔP_C(t) в момент t закачки α^-Й жидкости равна:
ΔP_C(t) = P_C(t)-P_ПЛ, (10)
где Р_ПЛ - пластовое давление, приведенное к глубине L башмака насосно-компрессорных труб, МПа.

Для определения коэффициента S скин-эффекта при проведении технологических операций на скважине замеряют и регистрируют устьевое давление, плотность и объемный расход α^-Й нагнетаемой жидкости с интервалом в 5÷60 с (т.е. с периодом опроса 5÷60 с). Для каждого замера в момент t по этим данным в реальном времени процесса рассчитывают потери давления Р_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление Р_Г(t) столба α^-Й нагнетаемой жидкости, динамическое забойное давление Р_С(t) по формуле (9), репрессию на пласт ΔP_C (t) по формуле (10), объемный расход α^-Й жидкости Q(t) в забойных условиях. Затем для каждого текущего момента времени t_N определяют величину функции репрессии Y(t_N), характеризующей работу единицы расхода α^-Й жидкости на нестационарное течение в призабойной зоне скважины, по формуле:

где N=2; 3; 4;... - номер текущего замера устьевого давления, плотности и объемного расхода α^-Й нагнетаемой жидкости;
n=0; 1; 2; 3;... N-1 - номера предыдущих замеров;
t₀ - время начала закачки (начальный замер n=0), с;
t₁;... t_n - время первого,... n замеров, с;
t_N - время текущего замера, с;
ΔP_C(t_O); ...ΔP_C (t_n)- репрессия на пласт в начале закачки и в момент времени t_n предшествующего n замера, Па;
ΔP_C(t_N)- репрессия на пласт в момент t_N текущего N замера, Па;
Q₀, . . . Q_n - объемный расход α^-Й жидкости в забойных условиях в начале закачки и в момент предшествующего n замера, м³/с;
Y(t_N) - функция репрессии, характеризующая работу, затраченную на нестационарное течение в призабойной зоне единицы расхода α^-Й жидкости, в текущий момент времени t_N с начала процесса, Па•с;
ε^- гидропроводность пласта, м²•м/Па•c:

k - проницаемость пласта для пластовой жидкости, м²;
h - эффективная толщина продуктивного пласта, принимающего нагнетаемую жидкость, м;
μ- вязкость пластовой жидкости, Па•с.

Одновременно с функцией репрессии Y(t_N) вычисляют накопленный объем жидкости в забойных условиях W(t_N), поступившего в пласт к моменту времени t_N с начала закачки, по формуле:

Полученные значения Y(t_N) и W(t_N) наносят на график.

На фиг. 5 представлен график зависимости функции репрессии Y(t_N) от накопленного объема жидкости W(t_N) на примере закачки в скважину гелеобразующего состава: водного раствора сополимера "Комета" и смолы "ДЭГ", где по оси абсцисс отложены величины накопленного объема W(t_N), м³, по оси ординат - величины функции репрессии Y(t_N), МПа•с. На фиг.5 приняты следующие условные обозначения: 1 - первый прямолинейный участок при закачке в пласт 6,7 м³ гелеобразующего состава, 2 - второй прямолинейный участок при дальнейшей закачке в пласт 1,2 м³ гелеобразующего состава, 3 - третий прямолинейный участок при дальнейшей закачке в пласт 4,4 м³ гелеобразующего состава, 4 - четвертый прямолинейный участок при продавке в пласт 4,2 м³ пластовой воды.

При наличии цифровой регистрации устьевых параметров и системы компьютерного анализа определение величин Y(t_N), W(t_N) и построение графика зависимости:
Y(t_N)=Y[W(t_N)] (14)
производят непосредственно в процессе технологического воздействия в реальном времени.

Производят аппроксимацию отдельных участков графика зависимости (14) прямолинейными отрезками. В интервале времени [t_j, t_j+1] линейной аппроксимации определяют наклон прямолинейного участка B_j. Величину коэффициента скин-эффекта S_j, отражающего состояние призабойной зоны скважины в интервале времени [t_j, t_j+1] технологической операции определяют по формуле:

где S_j - коэффициент скин-эффекта, отражающий дополнительные фильтрационные сопротивления потоку жидкости вследствие загрязнения и несовершенства вскрытия призабойной зоны скважины, в интервале времени [t_j, t_j+1] технологического воздействия;
r_C - радиус скважины, м;
χ- пьезопроводность продуктивного пласта, м²/с;
B_j - наклон графика зависимости (14) в интервале времени [t_j, t_j+1] технологического воздействия, Па•с/м³.

После достижения запланированной величины скин-эффекта изменяют режим закачки вплоть до ее прекращения.

При определении гидропроводности пласта ε проводят закачку пластовой жидкости в добывающую или нагнетательную скважину. До проведения операции на скважине задают произвольный ряд M значений гидропроводности пласта ε_m
ε₁ < ε₂ <...< ε_m< ...< ε_M, (16)
заведомо включающих истинную величину гидропроводности пласта ε_ИСТ
ε₁< ε_ИСТ< ε_M. (17)
Организуют процесс импульсной нестационарной закачки пластовой жидкости в скважину. При этом на устье выполняют замеры и регистрацию устьевых давлений, плотности и объемного расхода нагнетаемой пластовой жидкости. Для каждого замера в момент времени t по этим данным в реальном времени процесса рассчитывают потери давления P_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление P_Г(t) столба нагнетаемой жидкости, динамическое забойное давление Р_С(t) по формуле (9), репрессию на пласт ΔP_C(t) по формуле (10).

Рассчитывают объемный расход нагнетаемой жидкости Q(t) в забойных условиях. Затем для каждого из принятых значений ε_m гидропроводности пласта определяют значения производной ΔY_m/ΔX_m(t_N) функции репрессии на пласт по формуле

где N, N-1 - номера текущего и предыдущего замеров (N=2; 3; 4;...) устьевого давления, плотности и объемного расхода нагнетаемой жидкости;
i=0; 1; 2;... N-2 - номера предшествующих замеров;
t_N, t_N-1 - время текущего и предыдущего замеров, с;
t₀; t₁;... t_i - время предшествующих замеров, с;
∋P_C(t_N), ∋ P_C(t_N1 )- репрессия на пласт в текущем и предыдущем замерах, Па;
Q_N , Q_N-1 - объемный расход жидкости в забойных условиях в текущем и предыдущем замерах, м³/с;
Q_i - объемный расход жидкости в забойных условиях в предшествующих замерах, м³/с;
8_m(t_N) - функция репрессии, характеризующая работу, затраченную на нестационарное течение единицы расхода жидкости в пласте гидропроводностью H_m, в момент времени t_N с начала процесса, Па•с;
H_m - принятая в расчетах гидропроводность пласта, м²•м/Па• с:

k_m - принятая в расчетах проницаемость пласта для пластовой жидкости, м²;
h - эффективная толщина продуктивного пласта, принимающего нагнетаемую пластовую жидкость, м;
П - вязкость пластовой жидкости, Па•с.

Одновременно с ∋8_m/∋&_m(t_N) для того же момента времени t_N по формуле (13) определяют накопленный объем жидкости в забойных условиях W(t_N), поступившей в пласт к моменту времени t_N с начала закачки.

Полученные значения наносят на график.

На фиг. 6 представлен график зависимости производной ∋8_m/∋&_mt_N функции репрессии на пласт от накопленного объема жидкости в забойных условиях W(t_N) на примере закачки в скважину пластовой жидкости, где по оси абсцисс отложены величины накопленного объема жидкости, м³, по оси ординат - величины производной функции репрессии на пласт. На фиг.6 приняты следующие условные обозначения:

график производной, когда в расчетах принята гидропроводность пласта 5,1 мкм²•м/мПа•с;
-■- график производной, когда в расчетах принята гидропроводность пласта 20,4 мкм²•м/мПа•с;
-Δ- график производной, когда в расчетах принята гидропроводность пласта 10,3 мкм²• м/мПа•с.

Графики производной ∋8/∋& существенно зависят от принятой величины гидропроводности пласта H_m. Чем ближе значения H_m к истинной величине гидропроводности пласта H_ИСТ, тем ближе графики производной ∋8/∋ & к прямой, параллельной оси абсцисс. Если истинное значение H_ИСТ включено в диапазон (17), то среди полученных кривых зависимости

устанавливаются одна - две линии, ближе остальных отвечающих условию:
∋8/∋& >t, H_ИСТ≅ = const. (21)
Далее известным методом последовательного приближения находят величину гидропроводности пласта H, при которой производная ∋8/∋& может быть принята постоянной наилучшим образом. Выбор оптимального выполнения условия (21) достигают численными методами с использованием аппарата практической физики. Величина, обеспечивающая выполнение условия наилучшим образом, и является искомым значением гидропроводности пласта H.

До определения параметров призабойной зоны предлагаемым способом на скважине организуют предварительные исследования с целью уточнения гидропроводности пласта ε и установления факта существенного загрязнения призабойной зоны (S> 20÷30).

Если полученная тем или иным известным способом величина коэффициента скин-эффекта S≥20÷30, реализуют предлагаемый способ.

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

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

где P_Cмакс - максимально допустимое забойное давление в процессе закачки пластовой жидкости, МПа;
σ_раск- давление раскрытия трещин в призабойной зоне пласта, МПа.

Установлено, что для получения достоверных результатов необходимо провести закачку на нескольких (4÷6 и более) режимах нагнетания с резкой сменой расхода с большего на меньший и наоборот.

Продолжительность закачки Δθ на каждом режиме устанавливают опытным путем или приближенно оценивают так:

где S - значение коэффициента скин-эффекта, определенное при предварительных испытаниях скважины;
χ- пьезопроводность пласта, м²/с.

В формуле (23) коэффициент S безразмерен, а размерность продолжительности закачки на каждом режиме такова: |Δ θ| = c.
На основе выполненных оценок на устье организуют основной процесс импульсной нестационарной закачки пластовой жидкости так, чтобы кривая переменного расхода являлась ступенчатой функцией времени t нагнетания:

где t - текущее время от начала процесса основной закачки пластовой жидкости, с;
Z=1, 2,... - порядковый номер режима основной закачки;
время начала и окончания Z режима основной закачки (начало первого режима закачки ), с;
средний расход в течение Z режима основной закачки (приближенное равенство означает, что в процессе основной закачки пластовой жидкости допускаются колебания расхода до 20÷30% от средней величины расхода на данном режиме), м³/с.

В процессе закачки на скважине замеряют и регистрируют устьевое давление, плотность и объемный расход пластовой жидкости с интервалом в 5÷60 с (т. е. с периодом опроса 5÷60 с). Для каждого замера в момент времени t по этим данным в реальном времени процесса рассчитывают потери давления P_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление Р_C(t) столба нагнетаемой пластовой жидкости, динамическое забойное давление Р_C(t) по
формуле (9), репрессию на пласт ΔP_C(t) по формуле (10), объемный расход пластовой жидкости Q(t) в забойных условиях. Затем для каждого N замера, выполненного в текущем временном интервале

на Z режиме закачки, определяют функцию репрессии Ψ_Z(Δt_Z), характеризующую нестационарное течение нагнетаемой жидкости в призабойной зоне на данном режиме, по формуле:
(25)
где N=2; 3; 4... - номер текущего замера;
i=0; 1; 2;... N-1 - номера предшествующих замеров;
Δ P_C(t_N), ΔP_C(t_N-1)- репрессия на пласт в текущем и предыдущем замерах, Па;
Q_N, Q_N-1 - объемный расход жидкости в забойных условиях в текущем и предыдущем замерах, м³/с;
Q_i - объемный расход жидкости в забойных условиях в предшествующих замерах, м³/с;
t_N - время текущего замера, с;
t₀, t₁,... t_N-1 - время предшествующих замеров, с;
∋t_Z- временной интервал текущего режима основной закачки при

Ф - пьезопроводность пласта, м²/с;
r_c - радиус скважины, м.

Величина функции репрессии <_Z (∋t_Z) безразмерна.

Расчеты по формуле (25) выполняют последовательно для всех замеров устьевых параметров. Для каждого режима Z закачки на основе выполненных замеров устьевых параметров и расчетов по приведенным формулам строят график.

На фиг. 7 представлен график зависимости функции репрессии <_Z(∋t_Z) от временного интервала ∋ t_Z на примере закачки в скважину пластовой жидкости, где по оси абсцисс отложены величины ln∋t_Z, по оси ординат - величины функции репрессии <_Z(∋ t_Z), соответствующие данному временному интервалу ∋t_Z.
На фиг.7 приняты следующие условные обозначения: Z=1, 2... 10 - графики функции репрессии в период времени ∋t_Z закачки пластовой жидкости в нефтедобывающую скважину на Z режиме с расходом

При наличии цифровой регистрации устьевых параметров и системы компьютерного анализа определение величин ln∋t_Z, <_Z(∋ t_Z) и построение графиков зависимости Ψ_Z(Δt_Z) = Ф(lnΔt_Z) производят непосредственно в процессе закачки жидкости в реальном времени t_N текущего замера.

Таким образом, каждому из режимов основной закачки соответствует своя линия (фиг. 7). В общем случае на каждом из полученных графиков выделяют начальный наклонный прямолинейный участок (пример дан в табл. 1), отражающий нестационарное течение нагнетаемой жидкости в загрязненной призабойной зоне на данном Z режиме закачки и описываемый уравнением прямой линии:
Ψ_Z(Δt_Z) = a_Z+b_Z•lnΔt_Z. (27)
Известным математическим методом наименьших квадратов находят наклон b_Z и начальный участок а_Z выделенного прямолинейного отрезка на каждом Z режиме закачки. После этого определяют:
гидропроводность призабойной зоны скважины ε_ПЗС

пьезопроводность призабойной зоны скважины χ_ПЗС

Поскольку все прямые уравнения (9) имеют общую точку пересечения, используя найденные коэффициенты a_Z, b_Z, a_Z-1, b_Z-1 для двух смежных режимов (Z, Z-1) закачки, определяют коэффициент скин-эффекта S:

после чего рассчитывают радиус R_ПЗС зоны загрязнения:

В формулах (28)-(31) размерности величин таковы: [ε] = м² •м/Па•c; [χ] = м²/c; [R]=м, коэффициенты a_Z, b_Z, S безразмерны.

Примеры конкретного выполнения
Пример 1. На нагнетательной линии перед устьем скважины организуют измерительный участок согласно фиг.1 и 2 длиной 6 м, позволяющей фиксировать перепады давлений при течении жидких сред с минимальными гидравлическими сопротивлениями.

Измерительный участок выполняют в виде калиброванной трубы 1 с внутренним диаметром 62 мм, на которой монтируют два датчика расхода 2 марки ПЭА1, датчик давления 3 марки "МИДА" и дифференциальный манометр 4 типа "Сапфир" с импульсными трубками 5 и 6, соединенными с началом и концом измерительного участка. Замеры давления, расхода и перепада давлений при закачке рабочего агента проводят на измерительном участке 1. Измерительный участок 1 выполнен длиной, позволяющей фиксировать перепады давлений при течении жидких сред с минимальными гидравлическими сопротивлениями. Устройство соединено через нагнетательную линию 7 со скважиной 8.

Выходы датчика давления 3 марки "МИДА" и дифференциального манометра 4 типа "Сапфир" соединены электрическими кабелями 9 и 10 с блоками искрозащиты 11 типа "Корунд" и 12 типа "Взлет", размещенными в выносном блоке 13, и далее с блоком сбора информации 14. Выходы датчиков расхода 2 марки ПЭА1 соединены электрическими кабелями 15 и 16 с вторичными блоками расходомеров 17 типа "Взлет" БИИ и 18 "Днепр-7" БП и далее с блоком сбора информации 14. Блок сбора информации 14 соединен с компьютером 19 типа Ноутбук.

При нагнетании рабочего агента через измерительный участок 1 в скважину 8 аналоговые сигналы с датчика давления 3 и дифференциального манометра 4 посредством электрических кабелей 9 и 10 через блоки искрозащиты 11 и 12 поступают в выносной блок 13 и далее в блок сбора информации 14. В блоках искрозащиты 11 и 12 производится гальваническая развязка электрических цепей.

Частотные сигналы с датчиков расхода 2 посредством электрических кабелей 15 попадают на вторичные блоки расходомеров 17 и 18, откуда с помощью соединительных кабелей 16 поступают на вход соответствующих каналов блока сбора информации 14.

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

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

На устье скважины организуют процесс импульсной нестационарной закачки реагента с колебаниями расхода и давления с хаотичными изменениями величин по амплитуде и частоте. Расход изменяется по амплитуде от 0,084 до 7,6 л/с, по частоте - от 0,002 до 0,02 Гц. Давление изменяется по амплитуде от 1 до 10÷15 МПа при той же частоте.

Процесс обработки скважины заключается в закачке в призабойную зону пласта порций гелеобразующего состава (α = 1) и его продавке пластовой водой (α = 2). В качестве гелеобразующего состава используют водный раствор сополимера "Комета" и смолы "ДЭГ", образующих псевдопластичную систему. Закачку ведут с начальным расходом 5,3 л/с.

При закачке гелеобразующего состава на скважине выполняют замеры и регистрацию устьевого давления, перепада давлений на измерительном участке, плотности и объемного расхода нагнетаемого раствора с периодом опроса 5 с. Для каждого замера расхода нагнетаемой жидкости Q₁(t) и перепада давлений на измерительном участке ΔP_ИЗМ(t) в реальнoм времени процесса рассчитывают вспомогательные параметры G₁, U₁ по формулам (1). Так, для t=1150 с насосно-компрессорные трубы были полностью заполнены водным раствором сополимера "Комета" и смолы "ДЭГ", при этом устьевые замеры расхода, устьевого давления и перепада давлений на измерительном участке соответственно равны:
Q₁(t)=829,44 м³/сут;
Р_УСТ(t)=13, 614 МПа
ΔP_ИЗМ(t) = 0,01441 МПа.
Тогда из формулы (1) для t= 1150 с вспомогательные параметры G₁, U₁ равны:

где d_ИЗМ=0,062 м - внутренний диаметр измерительного участка;
L_ИЗМ= 6 м - длина измерительного участка (расстояние между осями импульсных трубок на измерительном участке).

Величины вспомогательных параметров G₁, U₁, рассчитанные по формулам (1) для момента t, наносят на график зависимости (фиг.3)
U₁=U₁(G₁ ),
где по оси абсцисс откладывают величины lgG, по оси ординат - величины lgU.

После получения первых 40 значений U₁=U₁(G₁) производят аппроксимацию полученного массива точек корреляционной зависимостью:
U₁=10^-13,981•G₁^1,5525
В дальнейшем по мере поступления новых данных (значений пар G_J и U_J) параметры функциональной зависимости U₁=U₁(G₁) практически не изменились.

После установления вида корреляционной зависимости (2) для каждого замера расхода Q₁(t) гелеобразующего состава в реальном времени процесса по формуле (3) вычисляют вспомогательный параметр

Так, для t=1150 с:

где d_НКТ=0,059 м - внутренний диаметр насосно-компрессорных труб.

Принимая

из корреляционной зависимости (2) U₁= U₁(G₁) определяют величину

Для каждого замера расхода Q(t) в реальном времени процесса по формуле (6) рассчитывают коэффициент λ[(fluid1),Δt] гидравлических сопротивлений потоку в насосно-компрессорных трубах гелеобразующего состава. Так, для t= 1150 с при расходе Q₁(t)=829,44 м³/сут:

где вспомогательный параметр; d_НКТ=0.059^м - внутренний диаметр насосно-компрессорных труб; ρ_УСT(fluid1) = 1000 кг/м³ - плотность гелеобразующего состава в устьевых условиях, кг/м³; Q₁(t)=829,44 м³/сут - расход жидкости в момент t=1150 с; коэффициент гидравлических сопротивлений потоку в насосно-компрессорных трубах гелеобразующего состава при расходе Q₁(t)=829,44 м³/сут.

Определенные из уравнения (6) значения коэффициента λ[(fluid1),Δt] гидравлических сопротивлений потоку в насосно-компрессорных трубах гелеобразующего состава наносят на график зависимости λ [(fluid1),Δt] = Ф(Q₁(t)) (фиг. 4). По оси абсцисс откладывают величины lg Q₁(t), по ординат - величины lgλ[(fluid1), Δt].
После получения первых 40 точек производят аппроксимацию полученного массива. Корреляционная зависимость λ[(fluid1),Δt] = Ф(Q₁(t)) имеет вид
λ[(fluid1),Δt] = 0,61873• Q.
В дальнейшем по мере поступления новых данных параметры функциональной зависимости U₁=U₁(G₁) практически не изменились.

Для каждого замера расхода Q₁(t) по корреляционной зависимости λ[(fluid1),Δt] = 0,61873•Q в реальном времени процесса рассчитывают коэффициент гидравлических сопротивлений в насосно-компрессорных трубах потока гелеобразующего состава (водного раствора сополимера "Комета" и смолы "ДЭГ", образующих псевдопластичную систему).

По этим данным для каждого замера расхода Q₁(t) в реальном времени процесса по формуле (8) рассчитывают потери давления Р_ТР(t) гелеобразующего состава на жидкостное трение в насосно-компрессорных трубах. Так, при расходе Q₁(t)= 829,44 м³/сут потери давления гелеобразующего состава на жидкостное трение равны:

где λ[(fluid1), Δt]- коэффициент гидравлических сопротивлений гелеобразующего состава в насосно-компрессорных трубах при расходе 829,44 м³/сут;

L=2230 м - длина насосно-компрессорных труб от устья до башмака.

Тогда динамическое забойное давление Р_С(t) в текущий момент времени t= 1150 с при расходе 829,44 м³/сут по формуле (9) равно:
P_C (t)=P_УСТ(t)+P_Г(t)-P_ТР(t)=13,614+21,876-7,124=28,366 МПа;
где: P_УСТ(t)= 13,614 МПа - устьевое давление в момент t=1150 с закачки гелеобразующего состава при расходе Q₁(t)=829,44 м³/сут;
P_Г (t)= 21,876 МПа - гидростатическое давление столба гелеобразующего состава;
P_ТР(t)=7,124 МПа - потери давления гелеобразующего состава на жидкостное трение при расходе Q₁(t)=829.44 м³/сут.

Отсюда репрессия на пласт ΔP_C(t) в момент t=1150 с закачки гелеобразующего состава при расходе Q₁(t)=829,44 м³/сут по формуле (10) равна:
ΔP_C (t)=P_C(t)-Р_ПЛ=28,366-14,952=14,414 МПа,
где Р_ПЛ=14,952 МПа - пластовое давление, приведенное к глубине L=2230 м башмака насосно-компрессорных труб.

Пример 2. На нагнетательной линии перед устьем скважины организуют измерительный участок, как в примере 1.

При эксплуатации скважины глубиной 2230 м проводят обработку призабойной зоны с целью изоляции водопритоков. Процесс обработки заключается в закачке в призабойную зону скважины порций гелеобразующего состава (α = 1) и его продавке пластовой водой (α = 2). В качестве гелеобразующего состава используют водный раствор сополимера "Комета" и смолы "ДЭГ", образующих псевдопластичную систему. Закачку ведут с начальным расходом 5,3 л/с.

На устье скважины организуют процесс импульсной нестационарной закачки реагента. Процесс импульсной нестационарной закачки реагента характеризуется значительными колебаниями расхода и давления с хаотичными изменениями величин по амплитуде и частоте. Расход изменяется по амплитуде от 0,084 до 13,6 л/с, по частоте - от 0, 002 до 0,02 Гц. Давление изменяется по амплитуде от 1 до 10÷15 МПа при той же частоте.

В качестве текущей проводимости используют величину, определенную по результатам кратковременного импульсного нестационарного испытания данной скважины на приемистость пластовой жидкостью. Предварительные испытания данной скважины на приемистость пластовой водой показали, что текущая проницаемость пласта k составляет 0,163 мкм², проводимость k•h равна 2,45 мкм²•м, коэффициент скин-эффекта оценен как 12,89. Вязкость пластовой жидкости составляет 1,02 мПа•с, таким образом, гидропроводность пласта ε по формуле (19) равна:

Пьезопроводность пласта χ составляет 0,05 м²/c, радиус скважины r_C равен 0, 084 м.

Процесс обработки заключается в закачке в призабойную зону скважины порций гелеобразующего состава и его продавке пластовой водой. При этом на скважине выполняют замеры и регистрацию устьевых давлений, плотности и объемного расхода нагнетаемых жидкостей с периодом опроса 5 с. Для каждого момента времени t по этим данным в реальном времени процесса рассчитывают потери давления Р_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление P_Г(t) столба нагнетаемого реагента, динамическое забойное давление Р_С(t) по формуле (9), репрессию на пласт ΔP_C(t) по формуле (10), объемный расход реагента Q(t) в забойных условиях. Затем для каждого N замера (в момент времени t_N) по формуле (11) рассчитывают величину функции репрессии Y, а по формуле (13) - накопленный объем реагента W в забойных условиях, поступившего в пласт к этому моменту времени.

Полученные значения Y и W наносят на график (фиг.5). В реальном времени проводят аппроксимацию отдельных участков полученного графика Y=Y(W) прямолинейными отрезками и определяют их наклон. Первый участок соответствует закачке в призабойную зону 6,7 м³ гелеобразующего состава, при этом его наклон B₁ составлял:
B₁=1167,9 МПа•с/м³=1167,9•10⁶ Пa•c/м³,
а коэффициент скин-эффекта S₁ по формуле (15):

Эта величина показывает, что проводимость призабойной зоны в результате закачки 6,7 м³ гелеобразующего состава несколько снизилась в процессе дальнейшего нагнетания реагента наклон второго прямолинейного участка, аппроксимирующего кривую Y=Y(W) в диапазоне 6,8≤W≤8,0 м³, возрос и составил:
B₂=1988,7 МПа•с/м³=1988, 7•10⁶ Па• с/м³.

Величина коэффициента скин-эффекта S₂, соответствующая второму участку с наклоном 1988,7•10⁶ Па• с/м³, равна S₂=28,605.

Полученная величина указывает на проходящий в призабойной зоне скважины процесс закупорки до проектной величины 28-30. В этой связи после закачки в пласт 8,0 м³ гелеобразующего состава на устье прекращают нагнетание реагента в насосно-компрессорные трубы и начинают качать продавочную жидкость.

На графике Y= Y(W) (фиг.5) этому соответствуют участки 3 и 4 с практически совпадающими наклонами 1958, 8 и 2022,2 МПа•с/м³. Участок 3 отражает продавку гелеобразующего состава из насосно-компрессорных труб в призабойную зону продавочной жидкостью (пластовой водой), коэффициент скин-эффекта на этом участке S₃ равен 28,154. Следовательно, проектное снижение проводимости призабойной зоны достигнуто, дальнейшей закупорки призабойной зоны скважины не происходит и закачка гелеобразующего состава прекращена правильно. Это подтверждает коэффициент скин-эффекта на 4 участке: S₄ =29,11, что соответствует проектному показателю.

Непосредственно перед проведением изоляционных работ гидродинамические исследования не проводились, поэтому в известных способах использовалась величина проводимости пласта по результатам ранее выполненных гидродинамических исследований k•h=4,59 мкм²•м. Вязкость нагнетаемого реагента составляла 15 мПа•с.

В результате известные способы показали, что закупорка призабойной зоны не достигнута, а коэффициент скин-эффекта находится в пределах [-0,5--0, 15].

Пример 3. На нагнетательной линии перед устьем скважины организуют измерительный участок как в примере 1 и определяют гидропроводность пласта.

Проводят закачку пластовой жидкости в нефтедобывающую скважину глубиной 2240 м. Для оценки точности определения гидропроводности пласта предлагаемым способом проводят предварительные исследования скважины методом восстановления давления, в соответствии с которым гидропроводность пласта составляет 10,2 мкм²•м/(мПа•с). Таким образом, для оценки точности определения гидропроводности пласта принято:
ε_ИСТ= 10,2 мкм²• м/(мПа•с).
До проведения операции на скважине задают произвольный ряд значений гидропроводности пласта ε_m
1 мкм²•м/(мПа•c) ≤ ε_m ≤ 30 мкм²•м/(мПа•c);
заведомо включающих истинную величину гидропроводности пласта
ε_ист= 10,2 мкм² •м/(мПа•c).
Процесс определения гидропроводности пласта заключается в закачке в пласт 3 м³ пластовой жидкости. Закачку ведут с начальным расходом 5,8 л/с. На устье скважины организуют процесс импульсной нестационарной закачки пластовой жидкости с колебаниями расхода от 5,2 до 6,4 л/с и частотой 0,02 Гц, аналогично изменяется давление нагнетания. При этом на скважине выполняют замеры и регистрацию с интервалом в 5 с устьевых давлений, плотности и объемного расхода нагнетаемой пластовой жидкости. Для каждого замера в момент времени t по этим данным в реальном времени процесса рассчитывают потери давления Р_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление P_Г(t) столба нагнетаемой жидкости, динамическое забойное давление Р_С(t) по формуле (9), репрессию на пласт ΔP_C(t) по формуле (10).

Рассчитывают объемный расход нагнетаемой жидкости Q(t) в забойных условиях. Затем для каждого из принятых значений ε_m гидропроводности пласта определяют значение производной Δ Y_m/ΔX_m(t_N) по формуле (16). Одновременно с ΔY_m/ΔX_m(t_N) для того же момента времени t_N определяют накопленный объем жидкости в забойных условиях W(t_N), поступившей в пласт к моменту времени t_N с начала закачки, по формуле (13). Полученные значения наносят на график (фиг. 6), где по оси абсцисс откладывают величины W, по оси ординат - соответствующие им значения ΔY_m/ΔX_m(t_N).
Графики производной ΔY_m /ΔX_m(t_N) существенно зависят от принятой величины гидропроводности пласта ε_m. Чем ближе значения ε_m к истинной величине гидропроводности пласта ε_ист, тем ближе графики производной ΔY_m/ΔX_m(t_N) к прямой, параллельной оси абсцисс. Среди полученных кривых зависимости ΔY_m/ΔX_m (t_N) = ΔY_m/ΔX_m[W(t_N)] установлены две линии, ближе остальных отвечающих условию (21):
ΔY/ΔX[t,ε_ИСT] = const.
Далее известным методом последовательного приближения находят величину гидропроводности пласта ε = 10,3 мкм²•м/(мПа•c), при которой производная ΔY/ΔX может быть принята постоянной наилучшим образом. Выбор оптимального выполнения условия (21) достигают численными методами с использованием аппарата практической физики. Величина ε, обеспечивающая выполнение условия (21) наилучшим образом, и является искомым значением гидропроводности пласта, а погрешность ее определения равна 1%.

Пример 4. На нагнетательной линии перед устьем скважины организуют измерительный участок как в примере 1 и определяют параметры призабойной зоны скважины.

Предлагаемый способ реализуют на нефтедобывающей скважине глубиной 2240 м.

Для оценки точности определения параметров призабойной зоны скважины предлагаемым способом до проведения работ по реализации предлагаемого способа на скважине дополнительно выполняют гидродинамические исследования методом восстановления забойного давления и гидропрослушивания, определяя при этом гидропроводность ε,ε_ПЗП и пьезопроводность χ,χ_ПЗП пласта и призабойной зоны, коэффициент скин-эффекта S и радиус зоны загрязнения R_ПЗП :
ε = 10,2 мкм² •м/(мПа•c);
ε_ПЗП= 0,51 мкм²•м/(мПа•c);
χ = 1410 см²/c;
χ_ПЗП= 70,6 см²/c;
S=56,9;
R_ПЗП=1,68 м.

До определения параметров призабойной зоны предлагаемым способом на скважине организуют предварительные исследования с целью уточнения гидропроводности пласта. Для этого проводят импульсную нестационарную закачку пластовой жидкости в пласт. Полученная при этом величина гидропроводности пласта совпадает с результатами гидродинамических исследований. Установлено также существенное загрязнение призабойной зоны (S>20÷30). Поэтому далее реализуют предлагаемый способ, при этом достоверное определение параметров призабойной зоны достигают регистрацией параметров процесса (расхода, плотности жидкости и давления нагнетания) на устье скважины.

Для реализации предлагаемого способа на устье нефтедобывающей скважины организуют основной процесс импульсной нестационарной закачки пластовой жидкости. Основной процесс импульсной нестационарной закачки характеризуется колебаниями расхода от минимальных величин 0,58 л/с, обеспечивающих стабильную закачку с противодавлением на устье скважины, до максимальных 5, 79 л/с, обеспечивающих при этом недопущение развития искусственной трещиноватости в призабойной зоне пласта, что достигается выполнением условия (22) для максимального забойного давления в процессе закачки пластовой жидкости: P_Cмакс<30,6 МПа,
где Р_Смакс - максимально допустимое забойное давление в процессе закачки пластовой жидкости, МПа.

Для получения достоверных результатов выполняют закачку пластовой жидкости на 10 режимах нагнетания с резкой сменой расхода с большего на меньший и наоборот (табл. 1). Продолжительность закачки Δθ на каждом режиме приближенно оценивают по формуле (23):

для реализации предлагаемого способа продолжительность закачки на каждом режиме принята Δθ = 200c (табл.1).

Таким образом, на основе выполненных оценок на устье организуют основной процесс импульсной нестационарной закачки пластовой жидкости с резкой сменой расходов от максимальных до минимальных значений и обратно (табл.1) каждые 200 с так, что кривая переменного расхода образует некоторую ступенчатую функцию (24) времени t нагнетания (табл. 1).

При закачке пластовой жидкости для определения параметров призабойной зоны на скважине замеряют и регистрируют с интервалом в 10 с (т.е. с периодом опроса 10 с) устьевое давление, плотность и объемный расход нагнетаемого реагента. Для каждого замера в момент времени t по этим данным в реальном времени процесса рассчитывают потери давления Р_ТР(t) на жидкостное трение в насосно-компрессорных трубах, гидростатическое давление Р_Г(t) столба нагнетаемой жидкости, динамическое забойное давление Р_С(t) по формуле (9). Далее рассчитывают репрессию на пласт ΔP_C(t) по формуле (10). Для каждого текущего N замера, сделанного в текущем временном интервале

на Z режиме закачки, рассчитывают значение функции репрессии Ψ_Z(Δ t_Z), характеризующей нестационарное течение нагнетаемой жидкости в призабойной зоне скважины на данном режиме.

Расчеты выполняют последовательно для всех выполняемых замеров устьевых параметров. Для каждого режима Z закачки на основе выполненных замеров устьевых параметров и расчетов строят график, где по оси абсцисс откладывают величины lnΔt_Z, по оси ординат - соответствующие данному временному интервалу Δt_Z величины функции репрессии Ψ_Z(Δt_Z). На фиг.7 представлены графики функции репрессии Ψ_Z(Δt_Z), характеризующей нестационарное течение нагнетаемой жидкости в загрязненной призабойной зоне скважины в течение Z(Z=1; 2;... 9; 10) режима закачки, последовательно для всех 10 режимов импульсной нестационарной закачки пластовой жидкости в нефтедобывающую скважину при определении параметров призабойной зоны пласта предлагаемым способом в примере 1. Z=1; 2; ... 9; 10 - графики зависимости функции Ψ_Z(Δt_Z) = Ф(lnΔt_Z) от времени Δt_Z закачки пластовой жидкости в нефтедобывающую скважину на Z режиме с расходом

(табл. 1).

Таким образом, каждому из 10 режимов основной закачки соответствует своя линия (фиг.7). На каждом из 10 полученных графиков в интервале 20c ≤ Δt_Z ≤ 140c выделяют начальный наклонный прямолинейный участок, отражающий нестационарное течение нагнетаемой жидкости в загрязненной призабойной зоне на данном Z режиме закачки и описываемый уравнением прямой линии (27) с весьма высокими (0,96÷0, 99) коэффициентами корреляции. Методом наименьших квадратов находят наклон b_Z и начальный участок a_Z выделенного прямолинейного отрезка на каждом Z режиме закачки (табл. 1). Далее определяют гидропроводность ε_ПЗП и пьезопроводность χ_ПЗП призабойной зоны пласта по формулам (28), (29):

Аналогично определяют для остальных режимов закачки (табл. 1). Используя найденные коэффициенты a_Z, b_Z, a_Z-1, b_Z-1 для двух нежных режимов Z, Z-1 закачки, по формулам (30):

после чего по формуле (31) рассчитывают радиус R_ПЗП зоны загрязнения:

Аналогично определяют S и R_ПЗП для остальных режимов нагнетания. Итоги определений параметров призабойной зоны на 10 режимах закачки приведены в табл. 1, здесь же даны средние величины параметров.

Из сопоставления с результатами гидродинамических исследований скважины известным методом восстановления давления следует, что точность предлагаемого способа вполне достаточна для использования в нефтепромысловой практике. В рассматриваемом случае погрешности определения гидропроводности и пьезопроводности призабойной зоны составляют 7,4%; коэффициента скин-эффекта 5,6%; радиуса зоны загрязнения 3,6%.

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

Реферат

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

Формула