CC BY
19
УДК 539.422
ПРИМЕНЕНИЕ НЕЛОКАЛЬНОЙ МОДЕЛИ
НЕУПРУГИХ ДЕФОРМАЦИЙ В ДИНАМИЧЕСКОЙ ЗАДАЧЕ
РАЗРУШЕНИЯ ГОРНОЙ ПОРОДЫ
Егор Владимирович Вторушин
Новосибирский технологический центр компании «Бейкер Хьюз», 630090, Россия, г. Новосибирск, ул. Академика Кутателадзе, 4а, кандидат физико-математических наук, научный сотрудник, тел. (383)332-94-43, e-mail: [email protected]
Виталий Николаевич Доровский
Новосибирский технологический центр компании «Бейкер Хьюз», 630090, Россия, г. Новосибирск, ул. Академика Кутателадзе, 4а, доктор физико-математических наук, советник по науке, тел. (383)332-94-43, e-mail: [email protected]
При бурении скважин режуще-скалывающими долотами основным фактором является динамическое взаимодействие резца и горной породы. В данной работе изучается процесс резания горной породы резцом, который описывается в рамках нестационарной нелокальной модели неупругих деформаций. Численный метод решения двумерной динамической модельной задачи основан на методе конечного элемента в смешанной постановке. Численное моделирование показало, что в рамках данной модели можно определить как характерный размер скола, так и оптимальную глубину резания, при условии, что известны неклассические модули породы, используемые для описания неупругой среды в рамках новой модели.
Ключевые слова: разрушение горной породы, резание горных пород, неупругость, неевклидовы деформации, метод конечного элемента.
NONLOCAL MODEL OF INELASTIC DEFORMATIONS APPLIED TO DYNAMIC PROBLEM OF ROCK FRACTURING
Egor V. Vtorushin
Novosibirsk Technology Center «Baker Hughes», 630090, Russia, Novosibirsk, 4a Kutateladze St., Ph. D., Scientists, tel. (383)332-94-43, e-mail: [email protected]
Vitaly N. Dorovsky
Novosibirsk Technology Center «Baker Hughes», 630090, Russia, Novosibirsk, 4a Kutateladze St., D. Sc., Technical Adviser, tel. (383)332-94-43, e-mail: [email protected]
The major factor of well drilling with PDC bit is a dynamic interaction of cutting tool and rock. In this paper non-stationary nonlocal model of inelastic deformations is applied to rock cutting with cutting tool problem. The suggested numerical solution of 2-dimensional dynamic problem is based on mixed finite element method. The nonlocal model leads to the appearance of two non-classical material constants that by analogy with classic model, can be called 'inelastic' moduli. Numerical modeling has proven that is possible to predict both chip size and depth of cut within the scope of the suggested model if only the non-classical moduli are defined.
Key words: Rock fracturing, rock cutting, inelasticity, non-Euclidean deformations, mixed finite element method.
The major factor of well drilling with PDC bit process is a dynamic interaction of cutting tool and rock i.e. rock cutting. The difficulties of mathematical description of the rock cutting process, as it is emphasized in Cherepanov [1], arise from the elongated plastic zones around the cutting edge of a cutter and the interaction of multiple fractures generated from cutting. Cutting the rock formation with an apparent shear plane typically occurs under the brittle mode with cuttings formed as rock particles. The size of the rock particles is characteristic for the cutting process. The main experimental results of the cutting process are well studied (e.g. see Borisov [2]) and are typically reduced to describing the key patterns of the transitional process leading to the steady mode, as well as to characteristic features of the steady mode. On the other hand, the steady cutting mode is characterized with the measureable force of resistance to cutting, including cases when it is a function of the linear velocity of the cutter (Cherepanov [1]). Any theoretical model of the rock cutting process under the brittle fracturing mode proposed must confirm the experimental results. Adding new characteristic features of the cutting process does not bring more clarity to the physics of this phenomenon. The main problem of the rock cutting process may be reduced to finding the characteristic scale of formation fracturing under the conditions of a load applied to the cutter (weight on bit and torque on bit), taking into account the features of the stressed state of the formation in the vicinity of the cutter tip. The presented paper is providing an answer to the question using numerical modeling of rock cutting based on non-stationary nonlocal model of inelastic deformations [3].
In the Guzev paper [4], the following fact was noted: accounting for plasticity in the classical nonlinear theory of elasticity takes the deformed field of the stressed state outside of the Euclidean space (in the general case, the deformation metrics tends to become Riemann metrics). For the isotropic systems, such deformations are characterized by nonzero invariants of the Ricci tensor. Under the conditions of the plastic flow, it is the curvature of the deformation space that is the linear invariant. Then Gusev used the curvature as an incompatibility parameter for building stationary solution of the zonal disintegration problem. As a result, a structure of periodical fractured zones observed around mining pits was achieved.
Using the measure of deformation incompatibility as the basis, the definition of R is established as a second-order invariant of the deformation tensor e:
R
-r = Asj/ -djdkeik. Let's introduce destruction y - a conjugated intensive pair
of the extensive incompatibility parameter R in the terms of the non-equilibrium thermodynamics as y = Pe^- + aR. Here a and p are parameters of state equation
e0 = const + T0(s - s0) + Xe2 /2p + \izikeik / p + PRe^ / p + aR2 / 2p,
where p is the density and X, ^ denote elastic moduli. Parameters a and p are considered as new material parameters referred to as 'inelastic' moduli. Let's introduce the stress tensor aik and velocity of the deformed medium v. Then motion equations for small strains with presence of destruction field are expressed as:
P^ - dk^ik = 0 > % - (dkvi + divk)!2 = -%k - tpjfiik / 3 , ^ik = te jfik + + P/?5;/,,
Vjj = £(ajj -4A^), Vik = £ [aik-ajAk/3) + 2(didky~^ik/3)
Here £ and £ are plastic (Maxwell relaxation) moduli. Note that as R and y are conjugated to each other only one of them can be used independently.
In this paper the rock medium' behavior in the vicinity of a rigid moving cutter is analyzed numerically in scope of the nonlocal model of inelastic deformations. The mixed finite element solver presented by E.V. Vtorushin in [5] is applied. The computational domain is a step of 1 cm in height on a rectangle of 5 cm in height and 24 cm in length (Fig. 1). The external boundaries are far enough from the rock/cutter interaction zone to exclude their significant effect. The step's base is fixed rigidly. The surface of the step's riser moved with the constant velocity u0 = 150cm/sec.
Such linear velocity is typical for the cutter at a distance from a center of a bit while drilling at 120 rpm.
Fig. 1. Stair step sample of rock with moving riser face
Figure 2 presents the value of the destruction y for a =0.04446, p =-0.5, £ =0.125, £ = 0.0625. By the time moment t = 37.5• 105 sec the field of y takes the characteristic periodic structure with a manifested failure scale. The destruction parameter reaches its maximum value at the first (closest to the cutter's tip) peak.
Figure 2a shows the destruction peak is formed near the cutter's tip as well as attenuating destruction waves. Fig. 2c shows possible cleavages formed in the direction of cutting, and Fig. 2d shows the loose layer forming in the depths. Figures 3 and 4 show the sensitivity domain of the solutions of interest relative to the "inelastic" modulus a. Set p = -1/3, £ = 0.125 and £ = 0.0625. The calculations were performed for four values of the modulus a: ( 0.15067, 0.06841, 0.04 and 0.03). Corresponding "weakened" modulus [5] X = X-p2 /a was : ( 0.9A, 0.8i, 0.62i and 0.5i).
Fig. 2. Three-dimensional surface y(x,y) at time moment t= 37.5 10-5 sec
Figure 3 demonstrates color charts of the volume plastic deformations at
moment t = 75 • 10-5 sec. Modulus a has a direct effect on the characteristic scale of the failure zone. There is spatial periodicity in the zone. Its period increases as a decreases (X "weakening"). The few first failure zones demonstrate there is a characteristic scale for the following rock defragmentation (chip extrusion).
Chip size scale. Hnri7ontal Peak count
-Ii
2 — n AT 2 f—
sw&UL
Fig. 3. Dependence of the volume plastic deformations on modulus a
_ _ 1/2
Figure 4 shows color charts of the shear plastic deformations (yik yik) at moment t = 75 • 10-5 sec. Note that the rock fracturing process can be considered as an accumulation of shear plastic deformations that reach their critical level followed by rock fragmentation (i.e., instant strain energy release and appearance of new surfaces). The possible fragmentation lines are marked in dash lines.
Possible Rock Fragmentation
Fig. 4. Dependence of the shear plastic deformations on modulus a
The four subfigures of Figures 3 and 4 demonstrate chip size dependence on modulus a .
CONCLUSION
In this paper, the nonlocal dynamic model of inelastic deformations was applied for solving dynamic problem of failure zone formation while rock cutting, and an example of non-stationary process lead to the establishment of disintegration field causes the rock failure was provided. To make it possible, the thermodynamically consistent dynamic theory of zonal disintegration was exploited. The numerical investigation of the model 2-dimensional rock cutting problem confirmed that: a zone of plastic deformation of a rock has a characteristic scale of rock failure and the non-classic a modulus is responsible for the failure scale caused by destruction.
REFERENCES
1. Cherepanov G.P. Theory of rock cutting // Strength of Materials. - 1986. - Vol. 18(8). -P.1103-1114.
2. Borisov K.I. A Scientific Approach to Evaluating Effectiveness of Dynamic Processes of Rock Failure While Drilling Wellbores Using Modern Shearing Tools: Diss. ... Dr of Science in Engineering. - Tomsk. - 2012.
3. Dorovsky V.N., Romensky E.I., Sinev A.V. Spatially non-local model of inelastic deformations: applications for rock failure problem // Geophysical Prospecting. - 2015. - Vol. 63(4). -P.1198-1212.
4. Guzev M.A., Paroshin A.A. NonEuclidean Model of the Zonal Disintegration of Rocks around an Underground Working // Journal of Applied Mechanics and Technical Physics. - 2001. -Vol. 42(1). - P. 131-139.
5. Vtorushin E.V. Application of mixed finite elements to spatially non-local model of inelastic deformations // International Journal on Geomathematics. - 2016. - Vol. 7(2). - P. 183-201.
© E. B. BmopywuH, B. H. ffopoecKUU, 2017
