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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 144

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Подставляя (3)—(7) в уравнение (2),

имеем:

 

 

 

д Г kh

_др_ "I ,

_д_ Г kh_

др ~|

C$Cr RiN, dp

(8)

дх \ \i

дх J

ду

L Ц

ду J

С{М2 ^

dt

 

Сопоставляя уравнения (1)

и

(8), получаем,

что подобие электрических

и фильтрационных процессов имеет место,

если

 

 

 

 

 

CtM*

 

 

 

(9)

 

 

 

 

 

 

Уравнение (9) — первое условие подобия протекания электрических и филь­

трационных процессов.

 

 

 

 

через некоторое сопротивле­

Сила электрического тока г, протекающего

ние R, согласно закону Ома, равна

 

 

 

 

 

 

 

и 1

ц 2

 

 

(10)

 

 

 

R

 

 

 

 

 

 

 

 

 

 

Здесь их и и2 — напряжения

на концах элемента

электрической сетки

с сопротивлением R.

Если в пористой среде выделить элементарный объем hAx Ау, то расход

жидкости через грань ft Ду в направлении оси х составит

к_

Pi Р2

ft Ду,

У= р

Ах

 

или при равенстве Ах = Ау

 

 

 

Г

(И)

 

 

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

 

 

i = Cgq.

 

( 12)

Подставляя (3), (5) и (12) в уравнение (10), получаем

 

 

 

kh

(Pl — Pi)-

(13)

 

CRCq V-

 

 

 

Необходимость

соблюдения

подобия

уравнений (11)

и (13) приводит

к условию

 

 

 

 

 

 

=

1.

(14)

 

 

CRCq

 

 

Соотношение (14) — второе

условие подобия протекания

фильтрационных

и электрических процессов в рассматриваемой сетке RC.

 

Таким образом,

электрическое моделирование процессов фильтрации упру­

гой жидкости в упругой пористой среде возможно при введении коэффициентов подобия согласно равенствам (3)— (7) и (12) и соблюдении условий подобия (9)

и(14).

Вусловия подобия (9) и (14) входят шесть коэффициентов пропорциональ­ ности. Это позволяет четыре коэффициента подобия выбирать произвольно, например, с учетом конструктивных особенностей модели и исходных данных,

адва других коэффициента вычислять с использованием условий (9) и (14). Известно, что, например, набор (комплект) номиналов сопротивлений и емкостей

вэлектрических моделях довольно ограничен. Поэтому коэффициенты CR и Cg

приходится выбирать исходя из наличных номиналов электрических сопроти­ влений, емкостей и значения параметров p/ftft и p*ft пласта.

128


При помощи метода электрогидродинамической аналогии можно решать задачи подземной газогидродинамики при сложных началь­ ном и граничных условиях. Решение задач на электрических моделях осуществляется достаточно просто и с большой наглядностью. При использовании электроинтеграторов можно проследить шаг за шагом весь процесс разработки месторождения. Не представляет труда в процессе расчетов на любом временном слое вносить изменения, коррективы, новые данные и т. д. К недостаткам метода электро­ аналогий можно отнести некоторую погрешность в получаемых ре­ шениях (в основном из-за невозможности точной аппроксимации карт \ i l k h и неточностей замеров и задания начальных и гра­

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

Перейдем к рассмотрению вопроса об' использовании метода электрогидродинамической аналогии в расчетах по нестационарной фильтрации идеального и реального газов в неоднородных пластах при произвольном расположении скважин.

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

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

д

Г к (х,

у) h (х,

у) _ дрЪ ~1

,__

д_ Г к (х ,

у) h (х, у) _

др*

~| _

 

дх\_

р

дх Jт ” ду

L

р

 

ду

J

 

 

 

 

к (х,

у) т (х ,

y)h (х,у)

др2

 

 

'

 

 

 

 

 

 

р

 

 

dt

 

 

'

при следующих условиях:

 

 

 

 

 

 

 

 

 

 

t =

0, р = р ( х ,

у) =

р я =

const,

(ж,

y ) £ G ;

(1 6 )

 

Я1 =

к (х, у) h (х, у)

 

ds = const;

(ж, у) 6 s,.;

(1 7 )

 

2|РРа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

2,

. . ., п;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 8 )

Здесь

G — область

газоносности;

— нормаль

к контуру

£-й

скважины s,-; — нормаль к внешней границе газоносного пласта Г ; п — число газовых скважин; условие (х , у) 6G означает, что ж и у принадлежат области G.

Начальное условие (16) означает, что в начальный момент вре­ мени (t = 0) пласт находился в невозмущенном состоянии, т. е. давление в каждой точке пласта равнялось начальному давлению р а.

9 Заказ 1013

129


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

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

ОВ,17 12,Зк 18,51 г щ 30,8537,02 $3,19 Ц 3655,53 51,70

£, Ю6с

Рис. 41. Зависимость т = т (t)

Впервые процесс неустановившейся фильтрации идеального газа удалось смоделировать П. М . Белашу на электроинтеграторе при линеаризации уравнения Л . С. Лейбензона на каждом временном слое. Однако при этом требуется значительный объем вычислитель­ ной работы.

Введение, согласно Л. С. Лейбензону [42], новой временной пере­ менной т

 

t

 

t = j ' p d t - { - C

или x — ^ p d t

(19)

 

о

 

позволило П. М. Белашу, А .

Л. Гофлину, А . М.

Кирилечеву,

Н . Г. Степанову, Ю. В. Фрумсону, Э. Чен-Син, И. В. Шипкову на разной основе и почти одновременно свести нелинейное дифферен­

циальное уравнение неустановившейся

фильтрации

идеального

газа (15) к линейному (линеаризованному) уравнению

 

д Г к (х, у) h (ж, у)

др* -] . д Г к (X, у) h (х, у)

др2 ~|

0х\_

р

дх J * ду L

р

‘ ду

J

 

= с ф , у ) m (х, у) h (х,

др%

 

( )

 

 

 

У) дх ’

 

20

 

 

 

 

 

130



Физический смысл введения новой временной переменной, со­ гласно (19), означает введение нового масштаба времени. В фор­ муле (19) давление р принимается независимым от координат и пред­

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

нием среднего пластового давления (при помощи уравнения мате­ риального баланса), не представляет труда путем численного или гра­ фического интегрирования найти зависимость т = т (t ). Эта зависи­

мость имеет вид, приведенный на рис. 41. Убедимся, что это так. Пусть отбор из месторождения постоянен во времени. Тогда

Рат(?£

Р (*) = Рн

а£2н

и для т имеем

t

QPat

dt p Ht

2aQH

В результате решения на электроинтеграторе уравнения (20) получаются показатели разработки газового месторождения в функ­ ции новой временной переменной т. Использование зависимости т == = т (t) позволяет вновь перейти к реальному времени t. Если усло­ вия по скважинам принимаются зависимыми от времени t, то при

решении задачи эти условия задаются в функции временной пере­

менной т. Для перехода от времени U r

также используется зависи­

мость т = % (t).

 

 

Рассуждая аналогично предыдущему случаю, введем коэффи­

циенты пропорциональности согласно

следующим соотношениям:

И =

Срр 2;

(2 1 )

t3 =

Cfi\

(2 2 )

 

 

(23)

г э

У э

(24)

 

С ~

Cmam h,

(25)

При подстановке выражений (21)— (25) в уравнение (2) получаем, что для достижения подобия процессов фильтрации газа в пористой среде и электрических процессов в сетках из сопротивлений и емко­ стей должно выполняться условие

с « Р п

1.

(26)

CtM*

 

 

Расход газа q вдоль оси х через элементарный объем h dx dy,

приведенный к атмосферному давлению и пластовой температуре, будет27

к

dp

Лр2

(27)

Я Р ат р

dx

2рат

9'

131