Application of Machine Learning Algorithms to Predict the Effectiveness of Radial Jet Drilling Technology in Various Geological Conditions

: This study presents a methodological approach to forecasting the efﬁciency of radial drilling technology under various geological and physical conditions. The approach is based upon the integration of mathematical statistical methods and building machine learning models to forecast the liquid production rate increment, as well as to forecast technological indexes using a hydrodynamic model. This paper reviewed the global practice of radial drilling and well intervention efﬁciency modeling. The efﬁciency of the technology in question was analyzed on the oil deposits of the Perm Territory. Mathematical statistical methods were used to determine the geological and technological parameters of the efﬁcient technology use. Based on the determined parameters, machine learning models were built, allowing us to forecast the oil and liquid production rate. A script was developed to integrate machine learning methods into a hydrodynamic simulator. When the method was tested, the deviations in the difference between the actual and the forecast cumulative oil production did not exceed 10%, which proves the reliability of the method. At the same time, the hydrodynamic model allows for taking into account the mutual inﬂuence of oil wells, the dynamics of water cut, and reservoir pressure.


Introduction
Many oil fields around the world face the growth of residual hard-to-extract oil deposits. Conditions such as the reduction of oil recovery efficiency, the growth of zones not covered with the water-oil displacement process, and a complex geological structure require specific geological and technical operations. However, the operations are not equally efficient and sometimes are not economically feasible. The search for technologies for the increasing of oil production efficiency, involving distant oil-saturated areas into production, reducing costs of hard-to-extract oil recovery, and increasing the current technology efficiency is one of the main priorities of the oil industry. The affordable and efficient technology of radial jet drilling (RJD) has been used globally, allowing for the recovery of non-drainable reserves by making highly-permeable channels up to 100 m long [1][2][3].
The analysis of the technology use reveals the variability of its efficiency under different geological and physical conditions. There is still no reliable tool for assessing the technology efficiency depending on the geological, physical, and technological properties of the oil site.
From the literature, it was noted that the radial drilling technology is a global practice, but the actual data that the previous studies are based upon appear insufficient, therefore, the influence of the geological and physical properties of the reservoir on the efficiency of the technology has not been properly studied. The accumulated experience of using the technology at the Perm Territory oil fields (Russia) allows for the use of statistical methods and machine learning methods for the analysis.

Materials and Methods
During the study, an analysis of the efficiency of radial formation drilling technology at the Perm Territory oil fields was carried out. It is important to note that the effect of radial jet drilling is different in various geological conditions. To identify the dominant properties of the reservoir and technological parameters, the methods of mathematical statistics were used-t-test, Pearson's criterion, and linear discriminant analysis. On the basis of the dominant parameters, models of neural networks and support vector machines are constructed for regression tasks. A program has been developed that allows for the combination of machine learning models and the reservoir simulator (Tempest More). A retrospective analysis was carried out to validate the new approach.

Analysis of the Efficiency of Radial Jet Drilling Technology at the Perm Territory Oil Fields
At the oil fields of the Perm Territory, the RJD technology has been used since 2006. By now, to date, 648 RJD operations have been completed ( Figure 1) [9].
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 15 In the articles [16][17][18], statistical methods were applied to compare RJD, acid treatment, and re-perforation in the Perm Territory. Radial jet drilling was more effective for the oil production rate increment.
From the literature, it was noted that the radial drilling technology is a global practice, but the actual data that the previous studies are based upon appear insufficient, therefore, the influence of the geological and physical properties of the reservoir on the efficiency of the technology has not been properly studied. The accumulated experience of using the technology at the Perm Territory oil fields (Russia) allows for the use of statistical methods and machine learning methods for the analysis.

Materials and Methods
During the study, an analysis of the efficiency of radial formation drilling technology at the Perm Territory oil fields was carried out. It is important to note that the effect of radial jet drilling is different in various geological conditions. To identify the dominant properties of the reservoir and technological parameters, the methods of mathematical statistics were used-t-test, Pearson's criterion, and linear discriminant analysis. On the basis of the dominant parameters, models of neural networks and support vector machines are constructed for regression tasks. A program has been developed that allows for the combination of machine learning models and the reservoir simulator (Tempest More). A retrospective analysis was carried out to validate the new approach.

Analysis of the Efficiency of Radial Jet Drilling Technology at the Perm Territory Oil Fields
At the oil fields of the Perm Territory, the RJD technology has been used since 2006. By now, to date, 648 RJD operations have been completed (Figure 1) [9]. The number of operations carried out on the carbonate reservoir (594) was much greater than that on the terrigenous one (54). The efficiency of the operations carried out on the carbonate reservoir was much higher than that on the terrigenous one. The average The number of operations carried out on the carbonate reservoir (594) was much greater than that on the terrigenous one (54). The efficiency of the operations carried out on the carbonate reservoir was much higher than that on the terrigenous one. The average incremental oil production was 33,291 barrels per well at the carbonate sites, while at the terrigenous sites it was 25,880 barrels. For one-third of all the operations, the effect of RJD at the terrigenous sites lasts for less than one year, and in 15% of cases the oil production increment was less than 14.56 barrels. At the carbonate sites, the number of well interventions with the effect lasting for less than a year was 7%, while the number of failed interventions was only 1%.
The main reason for the high efficiency of the RJD technology at carbonate reservoirs is the additional hydrochloric acid treatment (HAT) which is much more efficient for carbonate rocks than for terrigenous ones. There are several possible reasons for the low efficiency of RJD at terrigenous sites, such as clay swelling, canal colmation, the Jamin effect, and low residual oil reserves at the moment of RJD. The technological characteristics of RJD and the variability of its efficiency create the need for detailed special research of the conditions for using the RJD technology on carbonate and terrigenous sites [1]. Due to the rare use and low efficiency of RJD at terrigenous reservoirs, this study will only elaborate on its use on carbonate reservoirs.
At carbonate sites, 590 RJD operations were carried out. Figure 2 shows the distribution of sites by the cumulative incremental oil production due to RJD. terrigenous sites it was 25,880 barrels. For one-third of all the operations at the terrigenous sites lasts for less than one year, and in 15% of cases t increment was less than 14.56 barrels. At the carbonate sites, the number tions with the effect lasting for less than a year was 7%, while the numb ventions was only 1%.
The main reason for the high efficiency of the RJD technology at car is the additional hydrochloric acid treatment (HAT) which is much mor bonate rocks than for terrigenous ones. There are several possible reaso ficiency of RJD at terrigenous sites, such as clay swelling, canal colmation and low residual oil reserves at the moment of RJD. The technological RJD and the variability of its efficiency create the need for detailed spec conditions for using the RJD technology on carbonate and terrigenous si rare use and low efficiency of RJD at terrigenous reservoirs, this study w on its use on carbonate reservoirs.
At carbonate sites, 590 RJD operations were carried out. Figure 2 sh tion of sites by the cumulative incremental oil production due to RJD. The greatest amount of the incremental production was registered (42%-8386 thousand barrels); a comparable amount was collected fr (33%-6610 thousand barrels); they are followed by the V3V4 sites (25% barrels) and an insignificant contribution was made by the KV1 sites (0. barrels).
At the next stage, the results of well flow tests (WFTs) were analyzed not carried out for all the wells, the study only covers the data for 417 after RJD. Figure 3 shows the success rate, i.e., the percentage of positi hydrodynamic properties of the reservoir.
The success rate of the parameter change varies from 70 to 92% (on all the parameters except for skin factor (success rate 67 to 71%). The greatest amount of the incremental production was registered at the T-Fm sites (42%-8386 thousand barrels); a comparable amount was collected from the Bsh sites (33%-6610 thousand barrels); they are followed by the V3V4 sites (25%-5059 thousand barrels) and an insignificant contribution was made by the KV1 sites (0.5%-87 thousand barrels).
At the next stage, the results of well flow tests (WFTs) were analyzed. The WFTs were not carried out for all the wells, the study only covers the data for 417 wells before and after RJD. Figure 3 shows the success rate, i.e., the percentage of positive changes in the hydrodynamic properties of the reservoir.
The success rate of the parameter change varies from 70 to 92% (on average 80%) for all the parameters except for skin factor (success rate 67 to 71%).
In the last ten years, the RJD technology became one of the major technologies at the Perm Territory oil fields for wells with a decreasing productivity and a low production rate at the end of their service life. The technology efficiency varies depending on the geological and technological conditions, which requires a thorough approach to selecting candidate wells, in order to increase the general technological efficiency of the method. Appl. Sci. 2021, 11, x FOR PEER REVIEW 5 of 15 In the last ten years, the RJD technology became one of the major technologies at the Perm Territory oil fields for wells with a decreasing productivity and a low production rate at the end of their service life. The technology efficiency varies depending on the geological and technological conditions, which requires a thorough approach to selecting candidate wells, in order to increase the general technological efficiency of the method.

Analysis of the Current RJD Efficiency Forecasting Methods
This paper analyzed the main existing methods used to forecast the effectiveness of well interventions: geological and industrial analysis, statistical analysis, machine learning methods, reservoir simulation (RS) [4,19,20].
In the geological and industrial analysis, the results are limited to the assessment of specific productivity coefficients, while the totality of geological and technological parameters are not taken into account. The detailed "manual" well assessment, based on geological and industrial analysis using analytical and statistical methods, is quite subjective and time-consuming [1].
The statistical analysis requires finding and analyzing "emissions", which requires more time for data preparation. For this reason, in the last several years, forecasting based on machine learning methods has been developed, such as neural networks, decision trees, random forest algorithms, and cluster analysis. The main drawbacks of these methods are the absence of clear forecasting algorithms and physical proof, as well as low interpretability of final results.
When using RS, it is important to take into account the subjectivity of the well-intervention settings and adjustments, which can significantly affect the forecasting capacity of the model. Due to the time consumption and high cost of RS operations, they have to be used mainly for the development of high-cost well interventions (sidetracks, horizontal shafts) [19,21].
The leading research organizations are working on the development of methods for forecasting the oil and liquid production increment due to well intervention, using hydrodynamic and mathematical modeling technologies, machine learning methods, etc. The review of the previous studies revealed the need for a comprehensive assessment of geological and technological indexes and for increasing the reliability of well intervention efficiency assessment [2].

Analysis of the Current RJD Efficiency Forecasting Methods
This paper analyzed the main existing methods used to forecast the effectiveness of well interventions: geological and industrial analysis, statistical analysis, machine learning methods, reservoir simulation (RS) [4,19,20].
In the geological and industrial analysis, the results are limited to the assessment of specific productivity coefficients, while the totality of geological and technological parameters are not taken into account. The detailed "manual" well assessment, based on geological and industrial analysis using analytical and statistical methods, is quite subjective and time-consuming [1].
The statistical analysis requires finding and analyzing "emissions", which requires more time for data preparation. For this reason, in the last several years, forecasting based on machine learning methods has been developed, such as neural networks, decision trees, random forest algorithms, and cluster analysis. The main drawbacks of these methods are the absence of clear forecasting algorithms and physical proof, as well as low interpretability of final results.
When using RS, it is important to take into account the subjectivity of the wellintervention settings and adjustments, which can significantly affect the forecasting capacity of the model. Due to the time consumption and high cost of RS operations, they have to be used mainly for the development of high-cost well interventions (sidetracks, horizontal shafts) [19,21].
The leading research organizations are working on the development of methods for forecasting the oil and liquid production increment due to well intervention, using hydrodynamic and mathematical modeling technologies, machine learning methods, etc. The review of the previous studies revealed the need for a comprehensive assessment of geological and technological indexes and for increasing the reliability of well intervention efficiency assessment [2].

Brief Description of the Reservoirs
The objects of this study are the reservoirs of the Famenian, Turnesian, and Bashkirian age in the south-eastern territory of the Volga-Ural Province. The Famenian deposit rocks are accumulations of porous and cavernous reservoirs of the Solikamsk depression reef structures. The reservoir qualities of the Famenian deposits are mostly influenced by rock fissures. The Turnesian deposit reservoirs feature uneven geological sections and lesser thickness, therefore, lower well productivity [20,22,23]. The Bashkir reservoirs belong to the porous type; the porosity is created by intra-form and inter-form cavities. The best Appl. Sci. 2021, 11, 4487 6 of 14 reservoir properties are found in biomorphic limestones with a foraminiferal structure; in certain interlayers, they acquire additional density due to secondary calcitation [22,23].

Assessing the Influence of the Geological and Physical Properties of the Reservoir on the RJD Efficiency
At the present moment, it is not clear which geological conditions provide the highest efficiency of the RJD technology. An important task is to determine the influence of geological and geophysical parameters on the RJD technology efficiency. This will allow for the selection of candidate wells based on scientific data. In order to determine the most significant factors, an integrated database was made for all the wells where RJD was performed and for the relevant parameters. The database included the parameters of the geological and physical reservoir properties (GPRP) used for calculating oil deposits; the results of HDT carried out before RJD; the results of the geophysical well survey interpretation (GWSIR). In total, 31 parameters for 590 wells in 40 oil fields were included [23,24]. Table 1 shows all parameters. Oil density in the formation conditions ian age in the south-eastern territory of the Volga-Ural Province. The Famenian deposit rocks are accumulations of porous and cavernous reservoirs of the Solikamsk depression reef structures. The reservoir qualities of the Famenian deposits are mostly influenced by rock fissures. The Turnesian deposit reservoirs feature uneven geological sections and lesser thickness, therefore, lower well productivity [20,22,23]. The Bashkir reservoirs belong to the porous type; the porosity is created by intra-form and inter-form cavities. The best reservoir properties are found in biomorphic limestones with a foraminiferal structure; in certain interlayers, they acquire additional density due to secondary calcitation [22,23].

Assessing the Influence of the Geological and Physical Properties of the Reservoir on the RJD Efficiency
At the present moment, it is not clear which geological conditions provide the highest efficiency of the RJD technology. An important task is to determine the influence of geological and geophysical parameters on the RJD technology efficiency. This will allow for the selection of candidate wells based on scientific data. In order to determine the most significant factors, an integrated database was made for all the wells where RJD was performed and for the relevant parameters. The database included the parameters of the geological and physical reservoir properties (GPRP) used for calculating oil deposits; the results of HDT carried out before RJD; the results of the geophysical well survey interpretation (GWSIR). In total, 31 parameters for 590 wells in 40 oil fields were included [23,24]. Table 1 shows all parameters. Maximum oil production rate after RJD barrels/day Effect duration days Average daily oil production rate increment barrels/day Incremental oil production total barrels Table 2 demonstrates the results of t-test calculations for all wells. The values with a p-value under 0.05 (with the likelihood of differences in the compared data of over 95%) are highlighted in bold. It means that the studied parameter has a statistically relevant (not random) effect on the differences between the index values in the selected data sets. The values with a p-value ranging from 0.05 to 0.10 (also has some effect) are highlighted in italics. N1 and N2 mark the data selections for the first and the second group, respectively. R is the correlation coefficient between geological parameters and efficiency parameters. R was not high, but many parameters were significant at the p-level < 0.05, and for evaluating the influence of geological parameters on the efficiency parameters, in this case, we proposed the calculation of t-test between groups with different efficiencies. All the wells and the respective parameters were divided into two approximately even data sets selected by one of the efficiency indexes. Statistically irrelevant results were not included.
The t-test values were also calculated separately for different sites. During the research into the influence that these parameters have on the efficiency indexes, the following conclusions were made using the Student's t-test [22,25,26]:

1.
For the Fm site, the main influence on the efficiency is from the conductivity and permeability of the reservoir farfield, piezoconductivity, oil saturation, and oil-saturation thickness. At higher values of the said parameters, better efficiency of RJD should be expected. Moreover, if the bottom hole zone (BHZ) is deteriorated (positive skin factor value), the well intervention is more effective due to the use of the well potential during the BHZ cleaning.

2.
In the wells operating the T reservoir, the greatest efficiency is achieved at the higher viscosity and density of the oil and higher energy potential of the reservoir (reservoir and bottom-hole pressure). Poor BHZ condition improves the efficiency of the technology.

3.
For the Bsh site, the success of RJD was mostly due to the higher values of interval quantity, sand content and net pay reservoir thickness, the relatively contaminated BHZ (skin factor, BHZ permeability), the higher values of bottom-hole pressure, saturation pressure, piezoconductivity, BHZ, and farfield hydroconductivity, but the effect duration was longer with poorer oil viscosity and, consequently, poorer hydroconductivity.

4.
For all the sites, there is the influence of the skin factor value, oil production rate before RJD, reservoir oil saturation, and reservoir energy (bottom-hole and reservoir pressures).

Development of the Integrated RJD Efficiency Forecasting Method
When doing a long-term forecast of the technological indexes, one should consider the mutual influence of the wells, the water intrusion rate, the reservoir pressure changes as well as the geological and technological operations taking place on other wells. The integrated method is designed to resolve the mentioned problems.
The method includes the following stages: Stage 1. Forecasting the liquid production rate increment after RJD. The technological efficiency of the RJD technology may be assessed with a number of methods: statistics, machine learning, or hydrodynamic modeling. However, only a hydrodynamic simulator allows for long-term forecasting of the technological indexes taking into account the mutual influence of the wells. There are two ways of forecasting with a hydrodynamic simulator: by setting the liquid withdrawal or the bottom hole pressure limits. In the authors' practice, the method of setting the liquid rates is more common. For such estimation, it is necessary to know the changes of well operation mode after RJD to set the liquid production rate after RJD for the forecast period.
As a rule, it requires a geological analysis and calculation of the specific productivity factors that determine the production rate change after RJD [26]. The authors suggest forecasting the liquid production rate increment using machine learning methods and including these models into the hydrodynamic simulator. To develop the machine learning models, it is necessary to reduce the number of variables, in other words, to select the most relevant parameters that influence the liquid production increment. To reduce the scale, linear discriminant analysis was used for finding the set of the parameters that influence particular values.
For the Bsh site, the linear discriminant function was as follows (1): The result of the linear discriminant function classification was mostly influenced by, in descending order: oil-saturated thickness, net-to-gross sand ratio, oil viscosity, farfield hydroconductivity, water cut, and skin factor.
When the linear discriminant analysis was used on the Bsh site, 31 of 33 (94%) objects with the relative liquid production increment of less than 5 times were correctly recognized, and 30 of 35 (86%) objects with the increment of more than 5 times were correctly recognized.
The reliability of the models is due to the high training dataset classification correctness rate (83 to 97%) and the high values of the canonical correlation coefficients (0.71-0.79).
For the T site, the linear discriminant function was as follows (2): For the T site, 28 of 33 (85%) objects with the relative liquid production increment of less than 4 times were correctly recognized and 31 of 32 (97%) objects with an increment rate of more than 4 times were correctly recognized.
For the Fm site, the linear discriminant function was as follows (3): For the Fm site, 10 of 12 (83%) objects with the relative liquid production increment of less than 6 times were correctly recognized and 11 of 12 (92%) objects with an increment rate more than 6 times were correctly recognized.
Based on the revealed parameters, the machine learning models were developed to allow for the forecasting of the liquid production rate after radial jet drilling.
Initially, neural networks of different architecture were built, which quite reliably forecasted the liquid production rate after RJD (see Table 3).  Figure 4 shows the training results. The second selected forecasting method was the Support Vector Machine method. However, with the modification of the parameters influencing the learning productivity, the best productivity value achieved was 0.79 for the training selection and 0.78 for the test selection with parameters-kernel type = radial basis function, gamma = 0.16, C = 5, epsilon = 0.2. Table 4 shows the models and their corresponding statistical parameters. Where MAPE-mean absolute percentage error; STD-standard deviation; SEM-standard error mean, RMSE-root mean square error.
For further forecast, a multi-layer perceptron (MLP) was used, a neural network with a simpler architecture: 17 input layer neurons, 1 hidden layer with 5 neurons, and 1 neuron on the output layer, the neuron activation function is logistic, the error function is the sum of squares. The correlation coefficients achieved for this network were suitable for well stimulation forecasting tasks. Figure 5 shows the difference between the standard statistical method for forecasting well stimulation events efficiency in Perm region and the MLP.
A significant improvement in the quality of the forecast can be noted, which makes it possible to determine the potential of wells for RJD. Stage 2. Integrating the machine learning models with the hydrodynamic simulator. The combination of statistical and mathematical forecasting methods significantly increases the reliability of forecasting the effects of well interventions. This approach allows for the taking into account of the geological and technological parameters entered in the machine learning model, as well as the dynamic parameters from the hydrodynamic model at the moment of forecasting. The second selected forecasting method was the Support Vector Machine method. However, with the modification of the parameters influencing the learning productivity, the best productivity value achieved was 0.79 for the training selection and 0.78 for the test selection with parameters-kernel type = radial basis function, gamma = 0.16, C = 5, epsilon = 0.2. Table 4 shows the models and their corresponding statistical parameters. Where MAPE-mean absolute percentage error; STD-standard deviation; SEM-standard error mean, RMSE-root mean square error. For further forecast, a multi-layer perceptron (MLP) was used, a neural network with a simpler architecture: 17 input layer neurons, 1 hidden layer with 5 neurons, and 1 neuron on the output layer, the neuron activation function is logistic, the error function is the sum of squares. The correlation coefficients achieved for this network were suitable for well stimulation forecasting tasks. Figure 5 shows the difference between the standard statistical method for forecasting well stimulation events efficiency in Perm region and the MLP.   A significant improvement in the quality of the forecast can be noted, which makes it possible to determine the potential of wells for RJD. Stage 2. Integrating the machine learning models with the hydrodynamic simulator.
The combination of statistical and mathematical forecasting methods significantly increases the reliability of forecasting the effects of well interventions. This approach allows for the taking into account of the geological and technological parameters entered in the machine learning model, as well as the dynamic parameters from the hydrodynamic model at the moment of forecasting. The hydrodynamic simulator selected for the radial jet drilling modeling was Tempest More 8.6. The simulator is a universal hydrodynamic modeling pack used for oil, oil and gas, and gas condensate deposits.
The main parameter for the radial drilling efficiency assessment in the hydrodynamic modeling was the change of liquid production rate after RJD; skin factor and connectivity values are set by default (0, 1, respectively). Moreover, the line, length, and diameter of the radial canals are also entered into the simulator.
In order to integrate the hydrodynamic simulator and the machine learning models, a script in Python language was developed to collect the statistical parameters from the model (oil-saturated thickness, porosity, absolute permeability, number of intervals, etc.), as well as the dynamic parameters at the moment of radial drilling within the forecast period (reservoir pressure, current liquid production rate, water cut). Then, the neural network was used to calculate the liquid production rate after radial drilling.
This allows getting an automatic forecast of the liquid production increment due to drilling at any moment, which was used to assess the technological efficiency of drilling in both the long and short term.
Stage 3. Assessment of the oil production increment, water cut dynamics, reservoir pressure changes, as well as the incremental oil production due to drilling.
At the next stage, the hydrodynamic simulator was used to perform calculations taking into account the liquid production rate values calculated for the date of the planned RJD operation. Then, the technological development indexes were forecast, taking into account the potential liquid production increment after radial drilling, the dynamics of reservoir pressure and water cut, and the mutual influence of the wells.
The method makes it possible to predict the total production increment, oil production increment by years, the effect duration, changes to the water cut, and reservoir pressure dynamics, taking into account the mutual influence of all the wells and the cumulative effect of the geological and technological parameters.

Results
The developed integrated RJD efficiency forecasting method was applied to 15 wells located in the oil fields of the Perm Territory. Figure 5 shows a retrospective forecast of the oil production rate increment compared to the actual values for a number of wells in different oil fields. The oil production rate was compared as the basic parameter indicating the efficiency of well intervention. Figure 6 shows that the author's method was more accurate and closer to the actual effect of the operation, and the deviations from the incremental oil production rate did not exceed 11% (Table 5).
The developed method is universal and can be extrapolated to other improved oil recovery methods. different oil fields. The oil production rate was compared as the basic parameter indicating the efficiency of well intervention. Figure 6 shows that the author's method was more accurate and closer to the actual effect of the operation, and the deviations from the incremental oil production rate did not exceed 11% ( Table 5).
The developed method is universal and can be extrapolated to other improved oil recovery methods. Figure 6. Comparison of the forecast using the developed method to the actual average daily oil production increment data by years. Table 5. Statistical measures of the proposed forecasting method. Figure 6. Comparison of the forecast using the developed method to the actual average daily oil production increment data by years.

Discussion
The paper reviewed the global practice of radial drilling and well intervention efficiency modeling. It is noted that the RJD technology is a global practice, and its success depends on the geological and physical parameters of the reservoir.
The efficiency of the technology was analyzed on the oilfields of the Perm Territory, where a lot of experience in RJD use had been accumulated. It was revealed that the technology is much more efficient in carbonate reservoirs compared to terrigenous ones.
Mathematical statistics methods were used to determine the geological and technological parameters of the efficient technology use. For the sites of different geological age, it was revealed that various geological and technological parameters influence the efficiency indexes. For all the sites, there is the influence of the skin factor value, oil production rate before RJD, reservoir oil saturation, and reservoir energy (bottom-hole and reservoir pressures).
Based on the determined parameters, the machine learning models were built, allowing for the forecasting of the liquid production rate. A script was developed for integrating the hydrodynamic model with the machine learning models, in order to take into account a greater volume of statistical information along with the dynamic oil field parameters when forecasting.
When the method was tested, the deviation differences between the history and the forecast oil production total did not exceed 11%, which proves the reliability of the method. At the same time, the hydrodynamic model takes into account the mutual influence of oil wells, the dynamics of water cut, and reservoir pressure.

Conclusions
During the present study: 1.
The efficiency of radial jet drilling in Russian and foreign oil fields was analyzed.

2.
The industrial survey data on radial jet drilling in different wells were analyzed and systematized. For different carbonate sites of the Perm Territory oil fields, the geological, physical, and technological parameters, influencing the effectiveness of radial jet drilling, were determined. 3.
Machine learning models have been developed to predict post-RJD fluid rates for candidate wells. It can be noted that the quality of the forecast has improved according to the standard statistical method used in the fields of the Perm Territory. The Rsquared obtained for MLP is suitable for well stimulation prediction tasks.

4.
The integrated method of forecasting the incremental production rate due to RJD was developed and tested. It should be noted that with an accurate forecast of the liquid production rate after RJD, the forecast of increased oil production total on the reservoir simulation model turns out to be very accurate (relative error less than 11%).
The method makes it possible to predict the total production increment, oil production increment by years, the effect duration, changes of the water cut, and reservoir pressure dynamics, taking into account the mutual influence of all the wells and the cumulative effect of the geological and technological parameters. The method allows for both short-time and long-time forecasting of the technology's efficiency.