Free Resources

We know the secret of your success

Modelling and Simulation of the Fluid Catalytic Cracking Unit’s Riser-Reactor using a twelve-lump kinetic model

Introduction

Fluidized-Bed Catalytic Cracking (FCC) is needed to transform low-value heavy oil into more useful lighter chained products (LPG, fuel gas, gasoline etc.).  The FCC unit will give an approximate conversion of 75% of the feedstock into gasoline and lighter products.  The catalytic process has been shown to be more efficient than the thermal cracking process.  The FCC unit is the heart of the modern refinery (Paul et al., 2015). 

Although the FCC unit is the most effective unit in the refinery, its feed comes from other units.  The distillation unit separates the raw crude into several intermediates, namely, gasoline, kerosene, naphtha, diesel and gas oil (Sadeghbeigi, 2012).  The heaviest part of the crude, which cannot be refined in the barometrical tower, is warmed and sent to the vacuum tower, where separation into gas oil and residue occurs.  The residue is then sent units such as the de-asphalting unit, the residue cracker, the delayed coker for further processing.  The residue can also be used as asphalt for tarring roads.  The gas oil coming primarily from the atmospheric column, the vacuum tower, and the delayed coker is known to be the feed of the FCC unit (Lababidi and Al Humaidan, 2011).

A broad range of products is obtained from the FCC unit, namely, Liquefied Petroleum Gas (LPG), high-octane gasoline, fuel gases, etc. (Lababidi and Al Humaidan, 2011).

Due to its high sensitivity to various boundary conditions and several continuous disturbances from different sources, the complexity of the operation process of the FCC unit is very high.  The FCC unit is responsible for the production of almost 50% of all used gasoline in the world (Negrao and Baldessar, 2006).

The FCC unit consists of three main sub-units, namely;

  • The riser reactor
  • The stripping section/ dis-engager
  • The regenerator reactor

The riser-reactor is a key sub-unit in the Fluidized-bed Catalytic Cracking Unit (FCCU), as this is where the catalytic cracking of the gasoline occurs. The liquid feed enters the riser reactor through a nozzle at the bottom of the reactor to ensure optimum atomization. The feed droplets subsequently contact the hot catalyst particles and, as a result of heat transfer, are vaporized. The cracking reaction begins instantaneously. The resulting vapor-particles mixture moves upwards due to a pressure differential along the reactor length. As the motion occurs, a cracking reaction also occurs simultaneously. As a result of cracking, the velocity of the vapors increases along with the riser height. The reactions along the reactor length are known to be endothermic (Heydari et al., 2010).  The rate at which the feed vaporizes at the entry zone is a key performance indicator of the catalytic cracking reaction. At a point in the riser-reactor, the cracking reaction halts due to the deposition of coke on the catalyst surface, causing catalyst deactivation and a short contact time between the catalyst and the vaporized feed. The spent catalyst moves from the riser reactor to the top of the regenerator.  It is termed ‘deactivated’ as a result of the accumulation of coke and hydrocarbons on its surface. The accumulated hydrocarbons on the surface of the spent catalyst are separated from the catalyst before leaving the reactor by stripping with steam (Du et al., 2014). The hydrocarbon vapours generated are then sent to a fractionating column. In the regenerator, coke is burnt off from the catalyst with the passage of air for catalyst reactivation.  The flue gas leaving the regenerator contains a large quantity of carbon (II) oxide, which is burned to carbon (IV) oxide in a CO furnace or alternatively with a waste heat boiler so that available heat can be recovered. The catalyst leaves the regenerator through the bottom to be mixed with the feed stream. To reactivate the catalyst, the regenerator oxidizes the coke on the spent catalyst to form CO, CO2, and H2O. Compressed air enters the regenerator from the bottom of the regenerator through a grid distribution pattern.  The regenerated catalyst leaves the regenerator and is mixed with the feedstock at the riser reactor's base, ensuring a continuous cycle.  Fresh makeup catalyst could be added to maintain continuity and for provisions of withdrawal of used up catalyst over time (Ahmed et al., 2013). 

Modern FCC catalysts usually comprise four components: matrix, binder, filler, and the active component (Zeolite). Zeolite ranges from about 15-50 wt% of the catalyst. Industrially, zeolite-based FCC catalysts are usually termed faujasite or type Y (Kumar and Reddy, 2011). The catalytic sites in the zeolite are known to be strong acids providing most of the catalytic activity in the catalyst.  The nitrogen present usually restrains the catalytic activity of the catalyst in trace quantities in the feedstock; the Nitrogen content of a typical Vacuum Gas Oil (VGO) is in the range of 1000-2000 ppm (Dasila et al., 2012). Higher Nitrogen content would result in a significant reduction in conversion rate.  Amorphous alumina constitutes the matrix component of the FCC catalyst, enabling cracking of higher-boiling, larger feedstock molecules.  Mostly, silica sol is used as the binder and kaolin, the filler, providing the physical strength of the catalyst.  Contaminant metals, e.g., nickel, vanadium, iron, copper, etc., present in the feedstock also contribute to the catalyst deactivation (Nace et al., 1971; Paraskos et al., 1976; Shah et al., 1977; Weekman, 1968).

An essential advantage of the fluid catalytic cracking process is the fluidity of the catalyst between the reactor and regenerator, especially when fluidized with appropriate vapor phases. The fluidization of the catalyst ensures intimate interaction between the catalyst and hydrocarbons for better and more efficient cracking reactions (Amhammed, 2013). The feedstock quality, operating conditions of the riser reactor-regenerator sections, and the type of catalyst are factors that strongly determine the conversion and yield pattern of the FCC unit.  The interactions between the riser and regenerator reactors, the uncertainty in the cracking reactions, and the coke deposition and burning kinetics influences the complexity of the FCC unit (Dagde and Puyate, 2012).

Problem Statement

Several works have shown proposed models for the dynamic situation in the FCCU but have failed to imitate the industrial situation in reality for several reasons successfully. These reasons include the complexity of the reactions in the reactor and the inconsistency between the proposed models and the industrial situation (Ahmed et al., 2014).  This inconsistency can be attributed to the assumption of isothermal conditions in the reactor, the lumped model characteristics, the mode in which heat and mass transfer resistances are accounted for and in some cases not accounted for, the catalyst deactivation model, and the model accounting for the adsorption effects of aromatics (Bollas et al., 2007).

This work seeks to account for some of these ignored discrepancies. It imitates the steady-state situation in the FCCU successfully to analyze the optimal length of the riser reactor using the operating conditions of the Fawley refinery’s FCCU. (Berry et al., 2004; Das et al., 2003; Gupta et al., 2007; Gupta and Kumar, 2008).

Aim and Objectives

This study aims to model and simulate the riser-reactor of Fawley Refinery’s Fluid Catalytic Cracking Unit (FCCU) using a modified twelve-lump kinetics model. It has the following objectives;

  • Objective 1: Development of a feasible and comprehensive model for the riser reactor of Fawley’s FCCU, accounting for the mass and heat transfer resistances, catalyst deactivation function, the poisoning effect of basic Nitrogen, and the adsorption effect of aromatics, asphaltenes, and resins.
  • Objective 2:  Development of a detailed hydrodynamic model for the riser reactor of Fawley’s FCCU.
  • Objective 3: Obtaining numerical solutions of the developed models using MATLAB built-in ode23 function, needed to predict the respective yields of vacuum gas oil, gasoline, liquefied petroleum gas, fuel gas, and coke.
  • Objective 4:  Comparison of the simulated results with plant data gotten from Fawley refinery to reveal the accuracy and efficiency of the developed models to imitate the real-time situation in the FCCU.  This comparison would be done with respect to the yield of vacuum gas oil, the cracked products, the fluids pressure drop, solid-phase temperature, and other critical process parameters.
  • Objective 5:  Carrying out sensitivity analysis on the developed models.

Research Questions

The identified research questions for this project are provided below:

  • What energy models are incorporated into the developed model to take cognizance of the heat transfer resistance, instant vaporization of feed, and one-dimensional vaporization of feed?
  • What hydrodynamic model is incorporated into the developed model to account for the pressure drop along the reactor length?
  • Are models incorporated to account for the catalyst deactivation, poisoning effect of basic Nitrogen, and adsorption effects of aromatics?
  • How would the operating conditions and parameters of Fawley Refinery’s FCCU be gotten?
  • Are the representation of the viscosities and the specific heat capacities of vacuum gas oil and other lumped components a function of temperature in the developed model, since in reality, there is a temperature variation along the reactor length as a result of heat exchange? (Leon-Becerril et al. 2004; Sharma, 2011).

Deliverables

The deliverables of this project are a project report, mathematical models, and simulated results. The models should consistently mimic the real-time situation in the riser-reactor of Fawley refinery’s FCCU. Also, the report should contain complete documentation of how the mathematical models were arrived at and tables of the simulated results along the reactor length.

Relevance

This project mainly focuses on mimicking the real-time situation in the riser-reactor of Fawley Refinery’s FCCU.

Methodology

This project focuses on secondary research, development of mathematical models, and model simulation, and they are discussed below:

Secondary research

The secondary research in this project will utilize a systematic approach (Johnson et al., 2016) to review the works of literature. The steps involved in the systematic review of the literature are provided below:

  • Step 1: Identify the research questions that can be used for the project.
  • Step 2: Identify the keywords that should be used to research the works of literature.
  • Step 3: Extract the journals and books that are appropriate for this project.
  • Step 4: Write the literature review chapter.

Model development

The development of mathematical models for the riser-reactor situation in Fawley refinery’s FCCU and simulating the developed models. The development of the mathematical models are in stages:

  • Stage 1:  Model development for the vaporizer subsystem
  • Stage 2: Using a preferred and appropriate kinetics of catalytic cracking in the FCCU’s riser reactor.
  • Stage 3:  Developing an appropriate catalyst deactivation model (Dagde and Puyate, 2012).
  • Stage 4: Developing an appropriate hydrodynamic model along the reactor length (American Petroleum Institute, 1992).
  • Stage 5: Taking an appropriate mass balance around the riser reactor of Fawley refinery’s FCCU (Geankoplis, 2011).
  • Stage 6: Taking an appropriate energy balance around the riser-reactor of Fawley refinery’s FCCU.
  • Stage 7: Carrying out an appropriate force balance around the riser-reactor of Fawley refinery’s FCCU.

Model simulation

The developed models would be transformed into appropriate codes on the MATLAB platform, and the codes would also be run and simulated on the MATLAB platform.

The codes and simulated results would be documented both as MATLAB and Excel spreadsheets.

Evaluation

The risk assessment conducted for this project is provided in the table below:

Table 1:  Risk assessment

Risk

Impact

Mitigation Plan

Inability to meet the deadline

Low

Get an extension from the supervisor in due time

Inability to get plant data

High

Refer to journals and institutes to extrapolate plant data

Insufficient knowledge in developing and simulating mathematical models

High

Refer to journals, textbooks, online forums, and other capable colleagues for help.

Schedule

Table 2: Project Plan

Task Name

Start Date

End Date

Duration (Days)

Initial Research

23/09/2021

07/10/2021

14

Proposal

07/10/2021

28/10/2021

21

Secondary Research

28/10/2021

07/12/2021

40

Introduction Chapter

07/12/2021

12/12/2021

5

Literature Review Chapter

12/12/2021

05/01/2022

24

Methodology Chapter

05/01/2022

17/01/2022

12

Development of the Mathematical Models

17/01/2022

15/03/2022

60

Presentation 1

15/03/2022

23/03/2022

8

Simulation of the Mathematical Models

23/03/2022

06/04/2022

14

Evaluation of Simulated Results

06/04/2022

13/04/2022

7

Discussion Chapter

13/04/2022

23/04/2022

10

Evaluation Chapter

23/04/2022

28/04/2022

5

Conclusion Chapter

28/04/2022

30/04/2022

2

Project Management Chapter

30/04/2022

01/05/2022

2

Abstract and Report compilation

01/05/2022

03/05/2022

2

Report Proofreading

03/05/2022

13/05/2022

10

Presentation 2

13/05/2022

23/05/2022

10

References

Ahmed, A., Maulud, A., Ramasamy, M., Lau, K. K., and Mahadzir, S. (2014). 3D CFD Modelling and Simulation of RFCC Riser Hydrodynamics and Kinetics. J. Appl. Sc., Vol. 14, No. 23, 3172-3181.

Ahmed, H. S., Shaban, S. A., Menoufy, M. F., and El Kady F. Y. (2013). Effect of Catalyst Deactivation on Vacuum Residue Hydrocracking. Egypt. J. Pet.,Vol. 22., No. 3, 367-372.

American Petroleum Institute (API). (1992). Technical Data Book – Petroleum Refining (5th ed.). New York.

Amhammed, M. (2013). Applications of Lumping Kinetics Methodology to Complex Reactive Mixtures. Heriot-Watt University, Department of Chemical Engineering. Heriot-Watt University.

Berry, T. A., McKeen, T. R., Pugsley, T. S., and Dalai, A. K. (2004). Two- Dimensional Reaction Engineering Model of the Riser Section of a Fluid Catalytic Cracking Unit. Ind . Eng. Chem. Res., Vol. 43, 5571-5581.

Bollas, G. M., Lappasa, A. A., Iatridisa, D. K., and Vasalos, I. A. (2007). Five-Lump Kinetic Model with Selective Catalyst Deactivation for the Prediction of the Product Selectivity in the Fluid Catalytic Cracking Process. Catal. Today, Vol. 127, 31-43. doi:10.1016/j.cattod.2007.02.037

Dagde, K. K., and Puyate, Y. T. (2012). Modelling And Simulation Of Industrial FCC Unit: Analysis Based on Five-lump Kinetic Scheme for Gas-Oil Cracking. IJERA, Vol. 2 No. 5, 698-714.

Das, A. K., Baudrez, E., Martin, G. B., and Heynderickx, G. J. (2003). Three-Dimensional Simulation of a Fluid Catalytic Cracking Riser Reactor. Ind. Eng. Chem. Res., Vol. 42, 2602-2617.

Dasila, P. K., Choudhury I., Saraf D., Chopra S., and Dalai A. (2012). Parametric Sensitivity Studies in a Commercial FCC Unit. Advances in  Chemical Engineering and Science, Vol. 2, 136-149.

Du, Y. P., Yang, Q., Zhao, H., and Yang, C. H. (2014). An Integrated Methodology for the Modeling of Fluid Catalytic Cracking (FCC) Riser Reactor. Applied Petrochemical Research, Vol. 4, 423-433.

Geankoplis, C. J. (2011). Transport Processes and Separation Process Principles. New Jersey: Pearson Education Inc.

Gupta, R. K., and Kumar, V. (2008). Fluid Catalytic Cracking Riser Modeling in Heat Transfer Mode. CPPM, Vol. 3 No. 1, 1934-2659. doi:10.2202/1934-2659.1123.

Gupta, R. K., Kumar, V., and Srivastava, V. K. (2007). A New Generic Approach for the Modelling of Fluid Catalytic Cracking (FCC) Riser Reactor. Chem. Eng. Sci., Vol. 62 No. 17, 4510–4528.

Heydari, M., AleEbrahim, H., and Dabir, B. (2010). Study of Seven-lump Kinetic Model in the Fluid Catalytic Cracking Unit. American Journal of  Applied Sciences, Vol. 7 No. 1, 71-76.

Johnson, D., Deterding, S., Kuhn, K.A., Staneva, A., Stoyanov, S. and Hides, L., 2016. Gamification for health and wellbeing: A systematic review of the literature. Internet interventions, 6, pp.89-106.

Kumar, V., and Reddy, A. S. (2011). Why FCC Riser is Taller than Model Predictions. American Institute of Chemical Engineers Journal, Vol. 57 No. 10, 2917-2920. doi:10.1002/aic.12499

Lababidi, H., & Al Humaidan, F. (2011). Modeling the Hydrocracking Kinetics of Atmospheric Residue in Hydrotreating Processes by the Continuous Lumping Approach. Energy Fuel, 25(5), 1939-1949

Leon-Becerril, E., Maya-Yescas, R., and Salazar-Sotelo, D. (2004). Effect of Modelling Pressure Gradient in the Simulation of Industrial FCC Risers. Chem. Eng. J., Vol. 100, 181-186.

Nace, D. M., Voltz, S. E., and Weekman (Jr.), V. W. (1971). Application of a Kinetic Model for Catalytic Cracking. Effects of Charge Stocks. Ind. Eng. Chem., Process Des. Dev., Vol. 10, 530-537.

Negrão, C. O., and Baldessar, F. (2006). Simulation of Fluid Catalytic Cracking Risers – A Six Lump Model. The 11th Brazillian Congress of Thermal Sciences and Engineering. Curitiba, Brazil.

Paraskos, J. A., Shah, Y. T., Mckinney, J. D., and Carr, N. L. (1976). A Kinetic Model for Catalytic Cracking in a Transfer-Line Reactor. Ind. Eng. Chem. Process Des. Dev., Vol. 15, 165-169.

Paul, D., Chaudhari, P. K., and Thakur, R. S. (2015). Simulation of FCC Riser Reactor Based on Ten Lump Model. International Journal of  Engineering Research and Applications, Vol. 5 No. 7, 59-67.

Paul, D., Thakur, R. S., and Chaudhari, P. K. (2015). Simulation of FCC Riser Reactor using Five Lump Model. International Journal of  ChemTech Research, Vol. 8 No 6., 750-758.

Sadeghbeigi, R. (2012). Fluid Catalytic Cracking Handbook (3rd ed.). Oxford, UK: Elsevier Inc.

Shah, Y. T., Huling, G. T., Paraskos, J. A., and Mckinney, J. D. (1977). A Kinematic Model for an Adiabatic Transfer-Line Catalytic Cracking Reactor. Ind. Eng. Chem. Process Des. Dev., Vol. 16, 8994.

Sharma, H. (2011). Modeling and Simulation of Fluid Catalytic Cracking Riser. Patalia, India: Department of Chemical Engineering, Thapar University.

Weekman, V. (1968). A Model of Catalytic Cracking Conversion in Fixed, Moving and Fluid-Bed reactors. Ind. Eng. Chem. Proc. Des., Vol. 7, 90-95. doi:10.1021/i260025a018 

Last updated: Dec 02, 2021 12:51 PM

Do you need help with your academic work ? Let us help

AcademicianHelp

Your one-stop website for academic resources, tutoring, writing, editing, study abroad application, cv writing & proofreading needs.

Get Quote
TOP