Обработка результатов исследования
Как известно, задачей гидродинамических исследований является изучение коллекторских и фильтрационных свойств пласта. По результатам исследований скважин на неустановившихся режимах определяют значения различных коэффициентов таких как, коэффициент гидропроводности, коэффициент эффективной проницаемости, коэффициент пьезопроводности, а также значение комплексного параметра.
В [1] описана методика обработки КВД классическим способом – методом касательной, а в данных методических указаниях описана методика обработки исследований, проведенных на экспериментальной установке «Пласт».
Описанная ниже методика обработки КВД отличается от уже известных методик тем, что радиус контура соизмерим с радиусом скважины
(Rк=1,50 м, Rс=0,5 м), и КВД в данном случае обрабатывается по формуле [9]:
. (9)
где: ΔР(t) – значение депрессии после остановки скважины, Па;
χ – коэффициент пьезопроводности пласта, м2/с;
RС – радиус скважины, м;
t – время, с;
А – коэффициент фильтрации, рассчитываемый по формуле:
. (10)
где: α1 – часть кругового пласта – сегмент установки «Пласт»:
.
(11)
где α – угол сегмента, 10° (рис. 59).
Рисунок 59 – Схема пласта
Для получения значений коэффициентов необходимо выполнить
следующие действия:
перевести показания датчика Z8 в систему СИ (Па);
найти логарифм от разницы давления на контуре и на скважине. Для этого вводим в ячейку F2 формулу «=LN($G$218-G2)». Ячейка G218 – последняя ячейка в столбике и это значение давления на контуре питания (рис. 60);
копировать формулу во все ячейки столбика F;
построить график зависимости ln(PK-P(t)) от Т. Последовательность действий аналогична построению предыдущих графиков.
Рисунок 60 – Логарифм от разницы давления на контуре и на скважине
После получения графика нужно определить положение двух точек на графике, которые образуют прямую линию. Сначала визуально определяем этот отрезок на графике. Далее подвести курсор к крайним точкам этого отрезка и записать значения времени из всплывающих подсказок (рис. 61).
Рисунок 61 – Конечные точки прямой
В нашем случае это время 0,5 с и время 6,1 с. Теперь нужно провести через все точки, находящиеся в интервале [0,5 – 6,1], прямую линию.
Для этого:
подвести курсор к графику и выделить его щелчком правой кнопки мышки;
вызвать контекстное меню;
выбрать Исходные данные;
в открывшемся окне открыть вкладку Ряд и нажать на кнопку Добавить. В окне Ряд появляется дополнительно Ряд2;
нажать на кнопку «Значения Х:»;
открывается окно;
открыть лист «КВД»;
выделить курсором в столбике «Т,с» ячейки, начиная от времени 0,5 и кончая временем 6,1;
в строке «Исходные данные – Значения Х» появляются адреса ячеек;
нажать на кнопку строки «Исходные данные – Значения Х»;
oвозвращаемся в окно «Исходные данные»;
в открывшемся окне нажать на кнопку «Значения Y»
выделить на листе «КВД» курсором в столбике «ln(Pk-P(t))» ячейки, соответствующие отрезку времени от 0,5 до 6,1 с;
нажать на кнопку ОК;
на исходном графике строится дополнительный ряд (обычно выделяется розовым цветом) – требуемый отрезок:
Для того чтобы получить уравнение линии тренда необходимо выполнить следующее:
подвести курсор к новому отрезку графика;
вызвать контекстное меню;
выбрать «Добавить линию тренда»:
указать – Тип линии тренда «Линейная»;
в параметрах линии тренда установить: «Показывать уравнение на диаграмме» (рис. 62).
Рисунок 62 – Шаг «Установка параметров линии тренда»
После всего этого на графике будет проведена линия тренда и показано уравнение (рис. 63).
Рисунок 63 – Конечный вид графика
Итак, получено уравнение:
. (11)
Сначала необходимо рассчитать значение объемного расхода воды Q. Для этого в исходной таблице (см. рис. 2.2) на листе «КВД» надо записать первое значение показания датчика Z10. Оно равно 0,7968 л/мин. Перевести это значение в систему СИ.
Согласно формуле (9) ln(0,67·A)=13,047, откуда 0,67·A=4,6⋅105.
Распишем все известные значения, которые входят в состав фильтрационного коэффициента А:
И получаем
Отсюда находится значение коэффициента эффективной проницаемости:
k = 1,34·10-12 м2 = 1340 мкм2 = 1,34 дарси.
Далее оценим значение проницаемости из углового коэффициента.
В нашем случае угловой коэффициент равен 0,489.
Тогда, согласно формуле (9)
Известно, что коэффициент пьезопроводности рассчитывается по формуле:
, (12)
где: μ – динамическая вязкость жидкости, Па*с;
β* – коэффициент объемного упругого расширения:
, (13)
где: m – пористость пласта, доли единицы;
βВОДЫ – коэффициент объемного упругого расширения воды;
βСКЕЛЕТА – коэффициент объемного упругого расширения скелета породы.
Примем μ=10-3 Па*с, m=0,2, βВОДЫ=3*10-10 Па-1, βСКЕЛЕТА=10-7 Па-1.
Подставляем значения в формулу (13) и далее в формулу (12). Получаем
