CC BY
35
УДК 532.546:519.63
Х.М. ГАМЗАЕВ, канд. техн. наук, доц., Азербайджанская
государственная нефтяная академия, Баку, Азербайджан
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ
УПРУГОВОДОНАПОРНОГО РЕЖИМА ПЛАСТА
Для описания упруговодонапорного режима разработки пласта предлагается модель типа Стефана, которая характеризуется наличием неизвестной границы раздела жидкостей. Условие на границе раздела считается неизвестным и взамен этого задается дополнительное условие на фиксированной границе пласта. Для численного решения поставленной задачи предложен вычислительный алгоритм на основе методов выпрямления фронтов и разностной аппроксимации. Табл.: 1. Библиогр.: 8 назв.
Ключевые слова: упруговодонапорный режим, пласт, метод выпрямления фронтов, разностная аппроксимация.
Постанова задачи. Известно, что большинство нефтяных месторождений приурочено к водоносным пластам и разрабатывается в условиях упруговодонапорного режима. При проявлении этого режима движение нефти к скважинам происходит не только за счет потенциальной энергии упругой деформации пласта и нефти, а также в силу давления краевой воды [1]. В большинстве случаев вытеснение нефти краевой водой происходит полностью и в пласте образуется четкая граница раздела двух жидкостей, движущаяся по неизвестному заранее закону.
Предположим, что рассматривается недеформируемый, однородный горизонтально расположенный нефтеносный пласт протяженностью Ь , постоянной толщины и ширины. В сечении пласта х = 0 расположена галерея эксплуатационных скважин, а внешняя граница пласта окружена краевой водой, находящейся под давлением р!1 (?). Пусть в момент времени / = 0 впускается в эксплуатацию галерея эксплуатационных скважин и в пласте возникает упруговодонапорный режим. За счет потенциальной энергии упругой деформации нефти и за счет давления краевой воды происходит прямолинейно-параллельное течение нефти к скважинам. По мере отбора нефти через галерею краевая вода поступает в пласт, полностью замещая поры занятые нефтью, и образуется четкая граница раздела вода-нефть. Предполагается, что нефть является слабосжимаемой жидкостью и ее движения в пласте подчиняются закону Дарси. Тогда уравнение, описывающее нестационарное прямолинейнопараллельное течение нефти в пласте с учетом активного продвижения краевой воды можно представить в следующем виде
дР д2 Р I л
— = %—Т , (х, ?) ={0 < х < КО, 0 < ? < г}, С1)
д? дх
где Р(х, ?) - давление в пласте; % = к ; к - абсолютная проницаемость
цРт
пласта; ц - вязкость нефти; р - коэффициент объемного сжатия нефти; т - коэффициент пористости пласта; ¿О - положение границы раздела вода-нефть.
Пусть в начальный момент времени / = 0 распределение давления в нефтеносном пласте и положение границы раздела жидкостей известны, т.е. для уравнения (1) имеем следующие начальные условия
^(0) = Ь , (2)
Р?=0 =ф(х) , 0 < х < . (3)
Предположим, что изменение давления во времени на галерее эксплуатационных скважин описывается функцией /(?). Тогда на границе пласта х = 0 будем иметь следующее условие
Рх=0 = /(?) . (4)
На границе раздела вода-нефть х = s(t) давление нефти должно быть равно давлению краевой воды
Р|х=.*(0 = Ps ({) , (5)
и должно быть выполнено условие материального баланса
Ж к дР |
т17 =------*=*(0 . (6)
а? ц дх 1
Следует отметить, что прямая задача упруговодонапорного режима пласта, состоящая в нахождении функций Р( х), &•(?), удовлетворяющих уравнениям (1), (6) и заданным условиям (2) - (5), относится к задаче типа Стефана [2, 3].
Однако необходимо отметить, что в упруговодонапорном режиме пласта давление краевой воды не доступно для непосредственного измерения. Следовательно, функция рх (?) неизвестна и также подлежит определению. Очевидно, что для корректной постановки задачи необходимо задавать дополнительное условие. Пусть дополнительное условие задается на границе пласта х = 0
к£и = 9(0 . (7)
ц дх
Таким образом, требуется определить функции P( x, t), s(t), ps (t) из уравнения (1), (6) и условий (2) - (5), (7). Задачу (1) - (7) можно отнести к классу граничных обратных задач [4, 5]. Существенной особенностью данной задачи является наличие подвижной границы раздела жидкостей, закон перемещения которого определяется в ходе решения задачи. Однако граничные обратные задачи для Стефановской модели типа (1) -(7) мало исследованы [6, 7].
Метод решения. Используя метод выпрямления фронтов, преобразуем задачу (1) - (7). Путем замены переменных
область задания уравнения (1) ^ отобразим на область
О = {0 < у < 1, 0 < ї < Т}.
Тогда уравнения (1), (6) и условия (2) - (5), (7) принимают вид
Для перехода к разностной задаче введем равномерную разностную сетку в области О : гаИт = {(у,ї/): уі = іИ, ї/ = /с, И = Ь/Ы, т = Т/М,
і = 0,1,2,...,Ы, / = 0, 1, 2, ...,М}. Подставляя (8) в уравнение (7), разностные аналоги уравнений (7), (8) на сетке гаш запишем в виде
х
У = —, t = t, Р(х, t) = P(y, t) ,
s(t)
s 2(t) — = %-—p + ys (t)——, (y, t) ей= {0 < y < 1,0 < t < T}
() dt oy2 У () dt oy Vy’ ) 1 У ’ 1
(8)
pt=o = 9(y) , 0 < y < 1, s(0) = L,
py=0 = f (t) ,
(9)
(10)
(11)
-P і
— y=0 =Mt)q(t)/ к,
-y|y
py=1 = Ps (t) .
(12)
(13)
XpУі
h
h
і = 1, N -1, / = 0, М -1.
Полученную систему разностных уравнений преобразуем к виду
ар/Х - с,Р/+1 + Ьр/Х = -^ , і = 1,N -1, / = 0,М -1, (14)
рУ+І _ р .7+1
52/+1 = 52/ - 2тхР—--------— , (15)
И
где аі = х + хРУі (РЫ -РЫ-1); Ь =х; сі = аі + Ьі + 52/И2 /т; = 52] Рг]И2 /т.
К уравнениям (14), (15) добавим разностные аналоги начальных и граничных условий (9) - (13):
Р =Фі, (16)
Р0/+1 = //+1, (17)
Р+1 = Р+1 - Иц? V / к, (18)
РЫ+1 = /, (19)
50 = Ь. (20)
Для решения системы разностных уравнений (14), (16) - (19) используем подход, предложенный в [8]. Решение этой системы представляется в виде
Р/+1 = а^Р^1 + Рі+1, і = 0, 1, 2, ..., N -1, (21)
где аі+1 = Ьі /(сі -аіаі), Рі+1 = (аірі + )/(сі-а¡а¡), і = 1,2, ..., Ы-1,
а1 = 1, Р1 =-Ицц} 5 /к.
Для определения зависимости между Р/+ и Р+1 в явном виде соотношение (21) запишется при і = 0
Р0+1 = а1Р1/+1 +Р!.
Подставив сюда выражение для Р/ +1, т.е. Р/+1 = а 2 Р^741 + Р2, будем иметь
Р0+1 =а^ Р2+1 + а^ +Рі .
Далее, подставляя в последнее уравнение выражения для
Р+1,Р+1,...,Р/+1, получим формулу, в которой Р +1 выражается через
рУ+1
РЫ
N N і-1
Р0+1 = PN+1 Па і +ЕРі Па » +Р1.
і=1 і=2 и=1
ЮБМ 2079-0031 Вестник НТУ "ХПИ", 2012, № 62 (968)
29
Отсюда, учитывая (17), (19), получим
N г-1
I1 +1 -^Рг П“ п -^!
Р*1+1 =-------------------------------------------------^-• (22)
N
П “
г
г=1
Определив р1 по формуле (22), можно последовательно найти
Р^Х, Р^, Р/+1 по рекуррентной формуле (21). А положение границы раздела жидкостей можно определить как положительный корень уравнения (15). При переходе на следующий временной слой описанная процедура вычислений снова повторяется.
Таким образом, предложенный численный метод позволяет в каждом временном слое определить распределение давления в нефтеносном пласте и положение подвижной границы.
Результаты численных расчетов. Для выяснения эффективности предложенного вычислительного алгоритма были проведены численные эксперименты по следующей схеме: для заданных функций д(?), ¿•(Г) решалась прямая задача (7) - (9), (12). Найденная зависимость I(?) = Р(0, ?) принималась за точные данные для численного решения обратной задачи по восстановлению р8 (?). Первая серия расчетов выполнялась с использованием этих невозмущенных данных. Вторая серия расчетов проводились при наложении на I(?) некоторой функции, моделирующей погрешность экспериментальных данных I(?) = I(?) + 8(а(?) - 0.5), где ст(?) - случайный процесс, моделируемый с помощью датчика случайных чисел, 8 - уровень погрешности.
Расчеты выполнялись на пространственно-временной разностной сетке с шагами к = 0.04, т = 2 сут. Результаты численного эксперимента,
проведенного для случая р = 10-8 Па 1, к = 5 10 12 м2, ц = 0.003 Па -с,
Ь = 200 м, ¿(Г) = Ь - 2л/?, т = 0.3 , ф(х) = 200 атм., при использовании невозмущенных и возмущенных входных данных представлены в таблице; в ней ? - время; р/ - точные значения функции рх (?); рх -вычисленные значения рх (?) при невозмущенных данных; р -вычисленные значения ps (?) при возмущенных данных. Для возмущения входных данных использовалась погрешность уровня 8 = 5 атм.
Результаты численного эксперимента
t, сут. 10 20 30 40 50 60 70 80 90
P‘s 215.97 220.51 222.11 222.09 221.00 219.13 216.65 213.67 210.29
Ps 215.97 220.51 222.11 222.09 221.00 219.13 216.65 213.67 210.29
Ps 214.94 220.20 219.16 220.82 219.20 217.53 218.46 210.58 211.26
Как показывают результаты численного эксперимента, при использовании невозмущенных входных данных искомая функция ps (t) восстанавливается точно. При использовании возмущенных входных данных, в которых погрешность имеет флуктуационный характер, искомая функция p (t) восстанавливается с относительной погрешностью приблизительно в 2 %. При уменьшении уровня погрешности решение восстанавливается более точно.
Выводы. Предложена математическая модель упруговодонапорного режима разработки пласта и алгоритм ее численного решения.
Список литературы: 1. Басниев К.С. Нефтегазовая гидромеханика / К.С..Басниев, Н.М. Дмитриев, Г.Д. Розенберг. - М. Ижевск: Институт компьютерных исследований, 2005.
- 544 с. 2. Рубинштейн Л.И. Проблема Стефана / Л.И. Рубинштейн. — Рига, Звайгзне, 1967. -232 с. 3. Самарский А.А. Вычислительная теплопередача / А.А. Самарский, П.Н. Вабищевич.
— М.: Едиториал, 2003. - 784 с. 4. Самарский АА. Численные методы решения обратных задач математической физики / А.А. Самарский, П.Н. Вабищевич. - М.: Издательство ЛКИ, 2009. - 480 с. 5. Алифанов О.М. Экстремальные методы решения некорректных задач / О.М. Алифанов, Е.А. Артюхин, С.В. Румянцев. — М.: Наука, 1988. - 288 с. 6. Голдьман Н.Л. Классическое и обобщенное решение двухфазной граничной обратной задачи Стефана / Н.Л. Голдьман // Вычислительные методы и программирование. - 2002. - № 3. - С. 133143. 7. Лобанова В.В. Решение одной обратной задачи Стефана с линейной подвижной границей / В.В. Лобанова // Математический журнал. Алматы. - 2006. - № 2. - С. 76-80. 8. Гамзаев Х.М. Численное решение некорректной задачи однофазного течения в двумерном пласте / Х.М. Гамзаев // Проблемы управления и информатики. - Киев. - 2008. - № 5. -С. 74-84.
Статью представил д.т.н., проф., зав. кафедрой систем информации НТУ ’ХПИ" А.А. Серков.
УДК 532.546:519.63
Чисельне моделювання пруговодонапірного режиму пласта / Гамзаев Х.М.
// Вісник НТУ "ХПІ". Серія: Інформатика та моделювання. - Харків: НТУ "ХПІ". - 2012. -№ 62 (968). - С. 26 - 32.
Для опису пруговодонапірного режиму розробки пласта пропонується модель типу Стефана, яка характеризується наявністю невідомого кордону розділу рідин. Умова на межі розділу вважається невідомою і замість цього задається додаткова умова на фіксованому кордоні пласта. Для чисельного розв’язання поставленої задачі запропонований обчислювальний алгоритм на основі методів випрямлення фронтів і різницевої апроксимації. Табл. 1. Бібліогр.: 8 назв.
Ключові слова: пруговодонапірний режим, пласт, метод випрямлення фронтів, різницево апроксимація.
UDC 532.546:519.63
Numerical modeling of elastic water drive mode of layer / Gamzaev Kh.M. // Herald of the National Technical University "KhPI". Subject issue: Information Science and Modelling. -
Kharkov: NTU "KhPI". - 2012. - №. 62 (968). - P. 26 - 32.
For the description of water drive development mode of the reservoir proposed model of the type of Stephen, which is characterized by the presence of unknown border of section of liquids. The condition on the boundary is considered to be unknown, and instead set additional condition on the fixed boundary layer. For the numerical solution of the set task proposed numerical algorithm based on the methods of rectification of the fronts and difference approximation. Tabl.: 1. Refs.: 8 titles.
Keywords: elastic water drive mode, the layer, the method of straightening fronts, difference approximation.
