Материал: Разработка и эксплуатация нефтяных месторождений.

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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

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

Такая обработка данных исследований позволяет построить модель пласта. В зависимости от метода получают детерминиро­ ванную или вероятностную модель пласта.

§ 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

Источник: https://tut-files.ru/previewfile/77200