Здравствуйте
Моделирую залежь с двойной пористостью, в связи с чем возникли некоторые вопросы.
По порядку:
Всего в ремасштабированной геомодели 18 слоев.
Соответственно в dimens задаю количество слоев равное 36.
Куб NTG дублирую.
Далее загружаю кубы пористости и проницаемости.
Соотношение пористости матрицы и трещины составляет примерно 0,85/0,15, проницаемость примерно 1/8.
Все остальные необходимые кубы (насыщенность и т.д.) также дублируются.
Однако при попытке произвести пробный расчет, шаги сильно дробятся и в модели возникает отрицательное пластовое давление вплоть до -15000 атм.
Попытка задать поле давления явно, результатов не дала.
Однако, если задать кубы пористости и проницаемости матрицы и трещины одинаковыми, проблема исчезает.
Моделирую первый раз, мануалов никаких нету, помогите кто чем может.
Спасибо!
Что значит производите дублирование кубов? каким образом?
Появились новые данные
Я выяснил, что эта проблема связана с величиной сигма-фактора который характеризует связь между матрицой и трещиной.
Методом научного тыка я подобрал его величину так, чтобы модель нормально считалась, но в связи с этим
возникает другой вопрос: если сигма фактор определяется размером блока матрицы (не ячейки ГДМ а именно блока матрицы)
но как этот размер определить? Исследований никаких естественно нету. Может быть есть какие то усредненные, типичные значения?
Потому что не оч. хочется действовать наугад.
просто куб NTG матрицы идентичен кубу NTG трещин
Поскольку исследований никаких нет, то принимай "сигму" как параметр неопределенности при адаптации модели,
то есть подбирай так, чтобы не только хорошо сходилось решение, но и история воспроизводилась.
Тогда продолжим разговор.
Основной проблемой возникающей при разработке трещиноватых коллекторов является то, что при переводе скважин с фонтана на мех. добычу,
происходит резкий скачек обводненности, практически с 0 до 100%. Какие условия необходимо соблюсти в ГДМ и двойной пористостью/проницаемостью,
чтобы получить такой эффект? (У меня, в полученной модели в подавляющем большинстве скважин рост обводненности плавный)
Формула сигмы есть в мануале.
можно прикинуть ширину раскрытия трещины и частоту (густоту) трещин в породе.
Исходя из этого оценить сигму.
По поводу обводнения. Плавный рост обводнённости скорее всего связан с грубостью сетки.
Если у вас вертикальные ячейки около метра, то скачёк обводнённости, думаю, не получите.
Либо надо менять проводимость трещин и опять же подбирать сигму. 1/8 мне кажется маловато (или у вас проницаемость матрицы очень большая?).
SIGMA=4*{(1/Lx^2)+(1/LY^2)+(1/LZ^2)}
где Lx,y,z-размеры матрицы по осям X,Y,Z
PS: и кстати, посмотри фазухи по трещинам, какие они у тебя?
По поводу сигма-фактора, на мой взгляд, важно помнить, что это площадь взаимодействия между трещинами и матрицей на приходящаяся единицу объема коллектора. Я тоже думаю, что 1/8 маловато. Сейчас в моделях ставлю 0,5-2 - эта штука достаточно сильно влияет на выработку матрицы.
Основываясь на своём опыте моделирования DP, рекомендую посмотреть/подергать следующие опции, связанные с двойной пористостью (это сродни плясок с бубном, но тем не менее вроде помогает))):
1) NODPPM - если не поставить это слово, Эклипс будет домножать трещинную проницемость на трещинную пористость и использовать полученное значение проницаемости (может быть именно с этим связано то, что у тебя нет резкого скачка). Сам с этим накалывался.
2) попробуй переключить солвер на пре-97: 60я позиция в слове OPTIONS (ну и вообще другие позиции этого слова, которые касаются DP, например 48)
3) посмотри правильно ли ты прописал соединение скважин и аквиферов с сеткой - соединения должны быть к трещинной части
Удачи!
А какие они должны быть?
Мы, когда делали такую модель, сошлись, что они должны быть напоминать Х, т.е. тренды падения по нефти и роста по воде одинаковы с увеличением водонасыщенности.
А вообще, наверное, многое зависит от размера пор, какой именно там потом нефти с водой.