COMPUTER SIMULATION OF HYDRODEMETALLATION IN FIXED-BED REACTORS

Abstract

An isothermal model for hydrodemetallation (HDM) of crude oils in catalytic fixed-bed reactors is proposed. This model involves a consecutive reaction mechanism, which is capable of accounting for particle deposit profiles with interior maxima. Consistent with the fact that HDM catalysts are conglomerates formed by precipitation, the porous catalyst itself is modeled as randomly overlapping spheres of equal size. The metal is deposited as growing metal sulfide crystallites on the inner surface of the catalyst. These crystallites originate from a certain number of randomly scattered nuclei and increase in size as the deposition proceeds. The random sphere model for the catalyst and the deposit provides the changes in the catalyst pore structure—local porosity and surface area.

The mass transport within the domain of the particle is due to restricted liquid diffusion, since the diameter of the metal bearing compound (porphyrin) and the intermediate are comparable to the pore size. The diffusion restrictions taken into account are the enhanced drag imposed on a molecule by adjacent pore walls and steric partitioning.

Since the deposition process is much slower than diffusion and reaction, the pseudo-steady-state assumption can be justified. The equations of conservation for mass are solved by orthogonal collocation on finite elements. Based on this solution technique a computer simulation program of HDM is designed that allows two modes of operation: constant temperature and constant conversion. The simulation program “SIMULA” is highly flexible with regard to reaction kinetics, catalyst structure, reactor design, and operating conditions. In comparison to a base case with uniform activity, the effect of intraparticle (radial) and bed (axial) activity profiles on the conversion rate is discussed. For the case investigated, a radial distribution of activity higher at the center of the particle than at the edge can increase catalyst life by 25%, but axial distribution was less successful.

KEYWORDS:

• Hydrodemetallation
• Catalytic reactor
• Computer simulation
• Optimal catalyst distribution
• Random structure

Notes

Present Address: DEGUSSA AG, D-6450 Hanau, Fed. Rep. of Germany.

To whom correspondence should be addressed.