Постановка задачи (круговой контур)

 

Задача 2. Внешний и внутренний контуры нефтеносности одно-пластового нефтяного месторождения имеют форму, близкую к окружностям (рис.4). Площадь месторождения можно представить в виде круга радиусом R=2000м. Нефтяная залежь окружена обширной водоносной областью, из которой в нефтеносную часть пласта поступает вода при снижении пластового давление р0=20 М Па. По данным гидродинамических

Рис.4 Контур нефтеносности: 1-внешний, 2-внутренний.

и лабораторных исследований установлено, что средняя проницаемость как нефтеносной, так и водоносной частей пласта одинакова и составляет

0,5·10^{-12 м2}. Толщина пласта в среднем h=10 м,

средняя пористость пласта m=0,3,

начальная нефтенасыщенность S_н=0,45,

насыщенность пласта связанной водой S_св=0,05.

Вязкость нефти и воды в пластовых условиях равны соответственно: .

Коэффициент пьезопроводности .

Добыча жидкость из месторождения изменяется во времени следующим образом: .

где

время ввода месторождения в разработку .

Требуется определить в условиях разработки при упругом режиме в законтурной области пласта изменение пластового давления.

 

Решение задачи 2

 

 

Рис. 5

 

Для расчета изменения во времени давления на контуре нефтяной залежи, используя аппроксимацию соответствующих решений Карслоу и Егера, имеем

 

 

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

Однако, по условиям данной задачи в период разбурения месторождения объем воды, поступающей из законтурной области, и следовательно, отбираемой жидкости из пласта – переменные во времени. Поэтому для расчета давления на контуре нефтяного месторождения P_кон (t) необходимо использовать интеграл Дюамеля, согласно которому

 

 

В условиях задачи q_ж зависит от физического времени t. В интеграл необходимо поставить Поэтому найдем зависимость q_ж = q_ж(τ) или, что то же самое, q_ж = q_ж(λ). Имеем

 

 

Это формула применяется в том случае, если

 

 

Задание для лабораторной работы.

 

№ п/п	Известные параметры	Определить	Число разбивания
	μ, k, x, H, в, tk, α0, μ, k, x, H, в, tk, α0, μ, k, x, H, в, tk, α0, μ, k, x, H, в, tk, α0, μ, k, x, H, в, tk, α0, μ, k, x, H, R, tk, α0, ρ0, μ, k, x, H, R, tk, α0, ρ0, μ, k, x, H, R, tk, α0, ρ0, μ, k, x, H, R, tk, α0, ρ0,   μ, k, x, H, R, tk, α0, ρ0,	ΔP(t)   ΔP(t)     ΔP(t)   ΔP(t)   ΔP(t)   P(t)   P(t)     P(t)   P(t)     P(t)	n=100, 200   n=200, 400     n=300, 600   n=150, 300   n=250, 500   n=300, 600   n=250, 500     n=250, 500   n=200, 400     n=500, 1000

Алгоритм вычисления определенного интеграла

На рис.6 подпрограмма INT(a,b) вычисляет определенный интеграл из формулы (2.1). Интеграл вычисляется одним из методов: метод прямоугольника, трапеций или Симпсона. Параметры а,в указывают на пределы интегрирования. Функция

.

 

 

Структурная схема расчета.

 

 

- - - - - - - - - - ввод начальных данных

t = 0, t_max, Δt ΔP конец

 

Рис.6

Фильтрация жидкости и газа

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

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

 

 

Здесь -вектор скорости фильтрации; - плотность жидкости. Для того, чтобы получить замкнутую систему уравнений, нужно воспользоваться тем, что свойства жидкости (плотность ρ и вязкость μ), так же как и пористость и проницаемость пористой среды, являются функциями давления (мы предположим движение изотермическим). Исходя, из предположения, слабой сжимаемости жидкости и пористой среды получается уравнение [14-15]:

 

где коэффициент .

 

Здесь - коэффициент пьезопроводности. Полученное уравнение обычно называется уравнением упругого режима или уравнением пьезопроводности.

Рассмотрим постановку основных задач теории упругого режима. Определим распределение давления в некоторой замкнутой области пространства D на протяжении промежутка времени 0 ≤ t ≤ T. Из теории уравнения теплопроводности известно, что если задать на границе Г области D линейную комбинацию давления и его производной по нормали к границе области

 

и задать начальное распределение давления в области Д

 

,

 

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

Рассмотрим подробнее физический смысл тех или иных дополнительных условий. Область, в которой ищется распределения давления жидкости, обычно представляет собой пористый пласт, частично имеющий непроницаемые границы, а частично сообщающийся с другими пластами и вскрывающими его скважинами. На непроницаемых границах должно удовлетворяться очевидное условие отсутствие потока – равенство нормального компонента скорости фильтрации нулю:

 

U_n = 0,

 

Откуда используя закон Дарси, получаем.

 

 

На участках границы с областями, в которых перераспределение давления практически не происходит (область питания), давление можно считать постоянным, так что

 

.

 

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

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

Часть границы области фильтрации обычно образована стенками скважины и дренажных галерей. На этой части границы чаще всего задается либо давление жидкости, либо поток ее через стенки скважины. Выбор того или иного условия зависит от режима работы скважины или галерей. Могут быть и более сложные условия, когда задается связь давления с расходом жидкости. Задание потока жидкости согласно закону Дарси эквивалентна заданию нормальной производной от давления. Условия этого типа выполняются на тех участках границы, через которые может происходить обмен жидкости с соседними пластами через сравнительно слабопроницаемые перемычки. Если толщина перемычки Δ мала, а давление за ней можно считать постоянным, то расход вытекающей жидкости через участок перемычки площадью ds составит . Это количество жидкости должно быть равно

 

uде U_n – нормальная проекция скорости фильтрации на рассматриваемом участке границы. Отсюда имеем

 

 

т.е. условия третьего рода.

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

Задача 3. Одинарное прямолинейно-параллельное движение. Пусть скорость фильтрации жидкости параллельны оси х и не зависит от координат y и z. Давление при этом удовлетворяет уравнению Предположим, что в плоскости х = L давление сохраняет постоянное значение, равное контурному.

 

При х = 0 задан расход жидкости q(t). При заданных значениях μ, k, x и L определить давление пласта.

Решение. В работе [5], предполагая, что получено решение поставленной задачи. Она описывается в виде:

 

 

Полученная формула имеет двоякий смысл. С одной стороны, он описывает распределение давления в пласте конечной длины L при малом времени. С другой стороны он дает распределение давления в произвольный момент времени. В пласте «бесконечной» протяженностью L → ∞. Дело в том, что конечное (не бесконечно малое) изменения давления распространяется за заданное время лишь на конечное расстояние и, если рассматриваются малые времена, можно считать пласт бесконечным. Решение задачи для бесконечного пласта автомодельное: независимые переменные х и t входят в решение не порознь, а лишь в комбинации

Если , то [5]

 

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

 

P ₁ (0, t) = P₀, P(L, t) = P(L, 0) = const;

 

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

Оценим это время для систем различных режимов. При этом для χ характерное значение х = 10⁴см²/сек. В результате имеем: при L = 1м (переходный процесс в одном пуске-блике породы) t = 0,1сек; при L = 300м (порядка расстояния между скважинами) t = 10⁴сек»3г. L = 1км (порядка размеров месторождения) t = 10⁷сек £ 100 суткам; L = 100км (порядка размеров крупной водонапорной системы) t = 10⁹сек = 30 лет.

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

Задания для лабораторной работы.

 

№ п/п	Известные параметры	Определить	Число разбивания
	μ, k, χ, q(t) = α0 t, t £ T, P0 (x) μ, k, χ, q(t) = α0tsin , t £ T, P0 (x) μ, k, χ, q(t) = α0ttg ,t £ T, P0 (x) μ, k, χ, q(t) = α0t^2/(t^0 + 2) t £ T, P0 (x) μ, х, q(t) = α0t, t £ T, P (t1), P1 (t1+ Δt) k, х, q(t) = α0tsin P (t1), P (t1+ Δt) k, х, q(t) = α0t, ρ (t1), P (t1+ Δt) μ, х, q(t) = α0tsin , ρP(t1), P (t1+ Δt) μ, k, q(t) = α0t, ρ (t1), P(t1+ Δt) μ, k, q(t) = α0tsin , P(t1), P(t1+ Δt)	P (x, t) P (x, t)   P (x, t)   P (x, t) k, ΔP (x, t)   μ, ΔP (x, t)   μ, ΔP (x, t) k, ΔP (x, t)   x, ΔP (x, t) x, ΔP (x, t)	n=200, 400 n=300, 600   n=400, 800   n=500, 1000 n=250, 500   n=350, 700   n=400, 800 n=500, 1000   n=500, 1000 n=600, 1200

Задача 4. Рассмотрим теперь одномерное осесимметричное (плоско-радиальное) движение при упругом режиме. Распределение давления определяется при этом как решение уравнения теплопроводности в полярных координатах (r, j):

 

 

удовлетворяющие начальному условию

 

P (r, 0) = f(r)

 

а граничным условиям при r = ρ и r = R.

По прежнему основной интерес представляют решения, отвечающие стационарному начальному распределению давления

 

f(r)=C₁lnr+C₂.

 

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

 

 

 

где – постоянная размерность давления.

 

 

Рис. 7

 

Предполагая, что в работе [5] получена решение поставленной задачи в виде

(28)

 

Отметим важное обстоятельство: в соотношения для P (r, t) не входит радиус скважины . Это означает, что в области применяемости условия распределения не зависит от радиуса скважины.

(28) – является несобственным интегралом, с помощью замены переменной приведем его собственному интегралу. Для этого сделаем замену:

 

 

Первоначальное дифференциальное уравнение является линейным, поэтому здесь применим принцип суперпозиции. Если имеется две скважины расположенные друг от друга на расстоянии l, то давление пласта в точке А удаленной от первой скважины на расстоянии будет равен (рис. 7)

 

 

Задания для лабораторной работы

№ п/п	Известные параметры	Определить	Число разбивания
	k, χ, μн, q1, q2 , R1, R2, t £ T k, χ, μн, q1, q2 , 0 £ х^2+ у^2 £ R^2, t £ T k, χ, μн, q1, q2 , q3, 0 £ х^2+ у^2 £ R^2, t £ T k, χ, μн, q1, q2 , q3, q4, 0 £ х^2+ у^2 £ R^2, t £ T k, χ, μн, q1=q2, 0 £ х^2+ у^2 £ R^2, t £ T k, χ, μн, q1= q2 = q3, 0 £ х^2+ у^2 £ R^2, t £ T	ΔP (t) ΔP(x, у, t) ΔP(x, у, t) ΔP(x, у, t) ΔP(x, у, t) ΔP(x, у, t)	n=500 n=300 n=200 n=300 n=250 n=400