The Middle Jurassic Toutunhe Formation at depths of approximately 4000–6000 m has increasingly come into focus as a current deep reservoir target in the central Junggar Basin (NW China). Based on petrography, SEM, stable isotopes, and fluid inclusion analyses, the goals of this study were to investigate the effect of depositional lithofacies on sandstone diagenetic heterogeneity and to examine the relationship between diagenetic evolution and oil charge within a heterogeneous reservoir. Grain size controls the overall abundance of cement and porosity and reservoir properties through its effect on ductile lithic sand grains and hence on mechanical compaction. Early diagenetic calcite cement is an exception to this trend. Ductile lithic-rich, very fine-grained sandstones featured compaction of easily deformed, clay-rich grains, resulting in a very rapid loss of porosity during burial. In contrast, dissolution and cementation occurred as well as ductile compaction in the fine-grained sandstones. Two episodes of oil charge occurred in the relatively coarser-grained sandstone lithofacies. Diagenesis progressed alternately with oil emplacement, and some diagenetic alterations and oil charge occurred simultaneously. Ductile lithic-rich, highly compacted sandstones and tightly calcite-cemented sandstones can create permeability barriers embedded in permeable reservoir sandstones, probably resulting in heterogeneous flow.

1. Introduction

Reservoir heterogeneity refers to spatial variability of a system property (e.g., porosity, permeability, and pore structure) at a specified scale [14]. The intrinsic complexity of a reservoir is the result of various geological processes, including depositional environment, diagenesis, and tectonic imprint. The heterogeneity is a critical factor affecting fluid flow and distribution within a reservoir [5, 6]. Diagenetic studies that are pertinent to reservoir characterization have been confined almost entirely to the identification and quantification of mineral components and porosity types by combined use of scanning electron microscope (SEM) and thin section and surely do not characterize diagenetic heterogeneity, especially at scales larger than microscopic scale. The primary texture and detrital compositions of sandstones, sand body geometry, and architecture, ultimately controlled by depositional environment, are of decisive importance in determining the rate and pathways of fluid flow affecting the extent and distribution of early diagenetic alterations in sandstones [712]. Early diagenesis in turn exerts a substantial control on the deep-burial diagenetic evolution of the sandstones [9, 10, 13]. It seems plausible that a significant component of sandstone diagenetic heterogeneity arises from initial variations in texture and grain compositions in principle [8, 14]. Thus, linking diagenetic alterations to depositional environments and lithofacies has been an important concept for evaluating a degree of large-scale heterogeneity characterizing diagenetic features [9, 10]. This is essential to achieve a complete and adequate spatial representation of a heterogeneous reservoir.

Hydrodynamic fluctuations and varying depositional facies associations directly determine textural and compositional parameters of sandstone sediments, producing differential compaction and diagenetic processes, as well as reservoir quality evolutions [1517]. Some sandstone lithofacies have experienced a very rapid loss of porosity through compaction or cementation and have achieved a high capillary entry pressure during early diagenesis and have acted as low-permeability barriers to late fluid (hydrocarbons included) flow in a reservoir [1719]. However, studies of differential sandstone diagenesis and porosity evolutions at a reservoir scale level are still far from being fully explored in the literature.

In recent years, the Jurassic lithologic and stratigraphic traps have become a major focus of deep reservoir exploration in the central Junggar Basin (NW China). Some oil fields such as Mobei, Mosuowan, and Moxizhuang oil fields have been discovered in the Jurassic sandstones (Figure 1(a)), reflecting promising prospects for deep exploration. No significant discoveries have been obtained in the Fukang Sag. By the end of 2014, the Shengli Oilfield Company, SINOPIC, has completed seven wells in this area (Figure 1(b)) and has obtained substantial findings in the Lower-Middle Jurassic sandstones. The recoverable oil original reserves were estimated to be approximately 2109.35 × 104 t (155 million bbl). The focus of this study is the Middle Jurassic Toutunhe Formation (Figures 2 and 3). This formation is buried at depths of approximately 4000–6000 m. The sandstones have strongly heterogeneous petrophysics with porosities of 2.1–16.4% and permeabilities of 0.02–28.9 mD. The inhomogeneity of oil shows has been noted from core description, well logs to well tests. The aim of the study is to investigate the effect of depositional environments and lithofacies on sandstone diagenetic heterogeneity and to decipher the diagenetic evolution pathways of a heterogeneous reservoir and its impact on hydrocarbon migration.

2. Geological Setting

The Fukang Sag is located in the Junggar Basin, which covers an area of approximately 136,000 km2 (52,509 mi2) north of Xinjiang Uygur Autonomous Region, northwest China (Figure 1(a)). The Junggar Basin is triangle-shaped and is confined by the Qinggelidi and Kelameili Mountains to the northeast, by the Yilinheibiergen and Bogda Mountains to the south, and by the Zhayier and Halalate Mountains to the northwest. The basin lies at a three-block junction between Kazakhstan, Siberia, and Tarim and is a part of the Kazakhstan block [20]. The Junggar Basin is floored by Precambrian crystalline basement and part of Hercynian folded basement [21] and is an Upper Paleozoic, Mesozoic, and Cenozoic superimposed basin. The basin has experienced the Hercynian, Indosinian, Yanshanian, and Himalaya orogenies. Three main phases of deformation have been recognized in the basin: the Middle-Late Permian rift stage, Triassic to Paleogene depression stage, and Neogene to Quaternary foreland stage [22, 23] (Figure 3). The Junggar Basin contains six structural units: the Wulungu Depression, Luliang Uplift, Western Uplift, Central Depression, Southern Overthrust Belt, and Eastern Uplift [24]. Each structural unit has a distinct structural style and tectonic history. The Fukang Sag lies in the Central Depression (Figure 1(b)). The main structure within the sag is a northeast-southwest trending monocline (Figure 2), which formed as a result of the uplift of the Bogda Mountains at the southern margin of the basin since the Paleogene [24]. The Jurassic strata are relatively flat lying or dip gently to the southwest with a dip angle of 2-3° (Figure 2). Most of the observed faults are northeast-directed and east-west-directed normal faults confined to the Jurassic strata (Figures 1(b) and 2). Throw of the fault is generally tens to hundreds of meters and strike length is less than 2 km. These normal faults and deep reversed faults, in some localities, show a flower structure (Figure 2).

The central basin contains Carboniferous to Quaternary sediment fills up to 10 km thick (Figure 3). Rather than attempting an encyclopedic coverage of all the sediments, here, we only summarize general depositional evolution of the Jurassic succession. During the Early Jurassic and early Middle Jurassic period, the central basin was an intracontinental shallow lacustrine basin under the low-amplitude oscillating tectonic regime. Depositional facies identified in the central basin include alluvial fan and fluvial, deltaic, and lacustrine facies [25, 26]. The deposits are typical of a coal-bearing succession and consist of conglomerates, sandstones, and mudstones (Figure 3). Due to the Early Yanshanian orogeny, the lacustrine area began to shrink rapidly during the late Middle-Late Jurassic period. Fluvial, deltaic, and shallow lacustrine facies are included. The deposits are sandstones and mudstones (Figure 3). The contact between the Toutunhe Formation and the underlying Xishanyao Formation is unconformable. The Toutunhe Formation has local erosions, and the Upper Jurassic succession is entirely absent (Figure 3).

Previous exploration and geochemical studies indicate that the Fukang Sag contains four sets of potential source rocks (i.e., Carboniferous, Permian, Triassic, and Jurassic) (Figure 3). Among these source rocks, most discovered hydrocarbons in the Jurassic strata are generally believed to be originated from both the Permian and Jurassic source facies [2729]. As these source rocks are deeply buried and not easily penetrated by the drill, previous published organic geochemistry studies have been almost entirely restricted to the northern and eastern slope of the sag [29, 30]. The Lower Permian Pingdiquan Formation contains deep lacustrine black mudstones, oil shales, and dolomitic mudstones with a total thickness of 50–650 m [30]. The mudrocks have total organic carbon (TOC) contents ranging from 0.5 to 10.2%. The organic matter is of type I-II kerogen with vitrinite reflectance (Ro) values of 1.2–2.0% [29]. It can be inferred that the source rocks are mature to over mature towards the sag center [30]. The source rocks in the Lower Jurassic and lower Middle Jurassic succession are composed of lacustrine-swamp deep gray mudstones and carbargilites interbedded with thin coal-beds. The mudstones have a total thickness of up to 500 m [30]. The TOC contents range from 0.3 to 3.7%. The organic matter is of type II-III kerogen with Ro values of 0.6–1.0% [29]. The Badaowan Formation’s source rocks are of better quality than the Sangonghe and Xishanyao Formation. It is speculated that these source rocks are mature to highly mature towards the sag center [30].

In the Fukang Sag, the Toutunhe Formation’s deposits were derived from Kelameili Mountains in the northeast and Bogda Mountains in the south [3133]. During deposition of the Toutunhe Formation, the study area witnessed a variety of depositional environments ranging from shallow lacustrine through to delta to meandering fluvial system with a total thickness of 700–900 m [33]. The formation can be subdivided into three lithological units (Figure 4). The lower unit consists of grayish-green mudstones and siltstones interbedded with very fine- to fine-grained sandstones. The middle unit includes grayish-green very fine- to fine-grained sandstones, siltstones, and mudstones. The upper unit comprises mainly greyish-brown gravel-bearing medium-grained sandstones, fine-grained sandstones, very fine-grained sandstones, and interlaminated siltstones and mudstones. The Lower Cretaceous mudstones are widely distributed with a total thickness of 140–180 m and can act as the main regional seals for the Jurassic hydrocarbon plays.

3. Samples and Methods

Cores were logged and sampled from the D-1, D-2, D-3, D-6, D-7, and D-8 wells (Figure 4). The specific lithofacies and depositional facies associations of the sandstones were defined based on core description (proportion of sand to shale, bed thickness, grain size, degree of bioturbation and deformation, sedimentary structure, presence of bioclasts and mudstone intraclasts, carbonate cement distribution, etc.). This was combined with wireline logging interpretation to delineate the overall depositional environments for the Toutunhe Formation. All cores were examined under ultraviolet (UV) light for reservoir properties (oil saturation, oil immersion, oil patch, fluorescence, and no show).

Twenty-eight conventional core samples were selected to represent the various depositional facies and lithofacies and the different reservoir properties (Figure 4). Mineralogical, stable carbon and oxygen isotope, and fluid inclusion analyses were performed at the Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China.

All samples were impregnated with blue epoxy resin under vacuum to indicate porosity prior to preparing polished thin sections. Sandstone petrographic compositions were quantified by 300-pointed modal analysis. In addition, the thin sections were studied to identify growth occurrence and habits of interstitial diagenetic minerals and to define textural relationships between different cements. Cathodoluminescence (CL) observations were made using a Relion III cold-cathode device at beam voltages of up to 15 keV and 500 μA. To further determine their textural relationships, selected samples were examined using a high-resolution field-emission Nova NanoSEM 450 scanning electron microscope equipped with an integrated analysis system of the Oxford Aztec Energy Dispersive X-ray Spectrometry (EDS) and electron backscatter diffraction (EBSD). Carbon-coated thin sections were mounted on gold grids. All backscatter electron (BSE) images were acquired at 20 keV at a working distance of 5–8 mm. EDS was used for qualitative identification of minerals.

Fourteen samples containing a single generation of calcite cement (either early or late) were selected for carbon and oxygen isotope analysis. The samples were ground to >200 mesh and then soaked overnight in 3% sodium hypochlorite to remove organic matter. Subsequently, the powder samples were reacted with 100% phosphoric acid at 25°C for 4 hours, and CO2 was extracted and analyzed. The entire process was finished using a Finnigan Delta S mass spectrometer. Oxygen isotopic compositions were calculated using a phosphoric acid-CO2 fractionation factor of 1.0125 at 25°C. All results were represented in standard notation relative to PDB standard. The reproducibility of duplicate analyses generally was better than ±0.2%.

Nine samples were collected to prepare 100 μm thick doubly polished sections for petrographic observation and microthermometric measurement of fluid inclusions. Petrography of pore bitumen (free oils and solid bitumens) and fluid inclusion and fluid inclusion assemblages were firstly examined using a Nikon 80I microscope under both transmitted and UV light. The wavelength of the emission fluorescence is greater than 420 nm. Homogenization temperature measurement of fluid inclusions was carried out using a calibrated Linkam THMSG600 heating-freezing stage. A step-heating and temperature cycling technique [35] was adopted to record the homogenization temperatures. The precision is ±0.1°C. Fluids inclusions displaying textural or microthermometric evidence for stretching, necking down, leakage, or decrepitation were excluded from the study.

4. Results

4.1. Core Description

Sedimentological description of all cores from six wells with a total length of 53.5 m has been done (Figure 4). Five depositional facies associations of the sandstones have been distinguished: (1) single and stacked story meandering fluvial channel, (2) delta distributary channel, (3) distributary-mouth bar, (4) interdistributary bay, and (5) shore-shallow lacustrine facies association. These facies associations have been identified based on core sections and well log data and are briefly described as follows. Single and stacked story channels are up to 5 m thick, fining-up units ranging mainly from blocky gravel-bearing medium- to very fine-grained sandstones, with minor siltstone and mudstone laminas. The meandering channel facies association dominates the whole upper Toutunhe Formation in the study area. Delta distributary channels are fine-grained sandstones, grading upwards into siltstones. The individual channel units are up to 3 m and sometimes are erosional-based. Sedimentary structures, in some localities, are blocky lamination at the base, changing gradually upwards into planar cross bedding towards the top of the unit, where current ripple is common. Distributary-mouth bars comprise up to 3 m thick, coarsening-up units from siltstones to fine-grained sandstones. The very fine-grained sandstones often display current ripple. Interdistributary bay sediments are almost mudstones interlaminated with thin and lenticular very fine-grained sandstones and siltstones (10 to 20 cm thick). The sandstones often exhibit plane parallel lamination in which bioclasts are common. Delta distributary channel, distributary-mouth bar, and interdistributary bay facies associations are mostly limited to the middle and lower units of the Toutunhe Formation in the D-2, D-6, D-7, and D-8 wells. Shore-shallow lacustrine fills are bioturbated mudstones and interbedded very fine-grained sandstones and siltstones displaying current ripple and plane parallel lamination. Mud intraclast-bearing fine-grained sandstones are seen locally. Patchy pyrite and bioclasts are commonly found. Shore-shallow lacustrine facies association dominates the middle and lower Toutunhe Formation in the D-1 and D-3 wells. In the first three facies associations, calcite commonly forms discrete and tightly cemented concretions (up to 30 cm) embedded in weakly cemented or uncemented sandstones (Figure 6(e)), in particular, in the relatively coarse-grained part of graded sand beds (Figure 4).

4.2. Sandstone Petrology

The Toutunhe Formation’s sandstones are very fine- to fine-grained sandstones with the median grain size values of 0.08–0.25 μm. Sorting ranges from moderate to well-sorted. The roundness of detrital sand grains varies from subangular to subrounded.

The sandstones are mostly feldspathic litharenite and litharenite (Figure 5; [34]). The average framework composition is 20.4% quartz, 23.3% feldspar, and 56.4% rock fragments (Q20.4F23.3L56.3). The detrital components of the sandstones are dominated by volcanic rock fragments, averaging approximately 36% of the whole rock volume (WRV). The igneous textures recognized in thin sections include lathwork, microlithic, cryptocrystalline, felsitic, and vitric types. Most volcanic rock fragments are composed of euhedral lathwork and microlithic feldspar phenocrysts in the cryptocrystalline groundmass which displays characteristic alteration to clays (Figures 6(a)6(d)). Metamorphic rock fragments have an average content of about and comprise metamorphosed siltstone, mica schist, slate, and phyllite. Sedimentary rock fragment contents are generally less than . They are mostly shale and chert. Some of the lithic fragments are variably altered or deformed and have vague or irregular outlines (Figures 6(a)6(d)) and termed as “ductile” lithic sand grains. The percentage of the ductile grains is . Detrital quartz as a fraction of the whole rock averages and is mostly monocrystalline with minor amounts of polycrystalline grains. The polycrystalline quartz features sharp extinction and resorption embayments, which is interpreted as metamorphic origin. Total detrital feldspar averages . Plagioclase is the most abundant feldspar type observed, with an average ratio of plagioclase to total feldspar of about 0.7. In addition, detrital clay matrix is locally abundant, up to .

4.3. Diagenetic Heterogeneity

Three major sandstone petrofacies have been defined on the basis of texture and detrital compositions, diagenetic features, and pore space properties (i.e., pore abundance and type). Diagenetic features of the three defined sandstone petrofacies are described below.

4.3.1. Ductile Lithic-Rich Sandstone

Ductile lithic-rich sandstones are mainly ripple cross-laminated very fine-grained sandstones (median grain size: 0.08–0.13 μm) that occur in all facies associations (Figure 4). Ductile lithic grains constitute 22–40% of the whole rock in these samples; the contents of detrital matrix are (Table 1). Diagenesis of the sandstones is characterized by compaction of easily deformed, clay-rich grains. Ductile grains are severely deformed by mechanical compaction and are squeezed between rigid detrital grains, blocking pore throats (Figures 6(a)6(c)). In addition to clay minerals almost indiscernible in thin sections, pyrite is an important authigenic phase in the ductile lithic-rich sandstones, ranging from to . Pyrite occurs mainly as clusters of cubes in interstitial pores and is an early diagenetic product (Figures 8(f) and 9(h)). The sandstones have no visible thin porosity (Figures 6(a)6(c)).

4.3.2. Tightly Calcite-Cemented Sandstone

Tightly calcite-cemented sandstones occur in the relatively coarse-grained lithofacies of a graded sand bed (Figure 4), having the median grain size values of 0.17–0.19 μm (Table 1). Tightly calcite-cemented sandstones contain calcite in excess of . The calcite cements, which are composed mostly of anhedral calcite crystals (Figure 9(a)), represent a poikilotopic texture and local replacement of framework grains (Figures 6(f) and 6(g)). The calcite precipitated before significant compaction, thereby limiting further diagenetic modifications and resulting in a high intergranular volume. The contents of ductile lithic grains are in these samples. The early diagenetic calcite cements exhibit orange-red luminescence emission (Figure 6(g)). The EDS analysis indicates that the calcite has homogeneous almost 100% CaCO3 (Figure 8(g)). The calcite has δ13C values of −4.3 to and δ18O values of −13.4 to (Table 2). The sandstones have also no visible porosity in thin sections (Figure 6(f)).

4.3.3. Permeable Reservoir Sandstone

Permeable reservoir sandstones are blocky to planar cross-stratified, fine-grained sandstones deposited predominately in the fluvial channel, distributary channel, and distributary-mouth bar facies associations. Ductile grains are less abundant in these samples than in ductile lithic-rich sandstones, ranging from to ; the contents of detrital matrix are less than (Table 1). Authigenic minerals include calcite, clay, analcime, anhydrite, barite, quartz overgrowth, albite, and hematite. These cements display a patchy or sporadic texture. Primary intergranular porosity (Figures 6(h) and 6(i)) as counted in thin sections averages , while secondary porosity averages and is affiliated with feldspars and lithic fragments and even minor cements (Figures 6(i), 6(j), and 7(d)).

The calcite contents are low, generally less than , and display a rather patchy and local cementation (Figures 7(a), 7(b), and 9(b)). The calcite cements formed after significant compaction and postdated chlorite coats, quartz overgrowths, and feldspar overgrowths or authigenic albite (Figures 8(a) and 8(c)). The calcite displays yellow to orange luminescence (Figure 7(b)). The EDS analysis reveals that the calcite contains Mn in some samples (Figure 8(h)). They have δ13C values of −8.4 to and δ18O values of −19.0 to − (Table 2).

Chlorite has contents of less than . Chlorite cement seems to be facies-related and has been only observed in the distributary channel and distributary-mouth bar facies associations. The chlorite occurs as thick, continuous grain coats and not as pore-filling cement (Figures 6(i), 7(c)7(e), 8(a), 8(b), and 9(c)). Chlorite coats represent the growth outward from detrital grain surfaces and are absent at points of grain-to-grain contact. Chlorite coats underlie quartz overgrowths, analcime, authigenic albite, barite, and calcite (Figures 7(c)7(e), 8(a), 8(b) and 9(c)). In addition, minor amounts of illite/smectite are present (Figures 9(d) and 9(e)), particularly in the fluvial channel facies.

There is large variation in the volume of analcime, ranging from 0 to . Analcime cements are the most abundant in the fluvial channel facies. The analcime is almost completely extinct and isotropic under crossed-polarized light and shows mainly pore-filling cementation (Figures 7(d), 7(e), 8(e), and 9(e)). Analcime cements are underlain by chlorite coats and replaced by calcite (Figures 7(d)7(g)) and are intimately intergrown with authigenic albite and illite/smectite clays (Figures 9(d) and 9(e)).

Anhydrite contents are very variable with a maximum of . Barite is rare, however, in some samples; up to barite is present. Sulfate cements fill pores in patches and locally replace detrital grains (Figures 7(h), 8(b)8(d) and 9(f)). Both anhydrite and barite grow over and engulf chlorite coats, quartz overgrowths, feldspar overgrowths, or authigenic albite (Figures 7(i), 7(j), 8(b), 8(c), and 9(f)) and are commonly juxtaposed with some of calcite cements (Figures 7(k), 8(c), and 8(d)).

Authigenic quartz has low contents, generally less than . Quartz cements occur as overgrowths on detrital quartz grains (Figures 7(c), 7(i), 8(a), and 9(g)), underlain by chlorite coats (Figure 7(c)) and enclosed by calcite and sulfate cements (Figures 7(i) and 8(a)). The contents of albite cements are less than . Albite pore fills and overgrowths occur. Pore-filling albite crystals are euhedral tabular and prismatic-shaped and are underlain by chlorite coats and hematite rims and engulfed by calcite and sulfate cements (Figures 8(b), 8(c), 9(c), 9(d), and 9(f). The albite crystals, in some samples, have been observed to be closely associated with analcime and illite/smectite clays (Figure 9(d)).

Hematite occurs in trace amounts as very thin rims around detrital grains (Figures 7(l) and 8(c)). The grain-riming hematite is generally absent at points of grain-to-grain contact, indicating that it is not syndepositional precipitation. Precipitation of the iron oxides likely occurred from general phreatic oxidation during the late Jurassic uplift [36].

4.4. Pore Bitumens and Fluid Inclusions

Oil shows are localized in cores, indicating the heterogeneous nature of oil migration in a reservoir. Ductile lithic-rich sandstones and tightly calcite-cemented sandstones do not contain oil shows in cores, pore bitumens, or oil-filled fluid inclusions in thin sections. This suggests that oil charge never occurred in these sandstones. By contrast, permeable reservoir sandstones generally have oil shows in cores. Pore bitumens and oil-filled inclusions are present, indicating that oil charged the sandstones. Therefore, permeable reservoir sandstones have been targeted for study of pore bitumens and fluid inclusions.

4.4.1. Petrography

Abundant pore-free oils and solid bitumens have been found. In contrast to nonfluorescent solid bitumens, pore-free oils have yellow and blue-white fluorescence under UV light. The yellow fluorescence dominates. These pore bitumens mostly fill pores and rare feldspar cleavage planes (Figures 10(a)10(h)). Oil-stained chlorite coats display a yellow fluorescence (Figures 10(a) and 10(b)), indicating that oil charge postdated the precipitation of chlorite coats.

With respect to oil-filled inclusions, two fluorescence colors, yellow and blue-white, have been observed. Similar to pore bitumens, the dominant fluorescence color is yellow. Oil inclusions are subspherical and elongated in shape, and some irregularly shaped varieties also occur. Oil inclusions are mostly two-phase inclusions containing less than 30 vol.% vapor at 20°C. Three-phase inclusions that contain a hydrocarbon liquid phase (or nonfluorescent solid bitumen), an aqueous liquid phase, and a vapor phase are present occasionally. Most oil inclusions occur within healed microfractures confined to individual quartz grains (Figures 10(i)10(l)). These oil inclusions generally range in diameter from 3 to 5 μm. In such a case, it is readily easy to find the coexisting aqueous inclusions. A small proportion of oil inclusions also lie in cements (Figures 10(m) and 10(n)), indicating that they were trapped at the time of the precipitation or recrystallization. No oil inclusions have been observed in syntaxial quartz overgrowths, suggesting that quartz precipitated prior to oil arrival. Oil inclusions and associated aqueous inclusions in calcite cements are smaller than within healed microfractures, often less than 3 μm in diameter.

4.4.2. Fluid Inclusion Microthermometry

Oil inclusions and aqueous inclusions in calcite cements are too small to reliably measure their homogenization temperatures, and some of the inclusions appear to have stretched or decrepitated. The aqueous inclusions that are coeval with oil inclusions hosted within quartz microfractures were selected for homogenization temperature measurement. A total of seventy-two homogenization temperature measurements were obtained from twenty-four fluid inclusion assemblages (Table 3). The homogenization temperature values display a bimodal distribution pattern: one population ranges from 80 to 100°C, and the other falls between 110 and 120°C (Figure 11). Noteworthy is the fact that the homogenization temperature values of the aqueous inclusions that are coeval with blue-white fluorescing oil inclusions are marginally lower than those associated with yellow fluorescing ones.

5. Discussion

5.1. Differential Diagenetic Sequences of Sandstones

Fluid inclusion data indicate that two episodes of oil charge took place in permeable reservoir sandstones in the Toutunhe Formation. Different episodes of oil emplacement, here, can be used to distinguish between preoil minerals that formed prior to oil charge, oil-stage minerals whose formation either resulted directly from the charge process or overlapped at least partly in time, and postoil minerals that precipitated after oil arrival in the reservoir. In terms of hydrocarbon migration and accumulation, such a partition is of great significance for understanding the temporal and spatial characteristics of rock-fluid interactions and assessing reservoir quality during the critical period of hydrocarbon charge [17, 18].

Based on textural relationships between authigenic minerals themselves and pore bitumens or oil inclusions, a general paragenetic sequence is presented in Figure 12. The dominant diagenetic features before the first oil charge in permeable reservoir sandstones were mechanical compaction, dissolution of detrital grains, and the formation of hematite rims, chlorite coats, illite/smectite clay, analcime, authigenic albite, and quartz overgrowths. The earliest formed minerals were chlorite coats and hematite rims, which were followed by analcime, authigenic albite, and quartz overgrowths. Most of these diagenetic alterations were controlled by the interaction of subaerial meteoric water with unstable volcanic rock fragments and detrital feldspars. Synchronous with the first oil arrival, quartz cementation was retarded. Anhydrite, barite, and calcite precipitated. When late oil charge occurred, calcite cementation continued. These reactions were controlled by high temperatures and increasingly complicated pore water chemistry arising from the interactions between hydrocarbon fluids and rocks as burial progressed.

In contrast to permeable reservoir sandstones, the compaction of easily deformed, clay-rich grains and associated authigenic clays dominated ductile lithic-rich sandstones before the first oil emplacement. Trace amounts of pyrite precipitated at this period. The ductile compaction continued but reduced with increasing burial. Almost no fluids entered the sandstones later. For tightly calcite-cemented sandstones, calcite cements significantly filled and occluded intergranular pores prior to the first oil arrival. Similar to ductile lithic-rich sandstones, no other diagenetic reactions occurred within the sandstones later.

5.2. Texture-Dependent Mechanical Compaction

The mechanical compaction of lithic sandstones is controlled by the amount and type of lithic fragments [11, 37, 38]. Volcanic rock fragments are the dominant framework component of sandstones in the Toutunhe Formation and are commonly altered. Breakdown of volcanic rock fragments initiated probably during transportation or weathering [38, 39] and continued during burial. When altered by weathering or diagenesis, volcanic rock fragments are extremely ductile [37]. Therefore, the altered volcanic rock fragments are the most important contributor to the mechanical compaction of the sandstones. Upon burial, these ductile grains were deformed between rigid detrital grains. They became smeared over the surfaces of the rigid grains and were squeezed into pore throats leading to significant porosity loss (Figures 6(a)6(d)). A moderate inverse correlation between porosity and permeability and ductile grain content is presented after excluding the tightly calcite-cemented sandstone samples from the plots in Figure 13, indicating that the ductile compaction is a major cause of porosity loss in the sandstones. Ductile lithic-rich sandstones containing ductile grains in excess of have an average porosity of 5%, ranging from 3.5 to 6.3%, and an average permeability of 0.10 mD, ranging from 0.04 to 0.18 mD. In comparison, permeable reservoir sandstones with less than ductile grain content have an average porosity of 13%, ranging from 8.6 to 15.9%, and an average permeability of 2.76 mD, ranging from 0.57 to 11.50 mD. A threshold value of 20% abundance seems to apply to these ductile grains.

Notably, it is clear from Figure 14 that the ductile grain abundance has a rough negative correlation with grain size. Most of the finer-grained sandstones (ripple cross-laminated) tend to have greater amounts of ductile grains and hence lower porosity and poorer reservoir properties. This suggests that depositional facies and lithofacies were important for predicting ductile components [11, 15, 38].

5.3. Origins of Authigenic Minerals and Their Effects on Reservoir Quality
5.3.1. Chlorite

Chlorite precursor phases and source and transportation of Fe2+ irons, almost entirely controlled by paleoclimate, paleosoil, and river system development, significantly affect the formation of chlorite coats [40]. Volcanic rock fragments are the most likely source of Fe2+ and Mg2+ ions for chlorite coats in the Toutunhe Formation. Five out of seventeen permeable reservoir sandstone samples contain chlorite coats (Table 1), so, while not frequent, they indicate that chlorite coats nonuniformly precipitate in sandstone.

Chlorite coats on quartz grains can prevent pervasive cementation of quartz overgrowths by limiting nucleation sites and inhibit the occlusion of primary intergranular porosity in sandstones [41, 42]. Figure 7(c) illustrates that minor breaks in the continuity of chlorite coats on quartz grains allow precipitation of quartz overgrowth. This is a compelling argument for the preservation of primary porosity by chlorite coats. The sandstones that contain chlorite coats are expected to have relatively high porosity values with a given composition, texture, and burial/temperature. The role of chlorite coats in the preservation of primary porosity, however, is not that important in the Toutunhe Formation. Clay grain coats have no effect on the precipitation of cements other than quartz [42]. Because of the mineralogically immature nature of the sandstones, quartz precipitation is almost volumetrically negligible in the Toutunhe Formation. Therefore, chlorite coats have little impact on primary porosity preservation.

5.3.2. Analcime

The mechanisms reported in the literature which may cause the formation of authigenic analcime in sandstones include direct precipitation from pore water or lake water, reaction of saline and alkaline pore water with detrital silicates, formation from a gel, and alteration of precursor zeolites [43]. In the Toutunhe Formation, no evidence for a precursor zeolite other than analcime has been found. The precipitation of analcime from sodium aluminosilicate gels is generally considered to have been associated with hydrothermal alterations [44]. No hydrothermal gels or other types of gels have been observed. Analcime cements in sandstones have euhedral, pore-filling, or blocky morphologies, indicating direct precipitation [45]. However, the nonuniform, localized occurrence of the analcime suggests that such a mechanism was not the principal means by which they formed. The reaction of volcanic rock detritus and associated clay derivatives with saline and alkaline pore fluid possibly produced some of the cements [4648]. In some localities, the analcime cements are closely associated with poorly crystalline clays, possibly some illite/smectite clays (Figures 9(d) and 9(e)).

5.3.3. Calcite

The δ18O values of calcite in the tightly calcite-cemented sandstones are consistent with the precipitation of early calcite from meteoric water at low temperatures. The negative δ13C values are probably related to a significant carbon contribution from soil-derived bicarbonate or oxidation of organic matter during early diagenesis. The cementation of early calcite, especially in completely cemented sandstones, causes substantial porosity loss. As discussed previously, the calcite forms discrete and tightly cemented concretions in the relatively coarse-grained part of a graded sand bed. This is consistent with an outcrop study from Jiao et al. [49] which indicated that the early calcite concretions occur commonly in the relatively coarse-grained sandstones at the central part of a fining-upward delta distributary channel sequence. They indicated that these concretions, in some localities, are up to 4 m in size and are heterogeneously distributed.

With respect to late calcite cements, the negative δ13C signature suggests that they were related to the supply of CO2 derived from decarboxylation of organic matter and Ca2+ ions in the formation water. The lower δ18O values are related to increasing temperatures and fluid-rock interactions during burial diagenesis. The presence of oil inclusions (Figures 10(m) and 10(n)) suggests that, during the formation of the late calcite, the ratio of hydrocarbons to aqueous fluids was high. In addition, the late calcite cements show slight enrichment in Mn2+ ions (Figure 6(d)). Cao et al. [50] proposed that Mn2+ is most likely derived from volcanic materials in the Junggar Basin and that the enrichment of Mn2+ ions in the pore water is closely associated with the dissolution of volcanic materials by “hot” hydrocarbon fluids in the deep Permian source rocks rather than volcanic fragments themselves in shallow sandstones. Therefore, we conclude that some of the late calcite precipitated from a Mn2+-rich hydrocarbon fluid that was introduced from depths into the Toutunhe Formation. It is fair to say that an episode of oil charge, ultimately sourced from the Permian source facies in the deep parts of the Fukang Sag, occurred in the Toutunhe Formation. This is entirely consistent with the analyses of biomarker compositions and geochemical origin of crude oils in the Toutunhe Formation by Hao et al. [51] who indicated that the oils were derived from the Permian and Jurassic source facies in the Toutunhe Formation.

5.3.4. Anhydrite and Barite

A source of and Ba2+ ions is required for sulphate cements. These are considered to be outside the Toutunhe Formation. The Permian source rocks contain thin seams of gypsum [52]. These gypsum deposits can spatially provide ions for the sulphate cements. A possible source for Ba2+ ions is the underlying Middle-Lower Jurassic coals [53, 54]. Both anhydrite and barite precipitated during late diagenesis, possibly temporally responding to the northeast-southwest trending tilting of the Fukang Sag due to the uplift of the Bogda Mountains since the Paleogene [24]. We tentatively suggest that a great proportion of deep-seated hydrocarbons that were derived from the Permian source facies possibly migrated up along faults from depths and leaked into the Toutunhe Formation. Synchronous with the up-fault migration of deep hydrocarbons, and Ba2+ invaded the Toutunhe Formation and were precipitated as barite and anhydrite. Although no oil inclusions have been found in the sulphate cements because of limited thin sections for fluid inclusion study, petrographic evidence indicates that barite and anhydrite cements are possibly simultaneous with late calcite, which contain oil inclusions. Thus the precipitation of the sulphate cements may be lined to oil charge. This further confirms that there was once an episode of charge in the Toutunhe Formation, ultimately sourced from the Permian source facies. The preferential occurrence of the sulphate cements in relatively coarse-grained sandstones may be due to the easier movement of fluids through these sandstones. These sulphate cements generally reduce porosity and permeability but do not completely prevent late flow as they occur in localized patchy areas similar to late calcite.

5.4. Timing of Secondary Dissolution

Determining the timing of secondary dissolution is difficult. Detrital feldspars and volcanic rock fragments are very unstable during near-surface meteoric water flushing [55]. However, the early diagenetic analcime cements, in some samples, have been observed to be partly dissolved (Figures 7(d)7(g)). Additionally, relic chlorite coats represent delicate textures (Figure 9(c)), suggesting late dissolution. This is because such fabric, if formed relatively early, should have been destroyed or reoccupied by later minerals [56]. Therefore, secondary dissolution most likely continued throughout burial.

A geochemically nearly closed system for burial diagenesis has been accepted [42, 57]. In such a system below the reach of freshwater, the pore water flow is limited and the bulk chemical composition of the sediments does not change significantly. The idea that significant net increases in porosity may occur during burial diagenesis due to leaching of framework grains by organic acids has mostly been problematic [42, 57]. Furthermore, organic acid quantities are insufficient to significantly influence bulk rock-water equilibria, albeit being present extensively in oilfield formation water [58, 59]. Organic acids may also be neutralized during expelling from source rocks and migration in sandstones [57]. Therefore, we concluded that most of secondary porosity formed by meteoric water influx through the sandstones at shallow burials.

5.5. Effect of Diagenetic Heterogeneity on Fluid Flow

The two sandstone petrofacies, highly compacted, ductile lithic-rich, very fine-grained sandstones, and tightly calcite-cemented sandstones may create permeability barriers or baffles embedded in porous reservoir sandstones. For reasons discussed previously, these permeability barriers are of early diagenetic origin. Depositional environment and lithofacies and the pore water chemistry under near-surface conditions are all key parts in their formation and distribution. These barrier interbeds vary widely in the abundance and occurrence. If they are continuous over large areas in a reservoir, they have a significant effect on fluid flow [6, 19, 60].

Due to reservoir compartmentalization by the intersecting barrier interbeds, oil migration and accumulation in a reservoir are inevitably hindered. An invasion-percolation-based model of oil migration by Luo et al. [18, 19] indicated that oils migrate laterally through a larger fraction of a heterogeneous reservoir than a homogeneous one. Luo et al. [18, 19] thought that not all permeable reservoir compartments separated by the intersecting barrier interbeds act as the pathways for migration or trap for accumulation and a considerable portion of them are filled with formation water. Although the oil-water contact in individual permeable compartments is essentially horizontal, the blocking effect of the intersecting barrier interbeds causes the differences in the oil-water contact between compartments, resulting in an uneven oil-water contact at a reservoir scale. In the Toutunhe Formation, well testing shows that oil and water zones are commonly interbedded; there exists a complex oil-water relationship. Reservoir heterogeneity most likely gives a good explanation for the complex oil-water distribution in the Toutunhe Formation.

In addition, the intra-sand-body barrier interbeds may significantly increase the tortuosity of fluid flow and alter the sweep efficiency, affecting recovery during waterflooding [6, 60].

The study concentrated on reservoir heterogeneity at the core scale. The combination of outcrop study and subsurface data will be used to predict the regional distribution of both the ductile lithic-rich sandstones and tightly calcite-cemented sandstones and thus to characterize reservoir heterogeneity at the macroscale and higher hierarchies. This necessitates single-well heterogeneity characterization and interwell correlation of the area of interest, especially when cores are strictly limited for deep drilling. Therefore, well log responses to these barrier interbeds can be quantified at first. Our future work is to construct regional barrier interbed networks and modeling oil migration.

6. Conclusions

We emphasize that the evaluation of sandstone diagenetic heterogeneity requires systematic examinations of depositional environment and lithofacies that basically determine the primary texture and compositions of sandstones.

Five depositional facies associations of the sandstones have been distinguished based on core sections and well log data: single and stacked story meandering fluvial channel, delta distributary channel, distributary-mouth bar, interdistributary bay, and shore-shallow lacustrine facies association.

The Toutunhe Formation’s sandstones consist mostly of feldspathic litharenite and litharenite. Grain size and the abundance of ductile-lithic sand grains are the most important primary detrital controls on diagenesis, porosity, and reservoir properties. These controls operate largely through their effect on mechanical compaction.

Three major sandstone petrofacies have been defined in the Toutunhe Formation. Ductile lithic-rich, very fine-grained sandstones, which occur in all facies associations, are dominated by compaction of easily deformed, clay-rich grains, resulting in a very rapid loss of porosity during burial. Tightly calcite-cemented sandstones occur commonly in the relatively coarse-grained lithofacies of a graded sand bed. By contrast, fine-grained permeable reservoir sandstones, deposited predominately in the fluvial channel, distributary channel, and distributary-mouth bar facies associations, contain relatively less abundant ductile-lithic grains. Various authigenic minerals are present in the relatively coarser-grained sandstone lithofacies with predominately a patchy texture. Chlorite coats only occur in the distributary channel and distributary-mouth bar facies. Analcime and illite/smectite clays are commonly found in the fluvial channel facies. Diagenesis progressed alternately with oil charge, and some diagenetic alterations and oil charge occurred simultaneously. Late diagenetic calcite and sulphate cements may have an indication for oil migration.

Highly compacted, ductile lithic-rich sandstones and tightly calcite-cemented sandstones can create permeability barriers embedded in a porous reservoir, resulting in heterogeneous flow.


The authors accept full responsibility for the ideas and conclusions expressed in the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (no. XDA14010202) and by the National Natural Science Foundation of China (no. 41372151). The authors are indebted to Mr. Tao Xu at the Shengli Oilfield Company, SINOPIC, for helpful suggestions.