Methods of Predicting the E ﬀ ectiveness of Hydrochloric Acid Treatment Using Hydrodynamic Simulation

: The article describes the original technique of predicting the e ﬀ ectiveness of hydrochloric acid treatment of the bottom hole zone of a carbonate reservoir. The technique consists in determining the technological e ﬀ ectiveness of the oil recovery stimulation procedure at productive wells using hydrodynamic simulation based on the calculated value of skin factor change. In the course of the study, a number of parameters a ﬀ ecting the intensity of ﬂow coe ﬃ cient decline in the near-wellbore zone during acid treatment have been established. The paper presents a comparison of data on the actual change of the skin factor after acidizing jobs at the ﬁelds of Perm Krai (Russian Federation) and the calculated values obtained using the presented technique. This scientiﬁc research includes an example of practical application of the proposed technique for the target well of the Kokuyskoye oil ﬁeld in Perm Krai, its results with a minor deviation coinciding with the actual values of the ﬂow rate. to the study, it is noted that using the proposed technique, recommendations can be provided on the selection of acid composition and technology of its injection into formation for the preliminary evaluation of the cost e ﬀ ectiveness of the designed procedure.


Introduction
The bottom-hole zone (BHZ) condition deteriorates during the entire well operation due to damage of flow channels with different particles. Since the contaminated area near the well reduces the reservoir-to-well hydrodynamic connectivity, the well productivity does not reach the maximum possible value, while carbonate rocks are characterized by low flow properties per se. The necessity of detailed study of properties of carbonate reservoir and methods of their development is propelled by the presence in them of one-third to half of hydrocarbon reserves worldwide (Khizhnyak et al., 2015 [1]).
Hydrochloric acid treatment (HAT) is the most effective and cost-efficient method of oil production stimulation at the deteriorated state of near-wellbore zone (NWZ) as of today (Aidagulov et [5]; and Singh et al., 2019 [6]). However, acidizing does not always produce the planned results due to insufficient consideration of key factors such as flow properties of reservoir, bottom-hole zone structure, and its parameters before well stimulation, as well as their behaviour in the process of treatment with acid composition (AC).
The success of chemical treatment in the bottom-hole zone is demonstrated by a negative change in the skin factor, which is calculated from the Hawkins' formula (Equation (1)) (Hawkins, 1956 [7]): Equation (1) takes into account the difference in permeability of the near-wellbore and far field zones (Figure 1), based on which, the degree of the bottom-hole zone contamination is evaluated.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 2 of 13 The success of chemical treatment in the bottom-hole zone is demonstrated by a negative change in the skin factor, which is calculated from the Hawkins' formula (Equation 1) (Hawkins, 1956 [7]): Equation 1 takes into account the difference in permeability of the near-wellbore and far field zones (Figure 1), based on which, the degree of the bottom-hole zone contamination is evaluated. The acid treatment process is considered to be quite complex and inexplicable in the details of the process of dissolution of minerals from carbonate reservoir rocks by hydrochloric acid. Therefore, when forecasting the technological effect of the considered method of stimulation of hydrocarbons production, there are many issues (Ali et al., 2019 [8]) arising associated with the inability to forecast the nature of changes in the void space of the productive formation near the well. Conducting flow studies on core samples under reservoir conditions is one of the ways to address these issues (Abdrazakov et al., 2019 [9]; Hall-Thompson et al., 2020 [10]). Computer tomography allows to determine the change of pore space, namely, the size and position of the channel that has broken through the sample, which is called a "wormhole." However, in the article (Gurbatova et al., 2013 [11]), the authors point to the unreliability of the results obtained during the studies on standard samples of rock formation using a flow unit in relation to the real acid treatments at the oil field, which is referred to high anisotropy of reservoir properties of carbonate oil-and gas-saturated rocks (Farooq et al., 2019 [12]).
Most of scientific research works is aimed at studying the process of "wormholing" forming in rocks under the effect of acid compositions. The authors of the work (Kanevskaya and Novikov, 2018 [13]) have analysed the existing models describing this process in detail. The application of the studied models is considered in the publication (Trushin et al., 2019 [14]), where semiempirical models of Gong, Furui, and Buijse have been used, which allow to predict the direction of AC spreading in the pore space of BHZ and to estimate the increase in permeability coefficient in this area. The obtained flow characteristics allow to calculate additional fluid inflow using available formulas of plane radial flow of fluids to hydrodynamically imperfect wells due to reduction of additional flow resistance provided by change in the skin factor.
An example of initial evaluation of the hydrochloric acid treatment effectiveness on carbonate rocks is also given in the work (Khuzin et al., 2019 [15]); its authors have developed an acidizing effect simulator that allows to calculate changes in the structure of the rock near the well. The software is The acid treatment process is considered to be quite complex and inexplicable in the details of the process of dissolution of minerals from carbonate reservoir rocks by hydrochloric acid. Therefore, when forecasting the technological effect of the considered method of stimulation of hydrocarbons production, there are many issues (Ali et al., 2019 [8]) arising associated with the inability to forecast the nature of changes in the void space of the productive formation near the well. Conducting flow studies on core samples under reservoir conditions is one of the ways to address these issues (Abdrazakov et al., 2019 [9]; Hall-Thompson et al., 2020 [10]). Computer tomography allows to determine the change of pore space, namely, the size and position of the channel that has broken through the sample, which is called a "wormhole." However, in the article (Gurbatova et al., 2013 [11]), the authors point to the unreliability of the results obtained during the studies on standard samples of rock formation using a flow unit in relation to the real acid treatments at the oil field, which is referred to high anisotropy of reservoir properties of carbonate oil-and gas-saturated rocks (Farooq et al., 2019 [12]).
Most of scientific research works is aimed at studying the process of "wormholing" forming in rocks under the effect of acid compositions. The authors of the work (Kanevskaya and Novikov, 2018 [13]) have analysed the existing models describing this process in detail. The application of the studied models is considered in the publication (Trushin et al., 2019 [14]), where semiempirical models of Gong, Furui, and Buijse have been used, which allow to predict the direction of AC spreading in the pore space of BHZ and to estimate the increase in permeability coefficient in this area. The obtained flow characteristics allow to calculate additional fluid inflow using available formulas of plane radial flow of fluids to hydrodynamically imperfect wells due to reduction of additional flow resistance provided by change in the skin factor.
An example of initial evaluation of the hydrochloric acid treatment effectiveness on carbonate rocks is also given in the work (Khuzin et al., 2019 [15]); its authors have developed an acidizing effect simulator that allows to calculate changes in the structure of the rock near the well. The software is based on the development of a "wormhole" in the process of AC treatment, as a model for which a semiempirical model proposed by Buijse and Glasbergen was used (Buijse and Glasbergen, 2005 [16]). This invention integrates all permeability and porosity properties of oil-saturated interlayers into the model and simulates their change when simulating AC injection. The success of the HAT procedure is also evaluated by reducing flow resistance coefficient in the NWZ.
Proceeding from the analysis of publications, it is possible to draw a conclusion that the existing methods of acid treatment effectiveness prediction are largely based on the determination of the skin factor change, which, in turn, changes with the increase of NWZ permeability coefficient. For a more accurate evaluation of the effectiveness of the designed geotechnical procedures, the hydrodynamic simulation is applied, which allows to calculate and evaluate the increase in well productivity. In Perm Krai, when designing the acidizing effect on oil producing wells, nowadays in most cases, the field data from previously conducted procedures of oil recovery stimulation in the same field, or in analogue fields, are included as a basis. This accounts for the discrepancy between the actual and predicted values of oil flow rate due to the wrong approach to the creation of the HAT model carried out when interpreting the degree of change in the skin factor from previous chemical treatments of reservoirs.

Materials and Methods
To predict the effectiveness of hydrochloric acid treatment of carbonate deposits this work proposes hydrodynamic simulation to determine the gain in incremental oil production and duration of the effect, based on theoretical calculations of the decrease of the skin factor. Daccord and Lenormand's experiments in the study of the interaction of AC and rock have brought to an empirical equation (Equation (2)) allowing to evaluate the change in the skin factor after the acid effect on the bottom-hole zone (Daccord et al., 1989 [17]): When calculating the fractal dimension d, the quantitative content of limestone and dolomite in carbonate rock should be taken into account. Below is an equation for calculating the value of fractal dimension (Equation (3)): In order to determine the percentage of calcite and dolomite in the rock, it was necessary to carry out laboratory studies of the rock, using carbonate metering device, such as KM-04M that allowed to determine the mass content of calcite, dolomite, and insoluble mineral residue in the crushed rock sample. Figure 2 shows the image of the carbonate metering device KM-04M used in the present research for quantitative determination of calcite, dolomite, and insoluble mineral residue in the rock under study. The acid value was determined by the ratio given in Equation 4: The concentration value of the reservoir depends on the percentage ratio of limestone and The acid value was determined by the ratio given in Equation (4): The concentration value of the reservoir depends on the percentage ratio of limestone and dolomite in the rock under study. It was determined by the following equality (Equation (5)): To determine the concentrations from reference data, the values of molecular masses of the reacting rocks, i.e., straight limestone and dolomite, and hydrochloric acid solution, were used. The concentration of the corresponding substance was calculated using Equation (6): According to GOST 3118-77 standard, the density of hydrochloric acid should be 1150-1190 g/L, therefore, a mean value of 1170 g/L was used for the concentration calculation. According to the reference data, the density of calcite and dolomite were taken to be equal to 2710 and 2880 g/L, respectively, and their molecular masses were taken to be equal to 100.086 and 184.399 g/L, respectively.
According to the reaction equations, the stoichiometric coefficient for limestone was taken to be equal to 2, and 4 for dolomite, so the reservoir stoichiometric coefficient was determined like the fractal value (Equation (7)): This work proposes to estimate the acid diffusion coefficient using the Stokes-Einstein equation (Chordia and Trivedi, 2010 [18]) that is widely used to determine the D value in fluids (Equation (8)): One of the variants for determining the dynamic viscosity of AC was the method used as per GOST 33768-2015 standard, using VPZh-4 viscometer, shown in Figure 3. The radius of the diffusing particle can be determined by means of a binocular microscope.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 13 For the sake of further convenience, Equations 2-8 were grouped and transformed into Equation 9: For the sake of further convenience, Equations (2)-(8) were grouped and transformed into Equation (9): By substituting constants in Equation (9), Equation (10) was obtained: When analysing Equation (10), we can conclude that the change in the skin factor after the acidizing effect on the carbonate reservoir depends on a number of parameters given in Equation (11): Based on Equation (11), to predict the effectiveness of HAT, the content of calcite and dolomite in the reservoir rock, as well as its current absolute temperature, rheological properties of AC, and the radius of its diffusing particles were determined. It is also necessary to define the perforated formation thickness, the well radius, and the designed amount of AC to be injected into the formation, which can vary widely based on the BHZ condition and the expected effect. After defining these parameters and calculating the skin factor change by Equation (10), the obtained value of ∆S was inputted into the hydrodynamic model of the production facility. The standard isothermal flow model of the three-phase (oil, water, and gas) three-component (oil, water, and oil-dissolved gas) fluid in pore space was used when building the field hydrodynamic model. Well perforation interval is characterized by the skin factor and "well-to-reservoir connectivity" parameters. When simulating the acid treatment, the calculated change in the skin factor was taken into account and the key technological parameters of the well operation were calculated. As a result of the calculations conducted in the hydrodynamic simulator (e.g., IRAP Tempest or Eclipse), diagrams of changes in the annual oil and water production by the basic variant and after the designed job were built. Based on the change in the annual oil and water production, a conclusion was made on the predicted effectiveness of acid treatment of the target well operating an oil-saturated reservoir with deteriorated condition of BHZ.

Research Results and Obtained Factual Findings on the Predicted Effectiveness Evaluation of Hydrochloric Acid Treatment at the Fields of Perm Krai
Hydrochloric acid treatments are used in Perm Krai to stimulate oil production from carbonate formations. The work reviews the practice of 18 acidizing jobs at the oilfields in Perm Krai. Using multiple regression analysis, the dependency of additional oil production on a number of parameters has been determined, which can be used to promptly estimate the effectiveness of hydrochloric acid treatment (Equation (12)): As shown in Figure 4, the relationship can be characterised as tight, however, this equation is of little practical use (R 2 = 0.70).
As shown in Figure 4, the relationship can be characterised as tight, however, this equation is of little practical use (R 2 = 0.70). To determine the possibility of practical application of obtained Equation 10 for predicting the effectiveness of the designed acid treatment, statistical analysis has been conducted. For each of reviewed procedure, the expected change in the skin factor is determined, calculated based on the presented technique of prediction of change in the skin factor after the HAT job. Table 1 provides a comparison of the actual change in the skin factor for the wells under study and the predicted value obtained using the given technique.
As it can be seen, the relative errors presented in Table 1 vary in the range from 1.07% to 11.60%, which is an allowable deviation. For more accurate analysis and determination of the significance of the relationship between the two sets of values, i.e., actual and calculated changes in the skin factor, the main statistical indicators shall be used. As determined by calculations using STATISTICA software, the Student's t test and the achieved significance level constitute 0.079 and 0.937, respectively, while standard deviations do not exceed allowable value range, as well, i.e., the internal scatter about the mean is insignificant. The conducted analysis indicates the high significance of the relation between actual and design data, which testifies to the effective approach to the prediction of HAT results using the proposed technique.
Comparative analysis of actual and design results has afforded the authors of this work, the assurance in the idea of using Daccord and Lenormand's formula in the basis of predicting the success of acid treatment. An example of the practical application of the proposed technique on the existing well, which consists in determining the HAT effectiveness using hydrodynamic simulation, is given below in this paper.  To determine the possibility of practical application of obtained Equation (10) for predicting the effectiveness of the designed acid treatment, statistical analysis has been conducted. For each of reviewed procedure, the expected change in the skin factor is determined, calculated based on the presented technique of prediction of change in the skin factor after the HAT job. Table 1 provides a comparison of the actual change in the skin factor for the wells under study and the predicted value obtained using the given technique. As it can be seen, the relative errors presented in Table 1 vary in the range from 1.07% to 11.60%, which is an allowable deviation. For more accurate analysis and determination of the significance of the relationship between the two sets of values, i.e., actual and calculated changes in the skin factor, the main statistical indicators shall be used. As determined by calculations using STATISTICA software, the Student's t test and the achieved significance level constitute 0.079 and 0.937, respectively, while standard deviations do not exceed allowable value range, as well, i.e., the internal scatter about the mean is insignificant. The conducted analysis indicates the high significance of the relation between actual and design data, which testifies to the effective approach to the prediction of HAT results using the proposed technique.
Comparative analysis of actual and design results has afforded the authors of this work, the assurance in the idea of using Daccord and Lenormand's formula in the basis of predicting the success of acid treatment. An example of the practical application of the proposed technique on the existing well, which consists in determining the HAT effectiveness using hydrodynamic simulation, is given below in this paper.
In July 2019, an oil-producing well of the Kokuyskoye oil field in Perm Krai, entering the Tournaisian strata, was determined as the target facility for the hydrochloric acid treatment job. To evaluate the perfection of the proposed method of predicting the acidizing effectiveness, the authors have presented a comparison of actual and design production rates for the subsequent 6 months of operation of the well under study. Additionally, using Tempest MORE hydrodynamic simulator, a geotechnical procedure was designed, and the technological effect of planned acid treatment of bottom-hole zone for the following 10 years was analysed.
FLUXOCORE-210 (grade P) acid composition was selected as an agent, as it has proved to be excellent in the field experience of oil recovery stimulation methods and has shown good results in laboratory studies presented in the paper (Ravelev, 2019 [19]). According to the results of hydrodynamic studies of the target oil producing well of 15 June 2019, the value of the skin factor was determined as equal to 1.1. The authors calculated the reduction of flow resistance in the near-wellbore zone using the above technique, which allowed to further calculate the increase in the target well flow rate after acid treatment on the hydrodynamic model.
Before acid treatment of the productive formation, in order to eliminate undesirable consequences, it was expedient to carry out laboratory research on the compatibility of AC with fluids inflowing from the target facility. The study of their interaction plays an important role in the design of hydrochloric acid treatment, since in case of poor compatibility, there is a possibility of a high-viscosity emulsion, clots, flakes, and precipitation that result in damage of void space. For this purpose, testing with 15 different ratios of oil, water, and acid composition was performed under laboratory conditions. In the beginning, stable water-oil emulsion from the selected samples of oil and water from the target well under study was obtained; after that, the acid composition under test was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.
As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids. Table 2. Results of study of compatibility of formation fluids collected from the target well, with the acid composition under test.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Water-25%
Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the Water-0% Oil-50% AC-50% Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 13 was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 13 was added to the obtained emulsion. As a result, the state of the obtained mixture was assessed visually for precipitation and phase immiscibility at several-minute exposure; the obtained mixture was filtered through a 100-mesh sieve with registration of presence/absence of nonfilterable residue. Table 2 illustrates some of the laboratory findings on the compatibility of the acid composition under test with the fluids collected from the target well.

Prepared Mixture Composition of Interacting Products
Water-25% Oil-25% AC-50% Water-33.3% Oil-33.33% AC-33.34% Water-0% Oil-50% AC-50% As a result of laboratory research on compatibility of formation fluids with acid composition, it can be concluded that formation fluids and acid composition are well compatible, as in all cases, there were no signs of residue formation during filtration of the obtained mixtures. Therefore, FLUXOCORE-210 composition can be considered applicable to these conditions by the nature of interaction with the fluids.
To determine the predicted change in the skin factor from Equation 10, a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the To determine the predicted change in the skin factor from Equation (10), a number of laboratory studies to obtain the initial data for the calculation were conducted. The mass fractions of calcite and dolomite contained in the rock were found to be 89.7 and 1.2%, respectively, which are typical of the Tournaisian stage in Perm Krai. The viscosity value of FLUXOCORE-210 acid composition, determined as the average of several experiments, was 1.37 × 10 −3 Pa·s. The radius of the diffusing AC particle, as determined with Mikmed-5 microscope, was 1.0 × 10 −6 m. This radius is defined as the average of several radii of the particles under study, their values changing within a small range. Figure 5 shows the image of AC particles under microscope as obtained by the authors during the research in the scientific laboratory.
Tournaisian stage in Perm Krai. The viscosity value of FLUXOCORE-210 acid composition, determined as the average of several experiments, was 1.37 × 10 −3 Pa•s. The radius of the diffusing AC particle, as determined with Mikmed-5 microscope, was 1.0 × 10 −6 m. This radius is defined as the average of several radii of the particles under study, their values changing within a small range. Figure 5 shows the image of AC particles under microscope as obtained by the authors during the research in the scientific laboratory. Data on reservoir pressure provides an understanding of the required rate of acid injection into the formation. The HAT process was often carried out using AzINMASh-30 pumping unit. Its technical characteristics are given in Table 3 (Iolchuev and Savenyuk, 2018 [20]; Mishchenko, 2003 [21]). According to Table 3, the acid composition pumping rate was determined as 6.5 × 10 −3 m 3 /s. Bottom-hole temperature was 28 °C, as determined from telemetry data. The well diameter was 146 × 10 −3 m.
According to geophysical and hydrodynamic studies, the formation operated by the target oil production well is characterized by strong heterogeneity and low permeability. As a result, a largevolume hydrochloric acid treatment (LHAT) is recommended for the best results. In a work (Musabirov et al., 2019 [22]), it is reported that the volume of acid composition per 1 m of treatment interval varies in the range from 3 to 5 m 3 when the LHAT technology is applied for carbonate deposits. The perforated formation thickness in this target well reaches 8 m, therefore, it is necessary to inject acid in the volume of 40 m 3 to maximize the BHZ permeability.
When all of obtained values are substituted in Equation 10, the calculated predicted decrease of the skin factor is 4.0, which has a significant impact on changes in the amount of fluid recovered from Data on reservoir pressure provides an understanding of the required rate of acid injection into the formation. The HAT process was often carried out using AzINMASh-30 pumping unit. Its technical characteristics are given in Table 3 (Iolchuev and Savenyuk, 2018 [20]; Mishchenko, 2003 [21]). According to Table 3, the acid composition pumping rate was determined as 6.5 × 10 −3 m 3 /s. Bottom-hole temperature was 28 • C, as determined from telemetry data. The well diameter was 146 × 10 −3 m.
According to geophysical and hydrodynamic studies, the formation operated by the target oil production well is characterized by strong heterogeneity and low permeability. As a result, a large-volume hydrochloric acid treatment (LHAT) is recommended for the best results. In a work (Musabirov et al., 2019 [22]), it is reported that the volume of acid composition per 1 m of treatment interval varies in the range from 3 to 5 m 3 when the LHAT technology is applied for carbonate deposits. The perforated formation thickness in this target well reaches 8 m, therefore, it is necessary to inject acid in the volume of 40 m 3 to maximize the BHZ permeability.
When all of obtained values are substituted in Equation (10), the calculated predicted decrease of the skin factor is 4.0, which has a significant impact on changes in the amount of fluid recovered from this well. The conducted calculation of the predicted value of the skin factor after acid treatment of the target facility under study, which equals S = 1.1 − 4.0 = −2.9, allows to perform the comparative analysis between the actual values of the well oil and fluid flow rate obtained over 6 months and those predicted using hydrodynamic simulation, as given in Table 4. Figure 6 illustrates the hydrodynamic model. this well. The conducted calculation of the predicted value of the skin factor after acid treatment of the target facility under study, which equals S = 1.1 − 4.0 = −2.9, allows to perform the comparative analysis between the actual values of the well oil and fluid flow rate obtained over 6 months and those predicted using hydrodynamic simulation, as given in Table 4. Figure 6 illustrates the hydrodynamic model.  According to Table 4, the predicted oil and liquid flow rates almost fully correspond to their actual values, since the average deviations were 5.33% and 5.12%, respectively, while the standard deviations were 0.67 and 0.74 with variation coefficients of 12.58% and 14.47%, respectively. Therefore, it is possible to note high accuracy of predicting the flow rates according to the developed technique, which will allow evaluating effectiveness of acid treatment of productive deposits that are lithologically represented by carbonate varieties.
At further calculations using the hydrodynamic model, it is possible to analyse incremental recovery of hydrocarbons for several years to come, i.e., to make a conclusion on the success of the  According to Table 4, the predicted oil and liquid flow rates almost fully correspond to their actual values, since the average deviations were 5.33% and 5.12%, respectively, while the standard deviations were 0.67 and 0.74 with variation coefficients of 12.58% and 14.47%, respectively. Therefore, it is possible to note high accuracy of predicting the flow rates according to the developed technique, which will allow evaluating effectiveness of acid treatment of productive deposits that are lithologically represented by carbonate varieties.
At further calculations using the hydrodynamic model, it is possible to analyse incremental recovery of hydrocarbons for several years to come, i.e., to make a conclusion on the success of the oil recovery stimulation method. Below is the analysis of the effectiveness of the designed procedure at the target facility in comparison with the basic variant using hydrodynamic simulation; its results are given in Table 5. The calculation was conducted through 2029.
Analysing the data presented in Tables 4 and 5, we can conclude that following the acid treatment of the target oil producing well, a significant growth in oil production is expected, namely, a fourfold increase. Based on the conducted forecast, by the beginning of 2029 at the target oil producing well, the cumulative oil production will increase by 9.88 thousand m 3 , and the cumulative liquid production will increase by 50.25 thousand m 3 , following the performed HAT job.

Conclusions
The evaluation technique of the hydrochloric acid treatment effectiveness proposed herein allows to predict the technological effectiveness of the designed procedure. The technique consists in determining the technological effectiveness using hydrodynamic simulation based on the calculated value of change in the skin factor using the initial data obtained from laboratory studies. Application of the proposed approach will allow to raise accuracy of the technical and economic assessment carried out at decision-making on conducting a geotechnical procedure. The technique also includes parameters enabling to provide recommendations on the selection of an acid composition and technology of its injection into the formation for the best effect from acid treatment of carbonate reservoirs. This approach to predicting incremental oil production was used to compare expected and actual well oil flow rates measured over 6 months of continuous operation after hydrochloric acid treatment. According to the analysis results, there is a minor error (≈ 5%) in the predicted values and actual readings of the target well, i.e., the proposed technique of effectiveness evaluation has been successfully tested, which proves the potential of its integration in engineering activities of oil companies.
Funding: This research was funded by the state assignment of the Ministry of Science and Higher Education of the Russian Federation as part of a government assignment, grant number FSNM-2020-0027.

Conflicts of Interest:
The authors declare no conflicts of interest. Boltzmann constant, J/K T absolute reservoir temperature, K µ dynamic viscosity of acid composition, Pa·s r diffusing particle radius, m S 1 skin factor before HAT Q 1 well oil flow rate before HAT, m 3 /day ∆Q cum increase in cumulative oil production, t