Patent Application
20160139293
This disclosure relates to the field of laboratory whole core analysis as well as the field of well logging.
Level of organic maturity (“LOM”) is a parameter used in the assessment of the commercial value of a shale reservoir. LOM impacts interstitial fluids, production limits, pore pressure, and reservoir quality. One indicator of LOM is vitrinite reflectance, which is measured in a laboratory.
The presence of pyrite, graphite, and turbostatic carbon nanostructures (“TCN”) produces interfacial polarization that results in frequency-based dispersion of effective electrical properties of shale formations. Locating a concurrence of favorable geologic parameters such as thermal history, gas content, reservoir thickness, fractures, rock composition, total organic carbon (“TOC”), porosity, permeability, free and sorbed gases, water saturations, and thermal maturity assists in finding gas shale reservoirs. Gas is stored interstitially within pore spaces between rock grains, fractures, or kerogen, or it can be adsorbed to the surface of organic components contained within the shale. Kerogen is categorized into four broad groups, each of which has distinct bearing on what type of hydrocarbons will be produced. Type 2 kerogen is generated from the remains of planktons in reducing environments found in moderately deep marine settings. Sulfur is associated as pyrite, free sulfur, or in the organic structure of kerogen. Type 3 kerogen is primarily derived from terrestrial plant debris that has been deposited in shallow to deep marine or non-marine environments. Type 3 kerogen may generate dry gas, as it has lower hydrogen and oxygen content. Marine shale facies contain type 2 and type 3 blends. Heat and time convert organic material to hydrocarbons. Temperature and pressure increase during burial, resulting in organic material giving off oil and gas. Diagenesis refers to the low-temperature alteration of organic material due to oxidation, chemical processes, and biological processes. With increases in temperature and changes in pH, organic material is converted into kerogen and lesser amounts of bitumen. Pyritization occurs during the diagenesis phase when sulfate-reducing bacterial colonies sulfate in sea water to biodegrade organic material. Hydrogen sulfide and native sulfur released by these bacteria combine with iron in clays to form pyrites (FeS2). Catagenesis causes chemical bonds to break down in shale and the kerogen. Oil is produced from type 1 kerogen, waxy oil is produced from type 2 kerogen, and gas is produced from type 3 kerogen. Pyrite and marcasite are two sulfides in iron-rich sedimentary rocks: black shale and coal, respectively. Many shale reservoirs contain up to 10 wt % of pyrite and TOC that amount up to 7% pyrite and up to 20% kerogen by volume.
Shale formations may contain 12% to 15% kerogen volume. Spectral elemental analysis using geochemical logs is used to solve for the elements present in the formation. However, the volume of kerogen cannot be determined with spectral tools due to the presence of carbon in various other minerals. The Passey method is designed for use in the evaluation of organic content of hydrocarbon source rocks. The method uses knowledge of the maturity of the organic material, and is less accurate for sediments that are overmature, such as shale gas formations. Frequently, layers including a higher organic content are interspersed with layers of lower organic content or conductive minerals resulting in dielectric anisotropy that adversely affects the resistivity interpretation. Therefore, laboratory measurements of thermal maturation is another method and includes measurements of thermal indicators such as vitrinite reflectivity, LOM, spore color index, conodont alteration, sulfur content, hydrogen index, and pyrolysis byproducts. Vitrinite is a good indicator of the maximum temperature that a shale has been exposed to but some shales do not contain vitrinite.
Electrical conductivity values in the induction tool frequency range, and permittivity values in the propagation tool frequency range, are estimated using existing resistivity interpretation techniques to primarily evaluate water saturation in a conventional geological sample or the geological formation. The effects of pyrite, graphite, and other constituent minerals are not included in the resistivity interpretation for shale gas characterization. Source rocks are commonly shales and lime-mudstones that contain large amounts of organic matter. Richness and maturity of source rocks is evaluated through a variety of laboratory analyses like TOC analysis, pyrolysis, elemental analysis, vitrinite reflectance, thermal alteration index, gas chromatography, and visual kerogen description.
The proposed method of subsurface estimation of LOM as a continuous realtime well log may be used for directional drilling and well placement in shale, organic-rich mudrock, and source-rock reservoirs. A drilling engineer and well placement engineer may use the continuous assessment of LOM to ensure that the wellbore stays in the geological formation with the highest LOM, thereby improving exposure of the wellbore to hydrocarbon reserves.
The proposed method of subsurface estimation of LOM as a continuous realtime well log may enable directional drillers and well placement engineeers to steer away from overmature rocks containing graphite and turbostatic carbon and to stay within pyrite-rich hydrocarbon-bearing rocks. The proposed method of subsurface estimation of LOM as a continuous well log may also enable reservoir engineers to plan lateral wells into zones rich in mature organics and avoid overmature zones. The proposed method of subsurface estimation of LOM as a continuous realtime well log may also provide extra geological information for accurate well placement.
In one embodiment, a method disclosed herein includes characterizing the type of shale gas reservoirs. The method also includes estimating of volume fractions of electrically-active sites such as pyrite, graphite, pyrrhotite, and recently-discovered turbostatic carbon nanostructures that produce electrochemically-inactive interfacial polarization effects in the 1 kHz to 10 MHz frequency range. The method also includes identifying the LOM of a shale gas reservoir. The method also includes improving resistivity interpretations in shale gas reservoirs for improved TOC estimation.
A method to obtain a continuous estimation of LOM, to improve TOC estimation based on delta-log-R method, and to classify the shale gas formations based on the complex conductivity dispersive response based on multi-frequency complex conductivity measurement is disclosed. This method uses a characterization of TOC and LOM of source rock either as a standalone measurement or in conjunction with other wireline log or laboratory measurements such a nuclear spectroscopy tools. The method is based on quantification of various types of dielectric polarization phenomena occurring in a geological material within MHz to GHz frequency range. This includes surface conduction and membrane polarization (clay, <1 kHz), electrochemical spectral induced polarization (clay and other minerals, <1 kHz), metallic interfacial polarization (conductive minerals in absence of redox-active ions, 1 kHz to 10 MHz), Maxwell-Wagner polarization (interstitial pore structure, 500 KHz to 100 MHz), dipolar polarization (pore fluids, includes orientation and distortion polarization, >10 MHz), and ionic polarization. This method can be also implemented for gathering additional information for laboratory diagnosis of LOM. Existing resistivity tools may utilize excitation frequencies from near-DC like galvanic resistivity tools. For example, the frequency may range from 0 Hz (e.g., near-DC) to 100 Hz. Existing resistivity tools may utilize excitation frequencies from very high frequency dielectric tools operating in the GHz range. For example, the frequency may range from 10 MHz to 1 GHz or more. This may aid in determining the TOC.
In one embodiment, a method for determining a level of organic maturity of a shale gas formation includes inverting multifrequency complex conductivity data to estimate a volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The inversion of the multifrequency complex conductivity data may be validated using estimates of the volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The volume fraction of graphite and turbostatic carbon nanostructures may be correlated to a level of organic maturity log of the shale gas formation. The level of organic maturity log may be validated using sulfur content obtained from pyrolysis or vitrinite reflectance. The method may also include quantifying a variation of an electromagnetic response due to the volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The electromagnetic response may be modified by removing the quantified variation to obtain resistivity and permittivity values.
A computing system is also disclosed. The computing system includes a processor and a memory system including a non-transitory computer-readable medium storing instructions that, when executed by the processor, causes the computing system to perform operations. The operations include inverting multifrequency complex conductivity data to estimate a volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The inversion of the multifrequency complex conductivity data may be validated using estimates of the volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The volume fraction of graphite and turbostatic carbon nanostructures may be correlated to a level of organic maturity log of the shale gas formation. The level of organic maturity log may be validated using sulfur content obtained from pyrolysis or vitrinite reflectance. The operations may also include quantifying a variation of an electromagnetic response due to the volume fraction of graphite, turbostatic carbon nanostructures, and pyrite. The electromagnetic response may be modified by removing the quantified variation to obtain resistivity and permittivity values.
Other aspects and advantages of the claimed subject matter will be apparent from the following description and the appended claims.
One or more specific embodiments of the present disclosure are described below. These embodiments are merely examples of the presently disclosed techniques. Additionally, in an effort to provide a concise description of these embodiments, all features of an actual implementation may not be described in the specification. It should be appreciated that in the development of any such implementation, as in any engineering or design project, numerous implementation-specific decisions are made to achieve the developers' specific goals, such as compliance with system-related and business-related constraints, which may vary from one implementation to another. Moreover, it should be appreciated that such development efforts might be complex and time consuming, but would nonetheless be a routine undertaking of design, fabrication, and manufacture for those of ordinary skill having the benefit of this disclosure.
When introducing elements of various embodiments of the present disclosure, the articles “a,” “an,” and “the” are intended to mean that there are one or more of the elements. The embodiments discussed below are intended to be examples that are illustrative in nature and should not be construed to mean that the specific embodiments described herein are necessarily preferential in nature. Additionally, it should be understood that references to “one embodiment” or “an embodiment” within the present disclosure are not to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features.
A method based on multifrequency complex conductivity tensor measurements and modeling of polarization phenomena in shale gas formation is disclosed. For purposes of demonstration of this method, this discussion will focus on electrochemically-inactive interfacial polarization. The method relates to formation evaluation and petrophysical interpretation. More particularly, the method relates to formation evaluation in gas shale and oil-bearing shale reservoirs. The method also relates to techniques for evaluating a subterranean shale gas, shaly sand formations, and whole core samples to identify and quantify LOM and TOC.
This method is used to identify and quantify the volume fraction of minerals associated with organic matter, thereby indicating the organic richness and maturity. In clay-rich source rock, mineral grains align horizontally, and organic matter is distributed in sub-horizontal lamellae. In carbonate source rocks, matrix grains don't have a specific orientation, and such rocks are laminated with organic rich and organic-lean layers. In this methodology, the geological samples are assumed to be composed of a non-dispersive, non-dielectric, and non-conductive rock matrix, dispersive, dielectric, and conductive rock minerals, organic matter, and fluids filling the pore space. Immature source rocks are conductive due to water-filled pore space, and have non-dispersive R-signals. Such rocks have dielectric and dispersive X-signals when pyrite minerals, a product of diagenesis during early phase of maturation, are present. Over-matured rocks are conductive due to the presence of graphite due to excessive graphitization, and both the R-signal and X-signal are dispersive. While matured source rocks are non-conductive due to the presence of hydrocarbon-filled pore space, they can be slightly conductive due to presence of graphite, which is a product of metagenesis. As source rock matures, a portion of solid organic matter is transformed to liquid and part to graphite or other precursors which increase in conductivity and dispersive R-signal response. Increases in the depth of burial and temperature may result in maturation of shale gas source rock, which may result in the formation of pyrite due to the presence of organic material and graphite from graphitization of kerogen material. Pulsed-neutron spectral logs can provide carbon-oxygen neutron logs that help in evaluating TOC content.
This method may use a subsurface diagnostic tool for assessing the maturation of source rock. With increases in maturation, graphite concentration increases as bonds are broken and hydrogen and oxygen is expelled from the source rock with graphitic pyrobitumen. Thermogenic shale gas reservoir properties may be present in post-mature Marcellus shale gas prospects (e.g., the Laughrey Ruble Black shale diagenesis). Late-stage methane during metagenesis may be present. Since graphite is a paleothermometer, and pyrite is an indicator of the presence of initial organic material, the concentrations of these minerals can be correlated to hydrocarbon presence, source rock maturation, and hydrocarbon generation. A calibration with core analysis for LOM and TOC may then be conducted. This calibration may be used to provide a more accurate well placement and to enhance directional drilling.
Shale gas is produced from formations that are heterogeneous, but kerogen maturity follows a relationship with depth of burial, temperature, and pressure. As a general description, type 1 is algal, type 2 is liptinitic, type 3 is humic kerogen, while type 4 contains small amounts of H-atom due to aromaticity of carbon and higher amounts of oxygen compared to type 1. Graphite can form an interconnected thin film around the grains by metagenesis of C—O—H fluids. Geologic energy on organic deposits may produce many types of natural graphites, such as flake and amorphous.
A mechanistic model may be implemented with a conductivity tensor tool forward model to differentiate between metallic-luster minerals, graphites, and clays. Modeling of interfacial polarization can be done using circuit models, empirical models, phenomenological models, or mechanistic models. A physically consistent model satisfies the Kronig-Kramers relationship. Some of the analytical models are built for cases when dispersed phases do not interact with each other. As a result, a single element is analyzed, and an effective medium formulation to the entire mixture is applied.
Conductive minerals such as pyrite, pyrrhotite, and graphite commonly occur in source rocks and shale reservoirs. Also, isolated conductive turbostatic carbon nanostructures, whose volume fraction increases with increasing thermal maturity, have been found in organic-rich shales. These electrically-active phases in the absence of redox-active ions, possess electrochemically-inactive phases and boundaries, and the resulting geological mixture exhibits large dielectric polarization referred to as interfacial polarization. During this phenomenon, there may be spatiotemporal variation of dynamics and distribution of charge carriers around the phase boundary. Therefore, the electromagnetic (“EM”) response of shale reservoir rocks may be influenced by the electrical properties of these electrically-active phases and electrodiffusion processes. As a consequence, the EM response may provide an inaccurate hydrocarbon saturation estimation based on conventional resistivity interpretation techniques. Analytical modeling techniques have been developed and numerical modeling schemes implemented to understand and quantify this phenomenon. An EM whole core induction tool (“WCIT”) was developed to study interfacial polarization, in the 10 kHz to 300 kHz frequency range, of artificial whole cores containing either dispersed or layered phases of pyrite and/or graphite. An inversion scheme based on a Cole-Cole empirical model coupled with a TILTECC-forward model of the WCIT response was developed to quantify interfacial polarization effects on WCIT complex conductivity measurements on whole cores. The measurements on whole cores can be subsequently used to improve resistivity interpretation of subsurface conductivity and dielectric measurements in shale formations. A mechanistic model of charge accumulation and relaxation around the electrically-active phases uniformly dispersed in a medium was also developed.
Modified Petrophysical Model
Various mathematical models of accumulation and relaxation of space charges around interfaces have been developed for applications in colloidal, electrochemical, microfluidic, and geophysical sciences. A phase boundary has the tendency to acquire a surface charge due to ionic adsorption, protonation/deprotonation, dissociation of a surface chemical group, or external applied field. A net surface charge may lead to a counterion cloud and diffused-charge distribution around the interface. In the presence of an electric field, the equilibrium charge distribution around the interface may be perturbed. Conduction and diffusion processes may occur over a certain time (e.g., the relaxation time) to attain a new equilibrium state. The dynamics and distribution of charge carriers around the interface may influence the electromagnetic response of a geological mixture. Consequently, electromigration, electrodiffusion, and electrochemical processes occurring at the phase boundaries may be quantified for accurate resistivity interpretation in formations susceptible to large interfacial polarization phenomena. A mechanistic model of interfacial polarization of a geological mixture containing electrochemically-inactive phase boundaries has been developed. Therefore, in the unperturbed condition, the model disclosed herein assumes a negligible surface charge, counterion cloud, and diffused charge distribution around the interface. In the absence of redox-active ions, pyrite-bearing sedimentary rocks and graphite-bearing thermogenic source rocks may demonstrate interfacial polarization phenomena. The new model developed for petrophysical application is based on Poisson-Nernst-Planck (“PNP”) equations for a diluted solution in a weak electrical field regime. PNP equations may be more realistic compared to circuit models that assume invariant uniform bulk concentration of charge carriers. Using PNP equations, the dipole moment due to a single, isolated dispersed phase particle may be determined. A consistent effective medium formulation is employed to determine effective complex conductivity of the geological mixture. This new model may be valid for a dispersed phase that has a representative length<1 mm, conductivity<100 S/m, relative permittivity<10, and volume fraction<0.3. Analytical modeling results are demonstrated below.
Forward Modeling
An inversion scheme based on a Cole-Cole empirical model coupled with a TILTECC-forward model of a WCIT response to quantify EIP effects on WCIT complex conductivity measurements on whole cores was developed. These can be subsequently used to improve resistivity interpretations of subsurface conductivity and dielectric measurements in shale formations. This inversion scheme based on the interfacial polarization model, when coupled with a consistent forward model, may remove interfacial polarization effects from the tool's EM response and thus improve the resistivity interpretation. The level of organic maturity can be assessed by selectively estimating the volume fraction of graphite and its precursors using this inversion scheme. The proposed method facilitates the use of broadband conductivity tensor tools for geophysical applications. Without this method, geophysical applications of multifrequency conductivity tool measurements may be less reliable. This method may help completion engineers to design efficient completion plans that result in targeted fracturing of sweet spots.
The site of polarization in rocks implies the existence of space charge inhomogeneties and electrochemical potentials. The space charge inhomogeneties may cause a gradient of electrical field, and the electrochemical potentials result in local sources of electromotive forces. Colossal dielectric constant phenomena was observed with a constant low frequency value. There may be a loss peak in the imaginary part of the permittivity. Petrophysical interpretation using Cole-Cole empirical model is based on an understanding the relationship of EM response and the governing petrophysical parameters. Factors in the chargeability of earth material may include:
The grain size of particles in the sample;
The types of minerals present in matrix;
The distribution pattern of these minerals veins, disseminated, massive;
The type and mobility of ions within the pore fluids;
The microscopic interaction of solid surfaces and fluids;
The amount of surface area within specific volume; and
The total pore space and geometrical distribution of pores.
The time constant is relaxation time. The time constant may be directly proportional to square of distance between mineral grains. The frequency dispersion factor is the decay coefficient.
Qualitative Interpretation
Anomalies in thin beds and depth misalignment issues may be removed using a tensor tool and using curves with the same vertical resolution. Poor borehole conditions, hydrocarbon reservoirs, low-porosity intervals, volcanic, and evaporites and cemented formations can produce anomalous log (R/X)f values.
Quantitative Deliverables
A standalone continuous estimation of LOM of shale gas formations can be generated, which can be further interpreted with other supporting logs to obtain continuous TOC values. A distinct relationship of R-signal and X-signal frequency dispersions of pyrite and graphite allows a quick-look insight into the LOM of the shale formation, as the LOM is directly proportional to graphite concentration.
Core measurements may be used in the same well or field to correlate the volume fraction of graphite and carbon nanostructures to the LOM log of the shale rock. Laboratory measurements may be used to establish the relationship of graphite and carbon nanostructure concentration to the LOM log.
Graphite and carbon nanostructure volume fractions at each depth can be converted to an LOM log. The LOM log can be validated using S2 content obtained from pyrolysis or vitrinite reflectance of core samples. LOM log can also be validated using values of HI from the log response.
The alteration of EM response due to the volume fraction of electrically-conductive phases other than water may be quantified and used to correct the EM response to obtain accurate resistivity values (Rcorr.). The TOC calculation can now be performed using the Passey et al. TOC equation:
TOC=(Δ log Rcorr.)×10(2.297-0.1688×LOM
This TOC can be validated using measurements of pulsed neutron logging devices or from correlations using uranium logs or density logs.
Qualitative Deliverables
Complex conductivity dispersion can be used for rock typing in a shale gas formation. Experimental correspondence between the proposed electrical rock types and their complex conductivity behavior was demonstrated. Based on the laboratory results and modeling work, the quick-look shale gas electrical rock typing scheme plotted in
This quick-look electrical rock-typing method uses complex conductivity measurements, which are obtained at at least two different frequencies. One frequency may be 10 times greater than the other frequency.
A base line condition exists when the R-signal and (−)X-signal curves or the estimated resistivity and relative permittivity values for the frequencies overlay each other. This baseline may indicate non-source rock. The physics of complex conductivity measurements and polarization in non-source water-bearing formation may dictate the overlap of the curves.
The R- and (−)X-signals may be mathematically scaled for the purpose of interpretability so that the curves intersect in a clean water bearing rock. With the baseline established, in agreement with the physics of the effects of interfacial polarization, the relative separation and non-parallelism of the multifrequency resistivity and permittivity curves (or R- and (−)X-signal) can be qualitatively interpreted as organic-rich or organic-lean, as mature or non-mature, as water-bearing or hydrocarbon-bearing, and as pyrite-bearing or graphite-bearing or clean.
The log (R/X)f curve is a function of graphitic and pyritic composition, which is related to the maturity of the TOC based on the pyritization and graphitization of the source rock. So, log (R/X)f may be a subsurface indicator of Level of Organic Metamorphism (“LOM”). Log (R/X)f=log 10(R/Rbaseline)−log 10(X/Xbaseline) at a given frequency.
Using the LOM determination from a core analysis and a LOM-log (R/X)f relationship, variability of the LOM in subsurface can be characterized. A LOM of 7 corresponds to an onset of maturity, and a LOM of 12 corresponds to an onset of over-maturity.
Each separate interval of interest may have a distinct baseline for the R- and X-signals. A GR-cutoff may distinguish between organic-rich source rocks and hydrocarbon-bearing reservoir intervals for the TOC profile estimation. Baselines of R-signals at multiple frequencies may be used to remove porosity dependence. Once a baseline is established for the R-signals in a given lithology, the porosity variation may affect the responses of the R-signal curves such that a shift in one curve is reflected in a shift of comparable proportional magnitude in the other curve(s).
An appropriate scaling of the R-signal curves can remove porosity dependence, allowing for the calculation of LOM in wells. The ability to perform the log (R/X)f technique at a well site using properly scaled R-signal and X-signal curves provides an excellent method for identifying organic-rich intervals to be sampled with a side-wall coring tool.
A processor may include a microprocessor, microcontroller, processor module or subsystem, programmable integrated circuit, programmable gate array, or another control or computing device.
The storage media 1406 may be implemented as one or more computer-readable or machine-readable storage media. Note that while in the example embodiment of
In some embodiments, the computing system 1400 contains one or more LOM determination module(s) 1408. The LOM determination module 1408 may be used to perform at least a portion of the method 1300 herein.
It should be appreciated that computing system 1400 is only one example of a computing system, and that computing system 1400 may have more or fewer components than shown, may combine additional components not depicted in the example embodiment of
Further, the at least a portion of the processing methods described herein may be implemented by running one or more functional modules in information processing apparatus such as general purpose processors or application specific chips, such as ASICs, FPGAs, PLDs, or other appropriate devices. These modules, combinations of these modules, and/or their combination with general hardware are all included within the scope of protection of the invention.
While the claimed subject matter has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the scope of the claimed subject matter as disclosed herein. Accordingly, the scope of the claimed subject matter should be limited only by the attached claims.
This application claims the benefit of a related U.S. Provisional Patent Application having Ser. No. 62/081,812, filed Nov. 19, 2014, titled “Subsurface Estimation of Level of Organic Maturity Using Multifrequency Complex Conductivity Measurements In MHz to MHz Frequency Range,” to John Rasmus et al., the disclosure of which is incorporated by reference herein in its entirety.
| Number | Date | Country | |
|---|---|---|---|
| 62081812 | Nov 2014 | US |
