Analysis of Dynamical Heat Conductivity of the Reservoir and Fluid Evacuation Zone on the Gas Condensate Well Flow Rate

: This study shows that the thermal conductivity of the rock borehole adjacent to the wells varies depending on the operation of the well. This is due to the fact that the actual temperature and temperature di ﬀ erence a ﬀ ect the humidity and other thermal properties of the rocks, which in turn a ﬀ ect the heat transfer coe ﬃ cient across the section between the moving gas and the rocks. The static temperature ﬁeld of primitive geothermal gradients acquires changes in a dynamic form. Theoretical consideration of changes in the thermal conductivity of rocks near the face and the wells is proposed to improve the prediction of gas condensate wells production. The result is achieved by introducing the speciﬁed equations of the thermal energy balance in the radial ﬁltration and lifting of well products, which contain the coe ﬃ cients of heat exchange and throttling. The reﬁnement bias estimation of the 10%–15% level of gas condensate well extraction is shown using proposed methodological approach to relatively well-known (traditional in the ﬁeld development practice) methods for estimating the extraction of a “medium well” from a particular oil and gas ﬁeld evaluation. The results of this work demonstrate important scientiﬁc, applied, educational and methodological signiﬁcance of using the methodology presented by the authors. and K.F.A.; methodology, A.-S.M.B.; software, H.M.A.; validation, M.I.F. and K.F.A.; formal analysis, A.-S.M.B.; investigation, Y.A.S.; resources, Y.A.S.; data curation, K.F.A.; writing—original draft preparation, M.I.F., A.-S.M.B.; visualization, M.K.A.; supervision, M.I.F.; project administration, M.I.F. All


Introduction
Gas and condensate extraction from gas condensate wells is carried out after drafting development projects and preliminary estimation of extractive capacity. During the design of the field development, basic and advanced methods for forecasting the main development parameters are used. After commissioning, long-term operation and repair of wells, the estimation of possible and optimal operating modes is also required. This makes it necessary to continuously improve the methods of calculating thermal parameters and debit. In most cases, the forecasting of gas condensate wells production is influence by the thermal conductivity of the adjacent rocks. This may vary in the same time with, or after, thermal intensification of the hydrocarbons influx. From theoretical perspective, the connection of extraction and thermo physical parameters of rocks, the actual temperature distribution curving along the well, the heat exchange between a lot of extracted fluids and the reservoir rock, are also important to take in account the relevance of the problem in order to get a better gas condensate production estimate. Therefore, it is important to formalize the corresponding refinement problem on the extraction of gas condensate wells with essentially non-isothermal operating conditions and to develop new scientific approaches. In this paper we consider a simplified model of gas condensate well, working from one productive reservoir and has one pipe of tubing productive line, with a geothermal gradient remaining unchanged for the remote zones, and the heat exchange passes in the plane of the reservoir and near the trunk of the well. At the same time, in the thermo hydraulic calculations, it is proposed to take into account the drainage, heat dissipation and heat inflow of the "fluid-adjoining rocks", depending on the actual pressures and temperatures in the dynamics at different depths.

Research Review and Setting the Research Task
The equation of the inflow of the gas condensate mixture from the production formation to the wellbore is recorded using the filtration coefficients A, B and constant C, which generally determine the volume of natural gas extraction for a certain period [1]: where: A = z·µ·P at ·T rav In Equation (1), the system of units SI and the nomenclature of parameters are shown below in the table (except for A, B, C and k). Table 1 shows all the parameters of the study and their units. A comparison in the results of the heat flux values, using thermal conductivity [2], shows that the failure takes into account the influence of reservoir pressures and temperatures. The thermal conductivity of sedimentary rocks leads to an overestimation or undervaluation of the heat flow by 13%-30%.
During the filtration movement of the hydrocarbon mixture through the formation, there is a heat-mass-exchange process and a throttling effect, which can be simplified by the equation [3], T pl − T bh = D j · P pl − P bh + ∆T he . ( The constraint z, the isobaric heat capacity C p , the dynamic viscosity µ, and the density of natural hydrocarbons , as well as the Joule-Thomson D j coefficient from the formation to the face, will vary considerably depending on the actual pressure P and temperature T [4]. They can be determined for natural gas containing methane more than 90%, with a molar mass of M-const and a density at standard conditions ρ st -const based on the following empirical functional dependences of Latonov Gurevich, Starling Ellington, and others illustrated in [5,6]: (3) µ(P, T, M) = (9.41 + 0.02·M)(1.8·T) 1.5 (209 + 19·M + 1.8·T)·10 7 ·EXP 3.5 + 547.8 where: T pc (ρ st ) = 88.25·(0.9915 + 1.759·ρ st ); P pc (ρ st ) = 2.9585·(1.608 − 0.05994·ρ st )-pseudo-critical parameters of the gas-condensate mixture; and ρ(P, T, M) = P /(z(P, T, M).R.T)-the density of the mixture in working conditions.  After passing through a thick penetrating zone of the productive formation, the gas condensate mixture is evacuated through a lifting column, which also has heat exchange with the rocks adjacent to the well [7]. It can be seen from the typical schematic representation of the formation of the wellbore system in Figure 1. Two-way arrows show the heat exchange between the adjacent rocks and the moving extraction fluid. First, the fluid passes through the filtration zone in the reservoir (shown with the arrowheads in one direction), and then rises up through the tube or other product line of the well. After passing through a thick penetrating zone of the productive formation, the gas condensate mixture is evacuated through a lifting column, which also has heat exchange with the rocks adjacent to the well [7]. It can be seen from the typical schematic representation of the formation of the wellbore system in Figure 1. Two-way arrows show the heat exchange between the adjacent rocks and the moving extraction fluid. First, the fluid passes through the filtration zone in the reservoir (shown with the arrowheads in one direction), and then rises up through the tube or other product line of the well. In the elevator column, the heat exchange continues between adjoining rocks and overcoming its path to the extractive fluid. This can be described by the following non-isothermal vertical transport and longitudinal heat exchange equations [4]: • , , • • , ,  In the elevator column, the heat exchange continues between adjoining rocks and overcoming its path to the extractive fluid. This can be described by the following non-isothermal vertical transport and longitudinal heat exchange equations [4]: where:  -hydraulic resistance. The system of Equations (1)-(8) is closed and is used to calculate the gas flow rate of the well in cases of substantially non-isothermal processes in the radial inflow of the extracted fluid from the formation and removal onto the surface through the wellbore column [5]. The disadvantages mentioned in modern studies include, in particular: the assumption of the exchangeability of K t along the column; the equation to 0 of the parameter ∆T he ; the assumption of the invariance of D j from the edge of the formation to the mouth, based on change in the heat transfer efficiency at different depths.
Complete or partial elimination of these shortcomings of modern theories is an actual task. This will lead to an improvement in the accuracy of the gas condensate well simulation and the prediction of extraction in general. However, the developed technique, suitable for gas condensate wells, can be modified to a universal one as for light oil production (more in the liquid composition than gas). This leads to analyzing, further, the compressibility, the heat capacity and the viscosity of the light crude oil.

The Purpose and Task of the Research
The purpose of the work is to develop a mathematical model of non-isothermal radial inflow and lifting (evacuation on the inner cavity of the tubing) of a gaseous mixture in a gas condensate well. This takes into account the determined thermal coefficient of heat conductivity λ t (x), which also determined the geothermal gradient in depth and the temperature difference of rock-fluid in working conditions. To evaluate the thermal conductivity λ t , we used the resulting heat transfer coefficient K t , which are related by the calculated area during radial filtration in the reservoir and evacuation of the fluid-like coolant through the pipe.
To achieve the goal, the task is to fulfill comparative analysis of modeling results based on a developed theoretical improvement with modeling according to the indicated general effect

Research Methods
The research methodology consisted of a comparative analysis of two effects accounting combinations of gas condensate fluid extraction: the effect of throttling and alternating heat transfer along the way of mixing the mixture. In this case, the flow mode parameters of the extracting hydrocarbon mixture, the discharge of the wells were calculated.
To solve this problem, the authors used the mathematical method of introducing artificial coefficients of elimination effects (throttling and alternating heat transfer), which allowed us to make a comparative analysis under identical conditions by other technical parameters and indicators. Under the name coefficient of elimination, the authors proposed to understand a parameter with an absolute value from 0 to 1, which, when introduced into the formula, could eliminate (annihilate, if it is 0) or completely take into account (if it is 1) any effect.
The system of nonlinear equations of the developed mathematical model was solved by the Runge-Kutta fourth order method and the Quasi-Newtonian method.

Research Results
The equations development for pseudo critical parameters of a hydrocarbon gaseous extraction mixture was made on the assumption of a small amount (total percentage by mass) of nitrogen and carbon content. With a significant content of these gases in the extracted fluid, more complicated analytics are needed [8,9], but in this paper-investigated for the calculation example, the following was used: The calculation of the Joule-Thomson coefficient was made on the assumption that methane in the production well is more than 95%, and the relative density ∆ 0.6. After introducing the coefficient of elimination of the throttle effect K a (which is possible with increasing the humidity of natural gas), as well as the use of (5)-(6) and [6] the following working formula of the study is obtained: Equations (9)-(11) were used universally for calculations of both the inflow of the fluid to the face and the evacuation of the gas condensate mixture on the tubing of the well.
The hydraulic resistance in the pump-compressor or other lift pipe of the well λ is a function of R e , temperature and other regime and structural parameters [9]. Therefore, in the study of non-isothermal lifting on the surface of the hydrocarbon mixture, the working functional dependence λ (P, T, M, M q , D and k e ) was developed and used for the Equation (7) based on Colebrook-White equation, and the studies of S. Borisov and I. Hodanovich [10] (instead of the above functional dependence λ (k e , R e and D)): because R e P, T, M, M q , D = 4·M q µ(P,T,M)·π·D . The heat transfer coefficient K t in the basic traditional methods is not considered to be variable along the production casing of the well. So, the elimination coefficient K a led to the fixation of the value at one average level. However, in the developed and improved mathematical model, on the contrary, the explicit dependence of K t on longitudinal thermal conditions is used according to the Vlasov formula and the recommendations of [2,11] of the form: In the rocks along the depth of the axis, the temperature has a static distribution, but the temperature of the mobile gas condensate mixture in the well and the casing is characterized by a dynamic distribution of temperatures in depth. Accordingly, the thermal conductivity from the fluid to the rocks also receives a dynamic. Sometime after the well starts, the dynamic distribution of the thermal conductivity parameter will have a certain steady value, which is proposed to be taken into account in the calculations of the predicted well mass flow and thermal flow.
The temperature of the rocks at each separate depth of the well column will be different, so the Equation (8) turns into a formula: It follows from Equation (14) that losses (dissipation) or heat influx are possible at various depths. It depends on the temperature difference T o (x) − T(x).
Using the proposed additional functional dependences and the Equations (9)- (14) in the system of Equations (1)-(5), (7), (11) and (14), the authors show the influence of the actual distribution K t (T o (x)) and the value of D j (P av , T av and M) to the level of well extraction (relative to the simulation in the original system of Equations (1)-(8) of constant mean values K t = 1.5 W/(m 2 K)-const and D j = 2.5 K/MPa-const).  Figure 2 shows the mass flow rate dependence of well production and temperature at the mouth on the level of reservoir pressure in the productive deposit. It is seen that the influence of the throttling process in the actual thermobaric conditions, with the help of Equation (11), reflected more change in the debit than the parallel use of dependence Equation (13). Figure 2 shows the mass flow rate dependence of well production and temperature at the mouth on the level of reservoir pressure in the productive deposit. It is seen that the influence of the throttling process in the actual thermobaric conditions, with the help of Equation (11), reflected more change in the debit than the parallel use of dependence Equation (13).
It was also obvious that for a calculated example with the data in the table of parameters nomenclature, the effect of reducing the rocks heat conductivity with depth and the throttle effect were multi-directional on mass flow rate. For the analyzed variants of the wells regime parameters, the maximum change in well flow at average reservoir pressure reached 10%-15%. The last fact indicated led to a corresponding change in the heat flux by the coolant itself. The dynamic heat transfer coefficient depended on the throttle effect on the temperature distribution in the dynamics of the extraction along the wellbore. The dependence of Dj(P, T) is shown in Figure 3.
The rocks temperature, at different depths of the well hole, could be taken in accordance with the thermal gradients (established during drilling) only in the absence of heat exchange with the contents of the well. After a certain time with the extraction of a liquid heated by the lower horizons, a balance was established in the system of the wellbore as longitudinal dynamic pressures and longitudinal dynamic temperatures. Temperatures, in particular, will depend not only on the temperature of the fluid in the productive horizon, but on the heat transfer processes and throttling from the depths to the surface. The latter fact seemed to be trivial at the first stage of the research. However, the quantitative assessment using the refined mathematical model shows significant changes in relation to the basic methodology of the most up-to-date educational textbooks on the gas condensate field development. Thus, in industrial practice, with the precise purpose of estimating extraction, it is not appropriate to consider hydraulic parameters without using the specification and clarification of the dynamic distribution (established after the start of extraction) of rock temperatures along the column, at least with the help of the Equations (11)- (14). It was also obvious that for a calculated example with the data in the table of parameters nomenclature, the effect of reducing the rocks heat conductivity with depth and the throttle effect were multi-directional on mass flow rate. For the analyzed variants of the wells regime parameters, the maximum change in well flow at average reservoir pressure reached 10%-15%. The last fact indicated led to a corresponding change in the heat flux by the coolant itself.
The dynamic heat transfer coefficient depended on the throttle effect on the temperature distribution in the dynamics of the extraction along the wellbore. The dependence of D j (P, T) is shown in Figure 3.
The rocks temperature, at different depths of the well hole, could be taken in accordance with the thermal gradients (established during drilling) only in the absence of heat exchange with the contents of the well. After a certain time with the extraction of a liquid heated by the lower horizons, a balance was established in the system of the wellbore as longitudinal dynamic pressures and longitudinal dynamic temperatures. Temperatures, in particular, will depend not only on the temperature of the fluid in the productive horizon, but on the heat transfer processes and throttling from the depths to the surface. The latter fact seemed to be trivial at the first stage of the research. However, the quantitative assessment using the refined mathematical model shows significant changes in relation to the basic methodology of the most up-to-date educational textbooks on the gas condensate field development. Thus, in industrial practice, with the precise purpose of estimating extraction, it is not appropriate to consider hydraulic parameters without using the specification and clarification of the dynamic distribution (established after the start of extraction) of rock temperatures along the column, at least with the help of the Equations (11)- (14).  The practical value of the developed mathematical model could be especially noticeable in the case of using known thermal methods for the intensification of hydrocarbon production. The results of the first stage of such studies according to the last statement were made at the conference [12]. Figure 4 shows the dependence of the well flow rate on the formation temperature in the productive horizon for different values of the dynamic heat transfer coefficient from the fluid (inside the wellbore) to the adjacent rocks. The dynamic value of thermal conductivity was associated with the corresponding dynamic value of the heat transfer coefficient in the selected geometry of the heat engineering system (Figure 1). From Figure 4, it can be seen that the well flow significantly depended on the heat transfer coefficient. Changing the heat transfer coefficient from 1 to 2 W/(m 2 •K) led to a change in the flow rate to 5%-15%. The practical value of the developed mathematical model could be especially noticeable in the case of using known thermal methods for the intensification of hydrocarbon production. The results of the first stage of such studies according to the last statement were made at the conference [12]. Figure 4 shows the dependence of the well flow rate on the formation temperature in the productive horizon for different values of the dynamic heat transfer coefficient from the fluid (inside the wellbore) to the adjacent rocks. The dynamic value of thermal conductivity was associated with the corresponding dynamic value of the heat transfer coefficient in the selected geometry of the heat engineering system (Figure 1). From Figure 4, it can be seen that the well flow significantly depended on the heat transfer coefficient. Changing the heat transfer coefficient from 1 to 2 W/(m 2 ·K) led to a change in the flow rate to 5%-15%.  The practical value of the developed mathematical model could be especially noticeable in the case of using known thermal methods for the intensification of hydrocarbon production. The results of the first stage of such studies according to the last statement were made at the conference [12]. Figure 4 shows the dependence of the well flow rate on the formation temperature in the productive horizon for different values of the dynamic heat transfer coefficient from the fluid (inside the wellbore) to the adjacent rocks. The dynamic value of thermal conductivity was associated with the corresponding dynamic value of the heat transfer coefficient in the selected geometry of the heat engineering system (Figure 1). From Figure 4, it can be seen that the well flow significantly depended on the heat transfer coefficient. Changing the heat transfer coefficient from 1 to 2 W/(m 2 •K) led to a change in the flow rate to 5%-15%. From the use of the isothermal equation of Adamov and the binary equation of the gaseous fluid inflow to the wellbore, in thermobaric conditions methods, the intensification of borehole extraction could reach 15%-30%. This was due to the fact that the traditional adaptation of parameters A and B and the reduction of error require time and resources for wells research in different modes. That time and resources are usually lacking in real industrial production to ensure the highest quality of scientific research. Changing the reservoir temperature to 50-80 degrees led to a doubling of the heat transfer coefficient (see Figures 4 and 5), which could lead to a change in mass flow rate up to 30%. Figure 5 shows that for the smallest reservoir pressure (20 MPa), the effect on the flow rate was slightly lower, but for a temperature difference of 80 degrees, it was about 30%. Figure 5 shows the flow rate and temperature range of the fluids at the mouth of the gas condensate well with a total change in pressure in the reservoir up to 10 MPa. Changes in operating temperatures in the geothermal reservoir can also lead to large fluctuations in reservoir pressure. We had chosen the range of reservoir pressure changes from experience in operating real oil and gas fields obtained during intensification by injection of hot water and steam.
The simulation was carried out in the MathCAD 15 program, which allowed checking the results of calculations and short term forecasting in the office and field work with the help of compact "Gadgets" from Windows 10. The verification showed the suitability of the calculation method for practical use under industry conditions, which also excluded the need for using more demanding (than described above) computing resources and simulator software.
It should be noted that the solution of the system of nonlinear equations by the above-mentioned methods did not provide solutions at all intermediate points of the parameter ranges. The curves had significant gaps (Figures 2, 4 and 5), but in general tendencies and the possibility of interpolating and comparing the investigated functions. The next step in the above studies was to consider the heat exchange over the cross section and the total surface of the formation with ΔThe. That also confirmed the significant effect of this heat transfer on the prediction of a gas condensate well extraction, also provided that thermo chemical or thermo physical intensification methods were carried out. From the use of the isothermal equation of Adamov and the binary equation of the gaseous fluid inflow to the wellbore, in thermobaric conditions methods, the intensification of borehole extraction could reach 15%-30%. This was due to the fact that the traditional adaptation of parameters A and B and the reduction of error require time and resources for wells research in different modes. That time and resources are usually lacking in real industrial production to ensure the highest quality of scientific research. Changing the reservoir temperature to 50-80 degrees led to a doubling of the heat transfer coefficient (see Figures 4 and 5), which could lead to a change in mass flow rate up to 30%. Figure 5 shows that for the smallest reservoir pressure (20 MPa), the effect on the flow rate was slightly lower, but for a temperature difference of 80 degrees, it was about 30%. On the other side, light oil was less compressible than gas, so isothermal compressibility worked favorably in term of heat transfer in well extraction. It is also important to notice that the isothermal compressibility of light oil variation was depending on the impact of gravity changes and with the solution gas oil ratio (GOR). Figure 6 shows how isothermal compressibility changes with crude oil gravity. As oil gravity increases, isothermal compressibility should increase. Results predicted by De Ghetto [13,14] do not properly model the phenomena. De Ghetto proposed a method that uses several equations covering various American petroleum institute (API) gravity ranges. This technique results in discontinuities in predicted properties as the equations change.
A similar approach was used to develop a simple equation for the variation of the heat capacity or specific heat of crude oil as a function of API gravity and temperature. Data was regressed to obtain the algorithm presented as Equation (15) for foot pound second unit (FPS) units: Equation  Figure 5 shows the flow rate and temperature range of the fluids at the mouth of the gas condensate well with a total change in pressure in the reservoir up to 10 MPa. Changes in operating temperatures in the geothermal reservoir can also lead to large fluctuations in reservoir pressure. We had chosen the range of reservoir pressure changes from experience in operating real oil and gas fields obtained during intensification by injection of hot water and steam.
The simulation was carried out in the MathCAD 15 program, which allowed checking the results of calculations and short term forecasting in the office and field work with the help of compact "Gadgets" from Windows 10. The verification showed the suitability of the calculation method for practical use under industry conditions, which also excluded the need for using more demanding (than described above) computing resources and simulator software.
It should be noted that the solution of the system of nonlinear equations by the above-mentioned methods did not provide solutions at all intermediate points of the parameter ranges. The curves had significant gaps (Figures 2, 4 and 5), but in general tendencies and the possibility of interpolating and comparing the investigated functions. The next step in the above studies was to consider the heat exchange over the cross section and the total surface of the formation with ∆T he . That also confirmed the significant effect of this heat transfer on the prediction of a gas condensate well extraction, also provided that thermo chemical or thermo physical intensification methods were carried out.
On the other side, light oil was less compressible than gas, so isothermal compressibility worked favorably in term of heat transfer in well extraction. It is also important to notice that the isothermal compressibility of light oil variation was depending on the impact of gravity changes and with the solution gas oil ratio (GOR). Figure 6 shows how isothermal compressibility changes with crude oil gravity. As oil gravity increases, isothermal compressibility should increase. Results predicted by De Ghetto [13,14] do not properly model the phenomena. De Ghetto proposed a method that uses several equations covering various American petroleum institute (API) gravity ranges. This technique results in discontinuities in predicted properties as the equations change. From the practical results of the study, graphically presented in Figure 6, Equations (19) and (20), it could be determined that for various operating conditions and compositions of light oils, the heat capacity could change on average up to 1.5 times (T = 0-260 F), and compressibility on average changed up to 2.2 times. With an increase in heat capacity and a decrease in compressibility, heat transfer increased close to directly proportional, i.e., accounting for changes in the dynamic field was even more important (compared with gas-like mixtures, for accurate determination of the heat transfer intensity).
The recommended heat capacity calculation Equation (16) favorably complemented the developed Equations (10)- (14), which ultimately showed to us the best forecasting result for the dynamical heat properties of the reservoir and fluid evacuation zone.

1
Taking into account the actual change in the thermal conductivity of rocks due to depth, depending on their temperature (setting the value in dynamics), as well as the change in temperature of the extracting fluid from throttling in the well cavity led to a precision forecasting the flow of wells by 10%-15%.
A similar approach was used to develop a simple equation for the variation of the heat capacity or specific heat of crude oil as a function of API gravity and temperature. Data was regressed to obtain the algorithm presented as Equation (15) for foot pound second unit (FPS) units: Equation (16) for international system of units (SI). Note that the algorithm was developed for a crude oil with a universal oil products collection (UOP) index of 11.8 (indicating intermediate, paraffinic-naphthenic light crude oil). If the UOP index is known, the correction factor illustrated on the graph could be applied to the output from Equation (15) or (16).
The resulting Equation (14) is presented as Equation (15) and; From the practical results of the study, graphically presented in Figure 6, Equations (19) and (20), it could be determined that for various operating conditions and compositions of light oils, the heat capacity could change on average up to 1.5 times (T = 0-260 F), and compressibility on average changed up to 2.2 times. With an increase in heat capacity and a decrease in compressibility, heat transfer increased close to directly proportional, i.e., accounting for changes in the dynamic field was even more important (compared with gas-like mixtures, for accurate determination of the heat transfer intensity).
The recommended heat capacity calculation Equation (16) favorably complemented the developed Equations (10)- (14), which ultimately showed to us the best forecasting result for the dynamical heat properties of the reservoir and fluid evacuation zone.

1.
Taking into account the actual change in the thermal conductivity of rocks due to depth, depending on their temperature (setting the value in dynamics), as well as the change in temperature of the extracting fluid from throttling in the well cavity led to a precision forecasting the flow of wells by 10%-15%.

2.
Using developed mathematical models gave an opportunity to specify that 15%-30% of the debit in the conditions of implementation of thermobaric methods of intensification of the gas-condensate extraction (in relation to basic methods that do not consider the dynamic distribution of temperatures in adjacent well of rocks and the corresponding longitudinal thermal conductivity).

3.
Light oil well extraction with its parameters (compressibility, heat capacity and viscosity) required less heat exchange surface, and guarantee a more intense heat transfer during production. An increase in the percentage of hydrocarbon condensates (in a gas-liquid mixture-well production) required an even closer examination of the dynamical heat conductivity of the reservoir rock on the gas condensate flow rate. Funding: This research was partially funded by The International University of Grand Bassam, Ivory Coast (www.iugb.org).

Acknowledgments:
The authors are grateful to PJSC "Ukrgazvydobuvannya" and the UkrNDIgaz Institute for a grant approach in providing materials and data on gas condensate wells for the required period of time to the extent necessary to establish the adequacy of the developed mathematical model for comparison with real industrial research and the possibilities of adaptation in science-intensive production.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature
API American petroleum institute FPS foot pound second unit GOR gas oil ratio SI international system of units UOP universal oil products collection