Effects of Dolomitization on Porosity during Various Sedimentation-Diagenesis Processes in Carbonate Reservoirs

Carbonate reservoirs, especially dolomite reservoirs, contain large reserves of oil and gas. The complex diagenesis is quite challenging to document the dolomite reservoirs formation and evolution mechanism. Porosity development and evolution in dolomite reservoirs primarily reflect the comprehensive effect of mineral dissolution/precipitation during dolomitization. In this study, multicomponent multiphase flow and solute transport simulation was employed to investigate dolomitization in the deep carbonate strata of the Tarim Basin, Northwest China, where active exploration is currently under way. Oneand two-dimensional numerical models with various temperatures, fluid compositions and hydrodynamic characteristics were established to quantificationally study dolomitization and its effect on porosity. After determining the main control factors, detailed petrologic characteristics in the studied area were also analyzed to establish four corresponding diagenetic numerical models under different sedimentary environments. These models enabled a systematic analysis of mineral dissolution/precipitation and a quantitative recovery of porosity evolution during various sedimentation-diagenesis processes. The results allowed for a quantitative evaluation and prediction of reservoir porosity, which would provide a basis for further oil and gas exploration in deep carbonate reservoirs.


Introduction
Carbonate reservoirs are widely distributed over the world and crucial for oil and gas exploration because of their proven reserves. However, carbonate reservoirs are usually developed with complicated crack/pores during diagenesis, which make the reservoir prediction difficult and hinder the development of oil and gas exploration [1][2][3][4][5][6][7] It is of great importance to understand the effect of diagenesis on the reservoir porosity evolution as it is one of the key factors affecting the reservoir quality.
The formation and evolution of reservoir porosity are the topic of carbonate reservoir research and constitute the premise for the efficiency improvement of oil and gas exploration [5,[8][9][10]. In carbonate reservoirs, reservoir porosity evolution is mainly the result of carbonate minerals dissolution and precipitation, such as dolomitization [11][12][13][14]. Fluid-rock reactions are usually the main factor for the porosity evolution because it results in carbonates dissolving to form voids, although in some cases minerals may precipitate to fill such voids [2,15,16]. However, the reservoir condition although in some cases minerals may precipitate to fill such voids [2,15,16]. However, the reservoir condition (e.g., temperature, pressure, mineral species, formation water hydrochemistry and fluid dynamics) is crucial to the porosity evolution during the long process of diagenesis [9,[17][18][19][20][21][22].
As the most important carbonate reservoir, dolomite has received wide attention [13,[23][24][25]. With the continuous exploration and discovery of oil and gas in deep-buried dolomite formations, a string of new issues has arisen. The most important issue is how dolomitization affects the reservoir porosity during various sedimentation-diagenesis processes [26][27][28][29]. This, however, has not received notable attention as the majority of previous research works focused on qualitative testing and analysis of dolomitization, which did not provide a systematic description of carbonate minerals under different burial conditions during the whole diagenetic process [30].
In this study, the deep strata of the Tarim Basin in Northwestern China were selected to establish a number of one-and two-dimensional numerical models under various diagenetic conditions. This allowed us to analyze the effect of temperature, fluid and hydrodynamic conditions on dolomitization, and for the quantitative evaluation of the reservoir porosity change induced by dolomitization. Then, quantitative recovery of porosity evolution during various sedimentation-diagenesis processes was successfully implemented. The results allowed for a quantitative evaluation and prediction of reservoir porosity, which would provide a basis for further oil and gas exploration in deep carbonate reservoirs.

Tectonic Locations
Tarim Basin is one of the most important petroliferous basins in China, which contains thick marine carbonate deposits. The most well-developed carbonate sequence within the basin is the Cambrian-Ordovician carbonate, which is also the key intervals for oil and gas exploration [31][32][33]. The Ordovician carbonate reservoirs in the Shunnan (SN) area, located in the central uplift of the Tarim Basin, are selected in this study, as shown in Figure 1, based on [33,34].

Diagenesis
A good overview of the typical diagenesis in the studied area has been provided by various previous research, as shown in Figure 2 [32,33,35,36] A series of core and thin sections and other petrological analyses show that the diagenesis mainly includes micritization ( Figure.

Diagenesis
A good overview of the typical diagenesis in the studied area has been provided by various previous research, as shown in Figure 2 [32,33,35,36] A series of core and thin sections and other petrological analyses show that the diagenesis mainly includes micritization (Figure 2a reservoirs. Liu [39] theorized that the Middle-Lower Ordovician dolomites in the SN area are mainly resulted from shallow burial, burial (transitional environment) and deep burial dolomitizations. The fine-grained dolomite strata formed during the burial stage have the best reservoir properties, causing grain bank metasomatism and sustained dolomitization of the mud-sized dolostone formed under a shallow burial environment [37].
The diagenetic evolution sequence of the reservoirs in the studied area is summarized in Figure 3. According to the diagenetic history, the diagenetic evolution process in the studied reservoir can be categorized into six consecutive stages, including sedimentary-parasyngenetic, parasyngenetic-shallow burial, supergene, shallow burial, middle-deep burial and the deep burial stages.   The diagenetic evolution of the Ordovician reservoirs in the Tarim Basin is complex, which mainly experience submarine diagenesis, meteoric diagenesis and the deep burial diagenetic settings. It is characterized by an early weak cementation, a late compaction and pressure solution, formation of stylolite, dissolution and karstification, mixed dolomitization in the supergene stage and the rupture and formation of multiphase fractures, which are all favorable conditions for developing favorable reservoirs [28,37].
Several studies suggest that the syngenetic karst and inter-stratal karst of the upper highstand system under the Middle-Lower Ordovician III interface are controlled by paleotopography and paleo-fault [38]. Such eogenetic karst reservoirs have a "quasi-layered" development model, and are filled by the subsequent formation of carbonate cement. In the late diagenetic stage, acidic hydrothermal fluids with medium-low temperature, high salinity and rich SiO 2 penetrated along the fault zones, leading to hydrothermal alterations such as hydrothermal dissolution, silica replacement and quartz Minerals 2020, 10, 574 4 of 18 precipitation to become the main controlling factors for the formation of reservoirs. Liu [39] theorized that the Middle-Lower Ordovician dolomites in the SN area are mainly resulted from shallow burial, burial (transitional environment) and deep burial dolomitizations. The fine-grained dolomite strata formed during the burial stage have the best reservoir properties, causing grain bank metasomatism and sustained dolomitization of the mud-sized dolostone formed under a shallow burial environment [37].
The diagenetic evolution sequence of the reservoirs in the studied area is summarized in Figure 3. According to the diagenetic history, the diagenetic evolution process in the studied reservoir can be categorized into six consecutive stages, including sedimentary-parasyngenetic, parasyngenetic-shallow burial, supergene, shallow burial, middle-deep burial and the deep burial stages.

Sedimentary Environments
The diagenetic evolution of the Ordovician carbonate rocks in the SN area occurred under different sedimentary environments. According to the diagenesis and petrological characteristics, the sedimentary environment can be divided into four categories: (1) marine phreatic environment in a subtidal zone, (2) meteoric fresh water and fresh water phreatic environment in a supralittoral zone, (3) seawater and freshwater interaction environment in an intertidal zone and (4) seawater evaporation environment in a shoal setting.

Sedimentary Environments
The diagenetic evolution of the Ordovician carbonate rocks in the SN area occurred under different sedimentary environments. According to the diagenesis and petrological characteristics, the sedimentary environment can be divided into four categories: (1) marine phreatic environment in a subtidal zone, (2) meteoric fresh water and fresh water phreatic environment in a supralittoral zone, (3) seawater and freshwater interaction environment in an intertidal zone and (4) seawater evaporation environment in a shoal setting.

Model Tools
The model tool employed is a multicomponent multiphase flow and solute transport software TOUGHREACT, which introduces reactive chemistry into the multiphase fluid and heat flow code TOUGH. The numerical method for fluid flow and chemical transport simulations is based on the integral finite difference (IFD) method for space discretization. Thermal, hydraulic and chemical fields are coupled to model fluid flow, solute transport and reactions simultaneously [40].
At present, TOUGHREACT has been widely used in diagenesis-related research, enhanced geothermal system, CO 2 storage and other geological study. It can accommodate any number of chemical species present in the liquid, gas and solid phases. Several porosity-permeability models are available to calculate the reservoir porosity and permeability changed by mineral transformation in real time. Mineral dissolution and precipitation are simulated by equilibrium and kinetic equations, which are showed as follows [40].

Equilibrium Mineral Dissolution/Precipitation
The mineral saturation ratio can be expressed as where m is the equilibrium mineral index, K m is the corresponding equilibrium constant, c is concentration, γ is thermodynamic activity coefficient, and Nc is aqueous species. At equilibrium, we have SI m = log 10 Ω m = 0 where SI m is called the mineral saturation index.

Kinetic Mineral Dissolution/Precipitation
Kinetic rates could be functions of non-basis species as well. Usually the species appearing in rate laws happen to be basis species.
r n = f c 1 , c 2 , · · · , c N C = ±k n A n 1 − Ω θ n η n = 1, . . . , N q where positive values of r n indicate dissolution, and negative values precipitation, k n is the rate constant (moles per unit mineral surface area and unit time), which is temperature dependent, A n is the specific reactive surface area per kg H 2 O and n is the kinetic mineral saturation ratio.

One-Dimensional Flow
In order to clearly distinguish and understand the influence of any single factor on mineral transformation and porosity evolution, and to avoid multi-factor interference, a simple and typical one-dimensional model was established. The total length of this horizontal direction model was 150 m, divided into 50 equidistant grid cells. Fluid was injected from the left side at a constant rate, which promoted the flow of formation water from the left to the right.
In the Base Case, parameters were set according to the actual physical properties and chemical conditions of the carbonate reservoir studied. The initial temperature was set to 40 • C and the initial pressure was set to equal to the atmospheric pressure. The initial state of the reservoir was assumed to be homogeneous, and the initial porosity was 0.35. At the beginning of the model run, the minerals in the formation were all calcite. The flow rate of fluid in the formation was set to 4 m per year, while the pH of the inflowing external water was 8.5, and the Mg/Ca ratio was 5.25. Thirteen different models, Cases 1-13 were set up to compare the effects of different factors on dolomitization, such as temperature, flow rate, seawater concentration, Mg/Ca ratio, pH and SO 4 2− concentration. The specific conditions of these various experiments are listed in Table 1. The values not specified in Table 1 are shown with a dash (-) and indicate same conditions as the Base Case. Table 2 displays four seawater concentrations involved in the simulation.

Vertical Profile Flow
Shallow areas, close to the surface were not water-saturated, which inevitably led to unsaturated areas to exit. As different groundwater levels must affect the hydrodynamic conditions in the reservoir, two different groundwater levels (50 m and 100 m) were set to compare their effects on the dolomitization. Except for the different groundwater levels, the simulation conditions in both models were the same.
The initial mineral in the model was calcite, and the initial water type was set as shallow formation water, with a low concentration of K + , Ca 2+ , Na + , Mg 2+ , Cl − , SO 4 2− , HCO 3 − , etc.
During the simulation, external fluids were continuously injected into the reservoir from the upper side of the formation. The external fluid had the same properties as seawater with higher concentrations of Ca 2+ and Mg 2+ .

Diagenesis Evolution
Based on the six successive diagenetic stages categorized in Section 2, six corresponding sub-models were established. The parameters of temperature, pressure and fluid for each diagenetic stage were taken from measured data of the reservoirs studied (Table 3).
In four different sedimentary environments models, the parameter settings were as follows: the simulated temperature and pressure was 25 • C and 0.1 MPa respectively. The seawater component had a salinity of 35 ppt [41]. The initial mineral was calcite, and the mineral content of each sub-model was based on real-time data from the previous stage. 1 Marine phreatic environment in the subtidal zone. Under the influence of tides, seawater leaks from the surface. The infiltration rate corresponds to the water exchange rate between seawater and formation.
Minerals 2020, 10, 574 7 of 18 2 Meteoric fresh water and freshwater phreatic environment in the supralittoral zone. The diagenetic process consists of two sub-processes: (1) the atmospheric freshwater leaching process, with the following parameters: infiltration rate referenced to the infiltration rate of rainfall in equatorial regions, fluid composition corresponding to the equatorial region rainwater composition. (2) The shallow layer water flow process, with the fluid defined as mixed atmospheric fresh water and formation water. 3 Seawater evaporation environment in shoal. The infiltration rate corresponds to that at the surface. 4 Seawater and freshwater interaction environment in intertidal zones: The diagenetic process includes two sub-processes: (1) atmospheric freshwater leaching with an infiltration rate corresponding to the annual rainfall in the equatorial region and a fluid composition referring to the rainwater component in the equatorial region. (2) Seawater infiltration process after sea level rise. The infiltration rate corresponds to the water exchange rate between seawater and formation.

Results
The main chemical reactions occurring in carbonate reservoirs are shown in reaction formulas (4)- (6). With the intrusion of external fluid, the original balance was disrupted, which made the calcite dissolve and release Ca 2+ and CO 3 2− , which could then combine with Mg 2+ in the external fluid to generate dolomite. The released Ca 2+ could combine with SO 4 2− to generate gypsum under certain conditions in this process.

One-Dimensional Flow
Relative content of calcite, dolomite and porosity evolution in the Base Case are shown in Figure 4.  In Case 1, the concentration of SO4 2− in the injected water was 10 times higher than that in the Base Case. Results in Figure 5 show that the calcite nearest to the injection point dissolved first and gradually transformed to dolomite. No calcite was found within 70 m of the injection point after 1 My, when the dolomite content was up to 0.55. During the transformation of calcite to dolomite, the generated Ca 2+ moved continuously to the right area as the water flows. When Ca 2+ and SO4 2− concentrations in the system reached the conditions for gypsum precipitation, gypsum was formed.   In Case 1, the concentration of SO 4 2− in the injected water was 10 times higher than that in the Base Case. Results in Figure 5 show  In Case 1, the concentration of SO4 2− in the injected water was 10 times higher than that in the Base Case. Results in Figure 5 show that the calcite nearest to the injection point dissolved first and gradually transformed to dolomite. No calcite was found within 70 m of the injection point after 1 My, when the dolomite content was up to 0.55. During the transformation of calcite to dolomite, the generated Ca 2+ moved continuously to the right area as the water flows. When Ca 2+ and SO4 2− concentrations in the system reached the conditions for gypsum precipitation, gypsum was formed.   within 90 m and 75 m in Case 5 and Case 6, respectively.
The difference between Case 7, Case 8 and the Base Case was the Mg/Ca ratio in the injected fluid. The Mg/Ca ratio in the Base Case was 5.25, while that in Case 7 and Case 8 were 10 times and 2 times that of the Base Case. The porosity in the whole simulation range increased to more than 0.4 after 0.5 My in Case 7 while the range was only about 50 m in Case 8.
Case 9 and Case 10 differed from the Base Case in pH values of the injected fluid, which was 8.5 in the Base Case, against 6.5 and 9 for Case 9 and Case 10, respectively. Within 20 m of the injection point, the porosity in Case 9 was 44% at the injection point, whereas 41% in Case 10.
Case 11, Case 12 and Case 13 differed from the Base Case with regards to the reservoir temperature, which was 40 °C for the Base Case, against 60 °C, 80 °C and 100 °C, for Cases 11, 12 and 13, respectively. These porosity evolution curves were largely different.

Vertical Profile Flow
The calcite and dolomite content at three key times for groundwater levels of 100 m and 50 m are shown in Figure 7, respectively. Regardless of the level of the groundwater, the external fluid first infiltrated into the saturated zone, after which it infiltrated into the reservoir from the surface, initiating a reaction. The reaction was enhanced at the interface between the unsaturated and the saturated zone. In the unsaturated zone, the reaction was gradually developed in the fluid flow path as a result of the accumulation of ions. For a groundwater level of 100 m, the reaction extended to more than 200 m in the horizontal direction, while for a groundwater level of 50 m, the reaction expanded to 100 m in the horizontal direction after 0.2 My. The difference between Case 4, Case 5, Case 6 and the Base Case was the injected seawater. In the Base Case, Seawater 1# was injected, while Seawater 2#, 3# and 4# were injected to Case 4, Case 5 and Case 6, respectively ( Table 2). The salinities of the four types of seawater were 27 ppt, 15 ppt, 32 ppt and 42 ppt, respectively. The porosity in Case 4 was seen to increase significantly and reached a maximum of 46%. The porosity higher than 0.44 at 1 My was within 30 m in Case 4, while within 90 m and 75 m in Case 5 and Case 6, respectively.
The difference between Case 7, Case 8 and the Base Case was the Mg/Ca ratio in the injected fluid. The Mg/Ca ratio in the Base Case was 5.25, while that in Case 7 and Case 8 were 10 times and 2 times that of the Base Case. The porosity in the whole simulation range increased to more than 0.4 after 0.5 My in Case 7 while the range was only about 50 m in Case 8. Case 9 and Case 10 differed from the Base Case in pH values of the injected fluid, which was 8.5 in the Base Case, against 6.5 and 9 for Case 9 and Case 10, respectively. Within 20 m of the injection point, the porosity in Case 9 was 44% at the injection point, whereas 41% in Case 10.
Case 11, Case 12 and Case 13 differed from the Base Case with regards to the reservoir temperature, which was 40 • C for the Base Case, against 60 • C, 80 • C and 100 • C, for Cases 11, 12 and 13, respectively. These porosity evolution curves were largely different.

Vertical Profile Flow
The calcite and dolomite content at three key times for groundwater levels of 100 m and 50 m are shown in Figure 7, respectively. Regardless of the level of the groundwater, the external fluid first infiltrated into the saturated zone, after which it infiltrated into the reservoir from the surface, initiating a reaction. The reaction was enhanced at the interface between the unsaturated and the saturated zone. In the unsaturated zone, the reaction was gradually developed in the fluid flow path as a result of the accumulation of ions. For a groundwater level of 100 m, the reaction extended to more than 200 m in the horizontal direction, while for a groundwater level of 50 m, the reaction expanded to 100 m in the horizontal direction after 0.2 My.

Vertical Profile Flow
The calcite and dolomite content at three key times for groundwater levels of 100 m and 50 m are shown in Figure 7, respectively. Regardless of the level of the groundwater, the external fluid first infiltrated into the saturated zone, after which it infiltrated into the reservoir from the surface, initiating a reaction. The reaction was enhanced at the interface between the unsaturated and the saturated zone. In the unsaturated zone, the reaction was gradually developed in the fluid flow path as a result of the accumulation of ions. For a groundwater level of 100 m, the reaction extended to more than 200 m in the horizontal direction, while for a groundwater level of 50 m, the reaction expanded to 100 m in the horizontal direction after 0.2 My.

Marine Phreatic Environment in the Subtidal Zone
Porosity evolution under marine phreatic environment is shown in Figure 8. Strong seawater cementation resulted in a calcite cement content of 5.1% during the sedimentary-parasyngenetic stage. During the parasyngenetic-shallow burial stage, seawater with high concentration of Ca 2+ and Mg 2+ was filled, which led to a strong cementation and to a drop of porosity to 20%. Dissolution occurred with the development of numerous secondary pores during the supergene stage. During the shallow burial stage, a large amount of dolomite and calcite precipitated, which decreased porosity to 18%. As a considerable amount of calcite was replaced by dolomite, the porosity decreased to 4.8% during the middle-deep burial stage. Finally, the porosity decreased to 1.6% because of calcite cementation in the deep burial stage.

Marine Phreatic Environment in the Subtidal Zone
Porosity evolution under marine phreatic environment is shown in Figure 8. Strong seawater cementation resulted in a calcite cement content of 5.1% during the sedimentary-parasyngenetic stage. During the parasyngenetic-shallow burial stage, seawater with high concentration of Ca 2+ and Mg 2+ was filled, which led to a strong cementation and to a drop of porosity to 20%. Dissolution occurred with the development of numerous secondary pores during the supergene stage. During the shallow burial stage, a large amount of dolomite and calcite precipitated, which decreased porosity to 18%. As a considerable amount of calcite was replaced by dolomite, the porosity decreased to 4.8% during the middle-deep burial stage. Finally, the porosity decreased to 1.6% because of calcite cementation in the deep burial stage.

Meteoric Fresh Water and Freshwater Phreatic Environment in the Supralittoral Zone
Porosity evolution under meteoric fresh water and freshwater phreatic environment is shown in Figure 9. With the synergistic effect of fresh water leaching, freshwater cementation and mechanical compaction, the porosity decreased to 38% during the sedimentary-parasyngenetic stage. As sea levels rose, the overlying seawater seeped into the stratum, and the porosity decreased to 22.2% in the parasyngenetic-shallow burial stage. During the supergene stage, dissolution occurred in the freshwater leaching environment, and the porosity increased to 33.9%, while during the burial stage, calcite cementation and dolomitization made the porosity decrease to 14.5%. As high salinity acidic hydrothermal fluid acted as diagenetic fluid, contributing to cementation, dissolution and dolomitization, porosity was almost constant throughout the middle-deep buried stage. Poor fluidity of formation water led to the weak cementation of calcite and the porosity reduced to 4.2% during the deep burial stage.

Meteoric Fresh Water and Freshwater Phreatic Environment in the Supralittoral Zone
Porosity evolution under meteoric fresh water and freshwater phreatic environment is shown in Figure 9. With the synergistic effect of fresh water leaching, freshwater cementation and mechanical compaction, the porosity decreased to 38% during the sedimentary-parasyngenetic stage. As sea levels rose, the overlying seawater seeped into the stratum, and the porosity decreased to 22.2% in the parasyngenetic-shallow burial stage. During the supergene stage, dissolution occurred in the freshwater leaching environment, and the porosity increased to 33.9%, while during the burial stage, calcite cementation and dolomitization made the porosity decrease to 14.5%. As high salinity acidic hydrothermal fluid acted as diagenetic fluid, contributing to cementation, dissolution and dolomitization, porosity was almost constant throughout the middle-deep buried stage. Poor fluidity of formation water led to the weak cementation of calcite and the porosity reduced to 4.2% during the deep burial stage.

Seawater Evaporation Environment in Shoal
Porosity evolution under seawater evaporation environment is shown in Figure 10. During the sedimentary-parasyngenetic stage, cementation, dolomitization and mechanical compaction were the main reactions and porosity decreased to 36%. In the parasyngenetic-shallow burial stage, poor fluidity of the infiltration fluid caused a large amount of calcite cement, resulting in a decrease in porosity. The porosity increased to 30.5% during the supergene stage. Due to calcite cementation and dolomitization, porosity decreased to 18.2% during the shallow burial stage. During the middle-deep buried stage, a large amount of calcite was converted into dolomite and a small amount of siliceous cement was formed, resulting in porosity decrease. Lastly, porosity decreased to 10.8% during the deep burial stage.

Seawater Evaporation Environment in Shoal
Porosity evolution under seawater evaporation environment is shown in Figure 10. During the sedimentary-parasyngenetic stage, cementation, dolomitization and mechanical compaction were the main reactions and porosity decreased to 36%. In the parasyngenetic-shallow burial stage, poor fluidity of the infiltration fluid caused a large amount of calcite cement, resulting in a decrease in porosity. The porosity increased to 30.5% during the supergene stage. Due to calcite cementation and dolomitization, porosity decreased to 18.2% during the shallow burial stage. During the middle-deep buried stage, a large amount of calcite was converted into dolomite and a small amount of siliceous cement was formed, resulting in porosity decrease. Lastly, porosity decreased to 10.8% during the deep burial stage.

Seawater Evaporation Environment in Shoal
Porosity evolution under seawater evaporation environment is shown in Figure 10. During the sedimentary-parasyngenetic stage, cementation, dolomitization and mechanical compaction were the main reactions and porosity decreased to 36%. In the parasyngenetic-shallow burial stage, poor fluidity of the infiltration fluid caused a large amount of calcite cement, resulting in a decrease in porosity. The porosity increased to 30.5% during the supergene stage. Due to calcite cementation and dolomitization, porosity decreased to 18.2% during the shallow burial stage. During the middle-deep buried stage, a large amount of calcite was converted into dolomite and a small amount of siliceous cement was formed, resulting in porosity decrease. Lastly, porosity decreased to 10.8% during the deep burial stage. Porosity evolution under seawater and freshwater interaction environment is shown in Figure  11. Dissolution, cementation and mechanical compaction were the main diagenetic processes

Seawater and Freshwater Interaction Environment in Intertidal Zones
Porosity evolution under seawater and freshwater interaction environment is shown in Figure 11. Dissolution, cementation and mechanical compaction were the main diagenetic processes during the sedimentary-parasyngenetic stage, and the porosity decreased to 35%. The porosity decreased to 19.7% because of calcite cement during the shallow burial stage and increased to 29% under the action of atmospheric freshwater leaching during the surface stage. In the shallow burial stage, the pore water gradually evolved into dolomitized fluid, causing calcite cementation and dolomitization and the porosity decreased to 13.9%. During the middle-deep buried stage, acidic hydrothermal fluid with medium-low temperature, high salinity and rich SiO 2 penetrated along the fault zone, and the porosity decreased to 9.3%. However, the porosity increased to 19.0% after hydrothermal alteration. The poor fluidity of the formation water and low porosity-permeability resulted in a weak cementation and a porosity decrease of 17.1% during the deep burial stage.

Seawater and Freshwater Interaction Environment in Intertidal Zones
Porosity evolution under seawater and freshwater interaction environment is shown in Figure  11. Dissolution, cementation and mechanical compaction were the main diagenetic processes during the sedimentary-parasyngenetic stage, and the porosity decreased to 35%. The porosity decreased to 19.7% because of calcite cement during the shallow burial stage and increased to 29% under the action of atmospheric freshwater leaching during the surface stage. In the shallow burial stage, the pore water gradually evolved into dolomitized fluid, causing calcite cementation and dolomitization and the porosity decreased to 13.9%. During the middle-deep buried stage, acidic hydrothermal fluid with medium-low temperature, high salinity and rich SiO2 penetrated along the fault zone, and the porosity decreased to 9.3%. However, the porosity increased to 19.0% after hydrothermal alteration. The poor fluidity of the formation water and low porosity-permeability resulted in a weak cementation and a porosity decrease of 17.1% during the deep burial stage.

Effect of Geological Factors on Dolomitization
The effect of the temperature, flow rate, fluid properties and other geological factors on dolomitization varied [23,42,43]. Results of one-dimensional models indicated that if the flow rate of the external fluid, i.e., the hydrodynamic condition, was different, the degree of dolomitization varied. The porosity change extent was consistent under various flow rate conditions, but different amounts of time were needed for the porosity to reach the same level. In short, the better the hydrodynamic conditions, the faster the fluid migration, which then pushed the reaction and expands the scope of dolomitization.
In two-dimensional models, comparing different groundwater surfaces, it was found that when the groundwater level was low, a 100 m unsaturated zone exist in the shallow layer, which contained much more gas. When external fluid entered the formation, it took a longer path to reach saturation, while more pores need to be filled in the meantime. Therefore, the reaction occurred later, and the degree of reaction was also weak due to the large number of ions being diluted during the porous flow process.
Different solution properties also lead to various levels of dolomitization [13,43]. The higher the Mg/Ca ratio in the fluid, the better the precipitation of dolomite, due to the large amount of Mg needed in the dolomitization process. However, the pH value of the fluid only had a limited influence compared with other factors.
Dolomitization is a reaction that involves a volume decrease, i.e., the volume of minerals is reduced with the conversion of calcite to dolomite, so that the number of pores within the reservoir increases [13,23]. However, if a large amount of gypsum that can occupy some pores is formed in this process, then the porosity will decrease. The results of this study demonstrated that the amount and the degree of precipitation gypsum were various under different diagenetic temperatures, which led to a large disparity in reservoir porosity.

Comparison of Model and Test Results during Successive Diagenetic Stages
Diagenetic models' results were consistent with test data both in mineral transformation and porosity evolution. In models, calcite and dolomite contents increased and the porosity decreased from 45% to about 35% throughout the sedimentary-parasyngenetic stage (488-465 Ma), meanwhile, an observation also shows that the main diagenesis was carbonate cementation. It is worth noting that the porosity evolution varied under different sedimentary environments, which was determined by the complex fluid composition. The depth ranges in the parasyngenetic-shallow burial stage (465-460 Ma) was 50-600 m. As sea levels rose, the overlying seawater seeped into the formation and replenished Ca 2+ and HCO 3 − in the pore water. Both test and model results indicate that increased temperature and pressure with the burial depth was beneficial to calcite cementation. During the supergene stage (460-455 Ma), deeply buried rocks were uplifted below the diving surface due to tectonic movements and underwent short-term erosion and weathering. With the leaching and dissolution of meteoric fresh water, fresh water with low ion concentration infiltrated into the formation, thus dissolving calcite, which led to the development of a large number of secondary pores.
During the shallow burial stage, porosity and mineral contents evolution show that the main diagenetic reactions were cementation and dolomitization, associated with a porosity decrease. As the burial depth increased, a gradual increase of formation temperature favored the overcoming of the dynamic obstacles to dolomite formation, which were dominant at the beginning of the shallow burial stage [44]. High-salinity concentrated seawater is the diagenetic fluid at this stage, which can easily flow downward under the effect of gravity and cause convection with low-density seawater at the bottom. Following the penetration of the high-salinity seawater in the stratum, calcite would begin to cement, causing the Ca 2+ concentration in the fluid to decrease, and the Mg 2+ /Ca 2+ to increase, which is beneficial for the fluid to overcome the kinetic obstacles. Calcite can then be replaced by dolomite, hence resulting in the decrease in calcite content with burial time. Previous studies have also shown that massive dolomitization of sediments is the main diagenesis in the shallow burial stage [45,46].
During the middle-deep burial stage, the diagenetic fluids were medium-high salinity, SiO 2 -rich and CO 2 -rich deep hydrothermal fluids. SiO 2 had a relatively high saturation in the hot liquid, whose temperature was about 200 • C after entering the formation (the formation temperature was about 100 • C). As the fluid flowed upward, reactions became different because of the different crack and fracture system in the formation. Strong hydrothermal dissolution happened in the fault zone of SN4 and SN5, with a developed fracture system.
During the deep burial stage (252-0 Ma), the temperature and pressure were high, which was conducive to the cementation and dolomitization. Calcite cementation occurred when high Ca 2+ fluid entered the formation from the Lower Cambrian source rock, leading the content of calcite to increase and the porosity to decrease. With the high concentration of Ca 2+ , the Mg/Ca continuously diminished, which is unfavorable for the dolomite precipitation.

Implications for Reservoir Evolution under Various Sedimentary Environments
The carbonate rock in the studied area experienced complex diagenesis under different sedimentary environments [43]. The various diagenetic phenomena observed were in the present results from the comprehensive effects of the sedimentary environments [47].
In the marine phreatic environment, cementation and compaction were the main diagenesis. As the carbonate deposits were all located below the sea level, the diagenetic fluid was seawater with a high concentration of Ca 2+ and Mg 2+ , which favored the formation of calcium carbonate but not that of dolomite. A large number of thick ring-like cements and fibrous calcite cements could be seen in the thin slices of the SN6, which is a typical single well representative of sedimentary diagenesis in the seawater subsurface flow environment.
In the meteoric water and freshwater phreatic environment, dissolution, cementation and compaction were the main diagenesis. With the sea level decline, some carbonate sediments became exposed. Atmospheric fresh water leach causing calcite dissolution porosity increased. With continuous infiltration, the saturation increased and calcite cement was formed, leading the porosity to decrease. A large number of fenestrae secondary pores and few equiaxed-granular freshwater cements filling the slab hole could be seen in the thin slices of the SN7, which is a typical single well of sedimentary diagenesis in the meteoric fresh water and freshwater phreatic environment.
In the seawater evaporative environment, dolomitization and compaction were the main diagenesis. Under dry climatic conditions, strong evaporation caused seawater to concentrate. High salinity seawater flowed into the formation, resulting in an increase in the concentration of Mg 2+ and Ca 2+ in the pore water. This hypersaline water with a high concentration of Mg 2+ increased the probability of effective collision between Mg 2+ and CO 3 2− , increasing the dolomite content. SN3 and SN5 are typical wells of the sedimentary diagenesis in the seawater evaporation environment.
In the seawater and freshwater interaction environment, cementation, dissolution and compaction were the main diagenesis. Controlled by secondary sedimentary cycles and the sea level decline temporarily, carbonate deposits in the intertidal zone were exposed to the atmosphere. Afterwards, as the sea level rose, carbonate deposits were submerged again by seawater. A large number of intragranular and intergranular dissolved pores and casting pores were formed in SN4, and the pores were filled with equiaxed-granular freshwater cements and fibrous calcite seawater cements formed during the parasyngenetic stage.

Conclusions
A series of test analysis and numerical simulation were combined to study the effects of various factors (such as temperature, flow rate, seawater concentration, Mg/Ca ratio, pH and SO 4 s concentration) on dolomitization during diagenesis in carbonate reservoirs. The degree of dolomitization varied with the flow rate and other hydrodynamic conditions of the external fluid. The better the hydrodynamic conditions, the faster the fluid migration, which then pushed the reaction and expanded the scope of dolomitization. Different solution properties and minerals also led to various levels of dolomitization. During successive diagenetic stages under different sedimentary environments with various temperatures, minerals, solutions and other conditions, reservoir experienced complicated fluid-rock reactions. The diagenetic process and porosity evolution curves in four different sedimentary environments were re-established in the studied area. The main controlling factors of dolomitization were identified, which is helpful to clarify the genesis and distribution of carbonate reservoirs. As the fluid-rock interaction mechanism in carbonate reservoirs is a complex process, further experimental and numerical simulation studies are planned to elaborate and enhance our understanding.