Критерии подобия при моделировании процессов
Транспортировки частиц шлама в горизонтальных участках ствола скважины
П |
роцесс транспортировки твердых частиц потоком жидкости в общем виде не имеет решения и требует для изучения постановки экспериментов.
Проведение экспериментов требует выполнения критериев подобия, отражающих краевые условия, характер процессов и геометрию системы.
К краевым условиям при изучении потоков жидкости относится выделение начального и конечного участков трубопровода, на которых формируется (искажается) профиль скорости течения жидкости в канале, и процесс является неустановившимся. Все исследования необходимо проводить вне этих участков. Согласно многочисленным исследованиям длина начального и конечного участков при промывке ньютоновской жидкостью составляет:
для ламинарного течения в трубе
lнач>0,03∙D∙Re, (1)
для турбулентного течения
lнач=(25÷50)D, (2)
где D - гидравлический диаметр трубопровода.
Для неньютоновских жидкостей длина начального участка иногда составляла 80.D и более. При течении жидкости в кольцевом пространстве длина начального участка составляла не менее 50 внутренних диаметров внешней трубы.
В целом можно принять
. (3)
|
|
Непосредственно из-под долота выбуриваемая порода достаточно равномерно поступает по всей площади кольцевого пространства. Однако в процессе движения частицы твердой фазы стремятся к нижней стенке ГС, и на некотором расстоянии от долота создается установившееся распределение выбуренной породы по высоте кольцевого пространства. Длина этого начального участка, на котором распределение концентрации частиц шлама еще является неустановившимся, зависит от свойств и режима течения промывочной жидкости и размеров частиц шлама, нуждается в специальном изучении.
Поток жидкости, движущийся в канале, характеризуется критериями подобия, число которых равно числу реологических параметров жидкости. Для ньютоновской жидкости это число Рейнольдса
, (4)
где U - средняя скорость течения; ρ, μ, ν -соответственно плотность, динамическая и кинематическая вязкости жидкости.
Течение бингамовской жидкости характеризуется двумя критериями подобия – Рейнольдса и Сен-Венана
,
. (5)
Для степенных жидкостей, описываемых реологическим уравнением Оствальда (их часто называют псевдопластичными жидкостями – ППЖ),
|
|
, (6)
где к - постоянная вязкости; т - показатель нелинейности ; имеем также два критерия подобия – показатель нелинейности и параметр Ren:
; π3=m . (7)
Процесс движения твердых частиц в жидкости определяется параметром Архимеда и критериями, характеризующими режим обтекания.
Для ньютоновской жидкости
;
, (8)
где Δρ =ρТ – ρж - разность плотностей твердой частицы и жидкости; g- ускорение силы тяжести; d - эквивалентный размер частицы; UОС - скорость движения частицы относительно жидкости (скорость оседания, скорость витания).
Для жидкости Бингама
;
. (9)
Для жидкости Оствальда
;
. (10)
Геометрические критерии подобия включают соотношения размеров наружной и внутренней трубы, размеров частиц и канала, а также величину относительного эксцентриситета:
;
;
, (11)
где D1 - наружный диаметр внутренней трубы;D2 - внутренний диаметр внешней трубы; - расстояние между центрами внешней и внутренней трубы.
На искривленных участках ствола скважины на частицы шлама действует центробежное ускорение. Даже при малых радиусах искривления ствола скважины (R=20÷ 50 м) и скоростях движения частиц шлама (U=1÷2 м/с), это ускорение, равное
|
|
м/с2 , (12)
значительно меньше g=9,8 м/с2 и может не учитываться.
Наиболее простым с методической точки зрения было бы натурное моделирование, когда все коэффициенты моделирования были бы равны единице. При этом геометрические размеры трубопровода получаются достаточно большими. Имея внутренний диаметр внешней трубы D2=0,2 м, получим суммарную длину начального и конечного участков около 20 м, а всей установки не менее 30 м, что достаточно много.
Примем геометрический коэффициент моделирования
. (13)
Индексы «м» относятся к модели, а индексы «н» - к натуре. Чем меньше КL, тем меньше размеры установок.
Из равенства параметра Рейнольдса для модели и натуры Re=idemпри Кρ=1 и Кμ=1, то есть при использовании в экспериментах той же жидкости, что и в натуре, получим
. (14)
Для жидкости Бингама при Кρ=1 из Re=idem получаем
|
|
, (15)
а из Sen=idem
. (16)
Из (15) и (16) видно, что использование натурной бингамовской жидкости в модели невозможно, так как при Кτо= Кη=1 из (15) получаем, что КU=1/КL, а из (16) КU=КL.
Для соблюдения условия Sen=idem необходимо при Кη=1 иметь Кτо=1/ или Кη=КL при Кτо=1.
Очевидно, что изменять независимо друг от друга τо и η достаточно сложно.
Для жидкости Оствальда из Ren=idem при Кρ=1 и Кm=1 получаем
. (17)
Значения коэффициента КU при различных условиях моделирования для Re=idem представлены в табл.1.
Расчеты показывают, что наименьшее увеличение скорости при уменьшении масштабов геометрического моделирования для выполнения Re=idem достигается при использовании степенной жидкости Оствальда.
Из Ar=idem при Кρ=Кμ=1 имеем для жидкостей Ньютона и Бингама
, (18)
для жидкости Оствальда
. (19)
КL | Жидкость | |||||
Ньютона | Бингама | Оствальда | ||||
КU при Кr= Кm=1 | КU при Кr=1, Кτо=1 | КU при Кr=1 | ||||
Кη=1 | Кη=0,5 | m=0,8 | m=0,6 | m=0,4 | ||
0,5 | 2 | 2 | 1 | 1,6 | 1,34 | 1,19 |
0,33 | 3 | 3 | 1,5 | 2,08 | 1,6 | 1,32 |
0,25 | 4 | 4 | 2 | 2,52 | 1,81 | 1,41 |
0,2 | 5 | 5 | 2,5 | 2,92 | 1,41 | 1,5 |
Если принять, что в натуре Δρ=(1,2÷1,8)·103 кг/м3, то при моделировании необходимо иметь Δρм=(1,2÷1,8)-103·КΔρ.
Отсюда модель твердой частицы должна иметь плотность
ρmм=(1,2÷1,8)·103·КΔρ +ρж . (20)
Расчеты по формуле (20), приведенные в табл.2, показывают, что при использовании жидкостей Ньютона и Бингама Ar=idem можно обеспечить лишь при КL=0,5 , выполняя модели твердых частиц из сплавов свинца. Уже при КL=0,4 модели твердых частиц следует изготавливать из материалов плотностью (18,7÷28)·103 кг/м3 ,что почти невозможно.
2
КL | ρmм 10-3, кг/м3 | |||
Жидкость | ||||
Ньютона, Бингама | Оствальда | |||
m = 0,8 | m = 0,6 | m = 0,4 | ||
0,5 | 10,5÷15,5 | 6 ÷ 9 | 4,3÷6,5 | 3,4 ÷5 |
0,4 | 18,7÷28 | 10 ÷ 15 | 6,6÷10 | 4,7 ÷ 7,1 |
0,33 | - | - | 9,2÷14 | 6,2 ÷ 9,4 |
0,25 | - | - | - | 9,6 ÷ 14,4 |
0,2 | - | - | - | 13,4 ÷ 20 |
При использовании в эксперименте жидкости Оствальда моделирование твердых частиц можно обеспечить при КL=0,4 (если m=0,8) и КL=0,25 (если m=0,4).
Постановка исследований по транспорту шлама существенно упрощается, если промывка скважин осуществляется степенной жидкостью, а режим течения в кольцевом пространстве – ламинарный. При этом становится необязательным выполнение условия , так как все ламинарные течения при
автомодельны.
В данном случае важно соблюсти динамическое подобие. На некотором расстоянии от забоя большая часть выбуренной породы оказывается на нижней стенке горизонтального ствола и передвигается под воздействием потока жидкости, проскальзывая (перекатываясь). Поскольку осаждение частиц шлама при этом отсутствует, отпадает необходимость соблюдения условия ArН= ArМ.
Динамическое подобие означает равенство в натуре и модели отношения всех сходственных сил. На частицу шлама, передвигающуюся по нижней стенке скважины, действуют сила тяжести G, сила Архимеда Far и сила, с которой жидкость увлекает частицу Ft .
Процесс будет подобен при (21)
или ,
но ,
где tст - касательные напряжения в жидкости на стенке канала.
Отсюда для динамического подобия необходимо иметь
. (22)
Из (21) при имеем
. (23)
Соотношение (23) позволяет определить значение касательных напряжений, а следовательно, и расход промывочной жидкости, необходимый для транспортировки частиц шлама в скважине, по данным лабораторных исследований.
(11) Тема 2009-2-99 КП. Математическая модель эффекта Томса.
Манжай В.Н. Количественное описание эффекта Томса и применение его в трубопроводном транспорте нефти // Известия высших учебных заведений. Нефть и газ. 2009. № 2. С. 99 – 106.
УДК 532.517.4 : 532.135
Дата добавления: 2018-02-28; просмотров: 612; Мы поможем в написании вашей работы! |

Мы поможем в написании ваших работ!