Теории упругого режима фильтрации газа

      Формулы (36) пригодны как для бесконечного пласта, так и для конечного открытого и закрытого пласта радиусом R_k. В последнем случае они справедливы только для первой фазы движения, пока воронка депрессии не достигнет границы пласта, т.е. для

      

      Изменение давления во второй фазе зависит от граничных условий пласта. Если пласт  закрыт, то давление будет продолжать снижаться во всем пласте, включая границу. Этот случай будет рассмотрен в последующих параграфах. Если пласт открытый (р=р_к или r=R_k, т.е. режим водонапорный, то во второй фазе установится стационарный режим с постоянной депрессией p_k-p_c где

1.4Метод усреднения.

      Рассмотрим  еще один приближенный метод применительно  к задачам неустановившейся фильтрации газа - метод усреднения временной  производной по пространству.

      В качестве примера рассмотрим прямолинейно-параллельную фильтрацию реального газа. Соответствующее этому случаю точное дифференциальное уравнение имеет вид

      

      Сделаем допущение, что коэффициент сверхсжимаемости z(р) можно заменить на где p_ср некоторое среднее давление в области фильтрации. Введем обозначение р₁=р/z(р).Тогда уравнение примет вид

      

      Пусть имеется первоначально невозмущенный  газонасыщенный пласт шириной В, толщиной h, длиной L. С трех сторон пласт ограничен непроницаемыми поверхностями, а с четвертой стороны (х = 0) вскрыт галереей. В момент t = 0 через галерею начинает отбираться газ с постоянным массовым дебитом, который в соответствии с законом Дарси можно записать в виде:

      

      Требуется определить давление в пласте в любой  момент времени t> 0. Для этого нужно найти решение уравнения (6.43) в области изменения удовлетворяющее начальному и граничным условиям:

p₁=p₁₀ при t=0;

при x=0, где

при x=L.

      Как и в методе последовательной смены стационарных состояний. принимаем, что в каждый момент времени существует конечная возмущенная область l(t), на границе которой выполняются условия

      p₁²=p₁₀², при x=l(t).

      Центральным моментом в рассматриваемом методе усреднения является принятие условия

      

равносильного предположению, что во всей части  пласта, охваченной возмущением, давление изменяется с одинаковой скоростью; тогда уравнение (6.43) принимает вид

      Проинтегрировав это уравнение дважды по х получим:

      

      Использовав граничные условия на галерее (6.45) и на границе возмущенной области (6.47), в результате получим

      

      где

      Подставив выражение l(t) в формулу (6.51), получим явную зависимость давления, от координаты и времени.

      В момент Т, когда возмущенная зона достигнет непроницаемой границы  пласта (l=L) закончится первая фаза.

      Для определения ее продолжительности  положим в уравнении (6.54) l=L и выразим время T:

      

      Можно найти приближенное значение T из формулы (6.55) и убедиться, что погрешность не превышает 3-4%.

      В течение второй фазы давление на границе  x=L падает и выполняется условие (6.46). Соотношения для второй фазы истощения газового пласта строятся аналогичным образом. Проделав аналогичные выкладки, получим закон распределения давления по пласту

      

      и закон изменения давления на галерее 

        

 

      

    2. Расчетная часть

    2.1. Рассчитать депрессию на пласт по точной формуле и по приближенным формулам.

      Исходные данные к расчету:

    -коэффициент  проницаемости k=0.29ּ10^-12 м²

    -толщина  пласта h=6 м

    -радиус  скважины r_c=0.08 м

    -радиус  контура питания R_k=300м

    -коэффициент  вязкости газа η=1,2ּ10^-5Паּс

    -пластовое  давление P_k=13,8МПа

    - коэффициент пористости m₀=0.2

    Также для расчета зададимся временем работы скважины t=1час=3600с и зададимся дебитом Q_ат из условия, что λ=0,004994. Тогда из формулы определим Q_ат

    

    Подставив наши данные, получим:

    

    определим коэффициент 

    

    2.1.1.Точное  решение.

    Определим безразмерную величину ξ для r=r_c

    

    Сравнивая полученное значение ξ со значениями в таблице 1 для λ=0,004994 заключаем, что ξ< ξ* поэтому безразмерное давление F определим по формуле

     .

    Выразив давление P=P_c, получим

    P_c= F P_k = 0.996ּ13,8ּ10⁶ =13,7ּ10⁶ Па

    Депрессия на пласт через 1 час будет равна:

    ∆P=P_k-P_c=(13,8-13,7)ּ10⁶=0,1 МПа.

    2.1.2. Расчет по линеаризованной формуле.

    Подставляя  наши данные в формулу (22), определим  забойное давление P_c через 1 час

    

    Депрессия будет равна:

    ∆ P= (13,8-13,55)ּ10⁶=0,25 МПа.

    2.1.3. Расчет методом  последовательной  смены стационарных  состояний. 

    

    

    Депрессия будет равна:

    ∆ P= (13,8-13,54)ּ10⁶=0,26 МПа.

    2.2 Относительная погрешность  расчетов.

1) Расчет по линеаризованной формуле :

    2) Расчет методом последовательной смены стационарных состояний.

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

    Выводы.

    Линеаризованная формула эффективна только в тех  случаях, когда радиус скважины очень маленький, потому что в этом случае воронка депрессии очень крутая и давление по всему пласту в целом не сильно отличается от начального. Но при больших радиусах скважины эта формула будет давать большую погрешность, т.к. давление по пласту будет сильно отличаться от начального. В отличие от линеаризованной формулы, формула последовательной смены стационарных состояний эффективна для любых радиусов скважин, но только для первой фазы движения, т.е. пока воронка депрессии не достигнет радиуса контура. Как показали расчеты наиболее точной является линеаризованная формула. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

    Литература

1. К.С. Басниев и др. «Подземная гидромеханика». М. Недра. 1993г.