Woodford Case Study
Woodford Case Study
Woodford Case Study
Tom Arnold1, Andrew Sneddon1, Sheng Wu2, Andrei Deev2 , Yongchun Tang 2,3
1.
Paladin Geological Services, 13832 Santa Fe Crossing Dr., Edmond Ok 73013;
tom.arnold@pslogging.com;andrew.sneddon@pslogging.com
2.
Arrow Grand Corporation, 738 Arrow Grand Cir., Covina, Ca 91722;
sheng@arrowgrand.com; deev@arrowgrand.com
3.
Power Environmental Energy Research Institute, 738 Arrow Grand Cir., Covina, CA 91722;
tang@geoisochem.com
Abstract
The benefits to hydrocarbon exploration from geochemical and SEM analysis have been known
for many years. Yet the advent of smaller, cheaper and more specialized equipment is now
providing the ability to make highly precise and critical measurements during the drilling process
rather than having to wait for weeks or even months to obtain the results. When this data is
coupled with standard surface logging techniques of gas and formation analysis, they provide
valuable information about the reservoir that will help with current production of the subject well
and future development planning.
For this case study, several instruments were used in a mobile laboratory during the drilling of a
Woodford horizontal well on cutting samples taken at regular intervals.
Rock Eval
GCIR II for Gas Isotopes (proprietary)
Niton XRF
Hitachi TM-3030 SEM
The Rock Eval Pyrolysis provided data on producible hydrocarbons, thermal maturity and total
organic content in each sample. The GCIR II analyzed the gas isotopes from both isotubes and
isojars yielding data on c1, c2 and c3 isotopes. This data provided insight into permeability, porosity
and provided information on hydrocarbon compartmentalization and differentiation from OBM.
XRF analysis was able to quantify the elemental composition of the cutting samples, providing an
understanding of the best locations to frac based on silica content, changes in mineral composition
and an indication of TOC. Petrographic analysis was done using the SEM. This instrument gave
direct observation of both microstructures and available pore space. In turn this data gave
corroboration to the previous geochemical and elemental analysis.
Introduction
Identifying the zone of highest production potential, especially while drilling, is a critical issue.
Even with good initial understanding of the target zone, geosteering the lateral in order to remain
within that interval is often challenging due to inadequate mineralogy, geochemical or petrologic
data. With the advent of new techniques, instrumentation and near real-time analysis, the
uncertainty that often arises can be greatly reduced.
By virtue of a unique and proprietary analysis of drill cutting’s for gas isotope analysis of C1,2,3,4
and CO2 over time provides a dynamic measurement of porosity, permeability, and overall
production potential of a given zone. This procedure provides important data concerning the
identification of compartmentalization from immediate measurements on fresh cuttings. This
process incorporates data from the analysis of both isotubes, unaltered gas samples, and isojars,
unaltered mud samples, each providing data that completes the picture of the hydrocarbon
potential from a given sample.
Another important advantage of this instrumentation and its associated procedures is its ability
to differentiate OBM from potential oil reserves. This reduces the possibility of sample
contamination interfering with results.
For a direct measurement of hydrocarbon potential from the rock, we incorporate the Rock
Eval machine. This pyrolysis instrument drives hydrocarbons from the drill cuttings at specific
stages based on temperature. Typical data such as S1, the amount of free hydrocarbon present in a
rock, and Total Organic Carbon can be measured directly. This provides an excellent measure of the
‘sweet spot’ while drilling.
X-Ray Fluorescence is another technique that will provide quantitative data on mineral
composition, rock hardness and hydrocarbon potential. Providing corroboration to the geochemical
analysis, the presence of vanadium in a sample correlates directly to TOC as seen by the Rock Eval
Pyrolysis. Here are a few additional valuable measurements that can be found from the
XRF(Weedmark et al.).
Porosity is always a key element. Being able to directly measure it in a rock provides the linking
element for all the previous analytical tools. The Scanning Electron Microscope is the instrument of
choice for this. The SEM provides a means to measure and examine the pore space present. When
taken in real-time while drilling, the data provides a quick means to assist in locating the ‘sweetest
spot’ in a well. If we can determine from geochemistry that a zone has a given potential, the SEM
will identify if there is sufficient porosity present to make the zone a potential target.
Application of these instruments and procedures were provided on a horizontal Woodford play
well. This formation is well known for its hydrocarbon potential throughout Oklahoma. Although a
great deal of study has been given to the geologic setting and overall extent of the formation, there
is still much to be learned especially when looking at the formation from the perspective of a single
well.
Analysis
Samples were taken at regular intervals throughout the drilling process and subjected to a
specific regimen of preparation. For the Rock Eval, XRF and SEM, the following procedure was
followed:
1. Extracted and rinsed in Diesel (OBM), H2O (water-based) stored in plastic jars
2. Rinsed again with soap/water mixture (OBM), and water (water-based), removal of
diesel and any remaining LCM
3. Put into drying pans to dry
4. Weight samples:
-Rock Eval: .100 ± .05 mg into Crucibles
-XRF: 7.00 ± .50 mg into container cups
Sample preparation for the GCIR Gas Isotope analyzer was different.
1. Samples were collected at specific depth into Isojars, filled with water +10 drops
Zephiran Chloride to prevent any bacteria
2. Specific injection volumes carefully measured and inserted in GCIR for analysis
In order to ensure the quality of the data provided, it is necessary to follow a consistent and
regulated process for each sample analyzed.
Rock Eval Pyrolysis
The data from the Rock Eval was used to identify specific areas where the best recoverable
hydrocarbon potential was present. These measurements are based on the results provided by
heating the rock sample during several stages of pyrolysis (McCarthy et al., 2011).
Stage 3 (S2): 300 to 600 C⁰ - these are residual hydrocarbons and represent
the amount of hydrocarbon a source rock might produce if maturation
continues
Stage 4 (S3): 300 to 600 C⁰- CO2 from the thermal cracking of kerogen
Note: All depths are in Measured Depth due to directional data being unavailable.
Figure 1 shows the results of the analysis for the entire well plotted against the gamma ray and
total gas. Results from S0, S1, S1/S2 ratio and TOC are displayed versus depth. The TOC is important
because hydrocarbon production is related to its carbon content (McCarthy et al., 2011). The
sensitivity of the S1/S2 ratio to hydrocarbons and TOC provide the primary indicators for these
zones of interest. An additional indicator is seen in Figure 2 where the data for the presence of
vanadium, gleaned from the XRF analysis, is plotted versus the S1/S2 ratio and TOC. This was used
to corroborate the results from the Rock Eval. (Examples of this correlation are circled on the Figure
2.)
Kerogen is a key indicator of hydrocarbon type and maturity. By determining TMax, maximum
temperature for the generation of hydrocarbons (McCarthy et al., 2011), from the Rock Eval and
plotting it versus TOC, we are able to determine Kerogen type (Dellisanti et al, 2010). This leads to
an understanding of thermal maturity and the type of hydrocarbon present within that interval.
Figure 3 plots all the values for TMax measured throughout the lateral. From this plot we can
clearly see those depths that are the most mature and have the greatest potential to yield
productive hydrocarbons. In addition, this plot defines those zones that will produce oil, gas and
condensate based on the Hydrocarbon Index.
Figure 4 plots the values of S1 versus TOC. A linear relationship is expected due to the fact that
all these samples were taken within the lateral of a horizontal well. Had this been a vertical well,
the plot would have appeared much different.
Figure 5 plots S2 versus TOC. This defines the oil prone type II & gas prone type III kerogen
content. This is typically the result of a marine depositional environment. Figure 6 defines the
Production Index. This is derived from the pyrolysis of the first and second stage and defined as
S1/(S1 + S2) (McCarthy et al., 2011). This ratio relates to the evolution of organic matter over time
being subjected to increased temperature and pressure. Consequently it relates to source rock
maturity. From this data the oil generation potential is defined. In addition, those values that fall
outside a given set of standards can be defined as contamination.
Figure 7 plots TOC versus depth. Understanding the amount of oil and gas potential within a
source rock is critical in its evaluation. This figure indicates the progression of carbon by weight
percent as the lateral is drilled. Clearly the concentration of carbon content increases as the
wellbore penetrates deeper into the Woodford. In addition this could indicate different
hydrocarbon compartments within the source rock.
Gas Isotopes
The examination of mud gas isotopes adds a much greater understanding of the source rock.
Here is a list of some of the characteristics that can be learned from this analysis.
(1) Permeability assessment: Methane carbon isotope data for both mud gas and gas released
from cuttings from the same depth were examined. Built on the work of Xia and Tang (2012), it
was found that the isotope difference between these two types of gases may indicate the
degree of gas recovery and thus the permeability of the formation (Arnold, et al. 2014).
(2) Well compartmentalization: high-resolution isotope data for C1, C2 and C3 are good
indicators for well compartment (Arnold, et al. 2014).
(3) “Isotope reversal” indicates high production rates. With maturity increasing, gas isotope
ratios increase, that is they become less negative. However, “isotope reversal” is found in some
gas fields, where C2 and C3 isotope ratios stop increasing with maturation, but begin to decline
drastically at certain maturity. Field data indicates that isotope reversal of C2 and C3 gas
isotopes may indicate the high production rates for such wells (Arnold, et al. 2014).
(4) Delineation of gas sources: gases from different sources have their own fingerprints (shale,
coal, biogenic, thermogenic, primary gas, secondary gas, etc). The accurate gas source
correlation is necessary for a better understanding of gas in place, quality and producibility
(Arnold, et al. 2014).
(5) Gas maturity and formation temperature: This is a classic use of isotopes. With increasing
maturity, gas isotopes increase (for example, methane δ13C may increase from -60 per mil to -
30 per mil with maturation (Arnold, et al. 2014); ( Luwen, 1983; Schoell, 1983).
(6) With gas generation from shale, the gas isotope ratios follow a certain trend. Therefore gas
accumulated at different time intervals should show different gas isotope ratios. We can thus
derive the “cumulative”, ratios seen over time, vs. the “instantaneous”, ratio seen at a single
given interval. This data can then be used to qualify porosity and permeability within a given
zone based on the length of time and amount of gas released. (Proprietary data GOR-Isotopes2
Software Group, HeoIsoChem Corporation) (Arnold, et al. 2014),
(7) Nano-porosity and permeability assessment: A proprietary method has been established to
evaluate the nano-porosity and permeability based on the isotope fractionation through gas
releasing from the rock matrix—the stage-wise degassing. The basic idea is to reconstruct the
isotope fractionation variation trend through degassing, which is related to the “easiness” for
the gas moving out of the rock matrix. The stage-wise degassing procedure can be joined with
the GCIR II to obtain real-time data on site (Arnold et al. 2014).
(8) Production decline: A proprietary and novel method has been developed to predict the
production decline in a shale gas well. The scientific fundamentals are similar to those of
permeability assessment, as gas isotope changes through production (Xia and Tang, 2012, GCA).
The method depends on the instrument we construct and develop (The natural gas isotope
analyzer, NGIA), and has been successfully deployed for field trial in Barnett Shale (Arnold, et al.
2014).
Figure 8 displays all isotube gas isotopes: iC1, iC2, iC3. The decreased amount of fluctuation
of the iC3 data is the result of a much heavier molecule and closely relates to the stability of the GCIR II.
The changes visible in iC1 over the depth range indicate hydrocarbons from different compartments
within the source rock.
Figure 9 is a plot of isotubes versus isojars for iC3. This analysis provides an indication of
areas of maximum porosity. 14,600’ MD to 16,500’ MD appear to be the areas that are of greatest
promise.
In Figure 10 iC2 is plotted for both isotubes and isojars showing a smaller difference.
Figure 11 puts everything together for all the isotopes derived from both tubes and jars. A
large difference between Jars/Tubes and their dynamic change, over time, provides information on
porosity/permeability, underscored by the Jars’ concentration change. As seen previously, the 14,600’
MD to 16,500 MD zone appears to the area of greatest hydrocarbon production. All this data indicates
the potential for hydrocarbon recovery is high and especially favors liquids.
X-Ray Florescence
Figures 12 & 13 plot the elemental composition of all the rock samples collected. The increased
presence of silica is defined within several areas, marking zones where fracking would be much easier.
These areas range between 12,500’MD to about 14,000’MD. From 12,400’MD to 13,800’MD the values
of sulfur and iron are very high. This is would indicate an area with considerable pyrite which could
inhibit production. (A closer investigation with the SEM will be necessary to identify the degree of
degradation due to pyrite pore blockage.) The XRF analysis of vanadium indicates possible zones of
production that fall well within those areas prescribed from isotope analysis and Rock Eval. This level of
cross corroboration is important in developing a clear picture of those areas that are of particular
interest.
From all this data, we have found three zones that require greater investigation.
Figure 14 is an SEM image taken of a sample originating at 13,070’MD. The pyrite framboid seen
here is a good indication of the increased sulfur and iron content defined by the XRD. This image is at
2500X and represents an area of 30 micormeters across. The only porosity visible in this photo is intra-
granular within the pyrite. Production from this area would be challenging. The spectrogram seen in
Figure 15 shows the increased iron content in this sample, most likely from the pyrite present.
In Figure 16 taken from 13,430’MD, the porosity becomes clear as the amount of iron and sulfur
decrease. Here the porosity has increased to very good. The average pore size is 1.81 micrometers.
From the Rock Eval data we find that S1 has increased through this area indicating more free
hydrocarbons. Notice the increased calcium in the spectrogram seen in Figure 17. Obtained from the
XRF, this indicates the possibility of CaCO3 may be blocking pore throats or even present as calcite.
As an example of the second area of interest, Figure 18 taken from 15,050’ MD shows
considerable porosity. Numerous pores are present with an average size of 6.7 micrometers. (This is the
best visible porosity seen within the well.) Looking at the Rock Eval for S1 on Figure 4 at this depth, a
marked increase can be seen. The XRF spectrogram in Figure 19 shows significant decrease in both
calcium and iron thus removing the CaC03 and pyrite concern.
Our last area of concentration is represented by a sample taken at 16,130’ MD seen in Figure
20. There are an increasing number of large pores present with an average size of 2.43 micrometers.
Porosity appears very good.
When all the data is tabulated from the SEM analysis, we have provided two plots which
represent the porosity within the lateral. The first plot, seen in Figure 21, plots average pore size versus
depth. The second plot indicates the overall porosity quality. Figure 22 plots the quality estimate versus
depth. This value is a subjective analysis based on the ease of locating porosity in a sample, the size of
the pores found and number of pores present in the image. All this is analyzed under standard
magnification within a given area.
Conclusion
The combination of four techniques has provided an in-depth understanding of the hydrocarbon
potential of a Woodford well.
• GCIR tells C1-5+CO2 isotope & composition static (tube) & dynamics (cuttings in jars), relates to
gas & permeability
• Rock Eval tells S0,1,2, TOC (liquid and potential) of the cuttings at given depth
Putting the mosaic of data together as we have presented here will help geologists to decide where
the hydrocarbons are located, which zones provide the greatest potential for production, where to frack
the well and reveal information about the reservoir for future development. This type of information
will allow producers to provide the maximum production from a well and reduce cost from unnecessary
drilling practices.
Biographical Sketch
Tom Arnold
His major research interests are applying molecular modeling and experimental simulation
techniques for geological problems. He pioneered the molecular modeling technique to many fields
of organic geochemistry, surface chemistry, reaction kinetics and organic and inorganic reactions
under subsurface geological processes, real-time gas isotope technology, double substitute
isotopes and in-situ measurement of CO2 and nano-porosity of rocks.
He served as consultant for many energy firms including Chevron, Shell, Exxon, BP,
Conocophillips, Petrobras (Brazil), Saudi Aramco, ENI-Agip (Italy), Total (France), Devon, Petrochina
and China Petroleum Corporation (Taiwan). Hs is cofounder for several startup companies in the
USA and abroad.
Sheng Wu
Andrei Deev
VP of Research and Development, ArrowGrand Technologies
727 Arrow Grand Circle, Covina, CA 91722
Andrew Sneddon
References Cited
Arnold,T., Gao, L., Tang, Y, Advanced Techniques to Increase Production from Horizontal Wells and
Reservoirs, The Journal of the Oklahoma City Geological Society (Shale Shaker), July/August
2014, pp278-303
McCarthy, K., Niemann, M., Palmowski, D., Peters, K., Stankiewicz, A., (2011), Basic Petroleum
Geochemistry for Source Rock Evaluation, Oilfield Review Summer 2011: 23 no 2.
Dellisanti, F., Pini,G.A., Baudin, F., (2010), Use of Tmax As A Thermal Maturity Indicator in Organic
Successions and Comparison with Clay Mineral Evolution, Clay Minerals 2010: 45, pp115-130.
Weedmark, T.C., Spncer, R.J., Juss, D.P., Marechal, F., (2014) Applications of Portable XRF, Chemical
Stratigraphy and SEM in the Horn River Basin, CSOG/CSEG/CWLS GeoConvention 2012, Calgary
presentation.
Lewan, M.D. (1983) Effects of thermal maturation on stable organic carbon isotopes as determined by
hydrous pyrolysis of Woodford Shale, Geochimica et Cosmochimica Acta, 47(8), Pages 1471-
1479, ISSN 0016-7037, 10.1016/0016-7037(83)90306-X.
Schoell, M. 1983. Genetic characterization of natural gases, AAPG Bulletin, , v. 67, p. 2225-2238
Wu, S., Deev, A, Haught, M, Tang, Y. 2008. Hollow Waveguide Quantum Cascade Laser Spectrometer as
an Online Microliter Sensor for Gas Chromatography. Journal of Chromatography A, 1188(2): p.
327-330.
FIGURES
WOODFORD
Figure 1: Rock Eval Pyrolysis for the subject well. The top of the Woodford is at 12,120. The spikes in
the S1, S1/S2 and TOC are significant and relate to specific hydrocarbon zones within the lateral.
Figure 2: Vanadium is a good indicator of the presence of hydrocarbons. Here it is
plotted against the hydrocarbon sensitive S1/S2 ratio and TOC.
1000
Immature Mature Postmature
900
Type I Kerogen Oil
Oil Prone Window
800
indicates
Gas Zone
Condensate- Wet
HYDROGEN INDEX (HI, mg HC/g TOC
lacustrian
700 9500-12080
Type II Kerogen 12090-13370
600 Oil Prone-->
marine 13400-14900
500 14930-16580
16610-17240
400 Type II-III
Kerogen
Oil- Gas Prone
300
Figure 3: TMax vs Hydrogen Index provides identification of hydrocarbon maturity and the Kerogen
type.
TOC vs S1
40
35
30
S1 (mg HC/gRock)
25
9500-12080
20 12090-13370
13400-14900
15
14930-16580
10 16610-17240
0
0 1 2 3 4 5 6
TOC (wt %)
Figure 4: TOC vs S1. Note the linear relationship of the data as samples were collected throughout
the lateral.
TOC vs S2
25
15
9500-12080
12090-13370
13400-14900
10
14930-16580
16610-17240
0
0 2 4 6 8 10
TOC (wt%)
Figure 5: TOC vs S2. Seen here the majority of the production is within the oil/gas zone for Type II & II
Kerogen.
Tmax vs PI
2
1.8
Oil
Generation
1.6
1.4
Contamination
Production Index (PI)
9500-12080
1.2
12090-13370
1 13400-14900
14930-16580
0.8
16610-17240
0.6 Immature
0.2
0
250.0 300.0 350.0 400.0 450.0 500.0 550.0
Tmax (0C)
Figure 6: Tmax vs PI. This relates the maturity of the hydrocarbon to the type of production expected.
TOC vs Depth
TOC (wt %)
0 1 2 3 4 5 6
9000
10000
11000
9500-12080
12000
12090-13370
13400-14900
Depth (ft)
13000
14930-16580
14000 16610-17240
15000
16000
17000
18000
10000
iC1 iC2 iC3
11000
12000
13000
14000
15000
16000
17000
Figure 8: Plotted vs depth, the values for each isopote is provided. The fluctuations in iC1 indicates
hydrocarbons from different compartments within the Woodford.
Isotope Concentration Change
11000
iSotubes 12000
13000
iSojars
14000
15000
Most Productive Zone
16000
17000
Figure 9: Between the relationship of iC3 from isotubes and isojars, the most productive zones within
the well can be defined.
Isotope Concentration Change
11000
12000
13000
14000
15000
16000
17000
Figure 10: Plotting iC2 from data collected between iSotubes and iSojars, a smaller difference is
found than the plot of iC3 seen in Figure 9.
Isotope Concentration Change
11000
12000
13000
Jar…
14000
15000
16000
17000
Figure 11: When taken collectively, we find a large difference between isotubes and isojars. This data
indicates that the potential for hydrocarbon recovery is very high and tending toward liquids. A large
difference, green double arrow, between Jars /Tubes and their dynamic change over time provides
information on porosity/permeability, underscored by the Jars’ concentration change.
Figure 12: XRF measurements in the lateral within the Woodford. The relationship of sulfur to iron is an indication of
pyrite, while calcium helps identify calcification at a given interval. Potassium, K, is one of the radioactive elements that
make up the spectral gamma ray, the other two being Uranium and Thorium.
Figure 13: XRF measurements in the lateral within the Woodford, continued. The quantity of Chromium and
Molybdenum both relate to TOC in a rock along with Vanadium.
Figure 14: SEM image of a pyrite framboid from a sample depth of 13,070 MD. The magnification was 2500X.
XRF Spectrogram from 13,070 MD
Figure 15: XRF spectrogram taken at 13,070 MD. The peaks of iron are higher through this interval indicating the
increased potential for the presence of pyrite.
Figure 16: Now at 6000X, porosity becomes visible in this image from 13,430 MD.
XRF Spectrogram from 13,430 MD
Figure 17: XRF Spectrogram from the sample at 13,430 MD. The amount of iron has significantly decreased from the
sample at seen in Figure 16. Note the increase in Ca in this sample.
Figure 18: SEM image from 15,050 MD. The porosity has increased in this sample with an average pore diameter of 6.7
micrometers.
XRF Spectrogram from 15,050 MD
Figure 19: This XRF spectrogram shows a decrease in both iron and calcium indicating a good hydrocarbon recovery
potential.
Figure 20: SEM image from 16,130 MD. There are a good number of pore present with an average size of
2.43 micrometers. The porosity appears very good. The magnification in this image was 4000X.
Pore Size µm e S
0.1 2.1 4.1 6.1 8.1 10.1
11000
11170
12000
14000
15000 14930
1511015080 15050
15140
1520015170
15260
15320
15380
15470
15500
16000
1610016040 16130
16220
1628016310
Figure 21: From the SEM images, average pore sizes are plotted.
Porosity Estimate (1-5)e S
0.1 0.6 1.1 1.6 2.1 2.6 3.1 3.6 4.1 4.6
11000
11170
12000
14000
15000 14930
15080
15110 15050
15140
15200 15170
15260
15320
15380
15470 15500
16000 16040
16100 16130
16220 16280
16310
Figure 22: Based on evaluation of the SEM images pertaining to the ease of locating porosity, average
pore size and the number of pores present, a subjective measure of porosity was obtained.