Alteration effects of karstification and hydrothermalism on middle Permian Qixia formation at the Wulong section, South China

Middle Permian Qixia Formation in the southwestern region of Sichuan (SW China) has experienced multiphase fluidisation, resulting in an unclear understanding of the reservoir reconstruction effect. In this study, a systematic analysis of the Qi2 member in Wulong Town was carried out by combining field outcrop petrology and geochemistry. The results demonstrated that multiple sets of crystalline dolomite-bioclastic limestone cycles were stacked vertically in the Qi2 member, accompanied by the development of fractures and karst channels. The dolomite was mainly composed of silty-fine dolomite (D1) and recrystallised dolomite (D2). Furthermore, obvious multiphase dolomitic cements (Cd1-Cd2) were present in the fractures and pores. Early karst is known to have lithologic mutation surface development and karst channel development at the top of several secondary cycles. The vadose silt dolomites (Cd1) having karst channels developed dull luminescence under cathode luminescence (CL). Both the geochemical indicators of elements and rare earth element (REE) content indicated dysoxic-oxic environmental conditions. The hydrothermal solution displayed tectonic carniole characteristics in the strata burial stage. Fractures and pores were filled with hydrothermal minerals such as coarse dolomites-saddle dolomites (Cd2, with some caused by recrystallisation of the Cd1 hydrothermal solution) and fluorites. Coarse dolomites-saddle dolomites developed dull-red luminescence with a bright-red rim under CL and their δ18OVPDB values were more negative than those of middle Permian limestone samples. Both the geochemical indicators of elements and REE content indicated the suboxic-anoxic environmental conditions. Karstification had minor constructive impact on the reservoir of the Qi2 member in Baoxing in southwestern Sichuan. Most products of karstification were distributed as fillings in channels. Aside from creating certain networked fractures, the hydrothermal solution was mainly filled with hydrothermal minerals along the fractures, pores and early karst channels. Karst and the hydrothermal solution mainly damaged the middle and upper parts of the middle Permian Qixia Formation in Southwest Sichuan. The impact of episodic fluid on the restoration of the carbonate reservoir was mainly restricted by channels for fluid migration and thickness differences among the reservoir. However, certain thick-layered and massive crystalline dolomite may hold promise for exploration.

Tectonic and sea level changes may cause extensive fluid migration 1 and karstification in many sedimentary basins.This fluid can be composed of evaporated brine, seawater and hydrothermal solutions [2][3][4] , which may induce the dolomitisation of carbonate formations under specific physical and chemical conditions [5][6][7] .In the field of petroleum geology, considerable attention has been paid to the dolomite reservoirs produced by hydrothermal solutions, which are a common type of dolomitisation fluid 8,9 .Additionally, carbonatite dissolution (karstification) by surface unsaturated fluids may also produce karst-type reservoirs that are of significance to the development of carbonate reservoirs [10][11][12] .Karst reservoirs account for 20-30% of explored carbonatite petroleum and natural gas reservoirs 13 .They are extensively distributed globally, including South America 14 , Australia 15 and the Sichuan Basin in China 16 .However, certain carbonatites went through restoration by the simultaneous karstification of unsaturated fluid and alteration of middle-deep buried fluids 17,18 .The collective impact of karstification and alteration on reservoirs still requires further investigation.
The middle Permian dolomite reservoir in Sichuan Basin (SW, China) is an important natural gas exploration series.It has attracted considerable research attention because of its special structural framework, complicated sedimentary environment and reservoir-forming conditions 19,20 .The formation of middle Permian dolomite reservoirs is mainly related to dolomitisation.However, the reservoir performance is constrained by various factors such as the sedimentary phase, karstification and tectonism.The degree to which the bank facies sedimentation superimposed on the eogenetic karst is preserved during the burial process remains ambiguous.Moreover, the controversy persists regarding whether organic acids or thermal fluids are responsible for largescale reservoir reconstruction in areas with complex tectonic superposition and medium-deep burial conditions.It has been found that strata having karst systems have good potential for exploration and development 21,22 .However, it remains uncertain whether the modification effect is beneficial to reservoir development during the combined influence of late multi-phase tectonic/fluidic events on local strata.
Following the hydrothermal dolomitisation and initial in-depth studies, numerous studies have been conducted on reservoir modification by hydrothermal fluids [23][24][25] .However, different opinions could be found on the combined effects of thermal fluid on carbonate transformation.It has been reported that the rising thermal fluid controls the formation of high-porosity and high-permeability belts in fold or fault carbonate reservoirs via the strike-slip faults as channels 26,27 .Thermal fluid could block existing pores during migration, since CO 2 -rich thermal fluid may corrode surrounding rocks to form pores and precipitate saddle dolomite under specific conditions 28 .Hydrothermalism is a common phenomenon that has been observed in carbonate reservoirs in areas such as the Sichuan Basin 29 and Bohai Basin 30 .The fluorine-containing thermal fluid can be a major factor that controls carbonate reservoir transformation 31 and has established the paleokarst reservoir transformation pattern based on thermal fluid.This pattern summarises the process in which karstification provides fundamental pore frameworks for carbonate reservoir and thermal fluid reconstruction of the built reservoir space in the middledeep burial stage 32 .In general, the debate is on the degree of dissolution-filling of pores and the improvement of the percolation capacity after the reconstruction of early pores by late fluids.
Middle Permian carbonate reservoirs in the southwestern region of Sichuan are commonly reformed by complex diagenetic fluids.What are the action mechanisms of such diagenetic fluids?What are their transformation effects?What are the characteristics of their spatial and temporal coupling?Therefore, a series of studies were carried out on the rock-geochemical characteristics in the middle Permian Qixia Formation (Qi2 member) of Baoxing and the southwestern region of Sichuan.The objective of this study was to (1) determine the lithological and geochemical characteristics of the Qixia Formation; (2) implement the diagenetic evolutionary process and (3) understand the combined impact of karst-thermal fluids on carbonate formation.The findings of this research could support oil and gas exploration in middle Permian carbonate reservoirs in the Sichuan Basin.

Geological setting
The middle Permian Qixia Formation of the Wulong Section is situated in proximity to Wulong County in the southwestern region of the Sichuan Basin in southern China.From a tectonic perspective, it is very close to the Yingxiu-Wulong Fault (Figs. 1a,b and 2b) and is accompanied by the extensive occurrence of faults 33 .In terms of palaeogeography, the middle Permian Qixia Formation is a component of the north-western region of the Yangtze Platform and is close to the western edge of the Sichuan Basin.The Yangtze Platform was flooded due to transgression during middle Permian, resulting in its transformation into a shallow platform (Fig. 1c,d).According to previous studies, the overall water depth in the Yangtze Platform is about 50 m and could reach 200 m in certain deep-water areas 34 .The Longmenshan Fault Zone has undergone two major stages of tectonic evolution since the Longmen Mountains at the western edge of the Sichuan Basin formed a uniform base due to the influence of the Jinning Movement and Chengjiang Movement in the Proterozoic era.These stages are the tensional stage, from the Sinian to the Middle Triassic, and the extrusion stage, since the Late Triassic 35,36 .There has been strong rifting activity from the Late Silurian to the Middle Triassic, which reached its highest point during the Late Permian and was accompanied by large-scale magmation.This strong Permian rifting activity was named the 'Emei Taphrogenesis' 37 .The Yingxiu-Wulong Fault is a component of the Longmenshan Central Fault.Subsequently, the Indosinian Movement, Yanshan Movement, Himalayan Movement and other tectonic-related events occurred.The strata of Qixia Formation were possibly reformed by thermal fluid-karst due to frequent changes in sea level and tectonic-thermal events.
Middle Permian strata of the Sichuan Basin are unconformable to the underlying Carboniferous strata after exposure for a long time.They are composed of the bottom Liangshan, middle Qixia and overlying Maokou Formation.The Liangshan Formation covers the regional plane of unconformity.It mainly consists of carbonaceous shale, coal and sandstone, which reflects the coastal marsh environment.The Qixia Formation consists of a complete tertiary cycle of transgression and regression, which can be divided into two members.The Qi1 member is mainly composed of limestone that is rich in organic matter with flint lumps at the bottom.It is classified as part of the transgression system tract.The Qi2 member is mainly composed of light grey-off-white crystalline dolomite with light grey micrite-sparry bioclastic limestone in local areas.It is part of the highstand systems tract and generally belongs to the open platform sedimentary environment.In this study, the main focus was on the upper parts (~ 20 m thick) of the Qi2 member, which was the key section for oil and gas exploration in the region.The telogenetic karst of the Qixia Formation underwent partial alteration after sedimentation, following which the sea level rose sharply and the Maokou Formation was developed.The Maokou Formation mainly had thin layers of dark grey-cinereus argillaceous limestone containing organic matter, with the common development of a nodular structure (Fig. 2a).It belongs to the low-energy deep-water sedimentary environment.
In the strata sequence, the vertical lithofacies variation of the ~ 40 m thick section in the Qi2 member showed clear cyclicity (Fig. 2a).The basic stratigraphic unit (cycle) was composed of two parts: the lower bioclastic limestone and the upper crystalline dolomite.This formed a meter-level upward shallow cycle.Such stratigraphic

Samples and methods
A total of 32 rock samples were collected successively from the Qi2 member of the Wulong Section for sedimentological and geochemical studies.Alizarin red stain was prepared to distinguish calcite and dolomite.One hundred thin sections from 32 rock samples were processed by double-side polishing and immersion in blue epoxy resin to determine their lithological composition, porosity and diagenetic sequence.Images of the processed thin sections were captured under an optical microscope (Olympus BX53M) and cathode luminescence microscope (CL8200MK5) at Southwest Petroleum University (Figs. 4 and 5).Trace elements and rare earth elements (14 samples) were analysed by Analytical Chemistry and Testing Services (ALS) in Guangzhou.Quantitative analyses were performed using lithium borate fusion and inductively coupled plasma mass spectrometry.Carbon and oxygen isotope (14 samples) tests were performed following the microdrill sampling at Southwest Petroleum University.These samples were prepared by the phosphate method and analysed using an IsoPrimeGC5 isotope mass spectrometer (Elementar, Germany).Stable isotope data were converted to permil (‰) relative to Vienna Pee Dee Belemnite (V-PDB) and corrected by fractionation factors mentioned by Fairchild and Spiro 39 .The precision of the δ 18 O and δ 13 C ratio data was better than ± 0.1‰.The rare earth element (REE) data were normalised to post-Archean Australian shale (PAAS) 40 .The REE anomalies were calculated by the classical equation 41 : δCe = Ce N / (Pr N 2 /Nd N ), δEu = 2Eu N /(Sm N + Gd N ).

Petrography
In the Wulong Section close to the Yingxiu-Wulong Fault, the upper outcrop and the thin sections of the Qi2 member samples showed different lithological and structural features of carbonate rocks.The matrix limestone was found to be mainly grey-light grey thin-middle bioclastic limestone (L1).The thin-middle limestone and dolomite gradually came into contact along the vertical direction.The thick-massive dolomite mainly had abrupt vertical contact with limestone (Fig. 3a).The dolomite could be divided into three types, i.e., powder to fine dolostone (D1), recrystallised dolostone (D2) and dolomite cement (Cd).Additionally, calcite cement (CC), pyrite (Py) and fluorite (Fl) were present in residual fractures and pores.The types of lithofacies and petrologic characteristics are described below.It generally consisted of thin-middle strata and had a basic cycle with the overlying thick-massive crystalline dolomite strata.Numerous bioclastic particles were found to be present.Moldic and intragranular pores were completely filled by calcite.Whereas, some of them were dolomitised (vadose silt dolomite) and cut by late fractures (Fig. 4a).

Dolomite cement (Cd).
Dolomite cement only appeared as pore-filling minerals.Cd could be further divided into Cd1 (vadose silt dolomite) and Cd2 (saddle dolomite and transitional form) according to its genesis.Cd1 was mainly composed of subhedral and euhedral crystals (0.05-0.25 mm) and was mainly formed in karst channels.Cd1 silt dolomite might have been carried by fluid media which passed through karst channels.It was not reformed or slightly reformed to a relatively small degree by late diagenetic fluids.Thus, it was speculated that Cd1 was the result of karstification (Fig. 5).Cd2 was mainly composed of anhedral and euhedral crystals

Alteration of karstification and hydrothermalism characteristics
Characteristics and spatial distribution of dolomite.Strata in the Qi2 member of the Wulong Section were about 40 m thick.In this study, the focus was given to the upper part of the Qi2 member (20 m), which could be divided into seven fourth-order cycles (C3-C9).Except for C4 and C5, each cycle was composed of the underlying thin bioclastic limestone layer and the overlying middle crystalline dolomite layer.C4 and C5 were sequenced with variable dolomite grain sizes (Fig. 3a).With changes in the sedimentary cycle, limestone and dolomite strata were formed upward in an alternating manner.Karst channels were formed at the top of most cycles (C4, C5, C7, C8 and C9) in the dolomite strata, with fractures in local areas.There were abundant residual bioclastic particles in the limestone strata, mainly including Gastropoda, Barchiopoda and Echinodermata (Fig. 4a).

Characteristics of karstification.
In contrast to telogenetic karst, the formation of eogenetic karst is based on the conduction media of karst water and different karst materials 42 .This karstification mode is driven by high-frequency fluctuations in sea levels while taking place around the penecontemporaneous stage 43 .Most contemporaneous sediments have high pore permeability and can provide channels for the diffusion of karst water.The diffusion of karst water among particles in the early stage resulted in the gradual formation of channels in local areas.Hence, double-porosity systems formed, combining channel porosity and interparticle porosity 42 .The characteristics of both telogenetic karst and eogenetic karst were found in the upper cycle of Qi2 member.Karst characteristics were mainly evident in vertical lithologic sequence characteristics, lithologic fabric characteristics and pore characteristics.The lithologic sequence contained lithologic mutation surfaces, which are generally an important characteristic of eogenetic karst 44 .Lithologic mutation surfaces were found from C5 to C6 and from C7 to C8 in field outcrops.Specifically, the underlying bioclastic limestone strata changed abruptly into crystalline dolomite strata while having obvious mutation interfaces (Fig. 2a,b,c).The cycle top belonged to a high-frequency exposed surface.Many karst channels (0.1-0.3 mm) were found in outcrop and thin sections of the sample.They were mainly distributed at the top of cycles with harbour-shaped edges and internal spaces filled with Cd1 (Fig. 5a,b,c,d).Karst channels were found to be solution-enlarged matrix fractures with interparticle pores.They had a comprehensive manifestation of eogenetic karst and telogenetic karst characteristics 45 .The thin sections of the sample had considerable moldic and intraparticle porosity.This reflected selective corrosion and the dissolution of unstable aragonite in eogenetic bioclastic particles.Under CL, Cd1 in karst channels displayed dull luminescence (Fig. 5e,f).

Characteristics of hydrothermalism.
It is generally believed that hydrothermal dolomite reservoirs are formed by the reconstruction of precursor pores and the alteration of matrix limestone by the Mg-rich fluid, which is transported from the underlying strata to the overlying limestone strata through the fault system.This process has been verified in the Michigan Basin in the northeastern United States and the Appalachian Basin in eastern Canada 7,46,47 .Therefore, hydrothermalism characteristics are mainly recognised according to the formation of hydrothermal dolomite (HTD) and the dissolution of limestone.
Hydrothermalism characteristics in a region are mainly manifested by petrofabrics and mineral combinations.Observation of the thin sections under a microscope revealed D2 to have obvious crystal brecciation (Fig. 4c).The dolomite crystals on both sides of the karst channel were broken in an interconnected pattern.The size and orientation of the crystals had controlling effect on the fragmentation of the petrofabrics and hence, resulting in directional fragmentation.This fragmentation may have been caused by hydraulic fracturing induced by high pressure generated during the hydrothermal fluid migration.The filling of pores with saddle dolomite (Cd2) was the most direct evidence of a hydrothermal environment 7 (Fig. 6a,b,c).Other typical hydrothermal-associated materials such as Py and Fl.Py mainly existed in blocks and was distributed in some fractures (Fig. 4c).Fl was distributed in fractures and among breccia while having irregular crystal morphology and triangular solution formation (Fig. 6d).Under CL, Cd2 showed dull-red luminescence with a bright-red rim (Fig. 6e,f).

Geochemistry
Trace elements.According to trace elements data, the V/Cr and Ni/Co values of L1 in the Qi2 member were 0.5 and 0.29, respectively.The V/Cr value of D1 ranged from 0.12 to 1.  1).
V/Cr and Ni/Co values are often used to evaluate the redox environment of sedimentary water bodies during the formation of minerals 48 .The V/Cr and Ni/Co values of L1, D1 and Cd1 revealed dysoxic-oxic environmental conditions while the V/Cr and Ni/Co values of D2, Cd2 and CC revealed suboxic-anoxic environmental conditions (Fig. 8a).
Ce is extremely sensitive to redox conditions.In an oxidising environment, soluble Ce 3+ would be oxidised into insoluble Ce 4+ , which could first enter other particles and then precipitate while resulting in negative Ce anomalies.Positive Ce anomalies may be present in dysoxic environments 49      Vol  1, Fig. 8b).The mean δ 13 C VPDB of D1 and Cd1 (2.01‰, n = 6) was very close to that of the mean of middle Permian marine facies carbonate in the upper Yangtze Platform (2.86‰, n = 78) 53 .δ 18 O VPDB (3.8‰ to 1.9‰) was similar to that of middle Permian dolomite 54 while having low-temperature dolomite characteristics 52 .This demonstrated that the formation of D1 and Cd1 might be related to dolomitisation from the shallow burial stage to the penecontemporaneous stage under low temperatures 55,56 .D2 and Cd2 both had relatively low δ 18 O VPDB values.The similar C-O isotopic characteristics confirmed the identical diagenetic fluid composition of D2 and Cd2 and their high-temperature dolomite characteristics 52 .Under near-surface temperature, the δ 18 O VPDB of CC precipitated in atmospheric with higher than − 8‰ of water 57 .The relatively low δ 18 O VPDB value (11.31‰) of CC indicated  www.nature.com/scientificreports/ it might be influenced by the high temperatures from the middle-deep burial stages.It was speculated that the extremely low δ 13 C value (− 7.27‰) might be due to the influence of organic acid invasion on samples during the burial process 58 .From the perspective of carbon and oxygen isotopes, there were obvious products of the penecontemporaneous-shallow burial stage and the middle-deep burial stage in the study area.

Discussion
Paragenesis.Carbonate platform sediments were formed in the middle Permian Qixia Formation in the southwestern region of the Sichuan Basin.The Qi2 member mainly belongs to the particle beach facies sediments 59 .Multi-cycle superposition sequences of thin-middle layered dark grey bioclastic limestone and thick-layered light grey crystalline dolomite were present vertically.The entire sedimentation process of the Qi2 member involved a gradual decline in the relative sea level with a short-term rapid increase in the relative sea level.Once the sedimentation was complete, the Qi2 member experienced tectonic uplift movement 60 while exposing the sedimentary strata.These factors provided ideal conditions for the development of karst.Therefore, multiple dolomite facies were developed in the Qixia Formation due to the complicated burial process.A paragenetic sequence was established as shown in Fig. 9.
The paragenetic sequence (Fig. 9) was established based on the cross-cutting relationships obtained from the field and the thin sections.L1 was the earliest original rock (Fig. 4a) while D1 recrystallised from C1 during the penecontemporaneous period (Fig. 4b).In this stage, the strata experienced Eogenetic karstification and Cd1 was transported into karst channels by fluid media.Once the sedimentation was complete, the stratum was uplifted to the surface and exposed by the geotectonic movement while causing Telogenetic karstification.Some Cd1 was filled with karst channels.Under the influence of matrix activation in the middle and late diagenetic stages 36 , hydrothermal fluid entered the strata along the faults, causing D1 to recrystallise into medium-coarse crystalline dolomite (D2).Late dolomite cement (Cd2) rapidly got filled with cracks (Fig. 4c,d).Hydrothermal minerals (e.g., Py and Fl) filled the intercrystal pores of Cd2 and its formation was slightly later than that of Cd2.Calcite cement (CC) in microfractures intersected the hydrothermal minerals and Cd2 (Fig. 4c,d).
Karstification.There were typical karstification signs at the top of the Qi2 member such as lithologic mutation surfaces caused by sea level changes 44 and karst channels with irregular morphology 45 (Fig. 2).This karst system played a significant role in creating spaces for reservoirs and providing pathways for the migration of hydrothermal fluid.However, the formation of the karst system was controlled significantly by the karstification intensity, which was closely related to the vertical sequence of cycles.Karst water at the top of the sedimentary cycles could migrate along strata positions with high initial porosity and permeability.In the early stages of karstification, highly unsaturated karst water may have diffused through interparticle pores among loose limestone while forming solution pipes and vugs for fluid conduction.With the continuous migration of karst water, its dissolving capacity decreased due to the difference in structural selection and saturation.Therefore, the karst water could no longer form new solution pipes and could only diffuse via existing solution pipes.When the karst water supply increased (for example, rainfall increased), it might have exceeded the seepage capacity of the solution pipes.As a result, karst water seeped out via interparticle pores or fractures out of solution pipes, which further expanded into existing solution pipes and vugs to form a 'karst channel-pore' system 61 .Karst channels were found to be the result of an expanded solution.
The formation of the karst system was accompanied by sea level rising and dolomitisation in limestone strata under the collective action of active and latent reflux.Discontinuous dolomite strata formed, resulting in the alternating development of dolomite and undolomitised limestone strata in the area 62,63 .Therefore, the loose limestone strata (Fig. 10a) were replaced by powder and fine dolostone (Fig. 10b) due to penecontemporaneous dolomitisation.The geochemical characteristics of samples also proved the suboxic-oxic environmental www.nature.com/scientificreports/conditions of the shallow burial stage.Additionally, Cd1 was transported with karst water and entered karst channels and filled them (Fig. 10c).
Hydrothermalism.The strong Indosinian movement reactivated the major basement fault in the study area 36 .According to previous studies, a relatively high palaeo-heat flow value is present in the southwestern region of the Sichuan Basin 64 .The hydrothermal phenomenon is a manifestation of the thermal effect of magmatic activity.Deep Mg 2+ -rich fluid influenced by the high temperature, entered the Qixia Formation along the major Indosinian Fault (base fault in the Wulong Section; Figs.1b and 3b).The small tectonic fractures accompanied the formation of the fault and the early karst system provided favourable pathways for the migration of deep hydrothermal fluids 65,66 .These tectonic fractures were mainly distributed in dolomite strata in dendritic form.Dolomitisation was not discovered in relatively compact limestone strata and only a few fractures were developed and filled with CC (Fig. 4a).After hydrothermal fluid entered the dolomite strata, D1 was transformed locally by FR-1 and recrystallised into middle-coarse crystalline dolomite (D2).Some of the hydrothermal fluid that entered the karst system might have resulted in the solution-enlarged pores that developed as karst channel and the recrystallisation of vadose silt dolomite (Cd1) into middle crystalline dolomite (Cd2, transition form of saddle dolomite).Furthermore, Cd2 was precipitated directly in some pores (Fig. 10d).After Phase-I hydrothermalism, the recrystallised dolomite strata might have undergone Phase-II hydrothermalism alteration due to the influence of many factors such as base activation.The Phase-II hydrothermal Mg 2+ -rich fluid resulted in the recrystallisation of the middle crystalline dolomite (Cd2, transition form of saddle dolomite) in the karst system into coarse crystalline dolomite (Cd2, transition form of saddle dolomite) or saddle dolomite (Cd2).When Mg 2+ in the hydrothermal fluid was consumed, residual pores in fractures would have been filled with hydrothermal minerals and CC (Fig. 10e).Once the hydrothermal activity ceased, the reconstructed strata were influenced by tectonic movements in the Yanshanian or Himalayan Periods while generating a series of unfilled tectonic fractures (Fig. 10f).
Intercrystal porosity was developed in thick-massive dolomite strata with the largest porosity of 8% (Fig. 4b).Intercrystal porosity was well preserved, which could indicate changes in the solution pipes.Under the influence of Phase-I hydrothermalism, the protolith was recrystallised (compacted with few residual intercrystal www.nature.com/scientificreports/porosities).After recrystallisation, the relatively compact crystalline dolomite strata restricted the diffusion of deep hydrothermal fluid to some extent.Consequently, it was difficult for the late hydrothermal fluid to continuously fill previous residual intercrystal porosity.However, it could only migrate along residual fractures.Hydrothermal fluids dissolved and precipitated simultaneously, forming a series of harbour-shaped corrosion rims and CC.The reservoir space in the thin dolomite strata was barely developed and the porosity was found to be 1% (Fig. 4d).This could imply that the thin-layered dolomite was more influenced by the late hydrothermal alteration than the thick-massive dolomite strata.Additionally, the loose karst system had more favourable solution pipes than the surrounding matrix rocks having broken fractures.When deep hydrothermal fluid entered the Qixia Formation strata along the basement fault, it might have first filled the karst system due to its relatively good porosity and permeability, solution channels and fluid reservoir space.Consequently, the karst system was completely filled and only a few pores were formed.Therefore, hydrothermalism in the study area was dominated by the filling effect.Intercrystal porosity in dolomite was first formed during penecontemporaneous dolomitisation.Pores were further shrunk due to hydrothermal recrystallisation of matrix dolomite while subsequently forming residual pores transformed by hydrothermal fluid.
Regional comparisons and implications.In order to establish how representative, the data from the South China study area are, it is useful to compare them to other known karstification (hydrothermalism) -controlled dolomite reservoirs.These cases are distributed in the southern Sichuan Basin of China, offshore Spain and western Canada.Alteration mechanisms and controlling factors have been described in detail.However, the relationship between these mechanisms has not been systematically established.
Similarly located in the southwest of the Sichuan Basin, southeast of this study area, four types of middle Permian dolomites have been revealed in the field (Fig. 11a,b).Petrographic and geochemistry studies have described the fabric and geochemical characteristics of the four types of dolomites as well as the dolomitisation involved in their evolution process (seepage/reflux and hydrothermal dolomitization).It indicated that the distance between the dolomitisation process and fault zones might be related to the distribution of different phases of dolomite 67 .This situation is similar to what has been that observed in the Wulong section.However, it is still not understood how the porosity evolved.
Porosity evolution has been described from the well-exposed, partially dolomitized Cambrian carbonate platform of the Mount Whyte Formation and Whirlpool Point locality of western Canada (Fig. 11g).The dolomitisation fronts and halo zones have higher porosity than their core, which is associated with the progressive decrease in porosity inboard of the reaction front (Fig. 11e,f), caused by the increase in crystal size, and decrease in pore size during recrystallization in the core of the body 68 .The porosity evolution involved in this study is only explained from the perspective of genesis while the example of western Canada illustrates why porosity increases or decreases from the perspective of mechanism, which supports the research.
Compared to the south China data and Cambrian dolomite reservoirs of western Canada, the phenomenon described by the Spanish case has a better similarity with this study.Currently, dolomitized upper Jurassic limestones constitute the main reservoir of the Casablanca oil field, located in the western Mediterranean Sea.The dolomite reservoir was also affected by the alteration effects of karstification and hydrothermalism (Fig. 11c,d).It had two distinct characteristics: (i) the reservoir rock was uplifted and subaerially exposed between the Paleogene compression and the Neogene extension while significantly increasing their porosity; (ii) the porosity was obliterated by successive cement precipitations 69 .The outcomes were similar to the phenomenon investigated in South China having multistage hydrothermal fluids resulting in multiple recrystallizations of dolomite and a decrease in porosity.
The similarity in fabric development and paragenesis of global compound alteration effects of karstification and hydrothermalism was significant.This may imply a uniform carbonate rock alteration mechanism, regardless of geological age.This study attempted to establish the relationship between complex carbonate rock alteration and reservoir porosity evolution while emphasizing porous conditions in the karst system than the surrounding matrix.However, the preservation of porous conditions under fracture networks and hydrothermalism was challenging.Furthermore, the preservation of porous conditions was also influenced by the thickness of the reservoir.The intercrystal porosity of the thick-massive reservoir was well preserved after early dolomitization and hydrothermal alteration (Figs.2a and 4b).In contrast, the thin reservoir was more susceptible to late hydrothermal alteration than the thick-massive reservoir, which resulted in a decrease in porosity, the obvious filling of fractures as well as intense compaction of the reservoir (Fig. 4a,c).Therefore, the alteration effect of karst and hydrothermal fluid was mainly destructive to carbonate in the study area.Differences in spatiotemporal distribution in local porosity were controlled by the connectivity of solution pipes and the reservoir thickness.

Conclusions
This study included comprehensive research on the petrology and geochemistry of the middle Permian Qixia Formation in the Wulong Section in the southwest region of the Sichuan Basin.The field observations, petrological and geochemical analyses revealed multiple dolomite facies.The diagenetic sequence was found to result from the alternation of karstification and hydrothermalism effects on the carbonate strata.and had the following characteristics: depleted δ 18 O VPDB values, positive Eu anomalies and trace elements of dysoxic environments.b) The fault zone and karst system created pathways for the entrance and migration of hydrothermal fluid.
The fault zone provided strata with basement faults and multiphase fractures.The karst system contained numerous irregular channels and vugs under the exposed surface, which provided pathways for the entrance of hydrothermal fluid.Hydrothermal fluid entered from the basement fault into the Qixia Formation and migrated along the favorable pathways (e.g.karst system and fracture network).Hydrothermal fluid can migrate readily through the loose early karst system along the fracture network.This migration pattern revealed the heterogeneity of reservoir alteration.c) The combined effects of karstification and hydrothermalism mainly resulted in the degradation of the dolomite reservoirs in the Qixia Formation of the Wulong Section.Although karstification provided a basic porosity framework for the dolomite reservoir, the karstification-filling effect of multiphase hydrothermalism was detrimental to the karst reservoirs.The local residual intercrystal porosity in the matrix dolomite was caused by different connectivity among pathways in the same strata sequence.The thickness of the reservoir may also have impacted the integrity of the reservoir spaces.Following the same diagenetic alteration, the thick-massive dolomite in the region could preserve reservoir spaces better than the thin-layered dolomite.

Figure 1 .
Figure 1.Geographic position of the Sichuan Basin and lithofacies palaeogeography.(a) Geographic position of the Sichuan Basin.(b) Geographic position of the study area.(c), (d) Lithofacies palaeogeography map during the Chihsian in South China and the location of the studied section (modified from Wang and Jin 34 ).

Figure 2 .
Figure 2. Field photos of the Second member, Qixia Formation, Wulong section and Southwest Sichuan Basin.(a) Division of the top Qi2 member into 7 cycles (C3-C9).The yellow inverted triangles denote regression.The red dotted lines denote fractures and the white solid lines are boundary lines of the sedimentary cycle.The white dotted line is the plane of unconformity.The yellow dotted line area is enlarged in Figure b.The red solid line area is enlarged in Figure c.The yellow solid line area is enlarged in Figure d.The white dots are sampling points.(b) The white dotted line refers to the plane of unconformity between the Qixia and Maokou Formations.The blue dotted line shows the karst channel.(c) The white solid line is the boundary line between C6 and C7, where the upper is limestone and the lower is dolomite.(d) The white solid line is the boundary line between C5 and C6, where the upper is limestone and the lower is dolomite.

Figure 3 .
Figure 3. Lithological column and geological cross-section map.(a) Lithological column of Qixia Formation at the Wulong section.(b) Geological cross-section map of the Wulong section.Section line I-I' corresponds to Fig. 1b.

Figure 7 .
Figure 7. Post-Archean average shale (PAAS)-normalised rare earth element (REE) patterns.(a) REE distribution map of limestone and karst-related dolomite, showing relevant characteristics influenced by karstification; (b) REE distribution map of CC and hydrothermal-related dolomite, showing relevant characteristics influenced by hydrothermalism.The PAAS data were obtained from McLennan 40 .

Figure 8 .
Figure 8. Stable isotope and trace element distribution patterns of sample rocks in the Second Member, Qixia Formation and Wulong section.(a) Cross-plot of V/Cr and Ni/Co values.(b) Cross-plot of δ 13 C VPDB and δ 18 O VPDB values.Palaeoredox conditions range of trace element values from Bryn and Manning 48 .Typical isotopic values of dolomite formed under low temperatures and high temperatures from Allan 52 .
a) The different dolomite facies in the middle Permian Qixia Formation in the Wulong Section in the southwest region of the Sichuan Basin represented different redox environments.According to geochemical data, D1 and Cd1 were related to karstification and had the following characteristics: enriched δ18 O VPDB values, negative Ce anomalies and trace elements of oxic environments.D2 and Cd2 were related to hydrothermalism

Bioclastic limestone (L1). Bioclastic limestone was mainly formed at the upper part of the Qi2 member.
).

Table 1 .
Stable isotope and Trace element data of sample rocks in middle Permian in SW Sichuan Basin of WL-Wulong section.

Table 2 .
Rare earth element data of sample rocks in middle Permian in SW Sichuan Basin of WL-Wulong section.