ВЕСТНИК ПНИПУ. ГЕОЛОГИЯ. НЕФТЕГАЗОВОЕ И ГОРНОЕ ДЕЛО BULLETIN OF PNRPU. GEOLOGY. OIL & GAS ENGINEERING & MINING
ISSN 2224-9923
Том / Volume 15 №21 2016
http://vestnik.pstu.ru/geo/
УДК 553.982.2
Статья / Article © ПНИПУ / PNRPU, 2016
ОПРЕДЕЛЕНИЕ ЗАБОЙНОГО ДАВЛЕНИЯ С ПОМОЩЬЮ МНОГОМЕРНЫХ СТАТИСТИЧЕСКИХ МОДЕЛЕЙ (НА ПРИМЕРЕ ПЛАСТА ТЛ-ББ ЮРЧУКСКОГО МЕСТОРОЖДЕНИЯ)
И.А. Черных
ООО «ЛУКОЙЛ-ПЕРМЬ» (614990, Россия, г. Пермь, ул. Ленина, 62)
DETERMINATION OF BOTTOMHOLE PRESSURE BY USING MULTIVARIATE STATISTICAL MODELS (ON EXAMPLE OF FORMATION TL-BB YURCHUKSKOIE FIELD)
I.A. Chernykh
LLC LUKOIL-PERM (62 Lenina st., Perm, 614990, Russian Federation)
Получена / Received: 06.10.2016. Принята / Accepted: 15.10.2016. Опубликована / Published: 02.12.2016
Ключевые слова: показатели эксплуатации скважин, забойное давление, механизированные скважины, необорудованные глубинными измерительными системами, методика определения забойного давления, статистический анализ, факторы, влияющие на величину забойного давления, корреляция, пошаговый линейный дискриминантный анализ, линейная дискриминантная функция, многомерная математическая модель, разнонаправленность влияния факторов, регрессионный анализ, множественная регрессия, работоспособность модели.
Проблема определения забойных давлений в механизированных добывающих скважинах, не оборудованных глубинными измерительными системами, является актуальной для многих нефтяных месторождений Пермского края. На практике, при отсутствии специальных приборов под скважинным насосом, забойное давление определяют путем пересчета динамического уровня. При таком подходе основную сложность представляет расчет плотности газожидкостной смеси, достоверность которого бывает низкой из-за влияния многочисленных осложняющих факторов. В настоящей работе предложен принципиально другой подход к определению забойного давления, рассмотренный на примере одной из добывающих скважин, оснащенных высокоточной глубинной измерительной системой, эксплуатирующей тульско-бобриковскую залежь Юрчукского месторождения. В качестве исходных данных приняты результаты непосредственного измерения забойного давления, а также ряд других показателей ее эксплуатации (дебиты нефти и жидкости, обводненность, динамический уровень, глубина спуска насоса под динамический уровень, давление на устье в затрубном пространстве). Выполненный на первом этапе анализ данных позволил сделать вывод, что забойное давление в период наблюдения изменялось, причем разнонаправленно: сначала постепенно снижалось, затем - увеличивалось. В этой связи исследование влияния показателей эксплуатации на величину забойного давления проведено для трех случаев: для всего периода наблюдений, а также отдельно для периодов его снижения и увеличения. Статистический анализ средних значений и плотностей распределений позволил выделить параметры, оказывающее влияние на величину забойного давления, и установить, что влияние это разнонаправленно. На заключительном этапе построены многомерные статистические модели, учитывающие разнонаправленное влияние показателей эксплуатации на забойное давление. Проверка «работоспособности» разработанных моделей выполнена на примере трех других скважин этого же объекта разработки. Она подтвердила целесообразность использования разработанных моделей для определения величин забойного давления по известным значениям показателей эксплуатации скважин и всего предложенного подхода в целом.
Key words:
well operation indicators, bottomhole pressure, mechanized wells, unequipped depth measuring systems, method of determining bottomhole pressure, statistical analysis, factors affecting the value of bottomhole pressure; correlation, stepwise linear discriminant analysis, linear discriminant function, multidimensional mathematical model, different directions of influence factors, regression analysis, multiple regression, models functionality.
The problem of determining BHP in production wells not equipped with depth measuring systems, it is relevant for many of the oil fields of Perm Krai. In practice, the absence of special devices for downhole pump bottomhole pressure is determined by converting the dynamic level. With this approach, the main difficulty is the calculation of the density of the gas-liquid mixture, the accuracy of which is low due to the influence of numerous complicating factors. In this paper the fundamentally different approach to the definition of bottomhole pressure, considered on one of the well equipped with high-precision depth measurement system, of Tl-Bb layer of Yurchukskoe field. The initial data are direct measurements of downhole pressure, as well as a number of other indicators of its operation (flow rates of oil and liquid, water content, dynamic level, the pump depth below the dynamic level, the pressure of the annulus). The first stage analysis of the data led to the conclusion that the bottomhole pressure during the observation period varied, with different directions: the first is gradually reduced, then - increased. In this regard, the study of influence of operating parameters on the value of BHP held for three cases: for the entire period of observation, as well as separately for the period of its decline and the increase. Statistical analysis of the averages and distribution densities possible to identify the parameters that influence the bottomhole pressure, and found that the effect is mixed. At the final stage the multidimensional statistical models that take into account the effect of multidirectional operating indicators BHP have been built. "Functionality" verification of developed models was made on the example of three other wells of the same development object. This verification confirmed the feasibility of using developed models for determining the values of bottomhole pressure from known values of indicators for well operation and all the proposed approach in general.
Черных Ирина Александровна - начальник отдела геофизических и гидродинамических исследований (моб. тел.: +007 919 483 84 30, e-mail: irina.chemykh@lp.lukoil.com).
Irina A. Chernykh - Head of the Department of Geophysical and Hydrodynamic Studies (mob. tel.: +007 919 483 84 30, e-mail: irina.chernykh@lp.lukoil.com).
Введение
Забойное давление является важнейшим параметром, характеризующим технологический режим работы каждой конкретной скважины. И в этой связи контроль за величиной забойного давления является основой мониторинга эксплуатации скважин [1-4].
С целью качественного контроля за эксплуатацией скважин добывающий фонд нефтяных месторождений ООО «ЛУКОЙЛ-ПЕРМЬ» оснащается различного рода глубинными измерительными системами. Так, на рассматриваемой в настоящей работе тульско-бобриковской залежи Юрчукского месторождения примерно 30 % скважин добывающего фонда оснащены такими системами. В остальных скважинах значение забойного давления определяется путем пересчета величины непосредственно измеряемого динамического уровня [5-12]. При этом важнейшую роль играет соответствие выбранной для пересчета методики реальным геолого-физическим условиям. На сегодняшний день в ООО «ЛУКОЙЛ-ПЕРМЬ» применяется ряд методик, позволяющих определять забойное давление. В них используются модели течения флюида в скважине, не учитывающие всех процессов, возникающих в стволе скважины при перемещении по нему многофазной жидкости. Аналогичные выводы получены авторами в работе [13].
В этой связи актуальной представляется разработка принципиально другой методики, основанной на математической обработке накопленных для каждого объекта данных и создании многомерных статистических моделей [14-16].
С целью разработки такой методики в данной работе использованы промысловые данные по скважинам, оснащенным глубинными измерительными системами, в которых многократно выполнены непосредственные замеры давлений на приеме насоса. Эти замеры гораздо точнее пересчитываются в забойные, приведенные к поверхности водонефтяного контакта, нежели динамический уровень [17]. Из общего перечня промысловых данных выбраны те параметры, которые, вероятно, оказывают влияние на забойное давление:
- глубина динамического уровня Нд;
- обводненность В;
- дебит жидкости ^^
- дебит нефти Qн;
- глубина погружения скважинного насоса под динамический уровень Нпогр;
- температура Т;
- давление на устье скважины Ру.
Исследование влияния параметров работы скважин на величину забойного давления
Исследование влияния параметров эксплуатации на величину забойного давления выполнено на примере скважины 882, характеризующейся максимальным количеством измерений (п = 98) за период с 16 декабря 2011 г. по 17 июля 2015 г. Изменение забойного давления Рзаб за время t отражено на рис. 1.
I замера
Рис. 1. Изменение значений Рзаб во времени
Видно, что в период с 16 декабря 2011 г. по 30 ноября 2013 г. происходит достаточно закономерное снижение значений Рзаб, далее -его закономерное повышение, но с меньшей интенсивностью. Исследуем изменение значений Рзаб и других показателей за весь период анализа - выборка 1; за время, когда значения давления снижаются, - выборка 2; за время, когда значения Рзаб повышаются, - выборка 3. Средние значения вышеприведенных показателей для вариантов 2 и 3 представлены в табл. 1.
Для оценки различий в средних значениях показателей вычислим значение критерия t по следующей формуле:
1v =■
1 - *2
V
1 1
(»1 - 1)2 + ( - 1)S
2 2
(1)
n2 - 2
где Х1, Х2 - соответственно средние значения показателей в выборках 2 и 3; (51 , 52 - дисперсии показателей в этих выборках. Различие в средних значениях считается статистически значимым, если tp > 4 Значения tt определяются в зависимости от количества сравниваемых данных и уровня значимости (а = 0,05). Из табл. 1 видно,
что средние значения статистически различаются
по Рзаб, В, бж, бн, Т, Ру.
Таблица 1 Сравнение статистических характеристик показателей
Показатель Средние значения Статистические оценки*
выборка 2 выборка 3 критерий Стьюдента tp критерий Пирсона х2
Рзаб, МПа 7,6 ± 0,3 7,1 ± 0,1 10,64453 0,000000 74,88653 0,000000
Нж, м 827,3 ± 146,1 835,8 ± 175,1 -0,26233 0,793631 1,179766 0,554392
В, % 6,3 ± 5,2 3,2 ± 1,5 3,95906 0,000144 14,59630 0,000679
0ж, мз/сут 23,5 ± 2,5 26,0 ± 0,7 -6,91203 0,000000 42,73306 0,000000
0н, т/сут 18,3 ± 2,3 20,8 ± 0,4 -7,51304 0,000000 48,06495 0,000000
Нпогр, м 727,7 ± 146,1 719,2 ± 175,1 -0,26233 0,793631 1,179766 0,554392
Т, °К 293,8 ± 0,5 294,1 ± 0,1 -4,32005 0,000038 17,00090 0,000203
Ру, МПа 2,3 ± 1,4 1,7 ± 0,8 2,38870 0,018864 6,731925 0,034529
P(Z), отн. ед. 0,672 ± 0,328 0,182 ± 0,088 9,992116 0,000000 58,21086 0,000000
Примечание: * в числителе - значение в знаменателе - уровень значимости а.
критерия;
Для более глубокого статистического анализа изучаемых показателей исследуем их распределения. Оптимальные величины интервалов значений показателей вычисляются по формуле Стерджесса:
X — X
д X _ max min
1 + 3,32 • lg N'
(2)
где Хтах, Хтт - соответственно максимальное и минимальное значения показателя; N -количество данных.
В каждом интервале определяются частости
N
P (X)_ Nk-,
^ ' Nq
(3)
где Р(Х) - частость в k-м интервале для группы (q = 2, соответствует второй выборке); Nk -число случаев содержания показателя Р(Х) в k-м интервале; Nq - объем выборки для выборок 2 и 3. При сравнении плотностей распределений показателей в изучаемых выборках применялась статистика Пирсона (х2), использовалась формула
1
X2 _ -1N2X
( M,
1M1 + M2
M 2 ^
N N.
(4)
2 У
где N1, N - соответственно количество значений показателей в данных выборок 2 и 3; Ы\, М2 - количество значений, попавших в заданный интервал, соответственно для двух изучаемых выборок; е - количество интервалов. Значения критерия х2 приведены в табл. 1. По этому критерию статистически отличаются следующие показатели: Рзаб, В, бж, бн, Т, Ру.
Исследование влияния параметров работы скважин на величину забойного давления с учетом его изменения во времени
Количественную оценку изменения значений Рзаб во времени для трех ранее описанных выборок произведем с помощью вычисления коэффициента корреляции г (табл. 2).
Видно, что за весь период t анализа (выборка 1) наблюдается статистически значимая обратная корреляция с Рзаб (г = -0,76), для выборки 2 эта корреляция несколько усиливается (г = -0,87), для третьей корреляция остается статистически значимой, но уже прямой (г = 0,65). В предположении того, что и другие изучаемые показатели могут изменяться во времени, по этим выборкам были вычислены значения г. Анализ показал, что подобные корреляции по направлению имеются у В и Ру. Для В корреляция по всем данным обратная (г = -0,31), корреляции по выборкам 2 и 3 равны соответственно -0,17 и 0,83. Для Ру общая корреляция также обратная при г = -0,36; корреляции по выборкам 2 и 3 равна соответственно -0,27 и 0,44. Вышеприведенные данные показывают, что в изучаемом временном диапазоне происходит значительная перестройка работы данной скважины. Для того чтобы убедиться, что перестройка не является случайным процессом, выдвинем следующую гипотезу: выделенные временные диапазоны работы скважины статистически отличаются по комплексу показателей. Для этого применим пошаговый линейный дискриминантный анализ. Возможности его использования для решения аналогичных задач приведены в работах [18-20]. Для разработки моделей разделения на варианты 2, 3 использовалась вся изучаемая выборка. Критерий Рзаб при разделении данных на варианты 2, 3 не применялся. Возможность построения линейной дискриминантной функции (ЛДФ) сводится к следующему. Если обозначить через Ху значения показателей с номером I в точке наблюдения и с номером у для второй выборки, то в результате может
Таблица 2
Корреляционная матрица
Р заб t Нд В Qж Qh Нпогр Т Ру
Р 1 заб 1,00 1,00 1,00 -0,76* -0,87* 0,65* 0,01 0,20 -0,01 0,37* 0,10 0,48* -0,69* -0,62* 0,54* -0,69* -0,53* 0,33* -0,00 -0,20 0,01 -0,47* -0,33 0,13 -0,29* 0,27 -0,18
t 1,00 1,00 1,00 -0,10 -0,31* 0,77* 0,74* 0,10 0,40 0,20 0,38* -0,36*
-0,40 -0,20 -0,17 0,83* 0,76* 0,80* 0,68* 0,39* 0,05 0,20 -0,27 0,44
Нд 1,00 1,00 1,00 0,05 0,15 -0,14 -0,24* -0,39 -0,39* -0,22* -0,38 -0,39 -0,99* -0,99* -0,99* 0,15 0,18 0,22 0,33* 0,28 0,49
В 1,00 1,00 1,00 -0,31* -0,17 0,61* -0,64* -0,60 0,04 -0,01 -0,11 0,10 -0,26 -0,15 0,19 -0,02 -0,08 -0,41*
0ж 1,00 1,00 1,00 0,92* 0,89* 0,81* -0,24* 0,40* 0,41* 0,30* 0,10 -0,10 -0,21* -0,09 -0,52*
Qh 1,00 1,00 1,00 0,24* 0,41* 0,40* 0,34* 0,15 -0,27* -0,19 -0,03 -0,36*
Нпогр 1,00 1,00 1,00 -0,15 -0,18 -0,16 -0,33* -0,28* -0,49
Т 1,00 1,00 1,00 0,01 0,12 0,00
Ру 1,00 1,00 1,00
Примечание: числитель - выборка 1; средняя строка корреляционные связи.
быть получена матрица Ж1 порядка т и п1 результатов наблюдений над этой выборкой:
W =
xu X1
X2 X О'
X_ X.
X1 X 2
X
Обозначим через Х* значения коэффициентов при переменных с номером I в точке с номером ] из третьей выборки, получим матрицу Ж2 порядка т хп2:
W =
X,, X,
12
X
1.2
X1 X1
X1.
XI XI . XI.
где т - число показателей; п1, п2, - объем выборок.
Для построения ЛДФ составляются матрицы центрированных сумм квадратов и смешанных произведений, по ним вычисляется выборочная матрица. Далее для определения коэффициентов
выборка 2; нижняя строка - выборка 3;* - значимые
линейной дискриминантной функции находится обратная выборочная ковариационная матрица -матрица С. Затем вычисляют граничное значение дискриминантных функций (Я0), которое делит выборку на два подмножества. Надежность классификации вычисляют с помощью критерия Пирсона х2. В результате реализации данного метода получена следующая линейная дискри-минантная функция:
7 = О,5250н + 0,933Т - 0,327Ру -
- 0,002Нпогр - 282,200 (5)
при Я = 0,679, х2 = 58,21, р = 0,000.
Среднее значение 7 для варианта 2 равно -0,898, для варианта 3 - 0,935. По данной функции были вычислены значения вероятности принадлежности к данным варианта 2 - Р(7). Соотношение между 7 и Р(7) приведено на рис. 2.
Видно, что при увеличении значений 7 величина Р(7) закономерно снижается. Среднее значение Р(7) для выборки 2 равно (0,672 ± 0,328); для выборки 3 (-0,182 ± 0,088). Отметим, что по данной функции верное распознавание составило 80 %, при этом все данные по варианту 3 распознаются верно. Оценку сравнений средних значений Р(7) и плотностей распределений
выборок 2 и 3 выполним по критериям ( и х (см. табл. 1).
Рис. 2. Зависимость Р(1 от 1
Выполненный статистический анализ средних значений и плотностей распределений показал, что Р(1), вычисленное по значениям Qн, Т, Ру, Нпогр, статистически различается. Здесь необходимо отметить, что влияние изучаемых показателей в ряде случаев даже разнонаправленно (см. табл. 2). Данный факт подтверждает сложность протекающих в скважине процессов при многофазной фильтрации и обусловливает необходимость разработки специального подхода к определению забойного давления, учитывающего разнонаправленность влияния многих факторов.
Разработка многомерных моделей для определения забойного давления
Для учета многообразия различных, в ряде случаев разнонаправленных, влияний изучаемых показателей на Рзаб построим несколько вариантов многомерных моделей с помощью регрессионного анализа. Расчет регрессионных коэффициентов в разрабатываемой модели выполним при помощи метода наименьших квадратов. Под регрессионным анализом понимается статистический метод исследования зависимостей между зависимой переменной У и одной или несколькими независимыми переменными Х1, Х2, Хр. Зависимый признак в регрессионном анализе называется результирующим, независимый - факторным. Обычно на зависимую переменную действуют сразу несколько факторов. Совокупное влияние всех независимых факторов на зависимую переменную учитывается благодаря множественной регрессии.
В общем случае множественную регрессию оценивают параметры линейного уравнения вида
У = а + Ь1Х1 + Ь2Х2 +... + ЬрХр.
В данном уравнении регрессионные коэффициенты (Ь-коэффициенты) представляют независимые вклады каждой независимой переменной в предсказание зависимой переменной. Линия регрессии выражает наилучшее предсказание зависимой переменной У по независимым переменным Х. В нашем случае в качестве зависимого признака выступает Рзаб, а в качестве независимых факторов - значения
Ндч QЖ, НПOГр, Тч Ру.
Для решения задачи регрессионного анализа методом наименьших квадратов вводится понятие функции невязки:
1 M 2
)=2Z(y - Y).
(6)
k=1
Условие минимума функции невязки:
' М М N
Е -у-- =ЕЕ ьа > + ЬоМ,
d ст(Ь)-0 i = 0...N
i=1
M
i=1 j=1
M N
M
Z yx k =ZZ'bjX, jXk+ Mb0 Zx ,k i=1 i=1 j=1 i=1
k = 1...N.
(7)
Полученная система является системой N + 1) линейных уравнений с N + 1) неизвестными Ь0...Ь^
Если представить свободные члены левой части уравнений матрицей
M
B =
Z Уг
i=1 M
Z УгХг,1
Z У г
iXi, N
а коэффициенты при неизвестных в правой части матрицей
A =
M Z Хм i=1 Z X,2 i=1 ... Z X, N i=1
M Z xi,1 i=1 M Z xi,1 xi,1 i=1 M Z Xi,2Xi,1 i=1 M .. Z Xi,NXi,1 i=1
M Z Xi,2 i=1 M Z Xi,1Xi,2 i=1 M Z Xi,2Xi,2 i=1 Z .. Z Xi,NXi,2 i=1
M Z Xi,N i=1 M Z Xi,1 Xi,N i=1 M Z Xi,2Xi,N i=1 M . Z Xi, NXi, N i=1
i=1
то получаем матричное уравнение А х X = В, которое легко решается методом Гаусса. Полученная матрица содержит коэффициенты уравнения линии регрессии:
X =
Построение многомерных моделей произведем по трем выборкам. Первоначально построим многомерную модель по скважине 882 с использованием 50 % данных (вариант 1). Данные сформированы в обучающую выборку следующим образом: первыми используются значения данных от 16.12.2011 г. (п = 1); данные, полученные 17.12.2011 г. (п = 2), в построении многомерной модели не применяются. Они берутся при определении «работоспособности» построенной модели. Таким образом, по данным, имеющим нечетный номер, будут построены модели, а по данным, имеющим четный номер, будет определена их «работоспособность». В результате будет получена первая многомерная модель и выполнена оценка ее «работоспособности». При сопоставлении значений Рзаб и значений, полученных по построенной модели рабм, находим варианты 1 и 2. Модель, построенная с выполнением данного условия, имеет следующий вид:
Рзабм = 71,594 - 1,67465 & -- 0,22745Г + 0,04276Ру + 0,44088В + (8) + 2,13748 0н + 0,00016 ЯП0Гр
при Я = 0,862, р < 0,0000.
По данной модели были вычислены значения Р1абм, которые использованы при сравнении со значениями Рзаб по вариантам 1, 2 (рис. 3).
Далее была построена модель с учетом классов, выделенных с помощью метода ЛДФ (без учета деления скважин на две четные и нечетные выборки). По этим данным была построена модель, имеющая следующий вид:
Рзабм = 66,186 - 0,07541 Q.M -- 0,19377 Т + 0,06005PV
(9)
при Я = 0,722, р < 0,0000.
С помощью данной формулы вычислены значения по варианту 3 (см. рис. 3). Далее по этой модели были вычислены значения Рз;;бм для экзаменационной выборки. Сопоставление
Далее была построена модель с учетом классов, выделенных с помощью метода ЛДФ, и с учетом деления скважин на два класса по четности и нечетности выборок (вариант 5). С учетом этих ограничений была построена модель, имеющая следующий вид:
Рзабм = 3,736 + 0,208301 Q,K + + 0,000192Нд - 0,108676Qh
(10)
при Я = 0,799, р < 0,00000.
Затем по этой модели были вычислены значения Рз^бм для экзаменационной выборки.
Сопоставление
с ра
Рзаб С P
для экзаменационной выборки
приведено на рис. 3 (вариант 4).
1забм для экзаменационной выборки приведено на рис. 3 (вариант 6).
По данным моделям были определены значения Рз^м и Рз3бм и выполнено сопоставление со значениями Рзаб путем построения полей корреляции между ними (см. рис. 3).
Далее выполнена проверка работоспособности трех построенных моделей для других скважин этого же объекта разработки (скважины 206, 209 и 216) в двух вариантах. В первом случае для расчетов используется 1-я модель без учета деления с помощью ЛДФ (вариант 7), во втором случае расчеты выполнены с учетом деления на классы по ЛДФ по моделям 2 и 3 (вариант 8). Здесь для экзамена использовались только те данные, у которых значение Рзаб находилось в диапазоне 6,6-8,2 МПа. Это вызвано тем, что оценивается работоспособность построенных моделей только в интервале интерполяции. Выполним анализ созданных корреляционных полей. Корреляционное поле, построенное по данным обучающей выборки (вариант 1), даже на визуальном уровне состоит из двух подполей. Аналогичное поле наблюдается по данным экзаменационной выборки. Это свидетельствует о том, что данные достаточно близки. Корреляционное поле, построенное с учетом деления на классы с помощью ЛДФ (вариант 3), показывает, что визуального деления на классы уже не видно. Аналогичное поле наблюдается по данным экзаменационной выборки (вариант 4). В случае, когда в обучении используются не только значения Рзабм, вычисленные по моделям 2, 3, но и данные по обучению и экзамену и выполняется их сравнение с Рзаб (вариант 5, 6), то наблюдается разделение их на зоны (см. рис. 3). Для варианта 5 корреляционное поле располагается в верхнем правом углу графика, для варианта 6 - в нижнем левом углу. При этом необходимо отметить, что для варианта 5 г = 0,77, а для варианта 6 - 0,65. Эти связи статистически значимы. Выполним проверку работоспособности построенных моделей для данных скважин 206,
8,6 8,4 8,2 8,0 7,8 7,6 7,4 7,2 7,0 6,8 6,6 6,4 62
* • К
6,50 6,84 7,14 7,44 7,78 8,14 6,66 6,99 7,29 7,60 7,99 8,40 Ра6, МПа Вариант 1
шт
6,50 6,84 7,14 7,44 7,78 8,14 6,66 6,99 7,29 7,60 7,99 8,40 Л*, МПа Вариант 3
•Л**«
6,50 6,84 7,14 7,44 7,78 8,14 6,66 6,99 7,29 7,60 7,99 8,40 Л- МПа Вариант 5
• • * •
6,50 6,84 7,14 7,44 7,78 8,14 6,66 6,99 7,29 7,60 7,99 8,40 МПа Вариант 7
' 6,50 6,84 7,14 7,44 7,78 8,14 6,66 6,99 7,29 7,60 7,99 МПа Вариант 2
8,40
8,6 8 4 8,2 8,0 7,8 7,6 74 72 7,0 6,8 6,6 6,4 62
ъ
6,50 6,84 7,14 7,44 7,78 6,66 6,99 7,29 7,60 7.
Л.,, МПа Вариант 4
8,14 99
8,40
'S
о,'
8,6 84 8,2 8,0 7,8 76 7,4 7,2 7,0 6,8 6,6 6,4 6,2
8,6 8 4 8,2 8,0 7,8 7,6 74 72 7,0 6.8 6,6 64 62
,50 6,84 7,14 7,44 7,78 6,66 6,99 7,29 7,60 " К5, МПа Вариант 6
8,14 ,99
8,40
•
* • • ф
• *
*
i,50 6,84 7,14 7,44 7,78 6,66 6,99 7,29 7,60 7.
МПа Вариант 8
8,14 99
8,40
Рис. 3. Поля корреляции между Рзаб и Рзабм
209 и 216 (варианты 7, 8). В случае, когда используем модель 1, связи между Рзаб и Рзабм не наблюдается (г = -0,12). При построении модели с использованием ЛДФ и вычислении значений Рзабм по моделям 2 и 3 связь между Рзаб и Рзабм является прямой и статистически значимой (г = 0,87). Таким образом, можно констатировать, что использование изменения значений Рзаб во времени позволяет более точно прогнозировать значения Рзабм по косвенным признакам.
Выводы
По результатам выполненных исследований можно сделать следующие основные выводы:
1. Забойное давление в механизированных скважинах, не оборудованных глубинными измерительными системами, определяется пересчетом динамического уровня по известным методикам со значительными неточностями.
2. Математическая обработка промысловых данных (забойное давление, дебиты нефти и жидкости, обводненность и др.) по одной из скважин Юрчукского месторождения, где были проведены непосредственные измерения, позволила разработать методику построения многомерных моделей, которая может быть исполь-
зована и на других месторождениях для определения забойного давления по известным значениям других показателей эксплуатации скважин.
3. Достоверность получаемых результатов с использованием построенных многомерных моделей проверена на примере трех других скважин данного объекта разработки.
Список литературы
1. Лекомцев А.В., Мордвинов В.А., Турбаков М.С. Эксплуатация добывающих скважин электроцентробежными насосами на нефтяных месторождениях Верхнего Прикамья // Нефтяное хозяйство. - 2010. -№ 10. - С. 144-145.
2. Лекомцев А. В., Мордвинов В. А., Турбаков М. С. Оценка забойных давлений в добывающих скважинах Шершневского месторождения // Нефтяное хозяйство. -2011. - № 10. - С. 30-31.
3. Муравьев И.М., Халиков Г. А., Юрин И.Я. Приведение давлений, замеренных на глубинах выше глубины давления насыщения, к забойным условиям // Нефтяное хозяйство. - 1960. - № 8. - С. 26-29.
4. Carvalho P.M., Podio A.L., Sepehrnoori K. An electrical submersible jet pump for gassy oil wells // Journal of Petroleum Technology. - 1999. - May. - P. 34-35. DOI: 10.2118/0599-0034-JPT.
5. Временное методическое руководство по расчету режимов скважин, эксплуатируемых глубинными насосами (ЭЦН и ШГН). - Уфа: БашНИПИнефть, 1976. - 114 с.
6. Грон В. Г., Мищенко И. Т. Определение забойного давления в добывающих скважинах, оборудованных установками погружного насоса: учеб. пособие. - М.: Изд-во ГАНГ, 1993. - 128 с.
7. Мищенко И.Т. Расчеты при добыче нефти и газа. -М.: Нефть и газ, 2008. - 295 с.
8. Лекомцев А.В., Мордвинов В.А. К оценке забойных давлений при эксплуатации скважин электроцентробежными насосами // Научные исследования и инновации. -2011. - Т. 5, № 4. - С. 29-32.
9. Дроздов А.Н. Технология и техника добычи нефти погружными насосами в осложненных условиях: учеб. пособие для вузов. - М.: Изд-во РГУ нефти и газа им. И.М. Губкина, 2008. - 616 с.
10. Lea J.F., Minissale J.D. Beam pumps surpass ESP efficiency // Oil and Gas Journal. - 1992. - May 18. - P. 72.
11. McCoy J.N., Podio A.L., Huddleston K.L. Acoustic determination of producing bottomhole pressure // SPE Fomation Evaluation. - 1985. - August. - P. 617-621. DOI: 10.2118/14254-PA.
12. Podio A.L., McCoy J.N., Becker D. Integrated well performance and analysis // SPE Computer Applications. -1992. - June. - P. 43-48. DOI: 10.2118/24060-PA.
13. Лекомцев А.В., Мордвинов В.А. Определение давления у приема электроцентробежных насосов по данным исследований скважин // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. -2012. - № 4. - С. 84-90.
14. Гмурман В.Е. Теория вероятностей и математическая статистика: учеб. пособие для вузов. - 4-е изд., испр. и доп. - М.: Высшая школа, 1972. - 368 с.
15. Дементьев Л.Ф. Статистические методы обработки и анализа промыслово-геологических данных. -М.: Недра, 1966. - 206 с.
16. Многомерный анализ данных методами прикладной статистики / С.С. Барковский, В.М. Захаров, А.М. Лукашов, А.Р. Нурутдинова, С.В. Шалагин. -Казань, 2010. - 126 с.
17. Davaatseren B., Golovko A.K., Tuya M. A study of the ozonolysis and mechanochemical treatment on the properties of the high paraffinic Tamsagbulag Oil (Mongolia) // Scientific reports of MAS. - Ulan-bator, 2006. - Vol. 182, № 4. - P. 59-68.
18. Путилов И.С., Галкин В.И. Применение вероятностного статистического анализа для изучения фациальной зональности турне-фаменского карбонатного комплекса Сибирского месторождения // Нефтяное хозяйство. - 2007. - № 9. - С. 112-114.
19. Путилов И. С. Разработка технологий комплексного изучения геологического строения и размещения месторождений нефти и газа. - Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2014. - 285 с.
20. Галкин В.И., Козлова И.А. Разработка вероятностно-статистических регионально-зональных моделей прогноза нефтегазоносности по данным геохимических исследований верхнедевонских карбонатных отложений // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2016. - № 6. -С. 40-45.
References
1. Lekomtsev A.V., Mordvinov V.A., Turbakov M.S. Ekspluatatsiia dobyvaiushchikh skvazhin elektrotsentrobezhnymi nasosami na neftianykh mestorozhdeniiakh Verkhnego Prikam'ia [Submersible pumps operation in wells at oilfields in Verkhnee Prikamje]. Neftianoe khoziaistvo, 2010, no.10, pp.144-145.
2. Lekomtsev A.V., Mordvinov V.A., Turbakov M.S. Otsenka zaboinykh davlenii v dobyvaiushchikh skvazhinakh Shershnevskogo mestorozhdeniia [Estimation of bottom-hole
pressure in producing wells of Shershnevkoe oilfield]. Neftianoe khoziaistvo, 2011, no.10, pp.30-31.
3. Murav'ev I.M., Khalikov G.A., Iurin I.Ia. Privedenie davlenii, zamerennykh na glubinakh vyshe glubiny davleniia nasyshcheniia, k zaboinym usloviiam [Bringing pressures measured at depths above the saturation pressure's depths to downhole conditions]. Neftianoe khoziaistvo, 1960, no.8, pp.26-29.
4. Carvalho P.M., Podio A.L., Sepehrnoori K. An electrical submersible jet pump for gassy oil wells. Journal of Petroleum Technology, 1999, may, p.34-35. DOI: 10.2118/0599-0034-JPT.
5. Vremennoe metodicheskoe rukovodstvo po raschetu rezhimov skvazhin, ekspluatiruemykh glubinnymi nasosami (ETsN i ShGN) [Temporary methodological guidance on the modes calculation in wells exploited with borehole pumps (ESP and SRP)]. Ufa: BashNIPIneft', 1976, 114 p.
6. Gron V.G., Mishchenko I.T. Opredelenie zaboinogo davleniia v dobyvaiushchikh skvazhinakh, oborudovannykh ustanovkami pogruzhnogo nasosa [Determination of the bottomhole pressure in production wells equipped with submersible pump units]: ucheb. posobie. Moscow: Izd-vo GANG, 1993, 128 p.
7. Mishchenko I.T. Raschety pri dobyche nefti i gaza [Calculations during oil and gas production]. Moscow: Neft' i gaz, 2008, 295 p.
8. Lekomtsev A.V., Mordvinov V.A. K otsenke zaboinykh davlenii pri ekspluatatsii skvazhin elektrotsentrobezhnymi nasosami [On estimation of bottomhole pressures in wells operation with electric centrifugal pumps]. Nauchnye issledovaniia i innovatsii, 2011, vol.5, no.4, pp.29-32.
9. Drozdov A.N. Tekhnologiia i tekhnika dobychi nefti pogruzhnymi nasosami v oslozhnennykh usloviiakh [Technology and technics of oil production using submersible pumps under complicated conditions]: ucheb. posobie dlia vuzov. Moscow: RGU nefti i gaza im. I.M. Gubkina, 2008, 616 p.
10. Lea J.F., Minissale J.D. Beam pumps surpass ESP efficiency. Oil and Gas Journal, 1992, May 18, p.72.
11. McCoy J.N., Podio A.L., Huddleston K.L. Acoustic determination of producing bottomhole pressure. SPE Fomation Evaluation, 1985, august, p.617-621. DOI: 10.2118/14254-PA.
12. Podio A.L., McCoy J.N., Becker D. Integrated well performance and analysis. SPE Computer Applications, 1992, June, p.43-48. DOI: 10.2118/24060-PA.
13. Lekomtsev A.V., Mordvinov V.A. Opredelenie davleniia u priema elektrotsentrobezhnykh nasosov po dannym issledovanii skvazhin [Determination of pressure in receiving electrical pumps according to the data of the
researched well]. Vestnik Permskogo natsional'nogo issledo-vatel'skogo politekhnicheskogo universiteta. Geologiia. Neftegazovoe i gornoe delo, 2012, no.4, pp.84-90.
14. Gmurman V.E. Teoriia veroiatnostei i matemati-cheskaia statistika [Theory of probability and mathematical statistics]: ucheb. posobie dlia vuzov. 4 izd., ispr. i dop. Moscow: Vysshaia shkola, 1972, 368 p.
15. Dement'ev L.F. Statisticheskie metody obrabotki i analiza promyslovo-geologicheskikh dannykh [Statistical methods for processing and analysis of the field and geological data]. Moscow: Nedra, 1966, 206 p.
16. Barkovskii S.S., Zakharov V.M., Lukashov A.M., Nurutdinova A.R., Shalagin S.V. Mnogomernyi analiz dannykh metodami prikladnoi statistiki [Multivariate data analysis by methods of applied statistics]. Kazan', 2010, 126 p.
17. Davaatseren B., Golovko A.K., Tuya M. A study of the ozonolysis and mechanochemical treatment on the properties of the high paraffinic Tamsagbulag Oil (Mongolia). Scientific reports of MAS. Ulan-bator, 2006, vol.182, no.4, pp.59-68.
18. Putilov I.S., Galkin V.I. Primenenie veroiatnostnogo statisticheskogo analiza dlia izucheniia fatsial'noi zonal'nosti turne-famenskogo karbonatnogo kompleksa Sibirskogo mestorozhdeniia [The results of statistical analysis for study fades characterization of T-Fm stage of Sibirskoe oilfield]. Neftianoe khoziaistvo, 2007, no.9, pp.112-114.
19. Putilov I.S. Razrabotka tekhnologii kompleksnogo izucheniia geologicheskogo stroeniia i razmeshcheniia mestorozhdenii nefti i gaza [Development of technologies for a comprehensive study of the geological structure and location of oil and gas fields]. Perm': Izd-vo Perm. nats. issled. politekhn. un-ta, 2014, 285 p.
20. Galkin V.I., Kozlova I.A. Razrabotka veroiatnostno-statisticheskikh regional'no-zonal'nykh modelei prognoza neftegazonosnosti po dannym geokhimicheskikh issledovanii verkhnedevonskikh karbonatnykh otlozhenii [Development of probabilistic-statistical regional-zoning models of oil and gas potential prediction based on the data of geochemical studies of the Upper Devonian carbonate deposits]. Geologiia, geofizika i razrabotka neftianykh i gazovykh mestorozhdenii, 2016, no.6, pp.40-45.
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:
Черных И.А. Определение забойного давления с помощью многомерных статистических моделей (на примере пласта ТЛ-ББ Юрчукского месторождения) // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. - 2016. - Т.15, №21. - С.320-328. DOI: 10.15593/2224-9923/2016.21.3
Please cite this article in English as:
Chernykh I.A. Determination of bottomhole pressure by using multivariate statistical models (on example of formation TL-BB
Yurchukskoie field). Bulletin ofPNRPU. Geology. Oil & Gas Engineering & Mining, 2016, vol.15, no.21, рр. 320-328. DOI: 10.15593/2224-9923/2016.21.3