Hydrocarbon Generation Kinetics of a Heterogeneous Source Rock

Sinopec Shengli Oilfield Company, 3 Liaocheng Road, Dongying, Shandong 257015, China. Energy Fuels , 2017, 31 (12), pp 13291–13304. DOI: 10.1021/acs...
0 downloads 10 Views 3MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Hydrocarbon generation kinetics of a heterogeneous source rock system: - Example from the Eocene-Oligocene Shahejie Formation, Bohai Bay Basin in eastern China Zhuoheng Chen, Maowen Li, Tingting Cao, Xiaoxiao Ma, ZHiming Li, Qigui Jiang, Zheng Li, and Chunqing Jiang Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b02361 • Publication Date (Web): 03 Nov 2017 Downloaded from http://pubs.acs.org on November 3, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Hydrocarbon generation kinetics of a heterogeneous source rock system:

2

- Example from the Lacsutrine Eocene-Oligocene Shahejie Formation,

3

Bohai Bay Basin, China

4

Zhuoheng Chen1, Maowen Li2, Tingting Cao2, Xiaoxiao Ma3, Zhiming Li2, Qigui Jiang2, Zheng

5

Li4 and Chunqing Jiang1 1

6 7

Geological Survey of Canada, 3303-33 Street NW Calgary,Alberta T2L 2A7,Canada;

2

China State Key Laboratory of Shale Oil and Shale Gas Resources and Effective Development,

8

Sinopec Petroleum Exploration and Production Research Institute,31 Xueyuan Road,Haidian

9

District, Beijing 100083,China;

10

3

China University of Petroleum (Beijing), 18 Fuxue Road,Changping District, Beijng 102249,

11

China; 4Sinopec Shengli Oilfield Company, 3 Liaocheng Road, Dongying, Shandong 257015,

12

China

13

Abstract

14

Source rock heterogeneity is a common and important feature that requires

15

attention in shale oil resource evaluation. The Eocene-Oligocene Shahejie Formation

16

in the Jiyang Super-basin, Bohai Bay Basin, has been studied to investigate vertical

17

variations in source rock characteristics and their thermal decomposition behaviours

18

in a relation to lacustrine shale oil resource potential. Fifty-four core samples from

19

Luo 69 Well were analyzed using Rock-Eval 6 equipment and the resulting FID

20

(Flame Ionization Detector) hydrocarbon pyrograms were examined to study their

21

generation kinetics. Differences in bulk geochemical characteristics, petroleum

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22

generation capacities and petrophysical properties allow for sub-division of the

23

samples into three “homogenous” groups: a) a stratigraphically higher facies,

24

corresponding to samples with HI value around 600 mg HC/g TOC and high S1/TOC

25

(either due to exogenous hydrocarbon input or thermally less stable kerogen, or both;

26

b) stratigraphically a facies in the middle of the section, representing thermally more

27

stable facies group with similar generation potential, but much lower S1/TOC; and c)

28

a stratigraphically lower facies, corresponding to source rock with low TOC (100 mg HC/g TOC) (Fig. 3b). The

279

samples with the lowest Tmax values (Group III) occur in the deepest part of the well

280

and show the lowest HI values (Figs. 3c, 3b and 4). It appears that samples in Group I

281

and II align in a single straight line on a TOC-S2 cross plot (Fig. 3d), indicating

282

similar generation potential. Samples in Group III fall along another trend with lower

283

generation potential than those in Groups I and II, which are supported by the

284

Tmax-HI plot in Fig. 3c The samples with lowest Tmax values in Group III are

285

typical Type III kerogen by their low hydrocarbon generation potential (fig. 3c) and

286

low HI and relative high OI values (Fig. 3 b).

287

ACS Paragon Plus Environment

Page 16 of 46

Page 17 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

288 289

Figure 3. Cross-plots of various Rock-Eval parameters showing geochemical characteristics of the

290

Es3x samples from Luo 69 Well.

291

Figure 4a is a composite diagram displaying geochemical data and petrophysical

292

properties derived from well logs against depth to reveal the variation of these

293

attributes in vertical direction. Figure 4b links typical petrophysical log responses to

294

organic rich rock types, showing the diversity of mineral components and texture in

295

the source rock unit in the study. The stratigraphic interval can be divided into five

296

sub-lithological-geochemical units: a) organic-rich laminated argillaceous limestone

297

(sub-unit A: 2990 – 3020m); b) calcareous mudstone with variable clay mineral

298

contents (sub-unit B: 3020 - 3040 m); c) organic-rich laminated argillaceous

299

limestone (sub-unit C: 3040 – 3068 m); d) low porosity massive to thin bedded

300

limestone (sub-unit D: 3068 – 3123 m); and e) calcareous mudstone/shaly limestone

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

301

(sub-unit E: 3123-3130 m). The geochemical characteristics of samples from sub-unit

302

A and C show similarities and differences. Samples in sub-unit A have much higher

303

S1/TOC ratios than those in sub-unit C though TOC and S2/TOC are comparable;

304

while the carbonate content and Tmax in the upper part are lower. Petrophysical

305

differences between the A and C sub-units are observable, slightly higher Gamma Ray

306

as indicated by higher clay content, higher density and lower sonic transit time shown

307

by lower porosity estimates, and lower resistivity indicated by the separation of

308

measured and calculated sonic (the second last column of Fig. 4a) in sub-unit A than

309

those of sub-unit C.

310

The sub-units D and E contain low organic content with TOC lower than 2%

311

with variable remaining generation potentials and relatively high percentage of free

312

hydrocarbon content (S1). The lithology is different from the overlying sub-units

313

indicated by low estimated porosity and lack of obvious ∆T separation (the second

314

last column of Fig. 4a).

315

ACS Paragon Plus Environment

Page 18 of 46

Page 19 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

316 317

Figure 4a. Depth profiles of Rock-Eval and petrophysical parameters showing different

318

geochemical and petrophysical characteristics of the source rock at different depths. The ∆T

319

separation is the difference between measured sonic transit time and estimated sonic transit time

320

from resistivity log proposed by Passey et al.24. Large separation indicates potential good source

321

rock interval. Porosity curves are smoothed to show the general trends.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

322 323

Figure 4b Composite profile of petrophysical well logs with Photomicrographs showing

324

vertical source rock heterogeneity in the source rock interval of Sha3 Member (Es3) of

325

Shahejie Formation in this study. GR (column 1) indicates relative abundance in clay content;

326

SP (spontaneous potential) indicates relative permeability; PE (photoelectronic Index, red

327

solid in second column) show carbonate content (3-5 from dolomite to calcite); Resistivity:

328

oil and gas (organic matter) indicator; Porosity logs (column 4). Photomicrographs showing a)

329

organic-rich muddy limestone. Thin clay-rich and carbonate layers interbedded with

330

TOC>2%; b) organic-rich laminated argillaceous limestone containing solid bitumen in

331

inter-granular pores of calcite grains (3057.6 m); c) thin carbonate layer with organic-rich

332

laminae. Organic matter shown as solid bitumen in horizontal fractures (3 059.35 m); d)

333

interbedded clay and carbonate thin layers (3104.4m) and TOC is 600 mg HC/g TOC from Group I, which is typical of

539

immature Type I/II source rock from a lacustrine setting in China.

540

requires higher temperature to convert to hydrocarbons, and the generation

ACS Paragon Plus Environment

39

38

. While in

Type I kerogen

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

541

temperature is expected to be even higher under a high heating rate at 3°C/millions in

542

this study. The 50% conversion at about 175°C is reasonable for the Group I.

543

Composite Generation Kinetics

544

A single immature sample represents a specific facies and contains little

545

information on the diversity of organic matters from different facies. A hydrocarbon

546

generation kinetics model based on a single immature sample, therefore, may not be

547

able to capture the variation in kinetic properties due to source rock heterogeneity. In

548

contrast, multiple samples from either vertically stacked or laterally extended

549

stratigraphic sections offer better chance to reveal the diversity in organic matter and

550

allow for characterization of their variability spatially. The proposed composite

551

kinetics model addresses the issue of source rock heterogeneity by using a synthetic

552

pyrogram of aggregated samples from vertically stacked or/and laterally extended

553

sections that represent the diversity in organic matters due to facies change and

554

variable terrestrial inputs. Because samples in sub-units A and D+E are either in early

555

hydrocarbon windows or partly contaminated by exogenous hydrocarbons, they

556

cannot be used for construction of a kinetic model without corrections40, 41. As such,

557

the immature samples of 22 pyrograms from post-solvent extraction replicates in

558

stratigraphic sub-unit C (Group I) have been used to construct a composite kinetics

559

model to demonstrate its application to a heterogeneous source rock system.

560

Under laboratory constant heating rate (25 °C/min), the calculated hydrocarbon

561

generation rates are similar among the 22 samples in sub-unit C (Figs. 13a). The

562

transformation ratios show a moderate variation (maximum 25 °C) at onset of

ACS Paragon Plus Environment

Page 34 of 46

Page 35 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

563

massive hydrocarbon generation (TR=10%). The difference in temperature becomes

564

small (only 5 °C) when TR reaches 90% (Figs. 13b). However, when geological

565

heating rate is applied (3 °C/million years), the generation rates behaved differently

566

(Fig. 13c), and show two distinct groups (Figs. 13c), one with the peak rates of around

567

165 °C and the other about 175 °C. The transformation ratios also show two distinct

568

groups, with one displaying onset temperature of massive generation at about 150 °C

569

and offset at 175 °C, and the other group with onset at about 160°C and offset at

570

180 °C (Figs. 13d). The maximum differences in temperature at the onset and offset

571

of generation among the 22 samples are 30 °C (150-180 °C). The offset temperature

572

at geological heating rate does not seem to converge as those derived from laboratory

573

heating rate (Fig. 12b), and large difference remains (Fig. 13d). In the proposed

574

composite model, the overall average behaviour of thermal degradation is expressed

575

as the mean production (conversion) rate and mean transformation ratio as shown in

576

Figure 13, while the maximum temperature difference in generation is shown by a

577

temperature envelope defined by the lower and upper boundaries. The temperature

578

difference of onset (and offset) hydrocarbon generation among the samples in Group I

579

can be up to 15 °C, indicating a kinetically heterogeneous source rock (Fig. 13d). For

580

a geothermal gradient of 30 °C/km, this translates to 500 meters of difference in burial

581

depth for onset hydrocarbon generation for the same source rock unit in the extreme

582

case.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

583 584

585

Figure 13. Source rock kerogen conversion for the non-contaminated samples in this study; a) and

586

b) observed from laboratory pyrolysis at a heating rate of 25 °C/min; and c) and d) inferred for

587

geological conditions at a heating rate of 3 °C/million years. Kinetics estimation was made by

588

optimization of both frequency factor A and activation energy Ea.

589

Discussion

590

Implications for Facies Change and Maturity

591

In a study of source rock thermal maturity of the Bakken Formation in North

592

Dakota, Nordeng42 demonstrated that, a) samples from the same source rock with

593

similar thermal maturity in a single well follow a linear relationship between mean E

594

and log(A); b) samples at different thermal maturities fall into different trends; and c)

ACS Paragon Plus Environment

Page 36 of 46

Page 37 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

595

all the trends are roughly parallel showing a systematic shift from low to high

596

maturities (Fig. 14). Thermal maturation results in kinetic parameters shifting to high

597

activation energies; while the observed spread of samples from the same well along a

598

trend-line may indicate a combined effect from laboratory random errors and the

599

source rock heterogeneity due to seasonal changes, terrestrial inputs, water condition

600

and others.

601

by random laboratory error has been discussed by many authors

602

rule-of-thumb for this effect on the kinetic parameter estimation is that an error of 1

603

kcal/mol in E is compensated by a factor of 2 error in A and that a random 2 °C

604

temperature measurement error may cause a spread of 4 kcal/mol in E

605

geochemistry laboratory of Geological Survey of Canada, a finely powdered (< 212

606

µm) Cretaceous Second White Speckled shale sample with a 5.1% TOC and 441 oC

607

Tmax has been used as a standard rock sample for the purpose of QA/AC of the

608

Rock-Eval procedure. The statistics of 11542 standard samples analyzed from 1998 –

609

2016 show a mean of 441 °C and a standard deviation of 1.0 °C in Tmax. For a

610

normally distributed measurement errors with a standard deviation of 1 °C in Tmax,

611

the 99% confidence interval is 4 °C in temperature. This could translate to 8 kcal/mol

612

in E, meaning that there is a 99% possibility that the spread of data point along the

613

log(A)−E straight-line due to laboratory errors would not be greater than 8 kcal/mol in

614

E.

The spread of data points along a log(A)−E compensation line as caused 34, 43, 44

. A simple

43

. At the

615

In the Luo 69 well, the sampling depth interval spans about 130 metres, and

616

thermal maturity should be similar or slightly increase with increasing depth.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

617

However, Tmax values of the samples do not follow a systematic order (Fig. 3c and 4),

618

indicating different generation kinetics rather than exogenous input affecting the

619

thermal indicator of the samples as suggested by results from the extracted replicates

620

(Figs. 9 and 10). Each sub-unit shows slight difference in lithology and thus small

621

change in kinetics are expected. However, samples from the same stratigraphic

622

sub-units A and D+E do not align in straight lines in a semi-logarithmic cross-plot of

623

Ea and log10(A) (Figs 15a and c). Additional factor, such as maturity overprints due

624

to source rock heterogeneity, is apparent.

625 626

Figure 14. Comparison of sample thermal maturity expressed as kerogen decomposition path (a)

627

versus kerogen generation kinetics (b). Source of data is Table 1 of Nordeng.42

ACS Paragon Plus Environment

Page 38 of 46

Page 39 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

628 629

Figure 15. Estimated kinetic parameters for samples in different sub-stratigraphic units showing

630

variation in kinetic properties as well as maturity overprints, a) samples in sub-unit A; b) in

631

sub-unit C; c) in sub-units D+E; and d) comparison of all samples.

632

Lamination and Source Rock Heterogeneity

633

Lamination is a common feature and unique to shale. Two alternative functional

634

characters of self-sourcing and self-retaining make heterogeneity special in shale oil

635

reservoir. The two mutual competing (exclusive) capacities rely on alternations of

636

vertically stacked sedimentary rocks of different characteristics within a single source

637

rock unit. Organic-rich lamina are the primary sources of hydrocarbon generation in a

638

shale. In contrast, in a coarse grained segment embedded in the laminar, the rock has

639

more storage capacity with reduced generation capacity and serves as ideal storage for

640

expelled hydrocarbons from the organic-rich lamina. Oil generated from the lamina

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 46

641

likely migrate to and are stored in coarser, more porous, but organic-leaner thin beds

642

and lenses within the shale. The “contamination” by self-generated in the same shale

643

unit affects pyrolysis experiments, resulting in false high S1 peak with a retarded S1b

644

hump in default S2 temperature regime on their hydrocarbon pyrogram and

645

suppressed Tmax. 13, 445

646

On the other hand, the coarse-grained thin layer or lens is deposited in an

647

alternative environment or the same environment with more terrestrial input, leading

648

to a mixture of indigenous organic matter with more terrestrial organic matter

649

Alternation in lithology results in changes in not only organic matter mixing and bulk

650

kerogen components, but also porosity and permeability, which can lead to permeable

651

zones within source rock unit that are likely connected with more mature source rocks.

652

Therefore, an additional cause for high S1 and retarded S1b humps could be

653

exogenous oils generated from thermally mature source rock by up-dipping migration.

654

Effects of Free Hydrocarbons versus Varying Kinetic Property

655

Tmax value has been regarded as thermal indicator in source rock study, particularly

656

in marine source rock evaluation. However, caution should be exercised when applied

657

to source rock with mixed kerogen types. When comparing thermal maturity of a

658

heterogeneous source rock, absolute value of Tmax, Ro% or temperature (depth) may

659

not be appropriate and the rate of kerogen conversion to hydrocarbon could be a better

660

term to use. This is because thermally less stable source rock may start to generate,

661

while thermally stable source rock remains unchanged under the same thermal stress

ACS Paragon Plus Environment

46

.

Page 41 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

662

level. Thermally less table samples convert more kerogen to hydrocarbons, thus have

663

higher thermal maturity level than those of thermally stable samples.

664

For example, in Luo 69 well, the samples in Group III from deepest part of the well

665

show the lowest Tmax value among the samples from all three groups. Although the

666

lowest in Tmax, samples in Group III show the highest PI and S1/TOC values. The

667

low generation potential (HI600 mg HC/g TOC) and low OI suggest

710

that the samples in Group I contain primarily Type I kerogen. The narrow activation

711

energy spectrum indicate unitary kerogen (Figs 11) and more thermally stable source

712

rock (Fig. 12). Although Tmax values are the highest, source thermal maturity

713

remains lowest in Group I. The kinetic model (Fig. 12) represented by sample # 97 97

714

and 85 indicate an immature nature for the samples in Group I, consistent with the

715

Rock-Eval data interpretation.

716

Conclusions

717

Data analyses suggest considerable variations in geochemical and petrophysical

718

characteristics in the lower part of the Es3 Member of the Shahejie Formation in Luo

719

69 well, which led to a subdivision of the unit into three stratigraphic segments for

720

characterization of the heterogeneity. Source rock samples seem to have similar bulk

721

geochemical characteristics and petrophysical properties within each segment, but

722

differ considerably across segments. Source rock heterogeneity shows in not only

723

organic richness, type of organic matter and hydrocarbon generation potential, but

724

also in their kinetic properties.

725

Characterization of source rock heterogeneity in shale oil reservoir is challenging

726

because of the two mutually competing properties of self-sourcing and -retaining.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

727

Exogenous oil contamination from self-generated oil may blur source rock’s true

728

generation capacity and kinetic properties. Identification of the overprints from

729

contamination is an important step in characterization of the heterogeneity.

730

The composite generation kinetics study shows a temperature variation of about

731

25 °C from 10% to 90% transformation ratios in geological time scale for the

732

apparently “homogeneous” source rock in sub-unit C (Group I). The kinetic models

733

also indicates a more unitary and thermally resistant kerogen composition in samples

734

from Group I than those from the other two Groups.

735

Acknowledgements

736

The authors acknowledge the financial support of China Major Science 973

737

Project (No. 2014CB239101) and Sinopec Key Laboratory of Petroleum

738

Accumulation Mechanisms Open Funds, and Sinopec management for permission to

739

publish this work. This represents an output from Geoscience for New Energy Supply

740

Program of Natural Resources Canada and is partly supported by PERD funding.

741

References

742 743 744 745 746 747 748 749 750 751

(1) Maugeri, L. Oil and Gas Journal 2015, 3 (http://www.ogj.com/articles/uogr/print/volume-3/issue-1/comment-beware-of-break-even-and-margin al-cost-analyses.html). (2) Chen, Z.; Hannigan, P. Canadian Journal of Earth Sciences 2015, 53, 156-167. (3) Hill, R. J.; Zhang, E.; Katz, B. J.; Tang, Y.. AAPG bulletin 2007, 91, 501-521. (4) Passey, Q. R.; Bohacs, K.; Esch, W. L.; Klimentidis, R.; Sinha, S. From oil-prone source rock to gas-producing shale reservoir-geologic and petrophysical characterization of unconventional shale gas reservoirs:SPE 131350. 2010, (http://www.onepetro.org/mslib/servlet/onepetropreview?id=SPE-131350-MS) (5) Aplin, A. C.; Macquaker, J. H. AAPG Bulletin 2011, 95, 2031-2059.

ACS Paragon Plus Environment

Page 44 of 46

Page 45 of 46

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794

(6) Bohacs, K.; Passey, Q.; Rudnicki, M.; Esch, W.; Lazar, O. International Petroleum Technology Conference, paper IPTC 16676. 2013, (http://dx.doi.org/10.2523/16676-MS). (7) Lei, Y.; Luo, X.; Wang, X.; Zhang, L.; Jiang, C.; Yang, W.; Yu, Y.; Cheng, M.; Zhang, L. AAPG Bulletin 2015, 99, 661-687. (8) Chen, Z., Guo, Q., Jiang, C., Liu, X., Reyes, J. and Jia, Z. International Journal of Coal Geology, 2017, 182, 52-65. doi: 10.1016/j.coal.2017.08.017. (9) Glaser, K. S.; Miller, C. K.; Johnson, G. M.; Toelle, B.; Kleinberg, R. L.; Miller, P.; Pennington, W. D. Oilfield Review 2013, 25, 16-29. (10) Cudjoe, S.; Vinassa, M.; Gomes, J. H. B.; Barati, R. G. Journal of Natural Gas Science and Engineering 2016, 33, 1201-1218. (11) Chen, Z.; Jiang, C.; Lavoie, D.; Reyes, J. International Journal of Coal Geology 2016, 165, 290-302. (12) Chen, Z.; Lavoie, D.; Jiang, C.; Duchesne, M. J.; Malo, M., Geological Characteristics and Petroleum Resource Assessment of the Macasty Formation, Anticosti Island, Quebec, Canada, Geological Survey of Canada, Open File 8018, 2016, DOI: 10.4095/297865. (13) Li, M.; Chen, Z.; Cao, T.; Ma, X.; Liu, X.; Li, Z.; Jiang, Q.; Wu, S., Expelled oils and their impacts on Rock-Eval data interpretation and generation kinetics, Eocene Qianjiang Formation in Jianghan Basin, China, Part 1, Data analysis and identification of affected samples. Organic Geochemistry. (under review) (14) Li, M.; Li, Z.; Jiang, Q.; Cao, T. AAPG Search and Discovery Article #80342. Posted November 25, 2013, (http://www.searchanddiscovery.com/documents/2013/80342li/ndx_li.pdf). (15) Zhang, L.; Youshu, B.; Juyuan, L.; Zheng, L.; Rifang, Z.; ZHANG, J. Petroleum Exploration and Development 2014, 41, 703-711. (16) Jiang, Q.; Li, M.; Qian, M.; Li, Z.; Li, Z.; Huang, Z.; zhang, C.; Ma, Y. Petroleum Geology And Experiment 2016, 38, 842-849. (17) Wang, Y.; Wang, X.; Song, G.; Liu, H.; Zhu, D.; Zhu, D.; Ding, J.; Yang, W.; Yin, Y.; Zhang, S. Petroleum Exploration and Development 2016, 43, 759-768. (18) Bohacs, K.; Neal, J.; Carroll, A.; Reynolds, D. AAPG Search and Discovery Article #90914, accessed 2000, (http://www.searchanddiscovery.com/abstracts/html/2000/annual/abstracts/0062.htm). (19) Behar, F.; Lorant, F.; Budzinski, H.; Desavis, E. Energy and Fuels 2002, 16, 831-841. (20) Flynn, J. H. Thermochimica Acta 1997, 300, 83-92. (21) Burnham, A. K.; Braun, R. L. Energy and Fuels 1999, 13, 1-22. (22) Chen, Z.; Liu, X.; Guo, Q.; Jiang, C.; Mort, A. Fuel 2017, 194, 91-101. (23) Chen, Z.; Liu, X.; Jiang, C. Energy and Fuels 2017, 31, 1832-1841. (24)

Passey, Q.; Creaney, S.; Kulla, J.; Moretti, F.; Stroud, J. AAPG Bulletin 1990, 74, 1777-1794.

(25) Bowman, T. AAPG Search and Discovery Article #110128. Posted June 14, 2010, (http://www.searchanddiscovery.com/documents/2010/110128bowman/ndx_bowman.pdf). (26) Jarvie, D. M. AAPG Memoir 2012, 97, 69-87. (27) Jarvie, D. M. AAPG Memoir 2012, 97,89-119. (28) Mardia K.V., Kent J.T., Bibby J.M. Multivariate Analysis, 7th edn. 1989, Academic Press, New York, pp. 213–54. (29) Davis J. Statistics and Data Analysis in Geology, 2nd edn. 1986, John Wiley and Sons Inc., New York, pp. 527–45.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828

(30) Chen, Z.; Osadetz, K.; Gao, H.; and Hannigan, P. Bulletin of Canadian Petroleum Geology, 2001, v. 49, no. 3, p.367-375. (31) Chen, Z., Osadetz, K. G., Jiang, C. and Li, M., 2009, AAPG Bulletin, v.93, no. 6, p. 289-251. (32) Chen, Z. Grasby, S.E., Osadetz, K.; Geofluids, 2011, 11, p. 228–241 (33) Qian, M.; Jiang, Q.; Li, M.; Li, Z.; Liu, P.; Ma, Y.; Cao, T. Petroleum Geology and Experiments 2017, 39, 279-286. (34) Peters, K. E.; Burnham, A. K.; Walters, C. C. AAPG Bulletin 2015, 99, 591-616. (35) Burnham, A. K.

2014. Thermochimica Acta 2014, 597, 35–40.

(36) Chang, J., Qiu, N., Zhao, X., Xu, W.,Xu, Q,Jin, F., Han, C.,Ma, X.,Dong, X.,Liang, X., Chinese Journal Geophysics, 2016, 59(3): 1003-1016,doi: 10.6038/cjg20160322. (37) Liu, C. Li, K., Chen,Y. and Chen, J., 2016, Proceedings, 41st Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, February 22-24, 2016, SGP-TR-209 (38) Peters, K. E., Clifford C. Walters, and Mankiewicz, P.J., AAPG Bulletin, 2006, 90, 3, 387–403 (39) Chen, Z., C. Jiang, J. Chen, 2015. Comparison of Source Rock Kerogen Kinetics Using a Data-Driven Model and Based on Rock-Eval Pyrolysis Data. AAPG ICE, Melbourne Australia, Sept. 13-16, 2015. http://ice.aapg.org/2015/program/monday-oral. (40) Chen, Z. Orbermajer, M. Jiang, C. and Liu X.., Evaluating petroleum generation kinetics from Rock-Eval data , Example from Upper Ordovician Strata in southern Ontario. AAPG Bulletin. (under review) (41) Chen, Z.; Li, M.; Cao, T.; Ma, X.; Liu, X.; Li, Z.; Jiang, Q.; Wu, S., Impacts of exogenous hydrocarbons on Rock-Eval data interpretation and generation kinetics, Eocene Qianjiang Formation, Jianghan Basin, China, Part 2- Correction and kinetic models. Organic Geochemistry .(under review) (42) Nordeng, S. North Dakota Geological Survey 2013. (https://www.dmr.nd.gov/ndgs/documents/Publication_List/pdf/geoinv/GI_164.pdf). (43) Burnham, A. K., Thermochimica Acta 2014, 597, 35–40 (44) Stainforth, J. G., Marine and Petroleum Geology, 2009, 26, 552-572. (45) Snowdon, L. R. AAPG Bulletin 1995, 79, 1337-1348. (46) Chen, Z., Guo, Q., Jiang, C., Liu, X., Reyes, J. Mort, A. and Jia, Z., International Journal of Coal Geology, 2017, 182 52-65. (47) Raji, M., Gr€ocke, D.R., Greenwell, H.C., Gluyas, J.G. and Cornford C., Marine and Petroleum Geology, 2015, 67, 154-169. (48) Zhang, J., Jiang, Z, Jiang, X., Wang, S., Liang, C., and Wu M., Marine and Petroleum Geology, 2016, 71, 344-359.

ACS Paragon Plus Environment

Page 46 of 46