ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 88
Скачиваний: 0
Низкая и неравномерная фазовая проницаемость и упругость пород выбросоопасиых зон резко снижает эффективность исполь зуемых методов.
Более эффективными методами борьбы с внезапными выбросами будут: 1) предварительный прогрев выбросоопасных зон до темпера туры 55—60° С; 2) закачка в выбросоопасные зоны через опережа ющие скважины веществ, интенсифицирующих процессы выделения из пор связанного газа. В качестве таких веществ могут быть спирты
иэлектролиты; 3) закачка веществ, хорошо растворяющих метан;
4)закачка маловязких веществ в сопредельную зону залегания выбросоопасных линз с последующей их полимеризацией с целью
упрочнения (цементации) пород.
Для более рентабельного ведения предупреждения внезапных выбросов необходимо определять с максимальным приближением положение и размер линз выбросоопасных пород.-
В качестве одного из эффективных методов обнаружения гидрат содержащих выбросоопасных линз может быть использован сейсми ческий метод.
Глава V
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ РАЗРАБОТКИ ГАЗОГИДРАТНЫХ ЗАЛЕЖЕЙ
Открытий |
[12, 32, 39] |
газогидратных залежей, приурочен |
ных к зонам |
пониженных |
температур и высоких давлений, выдви |
нуло необходимость решения задач их разработки. Газогидратные залежи распространены на половине территории суши нашей страны и на 25% территории суши мира. Они могут встречаться практиче ски во всех придонных осадочных отложениях мирового океана. Запасы газа в таких залежах, вероятно, превышают запасы всех обычных газовых месторождений, встречающихся иа суше нашей планеты.
Разработка газогидратных залежей основана на одном общем принципе — газ в таких залежах должен быть переведен из твердого связанного состояния в свободное непосредственно в пласте.
В пластовых условиях перевод газа из твердого в свободное со стояние может быть осуществлен путем снижения давления ниже давления разложения гидрата при пластовой температуре; путем повышения температуры выше температуры разложения гидрата; путем ввода в пласт веществ, разлагающих гидрат при неизменных пластовых давлении и температуре.
В приведенной главе исследуются вопросы фильтрации газа при разложении гидрата в пористой среде. При решении поставлен ных задач используется аналогия рассматриваемого явления с тепло проводностью при изменении агрегатного состояния тел, известной в математической физике как задачи Стефана.
Во избежание усложнения постановки задач и математических выкладок в расчетах в первом варианте не учитываются выделение воды при разложении гидрата, изменение удельного объема при изменении растворимости газа в воде.
§ 1. Разложение гидратов газа понижением давления в залежіг
При разработке газогидратных залежей путем перевода ее в газо образное состояние в пласте образуется область разложения гидрата и движение свободного газа. В области, где еще гидрат не разложен, газ будет находиться в твердом состоянии.
Таким образом, в пласте можно выделить две области с разными коллекторскими свойствами.
172
Линейное течение
Рассмотрим условия перемещения границы разложения гидрата в пласте
&і д Р і _ |
а'2 д р 2 |
Н7ДгДо |
.ѵ |
,. |
дх |
дх |
рр ( Т ) dt ’ |
\ |
) |
где W — содержание газа в гидратпом состоянии в единице объема пористой среды; рр(Т) — давление разложения гидрата при пласто вой температуре Т\ к г и к 2 — коэффициенты проницаемости соот ветственно в областях свободного газа и гидрата; ррр2 — текущие давления; р 0 — атмосферное давление; рг — вязкость газа; | — координата границы разложения гидрата; х — текущие координаты;
t — время.
Изменения давления в залежи определяются системой дифферен циальных уравнений:
д-р\ |
_ 2/иіЦг |
дРі |
. |
/ Ѵ |
|
дхі |
к і |
dt |
’ |
<,v - 7 |
|
д*р\ |
2 л г 2 Ц г др2 |
|
(Ѵ.З) |
||
dxt |
ка |
dl |
’ |
||
|
здесь m lt m2 — коэффициенты динамической пористости в областях 1 и 2 соответственно.
Уравнения (Ѵ-2) и (Ѵ.З) можно линеаризовать по Л. С. Лейбен-- зону или И. А. Парному. Полученные после линеаризации уравне ния типа уравнения теплопроводности позволяют сформулировать проблему как задачу Стефана, которая имеет достаточное число решений. В частности, для постоянного давления на галерее реше
ние уравнения (Ѵ.1) |
будет: |
|
а |
|
||
|
|
|
|
|
||
M pJ - p?) |
|
|
|
4а» |
= 4 ^ г Р к /«■ |
|
V |
|
|
|
|||
/;і1Ф |
ttl‘2 Зіво |
1 - Ф - V |
||||
— |
||||||
|
||||||
|
|
|
2 |
V |
а2 |
Представленное уравнение является трансцендентным относитель но а и неудобным для практических расчетов. Поэтому остановимся только на методе Л. С. Лейбензона, позволяющем получить доста точно простое решение. Погрешность этого метода оценена в задачах Стефана и не превышает 6 %.
В области 1 распределение давления принимаем как в стацио
нарном потоке: |
|
|
|
• |
Л = ]/р г + (Pp — P r )j, |
(V.4) |
|
где р — постоянное |
давление |
на галерее. |
|
В области 2 примем закон |
Гаусса 1: |
|
|
Р2 = |
/■ |
A-sPk |
|
Рк + (Рк—P*)erf JünJL- |
|
||
|
|
2 ѴаЦі ’ |
|
|
City |
то[Ar |
(V.5) |
|
|
|
173
Решение уравнения (V.1) ищем в |
виде: |
§ = а]/7 . |
(V.6) |
Подставим выражения (V.4), (V.5) и (V.6) в уравнение (V.1).
После этого найдем корень |
полученного |
квадратного уравнения |
|
|
■(РІ-Рр |
ІІРІ-РІ? |
-/мИ'ЦгРо (Р \-Р І) |
а ■- |
V^rta2 |
PFprPo |
(V.7) |
|
|||
|
|
|
Отрицательный знак перед корнем не имеет физического смысла, поскольку а величина всегда положительная.
Из уравнений (V.4) и (V.5) определяем дебит галереи:
Вкф (р° —рі)
(V.8)
2 р гР о а V t
В частном случае, когда а = 0, из формулы (V.7) получается решение для газированной залежи, поровое пространство которой полностью заполнено гидратами.
Радиальное течение
Рассмотрим задачу при работе скважины в центре плоско-ради- альпого пласта.
В данном случае распределение давления описывается системой дифференциальных уравнений:
д*Рі |
|
, 1 |
|
2mi\lr |
дрг , |
dr2 |
|
1 г dr |
Л 1 ; 1 , |
dt ’ |
|
d'iPZ |
_| |
1 |
dpi |
2 m i p r |
dp2 |
dr°- |
' |
r |
dr |
/г 2 |
dt ’ |
(V.9)
(Ѵ.10)
здесь г — текущий радиус.
При постоянном отборе газа из центральной скважины точное решение уравнений (Ѵ-9) и (Ѵ.10), линеаризованных по Л. С. Лей-
бензону, приводит к соотношениям: |
|
|
|
|
а |
|
|
7РРат е 4а, |
к2(РІ-Рр е 4a= |
mp.raW |
R = y at. |
|
2Pp ; |
||
АШірр |
|
|
Для приближенных расчетов в области 1 закон распределения давления опишем как в стационарном потоке
Рі = |
(Ѵ.11) |
174
В области 2 примем автомодельное решение дифференциального уравнения (V.10):
„2 |
(Ѵ.12) |
р ; - ' 1 Е 1 ( - і £ г ) + в - |
|
Начальные и граничные условия примем следующие:
Рч (°°, t ) = p 2{r, 0) = р к;
ра (Л, ü)=lV
Из условии (Ѵ.13) получим
В = рІ.
Удовлетворяя решение (Ѵ.13) условию (Ѵ.14), найдем
А-. РІ-РІ
“ ( - т а - ) '
Подставим искомые постоянные коэффициенты в (Ѵ.12):
Е \ |
4*2t ) |
(Ѵ.13) (Ѵ.14)
(V.15)
(Ѵ.16)
(Ѵ.17)
На подвижной границе разложения гидрата соблюдается условие
д. |
dpi (Д ■0 ^ |
(Д. О |
^РгРо |
ад |
ГѴ 18> |
||
1 |
д/- |
2 |
<?г |
|
Р р ( Т ) |
d t |
\ ■ / |
§ 2. Термический |
метод |
разложения |
гидратов |
|
Плоско-параллельный пласт.
Гидраты |
можно разложить также и при повышении температуры. |
|
Рассмотрим |
задачу, когда тепло вводится в пласт через |
галерею. |
II в данном случае пласт можно разбить на две области. |
В первой |
области тепло распределяется в среде разложившегося гидрата. Во второй области происходит нагревание гидрата до температуры разложения. В такой постановке вопрос можно свести к классиче ской задаче Стефана, решения которой имеются как в точной поста новке, так и полученные приближенными методами. Здесь мы рас смотрим комбинированный метод, являющийся, как нам кажется, весьма простым и достаточно точным в инженерных задачах.
При точном решении задачи Стефана распределение' температуры в областях (Ѵ.19) и (Ѵ.20) соответственно описываются уравнениями
Т 1= Т Г |
|
ф |
|
OsSz s SZ-jl; |
(Ѵ.19) |
|
|
|
X |
а _ \ |
Li ^ X оо, |
|
|
Т2 = ------- |
ф 2 V a2t — Ф |
(Ѵ.20) |
||||
2 Ѵаг )\ ’ |
||||||
1 -ф |
( |
|
|
|
|
175
здесь |
Тр |
— температура на открытом конце полуограниченного |
||
параллелепипеда; Т 0 — начальная |
температура |
тела; а — коэф |
||
фициент |
температуропроводности; |
Ф — функция |
Крампа; а — по |
|
стоянный |
коэффициент. |
|
|
|
Длина распространения нового агрегатного состояния опреде |
||||
ляется |
уравнением: |
|
|
|
|
|
£ = а ]/7. |
(V.21) |
Рассмотрим распределение температуры в области 1 с новым агрегатным состоянием. Введем в функцию Крампа верхнее значе ние координаты х\
Ф |
а |
(Ѵ.22) |
|
2 КД |
|||
|
|
Так как практически х является величиной очень малой, то
ф Н М |
|
= |
- 7 |
^ |
= ^ ^ . |
(Ѵ.23) |
\ 2 Ѵах J |
|
Упах |
У7iaxt |
|
||
Учитывая, что х ^ Lx, |
можно |
принять |
|
|||
ФІ 2 Уaxt j |
|
Утіа-it |
(Ѵ.24) |
Подставив (V.24) в (V.19), получим:
Тг = Тг |
Г г |
X |
(Ѵ.25) |
|
Ф |
V 7ia.it |
|||
|
|
2 Ѵах
Определим из (Ѵ.25) температуру в точке плавления:
(Ѵ.26)
Решим уравнения (Ѵ.25) и (Ѵ.26) совместно:
Тх = Тт- ( Т т- Т р) ^ . |
(Ѵ.27) |
Формула (V.27) показывает, что в области 1 закон распростра нения температуры очень близок к линейному. Поэтому в прибли женных расчетах для области 1 примем линейный закон распреде ления температуры так же, как и при установившемся режиме.
В области 2 координата х сверху не ограничена. Поэтому закон распределения температуры здесь нельзя аппроксимировать прямой линией. Следовательно, в области 2 необходимо пользоваться более точными решениями.
176