Приток из численного аквифера

Последнее сообщение
oilman007 38 9
Ноя 10

Под численным аквифером я подразумеваю ячейки на краю модели с увеличенным объемом.
Т.е. увеличивая объем в 1000 раз я считаю эта ячейка на краю модели у меня больше в 1000 раз и соответственно ее объем больше в 1000 раз.
Т.е. ее линейные размеры больше.
Но чтобы не рисовать ее такую большую я просто увеличиваю ее объем(MULTPV например), и, стало быть, исходя из формулы для расчета проводимости:

должен скорректировать TRANX, TRANY, TRANZ (с помощью MULTX, MULTY, MULTZ например).
Хочу для этого прогу написать даже...

Вопросы:
1) кто-нибудь этим занимался;
2) нужно ли это по вашему мнению.

Спасибо! smile.gif

RomanK. 2158 11
Ноя 10 #1

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

oilman007 38 9
Ноя 10 #2

Немного предыстории об возникновении данного вопроса:
была задача: смоделировать участок с гор. скв. И рассчитать дебит скв. при различных депрессиях. При этом пласт считать бесконечным.
Залежь нефтяная. Все ЧНЗ.
Есть все данные.
Давление равно давлению насыщения нефти газом.
Как я смоделировал бесконечный пласт: сделал огромный множитель на поровый объем (~10e10).
Так вот от того какой будет TRANX, TRANY, TRANZ в этих ячейках будет зависеть приток флюида в залежь->уровень пластового давления->кол-во выделившегося газа->дебит скв.

Dorzhi 894 12
Ноя 10 #3

а чем не устраивает численный аквифер того же эклипса?

oilman007 38 9
Ноя 10 #4

Dorzhi пишет:

а чем не устраивает численный аквифер того же эклипса?

поподробней можно?

dr-ulm 4 8
Ноя 10 #5

oilman007 пишет:

Немного предыстории об возникновении данного вопроса:
была задача: смоделировать участок с гор. скв. И рассчитать дебит скв. при различных депрессиях. При этом пласт считать бесконечным.
Залежь нефтяная. Все ЧНЗ.
Есть все данные.
Давление равно давлению насыщения нефти газом.
Как я смоделировал бесконечный пласт: сделал огромный множитель на поровый объем (~10e10).
Так вот от того какой будет TRANX, TRANY, TRANZ в этих ячейках будет зависеть приток флюида в залежь->уровень пластового давления->кол-во выделившегося газа->дебит скв.


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

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

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

MULTPV на линейные размеры не влияет.

oilman007 38 9
Ноя 10 #6

dr-ulm пишет:

MULTPV на линейные размеры не влияет.

а идея с изменением TRANXYZ Вам не нравится? (см. выше)

dr-ulm 4 8
Ноя 10 #7

oilman007 пишет:

а идея с изменением TRANXYZ Вам не нравится? (см. выше)


А что с TRANXYZ? Если последняя (крайняя) ячейка работает, то это уже не бесконечный пласт.

RomanK. 2158 11
Ноя 10 #8

dr-ulm пишет:

А что с TRANXYZ? Если последняя (крайняя) ячейка работает, то это уже не бесконечный пласт.


Да, создать действительно бесконечный пласт нелегко. По теме, приток из аквайфера можно представить в виде закона Ома: I=U/R, то что говорите вы это изменение сопротивления R, для заданного I. Также этого можно достичь просто изменив U.

oilman007 38 9
Ноя 10 #9

dr-ulm пишет:

А что с TRANXYZ? Если последняя (крайняя) ячейка работает, то это уже не бесконечный пласт.

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

Представим, что у меня модель в виде последовательности ячеек.
Т.е. по i=X, по j=1 и по k=1.
допустим что dx=dy
в ячейке (1,1,1) находится скважина. Работает и что то добывает.
Так вот если X будет достаточно большим то можно считать что пласт бесконечный... [1]
Чем дальше от скважины тем меньше мне нужна измельченная сетка,
т.е. я могу апскелить ячейки начиная с какого то номера,
или просто объединить их в одну большую...[2]
а здесь вместо одной большой могу вставить маленькую с увеличенным поровым объемом и измененным TRANX...[3]

Так? smile.gif

RomanK. 2158 11
Ноя 10 #10

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

oilman007 38 9
Ноя 10 #11

RomanK. пишет:

Вместо увеличения проницаемости можешь увеличить поровый объем.

исходя из формулы при увеличении PV нужно уменьшать PERM( или TRAN) и наоборот

RomanK. 2158 11
Ноя 10 #12

oilman007 пишет:

исходя из формулы при увеличении PV нужно уменьшать PERM( или TRAN) и наоборот


исходя из закона Ома. для заданного потока I=U/R, сопротивление обратно пропорционально проницаемости R~1/K. следовательно I~U*K, где U~PVOL. I~PVOL*K.

oilman007 38 9
Ноя 10 #13

RomanK. пишет:

исходя из закона Ома. для заданного потока I=U/R, сопротивление обратно пропорционально проницаемости R~1/K. следовательно I~U*K, где U~PVOL. I~PVOL*K.

исходя из Ваших рассуждений видно что мы говорим об одном и том же smile.gif
Для сохранения неизменного I увеличивая одно нужно уменьшать другое.

Складывается впечатление что прога действительно не нужна...

dr-ulm 4 8
Ноя 10 #14

oilman007 пишет:

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


Если сделать серию моделей вида:

две ячейки: ячейка со скважиной, краевая ячейка с MULTPV;
три ячейки: ячейка со скважиной, обычная ячейка, краевая ячейка с MULTPV;
четыре ячейки: ячейка со скважиной, 2 обычных ячейки, краевая ячейка с MULTPV
пять ячеек - и т.д.

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

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

Когда ты делаешь модель небольших линейных размеров, сравнимых с длиной горизонтального ствола, и на краю ставишь MULTPV, то это очень далеко от модели бесконечного пласта.

Это модель конечного пласта со смешанными краевыми условиями (промежуточными между неперетоком и постоянным давлением).

Кривая падения забойного давления в такой модели будет отличаться от аналитического решения неустановившегося режима (кроме начального участка). Это можно проверить на радиальной модели с вертикальной скважиной, например.

dr-ulm 4 8
Ноя 10 #15

oilman007 пишет:

Представим, что у меня модель в виде последовательности ячеек.
Т.е. по i=X, по j=1 и по k=1.
допустим что dx=dy
в ячейке (1,1,1) находится скважина. Работает и что то добывает.
Так вот если X будет достаточно большим то можно считать что пласт бесконечный... [1]


При этом, если хотим считать пласт бесконечным, то есть нас интересуют только результаты, которые попадают в неустановившийся режим, то вообще говоря без разницы, что там на краю. Главное, чтобы пласт был длинный (линейно). Этого можно добиться увеличением линейных размеров ячеек. Можно вместо этого постепенно увеличивать MULTPV, но разницы нет. С линейными размерами нагляднее, по-моему.

Гоша 1180 12
Ноя 10 #16

oilman007 пишет:

Для сохранения неизменного I увеличивая одно нужно уменьшать другое.

Складывается впечатление что прога действительно не нужна...


Не нужна, Роман тебе верно сказал - примерно и получается: во сколько раз наращиваем размер (MULTPV), во столько раз и уменьшаем проводимость (MULTX/Y)
Менять линейный размер ячеек можно, но вряд ли так уж "наглядно" - при визуализации они будут намного больше ячеек со скважинами.

Наращивать размер можно до тех пор, пока изменение давления в увеличенной ячейке за расчетное время не станет достаточно малым.

Go to top