Расчет коэффициента сверхсжимаемости основных компонент природного газа методом молекулярной динамики

Abstract

Актуальность работы обусловлена необходимостью определения точного значения коэффициента сверхсжимаемости для подсчета запасов газа, прогнозирования технологических показателей разработки газовых и газоконденсатных месторождений, поведения углеводородной системы при движении в лифтовых трубах, а также подготовки скважинной продукции и ее транспортировки. Цель: разработка программного модуля расчета коэффициента сверхсжимаемости компонент природного газа для улучшения точности определения z-фактора в сравнении с классическим подходом, основанным на использовании уравнения состояния. Методы. Для определения z-фактора используется методика, основанная на использовании уравнения состояния Пенга–Робинсона, его модификация с шифт параметром, а также разработанный авторами программный модуль молекулярно-динамического моделирования, описывающий поведение системы на молекулярном уровне. В качестве модели межатомного взаимодействия частиц рассматривается потенциал Леннард-Джонса и NPT-ансамбль. Результаты. На основе предложенных методов моделирования показано, что метод молекулярной динамики позволяет с достаточной точностью рассчитывать коэффициент сверхсжимаемости для основных углеводородных компонентов природных газов. На основе уравнения состояния Пенга–Робинсона и его модификации с шифт параметром, а также с помощью моделирования методом молекулярной динамики для метана, этана, пропана, азота и двуокиси углерода определены зависимости коэффициента сверхсжимаемости от давления в диапазоне температур от 250 до 410 К. Произведено уточнение параметров потенциала Леннард–Джонса для этана, пропана, азота и двуокиси углерода, что позволило повысить точность при расчете коэффициента сверхсжимаемости. Показано, что средняя абсолютная ошибка метода молекулярной динамики практически во всех расчетах не превышает 3 %.

The relevance of the work is caused by the need to determine the exact value of the compressibility factor for calculating gas reserves, forecasting technological indicators for development of gas and gas condensate fields, the behavior of hydrocarbon systems flow in wells, and preparation of well production and its transportation. The aim of the research is to develop a program for computation of the compressibility factor for the main components of natural gas to achieve more accuracy in determining z-factor in comparison with classical approach, based on the equations of state. Methods. To determine the z-factor the authors have used the technique based on the Peng-Robinson equation of state, its modification with a shift parameter, and a software module for molecular dynamic simulation developed by the authors that describes the behavior of the system at the molecular level. The potential of Lennard-Jones and NPT-ensemble is considered as a model of interatomic interaction of particles. Results. Based on the proposed methods of modeling, it was shown that the molecular dynamics simulations make it possible to accurately calculate the compressibility factor for the main components of natural gases. Based of the Peng-Robinson equation of state and its modification with a shift parameter, as well as using the molecular dynamics simulation for methane, ethane, propane, nitrogen and carbon dioxide, the dependencies of the compressibility factor on pressure were determined in the temperature range from 250 to 410 K. Lennard-Jones potential parameters for ethane, propane, nitrogen and carbon dioxide were specified that made it possible to improve the accuracy in calculating the compressibility factor. It is shown that the average absolute error of the molecular dynamics method in almost all calculations does not exceed 3 %.