Abstract

Oil and gas have been found in the Upper Ordovician Lianglitag Formation carbonates in the Tahe Oilfield, Tarim Basin, NW China. This study documents the origin of diagenetic fluids by using a combination of petrology, SIMS, fluid inclusion, and radiogenic isotope analysis. Six stages of calcite cements were revealed. C1-C2 formed in marine to early burial environments. C3 has relatively low δ18 values (−8.45 to −6.50) and likely has a meteoric origin. Meteoric water probably fluxed into aquifers during the Early Paleozoic and Late Paleozoic uplift. C4 has δ18 values typically 3 higher than those of C3, and probably formed during shallow burial. C5 displays relatively negative δ18 values (−8.26 to −5.12), and the moderate-to-high fluid-inclusion temperatures imply that it precipitated in burial environments. C6 shows homogenization temperatures (up to 200°C) higher than the maximum burial and much lower salinities (<10.61 wt% NaCl), which may suggest that the fluid was deeply recycled meteoric water. The average 87Sr/86Sr ratios of fracture- and vug-filling calcite cements are much higher, indicative of incorporation of radiogenic Sr. Caves and fractures constitute the dominant reservoir spaces. A corresponding diagenesis-related reservoir evolution model was established that favors exploration and prediction.

1. Introduction

Fluid origin and flow in sedimentary basins have an important influence on various diagenetic alteration processes, including cementation, dissolution, recrystallization, and dolomitization, which have significant impacts on reservoir development [110]. Thus, constraining the origin and evolution of fluid and fluid flow events is of considerable significance for hydrocarbon exploration and reservoir development. Fluid is not only a component of the fluid-rock interaction but also a carrier of matter and energy in the fluid-rock interaction. Fluid is an active control factor of the fluid-rock interaction system [7, 10, 11]. However, for deeply buried reservoirs, due to their complicated burial and tectonic history, the fluid chemical composition and flow are very complicated; therefore, understanding and predicting fluid flow are still challenging.

The Tahe Oilfield, located within the Tabei tectonic block in the northern part of the basin, is the largest Paleozoic marine oilfield in China. For the Tahe Oilfield, the current main oil-gas production layer is the Middle Ordovician Yijianfang Formation and the Early Middle Ordovician Yingshan Formation, and considerable previous works have been done on these layers [1217]. By contrast, carbonate rocks of the Upper Ordovician Lianglitag Formation are potential exploration strata, which have received little attention. The carbonates of the Upper Ordovician Lianglitag Formation are deeply buried, with depths > 4500 m. Previous studies mainly focused on the biostratigraphy, lithostratigraphy [1820], and sedimentary environment [2123]. Lu et al. [24] analyzed the main reservoir characteristics of the Lianglitag Formation: the matrix of the Lianglitag Formation is characterized by low porosity and permeability, with strong reservoir heterogeneity.

In situ Secondary Ion Mass Spectrometry (SIMS) has provided a powerful technique for studying basin fluid activities [2530]. SIMS can solve the traditional test difficulties; for example, a small number of samples cannot meet the test requirements or the analysis results are often the average value of multiphase fluid mixing. In this paper, for the first time, in situ SIMS analysis has been used to characterize different fluids and discriminate their nature in study areas.

This paper takes the carbonate reservoirs of the Upper Ordovician Lianglitag Formation as a typical example and focuses on the activities of diagenetic fluids as well as their effects on reservoir development under a complex multistage tectonic setting. This work can provide useful insights into fluid evolution during diagenesis for typical deeply buried reservoirs. Specifically, we address the following questions:(1)What are the petrological and geochemical characteristics of diagenetic calcites in the Late Ordovician Lianglitag Formation?(2)What is the nature of the fluids and their impact on the Lianglitag Formation?(3)What mechanisms were responsible for the formation of the carbonate reservoir of the Lianglitag Formation in the Tahe Oilfield in the Tarim Basin?(4)What is the model for diagenetic fluid flow and reservoir performance in deeply buried carbonate successions?

2. Geological Background

The Tarim Basin, located in southern Xinjiang province, northwestern China, is surrounded by the Tian Shan, Kunlun, and Altyn mountains to the north, southwest, and southeast, respectively (Figure 1(a)). It is a multicycle composite basin with continental crust, composed of pre-Sinian basaltic basins, which are formed by prototypical basins with different properties in different periods [32, 33]. The basin can be divided into seven tectonic units: Kuqa Depression, Tabei Uplift, North Depression (Manjiaer-Awati Depressions), Central Uplift, Southwest Depression, Southeast Uplift, and Southeast Depression [32] (Figure 1(a)). The Tahe Oilfield, covering an area of 2400 km2, is located in the Akekule uplift in the eastern part of the Tabei Uplift (Figures 1(a) and 1(b)).

The warm and humid Ordovician climate provided favorable conditions for carbonate deposition on the Tarim Plate, which was situated at 20°S (Fang et al. 2001). From bottom to top, the Lower-Middle Ordovician strata are divided into the Penglaiba Formation (O1), Yingshan Fm. (), and Yijianfang Formation (O2), whereas the Upper Ordovician strata are divided into the Qiaerbak Formation (O3), Lianglitag Formation (O3), and Sangtamu Formation (O3) (Figure 2). The Early Ordovician paleogeographic framework of the study area was characterized by the widespread development of carbonate platforms that were open to the sea in an extensional setting [14]. The buildup of carbonate platforms of the Lianglitag Formation was much lower than that of typical Early to Middle Ordovician platforms as a result of the sea-level rise [34], with a maximum thickness of 145.5 m.

Core and detailed well log, seismic, and other exploration well data analyses indicate that, during the early stage, the Tahe region was characterized by a gentle slope with few high-energy depositions and characterized by grey to dark grey, thick-bedded micrites and microcrystalline marls interbedded with brown to grey-brown, thin- to medium-bedded bioclastic micrites and granular, microcrystalline limestones, whereas, at the late stage, it evolved into a weak-rimmed platform, characterized by the deposition of high-energy platform margin facies. The platform margin facies are distributed in TP5, S116, S110, and S107 wells, which are characterized by sparite oolitic limestone, calcarenite, and bioclastic limestone. Towards the south, the platform was mainly composed of dark grey, thinly laminated, nodular micrite, indicative of slope deposition (Figure 1(b)) [22, 23].

Previous studies [23, 35] have shown that there is an unconformity between the top of the Lianglitag Formation and the overlying Sangtamu Formation (EPU) (Figure 2), indicating a regional exposure surface and depositional break. The southward subduction of the Tianshan Ocean and formation of its associated suture at the end of Middle Devonian time resulted in thrusting and uplifting in the north of the Tahe Oilfield was significantly greater than in the south [31, 36] and led to the Upper Ordovician Lianglitag strata in the north being completely or partly removed by erosion, showing as the Lianglitag Fm. strata were overlain by the Carboniferous, the Devonian, or Silurian strata (LPU) (Figures 1(c) and 2), whereas the Lianglitag strata remained in the south (Figures 1(b) and 1(c)). In this paper, based on previous works on the basic geology of the Tarim Basin [31, 3739], the burial history of well S109 was rebuilt by the Basinmod 1D software and the Easy% Ro model. During the deposition of the Ordovician, the uplift of the Tahe area led to the denudation of the strata. After the Ordovician deposited, the burial history was roughly divided into two stages by the two periods of uplift at the Silurian-Carboniferous and Permian. There are many periods of weak erosion events, and, during the Neogene and Quaternary, the strata experienced rapid subsidence, leading to thick deposition. The Upper Ordovician strata have been buried to >4500 m, with paleotemperatures up to approximately 130°C (Figure 3).

3. Methods

This study is based on a detailed examination of 32 wells from the Tahe Oilfield. The description of cores focused on diagenetic features and also noted the lithology, depositional composition, and textures, as well as sedimentary structures. The petrographic examination of 120 thin sections of the diagenetic features and associated fabrics was completed using transmitted light and cathodoluminescence microscopy (CL, using a RELIOTRON, model RELION III; 5–8 kV and 300–400 μA gun current) at the Diagenesis Laboratory of Institute of Geology and Geophysics, Chinese Academy of Sciences (IGGCAS).

To study the fluid nature of vug-filling calcite cement, in situ δ13C and δ18O measurements were conducted using an in situ SIMS (Secondary Ion Mass Spectrometry) Cameca IMS-1280 ion microprobe at IGGCAS. The analytical procedures were similar to those described by Rollion-Bard et al. (2003). A 10–15 μm diameter Cs ion beam was used to ablate 10–15 μm diameter spots from the sample. The internal precision for each spot ranged between 0.009 and 0.015 (% Standard Error). The external precision was estimated to be 0.4 for oxygen and 0.6 for carbon as determined by consecutive analysis of a UWC (international calcite standard from University of Wisconsin) standard and OKA (Chinese primary calcite standard number GBW04481) standard.

44 core samples from 15 wells were collected for conventional oxygen and carbon isotope analyses. The samples analyzed in this study were whole rocks and fracture- and vug-filling calcite cement. Samples were taken using a dental drill. The analyses were carried out at the C-O Isotope Laboratory at IGGCAS. The carbon and oxygen isotope compositions were determined using a MAT-251 mass spectrometer (Vienna PeeDee belemnite [VPDB] standard), and the results were expressed in parts per thousand (delta units). Reproducibility was checked by replicate analysis of laboratory standards and was ±0.05 () for oxygen isotopes and ±0.02 () for carbon isotopes.

26 core samples from 7 wells were collected for 87Sr/86Sr analyses. Strontium isotope analyses were performed at the Rb-Sr Isotope Ultraclean Laboratory at IGGCAS. HAc was used to dissolve calcite to release Sr. The isotopic ratios were measured on a VG354 solid source mass spectrometer. The measured 87Sr/86Sr values were corrected using a quality fractionation standard of 86Sr/88Sr = 0.1194; the mean value of the measured strontium isotope NBS987 (standard sample) was .

Fluid inclusion petrography and microthermometry were conducted on 25 samples of intergranular calcite cement and vug-fill calcite cement to measure the homogenization temperature (, uncorrected for pressure) and salinity of the aqueous fluid inclusions. The salinity, in weight percent of NaCl equivalents, was calculated from the ice melting temperature ( ice) using the Bodnar [40] equation with a temperature step increase of 2°C/min. Special emphasis was placed on fluid inclusion assemblages (FIAs) [41], which represent the most finely discriminated petrographically associated groups of fluid inclusions, which were presumably trapped at approximately the same time. The measurements were performed using a Linkam THMGS600 heating-freezing stage with a 0.1°C precision and calibrated with synthetic pure water and CO2 inclusions. The analyses were conducted at the Diagenesis Laboratory at IGGCAS.

4. Results

4.1. Petrographic Description
4.1.1. Cement Sequences

Calcite Cement 1 (C1). C1 consists of isopachous fringe or fibrous (IFC) and early inclusion-poor calcite cement, constituting the first diagenetic phase to occlude the depositional pore space. IFC cement is noted for its blotchy to nonluminescent appearance under cathodoluminescence (CL) microscopy and is consistently present in grainstones and packstones (Figures 4(a) and 4(b)).

Calcite Cement 2 (C2). C2 is characterized by drusy calcite cement postdating calcite C1 and displays growth zones of dark red luminescence under CL (Figures 4(a) and 4(b)).

Calcite Cement 3 (C3). C3 is a blocky calcite and is characterized by nonluminescence (Figures 4(c)4(f)). C3 occurs not only in intraparticle and interparticle pores but also in vug. The volume of C3 varies with its stratigraphic position. It has been widely observed in cores within the Lianglitag Formation pinch-out boundary and Sangtamu Formation pinch-out boundary, shown in Figures 5(a) and 5(b). At the south of the Sangtamu Formation pinch-out, C3 mainly occurs in the wells that are close to the fracture.

Calcite Cement 4 (C4). C4 grows on top of C3 crystals and exhibits a bright orange luminescent zoned CL pattern. It has euhedral surfaces that show no signs of dissolution. The percentage of C3 and C4 can reach 80% in some wells (Figures 4(c) and 4(d)).

Calcite Cement 5 (C5). C5 is a clear calcite with brown to dull luminescence (Figures 4(e) and 4(f)). It occurs as isolated rhombic crystals or completely fills the pore space with blocky textures. Some dissolution likely occurred in C5 filled with oil and/or bitumen (Figures 4(e) and 4(f)).

Calcite Cement 6 (C6). The mineralogical characteristics of C6 are similar to those of C5, which has bright and thick crystals. Under luminescence, C6 mainly shows red luminescence (Figures 5(c) and 5(d)). C6 is closely associated with fractures that only present in wells where fractures were developed (Figures 5(e) and 5(f)).

4.1.2. Dissolution and Fracturing

Millimeter- to meter-scale dissolution features and fractures were observed in cores and suggestive of multiple dissolution events.

Early meteoric dissolution occurred locally in high-energy carbonate shoals, but this dissolution pore space was subsequently filled by late calcite cement (Figure 6(a)), with few effective fabric-selected pores reserved. Dissolution vugs and channels, with sizes up to several centimeters, were observed locally, and these vugs were mainly filled by calcite C3. The dissolution vugs were non-fabric-selective (Figure 6(b)) and were distributed along fractures (Figure 6(c)). Meter scale dissolution cavities were detected in some cores, with two types of cavity fillings. One was characterized by sands, silts, marls, and breccias (Figure 6(d)); and the other was characterized by calcite cement and marl, with little detrital quartz sand (Figure 6(e)). During drilling, lost cores, bit drops, and lost circulation occurred in nine wells (e.g., wells S106 and S110). The above observations suggest the development of large cavities in these carbonates.

Based on core measurements, the total fracture frequency is 8 fractures per meter, whereas the open fracture frequency is only 1.6 fractures per meter. Therefore, approximately 80% of the observed fractures are infilled, which limits fracture connectivity. Most of the fractures are partially or completely occluded by green marl (Figure 6(f)).

4.2. Isotope Geochemistry and Fluid Inclusion Microthermometry
4.2.1. In Situ Carbon and Oxygen Isotopes by SIMS

Carbon and oxygen isotopic data for pore-filling calcite cement from SIMS measurement were presented in Figure 7 and listed in Table 1.

C1. This calcite has δ18O values ranging from −3.62 to −0.10 VPDB, with a mean value of −1.32 (). The range of δ13C values on C1 fluctuates from 1.44 to 4.50 VPDB, with a mean value of 3.72 VPDB ().

C2. Compared with C1, the δ18O values of C2 are more negative, ranging between −6.35 and −2.56 VPDB, with a mean value of −4.24 VPDB (). The range of δ13C values on C2 fluctuates from −0.91 to 3.75, with a mean value of 1.55 VPDB ().

C3. The oxygen isotope values on C3 are more negative than C2, ranging from −8.45 to −6.50 VPDB, with a mean value of −7.59 VPDB (), and the δ13C values are from −0.07 to 4.07 VPDB, with a mean value of 2.32 VPDB ().

C4. The δ18O values are typically 3 higher than those of C3, lying between −5.37 and −2.29 VPDB, with a mean value of −4.01 VPDB (), and the δ13C values on C4 are close to those of C3, ranging from −0.32 to 4.90 VPDB, with a mean value of 2.21 VPDB ().

C5. C5 shows more negative δ18 values (−8.26 to −5.12 VPDB) than those of C4, with a mean value of −7.14 VPDB (), and similar δ13 values (0.45 to 2.38 VPDB).

C6. C6 displays low δ18O values ranging from −11.00 to −9.14 VPDB, with a mean value of −9.64 VPDB (), and δ13C values ranging from −0.88 to 1.70 VPDB (), except one low value (−3.88) of VPDB.

4.2.2. Traditional Carbon and Oxygen Isotopes

The micrite ranges between −8.05 and −3.78 VPDB for δ18O (mean −6.53 VPDB, ) and between 0.56 and 2.48 VPDB for δ13C (mean 1.92 VPDB, ). The oxygen isotopic compositions of grainstone are close to those of micrite, ranging from −8.45 to −3.23 VPDB (mean −5.98 VPDB, ), and the carbon isotopic compositions of grainstone varies from −0.64 to 2.27 VPDB (mean 0.97 VPDB, ). Fracture- and cave-filling calcite cement have δ18 values ranging from −11.02 to −8.16 VPDB (mean −9.71 VPDB, ) and δ13 values ranging from −2.16 to 3.16 VPDB (mean −0.02 VPDB, ), respectively. The analytical results are listed in Table 2 and shown in Figure 8.

4.2.3. Strontium Isotopes

The analytical results are listed in Table 2 and shown in Figure 9(a). The micrite spans a broad range from 0.709005 to 0.711927 (mean value of 0.710418, ). The grainstone yields 87Sr/86Sr ratios of 0.708136 to 0.709119 (mean value of 0.708659, ). The fracture and vug fillings present much higher 87Sr/86Sr ratios, ranging from 0.709103 to 0.712863 (mean value of 0.710377, ).

4.2.4. Fluid Inclusions (FI)

25 doubly polished detachable wafers of calcite samples were prepared and analyzed for FI. Petrography, both transmitted light and UV fluorescence, was performed on doubly polished wafers to identify inclusion types for further analysis and to determine their relationship to the host minerals. The large and/or elongate inclusions were disregarded, which could mean stretching and resetting since their formation [41].

The fluid inclusion homogenization temperatures and salinity (weight percent of NaCl equivalent) results of the cementing phases are shown in Figure 10. C1 is inclusion poor; therefore, no microthermometry was conducted in this cement.

C2 contains a few fluid inclusions, and the shape of the fluid inclusion is much more regular. The homogenization temperatures () range from 72°C to 91°C. The salinity, derived from the final melting temperature of ice using the equation from Bodnar [40], ranges from 4.96 wt% to 10.73 wt% NaCl equiv.

Fluid inclusions in C3 are more abundant and mainly occur as single-phase aqueous inclusions with sizes <4 μm in general, which are too small to be measured.

The characteristics of inclusions in C4 are very similar to those of C2, with fewer numbers and the shape of fluid inclusion being more regular. is less than 90°C, and salinity is moderate, ranging from 6.45 wt% to 9.98 wt% NaCl equiv.

C5 contains abundant aqueous fluid inclusions, sizes from 3 to 15 μm, and a vapor-to-liquid ratio of <15% in general. C5 shows a wide range of varying from 91°C to 129°C (). Water salinities at the time of C5 are variable, ranging from 0.71 wt% to 23.18 wt% NaCl equiv. (). Oil-filled fluid inclusions were observed and showed yellow or blue UV fluorescence.

Aqueous fluid inclusions in calcite cement C6 have a vapor-to-liquid ratio of <15% and sizes of 4 to 12 μm. The average homogenization temperature is high, ranging from 132.6°C to 200°C. Note that C6 with relatively high temperatures has relatively low-salinity values, less than 10.61 wt% NaCl equiv. Oil inclusions were observed with blue-white UV fluorescence emission colors.

5. Discussion

5.1. Origin of Digenetic Fluids
5.1.1. Stable Isotopes and Fluid Inclusion Constraints on Fluids

The combination of stable isotopic components and fluid inclusions can be used to estimate the formation temperature of diagenetic cement, which can further refine the diagenetic sequences and define the fluid properties.

Fibrous isopachous calcites are considered to be early marine cement that formed in the marine phreatic zone [42]. Therefore, C1 in this study formed in the early marine diagenesis environment. One type of blocky calcite cement surrounding bioclasts shows nonluminescence with several thin yellow bands (Figures 4(a) and 4(b)) based on similar δ18O values and δ13C values with fibrous isopachous calcites (C1) (Figure 7). Hence, this blocky calcite cement should precipitate in marine diagenesis environment as well. Shields et al. [43] have reported that δ18O values of −6 to −4 VPDB and δ13C of −2 to 2 VPDB characterize the global marine carbonates carbon and oxygen isotope values during the Late Ordovician. Assuming that the Late Ordovician limestones originally had a temperature of 25°C for subtropical seawater [44, 45] and using the O’Neil et al. [46] calcite-water oxygen isotope fractionation equation, the δ18O value of the water present during the growth of C1 was between −1.22 and 2.4 V-SMOW, higher than that of Late Ordovician seawater, that is −3 to −6 V-SMOW [43]. The slightly higher oxygen isotope values can be achieved by evaporation [47]. Therefore, the higher δ18O value suggests that evaporative concentration of seawater may have occurred during precipitation of C1. The higher δ18O value was mainly found in well S106, located at the southern part of the study area and characterized by deep slope facies (Figure 1(b)). Hence, environmental influence may be an alternative interpretation for oxygen isotopic changes in carbonate.

Compared with C1, the C2 shows slightly lower δ18O values (−6.35 to −2.56 VPDB; Figure 7), and homogenization temperature from fluid inclusion study indicates that C2 with coarser crystal was formed during burial diagenesis at approximately 60–90°C (Figure 10).

C3 shows lower oxygen values compared with C2. The more negative δ18O values could be due either to precipitation at a higher temperature or to precipitation from waters with lower δ18O values [48, 49]. The single-phase aqueous fluid inclusion petrographic evidence, which indicates a capture temperature of less than 50°C [41], suggests that C3 probably was in meteoric origin. Previous studies have indicated that the Lianglitag Formation carbonates are affected by both Caledonian and Hercynian orogenic events [5053]. Hence, C3 probably precipitated as a result of meteoric water diagenesis during Caledonian and/or Hercynian orogenic events.

Compared with C3, C4 shows much more positive δ18O values. Most of the δ18O values for C4 lie within the marine calcite ranges [43], indicating that C4 was most likely influenced by marine fluid. data from the fluid inclusions in C4 shows that it formed in temperatures up to 90°C; hence, C4 probably precipitated during the shallow buried stage after the tectonic uplift. Melim et al. [54] also showed that the carbon and oxygen isotopes changed from negative to positive from the meteoric phreatic zone to a seawater-shallow buried environment.

C5 has relatively negative δ18 values, with δ18 values ranging from −8.26 to −5.12 VPDB. The high fluid inclusion temperatures (91°C to 129°C) and wide range of salinity (0~23.18 wt%) imply that C5 blocky spars may have precipitated during burial environment from the mixing of multiple sources of fluids.

The δ18 values of C6 are obviously negative and could be as low as −11 VPDB. Fluid inclusions in C6 are shown to have relatively low salinities (from 0 to 10.61 wt% equiv. NaCl) but highest homogenization temperatures (132.6 to 200°C). Based on the δ18O values and measured temperatures, the isotopic signature of the parent fluid ranges from 5.57 to 11.96 (SMOW). Combined with the regional burial history curve (Figure 3), the homogeneous temperature of C6 exceeded the maximum burial experience of the strata. This type of calcite may have precipitated from a variety of processes: (1) the mixing of meteoric water with hydrothermal water, (2) thermochemical sulfate reduction, and (3) deeply recycled meteoric water. Hence, the high precipitation temperatures and relatively low salinity of the fluids indicate that C6 may have precipitated from meteoric water mixed with hydrothermal fluids. However, if this model applies, fluid inclusions with high temperatures are expected to have high salinities rather than low salinities, as was found in this area. The second possibility is that C6 with high homogenization temperatures and low salinities may have precipitated from thermochemical sulfate reduction (TSR) derived water [15, 55]. TSR is redox reactions between hydrocarbon and dissolved sulfates in temperatures >120°C [5658]. TSR is well documented in Tazhong area [57, 59]. However, this proposal is discounted by the nonisotopically light carbon in both vug fillings and intergranular cement. Alternatively, deeply recycled meteoric water could have flowed down along faults to deep strata and thus became significantly heated. Subsequently, the heated meteoric water underwent upward migration, leading to the precipitation of hydrothermal calcite. This model of water flow was previously proposed by Qing and Mountjoy [60] to explain the variation of homogenization temperatures as well as the 87Sr/86Sr and δ18O values from saddle dolomite cement from the west towards the east in the Western Canada sedimentary basin as a result of tectonic thrusting and compression, sedimentary loading, and tectonic uplift on the western margin. The structural history of the Western Canada basin is broadly analogous to that of the Tahe area. Thus, it is possible for a similar model of water flow to have occurred in the study area as previously proposed by Li et al. [15], Jiang et al. [61], and Jiang et al. [62].

Different from δ18O, δ13C is mainly controlled by carbon sources [63]. Most of the δ13C values from C1 to C6 lie between −2 and 4 (Figure 7), close to the δ13C range of contemporaneous seawater, indicating a local source of C.

Compared with the δ18O and δ13C analysis results of matrix pore-filling calcite cement by the SIMS method, the δ18O and δ13C data of fracture- and cave-filling calcite cement obtained by traditional methods are more negative (Figure 8), most likely due to the results of karst, which had been discussed by Li et al. [15] in detail. Most of the oxygen isotopes of micrite and grainstone are obviously more negative than those of the Late Ordovian marine carbonate (from −6 to −4 VPDB) [43], which indicates that the matrix limestone in the Lianglitag Formation is strongly affected by the late diagenetic fluid.

5.1.2. Strontium Isotopes

During the Late Ordovician, the carbonate strontium isotope composition of the global ocean ranged from 0.7078 to 0.7082 [64, 65]. The seawater 87Sr/86Sr record is controlled by two main strontium fluxes and their isotopic compositions: the hydrothermal brines at midocean ridges (MOR flux), characterized by 87Sr/86Sr ratios being close to 0.703, and the continental input (river flux), which has a modern global average 87Sr/86Sr ratio of 0.712 [6668].

Both micrite and fracture-vug-fill calcite present much higher 87Sr/86Sr ratio values than Late Ordovician seawater, indicating largely introduced radiogenic Sr (Figure 9(a)). The higher 87Sr/86Sr ratios of some micrite samples may be the result of recrystallization, which was commonly observed under the microscope (Figures 5(c) and 5(d)); it could also reflect the introduction of radiogenic strontium from terrigenous clastic internal to the Lianglitag Formation [23]. For radiogenic Sr in calcite cement, the stratigraphic hiatuses identified in Figure 2 and burial history diagrams shown in Figure 3 suggest that there are several uplift and erosion events which could have involved the weathering of radiogenic source terrains. Radiogenic strontium could also be derived from the Precambrian basement or Lower Cambrian clastic section via circulating fluids within conductive faults and fractures that penetrate the basement in the Tazhong area [69]. 87Sr/86Sr ratio shows increasing trends from lower to upper part of the formation. Therefore, the 87Sr-rich fluid is probably from meteoric water leaching of clastic material in the Tahe area during the periods of tectonic uplift, and this is supported by a recent study from formation water geochemistry in the Tabei area [70].

Previous studies have indicated that the Lianglitag Formation carbonates are likely affected by both Early Paleozoic Caledonian and Late Paleozoic Hercynian karstification [5053]. However, because of the relatively similar characteristics in the minerals and cathodoluminescence, it is difficult to distinguish these phases. According to the stratigraphic distribution characteristics (Figure 2), one of the distinguished features of the Late Paleozoic is that the clastics are more obvious. Hence, strontium isotope may be useful in distinction between the two karsts. The strontium isotopic composition could be broadly divided into two groups (Figure 9(b)): the first group represent the Early Paleozoic karst, characterized by 87Sr/86Sr ratios that range from 0.70712 to 0.7085, which is close to the strontium isotopic composition of the Late Ordovician seawater; the second group shows 87Sr/86Sr ratios ranging from 0.7085 to 0.71188, significantly higher than the Late Ordovician marine strontium isotopic composition. This group represents Late Paleozoic karst, which indicates that strontium may be added by continental weathering through the river. In addition, the strontium isotopic composition of the northern part of the Tahe Oilfield (the north of the Sangtamu Fm. pinch-out boundary) is much higher, which is mainly due to the Late Paleozoic karst, whereas the strontium isotopic composition of the middle and southern parts of the Tahe Oilfield (the south of the Sangtamu Fm. pinch-out boundary) are broad and are both influenced by the Early Paleozoic and the Late Paleozoic karst, such as well S108 and well T901, which are superimposed on the two karst interactions.

5.2. Diagenetic Fluid Evolution and Its Influence on Reservoir Properties

Based on petrographic, isotopic, and fluid inclusion data, the diagenetic sequence is identified and divided into three main episodes. The eogenetic, mesogenetic, and telogenetic stages were used in this study [71].

During the deposition of the Lianglitag Formation, the Tahe region was characterized by a ramp and/or weakly rimmed platform, with a low-gradient slope and small-scale reefs and shoals [23]. Marine calcite cement (C1) is commonly observed (Figures 4(a) and 4(b)), especially in the southern slope facies, indicating that the Lianglitag Formation carbonates experienced strong marine cementation in the eogenetic stage. Minor fabric-selective pores occur in high-energy shoals; however, they are more commonly partly or completely filled by calcite cement (Figure 6(a)). Therefore, the carbonates exhibit intensive marine cementation and weak syndepositional freshwater dissolution during the eogenetic diagenetic stage (Figure 11(a)). Thereafter, the strata enter a shallow burial diagenetic environment where C2 calcite cement was precipitated.

The large volume of C3 precipitated at the wells demonstrates that the telogenetic meteoric impact is ubiquitous. In the Caledonian, the arch developed as a monocline dipping gently to the SE; as the arch was uplifted several times, three episodes of karstification, designated episodes I, II, and III, occurred during this period [52, 53]. Episode I occurred between the deposition of the Yijianfang and Qiaerbake Formations. Episode II occurred between the deposition of the Lianglitag and Sangtamu Formations. Episode III occurred between the deposition of the Ordovician Sangtamu Formation and Silurian strata. The Lianglitag Formation is mainly affected by Episode II [17]. The Lianglitag Formation carbonates were exposed to the surface and infiltrated by ground water, leading to significant dissolution of the carbonates and formation of dissolution vugs and caves with sizes from centimeters to meters (Figure 11(b)). Some vugs and caverns have been filled by minerals and/or sediments during the following mesogenetic stage (Figure 11(c)).

At the end of the Middle Devonian, because of the proximity to the South Tianshan Ocean, the northern part of the Tahe Oilfield was uplifted and subjected to subaerial erosion once again [31, 36]. The Lianglitag subcrop shown in Figure 1(b) and erosional removal of the Middle Ordovician through Carboniferous section shown in Figure 2 are both the result of Hercynian uplift and erosion. Therefore, Hercynian telogenetic fluids could have laterally access to the Lianglitag Formation via the subcrop limit and result in the development of dissolution cavities (Figures 1(c) and 11(d)). Access would not necessarily require transit through the overlying Sangtamu section. Additionally, fractures are important diagenetic fluid channels for karstification. There are four periods of fault structures recognizable in the Tahe seismic profiles at present [14], which from the bottom up are as follows: (1) thrust faults developed in Cambro-Ordovician strata that formed during the Middle-Late Ordovician; (2) a strike-slip fault that developed under the reflecting layer (i.e., Carboniferous base boundary) and formed in the Middle-Late Devonian; (3) normal faults cutting Permian and older strata associated with early to middle Permian magmatism; and (4) thrust faults that developed in pre-Jurassic strata and formed during the latest Permian to earliest Triassic. The Late Ordovician strata in the Tahe region are mainly affected by mid to late Devonian strike-slip faulting. Other periods of faulting are of subsidiary importance. Therefore, strike-slip faults in the Middle-Late Devonian developed and freshwater along fractures could affect the Lianglitag Formation strata.

Subsequently, the sediments were again influenced by marine to shallow burial, from which bright orange CL calcite (Cement 4) precipitated. Burial dissolution pores are locally visible (Figures 4(e), 4(f), 5(e), and 5(f)). Although late burial dissolution does occur, it is volumetrically minor. Previous studies have shown that burial dissolution mainly includes organic acids, CO2, H2S [72], and deep hydrothermal dissolution [12, 73, 74]. According to Section 5.1.1, TSR and thermal fluid activity are not obvious in the Upper Ordovician Lianglitag Formation, which is consistent with the results of Li et al. [14]. Deeply recycled meteoric water may have entered the carbonate rocks during the late Permian favored by fracturing [75] and generated both dissolution and subsequent cementation C6 in vugs. As noted above, C6 was precipitated at temperatures from 130 to 200°C and showed significantly higher 87Sr/86Sr ratios and lighter δ18O and δ13C values. We have observed some C6 temperature and oxygen isotopic decreasing trend in distance, and the trend of oxygen isotope values and temperatures during its pathway suggest that C6 was likely precipitated from recycling of meteoric water. Unfortunately, the data are not enough. Combined with previous works in Tarim Basin [15, 61, 62], C6 was probably closely related to deeply recycled meteoric water. The reduced water coexisted with highly mature oils with blue-white UV fluorescence emission colors. However, because most fractures in the Tahe region are partially or completely occluded by green marl, the development of deeply recycled meteoric water is very local (Figure 11(e)).

Non-fabric-selective vug and cave were produced during near-surface karst processes and later burial diagenesis forms an important class of carbonate reservoir in the Lianglitag Formation in the Tahe Oilfield. Being near surface combined with burial diagenesis produced the complex reservoirs with several scales of heterogeneity. The main controlling factor of the reservoir distribution was multistage karstification and fracturing.

6. Conclusions

(1) Combined with isotopic geochemical and fluid inclusion data, the corresponding diagenetic fluid was identified. The early nonluminescent C1 to growth zones of dark red luminescent C2 is a syngenetic seawater-shallow burial diagenetic environment; the subsequent nonluminescent C3 is related to the meteoric fluid system resulting from tectonic uplift and exposure. Then, the strata were buried again and the bright orange luminescent C4 zone was formed during shallow burial. The C5 precipitated in a middle-deep burial environment and was characterized by brown to dull luminescence. In addition, an abnormally high temperature and low salinity were found in C6 along fractures, which may be the product of deeply recycled meteoric water.

(2) The abundance of C1 indicates that the Lianglitag Formation carbonate suffered strong marine cementation and weak meteoric dissolution and cementation during the eogenetic diagenetic stage. The paleokarst during Caledonian II and Hercynian telogenetic karstification closely related to C3 constitutes the most constructive diagenesis. In addition, the carbonate was locally infected by deeply recycled meteoric water during burial. Compared with the Early Paleozoic fracture- and vug-filling calcite, the Late Paleozoic fracture- and vug-filling calcite have a significantly higher 87Sr/86Sr ratio due to the presence of more terrigenous clastic material.

(3) The reservoir system in the Tahe carbonates is dominated by non-fabric-selective pores and remanent (unfilled) karst caves. The development of the reservoir is closely related to fractures. The vugs and caves may be the result of the infiltration and dissolution of meteoric fluid along the fractures when the strata are tectonically uplifted, of high-temperature and low-salinity meteoric water circulating along the deep basin, or of the two superimposed transformations. Multistage karstification and fracturing are the main controlling factors of reservoir development, which can provide important guidance for future deep-buried oil and gas exploration, and have important economic implications.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is financially supported by the National Natural Science Foundation of China (Grant no. 41302084), Strategic Priority Research Program Grant of the Chinese Academy of Sciences (XDA14010201-02), and National Science and Technology Major Project of China (2017ZX05008-003). The authors thank Drs. Xiaoxiao Ling, Guoqiang Tang, Yu Liu, Qiuli Li, and Xianhua Li for their help during SIMS analyses in IGGCAS. Dr. Lei Jiang gave constructive advice both in English language and in logic arrangement. Many thanks are also due to Professor Chunfang Cai for constructive comments that helped to improve the manuscript.