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

SIMULATION FEATURES OF FRACTURE-POROUS COLLECTORS IN THE FUNDAMENTAL PROBLEMS OF FLUID MECHANICS IN COMPLEX SYSTEMS

Grachev S.I. 1 Strekalov A.V. 1 Sarancha A.V. 1
1 Federal State Budget Educational Institution of Higher Education «Tyumen State Oil and Gas University»
Recently become topical methods of most wave impact on productive formations such as non-stationary injection of water, intermittent drainage, infrasound and hydro-pulse impact on the bottomhole formation zone. The effectiveness of wave processes caused by a number of factors: the compressibility of the skeleton, the formation of cracks and in General an increased level of energy impact on filtration processes. Due to the fact that the presence of a network of cracks in productive permeable reservoir causes significant structural heterogeneity, for the full simulation requires the use of mathematical apparatus of the theory of hydraulic circuits. Therefore, the aim of this article is to review fundamental issues in modeling hydraulic systems with a highly developed arbitrary structure, which in particular include a system of cracks in fractured porous reservoirs. Physical and computational experiments allowed to identify «weirdness» in the behavior of some hydraulic systems, which is in contradiction to the generalization of Maxwell is not in a stationary mode in the absence of external influence.
modeling
hydrovane
model
mathematical function
performance
software and computing complex
1. Gracheva S.K., Strekalov A.V., Husainov A.T. Modelirovanie obrazovanija seti treshhin / Neftegazovoe delo. 2013. no. 2. рр. 168–183.
2. Merenkov A.P., Hasilev V.Ja. Teorija gidravlicheskih cepej. N., 1985. 276 р.
3. Morozov V.Ju. Tehnologija regulirovanija gidrosistem neftjanyh promyslov na primere sistem podderzhanija plastovogo davlenija: dis. ... kand. tehn. nauk. Tjumen, 2010.
4. Strekalov A.V., Grachev S.I. Programmnyj kompleks gidrodinamicheskogo modelirovanija prirodnyh i tehnicheskih sistem «Nemezida Gidrasim 2014» (Nemesis Hydrasym 2014). Svidetelstvo o gosudarstvennoj registracii programmy JeVM no. 2014614505. Zajavka no. 2014612343. Data gos. registracii 28 aprelja 2014.
5. Strekalov A.V., Morozov V.Ju. Model nestacionarnyh processov v gidravlicheskih sistemah setevoj struktury / Neftjanoe hozjajstvo. 2010. no. 8. рр. 107–109.

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

Несмотря на достаточную изученность данных процессов в отдельности, к сожалению, при моделировании процессов фильтрации в поровых и трещинопоровых средах они в комплексе не нашли отражение в общей системе уравнений, решаемой в классических (принятых в качестве основных при проектировании разработки месторождений углеводородов) симуляторах – Эклипс, Темпест, Техсхема и т.п. Более того, вычислительная схема (полностью или частично неявная) не позволяет прогнозировать негармонические колебания, вызванные историей гидродинамического воздействия на продуктивные пласты и любые гидравлические цепи в принципе.

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

Авторами данной статьи разработан и используется вычислительный программный комплекс (ВПК) Немезида [4] (с. Гидрасим) (рис. 1), апробированный множеством научных работ и подтвержденный многочисленными промысловыми и лабораторными экспериментами [3, 5].

В частности, в 2008 г. была создана стендовая установка для испытания разветвленной гидросистемы в условиях нестационарных волновых процессов и в т.ч. гидроударов. Результаты стендовых испытаний подтверждают корректность разработанной модели [5] посредством сопоставления результатов вычислительных экспериментов на ВПК и физической модели (стендовой установки).

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

Вследствие наибольшей универсальности разработанная модель ВПК позволяет прогнозировать как гармонические, так и негармонические колебания, являющиеся естественными или следствием работы элементов управления, параметрами энергообмена и физическими свойствами текучих сред (рис. 2). Математический аппарат в отношении волновых процессов описан в работе [5].

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

Последние исследования авторов, как наследие более старых исследований [3], направлены на изучение нестационарных процессов в техногенных и природных гидравлических системах.

В ходе стендовых испытаний [3] в 2008 г. было обнаружено «странное свойство» некоторых гидравлических систем: при определенных условиях запуска гидросистема не выходит на стационарный режим.

Вместе с тем известно, что любая система стремится перейти от неравновесного состояния к равновесному, т.е. стационарному или квазистационарному состоянию, в котором во всех узлах системы устанавливаются определенные давления, а в участках системы определенные расходы. На этом принципе, например, основан часто применяемый метод исследования скважин – снятие и обработка «индикаторных линий и диаграмм». На данном же принципе основаны методы решения задачи потокораспределения [2].

pic_22.tif

Рис. 1. Моделирование волновых процессов в вычислительной стендовой установке в ВПК «Немезида»

pic_23.tif

pic_24.tif

Рис. 2. Демонстрация волнового процесса в системе трещин продуктивного пласта ВПК «Немезида»

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

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

Из такой ситуации гипотетически можно было бы предположить невозможное: нарушение закона сохранения энергии. Однако, разумеется, это не так.

Основным признаком такого – особого состояния гидросистемы, когда она не выходит на стационарный режим, является сохранение автоколебаний давления и скоростей в некоторых участках системы со сдвигом фаз и разными амплитудами. Иными словами, когда в одном из участков наблюдается снижение скорости течения и снижение диссипации гидравлической энергии в тепловую, при этом в другом элементе прослеживается рост скорости и перепада давления на вязкостное трение.

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

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

Экспериментально также установлено, что одна и та же гидросистема с различных начальных распределений давления переходит в различное условно стационарное состояние (рис. 3, б и в). То есть уровень энергии одинаков, при этом характеристики автоколебаний могут отличаться.

pic_25.tif

а

pic_26.tif

б

pic_27.tif

в

Рис. 3. Динамика давления узла 4 с различных начальных давлений: а – 11 МПа; б – 5 МПа; в – 20 МПа

Следует отметить, что, рассчитывая любые гидравлические системы, будь то наземные сети с развитой структурой или продуктивные пласты с изотропной структурой или трещиноватой структурой, и предполагая оптимальным и технологически номинальным (или даже возможным) состоянием системы исключительно стационарный режим (как на рис. 3, а), мы находимся в заблуждении. Иначе говоря, волновые процессы в природе гидросистем – есть одно из средств сохранения их состояния с наименьшим уровнем диссипации энергии.

Данный вывод тесно связан и с пониманием стационарного состояния в квантовой физике.

Колебательный волновой процесс формируется так, что потенциальная величина давления стационарного состояния находится ровно между максимумом и минимумом давления с периодическим профилем волны, как при субгармонических (рис. 3, в), так и не гармонических (рис. 3, б) колебаниях.

В таких неклассических условиях «квазистационарных процессов» материальный баланс в узловых соединениях не соблюдается (рис. 4), из-за чего и меняется давление вследствие сжимаемости среды.

pic_28.tif

Рис. 4. Динамика скорости в одном из участков гидросистемы для случая запуска с давления 5 МПа

Как уже отмечалось [4], волновые процессы в продуктивных коллекторах с развитой сетью трещин (рис. 2) требуют дальнейшего изучения и технологического использования с целью повышения нефтеотдачи пластов с низкой поровой проницаемостью.

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