Origin of Bedding-Parallel Calcite Veins from Lacustrine Shale in the Eocene Dongying Depression, Bohai Bay Basin, China

Calcite veins, which developed parallel to the bedding, are widespread in laminated source rocks in the Eocene Dongying Depression. However, there is a lack of systematic description and classification of the veins. This study presents a systematic characterization of the calcite veins, host rocks, and micritic carbonate laminae by applying petrographic and geochemical methods to understand vein-forming mechanisms. Antitaxial and syntaxial veins are examined. Antitaxial veins contain typical fibrous crystals with the most intense fluorescence, and the median zone of these veins is often the micritic carbonate. Calcite crystals in syntaxial veins develop a blocky morphology of various sizes, indicating obvious growth competition. Data of rare earth elements and trace elements obtained from the micritic laminae, host rocks, and calcite veins are very similar. This indicates that the vein-forming nutrients originated from the carbonate in the host rocks and micritic laminae. The minor difference in C and Sr isotopes between calcite veins and micritic carbonate within the host rock and the negative shift in O isotopes in the veins are caused by ion exchange and dehydration of swelling clay minerals in the burial environment. This further proves that the calcite veins are formed in a closed system. Geochemical analysis suggests that the rocks are in the oil window and have good hydrocarbon potential. Thermal evolution of the acidic fluids generated from organic matter (OM) resulted in the dissolution of carbonate and formed fluid overpressure in the rocks. Fluid overpressure induced the formation of fractures in the interlayer and expanded the veins with the force of crystallization due to fibrous calcite growth. Blocky crystals grow in the fractures from the margins toward the center. Hydrocarbon expulsed via OM maturation in the host rock fills the intercrystalline pores. Moreover, shale with bedding-parallel calcite has the characteristics of high-quality shale oil reservoirs. These characteristics will probably provide guidance for shale oil exploration.


Introduction
Calcite veins are common in low-permeability source rocks. In particular, bedding-parallel veins filled by fibrous calcite, in which the fibers are mutually parallel and form quasivertically, have been reported from more than 110 localities in the world [1][2][3]. The shape, orientation, and internal structure of calcite veins preserve abundant information on the diagenesis, palaeostress, and rock deformation [4,5]. According to various vein types and associated crystal growth direction, morphologies, and vein growth mechanisms, Durney and Ramsay [6] distinguished stretched, syntaxial, and antitaxial veins. They can be further subdivided into blocky, elongate-blocky, and fibrous veins based on their crystal forms and arrangements [7].
Syntaxial veins often contain blocky or elongate-blocky crystals, which are ascribed to crystal growth into open holes or rapidly opening veins and exhibit strong growth competition. Since the shape of crystals cannot define the opening trail, only the total comparative displacement of the vein walls can be determined [7]. Antitaxial veins normally consist of fibrous crystals with smooth crystal boundaries [5]. Curved fiber boundaries can be used to infer the opening trail during the growth of cracks [8]. In addition, both solid inclusions detached from the vein walls and different phases of vein growth can indicate the opening process or direction [9,10]. Ramsay [9] proposed that the crack-seal mechanism can explain the growth of veins with elongated to fibrous crystals. Numerical simulations of fiber growth based on the crack-seal mechanism revealed that crystal growth competition is suppressed when the fracture width is less than 10 μm [11]. Therefore, the genesis of veins has been attributed to fluid overpressure caused by ascending fluid migration, which results in tensile stress in low-permeability host rocks [12,13]. The force of crystallization and ongoing deformation can explain fibrous veins with continuous growth characteristics [14][15][16][17]. However, not all fibrous veins form by the crack-seal mechanism. Some may form without the loss of cohesion between vein fibers and host rock [15], such as solid inclusions formed by fluctuations of adhesion at the interface of the vein and host rock [17].
Secondary minerals in veins precipitate mainly from supersaturated aqueous solutions because of changes in physiochemical conditions [18]. Diffusion and advection are the fundamental material-transport mode that controls vein formation [19]. In an open system, advective flow carries the solutions over long distances until they precipitate minerals in veins [20]. In contrast, in a closed system, the solutions diffuse over shorter distances [21]. Nevertheless, diffusion always plays a vital role in the formation of relatively thicker veins [22]. The mode of transport can be inferred from the geometry of the microstructures in veins. Therefore, the mineralogical and geochemical analysis of veins can provide valuable clues regarding the origin and nature of fluids.
The Dongying Depression in eastern China is rich in oilbearing shales in which several bedding-parallel calcite veins are preserved [23]. Such veins are regarded to have favored the preservation of oil in the Dongying Depression [24,25]. However, there is a lack of systematic description and classification of the veins. Further, there is no consensus regarding the genesis of different calcite veins. The formation of calcite veins is linked to the timing of hydrocarbon generation, and fluid overpressure is considered a critical factor for the occurrence of natural fractures in shale [3,23]. Organic acids and by-products formed during oil generation promote the dissolution of micrite, providing ions for calcite precipitation in veins [3]. Wang et al. [26] suggested that the migration of hydrocarbon-bearing fluids and consequent mineral precipitation led to continuous vein growth and expansion.
In this study, the calcite veins in the Dongying Depression were systematically characterized using petrographic and geochemical methods. The mechanism and fluid sources responsible for vein formation were discussed. The results of this study provide a better understanding of the impact of vein formation on shale oil reservoirs.

Geological Setting
The Bohai Bay Basin, located in eastern China, is an important hydrocarbon-producing basin. It comprises several subbasins and covers approximately 200,000 km 2 (Figure 1(a)). The Jiyang Subbasin is one of the subbasins bounded by the Tan-Lu Fault to the east, Chengning Uplift to the northwest, and Luxi Uplift to the south. The study area, i.e., Dongying Depression, is a secondary tectonic unit of the Jiyang Subbasin and has an area of around 5,700 km 2 (Figure 1(b)) [27]. It is situated in the south-eastern part of the Jiyang Subbasin and is bordered by the Luxi Uplift to the south, Guangrao rise to the southeast, Chengjiazhuang rise to the north, and Binxian and Qingcheng rises to the west. The study region comprises the south-western Boxin, central Niuzhuang, northwestern Lijing, and north-eastern Minfeng sags and the central diapiric anticline (Figure 1(c)).
The tectonic evolution of the Dongying Depression in the initial syn-rifting stage can be divided into rifting I to IV ( Figure 2). During rifting II and III, extensive lacustrine expansion occurred, resulting in a broad semideep lacustrine and deposition of abundant fine-grained sedimentary rocks in the Es4s-Es3x interval [27]. The Dongying Depression is filled with Cenozoic sediments forming six formations, namely, the Paleogene Kongdian (Ek), Shahejie (Es), and Dongying (Ed) Formations; Neogene Guantao (Ng) and Minghuazhen (Nm) Formations; and the Quaternary Pingyuan (Qp) Formation ( Figure 2). The Es Formation is divided into four members as Es1, Es2, Es3, and Es4 from the top to the base. The Es4 and Es3 members can be further divided into Es4x, Es4s, Es3x, Es3z, and Es3s from the base to the top.
During the accumulation of the Es4x, sedimentation occurred in a dry and hot climate, and four sags were relatively isolated. The sags subsequently merged into a shallow lake in the Es4s. Under comparatively warm-damp conditions, salt lacustrine facies comprising gypsum, argillaceous carbonate, calcareous mudstone, and shale were deposited ( Figure 2) [28]. During the accumulation of the Es3x, the Dongying Depression was intensely faulted. The basement rocks rapidly subsided to form a relatively sustained deeplacustrine environment with the lacustrine sediments comprising mudstone, calcareous mudstone, and shale. The Es3z and Es3s lithofacies are similar to the Es3x lithofacies and are characterized by stable sedimentation in the lacustrine basin that resulted in calcareous mudstone, shale, etc. (Figure 2) [28]. The interval of interest for this research is Es4s and Es3x, in which calcareous mudstone and shale with many bedding-parallel calcite veins are preserved.

Materials and Methods
The Es3x-Es4s interval is penetrated by the NY1 and FY1 wells in the Dongying Depression. A total of 17 core samples were collected from the NY1 and FY1 wells for thin sections (Table 1). Samples with bedding-parallel calcite veins (n = 9) were trimmed into cubes and polished. The calcite veins and host rocks were then drilled from the cubes using a microscope-mounted drilling assembly (1 mm diameter drill bit; ×10 magnification). Similarly, samples of micritic carbonate laminae (n = 8) were collected and prepared.
Cathodoluminescence (CL) observation of thin sections was performed under a Leica DM2500P CL microscope with a voltage of 12 kV and current of 250 μA. The thin sections were coated by a Pt conducting layer and then observed under a ZEISS S-3400 N field emission scanning electron microscope (SEM) with ETH of 15-20 kV. All the 2 Geofluids microscopic analyses were performed at the State Key Laboratory of Geological Processes and Mineral Resources in Beijing. An X-ray diffraction (XRD) analysis of the powdered samples was conducted by using a D/max-2500 theta/theta rotating anode X-ray diffractometer. The total organic carbon (TOC) of bulk rocks was determined using a carbon sulfur analyzer, LECO CS600. The measurement technique used here is based on sample combustion in an O 2 atmosphere to convert TOC to CO 2 . Thermal maturities were determined    4 Geofluids using a Leitz MVP-3 microscope photometer, which helps measure vitrinite reflectance (Ro) values. Asphalt was extracted before pyrolysis. The Rock-Eval parameters measured include the free hydrocarbons (S1), hydrocarbons cracked from kerogen (S2), carbon dioxide relieved from organic matter (S3), and maximum yield temperature of pyrolysate (Tmax). The microanalysis of TOC was conducted using a LECO CS600 at the Laboratory of Geological Microbiology of the China University of Geosciences, Beijing. The microdrilling samples of micritic carbonate lami-nae, veins, and host rocks were analyzed for trace elements, rare earth elements (REEs), carbon and oxygen isotopes, organic carbon isotopes, and strontium isotopes. These analyses were conducted at the Analytical Laboratory of the Beijing Research Institute of Uranium Geology. Trace elements and REEs were measured using a NexION300D apparatus for inductively coupled plasma-mass spectrometry (ICP-MS; >10 ppm, error < 5%; <10 ppm, error < 10%). Carbon and oxygen isotopes were measured by adding concentrated phosphoric acids to produce the CO 2 which were tested by a MAT-253 gas isotope mass spectrometer. Carbon and oxygen isotope compositions are reported in the standard delta notation as permil (‰) relative to Vienna Pee Dee Belemnite (VPDB), with an analytical precision of ±0.2‰ for carbon isotopes and ±0.1‰ for oxygen isotopes. Organic carbon isotope was measured using a MAT-253 gas isotope mass spectrometer when burning a sample in flowing oxygen to produce the CO 2 . Strontium isotope was determined by using a phoenix thermal ionization mass spectrometer (TIMS). The powder samples (0.1 to 0.5 g) were digested in 5 mL HF and 2 mL HClO 4 . Strontium was separated in 5 mL HCl using 200 to 400 mesh cation exchange resin. During mass spectrometry measurement, the Sr isotope ratio was obtained at the 1250 ± 50°C and an ion current intensity of 1 × 10 −11 A in vacuum. The 87 Sr/ 86 Sr ratio was corrected for mass fractionation using 88 Sr/ 86 Sr = 8:37521.

Petrologic Features.
Observations at NY1 and FY1 wells suggest that the lithologies of Es4s and Es3x mainly consist of calcareous mudstone and shale, sandy mudstone, gypsum mudstone, and argillaceous limestone [25]. The carbonate minerals are interbedded in the mudstone or shale and occur    Figure 3. Collectively, they account for more than 90% of the bulk samples. The data of clay ( Table 2) show that the content in the NY1 well ranges from 4% to 59%, with an average of 24.1%. The average feldspar and quartz contents are 5% and 22.1%, respectively. The average calcite and dolomite contents are 32.8% and 13.3%, respectively. The clay mineral content in the FY1 well ranges from 2% to 57%, with an average of 17.8%. The average feldspar and quartz contents are 4.6% and 24.9%, respectively. The average calcite and dolomite contents are 38.3% and 11.9%, respectively.
Statistical analysis of the clay mineral content of the NY1 and FY1 wells and some previously published data from the study area indicate that the clay minerals include kaolinite, mixed illite/smectite, illite, and chlorite ( Figure 4) [30,31]. The illite content increases significantly with increasing depth, whereas the mixed illite/smectite content shows the opposite trend. The kaolinite content also displays a negative correlation with the depth. The chlorite content is less than 8%.  Fy1 data Chen's data Li's data Ny1 data Figure 4: Proportions of clay minerals in samples from Es3 to Es4. Data are from the NY1 and FY1 wells; the data were collected for this study and from some previously published work on the study area [30,31]. 6 Geofluids In the CL analysis, antitaxial veins with fibrous crystals exhibit a bright red and dully red luminescence (Figure 7(a)). The edge of the lensoid calcite veins exhibit slightly brighter luminescence (Figure 7(b)). Syntaxial calcite veins contain crystals with coarse-crystalline, blocky morphologies. The contact surface of the vein and host rock is smooth. The crystal size varies from several to a hundred microns. The dark solid inclusions often fill between the blocky crystal, creating irregular shapes (Figures 6(e), 6(g), and 6(h)). The syntaxial veins with blocky crystals show oscillatory CL zoning with the exterior of crystals being brighter than the interior. In addition, fibrous calcite veins are observed near the blocky calcite veins, ranging in size from tens to hundreds of microns (

Calcite
Veins. The solid inclusions in fibrous and blocky calcite veins were identified as clay minerals, quartz, pyrite framboids, and organic matter (OM) by SEM ( Figure 9). The pyrites are widely distributed in the shale and the solid inclusions between lenticular calcite veins (Figures 9(c) and 9(d)) and are also present in the center of the veins (Figure 9 (Figure 9(h)) and among irregularly shaped blocky crystals (Figure 9(a)). Barite and strontianite can be observed by a qualitative analysis by SEM-energy dispersive X-ray spectroscopy (EDS). Barites can be observed at the margin of fibrous calcite veins and are partially replaced by calcite (Figure 9(e)). Strontianites have irregular shapes and fill in calcite veins (Figures 9(g)) and are very common in the fibrous calcite vein.

Organic
Matter. The TOC contents of the bulk rock in the NY1 (188 samples) and FY1 (143 samples) wells are 0.15-9.32 wt. % and 0.11-9.05 wt. %, respectively, with average values of 3.00 wt. % and 2.41 wt. %, respectively. Samples with TOC greater than 1 wt. % account for 93.8% and 98.5% of the total samples analyzed from NY1 and FY1 well, respectively. The S1 + S2 values of the bulk samples are approximately 15 mg/g. The plot of HI versus the pyrolysis Tmax indicates principally type I kerogen [27]. The vitrinite reflectance (Ro) values for NY1 and FY1 are 0.52-0.80% and 0.51-0.91%, respectively ( Figure 10). Results of the microanalysis of the TOC content of the host rock are listed in Table 3

Formation Mechanism of Fractures.
Fractures provide the space for the precipitation of calcite in the veins and are thereby preconditions for the formation of the latter. The calcite veins normally show bedding-parallel features. This indicates that the internal stress of shale, rather than external stress (tectonic activity) [32,33], is the main factor that the precipitation of the calcite vein is accommodated. The presence of OM in the calcite veins indicates that hydrocarbonbearing fluids flowed through the veins and filled them since vein opening. The rocks in the NY1 and FY1 wells are characterized by the high carbonate content, low feldspar content, and a TOC of almost more than 2 wt. %. Moreover, the rocks are in the oil window (Ro > 0:5%) and have good hydrocarbon potential (type I kerogen and high S1 + S2 values) ( Figure 3; Figure 10). Wang et al. [26] reported that the fluid inclusions in calcite veins normally have yellow to blue-green colors under ultra-violet (UV) illumination. The trapping temperature exceeded 70°C because of the homogenization temperatures of inclusions. Thus, the thermal evolution of the OM probably generated the acidic fluids that dissolved the carbonate in the rocks. Previous research has also suggested that the fluids were generated as the OM became matured, resulting in an increase in the local volume with abnormal overpressure within the Es3 and Es4 intervals, rather than by disequilibrium compaction [27]. In addition, with increasing burial depth, the mixed illite/smectite content decreases, and it transforms into illite (Figure 4). The dehydration of clay minerals causes structural modifications, such as the opening of pores and/or horizontal fractures (Figure 9(c)) [34], which are likely responsible for the fluid overpressure [26]. The fibrous crystals that grow almost perpendicular to the fracture wall of host rocks and have smooth fiber boundaries indicate that overpressure is not likely caused by horizontal compression. According to morphology changes with respect to the wall of fibrous calcite cemented veins, Rodrigues et al. [18] identified the opening of veins owing to horizontal compression, which shows an obvious angle change from inner zones to outer zones in the fibrous calcite. Since the shear stresses acted at the vein margins, the previously consolidated vein resisted shortening, whereas the host rock did not. Thus, the calcite veins formed by horizontal compression are expected to curve with respect to the wall; however, such curved veins are not present in our samples (Figures 6(a) and 6(c)). Changes in lithology during sediment accumulation are another potential factor for the formation of beddingparallel fractures. The main changes are a result of the growth of micritic carbonates in the median zone, with flanking fibrous crystals, or the growth of fibrous crystals at the margin of the blocky vein (Figures 6(c), 6(d), and 6(f)). The major rocks in the Es4s-Es3x intervals are laminated organic-rich shale, consisting of micritic and clay laminae [25]. Interlayer fractures often grow at the interface between micritic and clay laminae [34,35]. This is because the transitional surface is the weakest zone in the rock, allowing for fracturing by fluid pressure. Therefore, the secondary fractures promote the migration of hydrocarbon-bearing fluids and the consequent calcite precipitation.
In summary, given increasing depth, OM can generate acidic fluids that dissolve carbonate minerals, and the mixed illite/smectite is transformed to illite, thereby promoting the formation of spaces and fluid overpressure. Besides, fluid overpressure has a vital impact on forming bedding-parallel fractures in the transitional surface of the rock and interlayer within rocks.

Source of Vein Fill.
The increase in pore pressure may affect the chemical equilibrium of the aqueous system by increasing the solubility of carbonate minerals [5]. The OM is thermally mature. Thus, it expels organic acids and large volumes of CO 2 [36], leading to the dissolution of the primary carbonate in the host rock. As the pore pressure increases, an increasing amount of CaCO 3 is dissolved into fluids until the maximum is reached [37]. However, as the pore pressure decreases drastically because of horizontal microfracturing and pore opening, the supersaturated pore fluid causes carbonate precipitation in the fracture system. Both strontianite and barite observed within the calcite veins indicate that they precipitated simultaneously with the calcite. Strontianite is a rare low-temperature hydrothermal mineral. Unlike strontianite, barite has multiple origins and could form through biogenic and hydrothermal processes [38][39][40]. The coexistence of strontianite and barite indicates the possible influence of magmatic hydrothermal fluid [39].

Trace Elements and Rare Earth Elements (REEs).
Data of trace elements and REEs can potentially offer information on the source of minerals in veins and the fluid pathways [7]. Chondrite-normalized REEs of micritic laminae, fibrous and blocky calcite, and host rocks show similar patterns ( Figure 11). The δEu values for the micritic laminae and host rocks are similar to those from most of the calcite veins. The NY1-10 dataset has an anomalous δEu value of 1.57, which may be associated with the preferential leaching of Eu-rich material (for example, feldspar), reductive vein formation conditions, or transformation into divalent Eu as temperature increases [41]. The crystal shape of the calcite (fibrous and blocky) in the veins cannot be explained by the REE data    15 Geofluids [7]. The REE patterns in the veins are consistent with those in host rocks and micritic carbonate laminae, indicating that the materials within the calcite veins may have been derived from the host rock carbonate and micritic laminae.
The trace elements in the micritic laminae, host rocks and veins also exhibit very similar patterns. The concentration of strontium is generally high not only in micritic laminae and host rocks but also in the veins ( Figure 11); this is attributed to the dissolved micritic laminae and the host rock with a high Sr ion concentration flowing into diagenetic fluids. During such transfers, precipitation results in the formation of Sr-rich calcite in horizontal fractures. This may explain the high Sr content in the calcite veins [25,42] and the source of Sr 2+ for strontianite precipitation. The Ba content is 10456 μg/g in sample FY1-28 owing to the occurrence of barite within the vein (Figure 9(e)). The Es4s in the Dongying Depression develops laminated gypsum mudstone, which is abundant calcium and sulfate ions [24]. Pyrite is packed in the calcite vein (Figure 9(b)), indicating that the sulfate was partially removed during the sulfate reduction reaction [43]. Thus, Liang et al. [44] suggested that gypsum dissolution in the diagenetic fluid increased the concentration of sulfate and resulted in barite precipitation in the fractures. Another possibility is that Ba is universally present in the OM and can be replaced by potassium [45]; thus, sufficient barium may have been liberated from the host rock to facilitate the precipitation of barite.
Therefore, the similar patterns of trace elements and REEs in the micritic laminae, host rocks, and veins indicate that ions that leached from adjacent host rocks and micritic laminae are concentrated in nascent veins because of the focused fluid expulsion along these pathways, after the maturation of OM and the associated generation of organic and carbonic acids. The formation of strontianite and barite is attributed to the contribution of ions from the host rocks and micritic laminae, rather than to the flow of external fluids into the host rock.

C and O Isotopes.
The carbon and oxygen isotope compositions of the calcite veins represent the original isotope compositions and specific temperature of the parent fluid precipitation, respectively [5,46]. However, because of the resetting of δ 18 O values during diagenesis in sediments, the initial δ 18 O composition of shale with veins cannot be determined [47]. Therefore, the δ 18 O values in veins can provide clues regarding the source of the diagenetic fluids and the associated physicochemical conditions [48]. The δ 18 O values of most calcite veins are lower than those of the micritic carbonate laminae and carbonates in host rocks ( Figure 12); this observation is similar to that made for the coarse-grain calcite veins by Liu et al. [49].
Water-rock interaction during the temperature rise in the rock is likely the major reason for oxygen isotope depletion in the fluids, and thus, it probably influenced the δ 18 O of authigenic calcite veins [50,51]. The oxygen isotopic composition of formation water could have been modified gradually because of oxygen exchange within the silicates of the host rock [52]. Meteoric water infiltration is another mechanism for the negative shift in δ 18 O values in the calcite veins. Given the low permeability of the shale and the lack of tectonic uplift [49], the influence of meteoric water infiltration is likely small. Thus, the negative shift in δ 18 O more likely results from the modification of the parent fluid via the dehydration of swelling clay minerals and the increase in temperature. The presence of pyrite in a calcite vein may imply that sulfate reductions in the decomposition of OM produced a minor negative shift in δ 18 O. However, if sulfate reduction occurred in the study area, it would have led to a large This possibly suggests that the blocky veins formed earlier than the fibrous veins. The majority of the δ 13 C values of the fibrous and blocky veins are similar to those of the micritic carbonate laminae, with an average of +3.5‰ VPDB, except for sample NY1-13 that has the lowest δ 13 C value ( Table 4). The enrichment of 13 C in the host rock is induced by exchanges with the dissolved organic carbon [53]. Hence, the host rock carbonate has a lower δ 13 C value than the calcite vein; if the carbon isotope composition of the veins is the same as that of host rock carbonate, the incorporation of depleted organic carbon into the carbonate in the host rock can lead to this difference. The depleted organic carbon can be exchanged with any existing veins, but a larger crystal size in veins could have helped maintain their carbon composition, whereas the fine- grained carbonate in the host rocks could have been more thoroughly exchanged with the depleted carbon from the OM [53]. The δ 13 C org results show that the calcite veins are lighter than the host rocks (Figure 13(a)), which too indicates the exchange of organic carbon with carbonate in host rocks in the burial environment. Blocky veins are related to hydrocarbons, and it would have been easier for OM to fill up blocky crystals (Figures 6(e) and 6(g)); therefore, the organic carbon isotope of blocky veins is lower than that of fibrous veins.
Thus, the carbonate content of the calcite veins is derived from the micritic carbonate dissolution-reprecipitation mechanism. The negative shift of the δ 18 O value in calcite veins is mainly caused by the dehydration of swelling clay   Figure 13: Organic carbon (a) and strontium (b) isotopes of calcite veins and host rocks. 18 Geofluids minerals during the increase in the burial temperature of the shale. The host rock carbonate with lower δ 13 C values results from exchanges with the carbon taken from the OM.

Sr
Isotopes. The Sr isotope of calcite veins is another factor that indicates whether the veins grew from local or regional fluids [19]. 87 Sr is obtained from the radioactive decay of 87 Rb, which is not incorporated into the calcite lattice. Therefore, radioactive decay over time does not affect the Sr isotope of a calcite vein. If the calcite veins were derived from the carbonates in the host rocks, the 87 Sr/ 86 Sr ratio of the veins and host rock carbonate would be similar. Therefore, by comparing the Sr isotope values in calcite veins and host rocks, it is possible to determine whether the calcite veins formed from local materials or from the advective transport of fluids over a longer distance [20]. The Sr isotopic ratio in the host rock is slightly greater than that in its vein ( Figure 13(b)), but still much lower than that in its noncarbonate sections. Perhaps, small amounts of Sr within detrital silicate were exchanged with the carbonate in the host rock, thereby increasing the 87 Sr/ 86 Sr ratio [54]. The average Sr isotopic composition of micritic carbonate is 0.710718, which may reflect the characteristics of the initial sedimentary environment. The Sr isotopic concentration of the fibrous calcite vein is similar to that of micritic carbonate, indicating that the fluids of the fibrous vein were derived from the micritic carbonate. The isotopic concentration of the blocky calcite veins is also close to those of micritic carbonate or slightly higher than those of fibrous veins, indicating the Sr isotope may be affected by the dehydration of clay minerals in the early stage because the blocky calcite forms earlier than the fibrous calcite ( Figure 14). Therefore, the materials in the calcite vein originate from micritic carbonate in the host rock and micritic laminae.

Formation
Mode of Calcite Veins. The rocks from well NY1 and FY1 are in the oil window (Ro > 0:5) and have good hydrocarbon potential. The maturation of the OM leads to the expulsion of organic acid and hydrocarbons and the dehydration of clay minerals, resulting in the dissolution of micritic carbonate and fluid overpressure within the host rocks [3,27]. These diagenetic fluids within the host rocks migrate to the interlayer via diffusion, forming overpressure fractures and providing space for the later precipitation of calcite veins. The syntaxial veins with blocky crystals are induced by the different growth rates of the individual crystals; they may form in a cracking event, thereby creating enough space for the dissolved carbonate within the host rock to precipitate [7,55]. The fibrous calcite grows in a narrow fracture, and crystal growth competition is suppressed [11]. The vein-wall interface represents an area with low adhesion and new material accumulation, causing the dilation of the vein because of the force of crystallization of calcite fibers [17]. Moreover, the petrological observation and differences in the oxygen and strontium isotopic concentrations in blocky and fibrous calcite probably indicate that the blocky calcite forms earlier than the fibrous calcite. Therefore, the migration of the nutrients of the calcite within the host rocks to the interlayer fracture promotes the growth of blocky calcite. Generally, the blocky crystals comprise the clay minerals that detached from the vein walls and the hydrocarbons expulsed upon OM maturation. After the wide fracture is filled with the blocky vein, the vertical migration of fluids is hindered causing the dissolved carbonate to gradually precipitate into fibrous calcite in the vein-wall interface (Figure 15(a)). Besides, the hydrocarbon-bearing fluid contributes to the formation of bedding-parallel fractures in the transitional surface of micritic carbonate lamina and organic-rich clay lamina and interlayer within rocks, and this in turn results in the growth and continuous dilation of fibrous crystals (Figure 15(b)) [56].

Significance of Calcite Veins for Shale
Oil. The Es4s-Es3x rocks were deposited in a semideep and deep lacustrine environment. Compared with the marine shale in China, the TOC and carbonate mineral contents are considerably higher in the Es4s-Es3x rocks [57][58][59][60]. Moreover, compared with some major shale in the United States, the Es4s-Es3x shale has deep burial depth, low thermal maturity, high carbonate content, and low quartz content. These differences affect the types of pores in the Es4s-Es3x shale [59,61,62]. Therefore, unlike marine gas shale with primarily OM pores, the carbonate recrystallization intercrystalline pores and dissolution pores are the main pores in the Es4s-Es3x shale [24]. Moreover, the pores are mainly concentrated in the laminated interface near the veins; further, there exists favorable horizontal connectivity along a given lamina because of the connection of microfractures [49].
Therefore, the shale with a bedding-parallel calcite vein is considered a high-quality reservoir, characterized by high TOC and hydrocarbon potential, abundant recrystallization intercrystalline pores and interlaminar fractures, and high content of brittle minerals (average greater than 70%) [61,63]. The shale with veins is the favorable lithofacies for shale oil reservoirs and a sweet spot for exploration [64].

Conclusions
(1) Bedding-parallel calcite veins are widespread in source rocks in the Dongying Depression. Antitaxial and syntaxial calcite veins can be distinguished through petrographic observations. Antitaxial veins typically contain fibrous crystals that grow almost perpendicular to the host rocks with smooth fiber boundaries. The median zone of antitaxial veins is often initial micritic carbonate. Calcite in syntaxial veins develops blocky morphologies with various sizes and clear crystal growth competition (2) The shale mineralogy of NY1 and FY1 has high carbonate and low feldspar content. The shale has been buried almost in the oil window and has good hydrocarbon potential. With increasing depth, the OM generates acidic fluids to dissolve carbonate. The mixed illite/smectite also transforms into illite. Consequently, the formation of spaces in the transitional lithology and interlayer within rocks and fluid overpressure is promoted (3) The similar trace element and REE patterns in the micritic laminae, host rocks, and veins indicate that the ions leached from the adjacent host rocks and micritic laminae are concentrated in the nascent veins. This concentration is attributed to focused fluid expulsion along these pathways after the OM matures and produces organic and carbonic acids. Further, the Sr and Ba for the precipitation of strontianite and barite, respectively, originated from the host rocks rather than from the external fluids (4) The host rock carbonate is enriched in 87 Sr and depleted in 13 C. This is caused by isotopic exchanges with the noncarbonate host rock components. The micritic carbonate lamina grows as the median zone along with the fibrous calcite, probably reflecting the nature of the vein-forming fluids. The carbon and strontium isotopic compositions of micritic laminae and veins are also very similar. The negative shift of the oxygen isotopes of the veins is induced by the dehydration of clay minerals and the increase in temperature. This indicates that the materials of the veins are derived from the dissolution of carbonate within the shale in a closed system (5) The migration of the diagenetic fluid is vertical and limited in the interlayer. Such migration promotes the growth of calcite fibers by crystallization. Blocky crystals rapidly grow from both sides of the fracture toward the center, allowing the hydrocarbons from the host rocks to occupy the intercrystalline pores. In addition, the shale with bedding-parallel calcite vein is the favorable lithofacies for shale oil reservoirs

Data Availability
The clay minerals data used to support the findings of this study have been deposited in Figure 4 repository (doi:10 .1016/j.marpetgeo.2016.02.022; doi:10.1016/j.marpetgeo .2016.12.004). The C-O isotopes of coarse-grained vein data used to support the findings of this study have been deposited in Figure 12 repository (doi:10.1007/s12583-016-0946-3). The TOC, δ13Ccarb, δ18O, δ13Corg, 87Sr/86Sr, trace elements and REEs data used to support the findings of this study are included within the article. The XRD analysis, TOC content, and Rock-Eval parameters of rocks data used to support the findings of this study were supplied by the Shengli Oilfield Company, SINOPEC, under license and so cannot be made freely available. Requests for access to these data should be made to Zhenhuan Shen (shenzh9311@163.com).

Conflicts of Interest
The authors declare that they have no conflicts of interest.   Figure 15: Pattern of hydrocarbon accumulation of blocky and fibrous veins (hy: hydrocarbon). 20 Geofluids