Моделирование установившегося циркуляционного движения жидкости в прифильтровой зоне скважины
Разработка математической модели установившегося циркуляционного движения жидкости в скважине, фильтр которой разделен на нагнетательные и всасывающие секторы. Расчет давления и скорости жидкости при циркуляции в любой точке прифильтровой зоны скважины.
Рубрика | Геология, гидрология и геодезия |
Вид | статья |
Язык | русский |
Дата добавления | 26.02.2019 |
Размер файла | 1,5 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Моделирование установившегося циркуляционного движения жидкости в прифильтровой зоне скважины
Ведение
В настоящее время для очистки фильтров водозаборных скважин применяют импульсные, реагентные и комбинированные методы. Совершенствованием различных способов регенерации фильтров скважин занимались Алексеев В. С., Гаврилко В. М., Гуринович А. Д., Коммунар Г. М., Тесля В. Г. И другие. Наиболее эффективным является комбинированный импульсно - реагентный метод восстановления производительности скважин [1, 2, 3]. Реагентный метод предполагает следующие способы регенерации: реагентная ванна, циклическое задавливание реагента и циркуляция реагента в прифильтровой зоне скважины. Циркуляционные способы регенерации фильтров скважин успешно применялись и применяются в России, Германии, США и других странах. Все существующие циркуляционные способы очистки скважин и устройства для их осуществления предполагают движение реагента в прифильтровой зоне за счет разделения фильтра одной или несколькими горизонтальными перегородками или пакерами на секции для закачки и откачки реагента [4, 5, 6, 7, 8]. При таком движении реагента давление в любой точке есть функция от двух координат и в любой момент времени [5]. Однако также возможно осуществление циркуляции реагента путем разделения фильтра одним либо несколькими вертикально расположенными пакерами на нагнетающие и всасывающие сектора. В этом случае давление в любой точке будет являться функцией от угла и радиуса в любой момент времени . Теоретические и экспериментальные исследования такой схемы циркуляции реагента в прифильтровой зоне скважины в литературе отсутствуют. Поэтому возникла необходимость в дальнейшем изучении циркуляционного движения реагента с целью исследования структуры потока жидкости и обоснования применения устройства циркуляционной регенерации (УЦР) для дальнейшего совершенствования технологии циркуляционной регенерации фильтров водозаборных скважин.
Математическое моделирование установившегося циркуляционного движения жидкости. Рассмотрим установившуюся напорную фильтрацию жидкости в прифильтровой зоне скважины, фильтр которой разделен вертикальными пакерами на ряд нагнетательных и всасывающих секторов (Рис. 1).
Скважина предполагается совершенной по степени и характеру вскрытия водоносного горизонта мощностью . Циркуляция жидкости осуществляется за счет того, что в один нагнетательный сектор подается жидкость с расходом , а из одного всасывающего сектора происходит отбор жидкости с расходом . Давление и соответственно скорость будут зависеть от двух переменных радиуса и угла , т.е. . В этом случае процесс установившиеся фильтрации описывается дифференциальным уравнением в частных производных и в цилиндрической системе координат будет иметь вид [9]:
, < < , , (1)
Согласно закону Дарси скорость будет равна:
, и (2)
(3)
Для нахождения распределения давления в прифильтровой зане скважины необходимо решить уравнение (1) при следующих граничных условиях:
Для нагнетательного сектора:
, (4)
Для всасывающего сектора:
, (5)
Условие на контуре питания пласта:
, (6)
(7)
Где - угол нагнетательного сектора; - угол всасывающего сектора; - число нагнетательных секторов равно числу всасывающих секторов, - общее количество секторов (необходимое условие); - радиус скважины; - радиус контура питания; - давление в любой точке пласта; - давление на забое скважины равное пластовому; - статический уровень; - проницаемость водовмещаемой породы; - динамический коэффициент вязкости жидкости; - радиальная составляющая вектора скорости; - угловая составляющая вектора скорости.
Решая уравнение (1) согласно граничным условиям (4)-(7) методом разделения переменных [10] получаем:
(8)
Учитывая, что и формула (8) имеет вид:
(9)
где - коэффициент фильтрации водовмещающей породы; - плотность жидкости; - ускорение силы тяжести; - напор в пласте;
Распределение скоростей, соответствующие распределению давления согласно (2) и (3) будут соответственно равны:
(10)
(11)
Большой практический интерес представляет исследование циркуляционного движения в прифильтровой зоне с меньшей проницаемостью, чем водоносный горизонт. Уменьшение проницаемости есть результат проявления химического, биологического либо механического кольматажа [11]. В зависимости от конкретных условий может наблюдаться и увеличение проницаемости прифильтровой зоны.. Это может быть связано с проведением декольматационных работ или с устройством гравийной обсыпки в процессе бурения скважины. В работе [12] приводится расчет установившиеся радиальной фильтрации воды к скважине с измененной проницаемостью в призабойной зоне при сохранении закона Дарси и отклонении от него.
Рассмотрим циркуляционное движение жидкости в прифильтровой закольматированной зоне скважины с постоянной проницаемостью и радиусом и проницаемостью водоносного горизонта (Рис.2). Тогда уравнение (1) для решения данной задачи примет вид:
, < < , , (12)
и , < < , , (13)
Граничные условия примут вид:
На контуре питания пласта:
, (14)
Для нагнетательного сектора:
, (15)
Для всасывающего сектора:
, (16)
Условие на границе смены проницаемости:
(17)
(18)
Решая уравнения (12) и (13) при вышеперечисленных граничных условиях (14) -(18) и учитывая, что и получаем следующее решение:
(19)
и
(20)
С помощью полученных решений (19) и (20) можно оценить зону распространения давления и скорости в прифильтровой закольматированной зоне проницаемостью и радиусом в процессе установившиеся циркуляции жидкости.
Физическое моделирование циркуляционного движения жидкости
Физическое моделирование осуществлялось на экспериментальной установке состоящей из фильтрационного лотка диаметром 1,22 м, внутри которого, устанавливался фрагмент фильтра водозаборной скважины. Фильтр представлял собой трубчатый полиэтиленовый каркас (по ТУ 51-644-74) диаметром 22,5 см обмотанный проволокой с шагом 1,5 мм. Внутренняя часть фильтра была разделена вертикальным пакером не нагнетательный и всасывающий сектора. Для равномерного питания фильтрационного лотка водой устанавливался кольцевой бьеф. Пространство межу кольцевым бьефом и фильтром послойно с трамбованием заполнялась гравием (ТУРБ 16844-241-2001). Давление в пласте регистрировалось при помощи пьезометров, соединенных с трубками в пласте. Жидкость подавалась при помощи насоса установленного в резервуаре с водой. Для поддержания постоянного статического уровня был установлен напорный бак с холостым сливом. При помощи нагнетательного трубопровода жидкость с заданным расходом от насоса подавалась в фильтр, откуда по отводящему трубопроводу поступала опять в насос, тем самым осуществлялась циркуляция жидкости в прифильтровой зоне скважины. Циркуляционный расход фиксировался при помощи счетчика воды СХВ 15 установленный на нагнетательной линии. Для обеспечения герметизации фильтрационного лотка служила крышка с уплотнителем.
Рис. 3 Схема экспериментальной циркуляционной установки
1 - напорный бак; 2 - холостой слив; 3 - фильтрационный лоток; 4 - фильтр скважины; 5 - крышка фильтрационного лотка; 6 -водовмещающий грунт; 7 - кольцевой бьеф; 8 - трубки к пьезометрам; 9- стойка; 10 - нагнетательный трубопровод; 11 - отводящий трубопровод; 12 - насос; 13 - нагнетательный сектор; 14 - всасывающий сектор; 15 - счетчик воды СХВ 15; 16 - разделительный пакер
Установка работала следующим образом. Вода насосом подавалась в нагнетательный сектор с расходом , который фиксировался счетчиком воды. Далее вода из фильтра поступала в прифильтровую зону скважины, где оттуда циркуляционный поток устремлялся во всасывающий сектор с расходом . Пренебрегая потерями расхода в насосе можно считать, что =. При помощи вентиля, установленного на нагнетательной линии, возможно регулирование циркуляционного расхода. Давление в пласте регистрировалось при помощи пьезометров, трубки которых были установлены в пласте в характерных точках (Рис. 4). Как показали лабораторные эксперименты давление в пласте, при пуске насоса, распределяется и устанавливается за непродолжительный промежуток времени (порядка 5-20 с) что подтверждается в работе [13].
С помощью напорного бака в пласте поддерживался постоянный статический уровень (пластовое давление) равный 1,5 м. Вода для экспериментов была предварительно отстояна при комнатной температуре.
Для проверки сходимости математической модели циркуляционного движения жидкости в прифильтровой зоне были проведено ряд лабораторных экспериментов с однородным грунтом и с кольматирующим кольцом (рис. 1 и 2). Параметры лабораторной установки и условия проведения экспериментов приведены в таблице 1.
Таблица 1 Условия проведения экспериментальных исследований
Параметры установки Усл-ия экс-та |
Коэф. фильтр м/сут |
Длина фильтра, м |
Расход =, м3/с |
Статич уровень, м |
R, м |
, м |
, м |
К - во секторов n, |
|
Весь пласт однородный грунт |
= 1236 |
0,51 и 0,24 |
3,14 х10-4 3,07 х10-4 |
1,5 |
0,5 |
- |
0,113 |
2 |
|
Прифильтровая зона в виде кольматирующ кольца |
= 864 = 1236 |
0,24 |
3,07 х10-4 |
1,5 |
0,5 |
0,213 |
0,113 |
2 |
Сопоставление результатов математического и физического моделирования приведено на рис. 5, где на расчетные графики нанесены опытные данные. Где - это понижение уровня при откачке , и повышение уровня при закачке , т. е. .
Хорошая сходимость результатов математического и физического моделирования свидетельствует о возможности использования данной математической модели при дальнейшем изучении структуры циркуляционного движения жидкости.
Рис. 5 Сопоставление результатов математического и физического моделирования: а) при длине фильтра 0,51 м; б) при длине фильтра 0,24 м; в) при длине фильтра 0,24 м с кольматирующим кольцом.
С этой целью при помощи математического пакета Maple 7 [14] были построены графики иллюстрирующие распределение давления и скоростей в однородной прифильтровой зоне скважины при длине фильтра 0,24 м (рис 6).
Рис. 6 График распределения давления и скорости: а) спектр давлений; б) радиальная составляющая вектора скорости; в) угловая составляющая вектора скорости.
Выводы
1. Разработана математическая модель установившегося циркуляционного движения жидкости в прифильтровой зоне совершенной скважины, фильтр которой разделен вертикальными пакерами на ряд нагнетательных и всасывающих секторов. Предложенная математическая модель позволяет определить давление и скорость жидкости в процессе циркуляции в любой точке прифильтровой зоны скважины с постоянной и с измененной проницаемостью.
2. Исследовано циркуляционное движение жидкости от нагнетательного к всасывающему сектору в лабораторных условиях на экспериментальной установке с однородной и меньшей проницаемостью прифильтровой зоны скважины, а также различной длиной фильтра и расходом.
3. При помощи математического пакета Maple изучен характер распределения давления, а также радиальной и угловой скорости в прифильтровой зоне скважины, что говорит о существовании циркуляционного потока жидкости необходимого для растворения и выноса кольматирующего осадка.
4. Хорошая сходимость результатов математического и физического моделирования свидетельствует о возможности использования математической модели при изучении структуры циркуляционного движения жидкости, а также дальнейшем использовании лабораторной установки для совершенствования технологии циркуляционной регенерации скважин, в частности, для определения конструктивных параметров устройства циркуляционной регенерации.
Список цитированных источников
циркуляционный скважина жидкость движение
1. Ивашечкин, В. В. Опыт применения комбинированных технологий восстановления дебита водозаборных скважин / В. В. Ивашечкин, А. Н. Кондратович, А. М. Шейко // Энерго- и метериалосберегающие экологически чистые технологии: Тезисы докладов 6-й Междунар. Междунар. науч.-технич. конф, Гродно, 1-2 ноябр. 2005 г. / Гроднен. гос. ун-т; редкол.: А. И. Свириденок [и др.]. - Гродно, 2005. - С. 78-79.
2. Ивашечкин, В. В. Газоимпульсная технология восстановления пропускной способности фильтров водозаборных скважин / В. В. Ивашечкин; под ред. А. Д. Гуриновича. - Минск: БНТУ, 2005. - 270 с.
3. Ивашечкин, В. В. Опыт применения импульсно-виброреагентной технологии восстановления дебита водозаборных скважин / В. В. Ивашечкин, А. Н. Кондратович // Инф. Бюллетень научно-произв. ассоциации «Аквабел». - Минск, 2004. -№7. - С. 9-10.
4. Алексеев, В. С. Опыт циркуляционной регенерации водозаборных скважин / В. С. Алексеев, Г. М. Коммунар, В. Г. Тесля // Водоснабжение и санитарная техника . - 1985. - № 9. - С. 9-10.
5. Тесля, В. Г. Технология циркуляционной регенерации скважин / В. Г. Тесля // Сб. науч. тр. / МДНТП им. Ф. Э. Дзержинского. - Москва, 1985. - Вып. 74: Повышение эффективности работы водозаборов из поверхностных и подземных источников. - С. 114- 121.
6. Houben, G. Regenerierung und sanierung von Brunnen / G. Houben, C. Treskatis. - Munchen: Oldenbourg industriever, 2003. - 280 s.
7. Устройство для циркуляционной обработки скважин на воду: а. с. 1182129 СССР, ЕО3В 3/18, Е21В 43/22 / В. С. Алексеев, Г. М. Коммунар, В. Г Тесля; Всесоюзный НИИ ВОДГЕО. - № 3721443/23-26; заявл. 02.02.84; опубл. 30.09.85 // Патент. - 1985. - № 36. - С. 5.
8. Устройство для реагентной обработки скважин: а. с. 1654473 СССР, Е О3В 3/15 / К. Г. Иванищев, В. Е. Воропанов, В. А Попов; Всесоюзный НИИ ВОДГЕО. - № 447788/33; заявл. 24.08.88; опубл. 07.06.91 // Патент. - 1991. - № 21. - С. 4.
9. Маскет, М. Течение однородных жидкостей в пористой среде / М. Маскет. - М. -Л: Гостоптехиздат, 1949. - 628 с.
10. Тихонов, А. Н. Уравнения математической физики: учебник для студ. физико-математических спец. ун-тов / А. Н. Тихонов, А. А. Самарский; Моск. гос. ун-т им. М. В. Ломоносова. - 7-е изд. - Москва: Изд-во Московского ун-та: Наука, 2004. - 798 с.
11. Шейко, А. М. Анализ долговечности водозаборных скважин г.Минска / А.М. Шейко, В.В. Ивашечкин, Н.В. Холодинская, Э.А. Макарова // Вестник БНТУ. - 2006. - № 1. - С. 27-32.
12. Муфтахов, А. Ж. К расчету дебита скважин с измененной проницаемостью призабойной зоны / А. Ж. Муфтахов // Сб. науч. тр. / НИИ ВОДГЕО. - Москва, 1966. - Вып. 13: Гидрогеология. - С. 181- 185.
13. Гандер, В. Решение задач в научных вычислениях с применением Maple и MatLab / В. Гандер, И. Гржебичек. - Минск: Вассамедиа, 2005. - 520 с.
Размещено на Allbest.ru
Подобные документы
Расчет параметров режима работы бурового насоса при прямой промывке нефтяной скважины роторного бурения. Схема циркуляции промывочной жидкости в скважине при прямой промывке. Основные геометрические характеристики участков движения промывочной жидкости.
курсовая работа [1,2 M], добавлен 23.12.2012Распределение давления в газовой части. Уравнение Бернулли для потока вязкой жидкости. Графики зависимости дебита скважины и затрубного давления от проницаемости внутренней кольцевой зоны. Формула Дюпюи для установившейся фильтрации в однородном пласте.
курсовая работа [398,4 K], добавлен 10.01.2015Геологическая характеристика разреза скважины, ее конструкция. Определение количества потребных материалов для приготовления промывочной жидкости с заданными свойствами. Анализ инженерно–геологических условий бурения скважины. Выбор буровой установки.
курсовая работа [124,5 K], добавлен 05.12.2017Напорный приток к дренажной галерее. Приток к совершенной скважине, расположенной в центре кругового пласта. Время движения частицы жидкости, движущейся по радиусу от контура питания к скважине. Стоки и источники. Фильтрация неньютоновских жидкостей.
курсовая работа [538,7 K], добавлен 03.04.2014Анализ работы газовой скважины в пористой среде при установившемся режиме фильтрации газа. Исследование зависимости дебита газовой скважины от ее координат внутри сектора. Диагностика газовой скважины по результатам гидродинамических исследований.
курсовая работа [741,1 K], добавлен 15.04.2015Дебит скважины при частично изолированном контуре питания кругового пласта. Эпюра скоростей вблизи скважины. Динамика фронта частиц, продвигающихся от контура к скважине, являющегося приближенным аналогом линии изосат или фронта воды, замещающей нефть.
контрольная работа [1,7 M], добавлен 25.07.2014Геологическая характеристика зоны дренируемой скважины. Цели и методы гидродинамических исследований пластов. Построение индикаторных диаграмм (зависимости дебита от депрессии) и анализ характера их выпуклости. Уравнение притока жидкости в скважину.
курсовая работа [247,7 K], добавлен 27.01.2016Геологические условия бурения. Расчет плотности растворов. Выбор конструкции скважины и способа бурения, гидравлической программы бурения скважины. Выбор типа промывочной жидкости. Расчет обсадных колонн на прочность. Характеристика бурильной установки.
курсовая работа [74,5 K], добавлен 20.01.2016Геологическая характеристика месторождения. Выбор конструкции, технологии бурения эксплуатационной скважины на Туймазинском месторождении. Расчет цементирования эксплуатационной колонны, расхода промывочной жидкости и программы промывки, потери давления.
курсовая работа [2,4 M], добавлен 14.09.2012Механические методы воздействия в твердых породах. Проведение оценки давления гидроразрыва пласта. Расчет потерь давления на трение в лифтовой колонне при движении рабочей жидкости. Расчет скорости закачивания рабочей жидкости при проведении ГРП.
курсовая работа [248,2 K], добавлен 11.11.2013