Ноя 10
Под численным аквифером я подразумеваю ячейки на краю модели с увеличенным объемом.
Т.е. увеличивая объем в 1000 раз я считаю эта ячейка на краю модели у меня больше в 1000 раз и соответственно ее объем больше в 1000 раз.
Т.е. ее линейные размеры больше.
Но чтобы не рисовать ее такую большую я просто увеличиваю ее объем(MULTPV например), и, стало быть, исходя из формулы для расчета проводимости:
должен скорректировать TRANX, TRANY, TRANZ (с помощью MULTX, MULTY, MULTZ например).
Хочу для этого прогу написать даже...
Вопросы:
1) кто-нибудь этим занимался;
2) нужно ли это по вашему мнению.
Спасибо! 
Опубликовано
27 Ноя 2010
Активность
16
ответов
4426
просмотров
5
участников
0
Рейтинг
Контекст
Eclipse x221
Нет, не нужно. 1000 это выдуманная величина, следовательно находить точные значения TRANX не имеет практического смысла. А для практики пожалуйста, если даже и программу решил написать в одно действие, любой опыт положителен.
Немного предыстории об возникновении данного вопроса:
была задача: смоделировать участок с гор. скв. И рассчитать дебит скв. при различных депрессиях. При этом пласт считать бесконечным.
Залежь нефтяная. Все ЧНЗ.
Есть все данные.
Давление равно давлению насыщения нефти газом.
Как я смоделировал бесконечный пласт: сделал огромный множитель на поровый объем (~10e10).
Так вот от того какой будет TRANX, TRANY, TRANZ в этих ячейках будет зависеть приток флюида в залежь->уровень пластового давления->кол-во выделившегося газа->дебит скв.
а чем не устраивает численный аквифер того же эклипса?
Это ты смоделировал не бесконечный пласт, а, скорее, границу постоянного давления.
Лучше сделать ряды ячеек с линейными размерами, растущими в геометрической прогрессии.
Смысл в том, чтобы линейное расстояние до края модели было большим, иначе у тебя наступит установившийся или псевдоустановившийся режим.
MULTPV на линейные размеры не влияет.
А что с TRANXYZ? Если последняя (крайняя) ячейка работает, то это уже не бесконечный пласт.
Да, создать действительно бесконечный пласт нелегко. По теме, приток из аквайфера можно представить в виде закона Ома: I=U/R, то что говорите вы это изменение сопротивления R, для заданного I. Также этого можно достичь просто изменив U.
Представим, что у меня модель в виде последовательности ячеек.
Т.е. по i=X, по j=1 и по k=1.
допустим что dx=dy
в ячейке (1,1,1) находится скважина. Работает и что то добывает.
Так вот если X будет достаточно большим то можно считать что пласт бесконечный... [1]
Чем дальше от скважины тем меньше мне нужна измельченная сетка,
т.е. я могу апскелить ячейки начиная с какого то номера,
или просто объединить их в одну большую...[2]
а здесь вместо одной большой могу вставить маленькую с увеличенным поровым объемом и измененным TRANX...[3]
Так?
Так и делают.
Поровый объем и проницаемость (из нее и вычисляется проводимость TRANX) напрямую зависят друг от друга. Вместо увеличения проницаемости можешь увеличить поровый объем. Так как точную величину все равно никто не знает, практическая значимость разделения мощности акфайфера на объем и проницаемость не имеет. Также никто не знает точно положение и конфигурацию А, поэтому на глазок рисуют линии, квадраты и проч. и подбирают объем. Хочешь в этой выдуманной линии и проницаемость в тысячу раз поднять - ради бога, кто мешает.
исходя из закона Ома. для заданного потока I=U/R, сопротивление обратно пропорционально проницаемости R~1/K. следовательно I~U*K, где U~PVOL. I~PVOL*K.
Для сохранения неизменного I увеличивая одно нужно уменьшать другое.
Складывается впечатление что прога действительно не нужна...
Если сделать серию моделей вида:
две ячейки: ячейка со скважиной, краевая ячейка с MULTPV;
три ячейки: ячейка со скважиной, обычная ячейка, краевая ячейка с MULTPV;
четыре ячейки: ячейка со скважиной, 2 обычных ячейки, краевая ячейка с MULTPV
пять ячеек - и т.д.
то между результатами прогонов этих моделей будет весьма существенная разница.
Что такое модель бесконечного пласта? Это модель, в которой существенные изменения давления не достигают края модели за расчетное время. Т.е., в расчетном периоде должен иметь место неустановшийся режим. Этого можно достичь, например, путем увеличения линейных размеров ячеек от скважины наружу в геометрической прогрессии.
Когда ты делаешь модель небольших линейных размеров, сравнимых с длиной горизонтального ствола, и на краю ставишь MULTPV, то это очень далеко от модели бесконечного пласта.
Это модель конечного пласта со смешанными краевыми условиями (промежуточными между неперетоком и постоянным давлением).
Кривая падения забойного давления в такой модели будет отличаться от аналитического решения неустановившегося режима (кроме начального участка). Это можно проверить на радиальной модели с вертикальной скважиной, например.
При этом, если хотим считать пласт бесконечным, то есть нас интересуют только результаты, которые попадают в неустановившийся режим, то вообще говоря без разницы, что там на краю. Главное, чтобы пласт был длинный (линейно). Этого можно добиться увеличением линейных размеров ячеек. Можно вместо этого постепенно увеличивать MULTPV, но разницы нет. С линейными размерами нагляднее, по-моему.
Не нужна, Роман тебе верно сказал - примерно и получается: во сколько раз наращиваем размер (MULTPV), во столько раз и уменьшаем проводимость (MULTX/Y)
Менять линейный размер ячеек можно, но вряд ли так уж "наглядно" - при визуализации они будут намного больше ячеек со скважинами.
Наращивать размер можно до тех пор, пока изменение давления в увеличенной ячейке за расчетное время не станет достаточно малым.