Анализ одномерных потоков при нелинейных законах фильтрации

Анализ

Анализ

Анализ

Dр_к

Рис. 3.4. Индикаторная диаграмма в случае плоскорадиального течения несжимаемой жидкости в недеформируемом пласте по закону Дарси

1. Дебит Q не зависит от r, а только от депрессии D р_к. График зависимости Q от D р (рис.3.4) называется индикаторной диаграммой, а сама зависимость – индикаторной. Отношение дебита к депрессии называется коэффициентом продуктивности скважины

. (3.11)

2. Градиент давления и, следовательно, скорость u обратно пропорциональны расстоянию (рис.3.5) и образуют гиперболу с резким возрастанием значений при приближении к забою.

Рис. 3.6. Распределение давления по радиусу

Рис. 3.5. Зависимость градиента давления и скорости от расстояния до центра скважины

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

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

5. Дебит слабо зависит от величины радиуса контура r_к для достаточно больших значений r_к /r_c, т.к. r_к /r_c входят в формулу под знаком логарифма.

Течение совершенного газа через недеформируемый пласт. Выражение для потенциала (2.5) запишется в виде

.

Выпишем соотношения для:

Ø распределения потенциала ;

Ø распределения градиента потенциала ;

Ø дебита ;

Ø средневзвешенного давления .

В вышеприведенных соотношениях: .

Для определения закона движения частиц жидкости проинтегрируем уравнение движения по времени от 0 до t и по расстоянию от r₀ до r, где r₀ – начальное положение частицы флюида.

Переходя в вышеприведенных соотношениях от потенциала к давлению, получим искомые выражения, позволяющие провести исследование в физических переменных (табл. 3.3).

Таблица 3.4

Закон фильтрации
Распределение давления Р =р^2
Градиент давления
Уравнение притока
Уравнение движения

Распределение давления. Если сравнить распределения давления в случае потока газа с соответствующим распределением для однородной несжимаемой жидкости (рис. 3.7), то увидим, что для газа давление вблизи стенок скважины изменяется более резко, чем для несжимаемой жидкости. Пьезометрическая кривая для газа имеет, следовательно, более пологий характер на большем своём протяжении, чем кривая несжимаемой жидкости; однако у неё более резкий изгиб у стенки скважины, чем у кривой несжимаемой жидкости.

Уравнение притока (уравнение индикаторной линии). Индикаторная зависимость для газа описывает параболическую зависимость дебита Q_ст от депрессии Dр_k (рис.3.8) и линейную зависимость дебита от разницы квадратов пластового и забойного давлений в отличие от индикаторной зависимости для несжимаемой жидкости, где устанавливается линейная связь дебита с депрессией. Уравнение притока устанавливает линейную связь между дебитом и разностью квадратов контурного и забойного давлений, поэтому для простоты исследований индикаторная диаграмма при фильтрации идеального газа по закону Дарси строится в координатах Q_ст ~(р_k²-р_с²). В этом случае имеем прямую линию (рис.3.9), проходящую через начало координат с угловым коэффициентом

. (3.12)

Рис. 3.9. Индикаторная завиимость при фильтрации газа по закону Дарси в переменных Q – Dp²

Рис. 3.8. Индикаторная зависимость при фильтрации газа по закону Дарси

Запишем уравнение притока в координатах Q_ст ~ (р_к-р_с). Так как Q_cт= a (р_к²-р_с²), а разность квадратов р_к²-р_с²=2р_к D р_с - (Dр_с)², где D р_с= р_к - р_с, то

.

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

Распределение градиента давления. Градиент давления вблизи забоя резко возрастает как за счёт уменьшения r, так и за счёт падения давления р, вызванного сжимаемостью газа.

Изменение скорости фильтрации получим из закона Дарси

. (3.13)

Из (3.13) видно, что скорость фильтрации слабо меняется вдали от скважины и резко возрастает в призабойной зоне.

Реальный газ и недеформируемый пласт. Следует использовать при давлении р_пл >10МПа и депрессии на пласт р_с/р_к <0.9.

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

р = zr R T. (2.30)

или для изотермического течения газа

, (3.14)

Потенциальная функция имеет вид

. (3.15)

где `z = (z<sub>c</sub>+z<sub>к</sub>) / 2; `μ = (μ_c+μ_к) / 2; z_с =z(p_с), μ_с =μ (p_с), z_к =z(p_к), μ_к =μ (p_к).

Подставив в (3.9) выражение потенциала (3.15) и перейдя от массового дебита к объёмному, приведённому к стандартным условиям, получим уравнение притока:

. (3.16)

Полученное выражение для дебита реального газа отличается от выражения для совершенного газа среднепластовыми множителями `h и ` z. Если сравнить расчётные значения, то можно заметить, что дебиты реального газа ниже дебитов совершенного при тех же условиях. Для тяжелых углеводородов дебит природного газа может составлять всего лишь 72% дебита совершенного.

Течение несжимаемой жидкости в трещиноватом (деформируемом) пласте. Для данных условий потенциал

(3.17)

и основные зависимости имеют вид

· распределение давления

(3.18)

· градиент давления

(3.19)

· объёмный дебит

, (3.20)

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

· скорость фильтрации

. (3.21)

При малых депрессиях на пласт из-за малости b* можно считать, что

и тогда зависимость для давления (3.18) переходит в вид, аналогичный распределению давления в недеформируемом пласте.

При b* =0, т.е. для недеформируемого трещиноватого пласта, после раскрытия неопределённости в формуле (3.20) получаем формулу Дюпюи.

Рис. 3.10. Кривые распределения давления: 1– недеформируемый пласт 2 – трещинный пласт

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

2. Из формулы для объёмного дебита (3.20) следует, что индикаторная кривая – парабола четвёртого порядка с координатами вершины:

. (3.21)

Рис. 3.11. Вид индикаторной кривой при фильтрации несжимаемой жидкости в трещиноватом пласте

Парабола проходит через начало координат, симметрична относительно оси, параллельной оси дебитов; вторая ветвь смысла не имеет (рис.3.11). Однако если учесть реальные пластовые условия (полного смыкания трещин не происходит, т.к. не учитываются факторы, связанные с изменением характеристик течения из–за изменения раскрытия трещин в направлении потока), то можно говорить только о приближённом выполнении экстремальных условий (3.21).

3. Комплексный параметр b * можно определить или графо-аналитически или непосредственно из (3.21), взяв по индикаторной кривой два известных значения дебита Q₁ и Q₂ при двух значениях депрессии D р_с1, D р_с2, т.е. из соотношения

. (3.22)

По найденному значению b * можно из уравнения (3.21) определить проницаемость k_0т.

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

. (3.23)

Подобно тому, как в случае однородной несжимаемой жидкости существует линейная зависимость между потенциалом j и давлением р, так и в установившимся потоке малосжимаемой жидкости существует линейная зависимость между j и плотностью r. Это означает, что для упругой жидкости зависимость между rи координатой r выражается точно теми же формулами, какими выражается зависимость между р и r при однородной несжимаемой жидкости. Чтобы найти зависимость для давления подставим в уравнения, связывающие переменные r и r, значения r, r_к и r_с, определяемые уравнением состояния (2.27). Тогда для плоскорадиального течения имеем

. (3.24)

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

Массовый дебит для упругой жидкости определяется из (3.5) при подстановке j из (3.23)

. (3.25)

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

. (3.26)

Пренебрегать сжимаемостью жидкости в установившемся потоке можно только при условии достаточно малой величины коэффициента b_f и не очень большого перепада давления D р_с = р_к - р_с. В этом случае можно, как для несжимаемой жидкости, считать постоянным вдоль потока не только массовый дебит, но и объёмный. В противном случае, вдоль потока: постоянен только массовый дебит; массовая скорость фильтрации изменяется по тому же закону, что скорость фильтрации для несжимаемой жидкости.

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

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

, (3.27)

где .

Несжимаемая жидкость в недеформируемом пласте. Выразим скорость фильтрации через дебит Q: u=Q / (2p rh)

и перепишем выражение (3.27) в виде

. (3.28)

Отсюда, разделяя переменные и интегрируя, в первом случае, по радиусу от r до R_к и по давлению от р до р_к, а, во втором случае, по радиусу от r_с до R_к и по давлению от р_с до р_к, получаем:

· распределение давления в пласте

; (3.29)

· дебит скважины

. (3.30)

Дебит находится как положительный корень квадратного уравнения (3.29). Из данного уравнения видно, что индикаторная линия – парабола. Кривая распределения давления (3.29) – гипербола и воронка депрессии – гипербола вращения. Крутизна воронки депрессии у стенки скважины будет больше, чем у чисто логарифмической кривой при течении по закону Дарси.

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

. (3.30)

Подставим выражение (3.30) в (3.27) и, заменив плотность по уравнению состояния (3.14), получим:

. (3.31)

Разделив переменные и проинтегрировав в пределах р – р_с и r – r_c получим:

. (3.32)

Распределение давления по (3.32) отличается от распределения давления по закону Дарси наличием последнего члена, что диктует более резкое изменение давления в призабойной зоне.

Интегрируя уравнение(3.31) в пределах р_к - р_с и R_к - r_c, получаем выражение для притока при пренебрежении 1/R_к по сравнению с 1/r_c:

, (3.33)

или в общепринятом виде

. (3.34)

Уравнение (3.34) – основное уравнение, используемое при разработке газовых и газоконденсатных месторождений, так как определяет приток газа к скважине. Коэффициенты А и В определяют по данным исследования газовых скважин при установившихся режимах.

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

, (1.46)

где ; l_бл – средний линейный размер блока.

Умножим все члены (1.46) на плотность r и вынесем за скобки вязкость m. Тогда применительно к плоскорадиальному потоку получим:

, (3.35)

где .

После разделения переменных и интегрирования (3.35) в пределах r_c - r_к; j_с - j_к получим

, (3.36)

Если в (3.36) подставим выражение для трещинной проницаемости и выразим массовый дебит через объёмный, то будем иметь окончательное выражение

. (3.37)

Как видно из (3.37), индикаторная кривая в этом случае определяется в результате сложения двух парабол – параболы четвёртого порядка, симметричной относительно оси, параллельной оси дебитов, и параболы второго порядка (относительно дебита Q) симметричной относительно оси, параллельной оси депрессий (Dр_с) и отстоящей от последней на расстоянии, равном

.

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

(3.38)

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

    а б Рис. 3.12. Аппроксимация индикаторной диаграммы различными уравнениями притока: . Q=0,0972∆p – линейный закон фильтрации, без скин-эффекта; . Q=0,132∆p -12,432 – линейный закон фильтрации, со скин-эффектом; . ∆p=0,0001Q²+0,04 Q – нелинейный закон фильтрации

Из приведенных рисунков видно, что все аппроксимации находятся в области точности, удовлетворяющей точности, принятой при обработке гидродинамических исследований. В то же время, в первом случае мы имеем расчетную проницаемость k = 0,25 дарси, во втором – 0,19 дарси, а в третьем – 0,61 дарси. Таким образом, получаем, что по одним и тем же промысловым данным мы, если не сделать предварительно анализ вида течения, получим проницаемости пласта отличающие в несколько раз. Следовательно, и в прогнозируемой продуктивности пласта мы ошибемся в несколько раз. Если же, в результате мероприятий по интенсификации притока изменится тип коллектора, то, считая его неизменным, можно получить результаты ещё более отличающие. Отсюда следует, что применение даже очень совершенных расчетных методик может привести к неправильным результатам без предварительной оценки вида течения и коллектора, так как любая программа подбирает необходимое уравнение притока по заданной точности, а часто отличия могут крыться в области, принятой за достаточно точную.

Дата добавления: 2014-11-25