СПБГУАП группа 4736 Контакты https://new.guap.ru/i03/contacts
1.Математические модели и математический аппарат теории фильтрации
1.1. Уравнение пьезопроводности
Для описания процессов фильтрации жидкости в пласте используют математическую модель, в основе которой лежат несколько предположений: 1) добыча происходит на упругом режиме фильтрации, то есть допустимо считать, что в обычных интервалах изменения пластового давления плотность жидкости и пористость среды линейно зависят от давления [1]; 2)скорость фильтрации прямо пропорциональна градиенту давления,то есть выполняется закон Дарси:
u
k |
p |
|
|
||
|
,
(1.1)
где |
k k x, y,z,t |
- |
проницаемость, |
x, y,z,t |
вязкость, |
x, y,z,t |
- |
пористость, p p x, y, z,t |
- давление в точке пласта x, y, z в момент времени t . |
|
|||||
При вышеизложенных предположениях распределение давления в пласте в любой момент времени описывается с помощью уравнения пьезопроводности, которое в случае
наличия анизотропии по проницаемости имеет вид [1]:
|
|
kx |
2 |
p |
|
k y |
2 |
p |
|
kz |
|
2 |
p |
* p |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
x |
2 |
|
y |
2 |
|
|
z |
2 |
|
t |
,. |
|
(1.2) |
|||
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где - вязкость, |
|
k |
- |
проницаемость в |
|
данном направлении, |
* |
- коэффициент |
|||||||||||
|
|
|
|||||||||||||||||
упругоемкости пласта, |
* |
|
|
|
s , |
f |
и s - |
коэффициенты сжимаемости жидкости и |
|||||||||||
|
f |
||||||||||||||||||
пористой среды.
В случае наличия распределенных стоков и источников в фильтрационном потоке уравнение пьезопроводности примет следующий вид:
kx 2 p |
|
k y 2 p |
|
kz 2 p |
|
* p |
|
q |
, |
|||||
|
|
|
|
|
|
|
|
|
||||||
|
x2 |
y2 |
z 2 |
t |
|
0 |
||||||||
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где q – массовая производительность стоков или источников,
(1.3)
кг .3м с
Рассмотрим скважину в центре однородного изотропного кругового пласта с отсутствием непрерывно распределенных стоков и источников (см. Рисунок 1.1), в случае наличия плоскорадиальной симметрии уравнение (1.3) имеет вид [1]:
1 p r,t |
|
1 p r,t |
|
2 p r,t |
, |
(1.4) |
|||
|
|
|
|
r 2 |
|||||
|
t |
r r |
|||||||
|
|
|
|
||||||
|
|
|
- 3 - |
|
|
|
|||
СПБГУАП группа 4736 Контакты https://new.guap.ru/i03/contacts
где
|
k |
|
|
|
* |
- коэффициент пьезопроводности.
re
Рисунок 1.1. Скважина в центре однородного изотропного кругового пласта
1.2. Начальные и граничные условия
Как уже было сказано, процесс фильтрации описывается уравнением пьезопроводности. Решение задачи фильтрации в однофазном случае в какой-либо области сводится к задаче решения соответствующего уравнения пьезопроводности, то
есть к отысканию необходимой функции |
p r,t . Процесс отыскания решения уравнения в |
частных производных также называется интегрированием этого уравнения. Для того чтобы проинтегрировать уравнение пьезопроводности, необходимо также знать граничные и начальные условия [2].
Начальное условие заключается в задании начального распределения давления в пласте (то есть распределение давления на момент пуска скважины в работу), обычно считают, что пласт вначале невозмущен – когда давление во всех точках пласта одинаково
и равно начальному среднему пластовому давлению:
p r,t p0 , r rw ,re , t 0 . |
(1.5) |
Граничные условия задаются на скважине и на границах пласта. На границах пласта чаще всего рассматриваются так называемые условия 1-го и 2-го родов [2]. Граничное условие 1-го рода – это условие поддержания постоянного давления на данной границе,
реализуются такие условия в пластах с мощным аквифером: |
|
|
p(r,t) |
p0 , t 0 . |
(1.6) |
|
r re |
|
- 4 -
СПБГУАП группа 4736 Контакты https://new.guap.ru/i03/contacts
Граничные условия 2-го рода – это условия поддержания постоянного значения производной давления, то есть поддержания постоянного потока через границу пласта.
Это условие используется, в частности, при описании пластов с отсутствием перетока через границу (значение производной в этом случае равно нулю):
pr
0 , t 0 r re
(1.7)
Также, как будет показано далее, данное условие может быть использовано при описании поведения скважины на псевдоустанвившемся и установившемся режимах в системе разработки.
Существуют еще граничные условия 3-го рода – условия связи давления и производной давления по координате (то есть связь между давлением и потоком через границу пласта)[5]:
pr
r,t p |
r r |
f r,t , t 0 |
|
|
|
r r |
Г |
|
|
|
|
Г |
|
|
(1.8)
В системе заводнения для корректного описания процесса взаимодействия скважин,
как на неустановившемся, так и на псевдоустановившемся и установившемся режимах необходимо использовать, как будет показано далее граничные условия 4-го рода [3],[4] –
условия сопряжения зоны закачки и зоны добычи. Конкретный вид данных граничных условий будет рассмотрен далее.
Также необходимо задать граничное условие на скважине. В зависимости от задания граничного условия на скважине, математические модели работы скважины классифицируются на модели постоянного давления и модели постоянного дебита. В
случае, когда на скважине поддерживается постоянное забойное давление (граничное условие 1-го рода), используется модель постоянного давления, а граничное условие выглядит следующим образом:
p r,t |
r r |
|
|
|
|
|
w |
|
p |
wf |
|
, t
0
,
(1.9)
где pwf - забойное давление.
Модель постоянного дебита используется, соответственно, когда на скважине задано
условие 2-го рода – поддерживается постоянный дебит и представляет собой закон фильтрации Дарси .
p( r,t ) |
|
|
Bo qs |
, t 0 , |
(1.10) |
|
r |
r r |
2 khrw |
||||
|
|
|
w
где qs - дебит в поверхностных условиях, h - мощность пласта, Bo - объемный коэффициент.
- 5 -
СПБГУАП группа 4736 Контакты https://new.guap.ru/i03/contacts
С учетом эффекта послепритока [6]:
p( r,t ) |
|
B q |
|
|
C |
|
p |
wf |
, t 0 |
|
o |
s |
|
s |
|
||||
|
|
|
|
|
|
||||
r |
r r |
2 khr |
|
2 khr |
t |
|
|||
|
|
w |
|
|
w |
|
|
|
|
|
w |
|
|
|
|
|
|
|
|
(1.11)
где Cs - коэффициент послепритока.
Модель постоянного дебита лучше применять, если имеется технические
ограничения по добыче: это может быть предельная пропускная способность поверхностного обустройства, законодательные ограничения на максимальный дебит.
Модель постоянного давления лучше применять в долгосрочном режиме работы скважины, когда дебит меняется сильно, а забойное давление остается приблизительно на одном уровне.
Необходимо упомянуть, что вообще говоря, на скважине меняется как забойное давление так и дебит, поэтому корректнее предположить, что дебит связан с давлением некоторой связью, например – в некоторых случаях с большой точностью можно аппроксимировать эту связь – линейной зависимостью (в первом приближении). Работа скважины с граничным условием 3-го рода, то есть с условием линейной связи забойного давления с дебитом рассматривается в статье [7]. Также есть ряд публикаций,
посвященных работе скважины при условии, что их дебит – есть некоторая функция
времени, например степенная: |
* |
s |
[1]. В [8] решена более общая задача о |
q qs t |
|
||
взаимодействии пласта со скважиной, |
дебит которой меняется следующим образом: |
||
q q0
at |
|
e |
t |
|
|
.
1.3. Способы интегрирования уравнения пьезопроводности
Существуют различные способы нахождения решений уравнения пьезопроводности.
Широкий обзор методов решения уравнения пьезопроводности можно найти в книгах [9], [1], [10], а также [2].
1.3.1. Фундаментальные решения
Рассмотрим одномерные потоки. В «одномерных» потоках все характеристики потока (давление, скорость фильтрации и т.д.) зависят не более чем от одной координаты.
Примерами одномерных потоков могут служить: 1) Прямолинейно-параллельный поток 2)
Плоскорадиальный поток (рисунок) 3) Сферический радиальный.
- 6 -
СПБГУАП группа 4736 Контакты https://new.guap.ru/i03/contacts
Изобары |
а) |
Изобары |
б) |
Линии тока
Лини тока
X
в)
Линии тока
Изобары
Рисунок 1.2. Одномерные потоки в пространстве: а) одного б) двух и в) трех измерений.
В этих простейших случаях интегрирование уравнения пьезопроводности для точечного источника в бесконечном пласте сводится к интегрированию по времени фундаментального решения, которое математически описывается формулой для случаев пространства одного, двух и трех измерений следующими формулами [1], [2], [5]:
|
|
|
|
|
|
|
x |
2 |
||
|
1 |
|
|
|
|
|
|
|||
p |
e |
4 t |
||||||||
|
|
|
||||||||
|
|
|
|
|
|
|||||
|
t |
|
|
|
|
|
|
|
||
|
|
|
|
|
r |
2 |
|
|
||
|
1 |
|
|
|
|
|
|
|
||
p |
e |
4 |
t |
|
|
|||||
|
|
|
|
|||||||
|
|
|
|
|
|
|||||
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
2 |
||
|
1 |
|
|
|
|
|
||||
p |
|
e |
4 t |
|||||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
||||
t |
|
t |
|
|
|
|
|
|
|
|
(1.12)
(1.13)
(1.14)
Физический смысл фундаментальных решений заключается в том, что они представляют собой распределение падения давления в бесконечном пласте в любой момент времени после мгновенного включения в пространстве соответственно бесконечного плоского, бесконечного линейного или точечного источника или стока [1].
Фундаментальные решения уравнения пьезопроводности были впервые получены
- 7 -