ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 9, 0 Studies of the Rheological Properties of Drillig Fluids Iurii Salyzhy, Mykhaylo Myslyuk Ivao-Frakivsk Natioal Techical Uiversity of Oil ad Gas, Ivao-Frakivsk, Ukraie ABSTRACT The methodology ad results of ivestigatios of rheological properties of drillig fluids o the basis of rotatioal viscosimetry data have bee described. Special attetio has bee paid to rheological properties of biviscosity fluids ad applied aspects of determiatio of state equatios. INTRODUCTION Techological processes i oil ad gas idustry are closely coected with the use of o-newtoia fluids. Efficiet maagemet of these processes requires precise determiatio of rheological parameters ad ifluece of ay factors upo them. For this purpose the rotatioal viscosimeters are most widely used but the simplified methodologies for viscosimetry data treatmet do ot allow usig all possibilities of the measurig equipmet. The Couette flow i a gap betwee two coaxial cyliders of rotatioal viscosimeters is described by equatio ω = τ τout γ ξ dξ ɺ () ξ where ω is agular speed of rotatio of outer cylider, τ ad τ out are shear stresses o ier ad outer cyliders, γɺ ( τ) is fluid rheological model, γɺ is shear rate gradiet. The depedece betwee stresses o the outer ad ier cyliders is determied by the followig relatio: α τ τ τ α, if 0 ; τ out = τ 0 τ 0 τ < τ 0 α, if, () where τ 0 is yield stress, α = R R, ad R ad R are the radii of the ier ad outer cyliders. I Eq. the rheologically statioary models are used which permits clear γ ɺ = γɺ τ. aalytical solutio of METHODOLOGY The proposed methodology for rotatioal viscosimetry data treatmet is based o exact solutio of the Eq. ad cosiders iformatioal relevace of the experimets. The aim of rotatioal viscosimetry data processig is formalized as the task of search of ˆ idex of rheological model ad its parameters â i some a priori kow class of ϑ models ( ω ) either A,a + ε;... ; τ = or A ( ω,a ) + ε, ϑ;..., (3)
where τ is vector of measured shear stresses at agular rates ω, A ( ω,a ) is operatio of direct task of rotatioal viscosimetry Eqs. ad, a is vector of rheological properties, is idex of rheological model; ε is vector of cetered ormal radom compoet, caused by measuremet iaccuracies. The solutio of the task Eq. 3 is built usig the priciple of maximum of likelihood fuctio ad is realized with the help of the followig procedures ( ( ω ) τ) ( ) mi C A,a a ˆ, ˆ ; (4) mi σ ˆ = N r N ω τ i= ( A ( i,aˆ ) i ) ˆ, (5) where N is umber of rates of cylider rotatio, r is umber of rheological parameters beig evaluated, C is matrix of covariatios of radom compoet i Eq. 3. Accordig to Eq. 4 at first for each model ϑ evaluatios â of rheological properties are built, ad the accordig to Eq. 5 rheological model ˆ is evaluated. Cotrary to other existig methods, methodology allows more exact selectio of the proper rheological model ad precise estimatio of its parameters. Methodology allows to select the most appropriate model i the class of the rheologicaly statioary models: Newto, Bigham, Ostwald, Herschel-Bulkley, Schulma-Casso ad biviscosity models, estimate their rheological properties ad build the estimatio error covariace matrix O. ( T ˆ ) ˆ ˆ ˆ O = A ω,a C A ω,a, (6) ˆ ˆ ω,a is matrix of derivatives of rheological parameters ad trasposed to it T is matrix A ( ω,aˆ ˆ ). where A The biviscosity model 3 is cosidered as ay combiatio of rheological models, which describes rheological properties of fluid i differet rages of shear rates. ( γ ) * ɺ ( τ,a ), τ τ ; γ ɺ = γ ( ) * ɺ ( τ,a ), τ > τ, (7) where a, a are rheological properties accordigly for the low ad high gradiets * of shear rates, τ is boudary shear stress, calculated from the equatio * * γɺ τ,a = γɺ τ,a. For the selectio of the most appropriate rheological model, which describes fluid i the etire techological process (uder the ifluece of parameters of state or ay additives) the procedure of batch processig is proposed 4. I this case Eqs. 4 ad 5 are geeralized as follows: ( ) mi A ω,a τ a ˆ, =, z ; (8) mi σ ˆ c = z( N r ) ω τ = i= ( A ( i,aˆ ) i ) ˆ z N ˆ, (9) where σ c is dispersio of adequacy for the whole umber of experimets z with differet values of ifluecig factors.
The polyomial ad splie models are used for the describig of ay factors ifluece o the rheological properties. The splie fuctios are represeted accordig to V.A. Vasyleko 5 z 3 α p( x) = b Gm,z x x + b z+ ( x ), = = where b, b z + are parameters of aalytical splie represetatio, z is umber of experimetal poits, G ( x x ) = m,z m z x x l x x, if z is eve; = m z x x, if z is odd, q z x x = ( xi x i ), i= q = z + m! m!z!, x 3 = ( x, x,, x ) T z is vector of ifluece factors i the experimet, α = α, α,, α is multiidex, z α α α αz = z ( x ) (x ) (x ) (x ) ad m is parameter of variatioal fuctioal. THE STUDY As it was metioed i methodology descriptio, biviscosity models have bee icluded ito the class of rheological models, amog which the most adequate oe is beig selected. First-priority task was to fid out whether fluids beig described by these models do exist. The search of such Table. Results of rotatioal viscometry data treatmet where most appropriate is biviscosity rhelogical model Studied fluid Device gap Most appropriate biviscosity rhelogical model Model parameters Boudary shear stress Dispersio σ ˆ, Pa Herschel- Dispersio for other models σ ˆ, Pa shear τ0, k, η, rates Pa Pa s Pa s τ, Pa Bulkley Ostwald Bigham Bigham low.04 0.994 0.79384 & Ostwald high 0.6555 0.767.08 0.008 0.0655 0.0573 0.9857 low.4680 0.489 0.9365 Ostwald & Newto high 0.0379.8 0.7674.5660 5.9580.7950 low 0.30 0.457 3 0.9365 Ostwald & Bigham high.775 0.0083 4.94 0.057 0.077 0.8 0.4 low.590 0.600 4 0.9365 Ostwald & Ostwald high 0.3690 0.560 6.0 0.0907 0.090 0.67 0.5939 Note. τ 0 is yield poit, k is cosistecy idex, η is plastic viscosity ad is flow-behavior idex. Figure. Rheogram for studied fluid Figure. Rheogram for studied fluid
Figure 3. Rheogram for studied fluid 3 Figure 4. Rheogram for studied fluid 4 fluids was carried out by treatmet of rheometry data, obtaied by our colleagues ad us ad take from literature. I table ad o figures 4 the most vivid results for differet biviscosity models, which are as follows: studied fluid is CMC water solutio (Sample 33 6 ), studied fluid is cemet slurry with alumia-silicate microspheres, studied fluid 3 ad 4 are betoite ligite suspesios (Sample 5 ad Sample 7 respectively 6 ). It should be metioed that biviscosity fluids amog drillig fluids occur rather frequetly, but the fractio of such fluids is difficult to evaluate while it depeds o the type of fluid. Aother importat issue for the evaluatio of parameters of biviscosity fluids is the umber of rates of viscosimeter, used to make measuremets. The rheological properties of drillig fluids are more sesitive to ifluece of temperature, pressure ad additives the other parameters. So, while selectig drillig fluids formulatio, the attetio must be previously paid o rheological properties especially whe the aim is to select the most thermostable prescriptio. I the case of selectio drillig mud formulatio the thermostability meas: mud properties before ad after heatig must be i specified boudaries ad the parameters chages must be miimal. I table there are rotatioal viscosimetry data, obtaied accordig to the pla of experimet, the aim of which was the selectio of optimal formulatio of humate-biopolimer drillig mud. Table. Rotatioal viscometry data for humate-biopolimer drillig mud Values of factors Agles of tur (degree) at rotatig rates (rpm) Experimet Polypac Duo- before heatig after heatig % UL, % vis, % 6 30 60 00 00 300 600 6 30 60 00 00 300 600 0. 5 50 55 60 69 75 85 5 4 0 6 35 4 5 0.6 0. 4 5 7 5 5 0 3 5 6 3 0. 0. 4 8 5 8 46 4 5 7 0 9 4 0.6 4 59 7 8 0 4 56 34 43 5 58 66 7 9 5 0. 0. 3 4 5 6 8 0 0 0 3 6 0.6 4 56 64 73 87 98 7 7 7 4 30 40 46 60 7 0. 50 65 75 85 05 7 48 40 5 58 65 74 8 9 8 0.6 0. 3 6 8 0 5 4 3 5 7 4 9 3.08 0.35 0.55 6 37 46 56 74 89 6 5 3 7 3 38 44 57 0 0.9 0.35 0.55 5 5 9 36 4 54 7 3 8 7 3 37 7 0.65 0.55 9 8 38 47 65 78 08 8 3 7 34 39 50 7 0.05 0.55 8 6 30 35 44 5 7 7 4 8 6 30 38 3 7 0.35.0 60 77 87 9 47 4 5 58 64 74 8 96 4 7 0.35 0 0 3 5 7 3 0 0 3 5 8 5 7 0.35 0.55 8 6 3 38 49 58 78 8 3 5 3 36 46 Note. CAR is coal-alkali reaget. Device costat is 0.4788 Pa/degree, device gap is 0.9365.
Thermostatig of drillig fluid i each experimet has bee carried out for 4 hours i autoclave (with the capacity of 400 cm 3 ) i roller ove at the temperature of 0 С. As a result of treatmet of preseted data, parameters of rheological models i each poit of experimet pla have bee evaluated ad the most adequate amog them have bee selected (Table 3). Exp. Table 3. Most appropriate rheological models at experimetal poits Most appropriate models ad dispersio σ ˆ, Pa before heatig after heatig Biviscosity Bigham & Ostwald 0.743 Ostwald 0.65 Herschel-Bulkley 0.0349 Ostwald 0.006 3 Schulma- Casso 0.0476 Ostwald 0.0463 4 Herschel-Bulkley.850 Herschel-Bulkley 0.6888 5 Herschel-Bulkley 0.038 Ostwald 0.0403 6 Herschel-Bulkley 0.37 Biviscosity Bigham & Ostwald 0.367 7 Herschel-Bulkley 0.6 Herschel-Bulkley 0.09 8 Herschel-Bulkley 0.08 Herschel-Bulkley 0.039 9 Herschel-Bulkley 0.4406 Herschel-Bulkley 0.097 Biviscosity 0 Herschel-Bulkley 0.000 Ostwald & Ostwald 0.704 Herschel-Bulkley 0.64 Herschel-Bulkley 0.009 Herschel-Bulkley 0.3087 Ostwald 0.0675 3 Herschel-Bulkley.040 Herschel-Bulkley 0.054 4 Ostwald 0.0078 Ostwald 0.0366 5 Herschel-Bulkley 0.004 Ostwald 0.36 But it is clear that it s practically impossible to use several rheological models for the aalysis of ifluece of reagets upo rheological properties of drillig fluid. Thus, accordig to Eqs. 8 ad 9 dispersio of adequacy for the whole umber of experimets was estimated (Table 4). Table 4 cotais estimatios of dispersio of adequacy ust for five rheological models, because for other models it was t possible to make such estimatio while data could ot be described with these models i some experimetal poits. From Table 4 oe ca make a coclusio that i geeral Herschel- Bulkley model is the most adequate for the descriptio of results of experimets, the evaluatios of parameters of which are give i Table 5. Table 4. Geeralized fuctioal for model selectio Dispersio of adequacy σ ˆ c, Pa Rheological model before heatig after heatig Herschel-Bulkley 0.605 0.363 Schulma-Casso.696 -* Ostwald.04 0.3885 Bigham 9.7 5.43 Newto 8 89.5 Note. * - data could ot be described with this model i some experimetal poits. Table 5. Results for Herschel-Bulkley model Model parameters before heatig after heatig Exp. τ0, k, τ0, k, a Pa Pa s Pa Pa s 0.64 0.350 0.575 0.90 0.460 00 0.60 0.0533 0.774 0 0.066 0.7597 47 3.33 0.55 0.738 0.053 0.60 0.576 80 4.65.450 0.4649.06.730 0.45 50 5 0.80 0.036 0.7045 0 5.6 0-7 5.4980 34 6 0.56 3.6700 0.3559 0.63 0.45 8 7 6.6.090 0.46 7.34 6.005 0.538 87 8 0.74 0.0989 0.7589 0.35 0.0834 0.6908 50 9 6.68.4740 0.5007 4.05.7 0.4300 79 0 4.5 0.7883 0.469 0.675 0.3386 34 4.53.0560 0.5439 0.97.745 0.3654 8 7.0 0.3866 0.607 0.66 0.3398 68 3 0.5.7670 0.4074 3.69.05 0.3779 44 4 0.8 0-9 3.36 0. 0-6 5.4470 5 5 5.48 0.845 0.53.6.79 0.3468 9 Rheological properties refer to vector quatities, specified chage of which ca be estimated o the basis of multidimesioal fuctio of probability desity 7 T a = a a O a a 0 0, (0) where a, a are vectors of rheological 0 properties of drillig fluid correspodigly before ad after thermal ifluece.
Figure 5. Ifluece of reagets o the rheological properties of humate-biopolimer drillig mud Accordig to Eq. 0 estimatios of chage of rheological properties i the course of thermostatig have bee built, which are give i Table 5. restrictios ad which allowed estimatig regio of acceptable formulatios. Regio of acceptable formulatios is quite broad, that s why selectio of optimal While experimet was carried out formulatio of drillig fluid requires accordig to orthogoal cetral compositio pla, for the descriptio of its parameters formalizatio of optimizatio criterio. I this case, the aim is to select the most stable secod-degree polyoms have bee used. formulatio i terms of rheological I figure 5 depedecies of rheological parameters chagig, which is esured by parameters upo cocetratios of reagets miimizatio of criterio of Eq. 0. i drillig fluid at the specified CAR value of 3.5 % are represeted by isolies. Aalogical depedecies were built for other properties of drillig fluid, upo Figure 5 represets depedece of criterio of Eq. 0 upo cocetratios of reagets ad optimal drillig fluid formulatio, obtaied due to its use. values of which were imposed iterval
ACKNOWLEDGMENTS A lot of thaks to Dr. Sami Hietala ad his colleagues for the chace to participate i the coferece. REFERENCES. Myslyuk, M.A. (988), "Determiig rheological parameters for a dispersio system by rotatioal viscometry", Izheero-Fizicheskii Zhural, Vol. 54, No. 6, pp. 975-979.. Myslyuk, M.A., Vasylcheko, A.A., Salyzhy, Yu.M., ad Kusturova, E.V. (006), Evaluatio of rheological properties of drillig muds from the rotatioal viscometry data, Stroit. Neft. Gaz. Skvazhi Sushe More, No., pp. 9-33. 3. Myslyuk, M.A., Salyzhy, Yu.M. (008), The evaluatio of biviscosity fluids rheological properties o the basis of rotatioal viscometry data, Neft. Khoz., No., pp. 40-4,. 4. Myslyuk, M.A., Salyzhy, Yu.M. (007), Rotatioal viscometry: ew approaches to data processig, Naft. Gaz. Prom., No. 6, pp. 7-. 5. Vasyleko, V.A. (983) Splie fuctios: theory, algorithms, programs, Nauka, Novosibirsk. 6. Kelessidis, V.C., Maglioe, R. (008) Shear rate correctios for Herschel-Bulkley fluids i Couette geometry, Appl. Rheol., Vol. 8, 3448. 7. Myslyuk, M.A., Vasylcheko, A.A., Salyzhy, Yu.M., Kusturova, E.V. (006), To the selectio of drillig mud coditioig prescriptio cosiderig the heat resistace, Stroit. Neft. Gaz. Skvazhi Sushe More, No. 8, pp. 47-5.