Status Report on Petroleum System Analysis Study of the Doig Formation, Western Sedimentary Basin, Northeastern

P.L. Silva, The University of British Columbia, Vancouver, BC, [email protected] R.M. Bustin, The University of British Columbia, Vancouver, BC

Silva, P.L. and Bustin, R.M. (2017): Status report on petroleum system analysis study of the Triassic Doig Formation, Western Canada Sedimentary Basin, northeastern British Columbia; in Geoscience BC Summary of Activities 2016, Geoscience BC, Report 2017-1, p. 25– 28.

Introduction • Basin modelling: modelling of the Doig petroleum sys- tem on a basin scale, incorporating structural, litholog- The Triassic section of the Western Canada Sedimentary ical, geochemical and burial-history data, in order to de- Basin (WCSB) is the richest interval in terms of volume of termine the timing of thermogenic hydrocarbon oil per volume of rock in the basin (Marshall et al., 1987). generation and retention An important part of the succession is the Lower and Mid- dle Triassic Doig Formation, which historically has been • Reservoir properties: characterizing the reservoir po- known for limited production from its conventional reser- tential of the Doig Formation in terms of storage capac- voirs in British Columbia (BC) and . More recently, ity, producibility and response to hydraulic stimulation the Doig has been recognized as an important unconven- tional reservoir for gas and natural-gas liquids (NGL), with Geological Background estimates of total gas-in-place ranging from 1.1 to 5.6 tril- lion m3 (Walsh et al., 2006). Little is currently known, how- The Doig Formation consists of mudstone, siltstone and ever, about the unconventional portion of the Doig succes- subordinate , bioclastic packstone and grain- sion: the spatial and stratigraphic distribution of litho- and stone, deposited under marine conditions in environments reservoir facies is poorly constrained and measurements of ranging from shoreface through offshore (Evoy and Mos- petrophysical properties are sparse and limited in scope. low, 1995). It has long been recognized that the Doig can be The distribution of unconventional reservoir properties broadly subdivided into a lower, more organic-rich and within the Doig succession and the geological controls on phosphatic zone, known as the Doig Phosphate Zone these properties and their distribution are therefore un- (DPZ); and an upper, relatively organic-lean siltstone inter- known. Furthermore, there is an uncertainty in the distribu- val. A more detailed subdivision has been proposed by tion of hydrocarbon phase and composition, due to varia- Chalmers and Bustin (2012), who recognized three units: a tions in thermal maturity, depth of burial and organic- basal Doig A, composed mainly of interbedded dark argil- matter type. The purpose of the ongoing research project laceous and calcareous siltstones, corresponding to the reported herein is to determine the unconventional hydro- DPZ and distinguishable in well logs by a gamma-ray sig- carbon potential of the Doig Formation through a petro- nature; an intermediate section, named Doig B, corre- leum system analysis (PSA). The principal goal of the study sponding to medium to dark grey argillaceous siltstone is to delineate the distribution of producible liquids from containing localized sand beds; and an upper Doig C, tight reservoirs in northeastern BC. composed of siltstone and argillaceous fine sandstone. The PSAapproach used in this study consists of three inter- Ongoing Work related elements (Figure 1): To date, 470 m of core from fifteen wells in BC and 190 m • Source rock evaluation: thorough quantification and from seven wells in Alberta (Table 1) have been logged for mapping of the source-rock properties across the basin, lithology, sedimentary structures, bioturbation, diagenetic such as thermal-maturity levels and organic-matter features and structural features. These wells comprise three abundance and type strike and six dip cross-sections (Figure 2), covering the en- tire lateral extent of the Doig Formation in BC and Alberta. The primary selection criterion was availability and length of core in the Doig, in an attempt to adequately represent its Keywords: British Columbia, gas, shale oil, unconventional reservoir, petroleum system analysis stratigraphic and spatial variability. Another factor taken into consideration when selecting wells was the availability This publication is also available, free of charge, as colour digital files in Adobe Acrobat® PDF format from the Geoscience BC of a complete well-log suite over the Doig interval. The website: http://www.geosciencebc.com/s/DataReleases.asp. minimum requirement for this study is a quad-combo

Geoscience BC Report 2017-1 25 Figure 1. Elements of the petroleum system analysis approach to reservoir evaluation and their respective inputs as implemented in this study.

(gamma ray, density/neutron, resistivity and compressional Analytical work in progress includes X-ray diffraction for slowness), and preference was given to wells with shear quantitative mineralogy, Rock-Eval pyrolysis for organic slowness, nuclear magnetic resonance, elemental capture geochemical properties, helium pycnometry and mercury spectroscopy and spectral gamma ray. These logs, although intrusion for porosity and pore size distribution, pressure not available for all wells in the study, will aid in the devel- pulse-decay permeameter measurements of permeability, opment of neural-network models and in the extrapolation and triaxial cell tests for geomechanical properties. These of petrophysical and organic geochemical properties properties will be extrapolated to non-cored intervals and measured on core and on cuttings to non-cored sections and wells using quad-combo well logs through multiresolution wells. graph-based clustering and neural-network techniques.

Table 1. Wells from which core of the Doig Formation was logged. See Figure 2 for location.

26 Geoscience BC Summary of Activities 2016 Figure 2. Location of the wells used in this study overlain on the isopach map of the Doig Formation–Halfway For- mation interval.

Thermal maturity from pyrolysis will be used in a basin- ity spanning four orders of magnitude, despite a narrow modelling platform to constrain the maturity of the Doig range in grain size (Chalmers and Bustin, 2012; Chalmers Formation across the basin. The maps of kerogen type and et al., 2012). It is further anticipated that the degree of bio- total organic carbon, also derived from pyrolysis, will con- turbation, which changes the fabric of the sediments, may strain the spatial distribution and timing of thermal also control porosity and permeability (Pemberton and Gin- hydrocarbon generation. gras, 2005; Baniak et al., 2015).

Preliminary core logging reveals a wide variety of lithofac- ies and complex stacking patterns that confer a high degree There are currently no publicly available regional studies of heterogeneity on the succession. Preliminary work has with a systematic petrophysical characterization of the lith- identified ten lithofacies, defined based on colour, lithol- ofacies found in the Doig Formation. This study aims to fill ogy and fabric. Diagenetic and structural features add an- this gap and will serve as the foundation for a full petro- other layer of heterogeneity within individual lithofacies. physical and organic-geochemical characterization of These facies have distinctly different reservoir properties, these facies and their spatial and stratigraphic distribution such as porosity, distribution of pore-throat size, perme- across the WCSB. Ultimately, the lithofacies will be ability, elastic moduli and organic content. It has been grouped into facies associations to enable work at the log shown that petrophysical properties vary significantly be- scale, and the facies associations with their petrophysical tween the different Doig lithofacies, with median pore- properties keyed to cored wells will be extrapolated throat size spanning one order of magnitude and permeabil- through logs to adjacent wells.

Geoscience BC Report 2017-1 27 Acknowledgments Evoy, R.W. and Moslow, T.F. (1995): Lithofacies associations and depositional environments in the Doig For- The authors acknowledge financial support from Geo- mation, Buick Creek Field, northeastern British Columbia; science BC, Trican Geological Solutions Ltd., EnCana Cor- Bulletin of Canadian Petroleum Geology, v. 43, no. 4, poration, Husky Energy Inc., Chevron Canada Limited and p. 461–475. Canadian Natural Resources Limited. Additionally, the Marshall, G.M., Noble, I.A. and Tang, R.C.W. (1987): Triassic/Ju- rassic fields; in Geophysical Atlas of Western Canadian Hy- thorough review by E. Letham, Ph.D. candidate at The Uni- drocarbon Pools, N.L. Anderson, L.V. Hills and D.A. versity of British Columbia, is greatly appreciated. Cederwall (ed.), Canadian Society of Exploration Geophys- icists–Canadian Society of Petroleum Geologists, p. 187– References 215, URL [No- Baniak, G.M., Gingras, M.K., Burns, B.A. and Pemberton, S.G. vember 2016]. (2015): Petrophysical characterization of bioturbated sand- stone reservoir facies in the Upper Ula Formation, Pemberton, S.G. and Gingras, M.K. (2005): Classification and Norwegian North Sea, Europe; Journal of Sedimentary Re- characterizations of biogenically enhanced permeability; search, v. 85, no. 1, p. 62–81. AAPG Bulletin, v. 89, no. 11, p. 1493–1517, URL [November 2016]. of Halfway-Doig-Montney hybrid gas shale–tight gas reser- voir, northeastern British Columbia; Marine and Petroleum Walsh, W., Adams, C., Kerr, B. and Korol, J. (2006): Regional Geology, v. 38, no. 1, p. 53–72. “shale gas” potential of the Triassic Doig and Montney for- mations, northeastern British Columbia; BC Ministry of En- Chalmers, G.R.L., Bustin, R.M. and Bustin, A.A.M. (2012): Geo- ergy and Mines, Oil and Gas Division, Resource Development logical controls on matrix permeability of the Doig- and Geoscience Branch, Petroleum Geology Open File Montney hybrid shale-gas–tight-gas reservoir, northeastern 2006-02, 19 p., URL [November 2016]. ofActivities2011/SoA2011_Chalmers.pdf> [November 2016].

28 Geoscience BC Summary of Activities 2016