Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

FINITE-ELEMENT MODEL OF FILTRATION LIQUID TO A WELL IN A DEFORMABLE FORMATION

Azhikhanov N.T. 1 Marasulov U.A. 1
1 International Kazakh-Turkish University A. Yasavi
1029 KB
At the present stage is not enough published works on theoretical and practical study of the tasks of filtration to horizontal wells in a deformable medium. In this paper, at the analysis of the effectiveness of a horizontal wells accounting geometric configuration of the plane stratification. The analyzed elastic static stress and strain state of differently oriented wellbore formation in transtropic layer depending elastic anisotropy, including the degree of discontinuity flaw clutch of small inclined layers at an angle φ, when the longitudinal axis of the cavities constitute an arbitrary angle ψ with the line extent plane of isotropy coinciding with the plane of the slits. Developed the method of calculation of a filtration of liquid to differently oriented horizontal well is for calculation of productive characteristics in the deformable transversal-isotropic medium with a tilted angle of the plane of isotropy. These calculated indicators can be used in engineering calculations when developing oil field in fractured anisotropic environment.
filtration
finite element
horizontal well
deformation
permeability
1. Erzhanov Zh.S., Aitaliev Sh.M., Masanov Zh.K., Ustoichivost’ gorizontal’nyx vyrobotok v naklonno-sloistom massive. Alma-Ata, Nauka, 1971, 160 р.
2. Zhumagulov B.T., Masanov Zh.K., Azhikhanov N.T., Calculation of production of horizontal well of type stretch in the anisotropic porous medium. Journal of Mathematics and Technology, ISSN: 2078-0257, February, 2010. рp. 86–90
3. Fadeev A.B. Metod konechnyx e’lementov v geomexanike. – M.: Nedra, 1987. – 221 р.
4. Shestakov V.M. Gidrogeomexanika. – M.: Izd-vo MGU. 1998. – 72 р.
5. Borja R., Chao Li., Finite Element Formulation of Poro-Elasticity Suitable for Large Deformation Dynamic Analysis.,Proceedings, Third UJNR Workshop on Soil-Structure Interaction, March 29–30, 2004. – 9 p.

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

Постановка задачи

Горное давление на пласт компенсируется как напряжениями твердого скелета массива, так и давлением жидкости. Изменение последнего возмущает напряженно-деформированное состояние пласта (НДС), т.е. изменение порогового давления в одной точке вызывает перестройку НДС всей системы и в том числе деформацию во всем пласте.

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

Рассмотрим упругое статическое состояние горизонтальной скважины, продольная ось которой составляет произвольный угол с линией простирания плоскости изотропии породного массива (рис. 1).

Введем прямоугольную декартову систему координат Оxyz таким образом, что ось Оz направлена вертикально вверх, горизонтальные оси Оx и Оy совпадают с линиями соответственно вкрест и вдоль простирания плоскости изотропии.

Упругое состояние трансверсально-изотропного массива описывается уравнением обобщенного закона Гука в системе координат Ox΄y΄z΄, полученной путем поворота Оxyz на угол φ вокруг вертикальной оси Oz, и имеет вид [2]

agim01.wmf

где

agim02.wmf

agim03.wmf

agim04.wmf

pic_1.tif

Рис. 1. Расчетная схема наклонного трансверсально-изотропного массива с горизонтальной скважиной

Здесь коэффициенты деформации dij определяются из [3].

Компоненты деформации также можно определить через перемещение u, v и w (по оси Ox, Oy и Oz соответственно) с помощью cоотношения Коши.

Граничные условия зададим в виде

u = v = 0, при

agim05.wmf

u = v = w = 0 при agim06.wmf.

Далее нестационарная фильтрация жидкости к горизонтальной скважине в трансверсально-изотропной пористой среде (рис.1) описывается следующим уравнением, определяющим давление при пространственной фильтрации:

agim07.wmf

с граничными условиями

agim08.wmf

где agim09.wmf – коэффициенты проницаемости анизотропного (трансверсально-изотропного) пласта.

При этом полное напряжение [4] трансверсально-изотропного пласта может быть выражено через эффективное напряжение и давление в виде

agim10.wmf

agim11.wmf

agim12.wmf

Характерной особенностью модели является предположение о том, что пористая матрица деформируется совершенно свободно до некоторого жесткого предела ε0.

Конечно-элементная модель

Рассмотрим трехмерный изопараметрический шестигранный конечный элемент первого порядка [1], для которого интерполяционный полином является линейной функцией локальных координат ξ, η, ζ (рис. 2). Пределы изменения локальных координат для всех элементов составляют agim13.wmf agim14.wmf agim15.wmf Функции формы такого элемента даются соотношениями вида

agim16.wmf

где ξi, ηi, ζi – локальные координаты узловой точки с номером agim17.wmf.

Интерполяционные соотношения для перемещений в трехмерном изопараметрическом элементе имеют вид

agim18.wmf

где uα, vα, wα – компоненты вектора перемещений в узлах «e»-го элемента.

agim19.wmf – единичная матрица.

Интерполяционные соотношения для конечного элемента определяются как

agim20.wmf

здесь xα, yα, zα – глобальные координаты узловой точки с номером α.

Эти соотношения можно представить в виде

agim21.wmf

где N(e) – матрица функции форм элемента; X(e), Y(e), Z(e) – векторы узловых значений глобальных координат.

Матрица жесткости конечного изопараметрического элемента имеет вид:

agim22.wmf

Здесь матрица градиентов agim23.wmf определяется из функции формы Nα.

pic_2.tif

Рис. 2. Трехмерный изопараметрический элемент

В итоге матрицу жесткости изопараметрического элемента получим из

agim24.wmf,

численное интегрирование которого осуществляется методом Гаусса, т.е.

agim25.wmf

Здесь ξi, ηj, ζk – точки интергрирования Гаусса; n – порядок интегрирования; Hi, Hj, Hk – весовые коэффициенты.

Перемещение точек элемента определяем u(e) = N(e)U. Матрица конечной деформации имеет вид ε(e) = B(e)u(e), в свою очередь компоненты напряжений вычисляются из σ(e) = Dε(e).

Искомые значения определяются из решения систем линейных алгебраических уравнений [5]

KU = F,

здесь

agim26.wmf

где agim27.wmf – векторы узловых сил конечного элемента статически эквивалентны действию конечной деформации, agim28.wmf поверхностных и agim29.wmf объемных сил.

agim30.wmf

Рассмотренную выше матрицу жесткости элемента можно записать в виде

agim31.wmf

Аналогичным путем применяем процедуру для шестигранного изопараметрического элемента (рис. 2) к решению задачи фильтрации.

Вычислительный эксперимент

Деформируемое состояние наклонного под углом φ трансверсально-изотропного массива приводится с применением закона Гука и коэффициентов деформации. Численный эксперимент проводился по следующим данным: в качестве пород наклонных слоем взяты [2], аргиллит, алевролит, песчаник, известняк, модуль упругости которых имеет соответственно значения E = 1,34; 0,62; 2,95; 5(10-5 кг/см2) постоянного Пуассона соответственно υ(k) = 0,3; 0,2; 0,35; 0,11. Для трансверсально-изотропных слоев округленные упругие характеристики определяются E1 = 1,54·10-5 кг/см2; E2 = 0,98·10-5 кг/см2; модуль сдвига G2 = 0,36·10-5 кг/см2; υ1 = 0,22; υ1 = 0,25. При этом граничные условия учитываются при определении напряжений в наклонном под ψ = 30° пласте по X и Z. Сравнительные изолинии показаны на рис. 3 в случаях φ = 30° (пунктирные линии) и φ = 60° (сплошные линий). Общий вид функции деформации представлен на рис. 4.

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

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

pic_3.wmfа

 pic_4.wmfб

Рис. 3. Изолинии нормальных напряжений: а – σx; б – σz

 pic_5.wmfа

 pic_6.wmfб

Рис. 4. Функции компонентов деформации: а – εx; б – εz

Заключение

В работе определены напряженно-деформируемое состояние анизотропного (трансверсально-изотропного) пласта при фильтрации в нем жидкости. Установлена связь между напряжениями пласта и давлением фильтрирующейся жидкости. Численное решение поставленной задачи получено на основе МКЭ с применением изопараметрического элемента первого порядка. Представлена конечно-элементная модель напряженного деформируемого состояния пласта с горизонтальной скважиной произвольного профиля. Проанализированы пути определения упругих и фильтрационных характеристик деформируемых неоднородных пластов.

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

Рецензенты:

Муратов А.С., д.т.н., профессор, заведующий кафедрой «Математика» Международного казахско-турецкого университета имени Х.А. Ясави Министерства образования и науки Республики Казахстан, г. Туркестан;

Ашибаев Н.К., д.ф.-м.н., профессор, заведующий кафедрой «Теория и методика преподавания математики» Южно Казахстанского государственного университета имени М. Ауезова Министерства образования и науки Республики Казахстан, г. Шымкент.

Работа поступила в редакцию 04.04.2014.