Прогрессивные технологии сооружения скважин
ДВИЖЕНИЕ ПОТОКА В ФИЛЬТРЕ
Фильтры устанавливают в скважину с целью очистки добываемого из пласта флюида от песка и других инородных примесей. Фильтр должен пропускать флюид и иметь при этом минимальное гидравлическое сопротивление, надежно предохранять скважину от проникновения твердой фазы, образования пробок и существенного снижения дебита.
В нашей стране и за рубежом разработано огромное число фильтров, что осложняет выработку единого подхода к гидравлике движения потока в каждой конструкции. Большинство зависимостей базируется на известной формуле истечения из затопленного отверстия
О = цпЩ., (1.48)
V Р
Где V - скорость истечения из отверстия; ц - коэффициент расхода; N - скважность фильтра; J5 - перепад давления на фильтре; р - плотность флюида.
И. Петерсен, С. Ровер, М. Альбертсон предложили для описания перепада давления на фильтре использовать закон гиперболических косинусов
ХТ +1
' v (1.49)
Ch J 5 =- Ch |
Хт _ 1
Где т, D^ - соответственно длина и диаметр фильтра; % - гидравлический параметр фильтра, характеризующий его конструктивные особенности.
По данным В. С. Алексеева параметр х равен произведению коэффициента расхода на скважность фильтра цп.
Определение цп по формуле (1.49) может дать удовлетворительные результаты только для одного отдельно взятого отверстия. В случае оценки гидравлического параметра фильтра влияние отверстий друг на друга весьма существенно, что вносит в расчеты погрешность. Э. А. Грикевич рекомендует определять 54
Гидравлический параметр фильтра в лабораторных условиях по следующей формуле:
(1.50) |
0,246- |
Q |
Цп = Дф arcth 566M
В представленной формуле предполагается, что расход Q через фильтр по его длине m остается постоянным. На самом деле это может быть справедливо только для оценки притока в фильтры малой длины. Входные скорости притока в фильтр и расход через фильтр увеличиваются от нижних к верхним отверстиям согласно формуле Э. А. Грикевича
О = p0ch566^m, (1.51)
Где v0 - скорость притока в фильтр ниже рассматриваемого сечения.
С увеличением гидравлического параметра сопротивление фильтра уменьшается. Как показывают многочисленные эксперименты, параметр % изменяется от 0,19 до 0,88. Такой широкий разброс данных для фильтров однотипной конструкции объясняется разнообразием условий работы фильтров в реальных условиях и сложностью их классификации и типизации.
Поэтому прогнозировать реальное значение % в скважинных условиях с достаточной достоверностью сложно. Существенно изменяются гидравлические характеристики фильтра в контакте с породой. В. С. Алексеев получил данные разброса значений гидравлического параметра для фильтров одинаковой конструкции при контакте с породой различного гранулометрического состава (табл. 1.7).
Степень ухудшения гидравлического параметра фильтра при наложении частиц породы на отверстия определяется формой частиц и отверстия, соотношением их размеров. Установлено, что с увеличением окатанности и сферичности частиц плотность их
Таблица 1.7
|
55 |
Сложения уменьшается и соответственно обеспечиваются улучшенные фильтрационные свойства слоя песка, контактирующего с фильтром. При вытянутой, продолговатой форме частиц происходит необратимая закупорка отверстий, уменьшение эффективной скважности фильтра и резкое увеличение его гидравлического сопротивления. Интенсивность закупорки зависит от формы отверстия фильтра, находящегося в контакте с породой. Наиболее сильно эффект закупорки проявляется в сетчатых и проволочных фильтрах, а также некоторых других конструкциях, предусматривающих расширение отверстий в направлении песка.
После установки фильтра в заданный интервал при подъеме инструмента возникает депрессия и порода продуктивного пласта обжимает фильтр; вокруг отверстий образуются своды естественного равновесия. При освоении скважины и стимуляции притока на частицы, которые непосредственно прилегают к отверстиям, действует гидродинамическая сила, способствующая уплотнению частиц и частичной суффозии. Вследствие расширяющейся к породе форме отверстий при уплотнении породы возникает эффект «клина», приводящий к интенсификации снижения пористости и проницаемости прифильтрового слоя породы. Часть частиц породы выносится через отверстие, а другая часть закупоривает поры и цементирует естественный скелет породы, контактирующей с отверстием.
Степень уплотнения породы и цементации отверстия зависит от количества расклинивающих граней. В случае проволочного фильтра таких граней две, а в случае сетчатого покрытия - четыре. Функцию расклинивающих граней может выполнять не только фильтрующее покрытие, но и поверхность уже заклинившихся в отверстиях частиц. При застревании частиц в отверстиях расклинивающий эффект в существенной степени ухудшает гидравлические параметры фильтра.
При каптировании продуктивных пластов, сложенных разно - зернистыми песками, избежать расклинивания породы в отверстиях фильтра очень сложно из-за самых разнообразных соотношений размеров отверстий и частиц. Крупные фракции застревают в отверстиях, мелкие в процессе откачки цементируют пористое пространство, снижая проницаемость контактной зоны. Некоторые конструкции фильтров предусматривают работу отверстий без контакта с породой. Это фильтры с коническими отверстиями, гравитационные и тарельчатые конструкции. Указанные типы фильтров предусматривают образование сводов равновесия из частиц в процессе устойчивой работы. В реальных условиях своды равновесия при изменении режимов эксплуатации периодически разрушаются, что приводит к проникновению 56
Внутрь фильтра песка или гравия и заклиниванию отверстий. В этой связи избежать закупорки отверстий породой сложно, хотя расклинивающий эффект, снижающий проницаемость зоны, возможно предотвратить.
Закупорка отверстий породой наиболее интенсивно происходит в интервалах фильтра, характеризующихся максимальными скоростями притока. Поэтому в процессе освоения скважины и стимуляции притока гидравлическое сопротивление фильтра за счет контактных потерь будет изменяться, увеличиваясь в наиболее обильных интервалах и оставаясь приблизительно постоянным в менее обильных. Процесс закупорки отверстий породой продуктивного пласта способствует выравниванию эпюры входных скоростей по длине фильтра. Однако выравнивание скоростей достигается за счет снижения притока в перспективных интервалах и не способствует получению больших дебитов.
Гидравлический параметр фильтров % = цп считается величиной переменной и изменяющейся по длине фильтра в зависимости от состава, свойств контактирующей породы и режимов притока.
Из гидравлики известно, что коэффициенты сопротивления с Меняются с изменением режима движения потока, характеризующегося числом Рейнольдса (рис. 1.14). С увеличением скорости движения потока коэффициент сопротивления уменьшается. Следовательно, по длине фильтра коэффициент сопротивления или гидравлический параметр будет меняться даже при одинаковом эффекте наложения частиц. Измеренный в лабораторных условиях гидравлический параметр %, рассчитанный по формуле (1.50), носит обобщенный характер и справедлив только для образца модели заданной длины в испытуемых режимах и только при контакте с аналогичной породой. Обеспечить такое подобие почти невозможно, поэтому переносить расчетный по данным эксперимента гидравлический параметр с достаточной степенью достоверности в реальные условия нельзя.
Рассмотрим обобщенную модель движения потока в фильтре в общем виде. Закон сопротивления при обтекании потоком любых тел описывается зависимостью
Fc = флц^ + фтр d 2v2, (1.52)
Где Fc - сила гидродинамического сопротивления; фл, фт - коэффициент соответственно ламинарного и турбулентного сопротивления обтеканию; D - размер структурного элемента фильтра; ц, р - соответственно вязкость и плотность пластового флюида; V - Скорость обтекания.
57
Рис. 1.14. Зависимость коэффициента сопротивления с от режима движения потока при обтекании различных поверхностей цп 1 - 0,125; 2 - 0,22; 3 - 0,6; 4 - 1; 5 - 1,6 |
Перепад давления на обтекаемом элементе фильтра можно определить из выражения (1.52)
Ф лц» 2
(154) |
J5 = + ф Tpv.
Решая квадратное уравнение (1.53), получаем общий вид зависимости скорости от перепада давления на фильтре
V = J^JU + j - i^JL. чфт 2pd J фтр 2 фт pd
При ламинарном потоке скорость истечения через фильтр определяется только линейным членом
V = dJ 5/ ф!^ (1.55)
^фтР |
А при турбулентном истечении - квадратичным членом
(1.56)
Нельзя рассматривать только ламинарный или только турбулентный режим истечения из отверстий в фильтре, так как обычно наблюдается и тот и другой режим в разных сечениях. Ламинарный режим характерен для областей с низкими скоростями фильтрации (в нижних интервалах фильтра), а турбулентный - для интервалов наиболее обильного притока.
(1.53) |
Сопоставляя выражение (1.48) с равенствами (1.55) и (1.56), можно сделать вывод, что % = цп в зависимости от режима 58
Фильтрации будет меняться от значений % = ^фт1 (турбулентный режим) до х = р J5 (ламинарный режим).
Установить значения цп для фильтра данной конструкции сложно. В этом случае оправдано пользоваться коэффициентами сопротивления для ламинарного и турбулентного режимов фильтрации фл фт, каждый из которых в своей области достоверно описывают процесс движения потока в фильтре.
(157) |
Рассмотрим потери напора в фильтре с учетом изменяющихся типовых скоростей потока по длине по закону гиперболических синусов или тангенсов [см. равенства (1.41) и (1.42)]. Уравнение силы сопротивления (1.52) с учетом изменяющихся скоростей имеет вид
Fc =Ф Л т2 v( z) + Ф т рт 4 v 2(z X
Где v(z) - функция скорости по длине фильтра z.
Дифференцирование выражения (1.57) дает следующие выражения:
(1.58) |
Для гиперсинусоидального закона
Dv v а і ( z І
V(z) = —= ——chІ а— J;
Dz m sh — I m
Для тангенциального закона
Dv v V(Z) = -T = - dz m |
1 |
(1.59) |
Tg Y cos2 I Y — M |
Подставляя уравнения (1.58) и (1.59) в равенство (1.57), получаем уравнение распределения силы сопротивления по длине элемента фильтра для типового характера распределения скоростей.
Ch І - — |
По гиперсинусоидальному закону
Ch І - —
Fc = v— |
(1.60) |
Sh а |
Sh а |
Фл цт + фт Pm2v — •
По тангенциальному закону
V y |
V y |
Fc =- |
(1.61) |
Ф лцт + ф трт2 |
Tg Y cos2 І Y — m |
TgYcos2 I Y — m |
59
Перепад давления на элементе фильтра с учетом изменяющихся скоростей притока выражается следующими зависимостями:
Chl а — |
Для гиперсинусоидального закона распределения скоростей по длине фильтра
Chl а —
Ц , — V m
Фл - + Фт Pv а ----------
(1.62) |
J5 = v а- |
Sh с |
M sh а
Для тангенциального закона распределения скоростей по длине фильтра
Ют pv y |
V у |
J 5 =- |
(1.63) |
Y л ^ +
Tg ycos2 [ y — M |
M
Tg ycos2 [ y — M
В интервале ламинарного притока распределения перепада давления на элементе фильтра по его длине удовлетворительно описываются следующими уравнениями:
Y |
(164) |
J5 =
Mtg y cos2 [ y — M
Или
Юл цю аЛІ а —
(1.65) |
(1.66) |
J5 = V m
MSh а
Для турбулентного притока уравнения (1.64) и (1.65) имеют
Вид
—2 2
Т = VXPv Y
J 5 = -------- 7Т~г
(167) |
Sh2а |
J5 = vTpv 2а2 |
Если принять, что в приведенных уравнениях v - средняя скорость притока, то можно выразить ее через расход Q и площадь поверхности фильтра 2nr4m. Заменяя v на расход и площадь поверхности фильтра и интегрируя выражение (1.62) и 60 |
Tg ycos l y— V m
Ch2 [ аz
(1.63) по радиусу для радиального притока к фильтру, получаем уравнения для определения перепада давления на фильтре любой конструкции
SochI—— |
Q |
J 5 = — |
(1.68) |
2nm2sh — |
Г5 |
YQo |
Ch I — —
ФлЦІП — + фтР—
2п shc
J 5 = —
2nm[1]Tg у cos2 I y —
Где r4, r5 - соответственно больший и меньший радиусы фильтров, разность которых соответствует толщине фильтра.
Основными параметрами, характеризующими сопротивление фильтра в уравнениях (1.68) и (1.69), считаются коэффициенты ламинарного и турбулентного сопротивлений обтеканию фл и фт. Значения параметров фл и фт, легко получить в лабораторных условиях, а затем использовать в реальных условиях. В этом случае применение двух коэффициентов позволяет охватить практически всю область режимов фильтрации в фильтрах и избежать погрешностей в расчетах. Проводя эксперименты на фильтрационных лотках при малых скоростях притока на всей длине образца при значениях числа Рейнольдса меньше критических и зная перепад давления на фильтрах J5, получаемый по разнице показаний пьезометров, с помощью формулы (1.55) определяют коэффициент ламинарного сопротивления обтеканию фл. При установлении на лотке на всей длине образца турбулентного режима фильтрации в фильтре по формуле (1.56) находят коэффициент турбулентного сопротивления фт. Параметры неравномерности притока по длине фильтра — и y определяются по методике Э. А. Грикевича на основании экспериментов.
Значения гидравлических сопротивлений фл и фт могут быть найдены из реальной расходограммы. Основной закон сопротивления (1.53), включающий члены, пропорциональные скорости и квадрату скорости, имеет два неизвестных коэффициента сопротивления фл и фт и скорость. Перепад давления на фильтре изменяется пьезометрами. Подставляя в уравнение (1.53) реальный
61
Закон распределения скоростей притока по длине филмра, получаем
|
Shl a — |
Shl a —
Фл^ + Фі Pv -
J5(z) = |
(1.70) |
Sh a |
Sh a |
AT4
(1.71) |
Устанавливая ламинарный режим фильїрации на модели на всей длине образца, при которой квадратачныш член в (1.70) близок к нулю, получаем уравнение для определения коэффициента ламинарного сопроіивления
Ф = _A_T4J5sh a
JovshI a—
Если иніенсифицироваїь приток на модели и усіановиїь іур - булен1ный режим фильїрации по всей длине образца фильїра, іо линейный член в уравнении (1.70) будеї близок к нулю. Тогда коэффициен! іурбуленіного сопроіивления фильїра
J5sh a |
Фі =- |
(172)
(1.73) |
(174) |
Фі =- |
J5tg a |
_2 21 — PV tg | Y— M ПарамеIры неравномерности притока y и a зависяі оі сопроіивления фильіра, і. е. чем больше сопроіивление, іем более равномерный наблюдается приток. Поэтому правильнее в качестве гидравлического парамеіра фильіров использоваіь произведение коэффициентов сопроіивления на функцию распределения скоростей притока на моделях 62 |
—2 і 2 l — pv sh I a—
В случае, когда приток в фильір лучше описывается іанген - циальным законом выражения (1.71) и (1.72), можно записаіь
Ф = 2t4J5tg Y .
Цс tg | y— M
Tg2 И,
X т =ф = J ; (1.75)
Tg2Y pv2
Tg I Y — I
L m I 2R4 J5 nc.
X ё = Фл ;g 1 = (1.76)
Tg Y Н»
При решении задач по определению перепада давления на фильтре неравномерность притока сложно учитывать введением в выражения (1.75) и (1.76) тангенциальных и гиперсинусоидальных сомножителей. В этом случае может быть применен параметр, равный произведению коэффициентов ламинарного и турбулентного сопротивлений на параметр неравномерности притока — или y. При ламинарном режиме притока коэффициент сопротивления умножается на параметр неравномерности — или Y в первой степени и при турбулентной фильтрации - на их квадрат. В выражениях для определения перепада давления на фильтре (1.62)-(1.69) присутствуют параметры хт = флY2 = фл—2 и Хт = ФлY = фл—, что упрощает расчеты.
Фильтр редко состоит из однотипных по длине и толщине элементов обтекания. Обычно конструкция включает фильтрующую оболочку, каркас и промежуточные звенья. Поэтому теоретически правильнее было бы выделить из обобщенных коэффициентов или гидравлических параметров составляющие сопротивлений фильтрующей оболочки, каркаса, промежуточных звеньев и учесть влияние одних элементов на другие, стесненные условия обтекания. Однако на практике это не приведет к существенному шагу вперед из-за сложности точного теоретического и экспериментального выделения разных составляющих сопротивления и оценки их реального вклада в общий результат.