US10860753B2ActiveUtilityA1

Characterization of fluids with drag reducing additives in a couette device

44
Assignee: SCHLUMBERGER TECHNOLOGY CORPPriority: Nov 7, 2013Filed: Nov 7, 2013Granted: Dec 8, 2020
Est. expiryNov 7, 2033(~7.3 yrs left)· nominal 20-yr term from priority
G06F 2113/14G06F 30/20G06F 30/28G01N 11/04G01N 2011/006G01F 1/28G06F 17/11
44
PatentIndex Score
0
Cited by
49
References
20
Claims

Abstract

A method is provided for characterizing fluid flow in a pipe where the fluid includes a drag reducing polymer of a particular type and particular concentration. A computational model is configured to model flow of a fluid in a pipe. The computational model utilizes an empirical parameter for a drag reducing polymer of the particular type and the particular concentration. The computational model can be used to derive information that characterizes the flow of the fluid in the pipe. The empirical parameter for the particular type and the particular concentration of the drag reducing polymer can be identified by solving another computational model that is configured to model turbulent Couette flow in a Couette device for a fluid that includes a drag reducing polymer of the particular type and the particular concentration. The empirical data needed for identification of the empirical parameter are obtained from Couette device experiments.

Claims

exact text as granted — not AI-modified
What is claimed is: 
     
       1. A method of characterizing fluid flow in a pipe where the fluid includes a drag reducing polymer of a particular type and particular concentration, the method comprising:
 i) providing a Couette device that includes an annulus disposed between an outer annular surface and an inner annular surface and a sensor disposed adjacent the outer annular surface, wherein the sensor is configured to measure shear stress of fluid at the outer annular surface; 
 ii) providing a first computational model associated with the Couette device of i) to model turbulent Couette flow of a fluid based upon a flow representation consisting of an inner sublayer, an outer sublayer, and a turbulent core between the inner sublayer and the outer sublayer, wherein the first computational model includes an empirical parameter for the particular type and the particular concentration of the drag reducing polymer that is based on the shear stress of the fluid at the outer annular surface measured by the sensor, and wherein the first computational model further includes a first drag reduction parameter associated with the outer sublayer and a second drag reduction parameter associated with the inner sublayer, wherein both the first and second drag reduction parameters are based on the empirical parameter; 
 iii) loading a fluid into the annulus of the Couette device of i), the fluid including the particular type and the particular concentration of the drag reducing polymer, and operating the Couette device of i) to induce turbulent flow of the fluid in the annulus of the Couette device and to obtain related experimental data including shear stress of the fluid at the outer annular surface measured by the sensor; 
 iv) processing the first computational model of ii) in conjunction with the experimental data obtained in iii) on a data processor to solve for the empirical parameter of the first computational model; 
 v) providing a second computational model that is configured to model flow of fluid in a pipe, wherein the second computational model includes a drag reduction parameter that is a function of the empirical parameter as solved for in iv), and wherein the second computational model further includes a friction factor that is a function of the drag reduction parameter, wherein the friction factor relates pressure loss due to friction along a given length of pipe to the mean flow velocity through the pipe; and 
 vi) processing the second computational model on a data processor to derive information that characterizes the flow of fluid in the pipe. 
 
     
     
       2. A method according to  claim 1 , wherein the drag reduction parameter is also a function of a dimensionless pipe radius. 
     
     
       3. A method according to  claim 2 , wherein the second computational model is configured to relate the drag reduction parameter to the empirical parameter by an equation of the form:
     D   * =1+α *   R   + 
 
 where
 D *  is the drag reduction parameter, 
 α *  is the empirical parameter, and 
 R +  is the dimensionless pipe radius. 
 
 
     
     
       4. A method according to  claim 1 , wherein the second computational model is configured to relate the friction factor to the drag reduction parameter by an equation of the form: 
       
         
           
             
               
                 1 
                 
                   f 
                   0.5 
                 
               
               = 
               
                 
                   4 
                   ⁢ 
                   
                     
                       log 
                       10 
                     
                     ⁡ 
                     
                       ( 
                       
                         Re 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           f 
                           0.5 
                         
                       
                       ) 
                     
                   
                 
                 + 
                 
                   8.2 
                   ⁢ 
                   
                     D 
                     * 
                     2 
                   
                 
                 - 
                 8.6 
                 - 
                 
                   12.2 
                   ⁢ 
                   
                     log 
                     10 
                   
                   ⁢ 
                   
                     D 
                     * 
                   
                 
               
             
           
         
         where
 f is the friction factor, 
 D *  is the drag reduction parameter, and 
 Re is the Reynolds number of the flow in the pipe. 
 
       
     
     
       5. A method according to  claim 4 , wherein the second computational model is further configured to relate the Reynolds number Re to a dimensionless pipe radius R + . 
     
     
       6. A method according to  claim 4 , wherein the information derived in iv) includes a solution for the friction factor f for given flow conditions. 
     
     
       7. A method according to  claim 6 , wherein the information derived in iv) includes a pressure drop over a given length of pipe based on the solution for the friction factor f. 
     
     
       8. A method according to  claim 1 , wherein the second computational model is based upon a representation of the flow as two layers consisting of a viscous outer sublayer that surrounds a turbulent core. 
     
     
       9. A method according to  claim 1 , wherein both the first and second drag reduction parameters are also based on a dimensionless torque G derived from the shear stress of the fluid at the outer annular surface measured by the sensor. 
     
     
       10. A method according to  claim 9 , wherein:
 the outer and inner annular surfaces of the Couette device of i) are concentric with respect to one another about a common center, wherein the outer annular surface is offset from the center by a first radius R and the inner annular surface is offset from the center by a second radius r 0 , wherein R is greater than r 0 ; and 
 the first computational model is configured to relate the first and second drag reduction parameters to the empirical parameter by equations of the following form: 
 
       
         
           
             
               
                 D 
                 
                   o 
                   * 
                 
               
               = 
               
                 1 
                 + 
                 
                   
                     
                       α 
                       * 
                     
                     2 
                   
                   ⁢ 
                   
                     ( 
                     
                       1 
                       - 
                       η 
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       G 
                       
                         2 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         π 
                       
                     
                   
                 
               
             
           
         
         
           
             
               
                 D 
                 
                   i 
                   * 
                 
               
               = 
               
                 1 
                 + 
                 
                   
                     
                       α 
                       * 
                     
                     2 
                   
                   ⁢ 
                   
                     
                       ( 
                       
                         1 
                         - 
                         η 
                       
                       ) 
                     
                     η 
                   
                   ⁢ 
                   
                     
                       G 
                       
                         2 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         π 
                       
                     
                   
                 
               
             
           
         
         where
 D 0*  is the drag reduction parameter associated with the outer sublayer, 
 D i*  is the drag reduction parameter associated with the inner sublayer, 
 α *  is the empirical parameter, 
 η is the ratio r 0 /R, and 
 G is the dimensionless torque derived from the shear stress of the fluid at the outer annular surface measured by the sensor. 
 
       
     
     
       11. A method according to  claim 1 , wherein the first computational model includes an equation that defines a fluid velocity at a boundary of a viscous sublayer adjacent one of the outer and inner annular surfaces. 
     
     
       12. A method according to  claim 11 , wherein the equation is derived by momentum conservation for a portion of a turbulent core adjacent the viscous sublayer. 
     
     
       13. A method according to  claim 1 , wherein:
 the sensor comprises a floating element, a mechanical cantilever beam and a fiber Bragg grating strain gauge. 
 
     
     
       14. A method according to  claim 9 , wherein:
 the dimensionless torque G is determined based upon the following relation: 
 
       
         
           
             
               G 
               = 
               
                 T 
                 
                   ρ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     v 
                     2 
                   
                   ⁢ 
                   L 
                 
               
             
           
         
         
           where
 T is a torque derived from a shear stress τ w  as T=τ w 2πR 2 L, where the shear stress τ w  is the shear stress of the fluid at the outer annular surface measured by the sensor, R is the radius of the outer annular surface, 
 ρ is the density of the fluid, 
 ν is the kinematic viscosity of the fluid, and 
 L is the height of the gap of the Couette device. 
 
         
       
     
     
       15. A method according to  claim 1 , wherein:
 the method is repeated for a number of different concentrations of a particular drag reducing polymer or over different drag reducing polymers to characterize the expected pipe flow for these different scenarios. 
 
     
     
       16. A method according to  claim 1 , wherein:
 the method is repeated for a number of different flow conditions to characterize the expected pipe flow for these different scenarios. 
 
     
     
       17. A method of characterizing fluid flow in a pipe where the fluid includes a drag reducing polymer of a particular type and particular concentration, the method comprising:
 i) providing a Couette device that includes an annulus disposed between an outer annular surface and an inner annular surface; 
 ii) providing a first computational model associated with the Couette device of i) to model turbulent Couette flow of a fluid, wherein the first computational model includes an empirical parameter for the particular type and the particular concentration of the drag reducing polymer; 
 iii) loading a fluid into the annulus of the Couette device of i), the fluid including the particular type and the particular concentration of the drag reducing polymer, and operating the Couette device of i) to induce turbulent flow of the fluid in the annulus of the Couette device and to obtain related experimental data; 
 iv) processing the first computational model of ii) in conjunction with the experimental data obtained in iii) on a data processor to solve for the empirical parameter; 
 v) providing a second computational model that is configured to model flow of fluid in a pipe, wherein the second computational model is configured to utilize the empirical parameter as solved for in iv); and 
 vi) processing the second computational model on a data processor to derive information that characterizes the flow of fluid in the pipe; 
 wherein the first computational model includes a first drag reduction parameter associated with the outer annular surface of the Couette device of i) and a second drag reduction parameter associated with the inner annular surface Couette device of i), wherein both the first and second drag reduction parameters are based on the empirical parameter and a dimensionless torque G applied to the Couette device; 
 wherein the outer and inner annular surfaces of the Couette device of i) are concentric with respect to one another about a common center, wherein the outer annular surface is offset from the center by a first radius R and the inner annular surface is offset from the center by a second radius r 0 , wherein R is greater than r 0 ; and 
 wherein the first computational model is configured to relate the first and second drag reduction parameters to the empirical parameter by equations of the following form: 
 
       
         
           
             
               
                 D 
                 
                   o 
                   * 
                 
               
               = 
               
                 1 
                 + 
                 
                   
                     
                       α 
                       * 
                     
                     2 
                   
                   ⁢ 
                   
                     ( 
                     
                       1 
                       - 
                       η 
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       G 
                       
                         2 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         π 
                       
                     
                   
                 
               
             
           
         
         
           
             
               
                 D 
                 
                   i 
                   * 
                 
               
               = 
               
                 1 
                 + 
                 
                   
                     
                       α 
                       * 
                     
                     2 
                   
                   ⁢ 
                   
                     
                       ( 
                       
                         1 
                         - 
                         η 
                       
                       ) 
                     
                     η 
                   
                   ⁢ 
                   
                     
                       G 
                       
                         2 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         π 
                       
                     
                   
                 
               
             
           
         
         where
 D 0*  is the drag reduction parameter associated with the outer annular surface, 
 D i*  is the drag reduction parameter associated with the inner annular surface, 
 α *  is the empirical parameter, 
 η is the ratio r 0 /R, and 
 G is the dimensionless torque applied to the Couette device. 
 
       
     
     
       18. A method according to  claim 17 , wherein:
 the Couette device of i) includes a sensor disposed adjacent the outer annular surface, wherein the sensor is configured to measure shear stress of fluid at the outer annular surface; 
 the empirical parameter of the first computation model is based on the shear stress of the fluid at the outer annular surface measured by the sensor; and 
 the experimental data obtained in iii) and processed in iv) includes the shear stress of the fluid at the outer annular surface measured by the sensor. 
 
     
     
       19. A method according to  claim 18 , wherein:
 the dimensionless torque G is determined based upon the following relation 
 
       
         
           
             
               
                 G 
                 = 
                 
                   T 
                   
                     ρ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       v 
                       2 
                     
                     ⁢ 
                     L 
                   
                 
               
               , 
             
           
         
         where T is a torque derived from a shear stress τ w  as T=τ w 2πR 2 L, where the shear stress τ w  is the shear stress of the fluid at the outer annular surface measured by the sensor, R is the radius of the outer annular surface, L is the height of the gap of the Couette device, ρ is the density of the fluid, and ν is the kinematic viscosity of the fluid. 
       
     
     
       20. A method according to  claim 1 , wherein:
 the sensor is located in a sensor enclosure disposed adjacent the outer annular surface.

Cited by (0)

No later patents cite this yet.

References (0)

No backward citations on record.