Расчет вязкости при заданном давлении.

Для определения вязкости газа при заданном давлении и температуре необходимо следующее:

1.    Определить вязкость газа при атмосферном давлении и температуре Т - mат.

2.    Найти критические и приведенные параметры газа.

3.    По графику, приведенному в [1, рис. 11.6] , для найденных Рпр и Т пр m*=m / mат.

4.    По известным mат и m* рассчитать искомую вязкость при заданных Р и Т - m=m* mат

 

Решение:

Запишем уравнение притока газа при нелинейном двучленном законе фильтрации газа к скважине:

Р2пл -Р2з=аQ+bQ2

где а и b – коэффициенты фильтрационного сопротивления, зависящие от несовершенства скважины, геометрии зоны дренирования, параметров пласта и свойств газа и определяемые по формулам  приведенным выше.

Для того, чтобы определить характер изменения коэффициентов фильтрационных сопротивлений a и b в процессе разработки месторождения «Ямбург» необходимо оперделить а0 , далее по формуле (1) видно что а  зависит только от h при неизменных остальных параметрах . Найдя коэффициент а из уравнения притока можно определить коэффициент b.

Для нахождения а0 необходимо построить индикаторную кривую в координатах DP2/Q от Q.Принимаем что пластовое давление постоянно и равно 7,97 Мпа. Индикаторная кривая позволит определить коэффициент а0 как отрезок, отсекаемый на оси  DP2/Q .

Полученные данные приведенны в табл.1.

 

Таблица 1.
























































































Pпл, МПа.



Pз, Мпа.



Q,тыс. м/сут



DP2



DP2/Q



7,97



7,83



391



2,212



0,006



7,97



7,76



502



3,3033



0,007



7,97



7,7



524



4,2309



0,008



7,97



7,56



627



6,3673



0,010



7,97



7,48



780



7,5705



0,010



 

Расчет вязкости при заданном давлении.


По данным приведенным в табл. 1 построим индикаторную кривую в координатах DP2 от Q (рис.3) и  DP2/Q от Q (рис.4).



















Расчет вязкости при заданном давлении.







По рис. 4 определим коэффициент а0=0,0015 МПа.сут/тыс.м3

Коэффициент b0 найдем из уравнения :

((DP2/Q)5 –(DP2/Q)1) /Q5-Q1=b0

(0,01-0,006)/(780-391)*103= 1,0283*10-5

МПа.сут/тыс.м3

b0=1,0283*10-5

МПа.сут/тыс.м3

Так как коэффициенты а и b зависят от hгвг без учета изменения Z(P) и m(P) которые изменяются во времени, а а0=0,0015 МПа.сут/тыс.м3 и h0=95.5 м. можно рассчитать а в зависимости от изменения hгвг. Коэффициент b зависит от h2 , b0=1,0283*10-5

МПа.сут/тыс.м3 и h0=95.5 м. Полученные данные приведены в табл.2.

 

Таблица 2. Зависимость а (h) и b(h).
















































































































Годы



hгвг, м.



а

МПа.сут/тыс.м3



b

МПа.сут/тыс.м3



1987



0



0,0015



1,028E-05



1990



2



0,001532



1,07276E-05



1992



5



0,001617



1,14506E-05



1995



13



0,001871



1,3779E-05



1998



20



0,002367



1,64525E-05



2001



34



0,003676



2,47957E-05



2004



46



0,007092



3,82751E-05



2010



63



0,02084



8,87891E-05



 

а = а0 h0 / h0 - hгвг =0,0015*95,5/95,5-2  = 0,001532 МПа.сут/тыс.м3

b = b0 h0 / (h0 – hгвг)2 = 1,0283*10-5

*95,52/(95,5-2)2 =1,07276*10-5 МПа.сут/тыс.м3

 

Коэффициенты фильтрационных сопротивлений а и b зависят также от коэффициента сверхсжимаемости и от вязкости газа которые в свою очередь зависят от давления.

 

ищи здесь, есть все, ну или почти все