Выбор численных методов решения систем обыкновенных дифференциальных уравнений для быстродействующей модели пиролиза углеводородного сырья

Abstract

Актуальность исследования обусловлена необходимостью увеличения производства легких олефинов. Применение систем усовершенствованного управления технологическим процессом и оптимизации в реальном времени позволяет повысить эффективность пиролизных производств, но требует быстродействующей математической модели процесса. Цель: выбор метода численного решения систем обыкновенных дифференциальных уравнений, обеспечивающего набольшее быстродействие при расчете реакционного змеевика печи пиролиза. Сокращение времени, затрачиваемого на расчет каждого сценария, позволит использовать предлагаемую модель для задач оптимизации технологического процесса в реальном времени. Объект: математическая модель пиролиза этана, методы численного решения систем обыкновенных дифференциальных уравнений. Методы: системный анализ, математическое моделирование. Для решения систем обыкновенных дифференциальных уравнений применялись различные численные методы в явном виде, отличающиеся способом выбора шага интегрирования. Результаты. Разработана и опробована стационарная модель пиролиза этановой фракции. На разработанной модели выполнено сравнение скорости расчета при использовании численных методов решения систем обыкновенных дифференциальных уравнений с фиксированным и с адаптивным шагом и показано, что использование адаптивного шага при интегрировании позволяет сократить время расчета более чем в 20 раз (с более чем 11 ч до 34 мин) при сохранении точности расчетов. Такие результаты связаны с различными скоростями реакций по профилю реакционного змеевика - в зонах высоких температур и высоких концентраций исходных веществ требуется сокращение шага интегрирования для получения требуемой точности, в то время как на участках, характеризующихся низкими скоростями реакций, допустимо увеличение шага и сокращение общего количества расчетных итераций.

Relevance. The need to increase production of light olefins. The use of advanced process control systems and Real- Time Optimization makes it possible to increase the efficiency of steam cracking plants, but requires a high-speed mathematical model of the process. Aim. To select a method for numerical solution of systems of ordinary differential equations, which provides the highest speed when calculating the reaction coil of a steam cracking furnace. Reducing the time spent on calculating each scenario will allow the proposed model to be used for real-time process optimization tasks. Object. Mathematical model of ethane steam cracking, numerical methods for ordinary differential equations systems solution. Methods. System analysis, mathematical modeling. To solve the ordinary differential equations systems, various explicit numerical methods were used, differing in approach to integration step determination. Results. The authors have developed and tested a steadystate model of ethane steam cracking. The developed model was used to compare the calculation time required for solving ordinary differential equations systems using different numerical methods. It was demonstrated, that the use of an adaptive integration step reduces calculation time by more than 20 times (from more than 11 hours to 34 minutes) while maintaining the accuracy of calculations. This is due to different reaction rates through the length of the reaction coil - in areas of high temperatures and high concentrations of reagents, a reduction in the integration step is required to obtain the desired accuracy. And in low reaction rates areas an increase in the step and reduction in the total calculated iterations are acceptable.