Главная  Переработка нефти и газа 

Скачать эту книгу

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 [ 124 ] 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199

лом значении "к можно свести к следующему:

(8.22>

откуда по аналогии с уравнением (7.31)

Аб2 = -г

- (й

Выражение для Лбг в размерном виде можно переписать так: [т == - Ei [- Wg/(I - т)\, (8.22>

где ЛРх снимается, как показано на рис. 8.2, а т = т12,Ъ.

б. Пример расчета

Между двумя скважинами (добывающей и наблюдательной), расстояние между которыми Z)=350 м, были проведены исследования интерференции. Дебит добывающей скван<;ины 1750 м/сут, в наблюдательной скважине измерялось давление во времени (рис. 8.3).

i=33if30 С 1

0,035 -

Рис. 8.3. Кривая зависимости падения давления от времени в наблюдательной скважине по регистрированным данным

Остальные данные следующие: (.Ih=1,2 мПа-с, 5н=1,27; h= = 109 м; Ф1 = 0,15 (по керну и геофизическим данным); Ci = 2,6x ХЮ-з 1/МПа, С2 = 1,3-10- 1/МПа; Гскв=10см.

Из графика АР-log t (см. рис. 8.3) в полулогарифмическом масштабе находим:

т = 0,08 МПа/цикл.

Трещинная проницаемость вычисляется непосредственно из уравнения (8.17):

-тр.пл = 1.15 » 0,65 мкм1

Комплексный коэффициент ньезопроводности определяется из уравнения (8.18):

а = 0,445 r-/to = 13 825 тУс.

Экстраполяция асимптоты до ЛР = 0 дает /о = 39430 с. Из уравнения (8.19)

ФС =/С/ра= 3,918-10- 1/МПа.

По уравнению (8.20) рассчитывается трещинная пустотность:

Фр да Фа = (ФС - Ф1С,)/С2 да 0,00138.

Уравнение (8.21) дает безразмерную характеристику вмещающей способности пласта:

(О = ф Со/ФС = 0,00457.

Расчет К. По кривой поведения давления на рис. 8.3 падение давления при t=to составляет ЛР = ЛРх =0,00118 МПа. Таким образом, из уравнения (8.22) получаем:

ЛРх cmEi

1-(й Л 0,0118 = (0,08/2,3) Ei

= 0,00118 МПа;

I-3825- здзр

0,995x100

0,719 =>.х 5,47 х10«; X да 1,312 X 10"

8.2.3.1. Метод решения для оценки источника

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

AWTp.n.n) = - -7-0 = 0, (8.23)

Рис. 8.4. Схема идеализированного пласта.

2-модель / - слоистая матрица; б - модель Р-матрица в виде блоков, / - повторяющиеся элементы; Af - число повторяющихся элементов; о - приведенная пьезопроводность; Г - гидропроводиость

Наджуриет [3] предложил использовать приближенное решение для начальных и граничных условий, выраженное решением для линейного источника:

тр.пл

4ас.д/

(8.24)

(8.25)-

(8.26) (8.27)

4т1ЛГтр

где йс.д - составная диффузивность,

khl Khl\l

(Фсн)тР + (ФсН)м*" ~ (ФС)с.д (ФС)с.д

При этом для моделей 1 и 2, представленных на рис. 8.4, получим соответственно: модель /

=]/){апЬ]/8У7;

модель 2

Ф = уТРг coth ]b/t - т. ).

Параметр б связан с временем реагирования матрицы:

В = hlla. (8.28)

При исследовании интерференции скважин с использованием уравнения (8.24) можно отметить следующее.

Для меньших значений времени отклонение от поведения пласта с одним видом пустотности является следствием резкого различия таких характеристик трещинно-матричной системы, как

(ФсЬ)тр и (Och)m.

Для больших значений времени в соответствии с геометрией типа модели результаты будут следующими:

для модели 1 (слои)

«c.A.=7/[(®ch),p-f (Och)J; (8.29)

для модели 2 (кубы)

ас.д 2 = ЗГ/[3 (Ф ch),p -Ь (Ф сЬ)м1. (8.30)