Статья опубликована в рамках: Научного журнала «Студенческий» № 22(234)
Рубрика журнала: Информационные технологии
Скачать книгу(-и): скачать журнал часть 1, скачать журнал часть 2, скачать журнал часть 3, скачать журнал часть 4, скачать журнал часть 5, скачать журнал часть 6, скачать журнал часть 7, скачать журнал часть 8, скачать журнал часть 9
СОЗДАНИЕ МЕТОДИКИ ПО ПРОГНОЗИРОВАНИЮ И ВЫЯВЛЕНИЮ АВАРИИ ГАЗОНЕФТЕВОДОПРОЯВЛЕНИЯ С ПОМОЩЬЮ ЦИФРОВЫХ ТЕХНОЛОГИЙ
CREATION OF A METHODOLOGY FOR PREDICTING AND DETECTING AN ACCIDENT OF GAS AND OIL PRODUCTION WITH THE HELP OF DIGITAL TECHNOLOGIES
Anastasia Sharafutdinova
master's student, Department of Automated Control Systems, Ufa state aviation technical University,
Russia, Ufa
Elena Dron
scientific supervisor, Candidate of Sciences in Technical, associate professor, Ufa state aviation technical University,
Russia, Ufa
АННОТАЦИЯ
Статья описывает разработку методики прогнозирования и выявления аварий газонефтеводопроявлений с использованием цифровых технологий, включая метод наименьших квадратов. Автор исследования представляет результаты анализа данных, полученных с помощью датчиков и других устройств, установленных на объектах нефтегазовой отрасли. Результаты исследования позволяют улучшить систему мониторинга и предотвращения аварий на объектах нефтегазовой отрасли. Статья будет полезна для специалистов в области нефтегазовой промышленности, а также для тех, кто интересуется применением цифровых технологий в промышленности.
ABSTRACT
The article describes the development of a methodology for predicting and detecting accidents of gas and oil occurrences using digital technologies, including the least squares method. The author of the study presents the results of the analysis of data obtained using sensors and other devices installed at oil and gas industry facilities. The results of the study make it possible to improve the system of monitoring and prevention of accidents at oil and gas industry facilities. The article will be useful for specialists in the oil and gas industry, as well as for those who are interested in the application of digital technologies in industry.
Ключевые слова: линейное регрессионное уравнение, газонефтеводопроявление, влияющие факторы, KPI, поглощение бурового раствора.
Keywords: linear regression equation, oil and gas occurrence, influencing factors, KPI, drilling mud absorption.
Введение
Нефтегазовая отрасль является одной из наиболее важных и перспективных отраслей промышленности. Однако, в связи с особенностями добычи и транспортировки нефти и газа, этот сектор также является одним из самых опасных и рискованных. Аварии на объектах нефтегазовой отрасли могут привести к серьезным последствиям для окружающей среды и здоровья людей. В данной статье рассматривается разработка методики прогнозирования и выявления аварий газонефтеводопроявлений с использованием цифровых технологий, которая может помочь предотвратить подобные происшествия.
Основная часть
Модель линейной множественной регрессии (Multiple Linear Regression)- это практичная статистическая модель для оценки связей между непрерывной зависимой переменной и переменными-предикторами.
Множественной называют линейную регрессию, в модели которой число независимых переменных две или более.
Уравнение множественной линейной регрессии имеет вид:
Yрасчетн = b0 + b1x1 + b2x2 + … + bnxn (1)
Y – исследуемый параметр;
x0 , x1, x2 – исследуемые факторы;
b0 , b1, b2 – неизвестные коэффициенты.
Рассмотрим построение линейной множественной регрессии на примере аварии газонефтеводопроявления (ГНВП). Будем считать, что исследуемый параметр Y – это проявление аварии ГНВП, а параметры X1, X2, X3, Xn – это влияющие факторы на проявление аварии.
Факторами X1, X2, X3, Xn будем считать следующие параметры:
- Общий объем бурового раствора Vобщ, л
- Расход бурового раствора Qбур
- Давление на выходе Рвых, атм
- Нагрузка на долото Gо, кН/мм
- Момент кручения буровой колонны, кН/м
- Механическая скорость (υ м), м/ч
- Объемное газосодержание в промывочной жидкости βг, % .
В таблице 1 приведены значения параметров по скважинам и с оценкой экспертов вероятности возникновения аварии ГНВП.
Таблица 1.
Таблица значений факторов
№ скв |
Авария ГНВП |
Общий объем бурового раствора Vобщ, л |
Расход БР Qбур, л/м.п. |
Давление на выходе Рвых, атм |
Нагрузка на долото Gо, кН/мм |
Момент кручения бур. колонны Т, кН/м |
Механическая скорость бурения (υм), м/ч |
Объемное газосодерж. в ПЖ βг, % |
Плотность бурового раствора, кг/м3 |
105 |
0,21 |
50 |
2,1 |
25 |
0,3 |
1262 |
24 |
5 |
1554 |
106 |
0,12 |
72 |
1,1 |
20 |
0,6 |
1596 |
48 |
12,3 |
1585 |
107 |
0,19 |
96 |
3 |
30 |
0,6 |
1451 |
46 |
19,2 |
1637 |
108 |
0,29 |
74 |
1,7 |
26 |
0,5 |
1120 |
40 |
18 |
1549 |
109 |
0,12 |
112 |
2,6 |
23 |
0,5 |
1471 |
28 |
10,4 |
1513 |
110 |
0,03 |
64 |
2,8 |
21 |
0,4 |
1746 |
21 |
17,9 |
1750 |
111 |
0,11 |
54 |
1,5 |
29 |
0,6 |
1630 |
39 |
7 |
1381 |
112 |
0,02 |
115 |
2,9 |
26 |
0,4 |
1738 |
23 |
15,5 |
1572 |
113 |
0,22 |
110 |
2,6 |
21 |
0,6 |
1340 |
22 |
11,7 |
1248 |
114 |
0,03 |
77 |
1 |
20 |
0,6 |
1478 |
39 |
21,6 |
1486 |
115 |
0,09 |
120 |
1,2 |
27 |
0,5 |
1321 |
18 |
14,1 |
1634 |
116 |
0,02 |
92 |
1 |
24 |
0,4 |
1434 |
42 |
9,9 |
1341 |
117 |
0,02 |
61 |
1,9 |
24 |
0,6 |
1589 |
43 |
16,7 |
1641 |
118 |
0,14 |
118 |
1,4 |
29 |
0,3 |
1277 |
15 |
14,6 |
1270 |
119 |
0,14 |
94 |
2 |
31 |
0,3 |
1344 |
28 |
6 |
1594 |
120 |
0,16 |
69 |
1,7 |
30 |
0,6 |
1766 |
25 |
19,1 |
1149 |
105 |
0,56 |
72 |
3 |
33 |
0,2 |
896 |
37 |
9,4 |
1197 |
106 |
0,49 |
113 |
1,8 |
27 |
0,4 |
862 |
73 |
22,3 |
1157 |
107 |
0,59 |
162 |
4,2 |
39,3 |
0,3 |
1088 |
70 |
30,9 |
1162 |
108 |
0,73 |
107 |
2,6 |
34,8 |
0,3 |
661 |
67 |
34 |
1146 |
109 |
0,76 |
175 |
3,9 |
31,1 |
0,4 |
765 |
45 |
18,2 |
1074 |
110 |
0,62 |
93 |
4,4 |
28,4 |
0,3 |
1083 |
35 |
32,8 |
1225 |
111 |
0,48 |
80 |
2,3 |
39,2 |
0,4 |
929 |
60 |
13,9 |
1105 |
112 |
0,51 |
193 |
4,2 |
34,1 |
0,3 |
1199 |
36 |
26 |
1179 |
113 |
0,68 |
167 |
3,7 |
27,3 |
0,4 |
871 |
35 |
22,7 |
949 |
114 |
0,47 |
126 |
1,6 |
26,2 |
0,4 |
902 |
67 |
39,1 |
1144 |
115 |
0,56 |
194 |
1,8 |
35,4 |
0,4 |
766 |
30 |
25,5 |
1144 |
116 |
0,39 |
129 |
1,5 |
31,4 |
0,2 |
946 |
70 |
17,5 |
1033 |
117 |
0,61 |
100 |
2,7 |
31,9 |
0,4 |
826 |
75 |
27,7 |
1264 |
118 |
0,67 |
197 |
2 |
38 |
0,2 |
639 |
25 |
28 |
953 |
119 |
0,49 |
138 |
2,9 |
40,6 |
0,2 |
995 |
47 |
10,3 |
1148 |
120 |
0,46 |
115 |
2,7 |
40,2 |
0,3 |
1060 |
38 |
31,3 |
885 |
До того как приступить к шагам построения уравнения Y= F(..Xn), необходимо убедиться, что каждый фактор Х существенно влияет на показатель Y. Если в процессе установки зависимости окажется, что какой-то из факторов не оказывает влияния на результирующий результат, то данный фактор включаться в итоговое уравнение не будет. Необходимо провести расчет всех возможных комбинаций между Y и факторами X.
В таблице 2 приведены зависимости между Y и факторами X.
Будем считаем, что связь между параметрами есть, если коэффициент корреляции превышает |0,5|. Все рассматриваемые параметры, кроме параметра «Механическая скорость бурения» (0,46) имеют связь с возникновением аварии ГНВП. Данный параметр не будет учитываться при построении регрессионного уравнения.
Таблица 2.
Зависимости между Y и факторами X.
|
Авария ГНВП |
Vобщ, л |
Qбур, л/м.п. |
Рвых, атм |
Gо, кН/мм |
Т, кН/м |
υм, м/ч |
βг, % |
ρ, кг/м3 |
Авария ГНВП |
- |
0,57 |
0,62 |
0,64 |
-0,56 |
-0,90 |
0,46 |
0,57 |
-0,67 |
Vобщ, л |
|
- |
0,39 |
0,5 |
-0,33 |
-0,57 |
0,07 |
0,31 |
-0,61 |
Qбур, л/м.п. |
|
|
- |
0,41 |
-0,47 |
-0,37 |
0,21 |
0,3 |
-0,47 |
Рвых, атм |
|
|
|
- |
-0,56 |
-0,56 |
0,29 |
0,4 |
-0,64 |
Gо, кН/мм |
|
|
|
|
- |
0,54 |
-0,27 |
-0,27 |
0,55 |
Т, кН/м |
|
|
|
|
|
- |
-0,46 |
-0,47 |
0,62 |
υм, м/ч |
|
|
|
|
|
|
- |
0,24 |
-0,33 |
βг, % |
|
|
|
|
|
|
|
- |
-0,4 |
ρ, кг/м3 |
|
|
|
|
|
|
|
|
- |
Некоторые факторы могут оказаться взаимосвязаны друг от друга. В таком случае нет необходимости включать их при построении модели. Если коэффициент корреляции между факторами меньше 0,5, то при построении модели их следует учитывать оба, так как они никак не связаны между собой. Если коэффициент корреляции между факторами больше 0,5, то необходимо провести выбор в пользу какого-то одного фактора, который будет использоваться в модели. В данном случае выбор будет в пользу того фактора, у которого связь с исследуемым параметром (аварией) выше.
Например, общий объем бурового раствора имеет коэффициент корреляции с расходом бурового раствора 0,39, что означает, что параметры между собой связь не имеют и в расчет регрессионного уравнения аварии ГНВП следует использовать оба параметра. Однако общий объем бурового раствора имеет связь с параметром плотности бурового раствора (-0,61), что говорит о том, что в регрессионное уравнение включать оба параметра не имеет смысла. Один из рассматриваемых параметров, который влияет меньше на вероятность возникновения аварии ГНВП, следует исключить, в данном случае плотность бурового раствора (-0,67) имеет более высокую связь с исследуемым фактором, чем общий объем бурового раствора (0,57). Получим результирующую таблицу (таблица 3).
Таблица 3.
Результирующая таблица параметров
Авария ГНВП |
Qбур, л/м.п. |
Т, кН/м |
βг, % |
|
Авария ГНВП |
- |
0,62 |
-0,9 |
0,57 |
Qбур, л/м.п. |
- |
-0,37 |
0,3 |
|
Т, кН/м |
- |
-0,47 |
||
βг, % |
- |
Произведем расчет множественной регрессии.по формуле (1). В рассматриваемом примере необходимо найти 4 неизвестных параметра b. Для этого необходимо рассчитать статистику для ряда с применением метода наименьших квадратов, чтобы вычислить прямую линию, которая наилучшим образом аппроксимирует имеющиеся данные и затем возвращает массив, который описывает полученную прямую. Массив матрицы значений, в которой количество столбцов равняется числу неизвестных коэффициентов b регрессионного уравнения:
Ниже в таблице 4 представлены неизвестные коэффициенты b:
Таблица 4.
Таблица коэффициентов исследуемых параметров
b3 |
b2 |
b1 |
b0 |
0,0031 |
-0,0005 |
0,0861 |
0,7263 |
Поскольку коэффициенты, вносимые в исследуемые факторы вес найдены, получим линейное регрессионное уравнение для выявления риска ГНВП (KPI):
KPI = 0,7263 + 0,0861× Qбур + (-0,0005) × Т + 0,0031 ×βг (2)
Qбур - расход бурового раствора;
Т - момент кручения буровой колонны
βг - Объемное газосодержание в промывочной жидкости.
На рисунке 1 изображен график сопоставимости вероятности возникновения ГНВП с вычисленным KPI по параметрам. Для прогноза возникновения аварии система будет вычислять KPI и по этому значению выдавать прогноз возникновения аварии.
Рисунок 1. График сопоставимости вероятности возникновения ГНВП с вычисленным KPI
При достижении KPI выше 0,5 следует принимать меры по предотвращению возникновения аварии.
ГНВП часто сопровождается поглощением БР, причем в этот момент буровой раствор заменяется газом. Для минимизации эффекта поглощения бурового раствора рекомендуется использовать разнофракционную кальматирующую пачку. Для ликвидации возникновения аварии необходимо закрыть дроссельную задвижку на устье скважины. Если эта задвижка не будет закрыта, то значительное количество бурового раствора продолжит выходить, и скважина будет практически полностью заполнена газом. Поскольку газ имеет высокую степень подвижности, что влечет к возникновению опасной ситуации.
Заключение
Разработка методики прогнозирования и выявления аварий газонефтеводопроявлений с использованием цифровых технологий является важным шагом в повышении безопасности нефтегазовой отрасли. Главным преимуществом методики является необходимость отслеживать только один параметр, а не несколько. Поскольку анализ каждого параметра в отдельности может свидетельствовать о проявлении сразу нескольких возможных аварий, однако при анализе параметров в совокупности (через расчет KPI) система достаточно точно диагностирует вид аварии на начальном этапе, что является ее главным преимуществом. Расчет KPI происходит для каждой скважины индивидуально в зависимости от геолого-геофизических характеристик разбуриваемого интервала. Недостатком данной методики является необходимость в выборке значений параметров, при которых была зарегистрирована авария.
Список литературы:
- Алали Валид. Еремин Н.А. О созданной в ИПНГ РАН интеллектуальной системе предупреждения осложнений при строительстве скважин на суше и на море // Экспозиция Нефть Газ. 2023. Том 1. С. 27-32.
- Гончаров А.С., Савельев А.О., Саклаков В.М. Алгоритмическое обеспечение интеграции данных о состоянии бурения скважин, получаемых из распределенных гетерогенных источников // XIV Международная научно-практическая конференция, Томск, 28–30 ноября 2018 г. Томск, 2018 С. 22-26.
- Еремин Н.А. Цифровые тренды в нефтегазовой отрасли // Нефть и газ. №12.2017
- Мазмишвили А.И. Теория ошибок и метод наименьших квадратов: учеб. пособие. Недра, Москва: Математические методы, 1978 г. – 200 с.
Оставить комментарий