Not applicable.
Not applicable.
The present disclosure is directed to a control system design for a mixing system having multiple inputs, and more particularly, but not by way of limitation, to a control system design for a cement mixing system used in well bore servicing applications.
A control system typically comprises one or more physical system components under some form of automated control that cooperate to achieve a set of common objectives. The control system may be designed to reliably control the physical system components in the presence of external disturbances, variations among physical components due to manufacturing tolerances, and changes in commanded input values for controlled output values, such as a cement mixture density, for example. The control system may also be designed to remain stable and avoid oscillations within a range of specific operating conditions.
In a well bore environment, a control system may be used when mixing materials to achieve a desired mixture output. For example, when drilling an oil or gas well, it is common to install a tubular casing into the well bore and cement the casing in place against the well bore wall. A cement mixing system that supports well bore servicing operations, such as cementing casing into a well bore, may be designed with a control system configured to provide a desired volumetric flow rate of mixed cement having a desired density. In particular, the cement mixing control system may control valves that allow the in-flow of dry cement material and water to obtain the desired cement mixture density and desired cement mixture volumetric flow rate. The control system may operate, for example, by monitoring the cement mixture flow rate and density, and by regulating an in-flow water valve and an in-flow dry cement material valve. However, because such systems conventionally control the output parameters, such as cement mixture flow rate and density, dependently, these systems tend to have long lag times in the response of one valve to changes in the position of the other valve. This can lead to unacceptable oscillations in the monitored parameters, and difficulty in stabilizing the system. Therefore, to make the system more stable, it would be desirable to control output parameters, such as a mixture flow rate and a mixture density, for example, independently of each other. Accordingly, a need exists for a mixing control system with multiple inputs that decouples the effects of changes in the commanded outputs.
Disclosed herein is a control system for mixing at least two materials in a physical system having two or more tanks comprising at least two actuators, each actuator being operable to introduce a material into a first tank to form a first mixture, the first mixture flowing into a second tank to form a second mixture and a controller operable, based on a commanded input, to control the at least two actuators to obtain a density of either the first mixture or the second mixture and a volume flow rate of the second mixture out of the second tank, wherein the density is controlled independently from the volume flow rate.
Further disclosed herein is a control system for mixing at least two materials in a tank comprising at least two actuators, each actuator being operable to introduce a material into the tank to form a mixture and a controller operable, based on a commanded input, to control the at least two actuators to obtain a density of the mixture and a volume flow rate of the mixture out of the tank, wherein the density is controlled independently from the volume flow rate.
Further disclosed herein is a control system for mixing at least two materials comprising a controller operable to control a physical system comprising a plurality of actuators, each actuator being operable to introduce a material into a first tank to form a first mixture, wherein the controller operates the actuators to control a desired characteristic of the first mixture.
Further disclosed herein is a method for mixing at least two materials comprising introducing the at least two materials into a physical system, combining the at least two materials to form a mixture thereof, and independently controlling a desired characteristic of the mixture of the at least two materials and a desired parameter of the physical system.
Further disclosed herein is a method for mixing at least two materials comprising determining a commanded combined mass flow rate of the at least two materials into a first tank [commanded dmin/dt], determining a commanded combined volume flow rate of the at least two materials into the first tank [commanded dvin/dt], and independently controlling, based on the commanded dmin/dt and the commanded dvin/dt, a density of a mixture of the at least two materials in the first tank [ρ12] and a volume flow rate of the mixture out of the first tank [dv12/dt].
These and other features and advantages will be more clearly understood from the following detailed description taken in conjunction with the accompanying drawings and claims.
For a more complete understanding of the present disclosure and the advantages thereof, reference is now made to the following brief description, taken in connection with the accompanying drawings and detailed description, wherein like reference numerals represent like parts.
It should be understood at the outset that the present disclosure describes various implementations of different embodiments of a control system having one or more inputs. However, the present control system may also be implemented using any number of other techniques, whether currently known or in existence. The present disclosure should in no way be limited to the descriptions, drawings, and techniques illustrated below, including the design and implementation illustrated and described herein.
In an embodiment shown in
In an embodiment, the physical plant 10 is a well bore servicing fluid mixing system, such as, for example, a cement mixer used to provide a continuous stream of a cement slurry for cementing a tubular casing against a well bore wall. In this embodiment, a dry cement material may be fluidized by the introduction of pressurized air, which promotes fluid flow of the dry cement through a first feed line 28, to be dispensed into the first tank 12 through the first actuator 18, for example, and a carrier fluid, such as water, for example, may flow through a second feed line 30 to be dispensed into the first tank 12 through the second actuator 20. In this embodiment, the first and second actuators 18, 20 may be valves, for example. These two materials, the dry cement material and the carrier fluid, are mixed by the first stirrer 24 to obtain the first mixture 13. Non-fluidized sand or other particulate matter may be dispensed through a third actuator (not shown), such as a screw feeder, for example, into the first tank 12 to mix with the cement slurry. In various embodiments, any of the actuators may comprise a valve, a screw feeder, an augur, an elevator, or other type of actuator known to those skilled in the art. The first mixture 13 is preferably substantially homogenous. The physical plant 10 should provide the cement slurry via the outflow pump 22 and the discharge line 32 at a volumetric flow rate sufficient to support the well bore servicing operation, and should mix the dry concrete material, the carrier fluid, and any particulate material in appropriate proportions so that the cement slurry dispensed thereby achieves a desired density. The physical plant 10 in other embodiments may support other mixing operations of other materials. For example, in another embodiment, proppants and a carrier fluid may be dispensed through the first and second actuators 18, 20 into the first tank 10 to form a portion of a fracturing fluid.
The first control system 100 disclosed hereinafter is expected to reduce control oscillations due to system time lags and to promote independent control of a mixture density and a mixture flow rate.
Turning now to
The first control system 100 comprises a controller 102 that receives input parameter values from an operator through an interface with the first control system 100, and also receives sensed parameter values from sensors coupled to or integral with the physical plant 10. The controller 102 distributes as output parameter values commands to the first actuator 18 and the second actuator 20. An input parameter value 104 for h2 provides to the controller 102 the desired height h2 of the second mixture 15 in the second tank 16, an input parameter value 106 for dvs/dt provides the desired volumetric flow rate dvs/dt of the second mixture 15 out of the second tank 16, and an input parameter value 108 for ρs provides the desired density ρs of the second mixture 15 out of the second tank 16. A h1 sensor 110 provides an indication of the height h1 of the first mixture 13 in the first tank 12, a h2 sensor 112 provides an indication of the height h2 of the second mixture 15 in the second tank 16, a ρ12 sensor 114 provides an indication of the density ρ12 of the first mixture 13, and a ρs sensor 116 provides an indication of the density ρs of the second mixture 15. These indications may be referred to as sensed parameter values.
In an embodiment, the first control system 100 controls the actuators 18, 20 to achieve a desired volumetric flow rate dv12/dt of the first mixture 13 over the weir 14 into the second tank 16, and a density ρs of the second mixture 15 out of the second tank 16, independently of each other. For example, changing the dvs/dt input parameter value 106 causes the controller 102 to control the actuators 18, 20 so that the actual volumetric flow rate dv12/dt of the first mixture 13 over the weir 14 into the second tank 16 changes until it substantially equals the dvs/dt input parameter value 106. Actual values may also be referred to as nominal values in some contexts by those skilled in the control systems art. However, the density ρs of the second mixture 15 as it leaves the second tank 16 remains substantially unchanged because the density and flow rate parameters are controlled independently. Similarly, changing the ρs input parameter value 108 causes the controller 102 to alter the control signals to actuators 18, 20 until the sensed density ρs of the second mixture 15 as read by the ρs sensor 116 substantially equals the ρs input parameter value 108. However, the volumetric flow rate dv12/dt of the first mixture 13 over the weir 14 into the second tank 16 remains substantially unchanged because the flow rate and density parameters are controlled independently.
Turning now to
In an embodiment, the flow modulator 150 generates an actuator1 signal to control the first actuator 18, and this actuator1 signal may be expressed mathematically as proportional to a function ƒ1 as follows:
actuator1 signal ∝ƒ1=[dmin/dt−(dvin/dt)ρm2]/(ρm1−ρm2) (1)
Similarly, the flow modulator 150 generates an actuator2 signal to control the second actuator 20, and the actuator2 signal may be expressed mathematically as proportional to a function ƒ2 as follows:
actuator2 signal ∝ƒ2=[(dvin/dt)ρm1−dmin/dt]/(ρm1−ρm2) (2)
where dmin/dt is the combined mass flow rate entering the tank, dvin/dt is the combined volume flow rate entering the tank, ρm1 is the density of the first material, for example dry cement, flowing into the first tank 12 from the first actuator 18, and where ρm2 is the density of the second material, for example water, flowing into the first tank 12 from the second actuator 20. One skilled in the art will readily be able to determine a suitable first constant of proportionality for equation (1) to drive the first actuator 18 and a suitable second constant of proportionality for equation (2) to drive the second actuator 20 in a particular embodiment. In an embodiment, the actuator1 and actuator2 signals may be conditioned by one or more components (not shown) between the flow modulator 150 and the actuators 18, 20 to conform the actuator1 and actuator2 signals to a non-linear response of one or more of the actuators 18, 20.
Turning now to
The second flow regulator 204 receives from another component of the controller 102 a commanded mass flow rate signal 203 for dm12/dt of the first mixture 13 over the weir 14 into the second tank 16, a sensed height h1 of the first mixture 13 in the first tank 12 from the h1 sensor 110, and a sensed density ρ12 of the first mixture 13 from sensor 114. The second flow regulator 202 may receive the commanded mass flow rate dm12/dt signal 203, for example, from a flow2 state feedback controller with command feed forward 252, to be discussed hereinafter, or from some other component of the controller 102. The second flow regulator 204 generates the commanded first mass flow rate dmin/dt signal 153 that is received as a commanded value by the flow modulator 150 or by the first and second flow modulators 152, 154 described above with reference to
In an embodiment, the function of the first flow regulator 202 may be expressed mathematically as:
commanded dvin/dt=F(h1)(1−Kv)+Kv(commanded dv12/dt) (3)
and the function of the second flow regulator 204 may be expressed mathematically as:
commanded dmin/dt=F(h1)ρ12(1−Km)+Km(commanded dm12/dt) (4)
where F(h1) is a non-linear function of the height h1 of the first mixture 13 in the first tank 12 and provides an estimate of a volumetric flow rate dv12/dt of the first mixture 13 over the weir 14 into the second tank 16; where commanded dv12/dt is the commanded volumetric flow rate dv12/dt signal 201; where commanded dm12/dt is the commanded mass flow rate dm12/dt signal 203; where ρ12 is an indication of density of the first mixture 13 in the first tank 12 based on the input from the ρ12 sensor 114; and where Kv and Km are constants of proportionality that are greater than zero. The values for Kv and Km may be chosen through a closed form solution and/or iteratively so as to minimize overall response time while maintaining stability and desired flow rate and density trajectories during transition phases. An exemplary non-linear function F(h1) for volumetric flow rate over a rectangular weir is given in Engineering Fluid Mechanics, 5th Edition, by Roberson and Crowe, published by Houghton Mifflin, 1993, and may be represented as:
dv12/dt=F(h1)=KL(h1−hw)(3/2) (5)
where L is the length of the weir 14 dividing the tanks 12, 16, K is a flow coefficient that may be empirically determined over a set of operating conditions for a specific weir geometry, and hw is a constant representing the height of the weir.
The volumetric outflow of the first mixture 13 dv12/dt and the mass outflow of the first mixture 13 dm12/dt from the first tank 12 to the second tank 16 may be modeled as negative state feedbacks in the physical plant 10. In equation (3), the effect of the 1×F(h1) term is to cancel or decouple the negative state feedback associated with dv12/dt, and in equation (4), the effect of the 1×ρ12F(h1) term is to decouple the negative state feedback associated with dm12/dt. The control system 102 may be more robust as a result of the state feedback decoupling in the first and second flow regulators 202, 204 because the control system 102 only has to correct for errors between desired and actual mass rate and desired and actual volumetric flow rate without having to also compensate for the first mixture 13 leaving the first tank 12.
Turning now to
The flow2 state feedback controller with command feed forward 252 receives the dvs/dt input parameter value 106 and the ρs input parameter value 108, from an operator interfacing with the first control system 100, and also receives an indication of the density ρs of the second mixture 15 in the second tank 16 from the ρs sensor 116. Based on these inputs, the flow2 state feedback controller with command feed forward 252 produces the commanded mass flow rate dm12/dt signal 203 that is received as a commanded value by the second flow regulator 204 described above with reference to
Turning now to
commanded dv12/dt=u1(t)=dvs/dt+Kp1e1(t)+Ki1∫e1(t)dt (6)
The outflow of the second mixture 15 from the second tank 16 may be modeled as negative state feedback in the physical plant 10. In equation (6), the effect of the dvs/dt term, which is the command feed forward term, is to decouple the negative state feedback associated with the outflow of the second mixture 15 from the second tank 16. While this may not be strict state feedback decoupling, the same benefits of robustness may be at least partially obtained using this technique.
Turning now to
commanded dm12/dt=u2(t)=inputρs{(dh2/dt)A2+dvs/dt}+Kp2e2(t)+Ki2∫e2(t)dt (7)
The outflow of mass from the second tank 16 may be modeled as negative state feedback in the physical plant 10. In equation (7), the effect of the ρs(dvs/dt) term, the command feed forward term, is to decouple the negative state feedback associated with the outflow of mass from the second tank 16. While this may not be strict state feedback decoupling, the same benefits of robustness may be at least partially obtained using this technique. An additional refinement of this relaxed state feedback decoupling technique may be obtained by decoupling the effect of a loss of mass associated with changes of height of the second mixture 15 in the second tank 16. In equation (7), the ρs(dh2/dt)A2 term also contributes to decoupling the negative state feedback associated with the outflow of mass from the second tank 16. The dh2/dt factor may be determined from a series of indications of the height of the second mixture 15 in the second tank 16 or by other means such as an estimate of dh2/dt produced by a height observer as discussed below.
One skilled in the art will recognize that the results of the analysis of the flow1 state feedback controller with command feed forward 250 and the flow2 state feedback controller with command feed forward 252 above may be applied to digital signals as well as analog signals. For example, analog parameters such as the indication of density ρs of the second mixture 15 in the second tank 16 from the ρs sensor 116 may be converted by an analog-to-digital converter (D/A converter) to a digital signal. Similarly, analog outputs may be produced by a digital-to-analog converter (A/D converter), optionally combined with an amplifier to provide sufficient power to drive an electromechanical device, converting a digital control signal to an analog control signal suitable for controlling the first actuator 18 and the second actuator 20.
Turning now to
Referring to the right side of the drawing, the first and second flow modulators 152, 154, as described above with reference to
The first flow regulator 202 receives the commanded volumetric flow rate dv12/dt signal 201 from the flow1 state feedback controller with command feed forward 250 and the h1 indication from the h1 sensor 110 to generate the commanded first volume flow rate dvin/dt signal 151 that feeds into the flow modulators 152,154. As discussed above with reference to
The second flow regulator 204 receives the commanded mass flow rate dm12/dt signal 203 from the flow2 state feedback controller with command feed forward 252, the indication from the h1 sensor 110, and the ρ12 indication from the ρ12 sensor 114 to generate the commanded first mass flow rate dmin/dt signal 153 that feeds into the flow modulators 152, 154. As discussed above with reference to
The flow1 state feedback controller with command feed forward 250 receives the h2 input parameter value 104 and the dvs/dt input parameter value 106 from an operator interfacing with the first control system 100, and also receives the indication of h2 from the h2 sensor 112 to generate the commanded volumetric flow rate dv12/dt signal 201 that feeds into the first flow regulator 202. As discussed above, the flow1 state feedback controller with command feed forward 250 may be referred to as a height controller, because it controls the height h2 of the second mixture 15 in the second tank 16. The command feed forward term, namely the dvs/dt input parameter value 106, may be considered to decouple the negative state feedback of volumetric flow out of the second tank 16 as discussed above with reference to
The flow2 state feedback controller with command feed forward 252 receives the dvs/dt input parameter value 106 and the ρs input parameter value 108 from an operator interfacing with the first control system 100, and also receives the indication of ρs from the ρs sensor 116 to generate the commanded mass flow rate dm12/dt signal 203 that feeds into the second flow regulator 204. As discussed above, the flow2 state feedback controller with command feed forward 252 may be referred to as a density controller, because it controls the density ρs of the second mixture 15 in the second tank 16. The command feed forward term, formed by the product of the ρs input parameter value 108 multiplied by the dvs/dt input parameter value 106, may be considered to decouple the negative state feedback of mass flow out of the second tank 16 as discussed above with reference to
The flow1 and flow2 state feedback controller with command feed forward 250, 252 may be said to provide yet another level of removal from the system parameters dvin/dt and dmin/dt directly associated with the states of the first and second actuators 18, 20, and provide the desired responsiveness to the parameters desired to be controlled.
It will be readily appreciated by one skilled in the art that portions of the control components described above may be combined, for example within a computer program implementing the functional blocks of the control components.
In operation, as set out in the logic flow diagram of
The process begins at block 400 in which the well bore servicing equipment, including the physical plant 10 and the first control system 100, is brought to the well bore site and assembled. The discharge line 32 of
Once the equipment has been set up at the well site, the process proceeds to block 402 where an operator provides input parameter values to the controller 102, for example through a console or lap top computer coupled to the controller 102, for example via a serial cable or a wireless link. The input parameter values may include the h2 input 104, the dvs/dt input 106, and the ρs input 108. In operation, the controller 102 will act to control the first and second actuators 18, 20, for example valves, such that the corresponding actual quantities of ρs, dvs/dt, and h2 approach or equal the input parameter values.
The process proceeds to block 404 where the operator engages the controller 102. The process continues hereinafter along two independent but at least partially coupled paths. Proceeding to block 406, the controller 102 actively controls the first and second actuators 18, 20 in accordance with the input parameter values and the sensed conditions of the physical plant 10. When the controller 102 is first engaged, it is likely that the first and second tanks 12, 16 will be empty. In this case, the controller 102 may open one of the first or second actuators 18, 20 fully open and then open the other of the first or the second actuators 18, 20 so as to achieve the desired cement slurry density as indicated by the ρs input parameter value 108. The controller 102 continually determines and updates the actuator control signals that are output to the first and second actuators 18, 20, and the controller 102 may be said to be operating within a control loop represented by blocks 406 and 408. The controller 102 remains in this control loop 406, 408 while compensating for changes of indications returned from the physical plant 10, for example the sensed density ρs of the cement slurry in the second tank 16 according to sensor 116, and for changes of input values, for example the h2 input parameter value 104.
The process also concurrently proceeds along a second branch to block 410 where the physical plant 10 begins to pump cement slurry via the outflow pump 22, through the discharge line 32, through the coupled pipes, and downhole into the well bore. As the flow rate dvs/dt of the cement slurry exiting the second tank 16 via the outflow pump 22 changes, for example due to fluctuations in electrical power driving the outflow pump 22 or due to fluctuations in well bore back pressure, the controller 102 adjusts and maintains the actual physical parameters of the physical plant 10 to approach or equal the input parameter values.
The process proceeds to block 412 where the operator may modify an input parameter value, for example the ρs input 108. This change causes the controller 102 to change the control signals that are output to the first and second actuators 18, 20 in the control loop 406, 408, for example by further opening the first actuator 18 and further closing the second actuator 20. The coupling between this action of modifying an input parameter value in block 412 and the control loop in blocks 406, 408 is indicated by a dotted line in
The process proceeds to block 414 where the operator stops the controller 102. This action in block 414 will affect the control loop in blocks 406, 408 as indicated by a dashed line in
In an embodiment, the indication of the height h1 of the first mixture 13 in the first tank 12, and the indication of the height h2 of the second mixture 15 in the second tank 16 are provided by two height observer components, which estimate rather than directly sense h1 and h2. In another embodiment, a single height observer may be employed to provide an indication of the height h2 of the second mixture 15 in the second tank 16. In yet another embodiment, a single height observer may be employed to provide an indication of the height h1 of the first mixture 13 in the first tank 12, for example in a physical plant 10 having only one tank.
Under field conditions, the height indications h1, h2 provided by the h1 sensor 110 and the h2 sensor 112 may be subject to various errors, for example height oscillations due to movement of the physical plant 10 onboard a floating platform or ship. Additionally, the stirring action of the first stirrer 24 and second stirrer 26 may significantly agitate the level surface of the first mixture 13 in the first tank 12 and of the second mixture 15 in the second tank 16, introducing variations in the height indications h1 and h2. The introduction of the first and second materials into the first tank 12, for example dry cement and/or water, may introduce further variations in the height indication h1 provided by the h1 sensor 110. All of these surface height variations may be analyzed as noise in the height signals. It may be desirable to employ estimated height indications of h1 and h2 rather than propagate the noise or oscillation that may be present in the indications of direct sensors, for example the h1 sensor 110 and the h2 sensor 112, into the controller 102.
Generally, a height observer is implemented as a dynamic control system to obtain an estimated height of the mixture in the tank in real time. First, this estimate of the mixture height is compared to the measured mixture height to obtain a height error. Then this mixture height error is used to drive the estimated mixture height to an actual mixture height through the use of a proportional-integral type controller, also referred to as a PI controller. By setting the gains of the PI controller, the noise and oscillations of the mixture in the tank can be substantially removed from the mixture height estimate while tracking the actual value of the mixture height. The height observer according to the present disclosure reduces the negative effects of noise and poor sensor performance due to environmental effects such as cement dust in the air or tank oscillations from the height readings. This height observer reflects the state of the actual mixture height with substantially no time lag.
Turning now to
u3(t)=Kp3e3(t)+Ki3∫e3(t)dt (8)
The output of the height PI controller component 272 may be employed as an estimate of system errors, for example the difference between the desired and actual volumetric flow rates as well as other parameter estimations, for example the cross sectional area A1 of the first tank 12. These errors may be collected and referred to as a disturbance, where the disturbance accounts for the discrepancies between what is demanded from the first control system 100 and what the first control system 100 actually produces, unmeasured quantities such as air flows into or out of the subject tank, and inaccuracies of estimates of system parameters. When the physical plant 10 mixes dry cement, there may be differences between the desired cement rate and the actual cement rate, because the cement delivery may be inconsistent. Differences between the desired cement rate and the actual cement rate may be accommodated by the disturbance. The output of the height PI controller component 272 may be referred to as a first disturbance estimate 281. The first disturbance estimate 281 may be fed back into the controller 102 to provide disturbance decoupling. Disturbance decoupling is described in more detail hereinafter.
The output of the height PI controller component 272 (i.e., the output from the sixth summation component 280) is positively summed by a seventh summation component 282 with one or more height feed forward inputs 283, for example volumetric flow rates such as the difference between the commanded volumetric flow rate dvin/dt into the first tank 12 and the volumetric flow rate dv12/dt of cement slurry out of the first tank 12 over the weir 14 into the second tank 16. The performance of the height observer 270 is expected to be improved by using height feed forward inputs 283 as compared to the performance of more traditional controllers, which only employ feedback terms with no feed forward inputs. Because the height of a mixture in a tank, for example the height h1 of the first mixture 13 in the first tank 12, will increase or decrease due to a net positive or negative volumetric flow rate into the tank, the estimate of the height of the mixture in the tank depends upon the net positive or negative volumetric flow rate into the tank: the height feed forward inputs 283. The height feed forward inputs 283 may be summed either negatively or positively at the seventh summation component 282. The output of the seventh summation component 282 conforms to a volumetric flow rate that may be represented generally as dv/dt.
The height tank model component 274 multiplies the volumetric flow rate output of the seventh summation component 282 by a fourth integral component 284 associated with a fourth integration factor 286. The fourth integral component 284 corresponds to the inverse of the cross-sectional area of the tank, represented by “1/A” in the block for the fourth integral component 284. Multiplying a volumetric flow rate term (dv/dt) by the inverse of an area term, for example the inverse of a cross-sectional tank area expressed as 1/A, results in a velocity term (dx/dt), or more particularly in the case of a height controller, an estimated height rate of change dh/dt 285. The intermediate result of the estimated height rate of change dh/dt 285 may be used by other components in the controller 102, for example by the flow2 state feedback controller with command feed forward 252.
Integrating the velocity term via the fourth integration factor 286 results in a displacement x, or in the present case a height h. Thus, the output of the height tank model component 274 is an estimated height 277 of the mixture in the tank. One skilled in the art will recognize that the results of the analysis of the height observer 270 above may be applied to digital signals as well as analog signals. For example, analog parameters such as the sensed height input 271 may be converted by an analog-to-digital converter (A/D converter) to a digital signal. Similarly, analog outputs may be produced by a digital-to-analog converter (D/A converter), optionally combined with an amplifier to provide sufficient power to drive an electromechanical device, converting a digital control signal to an analog control signal suitable for controlling the first actuator 18 and the second actuator 20. Height observers, such as the height observer 270, are further disclosed in related 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 Jan. 4, 2005, which is incorporated herein by reference for all purposes.
Turning now to
The h2 sensor 112 provides an indication of the height h2 of the second mixture 15 in the second tank 16 to the second height observer 270-b, which is negatively summed by a twelfth summation component 269-b with a second estimated height negative feedback term 273-b of the second mixture 15 in the second tank 16 to obtain an error term e3-b(t) that is fed into a second height PI controller component 272-b. The output of the second height PI controller component 272-b is positively summed by a thirteenth summation component 282-b with a second height feed forward input 283-b output by a first calculation component 268. The first calculation component 268 outputs the result of the function F(h1) based on the estimated height 277-a of the first mixture 13 in the first tank 12. In an embodiment, the value of F(h1) is calculated by the first flow regulator 202 and provided to the thirteenth summation component 282-b as the second height feed forward input 283-b. In this embodiment, the first calculation component 268 is not employed. The output of the thirteenth summation component 282-b is then processed by a second height tank model component 274-b to produce a second estimated height 277-b of the second mixture 15 in the second tank 16. The first proportional factor 284 (i.e., 1/A) for the second height tank model component 274-b employs the cross-sectional area A2 of the second tank 16. Note that the second estimated height 277-b and the second estimated height negative feedback term 273-b are the same signals but are identified by different labels to point out their different functions in the controller 102.
Within the second height observer 270-b, the output from the second height PI controller component 272-b provides the second disturbance estimate signal 281-b that is negatively summed by the eleventh summation component 279, as previously discussed. In an embodiment, the second disturbance estimate signal 281-b may provide a more accurate estimate of the volumetric rate disturbance in the system 100 as compared to the first disturbance estimate signal 281-a because the height of the second mixture 15 in the second tank 16 varies more than the height of the first mixture 13 in the first tank 12. Additionally, in the two tank system of
The advantages of the height observer 270, such as, for example, attenuation of noise and improvement of poor sensor performance, estimation of a disturbance term, and estimation of a parameter rate of change, may be obtained in a density observer having a structure related to the height observer 270. In an embodiment, the indication of the density ρ12 of the first mixture 13 in the first tank 12 and the indication of the density ρs of the second mixture 15 in the second tank 16 are provided by two density observer components which estimate rather than directly sense the densities ρ12, ρs of the mixtures 13, 15, respectively. In another embodiment, the density observer may be used to estimate the density of other mixtures or materials in systems other than the physical plant 10 depicted in
Turning now to
The output of the density PI controller component 292 may be employed as an estimate of system errors, for example the difference between the desired and actual mass flow rates as well as other parameter estimations. These errors may be collected and referred to as a disturbance, where the disturbance accounts for the discrepancies between what is demanded from the first control system 100 and what the first control system 100 actually produces, unmeasured quantities such as air flows into or out of the subject tank, and inaccuracies of estimates of system parameters. The output of the density PI controller component 292 may be referred to as a second disturbance estimate 301. The second disturbance estimate 301 may be fed back into the controller 102 to provide disturbance decoupling. Disturbance decoupling is described in more detail hereinafter.
The output of the density PI controller component 292, which conforms to a mass flow rate, is summed with one or more density feed forward inputs 293, for example a mass flow rate, such as the difference between the commanded mass flow rate into the first tank 12 and the mass flow rate of cement slurry out of the first tank 12 over the weir 14 into the second tank 16, by a third summation component 302. The output of a twelfth multiplication component 435 and the output of a thirteenth multiplication component 440 are also negatively summed by the third summation component 302. Generally, the density feed forward inputs 293 are associated with mass flow into the associated tank and the outputs of the twelfth and thirteenth multiplication components 435, 440 are associated with mass flow out of the associated tank. The output of the third summation component 302 is processed by the density tank model component 294. The density tank model component 294 multiplies the output of the third summation component 302 by a sixth integral component 304 associated with a sixth integration factor 306. The sixth integral component 304 is inversely proportional to the height of the mixture times the cross-sectional area of the tank, as represented by “1/hA” in the block for the sixth integral component 304. Note that dividing a mass flow rate by hA is substantially equivalent to dividing through by the volume of the tank resulting in an estimated density rate of change 307 dρ/dt. The height may be provided by a height sensor, for example the h1 sensor 110, or the height observer 270. Alternatively, the height may be a fixed constant determined by experimentation to provide a preferred response rate of the general density observer 290. Integrating this quotient with respect to time results in a density. The output of the density tank model component 294 is thus the estimated density 297 of the mixture in the tank. The estimated density negative feedback term 295 is the same signal as the estimated density 297, but the estimated density feedback term 295 is fed back into the eighth summation component 299 at the input to the density observer 290 to be processed through the PI controller 292 and tank model component 294 to yield a more accurate estimated density 297 each time through. The estimated density negative feedback term 295 is multiplied by a factor A(dh/dt) by the twelfth multiplication component 435 to produce a ρA(dh/dt) term that is negatively fed back to the third summation component 302 as described above. The ρA(dh/dt) term corresponds to a mass rate of change based on changes in the height of the mixture in the tank. The dh/dt factor may be determined from a height sensor or a height observer. Alternatively, in an embodiment of the system 100 that provides no indication or estimate of height of the mixture, the twelfth multiplication 435 may be absent from the general density observer 290. The estimated density negative feedback term 295 is also multiplied by a factor dv/dt by the thirteenth multiplier component 440 to produce a ρ(dv/dt) term that is negatively fed back to the third summation component 302 as described above. The ρ(dv/dt) term corresponds to a mass rate of change due to flow of the mixture out of the tank. In an embodiment, the sensed density input 291 may be provided by a density sensor, for example the ρ12 sensor 114, installed in a recirculation line or an outflow line, such as the discharge line 32, associated with the tank. Thus, the density of the mixture as measured in the recirculation line may lag several seconds behind the density of the mixture in the tank. In this embodiment, a time delay of several seconds, such as three seconds, for example, may be introduced in the estimated density negative feedback term 295 before negatively summing with the sensed density 291 at the eighth summation component 299 to determine the fourth error term e4(t). This allows the sensed and estimated density to be in the same time reference frame before determining the fourth error term e4(t). The appropriate time delay may be readily determined from experimentation by one skilled in the art and depends upon the viscosity of the mixture and the speed of mixing.
Because the structures of the height observer 270 and the density observer 290 are related, one skilled in the art need only determine the gains appropriate to the height observer 270, determine the gains appropriate to the density observer 290, and configure the two observer structures 270, 290 accordingly. One skilled in the art will recognize that the results of the analysis of the density observer 290 above may be applied to digital signals as well as analog signals. For example, analog parameters such as the indication of density ρs of the second mixture 15 in the second tank 16 from the ρs sensor 116 may be converted by an analog-to-digital converter (A/D converter) to a digital signal. Similarly, analog outputs may be produced by a digital-to-analog converter (D/A converter), optionally combined with an amplifier to provide sufficient power to drive an electromechanical device, converting a digital control signal to an analog control signal suitable for controlling the first actuator 18 and the second actuator 20.
Turning now to
The output from the first density PI controller component 292-a, as the second disturbance estimate 301, is negatively summed by the seventeenth summation component 305 with the output from the second flow regulator 204 plus the output generated by a fifth multiplier component 321 to determine the commanded first mass flow rate dmin/dt signal 153. The second disturbance estimate 301 provides disturbance decoupling to the determination of the commanded first mass flow rate dmin/dt signal 153, while the fifth multiplier component 321 provides state feedback decoupling. The fifth multiplier component 321 outputs the product formed by multiplying the first estimated density 297-a of the first mixture 13 in the first tank 12 by the height rate of change dh1/dt of the first mixture 13 in the first tank 12 and by the area A1 of the first tank 12, obtaining the height rate of change dh1/dt, for example, from the estimated height rate of change 285 output of the first height observer 270a. In another embodiment, the second disturbance estimate 301 may be provided by the second density observer 290-b.
To estimate the density ρs of the second mixture 15, the ρs sensor 116 provides an indication of the density ρs of the second mixture 15 in the second tank 16 to the second density observer 290-b. This indication of ρs is summed by an eighteenth summation component 299-b with a second estimated density negative feedback 295-b of the second mixture 15 in the second tank 16 to obtain an error term e4-b(t). As described above, in an embodiment where the ρs sensor 116 is located in a recirculation line or in an outflow line, such as the discharge line 32, the estimated density negative feedback term 295-b may be delayed several seconds before it is input to the eighteenth summation component 299-b. The error term e4-b(t) is fed into a second density PI controller component 292-b. A second density feed forward input 293-b, output by a sixth multiplier component 322, is positively summed by a nineteenth summation component 302-b with the output of the second density PI controller component 292-b. The sixth multiplier component 322 outputs the product formed by multiplying the first estimated density 297-a of the first mixture 13 in the first tank 12 by the volumetric flow rate dv12/dt of the first mixture 13 into the second tank 16. The output of a sixteenth multiplication component 435-b and a seventeenth multiplication component 440-b, corresponding to a mass rate of change due to changes in height in the second tank 16 and mass flow rate out of the second tank 16, respectively, are negatively summed by the nineteenth summation component 302-b. The output of the nineteenth summation component 302-b is processed by a second density tank model component 294-b. Note that the sixth integral component 304 for the second density tank model component 294-b employs the cross sectional area A2 of the second tank 16 and the height h2 of the second mixture 15 in the second tank 16. The second estimated density 297-b of the second mixture 15 in the second tank 16 is provided as an input to the flow2 state feedback controller with command feed forward 252 for use in generating the commanded mass flow rate dm12/dt, which is the commanded mass flow rate dm12/dt signal 203 provided as a commanded parameter value to the second flow regulator 204 as described above with reference to
One skilled in the art will readily appreciate that the components of the first control system 100 disclosed above are susceptible to many alternate embodiments. While several alternate embodiments are disclosed hereinafter, other embodiments are contemplated by the present disclosure.
Turning now to
Several components are combined to provide a first height controller 366 that is substantially similar to the flow1 state feedback controller with command feed forward 250 described above with reference to
Turning now to
One of ordinary skill in the art will readily appreciate that no integral processing is employed in the first height controller 366 described above with reference to
Turning now to
The density observer 290 determines the density estimate 297 of the density ρ12 of the first mixture 13 in the first tank 12 and the second disturbance estimate 362 based on the output of the ρ12 sensor 114 and density feed forward inputs 293. Because this embodiment provides no indication of height, the A(dh/dt) term in
A twelfth multiplication component 384 multiplies the density estimate 297 by the approximate actual volumetric flow dv12/dt. Because h1 is not measured, and the function F(h1) may not be used, the dvs/dt input 106 is used to approximate the actual volumetric flow dv12/dt over the weir 14. The output of the twelfth multiplication component 384 provides state feedback decoupling to the control system 376 to compensate for mass leaving the first tank 12.
The second disturbance estimate 362 is negatively summed with the output of the eleventh and twelfth multiplication components 382, 384 by a twenty-first summation component 374 to determine the commanded first mass flow rate dmin/dt signal 153. The commanded first mass flow rate dmin/dt signal 153 is provided to the first and second flow modulators 152, 154 to control the first and second actuators 18, 20 as described above. The output of the twenty-first summation component 374 also provides the density feed forward inputs 293 to the density observer 290.
Turning now to
While the physical plant 10 controlled by the several embodiments of the control systems 100, 364, 376, and 378 described above included either one or two tanks, in other embodiments the control systems may control mixing parameters of more than two tanks connected in series. In other embodiments, one or more sensors may be employed selected from the ρ12 sensor 114, the ρs sensor 116, the h1 sensor 110, the h2 sensor 112, and a flow rate sensor. Embodiments may use a mixture of sensors with no height observer 270 or density observer 290. Other embodiments may use one or more sensors with one or more height observers 270 and/or density observers 290. While the height observer 270 and the density observer 290 are discussed above, the observer concept may be extended to other sensors to reduce noise output from these and other sensors and to obtain the benefits associated with decoupling or removing disturbances and estimating intermediate rate terms. One skilled in the art will readily appreciate that coupling among control components may be altered to simplify processing, as for example combining two summation components in series by directing all the inputs to the two summation components to a single summation component or by replacing two components that calculate the same value by one component and routing the output of the one component to two destinations. Alternately, in an embodiment, a summation component summing more than two inputs may be expanded into a series of two or more summation components that collectively sum the several inputs.
The controller 102 used in the various control systems 100, 364, 376, and 378 described above may be implemented on any general-purpose computer with sufficient processing power, memory resources, and network throughput capability to handle the necessary workload placed upon it.
The secondary storage 584 typically comprises one or more disk drives or tape drives and is used for non-volatile storage of data and as an over-flow data storage device if RAM 588 is not large enough to hold all working data. Secondary storage 584 may be used to store programs that are loaded into RAM 588 when such programs are selected for execution. The ROM 586 is used to store instructions and perhaps data that are read during program execution. ROM 586 is a non-volatile memory device which typically has a small memory capacity relative to the larger memory capacity of secondary storage 584. The RAM 588 is used to store volatile data and perhaps to store instructions. Access to both ROM 586 and RAM 588 is typically faster than to secondary storage 584.
I/O devices 590 may include printers, video monitors, liquid crystal displays (LCDs), touch screen displays, keyboards, keypads, switches, dials, mice, track balls, voice recognizers, card readers, paper tape readers, or other well-known input devices. The network connectivity devices 592 may take the form of modems, modem banks, ethernet cards, universal serial bus (USB) interface cards, serial interfaces, token ring cards, fiber distributed data interface (FDDI) cards, wireless local area network (WLAN) cards, radio transceiver cards such as Global System for Mobile Communications (GSM) radio transceiver cards, and other well-known network devices. These network connectivity devices 592 may enable the CPU 582 to communicate with an Internet or one or more intranets. With such a network connection, it is contemplated that the CPU 582 might receive information from the network, or might output information to the network in the course of performing the above-described method steps. Such information, which is often represented as a sequence of instructions to be executed using processor 582, may be received from and outputted to the network, for example, in the form of a computer data signal embodied in a carrier wave.
Such information, which may include data or instructions to be executed using processor 582 for example, may be received from and outputted to the network, for example, in the form of a computer data baseband signal or signal embodied in a carrier wave. The baseband signal or signal embodied in the carrier wave generated by the network connectivity devices 592 may propagate in or on the surface of electrical conductors, in coaxial cables, in waveguides, in optical media, for example optical fiber, or in the air or free space. The information contained in the baseband signal or signal embedded in the carrier wave may be ordered according to different sequences, as may be desirable for either processing or generating the information or transmitting or receiving the information. The baseband signal or signal embedded in the carrier wave, or other types of signals currently used or hereafter developed, referred to herein as the transmission medium, may be generated according to several methods well known to one skilled in the art.
The processor 582 executes instructions, codes, computer programs, scripts which it accesses from hard disk, floppy disk, optical disk (these various disk based systems may all be considered secondary storage 584), ROM 586, RAM 588, or the network connectivity devices 592.
While several embodiments have been provided in the present disclosure, it should be understood that the disclosed systems and methods may be embodied in many other specific forms without departing from the spirit or scope of the present disclosure. The present examples are to be considered as illustrative and not restrictive, and the intention is not to be limited to the details given herein, but may be modified within the scope of the appended claims along with their full scope of equivalents. For example, the various elements or components may be combined or integrated in another system or certain features may be omitted, or not implemented.
Also, techniques, systems, subsystems and methods described and illustrated in the various embodiments as discrete or separate may be combined or integrated with other systems, modules, techniques, or methods without departing from the scope of the present disclosure. Other items shown or discussed as directly coupled or communicating with each other may be coupled through some interface or device, such that the items may no longer be considered directly coupled to each other but may still be indirectly coupled and in communication, whether electrically, mechanically, or otherwise with one another. Other examples of changes, substitutions, and alterations are ascertainable by one skilled in the art and could be made without departing from the spirit and scope disclosed herein.
The present application claims the benefit under 35 U.S.C. § 119(e) of U.S. Provisional Application Ser. No. 60/671,392 filed Apr. 14, 2005 and entitled “Implementation of Alternative Cement Mixing Control Schemes,” by Jason D. Dykstra, et al, which is incorporated herein by reference for all purposes.
Number | Name | Date | Kind |
---|---|---|---|
3379421 | Putman | Apr 1968 | A |
3591147 | Anderson et al. | Jul 1971 | A |
3605775 | Zaander et al. | Sep 1971 | A |
3886065 | Kappe et al. | May 1975 | A |
3933041 | Hyer | Jan 1976 | A |
4141656 | Mian | Feb 1979 | A |
4327759 | Millis | May 1982 | A |
4349435 | Ochiai | Sep 1982 | A |
4397561 | Strong et al. | Aug 1983 | A |
4421716 | Hench et al. | Dec 1983 | A |
4436431 | Strong et al. | Mar 1984 | A |
4654802 | Davis | Mar 1987 | A |
4764019 | Kaminski et al. | Aug 1988 | A |
4779186 | Handke et al. | Oct 1988 | A |
5027267 | Pitts et al. | Jun 1991 | A |
5038611 | Weldon et al. | Aug 1991 | A |
5098667 | Young et al. | Mar 1992 | A |
5103908 | Allen | Apr 1992 | A |
5114239 | Allen | May 1992 | A |
5281023 | Cedillo et al. | Jan 1994 | A |
5289877 | Naegele et al. | Mar 1994 | A |
5320425 | Stephenson et al. | Jun 1994 | A |
5365435 | Stephenson | Nov 1994 | A |
5382411 | Allen | Jan 1995 | A |
5441340 | Cedillo et al. | Aug 1995 | A |
5452954 | Handke et al. | Sep 1995 | A |
5503473 | Dearing, Sr. et al. | Apr 1996 | A |
5570743 | Padgett et al. | Nov 1996 | A |
5571281 | Allen | Nov 1996 | A |
5590958 | Dearing, Sr. et al. | Jan 1997 | A |
5624182 | Dearing, Sr. et al. | Apr 1997 | A |
5775803 | Montgomery et al. | Jul 1998 | A |
6007227 | Carlson | Dec 1999 | A |
6113256 | Bonissone et al. | Sep 2000 | A |
6120172 | Chen et al. | Sep 2000 | A |
6120173 | Bonissone et al. | Sep 2000 | A |
6253607 | Dau | Jul 2001 | B1 |
6491421 | Rondeau et al. | Dec 2002 | B2 |
7056008 | Rondeau et al. | Jun 2006 | B2 |
7284898 | Duell et al. | Oct 2007 | B2 |
7308379 | Dykstra et al. | Dec 2007 | B2 |
7353874 | Dykstra et al. | Apr 2008 | B2 |
20020093875 | Rondeau et al. | Jul 2002 | A1 |
20030072208 | Rondeau et al. | Apr 2003 | A1 |
20030161211 | Duell et al. | Aug 2003 | A1 |
20040016572 | Wylie et al. | Jan 2004 | A1 |
20050135185 | Duell et al. | Jun 2005 | A1 |
20050201197 | Duell et al. | Sep 2005 | A1 |
20060141107 | Schwimmer et al. | Jun 2006 | A1 |
20060161358 | Dykstra et al. | Jul 2006 | A1 |
20060231259 | Dykstra et al. | Oct 2006 | A1 |
20060233039 | Dykstra et al. | Oct 2006 | A1 |
20060235627 | Dykstra et al. | Oct 2006 | A1 |
20070153622 | Dykstra et al. | Jul 2007 | A1 |
20070153623 | Dykstra et al. | Jul 2007 | A1 |
20070153624 | Dykstra et al. | Jul 2007 | A1 |
20070171765 | Dykstra et al. | Jul 2007 | A1 |
Number | Date | Country |
---|---|---|
495098 | Jul 1992 | EP |
1356188 | Oct 2003 | EP |
9519221 | Jul 1995 | WO |
WO 9519221 | Jul 1995 | WO |
WO 0236959 | May 2002 | WO |
WO 0236929 | May 2002 | WO |
02044517 | Jun 2002 | WO |
WO 03065015 | Aug 2003 | WO |
Number | Date | Country | |
---|---|---|---|
20060233039 A1 | Oct 2006 | US |
Number | Date | Country | |
---|---|---|---|
60671392 | Apr 2005 | US |