строят карты распределения параметров пластов по площади (карты равных толщин пласта — изопахит, проницаемостей, по ристостей, вязкостей нефти; карты распространения зональных интервалов или слоев и др.) и схемы распределения по раз резу, строят блок-диаграммы. При этом важная роль отводится детальной корреляции продуктивного пласта, что позволяет точно определить толщины пластов, детально расчленить раз рез, выяснить прерывистость пластов по простиранию и посто янство свойств слагающих его пород, выделить отдельные слои (пропластки) и зоны (линзы).
Вероятностный метод обоснован следующим. Так как доступ в залежь имеется посредством скважин, то по результатам про водимых исследований скважин и пластов локально (в извест ных размерах области пласта) устанавливают параметры за лежи. Значения этих параметров изменяются по объему (пло щади, толщине) пласта в широких пределах. Изменения, можно полагать, носят случайный характер. Поэтому считают, что дан ные исследований — это выборка из генеральной совокупности параметров (залежи), что позволяет использовать аппарат ма тематической статистики, теории вероятностей и теории случай ных функций. В соответствии с законом больших чисел при уве личении объемов выборки выборочные характеристики сходятся по вероятности к генеральным, т. е. их надежность возрастает, погрешность уменьшается.
Такая обработка данных исследований позволяет построить модель пласта. В зависимости от метода получают детерминиро ванную или вероятностную модель пласта.
§ 2.2. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РАЗРАБОТКИ МЕСТОРОЖДЕНИИ
Процесс разработки конкретного нефтяного месторождения од нократен и непосредственному наблюдению «доступен» только в ограниченном числе скважнно-точек. Это обусловило разви тие методов его моделирования, как непрямых, опосредствован ных методов научного исследования. Моделирование, и как ре зультат, модель процесса обеспечивают возможность при срав нительно небольших затратах в короткие сроки многократно (многовариантно) «проиграть» медленно протекающие процессы разработки в различных технологических условиях и тем самым выбрать рациональную технологию. При создании моделей процесса разработки нефтяных месторождений моделируют ге олого-физические свойства пласта, его геометрическую форму, флюиды и процесс извлечения нефти и газа из недр.
Различают физическое и математическое моделирования. При физическом моделировании на модели, представляющей по су ществу натурный или масштабно уменьшенный образец ориги нала (лабораторную, пилотную установки), воспроизводят и ис следуют процессы, качественно одинаковые с процессами, про текающими в реальном объекте. В связи с трудностью создания полного подобия пласта и измерения параметров гидравличе ские модели нефтяных пластов не нашли применения, хотя фи зическое моделирование отдельных элементов процесса разра ботки незаменимо (например, вытеснение нефти водой).
Математическое моделирование заключается в исследовании процессов путем построения и решения системы математиче ских уравнений, относящихся к собственно процессу и краевым условиям. Математическая модель основана на упрощении ^идеализации) сложного реального процесса. Для ее создания природные условия соответствующим образом дифференцируют, выделяют среди них главные, определяющие факторы и пред ставляют HXJ в таком виде, который обеспечивает возможность достижения цели. Причем нефтегазоносный пласт рассматри вают как единую гидродинамически связанную систему не только во всей области нефтегазоносности, но и включая окру жающую водонапорную область. Перемещение флюидов внутри этой единой системы определяется начальными (до начала раз работки) и граничными (на поверхностях, ограничивающих пласт с внешних сторон, и на стенках скважин внутри пласта) условиями или в совокупности краевыми условиями.
Граничные условия задаются в виде искомой -функции (дав ление, расход жидкости), ее производной (градиент давления, скорость) или в смешанном виде (соответственно граничные ус ловия первого, второго или третьего рода).
Начальные условия характеризуют состояние пласта до на чала разработки, когда в процессе разработки неустановившееся (нестационарное) движение флюидов наряду с простран ственными координатами определяется еще временем.
Системы математических уравнений решают аналоговым и вычислительным методами. Аналоговый метод математического моделирования базируется на подобии явлений и процессов раз личной физической природы, т. е. на широкой физической ана логии. Можно назвать аналогии между полями фильтрации жидкости (закон Дарси), электрического тока в проводящей среде (закон Ома), электрическим в диэлектрике (закон индук ции), магнитным (закон магнитной индукции) и температурным (основное уравнение теплопроводности).
Электрическое моделирование процесса разработки основано на электрогидродинамической аналогии (ЭГДА), т. е. аналогии
Между .движением электрического тока в проводящей среде и фильтрацией жидкости в пористой среде.
Вычислительные методы подразделяются на аналитические, численные и статистические. Аналитические методы соответ ствуют классическому подходу к моделированию процессов, когда ставится исходная задача, вводятся упрощающие пред положения и на их основе формулируется новая задача, кото рая поддается решению в виде аналитического выражения, фор мулы, обеспечивающей получение значения функции для каж дого значения аргумента. Упрощающие предположения иногда приводят к существенным погрешностям в результатах проекти рования, а без них задача в аналитической форме не решается. К числу аналитических методов, дающих точные решения задач разработки нефтяных месторождений, т. е. в точности удовлет воряющих исходным уравнениям, начальным и граничным ус ловиям, относятся метод разделения переменных ■ (метод Фурье), методы теории функций комплексного переменного, ин тегральных преобразований и др. Приближенные решения по лучают с использованием методов эквивалентных фильтрацион ных сопротивлений, последовательной смены стационарных со стояний, интегральных соотношений и др.
Погрешность приближенных методов по сравнению с точ ными оценена для некоторых случаев в подземной гидрогазо динамике, и мы не будем останавливаться на ней. Отметим только, что с учетом точности исходной информации и сложно сти изучаемого процесса она в отдельных случаях вполне удов летворяет практику.
Более полный учет множества воздействующих на процесс разработки факторов можно выполнить с использованием чис ленных методов на базе применения ЭВМ как наиболее эффек тивных и универсальных моделирующих устройств. Основы при менения численных (разностных) методов математического моделирования для решения задач разработки нефтяных место рождений в настоящее время достаточно развиты. При проекти ровании разработки нефтяных месторождений чаще всего при
меняют численные методы.
Статистические методы моделирования базируются на ста тистических данных предшествующей разработки месторожде ний. Изучая фактические закономерности развития процесса в залежи за прошедший период (ретроспективу), они позволяют оперативно без больших‘затрат времени и труда сформулиро вать заключение о предстоящем развитии основных технологи ческих показателей разработки (перспективу). Эти методы из ложены в гл. 4.
Электромоделирование процесса разработки нефтяных место рождений осуществляется с использованием электрических мо делей-аналогов. Электрическая модель (электроинтегратор) мо-. жет быть: 1) сплошной средой — жидкой (в-виде электролити ческой ванны) или твердой (в виде листов электропроводящей бумаги или фольги различных металлов); 2) сеткой дискрет ных элементов — омических сопротивлений (R-сетки) или оми ческих сопротивлений и емкостей (RC-сетки); 3) различной ком бинацией первых двух. Предпочтение отдают RC-сеткам.
Вобщем случае для исследования фильтрационных потоков
взалежах (объемных пластах) электрические сетки должны быть трехмерными (пространственными). Толщина пласта и ее изменение очень малы по сравнению с его размерами в горизон тальной плоскости вдоль осей х и у, поэтому вертикальной со ставляющей скорости движения можно пренебречь. Это позво
ляет ограничиться созданием плоских сеток и, следовательно, решением приближенных двухмерных уравнений движения для весьма протяженных пластов с переменной толщиной h= h(x,y). Например, такое уравнение при упругом режиме (типа Фурье) применительно к неоднородному пласту можно записать
д |
Г k (х, у) h (х, у) |
др 1 |
.__ д |
Г k (х, |
у) h (х, у) |
др 1 _ |
|
дх |
L |
р |
дх J |
ду |
L |
Р |
ду J |
|
|
|
= h(x, |
У) ?>*-%-, |
|
(2.1) |
|
|
|
|
|
|
at |
|
|
где k, h — проницаемость и толщина пласта; р — абсолютная (динамическая) вязкость жидкости; р* —коэффициент упругоемкости насыщенного пласта; р — давление; / — время.
Допустим, требуется определить изменение пластовых и за бойных давлений в замкнутой залежи при следующих краевых (начальном и граничных) условиях:
|
|
р = р(х, |
у) — рк при / = 0; |
|
(2.2) |
||
Qt = |
$ |
k(x,y)h(x, y) |
др_ЛТ1> t. = |
1( |
2( |
. , |
п; (2.3) |
Р |
дпх |
|
|
||||
|
|
|
|||||
|
|
|
|
|
|
||
|
|
др/дп2= О, |
|
|
|
(2.4) |
|
где рк— начальное пластовое давление |
в |
залежи; |
qi — дебит |
||||
t-й скважины; п — число |
работающих |
скважин; пи |
п2— нор |
||||
мали к контурам Г, скважин и к внешней непроницаемой гра нице залежи.
Электрические процессы в сетке, образованной омическими сопротивлениями R и включающей емкости С0, в любой момент времени описываются аналогичным уравнением
дхэ ч Rx дхэ ) |
дуэ Ч Ry дуэ ) |
0 dt3 |
(2.5) |
|
где ха, Уэ— координаты сеточной области электрической модели; U — электрическое напряжение; Rx, Rv— сопротивление элемен тов электрической сети; Со— емкость конденсатора; /э— время протекания электрических процессов.
Конденсаторы одними концами присоединяются к узлам сетки, а другими — к одной общей точке. Конденсаторы перио дически заряжают на напряжение U0.
Для любого элементарного объема пласта AxAyh в направ лении оси л: и для соответствующего элемента электрической сетки можно записать объемный расход жидкости qx по закону Дарси и электрический ток ix по закону Ома:
q = k (х, у) h (х, |
у) Ду Ар |
_ |
k (х, у) h (х, у) д |
^ б) |
|Х |
Дх |
|
|А |
|
|
ix = AU/Rx, |
|
(2.7) |
|
где Ах=Ау.
Для пересчета электрических величин в фильтрационные вводят коэффициенты пропорциональности (подобия):
Cq = Uq-, |
CR =R /—£~; Ср— W p\ |
4 |
kh |
Се= С0//ф*; Ct - |
tslt\ Cx = Cy = x jx = yJy=\IM „, (2.8) |
где Мм— линейный масштаб модели.
Пересчитывая электрические величины в уравнениях (2.5) и (2.7) с помощью коэффициентов пропорциональности (2.8) и сопоставляя их соответственно с уравнениями (2.1) и (2.6), оп ределяют условия подобия протекания фильтрационных и элек трических процессов:
СЛСР/(С,М„) = 1; CPI{C„CR) = 1. |
(2.9) |
Таким образом, выполняя условия подобия (2.9), на элект рической модели с RC-сеткой решают задачу неустановившейся фильтрации упругой жидкости к скважинам. Задачи неустано вившейся фильтрации реального газа, встречающиеся при мо делировании разработки нефтегазовых залежей, можно решить по методике, предложенной Ю. П. Коротаевым, С. Н. Закиро
вым и А. Н. Тимашевым.
Основными недостатками аналоговых сеточных вычислитель ных машин являются необходимость трудоемкого перебора со противлений при изменении условий разработки (например,
45