# Patent

## Title
Methods for determining a volumetric ratio of a material to the total materials in a mixing vessel

## Abstract

Methods are provided for determining an estimated volumetric ratio of a material to total materials in a mixing vessel. In various embodiments, the methods may comprise: summing a commanded input flowrate of the material and a volumetric disturbance flowrate of the material being fed to the mixing vessel; estimating the output flowrate of the material exiting the mixing vessel; negatively feeding back the estimated output flowrate of the material to obtain an estimated volumetric rate of change of the material in the mixing vessel; and integrating the estimated volumetric rate of change of the material to compute the estimated volumetric ratio of the material to the total materials in the mixing vessel.

## Background
BACKGROUND AND SUMMARY OF THE INVENTIONS
The present invention generally relates to process control, and more particularly to methods of estimating the volumetric ratio of a material to the total materials in a mixing vessel.
The following applications filed concurrently herewith are not necessarily related to the present application, but are incorporated by reference herein in their entirety:
“Systems for Determining a Volumetric Ratio of a Material to the Total Materials in a Mixing Vessel” (U.S. application Ser. No. 11/323,323, filed simultaneously with the effective filing date of the present application);
“Systems for Volumetrically Controlling a Mixing Apparatus” (U.S. application Ser. No. 11/323,322, filed simultaneously with the effective filing date of the present application); and
“Methods of Volumetrically Controlling a Mixing Apparatus,” (U.S. application Ser. No. 11/323,324, filed simultaneously with the effective filing date of the present application).
Control systems are currently being employed to control processes for mixing together multiple components in a mixing vessel. An example of such a process is mixing together dry cement and water to form a cement slurry for use in well cementing. Well cementing is a process in which wells that penetrate subterranean formations are formed in the earth, allowing natural resources such as oil or gas to be recovered from those formations. In well cementing, a wellbore is drilled while a drilling fluid is circulated through the wellbore. The circulation of the drilling fluid is then terminated, and a string of pipe, e.g., casing, is run in the wellbore. Next, primary cementing is typically performed whereby a slurry of cement in water is placed in the annulus, which is located between the exterior of the pipe and the walls of the wellbore. Within the annulus, the cement slurry is allowed to set, i.e., harden into a solid mass, to thereby attach the string of pipe to the walls of the wellbore and seal the annulus. Subsequent secondary cementing operations, i.e., any cementing operation after the primary cementing operation, may also be performed. One example of a secondary cementing operation is squeeze cementing whereby a cement slurry is forced under pressure to areas of lost integrity in the annulus to seal off those areas.
Conventional control systems for such a cement mixing process often attempt to control the output flowrate and output density of the mixture exiting the mixing process by controlling the positions of input valves into the system. In the example in which the input valves are an input water valve and an input cement valve, an output slurry density measurement and a total output flowrate measurement are commonly used to control the process. A Proportional-Integral-Derivative (PID) controller may be used to calculate the commanded input water flowrate based on the total commanded input flowrate and the commanded slurry density. It may also be used to calculate the output water flowrate based on the total measured output flowrate and the measured slurry density. Further, a PID controller may be used to calculate the commanded input cement flowrate based on the commanded total input flowrate and the commanded slurry density. Moreover, it may be used to calculate the output cement flowrate based on the total measured output flowrate and the measured slurry density. However, this type of control system has a major drawback in that the response of the water and cement control loops are time lagged. Thus, a change in the water flowrate usually is not observed and corrected for by the cement control loop for some time and vice versa. As a result, oscillations in the density and flowrate may be experienced, especially during transitional phases such as an input disturbance or a commanded change. Another drawback of this control system is that often no densitometer is available to measure the output slurry density, or the output slurry density is ill-conditioned to be used as a control variable (i.e., the value of the density of one component being mixed is very close to the value of the density of the other component being mixed in a two-component system).
The physical system, i.e., the mixing process, being controlled typically exhibits some nonlinear behavior. Using a PI or a PID control system to overcome the physical system nonlinearity results in eigenvalue migration. That is, the eigenvalues, i.e., the parameters that define the control system, are dependent on the operating conditions such as the flowrate and thus experience relatively large shifts in value as the operating conditions change. Unfortunately, the system is a coupled system in that different portions of the system depend upon each other. Thus, fine tuning the control system is typically impossible to accomplish due to the differing time- or frequency-domain responses of the different portions of the system.
In addition to these limitations, the mixing process often experiences disturbances that can lead to inaccuracies in the measurements of the process. Such disturbances include oscillations in the height of the fluid in the mixing vessel, particularly when the mixing vessel is in motion such as in a ship-based mixing process. Another disturbance commonly encountered is that one material, e.g., the dry cement, may become plugged in the pipe being fed to the mixing vessel such that a significant amount of air is required to force the material into the mixing vessel. As such, the fluid in the mixing vessel may contain unaccounted for air.
A need therefore exists for a control system capable of controlling the output flowrate and composition of a mixing process without needing to control or measure the output density of the process. Further, it is desirable to reduce the lag-time of the control system, allowing the process to be monitored and controlled in real time with more accuracy and precision. It is also desirable that the control system be capable of more robustly accounting for disturbances, nonlinearities, and noise that may occur in the mixing process.
Methods of Determining a Volumetric Ratio of a Material to the Total Material in a Mixing Vessel
Some teachings and advantages found in the present application are summarized briefly below. However, note that the present application may disclose multiple embodiments, and not all of the statements in this section necessarily relate to all of those embodiments. Moreover, none of these statements limit the claims in any way.
The estimated volume of a material to total materials in a mixing vessel may be determined using a volumetric ratio observer comprising a feedback loop. The volumetric ratio observer advantageously provides for filtered, zero-lag estimations of the actual volumetric ratios within the mixing vessel in a manner that accounts for unwanted disturbances in the system. By way of example, the materials being combined in the mixing vessel may be dry cement and water, and the slurry formed therein may be pumped down a wellbore during a well cementing process. Knowing the relative volumes inside the mixing vessel at any time and thus the relative volumes of the cement and water being pumped downhole may be very useful.
The volumetric ratio observer may also be employed to estimate the volumetric ratios of the components in two or more mixing vessels in series that are separated by weirs or any other channeling devices that allow fluid to pass from one vessel to the next. The volumetric ratio observer desirably may be used in a control system of such a mixing process where the density of the slurry mixture is unavailable. It may also be employed to control the mixing process even if the densities of the materials being mixed are near the same value such that a densitometer cannot clearly differentiate between them. The volumetric ratio observer allows the mixing process to be controlled volumetrically, providing for tighter control over the relative volumes of the materials in the mixing vessels. As a result, the process may be optimized such that the overall cost of the process is minimized.
Systems for controlling such a mixing process may include multiple volumetric ratio observers (also referred to as volumetric estimators) for estimating the volumes of the respective components in the mixed product and a feedback block for combining at least one physical measurement of the mixed product with the estimators to provide a closed loop system. That is, the respective estimated volumes may be improved by feeding a correction based on one or more physical measurements of the system to the volumetric ratio observers. Examples of physical measurements include measurements of the height, density, total weight, and viscosity of the mixed product.
As described previously, the physical system, i.e., the mixing process, may be affected by nonlinearities. In particular, the feedback of the physical system may be nonlinear, making it difficult to control. To overcome this problem, a virtual system may be simulated in real-time, wherein the virtual system i.e., the system “seen” by the control system, represents the physical dynamics without the nonlinear physical system feedback. The physical system may be controlled with reference to this virtual system by simply using a proportional controller, resulting in more stable eigenvalue behavior than if the physical system were controlled by a PI or PID controller.
The outputs from the control loop employed to control the virtual system may be used to control at least one low-level control loop that controls low-level inputs to the physical system. The low-level inputs may be, for example, the positions of valves through which the respective components flow into the mixing process. In a first process, high-level commanded inputs (e.g., the height of the fluid or the volumetric ratio of one component to total components) may be converted into intermediate commanded targets (e.g., the total input flowrate and the input flowrate of one component). In at least one additional process, those intermediate commanded targets may be converted into the low-level control inputs or settings described above. In various embodiments, a disturbance value is fed back into the first process to decouple nonlinearities. In more embodiments, the desired or measured value of the total flowrate of the mixed product may be fed forward to decouple the effects of the output flow.



BRIEF DESCRIPTION OF THE DRAWING
FIG. 1 depicts a mixing apparatus comprising two mixing vessels separated by a weir.
FIG. 2 is a state block diagram of an embodiment of a physical system and a flow modulator being used to volumetrically mix components in the mixing vessels shown in FIG. 1.
FIG. 3A is a state block diagram of an embodiment of a portion of a volumetric ratio observer for use with a single mixing vessel.
FIG. 3B is a state block diagram of another embodiment of a portion of a volumetric ratio observer for use with two mixing vessels.
FIG. 4 is a state block diagram of an embodiment of a control system for controlling the mixing apparatus depicted in FIG. 1.
FIG. 5 is a state block diagram of another embodiment of a control system for controlling the mixing apparatus depicted in FIG. 1.
FIG. 6 is a state block diagram of yet another embodiment of a control system for controlling the mixing apparatus depicted in FIG. 1.
FIG. 7 is a state block diagram of yet another embodiment of a portion of a volumetric ratio observer for use with mixing three components utilizing a two-vessel mixing apparatus.
FIG. 8 is a state block diagram of still another embodiment of a volumetric ratio observer for use with mixing three components utilizing a two-vessel mixing apparatus.
FIG. 9 shows how a process for mixing multiple components in a mixing apparatus comprising a single vessel or tank can be controlled using a volumetric ratio mixing control scheme.
FIG. 10 shows results obtained from systems according to FIGS. 2-9.
FIG. 11 shows yet another embodiment, with a different implementation of the disclosed volumetric control ideas.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
Physical System Model
The physical system considered here is a mixing apparatus comprising two mixing vessels 10 and 12, e.g., tanks, separated by a weir 14 as shown in FIG. 1. It is understood that weir 14 may be replaced by other forms of channeling fluid from mixing vessel 10 to mixing vessel 12. The mixing process may be carried out through the action of rotating paddles 16 and 18 in respective mixing vessels 10 and 12. Two different materials may be separately added to mixing vessel 10 through pipes 20 and 24. Valves 22 and 26 may be disposed in respective pipes 20 and 24 for controlling the flow of the materials into mixing vessel 10. Within mixing vessel 10, the two materials are mixed together using rotating paddle 16. The mixture formed in mixing vessel 10 may then flow over weir 14 into mixing vessel 12 where the mixing process continues with second rotating paddle 18. The mixture in mixing vessel 12 is finally pumped out of the mixing apparatus through an output pipe 28 in which a pump 30 is disposed. The mixing system depicted in FIG. 1 may reside on the ground, on an oil platform, or on a ship.
In the embodiment depicted in FIG. 1, water and dry cement are the materials being subjected to the mixing process. It is understood that in other embodiments, liquids other than water and dry additives other than cement could be subjected to the mixing process. The volumetric flowrates of the water and the dry cement supplied to mixing vessel 10 are represented in FIG. 1 as {dot over (V)}w and {dot over (V)}c, respectively. The mixing apparatus may be capable of mixing the dry cement and the water to a desired density at a desired volumetric flowrate as required for use in oil well cementing applications. Additional parameters shown in FIG. 1 include the volumetric slurry flowrate {dot over (V)}12 over the weir from mixing vessel 10 to mixing vessel 12, the slurry height h1 in mixing vessel 10, the output slurry rate {dot over (V)}s from mixing vessel 12, and the slurry height h2 in mixing vessel 12. In various embodiments of the mixing apparatus, the approximate values for these parameters of the actual physical system are as follows:

    
    
        {dot over (V)}s ranges from about 1 bbl/min (barrels per minute) to about 15 bbl/min;
        ({dot over (V)}s/{dot over (V)}s) ranges from about 0.3 to about 0.90;
        h1 is approximately 4 ft. as defined by the weir height;
        h2 is approximately controlled to 3.5 ft.;
        h1A1 is approximately 220 gallons; and
        h2A2 is approximately controlled to 175 gallons.

In alternative embodiments, the mixing apparatus may be designed to run at a {dot over (V)}s of up to 100 bbl/min.

    
    


The physical system can be modeled mathematically using the law of mass conservation in a control volume, which is represented for mixing vessel 10 by the following equation:

ρw{dot over (V)}w+ρc{dot over (V)}c−ρ12{dot over (V)}12+{dot over (m)}D={dot over (ρ)}12h1A1+ρ12{dot over (h)}1A1   (1)

where ρw is the density of water, ρc is the dry cement density, ρ12 is the density of the slurry flowing over the weir, and A1 is the cross-sectional area of mixing vessel 10. The parameter {dot over (m)}D represents the sum of all disturbances accounting for unknown mass rate inputs into the system such as the input mass rate of air. The derivation of Equation 1 assumes instantaneous mixing such that any change in the relative proportions of {dot over (V)}w and {dot over (V)}c is instantaneously realized in the resulting value of the slurry density in mixing vessel 10. With this simplification ρ12 now represents the density of all the slurry in mixing vessel 10 at any given moment. The conservation of mass equation for mixing vessel 12 is given as follows:

ρ12{dot over (V)}12−ρs{dot over (V)}s={dot over (ρ)}sh2A2+ρs{dot over (h)}2A2   (2)

where ρs is the density of the output slurry and A2 is the cross-sectional area of mixing vessel 12. Equation 2 also assumes instantaneous mixing such that ρs represents the density of all the slurry in mixing vessel 12 at any given moment.

The physical system can also be modeled mathematically by volume conservation assuming that both the water and the cement added to the system are incompressible. This model is represented for mixing vessel 10 and mixing vessel 12 by the following respective equations:

{dot over (V)}w+{dot over (V)}c−{dot over (V)}12+{dot over (V)}D={dot over (h)}1A1   (3)

{dot over (V)}12−{dot over (V)}s={dot over (h)}2A2   (4)

The parameter {dot over (V)}D in Equation 3 represents the “volumetric disturbance flowrate,” which is herein defined as the sum of the flowrates of inputs, e.g., air, into the mixing process other than the primary materials being mixed. The term {dot over (V)}12, which represents the volumetric flowrate over the weir, is a non-linear function of the weir shape, fluid rheology and the height of fluid in mixing vessel 10. If the weir shape and the fluid rheology are assumed to be constant, {dot over (V)}12 is predominantly a function of h1 as indicated by the following equation:

{dot over (V)}12=F(h1)   (5)

It is understood that the equations herein could also be applied to other forms of channeling the slurry from one mixing vessel to the next. Thus, Equation 5 could also define the volumetric flow rate through other forms of channeling devices besides a weir.

FIG. 2 depicts the Laplace frequency domain state block diagram of the physical system 34 modeled by Equations 1 through 5, which will be described in more detail later. The inputs of water and dry cement to the system are shown to come from respective supply lines 36 and 46 that feed a physical water valve 38 and cement valve 48. These valves 38 and 48 are the control point for both slurry density and slurry flowrate through the system. The valves 38 and 48 also represent the boundary between the physical system and the control process.
Flow Modulator
A procedure known as the Flow Modulator 32 is also shown in FIG. 2 that incorporates the following Equations 6 through 13 by modulating from commanded volumetric flowrates to actual volumetric flow and mass rates through water and cement valves 38 and 48. The positions of valves 38 and 48 directly affect the rate of water and dry cement being input into the system. The resulting input volumetric rate and input mass rate may be represented by the following equations:

{dot over (V)}in={dot over (V)}w+{dot over (V)}c   (6)

{dot over (m)}in=ρin{dot over (V)}in=ρw{dot over (V)}w+ρc{dot over (V)}c   (7)

where ρin is the combined instantaneous density of both input water and dry cement. As can be seen from Equations 6 and 7, the input rates {dot over (V)}w and {dot over (V)}c are directly coupled with respect to volumetric rate and density of the slurry through the system. Designing separate control algorithms for the water valve and cement valve could produce a system in which {dot over (V)}w and {dot over (V)}c are competing to control density and flowrate simultaneously, resulting in undesirable behavior. As such, {dot over (V)}in and {dot over (m)}in may be chosen as the decoupled control variable. Through these control variables the density and volumetric flowrate can be controlled independently from each other. The desired input volumetric rate {dot over (V)}in* and the desired input mass rate {dot over (m)}in* can be modeled by the following equations:

{dot over (m)}in*={circumflex over (ρ)}w{dot over (V)}w*+{circumflex over (ρ)}c{dot over (V)}c*   (8)

{dot over (V)}in*={dot over (V)}w*+{dot over (V)}c*   (9)

where {dot over (V)}w* and {dot over (V)}c* represent the desired commanded rates of water and dry cement to each valve, respectively. The parameters {circumflex over (ρ)}w and {circumflex over (ρ)}c represent the predetermined estimated values of water density and dry cement density. Rearranging Equations 7 and 8 the commanded rates to the valves can be represented as follows:





  
    
      
        
          
            V
            .
          
          w
          *
        
        =
        
          
            (
            
              1
              
                
                  
                    ρ
                    ^
                  
                  c
                
                -
                
                  
                    ρ
                    ^
                  
                  w
                
              
            
            )
          
          ⁡
          
            [
            
              
                
                  
                    ρ
                    ^
                  
                  c
                
                ⁢
                
                  
                    V
                    .
                  
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                  *
                
              
              -
              
                
                  m
                  .
                
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
                *
              
            
            ]
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            V
            .
          
          c
          *
        
        =
        
          
            (
            
              1
              
                
                  
                    ρ
                    ^
                  
                  c
                
                -
                
                  
                    ρ
                    ^
                  
                  w
                
              
            
            )
          
          ⁡
          
            [
            
              
                
                  m
                  .
                
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
                *
              
              -
              
                
                  
                    ρ
                    ^
                  
                  w
                
                ⁢
                
                  
                    V
                    .
                  
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                  *
                
              
            
            ]
          
        
      
    
    
      
        (
        11
        )
      
    
  




In order to verify that {dot over (V)}in and {dot over (m)}in are actually decoupled, the output rate of each valve is assumed to closely approximate the commanded input rate to each valve as follows:

{dot over (V)}w≈{dot over (V)}w*   (12)

{dot over (V)}c≈{dot over (V)}c*   (13)

Combining Equations 6 through 13 results in the following set of equations:





  
    
      
        
          
            V
            .
          
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
        =
        
          
            V
            .
          
          
            i
            ⁢
            
              
            
            ⁢
            n
          
          *
        
      
    
    
      
        (
        14
        )
      
    
  
  
    
      
        
          
            m
            .
          
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
        =
        
          
            
              (
              
                
                  
                    ρ
                    c
                  
                  -
                  
                    ρ
                    w
                  
                
                
                  
                    
                      p
                      ^
                    
                    c
                  
                  -
                  
                    
                      ρ
                      ^
                    
                    w
                  
                
              
              )
            
            ⁢
            
              
                m
                .
              
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
              *
            
          
          +
          
            
              (
              
                
                  
                    
                      ρ
                      w
                    
                    ⁢
                    
                      
                        ρ
                        ^
                      
                      c
                    
                  
                  -
                  
                    
                      ρ
                      c
                    
                    ⁢
                    
                      
                        ρ
                        ^
                      
                      w
                    
                  
                
                
                  
                    
                      p
                      ^
                    
                    c
                  
                  -
                  
                    
                      ρ
                      ^
                    
                    w
                  
                
              
              )
            
            ⁢
            
              
                V
                .
              
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
              *
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




Equations 14 and 15 verify that the volumetric input rate is completely independent of the mass input rate. Additionally, if {circumflex over (ρ)}c≅ρc and {circumflex over (ρ)}w≅ρw, then Equation 15 reduces to

{dot over (m)}in≅{dot over (m)}in*   (16)

and the mass input rate becomes independent of the volumetric input flowrate. If the density estimations are incorrect or the valve delivery is not approximated exactly as assumed in Equations 12 and 13, these “errors” may be absorbed into the modeled disturbance terms {dot over (V)}D and {dot over (m)}D.

The density of the slurry mixture may be unavailable due to a lack of a density measuring device or to the density values of the dry cement and water being very similar (i.e., ρw≅ρc) such that density is a poorly conditioned variable for good control. A mixing system in which the input water rate {dot over (V)}w into the first mixing vessel, the fluid height h2 in the second mixing vessel, and the output slurry rate {dot over (V)}s from the second mixing vessel are available for measurement may be controlled using a so-called volumetric ratio mixing control approach. That is, the mixing process may be controlled volumetrically, and the chosen control variables may be the overall total flowrate of the slurry through the system and the percentage or ratio of the slurry which is water.
In an embodiment in which the density is no longer the variable by which the water and the cement are proportioned, {circumflex over (ρ)}w may be set equal to 1 ({circumflex over (ρ)}w=1) and {circumflex over (ρ)}c may be set equal to 0 ({circumflex over (ρ)}c=0). Turning back to FIG. 2, the inputs to physical system 34 now become the overall commanded input rate {dot over (V)}in* and the commanded input water rate {dot over (V)}w* into the first mixing vessel. A Flow Modulator 32 is shown that incorporates Equations 6 through 13 by modulating from commanded volumetric flowrates to actual volumetric flow and mass rates through water and cement valves 38 and 48. The Flow Modulator 32 may send {dot over (V)}w* directly to water valve 38 via signal 36. Further, it may send the overall commanded input rate {dot over (V)}n* via signal 40 and the commanded input water rate {dot over (V)}w* via signal 42 to a summation block 44 where {dot over (V)}w* is subtracted from {dot over (V)}in* to obtain the commanded input cement rate {dot over (V)}c*, which may then be sent to cement valve 48 via signal 46. The positions of valves 38 and 48 may be set according to those commanded input rates.
The resulting water flowrate {dot over (V)}w exiting water valve 38 and the resulting cement flowrate {dot over (V)}c exiting cement valve 48 may be measured. The total input mass flowrate {dot over (m)}in to the mixing process is the result of the summation (summation block 60) of the water mass flow rate {dot over (V)}w (signal 50) multiplied by ρw (gain element 52) and the cement mass flow rate {dot over (V)}c (signal 56) multiplied by ρc (gain element 58) as described in Equation 7.
Next, {dot over (m)}in may be sent to another summation block 67 to which the mass disturbance flowrate {dot over (m)}D, the total mass flowrate out of the first mixing vessel, and the total mass flowrate within the first mixing vessel also may be sent. At summation block 67, the total mass flowrate out of the first mixing vessel and the total mass flowrate within the first mixing vessel may be subtracted from the sum of {dot over (m)}in and {dot over (m)}D to obtain the total mass rate of change in the first mixing vessel. The total mass rate of change may then be sent via signal 72 to an Integral controller comprising gain element 74 for multiplying the total mass rate of change by 1/h1A1 to obtain the total density rate of change in the first mixing vessel. The Integral controller also comprises an integral element 76 for multiplying the total density rate of change by 1/s, which is the laplace transform representation of integration, to determine the density of the mixture flowing over the weir, ρ12. The Integral controller may then feed ρ12 back to summation block 67 via signals 78 and 82. On its way to summation block 67, signal 78 may pass through gain element 80 where it is multiplied by {dot over (h)}1A1 to obtain the total mass flowrate in the first mixing vessel. Also, signal 80 may pass through gain element 84 where it is multiplied by {dot over (V)}12 to obtain the total mass flowrate out of the first mixing vessel, i.e., over the weir. In this manner, the Integral controller may dynamically recompute ρ12. After being sent to integral element 76, signal 72 further may be sent to gain element 86 where it is multiplied by the total output volumetric flowrate from the first mixing vessel {dot over (V)}12 to obtain the total mass flowrate {dot over (m)}12 before being sent to yet another summation block 90.
At summation block 90, the total mass flowrate in the second mixing vessel, indicated by signal 96, and the total mass flowrate out of the second mixing vessel (a measured value), indicated by signal 88, may be subtracted from the {dot over (m)}12 to obtain the total mass rate of change in the second mixing vessel. The total mass rate of change may then be sent via signal 90 to an Integral controller comprising gain element 92 for multiplying the total mass rate of change by 1/h2A2 to obtain the total density rate of change in the first mixing vessel. The Integral controller also comprises an integral element 94 for determining the density of the slurry flowing out of the second mixing vessel, ρs. The Integral controller may then feed ρs back to summation block 90 via signal 96. On its way to summation block 67, signal 96 may pass through gain element 98 where it is multiplied by {dot over (h)}2A2 to obtain the total mass flowrate in the second mixing vessel. In this manner, the Integral controller may dynamically recompute ρs.
As further shown in FIG. 2, the cement flowrate {dot over (V)}c exiting cement valve 48 and the water flowrate {dot over (V)}w exiting water valve 38 may be sent via signals 54 and 62, respectively, to a summation block 64 to obtain the total volumetric input flowrate {dot over (V)}in. Then {dot over (V)}in and a total volumetric disturbance flowrate {dot over (V)}D may be sent to summation block 100 via signals 65 and 102, respectively. The volumetric flowrate within the first mixing vessel may also be fed back to summation block 100 where it is subtracted from the sum of {dot over (V)}in and {dot over (V)}D to obtain the total volumetric rate of change in the first mixing vessel. The volumetric mass rate of change may then be sent via signal 104 to an Integral controller comprising gain element 106 for multiplying the total volumetric rate of change by 1/A1 to obtain the total height rate of change in the first mixing vessel. The Integral controller also comprises an integral element 108 for determining the height of the mixture in the second mixing vessel, h1. The Integral controller may then feed h1 back to summation block 100 via signal 110. On its way to summation block 100, signal 110 may pass through gain element 112 where it is multiplied by F(h1) to obtain the total volumetric flowrate out of the first mixing vessel {dot over (V)}12. In this manner, the Integral controller may dynamically recompute h1.
Additionally, signal 104 may be sent to a gain element 114 for multiplying h1 by F(h1) to determine {dot over (V)}12 before it is sent to yet another summation block 115. The total flowrate of the slurry exiting the second mixing vessel, {dot over (V)}s, is also sent via signal 116 to summation block 115 where it is subtracted from {dot over (V)}12 to obtain the total volumetric rate of change in the second mixing vessel. The output of summation block 115 is further sent to gain elements 120 and 122 for multiplying the total volumetric rate of change by 1/A2 and 1/s, respectively, to thereby determine the height of the slurry in the second mixing vessel, h2.
Volumetric Ratio Observer
The volumetric ratio of one material relative to the total materials in one of the mixing vessels may be determined using a Volumetric Ratio Observer. This observer is based on the same physical dynamics described above and may be derived in a way that does not include density parameters. That is, a single mixing vessel with N number of components being mixed together therein can be modeled using the law of mass conservation by the following equation:

ρ1({dot over (V)}in)1+ρ2({dot over (V)}in)2+ . . . +ρN({dot over (V)}in)N−ρout{dot over (V)}out+{dot over (m)}D={dot over (ρ)}outVT+ρout{dot over (V)}T   (17)

where ρN is the density of the Nth component being mixed, ({dot over (V)}in)N is the volumetric flowrate at which the Nth component is being added to the mixing vessel, ρout is the density of the mixture flowing out of the mixing vessel, {dot over (V)}out is the output flowrate of the mixture from the mixing vessel, and VT is the volume of the mixture currently in the mixing vessel. The parameter {dot over (m)}D represents the sum of all disturbances accounting for unknown mass rate inputs into the system and is given as follows:

{dot over (m)}D=ρ1({dot over (V)}D)1+ρ2({dot over (V)}D)2+ . . . +ρN({dot over (V)}D)N   (18)

where ({dot over (V)}D)N represents the unknown volumetric flowrate disturbance of the Nth component. The total volumetric flowrate disturbance {dot over (V)}D is given as the sum of all the component disturbances as follows:

{dot over (V)}D=({dot over (V)}D)1+({dot over (V)}D)2+ . . . +({dot over (V)}D)N   (19)

Using instantaneous mixing as described before, the density ρout may be represented by the following equation:





  
    
      
        
          ρ
          out
        
        =
        
          (
          
            
              
                
                  
                    ρ
                    1
                  
                  ⁡
                  
                    (
                    
                      V
                      T
                    
                    )
                  
                
                1
              
              +
              
                
                  
                    ρ
                    2
                  
                  ⁡
                  
                    (
                    
                      V
                      T
                    
                    )
                  
                
                2
              
              +
              ...
              +
              
                
                  
                    ρ
                    N
                  
                  ⁡
                  
                    (
                    
                      V
                      T
                    
                    )
                  
                
                N
              
            
            
              
                
                  (
                  
                    V
                    T
                  
                  )
                
                1
              
              +
              
                
                  (
                  
                    V
                    T
                  
                  )
                
                2
              
              +
              ...
              +
              
                
                  (
                  
                    V
                    T
                  
                  )
                
                N
              
            
          
          )
        
      
    
    
      
        (
        20
        )
      
    
  




where (VT)N represents the volume of the Nth component currently in the mixing vessel. The total volume of the mixture in the mixing vessel VT may be represented by the following equation:

VT=(VT)1+(VT)2+ . . . +(VT)N   (21)

Eliminating ρout from Equations 17 through 21 and grouping terms with common density coefficients, the resulting volumetric equations describing the separate component flow through the mixing vessel are given as follows:




  
    
      
        
          
            
              
                
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        
                          i
                          ⁢
                          
                            
                          
                          ⁢
                          n
                        
                      
                      )
                    
                    1
                  
                  -
                  
                    
                      (
                      
                        
                          
                            (
                            
                              V
                              T
                            
                            )
                          
                          1
                        
                        
                          V
                          T
                        
                      
                      )
                    
                    ⁢
                    
                      
                        V
                        .
                      
                      out
                    
                  
                  +
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        D
                      
                      )
                    
                    1
                  
                
                =
                
                  
                    (
                    
                      
                        V
                        .
                      
                      T
                    
                    )
                  
                  1
                
              
            
          
          
            
              
                
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        
                          i
                          ⁢
                          
                            
                          
                          ⁢
                          n
                        
                      
                      )
                    
                    2
                  
                  -
                  
                    
                      (
                      
                        
                          
                            (
                            
                              V
                              T
                            
                            )
                          
                          2
                        
                        
                          V
                          T
                        
                      
                      )
                    
                    ⁢
                    
                      
                        V
                        .
                      
                      out
                    
                  
                  +
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        D
                      
                      )
                    
                    2
                  
                
                =
                
                  
                    (
                    
                      
                        V
                        .
                      
                      T
                    
                    )
                  
                  2
                
              
            
          
          
            
              ⋮
            
          
          
            
              
                
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        
                          i
                          ⁢
                          
                            
                          
                          ⁢
                          n
                        
                      
                      )
                    
                    N
                  
                  -
                  
                    
                      (
                      
                        
                          
                            (
                            
                              V
                              T
                            
                            )
                          
                          N
                        
                        
                          V
                          T
                        
                      
                      )
                    
                    ⁢
                    
                      
                        V
                        .
                      
                      out
                    
                  
                  +
                  
                    
                      (
                      
                        
                          V
                          .
                        
                        D
                      
                      )
                    
                    N
                  
                
                =
                
                  
                    (
                    
                      
                        V
                        .
                      
                      T
                    
                    )
                  
                  N
                
              
            
          
        
        }
      
    
    
      
        (
        22
        )
      
    
  




Here the volumetric ratio of the Nth component with respect to the overall volume of the mixture is given as





  
    
      
        
          
            (
            
              R
              out
            
            )
          
          N
        
        =
        
          (
          
            
              
                (
                
                  V
                  T
                
                )
              
              N
            
            
              V
              T
            
          
          )
        
      
    
    
      
        (
        23
        )
      
    
  




where the notation (Rout)N incorporates the instantaneous mixing assumption, indicating not only the volumetric ratio of the Nth component to the total materials in the mixing vessel but also the volumetric flowrate ratio of the Nth component to the total output flowrate {dot over (V)}out. Combining Equations 21 and 23 gives the relationship between all the component volumetric ratios as follows:

(Rout)1+(Rout)2+ . . . +(Rout)N=1   (24)

The state block diagrams of the primary components of Volumetric Ratio Observers (VRO's) for a single mixing vessel and for two mixing vessels are shown in FIG. 3A and FIG. 3B, respectively. Using the same notation from earlier, the symbol (^) indicates that a parameter has been estimated. The commanded or setpoint inputs for the VRO's represent the commanded rates to the actual physical system and are signified by the symbol (*). The VRO may be implemented using various arrangements of closed loops, as will be detailed later in particular embodiments of the Volumetric Mixing Control approach. The VRO may serve to decouple the effects of disturbances in the system.
As shown in FIG. 3A, the VRO for the Nth component being fed to a single mixing vessel may include a summation block 55 for subtracting an estimated output flowrate of the Nth component ({dot over ({circumflex over (V)}out)N from the sum of a volumetric disturbance flowrate of the Nth component, ({dot over ({circumflex over (V)}D)N, and a commanded input flowrate of the Nth component ({dot over (V)}in*)N. The ({dot over ({circumflex over (V)}out)N may be fed to summation block 55 via signal 65, the ({dot over ({circumflex over (V)}D)N may be fed to summation block 55 via signal 51, and the ({dot over (V)}in*)N may be fed to summation block 55 via signal 53. The output of summation block 55, as indicated by signal 57, may represent an estimated volumetric rate of change of the Nth component in the mixing vessel. The estimated volumetric rate of change of the Nth component may be fed to an Integral controller comprising an integral element 59 for computing the estimated volume of the Nth component in the mixing vessel. The Integral controller may also include a gain element 61 for multiplying the estimated volume of the Nth component in the mixing vessel by 1/(the estimated volume of the total materials in the mixing vessel) to obtain the estimated output ratio of the Nth component to the total materials in the mixing vessel, ({circumflex over (R)}out)N. It may further include gain element 63 for multiplying ({circumflex over (R)}out)N by the total estimated output flowrate from the mixing vessel, {dot over (V)}out, to estimate the output flowrate of the Nth component ({dot over ({circumflex over (V)}out)N, which may be negatively fed back to summation block 55. Thus, the Integral controller dynamically recomputes ({dot over ({circumflex over (V)}out)N.
As illustrated in FIG. 3B, the Volumetric Ratio Observer may be expanded to cover two mixing vessels. In this embodiment, the term ({dot over ({circumflex over (V)}12)N is the estimated flowrate of the Nth component out of the first mixing vessel and into the second mixing vessel. The total estimated flowrate between the two mixing vessels, i.e., over the weir, may be represented by

{dot over ({circumflex over (V)}12=({dot over ({circumflex over (V)}12)1+({dot over ({circumflex over (V)}12)2+ . . . +({dot over ({circumflex over (V)}12)N   (25)

The volumes of the mixture in the first mixing vessel and the second mixing vessel are given by {circumflex over (V)}1 and {circumflex over (V)}2, respectively. The portion of the state block diagram depicted in FIG. 3B that dynamically recomputes ({dot over ({circumflex over (V)}12)N using a first Integral controller is the same as the state block diagram shown in FIG. 3A with the exception that gain element 61 multiplies by 1/{circumflex over (V)}1 and gain element 63 multiplies by {dot over ({circumflex over (V)}12. The ({dot over ({circumflex over (V)}12)N computed by the first Integral controller may be sent to a summation block 67 via signal 79 where an estimated value of the output flowrate ({dot over ({circumflex over (V)}out)N from the second mixing vessel is subtracted from ({dot over ({circumflex over (V)}12)N to obtain the volumetric rate of change in the second mixing vessel. This volumetric rate of change is then sent to a second Integral controller via signal 69. The second Integral controller comprises an integral element 71 for determining the total volume of the Nth component in the second mixing vessel, ({circumflex over (V)}2)N. It also comprises gain element 73 for multiplying ({circumflex over (V)}2)N by 1/{circumflex over (V)}2 to determine ({circumflex over (R)}out)N and gain element 75 for multiplying ({circumflex over (R)}out )N by {dot over (V)}out, thereby determining the estimated output flowrate of the Nth component from the second mixing vessel, ({dot over ({circumflex over (V)}out)N. The ({dot over ({circumflex over (V)}out)N may then be negatively fed back to summation block 67 via signal 77 such that it may be dynamically recomputed. It is understood that the VRO is not limited to one or two mixing vessels but may be used for any number of mixing vessels by the addition of an Integral controller for each additional mixing vessel. Further, control schemes like those shown in FIGS. 3A and 3B may be implemented for any component being mixed in the one or more mixing vessels. There is no limit to the number of components that may be mixed together using the control system described herein.

Cement Mixing Control Scheme

FIG. 4 illustrates one embodiment of the volumetric ratio mixing control scheme mentioned earlier. The process being controlled comprises mixing cement and water together in a mixing apparatus containing two mixing vessels separated by a weir as shown in FIG. 1. FIG. 4 depicts a control system 130 that includes two Height Observers 132 and 134, a State Feedback Controller 136, a Flow Regulator 138, and a Volumetric Ratio Observer 140. The Flow Modulator 32 and the state block diagram of the physical system 34 that are depicted in FIG. 2 are also shown in FIG. 4. A detailed description of these parts of the control scheme may be found in the previous discussion of FIG. 2.
The first Height Observer 132 depicted in FIG. 4 takes as input the measured height h2 of fluid in the second mixing vessel, the measured output flowrate of the slurry {dot over (V)}s exiting the second mixing vessel, and the overall commanded volumetric input flowrate {dot over (V)}in*. This Height Observer 132 then estimates the fluid height in the second mixing vessel, which is used as feedback in State Feedback Controller 136. It further estimates the overall volumetric disturbance flowrate {dot over ({circumflex over (V)}D, which is used for disturbance input decoupling in Flow Regulator 138. The second Height Observer 134 depicted in FIG. 4, also known as the Weir Flow Observer, takes as input h2 and {dot over (V)}s. With only these two inputs, Height Observer 134 estimates the flowrate of the fluid {dot over ({circumflex over (V)}12 flowing over the weir from the first mixing vessel to the second mixing vessel.
Describing Height Observer 132 in more detail, h2 may be fed from physical system 34 to a summation block 146 via signal 142. The estimated height of the fluid ĥ2 in the second mixing vessel may also be sent via signal 144 to summation block 146 where it is subtracted from h2 to determine an estimation of a height error for the second mixing vessel. This estimation of height error may then be fed via signal 148 to a Proportional-Integral controller 152 comprising an integral element 154, an integral gain element 156 for multiplying it by a constant Nio1, and a proportional gain element 150 for multiplying it by a constant No1. The PI gains may be set to remove the noise and oscillations of the second mixing vessel from the height estimation. The output of integral gain element 156 and of proportional gain element 150 may then be summed at a summation block 158 to estimate the total volumetric disturbance flowrate {dot over ({circumflex over (V)}D. The {dot over ({circumflex over (V)}D may be sent via signal 160 to another summation block 166. In addition, both the {dot over (V)}in* and the {dot over (V)}s as measured by a sensor may be fed to summation block 166 via signals 162 and 164, respectively. At summation block 166, the {dot over (V)}s may be subtracted from the sum of {dot over (V)}in* and {dot over ({circumflex over (V)}D to obtain the volumetric rate of change in the second mixing vessel. The output of summation block 166 may be sent to an Integral controller comprising a gain element 170 where it is multiplied by 1/(the estimated cross-sectional area of the second mixing vessel) to convert the volumetric rate of change to the rate of height change in the second mixing vessel. This rate of height change may be sent to an integral element 172 to compute ĥ2. The Height Observer 132 may continue to dynamically recompute ĥ2 in this manner.
As shown in FIG. 4, the Weir Flow Observer 134 may be very similar to the Height Observer 132. That is, it may also include a summation block 178 to which h2 is fed via signal 142 and ĥ2 is negatively fed via signal 176. The output of summation block 178, i.e., an estimation of a height error for the second mixing vessel, may then be fed via signal 148 to a Proportional-Integral controller 184 comprising an integral element 186, an integral gain element 188 for multiplying it by a constant Nio2, and a proportional gain element 182 for multiplying it by a constant No2. The output of integral gain element 188 and of proportional gain element 182 may then be summed at a summation block 190 to estimate the total output flowrate {dot over ({circumflex over (V)}12 from the first mixing vessel. The {dot over ({circumflex over (V)}12 may then be sent via signal 192 to another summation block 194 to which the {dot over (V)}s may also be sent via signal 195. At summation block 194, the {dot over (V)}s may be subtracted from the {dot over ({circumflex over (V)}12 to obtain the volumetric rate of change in the second mixing vessel. The output of summation block 194 may be sent to an Integral controller comprising a gain element 198 where it is multiplied by 1/(the estimated cross-sectional area of the second mixing vessel) to convert the volumetric rate of change to the rate of height change in the second mixing vessel. This rate of height change may be sent to an integral element 200 to compute ĥ2. The Height Observer 134 may continue to dynamically recompute ĥ2 in this manner. Additional information related to height observers may be found in U.S. patent application Ser. No. 11/029,072, entitled “Methods and Systems for Estimating a Nominal Height or Quantity of a Fluid in a Mixing Tank While Reducing Noise,” filed on Jan. 4, 2005, which is incorporated by reference herein in its entirety.
In order to maintain enough fluid in the physical system to supply a desired output flowrate of the slurry {dot over (V)}s* from the second mixing vessel, a State Feedback Controller 136 may be implemented where ĥ2 is the state feedback. In particular, the ĥ2 determined by Height Observer 132 may be sent via signal 204 to a summation block 206 where it is subtracted from the commanded height of the fluid h2* in the second mixing vessel, indicated by signal 202, to estimate the height error for the second mixing vessel. The output of summation block 206 may be sent to a proportional gain element 210 via signal 216 where it is multiplied by the constant Np before being summed with {dot over (V)}s* at summation block 214. In this manner, State Feedback Controller 136 computes a commanded output flowrate {dot over (V)}12* of the total materials from the first mixing vessel. This desired output flowrate is then sent through Flow Regulator 138 and Flow Modulator 32 to adjust the water and cement valves as needed.
This implementation of Height Observer 132 and Weir Flow Observer 134 with full state feedback allows control system 130 to be fully enhanced. These height observers not only provide filtered, zero-lag estimations of the actual signals but also provide for disturbance estimation.
The estimated output flowrate {dot over ({circumflex over (V)}12 from the first mixing vessel determined by Weir Flow Observer 134 may be fed back to an upper portion of Flow Regulator 138 to “cancel” or decouple the negative state feedback that naturally occurs in the physical system. The estimated total volumetric disturbance flowrate {dot over ({circumflex over (V)}D determined by Height Observer 132 results from input air and errors between commanded volumetric rates and actual volumetric rates through the valves. This volumetric disturbance flowrate estimation may be negatively fed back to Flow Regulator 138 to decouple the effect of the disturbance in the system, thereby making the control system invariant to unmeasured volumetric flowrate disturbances.
Describing Flow Regulator 138 in more detail, {dot over ({circumflex over (V)}12 may be fed back to summation block 220 where it is subtracted from {dot over (V)}12*, which is fed to summation block 220 via signal 216. The output of summation block 220 may then be sent to a proportional gain element 224 where it may be multiplied by a constant KV before being sent to another summation block 230. The {dot over ({circumflex over (V)}D may be fed back to summation block 230 via signal 226, and {dot over ({circumflex over (V)}12 may also be fed to summation block 230 via signal 228 such that {dot over ({circumflex over (V)}D is subtracted from the sum of the output of gain element 224 and {dot over ({circumflex over (V)}12. The output of summation block 230 is the total commanded input flowrate {dot over (V)}in* to the mixing process, which may be fed to Flow Modulator 32 via signal 232. As described previously, Flow Modulator 32 may modulate from the commanded flowrate {dot over (V)}in* to the actual input flowrate {dot over (V)}in.
The Volumetric Ratio Observer 140 shown in FIG. 4 may be implemented to estimate the ratio of water to total materials in the first mixing vessel in accordance with the following equation:




  
    
      
        
          
            (
            
              
                R
                ^
              
              12
            
            )
          
          w
        
        =
        
          (
          
            
              
                (
                
                  
                    V
                    ^
                  
                  1
                
                )
              
              w
            
            
              
                V
                ^
              
              1
            
          
          )
        
      
    
    
      
        (
        26
        )
      
    
  




The inputs to Volumetric Ratio Observer 140 may include the commanded input water rate {dot over (V)}w and the measured input water rate {dot over (V)}w as well as the closed loop estimate of the volumetric disturbance {dot over ({circumflex over (V)}D from Height Observer 132. The measured and commanded input water rates may be used to estimate the input disturbance flowrate ({dot over ({circumflex over (V)}D)w in the water delivery. This disturbance may be used for disturbance input decoupling within Flow Regulator 138 and to determine the input disturbance flowrate ({dot over ({circumflex over (V)}D)c in the cement delivery within Volumetric Ratio Observer 140.

In the embodiment of Volumetric Ratio Observer 140 shown in FIG. 4, a comparator 240 is employed to determine the estimated volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w by comparing {dot over (V)}w to {dot over (V)}w*, which are fed to comparator 240 via signals 236 and 238, respectively. The (ė{dot over (,)}ċ{dot over (i)}{dot over (r)}{dot over ( )}{dot over (V)}D)w may then be fed to a summation block 266 to which {dot over (V)}w* is also fed via signal 262. Further, an estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel may be negatively fed to summation block 266. At summation block 266, the ({dot over ({circumflex over (V)}12)w may be subtracted from the summation of ({dot over ({circumflex over (V)}D)w and {dot over (V)}w* to determine an estimated volumetric rate of change of the water in the first mixing vessel. The output of summation block 266 may be fed via signal 268 to an Integral controller comprising an integral element 270 and a gain element 272 for multiplying it by 1/{circumflex over (V)}1, thereby determining the estimated volumetric ratio ({circumflex over (R)}12)w of the water to the total materials in the first mixing vessel. The Integral Controller further comprises another gain element 274 for multiplying ({circumflex over (R)}12)w by the total estimated output flowrate {dot over ({circumflex over (V)}12 from the first mixing vessel to estimate the output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel. This estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w may then be fed back to summation block 266 via signal 264. The Integral controller continues to dynamically recompute the estimated rate ({dot over ({circumflex over (V)}12)w in this manner.
The Volumetric Ratio Observer 140 may also determine the volumetric disturbance flowrate of the cement ({dot over ({circumflex over (V)}D)c through the use of another summation block 244 for subtracting the volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w from the total input volumetric disturbance flowrate {dot over ({circumflex over (V)}D determined by Height Observer 132. The ({dot over ({circumflex over (V)}D)w may be fed from the output of comparator 240 to summation block 244 via signal 242, and the {dot over ({circumflex over (V)}D may be fed to summation block 244 via signal 234. The volumetric disturbance flowrate of the cement ({dot over ({circumflex over (V)}D)c may then be sent to yet another summation block 252 via signal 246. Further, a commanded input cement flowrate {dot over (V)}c* and an estimated output flowrate of the cement ({dot over ({circumflex over (V)}12)c from the first mixing vessel may be fed to summation block 252 via signals 248 and 250, respectively. At summation block 252, the ({dot over ({circumflex over (V)}12)c may be subtracted from the summation of ({dot over ({circumflex over (V)}D)c and {dot over (V)}c* to determine an estimated volumetric rate of change of the cement in the first mixing vessel. The output of summation block 252 may be fed via signal 254 to an Integral controller comprising an integral element 256, a gain element 258 for multiplying it by 1/{circumflex over (V)}1, and another gain element 260 for multiplying it by the total estimated output flowrate {dot over ({circumflex over (V)}12 from the first mixing vessel. As a result, the estimated volumetric rate of change of the cement in the first mixing vessel may be converted to the estimated output flowrate of the cement ({dot over ({circumflex over (V)}12)c from the first mixing vessel. This estimated output flowrate of the cement ({dot over ({circumflex over (V)}12)c may then be fed back to summation block 252 via signal 250. The Integral controller continues to dynamically recompute the estimated rate ({dot over ({circumflex over (V)}12)c in this manner.
The estimated water ratio ({circumflex over (R)}12)w in the first mixing vessel may be fed back and compared to the desired water ratio (R12*)w in a Proportional controller within a lower portion of Flow Regulator 138. That is, the ({circumflex over (R)}12)w may be fed via signal 278 from Volumetric Ratio Observer 140 to a comparator 279 of Flow Regulator 138. Further, the (R12*)w may be fed via signal 276 to comparator 279. The output of comparator 279 may then be fed via signal 280 to a proportional gain element 282 for multiplying it by Km before being sent to a summation block 288 of Flow Regulator 138. The estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w exiting the first mixing vessel also may be fed back to Flow Regulator 138 for decoupling purposes. That is, the ({dot over ({circumflex over (V)}12)w may be fed via signal 286 to summation block 288. Further, the estimated volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w may be fed to summation block 288. At summation block 288, the ({dot over ({circumflex over (V)}D)w may be subtracted from the summation of the output of gain element 282 and ({dot over ({circumflex over (V)}12)w, thereby computing the commanded input flowrate of the water {dot over (V)}w*. The {dot over (V)}w* may be fed to Flow Modulator 32 via signal 290. As described previously, Flow Modulator 32 may modulate from the total commanded input volumetric flowrate {dot over (V)}in* and the commanded input volumetric flowrate of the water {dot over (V)}w* to the actual total input mass flowrate {dot over (m)}in. (See FIG. 2).
The foregoing implementation of Volumetric Ratio Observer 140 with state feedback allows control system 130 to be fully enhanced. The VRO provides for filtered, zero-lag estimations of actual signals.
FIG. 5 illustrates another embodiment of the volumetric ratio mixing control scheme in which the process being controlled comprises mixing cement and water together in a mixing apparatus containing two mixing vessels separated by a weir as shown in FIG. 1. FIG. 5 depicts a control system 291 that is the same as control system 130 of FIG. 4 except for some changes in the Volumetric Ratio Observer, the State Feedback Controller, and the Flow Regulator. In particular, this embodiment extends the VRO in FIG. 4 from a one vessel implementation to a two vessel implementation for estimating the ratio of water to total materials in the second mixing vessel rather than the first mixing vessel. This embodiment also provides for water ratio control within the State Feedback Controller.
The differences of FIG. 5 are described in more detail below, beginning with Volumetric Ratio Observer 141. In particular, the estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel may be further passed to another summation block 292. At summation block 292, an estimated output flowrate of the water ({dot over ({circumflex over (V)}s)w from the second mixing vessel may be subtracted from ({dot over ({circumflex over (V)}12)w to determine the volumetric rate of change in the second mixing vessel. The output of summation block 292 may then be sent via signal 294 to an Integral controller comprising an integral element 296 and a gain element 298 for multiplying it by 1/{circumflex over (V)}2 to determine the estimated volumetric ratio of the water to total materials {circumflex over (R)}w in the second mixing vessel. The Integral controller may further include a gain element 300 for multiplying {circumflex over (R)}w by the total output flowrate {dot over (V)}s from the second mixing vessel, which may be measured, to determine the estimated total output flowrate ({dot over ({circumflex over (V)}s)w of the water. The ({dot over ({circumflex over (V)}s)w may be fed back to summation block 292 via signal 302, allowing it to be dynamically recomputed.
Another difference between Volumetric Ratio Observer 141 and Volumetric Ratio Observer 140 is that the estimated output flowrate of the cement ({dot over ({circumflex over (V)}12)c from the second mixing vessel may be further passed to another summation block 303. At summation block 303, an estimated output flowrate of the cement ({dot over ({circumflex over (V)}s)c from the second mixing vessel may be subtracted from ({dot over ({circumflex over (V)}12)c to determine the volumetric rate of change in the second mixing vessel. The output of summation block 303 may then be sent via signal 304 to an Integral controller comprising an integral element 306 and a gain element 308 for multiplying it by 1/{circumflex over (V)}2 to determine the estimated volumetric ratio of the cement to total materials {circumflex over (R)}c in the second mixing vessel. The Integral controller may further include a gain element 310 for multiplying {circumflex over (R)}c by the total output flowrate {dot over (V)}s from the second mixing vessel, which may be measured, to determine the estimated total output flowrate ({dot over ({circumflex over (V)}s)c of the cement. Further, the ({dot over ({circumflex over (V)}s)c may be fed back to summation block 303 via signal 312, allowing it to be dynamically recomputed.
In this embodiment, State Feedback Controller 137 may be different in that it may engage in proportional control of the volumetric ratio of the water to the total materials in the second mixing vessel, comparing the desired water ratio Rw* to the estimated water ratio {circumflex over (R)}w determined by Volumetric Ratio Observer 141. The {circumflex over (R)}w may be calculated using the following equation:




  
    
      
        
          
            R
            ^
          
          w
        
        =
        
          (
          
            
              
                (
                
                  
                    V
                    ^
                  
                  2
                
                )
              
              w
            
            
              
                V
                ^
              
              2
            
          
          )
        
      
    
    
      
        (
        27
        )
      
    
  




Describing State Feedback Controller 137 in more detail, the {circumflex over (R)}w and the Rw* may be fed to comparator 318 via signals 314 and 316 respectively. The output of comparator 318 may be fed to a proportional gain element 322 for multiplying it by a constant Kp and a gain element 324 for multiplying it by the desired total output flowrate {dot over (V)}12* from the first mixing vessel before being positively sent to a summation block 330. The Rw* may also pass through a gain element 328 for multiplying it by the total desired output flowrate of the slurry {dot over (V)}s* exiting the second mixing vessel to determine the desired output flowrate of the water exiting the second vessel. This desired output flowrate of the water may be positively fed forward to summation block 330 via signal 326 to decouple the effect of the water exiting the second mixing vessel.

The output of State Feedback Controller 137 and the estimated flowrate of water ({dot over ({circumflex over (V)}12)w out of the first mixing vessel, as determined by Volumetric Ratio Observer 141, may be fed via respective signals 276 and 278 to a summation block 279 of Flow Regulator 139 where they may be compared. The Flow Regulator 193 is implemented in the same way as Flow Regulator 138 in FIG. 4 with the exception that the proportional control compares the estimated flowrate of water ({dot over ({circumflex over (V)}12)w out of the first mixing vessel with the commanded flowrate of water ({dot over (V)}12*)w from State Feedback Controller 137. In particular, the output of summation block 279 may be sent via signal 280 to a proportional gain element 282 for multiplying it by a constant Km before sending it to another summation block 288, where the estimated volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w may be subtracted from it and from the estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel. The ({dot over ({circumflex over (V)}D)w determined by Volumetric Ratio Observer 141 may be sent to summation block 288 via signal 284. Further, the estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel may be sent to summation block 288. The output of summation block 288 may be the commanded input water flowrate {dot over (V)}w*, which may be fed to Flow Modulator 32.
FIG. 6 depicts yet another embodiment of the volumetric ratio mixing control scheme in which the process being controlled comprises mixing cement and water together in a mixing apparatus containing two mixing vessels separated by a weir as shown in FIG. 1. FIG. 6 depicts a control system 331 that is similar to control system 130 of FIG. 4. Notably, control system 331 does not contain a Weir Flow Observer. Further, this embodiment extends the VRO in FIG. 4 from a one vessel implementation to a two vessel implementation for estimating the ratio of water to total materials in the second mixing vessel rather than the first mixing vessel. This two vessel VRO may also be used to estimate the total volumetric disturbance flowrate by applying an internal PI controller to the fluid height in the second mixing vessel. Moreover, within the VRO the commanded total volumetric flowrate {dot over (V)}12* out of the first mixing vessel may be used as an estimate of the actual flowrate out of the first mixing vessel to determine the state feedback decoupling term for the Flow Regulator.
In this embodiment, a PI control loop may act directly on the water valve within the Flow Modulator using the actual measured input water flowrate as feedback (not shown). Tuned for a faster response time than the rest of the system, the water valve thus may be driven to produce the desired input water flowrate, resulting in zero steady state error. Therefore, an assumption may be made that all resulting volumetric disturbances are a result of errors in the cement valve between the commanded input cement flowrate and the actual delivered input cement flowrate (({dot over ({circumflex over (V)}D)w=0; {dot over ({circumflex over (V)}D=({dot over ({circumflex over (V)}D)c). As mentioned earlier, the VRO may determine this disturbance by closing a loop on the estimated height of fluid in the second mixing vessel. The estimated height of fluid ĥ2 in the second mixing vessel may be found by assuming the estimated cross-sectional area Â2 of the second mixing vessel is known for a given volume of fluid in the vessel.
Since the estimated volumetric disturbance term {dot over ({circumflex over (V)}D is assumed to only contain errors due to the cement valve, it is only fed back into the upper portion of the Flow Regulator. When fed through the Flow Modulator this only makes adjustments to the cement command. In summary, valve errors in both valves are decoupled by the combined effects of the PI control on the water valve and the disturbance input decoupling on the cement valve.
The differences between control system 331 in FIG. 6 and control system 130 in FIG. 4 are described in more detail below. The volumetric disturbance flowrate of the cement being fed to summation block 248 via signal 246 may be determined by first feeding the height of the fluid h2 in the second mixing vessel to summation block 330 via signal 320. At summation block 330, the sum of the estimated height of the water (ĥ2)w and the estimated height of the cement (ĥ2)c in the second mixing vessel may be subtracted from h2, thereby estimating the height error for the second mixing vessel. This height error may be sent to a PI controller 332 via signal 338. The PI controller may comprise an integral element 334, an integral gain element 336 for multiplying the height error by Nio1, and a proportional gain element 340 for multiplying it by No1 before sending it to summation block 342. The output of summation block 342 is the estimated volumetric disturbance flowrate of the cement, which is equivalent to the estimated total volumetric disturbance flowrate {dot over ({circumflex over (V)}D as represented by signal 344. Also, no estimated volumetric disturbance flowrate of the water is fed to summation block 266 nor to summation block 288 of Flow Regulator 143 since this estimated rate is equivalent to zero.
Additionally, in Volumetric Ratio Observer 145, respective gain elements 274 and 260 may be replaced by respective gain elements 275 and 261, which may multiply the respective estimated ratios of the water and the cement in the first mixing vessel by the commanded total output flowrate {dot over (V)}12* from the first mixing vessel. Moreover, the estimated output flowrate of the water ({dot over ({circumflex over (V)}12)w from the first mixing vessel may be further passed to another summation block 292. At summation block 292, an estimated output flowrate of the water ({dot over ({circumflex over (V)}s)w from the second mixing vessel may be subtracted from ({dot over ({circumflex over (V)}12)w to determine the volumetric rate of change in the second mixing vessel. The output of summation block 292 may then be sent via signal 294 to an Integral controller comprising an integral element 296 and a gain element 299 for multiplying it by 1/Â2 to determine the estimated height of the water (ĥ2)w in the second mixing vessel. The Integral controller may further include a gain element 301 for multiplying (ĥ2)w by 1/h2 and a gain element 300 for multiplying (ĥ2)w by the measured total output flowrate {dot over (V)}s from the second mixing vessel to determine the estimated total output flowrate ({dot over ({circumflex over (V)}s)w of the water. The ({dot over ({circumflex over (V)}s)w may be fed back to summation block 292 via signal 302, allowing it to be dynamically recomputed.
Another difference between Volumetric Ratio Observer 145 and Volumetric Ratio Observer 140 is that the estimated output flowrate of the cement ({dot over ({circumflex over (V)}12)c from the second mixing vessel may be further passed to another summation block 303. At summation block 303, an estimated output flowrate of the cement ({dot over ({circumflex over (V)}s)c from the second mixing vessel may be subtracted from ({dot over ({circumflex over (V)}12)c to determine the volumetric rate of change in the second mixing vessel. The output of summation block 303 may then be sent via signal 304 to an Integral controller comprising an integral element 306 and a gain element 309 for multiplying it by 1/Â2 to determine the estimated height of the cement (ĥ2)w in the second mixing vessel. The Integral controller may further include a gain element 311 for multiplying (ĥ2)c by 1/h2 and a gain element 310 for multiplying (ĥ2)c by the measured total output flowrate {dot over (V)}s from the second mixing vessel to determine the estimated total output flowrate ({dot over ({circumflex over (V)}s)c of the water. Further, the ({dot over ({circumflex over (V)}s)c may be fed back to summation block 303 via signal 312, allowing it to be dynamically recomputed. The (ĥ2)w and (ĥ2)c may be fed to and added together at summation block 326 before being fed to comparator 330 via signal 328.
In addition, the total estimated volumetric disturbance flowrate {dot over ({circumflex over (V)}D determined by Volumetric Ratio Observer 145 may be negatively fed to a summation block 230 of Flow Regulator 143, which does not contain a proportional controller for the volumetric flowrate exiting the first mixing vessel as in FIG. 4. Instead, the {dot over ({circumflex over (V)}D may be subtracted from the commanded total output flowrate {dot over (V)}12* from the first mixing vessel, which may be fed to summation block 230 via signal 216. As in FIG. 4, the output of summation block 230 is the total commanded input flowrate {dot over (V)}in*. FIG. 6 also depicts {dot over (V)}in* being fed via signal 162 from Flow Regulator 143 to summation block 166 of Height Observer 132.
Modifications and Variations
For a mixing system in which the measured parameters include the input water flowrate {dot over (V)}w into the first mixing vessel, the slurry density ρ12 in the first mixing vessel, the fluid height h2 in the second mixing vessel, and the output slurry flowrate {dot over (V)}s from the second mixing vessel, any of the embodiments discussed previously may be employed to control the system. However, one of the inherent problems with the mixing system depicted in FIG. 1 is the introduction of air into the mixture. Air entrained in the mixture may cause the overall slurry volume in the mixing vessels to be larger than expected, resulting in an increased h2 value. Additionally, air entrained in the mixture may cause the measured density of the mixture to be lower than expected. For most applications it is ideal to be able to mix the water and the cement to a density and a volume that does not reflect the entrainment of air. Fortunately, for a system that includes four sensors for the four measured parameters mentioned above, the Volumetric Ratio Observer may be implemented to predict the ratio of not only the water and cement in the mixture but also the amount of air entrained therein. As such, the system can be controlled to mix exactly the desired proportions of water and cement.
FIG. 7 illustrates an embodiment of the primary components of a two vessel Volumetric Ratio Observer 350 for modeling a system in which three components, i.e., water, cement, and air, are mixed through the system. The Volumetric Ratio Observer 350 includes control schemes 352, 354, and 356 for the water, the cement, and the air, respectively. Those control schemes are very similar to the control scheme shown in FIG. 3B except that the gain element for multiplying 1/s by the estimated total output flowrate {dot over ({circumflex over (V)}12 from the first mixing vessel is replaced by a gain element for multiplying 1/s by a commanded total output flowrate from the first mixing vessel {dot over (V)}12*. In this embodiment, the commanded input water flowrate {dot over (V)}w * and the measured input water flowrate {dot over (V)}w are also known, allowing the disturbance in water flowrate to be calculated directly. That is, the {dot over (V)}w and the {dot over (V)}w* may be fed via respective signals 364 and 366 to a comparator 368 for comparing the two and thus determining the estimated volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w.
Disturbances due to the cement delivery and due to the entrained air may be provided from external observer controllers that may be implemented via hardware or software modules. The total mass disturbance flowrate {dot over ({circumflex over (m)}D may be estimated by a Density Observer and the total volumetric disturbance flowrate {dot over ({circumflex over (V)}D may be estimated by a Height Observer as described previously. Suitable density observers are described in U.S. patent application Ser. No. 11/121,278, filed on May 3, 2005. Using the estimated parameter values of water density and cement density, these disturbances may be converted into the estimated volumetric flowrate disturbance ({dot over ({circumflex over (V)}D)c of the cement and the estimated disturbance due to the volumetric flowrate of entrained air ({dot over ({circumflex over (V)}D)a. An assumption is made that the density of air is relatively insignificant (ρa≅0) compared to the density of water and cement.
More specifically, the ({dot over ({circumflex over (V)}D)w computed by comparator 368 may be multiplied by the estimated density of water by passing it to a gain element 372 before sending it to another comparator 376 via signal 370. Gain element 372 determines the estimated mass flowrate of the water. The {dot over ({circumflex over (m)}D is also sent to comparator 376 via signal 374 where it may be compared to the estimated mass flowrate of the water to determine the estimated mass flowrate of the cement. This estimated mass flowrate of the cement may be sent via signal 360 to gain element 378 where it is multiplied by 1/(the estimated density of cement) to determine ({dot over ({circumflex over (V)}D)c. Both ({dot over ({circumflex over (V)}D)w and ({dot over ({circumflex over (V)}D)c may be fed to a summation block 384 via respective signals 380 and 382 where they may be subtracted from {dot over ({circumflex over (V)}D, which is sent to element 384 via signal 364, to determine ({dot over ({circumflex over (V)}D)a. The disturbance flowrates ({dot over ({circumflex over (V)}D)w, ({dot over ({circumflex over (V)}D)c, and ({dot over ({circumflex over (V)}D)a for each component may then be sent to controllers via respective signals 358, 360, and 362 to implement respective control schemes 352, 354, and 356.
Using the foregoing implementation, the components may be separated and the densities of the water and cement mixture excluding entrained air for the first mixing vessel and the second mixing vessel may be calculated from estimated parameters within the VRO in accordance with the following equations:




  
    
      
        
          
            ρ
            ^
          
          12
        
        =
        
          (
          
            
              
                
                  
                    
                      ρ
                      ^
                    
                    w
                  
                  ⁡
                  
                    (
                    
                      
                        V
                        ^
                      
                      1
                    
                    )
                  
                
                w
              
              +
              
                
                  
                    
                      ρ
                      ^
                    
                    c
                  
                  ⁡
                  
                    (
                    
                      
                        V
                        ^
                      
                      1
                    
                    )
                  
                
                c
              
            
            
              
                
                  (
                  
                    
                      V
                      ^
                    
                    1
                  
                  )
                
                w
              
              +
              
                
                  (
                  
                    
                      V
                      ^
                    
                    1
                  
                  )
                
                c
              
            
          
          )
        
      
    
    
      
        (
        28
        )
      
    
  
  
    
      
        
          
            ρ
            ^
          
          s
        
        =
        
          (
          
            
              
                
                  
                    
                      ρ
                      ^
                    
                    w
                  
                  ⁡
                  
                    (
                    
                      
                        V
                        ^
                      
                      2
                    
                    )
                  
                
                w
              
              +
              
                
                  
                    
                      ρ
                      ^
                    
                    c
                  
                  ⁡
                  
                    (
                    
                      
                        V
                        ^
                      
                      2
                    
                    )
                  
                
                c
              
            
            
              
                
                  (
                  
                    
                      V
                      ^
                    
                    2
                  
                  )
                
                w
              
              +
              
                
                  (
                  
                    
                      V
                      ^
                    
                    2
                  
                  )
                
                c
              
            
          
          )
        
      
    
    
      
        (
        29
        )
      
    
  




FIG. 8 illustrates another embodiment of the primary components of Volumetric Ratio Observer 386 for modeling a system in which water, cement, and air are mixed in a two-vessel mixing apparatus. The Volumetric Ratio Observer 386 includes control schemes 388, 390, and 392 for the water, the cement, and the air, respectively. Those control schemes are very similar to the control scheme shown in FIG. 3B except that the {dot over ({circumflex over (V)}12 gain element may be replaced by a {dot over (V)}12* gain element. Again, with the commanded input water flowrate {dot over (V)}w* and the measured input water flowrate {dot over (V)}w being known, the disturbance in water flowrate may be calculated directly. That is, the {dot over (V)}w and the {dot over (V)}w* may be fed via respective signals 400 and 402 to a comparator 404 for comparing the two and thus determining the estimated volumetric disturbance flowrate of the water ({dot over ({circumflex over (V)}D)w.
In this embodiment, disturbances due to cement delivery and due to entrained air are provided from internal PI feedback loops on the slurry density in the first mixing vessel and the fluid height in the second mixing vessel as shown. The mass disturbance flowrate {dot over ({circumflex over (m)}D may be calculated through a PI controller that compares the measured slurry density in the first mixing vessel to the estimated density calculated from the combined water, cement, and air mixture in the first mixing vessel. More specifically, the volumetric ratio of each component to the total materials in the first mixing vessel may be calculated by the PI controller of each control scheme. Those volumetric ratios may then be sent via respective signals 410, 412, and 414 to respective gain elements 416, 418, and 420 for multiplying the volumetric ratios by the estimated densities of air {circumflex over (ρ)}a, of cement {circumflex over (ρ)}c, and of water {circumflex over (ρ)}w, respectively to determine the estimated fraction of the total density in the first mixing vessel for each component. Those estimated fractions may then be summed at summation block 422 to determine the estimated total density of the slurry {circumflex over (ρ)}12 in the first mixing vessel. The measured slurry density ρ12 and the estimated slurry density {circumflex over (ρ)}12 may be sent to a comparator 428 for calculating the difference between the two and then to a PI controller 430 for determining {dot over ({circumflex over (m)}D.
The ({dot over ({circumflex over (V)}D)w computed by comparator 404 may be multiplied by the estimated density of water by passing it to a gain element 408 for determining the estimated input mass flowrate of the water before sending it to another comparator 432 via signal 406 where it is subtracted from the {dot over ({circumflex over (m)}D. The output of comparator 432 thus may be the estimated mass flowrate of the cement. The estimated mass flowrate of the cement may be passed through gain element 434 where it may be multiplied by 1/(the estimated density of cement) to determine ({dot over ({circumflex over (V)}D)c.
The total volumetric disturbance flowrate {dot over ({circumflex over (V)}D may be calculated through a PI controller that compares the measured fluid height h2 to the estimated height ĥ2 in the second mixing vessel calculated from the combined water, cement, and air volumes in the second mixing vessel, assuming that its cross-sectional area is known. More specifically, ĥ2 may be calculated by sending the volumes of water, cement, and air in the second mixing vessel, as determined via control schemes 388, 390, and 392, to a summation block 442 via signals 436, 438, and 440, respectively. At summation block 442, those volumes may be summed together to determine the total volume of fluid in the second mixing vessel. The total volume may then be sent to a gain element 444 for multiplying it by 1/(the estimated cross-sectional area of the second mixing vessel) to determine ĥ2 before being sent to comparator 450. The comparator 450 may determine the difference between h2 and ĥ2, and that difference may be sent to a PI controller 452 via signal 451. The outputs of the integral portion and the proportional portion of PI controller 452 may then be summed at summation block 454 to determine the {dot over ({circumflex over (V)}D. Then the {dot over ({circumflex over (V)}D may be sent to another summation block 458 via signal 456. Both ({dot over ({circumflex over (V)}D)w and ({dot over ({circumflex over (V)}D)c may be fed to a summation block 458 via respective signals 460 and 462 where they may be subtracted from {dot over ({circumflex over (V)}D to determine the volumetric disturbance flowrate in the air ({dot over ({circumflex over (V)}D)a. The disturbance flowrates ({dot over ({circumflex over (V)}D)w, ({dot over ({circumflex over (V)}D)c, and ({dot over ({circumflex over (V)}D)a for each component may then be sent to controllers via respective signals 394, 396, and 398 to implement respective control schemes 388, 390, and 392. Additionally, Equations 28 and 29 may be implemented to estimate the mixture densities in the first and second mixing vessels due to water and cement but excluding entrained air.
As shown in FIG. 9, a process for mixing multiple components in a mixing apparatus comprising a single vessel or tank may also be controlled using a volumetric ratio mixing control scheme. In one embodiment, the components being combined in the mixing apparatus may be cement and water. However, it is understood that other materials may also be combined in the single vessel. FIG. 9 depicts a control system 500 that includes a Flow Regulator 502, a Height Observers 506, and a Volumetric Ratio Observer 530. The Flow Regulator 502 includes a Flow Modulator 32, shown in detail in FIG. 2, a State Feedback Controller 510, and a model of a physical system 508 similar to the physical system 34 shown in FIG. 2. The physical system 508 is different from physical system 34 of FIG. 2 in that it only models a single mixing vessel with the height and density of the mixture in the single mixing vessel given as outputs. That is, the volumetric rate of change in the mixing vessel of physical system 508 is converted to the rate of change of height in the mixing vessel, which when integrated results in the height h of the slurry in the mixing vessel. Further, the mass rate of change in the mixing vessel of physical system 508 is converted to the rate of change of the density in the mixing vessel, which when integrated results in the density ρ of the slurry in the mixing vessel.
The measured height h of the slurry in the mixing vessel as given by the model of physical system 508 may be sent to Height Observer 506, which contains the same components as the Height Observer 132 in FIG. 6. The Height Observer 506 may estimate the height ĥ of the fluid in the mixing vessel and feed that to Flow Regulator 502. The measured height h may also be fed to Volumetric Ratio Observer 530, which is similar to the Volumetric Ratio Observer 145 shown in FIG. 6 except that it only contains one feedback loop for estimating the volumetric flowrate ({dot over ({circumflex over (V)}s)w of the water exiting the mixing vessel and the ratio of water to total materials {circumflex over (R)}w in the mixing vessel and one feedback loop for estimating the volumetric flowrate ({dot over ({circumflex over (V)}s)c. The Volumetric Ratio Observer 530 may estimate the total volumetric disturbance flowrate {dot over ({circumflex over (V)}D in the same manner as does Volumetric Ratio Observer 145.
Turning to Flow Regulator 502, its upper portion includes a comparator 514 to which ĥ may be sent via signal 510 and a commanded height h* may be sent via signal 512. The comparator 514 may subtract ĥ from h*. The output of comparator 514 may then be sent to a proportional gain element 518 via signal 516 where it may be multiplied by a constant Kv before being sent to another comparator 524. A commanded volumetric flowrate {dot over (V)}s* of the slurry exiting the mixing vessel and {dot over ({circumflex over (V)}D as determined by Volumetric Ratio Observer 530 may be also be fed to comparator 524 via signals 520 and 522, respectively. The comparator 514 may subtract {dot over ({circumflex over (V)}D from the sum of the output of gain element 518 and {dot over (V)}s* to determine the total commanded input flowrate {dot over (V)}in* to the mixing vessel, which may be fed to Flow Modulator 32 via signal 526.
The lower portion of Flow Regulator 502 may include a comparator 538 to which a desired water ratio Rw* and the estimated water ratio ({circumflex over (R)}12)w in the first mixing vessel may be fed via signals 538 and 532, respectively. The comparator 538 subtracts ({circumflex over (R)}12)w from Rw*, and its output may then be fed via signal 540 to a proportional gain element 542 for multiplying the output by Km before being sent to a comparator 544. The estimated output flowrate of the water ({dot over ({circumflex over (V)}s)w from the mixing vessel also may be fed back to the Flow Regulator for decoupling purposes. That is, the ({dot over ({circumflex over (V)}s)w may be fed via signal 534 to comparator 544 where the commanded input flowrate of the water {dot over (V)}w* may be computed. The {dot over (V)}w* may be fed to Flow Modulator 32 via signal 546.
In the various embodiments shown in FIGS. 2-9, the control schemes may be implemented by hardware or by software via a computerized system. A person of ordinary skill in the art would know how to create and use such hardware or software to implement the control schemes.
EXAMPLES
The invention having been generally described, the following examples are given as particular embodiments of the invention and to demonstrate the practice and advantages thereof. It is understood that the examples are given by way of illustration and are not intended to limit the specification or the claims to follow in any manner.
The mixing apparatus shown in FIG. 1 was assembled and operated using the embodiment of the control scheme shown in FIG. 6. Various parameters of the mixing process were determined and plotted as a function of time in FIG. 10. More specifically, line 550, labeled as the slurry recirculation density, represents the change in the measured slurry density in the first vessel. Line 552, labeled as the Ve_density represents the change in the density as given by the volumetric ratio observer with active disturbance decoupling. Line 554, labeled as the tub level, represents the change in the height of the slurry in the second vessel. Line 556, labeled as the cement valve position, represents the change in the position of the valve for controlling the flowrate of the cement into the mixing apparatus. Line 558, labeled as h2_hat, represents the change in the estimated height of the slurry in the second vessel as determined by the height observer, which filters the height sensor without zero lag. Line 560, labeled as the water valve position, represents the change in the position of the valve for controlling the flowrate of the cement into the mixing apparatus.
The results shown in FIG. 10 illustrate that the system is capable of controlling the relative volumes of cement and water in the mixing tub. Having line 550 track line 552 indicates that the system is producing the desired density and therefore the desired relative volumes. Also, the tub level is maintained near a desired amount, showing that the flowrate is maintained near its desired amount. It should be noted that at time 14 hr.:17 min. the cement delivery system runs low and a new supply is initiated. This is a common occurrence and is not a problem with the control system.
FIG. 11 shows yet another embodiment, with a different implementation of the disclosed volumetric control ideas. Note the following two aspects of this embodiment:
1) Variable height control: The height setpoint is changed depending on the height observer error, or cement rate error. This is done to reduce effects of water/cement ratio problems if we have flow inconsistencies in the cement supply system. This normally occurs when we switch between cement supply bins or pods.
2) Ideal ratio control: Instead of a ratio observer outlined in the previous embodiments this system uses the idea case by only inputting the output flow rate and assuming all other values.
According to various embodiments, methods of determining an estimated volumetric ratio of a material to total materials in a mixing vessel comprise: summing a commanded input flowrate of the material and a volumetric disturbance flowrate of the material being fed to the mixing vessel; estimating the output flowrate of the material exiting the mixing vessel; negatively feeding back the estimated output flowrate of the material to obtain an estimated volumetric rate of change of the material in the mixing vessel; and integrating the estimated volumetric rate of change of the material to compute the estimated volumetric ratio of the material to the total materials in the mixing vessel.
In more embodiments, methods of determining an estimated volumetric ratio of a material to total materials in a second mixing vessel that is partially separated from a first mixing vessel comprise: summing a commanded input flowrate of the material and a volumetric disturbance flowrate of the material being fed to the first mixing vessel; estimating an output flowrate of the material exiting the first mixing vessel; negatively feeding back the estimated output flowrate of the material to obtain an estimated volumetric rate of change of the material in the first mixing vessel; integrating the estimated volumetric rate of change of the material in the first mixing vessel to dynamically recompute the estimated output flowrate of the material exiting the first mixing vessel; estimating an output flowrate of the material exiting the second mixing vessel; negatively feeding back the estimated output flowrate of the material exiting the second mixing vessel and summing it with the estimated output flowrate of the material exiting the first mixing vessel, thereby obtaining an estimation of a volumetric rate of change of the material in the second mixing vessel; and integrating the estimated volumetric rate of change of the material in the second mixing vessel to compute the estimated volumetric ratio of the material to the total materials in the second mixing vessel.
In additional embodiments, methods of determining an estimated volumetric ratio of a second material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: measuring a height of the total materials in the second mixing vessel; comparing the height of the total materials in the second mixing vessel to a summation of an estimated height of a first material in the second mixing vessel and an estimated height of the second material in the second mixing vessel to obtain an estimation of a height error for the second mixing vessel; feeding the estimation of the height error to a controller to compute an estimated total volumetric disturbance flowrate; computing a summation of (a) a commanded input flowrate of the second material to the first mixing vessel, (b) the estimated total volumetric disturbance flowrate, and (c) a negative value of an estimated output flowrate of the second material from the first mixing vessel, thereby obtaining an estimated volumetric rate of change of the second material in the first mixing vessel; and integrating the estimated volumetric rate of change of the second material to obtain the estimated volumetric ratio of the second material to total materials in the first mixing vessel.
According to other embodiments, systems for determining an estimated volumetric ratio of a material to total materials in a mixing vessel comprise: a summation block for determining an estimated volumetric rate of change of the material in the mixing vessel; an integration element for determining an estimated volume of the material in the mixing vessel based on the estimated volumetric rate of change of the material in the mixing vessel; a first gain element for converting the estimated volume of the material in the mixing vessel to the estimated volumetric ratio of the material to the total materials; and a second gain element for converting the estimated volumetric ratio of the material to the total materials to the output flowrate of the material from the mixing vessel.
In more embodiments, systems for determining an estimated volumetric ratio of a material to total materials in a second mixing vessel that is partially separated from a first mixing vessel comprise: a first summation block for determining an estimated volumetric rate of change of the material in the first mixing vessel; a first integration element for integrating the estimated volumetric rate of change of the material in the first mixing vessel to determine an estimated volume of the material in the first mixing vessel; a first gain element for converting the estimated volume of the material in the first mixing vessel to the estimated volumetric ratio of the material to the total materials in the first mixing vessel; a second gain element for converting the estimated volumetric ratio of the material to the total materials in the first mixing vessel to the output flowrate of the material from the first mixing vessel; a second summation block for determining an estimated volumetric rate of change of the material in the second mixing vessel based on the output flowrate of the material from the first mixing vessel; a second integration element for integrating the estimated volumetric rate of change of the material in the second mixing vessel to determine the estimated volume of the material in the second mixing vessel; a third gain element for converting the estimated volume of the material in the second mixing vessel to the estimated volumetric ratio of the material to the total materials in the second mixing vessel; and a fourth gain element for converting the estimated volumetric ratio of the material to the total materials in the second mixing vessel to the output flowrate of the material from the second mixing vessel.
In yet more embodiments, systems for determining an estimated volumetric ratio of a second material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: a sensor for measuring a height of the total materials in the second mixing vessel; a first summation block for determining an estimation of a height error for the second mixing vessel by comparing the height of the total materials in the second mixing vessel with a summation of an estimated height of a first material in the second mixing vessel and an estimated height of the second material in the second mixing vessel; a controller for determining an estimated volumetric disturbance flowrate of the second material based on the height error; a second summation block for determining an estimated volumetric rate of change of the second material in the first mixing vessel; an integration element for integrating the estimated volumetric rate of change of the second material in the first mixing vessel to determine the estimated volume of the second material in the first mixing vessel; a first gain element for converting the estimated volume of the second material in the first mixing vessel to the estimated volumetric ratio of the material to the total materials in the first mixing vessel; and a second gain element for converting the estimated volumetric ratio of the material to the total materials in the first mixing vessel to an output flowrate of the material from the first mixing vessel.
According to additional embodiments, methods of controlling a volumetric ratio of a material to total materials in a mixing vessel comprise: estimating the volumetric ratio of the material to the total materials in the mixing vessel and an output flowrate of the material from the mixing vessel using a volumetric ratio observer; dynamically recomputing the commanded input flowrate of the material based on outputs of the volumetric ratio observer using a flow regulator; and adjusting an input valve of the material based on the commanded input flowrate of the material using a flow modulator. In one embodiment, the mixing vessel comprises a first mixing vessel partially separated from a second mixing vessel. In this case, a height observer may be used to estimate the height of the total materials in the second mixing vessel, and the volumetric ratio observer may be used to estimate the volumetric ratio of the material to the total materials in the first mixing vessel and an output flowrate of the material from the first mixing vessel.
In additional embodiments, methods of controlling a volumetric ratio of a material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: estimating the volumetric ratio of the material to the total materials in the second mixing vessel, an output flowrate of the material from the first mixing vessel, and a volumetric disturbance flowrate of the material using a volumetric ratio observer having the following inputs: a commanded input flowrate of the material and a measured input flowrate of the material; computing a commanded output flowrate of the material from the first mixing vessel using a state feedback controller having the following inputs: a commanded volumetric ratio of the material to the total materials in the second mixing vessel and the estimated volumetric ratio of the material to the total materials in the second mixing vessel; dynamically recomputing the commanded input flowrate of the material using a flow regulator having the following inputs: the estimated input flowrate error of the material and the estimated output flowrate of the material from the first mixing vessel; and adjusting an input valve of the material based on the commanded input flowrate of the material using a flow modulator.
In yet more embodiments, methods of controlling a volumetric ratio of a material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: estimating a total volumetric disturbance flowrate, the volumetric ratio of the material to the total materials in the first mixing vessel, and an output flowrate of the material from the first mixing vessel using a volumetric ratio observer having the following inputs: a measured height of the total materials in the second mixing vessel; a commanded input flowrate of the material; and a commanded input flowrate of a second material that is also being fed to the first mixing vessel; dynamically recomputing the commanded input flowrate of the material using a flow regulator having the following inputs: a commanded volumetric ratio of the material to the total materials in the first mixing vessel; an estimated volumetric ratio of the material to the total materials in the first mixing vessel; and the estimated output flowrate of the material from the first mixing vessel; and adjusting an input valve of the material based on the commanded input flowrate of the material using a flow modulator.
According to additional embodiments, systems for controlling a volumetric ratio of a material to total materials in a mixing vessel comprise: a volumetric ratio observer for estimating the volumetric ratio of the material to the total materials in the mixing vessel and an output flowrate of the material from the mixing vessel; a flow regulator coupled to the volumetric ratio observer for dynamically recomputing a commanded input flowrate of the material based on outputs of the volumetric ratio observer; and a flow modulator coupled to the flow regulator for adjusting an input valve of the material based on the commanded input flowrate of the material. In one embodiment, the mixing vessel comprises a first mixing vessel partially separated from a second mixing vessel. In this case, a height observer may be used to estimate the height of the total materials in the second mixing vessel, and the volumetric ratio observer may be capable of estimating the volumetric ratio of the material to the total materials in the first mixing vessel and an output flowrate of the material from the first mixing vessel.
In more embodiments, systems for controlling a volumetric ratio of a material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: a volumetric ratio observer for estimating the volumetric ratio of the material to the total materials in the second mixing vessel, an output flowrate of the material from the first mixing vessel, and a volumetric disturbance flowrate of the material, the volumetric ratio observer having the following inputs: an estimated total volumetric disturbance flowrate and a commanded input flowrate of the material; a state feedback controller for computing a commanded output flowrate of the material from the first mixing vessel, the state feedback controller having the following inputs: a commanded volumetric ratio of the material to the total materials in the second mixing vessel and the estimated volumetric ratio of the material to the total materials in the second mixing vessel; a flow regulator coupled to the state feedback controller and to the volumetric ratio observer for dynamically recomputing the commanded input flowrate of the material, the flow regulator having the following inputs: the estimated volumetric disturbance flowrate of the material and the estimated output flowrate of the material from the first mixing vessel; and a flow modulator coupled to the flow regulator for adjusting an input valve of the material based on the commanded input flowrate of the material.
In still more embodiments, systems for controlling a volumetric ratio of a material to total materials in a first mixing vessel that is partially separated from a second mixing vessel comprise: a volumetric ratio observer for estimating a total volumetric disturbance flowrate, the volumetric ratio of the material to the total materials in the first mixing vessel, and an output flowrate of the material from the first mixing vessel, the volumetric ratio observer having the following inputs: a measured height of the total materials in the second mixing vessel; a commanded input flowrate of the material; and a commanded input flowrate of a second material that is also being fed to the first mixing vessel; a flow regulator coupled to the volumetric ratio observer for dynamically recomputing a commanded input flowrate of the material having the following inputs: a commanded volumetric ratio of the material to the total materials in the first mixing vessel; the estimated volumetric ratio of the material to the total materials in the first mixing vessel; and the estimated output flowrate of the material from the first mixing vessel; and a flow modulator coupled to the flow regulator for adjusting an input valve of the material based on the commanded input flowrate of the material.
In various embodiments, methods comprise: generating multiple estimates of multiple respective component volumes; generating a feedback correction using at least one physical measurement of a mixed product; and combining the feedback correction with at least one of the estimates to provide a closed loop system. In more embodiments, methods of controlling a mixing process comprise: in a first process, converting high-level commanded inputs (e.g., the height of the slurry in a mixing vessel and the volumetric ratio of one material to the total materials in a mixing vessel) into intermediate commanded targets (e.g. the desired total flowrate from a mixing vessel and the desired volumetric ratio of one material to the total materials in a mixing vessel); and in at least one additional process, converting the intermediate commanded targets into low-level control settings (e.g., the valve positions), wherein a disturbance value is fed back into the first process to decouple nonlinearities. In further embodiments, methods of controlling a mixing process comprise: in a first process, converting high-level commanded inputs into intermediate commanded targets; in at least one additional process, converting the intermediate commanded targets into low-level control settings; and using feedforward of a total flowrate to decouple effects of an output flow. In yet more embodiments, methods for controlling a mixing process which is affected by physical nonlinearities comprise: compensating for the nonlinearities to provide an equivalent virtual system having more stable eigenvalues; and controlling the mixing process with reference to the equivalent virtual system.
In various embodiments, systems comprise: multiple volumetric estimators for multiple respective components; and a feedback block for combining at least one physical measurement of a mixed product with the estimators to provide a closed loop system. In more embodiments, systems for controlling a mixing process comprise: a first process for converting high-level commanded inputs into intermediate commanded targets; and at least one additional process for converting the intermediate commanded targets into low-level control settings, wherein the system is capable of feeding a disturbance value back into the first process to decouple nonlinearities. In additional embodiments, systems for controlling a mixing process comprise: a first process for converting high-level commanded inputs into intermediate commanded targets; and at least one additional process for converting the intermediate commanded targets into low-level control settings, wherein the system is capable of using feedforward of a total flowrate to decouple effects of an output flow. In yet more embodiments, systems for controlling a mixing process which is affected by physical nonlinearities comprise: at least one low-level control loop for controlling inputs to the mixing process; a real-time simulation of an equivalent virtual system in which the physical nonlinearities of the mixing process are at least partially compensated to provide for more stable eigenvalue behavior; and an additional control loop for controlling the equivalent virtual system, wherein the system is capable of using outputs of the additional control loop to at least partially control the low-level control loop.
According to various disclosed embodiments, there is also provided: a method comprising the actions of: generating multiple estimates of multiple respective component volumes using open loop operations; generating a feedback correction using at least one physical measurement of a mixed product; and combining said correction with at least one said open loop estimate to provide a closed loop system.
Modifications and Variations
While preferred embodiments of the invention have been shown and described, modifications thereof can be made by one skilled in the art without departing from the spirit and teachings of the invention. The embodiments described herein are exemplary only and are not intended to be limiting. Many variations and modifications of the invention disclosed herein are possible and are within the scope of the invention.
Various modifications, alternatives and implementations are suggested above, but many others are possible. For example, the observers are not necessarily implemented as in any of the examples above, but can be modified in various ways. For another example, the embodiments described above do not stand in isolation, but can be combined in various ways.
The attached Appendix gives further implementation details which are believed not to be necessary for understanding the inventions.
Accordingly, the scope of protection is not limited by the description set out above but is only limited by the claims which follow, that scope including all equivalents of the subject matter of the claims. Each and every claim is incorporated into the specification as an embodiment of the present invention. Thus, the claims are a further description and are an addition to the preferred embodiments of the present invention. The discussion of a reference herein is not an admission that it is prior art to the present invention, especially any reference that may have a publication date after the priority date of this application. The disclosures of all patents, patent applications, and publications cited herein are hereby incorporated by reference, to the extent that they provide exemplary, procedural, or other details supplementary to those set forth herein.

## Claims

1. A method of operating a computer control system for a mixing vessel to determine an estimated volumetric ratio of a material to total materials in a mixing vessel, comprising:
summing a commanded input flowrate of the material and a volumetric disturbance flowrate of the material being fed to the mixing vessel;
estimating the output flowrate of the material exiting the mixing vessel based on the summed commanded input flowrate and the volumetric disturbance flowrate;
negatively feeding back the estimated output flowrate of the material to obtain an estimated volumetric rate of change of the material in the mixing vessel; and
integrating the estimated volumetric rate of change of the material to compute the estimated volumetric ratio of the material to the total materials in the mixing vessel;
wherein the mixing process computer control system controls the mixing process volumetrically utilizing the estimated volumetric ratio of the material to control at least one input valve.

2. The method of claim 1, further comprising dynamically recomputing the estimated output flowrate of the material.

3. The method of claim 1, wherein said integrating is performed by a feedback control loop.

4. The method of claim 3, wherein the feedback control loop comprises gain elements for multiplying the Laplace transform for integration, 1/s, by 1/(an estimated total volume in the mixing vessel) and by an estimated total output flowrate.

5. The method of claim 3, further comprising computing a difference between an estimated total volumetric disturbance flowrate and the volumetric disturbance flowrate of the material to obtain an estimated volumetric disturbance flowrate for other materials in the mixing vessel.

6. The method of claim 5, further comprising estimating the output flowrate of a second material exiting the mixing vessel.

7. The method of claim 6, further comprising negatively feeding back the estimated output flowrate of the second material.

8. The method of claim 7, further comprising computing a summation of the estimated volumetric disturbance flowrate of the other materials, a commanded input flowrate of the second material, and the negatively fed estimated output flowrate of the second material, thereby obtaining an estimation of a volumetric rate of change of the second material in the mixing vessel.

9. The method of claim 8, further comprising integrating the estimated volumetric rate of change of the second material to compute an estimated volumetric ratio of the second material to the total materials in the mixing vessel.

10. The method of claim 9, further comprising dynamically recomputing the estimated output flowrate of the second material.

11. The method of claim 9, wherein said integrating the estimated volumetric rate of change of the second material is performed by a second feedback control loop.

12. The method of claim 6, wherein the material comprises water and the second material comprises cement.

13. The method of claim 1, wherein the total materials comprise a liquid and a gas-transported dry material.

14. The method of claim 1, wherein the volumetric disturbance flowrate of the material is determined by comparing a measured input flowrate of the material to the commanded input flowrate of the material.

15. A method of operating a computer system for a first and second mixing vessel to determine an estimated volumetric ratio of a material to total materials in the second mixing vessel that is partially separated from the first mixing vessel, comprising:
summing a commanded input flowrate of the material and a volumetric disturbance flowrate of the materials being fed to the first mixing vessel;
estimating an output flowrate of the material exiting the first mixing vessel based on the summed commanded input flowrate and the volumetric disturbance flowrate;
negatively feeding back the estimated output flowrate of the material to obtain an estimated volumetric rate of change of the material in the first mixing vessel;
integrating the estimated volumetric rate of change of the material in the first mixing vessel to dynamically recompute the estimated output flowrate of the material exiting the first mixing vessel;
estimating an output flowrate of the material exiting the second mixing vessel based on the recomputed estimated output flowrate of the material exiting the first mixing vessel;
negatively feeding back the estimated output flowrate of the material exiting the second mixing vessel and summing it with the estimated output flowrate of the material exiting the first mixing vessel, thereby obtaining an estimated volumetric rate of change of the material in the second mixing vessel; and
integrating the estimated volumetric rate of change of the material in the second mixing vessel to compute the estimated volumetric ratio of the material to the total materials in the second mixing vessel;
wherein the mixing process computer system controls the mixing process within each respective mixing vessel volumetrically utilizing at least one of the respective estimated volumetric ratios of the materials to control at least one input valve.

16. The method of claim 15, further comprising measuring an output flowrate of the total materials exiting the second mixing vessel.

17. The method of claim 15, further comprising dynamically recomputing the estimated output flowrate of the material exiting the second mixing vessel.

18. The method of claim 15, wherein said integrating the estimated volumetric rate of change of the material in the first mixing vessel and said integrating the estimated volumetric rate of change of the material in the second mixing vessel are performed by respective first and the second feedback control loops.

19. A method of operating a computer control system to determine an estimated volumetric ratio of a second material to total materials in a first mixing vessel that is partially separated from a second mixing vessel, comprising:
measuring a height of the total materials in the second mixing vessel;
comparing the height of the total materials in the second mixing vessel to a summation of an estimated height of a first material in the second mixing vessel and an estimated height of the second material in the second mixing vessel to obtain an estimation of a height error for the second mixing vessel;
feeding the estimated height error to a controller to compute an estimated volumetric disturbance flowrate of the second material;
computing a summation of (a) a commanded input flowrate of the second material to the first mixing vessel, (b) the estimated volumetric disturbance flowrate of the second material, and (c) a negative value of an estimated output flowrate of the second material from the first mixing vessel, thereby obtaining an estimated volumetric rate of change of the second material in the first mixing vessel; and
integrating the estimated volumetric rate of change of the second material to obtain the estimated volumetric ratio of the second material to total materials in the first mixing vessel;
wherein the mixing process computer control system controls the mixing process within each respective mixing vessel volumetrically utilizing the estimated volumetric ratio of the second material to control at least one input valve.