УП Мусин 2
.pdf
|
|
|
Министерство образования и науки Республики Татарстан |
НИ |
|||||||||||||
|
|
|
|
||||||||||||||
|
|
|
|
Альметьевский государственный нефтяной институт |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
М.М. Мусин |
|
|
|
ка |
АГ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
о |
т е |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Разработка нефтяных |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
месторожден й |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
б |
л |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Часть II |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
Учебное пособие |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
для студентов, обучающихся по специальностии |
130503.65 «Разработка |
|||||||||||||||
|
и эксплуатация нефтяных и газовых месторождений», всех форм обучения |
||||||||||||||||
|
|
|
|
|
|
н |
н |
ая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
тр |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
е |
к |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Э |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
Альметьевск 2007 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
УДК 622.276 |
|
|
|
|
|
|
|
|
|
|
НИ |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Мусин М.М. |
|
|
|
|
|
|
|
|
|
АГ |
|
|||
|
|
|
|
Разработка нефтяных месторождений. Часть II: Учебное пособие |
||||||||||||
|
|
для студентов, обучающихся по специальности 130503.65 «Разработка и |
||||||||||||||
|
|
эксплуатация нефтяных и газовых месторождений. – Альметьевск: Аль- |
||||||||||||||
|
|
метьевский государственный нефтяной институт, 2007. – 140 с. |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ка |
|
|
|
|
Данное учебное пособие составлено на основе лекций, прочитанных ав- |
||||||||||||||
|
тором |
студентам Альметьевского государственного нефтяного института в |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
е |
|
|
|
|
2002-2007 г.г. В нем рассмотрены стратегические задачи разработки нефтяных |
|||||||||||||||
|
и нефтегазовых месторождений, модели нефтяных пластов и современные |
|||||||||||||||
|
представления механизма взаимного вытеснения неф и и воды в нефтяных пла- |
|||||||||||||||
|
стах. Изложены методики расчета технологических показателей, методы кон- |
|||||||||||||||
|
троля и регулирования, правовые основы проектир ваният |
и разработки нефтя- |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
ных месторождений. Рассмотрены наиболее широко применяемые на практике |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
гидродинамические, физико-химические и тепловыео |
методы повышения нефте- |
||||||||||||||
|
отдачи пластов. Даны основные понятия о чис енном методе решения плоской |
|||||||||||||||
|
задачи вытеснения нефти водой в системе скважин |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
и |
|
|
совета АГНИ. |
|
|
|
|
|
Печатается по решению учебно-методическогоб |
|
|
||||||||||||
|
Рецензенты: |
|
|
|
ая |
б |
|
|
|
|
|
|
|
|||
|
|
1. |
Липаев А.А. – зав. кафедрой «Разработки и эксплуатации нефтяных и га- |
|||||||||||||
|
|
|
зовых месторождений», д.т.н., профессор |
|
|
|
|
|
||||||||
|
|
2. |
Фазлыев Р.Т. - гл вный научный сотрудник ТатНИПИнефть, действи- |
|||||||||||||
|
|
|
|
|
|
|
н |
|
|
|
|
|
|
|
|
|
|
|
|
тельный член РАЕН, профессор, д.т.н. |
|
|
|
|
|
||||||||
|
|
|
|
тр |
о |
н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
е |
к |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|
© Альметьевский государственный |
|||||||
Э |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
нефтяной институт, 2007 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
Э
ВВЕДЕНИЕ
Программа курса «Разработка нефтяных месторождений» охватыва-
ет следующие основные разделы: |
|
|
|
|
АГ |
|
|
– |
схематизация залежи и создание ее геолого-физической модели; |
||||||
– |
выбор методов воздействия на нефтяные залежи; |
|
НИ |
||||
– |
выделение эксплуатационных объектов; |
|
|
|
|
||
– |
обоснование системы разработки; |
|
|
|
ка |
|
|
– |
расчет технологических показателей разработки; |
|
|
||||
– |
анализ состояния и регулирование процесса р зр ботки. |
|
|||||
В результате усвоения дисциплины студент должен знать: |
|
||||||
– |
|
|
|
е |
|
|
|
источники пластовой энергии и режимы работы нефтяных пла- |
|||||||
|
стов; |
|
|
|
|
|
|
– |
теоретические основы проектирования разработки нефтяных ме- |
||||||
|
сторождений; |
|
о |
|
|
|
|
– |
основные методы увеличения нефте тдачитпластов; |
|
|
||||
– |
|
и |
|
|
|
|
|
методы составления проектных документов и способы планиро- |
|||||||
|
вания мероприятий по регулирован ю выработки запасов нефти. |
Разработкой нефтяных месторождений называют осуществление на- |
||||
|
|
|
б |
|
учно обоснованного процесса извлечения из недр нефти и газа. Процесс |
||||
|
|
и |
|
и выработки запасов |
разработки состоит из разбуривания месторождениял |
||||
нефти и газа. |
б |
|
|
|
Курс “Разработка нефтяных месторождений“– область прикладной науки, как учебная дисциплина относится к инженерным дисциплинам. В
ней дается не только качественное описание месторождения, но и количе-
вания и проектиров нияаяр зработки, реализации проектных решений, контроля, анализа и регулированиян разработки нефтяных месторождений.
ственные характеристики процесса извлечения нефти и газа из пласта.
Она состоит из разделов о системах и технологиях разработки, планиро-
Теоретической ос овой разработки нефтяных месторождений являются фундаме тальные положения физики нефтяного пласта и подземной
гидромеханики. |
|
|
|
Разраб тку нефтяныхн |
месторождений трудно себе представить без |
|
тр |
|
геологии и ге физики, без изучения геологического строения месторож- |
||
дения. Она опиоается также на важные положения физики и химии, ма- |
||
тематической физики и технологии эксплуатации скважин. |
||
|
Прое тирование и осуществление разработки нефтяных месторож- |
|
|
е |
|
дений подчинены единой цели и требуют единой методики, которая по- |
||
зволя тксвязать все знания о нефтяной залежи и о процессах, происходя- |
||
л |
|
|
щих внутри нее при извлечении жидкости и газа. Отдельные положения науки о разработке нефтяных месторождений были обоснованы в трудах И.М Губкина, Л.С. Лейбензона, А.П. Крылова, С.А Христиановича, И.А.
3
Э
Чарного, В.Н. Щелкачева, М. Маскета, М. Леверетта и др. уже в 30-40х годах ХХ века.
Курс «Разработка нефтяных месторождений» как самостоятельная
область науки и учебная дисциплина возникла после выхода в свет |
рабо- |
|
|
АГ |
|
ты Крылова А.П. с соавторами «Научные основы разработки нефтяных |
||
месторождений» в 1948 г. |
|
НИ |
Наука о разработке нефтяных месторождений является одной из |
немногих прикладных наук, имеющих дело с системой, которую в целом нельзя ни увидеть, ни взвесить, ни измерить, ни исследовать. Даже на месторождениях, где каждая скважина пробурена с полным отбором керна, отобранные образцы породы (керн) составляет ничтожно малую долю породы пласта, которая может быть исследована напрямую. Так же ограниче-
тоды каротажа скважин сравнительно недороги и могут часто использо-
ны и пробы жидкостей, на которых производятся д тальные лабораторные |
|||
исследования. |
|
е |
ка |
|
|
||
|
|
|
|
Ввиду указанной ограниченности объема образцов породы, были |
|||
разработаны методы геофизических исследований тскважин (каротажа) и их |
|||
интерпретации для получения более детальной |
о |
нформации о пласте. Ме- |
|
и |
|
|
ваться и на старых скважинах, пробуренных много лет назад, для опреде- |
||||
|
|
|
б |
|
ления значений пористости, проницаемости, насыщенности и характера |
||||
литологической изменчивости. |
|
и |
|
л |
Методами каротажа измеряют электрические, электромагнитные, аку- |
||||
|
б |
|
|
|
стические или радиоактивные свойства пород в пределах радиуса, в несколько раз превышающего радиус скважины. Для определения характеристик пласта в более крупном масштабе используются исследования кривых изменения давления и их интерпретация. Еще более новыми являются методы комплексной обработки геологической информации, когда
данные, полученные по керну, каротажным диаграммам и шламам породы, |
||||
|
|
|
н |
|
откалываемым долотом, обр батывается статистически для получения де- |
||||
|
|
н |
|
|
тального описания измеаяе ия фильтрационно-емкостных свойств коллек- |
||||
тора по пласту. |
о |
|
|
|
|
Нефтедобывающее производство состоит из следующих основных |
|||
направлений: |
|
|
|
|
1. |
тр |
|
|
|
проведение разведки и обоснование запасов нефти и газа; |
||||
2. |
проекти ование разработки; |
|||
3. |
обус ройс во месторождения и эксплуатационное разбуривание место- |
|||
|
к |
|
|
|
|
рождения; |
|
|
|
4. |
осуществление мероприятий по обеспечению притока нефти и газа из |
|||
л |
пласта к забоям скважин; |
|||
5. |
подъ м продукции из забоя скважин на дневную поверхность: |
|||
6. |
внутрипромысловаяе |
транспортировка и подготовка добываемой про- |
дукции.
4
Э
В курсе «Разработка нефтяных месторождений» излагаются вопроНИ- сы, в основном, касающиеся второго и четвертого вышеприведенных направлений.
Данное пособие, которое состоит из двух частей, подготовлено на основе материалов лекций, прочитанных автором в АльметьевскомАГ государственном нефтяном институте в 2002-2007 годах студентам дневного и заочного отделений, обучающимся по специальности «Разработка и эксплуатация нефтяных и газовых месторождений».
Глава 1 |
|
|
е |
|
|
т |
|
ЧИСЛЕННАЯ МОДЕЛЬ МНОГОМЕРНОГО ТЕЧЕНИЯка |
|||
ДВУХ НЕСМЕШИВАЮЩИХСЯ ЖИДКОСТЕЙ |
|||
|
о |
|
|
В главе 5 первой части учебного пос бия была рассмотрена одно- |
мерная модель непоршневого вытеснен я нефти водой из пласта, для ко-
торой получено аналитическое решение. Однако область применения од- |
|
номерной модели для расчета процесса разработкии |
нефтяных месторож- |
дений с применением заводнения довольно ограничена. Эту модель мож- |
|
|
л |
но использовать лишь для приближенных расчетов применительно к ряд- |
|
ным системам заводнения однородныхб пластов. Вместе с тем известно, |
|
|
б |
что нефтяные пласты обладают с льной послойной и площадной неодно- |
|
родностью. В процессе разработкиичасто возникает необходимость уплот- |
|
ая |
|
нения сетки и ввода в разработку новых скважин, перевода отдельных добывающих скважин под закачку воды, отключения обводненных скважин. В целях регулирования процесса разработки пласта и повышения коэффициента нефтеотдачи широко распространена технология целенаправлен-
ного изменения аправле ия фильтрационных потоков нефти и воды в |
||||
|
|
|
|
н |
пласте. По этой причи е движение нефти и воды в пласте является суще- |
||||
ственно неодномер ойн, поэтому возникает необходимость использования |
||||
|
|
|
о |
|
в расчетах двухтрехмерных моделей совместного движения многофазной |
||||
|
|
тр |
|
|
(нефти, воды и газа) жидкости в пористой среде. |
||||
|
Основой этих моделей являются законы сохранения массы и Дарси. |
|||
|
к |
|
|
1.1. Закон Дарси |
|
е |
|
|
|
л |
Нефть и вода, являющиеся несмешивающимися жидкостями, отли- |
чаются друг от друга по вязкости и относительным фазовым проницаемостям. Совместное движение их в пористой среде подчиняется линейному закону Дарси, который в векторной форме записывается в виде:
5
|
|
|
|
|
|
|
|
|
|
|
|
в = − |
|
кhкв |
gradP , |
|
|
|
|
|
|
|
|
НИ |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ν |
|
|
|
|
|
|
(1.1) |
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
μв |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
н = − |
|
кhкн |
|
gradP , |
|
|
|
|
|
|
|
АГ |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
ν |
|
|
|
|
|
|
(1.2) |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
μн |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂P |
|
|
|
|
|
|
∂P |
|
|
|
|
|
∂P |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
gradP = |
|
|
|
+ |
|
|
|
|
+ |
|
|
. |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
где |
i |
j |
к |
|
|
|
(1.3) |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
∂у |
|
∂z |
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂х |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь к, h - абсолютная проницаемость и нефтен сыщенн я толщина |
|||||||||||||||||||||||||||||||||||||||||||||||
|
пласта соответственно; кв= кв(s), кн= кн(s) – относительные ф зовые про- |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
ницаемости воды и нефти, являющиеся функциями от насыщенности во- |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
дой s; µв, µн – динамическая вязкость воды и н фти; р- пластовое давле- |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ка |
|
|
|
ние; ν в , ν н - вектора объемной скорости воды и н ф и. |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
Уравнение (1.1) можно представить в виде: |
т |
е |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
о |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
ν в =vвхi+ vву j , |
|
|
|
|
|
|
(1.4) |
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
где vвх, vву, - проекция скорости воды на оси х и у. |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
∂p |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
khkв |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
Vв у = Vву (x, y,t) = − |
б |
|
∂у |
; |
|
|
|
(1.5) |
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
μв |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
khkв ∂p |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
vвх |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= − |
|
|
|
|
|
|
|
|
. |
|
(1.6) |
|
|
|||||||||||
|
|
|
|
|
|
|
|
= Vв х (x, y,t) |
|
μв |
|
∂x |
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
Аналогично можно записать и уравнение (1.2) для нефти. |
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
1.2.Уравнение неразрывности потока |
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
Уравнение |
|
|
|
|
ая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
еразрывности при совместном движении двух фаз без |
|||||||||||||||||||||||||||||||||||||||||||||||
|
учета сжимаемости породын |
и жидкостей записывается для каждой фазы в |
||||||||||||||||||||||||||||||||||||||||||||||||
|
отдельности: |
о |
н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
тр |
|
|
|
|
|
|
|
|
|
|
∂ν ву |
|
|
|
|
|
|
|
|
|
|
|
∂sв |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
∂ν вх |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
к |
|
|
∂x + ∂y |
|
|
= −mh ∂t , |
|
(1.7) |
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
е |
|
|
∂νнх |
+ |
|
∂ν ну |
= −mh |
∂sн |
, |
|
|
|
(1.8) |
|
|
|||||||||||||||||||||||||||||||||
|
л |
|
|
|
∂х |
|
|
|
|
|
|
∂у |
|
|
∂t |
|
|
|
|
|
||||||||||||||||||||||||||||||
Э |
|
|
|
|
|
|
|
|
|
|
sв+sн=1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.9) |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Э
|
|
|
mh |
|
∂t |
+ |
∂x (Vx F(s))+ |
|
∂y |
(Vy F(s)) |
= 0 |
|
(1.10) |
НИ |
||||||||||||||||||||||
|
Система уравнений (1.1), (1.2), (1.7)-(1.9) описывает совместное |
|
||||||||||||||||||||||||||||||||||
движение нефти и воды в пористой среде. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
Уравнение (1.7)-(1.9) можно привести к следующему виду: |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
∂s |
|
∂ |
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
АГ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
∂Vx (x, y,t) |
+ |
|
∂Vy |
(x, y,t) |
= 0 |
|
|
|
(1.11) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
ка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Здесь |
|
|
|
|
|
Vx (x, y, t) = −c(x, y, t) |
∂Р |
; |
|
е |
(1.11а) |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂х |
|
т |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
V у(x, y, t) = −c(x, y, t) |
∂Р |
; |
|
|
|
(1.11б) |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
∂у |
|
|
о |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
к |
(s) |
|
|
|
к |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
(s) |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
c(x,y,t)= kh( |
|
|
|
+ |
|
|
|
н |
|
|
|
) . |
|
(1.12) |
|
|
||||||||||||
|
|
|
|
|
|
|
|
μв |
|
|
|
|
μн |
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|||||||||||
|
Для постановки конкретной задачи фи ьтрации нефти и воды в по- |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ристой среде кроме уравнений (1.10)-(1.11) необходимо сформулировать |
||||||||||||||||||||||||||||||||||||
граничные и начальные условия. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
1.3. Гран чные условия |
|
|
|
|
|||||||||||||||||||||||
|
Если граница пласта непроницаемаи |
, то граничные условия для дав- |
||||||||||||||||||||||||||||||||||
ления имеют вид: |
|
|
|
|
|
ая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
∂Р |
|
= 0 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
(1.13) |
|
|
||||||
|
|
|
|
|
|
|
н |
|
|
|
∂n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где n- направление |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
ормали к границе. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
На проницаемой границе надо задавать либо давление |
|
|
|||||||||||||||||||||||||||||||||
|
|
тр |
о |
|
|
|
|
|
|
Р г = P(x, y,t), |
|
|
|
|
|
|
|
|
|
(1.14) |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
либо расход жидк сти |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
к |
|
|
|
|
|
|
|
|
v = Vn (x, y,t). |
|
|
|
|
|
|
(1.15) |
|
|
||||||||||||||||
|
Граничные условия по насыщенности задаются в зависимости от то- |
|||||||||||||||||||||||||||||||||||
|
е |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
го, вте ает или вытекает вода через эту границу. |
|
|
|
|
|
|||||||||||||||||||||||||||||||
л |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если вода вытекает через границу, то граничные условия по насыщенности не требуется задавать.
7
Э
|
|
|
|
|
|
|
НИ |
Если вода втекает через границу, то должно быть задано значение |
|||||||
насыщенности на этой границе. |
|
|
|
|
|
|
|
Кроме границ залежи необходимо задавать условия на скважинах. |
|||||||
Они могут иметь вид: |
|
|
|
|
|
|
|
∙ |
на скважине задано забойное давление Рскв |
|
|||||
∙ |
на скважине задан дебит жидкости qскв |
|
|||||
Для добывающей скважины задание насыщенности не требуется. |
|||||||
|
|
|
|
|
|
ка |
|
На нагнетательной скважине обычно задают максимальноеАГзначение |
|||||||
водонасыщенности, которое входит в относительную ф зовую проницае- |
|||||||
мость. |
|
|
|
|
|
|
|
|
1.4. Начальные условия |
|
|
|
|
||
При t=0 задаются начальное распределение давл ния и водонасы- |
|||||||
щенности во всех точках пласта. |
|
|
т |
е |
|
|
|
1.5. Формулировка постановкио |
|
|
|
||||
задачи |
|
||||||
Пусть однопластовая залежь нефти д иной Lx , шириной Ly разраба- |
|||||||
|
|
|
и |
|
|
|
|
тывается с Nс скважинами в режиме заводнения. Система расположения |
|||||||
нагнетательных и добывающих скважинлпроизвольная с известными ко- |
|||||||
|
и |
б |
|
|
|
|
|
ординатами; пласт по проницаемости, пористости и нефтенасыщенной |
|||||||
толщине неоднородный. Требуется определить распределение насыщен- |
ности водой (и нефтью) и давленияб в каждой точке пласта для любого момента времени, удовлетворяющие вышеприведенным уравнениям (1.10), (1.11), начальным и граничным условиям (1.13)-(1.15). По результатам расчета давления и насыщенности водой необходимо рассчитывать суточные, месячные и годовые технологические показатели разработки.
Nх* Nу.
|
Эта задача е имеет аналитического решения, поэтому применяют |
||||||
|
|
|
|
|
|
|
ая |
численный методметод конечных разностей /17/. |
|||||||
|
|
|
|
|
|
н |
|
|
|
|
|
1.6.Численный метод решения плоской |
|||
|
|
|
|
|
н |
|
|
|
|
|
|
задачи фильтрации двухфазной жидкости |
|||
|
|
тр |
о |
|
|
|
|
|
|
|
|
|
1.6.1. Построение сетки |
||
|
|
|
|
|
|
||
|
|
1.Пласт разбивается на прямоугольные ячейки. Центрам ячеек при- |
|||||
сваиваютсяк |
номера: i – номер столбца, j – номер строки. Количество ячеек |
||||||
л |
е |
|
|
|
|
|
|
|
|
|
|
|
|
|
Э
|
НИ |
3.В ячейки – скважины заносятся заданные параметры пласта: про- |
|
ницаемость, пористость, нефтенасыщенная толщина. |
|
4.В межскважинных ячейках эти параметры вычисляются методом
интерполирования. Таким образом, получают сеточные карты параметров |
||||||||||||||||
пласта. |
|
|
|
|
|
|
|
|
|
|
|
|
|
АГ |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
5. Функции Р(x,y,t) и S(x,y,t) от непрерывных аргументов заменяют- |
||||||||||||||||
ся на функции pi j и s i j |
от дискретных аргументов (i, j) |
ка |
|
|||||||||||||
Pi j n =P(xi, yj, tn), где n – номер временного слоя. |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i,j+1 |
|
|
|
|
|
|
е |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i-1,j |
|
|
i j |
|
|
|
|
о |
|
т |
|
|
|
|
|
|
|
|
|
|
|
i+1,j |
|
|
|
|
|
||||
|
|
|
|
|
i,j-1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i →л |
и |
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
↑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
Рис. 1.1. Шаблон для выч сления потоков жидкости |
|
||||||||||||||
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ая |
в сеточной области |
|
|
|
|
|
|||||||
|
|
|
|
|
|
(1.11)-(1.12) для каждой ячейки, |
||||||||||
6.Дифференциальные уравнения |
кроме ячеек со скважинами, заменяются разностными уравнениями. Таким образом, получ ют систему Nх* Nу – Nс алгебраических уравнений.
Решив эту систему, находят значения давления и водонасыщенности для |
||||||
каждой ячейки для зада ого момента времени t = t n. |
||||||
|
Затем для этогонмомента времени вычисляются дебиты жидкости, |
|||||
нефти, величина бводненности каждой скважины. |
||||||
|
|
|
|
|
|
н |
|
Далее все расчеты повторяются для следующего шага по времени t |
|||||
=tn+∆t. |
|
|
|
о |
|
|
|
|
|
|
|
||
|
|
|
1.6.2. Построение разностных уравнений для давления |
|||
|
|
|
тр |
|
|
|
|
Дифференциальное уравнение (1.10) представляет собой уравнение |
|||||
|
|
к |
|
|
|
|
баланса (сохранения) жидкости (нефть+вода). Поэтому для получения его |
||||||
разностногое |
аналога для ячейки (i j) достаточно вычислить потоки жидко- |
|||||
л |
|
|
|
|
|
|
|
|
|
|
|
|
9 |
Э
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
НИ |
сти через границы ячейки, затем их алгебраически просуммировать. По- |
||||||||||||||||||
лучим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cn |
( pn |
− pn |
) −cn |
i−1/ 2, j |
( pn |
− pn |
i−1, j |
) + |
|
|
|
|||||||
i+1/ 2, j |
|
i+1, j |
|
i, j |
|
|
|
i, j |
|
|
|
|
|
АГ |
|
|||
+ cn |
|
( pn |
|
− pn |
|
) −cn |
|
( pn |
|
− pn |
|
|
|
|
|
|||
|
|
i, j |
|
i, j |
i, j |
−1) = 0.(1.16) |
|
|||||||||||
i, j+1/ 2 |
i, j+1 |
|
|
i, j−1/ 2 |
|
|
|
|
|
|
|
|||||||
Здесь cn i+1/ 2, j =2 cn i+1, j |
cn i, j /( cn i+1, j |
+ cn i, j ) . (1.17) |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ка |
|
|
В разностном уравнении (1.16) первое слагаемое предст вляет собой поток жидкости через левую границу ячейки (i, j), второе сл г емое – через правую границу, третье слагаемое – через верхнюю границу, четвертое
слагаемое – через нижнюю границу. При записи этого уравнения условно |
|||
|
|
|
т |
принято, что горизонтальные потоки направлены сл ва направо, а верти- |
|||
|
|
о |
|
кальные потоки направлены сверху вниз. Резуль а ыерасчета от указанной |
|||
условности направления потоков не зависят. |
|
|
|
|
л |
|
|
1.6.3. Построение разностных уравнений |
|||
для насыщенностии |
|
|
|
Найдем разность входящих и выходящих потоков воды через грани- |
|||
и |
|
|
|
цы ячейки и умножим на ∆t, это удет величина изменения объема воды |
|||
∆Vв в данной ячейке. Если раздел ть ееб |
на объем пор ячейки, то получим |
прирост водонасыщенности в ячейке за данный шаг времени. В результате
получим алгебраическое уравнение, аппроксимирующее дифференциаль- |
|||||||||||||||||||||||||||||||
ное уравнение (1.11): |
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
sn+1ij - |
sn ij = |
|
|
|
|
|
t |
|
(vn |
i+1/ 2, j |
F(sn |
i+1/ 2, j |
) - vn |
i−1/ 2, j |
F(sn |
i−1/ 2, j |
) + |
|||||||||||
|
|
|
|
|
x2 (mh)ij |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
+ vn |
i, j+1/ 2 |
F(sn |
i, j+1/ 2 |
) - vn |
i, j−1/ 2 |
F(sn |
i, j−1/ 2 |
)), |
|
(6.18) |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
н |
|
ая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
где |
v |
+ |
|
|
|
= |
cn i+1/ 2, j |
|
×( pn i+1, j - pn i, j ) . |
|
|
|||||||||||||||||||
|
|
|
i |
1 / 2, j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
о |
|
н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Значения F(s) определяются с учетом направления потока: |
|
|||||||||||||||||||||||||||||
|
|
F(sn iтр+1/ 2, j |
) = |
F(sn |
i+1, j |
) при |
|
p |
i+1, j |
³ p |
i, j |
. |
|
(1.19) |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
F(s |
n |
|
) |
|
при |
|
p |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
е |
к |
|
|
|
|
|
|
|
|
i+1, j |
£ p |
i, j |
|
|
|
|
|
|
||||||||||||
л |
|
|
|
|
|
|
i, j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Разностные уравнения (1.18) позволяют рассчитывать все значения |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
насыщенности водой для (n+1) –го временного слоя, если уже определе-
10