Thermal Maturity of the Marcellus Shale and the Line of Death : NE Pennsylvania

Transcription

1 Thermal Maturity of the Marcellus Shale and the Line of Death : NE Pennsylvania Conducted by Ethan Shula Department of Science & Mathematics Cedarville University, Cedarville, Ohio eshula@cedarville.edu

2 Objective: The objective of the study was to identify characteristics of the Marcellus Shale along the Line of Death and determine the principal factors influencing production. Thermal maturity, porosity, clay mineralogy, structural geology, completion procedures play varying roles in the success or failure of unconventional operations. Resolving discrepancies in the location of the Line of Death has the potential for significant economic benefits for operators located in the area.

3 Study Area Study Area: Susquehanna Co., PA Wyoming Co., PA A total of four wells were included in the study. Chief Oil & Gas LLC : Noble Unit 1H Squier Unit B-2H Harvey Unit 6H Jerauld Unit 1H

4

5

6 Introduction: The Marcellus Shale has developed into the nations most productive unconventional shale play since its emergence in Drilling efforts have increased exponentially since 2007 and have been primarily located in NE and SW Pennsylvania 8,341 unconventional wells have been completed since 2007 The formation is found within the Appalachian Basin at depths between 4,000 and 9,000 ft. and ranges in thickness from 20 to 250 ft.

7 EIA Productivity Report Gas production Mmcf/day Region Apr-15 May-15 Bakken 1,549 1,528 Eagle Ford 7,532 7,487 Haynesville 7,118 7,162 Marcellus 16,706 16,716 Niobrara 4,680 4,630 Permian 6,437 6,441 Utica 1,972 2,007 Total 45,994 45,971 Data reported in EIA s April 2015 Drilling Productivity Report.

8 EIA Productivity Report

9 Figure 1. Map depicting the depth to the Marcellus Shale (Marcellus Center for Outreach and Research, 2010).

10 Figure 2. Isopach map of the Marcellus Shale.

11 Geologic Setting: Line of Death Line of Death (Line of Demarcation) Operators have approximated the location of the line which marks the southern most limit of drilling operations due to the lack of production. Operators point to the region s heightened thermal maturity as the key factor in poor production. High thermal maturation is believed to be the result of extreme thermal and tectonic conditions associated with the Pennsylvania anthracite region to the south. Lackawanna Synclinorium located directly to the southeast

12 Geologic Setting: Line of Death The Lackawanna Synclinorium is a 110 km long crescent-shaped anthracite basin located in portions of Columbia, Luzerne, Lackawanna, Wayne, and Susquehanna Counties. Anthracite basins threaten natural gas development because of the high thermal conditions necessary for their development and the corresponding degradation of the surrounding strata. Anthracite is the highest rated coal Over 87% carbon content Associated with increased metamorphism

13 Prepared by Bureau of Topographic and Geologic Survey. Third Edition, Revised, 2000; Third Printing, 2008.

14 Anthracite Deposition Four Major Anthracite Fields: Northern, Eastern Middle, Western Middle, and Southern Upper Mississippian/Pennsylvanian Strata Coal seems found in the Pottsville and Llewellyn Formations Dominant lithology is sandstone and conglomerate Pottsville includes up to 14 coal beds (many discontinuous) Llewellyn includes up to 40 mineable coal beds Pennsylvanian age strata are time equivalent to the Bituminous fields of Western Pennsylvania

15 Anthracite Deposition Three Major Orogenic Events: Taconic Orogeny ( Ma) Acadian Orogeny ( Ma) Allegheny Orogeny ( Ma) Allegheny Orogeny (Permian) responsible for the development of the anthracite fields Estimated coalification temperatures: o C Estimated geothermal gradients: o C/km Estimated burial depths: km

16 Anthracite Deposition Hydrothermal Alteration No significant igneous intrusions in close proximity to Lackawanna Synclynorium USGS Open Field Report # states: orogenic uplift on the eastern flank of the Appalachian basin during the Alleghanian orogeny created a gravity-driven flow system that flushed hot basinal fluids out of the core of the Alleghanian orogen toward the west and northwest Theory that the synclinorium were formed as a result of Upper Silurian Salina salt-collapse structures. Alleghanian fluid or foreland-directed salt migration Resistive sandstone ridges of the synclinal basins protected against erosion of Pottsville and Llewellyn Formations

17 Previous Work: Thermal maturity values of Devonian strata in the Appalachian Basin have been well documented. Vitrinite Reflectance (%R o ) and Conodont Alteration Indexing (CAI) %R o CAI Vitrinite reflectance is a measurement of the percent light reflectance off the vitrinite maceral at 500x magnification in oil immersion. Conodont alteration indexing utilizes the fossil remains of chordates included in the stratigraphy.

18 %R o vs CAI Figure 6. Table displaying the conodont alteration index and associated temperatures in comparison with %R values and % fixed carbon (McCarthy et al., 2011).

19 Figure 4. Correlation of Coal Rank, %Ro, TAI, etc. (Repetski et al., 2008).

20 %R o vs CAI %Ro vs CAI: Window %Ro CAI Oil Condensate Dry Gas Overmature > Study Area %Ro values: Study Area CAI values: (according to the USGS Professional Paper #1708 by Repetski et. al. 2008)

21 Figure 4. Devonian %R o(max) isograds of Pennsylvania and New York (Repetski et al., 2008).

22 Figure 5. Devonian CAI max isograds of Pennsylvania and New York (Repetski et al., 2008).

23 Clay Diagenesis The concept of shale diagenesis is a significant factor in determining the thermal maturity of the underlying strata. Quantitatively, the most important diagenetic clay reaction in shale is the progressive transformation of smectite into illite via mixed layer illite/smectite (I/S) This reaction is also commonly referred to as smectite diagenesis or smectite illitization, (Pollastro, 1993). Clay mineralogy indicative of thermal maturation

24 Illite vs %Ro Illite content of Illite-Smectite mixed layer clay minerals*: Ro = 0.5% approximately corresponds to 25% Illite Ro = % approximately corresponds to 25-50% Illite Ro 1.5% approximately corresponds to 75% Illite *From Jiang, S., 2012, Clay Minerals from the Perspective of Oil and Gas Exploration: Clay Minerals in Nature Their Characterization, Modification, and Application, Chapter 2. Figure 6. Correlation of Reynolds and Howard s I/S notation and the Hoffman and Hower model (Pollastro, 1993)

25 Methods: Each well s completion report was reviewed in order to determine the target zones and lithologic sequences encountered. Completion reports were used to determine the perfing records and to help identify specific intervals for sampling. Sample Summary: Sample # Noble Unit 1H 8 Squier Unit B-2H 8 Harvey Unit 6H 10 Jerauld Unit 1H 6

26 Well Summaries: API # Well County TD TVD Spud Date Noble 1H Susquehanna 10,385 ft 7,276 ft 12/16/2011 Formation TD TVD Upper Marcellus 6,886 ft 6,830 ft Purcell 7,086 ft 7,009 ft Hot Marcellus 7,254 ft 7,138 ft Sample Intervals Sample # Interval (TD) Weight (g) #1 9,490-9, #2 9,580-9, #3 9,670-9, #4 9,850-9, #5 9,970-10, #6 10,060-10, #7 10,150-10, #8 10,270-10,

27 API # Well County TD TVD Spud Date Squier B-2H Susquehanna 11,992 ft 7,408 ft 11/3/2011 Formation TD TVD Upper Marcellus 7,260 ft 7,215 ft Purcell 7,580 ft 7,392 ft Sample Intervals Sample # Interval (TD) Weight (g) #1 8,080-8, #2 8,590-8, #3 8,950-8, #4 9,550-9, #5 9,970-10, #6 10,480-10, #7 10,990-11, #8 11,800-11,

28 API # Well County TD TVD Spud Date Harvey 6H Wyoming 14,640 ft 7,896 ft 11/2/2013 Formation TD TVD Upper Marcellus 7,482 ft 7,395 ft Purcell 7,754 ft 7,608 ft Lower Marcellus 7,850 ft 7,661 ft Hot Marcellus 7,910 ft 7,696 ft Sample Intervals Sample # Interval (TD) Weight (g) #1 8,470-8, #2 9,100-9, #3 9,500-9, #4 10,200-10, #5 10,700-10, #6 11,200-11, #7 11,800-11, #8 12,300-12, #9 12,800-12, #10 13,300-13,

29 API # Well County TD TVD Spud Date Jerauld 1H Susquehanna 11,343 ft N/A 7/12/2011 Formation TD TVD Upper Marcellus 7,220 ft N/A Sample Intervals Sample # Interval (TD) Weight (g) #1 10,375-10, #2 10,475-10, #3 10,675-10, #4 10,775-10, #5 10,925-10, #6 11,025-11,

30 Production Trends

31 Noble 1H Year Month MCF / month , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,265 Squier B-2H Year Month MCF / month , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,397 Harvey 6H Year Month MCF / month , , , , , ,435 *Closest proximity to the Line *Initially considered to be to the south of the line

32 Interpreted Zones of Production Susquehanna Co. Wyoming Co.

33 Time Mcf/Month 1 77,088 82,806 83,365 94, , , , , , , , , , , , , , ,440 91,969 90,243 89,793 91,177 89, ,734 75,695 75,695 73,253 75,695 73, ,694 EUR = 12.7 Bcf b = 1.7 Zero-Time Graph: Line Zone Note: EUR evaluation conducted using Halliberton Aries software

34 Time Mcf/Month 1 234, , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,508 EUR = 24.6 Bcf b = 1.7 Zero Time Graph: Core Zone Note: EUR evaluation conducted using Halliberton Aries software

35 Methods

36 X-Ray Diffraction Individual samples were separated, hand ground, packed into micro powder holders, and placed in a goniometer which scanned the samples from 2-60 degrees (2θ). Clay separates, homogenized out of suspension, were placed in a second goniometer and scanned from 2-35 degrees (2θ). Clay separates included only those particle 5μm The same samples were then removed and placed in a container containing glycol vapors. This procedure is performed in order to identify expansive clays (smectite and layered illite/smectite)

37 X-Ray Diffraction: Results The findings reported by Calgary Rock & Materials Services Inc. were similar between the four sampled wells. Lithologies were comprised of quartz, plagioclase (albite), calcite, dolomite, pyrite, kaolinite, illite, chlorite, and insignificant/minor quantities of layered illite/smectite. The largest Weight % and Vol. Fraction were comprised of Quartz and Carbonate material (calcite/dolomite).

38 Noble 1H Mineral Mean Wgt. % Carbonate 48.0 Quartz 29.0 Harvey 6H Mineral Mean Wgt. % Carbonate 57.4 Quartz 26.1 Squier B-2H Mineral Mean Wgt. % Carbonate 71.7 Quartz 15.0 Jerald 1-H Mineral Mean Wgt. % Carbonate 33.0 Quartz 32.6 Phyllosilicate composition almost entirely of Illite clays Kaolinite and Chlorite also somewhat persistent

39 Bulk Sample: Air-dry Run

40 Clay Separates: Composite-run

41 Interpretation Peak shifting occurs in samples that contain expansive clays. Two curves are produced, the first (blue) represents the air-dried run, while the second (green) represents the glycolated run. The glycolated samples expand the interlayer distance of the phyllosilicate sheet structure due to the substitution of larger radii glycol molecules for the original H 2 O molecules. Where expansive clay mineralogy is present, the peaks of the blue and green curves will experience shifting.

42 Spectrum illustrating minor peak shifting occurring at 9θ of Noble Unit 1H #5 sample.

43 Spectrum illustrating major peak shifting signifying the presence of expansive clay mineralogy.

44 Photomicroscopy Photomicrograph analysis was conducted using both a Motic 2300 camera mounted on a Motic polarizing microscope and a Nikon Eclipse 50i Polarizing Microscope with NIS Element visual software. Sequential imagery was taken for individual samples and stacked using Helicon Focus Pro imaging software. Images were created from samples located in the treated lateral.

45 Photomicroscopy: Noble 1H (4x)

46 Noble 1H (10x)

47 Squier B-2H (4x)

48 Squier B-2H (10x)

49 Harvey 6H (4x)

50 Jerauld 1H (4x)

51 Jerauld 1H (10x)

52 Photomicrograph depiction of cuttings at 4x magnification. A= Squier Unit B-2H; B= Noble 1H; C= Jerauld Unit 1H; D= Harvey Unit 6H.

53 Conclusion: High carbonate and silicate mineral composition dominate the Lower/Hot Marcellus. Only two samples from the Noble Unit 1H included ordered illite/smectite content. XRD and photomicrograph analyses did not identify any significant differences in the composition of the samples. High % Illite in the clay species separate points to a uniformly mature shale devoid of expansive clays.

54 Future Consideration: Further evaluation should include: XRD, TOC%, RockEval, SEM Backscatter Imaging, etc. Addition of such data would reveal the type of hydrocarbons present, the quality of the organic matter, the internal pore structure controlling production, and the true thermal maturity of the unit. Significant funding and lab testing required to complete such an undertaking at the undergraduate level.

55 Acknowledgments: I would like to thank Mr. Terry Ward and Chief Oil & Gas Inc. for providing samples, drilling reports, and all geophysical logs included in the study. I would like to thank Mr. Ray Strom and Calgary Rock and Materials Services Inc. for advising and performing XRD analyses. Additionally, I would like to thank Mr. Keith Mangini and Huntley & Huntley Inc.. Without their help, this project would not have been possible. I would also like to thank Professor Rice and Dr. Whitmore for providing countless hours of work on the development, advisement and completion of the project.

56 References: Boswell, Z., 2010, A Study of Natural Gas Extraction in Marcellus: Masters Thesis, Massachusetts Institute of Technology. Cardott, B., 2012, Introduction to Vitrinite Reflectance as a Thermal Maturity Indicator: Oklahoma Geological Society. AAPG Search and Discovery Article. Curtis, M., et al., 2010, Microstructural Observations in Gas Shales: AAPG Hedberg Conference. AAPG Search and Discovery Article. Curtis, M. et al., 2012, Development of Organic Porosity in the Woodford Shale with Increasing Thermal Maturity: International Journal of Coal Geology. v pp Engelder, T., and Lash, G., 2008, Marcellus Shale Play s Vast Resource Potential Creating Stir in Appalachia: The American Oil and Gas Reporter. Epstein, A., et al., 1977, Conodont Color Alteration-an Index to Organic Metamorphism: USGS Professional Paper: 995. Harrison, M., et al., 2004, The Lackawanna synclinorium, Pennsylvania: A salt-collapse structure, partially modified by thin-skinned folding: GSA Bulletin. v no. 11/12. pp

57 Jiang, S., 2012, Clay Minerals from the Perspective of Oil and Gas Exploration: Clay Minerals in Nature Their Characterization, Modification, and Application, Chapter 2. InTech. Levine, J., and Eggleston, J., 1992, The Anthracite Basins of Eastern Pennsylvania: 1992 Joint Meeting of the International Committee for Coal and Organic Petrology & The Society for Organic Petrology. U.S. G.S Open-Rile Report # Marcellus Region, Drilling Productivity Report, September 2014: U.S. Energy Information Administration. McCarthy, K., et al., 2011, Basic Petroleum Geochemistry for Source Rock Evaluation: Oilfield Review. v. 23 no. 2. pp Pollastro, R., 1993, Considerations and Applications of the Illite/Smectite Geothermometer in Hydrocarbon-Bearing Rocks of the Miocene to Mississippian Age: Clays and Clay Minerals. v. 41. no. 2. pp Repetski, J., et al., 2008, Thermal Maturity Patterns (CAI% and %Ro) in Upper Ordovician and Devonian Rocks of the Appalachian Basin: A Major Revision of the USGS Map I-917-E Using New Subsurface Collections: U.S. Geological Survey, Reston, Virginia.

58 Ruppert, L., et al., 2014, Thermal maturity patterns in Pennsylvanian coal-bearing rocks in Alabama, Tennessee, Kentucky, Virginia, West Virginia, Ohio, Maryland, and Pennsylvania, chap. F.2 of Coal and petroleum resources in the Appalachian basin; Distribution, geologic framework, and geochemical character: USGS Professional Paper 1708, 13 p. Ruppert, L., et al., 2010, Geologic Controls on Thermal Maturity Patterns in Pennsylvania coal- bearing rocks in the Appalachian Basin: International Journal of Coal Geology. v. 81. pp Ryder, R., et al., 2012, Geologic Cross Section C-C Through the Appalachian Basin from Erie County, North-Central Ohio, to the Valley and Ridge Province Bedford County, South-Central Pennsylvania: USGS. Schieber, J., 2010, Common Themes in the Formation and Preservation of Intrinsic Porosity in Shales and Mudstones Illustrated with Examples Across the Phanerozoic: Society of Petroleum Engineers. Schieber, J., 2013, SEM Observations on Ion Milled Samples of Devonian Black Shales from Indiana and New York: The Petrographic Context of Multiple Pore Types: AAPG Memoir 102, pp Stephens, M., et al., 2009, Laboratory Methods to Assess Shale Reactivity with Drilling Fluids: American Association of Drilling Engineers National Technical Conference & Exhibition. Ver Straeten, C., 2009, The Classic Devonian of the Catskill Front: A Foreland Basin Record of Acadian Orogenesis: NYSGA 2009 Trip 7. Zagorski, W., et al., 2011, An Overview of Some Key Factors Controlling Well Productivity in Core Areas of the Appalachian Basin Marcellus Shale Play: Range Resources Corporation.