Relative sorption coefficient: Key to tracing petroleum migration and other subsurface fluids
L. Zhang1,2,
Y. Wang3,4,
M.-W. Li5,
Q.-Z. Yin ORCID: orcid.org/0000-0002-4445-50966 &
W. Zhang7
Scientific Reports volume 9, Article number: 16845 (2019) Cite this article
The accumulation and spatial distribution of economically important petroleum in sedimentary basins are primarily controlled by its migration from source rocks through permeable carrier beds to reservoirs. Tracing petroleum migration entails the use of molecular indices established according to sorption capacities of polar molecules in migrating petroleum. However, little is known about molecular sorption capacities in natural migration systems, rendering these indices unreliable. Here, we present a new approach based on a novel concept of relative sorption coefficient for quantitatively assessing sorption capacities of polar molecules during natural petroleum migration. Using this approach, we discovered previously unrecognized "stripping" and "impeding" effects that significantly reduce the sorption capacities of polar compounds. These discoveries provide new insights into the behaviors of polar compounds and can easily explain why traditional molecular indices yield incorrect information about petroleum migration. In light of these new findings, we established new molecular indices for tracing petroleum migration. We demonstrate via case studies that the newly established indices, unlike traditional molecular indices, are reliable and effective in tracing petroleum migration. Our approach can be applied to diverse basins around the world to reveal distribution patterns of petroleum, which would decrease environmental risks of exploration by reducing unsuccessful wells.
Petroleum is produced in source rocks through thermal alteration of organic matter buried in sedimentary basins. Its accumulation and distribution are mainly controlled by secondary petroleum migration (SPM) through permeable carrier beds following petroleum expulsion (primary migration) out of the source rocks. Tracing petroleum migration can reveal distribution patterns of petroleum reservoirs and thus increase exploration success rate. Meanwhile, environmental risks can be decreased by reducing unsuccessful wells. Biomarker hydrocarbon geochemistry can be applied to trace SPM1,2,3,4,5,6. The underlying principle of this approach is that polar molecules are preferentially removed from petroleum during secondary migration due to sorption onto immobile mineral surfaces, which causes their concentrations to decrease with increasing migration distance1,3,5,7. Thus, tracing SPM entails the use of molecular indices established on the basis of sorption capacities of polar compounds in migrating petroleum (Supplementary Text S-1.1). However, there have been no reports of quantitative research on the sorption capacities of polar molecules in naturally migrating petroleum, although it has been suggested that the sorption capacities of alkylcarbazoles are determined by the (partial) shielding effect related to alkylation at positions 1 and/or 8 (refs 3, 6) (Fig. 1). Sorption of trace polar compounds during lateral petroleum migration typically reaches equilibrium2,4,5,8 and thus their sorption capacities essentially represent equilibrium sorption capacities determined by both sorption and desorption. However, the shielding and partial shielding effects consider only the sorption, not the desorption, of alkylcarbazoles during petroleum migration. Therefore, the previously-proposed theory on molecular sorption capacities, which was based solely on the shielding and partial shielding effects, needs to be re-evaluated to avoid erroneous results when applying molecular indices to trace SPM, as demonstrated by the case studies herein.
Structure of 1,7-dimethylcarbazole. The bond length data and angles are from Kurahashi et al.32 and Song33; the Van der Waals radii of a methyl are from Chen et al.34.
In this paper, we propose a new parameter, the Relative Sorption Coefficient (RSC), which quantitatively describes equilibrium sorption capacities of polar compounds in migrating petroleum. We establish its computation method (see Methods), and then apply and test the validity of the new approach, by using natural petroleum samples collected from the Xifeng Oilfield in the Ordos Basin and the Rimbey-Meadowbrook Reef Trend in the Western Canada Sedimentary Basin (WCSB). Because the petroleum in the Xifeng Oilfield contains very little benzocarbazoles5, we only analyzed alkylcarbazoles in the samples from this oilfield, although both alkylcarbazoles and benzocarbazoles can serve as important tracers1,5,9,10. Using the RCS approach, we determined the sorption capacities of an important group of polar molecules, the alkylcarbazoles, and discovered two previously unrecognized effects (i.e., stripping and impeding effects) that strongly influence the equilibrium sorption capacities of polar compounds. In light of these new findings, we reclassify alkylcarbazoles into six subgroups, according to their equilibrium sorption capacities, and propose new ratios that consist of numerators with stronger sorption capacities than denominators. All the new ratios show a significantly decreasing trend with increasing migration distance, demonstrating that the new ratios are reliable and effective indices for migration and that the RSC provides the key to tracing SPM. Our new method is further validated through the analyses of benzocarbazoles in petroleum samples from the Rimbey-Meadowbrook reef trend in the WCSB. Because the RSC is established on the basis of a physicochemical sorption model, it should be widely applicable to assessing equilibrium sorption capacities of solutes in a diverse range of geofluids.
Modification of secondary migration fractionation index
We have previously investigated SPM in the Xifeng Oilfield in the southwest part of the Ordos Basin in China using the secondary migration fractionation index (SMFIs). The geological setting, samples and geochemical data are documented in Zhang et al.5. Besides sorption, this earlier study also examined other factors11,12,13,14,15,16,17,18,19,20 that may influence the concentrations of polar compounds in migrating petroleum. It illustrated that the effect of thermal maturation of source rocks on polar molecule can be eliminated in the derivation of SMFIs (see Zhang et al.5 for details). Other influences (i.e., organic facies of source rocks, biodegradation of petroleum and dissolution in water) can be neglected for alkylcarbazoles in the petroleum of this field5. In this study, we re-calculated the SMFIs, using the migration-sorption fractionation equation with a quadratic polynomial (Supplementary Equation (S3)), to improve the accuracy of the results (Supplementary Text, S-1.2-3). The results clearly show an exponential decrease of SMFIs with increasing relative migration distance (Supplementary Fig. S1A–C,G–I), suggesting that the Xifeng Oilfield was likely formed by SPM in the SW direction along the sand body from the source kitchen located in the NE of the reservoir (refer to Zhang et al.5).
However, the information about petroleum migration derived from the SMFIs needs to be verified by using the ratios of the SMFIs. Because the SMFI is affected by the relative rates of concentration variations of polar molecules at the starting point or at a reference point of SPM (Supplementary Text, S-1.2-3), it cannot be used to construct reliable ratios. To overcome this problem, we revised the SMFI (Supplementary Equation (S19)) and its related ratios (Supplementary Text, S-1-3). The amended indices are denoted by the subscript λ (e.g. SMFIλ). The values of SMFIλ (Supplementary Table S4) display similar distribution trends as the values of SMFI Supplementary Fig. S1), even though the powers of their regression equations are different, as shown by Supplementary Equation (S16). The ratios of the SMFIλs of alkylcarbazoles with stronger sorption capacities to those with weaker sorption capacities should decrease with increasing migration distance if the underlying assumptions about source facies, biodegradation and thermal maturity effects are valid. The ratios were initially established based on the previously proposed theory about sorption capacities that considered only the shielding and partial shielding effects3,6,21,22. Based on this theory, alkylcarbazole isomers can be divided into three groups3,6: N-H shielded (Group I), N-H partially shielded (Group II) and N-H exposed (Group III). Their sorption capacities are expected to decrease in the order of Group III> Group II> Group I, and the ratios of the SMFIλs of dimethylcarbazoles (DMCAs) in Group III to those in Group II, III to I and II to I would be predicted to decrease with increasing migration distance. However, many of these ratios (Fig. 2) do not display a decreasing trend with increasing migration distance, but instead they exhibit a clear increasing trend (Supplementary Text, S-1.5), which is completely opposite to the decreasing trend of SMFIs and SMFIλs (Supplementary Fig. S1), and inconsistent with the geological conditions (Fig. 1 in Zhang et al.5). From this, we can see that if the ratios of carbazoles, constructed on the basis of the existing theory on sorption capacities, are used to trace SPM, it would yield erroneous or misleading information about petroleum migration. Similarly, the ratios based on the current sorption capacity theory cannot be used to verify the information about SPM that is inferred from the SMFIs and SMFIλs. Therefore, their use should be discontinued.
Correlation diagrams showing the relationships between relative migration distances and the SMFIλ ratios of different groups of dimethylcarbazoles in the Xifeng Oilfield. SMFIλ: the amended Secondary Migration Fractionation Index; DMCA: dimethylcarbazole. The number of data points in each panel is nineteen; the grey error bars indicate one standard deviation (1σ) from the logarithmic values of the SMFIλ ratios. 1,4/1,8-DMCA SMFIλ: ratio of the SMFIλ of 1,4-DMCA to the SMFIλ of 1,8-DMCA (A); 3,4/1,8-DMCA SMFIλ: ratio of the SMFIλ of 3,4-DMCA to the SMFIλ of 1,8-DMCA (B); 2,4/1,5-DMCA SMFIλ: ratio of the SMFIλ of 2,4-DMCA to the SMFIλ of 1,5-DMCA (C); 2,4/1,6-DMCA SMFIλ: ratio of the SMFIλ of 2,4-DMCA to the SMFIλ of 1,6-DMCA (D); 2,5/1,7-DMCA SMFIλ: ratio of the SMFIλ of 2,5-DMCA to the SMFIλ of 1,7-DMCA (E); 2,7/1,2-DMCA SMFIλ: ratio of the SMFIλ of 2,7-DMCA to the SMFIλ of 1,2-DMCA (F); 3,4/1,7-DMCA SMFIλ: ratio of the SMFIλ of 3,4-DMCA to the SMFIλ of 1,7-DMCA (G); 2,5/1,2-DMCA SMFIλ: ratio of the SMFIλ of 2,5-DMCA to the SMFIλ of 1,2-DMCA (H); 2,4/1,2-DMCA SMFIλ: ratio of the SMFIλ of 2,4-DMCA to the SMFIλ of 1,2-DMCA (I). The SMFIλ ratio of 1 at the reference point (x = 0 km) is the model value, and was excluded from regression analyses shown in this figure.
Relative sorption coefficient (RSC)
We re-examined the sorption capacities of polar compounds in petroleum samples from the Xifeng Oilfield using our new approach described in the Methods and Supplementary Text S-1. We calculated the relative sorption coefficients - the Kr values of alkylcarbazoles in these petroleum samples (Supplementary Table S2). The Kr values vary widely. Some N-H partially shielded DMCAs (Group II) have higher Kr values than some of the N-H exposed DMCAs (Group III), which cannot be explained by the existing theory on sorption capacities that considered only the shielding and partial shielding effects. This suggests that there are other factors controlling equilibrium sorption capacities of alkylcarbazoles.
Through comparison of desorption of the adsorbed polar compounds under both static and dynamic conditions (Supplementary Text, S-1.6), a stripping effect was observed arising from petroleum migration that causes "tall" alkylcarbazoles to desorb more easily than the "short" ones (Supplementary Figs S4, S5). This stripping effect greatly reduces the equilibrium sorption capacities of alkylcarbazoles with the alkyl substituents at positions 4 and/or 5, as is demonstrated by the Kr values.
In Group II alkylcarbazoles, the molecular height of 1,4-DMCA is greater than for 1,5- and 1,3-DMCA (Supplementary Text, S-1.6), and the latter two are taller than the other DMCAs in this group. Consequently, 1,4-DMCA has a lower Kr value than 1,5- and 1,3-DMCA, which have lower Kr values than the rest of the DMCAs in this group (Supplementary Fig. S3). In light of the stripping effect, the Group II alkylcarbazoles are further divided into three subgroups (Fig. 3) with decreasing stripping effect and increasing sorption capacity in the following order: N-H partially shielded alkylcarbazole with the alkyl at position 4 (Subgroup II-1), N-H partially shielded alkylcarbazoles with the alkyl at positions 3 or 5 (Subgroup II-2), and N-H partially shielded without the alkyl at position 3, 4 or 5 (Subgroup II-3) (Figs 3 and S3).
The sequence of individual dimethylcarbazoles in different groups sorted by their relative sorption coefficients. DMCA = dimethylcarbazole; Kr = relative sorption coefficient (%). In order to plot the datum point of 3,4-DMCA (Kr = 0%) using the logarithmic coordinate, 0.01% were added to all the Kr values. The black bars indicate the errors of the relative sorption coefficient calculated with Supplementary Equation (S29) in Supplementary Text S-1.4.
In Group III alkylcarbazoles, 3,4-DMCA has two methyls sticking out, and is subject to a stronger stripping effect (two-methyl stripping; see Supplementary Text, S-1.6 for details) and thus has a lower Kr value than 2,4- or 2,5-DMCA (Supplementary Fig. S5A–C). The Kr value of 3,4-DMCA is even lower than that of 1,8-DMCA, which experiences the shielding effect (Fig. 3).
In addition, we discovered an impeding effect related to the alkyls at positions 2 and 7 (Supplementary Text, S-1.7). The impeding effect causes the equilibrium sorption capacity of 2,7-DMCA in the Group III alkylcarbazoles to become lower than those of 2,3- and 2,6-DMCA (Supplementary Figs. S5D–F). Due to the stripping and impeding effects, alkylcarbazoles in Group III show large variations in Kr (Supplementary Fig. S4) and are also divided into three subgroups (Fig. 3): N-H exposed alkylcarbazole with the alkyls at positions 3 and 4 (Subgroup III-1), N-H exposed alkylcarbazoles with one alkyl at position 2 and the other alkyl at positions 4, 5 or 7 (Subgroup III-2), and N-H exposed alkylcarbazoles without the alkyls at positions 4, 5 or 7 (Subgroup III-3).
The stripping and impeding effects, which control the sorption capacities of polar molecules in migrating petroleum, are also related to the molecular structures of the organic compounds, just like the shielding effect. Based on the three effects noted above and the Kr values (Supplementary Table S2), these subgroups and the Group I alkylcarbazole can be arranged in the following sequence with decreasing equilibrium sorption capacity (Fig. 3): Subgroup III-3 (Kr = 78.5–100%)> III-2 (Kr = 0.40–5.6%), II-2 and II-3 [II-3 (Kr = 0.97–12.0%)> II-2 (Kr = 0.38–0.46%)]> II-1 (Kr = 0.11%) and Group I (Kr = 0.10%)> III-1 (Kr = 0.0%).
From the above analyses of molecular structures and their relationships with Kr, we established the following sequence of various effects on reducing the equilibrium sorption capacities: two-methyl stripping (represented by Subgroup III-1)> shielding (Group I)> partial shielding plus one-methyl stripping (II-1 and -2)> partial shielding (II-3), one-methyl stripping and impeding (III-2)> partial impeding. The interplay of these three effects results in complex variations in equilibrium sorption capacities of the DMCAs within and among subgroups. The seemingly unreasonable relationships of SMFI ratios with relative migration distance (Fig. 2) can all be explained by these effects and their combination(Supplementary Text S-1.8).
The relative sorption coefficient is derived from the linear isotherm model that is the simplification of the Langmuir isotherm model of equilibrium sorption at low concentrations of adsorbents such as carbazoles (refer to the Methods Section, Supplementary Information and Zhang et al.5). Recent studies on sorption of asphaltenes onto minerals show that the Langmuir isotherm model can be used to describe the equilibrium adsorption of asphaltenes when interactions between the solute and the solvent as well as interactions that can occur at a non-ideal lattice of a mineral are negligible and that the sorption of asphaltenes is highly dependent on the heteroatoms (i.e. N, O, S) in their molecular structure23,24,25. These results confirm the validity of using the linear isotherm model to investigate the equilibrium sorption of polar heteroatom compounds such as carbazoles onto solid surfaces.
New ratios and their application
Given the sorption capacity sequence of the subgroups and Group I, eighteen SMFIλ ratios are established as indices for petroleum migration: alkylcarbazoles in Subgroup III-3 to those in III-2, III-3 to III-1, III-3 to II-3, III-3 to II-2, III-3 to II-1, III-3 to Group I, III-2 to III-1, III-2 to II-1, III-2 to Group I, II-3 to II-2, II-3 to II-1, II-3 to Group I, II-3 to III-1, II-2 to II-1, II-2 to Group I, II-2 to III-1, II-1 to III-1, and Group I to III-1. Since the equilibrium sorption capacities of the numerators are significantly higher than those of the denominators, these ratios decrease with increasing migration distance and thus can serve as odometers for SPM (Supplementary Equation (S23)). Similarly, the corresponding ratios of the geometric means of SMFIλs decrease with increasing migration distance and can also be used as indices for petroleum migration (Supplementary Equations (S21 and S22)). It is worth noting that the ratios of alkylcarbazoles within each group (except Group I with only one compound), which were not considered previously, can also be useful in the establishment of new indices (Supplementary Text S-1.9).
The new SMFIλ ratios for the Xifeng Oilfield fit the known data well, clearly showing exponential decreases with increasing migration distance with high correlation coefficients (Fig. 4). These are consistent with the migration fractionations inferred from SMFIλs and SMFIs, and geological conditions26,27,28,29 (Fig. 1 in Zhang et al.5). Thus, the new ratios confirm the validity of the influence elimination and migration information revealed by the SMFIs5 and SMFIλs, and demonstrate that the petroleum migrated along the sand body from the source kitchen into the Xifeng Oilfield in a SW direction (refer to Zhang et al.5).
Correlation diagrams showing the relationships between relative migration distances and SMFIλ ratios of different subgroups of dimethylcarbazoles in the Xifeng Oilfield. SMFIλ: the modified Secondary Migration Fractionation Index; DMCA: dimethylcarbazole. The number of the data points in each panel is nineteen; the grey error bars indicate 1σ from the logarithmic values of the SMFIλ ratios. 2,3/1,3-DMCA SMFIλ: ratio of the SMFIλ of 2,3-DMCA to the SMFIλ of 1,3-DMCA (A); 2,6/1,6-DMCA SMFIλ: ratio of the SMFIλ of 2,6-DMCA to the SMFIλ of 1,6-DMCA (B); 2,5/1,5-DMCA SMFIλ: ratio of the SMFIλ of 2,5-DMCA to the SMFIλ of 1,5-DMCA (C); 2,6/2,4-DMCA SMFIλ: ratio of the SMFIλ of 2,6-DMCA to the SMFIλ of 2,4-DMCA (D); 2,7/1,4-DMCA SMFIλ: ratio of the SMFIλ of 2,7-DMCA to the SMFIλ of 1,4-DMCA (E); 1,2/1,8-DMCA SMFIλ: ratio of the SMFIλ of 1,2-DMCA to the SMFIλ of 1,8-DMCA (F); 1,2/1,3-DMCA SMFIλ: ratio of the SMFIλ of 1,2-DMCA to the SMFIλ of 1,3-DMCA (G); 1,2/1,4-DMCA SMFIλ: ratio of the SMFIλ of 1,2-DMCA to the SMFIλ of 1,4-DMCA (H). GMλ(III-3), GMλ(II-3) and GMλ(II-2): the geometric means of the SMFIλs of subgroups III-3, II-3, II-2 and II-1, respectively. GMλ(I) and GMλ(II-1) are the SMFIλs of 1,8-DMCA and 1,4-DMCA, respectively (refer to Fig. 3). The SMFIλ ratio of 1 at the reference point (x = 0 km) is the model value, which was excluded in the regression analyses in this figure.
Molecular indices that are correlated with migration directions, pathways and distances have been sought based on sorption capacities of polar organic compounds in migrating petroleum for decades1,5, but with limited success5 because the sorption capacities of these polar compounds have been unclear. As a result, reliable indices have not been established and secondary petroleum migration still remains the least understood of the processes involved in petroleum accumulation5. The results of the application of the new ratios of alkylcarbazoles in the petroleum of the Xifeng Oilfield, however, demonstrate that the relative sorption coefficients (Kr) can be used to assess the sorption capacities and that the new ratios established on the basis of the new understanding of the sorption capacities can serve as effective indices for petroleum migration. These new indices provide a powerful tool for revealing migration directions, pathways and distances that control petroleum distribution patterns in reservoirs in basins, which would greatly facilitate future petroleum exploration and increase the success rate of wells.
Furthermore, the new observation of the stripping effect on equilibrium sorption capacities is supported by the analyses of benzocarbazoles in the petroleum samples from the Rimbey-Meadowbrook reef trend of central Alberta, Canada. In these petroleum samples, the "taller" benzo[c]carbazole has a lower Kr value than the "shorter" benzo[a]carbazole, consistent with predictions from the stripping effect (Supplementary Table S5).
The new concept of RSC overcomes the dependency of the sorption coefficient assessment on the migration velocity and rock characteristics of the carrier beds. As demonstrated in our case studies, the RSC provides a powerful tool with a sound scientific basis to quantitatively evaluate equilibrium sorption capacities of polar compounds during petroleum migration, and can help uncover factors controlling equilibrium sorption capacities. Without this tool, it would be impossible to quantitatively assess equilibrium sorption capacities of polar compounds in migrating petroleum and to establish reliable molecular indices for tracing petroleum migration. The lack of a quantitative assessment tool is also the primary reason why many of the previously proposed molecular indices failed to provide reliable information about SPM. Application of this approach to quantitative assessment of equilibrium sorption capacities of alkylcarbazoles has resulted in the discovery of the previously unrecognized stripping and impeding effects that significantly reduce the equilibrium sorption capacities of polar compounds. These findings have led to the reclassification of the polar compounds according to their sorption capacities. Based on the reclassification of the polar compounds, we established eighteen new ratios. As demonstrated in our case studies, these new indices provide reliable information about petroleum migration (i.e. migration directions, routes and distances). Therefore, this approach is the key to tracing secondary petroleum migration and can be applied to petroliferous basins around the world, to reveal distribution patterns of petroleum reservoirs, which would help to find more petroleum and decrease environmental risks of exploration by reducing unsuccessful wells.
Moreover, the concept of RSC and its evaluation method developed in this study should be applicable in hydrological and environmental studies (as well as other disciplines) to trace the movement of pollutants and water (and other geofluids) (Supplementary Text S-1.10).
The equilibrium sorption of a polar molecule or an adsorbable element in a natural migration system of petroleum or other geofluids can be described by the linear isotherm model if its concentration is sufficiently low5,8,30,31. In this physicochemical model, the sorption coefficient Kd (cm3/g) represents the sorption amount of a polar compound or an adsorbable element at a given concentration and saturation of petroleum or geofluid8,30,31 (see Supplementary Equation (S5) in Zhang et al.)5. This amount may describe the equilibrium sorption capacity of the compound or element, according to Delle Site8. However, the Kd values determined in laboratories are not necessarily applicable to natural migration systems, due to differences in size, time and distance between laboratory experiments and natural migration systems. Moreover, lab experimental studies for the determination of sorption coefficients are expensive and time consuming, and the results may not be accurate, especially when concentrations are low8. Above all, Kd is also controlled by many factors such as the porosity, density of carrier beds and the average velocity of migration (Supplementary Text S-1.2). Therefore, the sorption coefficient Kd cannot be used directly to describe the equilibrium sorption capacities of polar organic compounds or trace elements during lateral migration.
To evaluate equilibrium sorption capacities of polar compounds (or adsorbable elements) in natural migration systems, we introduce a new concept of relative sorption coefficient (RSC):
$${K}_{r}=\frac{{K}_{d}-{K}_{dmin}}{{K}_{dmax}-{K}_{dmin}}\times 100( \% )$$
where Kr is the RSC; Kd is the sorption coefficient (cm3/g); Kdmax is the maximum value in a series of Kd values of polar compounds in petroleum (or adsorbable elements in other geofluids); and Kdmin is the minimum value. The range of Kr values is 0–100%. Kr can be used quantitatively to evaluate equilibrium sorption capacities. High Kr values indicate strong equilibrium sorption capacities.
For the linear isotherm model of the equilibrium sorption in natural migration systems5,8,30,31, we can derive the following equation from Supplementary Equation (S8) in Zhang et al.5:
$${K}_{d}=({R}_{d}-1)\frac{n}{{n}_{s}\cdot {\rho }_{s}}\,$$
where Rd represents the retardation factor of a polar compound in migrating petroleum or an adsorbable trace element in migrating groundwater (a dimensionless constant), being related to the sorption of the compound or the element and the average migrating velocity of petroleum or groundwater (Supplementary Text S-1.2); n is the porosity of the carrier bed (%); ns = 100 − n (%); ρs is the density of the solids (g/cm3).
Migration of petroleum (or other geofluids) usually occurred in past geological times. Therefore, the current porosity and density of carrier beds do not represent the porosity and density during migration, as these lithological properties most likely have changed over time during diagenesis. Therefore, quantitative measurements of the porosity and density of carrier beds during migration can rarely be obtained. However, these parameters are the same for different compounds or for different elements in a migration system, and thus can be eliminated (Supplementary Text S-1.2) when Eq. (2) is substituted into Eq. (1):
$${K}_{r}=\frac{{R}_{d}-{R}_{dmin}}{{R}_{dmax}-{R}_{dmin}}\times 100( \% )$$
where Rdmax is the maximum value in a series of Rd values of polar compounds or elements; Rdmin is the minimum value. Rd is also controlled by the average velocity of migration and the difference in relative variation rates of concentrations with time at the starting point of a migration pathway between polar compounds. However, it is demonstrated that the RSC can also eliminate these two kinds of influences when Supplementary Equations (S11–S13) are substituted into Eq. (3) (see Supplementary Text S-1.2 for details):
$${K}_{r}=\frac{{a}_{\lambda max}-{a}_{\lambda }}{{a}_{\lambda max}-{a}_{\lambda min}}\times 100( \% )$$
where \({a}_{\lambda }\) is a constant controlling migration-sorption fractionation (km−1) and can be derived from Supplementary Equation (S9); \({a}_{\lambda max}\)is the maximum in a series of \({a}_{\lambda }\) values of polar compounds (km−1); \({a}_{\lambda min}\) is the minimum (km−1). Equation (4) provides a workable means to quantitatively evaluate sorption capacities of polar organic compounds or adsorbable trace elements.
To quantify equilibrium sorption capacities of polar organic compounds in migrating petroleum, we have established a new method for computing RSC (Kr values) of polar compounds in natural migration petroleum, on the basis of Eq. (4) (Supplementary Text S-1.2). The method for computing the relative sorption coefficient involves the following steps:
The 1st step is to conduct regression analysis using Supplementary Equation (S3) instead of Eq. (1) in Zhang et al.5, to obtain estimates of the values for the constants \({a}_{1}\), \({a}_{2}\), \({a}_{3}\) and \({a}_{4}\) that are more accurate than achievable with the previous equation in Zhang et al.5. The data preparation and the subsequent non-linear regression analyses are presented in Zhang et al.5. However, the non-linear regression analyses herein are conducted in an iterative manner (Supplementary Text S-1.2) to obtain more rational regression equations.
The 2nd step is to calculate the λ ratios (λ is the relative variation rate of the concentration at the reference point for a given polar compound) from Supplementary Equation (S8), the migration-sorption factor \({a}_{\lambda }\) (a constant controlling migration-sorption fractionation) and finally the relative sorption coefficient Kr with Supplementary Equation (S9) and Eq. (4), respectively.
The Kr values of the alkylcarbazoles in the petroleum in the Xifeng Oilfield were calculated and are listed in Supplementary Table S2(Supplementary Text S-1.2).
Larter, S. R. et al. Molecular indicators of secondary oil migration distances. Nature 383, 593–597 (1996).
Larter, S. R. et al. An experimental investigation of geochromatography during secondary migration of petroleum performed under subsurface conditions with a real rock. Geochem. Trans. 9, 1–7 (2000).
Li, M., Larter, S. R., Stoddart, D. & Bjorøy, M. Fractionation of pyrrolic nitrogen compounds in petroleum during migration: derivation of migration related geochemical parameters. Geol. Soci. London Spe. Pub. 86, 103–123 (1995).
Yang, Y. L., Aplin, A. C. & Larter, S. R. Mathematical models of the distribution of geotracers during oil migration and accumulation. Petrol. Geosci. 11(1), 67–78 (2005).
Zhang, L., Li, M., Wang, Y., Yin, Q.-Z. & Zhang, W. A novel molecular index for secondary oil migration distance. Sci. Rep. 3, 2487 (2013).
Larter, S. R. & Aplin, A. C. Reservoir geochemistry: methods, applications and opportunities. Geol. Soci. London Spe. Pub. 86, 5–32 (1995).
Li, M., Larter, S. R. & Frolov, Y. B. Adsorptive interactions between petroleum nitrogen compounds and organic/mineral phases in subsurface rocks as models for compositional fractionation of pyrrolic nitrogen compounds in petroleum during petroleum migration. J. High Res. Chrom. 17, 230–236 (1994).
Delle Site, A. Factors affecting sorption of organic compounds in natural sorbent/water systems and sorption coefficients for selected pollutants. A Review. J. Phys. Chem. Ref. Data 30(1), 187–439 (2001).
Larter, S. R. et al. Reservoir geochemistry: a link between reservoir geology and engineering? SPE Reser. Eng. 12(1), 12–17 (1997).
Li, M. Quantification of petroleum secondary migration distances: fundamentals and case histories. Petrol. Explor. Develop. 27, 11–19 (2000).
Clegg, H., Wilkes, H. & Horsfield, B. Carbazole distributions in carbonate and clastic source rocks. Geochim. Cosmochim. Acta 61, 5335–5345 (1997).
Clegg, H., Wilkes, H., Santamaria-Orozco, D. & Horsfield, B. Influence of maturity on carbazole and benzocarbazole distributions in crude oils and source rocks from the Sonda de Campeche, Gulf of Mexico. Org. Geochem. 29, 183–194 (1998).
Li, M., Yao, H., Stasiuk, L. D., Fowler, M. G. & Larter, S. R. Effect of maturity and petroleum expulsion on pyrrolic nitrogen compound yields and distributions in Duvernay Formation petroleum source rocks in central Alberta, Canada. Org. Geochem. 26, 731–744 (1997).
Galimberti, R., Ghiselli, C. & Chiaramonte, M. A. Acidic polar compounds in petroleum: a new analytical methodology and applications as molecular migration indices. Org. Geochem. 31, 1375–1386 (2000).
Bennett, B., Chen, M., Brincat, D., Gelin, F. J. P. & Larter, S. R. Fractionation of benzocarbazoles between source rocks and petroleums. Org. Geochem. 33, 545–559 (2002).
Hwang, R. J., Heidrick, T., Mertani, B. Q. & Li, M. Correlation and migration studies of North Central Sumatra oils. Org. Geochem. 33, 1361–1379 (2002).
Bakr, M. M. Y. & Wilkes, H. The influence of facies and depositional environment on the occurrence and distribution of carbazoles and benzocarbazoles in crude oils: a case study from the Gulf of Suez, Egypt. Org. Geochem. 33, 561–580 (2002).
Huang, H., Bowler, B. F. J., Zhang, Z., Oldenburg, T. B. P. & Larter, S. R. Influence of biodegradation on carbazole and benzocarbazole distributions in oil columns from the Liaohe Basin, NE China. Org. Geochem. 34, 951–969 (2003).
Lager, A., Russell, C. A., Love, G. D. & Larter, S. R. Hydropyrolysis of algae, bacteria, archaea and lake sediments: insights into the origin of nitrogen compounds in petroleum. Org. Geochem. 35, 1427–1439 (2004).
Bennett, B. & Olsen, S. D. The influence of source depositional conditions on the hydrocarbon and nitrogen compounds in petroleum from central Montana, USA. Org. Geochem. 38, 935–956 (2007).
Yamamoto, M., Taguchi, T. & Sasaki, K. Basic nitrogen compounds in bitumen and crude oils. Chem. Geol. 93, 193–206 (1991).
Yamamoto, M. Fractionation of azarenes during oil migration. Org. Geochem. 19, 389–402 (1992).
Joonaki, E., Buckman, J., Burgass, R. & Tohidi, B. Water versus asphaltenes; liquid–liquid and solid–liquid molecular interactions unravel the mechanisms behind an improved oil recovery methodology. Sci. Rep. 9, 11369 (2019).
Pradilla, D., Simon, S. & Sjöblom, J. Mixed interfaces of asphaltenes and model demulsifiers part I: Adsorption and desorption of single components. Colloids & Surfaces A: Physicochem. Eng. Aspects 466, 45–56 (2015).
Bai, Y. et al. Effects of the N, O, and S heteroatoms on the adsorption and desorption of asphaltenes on silica surface: A molecular dynamics simulation. Fuel 240, 252–261 (2019).
Yang, H. & Zhang, W. Leading effect of the seventh member high-quality source rock of the Yanchang Formation in the Ordos Basin during the enrichment of low-penetrating oil-gas accumulation: geology and geochemistry. Geochimica 34(2), 147–154 (2005).
Wu, M., Zhang, L., Luo, X., Mao, M. & Yang, Y. Analysis of hydrocarbon migration stages in the 8th member of the Yanchang Formation in the Xifeng Oilfield. Oil Gas. Geology 27(1), 33–36 (2006).
Zhang, L., Wu, M., Yang, W., Luo, X. & Chen, Z. A new model for the formation of the low-permeability oilfields in the Triassic Yanchang Formation of the Ordos Basin. Coexisting Mechanisms, Accumulation and Distribution Laws of Oil, Gas, Coal and Uranium in the Ordos Basin, (eds Liu, C. & Wu, B.) 601-631 (Science Press, Beijing, 2016).
Duan, Y., Wu, B., Zhang, H., Zheng, C. & Wang, C. Geochemistry and genesis of crude oils of the Xifeng Oilfield in the Ordos Basin. Acta Geol. Sinica. 80(2), 301–310 (2006).
Wu, Y. Contamination Transportation Dynamics in Porous Media 86–93 (Shanghai Jiaotong University Press, Shanghai, 2007).
Chen, G. Applied Physical Chemistry 86–90 (Chemical Industry Press, Beijing, 2008).
Kurahashi, M., Fukuyo, M. & Shimada, A. The crystal and molecular structure of carbazole. Bull. Chem. Soci. Japan 42, 2174–2179 (1969).
Song, H. Organic Chemistry 2–4 (China Medical Sci. & Tech. Press, Beijing, 2005).
Chen, C., Bao, J. & Zhu, C. Molecular structural parameters and thermal stability of 1-methylcarbazole. Geol. Sci. Tech. Infor. 25(2), 57–59 (2006).
We are grateful to Drs. Quan Shi and Minghui Wu who kindly provided assistance in oil sampling and analyses. Special thanks to Professor Lloyd Snowdon for his reviewing and English editing of the manuscript before submission. We also thank Professor Simon George, three anonymous reviewers and the journal editor for their insightful and constructive comments and suggestions. This research was supported by Strategic Priority Research Program of Chinese Academy of Sciences (No. XDB10020203), China National Major Science 973 Projects (Nos. 2003CB214608 and 2014CB239101) and Guangdong Province Higher Education Pearl River Scholar Program.
Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, 100029, China
L. Zhang
Innovation Academy for Earth Science, Chinese Academy of Sciences, Beijing, 100029, China
Department of Earth, Ocean and Atmospheric Science, Florida State University, and National High Magnetic Field Laboratory, Tallahassee, FL, 32306-4100, USA
Y. Wang
Institute of Groundwater and Earth Sciences, Jinan University, Guangzhou, P.R. China
Sinopec Key Laboratory of Petroleum Accumulation Mechanisms, Sinopec Research Institute of Petroleum Exploration and Production, 31 Xueyuan Road, Beijing, 100083, China
M.-W. Li
Department of Earth and Planetary Sciences, University of California at Davis, One Shields Avenue, Davis, CA, 95616, USA
Q.-Z. Yin
Changqing Oilfield Company, PetroChina, Xi'an, Shanxi, 710021, China
W. Zhang
Search for L. Zhang in:
Search for Y. Wang in:
Search for M.-W. Li in:
Search for Q.-Z. Yin in:
Search for W. Zhang in:
L.Z. designed research, originated the concept, developed its computation method. L.Z., Y.W., M. L. and Q.-Z. Y. analyzed the data and wrote the paper. W.Z. analyzed the data.
Correspondence to L. Zhang.
The authors declare no competing interests.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Zhang, L., Wang, Y., Li, M. et al. Relative sorption coefficient: Key to tracing petroleum migration and other subsurface fluids. Sci Rep 9, 16845 (2019) doi:10.1038/s41598-019-52259-6
Received: 30 November 2018
DOI: https://doi.org/10.1038/s41598-019-52259-6
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports menu
About Scientific Reports
Guest Edited Collections
Scientific Reports Top 100 2017
Scientific Reports Top 10 2018
Editorial Board Highlights
Author Highlights