При разработке, мониторинге и планировании эксплуатации месторождения возникает задача определения пространственного распределения коэффициента фильтрационного сопротивления (пьезопроводности) движению флюида, характеризующего пропускную способность продуктивного пласта [2]. Ее решение обеспечивает технология веерной гидродинамической томографии, которая состоит в обработке и анализе веерных измерений интервалов времени наступления реакции в скважинах-приемниках на изменение режима давления в скважинах-источниках [4]. Практическая реализация технологии гидродинамической томографии может быть основана на косвенных измерениях интервальных времен распространения характерных точек кривых восстановления давления по анализу истории разработки месторождения. Прогноз томографических данных реализуется вычислительным экспериментом над построенной математической моделью месторождения в рамках гипотезы о характере ее основных компонентов [3].
В данной работе представлен алгоритм решения прямой задачи гидродинамической томографии поиска кратчайшего пути между парами скважин модели месторождения на основе принципа Беллмана [1, 5]. Что обеспечивает решение обратной задачи, основанной на итеративной схеме последовательного уточнения [3].
Рассматривается область S, в пределах которой выполняются построения. Координата точки из этой области обозначается ξ. Область S покрыта сеткой, которой соответствуют пары индексов так, что каждая точка области в сеточном представлении однозначно представима своими индексами
. В области S задано стартовое распределение коэффициента пьезопроводности
. На сетке (рис. 1) расположены пары скважин
, где m – номер скважины-источника возмущения и n – номер скважины-приемника возмущения. Всего пар скважин Q = M*N.
Рис. 1. Исходные данные для расчета оптимального пути
Распространение возмущения для каждой пары происходит по траекториям , имеющим начало в точке
и заканчивающимся в точке
. Траектории нумерованы индексом
, время движения возмущения по Lq есть q. Количество пройденных точек сетки в каждой траектории есть
. Текущая координата вдоль траектории Lq характеризуется переменной
. Время распространения движения волны q по Lq рассчитывается по формуле в дискретной форме:
, (1)
где l(qk) – длина части траектории Lq до ячейки qk. Наличие этого члена проявляется в снижении скорости по мере увеличения длины пути:
. (2)
Траектория Lq, соединяющая точки ξm и ξn есть линия, служащая решением задачи:
, (3)
где G(ξm, ξn) – возможные траектории между парами скважин q(m, n).
В результате более длинный путь, обходящий зоны повышенных значений κ(qk), за счет снижения скорости по гиперболическому закону может оказаться проигрышным в сравнении с более коротким, проходящим через локальные максимумы κ(qk). Это специфика динамики распространения особой точки коэффициента пьезопроводности.
Алгоритм построения траектории Lq(m, n), вдоль которой будут определены значения l(qk), основан на последовательном нахождении окончания траектории начиная с текущего индекса zkq для траектории q до qn вместе с «остатком траектории» в предположении, что
, и, следовательно, l(qk) заданы своими нулевыми приближениями:
.
, (4)
где gq(qk) – подобранные длины второй половины траектории для .
Алгоритм расчета основан на последовательном решении и сравнении результатов решения задачи, начиная от до
. При этом для каждого qk сохраняются для последующего сравнения все значения конца траектории gq(zkq), служащие именем завершающего участка
.
После завершения счета найденная траектория Gq(qk), характеризующаяся длинами gq(qk), принимается за соответственно.
Рассмотрим скважины с координатами (im, jm) и (in, jn). Согласно принципу Беллмана, движение при расчете траекторий начинается с принимающей скважины, т.е. с точки ξn = ξ(in, jn) и заканчивается в ξm = ξ(im, jm)
Шаг 1 (рис. 2, А). На первом шаге zkq = qn – 1 ищем J путей из ξn = ξ(in, jn) в ячейки
и столько же путей
из (im, jm) в (in – 1, j1).
Получаем пути для первого шага:
(5)
и соответствующие им длины:
. (6)
Шаг 2 (рис. 2, Б). На этом шаге zkq = qn – 2. Далее из ξn = ξ(in, jn) ищем пути в ячейки (in – 2, j2), через ячейки, для которых были найдены пути на шаге 1
,
. И столько же путей из (im, jm) в (in – 2, j2). На рис. 2, Б показан путь для одной из J ячеек.
А Б
Рис. 2. А – Шаг 1; Б – Шаг 2
Шаг 3 (рис. 3, А). Получаем оптимальные траектории и длины
, рассчитанные согласно условию, представленному в формуле (4). На рис. 3, А показан найденный путь, проходящий через ячейку (in – 2, j2), из ячейки ξn = ξ(in, jn) в ячейку (im, jm). Пути
запоминаем, далее работаем с ними.
Шаг 4. Если in – 2 im и j2 / jm, перейти к шагу 5. В противном случае перейти к шагу 8.
Шаг 5 (рис. 3, Б). Допустим, мы уже на итерации zkq = qn – z. Далее из n = (in, jn) ищем пути в ячейки
, через ячейки, для которых были найдены пути на z – 1 шаге
. И столько же путей из (im, jm) в (in – z, jz).
Шаг 6 (рис. 3, Б). C учетом формулы (4), как условия оптимальности пути, находим и длины
. На рис. 3, Б показан путь для одной из J ячеек при z = 3.
Шаг 7. Если in – z / im и jz / jm, перейти к шагу 5. В противном случае перейти к шагу 8.
Шаг 8 (рис. 3, В). На последней итерации zkq = qm + 1 рассчитываем J значений . Далее находим минимальное по jz значение интервального времени распространения возмущения между скважинами, расположенными в точках n = (in, jn) и m = (im, jm), соответствующее значению индекса jz.
А Б
В
Рис. 3. А – Шаг 3; Б – Шаг 5 для z – 3; В – Шаг 8
Оптимальная траектория, удовлетворяющая (4) и соединяющая точки (in, jn) и (im, jm), есть
. (7)
Программная реализация данного алгоритма была выполнена на языке Си#.
Тестовый пример. Рассмотрим модель, состоящую из 6 скважин (30 пар скважин), расположенных на однородной сети с двумя зонами повышенного фильтрационного сопротивления. На рис. 4 показаны кратчайшие пути из одной скважины к остальным, а также соответствующие этим путям интервальные времена.
А Б
В Г
Д Е
Рис. 4. Результаты работы алгоритма: А – скважина-источник № 1; Б – скважина-источник № 2; В – скважина-источник № 3, Г – скважина-источник № 4, Д – скважина-источник № 5, Е – скважина-источник № 6
Получены ожидаемые результаты (пути и времена), что свидетельствует о корректной и эффективной работе алгоритма.
Таким образом, разработанный алгоритм и программа могут быть использованы для получения интервальных времен, обеспечивающих алгоритм технологии гидродинамической томографии. Синтезированные интервальные времена сопоставляются с реальными, а их разность составляет входные данные в итерационном алгоритме восстановления пространственного распределения фильтрационного сопротивления. Эта технология позволяет контролировать возникающие в процессе разработки месторождения зоны потери проницаемости пласта, ведущие к нарушению штатного режима работы месторождения. В результате локализации аномальных зон будут существенно снижены экономические затраты их ликвидации.
Библиографическая ссылка
Кобрунов А.И., Мотрюк Е.Н., Кунцев В.Е. АЛГОРИТМ ПОИСКА КРАТЧАЙШЕГО ПУТИ И ИНТЕРВАЛЬНОГО ВРЕМЕНИ МЕЖДУ СКВАЖИНАМИ НА ОСНОВЕ ПРИНЦИПА БЕЛЛМАНА // Современные наукоемкие технологии. – 2016. – № 8-1. – С. 51-55;URL: https://top-technologies.ru/ru/article/view?id=36101 (дата обращения: 17.08.2024).