4. Розрахунок показників розробки газового родовища при режимі

Р_роб =соnst

Проводимо розрахунок при режимі заданого тиску на усті свердловини (Р_роб =соnst_.=21,3 МПа).

„Контрольна точка” – станом на 01.01.2013 року.

Визначаємо приведені параметри пластового газу як відношення тиску і температури до їх критичних значень:

Р_пр= = = 7,43 МПа

Т_пр= = = 1,87 К

	1,8	1,9	1,87
7,4	1,005	1,021	1,016
7,45	1,007	1,023	1,018
7,43			1,017

За таблицями Катца знаходимо Z_сер.:

Z_сер.=1,017

Знаходимо приведений газонасичений об’єм за формулою:

Ω^* = = =6,16 млн.м³ /МПа

Визначаємо приведені параметри пластового газу як відношення тиску і температури до їх критичних значень:

Т_ср. = = (371-293)/ =330,5 К

Р_пр4.1= = 4,59 МПа

Т_пр= = = 1,67 К

	1,6	1,7	1,67
4,55	0,830	0,865	0,8545
4,6	0,832	0,867	0,8565
4,59			0,8561

За таблицями Катца знаходимо Z_сер.:

Z_сер.4.1= 0,8561

Уточнюємо значення θ:

= 0,0301

Уточнюємо значення S:

S_4.1= = 0,212

За уточненими значеннями S,θ,Z_сер визначають дебіт q(t):

q(t)_4.1= = 115,71тис.м³/добу

За знайденим значенням q(t) уточнюємо Р_виб.(t):

Р_виб.4.1(t)=√P_роб.²∙ e^2S+ θ ∙q²(t)= 33,11 МПа

Уточнюємо значення Z_сер. для цього знайдемо Р_ср.:

Р_ср= Р_виб.+ )= 27,63 МПа

Р_{пр 4.2} = = 5,95 МПа

Т_{пр 2} = = 1,67 К

За таблицями Катца знаходимо Z_{сер 2}:

	1,6	1,7	1,67
5,95	0,889	0,915	0,9072

Z_сер.4.2= 0,9072

Уточнюємо значення θ _4.2 :

θ _4.2= =0,0338

Уточнюємо значення S:

S_4.2= = 0,2

За уточненими значеннями S, θ,Z_сер визначають дебіт q(t):

q(t)_4.2 = = 111,94 тис.м³/добу

За знайденим значенням q(t) уточнюємо Р_виб.(t):

Р_виб4..2(t) =√P_роб.²∙ e^2S+Ɵ∙q²(t)= 33,17 МПа

∆=(( Р_{виб4. 1}- Р_{виб 4.2})/ Р_{виб 4.1})*100% =((33,17-33,11)/ 33,17)*100%=0,18%

Р_виб1=33,17МПа

Р_пл=34,49 МПа

∆Р1=34,49-33,17=1,32 МПа

Р_роб=21МПа

q(t)1= 115,71 тис.м³/добу

4.1. Розрахунок "контрольних точок"

"Kонтрольнa точкa"- на 31.12.2013 року

Виконуємо перше наближення.

1.Знайдемо сумарний видобуток газу на 31.12.2013 року:

Q_вид.=

де: ∆t=t-t_n-1

n-кількість свердловин на початку розрахунку;

Q_вид.(t_n-1)- сумарний видобуток газу на попередній момент часу t_n-1;

Q_вид=0 млн.м³;

q(t_n-1) – дебіт середньої свердловини на попередній момент часу t_n-1;

- коефіцієнт експлуатації свердловини (відношення числа днів роботи свердловини у році до календарного числа днів):

=330/365=0,9

Підставляючи дані в формулу та приймаючи в першому наближенні:

Q_вид.4.1=0+ ∙0.9∙1∙0,365=38,011 (млн.м³)

Знаходимо приведений газонасичений об’єм за формулою:

Ω^* = = = 6,16 млн.м³ /МПа

Визначаємо середній пластовий тиск методом ітерації за формулою:

P_пл.(t)=(P_n/Z_n- Q_вид(t)/Ω^*)∙Z(P_пл.)

Приймаючи Z(P_пл.) рівним його попередньому значенню:

P_пл4.1.1.(t) = (34,49/1,017-38,011/6,16)∙1,017= 28,21 МПа

За знайденим значенням пластового тиску уточнюємо Z(P_пл.) за таблицями Катца:

Р_пр4.1.1= = 6,08 МПа

Т_пр.= =371/197,9=1,87 К

	1,8	1,9	1,87
6,05	0,948	0,968	0,962
6,1	0,950	0,970	0,964
6,08			0,9632

Z_сер4.1.1. = 0,9632

P_пл.4.1.2(t) = (34,49/1,017-38,011/6,16)∙0,9632= 26,72МПа

∆=(( Р_пл.4.1.1- Р_пл4.1.2)/ Р_{виб 4.1.1})*100% =((28,21-26,72)/ 28,21)*100%=5,28%

За знайденим значенням пластового тиску уточнюємо Z(P_пл.) за таблицями Катца:

Р_пр4.12= = 5,76 МПа

Т_пр= =371/197,9=1,87 К

	1,8	1,9	1,87
5,75	0,937	0,957	0,951
5,8	0,939	0,959	0,953
5,76			0,9514

Z(P_пл.)_4.1.2= 0,9514

Приймаючи Z(P_пл.) рівним його попередньому значенню:

P_пл.4.1.3(t)= (34,49/1,017-38,011/6,16)∙0,9514= 26,39 МПа

∆=(( Р_пл4.1.2- Р_пл4.1.3)/ Р_пл4.1.3)*100% =((26,75-26,39)/ 26,75)*100%=1,35%

Знайдемо дебіт свердловини:

q(t)_4.1.1= =17,08 тис.м³/добу

Приймаючи в першому наближенні значення S,θ,Z_сер рівними їх значенням на попередній момент часу(початок розрахунків),визначаємо дебіт «середньої» свердловини.

Обчислюємо поточний вибійний тиск з формули:

Р_виб4.1.1.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 26,2 МПа

Уточнюємо значення Z_сер для цього знайдемо Р_ср для P_роб. та знайденого Р_виб.(t):

Р_ср4.1.1= Р_виб.+ )= 23,83 МПа

Р_пр2.1.1= = 5,14 МПа;

Т_пр2= =1,67 К;

	1,6	1,7	1,67
5,10	0,850	0,882	0,8724
5,15	0,851	0,884	0,8741
5,14			0,8738

За таблицями Катца знаходимо Z_сер.:

Z_сер4.1.1= 0,8738

Уточнюємо значення S:

S_4.1.1= 0,207

Уточнюємо значення Ѳ:

Θ_4.1.1= 0,0313

За уточненими значеннями S,та θ, обчислюємо q(t):

q(t)_4.1.2= = 11,13 тис.м³/добу

За знайденими значеннями q(t) уточнюємо Р_виб.(t):

Р_виб4.1.2.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 26,27 МПа

Виконуємо друге наближення:

Уточнюємо річний видобуток газу на кінець 2013 року:

Q_вид.4.2=0+ ∙0.9∙1∙0,365= 20,83 (млн.м³)

Визначаємо середній пластовий тиск методом ітерації за формулою:

P_пл.(t)=(P_n/Z_n- Q_вид(t)/Ω^*)∙Z(P_пл.)

Приймаючи Z(P_пл.) рівним його попередньому значенню

P_пл.4.2.1(t) = (34,49/1,017-20,83/6,16)∙0,9445= 29,05 МПа

Р_пр4.2.1.= = 6,26 МПа

Т_пр.= =371/197,9=1,87 К

	1,8	1,9	1,87
6,25	0,9560	0,975	0,9693
6,3	0,958	0,977	0,9713
6,26			0,9697

Z_сер. = 0,9697

P_пл.4.2.2(t) = (34,49/1,017-20,83 /6,16)∙0,9697= 29,61 МПа

∆=(( Р_пл4.2.2- Р_пл4.2.1)/ Р_пл4.2.2)*100% =((29,61-29,05)/ 29,61)*100%=1,89%

За знайденим значенням пластового тиску уточнюємо Z(P_пл.) за таблицями Катца:

Р_пр4..2.2= = 6,38 МПа

Т_пр= =371/197,9=1,87 К

	1,8	1,9	1,87
6,35	0,960	0,979	0,9733
6,4	0,962	0,980	0,9746
6,38			0,9741

Z(P_пл.)_4.2.2= 0,9741

Приймаючи Z(P_пл.) рівним його попередньому значенню

P_пл.4.2.3(t)= (34,49/1,017-20,83 /6,16)∙0,9741= 29,74 МПа

∆=(( Р_пл4.2.3- Р_пл4.2.2)/ Р_пл4.2.3)*100% =((29,74-29,61)/ 29,74)*100%=0,44%

Знайдемо дебіт свердловини з формули:

q(t)_4.2.1= = 69,22 тис.м³/добу

Приймаючи в першому наближенні значення S,θ,Z_сер рівними їх значенням на попередній момент часу(початок розрахунків),визначаємо дебіт «середньої» свердловини.

Обчислюємо поточний вибійний тиск з формули:

Р_виб4.2.1(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 28,92 МПа

Уточнюємо значення Z_сер для цього знайдемо Р_ср для P_роб. та знайденого Р_виб.(t)

Р_ср4.2.1= Р_виб.+ )= 25,3 МПа

Р_пр2= = = 5,45 МПа;

Т_пр2= =1,67 К;

	1,6	1,7	1,67
5,45	0,865	0,895	0,886

За таблицями Катца знаходимо Z_сер.

Z_сер4,2.2.= 0,886

Уточнюємо значення S

S_4.2.2= 0,2045

Уточнюємо значення θ

θ_4.2.2= 0,0308

За уточненими значеннями S,та θ, обчислюємо q(t):

q(t)_4.2.2= = 70,34 тис.м³/добу

За знайденими значеннями q(t) уточнюємо Р_виб.(t)

Р_виб2.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 28,9 МПа

∆=(( Р_{виб 1}- Р_{виб 2})/ Р_{виб 2})*100% =((28,9-26,27)/ 28,9)*100%=9,1%

Виконуємо третє наближення:

Уточнюємо річний видобуток газу на кінець 2013 року:

Q_вид.4.3=0+ ∙0.9∙1∙0,365= 30,559 (млн.м³)

Визначаємо середній пластовий тиск методом ітерації за формулою:

P_пл4.3.1(t)=(P_n/Z_n- Q_вид(t)/Ω^*)∙Z(P_пл.)

Приймаючи Z(P_пл.) рівним його попередньому значенню

P_пл.4.3.1(t) = (34,49/1,017-30,559/6,16)∙0,9741= 28,2МПа

Р_пр4.3.1= = 6,08 МПа

Т_пр.= =371/197,9=1,87 К

	1,8	1,9	1,87
6,05	0,948	0,968	0,962
6,1	0,950	0,970	0,964
6,08			0,9632

Z_4.3.1 = 0,9632

P_пл.4.3.2(t) = (34,49/1,017-30,559/6,16)∙0,9632= 27,89Мпа

∆=(( Р_пл4.3.1- Р_пл4.3.2)/ Р_пл4.3.1)*100% =((28,2-27,89)/ 28,2)*100%=1,1%

Знайдемо дебіт свердловини з формули:

q(t)_4.3.1= = 46,08 тис.м³/добу

Приймаючи в першому наближенні значення S, θ,Z_сер рівними їх значенням на попередній момент часу(початок розрахунків),визначаємо дебіт «середньої» свердловини.

Обчислюємо поточний вибійний тиск з формули:

Р_виб.4.3.1(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 27,36 МПа

Уточнюємо значення Z_сер для цього знайдемо Р_ср для P_роб. та знайденого Р_виб.(t)

Р_ср4.3.1= Р_виб.+ )= 24,46 МПа

Р_пр4.3.1= = 5,27 МПа;

Т_пр2= =1,67 К;

	1,6	1,7	1,67
5,25	0,855	0,887	0,8774
5,3	0,858	0,889	0,8797
5,27			0,8783

За таблицями Катца знаходимо Z_сер.

Z_сер4.3.1= 0,8783

Уточнюємо значення S

S_4.3= 0,206

Уточнюємо значення θ

θ₂= 0,0306

За уточненими значеннями S,та θ, обчислюємо q(t):

q(t)_4.3.2= = 45,62 тис.м³/добу

За знайденими значеннями q(t) уточнюємо Р_виб.(t)

Р_виб4.3.2(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 27,36 МПа

∆=(( Р_{виб 4.3.2}- Р_{виб 4.2.2})/ Р_{виб 4.3.2})*100% =((28,9-27,36)/ 28,9)*100%=5,32%

Виконуємо четверте наближення:

Уточнюємо річний видобуток газу на кінець 2013 року:

Q_вид4.4=0+ ∙0.9∙1∙0,365= 26,498 (млн.м³)

Визначаємо середній пластовий тиск методом ітерації за формулою:

P_пл.(t)=(P_n/Z_n- Q_вид(t)/Ω^*)∙Z(P_пл.)

Приймаючи Z(P_пл.) рівним його попередньому значенню

P_пл.4.4.1(t) = (34,49/1,017-26,498/6,16)∙0,9632= 28,52 МПа

Р_пр4.4.1.= = 6,15 МПа

Т_пр.2= =371/197,9=1,87 К

	1,8	1,9	1,87
6,15	0,952	0,971	0,9653

Z_сер. = 0,9653

P_пл.4.4.2 (t) = (34,49/1,017-26,498/6,16)∙0,9653= 28,58 МПа

∆=(( Р_пл4.4.2- Р_пл4.4.1)/ Р_пл4.4.2)*100% =((28,58-28,52)/ 28,58)*100%=0,2%

Знайдемо дебіт свердловини з формули:

q(t)_4.4.1= = 55,68 тис.м³/добу

Обчислюємо поточний вибійний тиск з формули:

Р_виб4.4.1(t)=√P_роб.²∙ e^2S+θ ∙q²(t)= 27,93 МПа

Уточнюємо значення Z_сер для цього знайдемо Р_ср для P_роб. та знайденого Р_виб.(t)

Р_ср4.4.1= Р_виб.+ )= 24,76 МПа

Р_пр4.4.1= = 5,35 МПа;

Т_пр2= =1,67 К;

За таблицями Катца знаходимо Z_сер.

	1,6	1,7	1,67
5,35	0,860	0,891	0,8817

Z_сер2.= 0,8817

Уточнюємо значення S_4.4

S_4.4= 0,2055

Уточнюємо значення θ

Θ_4.4= 0,0307

За уточненими значеннями S,та θ, обчислюємо q(t):

q(t)_4.4.2= = 55,77 тис.м³/добу

За знайденими значеннями q(t) уточнюємо Р_виб.(t)

Р_виб4.4.2.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 27,92 МПа

∆=(( Р_{виб 4.4.2}- Р_{виб 4.3.2})/ Р_{виб 4.4.2})*100% =((27,92-27,36)/ 27,92)*100%=1,21%

Виконуємо п’яте наближення:

Уточнюємо річний видобуток газу на кінець 2013 року:

Q_вид4.5=0+ ∙0.9∙1∙0,365= 28,166 (млн.м³)

Визначаємо середній пластовий тиск методом ітерації за формулою:

P_пл.4.5.1(t)=(P_n/Z_n- Q_вид(t)/Ω^*)∙Z(P_пл.)

Приймаючи Z(P_пл.) рівним його попередньому значенню

P_пл4.5.1.(t) = (34,49/1,017-28,166/6,16)∙0,9653= 28,32 МПа

Р_пр4.5.1.= = 6,1 МПа

Т_пр.4.5.1= =371/197,9=1,87 К

	1,8	1,9	1,87
6,1	0,950	0,970	0,964

Z_4.5.1. = 0,964

P_пл4.5.2(t) = (34,49/1,017-28,166/6,16)∙0,964= 28,28 МПа

∆=(( Р_пл4.5.2- Р_пл4.5.1)/ Р_пл4.5.2)*100% =((28,32-28,28)/ 28,32)*100%=0,14%

Знайдемо дебіт свердловини з формули:

q(t)_4.5.1= = 52,16 тис.м³/добу

Приймаючи в першому наближенні значення S,θ,Z_сер рівними їх значенням на попередній момент часу(початок розрахунків),визначаємо дебіт «середньої» свердловини.

Обчислюємо поточний вибійний тиск з формули:

Р_виб4.5.1.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 27,71 МПа

Уточнюємо значення Z_сер для цього знайдемо Р_ср для P_роб. та знайденого Р_виб.(t)

Р_ср4.5.1= Р_виб.+ )= 24,64 МПа

Р_пр4.5.1= = 5,31 МПа;

Т_пр= =1,67 К;

	1,6	1,7	1,67
5,3	0,858	0,889	0,8797
5,35	0,860	0,891	0,8817
5,31			0,8801

За таблицями Катца знаходимо Z_сер.

Z_сер.= 0,8801

Уточнюємо значення S

S_4.5= 0,2059

Уточнюємо значення θ

Θ_4.5= 0,03066

За уточненими значеннями S,та θ, обчислюємо q(t):

q(t)_4.5.2= = 51,48 тис.м³/добу

За знайденими значеннями q(t) уточнюємо Р_виб.(t)

Р_виб4.5.2.(t)=√P_роб.²∙ e^2S+Ɵ∙q²(t)= 27,68 МПа

∆=(( Р_{виб 4.5.2}- Р_{виб 4.5.1})/ Р_{виб 4.5.2})*100% =((27,71-27,68)/ 27,71)*100%=0,11%

Технологічний режим на кінець 2013 року:

Р_виб=27,68 МПа

Р_пл=28,28 МПа

∆Р=28,28-27,68=0,6 МПа

Р_роб=21,3МПа

q(t)= 51,48 тис.м³/добу

Перевірку розрахунків виконала програма Gas Drive. Результати перевірки зазначені в таблиці «Розрахунок показників розробки».