В упругой пористой среде

Неустановившиеся процессы возникают при пуске и остановке скважин, либо при изменении темпов отбора флюидов. Такие процессы характеризуются перераспределением давления, а также изменением скоростей фильтрационных потоков и дебитов во времени и зависят от упругих свойств пластов и насыщающих их жидкостей. Т.е. основной формой пластовой энергии, обеспечивающей приток жидкости к скважине, – энергия упругой деформации жидкости и материала пласта. При снижении пластового давления объём сжатой жидкости увеличивается, а объём порового пространства сокращается за счет расширения материала пласта, что определяет вытеснение жидкости из пласта в скважину. Хотя коэффициенты сжимаемости жидкости и пласта малы ( , ,

), но зато велики объемы пластов и за счет этого при упругом режиме, могут быть значительные притоки жидкости.

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

Подсчет упругого запаса жидкости в пласте.

Выделим V₀ – элемент объема пласта. Тогда, V_0ж – объем жидкости насыщающей этот элемент пласта при начальном давлении р₀ равен:

V_0ж = m V₀ (7.1)

В соответствии с законами Гука, изменение упругого запаса жидкости DV_з в объеме V₀ при изменении давления на Dр определяется как:

(7.2)

где - коэффициент упругоемкости пласта, численно равный изменению упругого запаса жидкости в единице объема при изменении пластового давления на единицу.

Продифференцировав (7.2) по времени и учитывая, что , получим:

(7.3)

Если формулы (7.1) - (7.3) относить к разрабатываемому в условиях замкнуто-упругого режима нефтяному месторождению, то под V₀ следует понимать объем пласта, в котором к данному моменту времени произошло изменение давления на величину Dр, при этом, , где - начальное пластовое давление; - средневзвешенное по объему возмущенной части пласта V₀ давление. Вычислить средневзвешенное пластовое давление можно, если известна геометрия возмущенной части пласта и конкретное распределение давления в ней, по формуле:

. (7.4)

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

1) уравнение неразрывности,

(7.5)

2) закон Дарси,

; (7.6)

3) уравнение состояния сжимаемой жидкости:

; (7.7)

4) зависимость пористости от давления:

. (7.8)

Подставив (7.6 – 7.8) в (7.5) и пренебрегая членами второго порядка малости, получим уравнение пьезопроводности:

- в дек. системе координат (7.9)

или

- в случае осевой симметрии (7.10)

где: - коэффициент пьезопроводности пласта, характеризующий темп перераспределения пластового давления в условиях упругого режима.

Некоторые точные решения уравнения пьезопроводности.

1. Плоско-параллельный случай, приток упругой жидкости в полубесконечном пласте к прямолинейной галерее скважин.

Для рассматриваемого одномерного движения жидкости уравнение пьезопроводности запишется в виде:

(7.11)

Уравнение (7.11) решается при следующих начальных и граничных условиях:

(7.12)

Условия (7.12) можно интерпретировать таким образом: в начальный момент времени t = 0 пластовое давление было всюду в пласте одинаковым и равным р_к. При пуске скважины в момент времени t > 0 на галерее при х = 0 давление мгновенно упало до величины р_г, при этом на бесконечности x = ¥ давление остается постоянным и равным начальному пластовому р_к.

Решение задачи (7.11) - (7.12) получено методом автомодельной переменной и имеет следующий вид, рис.24:

(7.13)

где: - автомодельная переменная, а

- интеграл вероятности или интеграл Гаусса, который табулирован и имеется в математических справочниках.

Дебит галереи Q при x = 0 выражается в виде:

, (7.14)

т.е. с течением времени Q убывает ~

2. Плоско-радиальный случай реализуется в задаче о притоке упругой жидкости к скважине (точечному стоку или источнику) на плоскости в неограниченном пласте с постоянной мощностью и абсолютной проницаемостью.

В этом случае уравнение пьезопроводности имеет вид:

, (7.15)

и решается при следующих граничных и начальных условиях:

(7.16)

Условия (7.16) интерпретируются иначе чем (7.12): в начальный момент времени t = 0 пластовое давление было всюду в пласте одинаковым и равным р_к. В момент времени t > 0 в точке r = 0 начинает работать добывающая скважина с постоянным объемным дебитом Q₀, на бесконечности r = ¥ давление остается неизменным и равным р_к.

Уравнение (7.15) при граничных условиях (7.16) также решается методом автомодельной переменной:

. (7.17)

Перераспределение давления в пласте выражается, рис.25:

, (7.18)

где: - интегрально показательная функция, которая табулирована и имеется в математических справочниках.

(7.19)

Перераспределение дебита Q(r,t) в пласте:

(7.20)

Определение коллекторских свойств пласта по данным исследования скважин при упругом режиме. Кривые восстановления (падения) давления.

1. Метод касательной. После пуска или остановки скважины при помощи скважинных манометров снимают зависимость забойного давления от времени и строят график такой зависимости: при остановке скважины - кривую восстановления давления, при пуске скважины кривую падения давления. Обработка таких кривых основана на полученном решении (7.18), которое преобразуют, разложив интегрально-показательную функцию в ряд и ограничиваясь первыми двумя членами разложения:

(7.21)

При условии, что значение давления замеряют на забое и r = r_c, формула (7.18) записывается в виде:

(7.22)

где: , (7.23)

(7.24)

Обычно кривую восстановле-ния (падения) давления строят в координатах Dр_с от lg(t). По прямому участку полученной кривой находится отрезок А, отсекаемый его продолжением на оси Dр_с и тангенс угла наклона j, равный i, рис.26. Затем по (7.23) и (7.24) определяют коэффициент гидропроводности пласта и коэффициент пьезопроводности пласта c, по которым в свою очередь можно оценить фильтра-ционноемкостные характеристики горных пород.

Если скважина была пущена в эксплуатацию с постоянным дебитом Q и через промежуток времени T остановлена, но продолжает работать с тем же дебитом Q, то, используя метод суперпозиции, можно получить обобщенную формулу для обработки КВД методом касательной:

(7.25)

В таком случае, удобнее кривую восстановления (падения) давления строить и анализировать в координатах Dр_с от ln(t×T/(t+T)).

2. Метод последовательной смены стационарных состояний (ПССС). Метод основан на предположении, что давление в пласте меняется во времени значительно медленнее, чем по координатам. Поэтому производную по времени можно в первом приближении отбросить, в результате чего для давления получается уравнение Лапласа, описывающее стационарный процесс.

В каждый момент времени вся область движения жидкости, в действительности охватывающая весь пласт, условно разделяется на две области - возмущенную и невозмущенную. При этом предполагается, что в возмущенной области, начинающейся от стенки скважины, давление распределяется так, как будто бы движениё жидкости в ней установившееся: внешняя граница этой области служит в данный момент контуром питания. В невозмущенной области пласта давление всюду постоянно и равно начальному.

В соответствии с методом ПССС принимаем, что через время t после пуска скважины вокруг нее образуется возмущенная область радиусом R(t), где давление будет распределяться по стационарному закону:

(7.26)

Закон движения подвижной границы раздела возмущенной и невозмущенной областей определяется при помощи уравнения .материального баланса (7.3) и граничных условий:

(7.27)

Формула для анализа КВД принимает вид:

(7.28)

А кривую восстановления (падения) давления удобней строить в координатах Dр_с от ln(t).

3. Метод Хорнера. В методе Хорнера реализован метод суперпозиции (наложения фильтрационных потоков) в теории упругого режима, что позволяет решать задачи, связанные с пуском, остановкой или с изменением темпа добычи скважины. Например, если скважина была пущена в эксплуатацию с постоянным дебитом Q и через промежуток времени T остановлена, но продолжает работать с тем же дебитом Q, то для обработки кривых восстановления давления скважин можно получить следующую формулу:

, (7.29)

где T – время работы скважины с расходом Q.

Тогда кривую восстановления давления строят в координатах Dр_с от , рис. 27. При этом, по углу наклона j определяется гидропроводность скважины:

, (7.30)

а, экстраполируя прямолинейный участок кривой восстановления давления до пересечения с прямой, параллельной оси ординат и проведенной из точки , находят начальное пластовое давление .

Если дебит исследуемой скважины с момента пуска ее в эксплуатацию до момента остановки не оставался постоянным, то промежуток времени T можно подсчитать по приближённой формуле:

, (7.31)

где V – накопленная добыча жидкости из скважины с момента пуска до момента остановки, Q – текущий установившийся дебит скважины непосредственно перед ее остановкой.