CC BY
26
UDC 622.245.4
Development of mathematical models to control the technological properties of cement slurries
Sergei E. CHERNYSHOV1», Vladislav I GALKIN1, Zoya V. ULYANOVA1, David Iain Macpherson MACDONALD2
1 Perm National Research Polytechnic University, Perm, Russia
2
University of Aberdeen, Aberdeen, Scotland, United Kingdom
Oil and gas producing enterprises are making increasingly high demands on well casing quality, including the actual process of injection and displacement of cement slurry, taking into account requirements for the annular cement level, eliminating possible hydraulic fracturing, with developing a hydraulic cementing program. It is necessary to prevent deep invasion of cement slurry filtrate into the formation to exclude bridging of productive layers. It is impossible to fulfill all these requirements at the same time without application of modifying additives; complex cement compositions are being developed and applied more often. Furthermore, need to adjust cement slurries recipes appears for almost every particular well.
In order to select and justify cement slurries recipes and their prompt adjustment, taking into account requirements of well construction project, as well as geological and technical conditions for cementing casing strings, mathematical models of the main technological properties of cement slurries for cementing production casing strings in the Perm Region were developed. Analysis of the effect of polycarboxylic plasticizer (Pl) and a filtration reducer (fluid loss additive) based on hydroxyethyl cellulose (FR) on plastic viscosity (V), spreadability (S) and filtration (F) of cement slurries is conducted. Development of mathematical models is performed according to more than 90 measurements.
Key words: production well; casing cementing; technological properties of cement slurry; modeling the properties of cement slurries; cementing technology; well casing quality
Acknowledgement. This work was financially supported by the Ministry of Education and Science of the Perm Region during the implementation of project MIG-21 (Agreement № C-26/785 dated 21.12.2017).
How to cite this article: Chernyshov S.E., Galkin V.I., Ulyanova Z.V., Macdonald David Iain Macpherson. Development of mathematical models to control the technological properties of cement slurries. Journal of Mining Institute. 2020. Vol. 242, p. 179-190. DOI: 10.31897/PMI.2020.2.179
Introduction. Recipes of cement slurries for casing of well are subject to frequent adjustments due to the need to simultaneously comply with following requirements: slurries must be easily pumped, have rheological properties that can prevent the occurrence of hydraulic fracturing during cementing, have reduced filtration rate, and meet all requirements for the main technological properties [1, 5-7].
Many researchers are studying the effect of various modifying additives on main technological properties of cement slurries. As a result, recipes are developed taking into account the specifics of cementing process conditions on the field and requirements well construction project, including fluid loss additives, plasticizers, defoamers, expanding additives, etc. [1, 4, 6, 7, 9, 12, 20]. Requirements for modern cement slurries are determined by applied technology for well casing, and precise control of their main technological properties determines the possibility of implementing selected cementing method [8, 10, 11, 29].
Much attention of Russian and foreign researchers is paid to rheological properties that determine hydraulic cementing program with regard to high bottomhole temperatures and pressures [5, 13, 15, 21, 22]. Therefore, in addition to actual development of cement slurries recipes, various methods for investigation of properties are proposed, which bring the conditions of laboratory tests closer to real field tests [14, 24, 27].
800
700
600
500
c3 400
Ph
e
> 300
200
100
0
-100
•
• •• • •
•
• •••• • •
• •• • • •
• • s
1 ■ • •
• ;>¡i • ■ :
1 • i
• •••• • •
•••• • •
•
■ i:
i •
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
Pl, %
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
FR, %
420 400 380 360 340 320 300 280 260 240 220 200 180
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
Pl, %
c 1600 1400
1200
1000
.¡3
^ 800
on
3 600
400 200 0
-200
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
Pl, %
Fig. 1. Correlation fields between V and Pl, FR (a); 180
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
FR, %
-0.05 0.05 0.15 0.25 0.35 0.45 0.00 0.10 0.20 0.30 0.40
FR, %
S and Pl, FR (b); F and Pl, FR (c)
Large amount of laboratory test results allows the application of mathematical modeling methods to reduce the number of experiments and to use resulting models for operational management of the main technological properties [17, 18, 23, 26, 30]. However, these models describe one of the technological properties indicators and influence of only one component on any basic composition and are incomplete.
Statement of the problem. The task was to develop mathematical models that describe the change in the parameters of cement slurries depending on the amount of modifying additives introduced. That would allow operational adjustment of cement slurries recipes taking into account the requirements of well construction projects for the main technological properties, as well as geological and technical conditions for casing string cementing.
The following cement slurry recipe was adopted as the basis: Portland cement of PCT-I-G-CC brand; defoamer (DF) based on polyols and esters of fatty acids -0.2 %; water-cement ratio - 0.46. In this recipe, modifying additives were introduced to obtain required indicators of the main technological properties [1, 11].
According to the results of laboratory study, it was noted that defoamer additive does not affect considered and simulated technological properties. Let us introduce restrictions on applicability of the developed mathematical models: temperature conditions 25 ± 4 °C; cement slurry density 1.87 ± 0.02 g/cm3.
Methodology. At the first stage of research, correlation fields
a
b
between V and Pl, FR; S and Pl, FR; F and Pl, FR were constructed (Fig.1) and pair correlation coefficients r between considered modifying additives and studied parameters of technological properties of cement slurries are calculated [3, 16].
Analysis of constructed correlation fields shows that the ratios and strength of the correlation between V, S, F and Pl, FR have different forms.
Configuration of the correlation fields indicates that the plasticizer has a very insignificant effect on the plastic viscosity, while the effect of the filtration reducer on same parameter is significant. Moreover, in concentration range for FR 0-0.225 %, plastic viscosity mainly depends on FR, and at a higher concentration (> 0.250 %), Pl also begins to make some contribution to value of V, without demonstrating a linear relationship and indicating the complexity and possible interdependence of Pl and FR influence on V.Layout of correlation fields indicates an antagonism of Pl and FR by their influence on S and joint unidirectional effect on F.
Calculation of regression coefficients in the developed models was performed with the least squares method. Possibilities of using statistical methods to solve similar problems are given in [3, 16]. Note that the models are built separately for F, S and V.
Let us compare average values of actual F and model FM in the selected intervals by calculating the criterion tp according to following formula
|Xi "X2|
tP =
1 + _l f (n -1)5? + (n - 1)S21
n2
v ni + n2 - 2 y
where X?, X2 - average values of F and FM respectively; S?2, S22 - dispersions.
Difference in average values is considered statistically significant if tp > tt. Values tt were determined depending on amount of compared data and significance level (a = 0.05).
Comparison of the average values showed that for F and FM they are not statistically different in both cases at p < 0.05. For a deeper statistical analysis of the values of F and FM, their distribution was studied. Optimal magnitudes of values intervals for indicators are calculated according to Sturgess formula:
AX = X™ -
1 + 3,32 lg N
where Xmax - maximum indicator value; Xmin - minimum indicator value; N - amount of data. In each interval, frequencies are determined:
P(X) = Nk / Ng,
where P(X) - frequency in k interval for F group (q = 1); q = 2 corresponds to FM group; Nk - number of cases for indicator content P(X) in k interval; Ng - sample size for classes 1 and 2.
When comparing the density of distribution for indicators in the studied classes, Pearson statistics were used:
x2 = nn2 x
1
fM1 M
=1 M1 + M2
v N1 N2 y
where N1, N2 - number of values of F and FM respectively (classes 1, 2); M1, M2 - number of values that fall within a given interval, respectively, for two classes studied; e - number of intervals [3, 16].
Discussion. To build filtration models in different ranges, studied sample was formed according to the principle from minimum (Fmin = 10 cm3/30 min) to maximum (Fmax = 1381 cm3/30 min)
n
values. Number of indicators studied is equal to two, so the first model is built on three data. Next model will be built at n = 4 and so on until n = 91.
As a result, three groups of values were distinguished according to influence of reagents Pl and FR on property F. For the first group, an updated nonlinear model (range F < 80 cm3/30 min):
FM1 = 388 - 1479.9FR - 920.1Pl + 1600.1Pl-FR + 1599.9FR2 + 1000.1Pl2,
at R = 0.999, p < 0 forecast error equals 0.002 cm3/30 min.
Following model is obtained for the second group of values (80 < F < 380 cm3/30 min):
FM2 = 440.81 - 1022.38FR - 757.82Pl,
at R = 0.982, p < 0 forecast error equals 34.015 cm3/30 min.
For the third group model is as follows (F > 380 cm3/30 min):
FM3 = 1221.67 - 5309.76 FR - 5149.04Pl,
at R = 0.893, p < 0 forecast error equals 149.87 cm3/30 min.
Using the above formulas, values FM1, FM2, FM3 were calculated and compared to actual values F. From a practical point of view, interval of F values corresponding to the first class (the smallest values of the filtration rate) is of more interest, since it reflects the requirements of vast majority of regulatory documents on casing cementing for the most critical sections of oil and gas wells. Forecast error is also minimal, which will allow very accurately calculate FM values by known values of Pl and FR. In this case, necessary range of F values can be achieved in the entire range of Pl concentration due to varying the content of FR.
To graphically solve the question on choice of cement slurries recipe, taking into account the requirements for the main technological properties, a graph of FMM dependence on Pl and FR is constructed (Fig.2). Similar graphical solutions were obtained for other technological properties of cement slurries. According to graph (Fig.2), it is possible to determine the concentration values of Pl and FR at which FMM values are minimal and have values that satisfy the requirements of well construction projects.
Analysis of the comparison of model and actual values showed that the non-linear model in selected range of Pl, FR values works much better than other developed models. It should be noted
that even in this case, five values are observed when the values of F and FMM do not significantly coincide. It was found that these discrepancies occur in cases where a zero concentration of FR is used in range of plasticizer concentrations of more than 0.06 %. Identified phenomenon is associated with process of measuring this parameter.
Filtration index is a direct reflection of system ability for filtrate loss additive, implemented in a specific instrument situation - procedure for measuring the parameter. Potential ability of the system for filtrate
0.42 0.40 0.38 0.36 0.34 It 0.32
Ph
0.30 0.28 0.26 0.24 0.22
A I I|N
\\ \ \ \
^ v A \
\\ l:;l )
III A
-
§§
¡¡I ^ a
Üi ÜI
0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 _ 40 Pl, %
— 100 cm3/ 30 min
90 80
— 70 60
— 50
Fig.2. Graph of F dependence on Pl and FR
loss directly depends on amount of free dispersion medium that can be sediment on a filter of a certain diameter at a given pressure drop. Dispersion medium is predominantly solvent molecules (in this case, water), free from participation in formation of solvate shells around cement grains. When measuring the parameter, cement slurry, loaded into the chamber of the filter press, which is a vertical cylinder, undergoes a pressure drop at which the dispersion medium begins to filter out through the filter located below. During the measurement, due to lack of mixing, the system is subject to sedimentation separation, as a result of which heavy particles of cement are lowered towards the filter and, being compressed, form a physical barrier that prevents free filtering of the dispersion medium. Plasticizer molecules adsorbed on surface of cement grains increase mobility of cement particles, especially in conditions of low viscosity of liquid phase, surrounding the particles. Such an increase in mobility directly depends on concentration of the plasticizer and, starting with a certain concentration (in this case, from 0.06 %), it has a significant effect on the process of measuring the filtration index. This effect is compensated if there is a minimum amount of filtration reducer in the system, input of which ensures thickening of aqueous phase in cement slurry that significantly slows down sedimentation process. Mentioned above shows that satisfactory sedimentation stability is achieved at a concentration of FR > 0.05 %.
Thus, the data obtained for a zero concentration of FR should not be used in construction of models, since they do not reflect real state of the system by its ability for filtrate loss and are the result of distortion due to noticeable influence of sedimentation.
Then, using the same methodology, models were built on technological properties of cement slurries S and V.
Comparison of S and SM ratios, in contrast to F and FM ratios, shows that significant non-linearity between them is not observed. According to the principle that was stated above, 83 multidimensional models were built. Analysis of S and SM ratios showed that there are three types of relationships. In the first section there is a significant increase in R values from 0.177 to 0.810 in the range of 190-250 mm. The second section is in the range of 250-300 mm, an insignificant change in the values of R occurs here. In the third section in the range of S > 300 mm, a slight increase in the values of R is observed. Analysis of spreadability models allowed to reveal three types of influence. At values of S < 240 mm, influence of Pl and FR values is almost the same, with S in the range of 240-340 mm, a higher influence of Pl is observed, then inversion of the influence for Pl and FR values is observed, i.e. the influence of FR becomes more significant than Pl.
Thus, according to the influence of Pl and FR reagents on property S, three groups of values were distinguished. A non-linear model is obtained for the first group (range S < 240 mm):
SMM = 220.079 + 594.23Pl - 54.7803FR - 661.3888Pl2 - 644.0422FR-Pl - 81.8865FR2,
at R = 0.885, p < 0 forecast error equals 7.29 mm.
Following model is obtained for the second group (240 < S < 350 mm):
SMM = 271.124 + 836.3447Pl - 464.8676FR - 1326.5561Pl2 - 905.5659FR-Pl + 803.4287FR2,
at R = 0.889, p < 0 forecast error equals 14.66 mm.
For the third group model is as follows (S > 350 mm):
SMM = 348.333 - 916.667FR + 175.000Pl,
at R = 0.998, p < 0 forecast error equals 2.04 mm.
From a practical point of view, interval of S values corresponding to the first class (S < 240 mm) is of more interest, since it reflects the requirements of well construction projects.
- 183
Journal of Mining Institute. 2020. Vol. 242. C. 179-190 • Oil and Gas
Forecast error for this interval of S values (180-240 mm) is 2-4 %, which will make it possible to accurately calculate SM values with known values of Pl and FR. In this case, necessary range of S values can be achieved over the entire range of FR concentrations by varying the content of Pl.
When developing models of plastic viscosity, it was found that formation of V values occurs depending on Pl and FR for two models. At values of V < 240 mPa-s, FR and Pl are involved in formation of values. In other cases, value of V is formed only by FR.
Total of 84 multidimensional models were built. According to influence of Pl and FR on V, two groups of values were distinguished, division into groups occurs mainly in terms of FR reagent concentration.
A model is obtained for the first group (V < 240 mPa-s)
VM1 = 29.294 + 819.8383FR - 32.8641Pl,
at R = 0.989, p < 0 forecast error equals 9.409 mPa-s.
For the second group of values in order to identify areas where non-linearities can be observed, a non-linear model of FR and Pl influence on plastic viscosity was constructed. Results of study are compared with the description of physical and chemical processes occurring during the mixing of cement slurry according to analyzed recipes.
It is known that molecules of cellulose ethers are moderately rigidly chained, and independent position of each molecule (their free orientation) is possible only in relatively dilute solutions. As the number of rigid macromolecules in a given volume of the solution increases, probability of their arbitrary orientation will be less, and when a certain critical concentration is reached, further increase in the number of macromolecules in this volume will be impossible without mutual ordering of their part. Therefore, starting from a certain critical polymer concentration, an additional phase appears in the system, which is represented by ordered macromolecules. With an even greater increase in polymer concentration, fraction of ordered phase will increase until all macromolecules are mutually ordered [2].
For solutions of rigid-chain polymers, following concentration dependence of viscosity is characteristic. First, with an increase in the polymer concentration, curve passes through a maximum corresponding to critical concentration at which the formation of anisotropic (liquid-crystal)
phase begins. Then the curve passes through a minimum corresponding to transition of the system to a fully ordered anisotropic state. Decrease in viscosity is due to layered flow of anisotropic solutions, which has a lower resistance to flow than the mesh network of molecules in isotropic solutions [25].
With a further increase in the polymer concentration, viscosity begins to increase again due to formation of a fluctuation network between the liquid-crystal domains of anisotropic phase [19], as well as a decrease in fraction of free solvent. This character of dependence is also present in cellulose derivatives solutions filled with fine-dispersion particles [28].
Thus, a refined non-linear model is obtained for the second group of values (V > 240 mPa-s):
>
700
600
500
400
300
200
100
-100
□ □
ŒP □ □
CD № nD
s
<S>d
&> aP°
-100 0 100 200
o Class 1
300 400 V, mPa-s
500 600 700 800
Class 2
Fig.3. Correlation field between VMM1, VM1 and V
VMM1 = -43.621 + 2811.9621P1 + 365.5409FR - 8351.5551P12 - 2275.8001FR-P1 + 2564.79FR2,
at R = 0.991, p < 0 forecast error equals 12.964 mPa-s.
Model values of plastic viscosity were calculated and compared with actual values (Fig.3) using mentioned above formulas. It can be seen that values of VMM1, VM1 and V coincide well over the entire range of values. Results of mathematical modeling are confirmed by peculiarities of rheological behavior of cement slurries with different concentrations of FR and Pl reagents.
Analysis of experimental data obtained during the study shows that rheological behavior is preserved in the range of plasticizer concentrations of 0-0.200 % (Fig.4). In this case, maximum, corresponding to beginning of anisotropic phase appearance in the HEC solution, and minimum, corre-
PV, mPa-s 600 500 400 300 200 100 0
0.06 % Pl
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 FR, %
PV, mPa-s 800 700 600 500 400 300 200 100 0
0.08 % Pl
■ f
/
/
r-*
/ r
, /
--—**" '
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 FR, %
PV, mPa-s 600 500 400 300 200 100 0
0.1 % Pl
PV, mPa-s 600 500 400 300 200 100 0
0.2 % Pl
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 FR, %
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 FR, %
PV, mPa-s 600 500 400 300 200 100 0
0.3 % Pl
If || Isotropic area J J Transition area I ^ Anisotropic aret
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 FR, %
Fig.4. Dependence of cement slurry plastic viscosity (PV) on FR concentration at various Pl concentrations Journal of Mining Institute. 2020. Vol. 242. C. 179-190 • Oil and Gas
sponding to complete anisotropization of the solution, shift from 0.200 and 0.225 % HEC at plasti-cizer concentrations <0.06 % towards increasing concentrations of HEC (0.250 and 0.275 %) at plasticizer concentrations > 0.06 %.
This feature is in good agreement with effecting mechanism of plasticizer, which is that when adsorbed on cement grains, it displaces a certain amount of water from their surface, which was bound in hydration shell of the grain, thereby increasing the concentration of free water molecules in the system, which can be involved in hydration and dissolution of HEC molecules. Given that specific surface of cement remains constant, it can be expected that amount of water that can be released from the hydration shell composition during chemisorption of plasticizer molecules on the grain surface is also finite. At plasticizer concentrations > 0.06 %, a certain saturation of cement grain surface is achieved, after which the amount of released water does not change so significantly.
For rheological model, three areas can be arbitrarily distinguished: a completely isotropic area (HEC concentrations of 0-0.200 % at any plasticizer concentration), a two-phase transition area (HEC concentrations from 0.200-0.250 to 0.225-0.275 % depending on plasticizer concentration) and a completely anisotropic area (concentration of HEC 0.275-0.400 % at any concentration of plasticizer), which clearly correlates with results of mathematical modeling.
At a plasticizer concentration of 0.3 %, characteristic maximum and minimum are not identified on the graph. This is probably the result of a noticeable decrease in amount of free solvent -free water becomes involved in hydration of an additional number of plasticizer molecules, which doesn't allow friction between layers to be lowered with increasing anisotropization of solution due to "lubricating" effect of the water interlayer. Poor distinction of maximum and minimum in the graphs for plasticizer concentrations of 0-0.04 % is caused by the same reason, although in this case, free water is involved in formation of hydration layers on the surface of cement grains free of plasticizer.
Noting this remark, in the range of plasticizer concentrations of 0.06-0.2 %, within which the number of free water molecules in the system capable of participating in hydration and dissolution of HEC molecules varies insignificantly, for rheological model, three areas with distinct boundaries can be arbitrary set: a completely isotropic area (HEC concentrations < 0.250 %), a two-phase transition area of appearance and accumulation of a liquid-crystal phase (0.250-0.275 %) and a completely anisotropic area (HEC concentrations > 0.275 %).
In the range of an isotropic HEC solution existence (< 0.250 %), plastic viscosity of the system does not show a dependence on concentration of plasticizer and remains approximately constant. In the range of HEC solution existence in an anisotropic form (> 0.275 %), effect of a plasticizer on viscosity of the system already exists and is difficult to describe.
Presence of free plasticizer molecules that are not associated with the surface of cement grains can cause strengthening in the mesh network between moving layers of anisotropic phase by forming a "bridge" bond mediated by the plasticizer molecule, connected simultaneously with both layers. To realize such a "bridge" bond, one segment of the plasticizer molecule must form a chain of hydrogen bonds with one layer of anisotropic phase, and its other segment with vacant region of other layer. Furthermore, for one plasticizer molecule to be associated simultaneously with two different layers of anisotropic phase, certain conditions must be met. In order to be able to realize the bridging, it is necessary that the region of potential fixation within another layer is not yet occupied by another similar plasticizer molecule and is at an accessible distance.
For example, with an increase in concentration of plasticizer, probability for a site, vacant for attachment within another layer, will be less. Therefore, bridging can be realized at plasticizer concentrations not exceeding a certain critical value. In the field of very low concentrations, such bonds will also be realized, but due to their small numbers, they cannot be identified by increasing the vis-
cosity of the system. On the other hand, with a decrease in the distance between layers of anisotropic phase (i.e., with an increase in HEC concentration), such a bond is easier to realize, and, therefore, this occurs at lower concentrations of plasticizer molecules. Starting at some critical distance, formation of such bonds becomes very unlikely.
It can be expected that described strengthening in the interaction between layers of anisotropic phase will be identified in the form of an increase in viscosity of the system in a rather narrow range of concentrations that are free from interaction of plasticizer molecules with cement. With an increase in HEC concentration, such an increase in viscosity will occur at lower concentrations of plasticizer.
Since the distance between anisotropic layers is a function of HEC concentration, region of bridge bonds' existence will be limited by a certain balance of HEC concentrations and free plasti-cizer molecules and can be described by their ratio, which will actually require additional study outside the scope of the objectives of this work.
From a practical point of view, range of V values corresponding to the first class (V < 240 mPa-s) is of more interest, since it meets the requirements of most hydraulic well cementing programs. Forecast error is not more than 9.5 mPa-s, which will make it possible to quite accurately calculate VM values by known Pl and FR values.
In this case, necessary range of V values can be achieved over the entire range of Pl concentrations at FR concentrations not exceeding 0.225 %. It is positive that in this case it is necessary to deal with the range of FR concentrations within which there is only an isotropic HEC solution, when plastic viscosity of the system does not show a sensitive dependence on Pl presence, which will allow varying the FR concentration relatively freely and independently.
With a higher FR content, achieved values of plastic viscosity are undesirable or unac-ceptably high.
Conclusion. Summarizing the simulation results from standpoint of the developed models' applicability, it should be noted that overlapping concentration range of Pl and FR in the area of working values for parameters S and F (class 1 in both cases) is quite wide, but when the area of required values for parameter V is applied, it decreases sharply.
Expansion of the overlapping field for Pl and FR concentration ranges in the area of operating values for all three parameters at once (S, F, and V) can be achieved by reducing the degree of HEC polymerization incorporated in the recipe. Despite the fact that interaction features for components of the system by objective reasons should be preserved, adjustment of models or, at least, refinement of the boundaries of the selected classes will be required.
Comparing the results of laboratory investigation and mathematical modeling of the main technological properties with a description of physical and chemical processes that occur during mixing of cement slurries compositions in accordance with the developed recipes, as well as proposed methods for analyzing available data used for correlation and regression analyzes, analysis of constructed multidimensional statistical models, both linear and non-linear, allowed quantifying the effect of various modifying additives' values of Pl, FR in recipes of cement slurries and their combinations on filtration index, spreadability and plastic viscosity.
Constructed models can be used in predicting the values of F, S, V taking into account requirements of well construction project, conditions for cementing, estimated values of hydraulic fracturing pressures and absorption when developing a hydraulic program for cementing of wells. Model values of the main technological properties will reduce the number of probable concentrations of modifying additives during development of recipes, which will lead to time and cost savings at the stage of preparatory work for well casing.
Developed mathematical models were used to optimize cement slurries recipes to develop a work program for cementing the production casing of an oil production well at one of the fields in
100 150 200 250 300 350 400
Pressure, kg/cm2
b 'ZZ "nZ 2ZZ r-Z -,ZZ -ZZ
Pressure, kg/cm2
-First batch of displacement fluid Second batch of displacement fluid
-Third batch of displacement fluid -Fourth batch of displacement fluid
-Formation pressure -----Hydraulic fracturing pressure with safety factor
-Hydraulic fracturing pressure
Fig.5. Calculation of open borehole pressures for injection of displacement fluid taking into account values of plastic viscosity for cement slurry made according to design recipe with probable hydraulic fracturing (a), adjusted model formulation without exceeding of fracturing pressure (b)
the Perm Region (Fig.5). Cement slurry was made on the basis of the design recipe, and it turned out that obtained rheological properties can lead to appearance of hydraulic fracturing during implementation of specified modes of cement slurry injection and jet-in (Fig.5, a). In addition, value of filtration rate was higher than required.
Using developed mathematical models, a correction of cement slurry was performed. In fact, defined filtration index of cement slurry, recipe of which was optimized based on obtained model values, did not exceed required design value, and rheological properties corresponded to requirements of the work program and excluded the likelihood of hydraulic fracturing during injection and displacement of cement slurry (Fig.5, b).
Recipe obtained on the basis of developed mathematical models is recommended for use. Laboratory analysis of cement slurry showed the correspondence of model and actual values of the main technological properties for cement slurry. Difference in values of the indicators for simulated technological properties did not exceed 3 %.
The obtained results made it possible to identify peculiarities of physical and chemical interaction of considered modifying additives and their mutual influence in composition of the base cement slurry and determine the limits of their concentrations, beyond which it is not possible or difficult to control the main technological properties of the resulting cement slurries.
REFERENCES
1. Agzamov F.A., Izmukhambetov B.S., Tokunova E.F. Chemistry of cementing and washing solutions. Moscow: Nedra, 2011, p. 245 (in Russian).
2. Bartenev G.M., Frenkel S.Ya. Polymer Physics. Leningrad: Khimiya, 1990, p. 432 (in Russian).
3. Galkin V.I., Kunitskikh A.A. Statistical modeling of expanding cementing composition. Vestnik Permskogo natsionalnogo issledovatelskogo politekhnicheskogo universiteta. Geologiya. Neftegazovoe i gornoe delo. 2017. Vol. 16. N 3, p. 215-244 (in Russian).
4. Litvinenko V.S., Dvoinikov M.V. Justification of the Technological Parameters Choice for Well Drilling by Rotary Steer-able Systems. Journal of Mining Institute. 2019. Vol. 235, p. 24-29. DOI: 10.31897/PMI.2019.1.24
5. Ovchinnikov V.P., Aksenova N.A., Ovchinnikov P.V. Physical and chemical hardening processes, well operation and cement stone corrosion. Neftegazovyi universitet. Tyumen, 2007, p. 369 (in Russian).
6. Pereima A.A., Minchenko Yu.S., Trusov S.G. On influence of chemical treatment of cement slurries on effectiveness of expanding additives. Stroitelstvo neftyanykh i gazovykh skvazhin na sushe i na more. 2011. N 5, p. 2-7 (in Russian).
7. Petrov V.S. Regulation of cement slurry properties - stone using additives of aminomethylene phosphonic complexones. Neftegazovoe delo. 2012. N 6, p. 46-59 (in Russian).
8. Chernyshov S.E. Improving the technology for construction of additional sidetracks from previously drilled wells. Neftyanoe khozyaistvo. 2010. N 6, p. 22-24 (in Russian).
9. Chernyshov S.E., Krapivina T.N. Effect of expanding additives on properties of cement slurry-stone. Vestnik Permskogo go-sudarstvennogo tekhnicheskogo universiteta. Geologiya, geoinformatsionnye sistemy, gorno-neftyanoe delo. 2010. Vol. 9. N 5, p. 3133 (in Russian).
10. Chernyshov S.E., Krysin N.I., Kunitskikh A.A. On improving the technology of production casing cementing at the fields of Western Siberia. Neftyanoe khozyaistvo. 2011. N 10, p. 20-21 (in Russian).
11. Chernyshov S.E., Krysin N.I., Kunitskikh A.A. Results of innovative technology introduction for cementing oil wells at the fields of Kazakhstan. Neftyanoe khozyaistvo. 2012. N 8, p. 108-110 (in Russian).
12. Chernyshov S.E., Kunitskikh A.A., Votinov M.V. Study of hydration dynamics and development of compositions of expanding additives to cement slurries. Neftyanoe khozyaistvo. 2015. N 8, p. 42-44 (in Russian).
13. Al Wahedi Yasser, Bassioni Ghada, Awayes Jafar. Influence of classical and modern superplasticisers on the chemical and rheological behaviour of oil well cement: a comparative study. Advances in Cement Research. 2011. Vol. 23. N 4, p. 175-184. DOI: 10.1680/adcr.2011.23.4.175
14. Zhang S., Li Z., Cheng X., Luo H., Guo X. A modified testing method for low density cement slurry. Drilling Fluid and Completion Fluid. 2014. N 31(3), p. 65-68. DOI: 10.3969/j.issn.1001-5620.2014.03.017
15. Le Saout G., Lecolier E., Rivereau A., Zanni H. Chemical structure of cement aged at normal and elevated temperatures and pressures: Part I. Class G oilwell cement. Cement and Concrete Research. 2006. Vol. 36, p. 71-78. DOI: 10.1016/j.cemconres.2004.09.018
16. Davis J.C. Statistics and data analysis in geology. John Wiley & Sons, 2002, p. 656.
17. Domanska U., Goskowska M. Experimental solid + liquid equilibria and excess molar volumes of alkanol + hexylamine mixtures: Analysis in terms of the ERAS, DISQUAC and Mod. UNIFAC models. Fluid Phase Equilibria. 2004. Vol. 216, p. 135-145.
18. Ghada Bassioni, Mohammed Mohammed Ali, Ali Almansoori, Gabriele Raudaschl-Sieber, Fritz E. Kühn. Rapid determination of complex oil well cement properties using mathematical models. RSCAdvances. 2017. Iss. 9. DOI: 10.1039/C6RA26045D
19. Hermans J. The viscosity of concentrated solutions of rigid rod-like molecules Poly (y-benzyl-L-glutamate) in m-cresol. Colloid Sci. 1962. Vol. 17. N 7, p. 638-648.
20. Abbas G., Irawan S., Kumar S., Elrayah A.A.I. Improving oil well cement slurry performance using hydroxypropylmethylcellu-lose polymer. Advanced Materials Research. 2013. Vol. 787, p. 222-227. DOI: 10.4028/wwwscientific.net/AMR.787.222
21. Memon K.R., Shuker M.T., Tunio S.Q., Lashari A.A., Abbass G. Investigating rheological properties of high performance cement system for Oil wells. Research Journal of Applied Sciences, Engineering and Technology. 2013. N 6(20), p. 3865-3870.
22. Zingg A., Winnefeld F., Holzer L., Pakusch J., Becker S., Figi R. Interaction of polycarboxylate-based superplasticizers with cements containing different C3A amounts. Cement and Concrete Composition. 2009. Vol. 31. Iss. 3, p. 153-162. DOI: 10.1016/j.cemconcomp.2009.01.005
23. Li Y., Nancy A. Bjorndalen, Ergun Kuru. Numerical Modelling of Cuttings Transport in Horizontal Wells Using Conventional Drilling Fluids. Journal of Canadian Petroleum Technology. 2007. Vol. 46. Iss. 7. DOI: 10.2118/07-07-TN
24. Plank J., Sachsenhauser B. Experimental determination of the effective anionic charge density of polycarboxylate superplasticizers in cement pore solution. Cement and Concrete Research. 2009. Vol. 39. Iss. 1, p. 1-5. DOI: 10.1016/j.cemconres.2008.09.001
25. Robinson C., Ward J.C., Beevers R.B. Liquid crystalline structure in polypeptide solutions. Disc. Faraday Soc. 1958. Vol. 25, p. 29-42.
26. Rocha S.A., Guirardello R. An approach to calculate solid-liquid phase equilibrium for binary mixtures. Fluid Phase Equilibria. 2009. Vol. 281, p. 12-21.
27. Roncero J., Valls S., Gettu R. Study of the influence of superplasticizers on the hydration of cement paste using nuclear magnetic resonance and X-ray diffraction techniques. Cement and Concrete Research. 2002. Vol. 32, p. 103-108. DOI: 10.1016/S0008-8846(01 )00636-6
28. Kulichikhin V.G., Makarova V.V., Tolstykh M.Yu., Picken S.J., Mendes E. Structural Evolution of Liquid Crystalline Solutions of Hydroxypropyl Cellulose and Hydroxypropy l Cellulose - Based Nanocomposites during Flow. Polym. Sci. A. 2011. Vol. 53. N 9, p. 748-764.
29. Jin J., Xie W., Zhang H., Wang G., Zhou Z. Technology for liner cementing in Well 1 of storage Bai-6, Bannan, Dagang. Drilling Fluid and Completion Fluid. 2014. Vol. 31. N 6, p. 58-61. DOI: 10.3969/j.issn.1001-5620.2014.06.016
30. Umeokafor C.V., Joel O.F. Modeling of Cement Thickening Time at High Temperatures with Different Retarder Concentrations. 2010. SPE-136973-MS. DOI: 10.2118/136973-MS
Authors: Sergei E. Chernyshov, Candidate of Engineering Sciences, Associate Professor, nirgnf@bk.ru (Perm National Research Polytechnic University, Perm, Russia), Vladislav I. Galkin, Doctor of Geological and Mineralogical Sciences, Professor, vgalkin@pstu.ru (Perm National Research Polytechnic University, Perm, Russia), Zoya V. Ulyanova, Senior Lecturer, z.v.uljanova@pstu.ru (Perm National Research Polytechnic University, Perm, Russia), David Iain Macpherson Macdonald, Professor, d.macdonald@abdn.ac.uk (University of Aberdeen, Aberdeen, Scotland, United Kingdom).
The paper was received on 9 July, 2019.
The paper was accepted for publication on 8 November, 2019.
