# The Dynamic Evolution Model of the Chemical and Carbon Isotopic Composition of C1–3 during the Hydrocarbon Generation Process

## TL;DR

This paper introduces a dynamic model to track how the chemical and isotopic composition of C1–3 gases evolve during hydrocarbon generation.

## Contribution

A novel dynamic model is proposed to quantify the isotopic and chemical evolution of C1–3 gases during hydrocarbon generation.

## Key findings

- The model successfully captures the isotopic evolution of C1–3 gases across different maturity levels.
- The model can be applied to evaluate natural gas resources and understand hydrocarbon origins.
- It provides insights into isotopic reversal mechanisms in hydrocarbon evolution.

## Abstract

A new approach is presented in this paper for the dynamic modeling of the chemical and isotopic evolution of C1–3 during the hydrocarbon generation process. Based on systematic data obtained from published papers for the pyrolysis of various hydrocarbon sources (type I kerogen/source rock, type II kerogen/source rock, type III kerogen/source rock, crude oil, and asphalt, etc.), the empirical evolution framework of the chemical and isotopic composition of C1–3 during the hydrocarbon generation process was built. Although the empirical framework was built only by fitting a large amount of pyrolysis data, the chemical and isotopic composition of C1–3 derived from the pyrolysis experiments all follow evolution laws, convincing us that it is applicable to the thermal evolution process of various hydrocarbon sources. Based on the simplified formula of the isotopic composition of mixed natural gas at different maturities (δ13Cmixed), δ13Cmixed = X×niA×δ13CiA+Y×niB×δ13CiBX×niA+Y×niB, it can be derived that the cumulative isotopic composition of alkane generated in a certain maturity interval can be expressed by the integral of the product of the instantaneous isotopic composition and instantaneous yield at a certain maturity point, and then divided by the cumulative yield of alkane generated in the corresponding maturity interval. Thus, the cumulative isotopic composition (A(X)), cumulative yield (B(X)), instantaneous isotope (C(X)), and instantaneous yield (D(x)) in the dynamic model, comply with the following formula during the maturity interval of (X0~X). A(X) = ∫X0XCX×DXdxB(X), where A(X) and B(X) can be obtained by the fitting of pyrolysis data, and D(x) can also be obtained from the derivation of B(X). The dynamic model was applied on the pyrolysis data of Pingliang Shale to illustrate the quantitative evolution of the cumulative yield, instantaneous yield, cumulative isotope, and instantaneous isotope of C1–3 with increasing maturity. The dynamic model can quantify the yield of methane, ethane, and propane, as well as δ13C1, δ13C2, and δ13C3, respectively, during the hydrocarbon generation process. This model is of great significance for evaluating the natural gas resources of hydrocarbon source rock of different maturities and for identifying the origin and evolutionary process of hydrocarbons by chemical and isotopic data. Moreover, this model provides an approach to study the dynamic evolution of the isotope series of C1–3 (including reversed isotopic series), which is promising for revealing the mechanism responsible for isotopic reversal when combined with post-generation studies.