METHOD OF DETERMINING SUPERFICIAL GAS VELOCITY IN FLUIDIZED BED REACTORS

Abstract

Systems and methods useful in determining the superficial gas velocity in fluidized bed reactors may utilize a pressure drop across a portion of the system but not associated with a flowmeter. For example, method may comprise: obtaining a pressure for each of two different locations within a fluidized bed reactor system that comprises a reactor capable of containing a fluidized bed and a cycle gas loop, wherein one or both of the two different locations is not at a flowmeter; calculating a pressure drop based on the two pressures; calculating a first superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; and operating the fluidized bed reactor system based at least in part on the SGV.sub.alt.

Claims

1. A method comprising: obtaining a pressure measurement for each of two different locations within a fluidized bed reactor system that comprises a reactor capable of containing a fluidized bed and a cycle gas loop, wherein one or both of the two different locations is not at a flowmeter; calculating a pressure drop based on the two obtained pressure measurements; calculating a first superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; obtaining a second superficial gas velocity (SGV.sub.fm) for the fluidized bed using the flowmeter; comparing the SGV.sub.alt and the SGV.sub.fm; and operating the fluidized bed reactor system based at least in part on said comparing of the SGV.sub.alt and the SGV.sub.fm.

2. The method of claim 1, wherein operating the fluidized bed reactor based at least in part on said comparing of the SGV.sub.alt and the SGV.sub.fm comprises: modifying an operating parameter of the fluidized bed reactor system when the comparison of the SGV.sub.alt and the SGV.sub.fm fits a threshold requirement.

3. The method of claim 2, wherein the operating parameter comprises a reactor temperature, a reactor pressure, a reactant feed rate, a purge rate, reactor composition, or any combination thereof.

4. The method of claim 1, further comprising taking an action if the comparison of the SGV.sub.alt and the SGV.sub.fm fits an action threshold requirement.

5. The method of claim 4, wherein the action comprises one or more of the following: performing remedial actions on the flowmeter; triggering an alarm; and adding a flow bias or correction factor to the determined SGV.sub.fm value so as to control reactor velocity in a desired operating window.

6. The method of claim 1, wherein the comparison is an absolute value of a difference between the SGV.sub.alt and the SGV.sub.fm, and wherein the threshold requirement is the absolute value of the difference between the SGV.sub.alt and the SGV.sub.fm being 0.05 feet per second or greater.

7. The method of claim 1, wherein the comparison is an absolute value of a difference between the SGV.sub.alt and the SGV.sub.fm, and wherein the threshold requirement is the absolute value of the difference between the SGV.sub.alt and the SGV.sub.fm being 0.20 feet per second.

8. The method of claim 1, further comprising: producing a polyolefin in the fluidized bed reactor system.

9. The method of claim 1, wherein the flowmeter comprises a Venturi flowmeter.

10. The method of claim 1, wherein the flowmeter comprises a Pitot tube.

11. A method comprising: measuring a pressure for each of two different locations within a fluidized bed reactor system that comprises a reactor capable of containing a fluidized bed and a cycle gas loop, wherein one or both of the two different locations is not at a flowmeter; calculating a pressure drop based on the two measured pressures; calculating a superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; operating the fluidized bed reactor system based at least in part on the SGV.sub.alt; and producing a polyolefin in the fluidized bed reactor system.

12. The method of claim 11, wherein the operating of the fluidized bed reactor comprises: modifying an operating parameter of the fluidized bed reactor system.

13. The method of claim 12, wherein the operating parameter comprises a reactor temperature, a reactor pressure, a reactant feed rate, a purge rate, reactor composition, or any combination thereof.

14. The method of claim 11, wherein (a) the flowmeter comprises a Venturi flowmeter or (b) the flowmeter comprises a Pitot tube.

15. The method of claim 1, wherein the two different locations within the fluidized bed reactor system comprise a location (P1) at a top portion of the reactor and a location (P3) along a bottom portion of the fluidized bed within the reactor, such that obtaining the pressure measurement for each of said two different locations comprises obtaining a pressure drop across the fluidized bed; and further wherein SGV.sub.alt is calculated using at least one of the following methods: (a) a modified Ergun equation having the form of EQ. 1 $\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left(C_1\frac{\varphi }{{\left(1-\varphi \right)}^2}+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\right)& \mathrm{EQ}.1\end{array}$ wherein ΔP is the pressure drop across the fluidized bed (Pascals, Pa), h is the height between said two different locations P1 and P3 from which the pressure measurements are taken (meters, m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, ρ is fluid density (kg/m.sup.3), and C.sub.1 and C.sub.2 are dimensionless coefficients calculated per EQ. 2
C.sub.n=a.sub.nε.sup.b.sup.n+c.sub.n  EQ. 2 where n is the number of the coefficient, and a.sub.n, b.sub.n, and c.sub.n are additional constants associated with a given reactor size and type, and may be determined using superficial gas velocity data obtained from a flowmeter; (b) an equation having the form of EQ. 3 $\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left({\left(1-\varphi \right)}^{-3.65}\left(C_1+C_2{\mathrm{Re}}^{0.687}\right)\right)& \mathrm{EQ}.3\end{array}$ wherein each variable is as described in connection with EQ. 1; (c) an equation having the form of EQ. 4 $\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)*d^2}\left(C_1\left[\frac{\varphi }{{\left(1-\varphi \right)}^2}+{\left(1-\varphi \right)}^2\left(1+1.5{\varphi }^{\frac{1}{2}}\right)\right]+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\left[\frac{{\left(1-\varphi \right)}^{-1}+3\varphi \left(1-\varphi \right)+8.4{\mathrm{Re}}^{-0.343}}{1+{10}^{3\varphi }{\mathrm{Re}}^{-\frac{\left(1+4\varphi \right)}{2}}}\right]\right)& \mathrm{EQ}.4\end{array}$ wherein each variable is as described in connection with EQ. 1; and (d) an equation having the form of EQ. 5 $\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)*d^2}*\left[{\left(1-\varphi \right)}^{-3.65}\left(C_1+C_2{\mathrm{Re}}^{0.687}\right)\right]*\left(1-\mathrm{fH}\right)& \mathrm{EQ}.5\end{array}$ wherein ΔP, h, μ, U, d, Re, C.sub.1 and C.sub.2 are each as described in connection with EQ. 1, f is determined from EQ. 6 $\begin{array}{cc}f=\frac{{\mathrm{Fr}}^{-1.6}}{{\mathrm{Fr}}^{-1.6}+0.4},& \mathrm{EQ}.6\end{array}$ wherein Fr is Froude number where Fr=v.sub.t.sup.2(gL), v.sub.t is single-particle terminal velocity (m/sec), g is gravitational acceleration (m/s.sup.2), and L is reactor diameter (m), and H is determined per EQ. 7 $\begin{array}{cc}H=\{\begin{array}{cc}2.7{\varphi }^{0.234},& \varphi <0.0012\\ -0.019{\varphi }^{-0.455}+0.963,& 0.0012\le \varphi <0.014\\ 0.868e^{-0.38\varphi }-0.176e^{-119.2\varphi },& 0.014\le \varphi <0.25\\ -\frac{4.59}{{10}^5}{e}^{19.75\varphi }+0.852e^{-0.268\varphi },& 0.25\le \varphi <0.455\\ \begin{array}{c}\left(\varphi -0.59\right)(-1501{\varphi }^3+2203{\varphi }^2-\\ 1054\varphi +162),\end{array}& 0.455\le \varphi <0.59\\ 0,& \varphi \ge 0.59\end{array}& \mathrm{EQ}.7\end{array}$

16. The method of claim 1, wherein SGV.sub.alt is calculated using an equation having the form of EQ. 8
dP=0.5ρU.sub.0.sup.2K.sub.overall  EQ. 8 where dP is the pressure drop between the two different locations from which pressure measurements are obtained (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U.sub.0 is the superficial gas velocity (m/sec), and K.sub.overall is a combination of the individual pressure drop correction factors needed to account for the pressure drop between said two different locations from which pressure measurements are obtained.

17. The method of claim 16, wherein K.sub.overall is calculated using an empirical correlation determined from data of past operation of the reactor system over a time period during which accuracy of the flow meter is not in question.

18. The method of claim 16, wherein obtaining the pressure measurement for each of the two different locations within the reactor system comprises obtaining a pressure drop across reactor inlet piping of the reactor system; and further wherein K.sub.overall is determined as K.sub.overall=ΣK.sub.i, wherein each Ki is a pressure drop correction factor associated with a location or feature along the reactor inlet piping that causes pressure drop, and further wherein each K.sub.i is selected from the following: K.sub.friction, K.sub.elbow, K.sub.exit, K.sub.disc, K.sub.statichead, and any combination thereof.

19. The method of claim 16, wherein obtaining the pressure measurement for each of the two different locations within the reactor system comprises obtaining a pressure drop across a distributor plate disposed in a bottom portion of the reactor, and further wherein K.sub.overall is determined as the K value of the distributor plate.

20. The method of claim 1, wherein the two different locations within the fluidized bed reactor system have a distributor plate therebetween, and wherein SGV.sub.alt is calculated using an equation having the form of EQ. 14 $\begin{array}{cc}\mathrm{dP}=0.5\rho U^2{K}_{\mathrm{overall}}+h*\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left(C_1\frac{\varphi }{{\left(1-\varphi \right)}^2}+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\right)& \mathrm{EQ}.14\end{array}$ where dP is the pressure drop between the two different locations from which pressure measurements are obtained (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U is the superficial gas velocity (m/sec), K.sub.overall is a combination of the individual pressure drop correction factors needed to account for the pressure drop between said two different locations from which pressure measurements are obtained, h is the height from the top of the distributor the location of the two locations that is downstream of the distributor plate, μ is fluid viscosity in (Pa*sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3)
C.sub.n=a.sub.nε.sup.b.sup.n+c.sub.n  EQ. 2 where n is the number of the coefficient, and a.sub.n, b.sub.n, and c.sub.n are additional constants associated with a given reactor size and type, and may be determined using superficial gas velocity data obtained from a flowmeter; and further wherein the two different locations within the fluidized bed reactor system also have therebetween one or both of: reactor inlet piping and a disc; and further wherein K.sub.overall is determined as K.sub.overall=ΣK.sub.i, wherein each Ki is a pressure drop correction factor associated with a location or feature along the reactor inlet piping that causes pressure drop, and further wherein each K.sub.i is selected from the following: K.sub.plate, K.sub.friction, K.sub.elbow, K.sub.exit, K.sub.disc, and any combination thereof.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0013] To assist those of ordinary skill in the relevant art in making and using the subject matter hereof, reference is made to the appended drawings. The following figures are included to illustrate certain aspects of the disclosure, and should not be viewed as exclusive configurations. The subject matter disclosed is capable of considerable modifications, alterations, combinations, and equivalents in form and function, as will occur to those skilled in the art and having the benefit of this disclosure.

[0014] FIG. 1 is a diagram of a nonlimiting portion of a fluidized bed reactor system of the present disclosure.

[0015] FIG. 2 is a diagram of a nonlimiting portion of a fluidized bed reactor system of the present disclosure

[0016] FIG. 3 illustrates a comparison between the superficial gas velocities determined via the modified Ergun equation and the Venturi flowmeter measurement for all metallocene grades.

[0017] FIG. 4 illustrates a comparison between the superficial gas velocities determined via the modified Ergun equation and the Venturi flowmeter measurement for a single metallocene grades.

[0018] FIG. 5 is plant data used for determining an empirical correlation for K.sub.overall.

[0019] FIG. 6 is a graph of SGV.sub.fm vs. SGV.sub.ΔP for a commercial polyolefin reactor over a 3.5 month period.

DETAILED DESCRIPTION

[0020] The present application relates to systems and methods for determining the superficial gas velocity in fluidized bed reactors, for example, polyolefin reactors. More specifically, the methods and systems described herein involve determining the superficial gas velocity in a fluidized bed using one or more pressure drops within the fluidized bed reactor system, including pressure drops in the fluidized bed, sections of the cycle gas piping, or both, but not pressure drops across a flowmeter (e.g., a Venturi flowmeter, a Pitot tube, or the like) or related instrument.

[0021] As used herein, the abbreviation “SGV.sub.alt” is used generally to refer to a superficial gas velocity calculated or otherwise obtained using a pressure drop within the fluidized bed reactor system where said pressure drop is not derived from a flowmeter (e.g., a Venturi flowmeter, a Pitot tube, or the like). As used herein, the abbreviation “SGV.sub.fm” is used generally to refer to a superficial gas velocity calculated or otherwise obtained based on flowmeter measurements (e.g., a Venturi flowmeter, a Pitot tube, or the like) from within the fluidized bed reactor system.

[0022] Using the pressure drop within the fluidized bed reactor system, the superficial gas velocity (SGV.sub.alt) may be calculated based on said pressure drop. The specific choice of equation used to estimate the superficial gas velocity depends on the two locations across which the pressure drop is measured or otherwise obtained.

[0023] The pressure drop measurements for the SGV.sub.alt rely on pressure sensors (e.g., differential pressure sensors) that are typically more robust and less prone to fouling, drift, and other sources of inaccuracy as compared to the Venturi flowmeters. Accordingly, the SGV.sub.alt may be more accurate and provide for a more confident monitoring of the fluidized bed reactor as compared to the SGV.sub.fm.

[0024] The methods and systems described herein provide multiple avenues for implementing the SGV.sub.alt with or without the SGV.sub.fm. For example, the SGV.sub.fm and the SGV.sub.alt may be used together, where comparisons between the superficial gas velocities may identify drift in the SGV.sub.fm. Identifying drift, or other inaccuracies, early may mitigate either system downtime or production losses caused by operating the system based on the inaccurate superficial gas velocity. Identification of SGV.sub.fm drift may, for example, be achieved by implementing deviation alarming between SGV.sub.fm and SGV.sub.alt in the plant's process control system.

[0025] FIG. 1 is a diagram of a nonlimiting fluidized bed reactor system 100 of the present disclosure. The fluidized bed reactor system 100 includes a fluidized bed reactor 102 (or simply reactor 102) for containing a fluidized bed 104 above a distributor plate 106 located in the reactor 102. The fluidized bed reactor system 100 also includes a cycle gas loop 108 where gas effluent from the reactor 102 is compressed via compressor 110 and cooled via heat exchanger 112. As illustrated, a line 114 is used for conveying the gas effluent from the reactor 102 to the compressor 110, a line 116 is used for conveying the compressed gas effluent to the heat exchanger 112, and a line 118 is used for conveying the resultant recycle gas (the cooled and compressed gas effluent) to the reactor 102. The lines 114, 116, 118 provide simple flow or conveyance illustrations and may, in reality, be composed of multiple lines and have equipment (e.g., sensors, pumps, and the like) located along said lines 114, 116, 118.

[0026] In the illustrated fluidized bed reactor system 100, feedstock is added to the cycle gas loop 108 via line 120 downstream of the heat exchanger 112 along line 118. Further, a Venturi flowmeter 124 or similar hardware is located upstream of the compressor 110 along line 114. The fluidized bed reactor system 100 further includes a disc 122 within the reactor 102 for distributing as the feed coming into the reactor 102 from line 118 at inlet 124.

[0027] The systems and methods of the present disclosure herein use pressure drop values across two locations in the system (e.g., a portion of the fluidized bed reactor where said portion may be being defined between two locations within the fluidized bed reactor system) for calculating SGV.sub.alt, as opposed to across a flowmeter as done with the Venturi 124. The two locations across which the pressure drop is defined (e.g., the portion of the fluidized bed reactor system from which the pressure drop is obtained) dictates on the equation used for calculating the SGV.sub.alt.

[0028] Referring back to FIG. 1, six locations P1-P6 in the fluidized bed reactor system 100 are illustrated as locations that can be used to define the portion of the fluidized bed reactor system from which the pressure drop is obtained. That is, the pressure drop may be for a portion of the fluidized bed reactor system defined between P1 and P2, between P1 and P3, between P1 and P4, between P1 and P5, between P1 and P6, between P2 and P3, between P2 and P4, between P2 and P5, between P2 and P6, between P3 and P4, between P3 and P5, between P3 and P6, between P4 and P5, between P4 and P6, or between P5 and P6. The first location P1 is the reactor pressure measured above the fluidized bed 104. The second location P2 is the reactor pressure between the distributor plate 106 and disc 122. The third location P3 is in the reactor above the distributor plate 106 and within the bottom 10 vol % of the fluidized bed portion 104 of the reactor 102. The fourth location P4 is the middle 80 vol % of the fluidized bed portion 104 of the reactor 102. The fifth location P5 is the top 10 vol % of the fluidized bed portion 104 of the reactor 102. The sixth location P6 is in the line 118 after feedstock is added via line 120 and before an elbow 126 closest to the inlet 124 of the reactor 102.

[0029] The pressure at each of these locations P1-P6 may be measured using (a) individual pressure sensor at each location or (b) a differential pressure sensor connected to the two locations. In either instance, taps (not illustrated) may be included in the fluidized bed reactor system 100 for connecting either type of pressure sensors. Examples of pressure sensors may include, but are not limited to, pressure transducers, viscometers, thermocouples, the like, and any combination thereof. The types of pressure sensors utilized may include, but are not limited to, differential, acoustic, electromagnetic, piezoelectric, optical, potentiometric, the like, or any combination thereof. Examples of commercially available pressure sensors may include, but are not limited to, the DYTRAN 2006M1 acoustic pressure sensor, the OMEGA PX3005 piezoelectric pressure sensor, the EMERSON ROSEMOUNT 3051 pressure sensors, and the HONEYWELL FP2000 potentiometric pressure sensor.

[0030] Table 1 provides a list of example equations for calculating the SGV.sub.alt and the corresponding pressure drop locations within a fluidized bed reactor system used as an input for said equation. Each of the equations is described in more detail herein.

TABLE-US-00001 TABLE 1 Examples of SGV.sub.alt equations and corresponding pressure drop locations within a fluidized bed reactor system Equation Pressure Drop Location SGV.sub.alt Method A Across the fluidized bed (e.g., P3 minus P1 of (modified Ergun model) FIG. 1) SGV.sub.alt Method B Across the fluidized bed (e.g., P3 minus P1 of (Wen-Yu model) FIG. 1) SGV.sub.alt Method C Across the fluidized bed (e.g., P3 minus P1 of (Beetstra model) FIG. 1) SGV.sub.alt Method D Across the fluidized bed (e.g., P3 minus P1 of (Igci model) FIG. 1) SGV.sub.alt Method E Across the reactor inlet piping (e.g., P6 minus P2 of FIG. 1) SGV.sub.alt Method F Across the distributor plate (e.g., P3 minus P2 of FIG. 1) SGV.sub.alt Method G Across the fluidized bed and distributor plate (e.g., P1 minus P2 of FIG. 1) SGV.sub.alt Method H From the reactor inlet piping to above the fluidized bed (e.g., P1 minus P6 of FIG. 1)

SGV.SUB.alt .Method a (Modified Ergun Equation)

[0031] A first example SGV.sub.alt method is a modified Ergun equation presented in EQ. 1.

[00001]$\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left(C_1\frac{\varphi }{{\left(1-\varphi \right)}^2}+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\right)& \mathrm{EQ}.1\end{array}$

where ΔP is the pressure drop across the fluidized bed (P3 minus P1 of FIG. 1) (Pascals, Pa), h is the height between the two locations across which the pressure drop is measured (e.g., the height between the pressure taps used for measuring the differential pressure between P3 and P1) (meters, m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3).

[0032] Values like gas density and viscosity may change during operation and can be accurately predicted from temperature, pressure, and composition data, which may each be measured or otherwise obtained in real time and/or periodically. Alternatively, such values that change during operation may be treated as constants, which reduces the accuracy and cost of the method. A combination of the foregoing may be implemented where some values are constants and other are measured or otherwise obtained in real time and/or periodically.

[0033] The coefficients C.sub.1 and C.sub.2 may have the form of EQ. 2.

C.sub.n=a.sub.nε.sup.b.sup.n+c.sub.n  EQ. 2

where n is the number of the coefficient (e.g., 1 for C.sub.1), and a.sub.n, b.sub.n, and c.sub.n are additional constants.

[0034] C.sub.1 and C.sub.2 of EQ. 1 and a, b, and c of EQ. 2 may be unique for each reactor size and type. The constants C.sub.1 and C.sub.2 of EQ. 1 and/or a, b, and c of EQ. 2 may be determined using experimental superficial gas velocity data obtained from a Venturi flowmeter or other suitable means. Preferably, the data used to derive the constants has a high accuracy (e.g., high certainty that the Venturi flowmeter was not drifting or fouled). As additional data is collected during the method and/or operation of the systems described herein, the constants C.sub.1 and C.sub.2 of EQ. 1 and/or a, b, and c of EQ. 2 may be updated using superficial gas velocity data acquired independently from the modified Ergun equation. In this way, the determination of the constants using the Venturi flowmeter or other suitable means may be considered a feedback loop that allows for constant refining of the modified Ergun equation to provide higher accuracy and confidence in the superficial gas velocity data.

[0035] Once the constants are determined and/or provided from a previous determination and the properties of the fluidized bed (e.g., length, porosity, and density) determined, EQ. 1 has two variables: pressure drop across the fluidized bed and superficial gas velocity.

[0036] The systems described herein may include pressure sensors for measuring the pressure drop across the fluidized bed so that solving EQ. 1 (the modified Ergun equation) yields the SGV.sub.alt. The methods described herein may include measuring the pressure drop across the fluidized bed of a fluidized bed reactor system and solving EQ. 1 to yield a SGV.sub.alt. Said SGV.sub.alt may be used, as described further herein, for operating the fluidized bed reactor system alone or in conjunction with the SGV.sub.fm.

SGV.SUB.alt .Method B

[0037] A second example SGV.sub.alt method is presented in EQ. 3.

[00002]$\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left({\left(1-\varphi \right)}^{-3.65}\left(C_1+C_2{\mathrm{Re}}^{0.687}\right)\right)& \mathrm{EQ}.3\end{array}$

where ΔP is the pressure drop across the fluidized bed (P3 minus P1 of FIG. 1) (Pa), h is the height between the two locations across which the pressure drop is measured (e.g., the height between the pressure taps used for measuring the differential pressure between P3 and P1) (m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3).

SGV.SUB.alt .Method C

[0038] A third example SGV.sub.alt method is presented in EQ. 4.

[00003]$\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)*d^2}\left(C_1\left[\frac{\varphi }{{\left(1-\varphi \right)}^2}+{\left(1-\varphi \right)}^2\left(1+1.5{\varphi }^{\frac{1}{2}}\right)\right]+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\left[\frac{{\left(1-\varphi \right)}^{-1}+3\varphi \left(1-\varphi \right)+8.4{\mathrm{Re}}^{-0.343}}{1+{10}^{3\varphi }{\mathrm{Re}}^{-\frac{\left(1+4\varphi \right)}{2}}}\right]\right)& \mathrm{EQ}.4\end{array}$

where ΔP is the pressure drop across the fluidized bed (P3 minus P1 of FIG. 1) (Pa), h is the height between the two locations across which the pressure drop is measured (e.g., the height between the pressure taps used for measuring the differential pressure between P3 and P1) (m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3).

SGV.SUB.alt .Method D

[0039] A fourth example SGV.sub.alt method is presented in EQS. 5-7.

[00004]$\begin{array}{cc}\frac{\Delta P}{h}=\frac{18\mu U\varphi }{\left(1-\varphi \right)*d^2}*\left[{\left(1-\varphi \right)}^{-3.65}\left(C_1+C_2{\mathrm{Re}}^{0.687}\right)\right]*\left(1-\mathrm{fH}\right)& \mathrm{EQ}.5\end{array}$$\begin{array}{cc}f=\frac{{\mathrm{Fr}}^{-1.6}}{{\mathrm{Fr}}^{-1.6}+0.4},& \mathrm{EQ}.6\end{array}$$\begin{array}{cc}H=\{\begin{array}{cc}2.7{\varphi }^{0.234},& \varphi <0.0012\\ -0.019{\varphi }^{-0.455}+0.963,& 0.0012\le \varphi <0.014\\ 0.868e^{-0.38\varphi }-0.176e^{-119.2\varphi },& 0.014\le \varphi <0.25\\ -\frac{4.59}{{10}^5}{e}^{19.75\varphi }+0.852e^{-0.268\varphi },& 0.25\le \varphi <0.455\\ \begin{array}{c}\left(\varphi -0.59\right)(-1501{\varphi }^3+2203{\varphi }^2-\\ 1054\varphi +162),\end{array}& 0.455\le \varphi <0.59\\ 0,& \varphi \ge 0.59\end{array}& \mathrm{EQ}.7\end{array}$

where ΔP is the pressure drop across the fluidized bed (P3 minus P1 of FIG. 1) (Pa), h is the height between the two locations across which the pressure drop is measured (e.g., the height between the pressure taps used for measuring the differential pressure between P3 and P1) (m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, ρ is fluid density (kg/m.sup.3), Fr is Froude number where Fr=v.sub.t.sup.2/(gL), v.sub.t is single-particle terminal velocity (m/sec), g is gravitational acceleration (m/s.sup.2), and L is resolution length scale of the model (also known as reactor diameter) (m).

SGV.SUB.alt .Method E

[0040] A sixth example SGV.sub.alt method is presented in EQS. 8-9.

dP=0.5ρU.sup.2K.sub.overall  EQ. 8

K.sub.overall=K.sub.elbow+K.sub.exit+K.sub.disc+K.sub.static  EQ. 9

where dP is the pressure drop across the reactor inlet piping (P6 minus P2 of FIG. 1) (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U is the superficial gas velocity (m/sec), and K.sub.overall or K factor is the sum of the pressure drop loss coefficients (ΣK.sub.i) between where each of the locations defining the pressure drop (e.g., between P6 and P2 of FIG. 1), which may include, but is not limited to, K.sub.friction, K.sub.elbow, K.sub.exit, K.sub.disc, K.sub.statichead, and any combination thereof. EQ. 9 for K.sub.overall is a nonlimiting example equation that can be changed to account for the configuration of the system.

[0041] The K factor depends on the configuration and dimensions of the reactor inlet piping (horizontal and vertical), number and type of pipefittings, or other internal components that would result in a pressure drop in the system between points P6 and P2. In summary, K.sub.overall is a combination of the individual pressure drop correction factors needed to account for the pressure drop across the reactor inlet piping line. FIG. 2 is a diagram of a nonlimiting example reactor inlet piping ΔP instrument connections illustrating where correction factors may be accounted for. In the illustrated portion of a fluidized bed reactor system, a line 206 provides gas feed (e.g., recycled gas optionally augmented with fresh reactant gas) to a fluidized bed reactor 202 capable of containing a fluidized bed 204 above a distributor plate 208, where the line 206 connects to the reactor 202 (e.g., at an inlet) and within the reactor 202 is illustrated a disc 210 that assists with distributing the gas and liquid feed into the reactor 202 to prevent liquid pooling. The system also includes tap 212 that connects the bottom portion of the reactor 202 between the distributor plate 208 and disc 210 to a differential pressure meter 216 and tap 214 that connects the line 206 before the elbow to the differential pressure meter 216.

[0042] For calculating SGV.sub.alt according to EQ. 8 for the illustrated configuration in FIG. 2 using a pressure drop of P6 minus P2, correction factors for friction pressure drop in the line 206 (K.sub.friction) between P6 and the elbow, for the elbow in the line 206 (K.sub.elbow), for the connection between line 206 and the reactor (K.sub.exit), for the disc 210 (K.sub.disc), and for the vertical distance of the bypass line (K.sub.statichead) should be considered. In this example, K.sub.overall is the summation of the five foregoing Ks. K.sub.overall may be determined using plant data to fit EQ. 8. Alternatively, K.sub.overall may be calculated.

[0043] For two-phase flow, known as condense mode or super-condense mode operation in a fluidized bed process, where the wt % cond frac or wt % condensing at the reactor inlet is >0 and less than 50%, but not limited to 50%, the K.sub.overall will vary as a function of wt % condensing. An empirical correlation (EQ. 11) for K.sub.overall can be developed using plant data (FIG. 5). The data used for the correlation needs to be taken over a period of time when the Venturi flow meter is not fouled or accuracy of the flow meter is not in question. K.sub.overall as shown by EQ. 10 can be calculated by re-arranging EQ. 11 and solving for K.sub.overall. FIG. 6 shows SGV.sub.fm vs. SGV.sub.alt (top traces) for a commercial polyolefin reactor over a 3.5 month period. The ratio of SGV.sub.fm to SGV.sub.alt (bottom trace) shows SGV.sub.alt to be within 3% of the Venturi SGV.

[00005]$\begin{array}{cc}{K}_{\mathrm{overall}}={\left(\frac{2g*144*\Delta P}{{\rho }_i}\right)\left[\frac{1}{v_0*\left(1-\mathrm{wt}\%\mathrm{cond}\mathrm{frac}\right)}*\left(\frac{{\rho }_i{A}_i}{{\rho }_o{A}_o}\right)\right]}^2)& \mathrm{EQ}.10\end{array}$$\begin{array}{cc}{K}_{\mathrm{Overall}}=\alpha +\beta *\mathrm{Wt}\%\mathrm{Cond}\mathrm{Fraction}& \mathrm{EQ}.11\end{array}$

SGV.SUB.alt .Method F

[0044] A fifth example SGV.sub.alt method is presented in EQS. 12-13.

dP=0.5ρU.sup.2K.sub.overall  EQ. 12

K.sub.overall=K.sub.plate  EQ. 13

where dP is the pressure drop across the distributor plate (P3 minus P2 of FIG. 1) (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U is the superficial gas velocity (m/sec), and K.sub.overall or K factor is the sum of the pressure drop loss coefficients (ΣK.sub.i) between where each of the locations defining the pressure drop (e.g., between P3 and P2 of FIG. 1), which may include, but is not limited to, K.sub.plate. EQ. 13 for K.sub.overall is a nonlimiting example equation that can be changed to account for the configuration of the system. K.sub.plate is determined as described for the K values described in Method E.

SGV.SUB.alt .Method G

[0045] A seventh example SGV.sub.alt method is presented in EQS. 14-15.

[00006]$\begin{array}{cc}\mathrm{dP}=0.5\rho U^2{K}_{\mathrm{overall}}+h*\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left(C_1\frac{\varphi }{{\left(1-\varphi \right)}^2}+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\right)& \mathrm{EQ}.14\end{array}$$\begin{array}{cc}{K}_{\mathrm{overall}}=K_{\mathrm{plate}}& \mathrm{EQ}.15\end{array}$

where dP is the pressure drop across the fluidized bed and distributor plate (P1 minus P2 of FIG. 1) (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U is the superficial gas velocity (m/sec), K.sub.overall or K factor is the sum of the pressure drop loss coefficients (ΣK.sub.i) between where each of the locations defining the pressure drop (e.g., between P1 and P2 of FIG. 1), which may include, but is not limited to, K.sub.plate, h is the height from the top of the distributor plate to P1 (m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3). EQ. 15 for K.sub.overall is a nonlimiting example equation that can be changed to account for the configuration of the system. K.sub.plate is determined as described for the K values described in Method E.

SGV.SUB.alt .Method H

[0046] An eighth example SGV.sub.alt method is presented in EQS. 16-17.

[00007]$\begin{array}{cc}\mathrm{dP}=0.5\rho U^2{K}_{\mathrm{overall}}+h*\frac{18\mu U\varphi }{\left(1-\varphi \right)d^2}*\left(C_1\frac{\varphi }{{\left(1-\varphi \right)}^2}+C_2\frac{\mathrm{Re}}{{\left(1-\varphi \right)}^2}\right)& \mathrm{EQ}.16\end{array}$$\begin{array}{cc}{K}_{\mathrm{overall}}=K_{\mathrm{plate}}+K_{\mathrm{elbow}}+K_{\mathrm{exit}}+K_{\mathrm{disc}}+K_{\mathrm{statichead}}& \mathrm{EQ}.17\end{array}$

where dP is the pressure drop across the fluidized bed and distributor plate (P1 minus P6 of FIG. 1) (Pa), ρ is the gas density at the reactor inlet (kg/m.sup.3), U is the superficial gas velocity (m/sec), and K.sub.overall or K factor is the sum of the pressure drop loss coefficients (ΣK.sub.i) between where each of the locations defining the pressure drop (e.g., between P1 and P6 of FIG. 1), which may include, but is not limited to, K.sub.plate, K.sub.friction, K.sub.elbow, K.sub.exit, K.sub.disc, K.sub.statichead, and any combination thereof, h is the height from the top of the distributor plate to P1 (m), μ is fluid viscosity in (Pa*sec), U is superficial gas velocity (m/sec), ϕ is solids volume fraction (dimensionless), d is particle size (m), C.sub.1 and C.sub.2 are coefficients (dimensionless, calculated per EQ. 2), Re is Reynolds number where Re=ρUd/μ, and ρ is fluid density (kg/m.sup.3). EQ. 17 for K overall is a nonlimiting example equation that can be changed to account for the configuration of the system. K.sub.plate is determined as described for the K values described in Method E.

Methods and Systems Implementing a SGV.SUB.alt

[0047] The methods and systems described herein may use a pressure drop within a portion of a fluidized bed reactor system for determining SGV.sub.alt and may optionally further include a Venturi flowmeter for determining SGV.sub.fm. The SGV.sub.alt optionally in conjunction with the SGV.sub.fm may be used during the operation of the fluidized bed reactor system.

[0048] For example, a method of the present disclosure may comprise: obtaining a pressure for each of two different locations within a fluidized bed reactor system that comprises a reactor capable of containing a fluidized bed and a cycle gas loop; calculating a pressure drop based on the two pressures; calculating a superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; and operating the fluidized bed reactor system based at least in part on the SGV.sub.alt. Operating the fluidized bed reactor system may include maintaining or changing one or more operating parameters like reactor temperature, reactor pressure, reactant feed rate, purge rate, reactor composition, the like, and any combination thereof. For example, in reactor composition, the actual composition and/or concentrations of reactants and/or induced condensing agents may be adjusted to effect a change in the reactor composition.

[0049] In another example, a method of the present disclosure may comprise: obtaining a pressure for each of two different locations within a fluidized bed reactor system that comprises a reactor capable of containing a fluidized bed and a cycle gas loop; calculating a pressure drop based on the two pressures; calculating a first superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; obtaining a second superficial gas velocity (SGV.sub.fm) for the fluidized bed using a Venturi flowmeter; comparing the SGV.sub.alt and the SGV.sub.fm; and operating the fluidized bed reactor system based at least in part on the comparison of the SGV.sub.alt and the SGV.sub.fm (for instance, when the comparison fits a threshold requirement, one or more operating parameters such as those just discussed may be changed, which may help adjust to an unexpected or unacceptable deviation between the SGV values). The method can further include taking an action if the comparison of the SGV.sub.alt and the SGV.sub.fm is beyond the threshold value.

[0050] Example of such actions may include, but are not limited to, triggering an alarm, performing remedial actions on the Venturi flowmeter (e.g., performing a blowback procedure on the Venturi flowmeter), modifying one or more operating parameters of the fluidized bed reactor system, the like, and any combination thereof. Another action could include adding a flow bias or correction factor to the indicated Venturi flow so that the reactor velocity can be controlled in the correct operating window until the issue with the Venturi flow deviation could be resolved by cleaning the leads or throat of the Venturi flow meter during a planned reactor outage. It will be appreciated that these “actions” as just discussed do not necessarily involve an actual change to operating parameters, but instead may only entail adjustments to control algorithms and/or to measurements, thereby achieving higher confidence in control parameters before problems in actual operation occur.

[0051] Regarding the threshold requirement, when using the SGV.sub.fm in tandem with the SGV.sub.alt, a difference between the two values as low as 0.05 ft/s may, for example, be cause for an action to be taken (similarly, a difference as low as 0.10, 0.125, 0.15, 0.175, or 0.2 ft/s may be cause for an action to be taken). Further, a difference of 0.20 ft/s or more, such as 0.22 ft/s or more or 0.25 ft/s or more, between the SGV.sub.alt and the SGV.sub.fm may, for example, indicate a partially fouled Venturi flowmeter or another issue that may require immediate attention. Of course, the skilled artisan with the benefit of this disclosure could set any desired threshold difference to trigger (or indicate the possible desirability of) a given action, and/or a change in process conditions. And, as just alluded to, it is contemplated that more than one thresholds could be set, each linked to different action and/or change in process conditions (following the example just noted, a first threshold could be a difference greater than or equal to 0.050 ft/s; and a second threshold could be a difference greater than or equal to 0.200 ft/s; and each threshold could be used to trigger a different action and/or change of operating parameters). Again, any number of thresholds at any difference levels may readily be employed by the skilled artisan with the benefit of this disclosure; all such threshold numbers and amounts would fall within the scope of the present invention.

[0052] The fluidized bed reactor systems described herein may be suitable for producing a polyolefin (e.g., polyethylene, polypropylene, polybutylene, and the like). For example, methods described herein may further comprise producing a polyolefin in the fluidized bed reactor system.

[0053] Various aspects of the systems and methods described herein utilize computer systems. Such systems and methods can include a non-transitory computer readable medium containing instructions that, when implemented, cause one or more processors to carry out the methods described herein.

[0054] “Computer-readable medium” or “non-transitory, computer-readable medium,” as used herein, refers to any non-transitory storage and/or transmission medium that participates in providing instructions to a processor for execution. Such a medium may include, but is not limited to, non-volatile media and volatile media. Non-volatile media includes, for example, NVRAM, or magnetic or optical disks. Volatile media includes dynamic memory, such as main memory. Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, a hard disk, an array of hard disks, a magnetic tape, or any other magnetic medium, magneto-optical medium, a CD-ROM, a holographic medium, any other optical medium, a RAM, a PROM, and EPROM, a FLASH-EPROM, a solid state medium like a memory card, any other memory chip or cartridge, or any other tangible medium from which a computer can read data or instructions. When the computer-readable media is configured as a database, it is to be understood that the database may be any type of database, such as relational, hierarchical, object-oriented, and/or the like. Accordingly, exemplary embodiments of the present systems and methods may be considered to include a tangible storage medium or tangible distribution medium and prior art-recognized equivalents and successor media, in which the software implementations embodying the present techniques are stored.

[0055] The methods described herein can be performed using computing devices or processor-based devices that include a processor; a memory coupled to the processor; and instructions provided to the memory, wherein the instructions are executable by the processor to perform the methods described herein. The instructions can be a portion of code on a non-transitory computer readable medium. Any suitable processor-based device may be utilized for implementing all or a portion of embodiments of the present techniques, including without limitation personal computers, networks of personal computers, laptop computers, computer workstations, mobile devices, multi-processor servers or workstations with (or without) shared memory, high performance computers, and the like. Moreover, embodiments may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits.

[0056] For example, a fluidized bed reactor system may comprise: a reactor capable of containing a fluidized bed; a cycle gas loop; two or more pressure sensors; and a computing device configured to receive pressure data from the two or more pressure sensors and comprising: a processor; a memory coupled to the processor; and instructions provided to the memory, wherein the instructions are executable by the processor to perform a method comprising: calculate a pressure drop based on the pressure data; and calculate a first superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop. The method may also include send a signal to a portion of the fluidized bed reactor system, wherein the signal is a change in an operating parameter of said portion of the fluidized bed reactor system.

[0057] In another example, a fluidized bed reactor system may comprise: a reactor capable of containing a fluidized bed; a cycle gas loop; two or more pressure sensors; a Venturi flowmeter; and a computing device (a) configured to receive pressure data from the two or more pressure sensors and received data from the Venturi flowmeter; and (b) comprising: a processor; a memory coupled to the processor; and instructions provided to the memory, wherein the instructions are executable by the processor to perform a method comprising: calculate a pressure drop based on the pressure data; calculate a first superficial gas velocity (SGV.sub.alt) for the fluidized bed based on the pressure drop; obtain a second superficial gas velocity (SGV.sub.fm) based on the data from the Venturi flowmeter; compare the SGV.sub.alt and the SGV.sub.fm; and cause an action to occur if the comparison of the SGV.sub.alt and the SGV.sub.fm fits a threshold requirement (e.g., by sending a signal to a portion of the fluidized bed reactor system and/or outputting an alarm).

[0058] Unless otherwise indicated, all numbers expressing quantities of ingredients, properties such as molecular weight, reaction conditions, and so forth used in the present specification and associated claims are to be understood as being modified in all instances by the term “about.” Accordingly, unless indicated to the contrary, the numerical parameters set forth in the following specification and attached claims are approximations that may vary depending upon the desired properties sought to be obtained by the incarnations of the present inventions. At the very least, and not as an attempt to limit the application of the doctrine of equivalents to the scope of the claim, each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.

[0059] One or more illustrative incarnations incorporating one or more invention elements are presented herein. Not all features of a physical implementation are described or shown in this application for the sake of clarity. It is understood that in the development of a physical embodiment incorporating one or more elements of the present invention, numerous implementation-specific decisions must be made to achieve the developer's goals, such as compliance with system-related, business-related, government-related and other constraints, which vary by implementation and from time to time. While a developer's efforts might be time-consuming, such efforts would be, nevertheless, a routine undertaking for those of ordinary skill in the art and having benefit of this disclosure.

[0060] While compositions and methods are described herein in terms of “comprising” various components or steps, the compositions and methods can also “consist essentially of” or “consist of” the various components and steps.

[0061] To facilitate a better understanding of the embodiments of the present invention, the following examples of preferred or representative embodiments are given. In no way should the following examples be read to limit, or to define, the scope of the invention.

Examples

[0062] Historical data was gathered from a fluidized bed polyethylene reactor that was run for several months. The data was gathered during operation for both the Venturi flowmeter and the sensor measurements required to determine the superficial gas velocity via a modified Ergun equation (EQ. 1). During this time, several metallocene grades were generated, and the data was then able to be analyzed across all metallocene grades. Across all metallocene grades, as illustrated in FIG. 3, the superficial gas velocity from the Ergun equation was in agreement with the superficial gas velocity obtained from the Venturi flowmeter, although the modified Ergun equation generated more noise around the accurate reading. Additionally, the data was analyzed for an individual metallocene grade to determine the trend of an individual metallocene grade. When analyzed on an individual level, as illustrated in FIG. 4, the noise was reduced. It may be seen that the modified Ergun equation generated a superficial gas velocity that was in agreement with the superficial gas velocity obtained from the Venturi flowmeter while having less noise around the accurate reading. The modified Ergun equation may be seen to generate the superficial gas velocity with a similar accuracy to the Venturi flowmeter.

[0063] Therefore, the present invention is well adapted to attain the ends and advantages mentioned as well as those that are inherent therein. The particular examples and configurations disclosed above are illustrative only, as the present invention may be modified and practiced in different but equivalent manners apparent to those skilled in the art having the benefit of the teachings herein. Furthermore, no limitations are intended to the details of construction or design herein shown, other than as described in the claims below. It is therefore evident that the particular illustrative examples disclosed above may be altered, combined, or modified and all such variations are considered within the scope and spirit of the present invention. The invention illustratively disclosed herein suitably may be practiced in the absence of any element that is not specifically disclosed herein and/or any optional element disclosed herein. While compositions and methods are described in terms of “comprising,” “containing,” or “including” various components or steps, the compositions and methods can also “consist essentially of” or “consist of” the various components and steps. All numbers and ranges disclosed above may vary by some amount. Whenever a numerical range with a lower limit and an upper limit is disclosed, any number and any included range falling within the range is specifically disclosed. In particular, every range of values (of the form, “from about a to about b,” or, equivalently, “from approximately a to b,” or, equivalently, “from approximately a-b”) disclosed herein is to be understood to set forth every number and range encompassed within the broader range of values. Also, the terms in the claims have their plain, ordinary meaning unless otherwise explicitly and clearly defined by the patentee. Moreover, the indefinite articles “a” or “an,” as used in the claims, are defined herein to mean one or more than one of the element that it introduces.