ИСТИНА

При численном моделировании течений в пористой среде необходимо определять теплофизические свойства флюидов, участвующих в фильтрации, в широком диапазоне давлений и температур. В некоторых случаях, например при моделировании течений в геотермальных системах или течений, возникающих при подземном захоронении углекислого газа, исследуемый диапазон термобарических условий содержит критические термодинамические параметры. Для устойчивости расчета гидродинамических процессов при критических условиях необходимо использовать неклассические независимые переменные моделирования, например переменные давление-энтальпия. Применение подобных переменных осложняется тем, что не существует быстрых методов определения теплофизических свойств веществ по заданному давлению и энтальпии, которые эффективно работают при до- и закритических условиях. В задачах фильтрации свойства флюидов часто определяются с помощью кубического уравнения состояния. Существуют алгоритмы расчета термодинамических равновесий в переменных давление-температура. Вычислительная эффективность данных алгоритмов основывается, в частности, на том, что для заданных давления и температуры термодинамические параметры фаз определяются явным образом. Алгоритмы имеют один итерационный процесс, необходимый только для определения многофазных термодинамических равновесий, а не для параметров фазы. Однако, переменные давление- температура и, следовательно, существующие алгоритмы для расчета свойств не могут использоваться в задачах численного моделирования фильтрации при околокритических условиях.