Skip to main content
  • Research article
  • Open access
  • Published:

Methane and carbon at equilibrium in source rocks

Abstract

Methane in source rocks may not exist exclusively as free gas. It could exist in equilibrium with carbon and higher hydrocarbons: CH4 + C < = > Hydrocarbon. Three lines of evidence support this possibility. 1) Shales ingest gas in amounts and selectivities consistent with gas-carbon equilibrium. There is a 50% increase in solid hydrocarbon mass when Fayetteville Shale is exposed to methane (450 psi) under moderate conditions (100°C): Rock-Eval S2 (mg g-1) 8.5 = > 12.5. All light hydrocarbons are ingested, but with high selectivity, consistent with competitive addition to receptor sites in a growing polymer. Mowry Shale ingests butane vigorously from argon, for example, but not from methane under the same conditions. 2) Production data for a well producing from Fayetteville Shale declines along the theoretical curve for withdrawing gas from higher hydrocarbons in equilibrium with carbon. 3) A new general gas-solid equilibrium model accounts for natural gas at thermodynamic equilibrium, and C6-C7 hydrocarbons constrained to invariant compositions. The results make a strong case for methane in equilibrium with carbon and higher hydrocarbons. If correct, the higher hydrocarbons in source rocks are gas reservoirs, raising the possibility of substantially more gas in shales than analytically apparent, and far more gas in shale deposits than currently recognized.

Introduction

Few questions in geoscience are as interesting and controversial as the origin of methane in natural gas. And fewer go as far back in time. Evans proposed thermal cracking in 1971 to explain oil trending to methane with depth in a Canadian basin [1]. It was assumed that hydrocarbons cracked more or less randomly to smaller hydrocarbon and ultimately to methane. This premise had broad appeal, but no empirical or theoretical foundation. McNab had attempted to replicate methane generation from oil cracking in long-term cracking experiments in 1952, but failed [2]. Thermal cracking theory was nevertheless elevated to text-book status [3, 4], but it was never to gain the empirical support McNab sought, and its predictive powers today are limited. It cannot, for example, explain the properties of natural gas, the thermal stability of its light hydrocarbons [5], their compositions in gas deposits [6], or their constraints to thermodynamic equilibrium [7].

Catalysis by transition metals explains these properties [57], and there is considerable experimental evidence supporting it. Source rocks release gas catalytically in laboratory experiments at ambient temperatures [810], and molecular probe experiments confirm natural catalytic activity as the source of gas in these experiments [11]. However, there is also experimental evidence against transition metal catalysis in natural gas generation [12].

Catalysis accounts for the composition of natural gas, but other questions remain troubling. Shales ingest and expel gas in the laboratory in ways contradicting classical physics. Non-classical behavior is apparent in chaotic generation curves [8], and in the generation of hydrocarbons in reverse-mass order [9]. It is also apparent in the molecular probe experiments [11]. Methane ingestion in Fayetteville and Mowry Shales reported here is another outstanding example.

The premise that methane should be independent of the solid hydrocarbons generating it could be flawed. If generation is catalytic and reversible, methane and its source would be more like carbon dioxide and calcite at equilibrium. Carbon dioxide is never independent of calcite. Ingested and expelled reversibly from CaCO3/CaO, its pressure increases and falls with temperature as a thermodynamic function of the equilibrium constant. Methane could similarly exist in two states, gaseous methane and solid-state methane, at equilibrium. It accounts for the results presented here and the non-classical behavior in earlier experiments.

An equilibrium hypothesis is presented in three parts. The general concept is laid out first (Theory). The second part is the supporting experimental evidence (Ingestion), and the third contains a kinetic model of the theoretical equilibrium (Gas-Solid Equilibrium Kinetic Simulations). A new general equilibrium model for light hydrocarbon generation summarizes the results.

Results and discussion

Theory - Solid hydrocarbons storing and generating gas is hypothetical. It considers two very different reservoirs of gas in shales. The first is open porosity storing and expelling gas as a passive container in accordance with classical theory. The second does not. It is comprised of solids like pyrobitumens, and similar high-molecular weight hydrocarbons found in organic-rich shales. It is not a passive container storing hydrocarbons in solution, however. Methane could only have a very limited solubility in solid hydrocarbons, perhaps no more than a few percent by my estimates. In the current hypothesis solid hydrocarbon reservoirs are co-polymers of methane and carbon, at equilibrium. The equilibrium in reaction (1) illustrates the concept in its simplest form.

n CH 4 + n C < = > CH 2 2 n .
(1)

The equilibrium hypothesis and supporting evidence are discussed below. Here, we consider the origin of solid hydrocarbons, also referred to as ‘carbon pools’. ‘Methane’ is used throughout to illustrate relationships for simplicity. All references to ‘methane’ or ‘CH4’ should apply to the higher light hydrocarbons as well, although methane is the primary focus here.

The existence of solid hydrocarbons in source rocks is not new. They have been recognized by organic geochemists for decades. Generating gas from higher hydrocarbons necessarily generates a carbon deposit to balance hydrogen. It is the text-book explanation for the origin of pyrobitumen, for example, a ubiquitous organic mineral in sources rocks with H/C ratios ranging from about .5 to 1.6 [4]. Pyrobitumens are typically rich in transition metals [13, 14] and have recently been cited for catalytic activity in methane generation [15]. Solid hydrocarbons have been proposed as catalysts resembling activated carbon in the decomposition of higher hydrocarbons to gas [16]. Solid hydrocarbon pools in this model can be pyrobitumen, or any other organic solid containing a catalyst that generate gas. They are both catalysts for generating methane and reservoirs for storing it.

Solid hydrocarbons are not uncommon in transition metal catalysis. The hydrogenation of carbon monoxide to methane over nickel, ‘methanation’, is an outstanding example [17]. The intermediates to methane are distributed in pools of solid hydrocarbons associated with nickel. The nature of the carbon polymer in methanation is unclear except that it is a co-catalyst, it is unsaturated in hydrogen and it is not graphite. It is like the carbon pool proposed here, some carbon polymer between saturate hydrocarbon CH2 and graphite CH0. The activated carbon suggested by Alexander et al. [16] in the catalytic generation of natural gas is another relevant example.

How the proposed equilibrium reaction might proceed mechanistically is discussed in Gas-Solid Equilibrium Kinetic Simulations, below. We are only concerned here with overall conversion. Methane generation feeds carbon into a pool of hydrogen-unsaturated carbon, into pyrobitumen, for example.

Since source rocks possess natural catalytic activity and generate catalytic gas in the laboratory [11], we shall assume that natural methane is catalytically generated and that the catalyst resides in pools of solid hydrocarbons. Because most catalytic reactions are reversible and approach equilibrium over time, we also assume that methane generation is reversible and should bring methane and solid hydrocarbons into equilibrium over time.

Reaction (2) illustrates a general equilibrium between methane and solid hydrocarbons.

x CH 4 + CH y < = > C x + 1 H y + 4 x
(2)

There is a 46% mass difference in methane equivalents between a solid hydrocarbons with compositions CH0.5 and CH1.6 , for example, a typical range for pyrobitumens [4]. Thus, CH0.5 can consume 46% of its mass in methane generating CH1.6 if methane exists in equilibrium with solid hydrocarbon.

However, there is nothing to suggest that methane actually reacts with solid hydrocarbons generating new compounds. It is not enough that methane merely dissolves in CHy. In that case, x in reaction (2) would be insignificant and the equilibrium in (2) meaningless.

The proposed equilibrium is purely hypothetical, but testable. Do solid hydrocarbons in source rocks consume extraordinary amounts of light hydrocarbons in ways distinguishable from simple adsorption and solution?

Light Hydrocarbon Ingestion- The experimental challenge is to distinguish gas addition to the solid hydrocarbons in source rocks from classical gas addition to ordinary solid hydrocarbons. There are several classical ways light hydrocarbons can add to heterogeneous materials like source rocks. They can be adsorbed on surfaces, go into liquid solutions, or into polymer solutions, for example. There is no chemical change in adsorption and solution, however. They are first-order reactions in which rates of addition are proportional to concentrations: Rate = k*(X), where (X) is the concentration of free hydrocarbon X and k is the first-order rate constant. Because hydrocarbons do not compete for surface sites or positions in solution, transfer rates and solubilities are typically independent of other hydrocarbons. Butane has about the same water solubility in helium as it has in methane, for example [18].

Assume that hydrocarbons add to the solid hydrocarbons in source rocks differently. Hydrocarbon X adds to some receptor [···] generating the adduct [X], a process referred to here as ‘ingestion’. X disappears by second-order kinetics at rates proportional to X and [···]: Rate = k*(X)*[···]. However, in contrast to adsorption and solution, hydrocarbons compete for [···]. Thus rates of ingestion are not independent of other hydrocarbons.

We can distinguish between ingestion and either adsorption or solution through competitive addition. If solid hydrocarbons are mere solvents, hydrocarbon X will go into solid-solution at rates proportional to concentrations of X, independent of some other hydrocarbon Y. If, however, X reacts with [···], and Y competes with X for [···], X will go into the solid at rates inversely proportional to Y.

Do hydrocarbons compete in their addition to solid hydrocarbons?

To put this question to experimental test, two vials were charged with identical amounts of Mowry Shale and an equal molar mixture of ethane, propane, iso-butane and n-butane (C2-C4). Both vials were heated at 75°C for 200 hours. One (Vial A) was diluted with argon by 50% five times over the 200 hours and the other (Vial B) was diluted with equal amounts of methane five times. The only difference in the two reactions was dilutions with argon in one (Vial A) and dilutions with methane in the other (Vial B). Thus, C2-C4 gas concentrations in each vial diminished equally with each argon and methane dilution. Concentrations over time should be about the same in A and B if Mowry Shale is passive and only removes C2-C4 hydrocarbons by adsorption or solution. If it removes them by ingestion, and C1-C4 hydrocarbons compete for [···], C2-C4 concentrations in Vial A should fall progressively below those in Vial B.

Another vial (C) was charged with beach sand with a thin coating of n-octadecane (~ 1%) to assess adsorption and solution. The C2-C4 hydrocarbon mixture was added and the vial was heated and diluted with methane as was Vial B.

If Mowry Shale consumes hydrocarbons classically through adsorption and solution exclusively, the ratio of sums ∑ (C2-C4)A / ∑ (C2-C4)B should remain constant over time. If, however, Mowry Shale ingests C2-C4, the ratio of slums ∑ (C2-C4)A / ∑ (C2-C4)B should progress to zero over time.

Hydrocarbon gas concentrations (C2-C4) over time in the three vials are shown in Table 1. The differences between Vials A and B over time are dramatic. Mowry Shale consumed 163 μg C2-C4 g-1 in argon (Vial A) compared to under 30 μg g-1 in methane (Vial B). Figure 1 shows the decline in n-butane gas concentrations over time in A and B. Concentrations of n-butane in B fall exponentially consistent with sequential 50% dilutions, but% composition of C2-C4 remains essentially constant over time (Table 1B). The blue line passing through the B data is not the regression line for that data. It is the dilution line indicating where the data should plot if only dilution were lowering n-butane gas concentrations. Hydrocarbon concentrations in A fall sharply below the blue dilution line progressing to zero over time signaling almost total ingestion. Moreover, n-butane was selectively withdrawn from Vial A. Percent n-butane (C2-C4) fell sharply in A and remained essentially constant in B (Figure 2). Concentrations at termination were 16 μg g-1 in B and 0.064 μg g-1 in A, a 250-fold difference.

Table 1 Hydrocarbon ingestion, Mowry Shale, 75°C, 200 hours
Figure 1
figure 1

n-Butane ingestion in Mowry Shale from argon and methane. Concentrations of hydrocarbons in vials over time (in arbitrary units) are in Table 1. Blue solid dots are gas concentrations of n-butane in methane (Vial B) and red solid dots are gas concentrations in argon (Vial A). The blue line traces the concentration of n-butane diminishing by 50% with each dilution. It is not the regression line for the Vial B data. The red line is the DED regression line for the Vial A data. Data points below the blue line reflect n-butane ingestion into Mowry Shale.

Figure 2
figure 2

% wt n-Butane in C2-C4 in Vials A and B. % Concentrations of hydrocarbons in vials over time (in arbitrary units) are in Table 1B. Red: Vial A (argon); Blue: Vial B (methane).

C2-C4 concentrations in Vial C fell over time proportional to their molecular weights consistent with solution in n-octadecane. The ratio of iso-butane/n-butane approached equilibrium (~ 3 at 75°C) consistent with acid-catalyzed isomerization promoted by the mild acidity of beach sand [19]. The equivalent reaction in Vial B showed no change in iso-butane/n-butane, or in the composition of C2-C4. Thus, Mowry Shale in B showed no evidence of isomerization activity or adsorption under the reaction conditions. We conclude from this that the solubility of hydrocarbons in Mowry Shale by solution in liquid and solid hydrocarbons is relatively insignificant under these conditions.

Mowry Shale in argon thus consumed C2-C4 hydrocarbons vigorously and selectively at 75°C over 200 hours. The same shale in methane totally rejected the same C2-C4 hydrocarbons under the same conditions.

Methane ingestion- Fayettville Shale was exposed to nine sequential cycles of methane pressurization and exhaustion at 100°C. Vessels filled with shale were pressurized (450 psi), sealed, and slowly vented. The shale released only trace amounts of C2 and higher hydrocarbons during exhaustion. There was no evidence of free methane in the shale after these experiments. Heating vented samples in Argon at 100°C produced only trace amounts of methane and higher hydrocarbons, substantially less than the original sample.

Methane was indeed consumed by Fayetteville Shale, but not as free methane. It appeared in the solid hydrocarbons, in the Rock-Eval S2 peak, not in the free hydrocarbon S1 peak (Figure 3 and Table 2). Remarkably, the shale returned to its original Rock-Eval composition when exposed to the same exhaustion experiments in 2% methane. Thus, the solid S2 hydrocarbons in Fayetteville Shale had increased in mass by about 50% on exposure to high methane partial pressures and decreased in mass by 50% when exposed to low methane partial pressures. Samples with increased solid hydrocarbon mass showed no evidence of free methane either at 300°C in Rock-Eval analysis, or at 100°C for one hour in our analysis.

Figure 3
figure 3

Rock-Eval signals for Fayetteville Shale before and after exposure to methane. Upper panel is the FID trace for the starting shale and the lower panel is the FID trace for the same shale exposed to 100% methane (Table 2).

Table 2 Rock-Eval analysis of Fayetteville Shale before and after exposure to methane 450 psi, 100°C

Fayetteville Shale thus consumes methane at 100°C, nearly doubling its S2 mass. Mowry Shale consumes n-butane selectively in preference to ethane, propane, and iso-butane. In excess methane, methane is consumed to the exclusion of butane. Shales therefore consume light hydrocarbons by competitive addition to receptors generating high molecular weight hydrocarbons consistent with the proposed equilibrium between methane, carbon, and higher hydrocarbons (Reaction 1).

Gas-Solid Equilibrium Kinetic Simulations - Consider a source rock with free gas of concentration (Cx), where Cx is a light hydrocarbon containing x carbon atoms. It is in communication with solid hydrocarbon with the capacity to generate (Cx) from some catalytic intermediate [Cx] where [···] denotes the concentration of open catalytic sites, the receptors noted above. [Cx] forms through reaction (3) from [Cn], a high molecular weight hydrocarbon bonded to an active site, for example (CH2)n-M where n > > x, an intermediate discussed elsewhere [10].

C n + [ ] < = > C x + C n x
(3)

We shall assume [Cn] and [Cn-x] in reaction 3 are indistinguishable when n is large (n ~ n-x), and that [Cn] is therefore a constant in the kinetics of [Cx] generation. Gas generation proceeds through reaction (4), a reversible reaction. [···] bonds selectively to various hydrocarbons, and only to hydrocarbons.

C x < = > C x + [ ]
(4)

A hydrogen balance is ignored throughout this scheme. All hydrocarbons in brackets are unsaturated, CnH2n, and those in parentheses are saturated, CnH2n+2. Thus, reaction (3) is hydrogen neutral. Reaction (4), however, is not, and should be ([Cx] + H2 = (Cx) + [···]) [10]. We assume for simplicity that hydrogen is in excess and that it moves easily from the pool to [Cx] generating (Cx) and from (Cx) to the pool generating [Cx].

Reactions (3) and (4) are within a closed system, a sealed source rock, for example. Gas escapes confinement through reaction (5).

C x = > C x
(5)

Reaction (6) summarizes the kinetic steps interconverting free gas and solid-state gas as symbolized by [Cx] and [Cn].

C n + [ ] < = > C x < = > C x + [ ]
(6)

The rate of [Cx] generation is proportional to [···]*[Cn]. Because [Cn] is constant and [···] becomes constant over time (steady-state), the rate of gas generation from solid hydrocarbon is constant over time. It will continue generating [Cx] and free gas (Cx) until their concentrations become sufficiently high to promote reverse reactions at rates equal to forward reactions. At that point, reaction (6) is at steady-state (equilibrium). (Cx) and [Cx] are then at their maximum concentrations, and [···] is at its minimum concentration. This is the state of a sealed (closed) source rock in the subsurface. We are interested in the dynamics of that system when the rock is opened, and reaction (5) is significant. How will (Cx) decline over time when the system is at equilibrium and the rate of reaction (5) exceeds the rate of generation (reaction 3)?

Three possibilities are considered in the following kinetic simulations. The first is classical first-order expulsion without generation or equilibrium of any kind. The kinetic model thus contains only free, in-place gas (Cx). [Cn], [Cx], and [···] are all zero. In the second, (Cx) and [Cx] interconvert at equilibrium, but [Cn] does not generate [Cx]; gas generation does not attend gas release, in other words. In the third, reaction (6) is fully operative, two reservoirs of gas, [Cx] and (Cx), at equilibrium, are depleted while [Cx] is generated at a constant rate from [Cn].

Each possibility is treated assuming equal rate constants and equivalent intermediate concentrations. Figure 4 displays a schematic of the reactions used in the kinetic simulations. A small fraction of gas in reservoir A was removed by first-order kinetics in each iteration, ∆(Cx) = (0.2*(Cx)), and all other concentrations altered accordingly as described in Figure 3. Thus, ∆(Cx) over iterations simulates expulsion rates over time. ∆(Cx) is henceforth referred to as ‘rate’ denoted R and iteration as ‘time’ (t). Figure 5 shows the decline curves, R vs. t, for three hypothetical rocks. In the first (no catalytic activity) R declines exponentially (ED) as expected, thus describing the straight line on log scale. In the equilibrium models, R declines by double exponential decay (DED) without generation and by double exponential decay with a constant (DED1) with generation, (equation 7), where Rt is the expulsion rate at time t, (R)i is the initial expulsion rate for free gas (Cx), [R]i is the initial expulsion rate for solid hydrocarbon gas [Cx], and C is the constant for rate of generation from [Cn].

R t = R i e at + R i e bt + C
(7)
Figure 4
figure 4

Schematic diagram of kinetic model simulating gas generation from gas-solid equilibrium. A and B are reservoirs charged with native gas, (Cx) is in A and [Cx] is in B (reaction 8). Reservoir C resides in B. Open catalytic sites [···] reside in C and B, unspecified in diagram. Gas passes reversibly between reservoirs as indicated by first-order kinetics (k1 = 0.07) and by second-order kinetics (k-1 = 0.007) for generation of [Cx]. Gas passes to the surface from A by first-order kinetics controlled by rate constant k2 (0.2). Catalytic gas is generated in C from [Cn] at a constant rate controlled by second-order rate constant k (0.0004), considered irreversible in this simulation. Generated gas passes directly into B. Each reservoir is filled with an initial quantity of gas, assumed non-depletable in C (100), and depletable by first-order kinetics in A ((Cx) = 200) and B ([Cx] = 200). The initial value of [···] was assumed 5, with all concentrations unitless. Reaction 8 was simulated by iterating kinetic steps by first or second-order changes until concentrations of [···], [Cx], and (Cx) were constant (steady-state). For example, the concentration of [Cx] at iteration ii was: [Cx] ii  = [Cx] i  + (k*[ ] i *[Cn]) + (k-1*(Cx) i *[···] i )-(k1*[Cx] i ). All concentrations were similarly calculated at each iteration.

Figure 5
figure 5

Kinetic simulation of gas-solid equilibrium. The kinetic scheme in Figure 3 was simulated in 50 iterations with the following restrictions. Black Curve: Reaction 8 where [Cn] = [Cx] = 0. Blue Curve : Reaction 8 where [Cn] = 0. Red Curve : Reaction 8. Starting concentrations: [···] = 5; [Cn] = 100 (constant); [Cx] = 200; (Cx) = 200 (400 in first curve). Rate constants: 0.0004: [Cn] = > [Cx]; 0.07: [Cx] = > (Cx); 0.007: (Cx) = > [Cx]; 0.2: (Cx) = > Cx. The curves represent 50 first and second-order iterations. For example, the concentration of (Cx) on iteration ii is: (Cx)ii = (Cx)i + (.07*[Cx]i) – (.007*(Cx)i) – (.2*(Cx)i). Equations for lines (t = time): black, (Cx) = 500 exp(-.223 t) blue, (Cx) = 224 exp(-.243 t) + 27.3 exp(-.029 t) red, (Cx) = 223 exp(-.243 t) + 15.5 exp(-.045 t) + 12.5.

Figure 6 shows DED1 curves with different values of C reflecting different levels of generation attending production. Therefore, rates of production from source rocks releasing gas in equilibrium (reaction 6) should decline by DED1 reflecting free gas declining exponentially, solid-state gas declining exponentially, and generated gas at a constant rate C.

Figure 6
figure 6

Effects of generation on double exponential decay curves from Gas-Solid Equilibrium Kinetic Model. The kinetic model in Figure 3 (modified units) generates the family of DED1 curves shown with zero-order rate constants for generation from k = 0 to k = 0.016 ([Cn] = > [Cx]).

Gas production Miller Heirs well, Fayetteville Shale - Many gas wells decline classically, by simple exponential decay, but there are notable exceptions mainly in low porosity-permeability wells including unconventional wells [2024]. Figure 7 shows typical non-classical decay for a well producing from a source rock, in this example the Miller Heirs well producing from Fayetteville Shale. The line passing through the data is the DED1 line calculated from the kinetic model in Figure 5. It is not a best-fit equation line to the data. Showing the theoretical curve superimposed on field data in this way illustrates the remarkable fit between theory and field data. Regressing Miller Heirs data to DED1 gives an exact fit (R2 = 0.99), while regressing the data to single exponential decay gives a substantially poorer fit, R2 = 0.84. The constant C in the DED1 regression equation was significant throughout production, accounting for 37% of produced gas at 25 months.

Figure 7
figure 7

Production decline data for Miller Heirs 1-10H well producing from Fayetteville Shale. The Millers Heirs well is in the Arkoma Basin, Arkansas, Lease # 383264348; January 2006 - August 2011. The line passing through the production data is the decline curve calculated from the kinetic model in Figure 3 with the following rate constants and initial concentrations (63 iterations): k = 0.00012; k1 = 0.1; k-1 = 0.005; k2 = 0.16; [Cn] = 100 (constant); [Cx] = 250; (Cx) = 90; [···] = 5. The rate constants and starting concentrations in Figure 3 were hand-adjusted to these values to approximate the general shape of the Miller Heirs data. The line is not a best-fit line or the DED1 regression line to the data. The respective DED1 regression equations are: Miller Heirs data (R2 = 0.99): R = 20.7 e-0.266t + 6.71 e-0.039t + 1.50. Model, Blue line: R = 16.0 e-0.197t + 6.60 e-0.0526t + 1.95.

Regressing the same data to double exponential decay without a constant (DED, eq. 7 where C = 0) gives an equally strong correlation (R2 = 0.99), however. It is therefore impossible to evaluate the constant C from the data fit to DED1 in Figure 7. If the data has genuine linearity – two sources of gas declining exponentially and one constant - then any DED equation that fits the data between time 0 and n will necessarily underestimate the data beyond n. The DED line will fall exponentially while the data approaches the constant C. Therefore, the test for linearity (C) lies beyond n, where the respective regression curves predict the future.

The Miller Heirs data were regressed between times 0 and 50 months by DED and DED1. The two regressions gave similar curves with high degrees of correlation to the data between these time limits (R2 ~ 0.99). However, only the DED1 equation predicts the Miller Heirs data beyond 50 months (Figure 8). The DED line fall sharply below the Miller Heirs data. It pojects exponential decline while the data describes a largely linear rate of decline beyond 50 months. Figure 8 leaves little doubt about the dimensions of decline in this well. There are clearly three, and DED1 (eq. 9) describes them very well. DED1 is thus a property of Miller Heirs production data and the gas-solid hydrocarbon equilibrium model as Figure 7 so clearly illustrates.

Figure 8
figure 8

Production data for Millers Heirs wells beyond 50 months with regression lines for DED and DED1 between 0 and 50 months. The data in Figure 6 were regressed to double exponential decay with a constant C (DED1) and without a constant (C = 0) (DED) between 0 and 50 months. Red line: R = 21.4 e-0.25t + 7.0 e-0.022t. Blue line: R = 20.8 e-0.263t + 6.68 e-0.0364t + 1.33.

A New light hydrocarbon generation model

Few question the biological pedigree of higher hydrocarbons (biomarkers) in petroleum [3, 4]. Their carbon structures are precise fits to the carbon skeletons of bio-precursors. Light hydrocarbons between C1 and about C10 are different. They do not resemble biological parents [25], and they display distinctive patterns in composition not seen in the higher hydrocarbons. Methane through butane (C1-C4) are constrained to thermodynamic equilibrium [7, 10] and the hexanes and heptanes (C6-C7) display invariant compositions [26, 27], for example. A proportionality between (n-C6*i-C7) and (n-C7*i-C6), perhaps the highest correlation yet reported for hydrocarbons in crude oils (R2 = 0.99), illustrates their extraordinary compositional order [27].

It is in this context that methane-solid hydrocarbon equilibrium must be weighed. With methane in equilibrium with ethane and propane [7], it should come as no surprise that it would be in equilibrium with the solid hydrocarbons generating it. In fact, there is a strong possibility that all light hydrocarbons form along the same path, with light hydrocarbons and solid hydrocarbons at equilibrium. It is illustrated in a new model for light hydrocarbon generation presented here. The model unifies otherwise disparate elements of light hydrocarbon chemistry: natural gas at thermodynamic equilibrium, invariance in the higher light hydrocarbons, ingestion, and gas generation during production. It does this by bringing light hydrocarbons and solid hydrocarbons into equilibrium through catalytic intermediates [Cx] and [···]:

Cn + [ ] < = > Cx < = > Cx + [ ]

[Cx] represents the intermediate shaping product structures and compositions. Metathesis of [Cx] brings C1 to C4 to equilibrium [10] and molecular rearrangements of [Cx] bring C5 to C7 isomers to invariance [26, 27]. Ingestion, which links methane to [Cx] and thus to solid hydrocarbons, is the critical step in the proposed methane-solid hydrocarbon equilibrium. Mowry and Fayetteville Shales ingest light hydrocarbons in substantial amounts and with high selectivity. Mowry Shale in argon with open receptors ingests C2-C4 hydrocarbons vigorously, and the same shale in methane with pacified receptors totally rejects the same hydrocarbons under the same conditions (Figures 1 and 2). Generation and storage is expressed in the model through: [Cn] + [···] < = > [Cx]. It accounts for the linearity in double exponential decay decline curves (Figures 7 & 8). Thus, each component of the general equilibrium model has empirical support.

Most catalytic reactions are reversible, and approach thermodynamic equilibrium over time (residence time). It therefore follows that C1-C4 hydrocarbons residing in closed source rocks over geologic time will be at molecular and isotopic thermodynamic equilibrium, and steady state with respect to compositional change. However, once the rock is opened and old hydrocarbons escape, new hydrocarbons will be generated and their residence times can be on the order of hours. C1-C4 could then be removed from equilibrium and the C6 and C7 hydrocarbons, typically constrained to constant compositions including metastable equilibria [2628], displaced from these compositions as well. In other words, molecular and isotopic compositions of hydrocarbons generated at steady state can be distinct from hydrocarbons generated at pre-steady state. Our research has focused on hydrocarbons generated at steady state over geologic time. There is less know about pre-steady state because there have been fewer opportunities to find and analyze pre-steady state products. However, we have encountered oils from conventional reservoirs with bizarre C6 and C7 distributions consistent with pre-steady state kinetics at the onset of oil generation [29]. Unconventional production from source rocks offers opportunities for finding similar pre-steady state hydrocarbons. Compositions could transition between states in early production, capturing molecular and isotopic biases only rarely seen in conventional reservoirs. The laboratory for finding evidence of that transition is in the field, from wells producing oil and gas from source rocks, where the transition might be captured and analyzed.

Conclusions

Production rates declining by DED1 is not in itself significant. But, the coincidence of theory and data in Figure 6 is. It implicates gas-solid hydrocarbon equilibria in the production of gas from Fayetteville Shale.

Ingestion gives the hypothesis additional and independent support. First, it is unprecedented. Light hydrocarbons react with solid hydrocarbons in source rocks under conditions where they do not react with ordinary solid hydrocarbons. Methane, perhaps the least reactive hydrocarbon known, disappears in Fayetteville Shale at 100°C leaving no trace of CH4. It does not emerge in RockEval analysis in the S1 peak at 300°C, but later as some higher hydrocarbon in the S2 peak at 443°C (Table 2). Methane thus becomes part of the solid hydrocarbon ingesting it: C1 + [···] = > [C1]. The reactions between C1-C4 hydrocarbons and Mowry Shale at 70°C are equally striking. Adamantanes, like methane in thermal stability and reactivity, become highly reactive on carbon surfaces [30]. Solid hydrocarbons in source rocks are not ordinary hydrocarbons. They possess receptors, perhaps catalytic, that react with light hydrocarbons at low temperatures generating solid hydrocarbons of greater mass. The Miller Heirs field data and the ingestion results reported here make a substantial case for methane and solid hydrocarbons in equilibrium in source rocks. The fact that a general equilibrium also accounts for other properties of light hydrocarbons - natural gas at thermodynamic equilibrium and higher hydrocarbons in constant compositions – makes that case even stronger. The capacity of solid hydrocarbons to ingest gas reported here raises the possibility of substantially more gas in shales than analytically apparent, and far more gas in shale deposits than currently recognized.

Only the question of reversibility and perhaps generation remain. Catalysis would seem a given since methane could not react with solid hydrocarbon at 70°C (Table 1) without catalytic assistance. It is, in my view, extremely unlikely that methane could be found in equilibrium with ethane and propane in reservoir rocks [7] and not have been in equilibrium with solid hydrocarbons in source rocks. The Miller Heirs production curve in Figure 7 suggests that it is and that gas generated from solid hydrocarbons sustains production over time (Figure 8). That possibility becomes near certainty should hydrocarbons produced from source rocks transition from steady state to pre-steady state and those produced from conventional reservoirs do not.

Acknowledgements

Author thanks Petroleum Habitats and Worldwide Geochemistry for their generous support.

References

  1. Evans CR, Rogers MA, Bailey NJL: Evolution and alteration of petroleum in western Canada. Chem Geol. 1971, 8: 147-170. 10.1016/0009-2541(71)90002-7.

    Article  Google Scholar 

  2. McNab JG, Smith PV, Betts RL: The evolution of petroleum. Petrol Eng Chem. 1952, 44: 2556-2563.

    Google Scholar 

  3. Tissot BP, Welte DH: Petroleum Formation and Occurrence. 1984, New York: Springer

    Book  Google Scholar 

  4. Hunt JM: Petroleum Geochemistry and Geology. 1995, New York: Freeman

    Google Scholar 

  5. Mango FD: The stability of hydrocarbons under time-temperature conditions of petroleum genesis. Nature. 1991, 352: 146-148. 10.1038/352146a0.

    Article  Google Scholar 

  6. Mango FD: Methane concentrations in natural gas: the genetic implications. Org Geochem. 2001, 32: 1283-1287. 10.1016/S0146-6380(01)00099-7.

    Article  Google Scholar 

  7. Mango FD, Jarvie DM, Herriman E: Natural gas at thermodynamic equilibrium: implications for the origin of natural gas. Geochem Trans. 2009, 10: 6-10.1186/1467-4866-10-6.

    Article  Google Scholar 

  8. Mango FD, Jarvie DM: Low-temperature gas generation from marine shales. Geochem Trans. 2009, 10: 3-10.1186/1467-4866-10-3.

    Article  Google Scholar 

  9. Mango FD, Jarvie DM: The low-temperature catalytic path to natural gas: wet gas to dry gas over experimental time. Geochem Trans. 2009, 10: 10-10.1186/1467-4866-10-10.

    Article  Google Scholar 

  10. Mango FD, Jarvie DM: Metathesis in the generation of low-temperature gas in marine shales. Geochem Trans. 2010, 11: 1-10.1186/1467-4866-11-1.

    Article  Google Scholar 

  11. Mango FD, Jarvie DM, Herriman E: Natural catalytic activity in a marine shale for generating natural gas. Proc Roy Soc A. 2010, 466: 3527-3537. 10.1098/rspa.2010.0032.

    Article  Google Scholar 

  12. Lewan MD, Korarba MJ, Wueckaw D, Iestrzynski A: Evaluating transition-metal catalysis in gas generation from Permian Kupferscheifer by hydrous pyrolysis. Geochica Cosmochim Acta. 2008, 72: 4069-4093. 10.1016/j.gca.2008.06.003.

    Article  Google Scholar 

  13. Parnell J: Native platinum in pyrobitumen from Fonda, New York. Am Mineral. 1988, 73: 1170-1171.

    Google Scholar 

  14. Yen TF: Chemical aspects of metals in native petroleum. The Role of Trace Metals in Petroleum. 1975, Ann Arbor MI: Ann Arbor Science Publishers Inc, 1-30.

    Google Scholar 

  15. Pan C, Jiang L, Liu J, Zhang S, Zhu G: The effects of pyrobitumen on oil cracking in confined pyrolysis experiments. Org Geochem. 2012, 45: 29-47.

    Article  Google Scholar 

  16. Alexander R, Berwick LJ, Pierce K: Single carbon surface reactions of 1-octadecene and 2,3,6-trimethylphenol on activated carbon: implications for methane formation in sediments. Org Geochem. 2011, 42: 540-547. 10.1016/j.orggeochem.2011.03.005.

    Article  Google Scholar 

  17. Somorjai GA: Introduction to Surface Chemistry and Catalysis. 1994, New York: John Wiley & Sons, 489-

    Google Scholar 

  18. Plambeck JA: Introductory University Chemistry I, Solutions and Solubility of Gases. 2013, http://www.chem.ualberta.ca/courses/plambeck/p101/pou82.htm,

    Google Scholar 

  19. Pines H: The Chemistry of Catalytic Hydrocarbon Conversions. 1981, New York: Academic, 6-

    Google Scholar 

  20. Arps JJ: Analysis of decline curves. Trans AIME. 1945, 160: 228-247.

    Article  Google Scholar 

  21. Raghavan R: Behavior of wells completed in multiple producing zones: SPEFE. 1989, 219-230.

    Google Scholar 

  22. Fetkovich MJ, Works AM, Thrasher TS, Bradley MD: Depletion performance of layered reservoirs without crossflow: SPEFE. 1990

    Google Scholar 

  23. Ilk D, Rushing JA, Perego AD, Blasingame TA: paper SPE 116731 presented 2008 Annual Technical Conference & Exhibition, Denver, CO. Exponential vs. hyperbolic decline in tight gas sands – Understanding the origin and implications for reserve estimates using Arps’ decline curves. 2008

    Chapter  Google Scholar 

  24. Valko PP, Lee J: paper SPE 134231 presented 2010 Annual Technical Conference & Exhibition, Florence, Italy. A better way to forecast production from unconventional gas wells. 2010

    Chapter  Google Scholar 

  25. Mango FD: The light hydrocarbons in petroleum: a critical review. Org Geochem. 1997, 26: 417-440. 10.1016/S0146-6380(97)00031-4.

    Article  Google Scholar 

  26. Mango FD: An invariance in the isoheptanes of petroleum. Science. 1987, 237: 514-517. 10.1126/science.237.4814.514.

    Article  Google Scholar 

  27. Mango FD: The origin of light hydrocarbons. Geochim Cosmochim Acta. 2000, 64: 1265-1277. 10.1016/S0016-7037(99)00389-0.

    Article  Google Scholar 

  28. Helgeson HC, Knox AM, Owens CE, Shock EI: Petroleum oil field waters, and authigenic mineral assemblages: are they in metastable equilibrium in hydrocarbon reservoirs. Geochim Cosmochim Acta. 1993, 57: 3295-3339. 10.1016/0016-7037(93)90541-4.

    Article  Google Scholar 

  29. Mango FD: Pre-steady state kinetics at the onset of oil generation. Org Geochem. 1990, 16: 41-48. 10.1016/0146-6380(90)90024-T.

    Article  Google Scholar 

  30. Berwick L, Alexander R, Pierce K: Formation and reactions of alkyl adamantanes in sediments: carbon surface reactions. Org Geochem. 2011, 42: 752-761. 10.1016/j.orggeochem.2011.05.008.

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Frank D Mango.

Additional information

Competing interests

The author declares that there are no competing interests.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mango, F.D. Methane and carbon at equilibrium in source rocks. Geochem Trans 14, 5 (2013). https://doi.org/10.1186/1467-4866-14-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1467-4866-14-5

Keywords