MONITORING THE ACTIVITY OF REFORMING CATALYTIC CONVERTERS

Description

BACKGROUND

[0001] The present invention relates generally to catalytic steam-hydrocarbon reforming and more specifically to determining changes in the catalytic activity of the reforming catalyst used in catalytic steam-hydrocarbon reformers. The determination of changes may be for decreases in the catalytic activity of the reforming catalyst as well as for increases due to refurbishment events or catalyst replacement.

[0002] Steam-hydrocarbon reforming catalysts are subject to numerous types of degradation, including sintering, carbon formation, and sulfur poisoning. Catalyst degradation results in a decline in the catalytic activity, which in turn reduces the efficiency and economic performance of the hydrogen or synthesis gas production facility. It is therefore important to be able to monitor the activity of the reforming catalyst and quantify the extent of reduced catalytic activity so that catalyst regeneration activities can be scheduled and failure of the reformer can be avoided.

[0003] Steam-hydrocarbon reforming catalysts may be refurbished in various ways, including replacement, partial replacement, and steaming. It is important to be able to receive an empirical confirmation that catalyst refurbishment activites have resulted in a significant improvement in catalytic activity, thereby allowing for continued operation of the reforming plant before refurbishment activities are repeated again.

[0004] EP 2 246 294 A1 discloses a control apparatus for controlling the supply of fuel and air to a fuel reformer. Degradation of a reforming catalyst of the fuel reformer is taken into account by detecting the temperature of the reforming catalyst and comparing the detected temperature with an estimated catalyst temperature. The estimated catalyst temperature is calculated on the basis of a thermal model of the catalyst. A target temperature for the catalyst, a detected temperature representative for fuel reformer heater and a correction coefficient calculated from the detected catalyst temperature and the estimated catalyst temperature are the inputs to the thermal model.

[0005] EP 1 835 141 A2 discloses a method of diagnosing a fuel injection system supplying fuel to a fuel reformer of an exhaust aftertreatment system. The method comprises measuring a temperature of a reformer exhaust stream and calculating an actual relationship between the measured temperature and a commanded fuel injection rate and, furthermore, comparing the actual relationship with an expected relationship between the commanded fuel injection rate and an estimated exhaust stream temperature derived from a thermal reformer model.

[0006] US 2012 / 0291424 A1 is concerned with diagnosing catalyst detoration in a fuel reforming system disposed in an exhaust gas recirculation apparatus of an internal combustion engine. An inlet temperature sensor is disposed upstream of the fuel reformer in order to detect an inlet temperature of the recycled exhaust gas upstream of the catalyst and an outlet temperature sensor is disposed downstream of the fuel reformer in order to detect an outlet temperature of the gas downstream of the catalyst. The catalyst is diagnosed on the basis of the difference between the detected inlet temperature and the detected outlet temperature.

[0007] In the prior art, a temperature approach to equilibrium has been used as an indicator for reduced activity of reforming catalyst. This approach requires collection of a sample of the reformate from the reformer furnace and offline measurement of the composition of the sample. The outlet temperature is estimated from the equilibrium constant for the water-gas shift reaction and the composition of the reformate at the outlet. The equilibrium temperature is calculated from the equilibrium constant for the steam reforming reaction at the measured composition for the reformate at the outlet. The temperature difference is the temperature approach to equilibrium.

[0008] For a fixed reformer operating condition, a small temperature approach to equilibrium (close to zero) means that the catalyst activity is high and that little degradation of the catalyst has occurred. On the other hand, a large temperature approach to equilibrium means that the catalyst activity has decreased and that the catalyst degradation has occurred. Although the temperature approach to equilibrium is widely used as a diagnostic tool for steam-hydrocarbon reformers, it has a major deficiency. In addition to being sensitive to catalyst degradation, the temperature approach to equilibrium also exhibits significant sensitivities to normal variations in furnace operating conditions, such as production rate, steam-to-carbon molar ratio, reformer operating temperature, and reformer operating pressure.

[0009] These variables all affect the conventionally calculated temperature approach to equilibrium, even for constant reforming catalyst conditions. The sensitivity of the conventionally calculated temperature approach to equilibrium to process conditions severely limits its utility as a diagnostic tool for catalyst degradation.

[0010] Industry desires new and / or improved methods for monitoring the catalytic activity of reforming catalyst for both increased catalytic activity and decreased catalytic activity.

[0011] Industry desires methods for monitoring the catalytic activity of reforming catalyst that is sensitive to changes in the catalyst activity but relatively insensitive to normal changes in process operating conditions.

[0012] Industry desires methods for monitoring catalytic activity of reforming catalyst using existing process data, and preferably process data which does not require a sample to be collected and analyzed offline.

[0013] Industry desires methods for enabling determination of changes in activity of reforming catalyst, in particular, for enabling determination of reduced activity of reforming catalyst, that is sensitive to changes in the catalyst activity but relatively insensitive to normal changes in process operating conditions.BRIEF SUMMARY

[0014] The present invention relates to a method for determining changes in the catalytic activity of reforming catalyst, either or both of increased catalytic activity and / or decreased catalytic activity as defined in claim 1.BRIEF DESCRIPTION OF SEVERAL VIEWS OF THE DRAWINGS

[0015] FIG. 1 is a process flow diagram of a hydrogen production process. FIG. 2 is a block flow diagram for an example using the method. FIG. 3 is a scatter plot of ΔT empricial versus ΔT approach , FIG. 4 is a graph of ΔT approach and of ΔT filtered versus time, overlaid with indications of detected catalyst degradation or refurbishment events. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0016] The ensuing detailed description provides preferred exemplary embodiments only, and is not intended to limit the scope, applicability, or configuration of the invention. Rather, the ensuing detailed description of the preferred exemplary embodiments will provide those skilled in the art with an enabling description for implementing the preferred exemplary embodiments of the invention, it being understood that various changes may be made in the function and arrangement of elements without departing from scope of the invention as defined by the claims.

[0017] The articles "a" and "an" as used herein mean one or more when applied to any feature in embodiments of the present invention described in the specification and claims. The use of "a" and "an" does not limit the meaning to a single feature unless such a limit is specifically stated. The article "the" preceding singular or plural nouns or noun phrases denotes a particular specified feature or particular specified features and may have a singular or plural connotation depending upon the context in which it is used.

[0018] The adjective "any" means one, some, or all indiscriminately of whatever quantity.

[0019] The term "and / or" placed between a first entity and a second entity includes any of the meanings of (1) only the first entity, (2) only the second entity, and (3) the first entity and the second entity. The term "and / or" placed between the last two entities of a list of 3 or more entities means at least one of the entities in the list including any specific combination of entities in this list. For example, "A, B and / or C" has the same meaning as "A and / or B and / or C" and comprises the following combinations of A, B and C: (1) only A, (2) only B, (3) only C, (4) A and B and not C, (5) A and C and not B, (6) B and C and not A, and (7) A and B and C.

[0020] The phrase "at least one of" preceding a list of features or entities means one or more of the features or entities in the list of entities, but not necessarily including at least one of each and every entity specifically listed within the list of entities and not excluding any combinations of entities in the list of entities. For example, "at least one of A, B, or C" (or equivalently "at least one of A, B, and C" or equivalently "at least one of A, B, and / or C") has the same meaning as "A and / or B and / or C" and comprises the following combinations of A, B and C: (1) only A, (2) only B, (3) only C, (4) A and B and not C, (5) A and C and not B, (6) B and C and not A, and (7) A and B and C.

[0021] The term "plurality" means "two or more than two."

[0022] For the purposes of simplicity and clarity, detailed descriptions of well-known devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.

[0023] The present invention relates to a method for determining changes in the catalytic activity of reforming catalyst in a catalytic reactor, such as a catalytic steam reforming furnace. The catalytic reactor may be used for generating hydrogen and / or synthesis gas. The method and system will be described with reference to FIG. 1 showing an exemplary process flow diagram of a production facility for generating hydrogen and / or synthesis gas.

[0024] While the method and system are described with reference to a catalytic steam reforming furnace such as shown in FIG. 1, the method and system are applicable to other types of reformers, for example, reformers with radiant and convective reforming as disclosed in U.S. Pat. No. 6,793,700B2, and tube-in-tube reformers as disclosed in U.S. Pat. Nos. 5,639,431, 5,938,800, 4,909,808, and catalyzed heat exchanger reformer as disclosed in U.S. 6,117,578.

[0025] Catalytic steam reforming, also called steam methane reforming (SMR) or steam reforming, is defined as any process used to convert reformer feedstock to synthesis gas by reaction with steam over a catalyst. Synthesis gas, commonly called syngas, is any mixture comprising hydrogen and carbon monoxide. The reforming reaction is an endothermic reaction and may be described generally as C n H m + n H 2 O → n CO + (m / 2 + n) H 2 . The water-gas shift reaction CO+H 2 O→CO 2 +H 2 , an exothermic reaction, also takes place in the reactor.

[0026] A production facility for generating hydrogen and / or synthesis gas is a reformer and associated equipment for the production of hydrogen and / or synthesis gas. Associated equipment may include adiabatic prereformers, heat exchangers, pumps, fans, water-gas shift reactors, pressure swing adsorbers, condensers, boilers, steam drums, desulphurizers, deaerators, headers, manifolds, piping, etc.

[0027] In the present method, a reformer feed gas mixture 15 is introduced into one or more catalytic reactors 20 in reformer 10, each of the one or more catalytic reactors 20 having an inlet and an outlet. The reformer feed gas mixture 15 comprises at least one hydrocarbon and steam and is reacted in a steam-methane reforming reaction and water-gas shift reaction under reaction conditions effective to form a reformate 25 comprising H 2 , CO, CH 4 , and H 2 O. The reformate 25 is withdrawn from the one or more catalytic reactors 20.

[0028] The reformer 10 may be a reformer furnace 10. Catalytic steam reforming takes place in reformer furnace 10. A reformer furnace, also called a catalytic steam reformer, or a steam methane reformer, is defined herein as any fired furnace used to convert feedstock containing elemental hydrogen and carbon to synthesis gas by a reaction with steam over a catalyst with heat provided by combustion of a fuel. The one or more catalytic reactors 20 may be catalyst-containing reformer tubes. Reformer furnaces with a plurality of catalyst-containing reformer tubes, i.e. tubular reformers, are well known in the art. Suitable materials and methods of construction are known. Catalyst in the catalyst-containing reformer tubes may be any suitable catalyst known in the art, for example, a supported catalyst comprising nickel.

[0029] The reformer feed gas mixture 15 comprises a hydrocarbon feedstock 75 and steam 151. Hydrocarbon feedstock may be natural gas, methane, naphtha, propane, refinery fuel gas, refinery off-gas, or other suitable reformer feedstock known in the art. Suitable operating temperatures range from 350°C to 650°C at the inlet and 750°C to 950°C at the outlet. Suitable pressures range from 1 to 50 atm. Preferred operating conditions for a catalytic steam reformer are known in the art. The hydrocarbon feedstock 75 may be heated and introduced along with hydrogen 106 into a hydrodesulphurization unit 155 to remove sulfur from the hydrocarbon feedstock for forming the reformer feed gas mixture 15. Hydrogen 106 for hydrodesulphurization may be added to the feedstock before or after heating the hydrocarbon feedstock 75. Hydrogen product 105 made by the hydrogen production facility may be used to provide hydrogen 106.

[0030] As used herein, a reformate is any mixture comprising hydrogen and carbon monoxide formed from the reforming reaction of a hydrocarbon and steam.

[0031] In a reformer furnace 10, a fuel 5 is combusted with an oxidant gas 7 in a combustion section 30 of the reformer furnace 10 external to the plurality of catalyst-containing catalytic reactors 20 (reformer tubes) under conditions effective to combust the fuel 5 to form a combustion product gas 35 comprising CO 2 and H 2 O and generate heat to supply energy for reacting the reformer feed gas mixture 15 inside the plurality of catalyst-containing reformer tubes 20. The combustion product gas 35 is then withdrawn from the combustion section 30 of the reformer furnace 10.

[0032] Any suitable burner may be used to introduce the fuel 5 and the oxidant gas 7 into the combustion section 30. Combustion of the fuel 5 with the oxidant gas 7 generates heat to supply energy for reacting the reformer feed gas mixture 15 inside the plurality of catalyst-containing reformer tubes 20. The combustion product gas 35 is withdrawn from the combustion section 30 of the reformer furnace 10 and passed to the convection section 45 of the reformer furnace 10 to supply heat to other process streams. The combustion section 30 (also called the radiant, radiation, or radiative section) of the reformer furnace is that part of the reformer furnace containing the plurality of catalyst-containing reformer tubes 20. The convection section 45 of the reformer furnace 10 is that part of the reformer furnace 10 containing heat exchangers other than the plurality of catalyst-containing reformer tubes 20. The heat exchangers in the convection section may be for heating process fluids other than reformate 25 from the plurality of catalyst-containing reformer tubes 20, such as water / steam, air 130, pressure swing adsorption unit by-product gas 115, reformer feed gas mixture 15 prior to introduction into the catalyst-containing reformer tubes 20, prereformed reformer feed gas (not shown), etc. The convection section 45 may also, if desired, contain a convective prereformer (not shown).

[0033] Furnace conditions effective to combust the fuel 5 may comprise a furnace temperature ranging from 600°C to 1500°C and a pressure ranging from 98 kPa to 101.4 kPa (absolute). Actual flame temperatures are generally higher. The temperature may be as measured by a thermocouple, an optical pyrometer, or any other calibrated temperature measurement device known in the art for measuring furnace temperatures. The pressure may be as measured by any suitable pressure sensor known in the art, for example a pressure gauge as available from Mensor.

[0034] The fuel 5 may be formed from a by-product gas 115 from a pressure swing adsorption unit 100 and a supplemental fuel 118. By-product gas 115 from a pressure swing adsorption unit 100 is often called pressure swing adsorber tail gas, and supplemental fuel 118 is often called trim fuel. The by-product gas 115 and / or the supplemental fuel 118 may be heated by indirect heat transfer with the combustion product gas 35 and / or the reformate 25 before being used as fuel. The by-product gas 115 and / or the supplemental fuel 118 may be heated to a temperature ranging from 100°C to 260°C. By-product gas 115 and supplemental fuel 118 may be blended and introduced together through a burner to the combustion section 30, or they may be introduced separately through different ports in the burner. Alternatively, the by-product gas 115 may be introduced through the primary burner and the supplemental fuel 118 may be introduced through lances near the burner.

[0035] The oxidant gas 7 is a gas containing oxygen and may be air, oxygen-enriched air, oxygen-depleted air such as gas turbine exhaust, industrial grade oxygen, or any other oxygen-containing gas known for use in a reformer furnace 10 for combustion. For example, as shown in the FIG. 1, air 130 may be compressed in forced draft fan 135, heated by the combustion product gas 35 in heat exchanger 4, and passed to the reformer furnace 10 as oxidant gas 7. Alternatively, (not shown) the air may be heated by a hot water / steam stream where the hot water / steam stream is heated by the combustion product gas 35.

[0036] Combustion product gas 35 may heat a number of different process streams in the convection section 45 of the reformer furnace 10. The combustion product gas 35 may heat the streams in various different configurations (order of heating). FIG. 1 shows an exemplary heat exchange network, although any suitable heat exchange network can be used with the present invention.

[0037] In FIG. 1, the combustion product gas 35 heats the reformer feed gas mixture 15 and then superheats steam stream 125 from steam drum 120. A portion 151 of the superheated steam may be used to form the reformer feed gas mixture 15 and another portion used to form a steam product 150 (i.e. export steam). After heating the steam, the combustion product gas then heats a portion of boiler feed water stream 127 from steam drum 120 in a heat exchanger to form a two-phase mixture of steam and water of which at least a portion is returned to the steam drum 120. The combustion product gas 35 then heats a water stream 123 from deaerator 110 after which the combustion product gas 35 heats compressed air stream 3 from compressor 135. Make up water stream 87 is heated by the combustion product gas 35 in heat exchanger 6 and passed to deaerator 110. After heating the make up water stream 87, the combustion product gas 35 is passed to an induced draft fan 140 and exhausted.

[0038] Reformate 25 withdrawn from the one or more catalytic reactors 20 may be passed to shift reactor 60 containing shift catalyst 61. The reformate 25 may exchange heat with any number of streams before being passed to the shift reactor 60. For example, the reformate 25 withdrawn from the one or more catalytic reactors 20 may be passed to heat exchanger 40 (a so-called waste heat boiler) where the reformate 25 heats a portion of boiler feed water stream 127 thereby forming a two-phase water and steam stream that is reintroduced into steam drum 120.

[0039] In the shift reactor 60, the reformate 25 is reacted in the presence of a shift catalyst 61 under reaction conditions sufficient to shift the reformate 25 to form additional H 2 in the reformate 25. Additional hydrogen gas may be obtained by the catalytic reaction of carbon monoxide and steam. This reaction is exothermic and is commonly referred to as the water-gas shift reaction or shift reaction: CO+H 2 O→CO 2 +H 2 . The reaction is affected by passing carbon monoxide and water through a bed of a suitable catalyst. The reaction conditions effective to form additional hydrogen in the reformate 25 may comprise a second temperature ranging from 190°C to 500°C and a second pressure ranging from 203 kPa to 5,066 kPa (absolute).

[0040] Any suitable shift catalyst may be used. The shift reactor may be a so-called high temperature shift (HTS), low temperature shift (LTS), medium temperature shift (MTS), or combination. One or more shift reactors may be used.

[0041] For high temperature shift, an inlet temperature in the range 310°C to 370°C, and an outlet temperature in the range 400°C to 500°C are typical. Usually an iron oxide / chromia catalyst is used for high temperature shift. High temperature shift may be preferred for the present process.

[0042] For low temperature shift, an inlet temperature in the range 190°C to 230°C, and an outlet temperature in the range 220°C to 250°C are typical. Usually a catalyst comprising metallic copper, zinc oxide, and one or more other difficulty reducible oxides such as alumina or chromia is used for low temperature shift

[0043] For medium temperature shift, an inlet temperature in the range 190°C to 230°C and an outlet temperature of up to 350°C are typical. A suitably formulated supported copper catalyst can be used for medium temperature shift.

[0044] A combination may include a sequence of high temperature shift, cooling by indirect heat exchange, and low temperature shift. If desired, either shift stage can be subdivided with interbed cooling.

[0045] In the exemplary embodiment shown in FIG. 1, after passing through the shift reactor 60, the reformate 25 is passed to heat exchanger 50 where the reformate 25 heats a stream of boiler feed water 123 that is withdrawn from deaerator 110. The reformate 25 is then passed to heat exchanger 70 where hydrocarbon feedstock 75 is heated and reformate 25 is cooled. Reformate is then passed to heat exchanger 80 to heat a second water feed stream 85 by indirect heat transfer the reformate 25, thereby cooling the reformate 25.

[0046] The second water feed stream 85, typically called make-up water, may be distilled water, treated water (decalcified, filtered, etc.) or other suitable water known in the art. The heated water stream 85 is passed to deaerator 110 and forms the boiler feed water 123. Low pressure steam 11 is introduced into deaerator 110 to assist with stripping dissolved gases from the condensate and make-up water. A vent stream 17 is withdrawn from the deaerator and comprises steam and gases formed from the dissolved gases.

[0047] To reduce VOC emissions from the hydrogen production facility, the deaerator vent streams from deaerator 110 may be injected into the reformer furnace 10 as described in the "Report on Emission Limits for Rule 1189 - Emissions from Hydrogen Plant Process Vents," South Coast Air Quality Management District, June 7, 2001 (http / / www3.aqmd.gov / hb / attachments / 2002 / 020620b.doc), and "Final Environmental Assessment: Proposed Rule 1189 - Emissions from Hydrogen Plant Process Vents" SCAQMD No. 1189JDN021199, South Coast Air Quality Management District December 17, 1999 (http: / / www.aqmd.gov / docs / default-source / ceqa / documents / aqmd-projects / 2000 / final-ea-for-proposed-amended-rule-1189.doc?sfvrsn=4).

[0048] The reformate 25 may be passed to a trim cooler (not shown) to condense at least a portion of the water in the reformate to form condensed water and a water-depleted reformate gas.

[0049] The condensed water is separated from the water-depleted reformate gas in separator 90 to form water condensate 97 from at least a portion of the condensed water. Water condensate 97 is blended with the second water feed stream 85. A slip stream of condensed water may be removed from the separator, if desired.

[0050] The water-depleted reformate is passed to a pressure swing adsorption unit 100 where the water-depleted reformate is separated to form a hydrogen-containing product gas 105 and a pressure swing adsorption unit by-product gas 115, also called PSA tail gas, and PSA purge gas. The fuel combusted in the combustion section 30 of the reformer furnace 10 may be formed from at least a portion of the pressure swing adsorption unit by-product gas 115.

[0051] The pressure swing adsorption unit 100 may be operated using any known pressure swing adsorption cycle. Pressure swing adsorption cycles are well-known in the art and the various steps are described for example in U.S. Pat. No. 9,381,460, incorporated herein by reference. Specific pressure swing adsorption cycles are provided in companion articles "Pressure Swing Adsorption cycles for 4 to 7 adsorption beds," IP.com number 000241449, April 29, 2015, and "Pressure Swing Adsorption cycles for 7 or more adsorption beds," IP.com number 000241619, May 18, 2015.

[0052] The method is characterized by repeatedly carrying out a number of steps over a period of time and for a plurality of times to monitor the activity of the reforming catalyst in the one or more catalytic reactors 20.

[0053] The method may be applied to monitor the activity of reforming catalyst in a single catalytic reactor (e.g. reformer tube), a group of catalytic reactors (e.g. group of reformer tubes), or all of the catalytic reactors (e.g. reformer tubes) in the reformer.

[0054] A first of the characterizing steps includes measuring an outlet temperature, T outlet , representative of a temperature at the outlet of the one or more catalytic reactors 20 for each time of the plurality of times. The outlet temperature representative of the temperature at the outlet of one or more catalytic reactors 20 may be measured at the outlet of one catalytic reactor, or in a subheader receiving reformate from more than one catalytic reactor feeding the subheader, or in the main header receiving reformate from all of the catalytic reactors in a reformer. In the case where the temperature is measured at the outlet of only a single catalytic reactor, the method would monitor the catalytic activity of the reforming catalyst in that single catalytic reactor. In the case where the temperature is measured in the subheader receiving reformate from a subset of catalytic reactors, the method would monitor the collective catalytic activity of the reforming catalyst in that subset of catalytic reactors. In the case where the temperature is measured in the main header receiving reformate from all of the catalytic reactors in the reformer furnace, the the method would monitor the collective catalytic activity of the reforming catalyst for all of the catalytic reactors in the reformer furnace.

[0055] It is well recognized that the measured outlet temperature is not a precise measure of the real outlet temperature for catalytic reactors due to the accuracy limitations of the temperature sensor and heat losses in the system. And while the measured outlet temperature is not a reliable measure for determining the species concentrations from equilibrium constant calculations, it is useful in the present method, since any bias is applied throughout.

[0056] A second of the characterizing steps includes determining a temperature approach to equilibrium for the steam-methane reforming reaction at the outlet of the one or more catalytic reactors 20 for each time of the plurality of times using the measured outlet temperature, T outlet . The temperature approach to equilibrium relates the measured outlet temperature, T outlet , to the equilibrium temperature, T equilibrium , corresponding to the composition of the reformate. The temperature approach to equilibrium, ΔT approach , is a measure of the difference between the outlet temperature, T outlet , and the equilibrium temperature, T equilibrium , and may conveniently be defined, ΔT approach = T outlet - T equilibrium , although other expressions are possible to describe the temperature approach to equilibrium.

[0057] The measure of the difference between the outlet temperature, T outlet , and the equilibrium temperature, T equilibrium , may be expressed as a function of (T outlet - T equilibrium ) and / or T outlet T equilibrium .

[0058] For example, the temperature approach to equilibrium, ΔT approach , may be expressed by any one of the following equations: ΔT approach = T outlet − T equilibrium , or ΔT approach = k * T outlet − T equilibrium , or ΔT approach = A + k * T outlet n − T equilibrium n , or Δ T approach = T outlet − T equilibrium T outlet , or Δ T approach = T outlet − T equilibrium T equilibrium , or Δ T approach = T outlet T equilibrium , or Δ T approach = T equlibrium T outlet , or Δ T approach = A + k T outlet n T equilibrium n , or Δ T approach = A + k T equilibrium n T outlet n , or Δ T approach = A + 1 − T outlet T equilibrium , or Δ T approach = A + k 1 − T equilibrium T outlet , or Δ T approach = A + k 1 − T outlet n T equilibrium n , or Δ T approach = A + k 1 − T equilibrium n T outlet n , where A is a constant value, k is another constant value, and n is another constant value.

[0059] It is clear that any equation form providing a measure of the difference between the outlet temperature, T outlet , and the equilibrium temperature, T equilibrium , may be used to describe the temperature approach to equilibrium, ΔT approach .

[0060] The equilibrium temperature, T equilibrium , for the reforming reaction may be calculated from any suitable equation for the equilibrium constant for the reforming reaction and a determined reformate composition representative of the reformate withdrawn from the one or more catalytic reactors 20 (i.e. the partial pressures, or mole fractions of H 2 , CO, CH 4 , and H 2 O).

[0061] From the definition of the equilibrium constant for the steam methane reforming reaction, it can be determined: K SR T equilibrium = n ˙ H 2 n ˙ total 3 n ˙ CO n ˙ total n ˙ CH 4 n ˙ total n ˙ H 2 O n ˙ total P outlet P o 2 , where ṅ is the molar flow rate of the specified species of the reformate, P outlet is the outlet pressure, P° is the standard state pressure, for example 100 kPa (1 bar), and ṅ total is the total molar flow rate of the reformate, n ˙ total = n ˙ CH 4 + n ˙ CO 2 + n ˙ CO + n ˙ N 2 + n ˙ H 2 + n ˙ H 2 O at the outlet of the one or more catalytic reactors 20.

[0062] The equilibrium temperature can be calculated from any known equation for the equilibrium constant for the steam methane reaction, for example, ln K SR = C 1 × ln T + C 2 T + C 3 + C 4 × T + C 5 T 2 + C 6 T 3 where C 1 = 8.611124, C 2 =-2.264801E4, C 3 (P in bar abs) = -2.898252E1, C 4 =-4.980062E-3, C 5 =3.977411E-7, and C 6 =2.013436E-11, as provided by Rostrup-Nielsen, Concepts in Syngas Manufacture, Imperial College Press, Catalytic Science Series, Vol. 10, p 317, 2011.

[0063] The reformate composition representative of the reformate withdrawn from the one or more catalytic reactors 20 may be measured in a chemical analyzer, for example as described in U.S. Pat. Appl. No. 2016 / 0131631, or determined from other reformer furnace measurements as described below.

[0064] The reformate composition representative of the reformate withdrawn from the one or more catalytic reactors 20 may be determined through a calculation using chemical element flow rate balances and a number of common furnace measurements, for example a hydrocarbon feedstock composition measurement, a flow rate measurement of the hydrocarbon feedstock, a flow rate measurement of the steam, and the measured outlet temperature.

[0065] Referring to FIG. 1, the composition of the hydrocarbon feedstock 75 may be measured using a chemical analyzer 76, such as a gas chromatograph, and the flow rate of hydrocarbon feedstock 75 may be measured using a flow meter 77. The flow rate of steam 151 forming the reformer feed gas mixture 15 may be measured using flow meter 152. The flow rate of hydrogen 106 to the hydrodesulphurization unit 155, if present, may be measured using flow meter 108.

[0066] The reformate composition representative of the reformate withdrawn from the one or more catalytic reactors 20 may be calculated using a carbon mass / mole balance for the reformer: n ˙ feedstock y CH 4 + 2 y C 2 H 6 + 3 y C 3 H 8 + 4 y C 4 H 10 + y CO 2 = n ˙ CO 2 + n ˙ CH 4 + n ˙ CO ; an oxygen mass / mole balance for the reformer: 2 n ˙ feedstock y CO 2 + n ˙ steam = 2 n ˙ CO 2 + n ˙ CO + n ˙ H 2 O ; a hydrogen mass / mole balance for the reformer: 2 n ˙ H 2 , feed + n ˙ feedstock 4 y CH 4 + 6 y C 2 H 6 + 8 y C 3 H 8 + 10 y C 4 H 10 + 2 n ˙ steam = 2 n ˙ H 2 O + 4 n ˙ CH 4 + 2 n ˙ H 2 ; and a nitrogen mass / mole balance for the reformer: 2 n ˙ feedstock y N 2 = 2 n ˙ N 2 , where ṅ H2,feed is the molar flow rate of the hydrogen feed 106 to the hydrodesulphurization unit 155, if present, ṅ feedstock is the molar flow rate of the hydrocarbon feedstock 75, ṅ steam is the molar flow rate of steam 151, ṅ CO2 is the molar flow rate of CO 2 in the reformate 25, ṅ CH4 is the molar flow rate of CH 4 in the reformate, ṅ CO is the molar flow rate of CO in the reformate, ṅ H2O is the molar flow rate of H 2 O in the reformate, and y i , is the mole fraction of the respective component, i, in the hydrocarbon feedstock.

[0067] Equations 4-7 are for an exemplary hydrocarbon feedstock 75 that includes methane through butane. The equations can be easily modified for fewer or additional species in the hydrocarbon feedstock.

[0068] The methane slip is another parameter often measured. The reformate composition at the outlet of the one or more catalytic reactors may be calculated using a measurement representative of the methane concentration, such as determined from the methane slip. The dry gas methane concentration after the separator (knock-out drum) 90 and prior to introducing the reformate into the pressure swing adsorption unit 100 may be measured using chemical analyzer 96, which may comprise a metal oxide semiconductor sensor, such as Model M2A 65-2649RK-CH4 available from RKI Instruments. The shift reactor doesn't react the methane in the reformate and the molar flow rate of methane doesn't change from the outlet of the reformer to the inlet of the pressure swing adsorption unit 100. Then, n ˙ CH 4 = n ˙ CH 4 , dry = y CH 4 , dry × n ˙ total , dry where ṅ total,dry is the total molar flow rate on a dry basis, and y CH4,dry is the measured dry basis methane mole fraction. The total dry basis molar flow rate, ṅ total,dry introduced into the pressure swing adsorption unit 100, may be measured by a flow meter 98 and the molar flow rate of methane, ṅ CH4 calculated from the total dry basis molar flow rate and the measured dry basis methane mole fraction. The calculation of the reformate composition may also use a measured flow rate representative of the flow rate of the reformate, e.g. the dry-basis flow rate.

[0069] Alternatively, a sample of the reformate may be taken before reacting the reformate in shift reactor 60, the water in the sample condensed and removed, and a dry-basis methane mole fraction, y CH4,outlet , measured, where y CH 4 , outlet = n ˙ CH 4 n ˙ CO 2 + n ˙ CH 4 + n ˙ CO + n ˙ H 2 + n ˙ N 2 .

[0070] A constraint based on the water-gas shift equilibrium constant can also be used for determining an effective reformate composition at the outlet of the one or more catalytic reactors: K shift T outlet = n ˙ H 2 × n ˙ CO 2 n ˙ CO × n ˙ H 2 O where the equilibrium constant is evaluated at the measured outlet temperature, T outlet . Any suitable equation for the equilibrium constant for the shift reaction may be used, for example, the equation by Moe "Design of water-gas shift reactors," Chemical Engineering Progress 58, pp 33-36, 1962, K shift (T outlet )=exp(4577.8 / T outlet - 4.33).

[0071] An algorithm for the solution of a system of nonlinear equations, such as the Levenberg-Marquardt or trust-region-reflective method may be used to calculate the molar flow rates, n, using the equations 4-9. The Levenberg-Marquardt method is described, for example, in Levenberg, "A Method for the Solution of Certain Problems in Least Squares," Quart. Appl. Math., Vol. 2, pp 164-168, 1944, and Marquardt, "An Algorithm for Least-Squares Estimation of Nonlinear Parameters," SIAM, J. Appl. Math., Vol. 11, pp 431-441, 1963.

[0072] With the molar flow rates of the reformate species, the equilibrium temperature for the steam-methane reforming reaction can be calculated from equations 1-3 and the reformate pressure, which may be measured or otherwise known. Subsequently the temperature approach to equilibrium may be calculated from the outlet temperature, T outlet , and the equilibrium temperature, T equilibrium .

[0073] The temperature approach to equilibrium calculated using the measured outlet temperature differs from the calculation of the conventional temperature approach to equilibrium since it is recognized that any measured value using a temperature sensor is inaccurate. In the calculation of the conventional temperature approach to equilibrium, the outlet temperature is calculated based on a measured composition of the reformate withdrawn from the catalyst-containing reformer tubes. A sample of the reformate is collected and the concentrations of H 2 , CO 2 , CO, and H 2 O are measured. At temperatures typical for the outlet temperature of the reformer tubes, the water-gas shift reaction is assumed to be at equilibrium. Using the equilibrium constant for the shift reaction, K shift T shift , equilibrium = y H 2 × y CO 2 y CO × y H 2 O , the equilibrium temperature, T shift, equilibrium , is calculated from the measured mole fractions, y, of H 2 , CO 2 , CO, and H 2 O. Any suitable equation for the equilibrium constant for the shift reaction may be used, for example, the equation by Moe "Design of water-gas shift reactors," Chemical Engineering Progress 58, pp 33-36, 1962, K shift (T shift, equilibrium )=exp(4577.8 / T- 4.33). So for example, if the mole fractions of H 2 , CO 2 , CO, and H 2 O are assumed to be 0.6, 0.05, 0.2, and 0.05, respectively, the equilibrium temperature for the shift reaction, T shift,equilibrium , is calculated to be 843K.

[0074] Likewise, using the measured concentrations of the reformate, the equilibrium temperature for the steam-methane reforming reaction is calculated using equations 1-3 and the outlet pressure. Equation 3 may be used or an alternative expression.

[0075] For the conventional temperature approach to equilibrium, the temperature approach to equilibrium is calculated as ΔT approach, conventional = T shift, equilibrium - T equilibrium ,

[0076] A third of the characterizing steps includes calculating an empirical model-based temperature approach to equilibrium using furnace operating data acquired at each time of the plurality of times in an empirical model based on historical operating data for the reformer or another reformer of the same type for each time of the plurality of times. The empirical model-based temperature approach to equilibrium is a "reference" value taking into account variations in operating conditions for the reformer.

[0077] The empirical model for calculating the empirical model-based temperature approach to equilibrium is preferably based on furnace operating data from the reformer being monitored for reduced activity of the reforming catalyst. In the case of a new reformer, historical data is not available and the empirical model may be based on a similar reformer of the same type. A reformer is of the same type as the monitored reformer if the geometry of the catalytic reactors is the same. For example, a reformer having a plurality of catalyst-containing reformer tubes heated by combustion of a fuel and oxidant gas external to the plurality of catalyst-containing reformer tubes is of the same type as another reformer having a plurality of catalyst-containing reformer tubes heated by combustion of a fuel and oxidant gas external to the plurality of catalyst-containing reformer tubes.

[0078] The empirical model may be developed from historical data from a reformer. Selected historical data may be used for the development of the empirical model, for example data when the catalyst is known to have high catalytic activity or any and all historical data may be used.

[0079] The empirical model relates an empirical temperature approach to equilibrium to several operating parameters for the reformer. Like the temperature approach to equilibrium in the first characterizing step, the empirical temperature approach to equilibrium is based on a measured temperature representative of a temperature at the outlet of the one or more catalytic reactors 20. The measured temperature for developing the empirical model is preferably at the same temperature sensor location as for measuring the outlet temperature in the first characterizing step.

[0080] The empirical model may be formulated from historical outlet temperatures, T outlet, historical , for the model reformer for a historical plurality of times, and calculated historical equilibrium temperatures, T equilibrium, historical , for the model reformer for the historical plurality of times. The empirical model may be formulated from the empirical model-based temperature approach to equilibrium, ΔT empirical, historical , for the historical plurality of times.

[0081] Consistent with the definition of the temperature approach to equilibrium, the empirical model-based temperature approach to equilibrium, ΔT empirical, historical , is a measure of a difference between historical outlet temperatures, T outlet, historical , for the model reformer for a historical plurality of times, and calculated historical equilibrium temperatures, T equilibrium, historical , for the model reformer for the historical plurality of times. The empirical model-based temperature approach to equilibrium, ΔT empirical, historical , may be most conveniently defined ΔT empirical, historical = T outlet, historical - T equilibrium, historical , although other expressions are possible to describe the empirical model-based temperature approach to equilibrium as a measure of the difference between the historical outlet temperatures, T outlet, historical , and the calculated historical equilibrium temperatures, T equilibrium, historical .

[0082] The measure of the difference between the historical outlet temperatures, T outlet, historical , and the calculated historical equilibrium temperatures, T equilibrium, historical , may be expressed as a function of (T outlet, historical - T equilibrium, historical ) and / or T outlet , historical T equilibrium , historical .

[0083] For example, the temperature approach to equilibrium, ΔT empirical, historical , may be expressed by any one of the following equations: ΔT empirical , historical = T outlet , historical − T equilibrium , historical , or ΔT empirical , historical = k * T outlet , historical − T equilibrium , historical , or ΔT empirical , historical = A + k * T outlet , historical n − T equilibrium , historical n , or Δ T empirical , historical = T outlet , historical − T equilibrium , historical T outlet , historical , or Δ T empirical , historical = T outlet , historical − T equilibrium , historical T equilibrium , historical , or Δ T empirical , historical = T outlet , historical T equilibrium , historical , or Δ T empirical , historical = T equlibrium , historical T outlet , historical , or Δ T empirical , historical = A + k T outlet , historical n T equilibrium , historical n , or Δ T empirical , historical = A + k T equilibrium , historical n T outlet , historical n , or Δ T empirical , historical = A + 1 − T outlet , historical T equilibrium , historical , or Δ T empirical , historical = A + k 1 − T equilibrium , historical T outlet , historical , or Δ T empirical , historical = A + k 1 − T outlet , historical n T equilibrium , historical n , or Δ T empirical , historical = A + k 1 − T equilibrium , historical n T outlet , historical n , where A is a constant value, k is another constant value, and n is another constant value, and wherein the equation selected for the empirical model-based temperature approach to equilibrium, ΔT empirical, historical , is consistent with the equation selected for the temperature approach to equilibrium, ΔT approach .

[0084] The historical equilibrium temperatures, T equilibrium, historical , is calculated in the same manner as described for the equilibrium temperature, T equilibrium , using any of equations 1-9.

[0085] The historical equilibrium temperatures, T equilibrium, historical , may be calculated from the equilibrium constant for the steam-methane reforming reaction at historical reformate compositions representative of reformate withdrawn from the one or more catalytic reactors of the model reformer for the historical plurality of times.

[0086] The empirical model is formulated from additional historical operating data including an historical reformate pressure representative of the pressure at the outlet of the one or more catalytic reactors, P outlet, historical , an historical flow rate of the steam, n steam , historical , and an historical flow rate of the feedstock, n feedstock, historical , for the historical plurality of times.

[0087] Coefficients for the empirical model is determined by regressing ΔT empirical, historical as a function of the historical outlet temperatures, T outlet , the historical reformate pressures representative of the pressure at the outlet of the one or more catalytic reactors, P outlet, historical , the historical flow rate of the steam, n steam , historical , preferably using a steam-to-carbon molar ratio, and the historical flow rate of the feedstock, n feedstock, historical .

[0088] The function for the empirical model is a linear or non-linear equation with the coefficients determined by regressing the historical operating data, where ΔT empirical, historical = T outlet, historical - T equilibrium, historical for the regression. If an alternative equation is used to describe the measure of the difference between historical outlet temperatures, T outlet, historical , and calculated historical equilibrium temperatures, T equilibrium, historical , the alternative form of the empirical model-based temperature approach to equilibrium may be regressed.

[0089] The empirical model function may include variables derived from additional historical operating data, for example, a steam-to-carbon molar ratio may be used as a variable representative of the steam flow rate. The variable used in the equation for the empirical model may include, for example, hydrocarbon feedstock flow rate, ṅ feedstock , reformer tube outlet pressure, P outlet , steam-to-carbon ratio, S / C, and the temperature representative of a temperature at the outlet of the one or more catalytic reactors, T outlet .

[0090] The steam-to-carbon molar ratio is a conventional term in the field of hydrogen production and syngas production. The steam-to-carbon ratio may be expressed as a mass ratio or more commonly, a molar ratio.

[0091] The steam-to-carbon molar ratio of the feed is the ratio of the molar flow rate of steam in the feed to the molar flow rate of hydrocarbon-based carbon in the feed at the mixing Tee. The molar flow rate of hydrocarbon-based carbon is the molar flow rate of carbon where the carbon is associated with hydrocarbons (i.e. excluding carbon associated with carbon monoxide and carbon dioxide). For example, if the total molar flow rate of the feed is 100 moles / h, and the mole fraction of methane is 0.35, the mole fraction of ethane is 0.1, and the mole fraction of carbon monoxide is 0.05, then the molar flow rate of hydrocarbon-based carbon is 55 moles / h. Methane contributes 35 moles / h of hydrocarbon based carbon. Ethane contributes 20 moles / h of hydrocarbon-based carbon. And carbon monoxide contributes zero moles / h of hydrocarbon based carbon.

[0092] An equation or function for the empirical model may be developed by regressing the historical data. In general terms, the empirical temperature approach to equilibrium, ΔT empirical = f(T outlet , P outlet, S / C, n feedstock ), where T outlet , is a measured outlet temperature, P outlet is a measured outlet pressure, S / C is the steam-to-carbon molar ratio, and is a flow rate of the hydrocarbon feedstock.

[0093] The functional form for the empirical model can be linear or nonlinear. A linear example is Δ T empirical = A + B × n ˙ feedstock + C × P outlet + D × S C + E × T outlet . A nonlinear example is ΔT empirical = A + B×(ṅ feedstock ) b< + C×(P outlet ) c< + D×(S / C) d< + E×(T outlet ) e< .

[0094] In the development of the empirical model, historical information is used in the functional form and the data regressed to determine the coefficients that best fit the data according to a selected criteria.

[0095] For example, for the linear example, Δ T empirical , historical = A + B × n ˙ feedstock , historical + C × P outlet , historical + D × S C historical + E × T outlet , historical .

[0096] Regression of the historical data to determine the coefficients for the linear and nonlinear equations is routine.

[0097] Having the empirical model, for each time, the operating parameters for the reformer are determined through measurements, and the empirical model-based temperature approach to equilibrium is calculated per the third characterizing step.

[0098] A fourth of the characterizing steps of the method includes comparing the temperature approach to equilibrium using the measured outlet temperature to the calculated empirical model-based temperature approach to equilibrium for each time of the plurality of times.

[0099] The comparison between the temperature approach to equilibrium, ΔT approach , using the measured outlet temperature and the empirical model-based temperature approach to equilibrium, ΔT empirical, is facilitated using a characteristic operational value ΔT residual derived from the temperature approach to equilibrium and the calculated empirical model-based temperature approach to equilibrium. The characteristic operational value ΔT residual is a residual temperature approach to equilibrium, ΔT residual , defined as a difference, ΔT residual = ΔT approach - ΔT empirical . Other equivalently functioning characteristic operational values may be derived from differences and / or ratios of the temperature approach to equilibrium and the calculated empirical model-based temperature approach to equilibrium to achieve the same result.

[0100] The characterizing steps are repeated for the plurality of times in order to monitor the catalytic activity of the reforming catalyst. The characterizing steps may be performed at least once per month, or at least once per week, or at least once per day, or at least once per hour, or at least once per minute during operation of the reformer 10.

[0101] The method allows to determine decreased and / or increased catalytic activity of the reforming catalyst. The method for determining decreased and / or increased catalytic activity of reforming catalyst comprises the characterizing steps where the characterizing steps are repeated during a time period where the temperature approach to equilibrium relative to the empirical model-based temperature approach to equilibrium for each time of the plurality of times differs by an amount for a selected period of time, the amount and selected period of time determined to indicate decreased and / or increased activity of the reforming catalyst.

[0102] As time progresses from a initial installation of fresh catalyst or a catalyst refurbishment event, the amount by which the temperature approach to equilibrium relative to the empirical model-based temperature approach to equilibrium differs, e.g. ΔT residual , will tend to increase and the increase will be readily apparent. The amount by which the temperature approach to equilibrium relative to the calculated empirical model-based temperature approach to equilibrium for each time differs will, in general, increase as the activity of the catalyst decreases and therefore is an indication of decreased catalytic activity of the reforming catalyst.

[0103] Decreased activity of the reforming catalyst is determined by an amount the temperature approach to equilibrium differs from the empirical model-based temperature approach to equilibrium and the period of time the amount exceeds an amount indicating decreased catalytic activity.

[0104] The criteria for indicating reduced activity of the reforming catalyst is established where ΔT residual increases by more than a selected amount for a selected period of time (i.e. ΔT residual increases by more than the selected amount from the initial installation of fresh catalyst or catalyst refurbishment event and the increase amount is maintained for a selected period of time).

[0105] The amount and selected period of time indicating decreased catalytic activity of the reforming catalyst and whether to act on it, may depend on economic considerations due to decreasing efficiency of the reformer resulting from the decreased catalytic activity. The amount and selected period of time indicating decreased activity of the reforming catalyst may be decided from a predetermined rule for the amount and time period. The selected period of time may be at least one or more days, or at least one or more weeks, or at least one or more months, and may be as long as 6 months or more depending on the severity of the catalyst degradation.

[0106] For example, a threshold rule may be applied such that reduced activity is deemed to have occurred the first time that ΔT residual has increased by a specified amount, such as 3 °C. Or, as described in the example, a filtered value of ΔT residual may be used in judging when reduced activity is deemed to have occurred. Various other rules for detecting a significant degradation in catalyst activity may be developed using the well-known literature on control charts, quickest detection, and linear and nonlinear filtering. These rules include Shewhart control charts, cumulative sum control charts, exponentially weighted moving average control charts, Bayesian quickest detectors, particle filters, and Kalman filters. See, for example: 1) Poor, H. Vincent, and Olympia Hadjiliadis, Quickest detection, Cambridge: Cambridge University Press, 2009; 2) Isom, Joshua D. Exact solution of Bayes and minimax change-detection problems, Dissertation, University of Illinios at Urbana-Champaign, 2009; 3) Basseville, Michèle, and Igor V. Nikiforov, Detection of abrupt changes: theory and application, Englewood Cliffs: Prentice Hall, 1993; 4) Arulampalam, M. Sanjeev, et al. "A tutorial on particle filters for online nonlinear / non-Gaussian Bayesian tracking," IEEE Transactions on signal processing 50.2 (2002): 174-188; and 5) Harvey, Andrew C. Forecasting, structural time series models and the Kalman filter, Cambridge University Press, 1990.

[0107] The amount determined to indicate decreased catalytic activity of the reformer catalyst may be an amount corresponding to an increase in a value of a residual temperature approach to equilibrium, ΔT residual , of 3°C or more during the selected period of time, where ΔT residual = ΔT approach - ΔT empirical, where ΔT approach = T outlet - T equilibrium , and corresponds to the temperature approach to equilibrium determined in the second characterizing step, where T outlet is the measured temperature representative of the temperature at the outlet of the one or more catalytic reactors (20), where T equilibrium is a temperature calculated from an equilibrium constant for the steam-methane reforming reaction at a reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20), and where ΔT empirical corresponds to the calculated empirical model-based temperature approach to equilibrium for the regression where ΔT empirical, historical = T outlet, historical - T equilibrium , historical for the regression.

[0108] Increased catalytic activity of the reforming catalyst may also be determined to observe the effect of a catalyst replacement or catalyst refurbishment event. In case of a refurbishment event, additional refurbishment could be applied in case the increase in the catalytic activity is not considered sufficient.

[0109] The amount determined to indicate increased catalytic activity of the reformer catalyst may be an amount corresponding to a decrease in a value of a residual temperature approach to equilibrium, ΔT residual , of 3°C or more during the selected period of time, where ΔT residual = ΔT approach - ΔT empirical , where ΔT approach = T outlet - T equilibrium , and corresponds to the temperature approach to equilibrium determined in the second characterizing step, where T outlet is the measured temperature representative of the temperature at the outlet of the one or more catalytic reactors (20), where T equilibrium is a temperature calculated from an equilibrium constant for the steam-methane reforming reaction at a reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20), and where ΔT empirical corresponds to the calculated empirical model-based temperature approach to equilibrium for the regression where ΔT empirical, historical = T outlet, historical - T equilibrium, historical for the regression.

[0110] With reference to FIG. 1, the system comprises a temperature sensor 26 which is operable to measure an outlet temperature, T outlet , where T outlet is representative of a temperature at the outlet of one or more catalytic reactors 20. Temperature sensor 26 is configured to transmit temperature information relating to the measured outlet temperature.

[0111] The system also comprises a computing device 200 operable to receive operating information from the reformer furnace 10 including the temperature information from the temperature sensor 26. Computing device 200 is operable to determine a temperature approach to equilibrium for the steam-methane reforming reaction at the outlet of the one or more catalytic reactors 20 over a period of time and for a plurality of times. Computing device 200 is operable to calculate an empirical model-based temperature approach to equilibrium from an empirical model based on historical operating data for the reformer furnace 10 or another reformer furnace for each time of the plurality of times. Computing device 200 is capable of providing an output suitable for comparing the temperature approach to equilibrium to the calculated empirical model-based temperature approach to equilibrium for each time over the period of time in order to monitor for reduced activity of the reforming catalyst.

[0112] The temperature approach to equilibrium, ΔT approach , is a measure of the difference between the outlet temperature, T outlet , and the equilibrium temperature, T equilibrium . The empirical model-based temperature approach to equilibrium, ΔT empirical, historical, is a measure of a difference between historical outlet temperatures, T outlet, historical , for the model reformer for a historical plurality of times, and calculated historical equilibrium temperatures, T equilibrium, historical , for the model reformer for the historical plurality of times.

[0113] Computing device 200 may be a programmable logic controller.

[0114] Computing device 200 is operable to perform any of the calculations of the present method.

[0115] The one or more catalytic reactors 20 may be a plurality of catalyst-containing reformer tubes in a reformer furnace.

[0116] The system may further comprise one or more sensors operable to obtain at least a subset of the operating information from the reformer furnace 10 and configured to transmit at least the subset of the operating information to the computing device 200.

[0117] The one or more sensors may include a pressure sensor 27 operatively disposed to determine an outlet pressure representative of the pressure at the outlet of the one or more catalytic reactors 20 where the pressure sensor 27 is operatively connected to the computing device 200. The connection between the pressure sensor 27 and the computing device 200 may be hardwired or wireless. The operating information received by the computing device 200 may then include the outlet pressure representative of the pressure at the outlet of the one or more catalytic reactors 20.

[0118] The system may further comprise a chemical analyzer 76 operatively disposed to measure a composition of a sample representative of a hydrocarbon feed 75 where the reformer feed gas mixture 15 comprises the hydrocarbon feed 75. If present, the chemical analyzer 76 is operatively connected to the computing device 200. The connection between the chemical analyzer 76 and the computing device may be hardwired or wireless. The operating information received by the computing device 200 may then include the composition of the sample.

[0119] The system may further comprise a flow meter 77 operatively disposed to measure a flow rate of the hydrocarbon feed 75 where the flow meter 77 is operatively connected to the computing device 200. The connection between the flow meter 77 and the computing device 200 may be hardwired or wireless. The operating information received by the computing device 200 may then include the flow rate of the hydrocarbon feed 75.

[0120] The system may further comprise a flow meter 152 operatively disposed to measure a flow rate of steam 151 where the reformer feed gas mixture 15 comprises the steam 151 where the flow meter 152 is operatively connected to the computing device 200. The connection between the flow meter 152 and the computing device 200 may be hardwired or wireless. The operating information received by the computing device 200 may then include the flow rate of the steam 151.

[0121] The system may further comprise a (reformate) chemical analyzer 96 operatively disposed to measure a concentration representative of at least one component in the reformate 25. If present, the chemical analyzer 96 is operatively connected to the computing device 200. The connection between the chemical analyzer 96 and the computing device 200 may be hardwired or wireless. The operating information received by the computing device 200 may then include the concentration representative of the at least one component in the reformate 25. The at least one component may be methane. Chemical analyzer 96 may have a metal oxide semiconductor to determine the concentration of methane. Chemical analyzer 96 may be Model M2A 65-2649RK-CH4 available from RKI Instruments.

[0122] The chemical analyzer 96 may be operatively disposed to measure the composition of the reformate 25 prior to the reformate 25 being passed to the shift reactor 60. Alternatively, the chemical analyzer 96 may be operatively disposed to measure a water-depleted reformate 95 withdrawn from separator 90.

[0123] The system may further comprise a flow meter 98 operatively disposed to measure a flow rate of the water-depleted reformate 95 where the flow meter 98 is operatively connected to the computing device 200. The connection between the flow meter 98 and the computing device may be hardwired or wireless. The operating information received by the computing device 200 may then include the flow rate of the water-depleted reformate 95.

[0124] The system may further comprise a flow meter 108 operatively disposed to measure a flow rate of a hydrogen feed 106 where the reformer feed gas mixture 15 comprises the hydrogen feed 106 where the flow meter 108 is operatively connected to the computing device 200. The connection between the flow meter 108 and the computing device may be hardwired or wireless. The operating information received by the computing device 200 may then include the flow rate of the hydrogen feed 106.EXAMPLE

[0125] The method and system was applied to a steam reformer facility with a catalytic steam-hydrocarbon reformer, like that shown in FIG. 1. The reformer was instrumented with a chemical analyzer 76 for measuring the feedstock composition, a flow meter 77 for measuring the flow rate of the feedstock 75, a flow meter 108 for measuring the flow rate of hydrogen 106 to the hydrodesulphurization unit 155, a flow meter 152 to measure the flow rate of process steam 151, a temperature sensor 26 to measure the collective temperature of the reformate from the plurality of reformer tubes 20, a chemical analyzer 96 for measuring the methane slip on a dry basis, and a pressure sensor 27 for measuring the pressure of the reformate exiting the catalyst-containing reformer tubes 20.

[0126] The frequency of obtaining measurements was once per minute during hydrogen production.

[0127] A block diagram for the method of the example is shown in FIG. 2. Five sets of inputs were used to determine the reformate composition. The five sets of inputs included 1) a measured composition of the hydrocarbon feed, 2) a measured flow rate of the hydrocarbon feed, 3) a measured steam flow rate, 4) a dry-basis methane concentration at the reformer outlet, and 5) measured reformate outlet temperature. The molar flow rates of CH 4 , H 2 O, CO, CO 2 , H 2 , and N 2 at the reformer outlet (a total of six unknowns) was determined by solving a set of six nonlinear equations: 1) hydrogen atom balance (eq. 6), 2) oxygen atom balance (Eq. 5), 3) carbon atom balance (Eq. 4), 4) nitrogen atom balance (Eq. 7), 5) dry basis methane concentration at the reformer exit (Eq. 8 alt.), and 6) water-gas shift equilibrium at the reformer exit (Eq. 9). The set of six nonlinear equations was solved using the MATLAB ®< function "Isqnonlin". The determined molar flow rates at the reformer outlet, in conjuction with the measured reformer outlet pressure P outlet , was used to calculate the reform equilibrium temperature T equilibrium using Equation (1). The temperature approach to equilibrium ΔT approach was calculated as ΔT approach = T outlet - T equilibrium .

[0128] Five years of historical data were used to calculate a linear regression model for ΔT approach as a function of T outlet , P outlet , the steam-to-carbon molar ratio (S / C), and the hydrocarbon feedstock flow rate (ṅ reedstock ). The linear regression model is: ΔT empricial = − 689 ∘ C − 0.3769 * n ˙ feedstock − 0.007154 * P outlet + 24.41 * S / C + 0.7540 * T outlet wth P outlet measured in units of kPa, T outlet measured in units of °C, ṅ feedstock measured in standard m 3< / s (standard conditions being 101.325 kPa and 15°C), and S / C a unitless molar ratio. A statistical test was used to reject, with high confidence, the null hypotheses that each of the regression coefficients is equal to zero. The R-squared statistic associated with the regression model is 0.615, and the root mean squared error of the regression model is 2.94 °C. These fit statistics for the regression model are acceptable, because it is known that changes in catalyst activity should affect ΔT approach but it is not intended that the regression model capture these effects of catalyst activity. As explained above, the purpose of the regression model is to model the normal sensitivity of ΔT approach to normal changes in process operating conditions. FIG. 3 is a scatter plot of ΔT empirical versus ΔT approach, which illustrates the normal variability of ΔT approach and the degree of correlation between ΔT empirical and ΔT approach .

[0129] A residual temperature approach to equilibrium was then calculated as ΔT residual = ΔT approach − ΔT empirical .

[0130] Six year's worth of data (from mid-2010 through mid-2016 were used to calculate the residual temperature approach to equilibrium, ΔT residual , on an hourly basis for the facility. FIG. 4 shows a plot of the residual temperature approach to equilibrium, ΔT residual , versus the date for the six year period. A blind test was then used to assess the utility of the method and system for identifying catalyst degradation events and catalyst refurbishment events.

[0131] In conducting the validation test, ΔT residual was used to detect catalyst degradation and refurbishment events using a filtered treand. ΔT filtered , which is defined recursively based on hourly calculated values of ΔT residual . The filtered value at any given hour is equal to 0.90 times the filtered value for the previous hour, plus 0.10 times the current value of ΔT residual : ΔT filtered t = 0.9 * ΔT filtered t − 1 + 0.1 * ΔT residual t .

[0132] Catalyst degradation was detected when ΔT filtered , increased by 6°C above the lowest value since the last detected refurbishment. Catalyst refurbishment was detected when ΔT filtered decreased by 4°C below the highest value since the last detected catalyst degradation.

[0133] The results of the validation test are shown in FIG. 4 along with Table 1. In 10 out of 12 results, the identified catalyst degradation events and catalyst refurbishment events identified using the ΔT filtered had corresponded to actual catalyst degradation or refurbishment events as determined using plant operating and maintenance logs.

[0134] The results of the validation test confirm the utility of the method and system for calculating an effective indicator of catalyst activity. TABLE 1ΔT filtered (°C)Start Date End Date Start Value End Value Change Detected Catalyst Event Actual Event / Condition 1 Jul 201019 Oct 2010-2.92.15.0DegradationCatalyst skimming was performed in Feb. 2011, hence this would have been preceded by a period of decreasing catalyst activity30 Jan 20111 Mar 20114.40.4-4.0RefurbishmentTop 40% of the catalyst was replaced in Feb. 201115 Mar 201123 Apr 2011-5.6-0.65.0DegradationA period of "break-in" degradation typically follows catalyst replacement28 Oct 201115 Mar 20123.7-0.3-4.0RefurbishmentSteaming of the catalyst was performed in Feb. 2012, because hot-band formation was observed in October 20114 Oct 20128 Nov 2012-8.9-3.95.0DegradationNo significant recorded observation17 Jun 20131 Jul 20130-4-4.0RefurbishmentNo significant recorded observation1 Jul 20131 Oct 2013-3.9715.0DegradationHot band formation observed in Oct. - Nov. 20136 Dec 20139 Mar 20141.7-2.3-4.0RefurbishmentA complete catalyst change-out was performed in Feb. 201419 Aug 20142 May 2015-3.21.85.0DegradationA complete catalyst change-out was performed in Feb. 2014, this would typically be followed by a period of degradation that one would expect after a new catalyst installation6 Jan 201618 Jan 20163-1-4.0RefurbishmentStarted operating at higher S / C ratios24 Jan 201627 Feb 2016-2.22.85.0DegradationSuspected hot band formation27 Feb 201610 May 20162.8-1.2-4.0RefurbishmentThe plant started operating at reduced rates and higher S / C in May 2016 due to fears of hot band degradation.

Claims

1. A method for determining changes in the catalytic activity of reforming catalyst, the method comprising: (a) introducing a reformer feed gas mixture (15) comprising at least one hydrocarbon and steam into one or more catalytic reactors (20) comprising the reforming catalyst in a reformer (10), each of the one or more catalytic reactors (20) having an inlet and an outlet, reacting the reformer feed gas mixture (15) in a steam-methane reforming reaction and water-gas shift reaction under reaction conditions effective to form a reformate (25) comprising H2, CO, CH4, and H2O, and withdrawing the reformate (25) from the one or more catalytic reactors (20); wherein for a plurality of times, the method comprises (b) measuring an outlet temperature, Toutlet, representative of a temperature at the outlet of the one or more catalytic reactors (20) for each time of the plurality of times; (c) determining a reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20), wherein this step comprises: - determining a hydrocarbon feed composition representative of a composition of the hydrocarbon feed via a composition measurement; - determining a flow rate of the hydrocarbon feed, ṅfeedstock, via a flow rate measurement of the hydrocarbon feed; - determining a flow rate of the steam, ṅsteam, via a flow rate measurement of the steam; - and calculating the reformate composition using chemical element flow rate balances, the representative hydrocarbon feed composition, the flow rate of the hydrocarbon feed, the flow rate of the steam, and a water-gas shift equilibrium constant evaluated at the measured outlet temperature, Toutlet; (d) measuring a reformate pressure, Poutlet, representative of a pressure at the outlet of the one or more catalytic reactors (20); (e) calculating an equilibrium temperature, Tequilibrium, from an equilibrium constant for the steam-methane reforming reaction at the reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20) and using the measured reformate pressure; (f) determining a temperature approach to equilibrium for the steam-methane reforming reaction at the outlet of the one or more catalytic reactors (20) for each time of the plurality of times using the measured outlet temperature, Toutlet, wherein the temperature approach, ΔTapproach, to equilibrium is a measure of a difference between the outlet temperature, Toutlet, and the equilibrium temperature, Tequilibrium; (g) calculating an empirical model-based temperature approach to equilibrium from reformer operating data in an empirical model based on historical operating data for a model reformer for each time of the plurality of times, wherein the model reformer is the reformer (10); (h) determining coefficients for the empirical model by regressing the empirical model-based temperature approach to equilibrium, ΔTempirical, historical, as a function of the historical outlet temperatures, Toutlet, historical, the historical reformate pressures representative of the pressure at the outlet of the one or more catalytic reactors, Poutlet, historical, the historical flow rate of the steam, ṅsteam, historical, preferably using a steam-to-carbon molar ratio, and the historical flow rate of the feedstock, ṅfeedstock, historical, the function being a linear or non-linear equation with the coefficients determined by regression of the historical operating data, wherein the empirical model-based temperature approach to equilibrium, ΔTempirical, historical, is a measure of a difference between historical outlet temperatures, Toutlet, historical, for the model reformer for a historical plurality of times, and calculated historical equilibrium temperatures, Tequilibrium, historical, for the model reformer for the historical plurality of times, ΔTempirical, historical being a measure consistent with ΔTapproach; and (i) comparing the temperature approach to equilibrium to the calculated empirical model-based temperature approach to equilibrium for each time of the plurality of times by calculating a characteristic operational value (ΔTresidual) from the temperature approach to equilibrium and the calculated empirical model-based temperature approach to equilibrium, the characteristic operational value including the difference and / or the ratio of the temperature approach to equilibrium and the calculated empirical model-based temperature approach to equilibrium; and (j) determining whether a change in the characteristic operational value (ΔTresidual) represents an objective increase or decrease in the catalytic activity of the reforming catalyst; wherein a reduced activity of the reforming catalyst is established when the characteristic operational value (ΔTresidual) increases by more than a selected amount for a selected period of time; wherein steps (b) through (j) are repeated for the plurality of times.

2. The method according to claim 1 further comprising: measuring a methane concentration representative of a concentration of methane in the reformate (25) withdrawn from the one or more catalytic reactors (20); and measuring a flow rate representative of a flow rate of the reformate (25); wherein the reformate composition is calculated also using the representative methane concentration; and wherein the reformate composition is calculated also using the representative flow rate of the reformate.

3. The method of claim 1 or claim 2 wherein steps (b) through (j) are repeated during a time period where the temperature approach to equilibrium relative to the empirical model-based temperature approach to equilibrium for each time of the plurality of times differs by an amount for a selected period of time, the amount and selected period of time determined to indicate decreased and / or increased activity of the reforming catalyst.

4. The method of the previous claim wherein the amount determined to indicate decreased activity of the reforming catalyst is an amount corresponding to an increase in a value of a residual temperature to equilibrium, ΔTresidual, of 3°C or more during the selected period of time, where ΔTresidual = ΔTapproach - ΔTempirical, where ΔTapproach = Toutlet - Tequilibrium, and corresponds to the temperature approach to equilibrium calculated in step (f), where Toutlet is the measured temperature representative of the temperature at the outlet of the one or more catalytic reactors (20), where Tequilibrium is a temperature calculated from an equilibrium constant for the steam-methane reforming reaction at a reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20), and where ΔTempirical corresponds to the calculated empirical model-based temperature approach to equilibrium for the regression where ΔTempirical, historical = Toutlet, historical - Tequilibrium, historical for the regression; and wherein the amount determined to indicate increased activity of the reforming catalyst is an amount corresponding to an decrease in a value of a residual temperature to equilibrium, ΔTresidual, of 3°C or more during the selected period of time, where ΔTresidual = ΔTapproach - ΔTempirical, where ΔTapproach = Toutlet - Tequilibrium, and corresponds to the temperature approach to equilibrium calculated in step (f), where Toutlet is the measured temperature representative of the temperature at the outlet of the one or more catalytic reactors (20), where Tequilibrium is a temperature calculated from an equilibrium constant for the steam-methane reforming reaction at a reformate composition representative of the reformate withdrawn from the one or more catalytic reactors (20), and where ΔTempirical corresponds to the calculated empirical model-based temperature approach to equilibrium for the regression where ΔTempirical, historical = Toutlet, historical - Tequilibrium, historical for the regression.