IMPROVED CATALYTIC REFORMING

CROSS-REFERENCE TO RELATED APPLICATIONS [0001] Not applicable.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT [0002] Not applicable.

FIELD OF THE INVENTION

[ΘΘ03] The present disclosure relates generally to arornatization of hydrocarbons with a reforming catalyst. Particularly, the present disclosure relates to a catalytic reforming process in terms of improving the benzene production over the run-length of the reforming catalyst.

BACKGROUND OF THE INVENTION

[0004] The catalytic conversion of hydrocarbons into aromatic compounds, referred to as arornatization or reforming, is an important industrial process. The arornatization reactions may include dehydrogenation, isomerization, and hydrocraeking the hydrocarbons, each of which produces specific aromatic compounds. These reactions are generally conducted in one or more arornatization reactors containing an arornatization catalyst. The catalyst may increase the reaction rates, production of desired aromatics, and/or the throughput rates for the desired aromatic compounds.

[0005] The economics of a chemical process are affected by many factors, including capital cost, plant efficiency, and margin of the product. In the catalytic reforming process, the cost of the catalyst, the run- length of the catalyst and the throughput of the plant play an important role in determining the economics. Given their commercial importance, an ongoing need exists for improved methods of catalytic reforming processes. SUMMARY OF THE INVENTION

[0006] Disclosed herein is a method of catalytic reforming comprising feeding a feedstream comprising C ₆ -convertibles to one or more reactors; contacting the feedstream with a reforming catalyst; selecting values for a LHSV, a H2/HC ratio, and a conversion of Ce-cotivertibies from a deactivation kinetic model so as to maximize a net present amount of benzene produced over a run-length of the reforming catalyst; operating the one or more reactors at the selected LHSV, the selected H2/HC ratio, and the selected conversion of Ce-convertibles; andrecovering an effluent from the reactor, wherein the effluent comprises at least about 40 wt % benzene.

[0007] Also disclosed herein is a method of catalytic reforming comprising kinetically modeling a deactivation of a reforming catalyst with a deactivation kinetic model to determine a mn-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of Ce-convertibles in the feedstream; calculating a net present amount of benzene that is produced over the mn- length of the reforming catalyst; and selecting values for the LHSV, the H2/HC ratio, and the conversion of ( ^" ^-convertibles that maximize the net present amount of benzene produced over the run-length of the reforming catalyst,

[0008] Also disclosed herein is a catalytic reactor comprising a feedstream and a reforming catalyst wherein the reactor is configured to produce a net present amount of benzene over a mn-length of the reforming catalyst; wherein the run-length is determined as a function of a LHSV of the feedstream, a H2/HC ratio of the feedstream, and a conversion of Ce-convertibles in the feedstream by kinetically modeling a deactivation of the reforming catalyst with a deactivation kinetic model; and wherein values for the LHSV of the feedstream, the H2/HC ratio of the feedstream, and the conversion of CVconvertibles in the feedstream are selected to maximize the net present amount of benzene produced over the run-length of the reforming catalyst.

[0009] Also disclosed herein is a computer system comprising one or more processors; a non-transitory memory coupled to the processor, wherein the non- transitory memory comprises executable instructions that, when executed, cause the one or more processors to: kinetically model a deactivation of a reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of Cg-convertibles in the feedstream; calculate a net present amount of benzene that is produced over a mn-length of the reforming catalyst; and select values for the LHSV, j

the H2/HC ratio, and the conversion of Ce-convertibles that maximize the net present amount of benzene prod u ced over a run-length of the reforming catalyst.

[0010] Also disclosed herein is a system, which may optionally be utilized in combination with any of embodiments 1 to 77, comprising a catalytic reactor comprising a feedstream and a reforming catalyst; and a computer system coupled to the catalytic reactor and having a deactivation kinetic model to: kinetically model a deactivation of the reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of Ce-convertibles in the feedstream; and calculate a net present amount of benzene that is produced over the run-length of the reforming catalyst; wherein the catalytic reactor operates at values for the LHSV, the H2/HC ratio, and the conversion of Ce-convertibles that maximize the net present amount of benzene produced over a run-l ength of the reforming catalyst.

[0011] Also disclosed herein is a non-transitory computer readable medium, which may be optionally be utilized in combination with any of embodiments 1 to 92, comprising a model unit to: kinetically model a deactivation of a reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of Ce-convertibles in the feedstream; and calculate a net present amount of benzene that is produced over a run-length of the reforming catalyst.

BRIEF DESCRIPTION OF THE DRAWINGS

[0012] For a more detailed description of the embodiments of the disclosure, reference will now be made to the accompanying drawings, wherein:

[0013] Figure 1 is a schematic of one embodiment of a catalytic reactor system, the operation of which may be improved in accordance with the di sclosure.

[0014] Figure 2 illustrates a computer system suitable for implementing one or more embodiments of the disclosure.

[0015] Figure 3 illustrates a computer system and catalytic reactor system suitable for implementing one or more embodiments of the disclosure.

[0016] Figure 4 is a plot of fouling curves predicted by a deactivation kinetic model and an analytical solution of Equation (12) fit to each of the fouling curves.

[0017] Figure 5 is a plot of Ad in units of hr ^~! as a function of H2/HC ratio and conversion of Ce-convertibles for an isothermal reactor.

[0018] Figure 6 is a plot showing k _d (T ₀ ) as a function of H2/HC ratio for an isothermal reactor, as given by Equations (9) and ( 13) for the average over the reactor.

[0019] Figure 7 is a plot of fouling curves predicted by a deactivation kinetic model and analytical solutions fit to the fouling curves under various process conditions for a series of adiabatic reactors.

[0020] Figure 8 is a plot of the linear correlation of t _max with t _max (iso) when the latter is obtained from Equation (18) under the same process conditions and when t _max is determined by fitting Equation (17) to fouling curves predicted from a deactivation kinetic model for a series of adiabatic reactors.

[0021] Figure 9 is a plot of T _a as a function of φ obtained when a second order polynomial represented by Equation (22), i.e., second order with respect to φ, is fitted to values of T _a determined by fitting Equation (17) to the fouling curves predicted by a deactivation kinetic model.

[0022] Figure 10 is a plot of t _max as a function of t _!nax (iso) found by repeating full calculations with a deactivation kinetic model using the correlation in Equation (23) and fitting Equation (17) to the results.

[0023] Figure 1 1 is a plot of T _a as a function of φ determined by repeating full calculations with a deactivation kinetic model using the correlation in Equation (23) and fitting Equation (17) to the results. „

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0024] Disclosed herein are methods, systems, and reactors for improving catalytic reforming processes. "Aromatization", "aromatizing", and "reforming" as used herein refer to treatment of a hydrocarbon-containing feedstream to provide an arornatics enriched product (i.e., a product whose arornatics content is greater than in the feedstream). Catalytic reforming may comprise the production of reforoiate in a catalytic reformer (e.g., a reactor), the reform ate having a high aromatic content and very low olefins content (alkenes). Various catalytic reforming processes and catalysts (e.g., silica-bound L-type zeolite supporting platinum, chloride, and fluoride, and alumina supported platinum, tin, and chloride) are described herein. The methods as disclosed herein improve the catalytic reforming process by maximizing the production of benzene per unit of catalyst. The broader economics can be evaluated by using this knowledge and factoring in product margin, unit cost of the catalyst and other commercial vari ables such as the economic value of the hydrogen produced.

[0025] The catalytic reforming processes which may be improved using the disclosed inventive subject matter are described in more detail below, including both isothermal and adiabatic reactor designs. In general, the disclosed inventive subject matter may comprise kineticaily modeling the deactivation of the reforming catalyst with a deactivation kinetic model. The run-length of the reforming catalyst as a function of the percentage of C ₆ hydrocarbons that are convertible under the employed catalytic reforming process conditions (herein referred to as "conversion of C ₆ - convertibles") in the feedstream, the hydrogen-to-hydrocarbon mole ratio (hereinafter referred to as "H2/HC ratio") of the feedstream, and the liquid hourly space velocity (hereinafter referred to as "LHSV") of the feedstream may be determined from the deactivation kinetic model. The amount of benzene produced per unit of catalyst during the life cycle (e.g., run-length) of the reforming catalyst may then be determined from the run-length and the benzene yield per month. A catalytic reforming process may be improved by selecting values for LHSV, H2/HC ratio, and conversion of C ₆ - convertibies that maximize a net present amount of benzene produced per unit of catalyst over the run-length of the catalyst,

[0026] The catalytic reforming processes which may be improved according to the inventive subject matter may comprise one or more reactors. A feedstream comprising Ce-convertibles may be fed to the one or more reactors. Within the one or more reactors, the feedstream may be contacted with a reforming catalyst under reaction conditions (e.g., temperature, flow, and pressure) and improved conditions as disclosed herein for reforming of one or more components in the feedstream. Particularly, the reactor may be operated at the selected LHSV, the selected H2/HC ratio, and the selected conversion of Ce-convertibles determined as disclosed herein. The effluent of the reactor may comprise at least about 40 wt % benzene by weight of the effluent.

[0027] "Reactor", "catalytic reactor", and "aromatization reactor" as used herein refer interchangeably to the reactor vessel, reactor internals, and associated processing equipment, including but not limited to the catalyst, inert packing materials, scallops, flow distributors, center pipes, reactor ports, catalyst transfer and distribution system, furnaces and other heating devices, heat transfer equipment, and piping. Any suitable reactor configuration and associated processing equipment may be used to carry out the improved catalytic reforming process as disclosed herein. In various embodiments, the catalytic reforming process may comprise a fixed bed system, a moving bed system, a fluidized bed system, or combinations thereof. The catalytic reforming process may be batch or continuous. In an embodiment, the catalytic reforming process may comprise a fixed bed system comprising one or more fixed bed aromatization reactors in series. In a fixed bed system, the feed may be preheated in furnace tubes and passed into at least one reactor that contains a fixed bed of the catalyst. The flow of the feed to such a system can be upward, downward, or radial through the reactors.

[0028] One or more embodiments of the inventive subject matter may be used for or with a catalytic reactor system comprising at least one catalytic reactor and its corresponding processing equipment. Alternatively, one or more embodiments of the inventive subject matter may be used for or with a catalytic reactor system comprising a plurality of catalytic reactors in series.

[0029] Figure 1 illustrates one embodiment of a catalytic reactor system 100 to be improved via the disclosed inventive subject matter. In the embodiment shown in Figure I , the catalytic reactor system 100 may comprise four catalytic reforming reactors in series: reactors 10, 20, 30, and 40. However, the catalytic reactor system 100 may comprise any suitable number and configuration of aromatization reactors, for example one, two, three, five, six, or more reactors in series and/or in parallel. As aromatization reactions are highly endothermic, large temperature drops occur across the reactors. Therefore, each reactor 10, 20, 30, and 40 in the series may comprise a corresponding furnace 11, 21 , 31, and 41, respectively, for preheating components to a desired temperature for maintaining a desired reaction rate. Alternatively, one or more reactors may share a common furnace where practical. In embodiments, oxygenate and/or nitrogenate may be added within reaction zone 120 to increase/maintain catalyst activity, as further discussed hereinbelow.

[0030} In Figure 1, hydrocarbon feed 101 may be combined with recycle stream 119 to form combined feedstream 102, which may be fed into purification process 80, The purification process 80 employs known processes to purify the hydrocarbon feed, which may include fractionation and/or treating the hydrocarbon feed. Treating as used herein includes removing impurities, such as oxygenates, sulfur, and metals, from the hydrocarbon feed. Fractionation as used herein includes removing heavy (e.g., C<T) hydrocarbons and light (e.g., C ₅ ^" ) hydrocarbons. The resulting purified feed 103 may be combined with a dry hydrogen recycle 116 to produce hydrogen rich purified feed 104, which may be combined with oxygenate and/or nitrogenate 105 to produce reactor feedstream 106, as will be described in more detail herein. Reactor feedstream 106 may be pre-heated in a first furnace 1 1, which heats the hydrocarbons to a desired temperature, thereby producing fmst reactor feed 107. First reactor feed 107 may be fed into reactor 10, where the hydrocarbons are contacted with an reforming catalyst under suitable reaction conditions (e.g., temperature and pressure) for aromatizing one or more components in the feed to increase the aromatics content thereof. First reactor effluent 108 comprising aromatics, unreacted feed, and other hydrocarbon compounds or byproducts is recovered from first reactor 10. Hydrocarbon feed 101 and any process piping or process vessels upstream that supply said feed, hydrocarbon recycle stream 119, combined feedstream 102, purification process 80, purified feed 103, dry- hydrogen recycle 1 16, reactor feedstream 106, the furnaces, and the reactor vessels are all within the reaction zone 120 to which oxygenate and/or nitrogenate may be added. In various embodiments, one or more reactor feedstreams of the type shown in Figure 1 are feedstreams upon entering the reactor, and such feedstreams may be further characterized as described herein.

[0031] First reactor effluent 108 may be pre-heated in second furnace 21, which heats the hydrocarbons to a desired temperature, thereby producing second reactor feed 109. Second reactor feed 109 may be then fed into reactor 20, where the hydrocarbons are contacted with a reforming catalyst under suitable reaction conditions for aromatizing one or more components in the feed to increase the aromatics content thereof. Second reactor effluent 1 10 comprising aromatics, unreacted feed, and other hydrocarbon compounds or byproducts is recovered from second reactor 20. [0032] Second reactor effluent 110 may be pre-heated in third furnace 31, which heats the hydrocarbons to a desired temperature, thereby producing third reactor feed 1 1 1. Third reactor feed 1 1 1 may be then fed into reactor 30, where the hydrocarbons are contacted with a reforming catalyst under suitable reaction conditions for aromatizing one or more components in the feed to increase the aromatics content thereof. Third reactor effluent 112 comprising aromatics, unreacted feed, and other hydrocarbon compounds or byproducts is recovered from third reactor 30.

[0033] Third reactor effluent 112 may be pre-heated in fourth furnace 41, which heats the hydrocarbons to a desired temperature, thereby producing fourth reactor feed 113. Fourth reactor feed 113 may be then fed into reactor 40, where the hydrocarbons are contacted with a reforming catalyst under suitable reaction conditions for aromatizing one or more components in the feed to increase the aromatics content thereof. Fourth reactor effluent 114 comprising aromatics, unreacted feed, and other hydrocarbon compounds or byproducts is recovered from fourth reactor 40.

[0034] Fourth reactor effluent 114 may be fed into hydrogen separation process 50 that uses one or more known processes to separate hydrogen recycle 115 from reformate 117. In some embodiments, pure hydrogen may be also produced as a product. Reformate 117 may comprise aromatization reaction products from reactors 10, 20, 30, and 40 (e.g., aromatic and non-aromatic compounds) in addition to any unreacted feed and other hydrocarbon compounds or byproducts, in embodiments, the effluent from the catalytic reforming process may comprise at least about 40 wt % benzene by total weight of the effluent, in other embodiments, the effluent from the catalytic reforming process may comprise from about 40 wt % to about 60 wt % benzene by total weight of the effluent.

[0035] Hydrogen recycle 1 15 may be dried in dryer 60, thereby forming dry hydrogen recycle 116, which may be recycled into purified feed 103. Hydrogen is preferred as diluent as it serves to not only lower the partial pressure of the cyclic hydrocarbon, but also suppresses the formation of hydrogen deficient, carbonaceous deposits, commonly referred to as "coke," on the catalytic composite.

[0036] Reformate 117 is fed to purification-extraction process 70, which separates raffinate recycle 1 19 (e.g., unreacted feed) and reactor byproducts (not shown) from aromatics 1 18. In some embodiments, raffinate recycle 119 may be not combined with hydrocarbon feed 101 to form combined feedstream 102, but may be instead produced as a product. Hydrogen separation processes 50 and purification-extraction processes 70 are well known in the art and are described in numerous patents, including U.S. Patent 5,401 ,386 to Morrison et. al. entitled "Reforming Process for Producing High- Purity Benzene," U.S. Patent 5,877,367 to Witte entitled "Dehydrocyciization Process with Downstream Dimethylbenzene Removal", and U.S. Patent 6,004,452 to Ash et al. entitled "Process for Converting Hydrocarbon Feed to High Purity Benzene and High Purity Paraxylene," each of which is incorporated herein by reference in its entirety. Raffmate recycle 1 19 may be then recycled into feed 101 and aromatics 1 18 are sold or otherwise used as desired. For the sake of simplicity. Figure 1 does not illustrate the byproduct streams that are removed from catalytic reactor system 100 at various points throughout the system. However, persons of ordinary skill in the art are aware of the composition and location of such byproduct streams.

[0037] The composition of the feedstream (e.g., any of streams 101, 102, 103, 104, 106, 107, 108, 109, 110, 1 11 , 112, and/or 1 13) may be a consideration when designing catalytic reforming systems. In an embodiment, the feedstream may comprise non- aromatic hydrocarbons containing at least six carbon atoms. In embodiments, the feedstream may feed to the catalytic reactor (e.g., aromatization reactor), and may comprise any combination of hydrocarbons having from about five to about nine carbon atoms. In embodiments, the feedstream to the catalytic reforming system may comprise a mixture of hydrocarbons comprising C ₆ to Cg hydrocarbons containing up to about 10 wt % of Cs and lighter hydrocarbons (C ₅ ^~ ) by weight of the feed; alternatively, up to about 15 wt % of C ₅ and lighter hydrocarbons (Cs ^" ), and containing up to about 10 wt % of Cg and heavier hydrocarbons (CV) by weight of the feed. Low levels of C9 and C ₅ hydrocarbons may maximize the yield of high value aromatics. Coke formation may be affected by the coke forming components in the feedstream such as cyclopentane and C9 ^r compounds. In some embodiments, the feedstream maximizes the percentage of C ₆ hydrocarbons. Such a feedstream can be achieved by separating a hydrocarbon feedstock, such as a full-range naphtha, into a light hydrocarbon feed fraction and a heavy hydrocarbon feed fraction, and using the light fraction in the feedstream. In embodiments, the feedstream may maximize the conversion of CVconvertibles. In embodiments, the feedstream may comprise up to about 50 wt % Ce-convertibles; alternatively, up to 60 wt % Ce-convertibles; alternatively, up to 70 wt % C ₆ - convertibles by weight of the feedstream.

[0038] In certain embodiments, the feedstream may comprise naphtha. The naphtha feed may comprise a light hydrocarbon, with a boiling range of about 70 °F to about 450 °F. The naphtha feed may contain aliphatic, naphthenic or paraffinic hydrocarbons. These aliphatic and naphthenic hydrocarbons may be converted, at least in part, to aromatics in the catalytic reactor system. While catalytic aromatization typically refers to the conversion of naphtha, other feedstocks can be treated as well to provide an aromatics enriched product. Therefore, while the conversion of naphtha may be one embodiment, the present disclosure can be useful for the conversion or aromatization of a variety of feedstocks such as paraffin hydrocarbons, olefin hydrocarbons, acetylene hydrocarbons, cyclic paraffin hydrocarbons, cyclic olefin hydrocarbons, and mixtures thereof, and particularly saturated hydrocarbons.

[0039] In an embodiment, the feedstock may be substantially free of sulfur, nitrogen, metals, and other known poisons for reforming catalysts. In an embodiment, the feedstock contains less than 100 ppb of sulfur. If present, such poisons can be removed using methods known to those skilled in the art. In some embodiments, the feed can be purified by first using conventional hydrofining techniques, then using sorbents to remove the remaining poisons, as is known to those of skill in the art. Such hydrofining techniques and sorbents may be paxt of the improved catalytic reforming process herein below.

[0040] The reforming reactions may occur under process conditions that thermodynamically favor the dehydrocyclization reaction and limit undesirable hydrocracking reactions. Operating ranges for a typical catalytic reforming process, such as an aromatization process, may include reactor inlet temperatures between 700 °F and 1050 °F; a system pressure between 0 pounds per square inch gauge (psig) and about 500 psig; a recycle hydrogen rate sufficient to yield a hydrogen to hydrocarbon mole ratio for the feed to the reforming reactor zone between 0.1 and 20; and a liquid hourly space velocity for the hydrocarbon feed over the reforming catalyst of between 0.1 and 10. The operating temperatures generally may include reactor inlet temperatures from about 700 °F to about 1050 °F. The subject matter disclosed herein may be used to kinetically model a reforming process (e.g., an aromatization process) in order to enable selection of improved operating conditions for benzene production.

[0041] The reforming catalyst may comprise an inorganic support, a group VIII metal such as platinum, and one or more halides such as fluorine, chlorine, iodine, bromine, or combinations thereof. In some embodiments, the catalyst may comprise Group VIII metals on an inorganic support such as platinum on alumina, Pt/Sn on alumina and Pt/Re on alumina. In other embodiments, the catalyst may comprise Noble Group VIII metals on a zeolitic support such as Pt, Pt/Sn and Pt/Re on zeolitic supports which may comprise a binder and zeolites such as L-zeolites, ZSM-5, silicaiite and beta; and Noble Group VIII metals on alkali- and aika line-earth exchanged L-zeolites. The catalyst may comprise a large-pore zeolite as the inorganic support, charged with at least one Group VIII metal. In embodiments, the Group VIII metal may comprise platinum, which may be more selective for dehydrocyclization and which may be more stable under reforming reaction conditions than other Group VIII metals. In other embodiments, the catalyst may comprise a Group VIII metal such as rhenium, or a Group IV metal or metalloid such as tin.

[0042] In an embodiment, the catalyst may comprise a non-acidic catalyst which may comprise a non-acidic zeolite support as the inorganic support, a group VI I I metal, and one or more halides. Suitable haiides include chloride, fluoride, bromide, iodide, or combinations thereof. Suitable Group VIII metals include iron, cobalt, nickel, ruthenium, rhodium, palladium, osmium, iridium, platinum, or combinations thereof.

[0043] Inorganic supports for reforming catalysts (e.g., aromatization catalysts) can generally include any inorganic oxide. These inorganic supports include bound large pore aluminosilicates (zeolites), amorphous inorganic oxides and mixtures thereof. Large pore aluminosilicates include, but are not limited to, L-zeolite, Y -zeolite, mordenite, omega zeolite, beta zeoiite and the like. Amorphous inorganic oxides include, but are not limited to, aluminum oxide, silicon oxide and titania. Suitable bonding agents for the inorganic supports include, but are not limited to, silica, alumina, clays, titania, magnesium oxide, and combinations thereof.

[0044] The inorganic support may be an aluminosilicate, for example a zeolite. In embodiments, the reforming catalyst may comprise a zeolitic catalyst. Zeolite materials, both natural and synthetic, are known to have catalytic properties for many hydrocarbon processes. The term "zeolite" generally refers to a particular group of hydrated, crystalline metal aluminosilicates. These zeolites exhibit a network of Si0 ₄ and A10 ₄ tetrahedra in which aluminum and silicon atoms are crosslinked in a three- dimensional framework by sharing oxygen atoms. The cavities and channels throughout the crystalline material of zeolites are interconnected and generally can be of a size to allow selective separation of hydrocarbons. In the three-dimensional framework, the ratio of oxygen atoms to the total of aluminum and silicon atoms may be equal to 2. The framework exhibits a negative electro vale ce that typically is balanced by the inclusion of cations within the crystal such as metals, alkali metals, alkaline earth metals, or hydrogen. Thus, zeolites are a group of natural or synthetic hydrated alummosilicate minerals that typically contain alkali and alkaline metals. Zeolites are characterized by a framework structure that encloses interconnected cavities occupied by ion-exchangeable large metal cations such as potassium and water molecules permitting reversible dehydration. The actual formula of the zeolite may vary without changing the crystalline structure. In an embodiment, the mole ratio of silicon to aluminum (Si/Al) in the zeolite may vary from about 1.0 to about 3.5.

[0045] In embodiments, the reforming catalyst may comprise a large pore zeolite. The term 'large-pore zeolite" may be defined as a zeolite having an effective pore diameter of from about 6 Angstroms (A) to about 15 A; alternatively, from about 7 A to about 9 A. Examples of large pore crystalline zeolites are type L-zeolite, zeolite X, zeolite Y, omega zeolite, beta zeolite, ZSM-4, ZSM-5, ZSM-1G, ZSM-12, ZSM-2G, REY, USY, RE-USY, LZ-210, LZ-210-A, LZ-210-M, LZ-210-T, 88Z-24, ZZA-26, SSZ-31, SSZ-33, SSZ-35, SSZ-37, S8Z-41, S8Z-42, SSZ-44, MCM-58, mordenite, mazzite, faujasite, and combinations thereof. In an embodiment, the large pore zeolite may comprise an isotypic framework structure. In an embodiment, the reforming catalyst may comprise L-zeolite.

[0046] L-zeolite, its X-ray diffraction pattern, its properties, and methods for its preparation are described in detail in, U.S. Pat. No. 3,216,789, the content of which is hereby incorporated by reference. Zeolite X is described in U.S. Pat. No. 2,882,244. Mazzite is described in LI.8. Pat. Nos. 4,503,023 and 4,021,447. Mordenite is described in U.S. Pat. No. 4,503,023. Zeolite Y is described in U.S. Pat. No. 3,130,007. U.S. Pat. Nos. 3,216,789; 2,882,244; 4,503,023; 4,021,447; and 3,130,007, are hereby incorporated by reference to show zeolites useful for reforming processes.

[0047] In embodiments, the reforming catalyst may comprise a large pore L-zeolite. L-type zeolite catalysts are a sub-group of zeolitic catalysts. Typical L-type zeolites contain mole ratios of oxides in accordance wi th the following formula:

M.2,¾0*AI.203*xSi02®yH. _? 0 wherein "M" designates at least one exchangeable cation such as barium, calcium, cerium, lithium, magnesium, potassium, sodium, strontium, and zinc as well as non- metallic cations like hydronium and ammonium ions which may be replaced by other exchangeable cations without causing a substantial alteration of the basic crystal structure of the L-type zeolite. The "n" in the formula represents the valence of "M", "x" may be 2 or greater; and "y" is the number of water molecules contained in the channels or interconnected voids with the zeolite.

[0048] Bound potassium L-type zeolites, or KL zeolites, have been found to be particularly desirable. The term KL zeolite as used herein refers to L-type zeolites in which the principal cation M incorporated in the zeolite may comprise potassium. A KL zeolite may be cation-exchanged or impregnated with another metal and one or more halides to produce a platinum-impregnated, halided zeolite or a KL supported Pt- halide zeolite catalyst. In an embodiment, the zeolite may comprise type L zeolite. In embodiments, the reforming catalyst may comprise a potassium L-type zeolite, hereafter referred to as KL-zeolite, which refers to L-type zeolites wherein the principal exchangeable cation M incorporated in the zeolite is potassium.

[0049] One or more Group VIII metals may be added to the catalyst support to form a metallized catalyst support. The metal may be added to the catalyst support by employing a variety of known and conventional techniques, e.g., ion-exchange, incipient wetness, pore fill, impregnation, vapor deposition etc. In embodiments, the platinum and optionally one or more halides may be added to the zeolite support by any suitable method, for example via impregnation with a solution of a platinum-containing compound and one or more balide-containing compounds. In an embodiment, the metal may be added to the catalyst support by impregnation with a metal-containing solution. The metal in the metal containing solution may comprise at least one Group VIII metal; alternatively, Pt, Pd, Rh, Ir, Ru, Os, or combinations thereof. In an embodiment, the metal may comprise platinum that may be added to the catalyst support via contact with a metal-containing solution containing at least one platinum- containing compound. Examples of suitable platinum-containing compounds for contact with the catalyst support include without limitation platinum compounds that form positively charged platinum complex ions in solution such as for example platinum salts such as chlorides and nitrates; platinum complexes with amines; or combinations thereof. For example, the platinum-containing compound can be any decomposable platinum-containing compound including, but not limited to, ammonium tetrachloroplatinate, chloroplatinic acid, diammineplatinum (II) nitrite, bis- (ethylenediarmne)platinum (II) chloride, platinum (II) acetylacetonate, dichiorodiammine platinum, platinum (II) chloride, tetraammineplatinum (II) hydroxide, tetraammineplatinum chloride, and tetraammineplatinum (II) nitrate. In an embodiment, the platinum source may comprise tetraamine platinum chloride (TAPC). The amount of platinum in the metallized catalyst support may range from about 0.1 to about 5 wt %; alternatively, from about 0.1 to about 3 wt %; alternatively, from about 0.3 to about 1.8 wt %.

[0050] In an embodiment, the catalyst may comprise a large pore zeolite support with a platmum-contammg compound and at least one halide. One or more hahdes may be added to the catalyst support by contact with a halide-containing compound to form a halided catalyst support. The halides may be added into the catalyst support separately; alternatively, the halides may be added to the catalyst support at the same time. Such halides may be incorporated during addition of a metal, alternatively, the halides may be incorporated in a separate step that may be pre- or post-addition of the metal, to form a halided, metallized catalyst support. Examples of suitable halides include without limitation fluoride, chloride, bromide, iodide, or combinations thereof. Such halides may be introduced as the ammonium halide compound.

[0051] in an embodiment, the catalyst may comprise a large pore zeolite support with a platinum-containing compound and at least one ammonium halide compound. The ammonium halide compound may comprise one or more compounds represented by the formula N(R)-}X, where X may comprise a halide and where R represents a hydrogen or a substituted or unsubstituted carbon chain molecule having 1-20 carbons wherein each R may be the same or different. In an embodiment, R may comprise methyl, ethyl, propyl, butyl, or combinations thereof. Examples of a suitable organic ammonium compound represented by the formula \{R ) ;X may include ammonium chloride, ammonium fluoride, and tetraalkylammonium halides such as tetram ethyl ammonium chloride (TMAC), tetramethylammonium fluoride (TMAF), teiraethyiarnmoni urn chlori.de, tetraeth.ylarnm.om um fluoride, tetrapropylammoni urn chloride, tetrapropylammoni um fluoride, tetrabutylammonium chloride, tetrabutylammonium fluoride, methyltriethyiammonium chloride, m ethyl triethyl ammonium fluoride, or combinations thereof.

[0052] In an embodiment, the ammonium halide compound may comprise at least one acid halide and at least one ammonium hydroxide represented by the formula N(R') ₄ 0H, where R' may comprise hydrogen or a substituted or unsubstituted carbon chain molecule having 1-20 carbon atoms wherein each R' may be the same or different. In an embodiment, R' may comprise methyl, ethyl, propyl, butyl, or combinations thereof. Examples of a suitable ammonium hydroxide represented by the formula N(R') ₄ OH may include ammonium hydroxide, tetraalkylammonium hydroxides such as tetramethylammonium hydroxide, tetraethylammonium hydroxide, tetrapropylammomum hydroxide, tetrabutylammomum hydroxide, and combinations thereof. Examples of suitable acid halides may include HC1, HF, HBr, HI, or combinations thereof,

[0053] In an embodiment the ammonium halide compound may comprise (a) a compound represented by the formula N(R)4X, where X may comprise a halide and where R represents a hydrogen or a substituted or unsubstituted carbon chain molecule having 1-20 carbons wherein each R may be the same or different, and (b) at least one acid halide and at least one ammonium hydroxide represented by the formula N(R') ₄ OH, where R' may comprise hydrogen or a substituted or unsubstituted carbon chain molecule having 1-20 carbon atoms wherein each R' may be the same or different.

[0054J The halide-containing compound may further comprise an ammonium halide such as ammonium chloride, ammonium fluoride, or both in various combinations with the ammonium halide compounds described previously. More specifically, ammonium chloride, ammonium fluoride, or both may be used with (a) as described previously, a compound represented by the formula \{ . ) ;X. where X may comprise a halide and where R represents a hydrogen or a substituted or unsubstituted carbon chain molecule having 1-20 carbons wherein each R may be the same or different and/or (b) as described previously, at least one acid halide and at least one organic ammonium hydroxide represented by the formula N(R ^" ) ;()! L where R' may comprise a substituted or unsubstituted carbon chain molecule having 1-20 carbon atoms wherein each R' may be the same or different. For example, a first fluoride- or chloride-containing compound can be introduced as a tetraalkylammonium halide with a second fluoride- or chloride-containing compound introduced as an ammonium halide. In an embodiment, tetraalkylammonium chloride may be used with ammonium fluoride.

[0055] In an embodiment, the reforming catalyst may comprise a metallized, halided support and the amount of halide in the catalyst ranges from about 0.05 wt % to about 5.0 wt %. In an embodiment, the halided catalyst support may comprise chloride present in an amount of from about 0.1 wt % to about 5 wt %; alternatively, from about 0.1 wt % to about 3 wt %; alternatively, from about 0.3 wt % to about 1.8 wt %. In an embodiment, the halided catalyst support may comprise fluoride present in an amount of from about 0.1 wt % to about 5 wt %; alternatively, from about 0.1 wt % to about 3 wt %; alternatively, from about 0.3 wt % to about 1.8 wt %. In an embodiment, the haKded catalyst support may comprise both chloride and fluoride, which may be present in a C1:F ratio of from about 1 :10 to 10: 1; alternatively, from about 1 :5 to 5: 1; alternatively, from about 1 :2 to 2:1 .

[0056] In embodiments, the catalyst may comprise a low-acidity silica-bound potassium L-type zeolite support, platinum, chloride, and fluoride. In embodiments, the reforming catalyst may readily convert near-linear C ₆ s (i.e., C ₆ hydrocarbons with no more than one branch), but may not readily convert Ce hydrocarbons having more than one branch, e.g., dimethylpentane. Such readily converted Ce hydrocarbons are referred to herein as "Ce-convertibles." In embodiments, the process to be improved via the disclosed method may comprise a low-acidity silica-bound potassium L-type zeolite support, platinum, chloride, and fluoride. In embodiments, the process to be improved via the disclosed method may comprise a conversion of Ce-convertibles greater than about 75 %; alternatively, greater than about 80 %; alternatively, greater than about 83 %.

[0057] Examples of catalysts and reactor systems suitable for use with the methods described herein are the AROMAX* Process and catalyst technologies available from the Chevron Phillips Chemical Company of The Woodlands, Texas. A specific example of a suitable commercially available reforming catalyst is the AROMAX* II Catalyst available from Chevron Phillips Chemical Company LP of The Woodlands, Texas. Suitable reforming catalysts are also disclosed in U.S. Patent 7,153,801 to Wu entitled "Aromatization Catalyst and Methods of Making and Using Same," and U.S. Patent 6,812,180 to Fukunaga entitled "Method for Preparing Catalyst," each of which is incorporated herein by reference as if reproduced in their entirety.

[0058] In embodiments, the reforming process to be improved via the disclosed method may comprise a series of adiabatic Aromax ^" Reactors with inter-stage heating. Alternatively, the process to be improved may comprise isothermal Aromax' ⁸ Reactors in which a bank of isothermal reactors is heated in a fire box. In an embodiment, the process may comprise six adiabatic Aromax* Reactors with interstage heating: alternative embodiments may comprise seven adiabatic Aromax ^® Reactors with inter-stage heating.

[0059] Figure 2 illustrates a computer system 200 suitable for implementing all, or a portion of, one or more embodiments disclosed herein. All of, or a portion of, the system 200 may be implemented on any particular machine, or machines, with sufficient processing power, memory resources, and throughput capability to handle the necessary workload placed upon the computer, or computers. The computer system 200 includes a processor 202 (which may be referred to as a central processor unit or CPU) that is in communication with memory devices including secondary storage 204, read only memory (ROM) 206, random access memory (RAM) 208, input/output (I/O) devices 210, and network connectivity devices 212. The processor 202 may be implemented as one or more CPU chips.

[0060] It is understood that by programming and or loading executable instructions onto the computer system 200, at least one of the CPU 202, the RAM 208, and the ROM 206 are changed, transforming the computer system 200 in part into a particular machine or apparatus having the novel functionality taught by the present disclosure. It is fundamental to the electrical engineering and software engineering arts that functionality that can be implemented by loading executable software into a computer can be converted to a hardware implementation by well-known design rules. Decisions between implementing a concept in software versus hardware typically hinge on considerations of stabili ty of the desi gn and numbers of units to be produced rath er than any issues involved in translating from the software domain to the hardware domain. Generally, a design that is still subject to frequent change may be preferred to be implemented in software, because re-spinning a hardware implementation is more expensive than re-spinning a software design. Generally, a design that is stable that wil l be produced in large volume may be preferred to be implemented in hardware, for example in an application specific integrated circuit (ASIC), because for large production runs the hardware implementation may be less expensive than the software implementation. Often a design may be developed and tested in a software form and later transformed, by well-known design rules, to an equivalent hardware implementation in an application specific integrated circuit that hardwires the instructions of the software. In the same manner as a machine controlled by a new ASIC is a particular machine or apparatus, likewise a computer that has been programmed and/or loaded with executable instructions may be viewed as a particular machine or apparatus.

[0061 J The secondary storage 204 is typically comprised of one or more disk dri ves or tape drives and is used for non-volatile storage of data and as an over-flow data storage device if RAM 208 is not large enough to hold all working data. Secondary storage 204 may be used to store programs which are loaded into RAM 208 when such programs are selected for execution. The ROM 206 is used to store instructions and perhaps data which are read during program execution. ROM 206 is a non-volatile memory device which typically has a small memory capacity relative to the larger memory capacity of secondary storage 204. The RAM 208 is used to store vol atile data and perhaps to store instructions. Access to both ROM 206 and RAM 208 is typically faster than to secondary storage 204. The secondary storage 204, the RAM 208, and/or the ROM 206 may be referred to in some contexts as computer readable storage media and/or non-transitory computer readable media.

[0062] 170 devices 210 may include printers, video monitors, liquid crystal displays (LCDs), touch screen displays, keyboards, keypads, switches, dials, mice, track balls, voice recognizers, card readers, paper tape readers, or other well-known input devices.

[0063] The network connectivity devices 212 may take the form of modems, modem banks, Ethernet cards, universal serial bus (USB) interface cards, serial interfaces, token ring cards, fiber distributed data interface (FDDI) cards, wireless local area network (WLAN) cards, radio transceiver cards such as code division multiple access (CDMA), global system for mobile communications (GSM), long-term evolution (LTE), worldwide interoperability for microwave access (WiMAX), and/or other air interface protocol radio transceiver cards, and other well-known network devices. These network connectivity devices 212 may enable the processor 202 to communicate with the Internet or one or more intranets. With such a network connection, it is contemplated that the processor 202 might receive information from the network, or might output information to the network in the course of performing the above-described method steps. Such information, which is often represented as a sequence of instructions to be executed using processor 202, may be received from and outputted to the network, for example, in the form of a computer data signal embodied in a carrier wave.

[0064] Such information, which may include data or instructions to be executed using processor 202 for example, may be received from and outputted to the network, for example, in the form of a computer data baseband signal or signal embodied in a carrier wave. The baseband signal or signal embedded in the carrier wave, or other types of signals currently used or hereafter developed, may be generated according to several methods well known to one skilled in the art. The baseband signal and'or signal embedded in the carrier wave may be referred to in some contexts as a transitory signal. [0065] The processor 202 executes instructions, codes, computer programs, scripts which it accesses from hard disk, floppy disk, optical disk (these various disk based systems may all be considered secondary storage 204), ROM 206, RAM 208, or the network connectivity devices 212, While only one processor 202 is shown, multiple processors may be present. Thus, while instructions may be discussed as executed by a processor, the instructions may be executed simultaneously, serially, or otherwise executed by one or multiple processors. Instructions, codes, computer programs, scripts, and/or data that may be accessed from the secondary storage 204, for example, hard drives, floppy disks, optical disks, and/or other device, the ROM 206, and/or the RAM 208 may be referred to in some contexts as non-transitory instructions and/or non-transitory information .

[0066] In an embodiment, the computer system 200 may comprise two or more computers in communication with each other that collaborate to perform a task. For example, but not by way of limitation, an application may be partitioned in such a way as to permit concurrent and/or parallel processing of the instructions of the application. Alternatively, the data processed by the application may be partitioned in such a way as to permit concurrent and/or parallel processing of different portions of a data set by the two or more computers. In an embodiment, virtual ization software may be employed by the computer system 200 to provide the functionality of a number of servers that is not directly bound to the number of computers in the computer system 200, For example, virtualization software may provide twenty virtual servers on four physical computers. In an embodiment, the functionality disclosed above may be provided by executing the application and/or applications in a cloud computing environment. Cloud computing may comprise providing computing sendees via a network connection using dynamically scalable computing resources. Cloud computing may be supported, at least in part, by virtualization software. A cloud computing environment may be established by an enterprise and/or may be hired on an as-needed basis from a third party provider. Some cloud computing environments may comprise cloud computing resources owned and operated by the enterprise as well as cloud computing resources hired and/or leased from a third party provider.

[0067J In a embodiment, some or all of the functionality disclosed above may be provided as a computer program product. The computer program product may comprise one or more computer readable storage medium having computer usable program code embodied therein to implement the functionality disclosed above. The computer program product may comprise data structures, executable instructions, and other computer usable program code. The computer program product may be embodied in removable computer storage media and'or non-removable computer storage media. The removable computer readable storage medium may comprise, without limitation, a paper tape, a magnetic tape, magnetic disk, an optical disk, a solid state memory chip, for example analog magnetic tape, compact disk read only memory (CD-ROM) disks, floppy disks, jump drives, digital cards, multimedia cards, and others. The computer program product may be suitable for loading, by the computer system 200, at least portions of the contents of the computer program product to the secondary storage 204, to the ROM 206, to the RAM 208, and/or to other non-volatile memory and volatile memory of the computer system 200. The processor 202 may process the executable instructions and/or data structures in part by directly accessing the computer program product, for example by reading from a CD-ROM disk inserted into a disk drive peripheral of the computer system 200. Alternatively, the processor 202 may process the executable instructions and/or data structures by remotely accessing the computer program product, for example by downloading the executable instructions and/or data structures from a remote server through the network connectivity devices 212. The computer program product may comprise instructions that promote the loading and ^' or copying of data, data structures, files, and/or executable instructions to the secondary storage 204, to the ROM 206, to the RA : 208, and'or to other non -volatile memory and vol atile memory of the computer system 200.

[0068] In some contexts, the secondary storage 204, the ROM 206, and the RAM 208 may be referred to as a non-transitory computer readable medium or a computer readable storage media. A dynamic RAM embodiment of the RAM 208, likewise, may be referred to as a non-transitory computer readable medium in that while the dynamic RAM receives electrical power and is operated in accordance with its design, for example during a period of time during which the computer 200 is turned on and operational, the dynamic RAM stores information that is written to it. Similarly, the processor 202 may comprise an internal RAM, an internal ROM, a cache memory, and/or other internal non-transitory storage blocks, sections, or components that may be referred to i some contexts as non-transitory computer readable media or computer readable storage media.

[0069] In an embodiment, the computer system 200 (through, e.g., the network connectivity devices 212) is capable of facilitating communications between the processor 202 and components of the catalytic reactor system 100 of Figure 1, including reactors 10, 20, 30, and 40; furnaces 1 1 , 21 , 31, and 41 ; hydrogen separation process 50; dryer 60; purification-extraction process 70; purification process 80; other components such as pumps and value, or combinations thereof. In an embodiment, the network connectivity devices 212 may be wired devices such as Ethernet cards, USB interface cards, etc., or combinations thereof, and the communications between the processor 202 and the components of the catalytic reactor system 100 may be wired network communications. Alternatively, the network connectivity devices 212 may be wireless devices such as WiFi, Bluetooth, etc. and the communications between the processor 202 and the components of the catalytic reactor system 100 may be wireless network communications.

[0070] In an embodiment, components of the catalytic reactor system 100 (e.g., reactors 10, 20, 30, and 40; furnaces 11, 21, 31, and 41 ; hydrogen separation process 50; dryer 60; purification-extraction process 70; and purification process 80; valves; pumps; etc.) may be coupled to the computer system 200 such that data obtained from one or more components may be stored and/or used by the computer system 200. In an embodiment, the run-length and net present amount of benzene recovered may be kinetically modeled as a function of the H2/HC ratio, LH8V, and/or conversion of Q- convertibles using the computer system 200. After kinetic modeling, the H2/HC ratio, LBSV, and conversion of Ce-eonvertibies which maximize the run-length and/or net present amount of benzene recovered may be selected by a user of the computer system 200 or by the computer system 200 itself (e.g., under system or user-defined guidelines), for example. The computer system 200 may communicate with one or more components of the catalytic reactor system 100 so as to set or adjust operating parameters to maximize the run-length and/or net present amount of benzene.

[0071] Figure 3 illustrates the deactivation kinetic model 310 as disclosed herein utilized in conjunction with a computer system 300 (for example, computer system 200 as described herein) and a catalytic reactor system 330 (for example, as shown in Figure 1 and described herein). All or a portion of the deactivation kinetic model 310 as disclosed herein may be located on computer system 300 or on a combination of computer systems 300 (e.g., a distributed computing system). A control program. 320 may also be located in whole or in part on the computer system 300. The deactivation kinetic model 310 and/or control program 320 may communicate with the catalytic reactor system 330 through the computer system 300 as described herein. [0072] in an embodiment, the deactivation kinetic model 310 may comprise an input, an output, a model unit 312, a selection unit 314, and a sensitivity unit 316. The input and output may move data into and out of the kinetic model 310, and may occur automatically (e.g., in computer system 300) or manually (e.g., through user input). The deactivation kinetic model 310 may predict run-length and/or net present amount of benzene of the catalytic reactor system 330 as a function of H2/HC ratio, LHSV, and/or conversion of Ce-convertibles in model unit 312. The selection unit 314 of the deactivation kinetic model 310 may then be used to select (either automatically or manually) H2/HC ratio, LHSV and/or conversion of C ^' e-convertibles whic maximize run-length and or net present amount of product (e.g., benzene). Sensitivity of the values predicted by the model unit 312 may be determined, in the sensitivity unit 316 (e.g., utilizing Equation 23), and the values calculated in the sensitivity unit 316 may be output (e.g., to a display of the computer system 300, to external storage 302, or combinations thereof). The values of H2/HC ratio, LHSV, and/or conversion of Cg- coiivertibles may be output to other components of computer system 300 (e.g., ROM, RAM, internal hard drive, external storage 302, a display, combinations thereof, etc.), or the values may be output through computer system 300 to the catalytic reactor system 330,

[0073] In one or more embodiments, the computer system 300 may have a control program 320 which comprises one or more control routines, such as a fuzzy logic control routine 324 and/or a PID control routine 322. The control program 320 may control the values of H2/HC ratio, LHSV, and/or conversion of Ce-eonvertibles selected by the deactivation kinetic model 310. The control program 320 may alternatively or additionally control other parameters of the catalytic control system 330 so as to maintain the system 330 with the values selected by or with the deactivation kinetic model 310. In Figure 3, the control program 320 may control (e.g., open, close, adjust a degree of operation, or combinations thereof) components 332 and 334 (e.g., pumps or valves) to control the operating conditions (e.g., pressure, temperature, flow rate) of the HC stream and H2 stream entering reactor 338 through feedstream 331. For example, the fuzzy logic control routine 324 and PID control routine 32.2 of the control program 320 may communicate analog input (AI) and analog output (AO) with components (e.g., components 332 and 334) of the catalytic reactor system 330. Data from the control program 320 may be stored in the computer system 300, the external storage 302, or combinations thereof. The effluent stream 336 of the catalytic reactor system 330 may contain the desired product produced according to the selected value of H2/HC ratio, LHSV, and/or conversion of C6-convertibl.es which maximize run-length and/or net present amount of product produced in catalytic reactor system 330.

[0074] In an embodiment, the computer system 300 may communicate with the catalytic reactor system 330 including reactor 338 and components 332 and 334. In an embodiment, components 332 and 334 and/or reactor 338 of the catalytic reactor system 330 may be coupled to the computer system 300 such that data obtained may be stored and/or used by the computer system 300.

[0075] In embodiments, components or conditions of the catalytic reactor system 330 may change or be modified. The run-length and net present amount of benzene may be again kinetically modeled to account for changes in the catalytic reactor system 330. For example, the catalytic reactor system 330 may change as a result of addition of fresh catalyst, regeneration of deactivated catalyst, or the addition or use of co- catalysts to the system 330. In such scenarios, coefficients and values in the equations of the disclosed deactivation kinetic model may be calibrated and/or adjusted to more accurately predi ct the run-length and/or net present amount of benzene for the catalytic reactor system 330. The computer system 300 may be utilized to adjust the equations, to kinetically model the changed system 330 as a function of H2/HC ratio, LHSV, and conversion of Ce-convertibles, and to predict new run-lengths and/or net present amount of benzene scenarios. A user or the computer system 300 may then select new values H2/HC ratio, LHSV, and conversion of-convertibles which maximize the production of benzene and/or run-length in the catalytic reactor system 330 under changed components or conditions. In alternative embodiments, none or only a portion of the disclosed embodiments may be utilized on the computer system 300 so as to account for changes in the catalytic reactor system 330.

[0076] In an embodiment, the computer system 300 may kmetically model different catalyst systems to predict which catalyst system has improved LHS V, H2 HC ratio, and conversion of Gs-coiivertibles for operating conditions which are available at a given site where a reactor system (e.g., catalytic reactor system 330) may be located. For example, catalyst systems which require high H2/HC ratios for maximum run- length and/or net present amount of benzene (as determined via the deactivation kinetic modeling disclosed herein) may not be suitable for sites which have limited access to required amounts of hydrogen; thus, other systems may be kinetically modeled (e.g., by adjusting coefficients and values in the equations of the disclosed deactivation kinetic model) to find a suitable catalytic reactor system which does not require unobtainable amounts of hydrogen. In alternative embodiments, none or only a portion of the kinetic modeling associatd with changes and'or modifications may be utilized on the computer system 300 so as to kineticaily model different catalyst systems.

[0077] The computer system 300 may be capable of receiving internal data and'or external data and generating and delivering signals to components of the catalytic reactor system 330. For example, the computer system 300 may receive automated instructions and ^' or manual instructions from a user input, and may send signals to, for example, the reactor 338, components 332 and 334, or combinations thereof of the catalytic reactor system 330 based on the kinetic modeling performed according to the disclosed subject matter. Thus, the computer system 300 may be coupled to the catalytic reactor system 330, and may display the selected and operating values for H2/HC ratio, LHSV, and ^' or conversion of Q-conv rtibles. As such, the computer system 300 may be capable of affecting various functions of catalytic reactor system 330. In alternative embodiments, the catalytic reactor system 330 may be controlled manually or by a system(s) in addition to or other than computer system 300,

[0078] As disclosed herein, improvement of the catalytic reforming process may comprise kineticaily modeling the catalytic reforming process. Typically, one or more components of the feedstream undergoes one or more reforming reactions to produce aromatics. Some of the hydrocarbon reactions that occur during the reforming operation include the dehydrogenation of cyclohexanes to aromatics, dehydroisomerizaiion of alkylcyclopentanes to aromatics, and dehydrocyclization of acyclic hydrocarbons to aromatics. A number of other reactions also occur, including the dealkylation of alkyibenzenes, isoraerization of paraffins, and hydrocracking reactions which produce light gaseous hydrocarbons, e.g., methane, ethane, propane, and butane. These reforming reactions may be kineticaily modeled in a deactivation kinetic model as described herein, for example, by using one or more aromatization reaction rate equations. While catalyzing the foregoing reforming reactions, the reforming catalyst itself may be subjected to a number of reactions that cause the catalyst to deactivate. These deactivation reactions of the reforming catalyst may be kineticaily modeled in a deactivation kinetic model as described herein, for example, by using one or more catalyst deactivation equations. Also, the overall process may be subject to heat balance, for example by using a heat balance equation. Kinetic models including catalyst deactivation equations are referred to herein as deactivation kinetic models. Accordingly, in an embodiment, the deactivation kinetic model may comprise one or more equations of at least three basic types, namely aromatization reaction rate equations, heat balance equations, and catalyst deactivation equations.

[0079] In an embodiment, a single rate equation may be used to represent the aromatization reaction. For example, an irreversible first order aromatization rate equation for a plug flow reactor may be written as:

k(t) - p(t)exp(-E RT) - -υ[(1 + δ) ln(l -J) + 5/] (1) wherein E _r is the activation energy of the reaction, p(t) is the pre-exponential factor of the Arrhenius rate constant (which is the intrinsic catalyst activity), T is the reaction temperature, υ is the gas hourly space velocity (GHSV),/is the fraction conversion of the convertibles, and 6 is a dimensioniess parameter that accounts for the change in the number of moles of gaseous species as a result of the aromatization reactions. Note that p is a function of time because of catalyst deactivation.

[0080] The conversions of paraffins, isoparaffins, and napthenes as a combined feedstream are all lumped into a single equation in Equation (1). In embodiments, the deactivation kinetic model may be more complex, including separate rate equations for the various components of the feedstream. For example, separate rate equations similar to Equation (1 ) may be written for the Co-convertibles as well as for the C ₇ and Cg- convertibles. The term "C _x -convertibles" as used herein refers to those components containing x carbons that are readily converted to aromatics via a catalytic reforming process. For example, when the reforming process utilizes a catalyst having platinum dispersed on a zeoiitic catalyst support, the Ce-convertibles can be all C ₆ hydrocarbons excluding the highly substituted hydrocarbons, such as the dimethylbutane isomers. By utilizing separate rate equations for the various components of the feedstream, the model can be run to target certain conversion of CVconvertibles as it may be often practiced in a plant. As the catalyst deactivates during the catalytic reforming processes, the inlet temperature of the reactor(s) may be raised to maintain the overall conversion at the process outlet, e.g., reactor effluent. By fitting the model to plant or laboratory data, kinetic parameters may be revised. Including hydrocraeking reactions and separate kinetic equations for the various components of the feedstream allows for a better description of the catalytic reactor temperature profile over the run-length of the catalyst. Accordingly, in another embodiment, a plurality of rate equations such as Equations (2), (3), and (4) may be used to represent the aromatization reaction. [0081] An irreversible first order aromatization reaction rate equation of the conversion of CVconvertibles in a plug flow reactor may be written as:

p ₆ (t) exp( - !·, ,, RT) = -υ[( 1 + δ _ό ) ln( 1 » + δ _{6 6} ] (2) where E _r , ₆ is the activation energy of the reaction, p ₆ is the pre-exponential factor, υ is the GHSV, fe is the fraction of conversion for the (^convertibles, and δ _ό accounts for the change in the number of moles of gaseous species as a result of the aromatization of the Ce-convertibles. ^ is a function of time due to catalyst deactivation.

[0082] A rate equation similar to Equation (2) for the C ₇ & Cg-convertibles may be written as:

p7&s(t) exp(-E _r>7&8 RT) = -υ[(1 + δ _7&8 ) Ml -./ _7& g) + δ _{7&8 7&8} ] (3) where E _r>7 & ₈ is the activation energy of the reaction, p ₇ &g is the pre-exponential factor, υ is the GHSV, is the fraction conversion of the C- and Cg-convertibles, and δ ₇ & ₈ accounts for the change in the number of moles of gaseous species as a result of the aromatization of the C ₇ and Cg-convertibles. p ₇ &g is a function of time due to catalyst deactivation.

[0083] In embodiments, the deactivation kinetic model also accounts for the exothermic hydrocracking reaction, which may be significant near end of run (EOR) conditions, and for accurate!)' determining the run-length of the catalyst. An equation describing the cracking of a paraffin molecule into two smaller molecules is:

p _hc it)exp(~E _ht /I T) = -υ1η(1 -J) (4) where E _hc is the activation energy of the hydrocracking reaction and p} _lc is the preexponential factor of the hydrocracking reaction. The hydrocracking of paraffins does not alter the number of molecules in the product. In embodiments, it may be assumed that hydrocracking does not alter the amount of convertibles and that all hydrocarbons in the feedstream can be hydrocracked.

[0084] In an embodiment, a single heat balance equation may be used for the reactor, for exam le:

C _P AT =./C ₀ H _r (5) wherein C _p is the heat capacity of the stream, ΔΤ is the change in the reactor temperature as the result of the reaction, Co is the inlet concentration of the convertibles in the stream, and H _r is the h eat of th e aromatization reaction.

[ 0085] in a similar manner to the use of separate kinetic equations for the various components of the feedstream, separate heat balance equations may be written. Accordingly, in an embodiment, a plurality of heat balance equations such as Equation (6). The heat balance for the aromatization reactions of C _x convertibles, wherein x = C ₆ , C ₇ , or Cg may be given by an equation analogous to Equation (5):

C AT, ./t \H, (6) where C _p is the heat capacity of the stream, ΔΤ _ν is the change in the reactor temperature as the result of the aromatization reactions, C _x is the inlet concentration of the C _x convertibles in the stream, and H _r is the heat of reaction for the aromatization reactions.

[0086] A similar equation may be written for the hydrocracking reaction:

cyvr,,., ./C\H _h , (7) wherein ΔΊΊ ₁₀ i s the change in the reactor temperature as the result of the hydrocracking reaction, C _x is the concentration of C _x hydrocarbons in the feedstream, and H;., _c is the heat of reaction for the hydrocracking reaction

[0087] Predictions of the kinetic models above (Equations (1) - (7) or Equations (2), (3), (4), (6), and (7)), i.e., kinetic models without catalyst deactivation, are restricted to those of start of run (SOR) conditions, and may be used, for example, to determine SOR temperature, hereinafter referred to as T ₀ . In order to project reactor and catalyst performance over time in plant conditions and to predict the run-length of the catalyst, catalyst deactivation may be included in the kinetic modeling herein disclosed for the determination of improved reforming process conditions, and use of such deactivation kinetic models is disclosed further herein.

[0088] In embodiments, the decay of reforming catalyst activity with on-stream time t (i.e., catalyst deactivation) may be described by a first order rate equation: da(t)/dt = -k _d (T)a(t) (8) wherein aft) is the normalized catalyst activity which is defined as p _x (t)/p _x (t=0) with x=6 or 7&8 and k _< j(T) is the deactivation rate constant. In embodiments, it may be assumed that Equation (8) applies to C ₆ , C and Cg aromatization and that the hydrocracking reaction can also be deactivated but at a much slower rate. The deactivation of the hydrocracking may be described as a(t) ^a where a is much smaller than unity and is one of the fitting parameters.

[0089] The deactivation rate constant of Equation (8), k.jf Π. can be expressed in the Arrhenius form:

k _d (T) = A _d exp(-E _d / R.T) (9) where T is the catalyst temperature, E _d is the activation energy of the reaction that leads to catalyst deactivation, and R is the gas constant. The deactivation amplitude Ad is a function of H2/HC ratio. In embodiments, Ad and H2/HC ratio may be linearly correlated.

[0090] Equations (1) to (9) are disclosed herein. In one or more embodiments, a deactivation kinetic model may comprise Equation (1 ), Equation (2), Equation (3), Equation (4), Equation (5), Equation (6), Equation (7), Equation (8), Equation (9), or combinations thereof. In an embodiment, a deactivation kinetic model may comprise Equations (1) to (9). In an additional or alternative embodiment, a deactivation kinetic model may comprise Equations (I), (4), (5), (7), and (8). The deactivation kinetic model may be represented by first order rate equations for the conversion of light naphtha to aromatics, and hydrocrackmg, two heat balance equations, and a catalyst deactivation equation represented by a first order decay of the catalyst activity, in an additional or alternative embodiment, a deactivation kinetic model may comprise Equations (2), (3), (4), (5), (7), and (8). The deactivation kinetic model may be represented by a rate equation for the conversion of CVconvertihles, a rate equation for the C ₇ and Cg-convertibles, a rate equation for exothermic hydrocracking, a heat balance for the aromatization reactions, a heat balance for the hydoeracking reaction, and a catalyst deactivation equation. Deactivation kinetic models, such as those disclosed herein, may be applied to a single reactor vessel, or may be applied to a plurality of reactor vessels, for example, a series of reactors. Individual parameters within a given equation (e.g., E _r , p, C _p , H _r , etc.) may be adjusted independently for each reactor in a plurality of reactors to provide a more accurate representation of the performance of each reactor.

[0091] A deactivation kinetic model may be used to calculate a fouling curve for a given reforming catalyst. A fouling curve may be defined as the catalyst temperature, T, for maintaining a given conversion of Ce-convertibles as a function of on-stream time, t. The intersection between the fouling curve and the prescribed end-of-run (EO ) catalyst temperature gives the catalyst run-length. The EOR catalyst temperature may be a function of catalyst and/or process equipment (e.g., reactor metallurgy) maximum operating temperatures. For example, the EOR temperature may range from about 900 °F to about 1100 °F; alternatively, from about 950 °F to about 1050 °F; alternatively, from about 1000 °F to about 1020 °F. In an embodiment, the EOR may be about 940 °F. The run-length of the catalyst may be further determined as a function of the conversion of Ce-convertibles in the feedstream, LHSV of the feedstream, and H2/HC ratio in the feedstream. [0092] After setting up the three dimensional process variable space spanned by conversion of Ce-convertibles, LHSV and H2/HC ratio, one may calculate the fouling curve for each point in that space with a deactivation kinetic model [e.g., a simplified model employing Equations (1), (4), (5), (7), and (8) or a more complex model employing Equations (2), (3), (4), (5), (7), and (8)]. For a given EOR temperature, the run-lengths of the reforming catalyst may be read from the various fouling curves. Accordingly, the run-length of the catalyst may be established as a function of the LHSV of the feedstream, the H2/HC ratio in the feedstream, and the conversion of C ₆ - convertib!es in the feedstream.

[0093] For a given reactor system (e.g., a plurality of reactors in series), the equations of a deactivation kmetic model may be coupled and sol ved numerically using methods recognized by those of skill in the art with the aid of this disclosure, for example, using a Visual Basic Language equation solver application wit an Excel spreadsheet interface. As one can see, just a few discrete values of each process variable (e.g., LHSV, H2/HC ratio, and conversion of Ce-convertibles) would mushroom into over hundreds of separate nins by a deactivation kinetic model. Such a brute force approach, while effective, is tedious and time consuming. Therefore, the fouling curve(s) generated by a deactivation kinetic model may be approximated with an analytical soiution(s) to the deactivation kinetic model, for example, using an approximated analytical equation, in embodiments, solving of a deactivation kinetic model may be simplified via an approximated analytical equation depending on the system reactor configuration, for example, an isothermal reactor system or a series of adiabatic reactors, as described in more detail below.

[0094] Once the run-length of the reforming catalyst has been determined, the amount of benzene produced per unit of catalyst (e.g., pounds of benzene per pound of catalyst) during the life cycle (e.g., run-length) may be calculated from the run-length of the catalyst and the benzene yield per month.

[0095] In an embodiment, the amount of benzene produced per unit of catalyst may comprise calculating the net present amount of benzene produced.

[0096] In an embodiment, calculating the net present amount of benzene produced may comprise using the run-length of the reforming catalyst and a benzene yield per month of the reforming catalyst to calculate an amount of benzene produced over the run-length, per unit of the reforming catalyst. Calculating the net present amount of benzene produced may further comprise discountmg the amount of benzene produced over the run-length of the reforming catalyst per unit of reforming catalyst to determine the net present amount of benzene produced per unit of reforming catalyst over the run- length of the reforming catalyst.

[0097] In an additional or alternative embodiment, calculating the net present amount of benzene produced may comprise calculating the number of pounds of benzene produced per pound of reforming catalyst per hour, N, which may be calculated by:

N = 0.907 f wee LHSV (d _if / 4)

(10)

where djf, d _c , f, and wee are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of Ce-eonvertibles, and the weight fraction of C ₆ -convertibles in the feedstream. The factor of 0.907 may represent the contraction in the molecular weight as a Ce paraffin converts to benzene. The total amount of benzene produced per pound of catalyst in the run-length of the catalyst may then be calculated by multiplying N by the catalyst run-length in hours. In an embodiment, N may be greater than about 1 ; alternatively, greater than about 1 .5; alternatively, greater than about 2; alternatively, greater than about 2.5; alternatively, greater than about 3; alternatively, greater than about 3.5; alternatively, greater than about 4; alternatively, greater than about 4.5; alternatively, greater than about 5; alternatively, greater than about 5,5; alternatively, greater than about 6. In an embodiment, N may be in a range from about 1.5 to about 6; alternatively, from about 2 to about 5.

[0098] in an additional or alternative embodiment, calculating the net present amount of benzene product may comprise calculating the net present pounds of benzene produced per pound of reforming catalyst over the run-length of the reforming catalyst, NPPB, which may be calculated by:

NPPB = (N/720) { 1 - (1 -!).,, V /D _m

(1 1)

where D _m is the monthly discount rate and n is the run-length of the reforming catalyst in units of months. The monthly discount rate D _m is related to the yearly discount rate by (I - D _y ) = (1 - yearly discount rate) ^' In an embodiment, NPPB may be greater than about 14,000; alternatively, greater than about 16,000; alternatively, greater than about 18,000. In an embodiment, NPPB may be in a range of about 14,000 to about 30,000; alternatively, about 16,000 to about 24,000. A time-value formula may be used to correct future benzene value to a present value. Because the value of one pound of benzene at the present is worth more than that a few years from now, the time value of benzene must be taken into account. Alternatively, one can consider that one pound of benzene produced in the future would have been actually less than one pound of benzene had it been produced today. The percent reduction may be the same as the discount rate in the value of the benzene.

[0099] In an embodiment, the deactivation kinetic mode! may be simplified for a catalytic reforming process comprising an isothermal reactor system. The deactivation kinetic model may be modified for isothermal calculations by setting all the heats of reactions to zero. Because at high conversions the H2 HC ratio in an isothermal reactor increases significantly as one moves toward the reactor outlet, A _d may be treated as an implicit function of the distance from the reactor inlet. Also, A _d at each location along the reactor axis may be a weak function of the on-stream time, t, because the H2/HC ratio at a given location changes as the catalyst deactivates. The simplification follows from the fact that the fouling curve calculated with A _d varying along the reactor can be represented by the analytical solution of an isothermal reactor with a single A _d for the whole catalyst bed. Accordingly, the change of the foulmg curve with on-stream time, t, for a fixed conversion may be given by the analytical solution represented by:

(12)

here k _d represents the average of 1¾ across the reactor; E _r is the average activation energy of the aromatization reaction of the CV, Cy and Cg-convertibles; and T ₀ is the start-of-run (SOR) temperature in Kelvin. T ₀ may depend on LHSV, H2/HC ratio, and conversion of CVconvertibles. For example, T ₀ may be calculated as a function of GHSV (υ) and conversion of convertibles (f) using Equation (1); alternatively, T ₀ may be calculated as a function of GHSV (υ) and conversion of Ce-convertibies ( ) using Equation (2); alternatively, T ₀ may be calculated as a function of LHSV (υ) and conversion of C ₇ - and Cg-convertibles if) using Equation (3); alternatively, T ₀ may be calculated as a function of GHSV ⁷ (υ) and conversion of convertibles if) using

Equation (4). According to Equation (9), k _d may depend on A _d , which is the average of A _d across the reactor and may be determined as a function of H2/HC ratio and conversion of CVconvertihles. For example, Aa may be calculated from the H2/HC ratio and conversion of Ce-convertibles by the following equation:

A _d =1.343x H) ^{1 1} -5.282X 10 ¹⁰ /-1.456x 10 ^{i 0} (H2/HC) _{! ! ei} (13) where f is the fraction representing conversion of Ce-convertibles and (H2/HC) _/ „ _/ei is the hydrogen-to-hydrocarbon mole ratio in the inlet of a reactor. Without being limited by theory, it is believed the dependence of Ad on conversion of Ce- convertibfes may be due to the increase of H2/HC ratio in the reactor when more aromatization reaction occurs,

[00100] In embodiments, the run-length (t) of the reforming catalyst may be calculated as a flmction of conversion of C6-convertibles /), LHSV, and H2/HC ratio using Equations (12) and (13) and the values of T ₀ from the isothermal reactor deactivation kinetic model. For example, for a given conversion of Ce-convertibles, H2/HC ratio, and LHSV, Equation (13) may be substituted for A _d (e.g., in Equation (9), where k _d may be calculated), EOR temperature may be selected for T, and the run- length of the reforming catalyst (t) may be calculated from Equation (12) when T ₀ is calculated from, for example, one or more of Equations (l)-(3). The process of determining the run- length of the catalyst for each point in the process variabl e space is reduced to the calculation of T ₀ at that point. Using a different correlation between A _d and H2 HC ratio may result in a different linear function for Equation (13). However, the deactivation kinetic model calculations for T ₀ need not be repeated in order to recalculate the run-length of the catalyst because, by definition, T _c . is independent of Ad or Ad.

[00101 ] As shown in Example 1 herembelow, several trends can be established. First, for a given H2/HC ratio and LHSV, a lower conversion of CVconvertibles may give a larger NPPB. When a reactor is operated at a lower conversion, the nm-length of the catalyst may be longer. The longer run-length may be more than sufficient to compensate the time-value discount. However, there may be a limit on how low the conversion can be attained because of other operational requirements such as mogas octane number or the loading of the separation units downstream of the reactor. The results in Example 1 herembelow suggest a greater NPPB may be attained by operating the reactor at the lowest conversion that is consistent with other plant requirements. In embodiments, improving the reforming process may comprise selecting the lowest conversion of Ce-eonvertibles consistent with other plant requirements. [00102] Second, for a given conversion, it is advantageous to run the catalytic reforming process at high LHSV. Although running at a higher LHSV means operating the reforming catalyst at a higher temperature, thus a shorter catalyst life, the amount of benzene produced per month is also proportionally higher. The end result is a larger NPPB. In practice, there is an upper limit on the LHSV because the aromatization reaction likely will become diffusion limited at high LHSV, More importantly, a high LHSV can lead to short run-length and frequent catalyst replacement. The cost of the associated down time as well as the lead time for the manufacture of the fresh catalyst may put a cap on high LHSV values.

[00103] In terms of logistics, it may be feasible to have a catalyst replacement cycle with a frequency of one year or longer. Using the run-length of the reforming catalyst of one year as a cut-off, according to Example 1 herembelow, a high NPPB may be achieved with a LHSV as high as 4, so long as mass transport does not become the rate limiting step at this LHSV. In embodiments, diffusion limitation occurs at a LHSV of about 5. In embodiments, the reforming process may be improved in terms of NPPB by running at a LHSV of about 4. In embodiments, the conversion of (^convertibles may range from about 83% to about 93% and selected values of LHSV of the feedstream and H2/HC ratio may be about 4.

[001Θ4] Surprisingly, as shown in Example 1 herembelow, in embodiments, there is little advantage to using a high H2/HC ratio unless it is larger than 3. Although intuitively a higher H2/HC ratio may be thought to lead to a lower fouling, it does not necessarily translate into a longer run-length. As shown in the results in Example 1 hereinbelow, run-lengths of the catalyst for H2/HC ratios between 0.5 and 2 are similar. For an isothermal reactor, a higher H2/HC ratio results in a higher overall gas-hourly-space-velocity in the reactor, which requires a higher reactor temperature to maintain the conversion. In embodiments, for a given conversion of C ₆ - convertibles and LHSV, a higher value of NPPB is attained by selecting a H2/HC ratio of greater than 2 or less than about 1. In embodiments, improving the reforming process may comprise selecting a H2/HC ratio of greater than about 3; alternatively, improving the reforming process may comprise selecting a H2/HC ratio of about 4.

[00105] Catalyst activity may be expressed as a function of A _d , Ed, E _T , and T ₀ . A very simple equation for normalized catalyst activity, a(f), can be derived using the fact that, at a fixed conversion; a(t) = exp{(E,/R)(l/T(t) - ^■ I/T ₀ )> .

(14)

Substituting Equation (14) into Equation (12), the normalized catalyst activity a(t) of an isothermal reactor over time is given by:

a{i ) { l-((&(T ₀ ) t)/a)} ^a

(15)

where a= E ₃ /E _d and k _d is given by an equation similar to Equation (9) except that A _tl is replaced by A _d . The run-length is, therefore, dictated by k _d (T ₀ ), That is, a larger k _d (I ^' _o ) leads to a shorter run-length. Referring to Equation (9), 1¾(Τ ₀ ) is a product of A _d and the Arrhenius factor exp(-E _d /RT ₀ ). An increase in the H2/HC ratio reduces A _d but increases the Arrhenius factor because of a higher T ₀ . These two opposing factors govern whether k _d (T _G ) increases or decreases with the H2/HC ratio, as further discussed in Example 1 hereinbelow.

[00106] in an embodiment, the deactivation kmetic model may be simplified for a catalytic reforming process comprising a series of adiabatic reactors to improve benzene production per catalyst usage. As for the isothermal reactor design discussed hereinabove, the deactivation kinetic model is used to generate the run-length of the catalyst as a function of conversion of C ₆ -convertibles, LHSV and H2/HC ratio. The run-length is in turn used to determine the amount of benzene produced per catalyst life cycle after factoring in the time value of the benzene.

[00107] The deactivation kinetic model as described above is used to generate the fouling curve, which for a series of adiabatic reactors is defined as the inlet temperature, T, of the first reactor for maintaining a given conversion of Ce-convertibles at the last reactor outlet as a function of run-length, t. In embodiments, equal inlet temperatures are assumed for all the reactors. The fol lowing discussion is for a six-reactor configuration with catalyst volume ratios of 10: 10: 10:20:20:30, although it is understood that the disclosed subject matter of kvnetically modeling the catalytic reforming process is equally applicable to alternative numbers of reactors having different catalyst volume ratios.

[001Θ8] To extend the deactivation model to run conditions with various H2/HC ratios, a correlation between the deactivation amplitude, A _d (in units of hr ^"1 ), of each of the six reactors with the H2/HC ratio to each of the reactors may be represented by:

Λ,; 1.231xl O ⁿ - 1 .507xl 0 ¹⁰ (H2/HC) _M et

(16) Because the heat capacity of the gas mixture in the reactor is also affected by the H2/HC ratio, it may be approximated by a linear combination of the heat capacities of the hydrogen and hydrocarbons.

[00109] As mentioned above, to determine the run-length of the reforming catalyst, one first calculates the fouling curve. The intersection between the fouling curve and the prescribed end-of-run (EQR) inlet temperature gives the run-length of the catalyst. However, instead of repeating the time consuming calculations of the fouling curve for various process conditions, one can seek an approximated analytical solution to the model describing the fouling curve, in an embodiment, an approximated analytical solution may be used to predict the fouling curve of a series of adiabatic reactors. As in the isothermal embodiment, the approxirnation greatly reduces the number of calculations required. Once the parameters in this analytical solution are established as a function of the process variables, determination of the various run-lengths and benzene yields is reduced to the determination of these parameters.

[00110] When comparing the fouling curves of an isothermal reactor to an adiabatic reactor under the same process conditions, one notices that the fouling curve of the latter has a smaller initial slope, namely a smaller initial fouling rate at the start of run (SOR). However, near EOR the fouling rate of the adiabatic reactors increases much more rapidly, and the fouling curve reaches EOR earlier. Without intending to be limited by theory, this difference can be explained as follows. Because of the large endothemi in an adiabatic reactor, only a small part of the catalyst bed is used at a given time. Initially, only the part of the catalyst bed near to the reactor inlet is utilized, and it is the first to be deactivated. The deactivation front moves down the bed with on- stream time. As long as the front is not close to the reactor outlet, one does not "sense" the catalyst deactivation and the catalyst fouling rate appears to be small. As the front moves close to the reactor outlet, there is not much catalyst activity left in the entire reactor, and rapid increases in the inlet temperature are necessary to maintain the conversion. An analytical solution describing this type of behavior may be represented by:

T(t) = T _a -p/(t ) ln(l~t/ )

(17)

where T(t) is the inlet temperature of at least one adiabatic reactor, t is on-stream time, tmax represents the theoretical maximum run-length of a catalyst (e.g., the time at which the catalyst has near zero activity and the reactor temperature approaches infinity), β/ ^' is the initial slope of the fouling curve and T _a is related to the SOR T by T _a = T(t = 0) -β. β is an adjustable parameter describing the fouling of the catalyst near start-of run (SOR), and β may be adjusted to fit the analytical solution (e.g.. Equation (17)) to the fouling curve(s) (see, e.g., Figure 7). Once SOR T is known, there are only two independent parameters in Equation (17). For an isothermal reactor, t _max (iso) is given by:

iso) ~ (E _r exp(E _d /RT ₀ ^iS0 ))/(E _d A _d )

(18)

where E _r and E _d are respectively the activation energy of the aromatization reaction (average of those of C ₆ , C ₇ and Cg) and the deactivation reaction. T ₀ ^1S0 is the SOR temperature of the isothermal reactor.

[00111] In order to establish the relationship between these three parameters and the process variables, full calculations with the deactivation kinetic model for the process variables that span a wide range of run-lengths and SOR temperatures are performed. Equation (17) is fitted to the resulting fouling curves to determine the t _ma>; and T _a . In embodiments, t _max is correlated linearly to t _max (iso) when the latter is obtained from Equation (18) under the same process conditions as follows:

max 0.9101 t _max (iso)

(19)

This relationship is expected because, under similar process conditions, the maximum run-lengths should be about the same between the isothermal and adiabatic reactors.

[00112] Because T _a or β reflects the slope of the fouling curve, such are expected to be a complicated function of LHSV, H2/HC ratio, conversion of Q-eonvertibl.es andmax- Accordingly, a new process variable φ is defined as:

φ = B t _max (1:12/ HC) ^a _iiS ie _t (LHSV ) ^b + C X (t = 0) ^c

(20)

and a second order polynomial of φ is fit to the T _a obtained from the sample calculations. In Equation (20), T(t = 0) and t _iriax are, respectively, in units of Kelvin and day. The coefficients B and C as well as the exponents a, b and c are varied to obtain the best fit. In embodiments, the best fit may be achieved with φ given by:

φ = -70 (H2 / HC) isu ** (LHSV) ^0'434 + T(t = O) ²

(21)

and T _a may be found to correlate with φ by the equation: T _a = 6.643x 10 ² -l .O B lO ^"3 φ +1.265x 10 ^"9 φ ²

(22)

Incorporation of the conversion of the C ₆ -convertibles _, / ^' into Equation (20) leads to a zero exponent for j which leads to the belief that T _a is not a function of/, ^'

[00113] In order to determine the sensitivity of Equations (19), (21) and (22) to the correlation between the deactivation amplitude and H2/HC ratio, a slightly different equation for A _d may be examined in another analytical solution:

Ad= 1.505x 10 ¹ - 2.01 Ox 10 ³ °(H2/HC) _avera ge

(23)

in which A4 of each reactor in the series is determined by the average H2 HC ratio in that reactor. The results are summarized in the Example 3. The linear correlation between t„, _ax and t _max (iso) remains valid although the proportionality constant is slightly different. Equation (21 ) remains unchanged and the coefficients in Equation (22) are slightly modified.

[00114] The SOR inlet temperatures T(t = 0) for various process conditions are calculated directly from the deactivation kinetic model of the adiabatic reactors. Since SOR temperatures are independent of the con-elation between A _d and H2/HC ratio, SOR T only needs to be calculated once and at t = 0.

[00115] With T(t = 0), t _max and T _a determined, the fouling curve according to Equation ( 17) for the adiabatic reactors is fully defined. The run-length of the reforming catalyst as a function of conversion of Ce-convertibles LHSV, and H2/HC ratio to the 1 ^st reactor is then determined from Equation (17) by setting T equal to the selected EOR temperature. Thus, the run-length for a series of adiabatic reactors is calculated using Equations (17), (19), (21) and (22) and the values of SOR T from the deactivation kinetic model of the adi abatic reactors.

[00116] The net present pounds of benzene produced per pound of catalyst during the run-length of the catalyst (N PPB) may be calculated from the Equations (10) and (1 1), as for the isothermal reactor embodiment.

EXAMPLES

[00117] The inventive subject matter having been generally described, the following examples are given as particular embodiments of the inventive subject matter and to demonstrate the practice and advantages thereof. It is understood that the examples are given by way of illustration and are not intended to limit the specification of the claims to follow in any manner. The data in the Examples is computer generated using Equations 1 to 9 as the embodiment for the deactivation kinetic model.

EXAMPLE 1

[00118] The above-described deactivation kinetic model was used to determine process parameters which improve a catalytic reforming process utilizing an isothermal reactor. The relevant deactivation kinetic parameters, EOR condition, and feedstream composition for an isothermal reactor embodiment are summarized below in Table 1.

TABLE 1

[00119] An end-of-run (EOR) catalyst temperature of 940 °F (777 K) was used, which translates to a reactor skin temperature of about 1000 °F. E _r in Equation (12), which is the average activation energy of the aromatization reaction of the C ₆ -, C ₇ - and Cg-convertibles, has a value of 71 ,5 kcal/mole. T _Q is the start-of-run (SOR) temperature, which depends on LHSV, H2/HC ratio, and conversion of C6-convertibl.es.

[00120] Figure 4 is a plot of fouling curves predicted by the deactivation kinetic model (i.e., the shaped data points) and the fit of Equation (12) to the deactivation kinetic model (i.e., the solid lines). The deactivation kinetic model calculations are based on a feedstream composition shown in Table 1 , a conversion of Ce-convertibles of 85% and a LHSV of 1 .6. It is clear that the fits are exceptionally good. Figure 5 shows A _d as a function of the inlet H2/HC ratio for different values of conversion of Ce-convertibles for an isothermal reactor. The results in Figure 5 can be reduced to a single equation. Equation (13) hereinabove, which was used in Equation (12). In Figure 4, Equation (12) was fitted to the fouling curves of the deactivation kinetic model, i.e., the curves made by the series of data-points generated with Equations (1) to (9).

[00121] With an EOR temperature of 940 °F, LHSV of 1 ,6 and conversion of CVeonvertibles of 85%, it can be seen in Figure 4 that the run-length for the inlet H2/HC ratios of 0.5, 1, 2, 3, and 4 are about 28.8, 27.7, 27.6, 30.6, and 38.8 months, respective!)'. Plots like Figure 4 may be made for different combinations of conversion of Ce-convertibles, LHSV, and inlet H2/HC ratio values. However, instead of reading run-length off numerous plots of Equation (12) (e.g., Figure 4) for numerous combinations of conversion of Ce-eonvertibles, LHSV, and inlet H2/HC ratio values, the run-length of the reforming catalyst was calculated as a function of conversion of Cs-convertibles, LHSV, and inlet H2/HC ratio as described above using Equations (12) and (13) and the values of T _c . from the deactivation kinetic model utilizing an isothermal reactor. The results are shown in Table 2, with run-length in units of months rounded to the nearest tenth.

TABLE 2

I Catalyst Life, months

i Conversion of LH SV

[ Convertibles = 83% 1 1.2 1.6 2 2.4 2.8 3 4 J 5 ί 0.5 I 37.8 30.4 25.5 22.1 19.5 18.5 14.6 12.0

1 36.3 29.1 24.4 21.1 18.6 17.6 13.8 1 1.3

H2/HG 36.3 28.9 24.2 20.8 18.3 17.2 13.4 10.9

3 40.2 31.9 26.5 22.7 19.9 18.7 14.4 1 1.6

1 4 I 50.6 40.0 33.1 28.3 24.6 23.1 17.6 14.0

I Conversion of LH SV

I Convertibles ~ 87% I 1.2 1.6 2 2.4 2.8 3 4 I 5

I 0.5 I 35.9 27.2 22.8 19.7 17.4 16.4 12.9 10.6

1 32.8 26.2 21.9 18.9 16.6 15.7 12.2 10.0

H2/HC 2 33.1 26.3 21.9 18.8 16.4 15.5 1 1.9 9.6

3 37.1 29.3 24.3 20.7 18.1 17.0 12.9 10.3

I 4 47.7 37.5 31.0 26.3 22.8 21.3 16.1 12.6

Conversion of Cg- LH SV ]

I Convertibles ~ 89% 1.2 1.6 2 2.4 2.8 3 4 5 I

0.5 I 32.0 25.6 21.4 18.5 16.3 15.4 12.0 9.8 I

1 30.9 24.7 20.6 17.7 15.6 14.7 1 1.4 9.3

H2/HC 2 31.4 24.9 20.7 17.7 15.5 14.5 1 1.2 9.0

3 35.4 27.9 23.1 19.6 17.1 16.0 12.1 9.6

I 4 I 46.1 36.1 29.7 25.2 21.8 20.4 15.2 1 1.8 I

Conversion of LH SV

Convertibles = 91 % I 1.2 1.6 2 2.4 2.8 3 4 j 5

0.5 29.8 23.8 19.9 17.1 15.1 14.2 1 1.1 9.0 I

1 29.0 23.0 19.2 16.5 14.5 13.6 10.5 8.5

H2/HC 2 29.5 23.3 19.4 16.5 14.4 13.5 10.3 8.2

3 33.5 26.3 21.7 18.4 16.0 15.0 1 1 ,2 8.8

44.1 34.5 28.3 23.9 20.6 19.2 14.2 j 1 1.0 j

Conversion of Ce- LH SV

Convertibles ~ 93% I 1 ,2 1.6 2 2.4 2.8 3 ! 5

0.5 I 27.5 21.9 18.3 15.7 13.8 13.0 10.1 8.2

1 26.8 21.3 17.7 15.2 13.3 12.5 9.6 7.7

H2/HC 2 27.4 21.6 17.9 15.2 13.2 12.4 9.4 7.4

3 31.3 24.6 20.2 17.1 14.7 13.8 10.2 j 8.0

4 I 41.8 32.6 26.6 22.4 19.2 17.9 13.0 9.9 [00122] NPPB was calculated from Equations (10) and (1 1) using the results from Table 2. The NPPB values as a function of conversion of C6-convertibl.es, LHSV, and inlet H2/HC ratio are shown in Table 3. The units of NPPB shown in Table 3 are expressed in pounds of benzene produced per pound of catalyst.

TABLE 3

Net Present pounds of benzene 7 pounds catalyst-catalytic life j

Conversion of C _e - LHSV i Convertsbies ~ 83% 1.2 1.6 j 2 j 2.4 2.8 3 4 f 5 I

I 0.5 10,177 11,076 12,203 I 12,602 13,480 13,110 14,756 I 14,947 j

1 9,722 10,752 11,360 12,081 12,861 13,110 13,831 13,761

H2/HC 2 9,722 10,752 11,360 j 12,081 12,236 12,435 12,899 j 13,761

3 10,624 11,716 12,619 13,119 13,480 13,780 13,831 14,947

4 12,949 14,166 15,039 15,637 16,497 16,399 17,481 17,289 |

Conversion of C ₆ - LHSV [ Convertibles = 85% 1.2 1.6 I 2 j 2.4 2.8 3 4 j 5 j

I 0.5 I 9,956 11,011 11,634 12,372 12,531 12,735 14,164 14,092 I

I 1 I 9,720 10,676 11,196 j 11,833 12,531 12,735 13,209 j 14,092

H2/HC 2 9,720 10,676 11,196 11,833 11 ,886 12,037 13,209 12,867

3 10,652 11,672 12,067 12,906 13,171 13,426 14,164 14,092

I 4 112,845 14,203 j 14,998 | 15,508 16,287 16,133 16,980 j 16,512 |

J Conversion of C ₆ - LHSV i I Convertibles ~ 87% 1.2 1.6 1 2 2.4 2.8 3 4 j 5 ]

I 0.5 I 9,956 10,582 j 11,460 12,111 12,165 12,320 13,520 14,424 I

1 9,459 10,233 11,008 11,555 12,165 12,320 12,534 13,170

H2/HC 2 9,459 10,233 11,008 11,555 11.499 11,600 12,534 13,170

3 10,430 11,270 11,907 12,663 12,826 13,034 13,520 13,170

I 4 I 12,931 14,223 14,933 15,349 16,044 15,829 16,427 16,900 |

J Conversion of C _e - LHSV I I Convertibles = 89% 1.2 1.6 I 2 j 2.4 2.8 3 4 j 5 j

I 0.5 I 9,422 10,468 I 10,795 11,246 11,763 11,867 12,822 13,472 I

1 9,166 10,108 10,795 11,246 11.763 11,867 11,805 12,178

H2/HC 2 I 9,166 10,108 j 10,795 11,246 11,076 11,867 11,805 12,178

3 10,177 11,179 11,723 12,390 12,445 12,604 12,822 13,472

1 4 112,782 13,899 14,846 15,162 15,766 15,487 15,822 16,028 |

Conversion of C _e - LHSV

Co!wertibies = 91% 1.2 1.6 I 2 j 2.4 2.8 3 4 j 5

0.5 I 9,108 9,964 10,557 10,907 11,324 11,373 12,070 12,451

1 8,841 9,589 10,072 10,309 10,615 11,373 12,070 12,451

H2/HC 2 9,108 9,589 10,072 10,907 10,615 11 ,373 11,020 11,116

3 9,893 10,703 11,514 11,499 12,028 12,133 12,070 12,451

4 I 12,604 13,874 14,287 14,946 15,453 15,108 15,164 15,087

Conversion of C _e - LHSV i Convertibies 93% 1.2 1.6 j 2 j 2.4 2.8 3 4 f 5 j

0.5 I 8,761 9,414 I 9,793 I 10,536 10,848 10,839 11,262 11,360 I

1 8,484 9,024 9,793 9,920 10,117 10,049 11,262 11,360

H2/HC 2 I 8,484 9,414 9,793 j 9,920 10,117 10,049 10,180 I 9,984

3 9,578 10,562 10,789 11,147 11 ,573 11,623 11,262 11,360

4 I 12,397 13,481 14,139 14,121 14,411 14,690 14,453 14,078 | [00123] From Tables 2 and 3, values for LHSV and inlet H2/HC ratio which improve the process are about 4, for a desired run-l ength of the reforming catalyst of at least 12 months and for a conversion of Ce-convertibles between about 83% and 93%. It can be seen there is little advantage to using a high H2/HC ratio unless it is greater than about 2. Higher values of NPPB are obtained when the catalyst is operated at lower conversion of C6-coiivertibles and high LHSV. For a given conversion of iVconvertibles and LHSV, higher NPPB is obtained with a H2/HC ratio greater than 2 or below 1.

[00124] Figure 6 shows an example of how k _d (T ₀ ) varies with the H2/H C ratio for isothermal reactor, as given by Equation (9). The curve in Figure 6 is for an LHSV of 1.6 and a conversion of Ce-convertibles of 85%. Because the deactivation kinetic model used in this Example 1 employs a linear correlation between A _f j and H2/HC ratio, the NPPB in Table 3 are sensitive to the uncertainties in the correlation. A change in the correlation can affect the shape of the hump in k _d (T ₀ ) shown in Figure 6 or the minima in the run-length in Table 2, making them more pronounced or shallower. The trends described above have been found to remain valid for different correlations.

EXAMPLE 2

[00125] The above-described deactivation kinetic model was used to determine process parameters which improve a catalytic reforming process comprising a series of adiabatic reactors. The process has a six reactor configuration with catalyst volume ratios of 10: 10:10:20:20:30. The relevant kinetic parameters are listed in Table 4.

TABLE 4

[00126] Figure 7 is a plot of fouling curves (i.e., the shaped data points), predicted by the deactivation kinetic model and fits of Equation (17) to the model (i.e., the solid lines). The four curves in Figure 7 are based on model calculations with four different combinations of values for conversion of CVeonvertibles, LHS V, and H2/HC ratio. Particularly, one of the curves represents a LHSV of 1.6, an H2/HC ratio of 3, and a conversion of CVeonvertibles of 0.87. Another of the curves represents a LHSV of 2, an H2/HC ratio of 2, and a conversion of Ce-convertibles of 0.85. Another of the fouling curves represents a LHSV of 1 .2, an H2/HC ratio of 4, and a conversion of Ce-convertibles of 0.83. Another of the curves represents a LHSV of 4, an H2/HC ratio of 4, and a conversion of Ce-convertibles of 0.9. The fits of the data generated by the deactivation kinetic model fit the data quite well, even for series of adiabatic reactors]

[00127] In Figure 7, with an EOR temperature of 540 °C, it can be seen that when the LHSV of 1.6, an H2/HC ratio of 3, and a conversion of Ce-convertibles of 0.87 is used, the run-length is about 26.6 months. When a LHSV of 2, an H2/HC ratio of 2, and a conversion of ^convertibles of 0.85 is used, the run-length is about 20.6 months. When a LHSV of 1.2, an H2/HC ratio of 4, and a conversion of Ce- converiibies of 0.83 is used, the run-length is about 46.5 months. When a LHSV of 4, an H2/HC ratio of 4, and a conversion of (^-convertibles of 0.91 is used, the run- length is about 12.3 months.

[00128] Plots like Figure 7 may be made for different combinations of conversion of Cs-convertibles, LHSV, and H2/HC ratio values. However, instead of reading run- length off numerous plots of Equation (17) (e.g., Figure 7) for numerous combinations of conversion of Ce-convertibles, LHSV, and H2/HC ratio values, the run-length of the reforming catalyst was calculated as a function of conversion of Ce-convertibles, LHS V, and H2/HC ratio as described above using Equations (16) to (19), (21) and (22). Catalyst fouling is accounted for by Equation (16). The results are shown in Table 5, with run-length in units of months rounded to the nearest tenth.

Conversion of C«- LH SV

Convertibles - 87% 1.2 ! 1.6 2 2.4 2.8 3 4 5

0.5 30.2 ! 22.2 18.2 15.4 13.2 12.3

1 28.0 I 22.1 18.2 15.4 13.3 12.5 3.1 0 8

1st reactor

2

inlet H2/HC 29.6 j 23.4 19.4 16.5 14.3 13.4 ^' ΊοΊ ? ^' ε

3 33.7 j 26.6 22.1 18.8 16.4 15.4 1 1 ft ¾ 1

4 43.0 ! 33.7 27.7 23.5 20.4 19.1 14.4 I ; 3

[00129] NPPB was calculated from Equations (10) and (I I) using the results from Table 5. The NPPB values as a function of conversion of Ce-convertibles, LHSV, and H2/HC ratio are shown in Table 6. The units of NPPB shown in Table 6 are expressed in pounds of benzene produced per pound of catalyst.

TABLE 6: Net present pounds of benzene /pound catalvst-catalvtic life

A _d = 1231x10 ¹¹ ~ 1507x10 ⁰ (H2/HC)iniei ^"

Conversion of C ₆ - LHSV

Convertifo!es ~ 83% 1.2 j 1.6 j 2 2.4 2.8 I 3 I 4 5

[ 0.5 8,687 j 9,319 j 9,953 10,369 10,845 j 10,9 40 ji; 11327

1st 1 8,687 † 9,319 T 9,953 10,369 10,845 j 10,9 40 11327 reactor

2 8,921 9,651 10,382 10,898 10,845 j 10,940 p : i .¾ : : i 227 inlet

H2/HC 3 9,837 10,629 11,230 11 ,943 12,097 12,293 12,752 : 0

4 11,784 j 13,116 [ 13,687 14,477 15,131 I 14,929 [ 15,493 15,939

Conversion of C ₆ - LHSV

ConvertifoSes = 85% 1.2 I 1.6 I 2 2.4 2.8 j 3 j 4 5

0.5 8,410 I 9,201 I 9,748 10,072 9,802 10,502 | 1 .17 102

1st 1 8,410 9,201 9,748 10,072 9,802 10,502 | 10.17 102 reactor

2 8,654 9,544 10,192 10,619 10,457 10,502

inlet i 1 ! ! .if.fi H2/HC 3 9,609 10,554 11,069 11 ,698 11,750 I 11,900 I 12,106 i .. ^v 2 ^< >

4 11,854 j 12,812 [ 13,607 14,316 14,885 j 14,623 j 14,939 15,133

Conversion of C ₆ - LHSV

Convertib!es = 91% 1.2 j 1.6 j 2 2.4 2.8 I 3 I 4 I 5

0.5 I 7,388 j 7,966 j 7,997 8,388 f. ;;¾ i S i 71 i ti 00 0 l

1st 1 I 7,388 T 7,966 T 7,997 8,388 f. ;;¾ i S i 71 j ? 7K- 0 l reactor

2 7,936 † 8,349 73

iniet I T 8,985 8,995 9,0 I s ^« wg i S 7 1 S3M H2/HC 3 9,004 9,850 9,957 10,192 10,494 j 10,485 p vGC

4 j 11,274 j 12,005 [ 12,772 13,094 13,263 j 13,478 ] 12,961 1 309

EXAMPLE 3

[00130] The sensitivity of the results to the Ad vs. H2/HC ratio correlation was estimated by using another analytical soluation for Ad, Equation (23). Run-lengths and NPPB were again calculated for the same values of conversion of Co-convertibles, LHSV, and H2/HC ratio, except that Equation (23) was used for Ad instead of Equation (16). Tables 7 and 8 display the results.

TABLE 71 Catalyst Life, months Aa= 1505x10 ¹¹ -2,010x10 ¹⁰ (H2/HC)average

TABLE 8: Net Present pounds of benzene /pounds catalyst-catalytic life

A _d = 1.505x1 G ¹ -2,010x10 ¹⁰ (H2/HC)average

Co r! vers so r! of C ₆ - υ ^■ 3SV

Cos ertibies = 85% 1.2 1.6 2 2.4 ! 2.8 ! 3 j 4 ! 5

0.5 7,413 8,154 8,393 8,963 9,141 9,080 1 3 ,:;

1 7,413 8,154 8,393 8,963 9,141 9,080 9.193

1st reactor

2 7,916 8,506 9,301 9,520 9,802 9,794 10

inlet H2/HC

3 9,136 9,884 10,632 11,161 11,107 11,204 12100 1 ; -·:ϋε

4 11,854 13,123 13,607 14,316 14,885 14,623 14,939 15,133

Conversion ί >f Cs- LH SV

I Convertibles = 87% 1.2 1.6 2 2.4 2.8 3 4 5 j 0.5 7,413 7,616 8,120 8,599 8,674 8 8.105 7.9-ίίί

1 7,326 7,616 8,120 8,599 8,674 ¾ G¾o y.-iU- ^' 7.- H

1st reactor

2 7,587 8,346 8,591 9,174 9,356 9,293 9,-M ">

3 8,858 9,768 10,432 10,868 10,703 10,749 11 -10* 1 70*

I 4 11,912 12,792 13,504 14,126 14,605 14,279 15,290 1420O

Con vers son ί rf C ₆ - LH SV

Convertibles = 89% 1.2 ! 1.6 2 2.4 2.8 3 4 5

0.5 6,954 7,412 7,821 8,204 £ ;C9 99 ; .£ ^"■ ¾ υ.ί-ϋό

1 6,954 7,412 7,821 8,204 ;C9 / 99 ; .£ ^"■ ¾ υ.ί-ϋό

! SSs. sw GlO

2 7,495 8,166 8,307 8,797 8,873 ? !- · 3.151 ssi!eiH2/HC

3 8,806 9,634 9,738 10,546 10,263 10,996 io κ>ζ 1074*

4 11,728 12,754 13,377 13,907 4,290 13,898 14,662 14 3.<

[00131] Figure 8 shows the linear correlation of t _max with t _max (iso) when the latter is obtained from Equation (18) under the same process conditions and when t,. _iax is determined by fitting Equation (17) to fouling curves from the full calculations with the deactivation kinetic model for a series of adiabatic reactors.

[00132] Equation (22) was plotted for the data in Example 2, and the results are shown in Figure 9 with an R ^z value of 0.9551. φ was given by Equation (21). Figure 9 shows T _a as a function of φ obtained when a second order polynomial represented by Equation (22), i.e., second order with respect to φ, is fitted to values of T _a determined by fitting Equation (17) to the full calculations with the deactivation kinetic model.

[00133] To compare the results of Example 2 in Figures 8 and 9 so as to determine sensitivity, full calculations with the deactivation kinetic model using the correlation in Equation (23) instead of Equation (16) were performed. Fitting Equation (17) to the results provides the following equation:

= 0.9591 t _max (iso) (24) a plot of which is sho wn in Figure 10 (R~ value is 0.9796), and the equation:

T _a = 6.959x 10 ² -l . lOlx 10 ^"3 φ +1.323x 10 ^~ V (25) a plot of which is shown in Figure 1 1 (R value is 0.9893). φ was again given by Equation (21 ). Figure 50 shows t _max as a function of t _max (iso) found by repeating full calculations with the deactivation kinetic model using the correlation in Equation (23) and fitting Equation (17) to the results. Figure 11 shows T _a as a function of φ determined by repeating full calculations with the deactivation kinetic model using the correlation in Equation (23) and fitting Equation (17) to the results. The higher R" values for Figures 10 and 1 1 indicate a better fit for the curves in Figures 10 and 1 1 than the curves in Figures 8 and 9.

[00134] As can be seen, Table 5 and Table 7 have different yalues as do Table 6 and Table 8. The different values between Tables 5 and 7 and between Tables 6 and 8 may serve as indicators of the magnitude of the uncertainty in the values given in the tables. The average difference in the benzene yields between Table 6 and 8 is about 7%, with a standard deviation of about 5%, This provides an estimate of the uncertainty and sensitivity in the numbers in these tables. [00135] The results in Tables 5-8 indicate more benzene may be produced per catalyst life cycle (e.g., run-length) when the catalyst is operated at lower conversion of C _o -convertibles. There is no significant advantage to using a high H2/HC ratio unless it is greater than about 2, At low to moderate conversions of C6-convertibles, there are maxima in benzene yield as a. function of LHSV. For a run-length greater than 12 months, LHSV values are from about 3 to about 4, For a given conversion of Ce-convertibles and LHSV, improved NPPB may be attained with H2/HC ratio of less than 1 or greater than 2.

[00136] Several trends and correlations have been observed. Specifically, more benzene can be produced in a catalyst, life cycle (e.g., run-length) when the catalyst is operated under conditions of low conversion of Ce-convertibles and high LHSV. Furthermore, more benzene can be produced in a catalyst life cycle when the catalyst is operated at a H2/HC ratio larger than about, 3. Some of the trends that were observed in the results of the isothermal reactor embodiment of Example I can also be found for the adiabatic reactors embodiment, as shown in Example 2. First, for a given LHSV and 1 12 I K ratio to the 1 ^st reactor, a lower conversion of C6-convertibles gives a larger NPPB. Second, a higher 1 12 I K ratio does not necessary lead to a longer run-length unless the ratio is 3 or larger. There are little differences in the run-lengths and the benzene yields in terms of NPPB when the H2/HC ratio is between 0.5 and 2, An increase in the H2/HC ratio to the 1 ^st reactor increases the heat capacity of the gas in all the reactors. For a given conversion of the feedstream, a higher H2/HC ratio will result in smaller endotherms in the reactors, which translates to higher average catalyst bed temperatures. The higher average bed temperature means a larger Arrhenius term in the deactivation rate constant. Because the deactivation rate constant is a product of A _d and the Arrhenius term, the higher H2/HC ratio has opposite effects on the deactivation rate, namely reducing the deactivation amplitude but increasing the Arrhenius factor. These opposing effects result in small changes in the overall deactivation rate of the catalyst when the H2/HC ratio to the 1 ^st reactor is between 0,5 and 2. At a H2/HC ratio of 3, the reduction in A _d is significantly large that it outweighs the increase in the Arrhenius factor. It is interesting to note that although an increase in the H2/HC ratio increases the average bed temperature, it actually lowers the reactor temperature. Thus, embodiments of the inventive subject matter may comprise selecting an LHSV greater than 3: alternatively, selecting an LHSV of about 4.

[00137] At low to moderate conversions of C ₆ -convertibles, there are maxima in the benzene yield as a function of the LHSV. In embodiments, there are maxima in the NPPB at a. LHSV of around 3, This may be particularly true when conversion is high. When designing a plant, it is important to match the LHSV with the H2 HC ratio and the required conversion of Q-convertibles in order to achieve the maximum benzene yield from the catalyst. Because it is difficult to have a catalyst replacement every year logistically and economically, run-lengths of less than one year are undesirable. In embodiments, for low to moderate conversions (83% to 89%), the LHSV may be between 2.8 and 3.

[00138] It should be noted that the occurrences of maxima are expected and are sensitive to the definition of the EOR inlet temperature. When EOR T is set to a high value, the operation window between the EOR T and SOR T can become so small for a high LHSV operation that the run-length of the reforming catalyst will drop precipitously.

[00139] Upon reviewing the data for Example 1 (isothermal reactor) and Example 2 (adiabatic reactors), it may be noted that capita! savings may be obtained for reforming processes using a reforming catalyst comprising large-pore zeolite support with a platinum containing compound and at least one ha!ide by going from an adiabatic design to an isothermal one.

ADDITIONAL DESCRIPTION

[00140] The following are nonlimiting, specific embodiments in accordance with the present disclosure.

[00141] Embodiment 1 is a method of catalytic reforming comprising feeding a feedstream comprising Ce-convertibles to one or more reactors: contacting the feedstream with a reforming catalyst selecting values for a LHSV, a H2/HC ratio, and a conversion of Cg-convertibles from a deactivation kinetic model so as to maximize a net present amount of benzene produced over a run-length of the reforming catalyst; operating the one or more reactors at the selected LHSV, the selected H2/HC ratio, and the selected conversion of Ce~ convertibles; and recovering an effluent from the reactor, wherein the effluent comprises at least about 40 wt % benzene.

[001 2] Embodiment 2 is the method of embodiment 1 wherein the selected H2/HC ratio is greater than 2 or less than 1.

[00143] Embodiment 3 is the method of any of embodiments 1 through 2 wherein the selected LHSV is from about 3 to about 4 and the selected H2/HC ratio is about 4.

[00144] Embodiment 4 is the method of any of embodiments 1 to 3 wherein the reforming catalyst comprises a silica-bound L-type a zeolite support, platinum, chloride, and fluoride.

[00145] Embodiment 5 is the method of any of embodiments 1 to 4 wherein selecting values comprises calculating a fouling curve and determining the run-length of the reforming catalyst at a time at which the fouling curve reaches a prescribed end of ran temperature.

[00146] Embodiment 6 is the method of embodiment 5 wherein calculating the fouling curve comprises numerically solving a plurality of equations representing three dimensional process variable space spanned by the conversion of CVconvertibles in the feedstream, the LHSV, and the H2/HC ratio.

[00147] Embodiment 7 is the method of any of embodiments 5 to 6 wherein selecting values further comprises approximating the fouling curve with an analytical solution to the deactivation kinetic model.

[00148] Embodiment 8 is the method of embodiment 7 wherein the analytical solution is represented by: r ) ·· ·.· ^ 1

1 4- (RT ₀ I E _d jlnjl - (E _d I E _r )k _d (T ₀ )t

where R is the gas constant; E _d is an activation energy of a deactivation reaction; E _r is an average activation energy of a conversion of C ₆ -, C ₇ - and Cg-eonvertibles; T ₀ is a start-of- run temperature: and k _d (To) is represented by A _d exp(-E _d / RT) with A _d being the average of i along at least one isothermal reactor and further represented by:

A _d =1.343x 10" -5.282X 10 ^{ϊ 0} /-1 ·456χ lG ^{! 0} (H2/HC) _¾/ei

where/is the fraction representing conversion of Ce-convertibles.

[00149] Embodiment 9 is the method of any of embodiments 7 to 8 wherein the analytical solution is represented by:

T T _i( p O ) ln(l-t/t _raax )

where T _a is the inlet temperature of at least one adiahatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst, near start-of-run, wherein β is adjusted to fit the analytical solution to the fouling cun'e, [00150] Embodiment 10 is the method of any of embodiments 1 to 9 further comprising calculating a net present amount, of benzene produced and a benzene yield per month of the reforming catalyst to calculate an amount of benzene produced over the run-length per unit of the reforming catalyst and then discounting the amount of benzene produced over the run- length of the reforming catalyst per unit of reforming catalyst to determine the net present amount of benzene produced per unit of reforming catalyst over the run-length of the reforming catalyst.

[00151] Embodiment 1 1 is the method of any of embodiments 1 to 10 further comprising calculating a net present amount of benzene comprises calculating a number of pounds of benzene produced per pound of the reforming catalyst per hour, N, which is given by:

N = 0.907/ wee, LHSV (¾ ^■ / d,)

where dif, d _c , f, and w _¾ are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of (Vconvertihles, and the weight fraction of C ₆ -convertibies in the feedstream. [00152] Embodiment 12 is the method of embodiment 1 1 wherein N is greater than about I ,

[00153] Embodiment 13 is the method of any of embodiments 1 to 12 further comprising calculating a net present amount of benzene comprises calculating a net present pounds of benzene produced per pound of the reforming catalyst over the run-length of the reforming catalyst, NPPB, which is given by:

NPPB = (N/720) { 1 - { l -Dj' l /D _m

where D,„ is the monthly discount rate and n is the run-length of the reforming catalyst in units of months,

[00154] Embodiment 14 is the method of embodiment wherein the NPPB is greater than about 14,000,

[00155] Embodiment 15 is the method of any of embodiments 11 to 14 wherein N is greater than about 1 or the NPPB is greater than about 14,000.

[00156] Embodiment 16 is a method of catalytic reforming, and optionally the method of any of embodiments 1 to 15, comprising kmetically modeling a deactivation of a reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of (^-convertibles in the feedstream ; calculating a net present amount of benzene that is produced over the run-length of the reforming catalyst; and selecting values for the LHSV, the H2/HC ratio, and the conversion of Cg-convertibles that maximize the net present amount of benzene produced over the run-length of the reforming catalyst,

[00157] Embodiment 17 is the method of any of embodiments 1 to 16 wherein the deactivation kinetic model comprises at least one aromatization reaction rate equation, at least one catalyst deactivation equation, and at least one heat balance equation.

[00158] Embodiment 18 is the method of embodiment 17 wherein the reforming catalyst comprises a Group VIII metal and at least one halide on an inorganic support and wherein the catalyst deactivation equation is a first order rate equation accounting for the deactivation of the reforming catalyst via coking,

[00159] Embodiment 19 is the method of any of embodiments 17 to 18 wherein a deactivation rate constant (A¾ (T) ) for the first order catalyst deactivation equation is: k _d (T) = A _d exp(-E _rf / RT)

where T is a. temperature of the reforming catalyst, J¾ is an activation energy of one or more catalyst deactivation reactions, and R is a gas constant.

[00160] Embodiment 20 is the method of any of embodiments 16 to 18 wherein kineticaily modeling the deactivation of the reforming catalyst with the deactivation kinetic model comprises calculating a fouling curve and determining the run-length of the reforming catalyst, at a time at which the fouling curve reaches a prescribed end of run temperature.

[00161] Embodiment 21 is the method of any of embodiments 5 and 20 wherein calculating the fouling curve comprises numerically solving a plurality of equations representing three dimensional process variable space spanned by the conversion of (^-convertibles in the feedstream, the LHSV, and the H2/HC ratio.

[00162] Embodiment 22 is the method of any of embodiments 16 to 21 wherein kineticaily modeling the deactivation of the reforming catalyst with the deactivation kinetic model further comprises approximating the fouling curve with an analytical solution to the deactivation kinetic model.

[00163] Embodiment 23 is the method of any of embodiments 7 and 22 wherein the analytical solution is represented b :

where R is the gas constant; E _f j is an activation energy of a deactivation reaction; E _f is an average activation energy of a conversion of C ₍ ,~, r and C ₈ -convertibies; T ₀ is a start-of- run temperature; and k _d (To) is represented by A _d exp(-E _< < / RT) with A _d being the average of A _d along at least one isothermal reactor and further represented by:

Ad 1 10 ⁿ -5.282X 10 ⁱ⁰ -1 .456x 10 ¹⁰ (H2/HC) _/ „ _fe/

where /is the fraction representing conversion of C6-convertibles.

[00164| Embodiment 24 is the method of any of embodiments 7 to 8 and 22 to 23 wherein the analytical solution is represented by:

where T _a is the inlet temperature of at least one adiabatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst near start-of-run, wherein β is adjusted to fit the analytical solution to the fouling curve.

[00165] Embodiment 25 the method of any of embodiments 1 to 24 further comprising determining a sensitivity of the analytical solution by calculating another analytical solution where Ad is represented by

Ad= 1.505x 10 ⁿ - 2.01 Ox 10 ^iC' (H2/HC) _aV e _rage where (H2/HC) _aV erage is the average H2/HC ratio in the at least one reactor.

[00166] Embodiment 26 is the method of any of embodiments 1 to 25 further comprising determining a sensitivity of the analytical solution by calculating another analytical solution where Ad is represented by

where (H2/!:iC)j _n j _et is the H2/HC ratio in the inlet of at least one reactor.

[00167] Embodiment 27 is the method of any of embodiments 10 to 26 wherein calculating the net present amount of benzene produced comprises using the run-length of the reforming catalyst and a benzene yield per month of the reforming catalyst to calculate an amount of benzene produced over the run-length per unit of the reforming catalyst and then discounting the amount of benzene produced over the run-length of the reforming catalyst per unit of reforming catalyst to determine the net present amount of benzene produced per unit of reforming catalyst o ver the run-length of the reforming catalyst.

[00168] Embodiment 28 is the method of any of embodiments 10 to 27 wherein calculating a net present amount of benzene comprises calculating a number of pounds of benzene produced per pound of the reforming catalyst per hour, I , which is given by:

N = 0.907/ wee, LHSV (¾ ^■ / d,)

where dif, dc, f, and w _¾ are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of (Vconvertihles, and the weight fraction of C ₆ -convertibies in the feedstream.

[00169] Embodiment 29 is the method of any of embodiments 1 1 to 28 wherein N is greater than about 1.

[00170] Embodiment 30 is the method of any of embodiments 10 to 29 wherein calculating a net present amount of benzene comprises calculating a net present pounds of benzene produ ced per pound of the reforming catalyst over the run-iength of the reforming catalyst, NPPB, which is given by:

PPB = (N/720) { 1 - (1-D _m )"} /D _m

where D _m is the monthly discount rate and n is the run-length of the reforming catalyst in units of months,

[00171] Embodiment 31 is the method of any of embodiments 13 to 30 wherein the NPPB is greater than about 14,000.

[00172] Embodiment 32 is the method of any of embodiments 28 to 31 wherein the NPPB is greater than about 14,000 or N is greater than about i .

[00173] Embodiment 33 is the method of any of embodiments 1 to 32 wherein, for a given conversion of the Ce-convertibles in the teedstream and the LHSV, a higher value for the net present amount of benzene produced is attained by selecting a H2/HC ratio of the feedstream of greater than about 2 or less than about i ,

[00174] Embodiment 34 is the method of any of embodiments 1 to 33 wherein, for a given conversion of (^-convertibles in the feedstream, the LHSV selected is about 4.

[00175] Embodiment 35 is the method of any of embodiments 1 to 34 wherein, for a given conversion of Cg-convertibles in the feedstream, the H2/HC ratio of the feedstream is about 4.

[00176] Embodiment 36 is the method of any of embodiments I to 35 wherein the conversion of C6-convertibles in the feedstream is from about 83% to about 93%, wherein the LHSV is about 4 and the H2/HC ratio of the feedstream is about 4.

[00177] Embodiment 37 is a catalytic reactor, which may optionally incorporate any of the methods of embodiments 1 to 36, comprising a feedstream and a reforming catalyst wherein the reactor is configured to produce a net present amount of benzene over a run- length of the reforming catalyst; wherein the run-length is determined as a function of a LHSV of the feedstream, a H2 HC ratio of the feedstream, and a conversion of CVconvertibles in the feedstream by kinetically modeling a deactivation of the reforming catalyst with a deactivation kinetic model: and wherein values for the LHSV of the feedstream, the H2/HC ratio of the feedstream, and the conversion of Q-convertibles in the feedstream are selected to maximize the net present amount of benzene produced over the m -length of the reforming catalyst. [00178] Embodiment 38 is the catalytic reactor of embodiment 37 wherein the deactivation kinetic model comprises at least one aromatizatioii reaction rate equatio , at least one catalyst deactivation equation, and at least one heat balance equation,

[00179] Embodiment 39 is the catalytic reactor of any of embodiments 37 to 38 wherein the reforming catalyst comprises a Group VIII metal and at least one halide on an inorganic support and wherein the catalyst deactivation equation comprises a first, order rate equation to account for the deactivation of the reforming catalyst via coking.

[00180] Embodiment 40 is the catalytic reactor of any of embodiments 37 to 39 wherein a deactivation rate constant {¾ (T)) for the catalyst, deactivation equation is:

k _d (T) = A _d exp(-E _d / RT)

where T is a temperature of the reforming catalyst, ¾ is an activation energy of one or more catalyst deactivation reactions, and R is a gas constant.

[00181] Embodiment 41 is the catalytic reactor of any of embodiments 37 to 40 wherein the deactivation of the reforming catalyst is kinetically modeled with the deactivation kinetic model by calculating a fouling curve and determining the run-length of the reforming catalyst at a time at which the fouling curve reaches a prescribed end of run temperature.

[00182] Embodiment 42 is the catalytic reactor of embodiment 45 wherein the fouling curve is calculated by numerically solving a plurality of equations representing three dimensional process variable space spanned by the conversion of Cg-convertibles in the feedstream, the LHSV of the feedstream, and the H2/HC ratio of the feedstream.

[00183] Embodiment 43 is the catalytic reactor of any of embodiments 37 to 42 wherein the deactivation of the reforming catalyst is kinetically modeled with the deactivation kinetic model by approximating the fouling curve with an analytical solution to the deactivation kinetic model.

[00184] Embodiment 44 is the catalytic reactor of embodiment 43 wherein the analytical solution is represented by:

where R is the gas constant; E _d is an activation energy of a deactivation reaction; E _r is an average activation energy of a conversion of C ₆ -, C ₇ - and Cg-eonvertibles; T ₀ is a start-of- ran temperature; and k _d (To) is represented by A _d exp(-E _< ; / RT) with A _d being the average of A,j along at least one isothermal reactor and further represented by:

A _d ==1.343x 10 ^{5 !} -5.282X l ⁽⁾ⁱ⁰ /-1.456x 10 ¹⁰ (H2/HC) _¾fe ,

where /is the fraction representing conversion of C ₆ -convertibles.

[00185] Embodiment 45 is the catalytic reactor of any of embodiments 43 to 44 wherein the analytical solution is represented by:

T = T, -p/(1 ) ln(l--t/t _mira )

where T _a is the inlet temperature of at least one adiabatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst near start-of-run, wherein β is adjusted to fit the analytical solution to the fouling curve.

[00186] Embodiment 46 is the catalytic reactor of any of embodiments 37 to 45 wherein the net present amount of benzene produced over the run-length of the reforming catalyst is calculated per unit of the reforming catalyst by using the run-length of the reforming catalyst and a benzene yield per month of the reforming catalyst and then discounting the amount, of benzene produced over the run-length of the reforming catalyst, per unit of reforming catalyst.

[00187] Embodiment 47 is the catalytic reactor of any of embodiments 37 to 46 wherein the net present amount of benzene produced over the run-length of the reforming catalyst is calculated as a number of pounds of benzene produced per pound of the reforming catalyst per hour, N, which is given by:

N = 0.907/ wee HSV (¾ _f / d,}

where ά _¾ d _c , and wee are, respecti vely, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of C^-convertibles, and the weight fraction of (^-convertibles in the feedstream.

[00188] Embodiment 48 is the catalytic reactor of embodiment 47 wherein N is greater than about 1.

[00189] Embodiment 49 is the catalytic reactor of any of embodiments 37 to 48 wherein the net present amount of benzene produced over the run-length of the reforming catalyst is calculated as net present pounds of benzene produced per pound of the reforming catalyst over the run-length of the reforming catalyst, NPPB, which is given by: NPPB = (N/720) { 1 - (1-D _m ) ⁿ } /D _m

where D _m is the monthly discount rate and n is the run-length of the reforming catalyst in units of months,

[00190] Embodiment 50 is the catalytic reactor of embodiment wherein the NPPB is greater than about 14,000.

[00191] Embodiment 51 is the catalytic reactor of any of embodiments 47 to 50 where the NPPB is greater than about 14,000 or N is greater than about I .

[00192] Embodiment 52 is the catalytic reactor of an of embodiments 37 to 51 wherein, for a given conversion of the Ce-convertibles in the feedstream and the LHSV, a higher value for the net present amount of benzene produced is attained by selecting a H2/HC ratio of the feedstream of greater than about 2 or less than about 1 .

[00193] Embodiment 53 is the catalytic reactor of any of embodiments 37 to 52 wherein, for a given conversion of Ce-convertibles in the feedstream, the LHSV selected is about 4.

[00194] Embodiment 54 is the catalytic reactor of any of embodiments 37 to 53 wherein, for a given conversion of (^-convertibles in the feedstream, the H2 HC ratio of the feedstream is about 4.

[00 95] Embodiment 55 is the catalytic reactor of any of embodiments 37 to 54 wherein the conversion of (^-convertibles in the feedstream is from about 83% to about 93%, wherein the LHSV is about 4 and the H2/HC ratio of the feedstream is about 4.

[00196] Embodiment 56 is the catalytic reactor of any of embodiments 1 to 55 wherein the reforming catalyst comprises a large pore zeolite, platinum, and optionally one or more halides.

[00197] Embodiment 57 is the catalytic reactor of any of embodiments 1 to 57 wherein the reforming catalyst comprises a silica-bound KL-zeolite support, platinum, chloride, and fluoride.

[00198] Embodiment 58 is the catalytic reactor of any of embodiments 37 to 57, further comprising an effluent having at least about 40 wt % benzene.

[00199] Embodiment 59 is a computer system, which may optionally be used in combination with any of embodiments 1 to 58, comprising one or more processors; a. non- transitory memory coupled to the processor, wherein the non-transitory memory comprises executable instructions that, when executed, cause the one or more processors to: kineticallv model a deactivation of a reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstream, a H2/HC ratio of the feedstream, and a conversion of Ce-convertibles in the feedstream; calculate a net present amount of benzene that is produced over a run-length of the reforming catalyst; and select values for the LHSV, the H2/HC ratio, and the conversion of Cg-convertibles that maximize the net present amount of benzene produced over a run- length of the reforming catalyst.

[00200] Embodiment 60 is the computer system of embodiment 59 wherein the deactivation kinetic model comprises at least one aromatization reaction rate equation, at, least one catalyst deactivation equation, and at least one heat, balance equation,

[00201 ] Embodiment 61 is the computer system of embodiments 59 to 60 wherein the reforming catalyst comprises a Group VIII metal and at least one halide on an inorganic support and wherein the catalyst deactivation equation is a first order rate equation.

[00202] Embodiment 62 is the computer system of embodiments 59 to 61 wherein a deactivation rate constant (k _d (TJ) for the first order catalyst deactivation equation is:

k _d (T) === A _d expi-Er _f / RT)

where T is a temperature of the reforming catalyst, E _d is an activation energy of one or more catalyst deactivation reactions, and R is a gas constant.

[00203] Embodiment 63 is the computer system of any of embodiments 59 to 62 wherein kinetically model the deactivation of the reforming catalyst with the deactivation kinetic model comprises calculate a fouling curve and determine the run-length of the reforming catalyst at a time at which the fouling curve reaches a prescribed end of ran temperature.

[00204] Embodiment 64 is the computer system of embodiment 63 wherein calculate the fouling curve comprises numerically solving a plurality of equations representing three dimensional process variable space spanned by the conversion of CVconvertibles in the feedstream, the LHSV, and the H2/HC ratio.

[00205] Embodiment 65 is the computer system of any of embodiments 59 to 64 wherein kinetically model the deactivation of the reforming catalyst with the deactivation kinetic model further comprises approximate the fouling curve with an analytical solution to the deactivation kinetic model.

[00206] Embodiment 66 is the computer system of embodiment 65 wherein the analytical solution is represented by:

7 i ^/ -

1 + (RT ₆ I E _d )ln[l -- ^} -——— 1

(E _d I E _r )k _d (Γ ₀ )t

where R is the gas constant; Ed is an activation energy of a deactivation reaction; E _r is an average activation energy of a conversion of C ₆ -, C ₇ - and (^-convertibles; T ₀ is a start-of- run temperature; and k _d (To) is represented by A., exp(-E, _f / RT) with A _d being the average of along at least one isothermal reactor and further represented by:

A _d =1.343x 10 ^{1 1} -5.282X 10 ¹⁰ /-1.456x 10 ^{i 0} (H2/HC) _¾ ; _ei

where /is the fraction representing conversion of Ce-convertibles.

[00207] Embodiment 67 is the computer system of any of embodiments 65 to 66 wherein the analytical solution is represented by:

T = T _a -pV(t t _max ) ln(l-1/ )

where T _a is the inlet temperature of at least one adiabatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst near start-of-rtm, wherein β is adjusted to fit the analytical solution to the fouling curve.

[00208] Embodiment 68 is the computer system of any of embodiments 59 to 67 wherein calculate the net present amount of benzene produced comprises use the run-length of the reforming catalyst and a benzene yield per month of the reforming catalyst to calculate an amount of benzene produced over the run-length per unit of the reforming catalyst and then discount the amount of benzene produced over the run-length of the reforming catalyst per unit of reforming catalyst to determine the net present amount of benzene produced per unit of reforming catalyst over the run-length of the reforming catalyst.

[00209] Embodiment 69 is the computer system of any of embodiments 59 to 68 wherein calculate a net present amount of benzene comprises calculate a number of pounds of benzene produced per pound of the reforming catalyst per hour, N, which is given by:

N = 0.907/ wee, LHSV fd _!f / d,) where d _c ,f, and wee are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of Ce-convertibles, and the weight fraction of Ce-convertibies in the feedstream.

[00210] Embodiment 70 is the computer system of embodiment 69 wherein N is greater than about 1.

[00211] Embodiment 71 is the computer system of any of embodiments 59 to 70 wherein calculate a net present, amount of benzene comprises cal culate a net, present pounds of benzene produced per pound of the reforming catalyst over the run-length of the reforming catalyst, NPPB, which is given by:

NPPB = (N/720) { 1 - (1 -D _m ) ⁿ } /¾

where D _m is the monthly discount rate and n is the run-length of the reforming catalyst, in units of months.

[00212] Embodiment 72 is the computer system of embodiment 71 wherein the NPPB is greater than about 14,000.

[00213] Embodiment 73 is the computer system of any of embodiments 69 to 72 wherein the NPPB is greater than about 14,000 or N is greater than about 1.

[00214] Embodiment 74 is the computer system of any of embodiments 59 to 73 wherein, for a given conversion of the (^-convertibles in the feedstream and the LHSV, a higher value for the net present amount of benzene produced is attained by selecting a

H2/HC ratio of the feedstream of greater than about 2 or less than about 1.

[00215] Embodiment 75 is the computer system of any of embodiments 59 to 74 wherein, for a given conversion of Ce-convertibles in the feedstream, the LHSV selected is about 4.

[00216] Embodiment 76 is the computer system of any of embodiments 59 to 75 wherein, for a given conversion of CVconvertibles in the feedstream, the H2/HC ratio of the feedstream is about 4.

[00217] Embodiment 77 is the computer system of any of embodiments 59 to 76 wherein the conversion of Ce-convertibles in the feedstream is from about 83% to about 93%, wherein the LHSV is about 4 and the H2/HC ratio of the feedstream is about 4.

[00218] Embodiment 78 is a. system, which may optionally be utilized in combination with any of embodiments 1 to 77, comprising a. catalytic reactor comprising a feedstream and a reforming catalyst; and a computer system coupled to the catalytic reactor and having a deactivation kinetic model to: kinetically model a deactivation of the reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstreani, a H2/HC ratio of the feedstreani, and a conversion of Co-convertibles in the feedstreani; and calculate a net present amount of benzene that is produced over the run-length of the reforming catalyst; wherein the catalytic reactor operates at values for the LHSV, the H2/HC ratio, and the conversion of Cg-convertibles that maximize the net present amount of benzene produced over a run-length of the reforming catalyst.

[00219] Embodiment 79 is the system of embodiment 78 wherein kinetically model the deactivation of the reforming catalyst with the deactivation kinetic model comprises calculate a fouling curve and determine the nin-!ength of the reforming catalyst at a time at which the fouling curve reaches a prescribed end of run temperature.

[00220] Embodiment 80 is the system of any of embodiments 78 to 79 wherein kinetically model the deactivation of the reforming catalyst with the deactivation kinetic model further comprises approximate the fouling curve with an analytical solution to the deactivation kinetic model.

[00221 [ Embodiment 81 is the system of embodiment 80 wherein the analytical solution is represented by:

where R is the gas constant; E _d is an activation energy of a deactivation reaction; Ε,· is an average activation energy of a conversion of C ₆ -, C ₇ - and Cg-convertibles; T ₀ is a start-of- run temperature; and k _d (To) is represented by A _d exp(~E _d / RT) with A _d being the average of A,j along at least one isothermal reactor and further represented by:

A _d ==1.343x 10 ^{5 !} -5.282X 10 ¹⁰ /-1.456x 1 ø' ^'' '( ! 12 I K ^' ).-, _. , _. ,

where /is the fraction representing conversion of C ₆ -convertibles.

[00222] Embodiment 82 is the system of any of embodiments 80 to 81 wherein the analytical solution is represented by:

T = T _a -β/it/ ) ½(l-t/t _max ) where T _a is the inlet temperature of at least one adiabatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst near start-of-run, wherein β is adjusted to fit the analytical solution to the fouling curve.

[00223] Embodiment 83 is the system of any of embodiments 78 to 82 wherein calculate the net present amount of benzene produced comprises use the run-length of the reforming catalyst and a benzene yield per month of the reforming catalyst to calculate an amount of benzene produced over the run-length per unit of the reforming catalyst and then discount the amount of benzene produced over the run-length of the reforming catalyst per unit of reforming catalyst to determine the net present amount of benzene produced per unit of reforming catalyst over the run-length of the reforming catalyst.

[00224] Embodiment 84 is the system of any of embodiments 78 to 83 wherein calculate a net present amount of benzene comprises calculate a number of pounds of benzene produced per pound of the reforming catalyst per hour, N, which is given by:

N = 0.907/ ^' wee LHSV i% / 4)

where dy, d _c , f, and wce are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of Cg-convertibles, and the weight fraction of (^-convertibles in the feedstream.

[00225] Embodiment 85 is the system of embodiment 84 wherein N is greater than about 1.

[00226] Embodiment 86 is the system of any of embodiments 78 to 85 wherein calculate a net present amount of benzene comprises calculate a net present pounds of benzene produced per pound of the reforming catalyst over the run-length of the reforming catalyst, NPPB, which is given by:

NPPB = (N/720) {1 - (1-D _m )"} /D _m

where∑)„, is the monthly discount rate and n is the run-length of the reforming catalyst in units of months.

[00227] Embodiment 87 is the system of embodiment 86 wherein the NPPB is greater than about 14,000.

[00228] Embodiment 88 is the system of any of embodiments 84 to 87 wherein the NPPB is greater than about 14,000 or N is greater than about 1. [00229] Embodiment 89 is the system of any of embodiments 78 to 88 wherei , for a given conversion of the Ce-convertibles in the feedstrearn and the LHSV, a higher value for the net present amount of benzene produced is attained by selecting a H2/HC ratio of the feedstrearn of greater than about 2 or less than about 1.

[00230] Embodiment 90 is the system of any of embodiments 78 to 89 wherein, for a given conversion of Ce-convertibles in the feedstrearn, the LHSV selected is about 4.

[00231] Embodiment 91 is the system, of any of embodiments 78 to 90 wherein, for a given conversion of Ce-convertibles in the feedstrearn, the H2 HC ratio of the feedstrearn is about 4.

[00232] Embodiment 92 is the system of any of embodiments 78 to 91 wherein the conversion of Ce-convertibles in the feedstrearn is from about 83% to about 93%, wherein the LHSV is about 4 and the H2/HC ratio of the feedstrearn is about 4.

[00233] Embodiment 93 is a non-transitory computer readable medium, which may be optionally be utilized in combination with any of embodiments 1 to 92, comprising a model unit to: kinetically model a deactivation of a reforming catalyst with a deactivation kinetic model to determine a run-length of the reforming catalyst as a function of a LHSV of a feedstrearn, a H2/HC ratio of the feedstrearn, and a conversion of (^-convertibles in the feedstrearn; and calculate a net present amount of benzene that is produced over a run- length of the reforming catalyst.

[00234] Embodiment 94 is the non -transitory computer readable medium of embodiment 93, further comprising a selection unit to select values for the LHSV, the H2/HC ratio, and the conversion of Ce-convertibles that maximize the net present amount of benzene produced over a run-length of the reforming catalyst.

[00235] Embodiment 95 is the non-transitory computer readable medium of any of embodiments 93 to 94, further comprising a sensitivity unit to determine a sensitivity of the values selected by the selection unit.

[00236] Embodiment 96 is the non-transitory computer readable medium of any of embodiments 93 to 95 wherein kinetically model the deactivation of the reforming catalyst with the deactivation kinetic model comprises calculate a fouling curve and determine the run-length of the reforming catalyst at a time at which the fouling curve reaches a prescribed end of run temperature. [00237] Embodiment 97 is the non-transitory computer readable medium of any of embodiments 93 to 96 wherein kineticaily model the deactivation of the reforming catalyst with the deactivation kinetic model further comprises approximate the fouling curve with an analytical solution to the deactivation kinetic model

[00238] Embodiment 98 is the non-transitory computer readable medium of embodiment, 97 wherein the analytical solution is represented by:

T(t) -— - r- -" — ]

1 + (RT ₆ I E _d jlnjl - (E _d I E _r )k _d (Γ ₀ )t\

where R is the gas constant; E _< j is an activation energy of a deactivation reaction; E _r is an average activation energy of a conversion of Ce-, C ₇ - and (^-convertibles; T ₀ is a start-of- run temperature; and k _d (To) is represented by A., exp(-E, _f / RT) with A _d being the average of along at least one isothermal reactor and further represented by:

A _d =1.343x 10 ^{3 i} -5.282X 10 ¹⁰ /-1.456x 10 ¹⁰ (H2/HC) _¾fcr

where /is the fraction representing conversion of Ce-convertibles.

[00239] Embodiment 99 is the non-transitory computer readable medium of any of embodiments 97 to 98 wherein the analytical solution is represented by:

T T _a ½(l-t/t _max )

where T _a is the inlet temperature of at least one adiabatic reactor, t is on-stream time, t _max is the time at which the reforming catalyst has near zero activity and the reactor temperature approaches infinity, and β is an adjustable parameter describing the fouling of the catalyst near start-of-run, wherein β is adjusted to fit the analytical solution to the fouling curve.

[00240] Embodiment 100 is the non-transitory computer readable medium of any of embodiments 93 to 99 wherein calculate a net present amount of benzene comprises calculate a number of pounds of benzene produced per pound of the reforming catalyst per hour, N, which is given by:

N = 0.907/ wee HSV {%/ d,)

where dj _f , d _c , f, and w¾6 are, respectively, the density of the feedstream, the density of the reforming catalyst, the fraction representing conversion of CVconvertibles, and the weight fraction of Ce-convertibles in the feedstream.

[00241] Embodiment 101 is the non-transitory computer readable medium of embodiment 100 wherein N is greater than about 1. [00242] Embodiment 102 is the non-transitory computer readable medium of any of embodiments 93 to 101 wherein calculate a net present amount of benzene comprises calculate a. net present pounds of benzene produced per pound of the reforming catalyst over the run-length of the reforming catalyst, NPPB, which is given by:

PPB = (N/720) { 1 - (1-D _m )"} /D _m

where D _m is the monthly discount rate and n is the run-length of the reforming catalyst in units of months,

[00243] Embodiment 103 is the non-transitory computer readable medium of embodiment, 102 wherein the NPPB is greater than about 14,000.

[00244] Embodiment 04 is the non- ransitory computer readable medium of any of embodiments 100 to 103 wherein the NPPB is greater than about 14,000 or N is greater than about 1 .

[00245] While several embodiments have been provided in the present disclosure, it should be understood that the disclosed systems and methods may be embodied in many other specific forms without departing from the spirit or scope of the present disclosure. The present examples are to be considered as illustrative and not restrictive, and the intention is not to be limited to the details given herein. For example, the various elements or components may be combined or integrated in another system or certain features may be omitted or not implemented.

[00246] Many variations and modifications of the inventive subject matter disclosed herein are possible and are within the scope of the inventive subject matter. Where numerical ranges or limitations are expressly stated, such express ranges or limitations should be understood to include iterative ranges or limitations of like magnitude falling within the expressly stated ranges or limitations (e.g., from about 1 to about 10 includes, 2, 3, 4, etc.; greater than 0.10 includes 0.1 1 , 0.12, 0.13, etc.). Use of the term "optionally" with respect to any element of a claim is intended to mean that the subject element is required, or alternatively, is not required. Both alternatives are intended to be within the scope of the claim. Use of broader terms such as comprises, includes, having, etc. should be understood to provide support for narrower terms such as consisting of, consisting essentially of, comprised substantially of, etc. [00247] The ordering of steps in the various methods are for illustration purposes and do not necessarily reflect the order that various steps must be performed. The steps may be rearranged in different orders in different embodiments to reflect the needs, desires and preferences of the entity impleme ting the systems. Furthermore, many steps may be performed simultaneously with other steps in some embodiments.

[00248] Also, techniques, systems, subsystems and methods described and illustrated in the various embodiments as discrete or separate may be combined or integrated with other systems, modules, techniques, or methods without departing from the scope of the present disclosure. Other items shown or discussed as directly coupled or communicating with each other may be coupled or communicating through some interface, device, or intermediate component such that the items may no longer be considered directly coupled to each other but may still be indirectly coupled and in communication, whether electrically, mechanically, or othenvise with one another. Other examples of changes, substitutions, and alterations are ascertainable by one skilled in the art and could be made without departing from the spirit and scope disclosed.

[00249] Accordingly, the scope of protection is not limited by the description set out above but is only limited by the claims which follow, that scope including all equivalents of the subject matter of the claims. Each and every claim is incorporated into the specification as an embodiment of the inventive subject matter. Thus, the claims are a further description and are an addition to the preferred embodiments of the inventive subject matter. The discussion of a reference in the disclosure is not an admission that it is prior art to the inventive subject matter, especially any reference that may have a publication date after the priority date of this application. The disclosures of ail patents, patent applications, and publications cited herein are hereby incorporated by reference, to the extent that they provide exemplary, procedural or other details supplementary to those set forth herein.