МИНОБРНАУКИ РОССИИРГУ нефти и газа (НИУ) имени И.М. ГубкинаФакультет разработки нефтяных и газовых месторожденийКафедра управление разработкой нефтяных месторождений
Курс: «Модели оптимальной разработки и обустройства месторождений нефти и газа»Домашнее задание №1 «Оптимизация режимов эксплуатации скважин»Вариант №1Выполнил:студент группы РНМ-22-04Баринов Н.В.Проверил:профессор, д.т.н.Ермолаев А. И.Москва 2023
ЗАДАНИЕРассматривается эксплуатация 2-х нефтяных добывающих скважин при естественном режиме работы пласта, i=1,2.
Исходные параметры (Вариант 1):с
1=0,15 – доля нефти в продукции 1-й скважины; с
2=0,93 – доля нефти в продукции 2-й скважины; ΔP
max1=1,5 МПа – максимально допустимый перепад между давлением на контуре питания залежи и забойным давлением 1-й скважины;ΔP
max2=1,5 МПа – максимально допустимый перепад между давлением на контуре питания залежи и забойным давлением 2-й скважины; R
к=4,8 км – радиус контура питания залежи; R
12=R
21=100 м - расстояние между скважинами; R
11=R
22=R
c=0,1 м – радиус скважин; h =15 м – толщина пласта; k=0,2 Д – проницаемость пласта; μ=0,8 сПз – вязкость пластового флюида. Свободные переменные: q
2, q
3.
Требуется найти: такие дебиты скважин по жидкости – q
1 и q
2 и забойные давления скважин – Р
1 и Р
2, которые обеспечат максимальный суммарный дебит по нефти при выполнении ограничения на величину перепадов между давлением на контуре питания залежи и забойными давлениями скважин.
Решение:Процесс эксплуатации группы взаимовлияющих скважин описывается системой линейных уравнений относительно q
i и P
i – параметров, характеризующих режимы работы скважин:
(1)
Математическая постановка задачи имеет вид модели линейного программирования:
(2)
(3)
(4)
Если найдены дебиты скважин, то при известном давлении на контуре питания по формулам (1) не трудно определить их забойные давления.
С учетом формул (2) - (4) математическая постановка задачи имеет вид:
(5)
(6)
(7)
(8)
Определим числовые значения параметров, входящих в задачу (5) - (8):
;
;
Теперь задача (5) - (8) принимает вид:
(9)
(10)
(11)
(12)
Дополняя задачу (9) - (12) искусственными переменными q3 и q4, которые
также должны подчиняться условию (12), ограничения-неравенства (10) и (11)
можно заменить ограничениями-равенствами. После этого симплекс-методом
подлежит решению задача:
(13)
(14)
(15)
(16)
Итерация № 1.
Шаг 1. Пусть свободными переменными являются q2 и q3, а базисными, соответственно, q1 и q4.
Шаг 2. Выразим базисные переменные через свободные, исходя из уравнений (14), (15):
-
из уравнения (14):
10,8q
1 = 2048 – 3,9q
2 – q
3 q
1 = 189,63 – 0,36q
2 – 0,09q
3; (17)
-
из уравнений (15) и (17):
q
4 = 2048 – 3,9q
1 – 10,8q
2 = 2048 – 3,9(189,63 – 0,36q
2 – 0,09q
3) – 10,8q
2 = 2048 – 739,56 + 1,4q
2 + 0,35q
3 – 10,8q
2;q
4 = 1308,44 – 9,4q
2 + 0,35q
3 (18)Проверим чему окажутся равны
q1 и
q4, если положить
q2=
q3=0:
q1 =189,63≥0,
q4=1308,44>0. Следовательно, можно переходить к шагу 3.Шаг 3. Выразим функцию цели (13) только через свободные: (19)Анализируя формулу (19), можно сделать вывод, что условие оптимальности выполнено: в функции цели (19) все коэффициенты при свободных переменных не имеют строго положительных значений. Поэтому оптимальное решение задачи (13)-(16) имеет вид: q
3=0; q
4=1308,44; q
1=189,63; q
2=0. При этом суммарный дебит по нефти составит 199,31 м
3/сут. Анализируя формулу (19), можно сделать вывод, что условие оптимальности выполнено: в функции цели (19) все коэффициенты при свободных переменных не имеют строго положительных значений. Поэтому оптимальное решение задачи (13)-(16) имеет вид: q
3=0; q
4=1308,44; q
1=189,63; q
2=0. При этом суммарный дебит по нефти составит 28,39 м
3/сут. Из полученных результатов и уравнений (1) следует, что ΔP
i* – разность между давлением на контуре питания залежи и оптимальным забойным давлением i-й скважины (i=1,2): (20) (21)Таким образом, согласно (20)-(21) разность между давлением на контуре питания залежи и оптимальным забойным давлением i-й скважины (i=1,2) равняется максимально допустимому перепаду: так как скважина 2 не работает, то ΔР
2* будет меньше, чем максимально допустимый перепад, ΔР
1*= =1 Мпа.