Май 24
Доброго времени суток. Пытаюсь смоделировать выпадение конденсата в CMG STARS. По непонятным мне причинам конденсат отказывается конденсироваться, преодолев давление конденсации. Композиционную можедь составил в WinProp, в 2-phase flash пишет около 20% жидких компонентов при давлении ниже давления конденсации. Переходя в STARS oil phase - 0.
Примечательно, что такую же модель, закинув в модуль GEM, конденсат выделился без проблем.
Вопрос знатокам в моделировании. Ключевая ошибка в моделировании конденсата допущена в STARS или же WinProp?
P.S. Такое ощущение, что в STARS есть дополнительная функция появления конденсата. Она вообще существует?
Опубликовано
25 мая 2024
Активность
5
ответов
738
просмотров
3
участника
0
Рейтинг
Это ты с K-values что то накосячил.
В том то и дело, что нет. К примеру NC6 всегда меньше нуля. Это говорит, что конкретно этот компонент должен быть в жидкой фазе изначально. К слову у метана k-value больше единицы, так как он не сжижаемый. Я правильно понимаю?
В смысле, меньше нуля? Отрицательный что ли? ))))))
Построй свои значения на лог-лог плоте K-value vs P. Когда компонент пересекает 1 - он в жидкуе фазу уходит, когда выше 1 - в газовой. Метан выше 1, праильно.
Ох, меньше единицы, конечно, переработал))
У меня кривые такие же по характеру, как и у вас.
В WinProp смотрю логи, там выделение на жидкую фазу есть
"Такое ощущение, что в STARS есть дополнительная функция появления конденсата. Она вообще существует?"
Нет, расчет фазового состояния осуществляется на основании заданных констант фазового равновесия (КФР, k-values).
"Ключевая ошибка в моделировании конденсата допущена в STARS или же WinProp?"
Возможно, была некорректно проведена конвертация композиционной модели в псевдокомпозиционную. Возможно, дело в неправильной инициализации модели (обратите внимание, мольные доли каких фаз Вы задаете для начального состояния). Наконец, попробуйте отобразить КФР всех компонентов в некоторой ячейке, глобальный компонентный состав и решить уравнение Речфорда-Райса вручную и убедиться, что получаете такие же значения мольных долей фаз, как и симулятор.
"К примеру NC6 всегда меньше единицы. Это говорит, что конкретно этот компонент должен быть в жидкой фазе изначально."
Нет, например, если взять Вашу PVT-модель и провести two-phase flash для 20000 кПа, 80 C и глобального компонентного состава (0.93, 0.04, 0.02, 0.01), то жидкой фазы не будет (мольная доля газовой фазы получится 1.0072, если включить опцию negative flash в WinProp), несмотря на то что КФР для трех из четырех компонентов меньше 1.