Research and application of acid fracturing stimulation mechanism in ultra-deep subsalt dolomite reservoir in Tarim Basin

Deep and ultra-deep carbonate reservoirs are the focus of exploration and development in future. However, the problems of high pressures in the treatment process, a limited effective etching distance of acid, great acid leak-off, and poor adaptability of the acid system are encountered in this type of oil and gas reservoir. The mechanism of acid fracturing stimulation under different processes and parameters is not clear. Aiming at these issues, the treatment schemes, process optimization, parameter optimization, and liquid system screening are studied in this paper, try to clarify the acid fracturing stimulation mechanism, and the following conclusions are drawn: The acid network fracturing could activate natural fracture to generate a complex fracture network to the greatest extent, and thereby a high output could be achieved; By using of weighted fracturing fluid, the wellhead injection pressure, as well as the performance of equipment required, could be effectively reduced; With 20% gelling acid and 20% retarded acid system, the non-uniform etching could be realized to improve the effective etching distance of acid liquid. The conclusions in this paper shed light on the acid fracturing treatment of deep and ultra-deep carbonate rocks.


List of Symbols C f
The hydraulic fracture compression coefficient (MPa −1 ) C L The compressibility coefficient of fluid (MPa −1 ) C m The comprehensive compressibility coefficient of the matrix (MPa −1 ) C n The compression coefficient of natural fracture (MPa −1 ) K f The hydraulic fracture permeability (D) K m The permeability of bedrock (D) K n The natural fracture permeability (D) L m The matrix rock block size (m) The pressure of hydraulic fracture system (MPa) p junc The intersection pressure (MPa) P m The pressure of bedrock pore system (MPa) P n The pressure of the natural fracture system (MPa) p net The net pressure in the fracture (MPa) p tip The fluid pressure at the tip of the natural fracture (MPa) q f The volumetric flow rate of liquid at the source divided by that at the sink per unit hydraulic fracture pore volume (s −1 ) q n The volumetric flow rate of liquid at the source divided by that at the sink per unit natural fracture pore volume (s −1 ) x o The distance between the restart and intersection points (m) The matrix shape factor (12/Lm 2 ) M − M � Is the delta function, of which the value equals 0 when M = M� Δp nf The pressure loss in the direction of the natural fracture (MPa) Δp open The pressure loss in the direction of the length of activated natural fracture l (MPa) Δp shear The pressure loss in the direction of shear natural fracture length a (MPa) The fluid viscosity (MPa s) σ hmin and σ Hmax The min and max far-field stress (MPa) σ n Normal stress (MPa) 0 The shear strength of the natural fracture (MPa) τ Tangential stress (MPa)

Introduction
In the Tarim oilfield, 23 subset ultra-deep dolomite wells have been exploited, and seven-time treatment measures have been implemented in 6 of these wells (Table 1). The acid fracturing process primarily applies to the reservoirs with acid corrosion rates higher than 90%; Sand fracturing technology is principally adopted for reservoirs with low acid corrosion rates or containing gypsum; Generally speaking, currently, the treatment technologies of subsalt ultradeep dolomite reservoirs do not match with the real situations, and the treatment performance is not good enough. The effective exploitation of subsalt ultra-deep dolomite reservoir is of great importance for the high stable output of the Tarim oilfield.
Research on dolomite reservoir treatment has been reported. Yan carried out acid-rock reaction kinetics experiments with carbonate rock cores containing different dolomite mass fractions (Yan et al. 2021b). The difference in acid corrosion mechanisms of limestone and dolomite was analyzed with the scanning electron microscope, and the models of dolomite acid-rock reaction kinetics determined by different factors were established. Dong et al. (2021) constructed a three-dimensional acid migration model to calculate the geometry of the fracture in acid fracturing treatment. Acid flow, reaction migration, and rock corrosion in fracture were taken into account in this model. The influence of dolomite heterogeneity of the acid corrosion process was analyzed, and the acid corrosion modes in carbonate reservoirs with different spatial distributions of calcite and dolomite minerals were compared. Nino-Penaloza conducted acid-injecting tests on fractured cores of the chalk rock strata in Austin and dolomite strata in San Andres to probe the effects of injection temperature, pumping time, and acid type on the conductivity of fracture. And, the influence of core fracture surface texture before acid injection was investigated with rough and smooth surface textures as the research objects for quantification (Nino-Penaloza et al. 2015). Qi recommended that the user of retarded acid hardly led to high conductivity in the acid fracturing of the dolomite reservoir (Qi et al. 2021). Therefore, according to the acid-etching morphology of fractures, the influencing factors and improvement methods of fracture conductivity were discussed. It was suggested that the etching morphology was not directly related to the determining mechanism of the acid-rock reaction, but dependent on the reaction rate. The reason for low conductivity was the slow reaction between dolomite and retarded acid and the difficult formation of irregular etching grooves on the fracture surface. A small amount of dolomite in the reservoir did not affect the conductivity, but > 80%  (Li et al. 2012). Wang studied the fine segmented design method and developed the high-temperature resistant authigenic acid and gelling acid for ultra-deep carbonate rocks, aiming at the problems of fast reaction between acid and rock, short action distance of acid and low fracture conductivity under high closure stress ). In addition, the use of encapsulated acid is also an important means to achieve high-efficiency treatment of deep carbonate rocks (Luo et al. 2019a, b). Generally speaking, the primary problems of dolomite reservoir treatment include great acid leak-off caused by the development of natural fractures (Fu et al. 2019), fast acidrock reaction caused by high reservoir temperature (Rabie et al. 2014;Sayed et al. 2013), difficult maintenance of conductivity in acid-corroded fractures caused by high Young's modulus, large Poisson's ratio, and high crustal stress (Desouky et al. 2020;Mehrjoo et al. 2022), poor conductivity caused by uniform etching of dolomite with high purity (Zhang et al. 2020), difficult vertical effective operation caused by numerous vertical small layers and thick reservoirs, and difficult fracture caused by great fracture pressure gradient of the reservoir (Yan et al. 2021a). Therefore, the corresponding treatment measures were proposed as follows: liquid thickening , leak-off suppression with solid particles (Pankov et al. 2012), pre-flush fluid cooling (Lei et al. 2018;Zhang et al. 2022), large-displacement treatment (Han et al. 2017), closed acidification (Aljawad et al. 2019), sand-acid fracturing (Zhang et al. 2021), diversion acid fracturing (Rabie et al. 2012), and so on.
However, the ultra-deep subsalt dolomite reservoir in Tarim Basin is significantly different from other dolomite reservoirs, with higher temperatures, higher formation pressures, probable secondary precipitation caused by gypsum, and greater treatment difficulty (Hu et al. 2021;Peng et al. 2018;Yan et al. 2021c). Therefore, it is necessary to propose targeted treatment strategies based on a deep analysis of the reservoir characteristics and treatment difficulties, to form appropriate treatment processes, and optimize treatment parameters, for the efficient development of the ultra-deep subsalt dolomite gas reservoir in Tarim.

Difficulties and measures of reservoir treatment
In the Tarim oilfield, 23 subsalt ultra-deep dolomite wells were exploited, and seven treatment measures were implemented for 6 of these wells. The treatment process included acid fracturing by multi-stage alternating injection, temporary plugging followed by sand fracturing, and acidizing. Three treatment measures were implemented for the dolomite reservoir, and another three measures were implemented for the gypsum-containing dolomite reservoir, with the treatment well depth range majority of about 8000 m, reservoir temperature range of 130-170 °C, displacement range of 2.5-7.1 m 3 min −1 , total liquid volume range of 140-1310 m 3 , and extension pressure gradient range of 1.83-2.46 MPa/100 m. After the treatment, the output of wells in different formations and different processes varied greatly, and the treatment performance was not satisfactory. The treatment difficulties are the following (Cai et al. 2009;Su et al. 2019;Yang et al. 2020a): 1. Treatment treatment pressure is high. The reservoir is deep, and the fracture pressure gradient and fracture extension pressure gradient are high. Therefore, the treatment is difficult and requires more rigid treatment string and injection equipment. 2. The effective working distance of acid is limited. The reservoir temperature is high in the range of 130-170 °C, and thereby the reactions between acid and reservoir rock are fast. A large proportion of acid is consumed at the fracture opening, and the length of acid-corroded fracture is small. 3. Reservoir fractures and cavities are developed, acid fluid leak-off is obvious during fracturing, and liquid efficiency is low. The leak-off balanced and microfracture leak-off zigzag types are dominant in the gypsumcontaining dolomite acid fracturing wells. The etching ability of acid for the fractured wall is limited, and the acid-corroded fracture is narrow, resulting in high net pressure in the fracture; The fracture-height-out-of-control, connected-natural-fracture, and leak-off balanced types are dominant in the acid fracturing dolomite reservoir. The net pressure in the fracture is constant or declines gradually. 4. The effluent is serious after acid fracturing treatment.
The buried depths of some wells exceed the oil-water/ gas-water interfaces, so water is easy to produce after treatment. In contrast, the gas output of some wells with the same depths exceeds 15 × 10 4 m 3 d −1 without water production. 5. The adaptability of the acid system is poor. The treatment success rate of gypsum-containing dolomite reservoir is low, and treatment is difficult because the existing treatment liquid or process is not suitable. The gypsumcontaining dolomite reservoir is acidified with ordinary gelling acid with a low corrosion rate, leading to the reinforcement of rock plasticity and difficult fracture initiation and propagation; the non-uniform etching of gypsum-containing dolomite is of low efficiency in the presence of the acid, and the corrosion rate, as well as conductivity of acid-etching fracture, are low, leading to difficult prop at a high closing pressure and rapid decline of output.
Targeting these treatment difficulties, the following treatment countermeasures are proposed: 1. CaCl 2 weighted fracturing fluid is selected as the preflush fluid for fracturing and opening the reservoir in the early stage, and KCl weighted fracturing fluid is selected in the later stage to improve the hydrostatic pressure of the fluid in the well and to reduce the risk of wellhead overpressure, so the requirements of pipe string and injection equipment will be less stringent. Adding CaCl 2 and KCl, the density of fracturing fluid can be increased and the liquid column pressure can be increased, thus avoiding increasing the bottom hole pressure by increasing the wellhead pressure (Mao et al. 2022;Yang et al. 2020b). 2. First, a large amount of pre-flush fluid is injected to reduce the reservoir temperature and to slow down the reactions between acid and reservoir rock. 3. The gel fracturing fluid with high viscosity is selected as the pre-flush fluid to reduce the leak-off of acid. The multi-stage alternating injection process is adopted for the injection course to reduce the leak-off of acid, improve the liquid efficiency, and increase the working distance of acid. 4. On the premise of ensuring the treatment safety and integrity of the whole life cycle of the well, the displacement is gradually increased to prevent the breaking of the reservoir and connection of aquifers. 5. Multi-stage alternating injection and closed acid fracturing treatment process are adopted. Cross-linked acids are preferable, and acids with high and low viscosities are alternately injected to realize the deep etching of near and far well zones. Closed acid fracturing is adopted for near wells to improve the fracture conductivity; In the early stage of acid injection, acid with low viscosity is used for etching to activate the natural fractures around the well.

Optimization of the treatment process
Based on the reservoir characteristics, treatment difficulties, and treatment countermeasures, the treatment process is optimized, and the fracture length and conductivity are optimized.

Treatment process optimization
The output in different treatment processes was simulated by using the dynamic analysis model of volume fracturing horizontal well. The seepage differential equation of singlephase compressible micro-fluid in the bedrock when starting pressure gradient can be expressed as follows Ren et al. 2015): where P m is the pressure of bedrock pore system, MPa; is the starting pressure gradient of bedrock, MPa m −1 ; C L represents the compressibility coefficient of fluid, MPa −1 ; m is matrix porosity, dimensionless; denotes fluid viscosity, mPa s; C m represents the comprehensive compressibility coefficient of the matrix, MPa −1 ; K m represents the permeability of bedrock, D; is matrix shape factor, 12/L m 2 ; L m is matrix rock block size, m; P n is the pressure of natural fracture system, MPa. The natural fracture system agreed with the Darcy seepage law. Similarly, the control equation of the natural fracture system can be expressed as follows: where C n represents the compression coefficient of natural fracture, MPa −1 ; K n represents the natural fracture permeability, D; q n represents the volumetric flow rate of liquid at the source divided by that at the sink per unit natural fracture pore volume, s −1 ; n represents the porosity of natural fracture, dimensionless; M − M � is the delta function, of which the value equals 0 when M = M′.
As an independent medium system, the seepage of fluid in hydraulic fracture obeys Darcy's law, so its control equation can be expressed as follows: where C f represents the hydraulic fracture compression coefficient, MPa −1 ; K f represents the hydraulic fracture permeability, D; q f represents the volumetric flow rate of liquid at the source divided by that at the sink per unit hydraulic fracture pore volume, s −1 ; f represents the porosity of natural fracture, dimensionless; P f represents the pressure of hydraulic fracture system, MPa.
The dimensions in Eqs.
(1)-(3) were eliminated by mathematical methods, and a finite-element integral equation was established by using the weighted residual method. The continuous infinite-freedom-degree solving elements were discretized into finite elements for the solution, and thereby the output under different fracture conditions can be simulated (Ren et al. 2015).
Targeting the maximizing of output, the treatment schemes of the Cambrian dolomite reservoir, including sand fracturing, acid network fracturing, and conventional acid fracturing, were screened (Zhao et al. 2014); The main parameters used are shown in Table 2, the maximum daily and cumulative output values in different processes are shown in Fig. 1.
The maximum daily and cumulative output in different processes varied greatly. The calculation results show that the reservoir treatment with the acid network fracturing technique could activate the natural fractures to the greatest extent to form a complex fracture network (Fig. 2), so the daily and cumulative output could be maximized. Therefore, the acid network fracturing technique is preferable. During treatment, an acid system with low viscosity should be injected to activate the natural fractures around the well in the early stage of injection; Then, different acid systems with low and high viscosities should be alternately injected to realize the deep etching of near and far wells. The alternate injection of high viscosity acid and low viscosity acid with two purposes: (1) Alternate injection of high and low viscosity acid fluid to form a finger phenomenon, improve the degree of non-uniform etching, and then improve the conductivity; (2) Etching can be realized in the whole injection process, avoiding the defect of low etching efficiency when fracturing fluid is used alternately.
The pressure conditions for the activation and extension of natural fractures were analyzed by using the interaction criterion between artificial and natural fractures (Luo et al. 2018). Warpinski and Teufel proposed the shear slip criterion for intersection interaction perpendicular to natural fractures, and this criterion can be expressed by Eq. (4): The intersection pressure p junc is replaced by the net pressure p net in the fracture, and the stress τ and σ n can be expressed by the far-field stress σ hmin and σ Hmax . The intersection criterion can be expressed as follows: where 0 is the shear strength of the natural fracture.
(4)  The intersection interaction criterion can be expressed in the form of net pressure Pnet: And, where Δp open is the pressure loss in the direction of the length of activated natural fracture l; Δp shear is the pressure loss in the direction of shear natural fracture length a, and x o represents the distance between the restart and intersection points. The criterion for the fracture restart at the tip of activated natural fracture can be described as follows: where p tip is the fluid pressure at the tip of the natural fracture, and Δp nf represents the pressure loss in the direction of the natural fracture.
Similarly, Eq. (11) can be expressed as follows:p net ≥ 1 2 H max − h min (1 − cos 2 ) + T o + Δp nf (12) The horizontal stress difference of Wusonger formation in Tarim subsalt dolomite was 29 MPa, and the azimuth angle between natural fracture and principal stress was around 15°; The horizontal stress difference of Shayirik formation was 33 MPa, and the azimuth angle between natural fracture and principal stress was about 20°; The pressure conditions for the activation and extension of natural fractures in Fig. 2 show that the net pressure required to activate natural fracture increased with the increase in the angle and horizontal stress difference. The maximum activated net pressure was almost 7 MPa, and the natural fractures with great stress difference and large angles were difficult to activate; The extension pressure after the opening and activation of natural fractures was in the range of 2-3.5 MPa, at a medium level, indicating that the activation of natural fracture was vital to the treatment of fracture network. Considering the limited pump pressure during treatment, chemical corrosion with acids could be synergistically combined with mechanical interaction to activate natural fractures. On the one hand, some carbonate fillings can be dissolved through the pretreatment of the acid, on the other hand, the mechanical properties of the reservoir rocks can be reduced, thereby reducing the fluid pressure that activates the natural fracture.

Fracture length
With the preferable technique of acid network fracturing, the output change after single-stage treatment and compaction was simulated. The fracture lengths of single-stage fractures were set to 100, 120, 150, and 180 m, and the effectively closed conductivity was set to 30 D cm. Twenty natural fractures were activated. The daily and cumulative output under the different conditions of fracture lengths was calculated, and the results are shown in Fig. 3.
When the fracture length is 100, 120, 150, and 180 m, respectively, the initial daily output of oil is about 105, 127,147, 163 t d −1 , the stable daily output is about 9, 13, 17, 20 t d −1 , and the cumulative oil is about 0.64 × 10 5 , 0.8 × 10 5 , 0.93 × 10 5 , 1.04 × 10 5 t. The increase in cumulative oil is 0.16 × 10 5 , 0.13 × 10 5 , and 0.11 × 10 5 t, respectively. The longer the fracture was, the larger the treatment range was, and the smaller the seepage resistance was. Therefore, the daily and cumulative output increased. When the effective fracture length exceeded 150 m, the output increased to a smaller extent. Therefore, the fracture length value of 150 m was preferable.

Conductivity capacity
The conductivity of artificial fractures is another important index for the evaluation of treatment performance, and it affects the production performance after treatment. It is necessary to optimize the conductivity of fracture after the acid network fracturing process. In the simulation and calculations, the length of the single-section fracture was set to be 150 m, twenty natural fractures were activated. The daily and cumulative output under the different conditions of conductivity capacity of 10, 20, 30, and 40 D·m was calculated, and the results are shown in Fig. 4.
When the conductivity is 10, 20, 30, and 40 D cm, the initial daily output of oil is about 90, 123, 147, and 158 t d −1 , respectively, and the cumulative oil is 0.66 × 10 5 , 0.8 × 10 5 , 0.93 × 10 5 , and 0.99 × 10 5 t, respectively, with an increase of 0.14, 0.13, and 0.07 t, respectively. The greater the conductivity capacity was, the smaller the flow resistance was, and the greater the daily and cumulative output were; the cumulative output with conductivity capacity in the range of 30-40 D cm was close to the single-section controlled output; when the conductivity capacity exceeded 30 D cm, the contribution of increased conductivity capacity to the output was weakened. Therefore, the conductivity capacity value of 30 D cm was preferable.

Evaluation of the performance of acid fracturing materials
According to the optimization results of process parameters, the performances of acid, fracturing fluid, and etching were evaluated separately to judge whether the materials met the treatment requirements.

Fracturing fluid system
For the reduction in wellhead treatment pressure and the harshness of requirements for the performance of injection string and equipment for the dolomite reservoir, the wellbore hydrostatic pressure should be increased (Dai et al. 2019;Liu et al. 2020). Therefore, CaCl 2 (first section) and KCl (second and third sections) weighted fracturing fluids were selected for alternate injection.
The viscosity of CaCl 2 and KCl weighted fracturing fluids was measured with a Huck MARSIII rotary rheometer with a shear rate of 170 s −1 and heating rate of 1 °C min −1 . The relationship between sample viscosity and temperature was recorded. The test results are shown in Figs. 5 and 6.
The initial viscosity of CaCl 2 weighted fracturing fluid was about 90 mPa s. With the increase in temperature, the viscosity first increased and then declined. 30 min later, the viscosity reached a maximum of about 500 mPa s at about 140 °C. After that, the viscosity gradually declined to about 200 mPa s, and the leak-off of the acid system to natural fractures could be impeded to a large extent. In contrast, the rheological profile of KCl weighted fracturing fluid was different from that of CaCl 2 weighted fracturing fluid. The initial viscosity of KCl weighted fracturing fluid was about 950 mPa s. With the increase in temperature, the viscosity gradually declined. 30 min later, the viscosity declined to a value of about 200 mPa s and then was kept constant. According to the rheological characteristics of CaCl 2 and   KCl weighted fracturing fluids, the former was selected as the fracturing fluid in the first stage. Its high viscosity suppressed leak-off and promoted the formation of high pressure at the well bottom to open the reservoir and generate fractures; the latter was selected as the fracturing fluid in the form of alternate injection in the following stages to reduce reservoir temperature and acid fluid leak-off.

Acid system
Gelling and retarded acid were the candidates for acid systems. According to the industry standard SY/T 5886-2012 "Evaluation approach of retarded acid performance," the corrosion rate, etching morphology, and conductivity capacity of acid-etching fracture of 20% gelling acid and 20% retarded acid (pH = 4-5) for gypsum-containing dolomite were tested. The results are shown in Table 3.
The corrosion rates of gelling acid and retarded acid for the gypsum-containing dolomite were similar; the corrosion rate of gelling acid was in the range of 44.42-61.34%, and the average corrosion rate was 51.72%; the corrosion rate of retarded acid was in the range of 45.19-55.93%, with a mean value of 50.56%.
The etching for gypsum-containing dolomite was carried out by using 20-min displacement experiments. The etching morphologies of gypsum-containing dolomite with gelling acid and retarded acid are exhibited in Figs. 7 and 8. The former showed poor etching performance for gypsumcontaining dolomite, and only local pitting corrosion was observed. The conductivity test was carried out after etching. When the closing pressure was in the range of 5-10 MPa, almost no outflow was observed at the outlet, and high pressure was held at the inlet; In contrast, the rock plate showed a good etching morphology with local etching grooves with retarded acid. Figure 9 shows the conductivity of the two acid liquids under different closing pressures. With the increase in closure pressure, the conductivity decreases rapidly. However, since retarded acid can form better non-uniform etching than that of gelling acid (Figs. 7 and 8), the conductivity of retarded acid etched fracture is still higher than that of gelling acid. Under the closing pressure of 5 MPa, the conductivity of retarded acid etched fracture is about 9 D cm, while that of gelled acid etched fracture is about 7.5 D cm. Under the closing pressure of 60 MPa, the conductivity of retarded acid etched fracture is about 1.5 D cm, while that of gelled acid etched fracture is about 0.5 D cm.

Field application
Taking well Z2 as an example, the field experiment of acid fracturing technology of ultra-deep dolomite reservoir was performed. Well Z2 was 8791 m deep, and the acid fracturing treatment was carried out in the 8376-8406-m well section, where the reservoir texture was granular and finegrained dolomite and contained fracture pores mostly semifilled with calcite and dolomite; This section consisted of 2  layers of the class-II reservoir with a thickness of 6.5 m and porosity of 3-3.2% (mean value is 3.1%) and 6 layers of the class-III reservoir with a thickness of 39.5 m and porosity of 1.3-1.8% (mean value is 1.61%); The well bottom was at about 157 °C; The treatment section was at 8648-8683 m, with high crustal stress. The minimum horizontal principal stress was in the range of 174.9-180.8 MPa, and the horizontal stress difference was in the range of 38.9-44.8 MPa.
The brittleness index was greater than 65, and the average compressibility index was 0.5, indicative of high compressibility (Yan et al. 2019); 30 natural fractures were developed in the treatment section, and most of the dip angles of fractures were in the range of 55°-85°, so these fractures were ready to activate. According to the results of reservoir characteristics, reservoir treatment difficulties and countermeasures, treatment process and fracture parameter optimization, and acid fracturing material performance evaluation, the injection procedure was designed and is shown in Table 4.
According to the pumping procedure shown in Table 4, the field implementation was carried out. The treatment profile is shown in Fig. 10. The injection commenced at 17:45, and the pressure and displacement gradually increased. At 18:7, under the condition of constant displacement, the oil pressure continued to decrease, confirming that the fracture was opened and the fracture was extended normally. Before the decrease in oil pressure, the maximum oil pressure reached 113 MPa, and the well bottom pressure was evaluated to be about 180.4 MPa; At 18:35, the displacement increased and oil pressure increased accordingly. In the subsequent injection process, the decrease in oil pressure was observed many times under the condition of constant displacement, indicating that multiple natural fractures were connected.
Based on the data of pressure and displacement injected during the treatment, the net pressure fit was conducted, and the geometry of the fracture was simulated through inversion. The simulated fracture was 145 m half-length and 49 m high, which are close to the values (140 and 55 m, respectively) required in the acid fracturing design. Therefore, the treatment was carried out smoothly and successfully.
As shown in Fig. 11, before the acid fracturing treatment (April 2, 2022), a 4-mm nozzle was used for the drainage of well Z2. The oil and casing pressures were 0 and 19.012-18.697 MPa, with a total drainage volume of 0.2 m 3 and without natural gas; The acid fracturing treatment was carried out on April 3. After the treatment, the 4 mm nozzle was still used. The oil pressure was in the range of 18. 524-19.613 MPa, equivalent to about 6000 m 3 of gas output per day, indicative of good treatment performance.

Conclusion
The problems of high pressures in the treatment process, a limited effective etching distance of acid, great acid leakoff, and poor adaptability of the acid system are encountered in deep and ultra-deep carbonate reservoirs. The mechanism of acid fracturing stimulation under different processes and parameters is not clear. Based on analyzing the difficulties of dolomite reservoir treatment and the mechanism of different stimulation methods, the corresponding optimal stimulation method is proposed, the mathematical model of ultra-deep reservoir productivity is established, and the optimization design of injecting process, parameters, and acid fracturing materials is carried out, the following conclusions are obtained: 1. The stimulation mechanism under different treatment methods, injection parameters, and acid fracturing materials is analyzed through numerical simulation and experimental tests, which provides a theoretical basis for the on-site implementation of acid fracturing of the dolomite reservoir. 2. The application of weighted fracturing fluid, pre-flush fluid cooling, leak-off reduction in high-viscosity fracturing fluid, and optimization of treatment processes are crucial to acid fracturing treatment. 3. The application of weighted fracturing fluid could effectively reduce the wellhead pressure, and reduce the harshness of equipment performance required. 4. By using 20% gelling acid and retarded acid, non-uniform etching could be realized to improve the effective action distance of acid and effective conductivity.
Funding No funding support.

Conflict of interest
On behalf of all the co-authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.