Миклина О.А., Корохонько О.М. Обработка результатов исследования скважины на неустановившихся режимах работы
Подождите немного. Документ загружается.
11
В файл добавить лист. Для этого навести курсор мышки к меню Вставка
и выбрать Лист (рис. 1.19).
Рисунок 1.19 – «Добавить» лист
Скопировать данные на новый лист, начиная с точки закрытия скважины
(с точки 74,4). Для этого нужно выделить данные из листа «Пласт_Протокол»
одним из следующих способов:
¾ выделить строку с временем 74,4. Прокрутить таблицу вниз до
конечной точки и выделить последнюю строку с нажатой клавишей
Shift;
¾ выделить строку
с временем 74,4. Нажать комбинацию клавиш
Crtl+Shift+End;
¾ выделить нужные строки мышкой.
Переименовать новый лист (Лист2) – КВД.
Добавить шапку таблицы (рис. 1.20).
Рисунок 1.20 – Готовая таблица
Для построения КВД из таблицы удалить столбцы с показаниями
датчиков от Z1 до Z7, т.к. показания этих датчиков не потребуются в
дальнейших расчетах. Для этого выделить столбцы B – H, нажать правую
кнопку мыши и выбрать пункт «Удалить» (рис. 1.21).
12
Рисунок 1.21 – Шаг «Удалить» столбики
Получается таблица:
Добавляем столбик « Т, с ». Вводим в ячейку В2 формулу «
2$$2 AA
−
=
»
($A$2 – абсолютная ссылка, она необходима для того, чтобы при копировании
формулы адрес ячейки не менялся, и она появляется после нажатия клавиши
F4) (рис. 1.22).
Рисунок 1.22 – Шаг «Добавить» столбик «Т,с»
Далее необходимо скопировать формулу во все ячейки столбика В для
чего:
o подвести курсор мыши к маркеру в правом нижнем углу ячейки В2
(он примет вид
);
o протянуть формулу до конца таблицы.
Датчиком Z8 регистрируется давление на забое скважины, поэтому
нужно переименовать ячейку С1 – Pзаб(t). Для построения КВД нужно
добавить еще один столбик (столбик D) для расчета разницы забойных
давлений ∆Pзаб (t) = Pзаб(t) - Pзаб(0). Перевести это значение в систему СИ.
Формула в ячейке будет иметь следующий
вид «=(C2-$C$2)*10^5» (рис.
1.23).
13
Рисунок 1.23 – Шаг «Добавить столбик» ∆Pзаб (t)
Копируем формулу и получаем столбик значений ∆Pзаб(t), Па.
Строим КВД. Для этого нужно выделить столбик B (Т,с) и, удерживая
нажатой клавишу Ctrl, выделить столбик D (∆Pзаб(t), Па). Нажать на кнопку
«Мастер диаграмм». Выбрать тип диаграммы – Точечная, со значениями,
соединенными отрезками
без маркеров (рис. 1.24).
Рисунок 1.24 – Шаг «Тип диаграммы»
На следующем шаге выбрать Ряды в столбцах.
На третьем шаге установить параметры графика и подписать оси
(рис. 1.25).
Рисунок 1.25 – Шаг «Установить параметры диаграммы»
На завершающем шаге определить место расположения графика – на
отдельном листе. После этого на экране будет диаграмма (рис. 1.26).
14
Рисунок 1.26 – Диаграмма
Необходимо определить положение точки, где начинается спад (точка 5
рис. 1.2). Для этого подвести курсор мыши к этой точке, всплывает подсказка,
которая показывает координату точки.
Рисунок 1.27 – Точка закрытия скважины
В нашем случае это Т=21,7 с. В исходной таблице, начиная с этой точки,
удалить все последующие строки и график примет вид, как на рис. 1.28.
0
10000
20000
30000
40000
50000
60000
0 5 10 15 20 25
t, c
∆Pзаб(t), Па
Рисунок 1.28 – График кривой восстановления давления
15
2 ОБРАБОТКА РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ
Как известно, задачей гидродинамических исследований является
изучение коллекторских и фильтрационных свойств пласта. По результатам
исследований скважин на неустановившихся режимах определяют значения
различных коэффициентов таких как, коэффициент гидропроводности,
коэффициент эффективной проницаемости, коэффициент пьезопроводности, а
также значение комплексного параметра.
В [1] описана методика обработки КВД классическим способом –
методом касательной, а в данных методических указаниях
описана методика
обработки исследований, проведенных на экспериментальной установке
«Пласт».
Описанная ниже методика обработки КВД отличается от уже известных
методик тем, что радиус контура соизмерим с радиусом скважины
(Rк = 1,50 м, Rс = 0,5 м), и КВД в данном случае обрабатывается по формуле [2]:
()
[]
()
,221,167,0lnln
2
t
R
Atp
c
χ
−⋅=∆
(2.1)
где: ∆Р(t) – значение депрессии после остановки скважины, Па;
χ – коэффициент пьезопроводности пласта, м
2
/с;
R
c
– радиус скважины, м;
t – время, с;
А – коэффициент фильтрации, рассчитываемый по формуле:
,
2
1
α
π
µ
⋅
⋅
⋅
⋅
⋅
=
h
k
Q
А
(2.2)
где: α
1
– часть кругового пласта – сегмент установки «Пласт»:
,028,0
360
10
360
1
===
α
α
где α – угол сегмента, 10
0
С (рис. 2.1).
Рисунок 2.1 – Схема пласта
α
Ось скважины
16
Для получения значений коэффициентов необходимо выполнить
следующие действия:
− перевести показания датчика Z8 в систему СИ (Па);
− найти логарифм от разницы давления на контуре и на скважине. Для
этого вводим в ячейку F2 формулу «=LN($G$218-G2)». Ячейка G218 –
последняя ячейка в столбике и это значение давления на контуре
питания (рис. 2.2);
Рисунок 2.2 – Логарифм от разницы давления на контуре и на скважине
− копировать формулу во все ячейки столбика F;
− построить график зависимости ln(Pk-P(t)) от Т. Последовательность
действий аналогична построению предыдущих графиков.
После получения графика нужно определить положение двух точек на
графике, которые образуют прямую линию. Сначала визуально определяем этот
отрезок на
графике. Далее подвести курсор к крайним точкам этого отрезка и
записать значения времени из всплывающих подсказок (рис. 2.2).
Рисунок 2.2 – Конечные точки прямой
В нашем случае это время 0,5 с и время 6,1 с. Теперь нужно провести
через все точки, находящиеся в интервале [0,5 – 6,1], прямую линию.
Для этого:
o подвести курсор к графику и выделить его щелчком правой кнопки
мышки;
17
o вызвать контекстное меню;
o выбрать Исходные данные;
o в открывшемся окне открыть вкладку Ряд и нажать на кнопку Добавить.
В окне Ряд появляется дополнительно Ряд2;
o нажать на кнопку
«Значения Х:» ,
o открывается окно
o открыть лист «КВД»;
o выделить курсором в столбике «Т,с» ячейки, начиная от времени 0,5 и
кончая временем 6,1;
o в строке «Исходные данные – Значения Х» появляются адреса ячеек
18
o нажать на кнопку строки «Исходные данные – Значения Х»;
o возвращаемся в окно «Исходные данные»;
o в открывшемся окне нажать на кнопку
«Значения Y»
o выделить на листе «КВД» курсором в столбике «ln(Pk-P(t))» ячейки,
соответствующие отрезку времени от 0,5 до 6,1 с;
o нажать на кнопку ОК;
o на исходном графике строится дополнительный ряд (обычно выделяется
розовым цветом) – требуемый отрезок.
Для того чтобы получить уравнение линии тренда необходимо выполнить
следующее:
o подвести курсор к новому отрезку графика;
o вызвать контекстное меню;
o выбрать «Добавить линию тренда»;
o указать – Тип линии тренда «Линейная»;
o в параметрах линии тренда установить: «Показывать уравнение на
диаграмме» (рис. 2.3).
19
Рисунок 2.3 – Шаг «Установка параметров линии тренда»
После всего этого на графике будет проведена линия тренда и показано
уравнение (рис. 2.4).
Рисунок 2.4 – Конечный вид графика
Итак, получено уравнение:
.047,13489,0
+
⋅
−
=
xy
Сначала необходимо рассчитать значение объемного расхода воды Q. Для
этого в исходной таблице (см. рис. 2.2) на листе «КВД» надо записать первое
значение показания датчика Z10. Оно равно 0,7968 л/мин. Перевести это
значение в систему СИ.
Согласно формуле (2.1)
(
)
047,1367,0ln
=
⋅
A
, откуда
5
106,467,0 ⋅=⋅ А .
Распишем все известные значения, которые входят в состав
фильтрационного коэффициента А:
kk
А
⋅⋅
⋅
=
⋅⋅⋅⋅⋅
⋅⋅
=
−
−
−
−−
2
8
2
35
104,1
103,1
108028,014,32
10103,1
.
И получаем
5
2
8
106,4
104,1
103,1
67,067,0 ⋅=
⋅⋅
⋅
⋅=⋅
−
−
k
А
Отсюда находится значение коэффициента эффективной проницаемости:
k = 1,34*10
-12
м
2
= 1340 мкм
2
= 1,34 дарси.
Далее оценим значение проницаемости из углового коэффициента.
20
В нашем случае угловой коэффициент равен 0,489.
Тогда, согласно формуле (2.1)
489,0221,1
2
=
c
R
χ
.
Известно, что коэффициент пьезопроводности рассчитывается по
формуле
*
βµ
χ
⋅
=
k
(2.3)
где:
µ
– динамическая вязкость жидкости, Па*с;
β
*
– коэффициент объемного упругого расширения.
скелетаводы
m
βββ
+= *
*
(2.4)
где: m – пористость пласта, доли единицы;
β
воды
– коэффициент объемного упругого расширения воды;
β
свелета
– коэффициент объемного упругого расширения скелета породы.
Примем
µ
=10
-3
Па*с, m = 0,2,
П
а
П
а
скелетаводы
1
10,
1
10*3
710 −−
==
ββ
.
Подставляем значения в формулу (2.4) и далее в формулу (2.3). Получаем
489,0
10105,0
221,1
732
=
⋅⋅
−−
k
.
Откуда коэффициент
дарсимкммk 10100121010
10*884,4
489,0
2212
10
≈≈⋅≈=
−
.