Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
Прогнозирование технологических показателей разработки (нефтеизвлечения) на основе цифровой модели процесса
. Под прогнозированием понимается установление заключения о предстоящем развитии, т. е. предсказание о течении технологического процесса разработки в будущем. Следовательно, к методам прогнозирования относят все методы моделирования процесса разработки. Моделирование, и как результат, модель процесса обеспечивают возможность при сравнительно небольших затратах в короткие сроки многократно (многовариантно) «проиграть» медленно протекающие процессы разработки в различных технологических условиях и тем самым выбрать рациональную технологию. Различают физическое и математическое моделирования. Математическое моделирование заключается в исследовании процессов путем построения и решения системы математических уравнений, относящихся к собственно процессу и краевым условиям. Математическая модель основана на упрощении (идеализации) сложного реального процесса. Системы математических уравнений решают аналоговым и вычислительным методами. Вычислительные методы подразделяются на аналитические, численные(цифровые) и статистические. Численные методы (цифровые) заключаются в определении с помощью ЭВМ численных значений функции в некоторых дискретных точках для заданных численных значений аргумента, т. е. решение получается в некоторых точках пространства. Для этого пространственная область фильтрации мысленно разделяется на ряд квадратов или блоков путем наложения сетки определенного типа (в большинстве равномерной квадратной сетки). Исследуемый интервал времени также разделяется на отдельные элементарные интервалы с постоянным шагом. Преобразование непрерывных дифференциальных уравнений к дискретному виду осуществляется с помощью метода конечных разностей. Получить конечно-разностные уравнения можно методом разложения функции в ряд Тейлора в заданной точке, решая уравнение относительно искомой производной. Разложение функции в ряд Тейлора с использованием разностей прямых (вперед) и обратных (назад) соответственно можно записать: Из этих уравнений вычитанием определяем первую производную: где O( )—погрешность усечения (остаток), связанная с аппроксимацией функции; имеет порядок . Сложив уравнения (2.10) и (2.11), получим вторую производную Таким образом, для дискретной системы точек, пренебрегая погрешностью усечения, имеем: Отсюда понятно, что численные методы всегда приближенные, так как замена производных отношением конечных приращений вносит погрешность. Она тем меньше, чем меньше приращения (шаг). Для перехода к конечно-разностному уравнению обозначим узловые точки вдоль оси х индексомi, вдоль оси у — индексом j, вдоль оси времениt — индексомk. Имеются два основных способа перехода от значений на прежнем уровне времени к значениям на новом уровне: явная схема, когда новые значения функции для каждой точки вычисляются по значениям соседних точек прежнего уровня; неявная схема, когда все неизвестные значения нового уровня определяются одновременно. Для решения двумерных задач применяется неявная схема. Использование ее дает конечно-разностный аналог, например, дифференциального уравнения упругого режима в однородном пласте в виде где - коэффициент пьезопроводности пласта. В данном уравнении пять неизвестных давлений Такие уравнения записываем для каждой узловой точки сеточной области интегрирования (фильтрации) на (k+1)-й момент времени. Получаем вместо дифференциального уравнения систему из N алгебраических уравнений с N неизвестными, решая которую, определяем с помощью ЭВМ искомые давления в каждой узловой точке. Выполняя аналогичные расчеты для других временных уровней, находим изменение давления во времени в каждой узловой точке, Для расчета при k=0 задается начальное условие. При значенияхi и j, соответствующих узлам на внешней границе, используются граничные условия. Внешняя граница аппроксимируется ломаной сеточной границей. Аппроксимировать контур скважины не представляется возможным, так как применяемый шаг сетки (100—2000 м) существенно больше радиуса скважины. Г. Г. Вахитов показал, что в узловых точках расположения скважин вычисляемые давления равны давлениям на забое некоторой фиктивной («точечной») скважины с радиусом rсф=0, 2Δ х (приΔ х = Δ у). Тогда для расчета забойного давления в реальной скважине требуется учесть фильтрационное сопротивление между контурами фиктивной и реальной (с приведенным радиусом) скважин.
|
Последнее изменение этой страницы: 2017-05-05; Просмотров: 761; Нарушение авторского права страницы