1275
.pdfпарциальное давление  | 
	кислорода, Па;  | 
	Е — энергия активации,  | 
||
кДж/( моль-К);  | 
	Rr—  | 
	универсальная  | 
	газовая  | 
	постоянная,  | 
кДж/(моль-К); Т — абсолютная температура, К; К — коэффициент, учитывающий влияние степени окисленности нефти на скорость
реакции; / — степень окисленности нефти, доли единицы; п — по казатель степени при парциальном давлении кислорода, доли еди ницы.
Энергию активации необходимо определять с использованием метода наименьших квадратов. Для этого зависимость (5.11) не обходимо привести к линейному виду. Прологарифмировав правую и левую части (5.11), получим
1пау=1пЛ — В (1/Г),  | 
	(5.12)  | 
где
Л = Л0р£2ехр ( — Kf)\
B = EIRr.
Так как в опытах f = 0 и все величины скорости реакции при ведены к одному значению парциального давления кислорода, то
Л = Л0ро2= const.
Тогда значение коэффициента В вычисляют по формуле
где кЭг/—  | 
	число замеров скорости  | 
	реакции;  | 
	winp — скорость ре  | 
||
акции  | 
	нефти с кислородом воздуха  | 
	в i-u  | 
	опыте, кг/(кг-с); 7\- —  | 
||
температура кернодержателя в i-м опыте,  | 
	К;  | 
||||
£ _  | 
	6-0,351678— 113,74071-0,0184258  | 
	_  | 
	_  | 
	8133  | 
|
_  | 
	6-56,87815-10-°— (18,4258)М0-С  | 
	~  | 
	
  | 
	
  | 
|
Е = B R r = 8133 -8,3736 = 68 102 кДж/моль.
З а д а ч а 5.24. Вычислить энергию активации Е по результа там лабораторных исследований кинетики окисления нефти.
Эксперименты проводились по методике, которая описана в за даче 5.20. Результаты опытов приведены в табл. 51. По результа там предыдущих исследований было определено, что показатель степени при парциальном давлении кислорода в формуле (5.11)
п = 1,0, а коэффициент, учитывающий влияние степени окислен ности нефти на скорость реакции, К = 5.
181
Таблица 51
1
Номер  | 
	Скорость реак  | 
	Температура  | 
	Среднее парциаль  | 
	Степень  | 
экспери  | 
	ции, кг/(кг-с)  | 
	кернодержателя.  | 
	ное давление  | 
	окисленности  | 
мента  | 
	
  | 
	К  | 
	кислорода, 105 Па  | 
	нефти  | 
1  | 
	3,34-10—9  | 
	294,1  | 
	4,98  | 
	0,003  | 
2  | 
	1,41 • 10- 8  | 
	306,8  | 
	5,01  | 
	0,008  | 
3  | 
	3,01-10-8  | 
	318,6  | 
	4,70  | 
	0,013  | 
4  | 
	1,49-10 - 7  | 
	337,1  | 
	4,93  | 
	0,029  | 
5  | 
	3,35-10—7  | 
	352,0  | 
	4,80  | 
	0,056  | 
6  | 
	3,41 • 10—7  | 
	357,3  | 
	4,83  | 
	0,098  | 
При решении задачи необходимо замеренные скорости реакции привести к средним значениям парциального давления кислорода и степени окисленности нефти.
О т в е т :  | 
	£ = 64 510 кДж/(моль-К).  | 
З а д а ч а  | 
	5.25. В нагнетательную скважину с целью прогрева  | 
призабойной зоны пласта при инициировании процесса внутрипла-
стового горения закачивается горячий воздух.  | 
	температура  | 
|||||||
Темп  | 
	нагнетания  | 
	воздуха  | 
	^Воз = 6-104  | 
	м3/сут;  | 
||||
воздуха Т =  | 
	300 °С;  | 
	теплоемкость  | 
	воздуха  | 
	с^з = 1  | 
	кДж/(кг-°С);  | 
|||
плотность  | 
	воздуха  | 
	рв03 = 1,3  | 
	кг/м3; толщина пласта  | 
	h = 10 м;  | 
||||
коэффициент  | 
	охвата  | 
	пласта  | 
	по  | 
	толщине  | 
	т|о = 0,8;  | 
	начальная  | 
||
пластовая температура Т 0 = 20 °С; теплопроводность пласта и ок
ружающих  | 
	его  | 
	пород  | 
	А,пл = А,оп =  | 
	8 Вт/(м-К);  | 
	теплоемкость  | 
|
пласта и  | 
	окружающих  | 
	его  | 
	пород  | 
	спл = соп = 1,1  | 
	кДж/(кг-К);  | 
|
плотность  | 
	пласта  | 
	рпл =  | 
	2600  | 
	кг/м3.  | 
	
  | 
	
  | 
Определить температуру пласта Т на расстоянии 0,5 м от нагне
тательной скважины через 1 день после начала прогрева.  | 
	
  | 
|||||
Р е ш е н и е .  | 
	Температуру пласта можно рассчитать по за  | 
|||||
висимости, полученной А. Б. Золотухиным,  | 
	
  | 
|||||
0 =■  | 
	О*)2  | 
	(5.13)  | 
||||
12v(l — V и)  | 
||||||
  | 
	
  | 
	
  | 
||||
где  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	hr\2  | 
	
  | 
	
  | 
	
  | 
|
v = ? В О З Р в о з ^ в о з  | 
	
  | 
|||||
  | 
	
  | 
	4л/АплП2  | 
	
  | 
|||
и  | 
	=  | 
	Т -  | 
	т0  | 
	
  | 
||
т Б О З  | 
	
  | 
	T Q  | 
	
  | 
|||
  | 
	
  | 
	------  | 
	
  | 
|||
0  | 
	
  | 
	^-ПЛt  | 
	.  | 
	
  | 
||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||
СплРпл^ ^2
г — расстояние от оси нагнетательной скважины, м; t — время нагнетания агента, с.
182
Для условий задачи преобразуем формулу (5.13) к виду
Г - Г „ + ( Г „ —  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
Подставив  | 
	исходные данные, получим  | 
	
  | 
	
  | 
	
  | 
|||||
г* = 0,5/(10-0,8) = 0,0625;  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
v =  | 
	0,07-1,3-1  | 
	= 0,1132;  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
  | 
	4-3,14-Ю -0,88 -10—3  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	8 -10~3-8,64-104  | 
	0,003776.  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	1,1-2600- 102 (0,8)2  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
Тогда  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||
Т = 2 0 + (3000- 20) [ 1-  | 
	(0,0625)2  | 
	
  | 
	88  | 
	°С.  | 
|||||
12-0,1132-3,776-10= rj  | 
|||||||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
З а д а ч а  | 
	5.26. Вычислить время t, через которое температура  | 
||||||||
пласта достигнет 329,6 °С на расстоянии  | 
	1 м от нагнетательной сква  | 
||||||||
жины при закачке в пласт воздуха.  | 
	
  | 
	температура  | 
	нагнетае  | 
||||||
Расход воздуха  | 
	qB03 =  | 
	1000 м3/сут;  | 
|||||||
мого воздуха  | 
	Т воз = 500 сС;  | 
	объемная  | 
	теплоемкость  | 
	нагнетаемого  | 
|||||
воздуха Своз =  | 
	1,507  | 
	кДж/(м3-К); толщина  | 
	пласта  | 
	Д =  | 
	10 м; на  | 
||||
чальная пластовая температура Т 0 = 20 °С; объемная теплоемкость
пласта  | 
	с'пл = 2860  | 
	кДж/(м3-К);  | 
	теплопроводность  | 
	пласта  | 
	Хпл =  | 
|||||
= 2,4 Вт/(м- К).  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	О т в е т :  | 
	t = 8,64-105 с = 10 сут.  | 
	qB03,  | 
	при котором тем  | 
||||||
  | 
	З а д а ч а  | 
	5.27. Найти расход  | 
	воздуха  | 
|||||||
пература пласта будет равна 200 °С на расстоянии  | 
	3,3 м  | 
	через  | 
||||||||
100 сут после начала закачки горячего воздуха.  | 
	
  | 
	
  | 
	
  | 
|||||||
  | 
	Толщина  | 
	пласта  | 
	h = 10 м; начальная  | 
	пластовая температура  | 
||||||
Т = 20 °С; плотность пласта рпл =  | 
	2600 кг/м3; теплоемкость пласта  | 
|||||||||
сПл=1,1  | 
	кДж/(кг-К); теплопроводность пласта Хпл =  | 
	2,4 Вт/(м-К);  | 
||||||||
плотность воздуха  | 
	рв03 = 1,3 кг/м3; теплоемкость  | 
	воздуха С =  | 
||||||||
=  | 
	1,1  | 
	кДж/(кг-К);  | 
	температура  | 
	нагнетаемого  | 
	воздуха  | 
	7+3 =  | 
||||
=  | 
	500  | 
	°С.  | 
	qB03 =  | 
	5,787-10~2  | 
	м3/с = 5000 м3/сут.  | 
	
  | 
	
  | 
|||
  | 
	О т в е т :  | 
	
  | 
	
  | 
|||||||
§ 3. РАСЧЕТ ТЕХНОЛОГИЧЕСКИХ ПОКАЗАТЕЛЕЙ РАЗРАБОТКИ НЕФТЯНЫХ МЕСТОРОЖДЕНИЙ С ПРИМЕНЕНИЕМ ВЛАЖНОГО ВНУТРИПЛАСТОВОГО ГОРЕНИЯ И ЗАКАЧКИ В ПЛАСТ ПАРА.1 ЗАДАЧИ ДЛЯ УЧЕБНОЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКОЙ РАБОТЫ
При разработке нефтяных месторождений с применением] влаж ного внутрипластового горения (ВВГ) в пласт вместе с воздухом закачивается вода. В результате этого в пласте образуется высоко температурная зона, заполненная продуктами горения, азотом и во дяным паром. В передней части этой зоны выпариваются легкие фракции нефти и остается кокс, который затем сгорает.
183
Изменение отношения количества закачиваемой в пласт воды к количеству закачиваемого воздуха (водо-воздушного отношения) может приводить к существенному изменению температурного поля и скорости окислительной реакции. При проектировании ВВГ необходимо уметь рассчитывать его технологические параметры.
Срасчетом этих параметров связаны задачи 5.28К—5.31 К. Приближенный способ оценки показателей разработки нефтя
ного месторождения при нагнетании в него водяного пара дан в за даче 5.32К, вариант которой — задача 5.33К.
Задачи 5.34Н—5.38Н предназначены для учебной научно-иссле- довательной работы.
З а д а ч а 5.28К. Нефтяное месторождение, толщина продук тивного пласта которого составляет 17 м, насыщенное нефтью вяз костью р,н = ЮО мПа с, разрабатывается с самого начала, т. е. после ввода из промышленной разведки в разработку, с примене нием метода влажного внутри пластового горения (ВВГ).
Система разработки — однорядная,  | 
	при этом длина одного эле  | 
|||||||||
мента  | 
	(«полуэлемента») I =  | 
	400 м,  | 
	а  | 
	ширина Ь0 =  | 
	200 м (общий  | 
|||||
sc = 8-104 м2/скв).  | 
	Коэффициент  | 
	охвата  | 
	пласта  | 
	процессом  | 
	ВВГ  | 
|||||
по толщине г)2т = 0,882, а по площади т]2п =0,9.  | 
	
  | 
	
  | 
	
  | 
|||||||
Абсолютная проницаемость пласта k =  | 
	0,8 -10~12 м2; пористость  | 
|||||||||
т =  | 
	0,22; плотность  | 
	нефти  | 
	р„ =  | 
	0,85-103  | 
	кг/м3;  | 
	плотность  | 
	гор  | 
|||
ных  | 
	пород  | 
	рт =  | 
	2,5-103  | 
	кг/м3;  | 
	теплоемкость  | 
	пород  | 
	сТ =  | 
|||
=1,257 кДж/(кг-°С); коэффициент теплопроводности пород кровли
иподошвы пласта Хт = 3-(Вт/(м-°С); температуропроводность пород хт = 9-10-7 м2-с = 0,07776 м2-сут. Среднее пластовое
давление в элементе р = 102 Па.
При осуществлении влажного горения в нагнетательную сква жину закачивается воздух с расходом 0,926 м3/с = 80-103 м3/сут при Т = 30 °С, так что в каждый элемент поступает 40-103 м3/сут воздуха. По данным лабораторных исследований процесса ВВГ
содержание кокса в породе пласта составляет  | 
	zT = 28  | 
	кг/м3, теп  | 
|||||
лота  | 
	сгорания  | 
	кокса  | 
	А = 25,14-103 кДж/кг;  | 
	стехиометриче  | 
|||
ский  | 
	коэффициент а =  | 
	11,2 м3/кг. Пластовая  | 
	температура Тпл =  | 
||||
- 30 °С.  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
Начальная нефтенасыщенность  | 
	пласта s„0 =  | 
	0,93;  | 
	насыщен  | 
||||
ность  | 
	связанной  | 
	водой  | 
	sCB= 0,07.  | 
	Относительные  | 
	проницаемости  | 
||
пласта для нефти, воды и газов зависят линейно от соответствую щих насыщенностей. Остальные исходные данные, необходимые для расчета ВВГ, указаны в решении задачи.
Технология осуществления ВВГ в элементе системы разработки такова, что в начальной стадии ВВГ в пласте создаются высоко температурные зоны (паровое плато, зона горения) таким образом, что фронт конвекции хт(рис. 67) движется быстрее, чем фронт го рения, в 1,5 раза. Начальная стадия ВВГ продолжается до момента
времени  | 
	t = tlt когда фронт горения достигнет расстояния Хф =  | 
||
= 50 м.  | 
	Затем наступает стадия установившегося ВВГ,  | 
	когда ско  | 
|
рость  | 
	фронта конвекции vT = dxT/dt равна скорости  | 
	фронта го  | 
|
рения  | 
	Уф = d ^/d t. Эта последняя стадия продолжается до конца  | 
||
184
1 Z  | 
	3 4  | 
Рис. 67. Схема влажного внутрипластового горения (ВВГ). Зоны:
1 — /; 2 — горения / / ; 3 — парового плато I I I ; 4 — I V
разработки элемента пласта методом ВВГ, т. е. до момента времени, когда х? = I.
Требуется определить для рассматриваемого элемента пласта
температуру Т ш в паровом  | 
	плато в зоне / / / , водо- и газонасыщен-  | 
ности в зонах I, III и IV,  | 
	водовоздушное отношение Хвв на входе  | 
в пласт, дебиты нефти и воды, обводненности продукции в среднем за начальную стадию и в течение всего остального времени ВВГ.
На основе экспериментальных и теоретических исследований процесса принимается, что средняя температура в зоне/ / постоян ная и составляет Т% = 330 °С, а длина зоны / / Д£ = 1 м.
Р е ш е н и е . 1. Расчет начальной стадии создания высокотем пературной зоны в пласте.
При определении технологических показателей разработки пласта методом ВВГ, а также параметров самого процесса ВВГ применим одну из разновидностей так называемых зонных моделей, описание которых можно найти, например, в книге [5].
Согласно модели, используемой при решении настоящей задачи, поле температуры и насыщенностей в пласте при ВВГ разбивается на несколько характерных зон. На рис. 67 сплошной линией по казано распределение температуры в пласте, наблюдаемое, напри мер, в лабораторных экспериментах, а пунктирной — расчетное распределение температуры, принимаемое в «зонной» модели про-
185
цесса, на основе которой будут рассчитываться технологические показатели ВВГ
Всоответствии со схемой распределения веществ и температуры
впласте, показанной на рис. 67, в зоне I в пласте совместно дви жутся воздух и вода, в узкой зоне II происходит основная часть внутрипластовой окислительной реакции. В этой зоне находятся водяной пар, кислород, реагирующий с нефтяным остатком — кок сом, газообразные вещества, остающиеся или получающиеся в ре зультате реакции горения (азот, углекислый газ, окись углерода),
именуемые в дальнейшем как «газы горения». Водяной пар в зоне
/ / во многих случаях, по-видимому, бывает перегретым.  | 
	
  | 
||||||||
  | 
	Зона I I I (паровое плато) содержит в основной своей части газы  | 
||||||||
горения и водяной пар, а в части, прилегающей к зоне IV,  | 
	где про  | 
||||||||
исходят сложные  | 
	процессы дистилляции  | 
	легких  | 
	фракций нефти  | 
||||||
и образование кокса,— смесь всех движущихся  | 
	в  | 
	пласте веществ.  | 
|||||||
Однако приближенно принимается, что  | 
	в зоне  | 
	I I I  | 
	содержатся  | 
||||||
только газы горения и водяной пар. В зоне IV происходит совмест  | 
|||||||||
ное движение газов горения, нефти и воды.  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||||
  | 
	Граница между III и IV зонами считается фронтом конвекции  | 
||||||||
с координатой хт,  | 
	а за  | 
	координату фронта горения  | 
	Хф,  | 
	условно  | 
|||||
для  | 
	удобства расчетов,  | 
	принимается  | 
	граница  | 
	между зонами II  | 
|||||
и / / / .  | 
	
  | 
	
  | 
	оф, исходя из стехиометрии  | 
||||||
  | 
	Определим скорость фронта горения  | 
||||||||
горения.  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	Имеем следующую формулу:  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||
  | 
	Уф — Ц0г/-^В03>  | 
	Уот = Явоз^г  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	(5»14)  | 
|
где  | 
	<7во3 — расход  | 
	закачиваемого в пласт воздуха  | 
	в  | 
	нормальных  | 
|||||
условиях, м3; 5 — площадь поперечного  | 
	сечения  | 
	пласта;  | 
	3 —  | 
||||||
объем воздуха, требующийся для выжигания кокса из 1 м3 пласта в целом, м3.
Соответственно для # в03  | 
	имеем  | 
	формулу  | 
R воз —  | 
	
  | 
	(5.15)  | 
где zx — содержание кокса  | 
	в 1 м3  | 
	пласта, кг/м3; а — стехиомет  | 
рический коэффициент.  | 
	процесса ВВГ, проведенным примени  | 
|
Согласно исследованиям  | 
||
тельно к конкретным условиям рассматриваемого месторождения, а = 11,2 м3/кг. По условию задачи zr = 28 кг/м3. Следовательно, Яв03 = 11,2-28 = 313,6 м3/м3.
Поскольку по условию задачи коэффициент охвата пласта про
цессом  | 
	по  | 
	площади т)2П =  | 
	0,9, а по  | 
	толщине  | 
	т)2Т =  | 
	0,882,  | 
	то за  | 
|||
эффективную  | 
	ширину пласта  | 
	примем  | 
	b =  | 
	Ь0 0,9 = 200 0,9 =  | 
||||||
— 180  | 
	м,  | 
	а  | 
	за эффективную толщину пласта h = /г0 -0,882 =  | 
	17X  | 
||||||
X0,882 =  | 
	15 м. Тогда 5 =  | 
	bh =  | 
	180-15 =  | 
	2700 м2.  | 
	м3/с =  | 
	40X  | 
||||
По  | 
	формуле (5.14) с учетом того,  | 
	что  | 
	<7воэ= 0,46  | 
|||||||
Х103 м3/сут  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||
  | 
	
  | 
	
  | 
	40 -10s  | 
	= 0,547-10~6 м/с = 0,04724  | 
	м/сут-  | 
	
  | 
||||
  | 
	2700-313,6 0,864-10»  | 
	
  | 
||||||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
1S6
По условию задачи в начальной стадии ВВГ скорость фронта конвекции ут = 1,5 Уф = 1,5 0,04724 = 0,07086 м/сут.
Перейдем к определению температуры Т щ в зоне III, учитывая, что по условию задачи температура в зоне окислительной реакции II задана и остается в течение всего процесса ВВГ постоянной, равной Т* = 330 °С.
Вначале получим формулу для текущего ухода тепла из зоны I I I в кровлю и подошву пласта. В соответствии с условием задачи
и ее схематизацией будем считать,  | 
	что температуры  | 
	в зонах I и  | 
||
IV равны  | 
	пластовой, так что  | 
	=  | 
	7 IV = Тпл. Будем  | 
	также при  | 
ближенно  | 
	принимать, что средняя температура в зоне I I I остается  | 
|||
неизменной в течение всей начальной стадии процесса ВВГ. Если бы температура в зоне / была равной температуре в зоне III, т. е. при 0 < х < хтзначение Т = Тщ, то тепловые потери в кровлю и подошву пласта за счет теплопроводности можно определить по формуле
t
1  | 
	
  | 
	(* 2Хт (Гщ  | 
	Т пл) vTbdx  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	j  | 
|||||
* Ш“  | 
	1 ----- [п*, (< -*)]'2  | 
	•  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	(5Лб)  | 
|||||||
  | 
	
  | 
	о  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
Однако согласно схеме, представленной на рис. 67, температура  | 
|||||||||||||||
при х <  | 
	Хф снижается  | 
	скачком до пластового значения. Если же  | 
|||||||||||||
рассматривать только уход тепла из зоны III, то необходимо при  | 
|||||||||||||||
нимать, что уже при х <  | 
	
  | 
	Хф температура снижается до пластовой.  | 
|||||||||||||
При этом нужно учитывать, что граница  | 
	зон I I I  | 
	и IV движется  | 
|||||||||||||
со скоростью  | 
	ут,  | 
	а  | 
	граница  | 
	зон II  | 
	и  | 
	I I I — со скоростью Уф.  | 
|||||||||
Учитывая  | 
	сказанное,  | 
	получаем  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||||||
  | 
	
  | 
	t 2А»т  | 
	( 7*H  | 
	I  | 
	ТПл) vTbdx  | 
	
  | 
	i  | 
	2ХТ (Т HI — Тпл) ифbdx  | 
|||||||
?тШ = ^  | 
	[яхт (t —т)]12  | 
	
  | 
	
  | 
	[лхт (f — т)]1,2  | 
|||||||||||
  | 
	
  | 
	о  | 
	
  | 
	
  | 
|||||||||||
  | 
	
  | 
	
  | 
	г  | 
	(^т  | 
	
  | 
	Уф) bdx  | 
	
  | 
	л т",  | 
	'г  | 
	
  | 
	т  | 
	
  | 
|||
=  | 
	2Arj-АТ п I f  | 
	
  | 
	’  | 
	
  | 
	(5.17)  | 
||||||||||
(п*т)"2  | 
	'1  | 
	(f- т ) ' 2  | 
	4 Г - Г ш  | 
	- Г п -  | 
|||||||||||
По  | 
	условию  | 
	задачи  | 
	
  | 
	ут =  | 
	1,5,  | 
	
  | 
	Уф =  | 
	const.  | 
	Следовательно,  | 
||||||
ут— Уф =  | 
	0,5 уф. В этом случае  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||||||
йтП1 =  | 
	2ХтД'Гц10,5иф6  | 
	Г  | 
	dx  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||||
  | 
	
  | 
	(лхт) ь  | 
	
  | 
	оS  | 
	( t -  | 
	X)12  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	t  | 
	\1:2  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	(5.18)  | 
||
= 2ЯтДТи1&Цф ( ------^  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||||||
В более общем же случае, когда существует другое, но постоян  | 
|||||||||||||||
ное отношение  | 
	ут/Уф,  | 
	имеем  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||||
? Т Ш = г ^ Д Г ш Н (  | 
	V ф  | 
	— 1 ) ( — — )'я  | 
	
  | 
	( 5. 19)  | 
|||||||||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	\  | 
	
  | 
	
  | 
	/  | 
	\ ЯХт  | 
	/  | 
	
  | 
	
  | 
	
  | 
||
187
Всоответствии с условием задачи зона окислительной реакции
сне изменяющейся во времени шириной Д£ и температурой Т * движется в пласте также со скоростью Оф.
в  | 
	Обозначим через  | 
	t* время, за которое фронт горения проходит  | 
|
пласте расстояния  | 
	Д£. В этом случае t%= А|/«ф.  | 
	пласта  | 
|
в  | 
	Выражение для текущих тепловых потерь из зоны II  | 
||
кровлю и подошву  | 
	получаем аналогично формуле (5.17).  | 
	Имеем  | 
|
„  | 
	_ 2Хт6иф (Т%  | 
	Тпл)  | 
|
^11  | 
	/  | 
	N|/2  | 
	
  | 
  | 
	(лхт)1'2  | 
	
  | 
|
_  | 
	2\тЬифАТ*  | 
	г'  | 
	
  | 
  | 
	(лхт)Ь2  | 
	Lo  | 
	
  | 
-  | 
	(п*т)  | 
	[2  | 
	+  | 
  | 
	
  | 
	
  | 
|
г t  | 
	
  | 
	dx  | 
	t-U  | 
	dx  | 
Г  | 
	
  | 
	Г  | 
||
)  | 
	V - ' ) 112 )  | 
	(t-т;)112  | 
||
Lo  | 
	
  | 
	
  | 
	t*  | 
	1 Н и1  | 
------  | 
	—  | 
	1 ьо  | 
||
- т  | 
	
  | 
	1  | 
	
  | 
|
  | 
	
  | 
	U  | 
	1  | 
|
  | 
	
  | 
	- 2  | 
	V t - U  | 
	]•  | 
При t > t* имеем
ЧП1 = 4 )ф > * Ь Т ,( - ± - ') .  | 
	(5.20)  | 
V лхт /
Поскольку согласно условию задачи из зон I и IV тепло не уходит в кровлю и подошву пласта (TY = T1V = Тпл), то при / > полные тепловые потери в кровлю и подошву qT можно оп ределить следующим образом:
= ^7тц+ <7т1ц =  | 
	2Я,тЬиф j^A r in  | 
	^ -f-  | 
+ 2ДТ, ( ———У 2] .  | 
	(5.21)  | 
|
V лхт /  | 
	J  | 
	
  | 
Тепловые потери в кровлю и подошву пласта QT, накопленные ко времени t , получают интегрированием (5.21). Имеем
Q T = Q T II + Q T I I I —  | 
	4Хт6Д7'.1.Оф ( — 5— ^  | 
	t +  | 
  | 
	V лхт /  | 
	
  | 
+ -?-М щ ,А Г „,-  | 
	^ v .  | 
	(5.22)  | 
3(лхт)
Всоответствии со схемой процесса ВВГ, представленной на рис. 67, для количества тепла Qnjl, накопленного к моменту вре мени t в самом пласте, имеем следующую формулу:
<2пл = Стрт(1 —m) bhA Т ш (от—оф) t +
н - m C epaSni ATu l(vT — Оф) bht - f
+  | 
	m (срРп + сгРг) (1 — 5 щ ) Ы гА Т щ (ит — Цф) t -}-  | 
	
  | 
+  | 
	СтРт (1 — т) Д Т *& Л Д £ ,  | 
	(5.23)  | 
188
где sn  | 
	и sm — соответственно насыщенности  | 
	водой  | 
	в  | 
	зонах II  | 
и III.  | 
	Количество тепла Q*, выделившегося в  | 
	пласте  | 
	к  | 
	моменту  | 
времени t в результате реакции горения, выражается следующим образом:
Q* ~ z TAbh(v^t — А£).  | 
	(5.24)  | 
Согласно балансу тепла имеем  | 
	
  | 
Q* = QT+ Q™.  | 
	(5.25)  | 
Как видно из структуры формул, с помощью которых опреде  | 
|
ляли Q*, QT и <2пл, члены, входящие в указанные выше  | 
	величины,  | 
по-разному зависят от времени t. Это означает, что либо в случае постоянства АТ* и ATlU в зоне I I I будет изменяться насыщенность sIH пласта жидкой водной фазой, либо в случае существования в зоне I I I только перегретого пара АТщ не останется постоянной. Однако при сравнительно небольшом продвижении в пласте фронта горения АТ 1п может изменяться в сравнительно небольших преде лах, так что приближенно в этих случаях можно считать, что АТ т
равна среднему значению ATjn в зоне I I I и постоянна в течение рас
сматриваемого времени протекания ВВГ.  | 
	с зоны  | 
	IV,  | 
|
Расчет  | 
	технологических показателей ВВГ начнем  | 
||
в которой  | 
	необходимо определить насыщенности зоны  | 
	нефтью,  | 
	га  | 
зами горения и водой и, как наиболее важный показатель, скорость
фильтрации воды uBiv в зоне 1У-
Для этого рассмотрим совместное движение газов горения, нефти и воды в зоне IV
Обозначая насыщенность зоны IV газами горения srIV водой
sIV,  | 
	нефтью sHiv, имеем  | 
	
  | 
	
  | 
	
  | 
||
  | 
	•Sriv _1_SHIV “Г siv =  | 
	1*  | 
	
  | 
	
  | 
	(5.26)  | 
|
Для скоростей фильтрации соответствующих веществ известны  | 
||||||
следующие выражения:  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	Vr —  | 
	kkv (sr)  | 
	dp  | 
	_  | 
	kkH(sH)  | 
	dp  | 
  | 
	-------------------------1 ,  | 
	VH —  | 
	■  | 
	,  | 
||
  | 
	
  | 
	(Xr  | 
	dx  | 
	
  | 
	|i„  | 
	dx  | 
  | 
	Vв  | 
	kkB(sB)  | 
	dp  | 
	
  | 
	
  | 
	(5.27)  | 
  | 
	Ив  | 
	dx  | 
	
  | 
	
  | 
||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||
где  | 
	kr,  | 
	kn и /гв — относительные  | 
	проницаемости соответственно  | 
|||
для газа, нефти и воды.
Примем, исходя из условия задачи, что относительные прони цаемости в зоне IV для газов, нефти и воды линейно зависят от со ответствующих насыщенностей, а именно
К  | 
	Sp — Spo  | 
	kB  | 
	S —  | 
	S e e  | 
	(5.28)  | 
||
1  | 
	
  | 
||||||
1  | 
	SCB  | 
	SCB  | 
|||||
  | 
	
  | 
	
  | 
|||||
где sro — насыщенность пористой  | 
	среды газом, при которой про  | 
||||||
ницаемость для газа равна нулю;  | 
	sCB— насыщенность пласта свя  | 
||||||
занной водой; s* — водогазонасыщенность, при  | 
	которой проницае  | 
||||||
мость для  | 
	нефти становится равной нулю.  | 
	
  | 
	
  | 
	
  | 
|||
189
Из зоны I I I  | 
	в зону IV входит количество нефти,  | 
	пропорцио  | 
|||||
нальное скорости  | 
	ыт> а также первоначальному содержанию нефти  | 
||||||
в пласте за вычетом количества сгоревшего кокса, так что  | 
|||||||
и„ = ^ms„0---- vTt  | 
	
  | 
	
  | 
	
  | 
	(5.29)  | 
|||
где рк — плотность  | 
	кокса.  | 
	0,95* 103 кг/м3  | 
	и  | 
	учитывая,  | 
|||
Вычислим  | 
	UHIV-  | 
	Принимая рк =  | 
|||||
что т = 0,22; sH0 =  | 
	0,93; sCB— 0,07;  | 
	vr = 1,5 иф =  | 
	0,07086 м/сут,  | 
||||
получим  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
uHiV =  | 
	Г0,22 -0,93---------------- 0,07086 - 0,0124 м/сут.  | 
	
  | 
|||||
  | 
	V  | 
	
  | 
	
  | 
	0,95-10_3 J  | 
	
  | 
	J  | 
	
  | 
Составим  | 
	отношение  | 
	
  | 
	
  | 
	
  | 
|||
и г  | 
	________ Р н  | 
	(s rIV — s r o )  | 
	
  | 
	
  | 
	(5.30)  | 
||
Z>HIV  | 
	Н-г (s* — SIV — Sriv)  | 
	
  | 
	
  | 
||||
  | 
	
  | 
	
  | 
|||||
Будем полагать, согласно закону идеальных газов, при среднем дав лении в пласте р = 107 Па, что
yriv — ^го Ро 40-10М05 = 0,1481 м/сут.
р180-15-107
При sr0 = 0,05  | 
	значение  | 
	s* = 0,95.  | 
	Учитывая, что в зоне IV  | 
||
рн =  | 
	ЮО  | 
	мПа-с,  | 
	рг = 0,02  | 
	мПа-с, в  | 
	соответствии с формулой  | 
(5.30)  | 
	получаем  | 
	
  | 
	
  | 
	
  | 
|
0,1481_________ 100 (SrIV — 0,05)  | 
	
  | 
||||
0,0124  | 
	0,02(0,95 — sjy — sriy)  | 
	
  | 
|||
Выпишем отношение  | 
	
  | 
	
  | 
|||
a BIV _  | 
	Рн (SIV — $св)  | 
	
  | 
	(5.31)  | 
||
«HIV  | 
	MB (S* — «IV — Sriv)  | 
	
  | 
|||
  | 
	
  | 
||||
Из этого отношения получаем выражение для uBlv:
l>uiv  | 
	1,24 (Siv — 0,07)  | 
	(5.32)  | 
|
0,9 — siv  | 
|||
  | 
	
  | 
Таким образом, на основе теории установившегося трехфазного течения в зоне IV получены уравнения (5.30) и (5.31) для опреде ления трех неизвестных uBlv, sIV и sriv. Недостающее третье уравнение получаем на основе выражения для скорости конвек
тивного переноса тепла на границе зон II I  | 
	и IV в виде  | 
	
  | 
||
vT  | 
	_______________СгРг^Г ~Т Снр нЦц 4~ СврвРв_______________  | 
	(5.33)  | 
||
стрт (1 — т) + т [Сцр„ (1 — S — sr) + cBpDs + crprsr]  | 
||||
  | 
	
  | 
|||
Систему уравнений (5.30), (5.31) и (5.33) будем решать методом последовательных приближений.
190
