Файл: Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие.pdf

Установим подобие между электрическими и фильтрационными параметрами при помощи следующих соотношений:

(индекс «э» относится к электрическим параметрам).

Тогда условия подобия протекания электрических и фильтра­ ционных процессов запишутся в виде:

Излагаемый подход основывается на решении исходной задачи последовательно шаг за шагом по времени. Перед отысканием реше­ ния задачи в момент времени t + A t предварительно пересчиты­

ваются электрические емкости по формуле

Здесь за р принимаются значения давления в узловых точках, найденных на момент времени t (или соответственно t3).

Затем производится перенабор электрических емкостей. Для отыскания решения задачи на момент времени t -f- A t в качестве

начального условия задается распределение напряжений на сетке, полученное на момент времени t. Точность решения задачи возра­ стает при уточнении емкостей С согласно (8) в итерационном цикле

на каждом временном слое.

При решении задач неустановившейся фильтрации газа по шагам могут быть применены и иные подходы.

Запишем уравнение (1) в виде:

Введем коэффициенты пропорциональности при помощи следу­ ющих соотношений:

Если на каждом временном слое производить перенабор электри­ ческих сопротивлений, то можно получить подобие фильтрационных и электрических процессов на сетках R C . Тогда условия подобия

протекания электрических и фильтрационных процессов записы­ ваются в виде:

146

Таким образом, по известному решению задачи на момент вре­ мени t предлагается рассчитывать электрические сопротивления

по формуле

и перенабирать их для отыскания решения задачи (1)— (4) на момент времени £ + Af. Это означает, что линеаризация уравнения (1) осуществляется на каждом временном слое. При достаточно малом шаге по времени A t данная линеаризация обеспечивает большую

точность, чем методы, рассмотренные в § 7. При необходимости (при достаточно большом шаге Д t) на каждом временном слое реа­

лизуется итерационный процесс.

После перенабора сетки сопротивлений в качестве начального условия задается поле напряжений, полученное для момента вре­ мени t, и отыскивается решение задачи на момент времени t + A t .

Итак, решение на электрических моделях задач неустановившейся фильтрации идеального газа возможно осуществлять по шагам при введении коэффициентов пропорциональности согласно (10), удо­ влетворении условиям подобия (11) и (12), пересчете и перенаборе на каждом временном шаге электрических сопротивлений.

Дифференциальное уравнение неустановившейся фильтрации реального газа в реальной, неоднородной по коллекторским свой­

Краевые условия относительно функции ф переписываются ана­ логично тому, как показано в § 7.

Уравнение (16) в отличие от (14) и в левой и в правой частях за­ писано относительно одной неизвестной функции — ср. Кроме того, нелинейность уравнения (16) определяется лишь членом перед производной по времени.

Нетрудно видеть, что если при отыскании решения на момент времени t + A t нелинейный член в уравнении (16) вычислять по из­ вестному решению на момент времени t, то решение уравнения (16)

при соответствующих краевых условиях может быть получено на сетках R C . В этом случае значение емкости в каждой узловой точке

будет вычисляться (а затем уточняться в итерационном цикле) по формуле

Связь между электрическим и натурным временами задается соотношенпем t3 — Ct t. Остальные коэффициенты пропорциональ­

ности вводятся так же, как в § 7, при соблюдении соответствующих условий подобия.

Перед отысканием решения задачи на момент времени t + A t

в качестве начального задается поле напряжений, соответствующее моменту времени t.

Отметим, что при решении уравнения (16) может быть использован также подход, основанный на методе корректирующих токов (Н. Г. Степанов).

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

На электроинтеграторе УСМ-1 была решена задача о неустановившемся притоке идеального газа к галерее при безразмерном дебите д* = 0,5. Для со­ поставления использовано практически точное решение соответствующей задачи, полученное на ЭВМ [39].

Результаты расчетов показывают, что линеаризация дифференциального уравнения неустановившейся фильтрации газа с использованием среднего пластового давления при больших процентах отбора от запасов газа (или при значительной «глубине» общей депресспонной «воронки») приводит к погреш­ ностям в десятки процентов и более. Пошаговое решение соответствующей задачи дает погрешность около 5—6% (к концу разработки газовой залежи).

§ 10. Численные методы определения показателей разработки газовой залежи

при неравномерном расположении скважин

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

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

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

148

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

Применительно к теории и практике разработки газовых (и нефтяных) месторождений достаточно рассмотреть численные методы решения двумерных (по х и у) уравнений параболического типа. Строго говоря, процессы фильтра­ ции, происходящие в пласте при разработке газового месторождения, протекают во времени в трехмерном эвклидовом пространстве. Однако для рассмотрения трехмерных фильтрационных потоков требуется огромная геолого-физическая информация. Получение по ограниченному числу скважин исчерпывающей и достоверной исходной геолого-физической информации для решения трех­ мерных фильтрационных задач представляет сложную проблему. Этим в извест­ ной мере определяется рассмотрение здесь лишь двух методов (Д. Дугласа и А. А. Самарского) численного интегрирования двумерных уравнений пара­ болического типа. Этп методы применительно к фильтрационным задачам апро­ бированы в работах кафедры и проблемной лаборатории по газу МИНХ и ГП им. И. М. Губкина.

Идея методов Дугласа и Самарского наиболее наглядно иллюстрируется

Здесь / —■плотность источника (стока).

Рассмотрпм пока алгоритм решения уравнения (1) в квадратной области. На формулировке и записи начального и граничных условии останавливаться

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

Сущность методов численного интегрирования двумерных уравнений пара­ болического типа состоит в таком «расщеплении» исходного уравнения, что решение задачи получается в результате последовательного решения одномер­ ных разностных задач.

Согласно методу Д . Дугласа, при известном решении задачи на к-м времен­ ном слое решение на (к + 1)-м временном слое получается в результате после­ довательного решения следующих двух систем копечпо-разностпых уравнений:

сматриваемой квадратной области интегрирования; при к = 0 используется

начальное условие.

Системы уравнений (2) или (3) представляют собой системы из (п — I)2 алгебраических уравнений с (re — I )2 неизвестными. С учетом граничных усло­ вий системы (2) и (3) в общем случае будут системамп из (п + I)2 алгебраических уравнений с (п + I)2 неизвестными. Характерным для систем (2) и (3) является то, что они имеют трехдиагональную матрицу (на каждой строке или столбце сеточной области). Это обстоятельство позволяет при нахождении решения на

149

промежуточном временном шаге и к + 1-м временных шагах использовать метод прогонки.

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

на к-м временном слое. Получается промежуточное решение задачи на & + — -м

временном слое. Так же решается система уравнений (3) (с учетом граничных условий). Получаемое решение является искомым и соответствует к + 1-му моменту времени.

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

По методу А. А. Самарского уравнение (1) аппроксимируется следующими

Каждая из систем уравнений (4) и (5) вместе с граничными условиями содер­ жит в общем случае по (п + I)2 алгебраическому уравнению с (п 4- I)2 неиз­ вестными. Системы (4) и (5) имеют трехдиагональную матрицу (на каждой строке или столбце сеточной области). Поэтому, как и в предыдущем случае, для их решения можно использовать эффективный метод — прогонку.

Последовательность решения систем уравнений (4) и (5) ничем не отли­ чается от порядка решения систем (2) и (3). Обращаем лишь внимание на то,

Конечно-разностные уравнения (2) и (3); (4) и (5) являются неявными, абсолютно-устойчивым. Они аппроксимируют исходное уравнение с погреш­

то, согласно А. А. Самарскому, уравнение (6) может быть аппроксимировано следующими системами конечно-разностных уравнений:

150