US20120134749A1 - Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales - Google Patents
Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales Download PDFInfo
- Publication number
- US20120134749A1 US20120134749A1 US13/297,263 US201113297263A US2012134749A1 US 20120134749 A1 US20120134749 A1 US 20120134749A1 US 201113297263 A US201113297263 A US 201113297263A US 2012134749 A1 US2012134749 A1 US 2012134749A1
- Authority
- US
- United States
- Prior art keywords
- fractures
- fracture
- fluid flow
- shale
- hydrocarbon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 134
- 229910052756 noble gas Inorganic materials 0.000 title claims description 36
- 230000005012 migration Effects 0.000 title claims description 21
- 238000013508 migration Methods 0.000 title claims description 21
- 229930195733 hydrocarbon Natural products 0.000 title abstract description 61
- 239000004215 Carbon black (E152) Substances 0.000 title abstract description 38
- 235000015076 Shorea robusta Nutrition 0.000 title abstract description 18
- 244000166071 Shorea robusta Species 0.000 title abstract description 18
- 125000001183 hydrocarbyl group Chemical group 0.000 title 1
- 238000000034 method Methods 0.000 claims abstract description 29
- 239000011573 trace mineral Substances 0.000 claims description 18
- 235000013619 trace mineral Nutrition 0.000 claims description 18
- 150000002835 noble gases Chemical class 0.000 claims description 9
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 60
- 239000011435 rock Substances 0.000 abstract description 39
- 238000000605 extraction Methods 0.000 abstract description 16
- 230000001965 increasing effect Effects 0.000 abstract description 12
- 238000005553 drilling Methods 0.000 abstract description 11
- 235000009508 confectionery Nutrition 0.000 abstract description 10
- 239000003673 groundwater Substances 0.000 abstract description 10
- 238000011084 recovery Methods 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 4
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 230000004888 barrier function Effects 0.000 abstract 1
- 206010017076 Fracture Diseases 0.000 description 193
- 208000010392 Bone Fractures Diseases 0.000 description 104
- 210000003462 vein Anatomy 0.000 description 31
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 30
- 239000007789 gas Substances 0.000 description 24
- 238000004519 manufacturing process Methods 0.000 description 21
- 230000003993 interaction Effects 0.000 description 19
- 230000035699 permeability Effects 0.000 description 16
- 229910052500 inorganic mineral Inorganic materials 0.000 description 14
- 235000010755 mineral Nutrition 0.000 description 14
- 239000011707 mineral Substances 0.000 description 14
- 239000000203 mixture Substances 0.000 description 13
- 239000010453 quartz Substances 0.000 description 13
- 229910052761 rare earth metal Inorganic materials 0.000 description 13
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 13
- 230000015572 biosynthetic process Effects 0.000 description 12
- 238000009792 diffusion process Methods 0.000 description 12
- 238000005755 formation reaction Methods 0.000 description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 12
- 229910001868 water Inorganic materials 0.000 description 12
- 239000003345 natural gas Substances 0.000 description 11
- 150000003839 salts Chemical class 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 10
- 230000033558 biomineral tissue development Effects 0.000 description 10
- 229910021532 Calcite Inorganic materials 0.000 description 9
- 229910052734 helium Inorganic materials 0.000 description 9
- 239000011148 porous material Substances 0.000 description 9
- 239000013078 crystal Substances 0.000 description 8
- 229910052768 actinide Inorganic materials 0.000 description 7
- 150000001255 actinides Chemical class 0.000 description 7
- 229910052748 manganese Inorganic materials 0.000 description 7
- 239000000126 substance Substances 0.000 description 7
- 238000011065 in-situ storage Methods 0.000 description 6
- 238000013507 mapping Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000004904 shortening Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- 229910052723 transition metal Inorganic materials 0.000 description 6
- 150000003624 transition metals Chemical class 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000001307 helium Substances 0.000 description 5
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 5
- 229910052742 iron Inorganic materials 0.000 description 5
- 238000000095 laser ablation inductively coupled plasma mass spectrometry Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 229910052770 Uranium Inorganic materials 0.000 description 4
- 238000003776 cleavage reaction Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 238000000608 laser ablation Methods 0.000 description 4
- 230000037361 pathway Effects 0.000 description 4
- 230000007017 scission Effects 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 3
- 235000019738 Limestone Nutrition 0.000 description 3
- 208000002565 Open Fractures Diseases 0.000 description 3
- 229910052776 Thorium Inorganic materials 0.000 description 3
- 230000003466 anti-cipated effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 229910052785 arsenic Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000155 isotopic effect Effects 0.000 description 3
- 229910052754 neon Inorganic materials 0.000 description 3
- 238000003012 network analysis Methods 0.000 description 3
- 230000002441 reversible effect Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000000476 thermogenic effect Effects 0.000 description 3
- 239000011486 Onondaga limestone Substances 0.000 description 2
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 description 2
- 241001466077 Salina Species 0.000 description 2
- 230000000035 biogenic effect Effects 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 230000003197 catalytic effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000007872 degassing Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011049 filling Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000008719 thickening Effects 0.000 description 2
- 229910052719 titanium Inorganic materials 0.000 description 2
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 2
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 description 1
- 241000167854 Bourreria succulenta Species 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 229910052684 Cerium Inorganic materials 0.000 description 1
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 1
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 description 1
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- 208000006670 Multiple fractures Diseases 0.000 description 1
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 150000001218 Thorium Chemical class 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000002378 acidificating effect Effects 0.000 description 1
- 230000005262 alpha decay Effects 0.000 description 1
- LBDSXVIYZYSRII-IGMARMGPSA-N alpha-particle Chemical compound [4He+2] LBDSXVIYZYSRII-IGMARMGPSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 235000019693 cherries Nutrition 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 229910000037 hydrogen sulfide Inorganic materials 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010249 in-situ analysis Methods 0.000 description 1
- 238000001095 inductively coupled plasma mass spectrometry Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 229910052746 lanthanum Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000004452 microanalysis Methods 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- IJDNQMDRQITEOD-UHFFFAOYSA-N n-butane Chemical compound CCCC IJDNQMDRQITEOD-UHFFFAOYSA-N 0.000 description 1
- 238000005380 natural gas recovery Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000000399 optical microscopy Methods 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052700 potassium Inorganic materials 0.000 description 1
- 239000011591 potassium Substances 0.000 description 1
- 239000002244 precipitate Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000001294 propane Substances 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000007420 reactivation Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 229910052725 zinc Inorganic materials 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
- E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
- E21B47/11—Locating fluid leaks, intrusions or movements using tracers; using radioactivity
Definitions
- the following description relates generally to natural gas recovery, and more particularly to a method for identifying highly productive locations for hydrocarbon extraction from black shale.
- Hydraulic fracturing involves the injection and extraction of fluids and propping agents in the subsurface to stimulate fluid flow through natural fractures and increase fracture related permeability (e.g. increased fracture aperture, fracture size, and fracture network connectivity) thus enhancing hydrocarbon production.
- hydraulic fracturing is limited by: 1) inefficient resource recovery, 2) the potential for groundwater contamination from drilling fluids or mobilized hydrocarbons that can migrate through fractures and interact with groundwater, and 3) the current inability to develop accurate models for fracture fluid flow in the exceedingly complex fracture network present in black shale and other fractured lithologies. Therefore successful, economically viable, and environmentally safe application of these techniques requires a detailed understanding of fluid transport within the subsurface, specifically fluid flow within fracture networks.
- Black shales are of great interest as domestic hydrocarbon plays because of their high organic content and tight gas-retaining nature.
- market demand is increasing interest in shale hydrocarbon extraction, the present inventor is unaware of any comprehensive methodology capable of identifying highly productive areas known as “sweet spots” because of a paucity of information about geological fracture-related fluid flow.
- sweet spots Because of the immense costs associated with horizontal drilling and hydraulic fracturing, as well as the risk of environmental impact, hydrocarbon recovery from unconventional sources, such as black shale, is not economically favorable unless wells hit a sweet spot.
- the present invention By directly evaluating natural fluid migration in situ using conservative noble gas diffusion profiles, trace element proxies for geological fluid flow, and the characteristics of the fracture network, the presently disclosed predictive model for natural gas migration and determination of the location of sweet spots in shale has been developed.
- the present invention will directly contribute to a comprehensive strategy for hydrocarbon extraction, specifically in regions which lack deformed structures that lead to bed thickening shale and increased fracturing.
- the present invention allows for the evaluation of risk for aquifer contamination from hydraulic fluid and shale gas.
- the present invention comprises a combination of methods.
- Noble gas abundances and isotopic ratios (He, Ne) and trace element geochemistry (transition metals (Ti, Mn, Fe), rare earth elements (La—Lu), actinide chemistry (Th/U)) with advanced techniques in fracture network analysis are integrated, in order to: 1) determine the multi-stage fluid flow history through individual fractures; 2) quantify gas diffusion from relatively impermeable hydrocarbon host rock (i.e. black shale) to fracture sets; 3) develop a 3-D geospatial map of the regional fracture network to determine pathways of fluid migration; and 4) map chemical changes for the development of a regional hydrocarbon “sweet spot” map.
- the Marcellus shale in New York is located on the Appalachian plateau and is in the foreland of the Appalachian Fold Thrust Belt. It lies to the North and West of the Valley and Ridge province and is characterized by salt-detachment tectonics.
- the first order deformation structures of the Appalachian plateau are detachment folds, which have a basal decollement in the Silurian-aged Salina Salt.
- This salt layer acts as the glide plane for the Appalachian plateau detachment sheet and causes surface features characteristic of salt tectonics such as broad, gentle folds and abrupt changes in deformation style at the lateral and frontal termination of the salt.
- Salt also plays a role in faulting. Salt is thickened in the hinges of detachment anticlines and provides a weak zone that is easily faulted, resulting in blind reverse splays that cut up from the decollement in the salt and slice through this weak core.
- the stratigraphy of the Appalachian plateau varies, but in New York the Marcellus shale is the stratigraphically lowest subgroup of the siliciclastic Devonian Hamilton group.
- the Hamilton group is sourced from the Acadian orogeny and Rb-Sr dating of the lower Devonian black shales in the Hamilton groups puts their age at 384 ⁇ 9 to 377 ⁇ 11 Ma.
- the Hamilton group is sandwiched between the overlying late Devonian Catskill deltaic sequence and the underlying clean carbonates of the Devonian Onondaga limestone.
- the Marcellus subgroup is comprised of three formations: 1) the Union Springs-the lowest stratigraphic black shale unit; 2) the Oatka Creek Formation-the highest stratigraphic black shale; and 3) its lateral stratigraphic equivalent, the Mount Marion Formation.
- the Union Springs formation is made of black shales and dark grey limestones, and is separated from the black shales of the Oatka Creek and Mount Marion formations by the Cherry Valley limestone.
- the Appalachian plateau was deformed during the Pennsylvanian-Permian Alleghany orogeny. During this event, deformation progressed from the hinterland to the foreland along the basal salt layer.
- the Appalachian Plateau detachment sheet progressed to the northwest during this orogeny and did not interact with the underlying basement. Deformation within the detachment sheet varies with stratigraphy; in the lowermost strata shortening was accommodated by low-angle thrust faulting, while at higher levels shortening was first taken up by layer-parallel shortening and then accommodated by broad folding. In all, Alleghanian deformation occurred during two progressive stages: layer parallel shortening occurred first and was followed by a period of detachment folding and reverse faulting.
- the first stage of deformation, layer parallel shortening, is expressed in surficial geology by the deformed fossils and solution cleavage. Strains of up to 20% are observed in fossils and both solution cleavage and fossil shortening indicates a strain ellipsoid that has its short axis perpendicular to the regional structural trend. Solution cleavage maintains its bed-perpendicular orientation around folds and confirms that cleavage formed before folding.
- detachment folds formed by buckling above the Salina salt. These folds are characterized by comparatively tight anticlines cored by salt and broad synclines (FIG. 3). The cores of folds are significantly tighter in the salt horizon, but at higher stratigraphic levels folds are open with very gently dipping limbs ( ⁇ 5 degree limb dips). Folds are slightly asymmetric with steeper southeastern limbs. Thrust splays cut up from the decollement through the weak anticlinal fold-cores but do not make it to the surface and both antithetic and synthetic faults are common.
- Fracturing is a pervasive and complex feature of the Appalachian plateau that developed during successive phases of the Alleghanian orogeny and is associated with the first order structures. Fractures can be grouped into sets by their orientation, but the relationships between these sets are not completely understood. N-S striking fractures are interpreted as extension fractures due to early E-W extension in the forebulge of the Alleghanian orogeny.
- Cross-fold fractures (strikes ranging from 012° to 327° are difficult to interpret; explanations vary from fracturing due to multiple phases of the Alleghanian orogeny, to fracturing due to the stress field before the Alleghanian orogeny or tectonic unloading after the orogeny, or a combination of these mechanisms that invokes reactivation of fractures.
- ENE-striking fractures( ⁇ 071°) are interpreted as neotectonic and due to overpressure caused by hydrocarbon generation.
- Hydrocarbons can be created through biogenic and thermogenic means, with biogenic processes being significant at shallow depths and thermogenic production dominating at deeper levels. Significant hydrocarbon generation is usually attributed to a thermogenic process; hydrocarbons are formed at depth from the thermal degradation of kerogen. As rock is buried, temperature and pressure increase and the structure of kerogen becomes unstable. Kerogen progressively adjusts to this increasing temperature and pressure by eliminating functional groups and the linkages between nuclei, thus generating a wide range of compounds including hydrocarbons, CO 2 , Water, and hydrogen sulfide. Additionally, natural gas comprised of methane, ethane, propane, and n-butane (C1, C2, C3, and C4, respectively) can be generated through the mechanism of transition-metal catalysis.
- the present invention can develop reservoir-, field-, and regional-scale interpretations of fluid flow without the need for extrapolation; this makes the present predictive technology especially useful to local inhabitants, state governments, and hydrocarbon extractions corporations when implemented in the early stages of exploration.
- the Marcellus shale offers logistical advantages that make it an excellent case study. For example, it is exposed in many quarries across New York State based on its stratigraphic position above the quarried Onondaga limestone. This coincidence allows study of the three-dimensional relationships of fractures with great accuracy, and sampling of fresh, unaltered, outcrops that are revealed through quarrying activities. These outcrops are not weathered and their geochemistry is preserved, making them a useful analogue for more deeply buried rocks.
- Some economically important geological advantages to studying the Marcellus shale in New York are the type of deformation, the amount of deformation, the fracture pattern, and the regional fluid flow that this area experienced. Compared with the intensely deformed sections of the Marcellus shale in the Valley and Ridge province of Pennsylvania, the Marcellus shale in New York is relatively undeformed and fracture patterns have not been overprinted by larger structures. This lack of deformation in Appalachian plateau region of the Marcellus prevented the thickening and increased fracturing of shale beds, making it very difficult to determine the best location for extraction wells and predict the location of highly fractured “sweet spots”. This necessitates the understanding of fracture related flow in the plateau section of the Marcellus shale for any drilling program.
- Fractures are surfaces in rock along which mechanical failure has occurred and the rock has lost cohesion. They can form in tension (mode 1) or shear (mode 2 or 3) and often form sets of similarly oriented members. Fractures that accommodate some degree of slip along their surfaces are faults, while fractures that have no observable slip are joints. A grouping of fractures with sub-parallel orientations is a fracture set, while all of the fractures regardless of orientation form the fracture network.
- Hydraulic conductivity is a function of a rock's bulk porosity and permeability and describes the ease with which a fluid is transported through pore spaces or fractures.
- the extractable volume of hydrocarbons is directly proportional to the system's hydraulic conductivity (hereafter: K), which is a factor of 1) the hydraulic conductivity of the country rock (K CR ) and 2) the fracture network (hereafter: K FN ).
- K hydraulic conductivity
- K FN the hydraulic conductivity of the fracture network
- Size, location, termination style, aperture, planarity and roughness are key characteristics to determine flow within individual fractures.
- the size of a fracture refers to the three-dimensional surface area of the fracture, while aperture is the openness of fracture planes. Planarity (a measure of a fracture's deviation from a plane) and roughness (planar tortuosity) are important factors in determining permeability.
- the termination style defines the geometric orientation of the end of a fracture. There are four types of terminations including: T—a perpendicular intersection between fractures; J—an intersection in which one fractures curves into the other; I—a fracture that ends at its tip line without intersection and X—cross-cutting fractures.
- Combinations of these termination geometries and the location of individual fractures define the 3-D geometry of a fracture network. All of the above parameters influence hydraulic conductivity (i.e. the amount of hydrocarbons that can migrate through the fracture) and each other. For example, permeability and porosity concomitantly increase with increasing country rock grain size and fracture size, and the statistical probability of fracture intersection (connectivity) also increases with larger fracture size. Fracture connectivity, porosity, and permeability all affect hydraulic conductivity implicating a complex relationship between K FN and fracture properties.
- Fractures are grouped by their geometry into fracture sets or groups of sub-parallel oriented fractures. Orientation and spacing of fracture sets in a network are characteristics that affect fluid flow in different ways. For dense and homogeneous fracture networks fluid flow can be treated as flow through a porous medium, while in sparsely fractured areas a few large fractures may dominate flow. Fracture network hydraulic properties depend on fracture intensity (surface area of fractures per unit volume), connectivity (number of fracture intersections per unit volume), hierarchy, and chronology.
- Modeling the migration of hydrocarbons in fractured black shales is exceedingly complex due in part to the complex nature of hydraulic conductivity in a fractured medium, but also to the many dynamic processes of the earth. For example, dynamic changes in parameters such as regional stress field, in fracture mineralization, fluid pressure, climatic changes (wetness/dryness), fluid gradient, anthropogenic water use, and tectonic processes reduce the accuracy of model inputs significantly and retard the understanding of fluid transport through fractured media.
- fractures must be open in order to accommodate fluid flow. How fractures remain open (aperture>0) and the relative importance of mechanical and diagenetic characteristics in keeping fractures open is still contentious. Some authors argue for the role of the in situ stress field and suggest that only fractures oriented parallel to the maximum compressive stress will stay open and accommodate fluid flow. Fractures may never completely close if there is a sufficient hydraulic gradient, even though permeability decreases significantly as stress normal to the fracture increases. In addition, some component of shearing can keep fractures open, causing asperities on opposite faces of the fracture to ride up over one another and prop open the fracture.
- hydrocarbon-rich Barnett Shale of Northern Texas has fault induced fractures, but drill cores show pervasive calcite veining which correlates with low hydrocarbon production in heavily fractured areas and suggests that fractures can be completely sealed by mineralization.
- Past research indicates that partially mineralized fractures have the greatest potential to stay open, but fracture type and geometry as well as hydraulic gradient can play a role.
- the present invention develops a regional “sweet spot” model, by first understanding micro-scale fluid flow and meso-scale gas diffusion and flow. Therefore, by determining fracture flow rates, fracture flow direction, and the geometry and properties of the fracture network before horizontal drilling and hydraulic fracturing the present invention improves the success rate of drilled wells.
- the present invention first considers the appropriate geochemical tracers for evaluating micro-scale and macro-scale fluid flow in fractures.
- Two tools are chosen for analyzing fluid flow in fractures: (1) Noble Gas Geochemistry (NG): He, Ne, and Ar (useful for directly quantifying fluid migration through the complex fracture network on the meso-scale and macro-scale) and ( 2 ) trace element (TE) microchemistry by Cryogenic Laser Ablation Inductively Coupled Plasma Mass Spectrometry (CLA-ICP-MS): transition metals (Mn, Fe, Ti), rare earth elements (La—Lu), and actinides (Th/U) (used to evaluate microchemistry changes ( ⁇ 5 ⁇ m scale) providing a geological record of fluid through fractures).
- CLA-ICP-MS Cryogenic Laser Ablation Inductively Coupled Plasma Mass Spectrometry
- Mn, Fe, Ti transition metals
- La—Lu rare earth elements
- Th/U actinides
- noble gases are derived from three main sources including mantle (M), crust (C), and atmosphere (A).
- M mantle
- C crust
- A atmosphere
- mantle-sourced noble gases do not play a significant role and are therefore excluded for brevity.
- Crustal (C) and atmospheric noble gases do have significant sources in such organic-rich shales, while each respective reservoir has a unique noble gas elemental and isotopic composition.
- the changes in the noble gas composition that occur as fluids migrate along fractures and interact with crustal fluids primarily relate to the radiogenic nature of the rock protolith and its geologic history.
- Uranium (U) and thorium (Th) both of which are present at relatively high concentrations in most black shales decay to 4 He (alpha-particle: ⁇ ) (i.e. (235 or 238) U and 232 Th 4 He) simultaneously producing an array of minor nuclear reactions.
- ⁇ alpha-particle
- ⁇ alpha-particle
- U and 232 Th 4 He simultaneously producing an array of minor nuclear reactions.
- Ne-21 when the alpha particle strikes an O-18 nucleus [ 18 O( ⁇ ,n) ⁇ 21 Ne].
- Other various reactions that produce Ne isotopes i.e. 24 Mg (n, ⁇ ) ⁇ 21 Ne and 3 He, 25 Mg (n, ⁇ ) ⁇ 22 Ne and 3 He, and 23 Na(n, ⁇ ) ⁇ 20 Ne and 3 He or are not significant in most crustal settings with the exception of fluorine-rich rocks that produce Ne-22.
- Black shales also contain significant amounts of potassium ( 40 K) which decays to ( 40 Ar) ( 40 K ⁇ 40 Ar) that ultimately ends up in many crustal natural gases.
- the above interactions lead to significant increases in [ 4 He] (i.e. low radiogenic or crustal 3 He/ 4 He (e.g. 1 ⁇ 10 ⁇ 8 or 0.01Ra, where Ra: 1.39 ⁇ 10 ⁇ 6 )), enriched 21 Ne/ 22 Ne (e.g. 0.035-0.050 elevated from the air value of 0.029 by nucleogenic production), and drastically increased 4 He/ 21 Ne (excess) (e.g. 20 ⁇ 10 6 ) as hydrocarbon and groundwater fluids interact with fracture surfaces.
- Atmospheric noble gases (ANG) are incorporated into crustal fluids (i.e.
- gas (and rock) samples with dominantly ASW composition may have witnessed major fracture flow that led to extensive hydrocarbon loss.
- the origin and subsurface interaction of hydrocarbons in the subsurface can be constrained, enabling the quantifying of the migration with a fractured network.
- NG Noble gas
- vein filling minerals such as calcite
- vein filling minerals incorporate the chemical composition of pore fluids during mineralization providing a geochemical archive of pore fluid chemistry throughout various flow events during vein formation.
- vein minerals record the micro-scale interactions of fluids with fracture surfaces.
- Trace element concentrations in calcite veins specifically rare earth elements (REEs), oxyanion forming trace metals (OFTM) (e.g. Mn, Fe, As), and actinides (Th and U) may be used to estimate the volume of fluid migration, number of pulses of fluid migration, the source chemistry with which fluids interact, changes in fluid chemistry (e.g. pH, redox potential) and the time of vein filling.
- REEs rare earth elements
- OFTM oxyanion forming trace metals
- Th and U actinides
- This suite of TEs i.e. oxyanions forming transition metals, REEs, and actinides
- REEs transition metals
- actinides i.e. oxyanions forming transition metals, REEs, and actinides
- Mn and As oxyanions complexes are conservative and highly mobile at a neutral to basic pH across a range of E H conditions, while REEs are only mobile in acidic and highly saline conditions and travel primarily with dissolved organic carbon (DOC) all of which lead to a measurable fraction of these elements in vein calcites.
- U and Th in country rock are fractionated when interacting with crustal fluids because U has the ability to complex with DOC or form oxyanions, leading to low relative Th/U in fluids transported along fractures as compared to fluids directly interacting with country rock.
- Comparatively immobile trace elements such as REEs and Th will increase during low fluid transport (i.e. greater interaction with country rock) and decrease during periods of high rates of fluid transport as has been observed for some transition metals including Fe and Mn.
- LA-ICP-MS enables the in situ analysis of these selected trace elements within individual veins to a resolution of several dozen microns ( ⁇ 20 ⁇ m). This capability allows microchemical spatial determination of calcite vein chemistry, a proxy for pore fluid evolution through fractures.
- a high sensitivity, fast washout cryogenic laser ablation system (GMA 4200Volante CryoCell) is used.
- This cryogenic laser ablation system cryogenically freezes the organic material enabling robust and reproducible analysis of organic-rich samples.
- the GMA Volante laser ablation cell when combined with cryogenic capability, improves analytical sensitivity ⁇ 10 ⁇ enabling analysis of REEs and actinides (Th/U) and providing sufficient spatial resolution to monitor the record of fluid flow fluctuations throughout geological time.
- the claimed combination of trace element microanalysis and noble gas techniques provides a framework to understand the characteristics, cycles, and timing of fluid flow and the potential for hydrocarbon fluid extraction.
- These claimed methodologies enable the development of basin scale maps depicting “sweet spots” (hydrocarbon-bearing and conductive along natural fracture network), “dead spots” (and highly fractured and extensively diffused without a history of crustal interaction), “potential spots” (where “tight gas” extraction will be expensive despite the presence of hydrocarbons), and ground truthing well selection procedures within hydrocarbon producing black shales.
- the present invention takes a uniquely integrated approach, using a combination of fracture network analysis and applying geochemical tools to evaluate the flow of hydrocarbons through the fracture network within hydrocarbon producing black shales.
- the present invention addresses the role of natural fractures on fluid distribution and flow in shale in four ways: (1) mapping the physical characteristics of fracture sets such as the 3-D fracture network geometry, and interpreting fracture history based on cross cutting relationships (geostructural analysis); (2) analyzing the microchemistry of individual types of open and vein-filled fractures to determine fracture interaction with fluids (cryogenic LA-ICPMS (CLA-ICP-MS) TE geochemistry) ; (3) measuring the percent change in the bulk permeability of the shale due to fracturing (NG and (CLA-ICP-MS) TE geochemistry); and (4) assessing the current basin scale variations in the NG chemistry to determine system scale gas diffusional loss.
- the present invention enables the identifying of the contributions of different variables, whether regional-scale or microscale, and the understanding of the interplay between them. This is of fundamental importance in understanding geologic fluid flow through a fractured medium and something that previous studies have failed to capture.
- the present invention comprises a suite of methods that can accurately assess modern day fluid flow through fractured rock providing near real-time in situ measures of fluid migration.
- Geostructural analysis evaluates the role of individual fracture characteristics (i.e. aperture, size, roughness, mineralization, and chemical characteristics) on fluid flow. Preliminary data finds un-mineralized, partially mineralized, and completely healed fractures in some shale. When fluids migrating along fractures become supersaturated with respect to dissolved elemental concentrations they grow crystals into open fractures, which incorporate and record elemental composition of the fluids. By examining elemental crystal microchemistry of minerals across the aperture of a fracture fluid composition changes over time can be evaluated to determine the fracture characteristics that best accommodate fluid flow and maximize extraction efficiency.
- Preliminary optical microscopy shows euhedral calcite crystals, which indicate growth into open and fluid filled fractures. Nascent crystal growth proceeds from fracture walls progressively towards the fracture opening until the fracture heals or fluid flow stops.
- Preliminary CLA-ICP-MS analyses of vein chemistry show an exciting pattern of Mn and REE concentrations in calcite crystals.
- Microsampling of an individual vein ( ⁇ 4.5 mm thick) at 10-micron resolution shows cyclic variations in REE concentrations with three peaks per cross section showing a 6-8 times increase in
- CLA-ICP-MS mapping multiple stitched lines at 5 ⁇ m resolution (spot size) to produce a map ⁇ 6 mm ⁇ 10 mm
- CLA-ICP-MS provides accurate spatial micro-sampling at geochemically relevant intervals in organic-rich samples with sufficient analytical sensitivity and specificity to accurately determine actinides and REEs in calcite veins OFTM (e.g. Mn, As), REEs, and actinides (Th/U) in vein minerals are analyzed because they are ideal tracers for the evolution of fluid chemistry and the interaction of migration fluids with fracture surfaces within country rock.
- noble gas chemistry of vein fluid inclusions is analyzed to provide a snapshot of basin chemistry at the time of vein formation.
- Noble gas composition from these fluid inclusions is compared with trace element and noble gas chemistry of vein minerals to evaluate gas diffusion as the system changes.
- the combination of these tracers can then be used to develop a gas diffusion/migration model for fluids on the microscale.
- the findings are integrated to determine the fluid flow properties of a fracture system as a whole to determine the most efficient fluid flow pathways.
- mesoscale factors stress field, fracture geometry
- the history of fracture sets in the target shale and fracture fluid flow relationships are determined. Two techniques are combined: detailed three-dimensional, sub-meter to km-scale mapping of fractures and sampling for NG and vein TE signatures to correlate fracture patterns with fluid flow.
- geostructural mapping is conducted in quarries that cut through the target shale into the underlying limestone.
- Quarrying leaves a terrace at the base of the shale providing optimal opportunity to observe 3-D exposure. Corners of a quarry have been found to allow examination of the intersecting walls and the quarry floor, providing a detailed picture of the 3-D fracture network.
- the quarry opening also allows opposite walls to be compared at ⁇ 1 km distance and large scale structures (faults, folds) to be accurately mapped and correlated with variations in fracture patterns.
- Detailed mapping is done using a meter-square grid with 10-cm subdivisions, while mapping at the 10-1000 m scale is carried out with GPS-based laser rangefinder techniques that allow positional accuracy of ⁇ 10 cm at a distance of 100 m.
- drill-cores will be available from a few select test wells in the area, so that variations between the shale at-depth and newly exposed shale at the surface in quarries can be tracked.
- Basin wide changes due to fracturing of the target shale are evaluated by combining [U] and [Th] data with NG-MS analyses.
- the anticipated 4 He/2 21 Ne production ratio is calculated in comparison to measured 4 He/ 21 Ne.
- the mechanically calculated permeability of the shale is then used to determine the partitioning of He between pore space and crystal, which values are compared with the observations to gauge the effect of different degrees of fracturing in different areas. Assuming that the 21 Ne remains in the mineral phases, the ratio of 4 He/ 21 Ne relative to production will provide evidence for the amount (volume) of fracturing and hydrocarbon loss in the shale.
Abstract
American energy costs steadily increase leading to increased unconventional energy use. However, scientific barriers prevent widespread and economic development of these resources. In gas shale plays, resource utilization is limited by the understanding of how hydrocarbons and other fluids migrate in fractured rock. This limits industry's ability to extract hydrocarbons with enhanced recovery methods; prevents exploration in areas where hydrocarbons do not collect in traditional traps; and could unfortunately lead to dangerous environmental consequences such as contaminated groundwater. To address these concerns, the present invention integrates geochemical and geostructural techniques in a novel method for evaluating and optimizing the placement and drilling strategies of extraction wells in hydrocarbon rich black shales. By correctly choosing hydrocarbon “sweet spots” companies can reduce the number of unprofitable wells, choose directional drilling and completion strategies to accurately reflect the subsurface, and better select prime small- and full-field reservoirs.
Description
- This application claims the benefit of U.S. Provisional Application No. 61/413,967, 2010 Filed Nov. 15, 2010.
- The following description relates generally to natural gas recovery, and more particularly to a method for identifying highly productive locations for hydrocarbon extraction from black shale.
- Driven by geopolitical and hydrocarbon reserve uncertainties and the continuously increasing real cost of energy, domestic energy production is necessary to ensure the energy security and independence of the United States. In order to meet these increasing energy demands with domestic production, unconventional energy resources such as shale gas, shale oil and alternative coal technologies must be increasingly utilized. Economically viable production of unconventional energy resources requires enhanced methods of oil and gas recovery such as hydraulic fracturing and horizontal drilling.
- The combination of these techniques has had mixed success at extracting economic quantities of natural gas from low permeability shale deposits, which have poor country rock permeability and transmissivity, but can contain natural fractures. Hydraulic fracturing involves the injection and extraction of fluids and propping agents in the subsurface to stimulate fluid flow through natural fractures and increase fracture related permeability (e.g. increased fracture aperture, fracture size, and fracture network connectivity) thus enhancing hydrocarbon production. The use of hydraulic fracturing is limited by: 1) inefficient resource recovery, 2) the potential for groundwater contamination from drilling fluids or mobilized hydrocarbons that can migrate through fractures and interact with groundwater, and 3) the current inability to develop accurate models for fracture fluid flow in the exceedingly complex fracture network present in black shale and other fractured lithologies. Therefore successful, economically viable, and environmentally safe application of these techniques requires a detailed understanding of fluid transport within the subsurface, specifically fluid flow within fracture networks.
- Accordingly, there is a need for improved strategies for hydrocarbon extraction. The embodiments of a predictive methodology for directly evaluating fluid flow through natural and stimulated fractures in situ by integrating noble gas geochemistry, trace element geochemistry, and fracture analysis, disclosed below, satisfies this need.
- The following simplified summary is provided in order to provide a basic understanding of some aspects of the claimed subject matter. This summary is not an extensive overview, and is not intended to identify key/critical elements or to delineate the scope of the claimed subject matter. Its purpose is to present some concepts in a simplified form as a prelude to the more detailed description that is presented later.
- Black shales are of great interest as domestic hydrocarbon plays because of their high organic content and tight gas-retaining nature. Although market demand is increasing interest in shale hydrocarbon extraction, the present inventor is unaware of any comprehensive methodology capable of identifying highly productive areas known as “sweet spots” because of a paucity of information about geological fracture-related fluid flow. Because of the immense costs associated with horizontal drilling and hydraulic fracturing, as well as the risk of environmental impact, hydrocarbon recovery from unconventional sources, such as black shale, is not economically favorable unless wells hit a sweet spot. By directly evaluating natural fluid migration in situ using conservative noble gas diffusion profiles, trace element proxies for geological fluid flow, and the characteristics of the fracture network, the presently disclosed predictive model for natural gas migration and determination of the location of sweet spots in shale has been developed. The present invention will directly contribute to a comprehensive strategy for hydrocarbon extraction, specifically in regions which lack deformed structures that lead to bed thickening shale and increased fracturing. In addition, because its techniques directly monitor fracture fluid migration on multiple geologic scales, the present invention allows for the evaluation of risk for aquifer contamination from hydraulic fluid and shale gas.
- The present invention comprises a combination of methods. Noble gas abundances and isotopic ratios (He, Ne) and trace element geochemistry (transition metals (Ti, Mn, Fe), rare earth elements (La—Lu), actinide chemistry (Th/U)) with advanced techniques in fracture network analysis are integrated, in order to: 1) determine the multi-stage fluid flow history through individual fractures; 2) quantify gas diffusion from relatively impermeable hydrocarbon host rock (i.e. black shale) to fracture sets; 3) develop a 3-D geospatial map of the regional fracture network to determine pathways of fluid migration; and 4) map chemical changes for the development of a regional hydrocarbon “sweet spot” map.
- Preliminary research focused on the New York State portion of the Marcellus shale because of its potential for hydrocarbon production and because the naturally fractured black shale provides an optimal location for describing in detail the claimed methodologies. A short synopsis of the current state of hydrocarbon production in the Marcellus shale, a review of fracture-related fluid flow, fracture network analysis techniques, and relevant geochemical proxies are presented below.
- The Marcellus shale in New York is located on the Appalachian plateau and is in the foreland of the Appalachian Fold Thrust Belt. It lies to the North and West of the Valley and Ridge province and is characterized by salt-detachment tectonics. The first order deformation structures of the Appalachian plateau are detachment folds, which have a basal decollement in the Silurian-aged Salina Salt. This salt layer acts as the glide plane for the Appalachian plateau detachment sheet and causes surface features characteristic of salt tectonics such as broad, gentle folds and abrupt changes in deformation style at the lateral and frontal termination of the salt. Salt also plays a role in faulting. Salt is thickened in the hinges of detachment anticlines and provides a weak zone that is easily faulted, resulting in blind reverse splays that cut up from the decollement in the salt and slice through this weak core.
- The stratigraphy of the Appalachian plateau varies, but in New York the Marcellus shale is the stratigraphically lowest subgroup of the siliciclastic Devonian Hamilton group. The Hamilton group is sourced from the Acadian orogeny and Rb-Sr dating of the lower Devonian black shales in the Hamilton groups puts their age at 384±9 to 377±11 Ma. The Hamilton group is sandwiched between the overlying late Devonian Catskill deltaic sequence and the underlying clean carbonates of the Devonian Onondaga limestone. The Marcellus subgroup is comprised of three formations: 1) the Union Springs-the lowest stratigraphic black shale unit; 2) the Oatka Creek Formation-the highest stratigraphic black shale; and 3) its lateral stratigraphic equivalent, the Mount Marion Formation. The Union Springs formation is made of black shales and dark grey limestones, and is separated from the black shales of the Oatka Creek and Mount Marion formations by the Cherry Valley limestone.
- The Appalachian plateau was deformed during the Pennsylvanian-Permian Alleghany orogeny. During this event, deformation progressed from the hinterland to the foreland along the basal salt layer. The Appalachian Plateau detachment sheet progressed to the northwest during this orogeny and did not interact with the underlying basement. Deformation within the detachment sheet varies with stratigraphy; in the lowermost strata shortening was accommodated by low-angle thrust faulting, while at higher levels shortening was first taken up by layer-parallel shortening and then accommodated by broad folding. In all, Alleghanian deformation occurred during two progressive stages: layer parallel shortening occurred first and was followed by a period of detachment folding and reverse faulting.
- The first stage of deformation, layer parallel shortening, is expressed in surficial geology by the deformed fossils and solution cleavage. Strains of up to 20% are observed in fossils and both solution cleavage and fossil shortening indicates a strain ellipsoid that has its short axis perpendicular to the regional structural trend. Solution cleavage maintains its bed-perpendicular orientation around folds and confirms that cleavage formed before folding.
- During the second stage of deformation, detachment folds formed by buckling above the Salina salt. These folds are characterized by comparatively tight anticlines cored by salt and broad synclines (FIG. 3). The cores of folds are significantly tighter in the salt horizon, but at higher stratigraphic levels folds are open with very gently dipping limbs (<5 degree limb dips). Folds are slightly asymmetric with steeper southeastern limbs. Thrust splays cut up from the decollement through the weak anticlinal fold-cores but do not make it to the surface and both antithetic and synthetic faults are common.
- Fracturing is a pervasive and complex feature of the Appalachian plateau that developed during successive phases of the Alleghanian orogeny and is associated with the first order structures. Fractures can be grouped into sets by their orientation, but the relationships between these sets are not completely understood. N-S striking fractures are interpreted as extension fractures due to early E-W extension in the forebulge of the Alleghanian orogeny. Cross-fold fractures (strikes ranging from 012° to 327° are difficult to interpret; explanations vary from fracturing due to multiple phases of the Alleghanian orogeny, to fracturing due to the stress field before the Alleghanian orogeny or tectonic unloading after the orogeny, or a combination of these mechanisms that invokes reactivation of fractures. ENE-striking fractures(˜071°) are interpreted as neotectonic and due to overpressure caused by hydrocarbon generation.
- During Alleghanian tectonics, fluids migrated to the Appalachian plateau from the Appalachian fold-thrust belt, with the hypotheses for the driving force ranging from a mechanical “squeegee” to a thermally driven mechanism. This fluid migration caused a widespread resetting of the magnetic signatures in the rocks, suggesting regional-scale fluid flow. This regional fluid flow utilized fractures within the Marcellus shale and is recorded in the geochemistry of the country rock and veins. Although the natural gas found in the Marcellus shale formed in place, this regional fluid flow caused the transport of fluids through fractures and is evidence of past fluid flow through the fracture network. This fluid flow may have altered gas concentrations in the Marcellus, and quantifying these changes can serve to define a model for understanding gas recovery through induced fracturing (eg. hydraulic fracturing)
- Hydrocarbons can be created through biogenic and thermogenic means, with biogenic processes being significant at shallow depths and thermogenic production dominating at deeper levels. Significant hydrocarbon generation is usually attributed to a thermogenic process; hydrocarbons are formed at depth from the thermal degradation of kerogen. As rock is buried, temperature and pressure increase and the structure of kerogen becomes unstable. Kerogen progressively adjusts to this increasing temperature and pressure by eliminating functional groups and the linkages between nuclei, thus generating a wide range of compounds including hydrocarbons, CO2, Water, and hydrogen sulfide. Additionally, natural gas comprised of methane, ethane, propane, and n-butane (C1, C2, C3, and C4, respectively) can be generated through the mechanism of transition-metal catalysis. Laboratory efforts to generate gas by purely thermal mechanisms show higher formation temperatures than observed in nature (up to 400° C.) or result in a higher fraction of heavy gasses (C2-C4) than seen in nature. However, a catalytic mechanism produces gas at low temperatures (□200° C.) with C1—C4 fractions that mimic natural gas. In particular, marine shales (such as the Marcellus shale) show an increase in released light hydrocarbons over time, indicating the opposite effect of traditional desorption. These shales generally contain the necessary transition metals for catalytic gas generation and the Marcellus shale is no exception.
- Within the last 5 years interest in natural gas production from the Marcellus shale has spiked because of the development of enhanced recovery technologies. Although there is current drilling for natural gas in the Marcellus shale in Pennsylvania, permitting issues have stalled work in New York. Current drilling efforts in Pennsylvania are focused on the hinges of anticlines where the Marcellus is thickened and where fracturing is most intense, but this strategy is not viable in New York because deformation towards the east is weaker and large scale geologic structures (e.g. folds) are more subtle. Despite political delays and geostructural challenges, there is interest in drilling in New York. The lack of structural controls and an insufficient understanding of fracture fluid flow necessitate more research in order to ensure efficient and safe production of natural gas. The lack of recent exploration in the NY Marcellus shale makes the present invention both timely and useful. With sufficient data, the present invention can develop reservoir-, field-, and regional-scale interpretations of fluid flow without the need for extrapolation; this makes the present predictive technology especially useful to local inhabitants, state governments, and hydrocarbon extractions corporations when implemented in the early stages of exploration.
- In addition to the appropriate timing for evaluation and increasing resource demands, the Marcellus shale offers logistical advantages that make it an excellent case study. For example, it is exposed in many quarries across New York State based on its stratigraphic position above the quarried Onondaga limestone. This coincidence allows study of the three-dimensional relationships of fractures with great accuracy, and sampling of fresh, unaltered, outcrops that are revealed through quarrying activities. These outcrops are not weathered and their geochemistry is preserved, making them a useful analogue for more deeply buried rocks.
- Some economically important geological advantages to studying the Marcellus shale in New York are the type of deformation, the amount of deformation, the fracture pattern, and the regional fluid flow that this area experienced. Compared with the intensely deformed sections of the Marcellus shale in the Valley and Ridge province of Pennsylvania, the Marcellus shale in New York is relatively undeformed and fracture patterns have not been overprinted by larger structures. This lack of deformation in Appalachian plateau region of the Marcellus prevented the thickening and increased fracturing of shale beds, making it very difficult to determine the best location for extraction wells and predict the location of highly fractured “sweet spots”. This necessitates the understanding of fracture related flow in the plateau section of the Marcellus shale for any drilling program.
- Fractures are surfaces in rock along which mechanical failure has occurred and the rock has lost cohesion. They can form in tension (mode 1) or shear (mode 2 or 3) and often form sets of similarly oriented members. Fractures that accommodate some degree of slip along their surfaces are faults, while fractures that have no observable slip are joints. A grouping of fractures with sub-parallel orientations is a fracture set, while all of the fractures regardless of orientation form the fracture network.
- Because black shale (country rock) is the ultimate source of hydrocarbons, the rate at which hydrocarbons diffuse into fractures and flows through the fracture network limits hydrocarbon production. The rate of hydrocarbon diffusion into fractures and through the fracture network is controlled by the hydraulic conductivity of the system (K). Hydraulic conductivity is a function of a rock's bulk porosity and permeability and describes the ease with which a fluid is transported through pore spaces or fractures.
- In black shale, the extractable volume of hydrocarbons is directly proportional to the system's hydraulic conductivity (hereafter: K), which is a factor of 1) the hydraulic conductivity of the country rock (KCR) and 2) the fracture network (hereafter: KFN). Thus fluid flow is simultaneously affected by KCR and KFN and their interaction (Eaton, 2006). In areas that have low country rock permeability (i.e. shales), flow properties are dominated by fractures, while fractures are less important in areas with higher country rock permeability (e.g. sandstones). As hydrocarbons diffuse from the country rock into fractures, the hydraulic conductivity of the fracture network (KFN) is directly related to the characteristics of the individual fractures, fracture sets, and the entire fracture network.
- Size, location, termination style, aperture, planarity and roughness are key characteristics to determine flow within individual fractures. The size of a fracture refers to the three-dimensional surface area of the fracture, while aperture is the openness of fracture planes. Planarity (a measure of a fracture's deviation from a plane) and roughness (planar tortuosity) are important factors in determining permeability. The termination style defines the geometric orientation of the end of a fracture. There are four types of terminations including: T—a perpendicular intersection between fractures; J—an intersection in which one fractures curves into the other; I—a fracture that ends at its tip line without intersection and X—cross-cutting fractures. Combinations of these termination geometries and the location of individual fractures define the 3-D geometry of a fracture network. All of the above parameters influence hydraulic conductivity (i.e. the amount of hydrocarbons that can migrate through the fracture) and each other. For example, permeability and porosity concomitantly increase with increasing country rock grain size and fracture size, and the statistical probability of fracture intersection (connectivity) also increases with larger fracture size. Fracture connectivity, porosity, and permeability all affect hydraulic conductivity implicating a complex relationship between KFN and fracture properties.
- Fractures are grouped by their geometry into fracture sets or groups of sub-parallel oriented fractures. Orientation and spacing of fracture sets in a network are characteristics that affect fluid flow in different ways. For dense and homogeneous fracture networks fluid flow can be treated as flow through a porous medium, while in sparsely fractured areas a few large fractures may dominate flow. Fracture network hydraulic properties depend on fracture intensity (surface area of fractures per unit volume), connectivity (number of fracture intersections per unit volume), hierarchy, and chronology.
- Modeling the migration of hydrocarbons in fractured black shales is exceedingly complex due in part to the complex nature of hydraulic conductivity in a fractured medium, but also to the many dynamic processes of the earth. For example, dynamic changes in parameters such as regional stress field, in fracture mineralization, fluid pressure, climatic changes (wetness/dryness), fluid gradient, anthropogenic water use, and tectonic processes reduce the accuracy of model inputs significantly and retard the understanding of fluid transport through fractured media.
- Most importantly, even at great depths under high overburden pressures, fractures must be open in order to accommodate fluid flow. How fractures remain open (aperture>0) and the relative importance of mechanical and diagenetic characteristics in keeping fractures open is still contentious. Some authors argue for the role of the in situ stress field and suggest that only fractures oriented parallel to the maximum compressive stress will stay open and accommodate fluid flow. Fractures may never completely close if there is a sufficient hydraulic gradient, even though permeability decreases significantly as stress normal to the fracture increases. In addition, some component of shearing can keep fractures open, causing asperities on opposite faces of the fracture to ride up over one another and prop open the fracture.
- The diagenetic approach to finding open fractures focuses on mineralization within open fractures and country rock stiffening due to cementation. Mineral bridges can form in fractures and cement can precipitate in the host rock holding fractures open regardless of fluid pressure or stress field changes. In the Travis Peak formation in East Texas, fluid inclusions were used to reconstruct the temperature and pressure of vein formation. Burial models suggest a 48 Myr history of vein growth, indicating that fractures were open and slowly grew minerals for an extended period of time. However, mineralization does not necessarily lead to increased permeability through fractures since complete mineralization can cause fractures to close. For example, the hydrocarbon-rich Barnett Shale of Northern Texas has fault induced fractures, but drill cores show pervasive calcite veining which correlates with low hydrocarbon production in heavily fractured areas and suggests that fractures can be completely sealed by mineralization. Past research indicates that partially mineralized fractures have the greatest potential to stay open, but fracture type and geometry as well as hydraulic gradient can play a role.
- Theoretical models and field observations suggest that, local fluid pressures can exceed lithostatic pressure, generating large hydrofractures that are capable of cutting up from a reservoir and through impermeable cap rock. Although only some large fractures cut through many stratigraphic layers, they interact with the entire fracture network by crosscutting smaller features and are capable of transporting fluids over large distances as observed in the Uinta basin, where field observations have identified natural hydrofractures that transported fluids for several kilometers vertically and tens of kilometers horizontally.
- The exhaustive list of considerations included above provides an example of the varied and complex manner in which fractures can influence fluid migration and the numerous dynamic fracture-related processes that can change both geospatially and over time. These considerations depict the difficulty of developing a theoretical model for fracture fluid flow and hydrocarbon extraction from shales. The level of current modeling capabilities, varied geological structure throughout hydrocarbon lithology, and dynamic changes within the fracture network lead to expensive and economically imprudent drilling of many failed, non-productive wells. By placing direct empirical measurement of conservative, in situ, and natural tracers for methane diffusion and fluid flow on the micro-, meso-, and macro-scale in its geostructural framework, the present invention provides a cost effective solution. The present invention develops a regional “sweet spot” model, by first understanding micro-scale fluid flow and meso-scale gas diffusion and flow. Therefore, by determining fracture flow rates, fracture flow direction, and the geometry and properties of the fracture network before horizontal drilling and hydraulic fracturing the present invention improves the success rate of drilled wells.
- The present invention first considers the appropriate geochemical tracers for evaluating micro-scale and macro-scale fluid flow in fractures. Two tools are chosen for analyzing fluid flow in fractures: (1) Noble Gas Geochemistry (NG): He, Ne, and Ar (useful for directly quantifying fluid migration through the complex fracture network on the meso-scale and macro-scale) and (2) trace element (TE) microchemistry by Cryogenic Laser Ablation Inductively Coupled Plasma Mass Spectrometry (CLA-ICP-MS): transition metals (Mn, Fe, Ti), rare earth elements (La—Lu), and actinides (Th/U) (used to evaluate microchemistry changes (˜5 μm scale) providing a geological record of fluid through fractures).
- The inert chemical nature of noble gases makes them ideal tracers of fluid origin, fluid diffusion, fluid-rock interaction, and fluid flow mass balance in the Earth's crust. In crustal fluids, including hydrocarbons, noble gases are derived from three main sources including mantle (M), crust (C), and atmosphere (A). In most organic-rich shales, mantle-sourced noble gases do not play a significant role and are therefore excluded for brevity. Crustal (C) and atmospheric noble gases, however, do have significant sources in such organic-rich shales, while each respective reservoir has a unique noble gas elemental and isotopic composition. The changes in the noble gas composition that occur as fluids migrate along fractures and interact with crustal fluids primarily relate to the radiogenic nature of the rock protolith and its geologic history. Uranium (U) and thorium (Th) (both of which are present at relatively high concentrations in most black shales decay to 4He (alpha-particle: α) (i.e. (235 or 238)U and232Th 4He) simultaneously producing an array of minor nuclear reactions. For this study, an important interaction produces Ne-21 when the alpha particle strikes an O-18 nucleus [18O(α,n)→21Ne]. Other various reactions that produce Ne isotopes (i.e. 24Mg (n,α)→21Ne and 3He, 25Mg (n,α)→22Ne and 3He, and 23Na(n,α)→20Ne and 3He or are not significant in most crustal settings with the exception of fluorine-rich rocks that produce Ne-22. Black shales also contain significant amounts of potassium (40K) which decays to (40Ar) (40K→40Ar) that ultimately ends up in many crustal natural gases. The above interactions lead to significant increases in [4He] (i.e. low radiogenic or crustal 3He/4He (e.g. 1×10−8 or 0.01Ra, where Ra: 1.39×10−6)), enriched 21 Ne/ 22Ne (e.g. 0.035-0.050 elevated from the air value of 0.029 by nucleogenic production), and drastically increased 4He/21Ne (excess) (e.g. 20×106) as hydrocarbon and groundwater fluids interact with fracture surfaces. Atmospheric noble gases (ANG) are incorporated into crustal fluids (i.e. mainly groundwater) either when water equilibrates with atmospheric gases prior to being recharged into the subsurface (termed air saturated water (ASW) or as sedimentation pore water at the time of sediment deposition. The relevant concentrations of ANG in groundwater are dependent upon temperature equilibrium at the time of recharge and the Henry's Law solubility of each noble gas where the Henry's Law constant increases in the heavier noble gases (i.e. solubility: He<Ne<Ar<Kr<Xe). In comparison to crustal gas interaction, circulating fluids with ASW composition have low [4He] (but higher 3He/4He (e.g. 1.36×10−6 or ˜0.985Ra, where Ra: 1.39×10−6)), atmospheric 21Ne/22Ne (e.g. 0.0289), and low solubility controlled 4He/21Ne (e.g. 85). Noble gas compositions with ASW composition would indicate fluid flow through permeability, highly fractured fracture network. Thus, the amount of ASW gas in a natural gas deposit (e.g. 36Ar content-ppm) is often a function of the amount of fluid flow or residual pore water. Ballentine et al. (2008) and Gilfillian (2009) have modeled these interactions in a series of papers on the major carbon dioxide rich gases of the western US. It is herein proposed that gas (and rock) samples with dominantly ASW composition may have witnessed major fracture flow that led to extensive hydrocarbon loss. By measuring the noble gas composition in pore fluids and retained in mineral lattices (country rock and vein minerals) the origin and subsurface interaction of hydrocarbons in the subsurface can be constrained, enabling the quantifying of the migration with a fractured network.
- Noble gas (NG) geochemistry has been used to constrain the permeability, effective porosity and the interaction of sedimentary basins with groundwater and to trace basin-wide migration of methane and other hydrocarbons. In addition, NG studies have identified when groundwater flow is dominated by advection through fractures and quantified interaction of fracture network fluids with surrounding rock.
- In shale, the production of 4He and 21Ne by radioactive decay of the uranium and thorium series produces an alpha particle that travels 6 to 8 microns and can either embed in a quartz grains as a He atom, or, interact with an 18O atom within the quartz to produce 21Ne. The 4He/21Ne production ratio in quartz is 2.2×107, that is one out of every 22 million alpha decays produces 21Ne in quartz. These two decay products (4He, 21Ne) are useful for tracing fluid flow because they interact with quartz crystals differently. 4He has a small atomic radius that can diffuse through quartz over geologic time scales. Over millions of years, the helium in the pore space (freely available to interact with circulating fluids) and the helium concentration in the quartz crystal reach equilibrium. 21Ne formed within the quartz grain has a larger atomic radius and has limited diffusion in quartz at room temperature but is only re-released at higher temperature or as the result of quartz breakdown. Thus, fluid flow along fractures in sedimentary basins reduces 4He concentrations in the quartz (as gas is removed with circulating crustal fluids) while 21Ne remains trapped in the quartz. Because 1He and 21Ne are produced throughout the lifetime of the shale they give an estimate of the total flow of gases since lithification and when measuring spatially within and/or near fractures can provide an estimate of volume of shale that has been degassed. Helium, which is more diffusive than methane, effectively provides a tracer of gas diffusion from the country rock into fractures with fluid mobility.
- This hypothesized behavior of the He and Ne in quartz suggests that in areas with widely spaced and flow accommodating fractures, draining of noble gasses along fractures would result in a gradational decrease in 4He/21Ne concentrations in rock as the distance to the fracture decreases (i.e. closer to fractures more degassed). In areas that have seen little fluid flow, the He/Ne ratio will approach the anticipated production ratio (e.g. 4He/21Ne: 2.2×107) as determined by measuring [U] and [Th]. Conversely, in regions close to very conductive fractures more than 95% of the helium will be lost. By contrast, a much lower 4He/21Ne ratio can be expected in areas with a high density of conductive fractures as has been observed previously (Cook et al., 1996). Measuring the 4He and 21Ne concentrations and constructing a degassing/diffusion profile is useful for identifying areas where extensive fluid flow has occurred. Testing the retained 4He (i.e. highest diffusion coefficient) in a fractured but relatively impermeable rock provides an estimate of total permeability of the formation. If noble gas ratios, specifically at depth, show an ASW profile without significant interaction with crustal fluids (e.g. 3He/4He: 0.51.0 Ra, and ASW 21Ne/22Ne: 0.029), there is a potential for extensive noble gas and hydrocarbon loss from previous fluid flow along fractures throughout geological time. These areas can be avoided when choosing where to drill. By contrast, if noble gas compositions show extensive interaction with crustal fluids and a diffused 4He/21Ne profile (much below production ratio), then we anticipate high fracture network hydraulic conductivity (KFN) without prior loss of crustal noble gases and hydrocarbons is anticipated. These locations would define sweet spots because of their high KFN will enable fluid extraction along the naturally occurring fracture networks. Alternatively, if an area has a 4He/21Ne approaching production ratios it would imply true “tight gas” country rock with a poor fracture network. While these potential plays would still retain hydrocarbons, economically viable extraction along natural fracture networks is unlikely. These areas would require more costly horizontal drilling and hydraulic fracturing, but still lead to less production making less than optimal hydrocarbon plays.
- Past research has shown that hydrocarbons, groundwater, and radiogenically produced gases interact in the subsurface with circulating fluids imprinting their chemical signature on the immobile fraction. Indeed, the movement of water in the subsurface has profound implications for collection, migration, and entrapment of natural gas and oil. Helium, methane, and water all migrate along the same fracture pathways, while 4He/21Ne and 4He record the relative amount of fluid degassing and the pathways along which water and mobilized hydrocarbon fluids travel. However, noble gas methodologies alone do not preserve a record of the volume of water flux, timing, or cycles of fluid migration in black shales.
- Conversely, vein filling minerals (such as calcite) incorporate the chemical composition of pore fluids during mineralization providing a geochemical archive of pore fluid chemistry throughout various flow events during vein formation. As a result, vein minerals record the micro-scale interactions of fluids with fracture surfaces. Trace element concentrations in calcite veins, specifically rare earth elements (REEs), oxyanion forming trace metals (OFTM) (e.g. Mn, Fe, As), and actinides (Th and U) may be used to estimate the volume of fluid migration, number of pulses of fluid migration, the source chemistry with which fluids interact, changes in fluid chemistry (e.g. pH, redox potential) and the time of vein filling.
- This suite of TEs (i.e. oxyanions forming transition metals, REEs, and actinides) is selected in order to evaluate the interaction of migrating fluids and fractured country rock during fluid migration because they have a high degree of interaction with fracture surfaces and a preference to precipitate from water and incorporate into vein forming minerals in a predictable pattern dependent upon each of their individual chemical affinities. These characteristics result in their ability to accurately record relevant changes in pH, EH (oxygen fugacity), and saturation conditions (i.e. relative volume of water flux). For example, Mn and As oxyanions complexes are conservative and highly mobile at a neutral to basic pH across a range of EH conditions, while REEs are only mobile in acidic and highly saline conditions and travel primarily with dissolved organic carbon (DOC) all of which lead to a measurable fraction of these elements in vein calcites. U and Th in country rock are fractionated when interacting with crustal fluids because U has the ability to complex with DOC or form oxyanions, leading to low relative Th/U in fluids transported along fractures as compared to fluids directly interacting with country rock. Comparatively immobile trace elements such as REEs and Th will increase during low fluid transport (i.e. greater interaction with country rock) and decrease during periods of high rates of fluid transport as has been observed for some transition metals including Fe and Mn.
- Because vein mineralization occurs slowly over geological time, exceedingly small analytical resolution is needed in order to evaluate spatial changes in chemical composition within vein minerals. The advent of high resolution LA-ICP-MS enables the in situ analysis of these selected trace elements within individual veins to a resolution of several dozen microns (˜20 μm). This capability allows microchemical spatial determination of calcite vein chemistry, a proxy for pore fluid evolution through fractures.
- However, the current state of LA-ICP-MS capabilities poses two potential problems for the analysis of vein mineralization, which include significantly lower analytical sensitivity as compared to solution based-ICP-MS and poor laser coupling with organic-rich country rock and vein minerals.
- While some trace elements are present at easily detectable concentration in vein minerals (Mn, Fe, Zn, La, Ce), these analytes can only be reliably measured at a resolution of ˜20 μm (20 μm spot size). Although this spatial resolution is markedly better than in solution analysis, optical mineralogical observations show precipitations fronts on the scale of a few microns (˜10 μm). A current laser ablation system can ablate to a spot size approaching 2 microns, but sensitivity is decreased at a smaller spot size. Additionally, even if smaller spot size reaches sufficient spatial resolution, it still does not permit robust ablation of organic-rich materials. To overcome this analytical hurdle, a high sensitivity, fast washout cryogenic laser ablation system (GMA 4200Volante CryoCell) is used. This cryogenic laser ablation system cryogenically freezes the organic material enabling robust and reproducible analysis of organic-rich samples. Additionally, the GMA Volante laser ablation cell, when combined with cryogenic capability, improves analytical sensitivity ˜10× enabling analysis of REEs and actinides (Th/U) and providing sufficient spatial resolution to monitor the record of fluid flow fluctuations throughout geological time.
- Therefore, the claimed combination of trace element microanalysis and noble gas techniques provides a framework to understand the characteristics, cycles, and timing of fluid flow and the potential for hydrocarbon fluid extraction. These claimed methodologies, enable the development of basin scale maps depicting “sweet spots” (hydrocarbon-bearing and conductive along natural fracture network), “dead spots” (and highly fractured and extensively diffused without a history of crustal interaction), “potential spots” (where “tight gas” extraction will be expensive despite the presence of hydrocarbons), and ground truthing well selection procedures within hydrocarbon producing black shales.
- Although several studies have examined the chemistry of fractured rocks in the Appalachian basin by using fluid inclusions in veins, little work has been done on the partially mineralized fractures that likely accommodate fluid flow. In addition, past research has focused on either the mechanical properties of individual fractures, or modeling flow through fracture networks based on fracture orientations. The present invention takes a uniquely integrated approach, using a combination of fracture network analysis and applying geochemical tools to evaluate the flow of hydrocarbons through the fracture network within hydrocarbon producing black shales.
- The present invention addresses the role of natural fractures on fluid distribution and flow in shale in four ways: (1) mapping the physical characteristics of fracture sets such as the 3-D fracture network geometry, and interpreting fracture history based on cross cutting relationships (geostructural analysis); (2) analyzing the microchemistry of individual types of open and vein-filled fractures to determine fracture interaction with fluids (cryogenic LA-ICPMS (CLA-ICP-MS) TE geochemistry) ; (3) measuring the percent change in the bulk permeability of the shale due to fracturing (NG and (CLA-ICP-MS) TE geochemistry); and (4) assessing the current basin scale variations in the NG chemistry to determine system scale gas diffusional loss. By looking at what controls fluid flow at different scales, the present invention enables the identifying of the contributions of different variables, whether regional-scale or microscale, and the understanding of the interplay between them. This is of fundamental importance in understanding geologic fluid flow through a fractured medium and something that previous studies have failed to capture. In addition, the present invention comprises a suite of methods that can accurately assess modern day fluid flow through fractured rock providing near real-time in situ measures of fluid migration.
- In order to understand the effect of fracture systems on fluid flow, the role of individual fractures in the process must first be understood. Geostructural analysis evaluates the role of individual fracture characteristics (i.e. aperture, size, roughness, mineralization, and chemical characteristics) on fluid flow. Preliminary data finds un-mineralized, partially mineralized, and completely healed fractures in some shale. When fluids migrating along fractures become supersaturated with respect to dissolved elemental concentrations they grow crystals into open fractures, which incorporate and record elemental composition of the fluids. By examining elemental crystal microchemistry of minerals across the aperture of a fracture fluid composition changes over time can be evaluated to determine the fracture characteristics that best accommodate fluid flow and maximize extraction efficiency.
- Preliminary optical microscopy shows euhedral calcite crystals, which indicate growth into open and fluid filled fractures. Nascent crystal growth proceeds from fracture walls progressively towards the fracture opening until the fracture heals or fluid flow stops. Preliminary CLA-ICP-MS analyses of vein chemistry show an exciting pattern of Mn and REE concentrations in calcite crystals. Microsampling of an individual vein (˜4.5 mm thick) at 10-micron resolution shows cyclic variations in REE concentrations with three peaks per cross section showing a 6-8 times increase in
- REE levels. The average wavelength of these cycles is approximately 1.3 mm and suggests the occurrence of at least three distinct fluid flow events during vein formation. CLA-ICP-MS mapping (multiple stitched lines at 5 μm resolution (spot size) to produce a map ˜6 mm×10 mm) on partially mineralized fractures and healed fractures (veins) is conducted to examine the spatial changes in trace element chemistry across individual fractures as a proxy for cycles of fluid flow and mineralization. CLA-ICP-MS provides accurate spatial micro-sampling at geochemically relevant intervals in organic-rich samples with sufficient analytical sensitivity and specificity to accurately determine actinides and REEs in calcite veins OFTM (e.g. Mn, As), REEs, and actinides (Th/U) in vein minerals are analyzed because they are ideal tracers for the evolution of fluid chemistry and the interaction of migration fluids with fracture surfaces within country rock.
- In addition to studying the cycles of fluid flow, noble gas chemistry of vein fluid inclusions is analyzed to provide a snapshot of basin chemistry at the time of vein formation. Noble gas composition from these fluid inclusions is compared with trace element and noble gas chemistry of vein minerals to evaluate gas diffusion as the system changes. The combination of these tracers can then be used to develop a gas diffusion/migration model for fluids on the microscale.
- After determining the fluid flow characteristics of individual fractures (microscale), the findings are integrated to determine the fluid flow properties of a fracture system as a whole to determine the most efficient fluid flow pathways. To understand the role of the mesoscale factors (stress field, fracture geometry) on fluid flow, the history of fracture sets in the target shale and fracture fluid flow relationships are determined. Two techniques are combined: detailed three-dimensional, sub-meter to km-scale mapping of fractures and sampling for NG and vein TE signatures to correlate fracture patterns with fluid flow.
- In one embodiment of the present invention, geostructural mapping is conducted in quarries that cut through the target shale into the underlying limestone. Quarrying leaves a terrace at the base of the shale providing optimal opportunity to observe 3-D exposure. Corners of a quarry have been found to allow examination of the intersecting walls and the quarry floor, providing a detailed picture of the 3-D fracture network. The quarry opening also allows opposite walls to be compared at ˜1 km distance and large scale structures (faults, folds) to be accurately mapped and correlated with variations in fracture patterns. Detailed mapping is done using a meter-square grid with 10-cm subdivisions, while mapping at the 10-1000 m scale is carried out with GPS-based laser rangefinder techniques that allow positional accuracy of <10 cm at a distance of 100 m.
- Preliminary research has shown multiple fracture sets of different types and orientations in the target shale formations. In addition to joints, conjugate sets of transverse shear fractures, extensional joints and low angle fractures with reverse shear are shown. This allows observation of calcite veins and partially mineralized fractures with calcite bridges that record several generations of fluid flow as described earlier.
- When fracture network geometry from field studies is combined with NG isotopic data, it can determine which fracture sets accommodate fluid flow and identify sweet spots, as described earlier from work on groundwater systems in fractured rock. Samples for NG-MS will be collected from different sets of fractures and host rocks in a variety of locations using a 1-inch core drill. Quantitative fluid flow data is compared with fracture patterns to constrain fracture network effects, and group fractures so that the properties of individual fractures can be used to further understand the flow of hydrocarbons in the target shale.
- Quarries from different parts of the target outcrop belt are studied in this way so that the data can be compared to understand regional variations. In one embodiment of the present invention, drill-cores will be available from a few select test wells in the area, so that variations between the shale at-depth and newly exposed shale at the surface in quarries can be tracked.
- Basin wide changes due to fracturing of the target shale are evaluated by combining [U] and [Th] data with NG-MS analyses. By measuring shale [U] and [Th], the anticipated 4He/221Ne production ratio is calculated in comparison to measured 4He/21Ne. The mechanically calculated permeability of the shale is then used to determine the partitioning of He between pore space and crystal, which values are compared with the observations to gauge the effect of different degrees of fracturing in different areas. Assuming that the 21Ne remains in the mineral phases, the ratio of 4He/21Ne relative to production will provide evidence for the amount (volume) of fracturing and hydrocarbon loss in the shale.
- It can thus clearly be seen that the predictive methodology for directly evaluating fluid flow through natural and stimulated fractures in situ by integrating noble gas geochemistry, trace element geochemistry, fracture analysis, and regional structural geology is a significant improvement over the extant shale hydrocarbon extraction methods. Not only is well-selection success rate improved and the use of geologic features maximized, but hydrocarbon extraction is improved and recovery costs are reduced dramatically
- What has been described above includes examples of one or more embodiments. It is, of course, not possible to describe every conceivable combination of components or methodologies for purposes of describing the aforementioned embodiments, but one of ordinary skill in the art may recognize that many further combinations and permutations of various embodiments are possible. Accordingly, the described embodiments are intended to embrace all such alterations, modifications and variations that fall within the spirit and scope of the appended claims. Furthermore, to the extent that the term “includes” is used in either the detailed description or the claims, such term is intended to be inclusive in a manner similar to the term “comprising” as “comprising” is interpreted when employed as a transitional word in a claim.
Claims (1)
1. A method comprising:
a. Analyzing noble gases to determine the style of fluid migration in the sub-surface;
b. Analyzing noble gases to distinguish fracture density; and
c. Optimizing the direction and orientation of fluid migration with noble gas and trace element chemistry
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/297,263 US20120134749A1 (en) | 2010-11-15 | 2011-11-15 | Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US41396710P | 2010-11-15 | 2010-11-15 | |
US13/297,263 US20120134749A1 (en) | 2010-11-15 | 2011-11-15 | Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales |
Publications (1)
Publication Number | Publication Date |
---|---|
US20120134749A1 true US20120134749A1 (en) | 2012-05-31 |
Family
ID=46126764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/297,263 Abandoned US20120134749A1 (en) | 2010-11-15 | 2011-11-15 | Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales |
Country Status (1)
Country | Link |
---|---|
US (1) | US20120134749A1 (en) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140256055A1 (en) * | 2011-11-11 | 2014-09-11 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons |
US20150032425A1 (en) * | 2013-07-25 | 2015-01-29 | Halliburton Energy Services, Inc. | Determining Flow Through A Fracture Junction In A Complex Fracture Network |
US20150361792A1 (en) * | 2013-03-08 | 2015-12-17 | Halliburton Energy Services, Inc. | Systems and methods for optimizing analysis of subterranean well bores and fluids using noble gases |
US9891331B2 (en) | 2014-03-07 | 2018-02-13 | Scott C. Hornbostel | Exploration method and system for detection of hydrocarbons from the water column |
US9890617B2 (en) | 2014-09-18 | 2018-02-13 | Exxonmobil Upstream Research Company | Method to determine the presence of source rocks and the timing and extent of hydrocarbon generation for exploration, production and development of hydrocarbons |
US10094815B2 (en) | 2014-09-18 | 2018-10-09 | Exxonmobil Upstream Research Company | Method to enhance exploration, development and production of hydrocarbons using multiply substituted isotopologue geochemistry, basin modeling and molecular kinetics |
US10132144B2 (en) | 2016-09-02 | 2018-11-20 | Exxonmobil Upstream Research Company | Geochemical methods for monitoring and evaluating microbial enhanced recovery operations |
US10309217B2 (en) | 2011-11-11 | 2019-06-04 | Exxonmobil Upstream Research Company | Method and system for reservoir surveillance utilizing a clumped isotope and/or noble gas data |
US10400596B2 (en) | 2014-09-18 | 2019-09-03 | Exxonmobil Upstream Research Company | Method to enhance exploration, development and production of hydrocarbons using multiply substituted isotopologue geochemistry, basin modeling and molecular kinetics |
US10415379B2 (en) | 2015-02-03 | 2019-09-17 | Exxonmobil Upstream Research Company | Applications of advanced isotope geochemistry of hydrocarbons and inert gases to petroleum production engineering |
US10494923B2 (en) | 2014-09-18 | 2019-12-03 | Exxonmobil Upstream Research Company | Method to perform hydrocarbon system analysis for exploration, production and development of hydrocarbons |
US10533414B2 (en) | 2015-02-03 | 2020-01-14 | Michael Lawson | Applications of advanced isotope geochemistry of hydrocarbons and inert gases to petroleum production engineering |
US11041384B2 (en) | 2017-02-28 | 2021-06-22 | Exxonmobil Upstream Research Company | Metal isotope applications in hydrocarbon exploration, development, and production |
US11237146B2 (en) | 2015-03-02 | 2022-02-01 | Exxonmobil Upstream Research Company | Field deployable system to measure clumped isotopes |
WO2022060393A1 (en) * | 2020-09-21 | 2022-03-24 | Saudi Arabian Oil Company | Systems and methods for identifying gas migration using helium |
CN114264680A (en) * | 2021-11-15 | 2022-04-01 | 中煤科工集团西安研究院有限公司 | Method for predicting fluoride concentration in mine water based on analogy method |
WO2023239798A1 (en) * | 2022-06-07 | 2023-12-14 | Koloma, Inc. | Systems and methods for monitoring, quantitative assessment, and certification of low-carbon hydrogen and derivative products |
US11867682B2 (en) | 2020-09-21 | 2024-01-09 | Baker Hughes Oilfield Operations Llc | System and method for determining natural hydrocarbon concentration utilizing isotope data |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2318689A (en) * | 1943-05-11 | Tracing gas through underground | ||
US3835710A (en) * | 1971-05-03 | 1974-09-17 | L Pogorski | Geo-chemical exploration method |
US4378055A (en) * | 1977-08-05 | 1983-03-29 | Phillips Petroleum Company | Analyzing for helium in drilling muds to locate geothermal reservoirs |
US4537062A (en) * | 1982-10-19 | 1985-08-27 | Kohlensaure-Werke Rudolf Buse Sohn Gmbh & Co. | Method and apparatus for investigating the structure and porosity of earth and stony regions |
US5286651A (en) * | 1989-08-24 | 1994-02-15 | Amoco Corporation | Determining collective fluid inclusion volatiles compositions for inclusion composition mapping of earth's subsurface |
US5328849A (en) * | 1989-08-24 | 1994-07-12 | Amoco Corporation | Inclusion composition mapping of earth's subsurface using collective fluid inclusion volatile compositions |
US5416024A (en) * | 1989-08-24 | 1995-05-16 | Amoco Corporation | Obtaining collective fluid inclusion volatiles for inclusion composition mapping of earth's subsurface |
US20070169540A1 (en) * | 2001-06-02 | 2007-07-26 | Sterner Steven M | Method and apparatus for determining gas content of subsurface fluids for oil and gas exploration |
US20110030465A1 (en) * | 2008-04-09 | 2011-02-10 | Philip Craig Smalley | Geochemical surveillance of gas production from tight gas fields |
US20130091925A1 (en) * | 2011-10-12 | 2013-04-18 | Thomas Henry Darrah | Method for Determining the Genetic Fingerpring of Hydrocarbons and Other Geological Fluids using Noble Gas Geochemistry |
-
2011
- 2011-11-15 US US13/297,263 patent/US20120134749A1/en not_active Abandoned
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2318689A (en) * | 1943-05-11 | Tracing gas through underground | ||
US3835710A (en) * | 1971-05-03 | 1974-09-17 | L Pogorski | Geo-chemical exploration method |
US4378055A (en) * | 1977-08-05 | 1983-03-29 | Phillips Petroleum Company | Analyzing for helium in drilling muds to locate geothermal reservoirs |
US4537062A (en) * | 1982-10-19 | 1985-08-27 | Kohlensaure-Werke Rudolf Buse Sohn Gmbh & Co. | Method and apparatus for investigating the structure and porosity of earth and stony regions |
US5286651A (en) * | 1989-08-24 | 1994-02-15 | Amoco Corporation | Determining collective fluid inclusion volatiles compositions for inclusion composition mapping of earth's subsurface |
US5328849A (en) * | 1989-08-24 | 1994-07-12 | Amoco Corporation | Inclusion composition mapping of earth's subsurface using collective fluid inclusion volatile compositions |
US5416024A (en) * | 1989-08-24 | 1995-05-16 | Amoco Corporation | Obtaining collective fluid inclusion volatiles for inclusion composition mapping of earth's subsurface |
US20070169540A1 (en) * | 2001-06-02 | 2007-07-26 | Sterner Steven M | Method and apparatus for determining gas content of subsurface fluids for oil and gas exploration |
US20110030465A1 (en) * | 2008-04-09 | 2011-02-10 | Philip Craig Smalley | Geochemical surveillance of gas production from tight gas fields |
US8505375B2 (en) * | 2008-04-09 | 2013-08-13 | Bp Exploration Operating Company Limited | Geochemical surveillance of gas production from tight gas fields |
US20130091925A1 (en) * | 2011-10-12 | 2013-04-18 | Thomas Henry Darrah | Method for Determining the Genetic Fingerpring of Hydrocarbons and Other Geological Fluids using Noble Gas Geochemistry |
Non-Patent Citations (14)
Title |
---|
Bhullar, A. G. et al, Organic Geochemistry 1998, 29, 735-768. * |
Etiope, G. et al, Physics of the Earth and Planetary Interiors 2002, 129, 185-204. * |
Fengjun, N. et al, Marine and Petroleum Geology 2001, 18, 561-575. * |
Gratier, J.-P. et al, Journal of Geophysical Research 2003, 108, B2, 2104, 25 pages. * |
Hao, F. et al, Marine and Petroleum Geology 2007, 24, 1-13. * |
Heath, J. et al, Energy Procedia 2009, 1, 2903-2910. * |
Huang, B. et al, Organic Geochemistry 2003, 34, 1009-1025. * |
Knies, J. et al, Geochemistry Geophysics Geosystems 2004, 5, 14 pages. * |
Mazurek, M. Applied Geochemistry 1999, 15, 211-234. * |
Shi, D. et al, Organic Geochemistry 2005, 36, 2003-223. * |
Torgersen, T. et al, Earth and Planetary Science Letters 2004, 226, 477-489. * |
Zeng, J. et al, Journal of Geochenical Exploration 2006, 89, 455-459. * |
Zhang, S. et al, Marine and Petroleum Geology 2004, 21, 651-668. * |
Zhou, Z. et al, Geochimica et Cosmochimica Acta 2005, 69, 5413-5428. * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10309217B2 (en) | 2011-11-11 | 2019-06-04 | Exxonmobil Upstream Research Company | Method and system for reservoir surveillance utilizing a clumped isotope and/or noble gas data |
US9146225B2 (en) | 2011-11-11 | 2015-09-29 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons with an underwater vehicle |
US9612231B2 (en) * | 2011-11-11 | 2017-04-04 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons |
US20140256055A1 (en) * | 2011-11-11 | 2014-09-11 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons |
US10527601B2 (en) | 2011-11-11 | 2020-01-07 | Exxonmobil Upstream Research Company | Method for determining the location, size, and fluid composition of a subsurface hydrocarbon accumulation |
US10330659B2 (en) | 2011-11-11 | 2019-06-25 | Exxonmobil Upstream Research Company | Method for determining the location, size, and fluid composition of a subsurface hydrocarbon accumulation |
US20150361792A1 (en) * | 2013-03-08 | 2015-12-17 | Halliburton Energy Services, Inc. | Systems and methods for optimizing analysis of subterranean well bores and fluids using noble gases |
US10060258B2 (en) * | 2013-03-08 | 2018-08-28 | Halliburton Energy Services, Inc. | Systems and methods for optimizing analysis of subterranean well bores and fluids using noble gases |
US20150032425A1 (en) * | 2013-07-25 | 2015-01-29 | Halliburton Energy Services, Inc. | Determining Flow Through A Fracture Junction In A Complex Fracture Network |
US9418184B2 (en) * | 2013-07-25 | 2016-08-16 | Halliburton Energy Services, Inc. | Determining flow through a fracture junction in a complex fracture network |
US9891331B2 (en) | 2014-03-07 | 2018-02-13 | Scott C. Hornbostel | Exploration method and system for detection of hydrocarbons from the water column |
US10494923B2 (en) | 2014-09-18 | 2019-12-03 | Exxonmobil Upstream Research Company | Method to perform hydrocarbon system analysis for exploration, production and development of hydrocarbons |
US10094815B2 (en) | 2014-09-18 | 2018-10-09 | Exxonmobil Upstream Research Company | Method to enhance exploration, development and production of hydrocarbons using multiply substituted isotopologue geochemistry, basin modeling and molecular kinetics |
US10400596B2 (en) | 2014-09-18 | 2019-09-03 | Exxonmobil Upstream Research Company | Method to enhance exploration, development and production of hydrocarbons using multiply substituted isotopologue geochemistry, basin modeling and molecular kinetics |
US9890617B2 (en) | 2014-09-18 | 2018-02-13 | Exxonmobil Upstream Research Company | Method to determine the presence of source rocks and the timing and extent of hydrocarbon generation for exploration, production and development of hydrocarbons |
US10533414B2 (en) | 2015-02-03 | 2020-01-14 | Michael Lawson | Applications of advanced isotope geochemistry of hydrocarbons and inert gases to petroleum production engineering |
US10415379B2 (en) | 2015-02-03 | 2019-09-17 | Exxonmobil Upstream Research Company | Applications of advanced isotope geochemistry of hydrocarbons and inert gases to petroleum production engineering |
US11237146B2 (en) | 2015-03-02 | 2022-02-01 | Exxonmobil Upstream Research Company | Field deployable system to measure clumped isotopes |
US10132144B2 (en) | 2016-09-02 | 2018-11-20 | Exxonmobil Upstream Research Company | Geochemical methods for monitoring and evaluating microbial enhanced recovery operations |
US11041384B2 (en) | 2017-02-28 | 2021-06-22 | Exxonmobil Upstream Research Company | Metal isotope applications in hydrocarbon exploration, development, and production |
WO2022060393A1 (en) * | 2020-09-21 | 2022-03-24 | Saudi Arabian Oil Company | Systems and methods for identifying gas migration using helium |
US11656211B2 (en) | 2020-09-21 | 2023-05-23 | Saudi Arabian Oil Company | Systems and methods for identifying gas migration using helium |
US11867682B2 (en) | 2020-09-21 | 2024-01-09 | Baker Hughes Oilfield Operations Llc | System and method for determining natural hydrocarbon concentration utilizing isotope data |
CN114264680A (en) * | 2021-11-15 | 2022-04-01 | 中煤科工集团西安研究院有限公司 | Method for predicting fluoride concentration in mine water based on analogy method |
WO2023239798A1 (en) * | 2022-06-07 | 2023-12-14 | Koloma, Inc. | Systems and methods for monitoring, quantitative assessment, and certification of low-carbon hydrogen and derivative products |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20120134749A1 (en) | Using noble gas geochemistry to evaluate fluid migration in hydrocarbon bearing black shales | |
Mazurek et al. | Natural tracer profiles across argillaceous formations | |
Rutqvist et al. | Coupled reservoir-geomechanical analysis of CO2 injection at In Salah, Algeria | |
LeBlanc | High resolution sequence stratigraphy and reservoir characterization of the" Mississippian Limestone" in north-central Oklahoma | |
Sundal et al. | Modelling CO2 migration in aquifers; considering 3D seismic property data and the effect of site-typical depositional heterogeneities | |
Beaudoin et al. | Regional-scale paleofluid system across the Tuscan Nappe–Umbria Marche Arcuate Ridge (northern Apennines) as revealed by mesostructural and isotopic analyses of stylolite-vein networks | |
Carminati et al. | Control of Cambrian evaporites on fracturing in fault-related anticlines in the Zagros fold-and-thrust belt | |
Gabellone et al. | Fluid channeling along thrust zones: the Lagonegro case history, southern Apennines, Italy | |
Riestenberg et al. | Project ECO2S: Characterization of a world class carbon dioxide storage complex | |
Mazurek et al. | Derivation and application of a geologic dataset for flow modelling by discrete fracture networks in low-permeability argillaceous rocks | |
Bruhn et al. | Tectonics, fluid migration, and fluid pressure in a deformed forearc basin, Cook Inlet, Alaska | |
Yousef et al. | Microfracture Characterization in Sandstone Reservoirs: A Case Study from the Upper Triassic of Syria’s Euphrates Graben | |
Dick et al. | The internal architecture and permeability structures of faults in shale formations | |
Ryerson et al. | Natural CO2 accumulations in the western Williston Basin: A mineralogical analog for CO2 injection at the Weyburn site | |
Collins et al. | Natural Fractures in the Niobrara Formation, Boulder to Lyons, Colorado | |
Feng et al. | Mineral filling pattern in complex fracture system of carbonate reservoirs: implications from geochemical modeling of water-rock interaction | |
Dembicki Jr et al. | Lessons learned from the Floyd shale play | |
Roure | Crustal architecture, thermal evolution and energy resources of compressional basins (André Dumont medallist lecture 2013) | |
Şen | Producible fluid oil saturations of the Upper Cretaceous unconventional carbonate plays, northern Arabian plate | |
Nicolas et al. | Introduction to the GEM-2 Hudson–Ungava Project, Hudson Bay Lowland, northeastern Manitoba | |
Lotfiyar et al. | Geochemical, geological, and petrophysical evaluation of Garau Formation in Lurestan basin (west of Iran) as a shale gas prospect | |
Löfgren et al. | Tracer tests-possibilities and limitations. Experience from SKB fieldwork: 1977-2007 | |
Murray | Mechanical stratigraphy and sonic log relationships using the Proceq Bambino in the Niobrara Formation, Denver Basin | |
Brosse et al. | CO2 storage in the struggle against climate change | |
Fitz-Díaz | Progressive deformation, fluid flow and water-rock interaction in the Mexican Fold-Thrust Belt, Central Mexico |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |