Geological and geochemical characterization of lacustrine shale: A case study of the Jurassic Da'anzhai member shale in the central Sichuan Basin, southwest China

Geological and geochemical characterization of lacustrine shale: A case study of the Jurassic Da'anzhai member shale in the central Sichuan Basin, southwest China

Accepted Manuscript Geological and geochemical characterization of lacustrine shale: A case study of the Jurassic Da'anzhai member shale in the centra...

6MB Sizes 1 Downloads 76 Views

Accepted Manuscript Geological and geochemical characterization of lacustrine shale: A case study of the Jurassic Da'anzhai member shale in the central Sichuan Basin, southwest China Qilu Xu, Bo Liu, Yongsheng Ma, Xinmin Song, Yongjun Wang, Zhangxin Chen PII:

S1875-5100(17)30380-3

DOI:

10.1016/j.jngse.2017.09.008

Reference:

JNGSE 2314

To appear in:

Journal of Natural Gas Science and Engineering

Received Date: 21 June 2017 Revised Date:

15 September 2017

Accepted Date: 27 September 2017

Please cite this article as: Xu, Q., Liu, B., Ma, Y., Song, X., Wang, Y., Chen, Z., Geological and geochemical characterization of lacustrine shale: A case study of the Jurassic Da'anzhai member shale in the central Sichuan Basin, southwest China, Journal of Natural Gas Science & Engineering (2017), doi: 10.1016/j.jngse.2017.09.008. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT

2 3

4

First: Manuscript with track changes

Geological and geochemical characterization of lacustrine shale: a case study of the Jurassic Da'anzhai Member shale in the central Sichuan Basin, Southwest China

RI PT

1

Qilu Xua,b,c, Bo Liuc,*, Yongsheng Maa, Xinmin Songd, Yongjun Wangd, Zhangxin Chenb a

School of Earth Sciences and Resources, China University of Geosciences(Beijing), Beijing 100083, China

6

b

Chemical and Petroleum Engineering, University of Calgary, Calgary, Alberta T2N1N4, Canada

7

c

Oil and Gas Research Center, Peking University, Beijing 100871, China

8

d

PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China

10

Abstract: :

M AN U

9

SC

5

Although numerous progresses have been achieved in characterizing marine shales, studies that are associated with lacustrine shales are limited.

To study lacustrine shales, samples were collected from

12

the Jurassic Da'anzhai Member in the central Sichuan Basin of China, and their mineralogical, reservoir,

13

OM, and paleoenvironmental characteristics were determined, as well as the relationships between them.

14

Analysis of trace elements reveals that the shales formed in paleoenvironments that were oxic to suboxic,

15

dry to humid, had moderate to strong weathering, and were characterized by fresh to salt water

16

conditions. These environments are more variable than those of marine shales. The paleoenvironmental

17

conditions and mineralogy of the shales, particularly the oxic to suboxic paleo-redox conditions, resulted

18

in the relatively low levels (0.11-2.18%, average 0.97%). However, based on evaluation criteria for

19

continental source rocks, the OM is of high quality because of its high level of maturity (Ro: 0.95-1.43;

20

Tmax: 428-500°C) and favorable kerogen type (II2). There are well-developed intraparticle pores,

AC C

EP

TE D

11

1

ACCEPTED MANUSCRIPT

interparticle pores and microcracks. The SBET (5.42-10.69 m2/g, average 6.92 m2/g) and VBJH (1.19-4.10

22

mL/100 g, average 2.97 mL/100 g) values also indicate good nanoscale storage space. Terrestrial

23

minerals (i.e., quartz and clay) and authigenic carbonate minerals (i.e., calcite) are, respectively,

24

positively and negatively correlated with the nanoscale storage space, .. Small pores (3-5 nm) dominate

25

the nanoscale storage space. The isotherms and hysteresis loops are of Type П and Type H3, respectively,

26

which indicates wedge-shaped pores. However, the hysteresis loops indicate that the lacustrine shale has

27

more dead-end pores and larger pores with more complex microstructures than other lacustrine

28

reservoirs. In general, the Da’anzhai lacustrine shale has the potential for unconventional oil and gas

29

exploration. The lacustrine shale’s mineralogy, reservoirs, and OM are closely related to each other, and

30

their differences are mainly caused by the paleoenvironmental conditions in which the shale formed.

32

SC

M AN U

Keywords: Lacustrine shale; Paleoenvironment; Nanoscale reservoir; Organic matter; Minerals; Mechanisms.

TE D

31

RI PT

21

1. Introduction

34

The high commercial production of shale oil and gas in North America has made shale a focus of

EP

33

exploration in many countries and regions. However, the most widely developed shales were formed in

36

marine systems; few studies have focused on the characteristics of lacustrine shales. Additionally,

37

lacustrine shales are widely distributed in many areas, such as Africa, South America, Southeast Asia

38

and China (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et al., 2015;

39

Jiang et al., 2016a; Yang et al., 2015).

AC C

35

2

ACCEPTED MANUSCRIPT

40

The most obvious characteristic of lacustrine shales is that they are subject to more variable paleoenvironments than marine shales. The formation environments of lacustrine and marine shales are

42

very different in terms of the water depth, energy, provenance, sedimentary type and biological effects

43

(Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002). These features lead to

44

differences in their basic characteristics including mineralogy, organic matter (OM), and reservoir type

45

(Table 1). Minerals are the material basis of shales, and they affect not only OM generation but also

46

reservoir development (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et

47

al., 2015; Yang et al., 2015; Jiang et al., 2016a). A major breakthrough in marine shale oil and gas

48

exploration was based on the recognition that shale acts not only as a source rock but also as a reservoir.

49

Unlike studies of conventional reservoirs, research on unconventional reservoirs has primarily focused

50

on nanoscale reservoirs. Previous studies of marine and lacustrine shales suggest that shale gas is mainly

51

controlled by a nanoscale pore system that is closely related to the mineral characteristics (Zou et al.,

52

2011; Morad et al., 2010; He et al., 2012; Luo et al., 2013; Dang et al., 2015; Chalmers et al., 2012;

53

Curtis et al., 2012; Jiang et al., 2015; Li et al., 2016; Jiang et al., 2017; Zhang et al., 2017). OM is the

54

basic material that generates hydrocarbons, and its type and maturity dominate the formation and

55

enrichment of shale oil and gas. The preservation and enrichment of OM in sediments is controlled by

56

many factors, including the primary productivity, paleo-redox conditions, nutrient availability, clastic

57

influx and minerals (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et al.,

58

2015; Yang et al., 2015; Jiang et al., 2016a; Xu et al., 2017). In summary, the paleoenvironment,

59

minerals, reservoir and OM, which are the basic characteristics of shales, have been the focuses of shale

AC C

EP

TE D

M AN U

SC

RI PT

41

3

ACCEPTED MANUSCRIPT

oil and gas studies, and these characteristics are often interrelated. Additionally, studies of the

61

characteristics and mechanisms have mainly focused on marine shales; fewer studies have examined

62

lacustrine shales.The lacustrine shale of the Jurassic Da’anzhai Member from the central Sichuan Basin

63

of China has been interpreted as fractured reservoirs, low permeability fractured reservoirs and tight

64

reservoirs. However, most studies have only focused on the conventional carbonate reservoirs, and the

65

shales have only been considered as source rocks for the carbonate reservoirs (Lu et al., 2014; Chen et

66

al., 2015). Encouragingly, the Da’anzhai shale of the Yuanba Block in the northwestern Sichuan Basin,

67

which has similar geological conditions to those of the central Sichuan Basin, hosts high-yield wells (He

68

et al. 2012; Zhu, et al. 2016). In addition, some lacustrine shales have yielded oil or gas, such as the

69

Triassic Chang 7 and Chang 9 Member in the Ordos Basin, the Paleogene strata (Es3 and Es4) in the

70

Bohai Bay Basin and the Cretaceous strata in the Songliao Basin (Qing-1) (Zou et al., 2011; Luo et al.,

71

2013; Dang et al., 2015; Li et al., 2015; Yang et al., 2015; Jiang et al., 2016b). These successful

72

developments strongly indicate that the Da’anzhai shale in the central Sichuan Basin is worth studying.

SC

M AN U

TE D

In this study, we assessed the Da’anzhai lacustrine shale’s characteristics. We performed

EP

73

RI PT

60

paleoenvironmental evaluation using TE, mineralogical evaluation usingXRD analysis, reservoir

75

evaluation using FE-SEMand LTNA analyses, and OM evaluation using analyses of the TOC,

76

chloroform bitumen ‘A’, Rock-Eval, carbon isotopes, kerogen macerals and organic elements. We also

77

compared the Da’anzhai shale to other lacustrine shales which have been successfully developed and

78

assessed the similarities and differences of lacustrine shales. In addition, we summarized the

AC C

74

4

ACCEPTED MANUSCRIPT

relationships between these characteristics and mechanisms affecting them. Finally, we judged the

80

development potential of the Da’anzhai shale.

AC C

EP

TE D

M AN U

SC

RI PT

79

5

ACCEPTED MANUSCRIPT

81

Table 1 Comparisons between basic parameters in the Da’anzhai shale and other successfully developed shales

RI PT

82 83 84

lacustrine shale

Central Sichuan Basin, China

layer

Barnett

Ordos Basin,

Ordos Basin,

Bohai Bay

Bohai Bay

Songliao

Fort Worth Basin,

China

China

Basin, China

Basin, China

Basin, China

Texas, US

downfaulted

downfaulted

downfaulted

basin

basin

basin

Paleogene

Paleogene

Cretaceous

Carboniferous

shale with

shale with

shale with

siliceous and

siltstone

siltstone

argillaceous

carbonaceous

limestone

shale, limestone,

interlayers

dolomite

Ohio Appalachian Basin, Kentucky, US

polycyclic superimposed

basin

basin

Triassic

Triassic

shale with

shale with

siltstone

siltstone

interlayers

interlayers

interlayers

interlayers

1900~3200

1500~2500

1600~2600

1500~4200

1500~5200

1000~2500

1980~2590

510~1800

10~70

30~70

10~15

100~500

250~350

70~150

15~60

9~31

(20~80) /54

(37~53) /43

(29~56) /45

(25~80) /36

(30~80) /36

(30~80) /36

35~80

45~60

Jurassic

shells and

combination

carbonaceous interlayers

effective

(%) brittle minerals

Qing-1

polycyclic

lithological

thickness

Es4

superimposed

foreland basin

shale with

depth (m)

Es3

TE D

basin type

Chang 9

EP

basin

marine shale

Chang 7

M AN U

Da'anzhai

AC C

parameters

SC

85

(%)

6

foreland basin

piedmont depression Devonian interbed of black shale and gray siliceous shale

ACCEPTED MANUSCRIPT

0.20~5.92

1~12

10

2~8

1.3~9.3

6~12

4~5

4~7

permeability

0.002~0.9

0.01~0.3

(0.3~5.0) /2.97

<0.50

<0.50

<0.15

0.01

<0.1

TOC (%)

(0.1~2.2) /1.0

(0.3~36) /8.3

(0.3~11.3) /3.1

(0.5~13.8) /3.5

(0.8~16.7) /3.2

(0.4~4.5) /2.2

2.0~7.0

0.5~4.7

Ro (%)

0.9~1.5

0.8~1.2

0.9~1.3

0.4~1.2

0.4~2.0

0.5~1.5

1.1~2.2

0.3~1.3

pyrolysis gas

pyrolysis gas

pyrolysis gas

Ⅰ-Ⅱ





8.50~9.91

1.69~2.83

genetic type

biogenic and

pyrolysis gas

pyrolysis gas

pyrolysis gas



Ⅰ-Ⅱ1

Ⅰ-Ⅱ1

Ⅰ-Ⅱ1

0.87~1.98

0.92~2.91

(1.68~4.25)

(0.60~3.70)

/2.92

/2.00

pyrolysis gas

pyrolysis gas

kerogen

2011; Luo et references

al., 2013; Lu et al., 2014; Chen et al., 2015;

2011; Luo et al., 2013; Yang et al., 2015; Jiang et al., 2016a.

SC

Zou et al.,

Zou et al.,

2011; Luo et

Zou et al., 2011;

Zou et al.,

al., 2013; Shi et

Luo et al., 2013;

2011; Luo et

al., 2013; Yang

Dang et al.,

al., 2013; Li,

et al., 2015;

2015; Jiang et

2014; Jiang et

Jiang et al.,

al., 2016a

al., 2016a

2016a.

M AN U

Zou et al.,

0.31~1.24

TE D

(m3/t)

EP

gas content

Ⅰ-Ⅱ1

AC C

type

RI PT

porosity (%)

7

Schmoker, 2002;

Schmoker, 2002;

Zou et al.,

Bowker et al.,

Bowker et al.,

2011; Luo et

2007; Jarvie et al.,

2007; Jarvie et al.,

al., 2013; Jiang

2007; Martini et

2007; Martini et

et al., 2016a

al., 2008; Lu et al.,

al., 2008; Lu et al.,

2014

2014.

ACCEPTED MANUSCRIPT

86

2. Geological setting

89 90 91

Fig. 1. Study location, geological map (Chen et al., 2015) and simplified stratigraphic column of the Jurassic Da’anzhai

AC C

88

EP

TE D

M AN U

SC

RI PT

87

Member in the central Sichuan Basin.

The study area is located in the central Sichuan Basin, which includes the Nanchong and Suining

92

areas. It is bounded by the Huafu Mountains to the east and the Longquanshan basement faults to the

93

west, and it covers an area of 6×104 km2 (Fig. 1) (Lu et al., 2014; Chen et al., 2015). A lacustrine

94

environment dominated the Sichuan Basin during the Jurassic Period, and the Da’anzhai Member 8

ACCEPTED MANUSCRIPT

represents the largest and deepest stage of this environment. The sedimentary facies of the Da’anzhai

96

Member can be divided into shore-shallow lacustrine, shallow lacustrine and semi-deep lacustrine

97

subfacies, which include shell beach, calcareous flat, mud flat, and semi-deep lake microfacies. These

98

facies are distributed in concentric belts around the lake basin’s center in the Yilong-Pingchang areas

99

(Fig. 1) (Lu et al., 2014; Chen et al., 2015). From top to bottom, the Da’anzhai Member in the central

100

Sichuan Basin can be divided into main three sub-members: Da1, Da13 and Da3 (Fig. 1). The thickest

101

Da13 sub-member (30-50 m) is a set of widely distributed dark mudstones with lamellated fractures and

102

rich OM; this sub-member is considered the most important source rock for hydrocarbons. The Da3

103

sub-member deposits (5-14 m), which are initial deposition products, are primarily (mud) shell

104

limestones with dark gray mudstone interlayers. The Da1 sub-member (30-45 m thick) primarily

105

consists of shales and (mud) shell limestones; this sub-member contains the sedimentary products of the

106

period of basin uplift (Lu et al., 2014; Chen et al., 2015).

TE D

M AN U

SC

RI PT

95

3. Samples and experimental methods

108

3.1. Samples

109

Thirty-nine representative core samples were obtained from the Da’anzhai Member shale in the

AC C

EP

107

110

central Sichuan Basin based on drilling, core observations, stratigraphic data, and log responses (Fig. 1

111

and Table 3).

112

3.2. Methods

113

All 39 shale samples were analyzed to determine their bulk mineralogy viaXRD

114

using a X'Pert

Pro MPD device with a working voltage of 40 kV and a current of 40 mA in the Key Laboratory of

9

ACCEPTED MANUSCRIPT

Orogenic Belts and Crustal Evolution at Peking University. To investigate the clay mineral fraction in

116

detail, textured mounts were measured four times at a goniometer rate of 0.5θ/minute with a registration

117

range from 2 to 42°2θ at the Beijing Research Institute of Uranium Geology.

118

RI PT

115

The micromorphology and oiliness of the samples were determined at the Key Laboratory of Marine Reservoir Evolution and Hydrocarbon Accumulation of the China University of Geosciences

120

(Beijing) using an Axio Scope A1 pol Fluorescence-Polarization Microscope. An INCA Synergy system

121

at the Key Laboratory of Orogenic Belts and Crustal Evolution at Peking University was used to analyze

122

the microstructure of each sample via FE-SEM and the mineral composition of each sample via EDS.

123

The samples were subjected to argon-ion polishing to obtain higher resolution images.

M AN U

124

SC

119

LTNA analysis is capable of characterizing nanoscale pores in porous media, and adsorption tests were performed with dried samples (60-80 mesh) on an ASAP2020M device at the College of

126

Chemistry and Chemical Engineering of the China University of Petroleum (Beijing). N2 adsorption

127

isotherms were obtained at -196.15°C within a relative pressure range of 0.005 to 0.998. The BET

128

method for the equivalent specific surface areas and theBJH method for the total pore volumes were

129

applied in the calculations (Brunauer et al., 1938; Barrett et al., 1951).

AC C

EP

TE D

125

130

The OM richness was characterized by measuring the TOC values of 25 samples using a GHM-02

131

analyzer; 10% hydrochloric acid was used to remove carbonates, and the test sensitivity was 10-13 mg/g.

132

Chloroform bitumen ‘A’ was extracted from the crushed shale powders (<0.09 mm), and the weight was

133

calculated at the China University of Petroleum (Beijing). A Rock-Eval 6 instrument was used to

134

perform a pyrolysis analysis, which yielded the temperature of the Tmax, S0, S1 and S2. The OM maturity

10

ACCEPTED MANUSCRIPT

was determined based on the Ro value, which was obtained using a reflected light microscope under oil

136

immersion and calculated as the mean of 10-30 measurements taken for a single sample. Additionally,

137

measurements were taken to determine the isolated kerogen content and identify the organic kerogen

138

macerals, carbon isotope content and organic element content according to the (GB/T) 19144-2010,

139

(SY/T) 5125-2014, (GB/T) 18340.2-2010, (GB/T) 19143-2003 standards, respectively. The kerogen

140

measurements and the rock pyrolysis analysis were performed at the Keyuan Engineering Technology

141

Testing Center.

SC

M AN U

142

RI PT

135

The TE concentrations of 17 samples were determined using ICP-MS, which was performed at the Beijing Research Institute of Uranium Geology. The detailed sample processing procedure as well as the

144

analytical precision and accuracy followed the methods described for the (GB/T) 14506.28-2010

145

standard.

TE D

143

4. Results and discussions

147

4.1. Paleoenvironmental properties and differences

148

Lacustrine shales are different from marine shales mainly due to their formation paleoenvironments.

EP

146

Lakes are complex and dynamic systems, and the development of lacustrine shales is controlled more

150

significantly by the sedimentary environment (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002;

AC C

149

11

ACCEPTED MANUSCRIPT

Table 2 Trace elements of the Da’anzhai lacustrine shale. G10-2

G4-3

G6-2

J45-1

J45-2

J53-1

J61-1

J61-2

J61-3

J61-4

Li-1

Li-2

P1-1

S2-1

S2-3

S2-4

X44-1

Th(ppm)

11.70

10.40

5.43

4.32

12.00

10.60

1.72

4.13

7.74

10.20

11.30

9.50

11.60

2.31

10.50

6.21

13.70

Ni/Co

2.68

2.98

3.26

4.10

2.87

2.72

5.49

3.87

3.50

2.51

2.65

2.97

3.07

3.22

2.32

3.37

2.78

RI PT

code

1.76

1.58

1.44

1.50

1.64

1.64

2.47

1.56

1.33

1.69

1.60

1.50

1.61

1.77

1.30

1.41

1.45

0.78

0.80

0.63

0.63

0.78

0.78

0.67

0.68

0.70

0.77

0.78

0.74

0.78

0.61

0.68

0.68

0.77

V/Sc

9.78

8.70

6.88

7.73

9.36

8.91

13.81

8.41

7.94

9.03

7.93

7.21

9.22

9.80

8.54

7.14

7.81

Th/U

2.77

4.24

4.45

2.79

4.07

2.95

1.28

3.53

4.63

3.81

4.25

4.63

4.20

2.74

5.25

4.60

4.72

Sr/Cu

3.75

4.60

42.98

77.10

3.50

3.62

58.70

35.54

53.09

8.46

7.84

19.01

3.35

7.47

24.69

6.05

Sr/Ba

0.29

0.26

3.66

3.08

0.23

0.24

4.77

2.78

3.45

0.52

0.46

1.28

0.18

3.77

0.44

1.84

0.36

118.7

115.3

109.8

131.7

110.3

187.5

141.2

128.9

194.3

112.9

60.7

130.6

112.6

153.9

3

6

4

1

7

7

6

9

0

9

6

7

2

3

7.79

9.03

6.39

6.65

9.34

9.58

8.46

8.67

8.49

8.04

9.44

4.19

9.93

9.86

8.88

7.88

10.03

(La/Yb)N

6.65

7.75

6.57

8.37

6.82

8.04

9.85

9.82

11.14

8.00

7.50

4.35

7.49

8.62

7.67

9.05

(La/Sm)N

4.66

4.74

3.04

3.32

5.58

5.42

3.19

3.06

2.37

3.88

5.71

2.59

5.90

3.07

4.20

3.90

5.01

ΣLREE/ ΣHREE

88.73

38.15

M AN U

m)

130.97

TE D

ΣREE(pp

SC

V/Cr V/(V+Ni)

80.5 4

10.8 8

0.99

1.04

1.52

1.97

0.77

0.97

2.04

2.01

2.92

1.45

0.84

1.41

0.79

2.10

1.35

1.38

1.17

(Ce/Yb)N

4.48

5.38

4.98

5.83

4.54

5.42

7.21

7.33

9.24

5.58

4.97

3.22

5.00

8.49

6.10

5.59

6.15

δCe

0.92

0.94

0.97

0.95

0.92

0.93

0.94

0.93

0.99

0.92

0.92

0.92

0.93

0.96

0.95

0.97

0.92

CSr(‰)

0.19

0.18

0.77

1.12

0.17

0.16

0.63

0.56

0.88

0.36

0.27

0.61

0.15

0.72

0.12

0.35

0.23

AC C

EP

(Gd/Yb)N

XN, where N refers to the chondrite-normalized value (Gromet et al., 1984); δCe=CeN/(LaN×PrN)1/2 (Wright et al.,1984; Murray et al.,1990).

12

ACCEPTED MANUSCRIPT

Schmoker, 2002). TEs can be used as effective quantitative indicators of the paleoenvironment by

2

comparing their type and abundance with standard values (Wright et al., 1984; Murray et al., 1990;

3

Manning, 1994; Ross et al., 1995; Kimura and Watanabe, 2001; Jarvis et al., 2001; Rimmer, 2004; Jones

4

and Roy et al., 2007).

5 6

TE D

M AN U

SC

RI PT

1

Fig. 2. Identification of paleoenvironments using TEs (a-c) and clay (d).

4.1.1. Paleo-redox conditions

8

Due to their responses to changing redox conditions, the ratios of redox-sensitive TEs, such as

AC C

9

EP

7

Ni/Co, V/Cr, V/Sc, Th/U, and δCe, are usually used to evaluate redox conditions. According to Jones

10

and Manning (1994), V/Cr ratios <2 indicate oxic conditions, V/Cr ratios from 2 to 4.25 indicate dysoxic

11

conditions, and V/Cr ratios >4.25 indicate suboxic to anoxic conditions. Kimura and Watanabe (2001)

12

found that the V/Sc ratio is positively related to oxic conditions. Jones and Manning (1994) showed that

13

Ni and Co are enriched in pyrite, and high Ni/Co ratios indicate anoxic conditions, Ni/Co ratios <5 13

ACCEPTED MANUSCRIPT

indicate oxic conditions, Ni/Co ratios between 5 and 7 indicate dysoxic conditions, and Ni/Co ratios >7

15

indicate suboxic to anoxic conditions. Similarly, Th/U ratios <1.5 represent anoxic conditions, Th/U

16

ratios between 1.5 and 3 represent dysoxic to suboxic conditions, and Th/U ratios >3 represent oxic

17

conditions (Jones and Manning, 1994). In addition, δCe values >1 (positive anomaly) and <0.95

18

(negative anomaly) suggest anoxic and oxic to suboxic conditions, respectively (Wright et al., 1984;

19

Murray et al., 1990).

SC

For the Da’anzhai shale in the central Sichuan Basin (Table 2), the ranges of the Ni/Co, V/Cr, V/Sc,

M AN U

20

RI PT

14

and Th/U ratios are relatively wide, from 2.32-5.49, 1.30-2.47, 6.88-13.81, and 1.28-5.25, respectively.

22

The values of δCe vary from 0.92 to 0.99 (average 0.94), which indicate a weakly negative anomaly

23

(Wright et al., 1984; Murray et al., 1990).Collectively, these proxies reveal that the shale was deposited

24

under variable paleo-redox conditions that were mainly oxic to suboxic (Fig. 2).

TE D

21

4.1.2. Paleo-salinity

26

Sr/Ba values vary directly with the distance to a sea or lake coast; Sr/Ba ratios <1 suggest

EP

25

freshwater conditions, and Sr/Ba ratios >1 suggest salt water conditions (Jarvis et al., 2001; Kimura and

28

Watanabe, 2001). In addition, low Sr concentrations (0.1-0.3‰) indicate a freshwater setting, and high

29

values (0.8-1‰) represent a salt water setting (Jarvis et al., 2001; Kimura and Watanabe, 2001). Th

30

easily dissolves in acidic water and hydrolyzes to oxide or hydroxide sediments under alkaline

31

conditions; therefore, high Th values indicate a salt water environment (Jones and Manning, 1994).

32 33

AC C

27

In the Da’anzhai shale samples (Table 2), the Sr/Ba ratios vary widely (0.18-4.77) with an average of 1.62. The Th content also varies widely (1.72-13.70 ppm) and averages 8.43 ppm, and the Sr 14

ACCEPTED MANUSCRIPT

34

concentrations vary from 0.12 to 1.12‰ (average 0.44‰). Generally, these proxies indicate variable

35

paleo-salinity conditions including both fresh and salt water settings (Fig. 2).

RI PT

4.1.3. Paleo-weathering

M AN U

SC

36

37

39

Fig. 3. Chondrite normalized REE distribution curves for REE minerals. Chondrite values are from Gromet (1984).

When rocks are subjected to strong weathering, heavy rare earth elements (HREEs) dissolve and

TE D

38

migrate much more easily than light rare earth elements (LREEs). Therefore, high ΣLREE/ΣHREE

41

values indicate strong paleo-weathering conditions (Kimoto et al., 2006; Stevens and Quinton, 2008). In

42

addition, the (La/Yb)N, (La/Sm)N, (Gd/Yb)N and (Ce/Yb)N ratios, which are positively related to LREE

43

enrichment, are also positively correlated with paleo-weathering (Table 2) (Wright et al., 1984; Murray

44

et al., 1990; Ross et al., 1995; Roy et al., 2007).

AC C

45

EP

40

The ΣLREE/ΣHREE rates are between 4.19 and 10.03 (average 8.39), which indicate LREE

46

enrichment. The (La/Yb)N ratio is represented by the slope of the chondrite-normalized REE distribution

47

curve (Fig. 3), which indicates the degree of the graph’s inclination (Ross et al., 1995; Roy et al., 2007).

48

The (La/Yb)N values vary from 4.35 to 11.14 (average 8.15), which indicate that the curves are tilted to 15

ACCEPTED MANUSCRIPT

the right and that the samples are rich in LREE-acidic rocks (Fig. 3). The (La/Sm)N ratio, which

50

represents the degree of LREE fractionation, is positively related to the ΣLREE value (Ross et al., 1995;

51

Roy et al., 2007). The (La/Sm)N ratios range from 2.37 to 5.90 (average 4.10), which indicates a P(E)

52

model with abundant LREEs. The (Gd/Yb)N ratio, which indicates the degree of HREE fractionation, is

53

negatively correlated with the ΣHREE value (Ross et al., 1995; Roy et al., 2007). The (Gd/Yb)N values

54

vary from 0.77 to 2.92 (average 1.46), and the (Ce/Yb)N ratios vary from 3.22 to 9.24 (average of 5.85),

55

which indicate that the Da’anzhai shale is depleted in HREEs (Fig. 3). In summary, the REE results

56

(Table 2) indicate that the shale was deposited under moderate weathering conditions with periods of

57

strong weathering.

M AN U

SC

RI PT

49

4.1.4. Paleoclimate

59

∑REE values are closely related to paleoclimate conditions; that is, ∑REE values are higher in

TE D

58

warm and humid climates and lower in cold and dry climates (Wright et al., 1984; Murray et al., 1990;

61

Ross et al., 1995; Roy et al., 2007). The ΣREE values of the Da’anzhai shale vary widely from 38.15 to

62

194.3 µg/g (average 121.59 µg/g) (Table 2). The Sr/Cu ratio is another useful indicator of lacustrine

63

paleoclimate, with high values reflecting humid climates (Reheis, 1990; Chen et al., 2009).The Sr/Cu

64

ratios in the Da’anzhai shale also vary widely from 3.35 to 80.54%. In addition, illite usually forms in

65

dry conditions by the weathering depotassication of feldspar, mica and other aluminosilicate minerals.

66

Conversely, kaolinite forms in humid climates by the eluviation of feldspar, mica and pyroxene.

67

Therefore, high illite/clay ratios indicate dry, saline waters with high K+ content, and high kaolinite/clay

68

ratios suggest humid climates (Vanderaveroet, 2000; Gingele et al., 2001). The illite/clay ratios of the

AC C

EP

60

16

ACCEPTED MANUSCRIPT

Da’anzhai shale are relatively high (between 26 and 86%, average 52%). However, the kaolinite/clay

70

ratios are relatively low (between 3 and 19%, average 11%) (Fig. 2d). These results indicate a variable

71

climate with both dry and humid periods as well as a drying trend from humid climates.

72

RI PT

69

In general, the TE contents and types t in lacustrine shales are more complex and varied than those in marine shales, which indicates that lacustrine shales tend to experience more variable

74

paleoenvironmental conditions (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002).

75

The lacustrine Da’anzhai shale was formed under relatively dry climatic conditions that included both

76

arid and humid stages. In addition, the shale experienced a moderate degree of paleo-weathering with

77

stages of strong weathering, variable oxic to suboxic paleo-redox conditions that were weakly reducing,

78

and a predominantly freshwater setting interspersed with periods of saltwater.

M AN U

SC

73

4.2. Mineralogy properties and differences

80

As shown in Table 3, the Da’anzhai shale primarily contains clay (42.4%), quartz (27.7%) and

81

calcite (21.1%). The minerals in these lacustrine shales vary widely. Clay comprises the majority of the

82

shale and ranges from 16.3 to 68.4%, while the quartz content ranges from 5.6 to 53.4%. Carbonate

83

minerals, which mainly include calcite, dolomite and aragonite, range from 3.0 to 46.6% with an

84

average of 27.4%. The Da’anzhai shale contains more brittle minerals than the other lacustrine shales

85

and marine shales (Table 1 and Fig. 4). For example, the brittle mineral content of the Da’anzhai shale

86

ranges from 24.2 to 80.3%, with an average of 53.7%, indicating medium-high brittleness that is

87

favorable for natural fractures and hydraulic fracturing. Similar to other lacustrine shales, lacustrine

88

shales have more siliceous or calcareous interlayers than marine shales, which can cause the wider

AC C

EP

TE D

79

17

ACCEPTED MANUSCRIPT

ranges of carbonate and quartz minerals and improve the brittleness index and fracturing (Figs1 and

90

5I).Ternary diagrams were constructed to compare the mineral compositions of hot shales that have been

91

successfully exploited including lacustrine shales (Fig. 4a) and marine shales (Fig. 4b). Generally, the

92

comparison shows that the mineral composition of the lacustrine shales is more variable and complex

93

than in marine shales (Fig. 4). For example, the Da’anzhai shale has relatively wide ranges of carbonate

94

(1-70%), clay (15-70%) and quartz (5-50%) (Fig. 4a). In addition, the Da’anzhai shale has the wider

95

mineral ranges than other lacustrine shales, especially for the range of carbonate minerals.

M AN U

SC

RI PT

89

EP

97

TE D

96

Fig. 4. Ternary diagrams of mineral compositions. (a) Lacustrine shales: 1=Da’anzhai shale, 2=Jurassic shale in Sichuan

99

Basin (Jiang et al., 2017), 3=Triassic shale in Sichuan Basin(Jiang et al., 2017), 4=ES3 shale (Wang et al., 2015a). (b) Marine

AC C

98

100

shales: 5 Woodford shale (Jarvie et al. 2007 and Han et al. 2013), 6=Barnett shale (Jarvie et al. 2007 and Han et al. 2013),

101

7= Wufeng-Longmaxi shale (Wu et al., 2014), 8=Ohio shale (Jarvie et al. 2007 and Han et al. 2013).

102

Table 3 Minerals and TOC results for the Da’anzhai lacustrine shale.

103

18

ACCEPTED MANUSCRIPT

G10-3 2700.7 G10-4 2646.6 G4-1 2385.3 G4-2 2395.8

subquartz feldspar calcite dolomit aragonite siderite pyrite hematite total clay others carbonate brittleness member (%) (%) (%) e (%) (%) (%) (%) (%) (%) (%) (%) minerals (%) 44.2 3.0 52.8 3.0 47.2 Da1 22.1 1.1 0.9 1.2 6.3 68.4 2.1 24.2 Da13 37.5 2.0 2.9 2.6 55.0 2.9 40.4 Da13 1

4.2 2.1

G4-3 2405.7

Da Da1 Da1 Da13

28.1 28.2 35.2 25.8

2.5

G6-1 2537.1

Da13

21.4

2.6

G6-2 2580.1

8.7

1.2

47.6

19.5

1.1

0.8

J45-1 2632.4 J45-2 2656.0

Da3 Da1 Da13

12.5 23.4

0.6 6.1

27.9 3.9

36.3

0.4 0.4

0.7 1.1

J53-1 2827.3

Da13

31.9

4.1

1

14.4 55.3 14.3 17.5 69.4

6.4

17.0 13.4 12.5

4.5

4.6

0.7

1.6

1.2

Da13

8.5

J61-6 Li-1 Li-2 Li-3

2665.9 1713.8 1726.1 1734.3

Da1 Da1 Da1 Da13

25.1 30.4 24.2 40.7

1.8 2.6 1.0 2.3

15.6 9.8 28.9 3.1

1.0

P1-1 3215.9

Da13

23.3

0.9

0.6

1.0

P1-2 3221.2

Da13

34.8

4.1

P1-3 3239.2

Da13

30.3

14.9

P1-4 3258.6 S2-1 2836.1 S2-2 2849.9

3

Da Da1 Da1 Da13

23.2 19.8 40.3

51.8 48.2 10.7

Da3 Da3

18.9

S2-4

2889.0

S2-5 2891.6 X20-1 1844.4 X28-1 1952.2 X44-1 2089.4 X44-2 2093.0

Da1 Da13 Da13 Da13 Da13 Da13 1

53.4 49.4

3.6

LQ-1 3548.1 X8-1 1780.7

W8-1 2065.7

Da Da1

1.3

0.2 0.8 0.0 4.2

3.2

6.3 0.0

5.3 2.2

69.2 41.6 54.1

0.40 0.98

5.3

31.1

1.43

4.2

7.4

28.8

1.18

20.1

1.0

67.1

75.8

20.5 64.0

1.1 1.1

64.2 3.9

76.7 27.3

55.5

5.4

22.6 26.7 27.7 33.1

10.3

1.1

0.81

41.1 13.4 18.9

1.0 1.1 0.9 2.3

20.0 46.2 53.3 38.5 53.9

0.7

31.9

2.18

59.3 64.8 57.0 24.3

76.2 70.4 71.4 59.1

1.49

69.4

77.9

0.30

26.9 9.8 35.2 3.1

52.0 40.2 59.4 43.8

1.17 1.07 0.69

6.4

67.8

1.6

24.9

1.59

0.9

60.2

4.1

38.9

1.51

3.6

51.2

14.9

45.2

1.55

1.5 0.2 0.9

16.3 29.6 42.5

57.1 50.4 10.7

80.3 70.2 51.0

0.61 0.11 1.18

1.9 2.0

47.7

1.3

33.4

49.0

66.6

0.33

48.9

1.5

30.7

50.4

69.3

0.18 0.76

8.1

24.9 35.9

1.4

27.9

3.0

4.0

23.8

3.0

43.2

AC C

X44-3 2116.7

Da13

17.6

1.0

2.1

TE D

2867.1

44.9 9.5 42.7 6.8

0.9

TOC (%) 0.82 1.33

60.9

SC

J61-5 2706.8

S2-3

3.5

M AN U

Da1 Da13

16.9 5.6 14.4 34.8

26.6 54.9 38.4 64.8

3.1

2662.6 2664.4 2667.6 2698.1

Da Da1

3.5

7.4

25.6 8.5

EP

J61-1 J61-2 J61-3 J61-4

22.9

RI PT

depth (m) G10-1 2650.8 G10-2 2668.4 code

46.6 38.7

2.3

3.1 3.8 3.7 9.5

1.5

8.1

53.4 57.5

44.3 51.8

0.7

28.7 8.5

53.6 44.4

0.53

60.8

0.6

4.0

31.9

1.15

52.7

76.5

20.5

44.1

2.9

18.7

40.5

2.6

2.5

34.3

28.2

2.2

18.5

7.8

43.3

23.9

1.9

35.2

11.4

26.6

2.5

50.8

1.1

104

W8-2 2069.5

105

4.3. Characterization and controlling factors of lacustrine shale reservoirs

106

4.3.1. Microscopic characteristics

19

1.0

18.7

62.8

2.5

43.0

26.3

54.5

46.6

70.5

0.84

ACCEPTED MANUSCRIPT

108

AC C

EP

TE D

M AN U

SC

RI PT

107

109

Fig.5. Typical reservoirs images of lacustrine shale A:

Shells were stained using alizarin red to highlight the carbonate

110

minerals; the shells are well developed and preserved in terrestrial inputs and matrices. B: Greenish-yellow fluorescence.

111

Matrices with large amounts of terrestrial inputs between the shell grains have strong fluorescence intensity, indicating

112

effective reservoirs with good oiliness in the matrices. C: Lacustrine shale core sample with well-developed calcareous

113

interlayers, indicating a changeable sedimentary environment. D: Quartz particles with well-developed interparticle pores

20

ACCEPTED MANUSCRIPT

filled with clay minerals; quartz and clay minerals have well-developed intraparticle pores. E: Quartz particles have

115

well-developed intraparticle pores and interparticle pores; calcite minerals fill the storage spaces between quartz particles

116

(calcite cementation). F: Intraparticle pores located along the cleavage planes of clay particles; interparticle pores between

117

clay and carbonate minerals. G: OM with well-developed pores is accompanied by clay minerals. H: Calcite minerals fill the

118

spaces between particles (calcite cementation); quartz has well-developed intraparticle pores. I: Brittle minerals with

119

well-developed micro-cracks that are mainly filled by clay minerals. J: Model of different types of pores. K: Model of

120

carbonate mineral (e.g., calcite) diagenesis that results in reduced storage spaces.

SC

M AN U

121

RI PT

114

According to the International Union of Pure and Applied Chemistry (IUPAC) classification scheme, the pore types in unconventional reservoirs are divided by pore size into micropores, mesopores,

123

and macropores (<2 nm, 2-50 nm, and >50 nm, respectively) (IUPAC, 1994; Thommes et al., 2015).

124

Thin section observations of the shale indicate that the matrices fluoresce intensely (Fig. 5B). Nanoscale

125

reservoirs were examined using FE-SEM with argon-ion milling technology, and the minerals were

126

identified using EDS. Similar to other lacustrine shales, the pores and fractures in the Da’anzhai shale

127

are well developed and fall into four categories.

EP

(1) Interparticle pores: these pores represent the largest number of pores, and they have good

AC C

128

TE D

122

129

connectivity and form effective pore networks. The diameters generally range between 5 nm and 90 µm.

130

Many types of interparticle pores were observed, including pores between grains (Fig. 5D,E,H,J),

131

crystals and clay platelets (Fig. 5D,F,G,J), and pores at the edges of rigid grains (Fig. 5E,J).

132 133

(2) Intraparticle pores: these pores primarily formed through diagenesis, and some are primary. The diameters mainly range from 2 nm to 10 µm. The primary sub-types in the Da’anzhai shale are

21

ACCEPTED MANUSCRIPT

134

intraplatelet pores within clay aggregates (Fig. 5D,F,J), intercrystalline pores within pyrite framboids,

135

pores within fossil bodies and moldic pores after fossils (Fig. 5D,J). (3) Fracture pores: these pores are not controlled by individual particles and can be regarded as a

RI PT

136 137

type of pore that formed by rock deformation because of tectonic movements, sedimentation, diagenesis,

138

OM hydrocarbon generation and other geological effects (Fig. 5I,J).

(4) OM pores: these pores are found within OM (Fig. 5L) and are usually irregular, bubble-like, and

SC

139

oval. They generally range in diameter from 2 to 500 nm. Although most OM pores are incorrectly

141

regarded as isolated pores in two-dimensional planes, they are interconnected in three-dimensional space

142

and facilitate the development of reservoirs.

143

M AN U

140

Similar to other lacustrine and marine shale reservoirs, the Da’anzhai shale commonly contains well-developed pores and fractures, particularly intraparticle and interparticle pores (Curtis et al., 2002;

145

Chalmers et al., 2012; Milliken et al., 2013; Nie et al., 2015; Wang et al., 2015a; Zhou and Kang, 2016;

146

Jiang et al., 2017; Zhang et al., 2017). Especially, due to their rich terrestrial clastic mineral content, the

147

interparticle pores of lacustrine shales are better developed than those of marine shales. Most of the

148

intraparticle pores are nanoscale pores (<200 nm) (Fig. 5D,E,F) and the OM also contains nanopores

149

(Fig.5G). Microcracks can develop due to the high brittle mineral

150

However, these cracks are often filled by calcite and clay minerals due to the high clay content (42.4%)

151

(Fig.5I,J).

AC C

EP

TE D

144

22

content for the Da’anzhai shale.

ACCEPTED MANUSCRIPT

4.3.2. Nanoscale reservoirs characterized by LTNA

153

The LTNA test is the most widely used method of characterizing several parameters of nanoscale

154

pores in porous media (micropores and mesopores), includingVBJH, SBET, DA, ad/desorption isotherms

155

and PSD (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015; Wang et al., 2015c; Zhou et al., 2016;

156

Zhang et al., 2017).

157

Table 4 Comparisons of primary pore structure parameters measured by LTNA between the Da’anzhai shale and other

158

lacustrine shales

159

3

4

5

6

7

160

SBET (m2/g)

Sichuan

Jurassic Da'anzhai

1.75-10.69

Basin

Member

(6.92)

Upper Triassic

0.25-4.39

Chang 7 Member Xujiahe

Ordos Basin Sichuan Basin

3

(cm /100g)

D (nm)

1.19-4.10

7.21-24.71

(2.97)

(14.11)

0.02-0.90

6.67-18.40

(2.62)

(0.49)

(11.23)

0.62-14.70

0.22-2.97

4.91-13.93

(7.76)

(1.41)

(6.97)

Upper Triassic

1.46-5.24

0.72-1.40

9.40-25.00

Chang 7 Member

(2.94)

(1.01)

(16.40)

Bohai Bay

Shahejie

0.40-12.45

0.08-1.60

4.58-41.36

Basin

Formation

(3.62)

(0.62)

(16.84)

0.39-34.06

0.13-4.18

3.16-16.48

(10.92)

(1.59)

(9.56)

Upper Triassic

1.10-1.90

0.69-1.09

15.4-23.4

Chang 7 Member

(1.57)

(0.83)

(19.27)

Ordos Basin

Songliao

Formation

VBJH

TE D

2

strata

EP

1

location

Upper Cretaceous Qingshankou

AC C

NO.

M AN U

SC

RI PT

152

Basin

Ordos Basin

Formation

161

23

hysteresis

loops types

references

H3

H3 H3 H3 H3,H4

H2,H3

H3,H4

Jiang et al., 2016b Peng et al., 2016 Fu et al., 2015 Zhang et al., 2016 Wang et al., 2015b Yang et al., 2017

RI PT

ACCEPTED MANUSCRIPT

162

164

Fig. 6. PSD curves obtained from the low temperature N2 adsorption isothermswith BJH pore sizes.

SC

163

The VBJH value of the Da’anzhai shale (2.97 cm3/100 g) is slightly higher than the values in other shales, and the SBET value (6.92 m2/g) is nearly equal to those of the other shales (Table 4). These values

166

indicate that the Da’anzhai shale can also have significant nanoscale storage space.

M AN U

165

To investigate the distribution and contributions of each pore type, the VBJH with respect to the PSD

168

was determined from the adsorption branches of the BJH model (Wang et al., 2015b; Jiang et al., 2016b;

169

Zhou et al., 2016; Zhang et al., 2017). The PSD plots show that pores with diameters between 3 and 5

170

nm make up the largest proportion of the total pore volume and that the proportions decrease with

171

increasing pore size, which is similar to

EP

some lacustrine shales (Fig. 6) (Jiang et al., 2016b).

AC C

172

TE D

167

24

173 174

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Fig. 7. Comparisons between the nanopore structures of the Da’anzhai shale and other lacustrine shales. A, B, C: the Da’anzhai shale. D: Lacustrine shale of the Chang 7 Member in the Ordos Basin, China (Fu et al., 2015). E: Lacustrine shale

176

of the Songliao Basin, China (Wang et al., 2015b). F: Lacustrine shale of the Yanchang Formation in the Ordos Basin, China

177

(Jiang et al., 2016b; Wang et al., 2015c). G: Lacustrine shale of the Shahejie Formation in the Dongying Sag, China (Liu et

178

al., 2017). H: Lacustrine shale of the Shahejie Formation in the Zhanhua Sag, China. I: Lacustrine shale of the 5

179

Member of the Xujiahe Formation of Upper Triassic in Sichuan Basin, China (Peng, 2016).

180

AC C

EP

175

th

LTNA isotherms and hysteresis patterns can be used to investigate the physisorption mechanisms

181

and structural characteristics of shale pores (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015;

182

Zhou et al., 2016; Zhang et al., 2017). Adsorption isotherms are divided into five types (Types I-VI), and 25

ACCEPTED MANUSCRIPT

hysteresis loops are divided into four types (Types H1-H4) (Sing et al., 1985; IUPAC, 1994; Thommes et

184

al., 2015). Following this classification, the isotherms of the Da’anzhai shale samples are similar and

185

belong to Type П, and the hysteresis loops are characterized by Type H3 hysteresis loops that indicate

186

wedge-shaped pores (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015). However, compared to

187

other lacustrine shales, the Da’anzhai shale has relatively weak hysteresis loops and wide overlapping

188

ranges of adsorption and desorption isotherms (Fig. 7), which indicate many dead-end pores with highly

189

complex structures (De Boer, 1958; Nie et al., 2015).

M AN U

SC

RI PT

183

190

4.3.3. Factors that influence pore structure

191

SBET and VBJH are the two most basic LTNA parameters, and they can directly indicate the storage spaces of nanoscale reservoirs. Fig. 8a shows that the TOC values are positively correlated with SBET

193

and VBJH. However, the correlation is not strong because of the relatively low TOC values (average

194

0.97%). Fig. 8b shows that the lacustrine clay content is positively correlated with SBET and VBJH; this

195

correlation is mainly caused by the high clay content (42.4%) and well-developed intraparticle pores

196

(e.g., intraplatelet pores within clay aggregates) and interparticle pores (e.g., pores at the edges of

197

particles, crystals and clay platelets) (Figs. 5F, L). The quartz content is also positively correlated with

198

SBET and VBJH (Fig. 8c). Quartz particles in lacustrine shales usually occur within mineral matrices and

199

have well-developed interparticle pores at the edges of other minerals (e.g., clay and carbonate minerals).

200

The quartz particles also have well-developed nanoscale intraparticle pores. Additionally, the quartz

201

particles have no obvious cementation (Fig.5D). Conversely, the carbonate content is negatively

202

correlated with SBET and VBJH (Fig. 8d). Carbonate minerals, particularly calcite, can significantly affect

AC C

EP

TE D

192

26

ACCEPTED MANUSCRIPT

reservoir quality and generally has adverse effects (Wang et al., 2016).

Carbonate minerals (e.g.,

204

calcite) are easily dissolved and redistributed and can be formed in primary pores by various

205

mechanisms during the diagenesis processes (Figs. 5E, H, K) (Huang et al., 1984; Wang et al., 2016).

206

Based on qualitative observations using SEM, polarized light microscopy and fluorescence microscopy,

207

lacustrine carbonate minerals can effectively reduce the available storage space by diagenesis processes

208

including cementation, compaction, pressure-solution, recrystallization and replacement (Figs. 5E, H, K)

209

(Huang et al., 1984; Wang et al.,

210

2016).

AC C

211

EP

TE D

M AN U

SC

RI PT

203

Fig. 8. Correlation plots between nanopore storage space indicators (SBET and VBJH) and values of TOC and minerals.

27

ACCEPTED MANUSCRIPT

4.4. Characterization and controlling factors of OM in lacustrine shale

213

4.4.1. Characterization of the OM content

214

Water conditions play an important role in the formation of source rocks. Lacustrine shales

RI PT

212

primarily form in semi-deep to deep lakes with fresh to low salinity water, whereas marine shales

216

primarily form in reductive, salty, weakly alkaline and low energy water. These features give rise to

217

different source rock characteristics (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker,

218

2002; Chen et al., 2015; Yang et al., 2015; Jiang et al., 2016a). Due to these differences, we evaluate the

219

lacustrine OM using a terrestrial oil generation theory that was derived from studies of continental

220

petroliferous basins in China and plays a special role in global petroleum theories (Huang et al., 1984;

221

Lu and Zhang, 2008) (Table 5).

222

Table 5 Factors used to evaluate OM in lacustrine source rock (Huang et al., 1984; Lu and Zhang, 2008).

M AN U

TE D

lake types

TOC(%)

freash~brackish water salt water~high salt water

chloroform bitumen ‘A’ (%) S1 +S2 (mg/g)

non-source rock <0.4

EP

index

AC C

223

SC

215

<0.2

<0.015

poor

source rock medium good

very good

0.4~0.6

>0.6~1.0

>1.0~2.0

>2.0

0.2~0.4

>0.4~0.6

>0.6~0.8

>0.8

0.015~0.050 <0.050~0.100 >0.100~0.200 <2

2~6

>6~20

>0.200 >20

The mean

224

values (and ranges) of the shale’s kerogen macerals consist of 32% exinites (29-35%), 27%

225

sapropelinites (23-32%) and 24% vitrinites (18-26%), which indicates that the OM is primarily derived

226

from plankton and bacteria. The δ13CPDB values (-29.1 to -29.5‰, average -29.3‰) support this finding

227

(Table 7) (Huang et al., 1984; Lu and Zhang, 2008). KI is an effective indictor of OM type. The OM

228

types (I, II1, II2, and III) can be identified using the KI ranges >>80-100, 80-40, 40-0 and <<0 to -100, 28

ACCEPTED MANUSCRIPT

respectively (Huang et al., 1984; Lu and Zhang, 2008). The KI values for the Da’anzhai shale vary from

230

2 to 15 with an average of 8.7, which indicates that the OM is chiefly sapropelic-humic type II2 and has

231

significant potential to generate hydrocarbons. The mean values (and ranges) of the organic elements C,

232

H, and O make up 46.59% (31.12-61.41%), 3.53% (2.09-4.66%) and 1.56% (0.35-1.88%) of the shale’s

233

OM, respectively. The high H/C ratios (0.79-1.05, average 0.90) and low O/C ratios (0.10-0.30, average

234

0.17) (Fig. 9a) indicate that the OM is type II and has relatively high hydrocarbon potential (Huang et al.,

235

1984; Lu and Zhang, 2008). In addition, the Tmax-HI index (Fig. 9b), which is an important evaluation

236

index for the OM type based on Rock-Eval pyrolysis, gives a similar result (Espitaliéet al., 1985; Lu and

237

Zhang, 2008; Misch et al., 2016). Generally, the Da’anzhai shale has similar OM types with other

238

lacustrine and marine shales and the dominant type is of type II (Table 1), which can support the

239

Da’anzhai shale hydrocarbon generation capacity.

240

The TOC values of the Da’anzhai shale vary from 0.11 to 2.18% (average 0.97%), these values are

241

generally lower than those measured in the world's major lacustrine and marine shales (Table 1 and Fig.

242

9). The Ro values of the Da’anzhai shale are similar or even higher than those of other lacustrine and

243

marine shales, which indicate that the OM is highly mature. This could also be deduced from the Tmax

244

values, which vary from 428 to 500 ℃ with an average of 449 ℃ (Table 6). Table 6 Analysis results for

245

Rock-Eval, Chloroform Bitumen ‘A’ and TOC.

AC C

EP

TE D

M AN U

SC

RI PT

229

29

ACCEPTED MANUSCRIPT

1.15

0.50 3.97 0.16 1.03 0.71 2.12 1.09 1.28 0.55 1.85 0.94 0.26 1.32 0.05 1.83

OPI 0.14 0.34 0.12 0.12 0.24 0.07 0.21 0.10 0.11 0.48 0.14 0.29 0.10 0.35 0.04

G10-2 G4-3 J53-1 Li-1 Li-2 P1-1

maceral (%) 13 exinite vitrinite inertinite δ CP DB Ro (%) KI (‰) sporo normal amorphinite resinite sporinite subtotal fusinite pollinite vitrinite 0.73-1.17 27 4 4 26 34 24 15 11 -29.2 (0.95) 1.13-1.65 27 3 2 30 35 18 20 11 -29.5 (1.43) 1.19-1.47 32 4 4 22 30 22 16 15 -29.2 (1.32) 0.78-1.08 28 6 6 19 31 26 15 9 -29.1 (0.95) 0.95-1.25 27 6 2 21 29 26 18 4 -29.5 (1.17) 0.88-1.28 23 5 6 22 33 26 18 2 -29.2 (1.02) sapropelinite

AC C 249 250

HI (mg/g) 60.60 298.77 40.65 117.53 49.69 142.48 93.29 119.80 79.33 116.50 61.96 43.18 112.14 97.82 159.28

PC (%) 0.05 0.50 0.02 0.10 0.08 0.19 0.12 0.12 0.05 0.30 0.09 0.03 0.12 0.01 0.16

D (%) 5.95 37.51 3.98 11.13 5.48 12.79 9.86 11.17 7.48 18.58 6.05 5.18 10.34 14.22 13.91

RI PT

0.08 2.03 0.02 0.14 0.22 0.17 0.29 0.15 0.07 1.70 0.16 0.11 0.14 0.03 0.09

S1 +S2 (mg/g) 0.58 6.00 0.18 1.18 0.94 2.29 1.38 1.43 0.62 3.55 1.09 0.38 1.46 0.09 1.92

SC

0.009 0.006 0.007 0.004 0.009 0.009 0.008 0.009 0.006 0.006 0.008 0.006 0.006 0.008 0.009

Tmax (℃) 456 444 461 439 428 443 441 446 445 452 453 448 444 500 439

HCI (mg/g) 11.04 153.11 7.30 16.61 16.33 11.67 25.49 14.82 10.77 107.39 10.88 19.28 12.41 73.45 8.27

Table 7 Results for Ro, macerals, carbon isotopes and organic elements.

code

248

0.24 0.11 0.27 0.01 0.40

S0 (mg/g) S1 (mg/g) S2 (mg/g)

TE D

247

0.02 0.25 0.42 0.40 0.35 0.31 0.16

TOC (%) 0.82 1.33 0.40 0.88 1.43 1.49 1.17 1.07 0.69 1.59 1.51 0.61 1.18

EP

246

G10-1 G10-2 G4-1 G4-2 G4-3 J61-4 Li-1 Li-2 Li-3 P1-1 P1-2 P1-4 S2-2 S2-5 XI44-2

chloroform bitumen A(%) 0.03

M AN U

code

Fig.9. Discrimination diagrams for the OM type.

30

organic elements C (%) H (%) O (%)

H/C

O/C

39.89

3.07

0.75

0.92

0.1

46.31

4.07

1.88

1.05

0.3

49.26

3.87

1.37

0.94

0.2

31.12

2.09

0.35

0.81

0.1

61.41

4.66

1.62

0.91

0.2

51.54

3.39

0.96

0.79

0.1

ACCEPTED MANUSCRIPT

251

Terrestrial oil-generation theory (Huang et al., 1984; Lu and Zhang, 2008) suggests the following conclusions regarding the lacustrine shale:

253

medium to very good, even though the values of S1+S2 (genetic potential) are poor to medium (Tables 5

254

and 6). (2) Because of the strong hydrocarbon expulsion of the high maturity OM (Ro>0.5-0.7%), the

255

evaluation criteria for high maturity OM should be reduced. The shale with high Ro values (0.95-1.43%)

256

is highly mature. (3) Lacustrine carbonate can also generate hydrocarbons, which is an important feature

257

but easy to overlook. The petroleum generation threshold of carbonate (TOC: 0.12-0.4%) is lower than

258

that of mudstone (TOC: 0.4-0.5%), which means that carbonates can produce hydrocarbons at lower

259

TOC values. Thus, the low TOC values of the high carbonate shale may not indicate a low hydrocarbon

260

generation capacity. Additionally, carbonate interlayers in shale can play a role similar to that of seal

261

rocks, which are conducive to the preservation of oil and gas. In summary, the results of this study

262

indicate that the Da’anzhai shale has some hydrocarbon generation capacity, which is also supported by

263

the clear fluorescence of the thin sections (Figs. 5B, C).

SC

M AN U

TE D

Based on microscopic observations of the Da’anzhai and lacustrine shales, the OM occurrence

EP

264

(1) The TOC and chloroform bitumen ‘A’ values are

RI PT

252

forms are generally similar to other lacustrine shales; they are mainly dispergated in shale matrices

266

(Fig.5) (Morad et al., 2010; Zou et al., 2011; He et al., 2012; Luo et al., 2013; Tian et al., 2014; Dang et

267

al., 2015; Xu et al., 2017). However, these lacustrine shales have more amorphous forms than marine

268

shales because of their different water conditions and biodegradation (Huang et al., 1984; Lu and Zhang,

269

2008; Tian et al., 2014; Xu et al., 2017). In addition, lacustrine OM is distributed in concentric belts

270

around the lake basin’s center, and OM abundance increases with lake depth. Thus, OM is generally

AC C

265

31

ACCEPTED MANUSCRIPT

more abundant in moderately deep to deep lake facies than in shallow lake facies, which can be

272

supported by the wide TOC range of lacustrine shales (Table 1 and Fig. 9) (Huang et al., 1984; Lu and

273

Zhang, 2008; Tian et al., 2014; Xu et al., 2017).

4.4.2. Factors that influence the OM content

Fig. 10. Correlation plots between TOC and V/Cr, V/Sc, Ni/Co and V/(V+Ni) values.

EP

276

AC C

275

TE D

M AN U

SC

274

RI PT

271

32

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

277 278

Fig. 11. Correlation plots between TOC values and minerals.

279

TOC is commonly used to represent OM. As discussed in Section 4.1, the reducibility is positively correlated with the V/Cr, V/(V+Ni), and V/Sc ratios and negatively with the Ni/Co ratio. In addition,

281

pyrite is an important mineral in organic-rich sediments and is useful for reconstructing

282

paleoenvironmental conditions. A high pyrite content reflects a stable, high-salinity environment with

283

strong reducibility (Leventhal, 1983). The TOC values of the Da’anzhai shale are positively correlated

284

with the V/Cr, V/(V+Ni) and V/Sc ratios (Figs. 10a, b, c) and the pyrite content (Fig. 11a) and correlated

285

negatively with the Ni/Co ratio (Fig. 10d). These values indicate that paleo-redox conditions had a

286

significant effect on the formation and preservation of OM and that the strongly reducing conditions

287

were conducive to the preservation of OM (Zeng et al., 2015; Chen et al., 2016; Xu et al., 2017).

288 289

AC C

EP

TE D

280

Similar to other lacustrine shales, the TOC values of the Da’anzhai shale are positively correlated with the clay mineral content (Fig. 11b), which is due to the OM’s strong sorption ability of clay and its 33

ACCEPTED MANUSCRIPT

relatively high content for the lacustrine shales (Huang et al., 1984; Gingele et al., 2001; Dang et al.,

291

2015; Xu et al. 2017). The lacustrine clay minerals are small and rich in intraparticle and interparticle

292

pores, which creates storage space for OM (Figs.5D, F, G and J). The OM observed within the clay

293

minerals primarily has complex, dispergated and residual biological shapes, which indicate that the OM

294

and clay minerals are in close contact (Figs. 5G and J) (Huang et al., 1984; Lu and Zhang, 2008; Tian et

295

al., 2014; Xu et al., 2017). Conversely, the carbonate content is negatively associated with OM (Fig.

296

11c), which is similar to some lacustrine shales (Huang et al., 1984; Jarvis et al., 2001; Tian et al., 2014;

297

Wang et al., 2016). Carbonate minerals fill the pores via cementation, compaction, pressure-solution,

298

recrystallization and replacement , which effectively reduces the storage space for OM (Figs. 5E, H, K)

299

(Huang et al., 1984; Wang et al., 2016).. Furthermore, the high carbonate mineral content in the

300

lacustrine facies reflects relatively shallow water environments, such as the shore-shallow lake face and

301

shallow lacustrine facies, which are not favorable for the generation or storage of OM (Huang et al.,

302

1984; Wang et al., 2016; Xu et al., 2017). There is a weak correlation between the quartz content and

303

TOC values (Fig. 11d). The correlations between quartz and OM are not uniform; differences in the

304

correlation between quartz and TOC contents is related to the depositional environment of shale

305

reservoirs (Zeng et al., 2014; Liu et al., 2015; Ross and Bustin, 2007; Chalmers et al., 2012b; Tian et al.,

306

2013). Negative correlations are caused by the terrestrial quartz inputs (Zeng et al., 2014; Liu et al.,

307

2015), while the positive correlations are related to biogenic quartz (Ross and Bustin, 2007; Chalmers et

308

al., 2012b; Tian et al., 2013). Therefore, the quartz in the lacustrine shales may have both related to

309

biogenic and terrestrial origins. This could lead to the weak correlations between quartz and OM.

AC C

EP

TE D

M AN U

SC

RI PT

290

34

ACCEPTED MANUSCRIPT

310 311

Based on the analyses in Figs. 10 and 11, the low TOC values of the Da’anzhai lacustrine shales were primarily caused by the oxygen paleoenvironment and carbonate minerals. The relatively shallow

313

and turbulent lake water formed a high oxygen paleoenvironment with greater inputs of terrestrial

314

nutrients, which was conducive to the breeding of shell organisms and the formation of carbonate

315

minerals but not to the preservation of OM. In summary, the oxygen paleoenvironment and high

316

carbonate mineral content are primarily related to the unique lacustrine facies, which was not conducive

317

to the formation of OM.

320

SC

M AN U

319

4.5. Relationships between the mineralogy, reservoirs, OM and paleoenvironmental conditions of lacustrine shales

The mineral content of lacustrine shales is mainly related to the paleoenvironment. Strong

TE D

318

RI PT

312

paleo-weathering can lead to a greater number of weathering products, which is favorable to the

322

development of terrestrial minerals in lacustrine facies. Humid (rainy) paleoclimates with more water

323

can improve river transportation and reduce lake salinity, and low salinity paleoenvironments inhibit the

324

formation of authigenic carbonate. In addition, fine-grained terrestrial minerals, especially clay minerals,

325

inhibit the formation of authigenic carbonate (Leventhal, 1983; Huang et al., 1984; Bohacs et al., 2000;

326

Chi et al., 2003; Fu et al., 2015). High hydrodynamic conditions can effectively remove fine-grained

327

terrestrial minerals, which favors the formation of authigenic carbonate.

328 329

AC C

EP

321

The lacustrine sedimentary environment is complex and the lacustrine shale formation is easily influenced by the paleoenvironment. The two most important differences are as follows: first, the

35

ACCEPTED MANUSCRIPT

lacustrine paleoenvironments are complex and varied with the short detention time and strong variation

331

of water level (Huang et al., 1984; Chen et al., 2015). (Figs. 1 and 5C, Tables 2 and 6). Second, lakes are

332

close to the terrestrial sources with less water, which lead to rich terrestrial inputs with efficient impact

333

on shale (Figs. 1 and 5A,B, Tables 3 and 6) (Huang et al., 1984; Chen et al., 2015) (Huang et al., 1984;

334

Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002). Lakes are generally located near paleo-uplift

335

environments with large topographical differences, and terrestrial debris can be transported directly into

336

a lake basin over short distances. These features result in rich, proximal and variable sources for

337

lacustrine facies. (Figs. 1, 5A,5B; Table 3). The above two characteristics can effectively lead to the

338

complex and changeable mineral contents of lacustrine shales (Table 3; Fig. 4). In addition, there are

339

some other common points of the lacustrine shales including relatively small thickness of monolayer,

340

more siliceous or calcareous intercalations, and poor comparability with strong heterogeneity (Table 1),

341

which are also the results of the two characteristics of lacustrine sedimentation environments (Tian et al.,

342

2014; Wang et al., 2015a; Zhu et al., 2016; Xu et al., 2017).

SC

M AN U

TE D

Similar to other lacustrine shales, the storage space of the Da’anzhai lacustrine shale is directly

EP

343

RI PT

330

controlled by minerals (Fig. 12). Because of the large amount of terrestrial clastic minerals, the

345

interparticle pores of lacustrine shales are well-developed. Developed interparticle pores and intraparicle

346

pores can guarantee a considerable amount of storage space. In addition, numerous siliceous or

347

calcareous interlayers in lacustrine shale formations can result from varying changeable

348

paleoenvironments with rich terrestrial inputs. These interlayers can increase the brittleness index and

349

facilitate the development of cracks.

AC C

344

36

ACCEPTED MANUSCRIPT

350

Compared to other shales, the Da’anzhai lacustrine shale has a relatively low OM content (Table 6), which is mainly the result of more variable paleoenvironments with high oxygen content (Fig. 12)

352

(Huang et al., 1984; Lu and Zhang, 2008; Tian et al., 2014; Xu et al., 2017). Different minerals have

353

different effects on storage spaces of OM, and minerals in lacustrine shales are mainly controlled by

354

paleoenvironments (Fig. 12). In addition, the relatively low OM content of the Da’anzhai lacustrine

355

shale is the reason for OM weak effects on storage spaces (Fig. 8a).

SC

RI PT

351

In summary, minerals, reservoirs, and OM in lacustrine shales are closely associated with each

M AN U

356

other, and they were all controlled by the paleoenvironment which was the dominant factor and the link

358

between these parameters (Fig.

12).

359

Fig.12 Relationships between mineralogy, reservoirs, OM and paleoenvironment for lacustrine shale

362

4. Conclusions

363 364

AC C

361

EP

360

TE D

357

1. There are some similarities between the Da’anzhai lacustrine shale and other lacustrine shales. (1)

365

The two most important differences are that lacustrine shales formed in more complex and variable

366

paleoenvironments and have more abundant terrestrial inputs than marine shales. (2) The mineral

367

compositions of lacustrine shales are variable and complex, and they have abundant terrestrial minerals. 37

ACCEPTED MANUSCRIPT

In addition, lacustrine shales have numerous siliceous or calcareous interlayers. (3) The interparticle

369

pores and interlayer pores of lacustrine shales are well-developed. Unlike other lacustrine shales, the

370

Da’anzhai shale has relatively weak hysteresis loops and wide overlapping ranges of adsorption and

371

desorption isotherms, indicating a greater number of dead-end pores with more complex structures. In

372

addition, the Da’anzhai shale has lower TOC values with more amorphous OM forms.

2. The minerals, reservoirs, and OM for lacustrine shales are closely associated with each other;

SC

373

RI PT

368

these parameters are all controlled by the paleoenvironmental conditions in which the shales formed. (1)

375

Highly variable lacustrine paleoenvironment lead to variable mineral contents and abundant terrestrial

376

inputs give rise to a high terrestrial mineral content in lacustrine shales, which can also create a higher

377

number of siliceous or calcareous intercalations relative to marine shales. (2) Terrestrial minerals (e.g.,

378

clay and quartz) can effectively improve the storage spaces, but the authigenic carbonate minerals (e.g.,

379

calcite) have the opposite effect. Well-developed interparticle pores in lacustrine shales are caused by

380

abundant particles. (3) Clay content has a positive effect on the formation and preservation of OM, while

381

higher oxygenation and carbonate contents have negative effects. Lakes generally have relatively

382

shallow water and a variable paleoenvironment, which leads to a higher oxygen content and reduces the

383

OM abundance of lacustrine shales. The weak influence of OM on nanoscale storage spaces is due to the

384

low TOC values of the lacustrine shales.

TE D

EP

AC C

385

M AN U

374

3. Compared to other lacustrine shales, the Da’anzhai shale also demonstrates the potential for

386

unconventional oil and gas exploration. (1) The lacustrine shale has well-developed storage spaces due

387

to the high number of interparticle pores, intraparticle pores, and fractures. Although the lacustrine shale

38

ACCEPTED MANUSCRIPT

has more dead-end pores and larger pores with more complex microstructures, the types of isotherms

389

(Type П) and hysteresis loops (Type H3) are generally similar to those of lacustrine shales. The values of

390

VBJH (1.19-4.10 m3/100 g, average 2.97 m3/100 g) and SBET (1.75-10.69 m2/g, average 6.92 m2/g) are

391

similar to other shales, indicating a high amount of storage space in nanoscale reservoirs. (2) Abundant

392

terrestrial minerals in lacustrine facies can support a greater number of interparticle pores and

393

intraparticle pores. Numerous siliceous or calcareous interlayers in lacustrine shales can also improve

394

the brittleness index, facilitating the development of fractures and creating a seal for oil and gas. (3)

395

Although the TOC values of the Da’anzhai shale are lower than other lacustrine shales, the OM still has

396

several advantages according to the evaluation criteria for continental source rocks, including moderate

397

values of chloroform bitumen ‘A’ and TOC, a high level of OM maturity and a favorable kerogen type

398

(II).

399

Acknowledgments

TE D

M AN U

SC

RI PT

388

Our study is supported by the Institute Program of PetroChina Research Institute of Petroleum

401

Exploration and Development (Grant 2016yj01), and the National Natural Science Foundation of China

402

(Grants 41272137 and 41572117). The authors also sincerely appreciate the support from the China

403

Scholarship Council (CSC).

404

Nomenclature

405

TOC

406

Ro

Vitrinite reflectance

407

Tmax

Pyrolysis Temperature at Maximum Hydrocarbon Generation

AC C

EP

400

Total organic carbon, %

39

ACCEPTED MANUSCRIPT

SBET

Brunauer-Emmett-Teller specific surface area

409

VBJH

Barrett-Joyner-Halenda total pore volume

410

DA

Average pore diameter

411

OM

Organic matter

412

TE

Trace element

413

XRD

X-ray diffraction

414

FE-SEM

Field emission scanning electron microscopy

415

LTNA

Low temperature N2 adsorption

416

EDS

Energy dispersive spectrometry

417

S0

Gaseous hydrocarbons, mg/g

418

S1

Free hydrocarbons, mg/g

419

S2

Pyrolysis of hydrocarbons, mg/g

420

ICP-MS

Inductively coupled plasma-mass spectrometry

421

REE

Rare earth element, ppm

422

ΣHREE

Total content of heavy rare earth elements

423

ΣLREE

Total content of light rare earth elements

424

PSD

Pore size distribution

425

OPI

426

HI

Hydrogen index; HI=(S2/TOC)×100%, mg/g

427

PC

Effective carbon content; PC=0.083×(S0+S1+S2), %

AC C

EP

TE D

M AN U

SC

RI PT

408

Oil production index; OPI=S1/(S0+S1+S2)

40

ACCEPTED MANUSCRIPT

D

Degradation rate; D=(PC/TOC) ×100%, %

429

HCI

Hydrocarbon index; HCI=(S0+S1)/ TOC×100, mg/g

430

KI

Kerogen type index; KI=Sapropelinite (%)×1+Exinite (%)×0.5+Vitrinite

RI PT

428

(%)×(-0.75)+Inertinite (%)×(-1) (Huang et al., 1984; Lu and Zhang, 2008)

432

References

433

Barrett, E.P., Joyner, L.G., Halenda, P. P., 1951. The determination of pore volume and area distributions in porous

435

substances. I. Computations from nitrogen isotherms. J. Am. Chem. Soc. 73(1), 373-380.

M AN U

434

SC

431

Bohacs, K.M., Carroll, A.R., Neal, J.E., Mankiewicz, P.J., 2000.Lake-basin type, source potential, and

436

hydrocarbon character: an integrated sequence-stratigraphic-geochemical framework. Lake basins through

437

space and time: AAPG Studies in Geology 6, 3-34.

441 442

TE D

440

309-319.

Bowker, K.A., 2007. Barnett shale gas production, Fort Worth Basin: issues and discussion. AAPG Bull. 91(4), 523-533.

EP

439

Brunauer, S., Emmett, P.H., Teller, E.J., 1938. Adsorption of gases in multimolecular layers. Am. Chem. Soc. 60,

Chalmers, G.R., Bustin, R.M., Power, I.M., 2012. Characterization of gas shale pore systems by porosimetry,

AC C

438

443

pycnometry, surface area, and field emission scanning electron microscopy/transmission electron

444

microscopy image analyses: Examples from the Barnett, Woodford, Haynesville, Marcellus, and Doig units.

445

AAPG Bull. 96(6), 1099-1119.

41

ACCEPTED MANUSCRIPT

Chi, G., Giles, P.S., Williamson, M.A., Lavoie, D., Bertrand, R., 2003. Diagenetic history and porosity evolution

447

of Upper Carboniferous sandstones from the Spring Valley #1 well, Maritimes Basin, Canada-implications

448

for reservoir development. J. Geochem. Explor. 80, 171-191.

449

RI PT

446

Chen, C., Mu, C.L., Zhou, K.K., Liang, W., Ge, X.Y., Wang, X.P., Wang, Q.Y., Zheng, B.S., 2016. The geochemical characteristics and factors controlling the organic matter accumulation of the Late

451

Ordovician-Early Silurian black shale in the Upper Yangtze Basin, South China. Mar. Pet. Geol. 76, 159-175.

453 454

Chen, L., Liu, C.L., Zhuang, C., Che, X.G., Wu, J., 2009. Rare Earth Element records of the lower Paleogene

M AN U

452

SC

450

sediments in the Sanshui Basin and their paleoclimate implications. Acta Sedimentol. Sin. 27 (6), 1155-1162. Chen, S.J.; Zhang, H.X., Lu, J.G, Yang, Y.M., Liu, C.W., 2015. Controlling factors of Jurassic Da’anzhai Member tight oil accumulation and high production in central Sichuan Basin, SW China. Petrol. Explor. Dev. 42(2),

456

206-214.

TE D

455

457 Curtis, J. B., 2002. Fractured shale-gas systems. AAPG Bull. 86(11), 1921-1938. 458 Dang, W., Zhang, J.C., Huang, X., Li, X.G., Chen, Q., Sun, R., Xue, B., Han, S.B., 2015. Main-controlling

geological factors of gas-bearing property of continental shale gas: a case study of Member 3rd of Shahejie

460

Formation in western Liaohe sag. Acta Petrolei Sin. 36(12), 1516-1530.

AC C

EP

459

461 De Boer, J.H., 1958. The structure and properties of porous materials. In: Proceedings of the tenth symposium of the 462

colston research society held in the University of Bristol, Butterworths, London, pp. 68-94.

463

Espitalié, J., Deroo, G., Marquis, F., 1985. La pyrolyse Rock-Eval et ses applications. Deuxième partie. Revue de

464

l'Institut français du Pétrole 40(6), 755-784.

42

ACCEPTED MANUSCRIPT

Fu, H.J., Wang, X.Z., Zhang, L.X., Gao, R.M., Li, Z.T., Xu, T., Zhu, X.L., Xu, W., Li, Q., 2015. Investigation of

466

the factors that control the development of pore structure in lacustrine shale: A case study of block X in the

467

Ordos Basin, China. J. Nat. Gas. Sci Eng. 26, 1422-1432.

RI PT

465

Gingele, F.X., De Deckker, P., Hillenbrand, C.D., 2001. Late Quaternary fluctuations of the Leeuwin Current and

469

palaeoclimates on the adjacent land masses: clay mineral evidence. Aust. J. Earth Sci. 48(6), 867-874.

470

Gromet, L.P., Haskin, L.A., Korotev, R.L., Dymek, R.F., 1984. The “North American shale composite”: its

471

compilation, major and trace element characteristics. Geochim. Cosmochim. Ac.

SC

468

M AN U

48(12), 2469-2482.

472 Han, S.B., Zhang, J.C., Li, Y.X., Hosrfield, B., Tang, X., Jiang, W.L., Chen, Q., 2013. Evaluation of Lower

Cambrian shale in northern Guizhou Province, South China: Implications for shale gas potential. Energ. Fuel 27,

474

2933-2914.

475

He, F.Q., Zhu, T., 2012. Favorable targets of breakthrough and built-up of shale gas in continental facies in Lower

TE D

473

Jurassic, Sichuan Basin. Petrol. Geol. Exper. 34(3), 246-251. Huang, D.F., Li, J.C., Zhang, D.J., 1984.

477

Evolution and hydrocarbon generation mechanisms of Terrestrial organic matter. Petroleum Industry Press:

478

Beijing, pp. 1-228.

EP

476

AC C

479 IUPAC (International Union of Pure and Applied Chemistry), 1994. Physical Chemistry Division Commission on 480

Colloid and Surface Chemistry, Subcommittee on Characterization of porous solids. Recommendations for the

481

characterization of porous solids (Technical Report). Pure Appl. Chem. 66 (8), 1739-1758.

482

Jarvie, D.M., Hill, R.J., Ruble, T.E., 2007. Unconventional shale-gas systems: The Mississippian Barnett Shale of

483

north-central Texas as one model for thermogenic shale-gas assessment. AAPG Bull. 91, 475-499.

43

ACCEPTED MANUSCRIPT

484 485

Jarvis, I.A.N., Murphy, A.M., Gale, A.S., 2001. Geochemistry of pelagic and hemipelagic carbonates: criteria for identifying systems tracts and sea-level change. J. Geol. Soc. 158(4), 685-696. Jiang, S., Xu, Z.Y., Feng, Y.L., Zhang, J.C., Cai, D.S., Chen, L., Wu,Y., Zhou, Z.S., Bao, S.J., Long, S.X., 2016a.

487

Geologic characteristics of hydrocarbon-bearing marine, transitional and lacustrine shales in China. J. Asian

488

Earth Sci. 115, 404-418.

RI PT

486

Jiang, F.J., Chen, D., Wang, Z.F., Xu, Z.Y., Chen, J., Liu, L., Huyan, Y.Y., Liu, Y., 2016b. Pore characteristic

490

analysis of a lacustrine shale: A case study in the Ordos Basin, NW China. Mar. Pet. Geol. 73, 554-571.

492

M AN U

491

SC

489

Jiang, C., Bryant, S., Daigle, H., 2015. A bundle of short conduits model of the pore structure of gas shale. Unconventional Resources Technology Conference (URTEC).

Jiang, S., Tang, X., Cai, D., Xue, G., He, Z., Long, S., Peng, Y., Gao, B., Xu, Z., Dahdah, N., 2017. Comparison of

494

marine, transitional, and lacustrine shales: A case study from the Sichuan Basin in China. J. Petrol. Sci. Eng.,

495

150, 334-347.

498 499 500 501 502 503

conditions in ancient mudstones. Chem. Geol. 111(1-4), 111-129.

EP

497

Jones, B., Manning, D.A.C., 1994. Comparison of geochemical indices used for the interpretation of palaeoredox

Kimoto, A., Nearing, M.A., Shipitalo, M.J., Polyakov, V.O., 2006. Multi℃year tracking of sediment sources in a

AC C

496

TE D

493

small agricultural watershed using rare earth elements. Earth Surf. Proc. Land. 31(14) 1763-1774. Kimura, H., Watanabe, Y., 2001. Oceanic anoxia at the Precambrian-Cambrian boundary. Geology 29(11), 995-998.

Leventhal, J.S., 1983. An interpretation of carbon and sulfur relationships in Black Sea sediments as indicators of environments of deposition. Geochim. Cosmochim. Ac. 47(1), 133-137.

44

ACCEPTED MANUSCRIPT

504

Li, A., Ding, W.L., He, J.H., Dai, P., Yin, S., Xie, F., 2016. Investigation of pore structure and fractal

505

characteristics of organic-rich shale reservoirs: A case study of Lower Cambrian Qiongzhusi Formation in

506

Malong block of eastern Yunnan Province, South China. Mar.

508

RI PT

507

Pet. Geol. 70, 46-57.

Li, J.Y., 2014. Oil and gas contents and movable oil amounts of shales in 4th member of Shahejie Formation, Lijin sub Sag, Dongying Sag. Petrol. Geol. Exp. 36(3), 365-369.

SC

509 Li, P.P., Hao, F., Guo, X.S., Zou, H.Y., Yu, X.Y., Wang, G.W., 2015. Processes involved in the origin and

accumulation of hydrocarbon gases in the Yuanba gas field, Sichuan Basin, southwest China. Marine and

511

Petroleum Geology 59, 150-165.

512

Lu, S.F., Zhang, M., 2008. Petroleum Geochemistry. Petroleum Industry Press: Beijing, pp. 200-246.

513

Lu, B.X., Zheng, R.C., Liang, X.W., Xu, W.L., 2014. Characteristics analysis of Da’anzhai shale gas (oil)

517 518 519 520

TE D

516

Luo, P., Ji, L.M., 2013. Reservoir characteristics and potential evaluation of continental shale gas. Nat. Gas Geosci. 24(5), 1060-1068.

Liu, Y., Lu, Z.Y., Feng, M.S., Wang, J., Tian, T.H., Chao, J., 2017. Micro-pore characteristics of shale oil

EP

515

reservoirs in eastern Sichuan Basin. Geology in China 41, 1387-1398.

reservoirs of the Shahejie Formation in the Dongying Sag, Bohai Bay Basin. Acta Geol. Sin. 91(3), 629-644.

AC C

514

M AN U

510

Martini, A. M., Walter, L.M., McIntosh, J.C., 2008. Identification of microbial and thermogenic gas components from Upper Devonian black shale cores, Illinois and Michigan basins. AAPG Bull. 92(3), 327-339.

521

Morad, S., Al-Ramadan, K., Ketzer, J.M., De Ros, L.F., 2010. The impact of diagenesis on the heterogeneity of

522

sandstone reservoirs: a review of the role of depositional facies and sequence stratigraphy. AAPG Bull. 94,

523

1267-1309.

45

ACCEPTED MANUSCRIPT

527 528 529 530 531

Milliken, K.L., Rudnicki, M., Awwiller, D.N., Zhang, T., 2013. Organic matter-hosted pore system, Marcellus Formation (Devonian), Pennsylvania. AAPG Bull. 97 (2), 177-200.

RI PT

526

Upper Visean Black Shales in the Dniepr-Donets Basin (Ukraine). Mar. Pet. Geol. 75, 203-219.

Murray, R.W., ten Brink, M.R.B., Jones, D.L., Gerlach, D.C., Russ, G.P., 1990. Rare earth elements as indicators of different marine depositional environments in chert and shale. Geology 18(3), 268-271.

SC

525

Misch, D., Gross, D., Mahlstedt, N., Makogon, V., Sachsenhofer, R.F., 2016. Shale gas/shale oil potential of

Nie, B.S., Liu, X.F., Yang, L.L., Meng, J.Q., Li, X.C., 2015. Pore structure characterization of different rank coals

M AN U

524

using gas adsorption and scanning electron microscopy. Fuel 158, 908-917.

532 Potter, P., Maynard, J., Pryor, W., 1980. Sedimentology of Shale. Springer: New York, pp. 98. 533 Peng, J.L., Zhang, J.C., Wei, P.F., Long, S., Huang, X., Huang, Z.L., Ma, Y.L., Zhang, M.Q., 2016. Characteristics

and Controlling Factors of Micro-pore Structures of the 5th Shale Member of the Xujiahe Formation of Upper

535

Triassic in Sichuan Basin. Geoscience 30 (4), 896-904.

536

Reheis, M.C., 1990. Influence of climate and eolian dust on the major-element chemistry and clay mineralogy of

539 540 541 542 543

EP

538

soils in the northern Bighorn Basin, USA. Catena 17(3), 219-248. Rimmer, S.M., 2004. Geochemical paleoredox indicators in Devonian-Mississippian black shales, central

AC C

537

TE D

534

Appalachian Basin (USA). Chem. Geol. 206(3), 373-391. Ross, G. R., Guevara, S. R., Arribere, M. A., 1995. Rare earth geochemistry in sediments of the Upper Manso River basin, Rio Negro, Argentina. Earth Planet. Sc. Lett. 133(1-2), 47-57. Roy, P. D., Smykatz-Kloss, W., 2007. REE geochemistry of the recent playa sediments from the Thar Desert, India: an implication to playa sediment provenance. Chem. Erde-Geochem. 67(1), 55-68.

46

ACCEPTED MANUSCRIPT

545 546

Schmoker, J. W., 2002. Resource-assessment perspectives for unconventional gas systems. AAPG Bull. 86(11), 1993-1999. Shi, B.H., Yao, J.L., Zhang, Y., Zhang, L., Yang, Y.J., Li, H., 2013. Geologic conditions for hydrocarbon

RI PT

544

547

accumulation in Chang-9 oil reservoir of the Yanchang Formation, Ordos Basin. Oil gas Geol. 34(3),

548

294-300.

Sing, K.S., Everett, D.H., Haul, R.A.W., Moscou, L., Pierotti, R.A., Rouquerol, J., Siemieniewsha, T., 1985.

SC

549

Reporting physisorption data for gas/solid systems with special reference to the determination of surface area

551

and porosity. Pure Appl. Chem. 57, 603-619.

552 553

M AN U

550

Stevens, C.J., Quinton, J.N., 2008. Investigating source areas of eroded sediments transported in concentrated overland flow using rare earth element tracers. Catena 74(1), 31-36.

Thommes, M., Kaneko, K., Neimark, A. V., Olivier, J. P., Rodriguez-Reinoso, F., Rouquerol, J., Sing, K. S., 2015.

555

Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution

556

(IUPAC Technical Report). Pure Appl. Chem. 87(9-10), 1051-1069. Tian, C.T., Ma, S.P., Yang, Y., Meng, J., Xia, Y.Q., 2014. Comparison of hydrocarbon generation conditions

EP

557

TE D

554

between lacustrine and marine carbonate source rocks. Oil Gas Geol. 35(3), 336-341. Vanderaveroet, P.,

559

2000.Miocene to Pleistocene clay mineral sedimentation on the New Jersey Shelf. Oceanol. Acta 23(1),

560

25-36.

561

AC C

558

Veizer, J., Demovic, R., 1974. Strontium as a tool in facies analysis. J.

47

Sediment. Res. 44(1), 93-115.

ACCEPTED MANUSCRIPT

562

Wang, J., Cao, Y.C., Liu, K.Y., Liu, J., Xue, X.J., Xu, Q.S., 2016. Pore fluid evolution, distribution and water-rock interactions of carbonate cements in red-bed sandstone reservoirs in the Dongying Depression, China. Mar.

564

Pet. Geol. 72, 279-294.

565

RI PT

563

Wang, M., Wilkins, R. W., Song, G., Zhang, L., Xu, X., Li, Z., Chen, G., 2015a. Geochemical and geological characteristics of the Es3 L lacustrine shale in the Bonan sag, Bohai Bay Basin, China. In. J.Coal Geol. 138,

567

16-29. Wang, M., Xue, H.T., Tian, S.S., Wilkins, R.W.T., Wang, Z.W., 2015b. Fractal characteristics of

568

Upper Cretaceous lacustrine shale from the Songliao Basin, NE China. Mar. Pet. Geol. 67, 144-153.

569

Wang, Y., Zhu, Y.M., Wang, H.Y., Feng, G.J., 2015c. Nanoscale pore morphology and distribution of lacustrine

570

shale reservoirs: Examples from the Upper Triassic Yanchang Formation, Ordos Basin. J. Energy Chem.

571

24(4), 512-519.

M AN U

SC

566

Wu, Y., Fan, T., Zhang, J., Jiang, S., Li, Y., Zhang, J., Xie, C., 2014. Characterization of the upper Ordovician and

573

lower Silurian marine shale in northwestern Guizhou province of the upper Yangtze block, South China:

574

Implication for shale gas potential. Energy Fuel 28(6), 3679-3687.

577

EP

576

Wright, J., Seymour, R.S., Shaw, H.F., 1984. REE and Nd isotopes in conodont apatite: variations with geological age and depositional environment. Geol. Soc.

Am. Spec. Pap. 196, 325-340.

AC C

575

TE D

572

Xu, Q., Liu, B., Ma, Y., Song, X., Wang, Y., Xin, X., Chen, Z., 2017. Controlling factors and dynamical formation

578

models of lacustrine organic matter accumulation for the Jurassic Da'anzhai Member in the central Sichuan

579

Basin, southwestern China. Mar. Pet. Geol. 86, 1391-1405.

580

Yang, C., Zhang, J., Tang, X., Ding, J., Zhao, Q., Dang, W., Chen, H., Su, Y., Li, B., Lu, D., 2017. Comparative

581

study on micro-pore structure of marine, terrestrial, and transitional shales in key areas, China. Int. J. Coal

48

ACCEPTED MANUSCRIPT

582

Geol. 171, 76-92. Yang, Z., Hou, L.H., Tao, S.Z., Cui, J.W., Wu, S.T., Lin, S.H., Pan, S.Q., 2015. Formation

583

and “sweet area” evaluation of liquid-rich hydrocarbons in shale strata. Petrol. Explor. Dev. 42(5), 609-620. Zeng, S., Wang, J., Fu, X., Chen, W., Feng, X., Wang, D., Song, C., Wang, Z., 2015. Geochemical characteristics,

RI PT

584 585

redox conditions, and organic matter accumulation of marine oil shale from the Changliang Mountain area,

586

northern Tibet, China. Mar. Pet. Geol. 64, 203-221.

Zhang, J.Z., Li, X.Q., Wei, Q., Cao, W.J., Liang, W.L., Wang, Z., Wang, F.Y., 2017. Quantitative characterization

588

of pore-fracture system of organic-rich marine-continental shale reservoirs: A case study of the Upper

589

Permian Longtan Formation, Southern Sichuan Basin, China. Fuel 200, 272-281.

M AN U

590

SC

587

Zhang, Q., Zhu, X.M., Li, C.X., Qiao, L.J.Y., Liu, C., Mei, X.H., Du, H.Y., Lu, J.J., 2016. Classification and quantitative characterization of microscopic pores in organic-rich shale of the Shahejie Formation in the

592

Zhanhua Sag, Bohai Bay Basin. Oil Gas Geol. 37(3), 422-438.

593

TE D

591

Zhou, L., Kang, Z.H., 2016. Fractal characterization of pores in shales using NMR: A case study from the Lower Cambrian Niutitang Formation in the Middle Yangtze Platform, Southwest China. J. Nat. Gas Sci. Eng. 35,

595

860-872.

597 598 599

Zou, C.N., Dong, D.Z., Yang, H., Wang, Y.M., Huang, J.L., Wang, S.F., Fu, C.X., 2011. Conditions of shale gas

AC C

596

EP

594

accumulation and exploration practices in China. Natur. Gas Ind. 31, 26-39. Zhu, T., Long, S.X., Wang F., Peng, Y.M., 2016. Sedimentary models and lithofacies types of lacustrine mud shale in the Sichuan Basin. Natur. Gas Ind. 36(8), 22-28.

600 601

49

ACCEPTED MANUSCRIPT

603 604

605

Second: Manuscript (clean copy)

Geological and geochemical characterization of lacustrine shale: a case study of the Jurassic Da'anzhai Member shale in the central Sichuan Basin, Southwest China

RI PT

602

Qilu Xua,b,c, Bo Liuc,*, Yongsheng Maa, Xinmin Songd, Yongjun Wangd, Zhangxin Chenb a

School of Earth Sciences and Resources, China University of Geosciences(Beijing), Beijing 100083, China

607

b

Chemical and Petroleum Engineering, University of Calgary, Calgary, Alberta T2N1N4, Canada

608

c

Oil and Gas Research Center, Peking University, Beijing 100871, China

609

d

PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China

M AN U

SC

606

610

Abstract:

611

Although numerous progresses have been achieved in characterizing marine shales, studies that are associated with lacustrine shales are limited.

To study lacustrine shales, samples were collected from

613

the Jurassic Da'anzhai Member in the central Sichuan Basin of China, and their mineralogical, reservoir,

614

OM, and paleoenvironmental characteristics were determined, as well as the relationships between them.

615

Analysis of trace elements reveals that the shales formed in paleoenvironments that were oxic to suboxic,

616

dry to humid, had moderate to strong weathering, and were characterized by fresh to salt water

617

conditions. These environments are more variable than those of marine shales. The paleoenvironmental

618

conditions and mineralogy of the shales, particularly the oxic to suboxic paleo-redox conditions, resulted

619

in the relatively low levels (0.11-2.18%, average 0.97%). However, based on evaluation criteria for

620

continental source rocks, the OM is of high quality because of its high level of maturity (Ro: 0.95-1.43;

621

Tmax: 428-500°C) and favorable kerogen type (II2). There are well-developed intraparticle pores,

AC C

EP

TE D

612

50

ACCEPTED MANUSCRIPT

interparticle pores, and microcracks. The SBET (5.42-10.69 m2/g, average 6.92 m2/g) and VBJH (1.19-4.10

623

mL/100 g, average 2.97 mL/100 g) values also indicate good nanoscale storage space. Terrestrial

624

minerals (i.e., quartz and clay) and authigenic carbonate minerals (i.e., calcite) are, respectively,

625

positively and negatively correlated with the nanoscale storage space, .. Small pores (3-5 nm) dominate

626

the nanoscale storage space. The isotherms and hysteresis loops are of Type П and Type H3, respectively,

627

which indicates wedge-shaped pores. However, the hysteresis loops indicate that the lacustrine shale has

628

more dead-end pores and larger pores with more complex microstructures than other lacustrine

629

reservoirs. In general, the Da’anzhai lacustrine shale has the potential for unconventional oil and gas

630

exploration. The lacustrine shale’s mineralogy, reservoirs, and OM are closely related to each other, and

631

their differences are mainly caused by the paleoenvironmental conditions in which the shale formed.

633

SC

M AN U

Keywords: Lacustrine shale; Paleoenvironment; Nanoscale reservoir; Organic matter; Minerals; Mechanisms.

TE D

632

RI PT

622

1. Introduction

635

The high commercial production of shale oil and gas in North America has made shale a focus of

EP

634

exploration in many countries and regions. However, the most widely developed shales were formed in

637

marine systems; few studies have focused on the characteristics of lacustrine shales. Additionally,

638

lacustrine shales are widely distributed in many areas, such as Africa, South America, Southeast Asia

639

and China (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et al., 2015;

640

Jiang et al., 2016a; Yang et al., 2015).

AC C

636

51

ACCEPTED MANUSCRIPT

641

The most obvious characteristic of lacustrine shales is that they are subject to more variable paleoenvironments than marine shales. The formation environments of lacustrine and marine shales are

643

very different in terms of the water depth, energy, provenance, sedimentary type and biological effects

644

(Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002). These features lead to

645

differences in their basic characteristics including mineralogy, organic matter (OM), and reservoir type

646

(Table 1). Minerals are the material basis of shales, and they affect not only OM generation but also

647

reservoir development (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et

648

al., 2015; Yang et al., 2015; Jiang et al., 2016a). A major breakthrough in marine shale oil and gas

649

exploration was based on the recognition that shale acts not only as a source rock but also as a reservoir.

650

Unlike studies of conventional reservoirs, research on unconventional reservoirs has primarily focused

651

on nanoscale reservoirs. Previous studies of marine and lacustrine shales suggest that shale gas is mainly

652

controlled by a nanoscale pore system that is closely related to the mineral characteristics (Zou et al.,

653

2011; Morad et al., 2010; He et al., 2012; Luo et al., 2013; Dang et al., 2015; Chalmers et al., 2012;

654

Curtis et al., 2012; Jiang et al., 2015; Li et al., 2016; Jiang et al., 2017; Zhang et al., 2017). OM is the

655

basic material that generates hydrocarbons, and its type and maturity dominate the formation and

656

enrichment of shale oil and gas. The preservation and enrichment of OM in sediments is controlled by

657

many factors, including the primary productivity, paleo-redox conditions, nutrient availability, clastic

658

influx and minerals (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002; Chen et al.,

659

2015; Yang et al., 2015; Jiang et al., 2016a; Xu et al., 2017). In summary, the paleoenvironment,

660

minerals, reservoir, and OM, which are the basic characteristics of shales, have been the focuses of shale

AC C

EP

TE D

M AN U

SC

RI PT

642

52

ACCEPTED MANUSCRIPT

oil and gas studies, and these characteristics are often interrelated. Additionally, studies of the

662

characteristics and mechanisms have mainly focused on marine shales; fewer studies have examined

663

lacustrine shales.The lacustrine shale of the Jurassic Da’anzhai Member from the central Sichuan Basin

664

of China has been interpreted as fractured reservoirs, low permeability fractured reservoirs and tight

665

reservoirs. However, most studies have only focused on the conventional carbonate reservoirs, and the

666

shales have only been considered as source rocks for the carbonate reservoirs (Lu et al., 2014; Chen et

667

al., 2015). Encouragingly, the Da’anzhai shale of the Yuanba Block in the northwestern Sichuan Basin,

668

which has similar geological conditions to those of the central Sichuan Basin, hosts high-yield wells (He

669

et al. 2012; Zhu, et al. 2016). In addition, some lacustrine shales have yielded oil or gas, such as the

670

Triassic Chang 7 and Chang 9 Member in the Ordos Basin, the Paleogene strata (Es3 and Es4) in the

671

Bohai Bay Basin and the Cretaceous strata in the Songliao Basin (Qing-1) (Zou et al., 2011; Luo et al.,

672

2013; Dang et al., 2015; Li et al., 2015; Yang et al., 2015; Jiang et al., 2016b). These successful

673

developments strongly indicate that the Da’anzhai shale in the central Sichuan Basin is worth studying.

TE D

M AN U

SC

RI PT

661

In this study, we assessed the Da’anzhai lacustrine shale’s characteristics. We performed

675

paleoenvironmental evaluation using TE, mineralogical evaluation using XRD analysis, reservoir

676

evaluation using FE-SEMand LTNA analyses, and OM evaluation using analyses of the TOC,

677

chloroform bitumen ‘A’, Rock-Eval, carbon isotopes, kerogen macerals and organic elements. We also

678

compared the Da’anzhai shale to other lacustrine shales which have been successfully developed and

679

assessed the similarities and differences of lacustrine shales. In addition, we summarized the

AC C

EP

674

53

ACCEPTED MANUSCRIPT

relationships between these characteristics and mechanisms affecting them. Finally, we judged the

681

development potential of the Da’anzhai shale.

AC C

EP

TE D

M AN U

SC

RI PT

680

54

ACCEPTED MANUSCRIPT

Table 1 Comparisons between basic parameters in the Da’anzhai shale and other successfully developed shales lacustrine shale

basin

Sichuan Basin, China

basin type

layer

foreland basin

Jurassic shale with

Chang 9

Es3

Es4

Qing-1

Barnett

Ordos Basin,

Ordos Basin,

Bohai Bay

Bohai Bay

China

China

Basin, China

Basin, China

polycyclic

polycyclic

superimposed

superimposed

downfaulted

downfaulted

basin

basin

basin

basin

basin

Triassic

Triassic

Paleogene

Paleogene

Cretaceous

Carboniferous

shale with

siliceous and

argillaceous

carbonaceous

limestone

shale, limestone,

interlayers

dolomite

RI PT

Central

Chang 7

Songliao

Fort Worth Basin,

Basin, China

Texas, US

downfaulted

SC

Da'anzhai

marine shale

M AN U

parameters

foreland basin

Ohio Appalachian Basin, Kentucky, US piedmont depression Devonian

shale with

shale with

shale with

shale with

siltstone

siltstone

siltstone

siltstone

interlayers

interlayers

interlayers

interlayers

1900~3200

1500~2500

1600~2600

1500~4200

1500~5200

1000~2500

1980~2590

510~1800

10~70

30~70

10~15

100~500

250~350

70~150

15~60

9~31

(20~80) /54

(37~53) /43

(29~56) /45

(25~80) /36

(30~80) /36

(30~80) /36

35~80

45~60

porosity (%)

0.20~5.92

1~12

permeability

0.002~0.9

TOC (%) Ro (%)

lithological

shells and

combination

carbonaceous interlayers

depth (m)

thickness (%)

minerals (%)

genetic type

EP

brittle

TE D

effective

interbed of black shale and gray siliceous shale

10

2~8

1.3~9.3

6~12

4~5

4~7

0.01~0.3

(0.3~5.0) /2.97

<0.50

<0.50

<0.15

0.01

<0.1

(0.1~2.2) /1.0

(0.3~36) /8.3

(0.3~11.3) /3.1

(0.5~13.8) /3.5

(0.8~16.7) /3.2

(0.4~4.5) /2.2

2.0~7.0

0.5~4.7

0.9~1.5

0.8~1.2

0.9~1.3

0.4~1.2

0.4~2.0

0.5~1.5

1.1~2.2

0.3~1.3

pyrolysis gas

pyrolysis gas

pyrolysis gas

pyrolysis gas

Ⅰ-Ⅱ1

Ⅰ-Ⅱ





AC C

682

pyrolysis gas

pyrolysis gas

pyrolysis gas



Ⅰ-Ⅱ1

Ⅰ-Ⅱ1

biogenic and pyrolysis gas

kerogen type

Ⅰ-Ⅱ1

55

ACCEPTED MANUSCRIPT

al., 2013; Lu et al., 2014; Chen et al., 2015;

/2.00

0.31~1.24

8.50~9.91

Zou et al.,

Zou et al., 2011; Luo et al., 2013; Yang et al., 2015; Jiang et al., 2016a.

/2.92

2011; Luo et

Zou et al., 2011;

Zou et al.,

al., 2013; Shi et

Luo et al., 2013;

2011; Luo et

al., 2013; Yang

Dang et al.,

al., 2013; Li,

et al., 2015;

2015; Jiang et

2014; Jiang et

Jiang et al.,

al., 2016a

al., 2016a

2016a.

Zou et al.,

56

1.69~2.83

Schmoker, 2002;

Schmoker, 2002;

Bowker et al.,

Bowker et al.,

RI PT

references

(0.60~3.70)

2011; Luo et

2007; Jarvie et al.,

2007; Jarvie et al.,

al., 2013; Jiang

2007; Martini et

2007; Martini et

et al., 2016a

al., 2008; Lu et al.,

al., 2008; Lu et al.,

2014

2014.

SC

2011; Luo et

(1.68~4.25)

M AN U

Zou et al.,

0.92~2.91

TE D

(m /t)

0.87~1.98

EP

3

AC C

gas content

ACCEPTED MANUSCRIPT

2. Geological setting

686 687

EP

685

Fig. 1. Study location, geological map (Chen et al., 2015) and simplified stratigraphic column of the Jurassic Da’anzhai

Member in the central Sichuan Basin.

AC C

684

TE D

M AN U

SC

RI PT

683

The study area is located in the central Sichuan Basin, which includes the Nanchong and Suining

688

areas. It is bounded by the Huafu Mountains to the east and the Longquanshan basement faults to the

689

west, and it covers an area of 6×104 km2 (Fig. 1) (Lu et al., 2014; Chen et al., 2015). A lacustrine

690

environment dominated the Sichuan Basin during the Jurassic Period, and the Da’anzhai Member

691

represents the largest and deepest stage of this environment. The sedimentary facies of the Da’anzhai 57

ACCEPTED MANUSCRIPT

Member can be divided into shore-shallow lacustrine, shallow lacustrine and semi-deep lacustrine

693

subfacies, which include shell beach, calcareous flat, mud flat, and semi-deep lake microfacies. These

694

facies are distributed in concentric belts around the lake basin’s center in the Yilong-Pingchang areas

695

(Fig. 1) (Lu et al., 2014; Chen et al., 2015). From top to bottom, the Da’anzhai Member in the central

696

Sichuan Basin can be divided into main three sub-members: Da1, Da13 and Da3 (Fig. 1). The thickest

697

Da13 sub-member (30-50 m) is a set of widely distributed dark mudstones with lamellated fractures and

698

rich OM; this sub-member is considered the most important source rock for hydrocarbons. The Da3

699

sub-member deposits (5-14 m), which are initial deposition products, are primarily (mud) shell

700

limestones with dark gray mudstone interlayers. The Da1 sub-member (30-45 m thick) primarily

701

consists of shales and (mud) shell limestones; this sub-member contains the sedimentary products of the

702

period of basin uplift (Lu et al., 2014; Chen et al., 2015).

TE D

M AN U

SC

RI PT

692

3. Samples and experimental methods

704

3.1. Samples

705

Thirty-nine representative core samples were obtained from the Da’anzhai Member shale in the

EP

703

central Sichuan Basin based on drilling, core observations, stratigraphic data, and log responses (Fig. 1

707

and Table 3).

AC C

706

708

3.2. Methods

709

All 39 shale samples were analyzed to determine their bulk mineralogy via XRD using a X'Pert Pro

710

MPD device with a working voltage of 40 kV and a current of 40 mA in the Key Laboratory of Orogenic

711

Belts and Crustal Evolution at Peking University. To investigate the clay mineral fraction in detail,

58

ACCEPTED MANUSCRIPT

712

textured mounts were measured four times at a goniometer rate of 0.5θ/minute with a registration range

713

from 2 to 42°2θ at the Beijing Research Institute of Uranium Geology. The micromorphology and oiliness of the samples were determined at the Key Laboratory of

RI PT

714

Marine Reservoir Evolution and Hydrocarbon Accumulation of the China University of Geosciences

716

(Beijing) using an Axio Scope A1 pol Fluorescence-Polarization Microscope. An INCA Synergy system

717

at the Key Laboratory of Orogenic Belts and Crustal Evolution at Peking University was used to analyze

718

the microstructure of each sample via FE-SEM and the mineral composition of each sample via EDS.

719

The samples were subjected to argon-ion polishing to obtain higher resolution images.

M AN U

720

SC

715

LTNA analysis is capable of characterizing nanoscale pores in porous media, and adsorption tests were performed with dried samples (60-80 mesh) on an ASAP2020M device at the College of

722

Chemistry and Chemical Engineering of the China University of Petroleum (Beijing). N2 adsorption

723

isotherms were obtained at -196.15°C within a relative pressure range of 0.005 to 0.998. The BET

724

method for the equivalent specific surface areas and the BJH method for the total pore volumes were

725

applied in the calculations (Brunauer et al., 1938; Barrett et al., 1951).

EP

TE D

721

The OM richness was characterized by measuring the TOC values of 25 samples using a GHM-02

727

analyzer; 10% hydrochloric acid was used to remove carbonates, and the test sensitivity was 10-13 mg/g.

728

Chloroform bitumen ‘A’ was extracted from the crushed shale powders (<0.09 mm), and the weight was

729

calculated at the China University of Petroleum (Beijing). A Rock-Eval 6 instrument was used to

730

perform a pyrolysis analysis, which yielded the temperature of the Tmax, S0, S1, and S2. The OM maturity

731

was determined based on the Ro value, which was obtained using a reflected light microscope under oil

AC C

726

59

ACCEPTED MANUSCRIPT

immersion and calculated as the mean of 10-30 measurements taken for a single sample. Additionally,

733

measurements were taken to determine the isolated kerogen content and identify the organic kerogen

734

macerals, carbon isotope content and organic element content according to the (GB/T) 19144-2010,

735

(SY/T) 5125-2014, (GB/T) 18340.2-2010, (GB/T) 19143-2003 standards, respectively. The kerogen

736

measurements and the rock pyrolysis analysis were performed at the Keyuan Engineering Technology

737

Testing Center.

SC

The TE concentrations of 17 samples were determined using ICP-MS, which was performed at the

M AN U

738

RI PT

732

739

Beijing Research Institute of Uranium Geology. The detailed sample processing procedure as well as the

740

analytical precision and accuracy followed the methods described for the (GB/T) 14506.28-2010

741

standard.

4. Results and discussions

743

4.1. Paleoenvironmental properties and differences

744

Lacustrine shales are different from marine shales mainly due to their formation paleoenvironments.

EP

TE D

742

Lakes are complex and dynamic systems, and the development of lacustrine shales is controlled more

746

significantly by the sedimentary environment (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002;

AC C

745

60

ACCEPTED MANUSCRIPT

Table 2 Trace elements of the Da’anzhai lacustrine shale. G10-2

G4-3

G6-2

J45-1

J45-2

J53-1

J61-1

J61-2

J61-3

J61-4

Li-1

Li-2

P1-1

S2-1

S2-3

S2-4

X44-1

Th(ppm)

11.70

10.40

5.43

4.32

12.00

10.60

1.72

4.13

7.74

10.20

11.30

9.50

11.60

2.31

10.50

6.21

13.70

Ni/Co

2.68

2.98

3.26

4.10

2.87

2.72

5.49

3.87

3.50

2.51

2.65

2.97

3.07

3.22

2.32

3.37

2.78

RI PT

code

1.76

1.58

1.44

1.50

1.64

1.64

2.47

1.56

1.33

1.69

1.60

1.50

1.61

1.77

1.30

1.41

1.45

0.78

0.80

0.63

0.63

0.78

0.78

0.67

0.68

0.70

0.77

0.78

0.74

0.78

0.61

0.68

0.68

0.77

V/Sc

9.78

8.70

6.88

7.73

9.36

8.91

13.81

8.41

7.94

9.03

7.93

7.21

9.22

9.80

8.54

7.14

7.81

Th/U

2.77

4.24

4.45

2.79

4.07

2.95

1.28

3.53

4.63

3.81

4.25

4.63

4.20

2.74

5.25

4.60

4.72

Sr/Cu

3.75

4.60

42.98

77.10

3.50

3.62

58.70

35.54

53.09

8.46

7.84

19.01

3.35

80.54

7.47

24.69

6.05

Sr/Ba

0.29

0.26

3.66

3.08

0.23

0.24

4.77

2.78

3.45

0.52

0.46

1.28

0.18

3.77

0.44

1.84

0.36

118.7

115.3

109.8

131.7

110.3

187.5

141.2

128.9

194.3

112.9

130.6

112.6

3

6

4

1

7

7

6

9

0

9

7

2

7.79

9.03

6.39

6.65

9.34

9.58

8.46

8.67

8.49

8.04

9.44

4.19

9.93

9.86

8.88

7.88

10.03

(La/Yb)N

6.65

7.75

6.57

8.37

6.82

8.04

9.85

9.82

11.14

8.00

7.50

4.35

7.49

10.88

8.62

7.67

9.05

(La/Sm)N

4.66

4.74

3.04

3.32

5.58

5.42

3.19

3.06

2.37

3.88

5.71

2.59

5.90

3.07

4.20

3.90

5.01

(Gd/Yb)N

0.99

1.04

1.52

1.97

0.77

0.97

2.04

2.01

2.92

1.45

0.84

1.41

0.79

2.10

1.35

1.38

1.17

(Ce/Yb)N

4.48

5.38

4.98

5.83

4.54

5.42

7.21

7.33

9.24

5.58

4.97

3.22

5.00

8.49

6.10

5.59

6.15

ΣLREE/ ΣHREE

88.73

M AN U

)

130.97

38.15

TE D

ΣREE(ppm

SC

V/Cr V/(V+Ni)

60.76

153.93

0.92

0.94

0.97

0.95

0.92

0.93

0.94

0.93

0.99

0.92

0.92

0.92

0.93

0.96

0.95

0.97

0.92

0.19

0.18

0.77

1.12

0.17

0.16

0.63

0.56

0.88

0.36

0.27

0.61

0.15

0.72

0.12

0.35

0.23

EP

δCe CSr(‰)

AC C

XN, where N refers to the chondrite-normalized value (Gromet et al., 1984); δCe=CeN/(LaN×PrN)1/2 (Wright et al.,1984; Murray et al.,1990).

61

ACCEPTED MANUSCRIPT

Schmoker, 2002). TEs can be used as effective quantitative indicators of the paleoenvironment by

2

comparing their type and abundance with standard values (Wright et al., 1984; Murray et al., 1990;

3

Manning, 1994; Ross et al., 1995; Kimura and Watanabe, 2001; Jarvis et al., 2001; Rimmer, 2004; Jones

4

and Roy et al., 2007).

5 6

TE D

M AN U

SC

RI PT

1

Fig. 2. Identification of paleoenvironments using TEs (a-c) and clay (d).

4.1.1. Paleo-redox conditions

8

Due to their responses to changing redox conditions, the ratios of redox-sensitive TEs, such as

AC C

9

EP

7

Ni/Co, V/Cr, V/Sc, Th/U, and δCe, are usually used to evaluate redox conditions. According to Jones

10

and Manning (1994), V/Cr ratios <2 indicate oxic conditions, V/Cr ratios from 2 to 4.25 indicate dysoxic

11

conditions, and V/Cr ratios >4.25 indicate suboxic to anoxic conditions. Kimura and Watanabe (2001)

12

found that the V/Sc ratio is positively related to oxic conditions. Jones and Manning (1994) showed that

13

Ni and Co are enriched in pyrite, and high Ni/Co ratios indicate anoxic conditions, Ni/Co ratios <5 62

ACCEPTED MANUSCRIPT

indicate oxic conditions, Ni/Co ratios between 5 and 7 indicate dysoxic conditions, and Ni/Co ratios >7

15

indicate suboxic to anoxic conditions. Similarly, Th/U ratios <1.5 represent anoxic conditions, Th/U

16

ratios between 1.5 and 3 represent dysoxic to suboxic conditions, and Th/U ratios >3 represent oxic

17

conditions (Jones and Manning, 1994). In addition, δCe values >1 (positive anomaly) and <0.95

18

(negative anomaly) suggest anoxic and oxic to suboxic conditions, respectively (Wright et al., 1984;

19

Murray et al., 1990).

SC

For the Da’anzhai shale in the central Sichuan Basin (Table 2), the ranges of the Ni/Co, V/Cr, V/Sc,

M AN U

20

RI PT

14

and Th/U ratios are relatively wide, from 2.32-5.49, 1.30-2.47, 6.88-13.81, and 1.28-5.25, respectively.

22

The values of δCe vary from 0.92 to 0.99 (average 0.94), which indicate a weakly negative anomaly

23

(Wright et al., 1984; Murray et al., 1990).Collectively, these proxies reveal that the shale was deposited

24

under variable paleo-redox conditions that were mainly oxic to suboxic (Fig. 2).

TE D

21

4.1.2. Paleo-salinity

26

Sr/Ba values vary directly with the distance to a sea or lake coast; Sr/Ba ratios <1 suggest

EP

25

freshwater conditions, and Sr/Ba ratios >1 suggest salt water conditions (Jarvis et al., 2001; Kimura and

28

Watanabe, 2001). In addition, low Sr concentrations (0.1-0.3‰) indicate a freshwater setting, and high

29

values (0.8-1‰) represent a salt water setting (Jarvis et al., 2001; Kimura and Watanabe, 2001). Th

30

easily dissolves in acidic water and hydrolyzes to oxide or hydroxide sediments under alkaline

31

conditions; therefore, high Th values indicate a salt water environment (Jones and Manning, 1994).

32 33

AC C

27

In the Da’anzhai shale samples (Table 2), the Sr/Ba ratios vary widely (0.18-4.77) with an average of 1.62. The Th content also varies widely (1.72-13.70 ppm) and averages 8.43 ppm, and the Sr 63

ACCEPTED MANUSCRIPT

34

concentrations vary from 0.12 to 1.12‰ (average 0.44‰). Generally, these proxies indicate variable

35

paleo-salinity conditions including both fresh and salt water settings (Fig. 2).

RI PT

4.1.3. Paleo-weathering

M AN U

SC

36

37

39

Fig. 3. Chondrite normalized REE distribution curves for REE minerals. Chondrite values are from Gromet (1984).

When rocks are subjected to strong weathering, heavy rare earth elements (HREEs) dissolve and

TE D

38

migrate much more easily than light rare earth elements (LREEs). Therefore, high ΣLREE/ΣHREE

41

values indicate strong paleo-weathering conditions (Kimoto et al., 2006; Stevens and Quinton, 2008). In

42

addition, the (La/Yb)N, (La/Sm)N, (Gd/Yb)N and (Ce/Yb)N ratios, which are positively related to LREE

43

enrichment, are also positively correlated with paleo-weathering (Table 2) (Wright et al., 1984; Murray

44

et al., 1990; Ross et al., 1995; Roy et al., 2007).

AC C

45

EP

40

The ΣLREE/ΣHREE rates are between 4.19 and 10.03 (average 8.39), which indicate LREE

46

enrichment. The (La/Yb)N ratio is represented by the slope of the chondrite-normalized REE distribution

47

curve (Fig. 3), which indicates the degree of the graph’s inclination (Ross et al., 1995; Roy et al., 2007).

48

The (La/Yb)N values vary from 4.35 to 11.14 (average 8.15), which indicate that the curves are tilted to 64

ACCEPTED MANUSCRIPT

the right and that the samples are rich in LREE-acidic rocks (Fig. 3). The (La/Sm)N ratio, which

50

represents the degree of LREE fractionation, is positively related to the ΣLREE value (Ross et al., 1995;

51

Roy et al., 2007). The (La/Sm)N ratios range from 2.37 to 5.90 (average 4.10), which indicates a P(E)

52

model with abundant LREEs. The (Gd/Yb)N ratio, which indicates the degree of HREE fractionation, is

53

negatively correlated with the ΣHREE value (Ross et al., 1995; Roy et al., 2007). The (Gd/Yb)N values

54

vary from 0.77 to 2.92 (average 1.46), and the (Ce/Yb)N ratios vary from 3.22 to 9.24 (average of 5.85),

55

which indicate that the Da’anzhai shale is depleted in HREEs (Fig. 3). In summary, the REE results

56

(Table 2) indicate that the shale was deposited under moderate weathering conditions with periods of

57

strong weathering.

M AN U

SC

RI PT

49

4.1.4. Paleoclimate

59

∑REE values are closely related to paleoclimate conditions; that is, ∑REE values are higher in

TE D

58

warm and humid climates and lower in cold and dry climates (Wright et al., 1984; Murray et al., 1990;

61

Ross et al., 1995; Roy et al., 2007). The ΣREE values of the Da’anzhai shale vary widely from 38.15 to

62

194.3 µg/g (average 121.59 µg/g) (Table 2). The Sr/Cu ratio is another useful indicator of lacustrine

63

paleoclimate, with high values reflecting humid climates (Reheis, 1990; Chen et al., 2009).The Sr/Cu

64

ratios in the Da’anzhai shale also vary widely from 3.35 to 80.54%. In addition, illite usually forms in

65

dry conditions by the weathering depotassication of feldspar, mica and other aluminosilicate minerals.

66

Conversely, kaolinite forms in humid climates by the eluviation of feldspar, mica, and pyroxene.

67

Therefore, high illite/clay ratios indicate dry, saline waters with high K+ content, and high kaolinite/clay

68

ratios suggest humid climates (Vanderaveroet, 2000; Gingele et al., 2001). The illite/clay ratios of the

AC C

EP

60

65

ACCEPTED MANUSCRIPT

Da’anzhai shale are relatively high (between 26 and 86%, average 52%). However, the kaolinite/clay

70

ratios are relatively low (between 3 and 19%, average 11%) (Fig. 2d). These results indicate a variable

71

climate with both dry and humid periods as well as a drying trend from humid climates.

72

RI PT

69

In general, the TE contents and types t in lacustrine shales are more complex and varied than those in marine shales, which indicates that lacustrine shales tend to experience more variable

74

paleoenvironmental conditions (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002).

75

The lacustrine Da’anzhai shale was formed under relatively dry climatic conditions that included both

76

arid and humid stages. In addition, the shale experienced a moderate degree of paleo-weathering with

77

stages of strong weathering, variable oxic to suboxic paleo-redox conditions that were weakly reducing,

78

and a predominantly freshwater setting interspersed with periods of saltwater.

M AN U

4.2. Mineralogy properties and differences

TE D

79

SC

73

As shown in Table 3, the Da’anzhai shale primarily contains clay (42.4%), quartz (27.7%) and calcite

81

(21.1%). The minerals in these lacustrine shales vary widely. Clay comprises the majority of the shale

82

and ranges from 16.3 to 68.4%, while the quartz content ranges from 5.6 to 53.4%. Carbonate minerals,

83

which mainly include calcite, dolomite, and aragonite, range from 3.0 to 46.6% with an average of

84

27.4%. The Da’anzhai shale contains more brittle minerals than the other lacustrine shales and marine

85

shales (Table 1 and Fig. 4). For example, the brittle mineral content of the Da’anzhai shale ranges from

86

24.2 to 80.3%, with an average of 53.7%, indicating medium-high brittleness that is favorable for natural

87

fractures and hydraulic fracturing. Similar to other lacustrine shales, lacustrine shales have more

88

siliceous or calcareous interlayers than marine shales, which can cause the wider ranges of carbonate and

AC C

EP

80

66

ACCEPTED MANUSCRIPT

quartz minerals and improve the brittleness index and fracturing (Figs1 and 5I).Ternary diagrams were

90

constructed to compare the mineral compositions of hot shales that have been successfully exploited

91

including lacustrine shales (Fig. 4a) and marine shales (Fig. 4b). Generally, the comparison shows that

92

the mineral composition of the lacustrine shales is more variable and complex than in marine shales (Fig.

93

4). For example, the Da’anzhai shale has relatively wide ranges of carbonate (1-70%), clay (15-70%)

94

and quartz (5-50%) (Fig. 4a). In addition, the Da’anzhai shale has the wider mineral ranges than other

95

lacustrine shales, especially for the range of carbonate minerals.

M AN U

SC

RI PT

89

EP

97

TE D

96

Fig. 4. Ternary diagrams of mineral compositions. (a) Lacustrine shales: 1=Da’anzhai shale, 2=Jurassic shale in Sichuan

99

Basin (Jiang et al., 2017), 3=Triassic shale in Sichuan Basin(Jiang et al., 2017), 4=ES3 shale (Wang et al., 2015a). (b) Marine

AC C

98

100

shales: 5 Woodford shale (Jarvie et al. 2007 and Han et al. 2013), 6=Barnett shale (Jarvie et al. 2007 and Han et al. 2013),

101

7= Wufeng-Longmaxi shale (Wu et al., 2014), 8=Ohio shale (Jarvie et al. 2007 and Han et al. 2013).

67

ACCEPTED MANUSCRIPT

Table 3 Minerals and TOC results for the Da’anzhai lacustrine shale. feldspar

calcite

dolomite

aragonite

siderite

pyrite

hematite

(%)

(%)

(%)

(%)

(%)

(%)

(%)

(%)

G10-1

2650.8

Da1

44.2

G10-2

2668.4

Da13

22.1

1.1

0.9

2700.7

Da13

37.5

2.0

2.9

2646.6

1

28.1

4.2

22.9

1

G10-3 G10-4

Da

3.0 1.2

6.3 2.6

1.2

17.0

2385.3

Da

28.2

G4-2

2395.8

Da1

35.2

2.1

6.4

12.5

G4-3

2405.7

Da13

25.8

2.5

4.6

0.7

2537.1

Da13

21.4

2.6

2580.1

3

G6-2

Da

7.4

0.82

24.2

1.33

55.0

2.9

40.4

0.81

26.6

41.1

69.2

54.9

4.5

38.4

1.6

64.8

3.5

60.9

13.4

41.6

0.40

18.9

54.1

0.98

5.3

31.1

1.43

4.2

7.4

28.8

1.18

0.9

36.3

0.4

0.7

20.5

1.1

64.2

76.7

0.4

1.1

64.0

1.1

3.9

27.3

3.1

55.5

5.4

0.2

22.6

1.0

59.3

76.2

1.0

0.8

26.7

1.1

64.8

70.4

0.0

27.7

0.9

57.0

71.4

1.3

4.2

33.1

2.3

24.3

59.1

1.49

2.1

20.0

69.4

77.9

0.30

46.2

26.9

52.0

9.8

40.2

1.17

38.5

35.2

59.4

1.07

53.9

3.1

43.8

0.69

6.4

67.8

1.6

24.9

1.59

4.1

0.9

60.2

4.1

38.9

1.51

3.6

51.2

14.9

45.2

1.55

1.5

16.3

57.1

80.3

0.61

3.9

2827.3

Da13

31.9

4.1

2662.6

1

16.9

14.4

1

Da Da

5.6

55.3

1

44.9 9.5

J61-3

2667.6

Da

14.4

14.3

J61-4

2698.1

Da13

34.8

17.5

2706.8

Da13

8.5

69.4

2665.9

1

25.1

1.8

15.6

1

30.4

2.6

9.8

1.0

28.9

6.3

2.3

3.1

0.0

0.9

0.6

Da

1

1726.1

Da

24.2

Li-3

1734.3

Da13

40.7

3215.9

Da13

23.3

3221.2

Da13

34.8

3239.2

Da13

30.3

14.9

3258.6

3

23.2

51.8

Da

42.7 6.8

EP

Da

Li-2

P1-4

47.2

2.1

75.8

27.9

P1-3

3.0

68.4

67.1

6.1

P1-2

52.8

1.0

0.6

P1-1

(%)

20.1

12.5

1713.8

minerals

0.8

23.4

Li-1

(%)

1.1

Da

J61-6

(%)

(%)

19.5

Da13

J61-5

TOC

47.6

2656.0

2664.4

brittleness

1.2

2632.4

J61-2

carbonate

8.7

J45-2

J61-1

others

clay

1

J45-1

J53-1

3.5

TE D

G6-1

13.4

M AN U

G4-1

total

RI PT

(m)

sub-member

quartz

SC

code

depth

AC C

102

1.0

10.3 3.2

53.3 1.1

1.0

5.3

68

0.7

1.9

31.9

2.18

ACCEPTED MANUSCRIPT

2836.1

Da1

19.8

2849.9

1

Da

40.3

2867.1

Da13

17.6

47.7

1.3

S2-4

2889.0

3

Da

18.9

48.9

1.5

S2-5

2891.6

Da3

53.4

X20-1

1844.4

Da13

49.4

1952.2

1

Da

24.9

2089.4

Da13

35.9

X44-2

2093.0

Da13

27.9

3.0

4.0

X44-3

2116.7

Da13

23.8

3.0

43.2

LQ-1

3548.1

Da13

44.1

2.9

18.7

1780.7

Da13

40.5

2.6

2.5

2065.7

1

28.2

2.2

18.5

7.8

43.3

1

23.9

1.9

35.2

11.4

26.6

X8-1 W8-1 W8-2

2069.5

Da Da

29.6

0.9

42.5

50.4

70.2

0.11

10.7

51.0

1.18

33.4

49.0

66.6

0.33

30.7

50.4

69.3

0.18

2.0

8.1 1.4

2.3

25.6

3.1

8.5

RI PT

46.6

3.8

3.7

SC

X44-1

0.2

9.5

2.5

M AN U

X28-1

10.7

TE D

S2-3

3.6

2.2

EP

S2-2

48.2

AC C

S2-1

69

53.4

38.7

1.5

8.1

57.5

44.3

0.7

28.7

53.6

8.5

44.4

0.53 1.15

51.8

60.8

4.0

31.9

20.5

52.7

76.5

34.3

18.7

62.8

2.5

43.0

26.3

54.5

46.6

70.5

50.8

0.6

1.1

1.0

0.76

0.84

ACCEPTED MANUSCRIPT

4.3. Characterization and controlling factors of lacustrine shale reservoirs

104

4.3.1. Microscopic characteristics

105

AC C

EP

TE D

M AN U

SC

RI PT

103

106

Fig.5. Typical reservoirs images of lacustrine shale. A: Shells were stained using alizarin red to highlight the carbonate

107

minerals; the shells are well developed and preserved in terrestrial inputs and matrices. B: Greenish-yellow fluorescence.

108

Matrices with large amounts of terrestrial inputs between the shell grains have strong fluorescence intensity, indicating

70

ACCEPTED MANUSCRIPT

effective reservoirs with good oiliness in the matrices. C: Lacustrine shale core sample with well-developed calcareous

110

interlayers, indicating a changeable sedimentary environment. D: Quartz particles with well-developed interparticle pores

111

filled with clay minerals; quartz and clay minerals have well-developed intraparticle pores. E: Quartz particles have

112

well-developed intraparticle pores and interparticle pores; calcite minerals fill the storage spaces between quartz particles

113

(calcite cementation). F: Intraparticle pores located along the cleavage planes of clay particles; interparticle pores between

114

clay and carbonate minerals. G: OM with well-developed pores is accompanied by clay minerals. H: Calcite minerals fill the

115

spaces between particles (calcite cementation); quartz has well-developed intraparticle pores. I: Brittle minerals with

116

well-developed micro-cracks that are mainly filled by clay minerals. J: Model of different types of pores. K: Model of

117

carbonate mineral (e.g., calcite) diagenesis that results in reduced storage spaces.

SC

M AN U

118

RI PT

109

According to the International Union of Pure and Applied Chemistry (IUPAC) classification scheme, the pore types in unconventional reservoirs are divided by pore size into micropores, mesopores,

120

and macropores (<2 nm, 2-50 nm, and >50 nm, respectively) (IUPAC, 1994; Thommes et al., 2015).

121

Thin section observations of the shale indicate that the matrices fluoresce intensely (Fig. 5B). Nanoscale

122

reservoirs were examined using FE-SEM with argon-ion milling technology, and the minerals were

123

identified using EDS. Similar to other lacustrine shales, the pores and fractures in the Da’anzhai shale

124

are well developed and fall into four categories.

EP

AC C

125

TE D

119

(1) Interparticle pores: these pores represent the largest number of pores, and they have good

126

connectivity and form effective pore networks. The diameters generally range between 5 nm and 90 µm.

127

Many types of interparticle pores were observed, including pores between grains (Figs. 5D,E,H,J),

128

crystals and clay platelets (Figs. 5D,F,G,J), and pores at the edges of rigid grains (Figs. 5E,J).

71

ACCEPTED MANUSCRIPT

129

(2) Intraparticle pores: these pores primarily formed through diagenesis, and some are primary. The diameters mainly range from 2 nm to 10 µm. The primary sub-types in the Da’anzhai shale are

131

intraplatelet pores within clay aggregates (Figs. 5D,F,J), intercrystalline pores within pyrite framboids,

132

pores within fossil bodies and moldic pores after fossils (Figs. 5D,J).

133

RI PT

130

(3) Fracture pores: these pores are not controlled by individual particles and can be regarded as a type of pore that formed by rock deformation because of tectonic movements, sedimentation, diagenesis,

135

OM hydrocarbon generation and other geological effects (Figs. 5I,J).

M AN U

136

SC

134

(4) OM pores: these pores are found within OM (Fig. 5L) and are usually irregular, bubble-like, and oval. They generally range in diameter from 2 to 500 nm. Although most OM pores are incorrectly

138

regarded as isolated pores in two-dimensional planes, they are interconnected in three-dimensional space

139

and facilitate the development of reservoirs.

140

TE D

137

Similar to other lacustrine and marine shale reservoirs, the Da’anzhai shale commonly contains well-developed pores and fractures, particularly intraparticle and interparticle pores (Curtis et al., 2002;

142

Chalmers et al., 2012; Milliken et al., 2013; Nie et al., 2015; Wang et al., 2015a; Zhou and Kang, 2016;

143

Jiang et al., 2017; Zhang et al., 2017). Especially, due to their rich terrestrial clastic mineral content, the

144

interparticle pores of lacustrine shales are better developed than those of marine shales. Most of the

145

intraparticle pores are nanoscale pores (<200 nm) (Figs. 5D,E,F) and the OM also contains nanopores

146

(Fig.5G). Microcracks can develop due to the high brittle mineral content for the Da’anzhai shale.

147

However, these cracks are often filled by calcite and clay minerals due to the high clay content (42.4%)

148

(Figs. 5I,J).

AC C

EP

141

72

ACCEPTED MANUSCRIPT

4.3.2. Nanoscale reservoirs characterized by LTNA

150

The LTNA test is the most widely used method of characterizing several parameters of nanoscale

151

pores in porous media (micropores and mesopores), including VBJH, SBET, DA, ad/desorption isotherms

152

and PSD (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015; Wang et al., 2015c; Zhou et al., 2016;

153

Zhang et al., 2017).

154

Table 4 Comparisons of primary pore structure parameters measured by LTNA between the Da’anzhai shale and other

155

lacustrine shales

4

5

6

7

156 157

SC

Sichuan

Jurassic Da'anzhai

1.75-10.69

Basin

Member

(6.92)

M AN U

3

SBET (m2/g)

VBJH 3

(cm /100g)

D (nm)

1.19-4.10

7.21-24.71

(2.97)

(14.11)

0.02-0.90

6.67-18.40

(0.49)

(11.23)

0.22-2.97

4.91-13.93

Upper Triassic

0.25-4.39

Chang 7 Member

(2.62)

Sichuan

Xujiahe

0.62-14.70

Basin

Formation

(7.76)

(1.41)

(6.97)

1.46-5.24

0.72-1.40

9.40-25.00

Ordos Basin

Ordos Basin

TE D

2

strata

Upper Triassic

Chang 7 Member

(2.94)

(1.01)

(16.40)

Bohai Bay

Shahejie

0.40-12.45

0.08-1.60

4.58-41.36

Basin

Formation

(3.62)

(0.62)

(16.84)

0.39-34.06

0.13-4.18

3.16-16.48

(10.92)

(1.59)

(9.56)

Upper Triassic

1.10-1.90

0.69-1.09

15.4-23.4

Chang 7 Member

(1.57)

(0.83)

(19.27)

Songliao Basin

Upper Cretaceous

EP

1

location

Qingshankou Formation

AC C

NO.

RI PT

149

Ordos Basin

73

hysteresis

loops types

references

H3

H3 H3 H3 H3,H4

H2,H3

H3,H4

Jiang et al., 2016b Peng et al., 2016 Fu et al., 2015 Zhang et al., 2016 Wang et al., 2015b Yang et al., 2017

RI PT

ACCEPTED MANUSCRIPT

158

160

Fig. 6. PSD curves obtained from the low temperature N2 adsorption isotherms with BJH pore sizes.

SC

159

The VBJH value of the Da’anzhai shale (2.97 cm3/100 g) is slightly higher than the values in other shales, and the SBET value (6.92 m2/g) is nearly equal to those of the other shales (Table 4). These values

162

indicate that the Da’anzhai shale can also have significant nanoscale storage space.

M AN U

161

To investigate the distribution and contributions of each pore type, the VBJH with respect to the PSD

164

was determined from the adsorption branches of the BJH model (Wang et al., 2015b; Jiang et al., 2016b;

165

Zhou et al., 2016; Zhang et al., 2017). The PSD plots show that pores with diameters between 3 and 5

166

nm make up the largest proportion of the total pore volume and that the proportions decrease with

167

increasing pore size, which is similar to some lacustrine shales (Fig. 6) (Jiang et al., 2016b).

EP AC C

168

TE D

163

74

169 170

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Fig. 7. Comparisons between the nanopore structures of the Da’anzhai shale and other lacustrine shales. A, B, C: the

Da’anzhai shale. D: Lacustrine shale of the Chang 7 Member in the Ordos Basin, China (Fu et al., 2015). E: Lacustrine shale

172

of the Songliao Basin, China (Wang et al., 2015b). F: Lacustrine shale of the Yanchang Formation in the Ordos Basin, China

173

(Jiang et al., 2016b; Wang et al., 2015c). G: Lacustrine shale of the Shahejie Formation in the Dongying Sag, China (Liu et

174

al., 2017). H: Lacustrine shale of the Shahejie Formation in the Zhanhua Sag, China. I: Lacustrine shale of the 5

175

Member of the Xujiahe Formation of Upper Triassic in Sichuan Basin, China (Peng, 2016).

176

AC C

EP

171

th

LTNA isotherms and hysteresis patterns can be used to investigate the physisorption mechanisms

177

and structural characteristics of shale pores (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015;

178

Zhou et al., 2016; Zhang et al., 2017). Adsorption isotherms are divided into five types (Types I-VI), and 75

ACCEPTED MANUSCRIPT

hysteresis loops are divided into four types (Types H1-H4) (Sing et al., 1985; IUPAC, 1994; Thommes et

180

al., 2015). Following this classification, the isotherms of the Da’anzhai shale samples are similar and

181

belong to Type П, and the hysteresis loops are characterized by Type H3 hysteresis loops that indicate

182

wedge-shaped pores (Sing et al., 1985; IUPAC, 1994; Thommes et al., 2015). However, compared to

183

other lacustrine shales, the Da’anzhai shale has relatively weak hysteresis loops and wide overlapping

184

ranges of adsorption and desorption isotherms (Fig. 7), which indicate many dead-end pores with highly

185

complex structures (De Boer, 1958; Nie et al., 2015).

M AN U

SC

RI PT

179

186

4.3.3. Factors that influence pore structure

187

SBET and VBJH are the two most basic LTNA parameters, and they can directly indicate the storage spaces of nanoscale reservoirs. Fig. 8a shows that the TOC values are positively correlated with SBET

189

and VBJH. However, the correlation is not strong because of the relatively low TOC values (average

190

0.97%). Fig. 8b shows that the lacustrine clay content is positively correlated with SBET and VBJH; this

191

correlation is mainly caused by the high clay content (42.4%) and well-developed intraparticle pores

192

(e.g., intraplatelet pores within clay aggregates) and interparticle pores (e.g., pores at the edges of

193

particles, crystals and clay platelets) (Figs. 5F, L). The quartz content is also positively correlated with

194

SBET and VBJH (Fig. 8c). Quartz particles in lacustrine shales usually occur within mineral matrices and

195

have well-developed interparticle pores at the edges of other minerals (e.g., clay and carbonate minerals).

196

The quartz particles also have well-developed nanoscale intraparticle pores. Additionally, the quartz

197

particles have no obvious cementation (Fig.5D). Conversely, the carbonate content is negatively

198

correlated with SBET and VBJH (Fig. 8d). Carbonate minerals, particularly calcite, can significantly affect

AC C

EP

TE D

188

76

ACCEPTED MANUSCRIPT

reservoir quality and generally has adverse effects (Wang et al., 2016). Carbonate minerals (e.g., calcite)

200

are easily dissolved and redistributed and can be formed in primary pores by various mechanisms during

201

the diagenesis processes (Figs. 5E, H, K) (Huang et al., 1984; Wang et al., 2016). Based on qualitative

202

observations using SEM, polarized light microscopy and fluorescence microscopy, lacustrine carbonate

203

minerals can effectively reduce the available storage space by diagenesis processes including

204

cementation, compaction, pressure-solution, recrystallization and replacement (Figs. 5E, H, K) (Huang

205

et al., 1984; Wang et al., 2016).

207

AC C

206

EP

TE D

M AN U

SC

RI PT

199

Fig. 8. Correlation plots between nanopore storage space indicators (SBET and VBJH) and values of TOC and minerals.

77

ACCEPTED MANUSCRIPT

4.4. Characterization and controlling factors of OM in lacustrine shale

209

4.4.1. Characterization of the OM content

210

Water conditions play an important role in the formation of source rocks. Lacustrine shales

RI PT

208

primarily form in semi-deep to deep lakes with fresh to low salinity water, whereas marine shales

212

primarily form in reductive, salty, weakly alkaline and low energy water. These features give rise to

213

different source rock characteristics (Huang et al., 1984; Bohacs et al., 2000; Curtis, 2002; Schmoker,

214

2002; Chen et al., 2015; Yang et al., 2015; Jiang et al., 2016a). Due to these differences, we evaluate the

215

lacustrine OM using a terrestrial oil generation theory that was derived from studies of continental

216

petroliferous basins in China and plays a special role in global petroleum theories (Huang et al., 1984;

217

Lu and Zhang, 2008) (Table 5).

218

Table 5 Factors used to evaluate OM in lacustrine source rock (Huang et al., 1984; Lu and Zhang, 2008). lake types freash~brackish

salt water~high salt water

AC C

chloroform bitumen ‘A’ (%)

non-source rock

source rock poor

medium

good

very good

0.4~0.6

>0.6~1.0

>1.0~2.0

>2.0

<0.2

0.2~0.4

>0.4~0.6

>0.6~0.8

>0.8

<0.015

0.015~0.050

<0.050~0.100

>0.100~0.200

>0.200

<2

2~6

>6~20

>20

<0.4

EP

TOC(%)

water

TE D

index

M AN U

SC

211

S1+S2(mg/g)

219

The mean values (and ranges) of the shale’s kerogen macerals consist of 32% exinites (29-35%),

220

27% sapropelinites (23-32%) and 24% vitrinites (18-26%), which indicates that the OM is primarily

221

derived from plankton and bacteria. The δ13CPDB values (-29.1 to -29.5‰, average -29.3‰) support this

222

finding (Table 7) (Huang et al., 1984; Lu and Zhang, 2008). KI is an effective indictor of OM type. The 78

ACCEPTED MANUSCRIPT

OM types (I, II1, II2, and III) can be identified using the KI ranges >>80-100, 80-40, 40-0 and <<0 to

224

-100, respectively (Huang et al., 1984; Lu and Zhang, 2008). The KI values for the Da’anzhai shale vary

225

from 2 to 15 with an average of 8.7, which indicates that the OM is chiefly sapropelic-humic type II2

226

and has significant potential to generate hydrocarbons. The mean values (and ranges) of the organic

227

elements C, H, and O make up 46.59% (31.12-61.41%), 3.53% (2.09-4.66%) and 1.56% (0.35-1.88%)

228

of the shale’s OM, respectively. The high H/C ratios (0.79-1.05, average 0.90) and low O/C ratios

229

(0.10-0.30, average 0.17) (Fig. 9a) indicate that the OM is type II and has relatively high hydrocarbon

230

potential (Huang et al., 1984; Lu and Zhang, 2008). In addition, the Tmax-HI index (Fig. 9b), which is an

231

important evaluation index for the OM type based on Rock-Eval pyrolysis, gives a similar result

232

(Espitaliéet al., 1985; Lu and Zhang, 2008; Misch et al., 2016). Generally, the Da’anzhai shale has

233

similar OM types with other lacustrine and marine shales and the dominant type is of type II (Table 1),

234

which can support the Da’anzhai shale hydrocarbon generation capacity.

TE D

M AN U

SC

RI PT

223

The TOC values of the Da’anzhai shale vary from 0.11 to 2.18% (average 0.97%), these values are

236

generally lower than those measured in the world's major lacustrine and marine shales (Table 1 and Fig.

237

9). The Ro values of the Da’anzhai shale are similar or even higher than those of other lacustrine and

238

marine shales, which indicate that the OM is highly mature. This could also be deduced from the Tmax

239

values, which vary from 428 to 500 ℃ with an average of 449 ℃ (Table 6).

AC C

EP

235

79

ACCEPTED MANUSCRIPT

Table 6 Analysis results for Rock-Eval, Chloroform Bitumen ‘A’ and TOC.

code

bitumen A

TOC

Tmax

S1+S2

(℃)

(mg/g)

S0 (mg/g)

S1 (mg/g)

S2 (mg/g)

0.82

0.009

0.08

0.50

456

1.33

0.006

2.03

3.97

444

(%)

(%) G10-1

0.03

G10-2 G4-1

0.02

0.40

0.007

0.02

0.16

461

G4-2

0.25

0.88

0.004

0.14

1.03

439

0.42

1.43

0.009

0.22

0.71

428

0.40

1.49

0.009

0.17

2.12

443

Li-1

0.35

1.17

0.008

0.29

1.09

441

Li-2

0.31

1.07

0.009

0.15

1.28

Li-3

0.16

0.69

0.006

0.07

D

HCI

(mg/g)

(%)

(%)

(mg/g)

0.58

0.14

60.60

0.05

5.95

11.04

6.00

0.34

298.77

0.50

37.51

153.11

0.18

0.12

40.65

0.02

3.98

7.30

1.18

0.12

117.53

0.10

11.13

16.61

0.24

49.69

0.08

5.48

16.33

0.07

142.48

0.19

12.79

11.67

1.38

0.21

93.29

0.12

9.86

25.49

446

1.43

0.10

119.80

0.12

11.17

14.82

0.55

445

0.62

0.11

79.33

0.05

7.48

10.77

0.006

1.70

1.85

452

3.55

0.48

116.50

0.30

18.58

107.39

0.24

1.51

0.008

0.16

0.94

453

1.09

0.14

61.96

0.09

6.05

10.88

P1-4

0.11

0.61

0.006

0.11

0.26

448

0.38

0.29

43.18

0.03

5.18

19.28

S2-2

0.27

1.18

0.006

0.14

1.32

444

1.46

0.10

112.14

0.12

10.34

12.41

S2-5

0.01

0.008

0.03

0.05

500

0.09

0.35

97.82

0.01

14.22

73.45

XI44-2

0.40

0.009

0.09

1.83

439

1.92

0.04

159.28

0.16

13.91

8.27

1.15

TE D

1.59

P1-2

EP

241 242

PC

2.29

AC C

P1-1

HI

0.94

M AN U

G4-3 J61-4

OPI

RI PT

chloroform

SC

240

80

ACCEPTED MANUSCRIPT

Table 7 Results for Ro, macerals, carbon isotopes and organic elements. maceral (%)

Li-1

Li-2

P1-1

(1.43) 1.19-1.47 (1.32) 0.78-1.08 (0.95) 0.95-1.25 (1.17) 0.88-1.28 (1.02)

27

4

4

26

34

24

27

3

2

30

35

18

32

4

4

22

30

22

28

6

6

19

31

27

6

2

21

29

23

5

6

22

33

vitrinite

KI

fusinite

15

RI PT

1.13-1.65

sporinite

pollinite

subtotal

normal

resinite

(‰)

11

SC

J53-1

(0.95)

sporo

amorphinite

organic elements

δ13CPDB

inertinite

C (%)

H (%)

O (%)

H/C

O/C

-29.2

39.89

3.07

0.75

0.92

0.1

20

11

-29.5

46.31

4.07

1.88

1.05

0.3

16

15

-29.2

49.26

3.87

1.37

0.94

0.2

26

15

9

-29.1

31.12

2.09

0.35

0.81

0.1

26

18

4

-29.5

61.41

4.66

1.62

0.91

0.2

26

18

2

-29.2

51.54

3.39

0.96

0.79

0.1

M AN U

G4-3

0.73-1.17

vitrinite

TE D

G10-2

Ro(%)

exinite

EP

code

sapropelinite

AC C

243 244

81

ACCEPTED MANUSCRIPT

RI PT

245

SC

246

Fig.9. Discrimination diagrams for the OM type.

247

Terrestrial oil-generation theory (Huang et al., 1984; Lu and Zhang, 2008) suggests the following

M AN U

248

conclusions regarding the lacustrine shale: (1) The TOC and chloroform bitumen ‘A’ values are medium

250

to very good, even though the values of S1+S2 (genetic potential) are poor to medium (Tables 5 and 6).

251

(2) Because of the strong hydrocarbon expulsion of the high maturity OM (Ro>0.5-0.7%), the

252

evaluation criteria for high maturity OM should be reduced. The shale with high Ro values (0.95-1.43%)

253

is highly mature. (3) Lacustrine carbonate can also generate hydrocarbons, which is an important feature

254

but easy to overlook. The petroleum generation threshold of carbonate (TOC: 0.12-0.4%) is lower than

255

that of mudstone (TOC: 0.4-0.5%), which means that carbonates can produce hydrocarbons at lower

256

TOC values. Thus, the low TOC values of the high carbonate shale may not indicate a low hydrocarbon

257

generation capacity. Additionally, carbonate interlayers in shale can play a role similar to that of seal

258

rocks, which are conducive to the preservation of oil and gas. In summary, the results of this study

259

indicate that the Da’anzhai shale has some hydrocarbon generation capacity, which is also supported by

260

the clear fluorescence of the thin sections (Figs. 5B, C).

AC C

EP

TE D

249

82

ACCEPTED MANUSCRIPT

261

Based on microscopic observations of the Da’anzhai and lacustrine shales, the OM occurrence forms are generally similar to other lacustrine shales; they are mainly dispergated in shale matrices

263

(Fig.5) (Morad et al., 2010; Zou et al., 2011; He et al., 2012; Luo et al., 2013; Tian et al., 2014; Dang et

264

al., 2015; Xu et al., 2017). However, these lacustrine shales have more amorphous forms than marine

265

shales because of their different water conditions and biodegradation (Huang et al., 1984; Lu and Zhang,

266

2008; Tian et al., 2014; Xu et al., 2017). In addition, lacustrine OM is distributed in concentric belts

267

around the lake basin’s center, and OM abundance increases with lake depth. Thus, OM is generally

268

more abundant in moderately deep to deep lake facies than in shallow lake facies, which can be

269

supported by the wide TOC range of lacustrine shales (Table 1 and Fig. 9) (Huang et al., 1984; Lu and

270

Zhang, 2008; Tian et al., 2014; Xu et al., 2017).

SC

M AN U

4.4.2. Factors that influence the OM content

AC C

EP

TE D

271

RI PT

262

272

83

ACCEPTED MANUSCRIPT

Fig. 10. Correlation plots between TOC and V/Cr, V/Sc, Ni/Co and V/(V+Ni) values.

M AN U

SC

RI PT

273

274

Fig. 11. Correlation plots between TOC values and minerals.

276

TOC is commonly used to represent OM. As discussed in Section 4.1, the reducibility is positively

TE D

275

correlated with the V/Cr, V/(V+Ni), and V/Sc ratios and negatively with the Ni/Co ratio. In addition,

278

pyrite is an important mineral in organic-rich sediments and is useful for reconstructing

279

paleoenvironmental conditions. A high pyrite content reflects a stable, high-salinity environment with

280

strong reducibility (Leventhal, 1983). The TOC values of the Da’anzhai shale are positively correlated

281

with the V/Cr, V/(V+Ni) and V/Sc ratios (Figs. 10a, b, c) and the pyrite content (Fig. 11a) and correlated

282

negatively with the Ni/Co ratio (Fig. 10d). These values indicate that paleo-redox conditions had a

283

significant effect on the formation and preservation of OM and that the strongly reducing conditions

284

were conducive to the preservation of OM (Zeng et al., 2015; Chen et al., 2016; Xu et al., 2017).

AC C

EP

277

84

ACCEPTED MANUSCRIPT

285

Similar to other lacustrine shales, the TOC values of the Da’anzhai shale are positively correlated with the clay mineral content (Fig. 11b), which is due to the OM’s strong sorption ability of clay and its

287

relatively high content for the lacustrine shales (Huang et al., 1984; Gingele et al., 2001; Dang et al.,

288

2015; Xu et al. 2017). The lacustrine clay minerals are small and rich in intraparticle and interparticle

289

pores, which creates storage space for OM (Figs. 5D,F,G,J). The OM observed within the clay minerals

290

primarily has complex, dispergated and residual biological shapes, which indicate that the OM and clay

291

minerals are in close contact (Figs. 5G and J) (Huang et al., 1984; Lu and Zhang, 2008; Tian et al., 2014;

292

Xu et al., 2017). Conversely, the carbonate content is negatively associated with OM (Fig. 11c), which is

293

similar to some lacustrine shales (Huang et al., 1984; Jarvis et al., 2001; Tian et al., 2014; Wang et al.,

294

2016). Carbonate minerals fill the pores via cementation, compaction, pressure-solution,

295

recrystallization and replacement, which effectively reduces the storage space for OM (Figs. 5E,H,K)

296

(Huang et al., 1984; Wang et al., 2016). Furthermore, the high carbonate mineral content in the

297

lacustrine facies reflects relatively shallow water environments, such as the shore-shallow lake face and

298

shallow lacustrine facies, which are not favorable for the generation or storage of OM (Huang et al.,

299

1984; Wang et al., 2016; Xu et al., 2017). There is a weak correlation between the quartz content and

300

TOC values (Fig. 11d). The correlations between quartz and OM are not uniform; differences in the

301

correlation between quartz and TOC contents are related to the depositional environment of shale

302

reservoirs (Zeng et al., 2014; Liu et al., 2015; Ross and Bustin, 2007; Chalmers et al., 2012b; Tian et al.,

303

2013). Negative correlations are caused by the terrestrial quartz inputs (Zeng et al., 2014; Liu et al.,

304

2015), while the positive correlations are related to biogenic quartz (Ross and Bustin, 2007; Chalmers et

AC C

EP

TE D

M AN U

SC

RI PT

286

85

ACCEPTED MANUSCRIPT

305

al., 2012b; Tian et al., 2013). Therefore, the quartz in the lacustrine shales may have both related to

306

biogenic and terrestrial origins. This could lead to the weak correlations between quartz and OM. Based on the analyses in Figs. 10 and 11, the low TOC values of the Da’anzhai lacustrine shales

RI PT

307

were primarily caused by the oxygen paleoenvironment and carbonate minerals. The relatively shallow

309

and turbulent lake water formed a high oxygen paleoenvironment with greater inputs of terrestrial

310

nutrients, which was conducive to the breeding of shell organisms and the formation of carbonate

311

minerals but not to the preservation of OM. In summary, the oxygen paleoenvironment and high

312

carbonate mineral content are primarily related to the unique lacustrine facies, which was not conducive

313

to the formation of OM.

316

M AN U

315

4.5. Relationships between the mineralogy, reservoirs, OM and paleoenvironmental conditions of lacustrine shales

TE D

314

SC

308

The mineral content of lacustrine shales is mainly related to the paleoenvironment. Strong paleo-weathering can lead to a greater number of weathering products, which is favorable to the

318

development of terrestrial minerals in lacustrine facies. Humid (rainy) paleoclimates with more water

319

can improve river transportation and reduce lake salinity, and low salinity paleoenvironments inhibit the

320

formation of authigenic carbonate. In addition, fine-grained terrestrial minerals, especially clay minerals,

321

inhibit the formation of authigenic carbonate (Leventhal, 1983; Huang et al., 1984; Bohacs et al., 2000;

322

Chi et al., 2003; Fu et al., 2015). High hydrodynamic conditions can effectively remove fine-grained

323

terrestrial minerals, which favors the formation of authigenic carbonate.

AC C

EP

317

86

ACCEPTED MANUSCRIPT

324

The lacustrine sedimentary environment is complex and the lacustrine shale formation is easily influenced by the paleoenvironment. The two most important differences are as follows: first, the

326

lacustrine paleoenvironments are complex and varied with the short detention time and strong variation

327

of water level (Huang et al., 1984; Chen et al., 2015) (Figs. 1 and 5C, Tables 2 and 6). Second, lakes are

328

close to the terrestrial sources with less water, which lead to rich terrestrial inputs with efficient impact

329

on shale (Figs. 1 and 5A,B, Tables 3 and 6) (Huang et al., 1984; Chen et al., 2015) (Huang et al., 1984;

330

Bohacs et al., 2000; Curtis, 2002; Schmoker, 2002). Lakes are generally located near paleo-uplift

331

environments with large topographical differences, and terrestrial debris can be transported directly into

332

a lake basin over short distances. These features result in rich, proximal and variable sources for

333

lacustrine facies. (Figs. 1, 5A, 5B; Table 3). The above two characteristics can effectively lead to the

334

complex and changeable mineral contents of lacustrine shales (Table 3; Fig. 4). In addition, there are

335

some other common points of the lacustrine shales including relatively small thickness of monolayer,

336

more siliceous or calcareous intercalations, and poor comparability with strong heterogeneity (Table 1),

337

which are also the results of the two characteristics of lacustrine sedimentation environments (Tian et al.,

338

2014; Wang et al., 2015a; Zhu et al., 2016; Xu et al., 2017).

SC

M AN U

TE D

EP

AC C

339

RI PT

325

Similar to other lacustrine shales, the storage space of the Da’anzhai lacustrine shale is directly

340

controlled by minerals (Fig. 12). Because of the large amount of terrestrial clastic minerals, the

341

interparticle pores of lacustrine shales are well-developed. Developed interparticle pores and intraparicle

342

pores can guarantee a considerable amount of storage space. In addition, numerous siliceous or

343

calcareous interlayers in lacustrine shale formations can result from varying changeable

87

ACCEPTED MANUSCRIPT

344

paleoenvironments with rich terrestrial inputs. These interlayers can increase the brittleness index and

345

facilitate the development of cracks. Compared to other shales, the Da’anzhai lacustrine shale has a relatively low OM content (Table 6),

RI PT

346

which is mainly the result of more variable paleoenvironments with high oxygen content (Fig. 12)

348

(Huang et al., 1984; Lu and Zhang, 2008; Tian et al., 2014; Xu et al., 2017). Different minerals have

349

different effects on storage spaces of OM, and minerals in lacustrine shales are mainly controlled by

350

paleoenvironments (Fig. 12). In addition, the relatively low OM content of the Da’anzhai lacustrine

351

shale is the reason for OM weak effects on storage spaces (Fig. 8a).

M AN U

SC

347

In summary, minerals, reservoirs, and OM in lacustrine shales are closely associated with each

353

other, and they were all controlled by the paleoenvironment which was the dominant factor and the link

354

between these parameters (Fig. 12).

EP

TE D

352

355

357 358

Fig.12 Relationships between mineralogy, reservoirs, OM and paleoenvironment for lacustrine shale

AC C

356

4. Conclusions

1. There are some similarities between the Da’anzhai lacustrine shale and other lacustrine shales. (1)

359

The two most important differences are that lacustrine shales formed in more complex and variable

360

paleoenvironments and have more abundant terrestrial inputs than marine shales. (2) The mineral

361

compositions of lacustrine shales are variable and complex, and they have abundant terrestrial minerals. 88

ACCEPTED MANUSCRIPT

In addition, lacustrine shales have numerous siliceous or calcareous interlayers. (3) The interparticle

363

pores and interlayer pores of lacustrine shales are well-developed. Unlike other lacustrine shales, the

364

Da’anzhai shale has relatively weak hysteresis loops and wide overlapping ranges of adsorption and

365

desorption isotherms, indicating a greater number of dead-end pores with more complex structures. In

366

addition, the Da’anzhai shale has lower TOC values with more amorphous OM forms.

2. The minerals, reservoirs, and OM for lacustrine shales are closely associated with each other;

SC

367

RI PT

362

these parameters are all controlled by the paleoenvironmental conditions in which the shales formed. (1)

369

Highly variable lacustrine paleoenvironment lead to variable mineral contents and abundant terrestrial

370

inputs give rise to a high terrestrial mineral content in lacustrine shales, which can also create a higher

371

number of siliceous or calcareous intercalations relative to marine shales. (2) Terrestrial minerals (e.g.,

372

clay and quartz) can effectively improve the storage spaces, but the authigenic carbonate minerals (e.g.,

373

calcite) have the opposite effect. Well-developed interparticle pores in lacustrine shales are caused by

374

abundant particles. (3) Clay content has a positive effect on the formation and preservation of OM, while

375

higher oxygenation and carbonate contents have negative effects. Lakes generally have relatively

376

shallow water and a variable paleoenvironment, which leads to a higher oxygen content and reduces the

377

OM abundance of lacustrine shales. The weak influence of OM on nanoscale storage spaces is due to the

378

low TOC values of the lacustrine shales.

TE D

EP

AC C

379

M AN U

368

3. Compared to other lacustrine shales, the Da’anzhai shale also demonstrates the potential for

380

unconventional oil and gas exploration. (1) The lacustrine shale has well-developed storage spaces due

381

to the high number of interparticle pores, intraparticle pores, and fractures. Although the lacustrine shale

89

ACCEPTED MANUSCRIPT

has more dead-end pores and larger pores with more complex microstructures, the types of isotherms

383

(Type П) and hysteresis loops (Type H3) are generally similar to those of lacustrine shales. The values of

384

VBJH (1.19-4.10 m3/100 g, average 2.97 m3/100 g) and SBET (1.75-10.69 m2/g, average 6.92 m2/g) are

385

similar to other shales, indicating a high amount of storage space in nanoscale reservoirs. (2) Abundant

386

terrestrial minerals in lacustrine facies can support a greater number of interparticle pores and

387

intraparticle pores. Numerous siliceous or calcareous interlayers in lacustrine shales can also improve

388

the brittleness index, facilitating the development of fractures and creating a seal for oil and gas. (3)

389

Although the TOC values of the Da’anzhai shale are lower than other lacustrine shales, the OM still has

390

several advantages according to the evaluation criteria for continental source rocks, including moderate

391

values of chloroform bitumen ‘A’ and TOC, a high level of OM maturity and a favorable kerogen type

392

(II).

393

Acknowledgments

TE D

M AN U

SC

RI PT

382

Our study is supported by the Institute Program of PetroChina Research Institute of Petroleum

395

Exploration and Development (Grant 2016yj01), and the National Natural Science Foundation of China

396

(Grants 41272137 and 41572117). The authors also sincerely appreciate the support from the China

397

Scholarship Council (CSC).

398

Nomenclature

399

TOC

400

Ro

Vitrinite reflectance

401

Tmax

Pyrolysis Temperature at Maximum Hydrocarbon Generation

AC C

EP

394

Total organic carbon, %

90

ACCEPTED MANUSCRIPT

SBET

Brunauer-Emmett-Teller specific surface area

403

VBJH

Barrett-Joyner-Halenda total pore volume

404

DA

Average pore diameter

405

OM

Organic matter

406

TE

Trace element

407

XRD

X-ray diffraction

408

FE-SEM

Field emission scanning electron microscopy

409

LTNA

Low temperature N2 adsorption

410

EDS

Energy dispersive spectrometry

411

S0

Gaseous hydrocarbons, mg/g

412

S1

Free hydrocarbons, mg/g

413

S2

Pyrolysis of hydrocarbons, mg/g

414

ICP-MS

Inductively coupled plasma-mass spectrometry

415

REE

Rare earth element, ppm

416

ΣHREE

Total content of heavy rare earth elements

417

ΣLREE

Total content of light rare earth elements

418

PSD

Pore size distribution

419

OPI

420

HI

Hydrogen index; HI=(S2/TOC)×100%, mg/g

421

PC

Effective carbon content; PC=0.083×(S0+S1+S2), %

AC C

EP

TE D

M AN U

SC

RI PT

402

Oil production index; OPI=S1/(S0+S1+S2)

91

ACCEPTED MANUSCRIPT

D

Degradation rate; D=(PC/TOC) ×100%, %

423

HCI

Hydrocarbon index; HCI=(S0+S1)/ TOC×100, mg/g

424

KI

Kerogen type index; KI=Sapropelinite (%)×1+Exinite (%)×0.5+Vitrinite

RI PT

422

(%)×(-0.75)+Inertinite (%)×(-1) (Huang et al., 1984; Lu and Zhang, 2008)

426

References

427

Barrett, E.P., Joyner, L.G., Halenda, P. P., 1951. The determination of pore volume and area distributions in porous

429

substances. I. Computations from nitrogen isotherms. J. Am. Chem. Soc. 73(1), 373-380.

M AN U

428

SC

425

Bohacs, K.M., Carroll, A.R., Neal, J.E., Mankiewicz, P.J., 2000.Lake-basin type, source potential, and

430

hydrocarbon character: an integrated sequence-stratigraphic-geochemical framework. Lake basins through

431

space and time: AAPG Studies in Geology 6, 3-34.

435 436

TE D

434

309-319.

Bowker, K.A., 2007. Barnett shale gas production, Fort Worth Basin: issues and discussion. AAPG Bull. 91(4), 523-533.

EP

433

Brunauer, S., Emmett, P.H., Teller, E.J., 1938. Adsorption of gases in multimolecular layers. Am. Chem. Soc. 60,

Chalmers, G.R., Bustin, R.M., Power, I.M., 2012. Characterization of gas shale pore systems by porosimetry,

AC C

432

437

pycnometry, surface area, and field emission scanning electron microscopy/transmission electron

438

microscopy image analyses: Examples from the Barnett, Woodford, Haynesville, Marcellus, and Doig units.

439

AAPG Bull. 96(6), 1099-1119.

92

ACCEPTED MANUSCRIPT

Chi, G., Giles, P.S., Williamson, M.A., Lavoie, D., Bertrand, R., 2003. Diagenetic history and porosity evolution

441

of Upper Carboniferous sandstones from the Spring Valley #1 well, Maritimes Basin, Canada-implications

442

for reservoir development. J. Geochem. Explor. 80, 171-191.

443

RI PT

440

Chen, C., Mu, C.L., Zhou, K.K., Liang, W., Ge, X.Y., Wang, X.P., Wang, Q.Y., Zheng, B.S., 2016. The geochemical characteristics and factors controlling the organic matter accumulation of the Late

445

Ordovician-Early Silurian black shale in the Upper Yangtze Basin, South China. Mar. Pet. Geol. 76, 159-175.

447 448

Chen, L., Liu, C.L., Zhuang, C., Che, X.G., Wu, J., 2009. Rare Earth Element records of the lower Paleogene

M AN U

446

SC

444

sediments in the Sanshui Basin and their paleoclimate implications. Acta Sedimentol. Sin. 27 (6), 1155-1162. Chen, S.J.; Zhang, H.X., Lu, J.G, Yang, Y.M., Liu, C.W., 2015. Controlling factors of Jurassic Da’anzhai Member tight oil accumulation and high production in central Sichuan Basin, SW China. Petrol. Explor. Dev. 42(2),

450

206-214.

TE D

449

451 Curtis, J. B., 2002. Fractured shale-gas systems. AAPG Bull. 86(11), 1921-1938. 452 Dang, W., Zhang, J.C., Huang, X., Li, X.G., Chen, Q., Sun, R., Xue, B., Han, S.B., 2015. Main-controlling

geological factors of gas-bearing property of continental shale gas: a case study of Member 3rd of Shahejie

454

Formation in western Liaohe sag. Acta Petrolei Sin. 36(12), 1516-1530.

AC C

EP

453

455 De Boer, J.H., 1958. The structure and properties of porous materials. In: Proceedings of the tenth symposium of the 456

colston research society held in the University of Bristol, Butterworths, London, pp. 68-94.

457

Espitalié, J., Deroo, G., Marquis, F., 1985. La pyrolyse Rock-Eval et ses applications. Deuxième partie. Revue de

458

l'Institut français du Pétrole 40(6), 755-784.

93

ACCEPTED MANUSCRIPT

Fu, H.J., Wang, X.Z., Zhang, L.X., Gao, R.M., Li, Z.T., Xu, T., Zhu, X.L., Xu, W., Li, Q., 2015. Investigation of

460

the factors that control the development of pore structure in lacustrine shale: A case study of block X in the

461

Ordos Basin, China. J. Nat. Gas. Sci Eng. 26, 1422-1432.

RI PT

459

Gingele, F.X., De Deckker, P., Hillenbrand, C.D., 2001. Late Quaternary fluctuations of the Leeuwin Current and

463

palaeoclimates on the adjacent land masses: clay mineral evidence. Aust. J. Earth Sci. 48(6), 867-874.

464

Gromet, L.P., Haskin, L.A., Korotev, R.L., Dymek, R.F., 1984. The “North American shale composite”: its compilation, major and trace element characteristics. Geochim. Cosmochim. Ac.

48(12), 2469-2482.

M AN U

465

SC

462

466 Han, S.B., Zhang, J.C., Li, Y.X., Hosrfield, B., Tang, X., Jiang, W.L., Chen, Q., 2013. Evaluation of Lower

Cambrian shale in northern Guizhou Province, South China: Implications for shale gas potential. Energ. Fuel 27,

468

2933-2914.

469

He, F.Q., Zhu, T., 2012. Favorable targets of breakthrough and built-up of shale gas in continental facies in Lower

470 471

Jurassic, Sichuan Basin. Petrol. Geol. Exper. 34(3), 246-251. Huang, D.F., Li, J.C., Zhang, D.J., 1984. Evolution and hydrocarbon generation mechanisms of Terrestrial organic matter. Petroleum Industry Press: Beijing, pp. 1-228.

EP

472

TE D

467

AC C

473 IUPAC (International Union of Pure and Applied Chemistry), 1994. Physical Chemistry Division Commission on 474

Colloid and Surface Chemistry, Subcommittee on Characterization of porous solids. Recommendations for the

475

characterization of porous solids (Technical Report). Pure Appl. Chem. 66 (8), 1739-1758.

476

Jarvie, D.M., Hill, R.J., Ruble, T.E., 2007. Unconventional shale-gas systems: The Mississippian Barnett Shale of

477

north-central Texas as one model for thermogenic shale-gas assessment. AAPG Bull. 91, 475-499.

94

ACCEPTED MANUSCRIPT

478 479

Jarvis, I.A.N., Murphy, A.M., Gale, A.S., 2001. Geochemistry of pelagic and hemipelagic carbonates: criteria for identifying systems tracts and sea-level change. J. Geol. Soc. 158(4), 685-696. Jiang, S., Xu, Z.Y., Feng, Y.L., Zhang, J.C., Cai, D.S., Chen, L., Wu, Y., Zhou, Z.S., Bao, S.J., Long, S.X., 2016a.

481

Geologic characteristics of hydrocarbon-bearing marine, transitional and lacustrine shales in China. J. Asian

482

Earth Sci. 115, 404-418.

RI PT

480

Jiang, F.J., Chen, D., Wang, Z.F., Xu, Z.Y., Chen, J., Liu, L., Huyan, Y.Y., Liu, Y., 2016b. Pore characteristic

484

analysis of a lacustrine shale: A case study in the Ordos Basin, NW China. Mar. Pet. Geol. 73, 554-571.

486

M AN U

485

SC

483

Jiang, C., Bryant, S., Daigle, H., 2015. A bundle of short conduits model of the pore structure of gas shale. Unconventional Resources Technology Conference (URTEC).

Jiang, S., Tang, X., Cai, D., Xue, G., He, Z., Long, S., Peng, Y., Gao, B., Xu, Z., Dahdah, N., 2017. Comparison of

488

marine, transitional, and lacustrine shales: A case study from the Sichuan Basin in China. J. Petrol. Sci. Eng.,

489

150, 334-347.

492 493 494 495 496 497

conditions in ancient mudstones. Chem. Geol. 111(1-4), 111-129.

EP

491

Jones, B., Manning, D.A.C., 1994. Comparison of geochemical indices used for the interpretation of palaeoredox

Kimoto, A., Nearing, M.A., Shipitalo, M.J., Polyakov, V.O., 2006. Multi℃year tracking of sediment sources in a

AC C

490

TE D

487

small agricultural watershed using rare earth elements. Earth Surf. Proc. Land. 31(14) 1763-1774. Kimura, H., Watanabe, Y., 2001. Oceanic anoxia at the Precambrian-Cambrian boundary. Geology 29(11), 995-998.

Leventhal, J.S., 1983. An interpretation of carbon and sulfur relationships in Black Sea sediments as indicators of environments of deposition. Geochim. Cosmochim. Ac. 47(1), 133-137.

95

ACCEPTED MANUSCRIPT

498

Li, A., Ding, W.L., He, J.H., Dai, P., Yin, S., Xie, F., 2016. Investigation of pore structure and fractal

499

characteristics of organic-rich shale reservoirs: A case study of Lower Cambrian Qiongzhusi Formation in

500

Malong block of eastern Yunnan Province, South China. Mar.

502

RI PT

501

Pet. Geol. 70, 46-57.

Li, J.Y., 2014. Oil and gas contents and movable oil amounts of shales in 4th member of Shahejie Formation, Lijin sub Sag, Dongying Sag. Petrol. Geol. Exp. 36(3), 365-369.

SC

503 Li, P.P., Hao, F., Guo, X.S., Zou, H.Y., Yu, X.Y., Wang, G.W., 2015. Processes involved in the origin and

accumulation of hydrocarbon gases in the Yuanba gas field, Sichuan Basin, southwest China. Marine and

505

Petroleum Geology 59, 150-165.

506

Lu, S.F., Zhang, M., 2008. Petroleum Geochemistry. Petroleum Industry Press: Beijing, pp. 200-246.

507

Lu, B.X., Zheng, R.C., Liang, X.W., Xu, W.L., 2014. Characteristics analysis of Da’anzhai shale gas (oil)

511 512 513 514

TE D

510

Luo, P., Ji, L.M., 2013. Reservoir characteristics and potential evaluation of continental shale gas. Nat. Gas Geosci. 24(5), 1060-1068.

Liu, Y., Lu, Z.Y., Feng, M.S., Wang, J., Tian, T.H., Chao, J., 2017. Micro-pore characteristics of shale oil

EP

509

reservoirs in eastern Sichuan Basin. Geology in China 41, 1387-1398.

reservoirs of the Shahejie Formation in the Dongying Sag, Bohai Bay Basin. Acta Geol. Sin. 91(3), 629-644.

AC C

508

M AN U

504

Martini, A. M., Walter, L.M., McIntosh, J.C., 2008. Identification of microbial and thermogenic gas components from Upper Devonian black shale cores, Illinois and Michigan basins. AAPG Bull. 92(3), 327-339.

515

Morad, S., Al-Ramadan, K., Ketzer, J.M., De Ros, L.F., 2010. The impact of diagenesis on the heterogeneity of

516

sandstone reservoirs: a review of the role of depositional facies and sequence stratigraphy. AAPG Bull. 94,

517

1267-1309.

96

ACCEPTED MANUSCRIPT

521 522 523 524 525

Milliken, K.L., Rudnicki, M., Awwiller, D.N., Zhang, T., 2013. Organic matter-hosted pore system, Marcellus Formation (Devonian), Pennsylvania. AAPG Bull. 97 (2), 177-200.

RI PT

520

Upper Visean Black Shales in the Dniepr-Donets Basin (Ukraine). Mar. Pet. Geol. 75, 203-219.

Murray, R.W., ten Brink, M.R.B., Jones, D.L., Gerlach, D.C., Russ, G.P., 1990. Rare earth elements as indicators of different marine depositional environments in chert and shale. Geology 18(3), 268-271.

SC

519

Misch, D., Gross, D., Mahlstedt, N., Makogon, V., Sachsenhofer, R.F., 2016. Shale gas/shale oil potential of

Nie, B.S., Liu, X.F., Yang, L.L., Meng, J.Q., Li, X.C., 2015. Pore structure characterization of different rank coals

M AN U

518

using gas adsorption and scanning electron microscopy. Fuel 158, 908-917.

526 Potter, P., Maynard, J., Pryor, W., 1980. Sedimentology of Shale. Springer: New York, pp. 98. 527 Peng, J.L., Zhang, J.C., Wei, P.F., Long, S., Huang, X., Huang, Z.L., Ma, Y.L., Zhang, M.Q., 2016. Characteristics

and Controlling Factors of Micro-pore Structures of the 5th Shale Member of the Xujiahe Formation of Upper

529

Triassic in Sichuan Basin. Geoscience 30 (4), 896-904.

530

Reheis, M.C., 1990. Influence of climate and eolian dust on the major-element chemistry and clay mineralogy of

533 534 535 536 537

EP

532

soils in the northern Bighorn Basin, USA. Catena 17(3), 219-248. Rimmer, S.M., 2004. Geochemical paleoredox indicators in Devonian-Mississippian black shales, central

AC C

531

TE D

528

Appalachian Basin (USA). Chem. Geol. 206(3), 373-391. Ross, G. R., Guevara, S. R., Arribere, M. A., 1995. Rare earth geochemistry in sediments of the Upper Manso River basin, Rio Negro, Argentina. Earth Planet. Sc. Lett. 133(1-2), 47-57. Roy, P. D., Smykatz-Kloss, W., 2007. REE geochemistry of the recent playa sediments from the Thar Desert, India: an implication to playa sediment provenance. Chem. Erde-Geochem. 67(1), 55-68.

97

ACCEPTED MANUSCRIPT

539 540

Schmoker, J. W., 2002. Resource-assessment perspectives for unconventional gas systems. AAPG Bull. 86(11), 1993-1999. Shi, B.H., Yao, J.L., Zhang, Y., Zhang, L., Yang, Y.J., Li, H., 2013. Geologic conditions for hydrocarbon

RI PT

538

541

accumulation in Chang-9 oil reservoir of the Yanchang Formation, Ordos Basin. Oil gas Geol. 34(3),

542

294-300.

Sing, K.S., Everett, D.H., Haul, R.A.W., Moscou, L., Pierotti, R.A., Rouquerol, J., Siemieniewsha, T., 1985.

SC

543

Reporting physisorption data for gas/solid systems with special reference to the determination of surface area

545

and porosity. Pure Appl. Chem. 57, 603-619.

546 547

M AN U

544

Stevens, C.J., Quinton, J.N., 2008. Investigating source areas of eroded sediments transported in concentrated overland flow using rare earth element tracers. Catena 74(1), 31-36.

Thommes, M., Kaneko, K., Neimark, A. V., Olivier, J. P., Rodriguez-Reinoso, F., Rouquerol, J., Sing, K. S., 2015.

549

Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution

550

(IUPAC Technical Report). Pure Appl. Chem. 87(9-10), 1051-1069. Tian, C.T., Ma, S.P., Yang, Y., Meng, J., Xia, Y.Q., 2014. Comparison of hydrocarbon generation conditions

EP

551

TE D

548

between lacustrine and marine carbonate source rocks. Oil Gas Geol. 35(3), 336-341. Vanderaveroet, P.,

553

2000.Miocene to Pleistocene clay mineral sedimentation on the New Jersey Shelf. Oceanol. Acta 23(1),

554

25-36.

555

AC C

552

Veizer, J., Demovic, R., 1974. Strontium as a tool in facies analysis. J.

98

Sediment. Res. 44(1), 93-115.

ACCEPTED MANUSCRIPT

556

Wang, J., Cao, Y.C., Liu, K.Y., Liu, J., Xue, X.J., Xu, Q.S., 2016. Pore fluid evolution, distribution and water-rock interactions of carbonate cements in red-bed sandstone reservoirs in the Dongying Depression, China. Mar.

558

Pet. Geol. 72, 279-294.

559

RI PT

557

Wang, M., Wilkins, R. W., Song, G., Zhang, L., Xu, X., Li, Z., Chen, G., 2015a. Geochemical and geological characteristics of the Es3 L lacustrine shale in the Bonan sag, Bohai Bay Basin, China. In. J.Coal Geol. 138,

561

16-29. Wang, M., Xue, H.T., Tian, S.S., Wilkins, R.W.T., Wang, Z.W., 2015b. Fractal characteristics of

562

Upper Cretaceous lacustrine shale from the Songliao Basin, NE China. Mar. Pet. Geol. 67, 144-153.

563

Wang, Y., Zhu, Y.M., Wang, H.Y., Feng, G.J., 2015c. Nanoscale pore morphology and distribution of lacustrine

564

shale reservoirs: Examples from the Upper Triassic Yanchang Formation, Ordos Basin. J. Energy Chem.

565

24(4), 512-519.

M AN U

SC

560

Wu, Y., Fan, T., Zhang, J., Jiang, S., Li, Y., Zhang, J., Xie, C., 2014. Characterization of the upper Ordovician and

567

lower Silurian marine shale in northwestern Guizhou province of the upper Yangtze block, South China:

568

Implication for shale gas potential. Energy Fuel 28(6), 3679-3687.

571

EP

570

Wright, J., Seymour, R.S., Shaw, H.F., 1984. REE and Nd isotopes in conodont apatite: variations with geological age and depositional environment. Geol. Soc.

Am. Spec. Pap. 196, 325-340.

AC C

569

TE D

566

Xu, Q., Liu, B., Ma, Y., Song, X., Wang, Y., Xin, X., Chen, Z., 2017. Controlling factors and dynamical formation

572

models of lacustrine organic matter accumulation for the Jurassic Da'anzhai Member in the central Sichuan

573

Basin, southwestern China. Mar. Pet. Geol. 86, 1391-1405.

574

Yang, C., Zhang, J., Tang, X., Ding, J., Zhao, Q., Dang, W., Chen, H., Su, Y., Li, B., Lu, D., 2017. Comparative

575

study on micro-pore structure of marine, terrestrial, and transitional shales in key areas, China. Int. J. Coal

99

ACCEPTED MANUSCRIPT

Geol. 171, 76-92. Yang, Z., Hou, L.H., Tao, S.Z., Cui, J.W., Wu, S.T., Lin, S.H., Pan, S.Q., 2015. Formation

577

and “sweet area” evaluation of liquid-rich hydrocarbons in shale strata. Petrol. Explor. Dev. 42(5), 609-620.

578

Zeng, S., Wang, J., Fu, X., Chen, W., Feng, X., Wang, D., Song, C., Wang, Z., 2015. Geochemical characteristics,

RI PT

576

579

redox conditions, and organic matter accumulation of marine oil shale from the Changliang Mountain area,

580

northern Tibet, China. Mar. Pet. Geol. 64, 203-221.

Zhang, J.Z., Li, X.Q., Wei, Q., Cao, W.J., Liang, W.L., Wang, Z., Wang, F.Y., 2017. Quantitative characterization

582

of pore-fracture system of organic-rich marine-continental shale reservoirs: A case study of the Upper

583

Permian Longtan Formation, Southern Sichuan Basin, China. Fuel 200, 272-281.

M AN U

584

SC

581

Zhang, Q., Zhu, X.M., Li, C.X., Qiao, L.J.Y., Liu, C., Mei, X.H., Du, H.Y., Lu, J.J., 2016. Classification and quantitative characterization of microscopic pores in organic-rich shale of the Shahejie Formation in the

586

Zhanhua Sag, Bohai Bay Basin. Oil Gas Geol. 37(3), 422-438.

587

TE D

585

Zhou, L., Kang, Z.H., 2016. Fractal characterization of pores in shales using NMR: A case study from the Lower Cambrian Niutitang Formation in the Middle Yangtze Platform, Southwest China. J. Nat. Gas Sci. Eng. 35,

589

860-872.

591 592 593

Zou, C.N., Dong, D.Z., Yang, H., Wang, Y.M., Huang, J.L., Wang, S.F., Fu, C.X., 2011. Conditions of shale gas

AC C

590

EP

588

accumulation and exploration practices in China. Natur. Gas Ind. 31, 26-39.

Zhu, T., Long, S.X., Wang F., Peng, Y.M., 2016. Sedimentary models and lithofacies types of lacustrine mud shale in the Sichuan Basin. Natur. Gas Ind. 36(8), 22-28.

594

100

ACCEPTED MANUSCRIPT Highlights

1. The Da'anzhai Member is systematically studied as a shale system for the first time.

RI PT

2. The mechanisms of lacustrine shale paleoenvironment and mineral are discussed.

3. The mechanisms of lacustrine shale reservoirs and OM are discussed.

SC

4. The relationships between paleoenvironment, mineral, reservoir and OM are

AC C

EP

TE D

M AN U

discussed.