Abstract

The prediction of oil cracked gas resources is necessary and urgent in the gas exploration of these basins at high to over stage in China. A marine crude oil sample was pyrolyzed using sealed gold tubes system in our study. The pyrolysates including gas, liquid, and solid were quantitatively analyzed. Based on the pyrolysis data and kinetic calculation, the yield correlativity among gas, liquid, and solid products was regressed with high correlative coefficients to establish a prediction model suitable for the resource estimation of oil cracked gas. The verification formula for this model was also established on the principle of mass conservation. The affecting factors and the application preconditions of this model were discussed. Finally the model was extrapolated to the prediction of oil cracked gas resources of Feixianguan formation in NE Sichuan basin, SW China. The prediction value of oil cracked resources is about  m3, and generation intensity of oil cracked gas is about  m3/km2, and the paleo-oil reserves is about  t. The verifying value for this prediction is approximately equal to 1, indicating the model is reliable in the resource estimation of oil cracked gas.

1. Introduction

The presentation of oil cracked gas and the research of kinetic model of oil cracking have made contribution to the theory of natural gas generation [16]. It is shown that oil cracking is an important passway for gas generation from marine source rocks [7]. Oil cracked is a pyrolysis process to generate gaseous hydrocarbons and residues (solid bitumen) under certain temperature and pressure, which is similar to the generation of oil and gas from Kerogen. This process can be described by chemical kinetic equations. Kinetic pyrolysis experiment is an effective method to simulate it [8]. Natural gas should be mainly from oil cracked at the high-over mature stage of marine source rocks according to normal kinetic model [9]. The exploration of oil cracked gas may show a bright prospect because of the great thickness and wide distribution of marine sediment in China. It is important to predict resources of oil cracked gas accurately and effectively.

Conventional oil and gas resource assessment is normally based on the identification and evaluation of effective source rocks and the estimation of the potential of hydrocarbon generation [10, 11]. However, it is difficult to draw reliable conclusions by applying traditional resource assessment parameters to the prediction of oil cracked gas resource since it is generally distributed in the lower part of superposition basins. The mechanism of the hydrocarbon distribution and accumulation in superposition basin is quite complicated [11], and the source rocks have evolved to high and over mature stage [12, 13]. It is indicated in the exploration that the solid bitumen in reservoir in Sichuan Basin, Tarim Basin, and the extensive carbonate areas in south China is pyrobitumen formed in the process of oil cracked in paleo-oil reservoir [14, 15]. The pyrobitumen is an important indicator of oil cracked gas. Meanwhile, it can be a new key to the resource assessment of oil cracked gas as an alternative of hydrocarbon source rock.

Previous studies on crude oil pyrolysis mainly focused on the mechanism of gas generation, the evolution of gaseous hydrocarbons and isotope composition [16, 8, 16, 17]. Less attention has been paid to the quantity changes of gas, liquid, and solid pyrolysates and the correlativity among them in pyrolysis. The correlation was not yet extrapolated to resource assessment at present.

In the process of oil pyrolysis in a closed system, it was noted that the yields of solid pyrobitumen and gaseous hydrocarbons continued to increase (a positive correlation between them), while the yield of liquid hydrocarbons decreased. Inspired from that, a prediction model for the resource assessment of oil cracked gas based on solid bitumen in reservoir was established according to the pyrolysis data and mass balance calculation [18]. Approach was provided for the estimation of natural gas resources and the restoration of paleo-oil accumulation in the areas with high-over mature marine source rocks in south China, furthermore, its application in the gas pools of Feixianguan formation is discussed here. It is significant for the decision-making of natural gas exploration in China.

2. Sample and Experiment

2.1. Sample

A normal marine crude oil sample was selected for pyrolysis experiment from Ordovician interior reservoirs in the depth of 5940–5983 m in Well YM2 in Tarim Basin. The oil sample contains 47.3% saturates, 26.4% aromatics, 13% nonhydrocarbons, and 5.95% asphaltenes. The ratio of saturates to aromatics is 1.79.

2.2. Experiment

Oil cracked generally occurred in a relatively closed paleo-oil reservoir in geological condition, so a pressurized confined gold tube system was adopted in our experiment. The equipment and principle were presented elsewhere in detail [12, 17, 19]. Here the experimental procedures are briefly described as follows. The prepared oil sample (5–40 mg) was sealed in gold tubes (40 mm × 4.2 mm i.d.) and charged with argon. All of the gold tubes were placed in stainless steel autoclaves quickly placed in the oven. The pressure of autoclaves was kept at 50 MPa. The temperature were programmed to increase at a rate of 2°C/h and 20°C/h, respectively, from 350°C to 600°C. For each of the two heating series, oil samples were loaded into gold tubes under anhydrous conditions. All autoclaves were heated in a uniform oven. At last, the autoclaves were removed from the oven one by one and then quenched in a cold water bath when the desired temperatures were reached. During the experimental heating process, the temperature error is less than 1°C, and the pressure error is less than 5 MPa.

2.2.1. Gas Analysis

Cleaned gold tubes were placed in a vacuum line connected to an Agilent 6890N gas chromatography. After isolating the extraction line from the vacuum pump, the tube was pierced with a needle, allowing the gas to be volatilized into the line. Then, molecular characterization and quantification of the total gas fraction were performed by gas chromatography. External standard method was applied to the quantification of gas components. The GC conditions were as follows. The column was kept at 40°C for 6 min, programmed to 180°C at a rate of 25°C/min, and finally kept at 180°C for 4 min.

Residual liquid hydrocarbons comprise light hydrocarbons (C6–C13) and heavy hydrocarbons (C13+). Following gas analysis, the light hydrocarbons were frozen in a vial (4 mL) by liquid nitrogen and quickly added into pentane solution. The punctured gold tubes were removed from the vacuum device, cut open, and added to the vial. Then it was extracted three times by sonication in pentane and taken out. The internal standard of deuterated-nC18 was added into the solution in vial which was placed for hours after constant volume. A part of the solution (200 μL) was taken out for the quantitative analysis of light hydrocarbons. The remaining part was filtered using repeated washing with dichloromethane. The filtrates were evaporated and weighted. The C13+ fraction was then measured. The light hydrocarbons analysis was carried out by Finnigan Trace GC Ultra gas chromatography and Quadax column (50 m × 0.25 mm). The GC conditions were as follows. The column was kept at 40°C for 6 min, programmed to 180°C at a rate of 25°C/min, and finally kept at 180°C for 4 min.

Residual liquid hydrocarbons comprise light hydrocarbons (C6–C13) and heavy hydrocarbons (C13+). Following gas analysis, the light hydrocarbons were frozen in a vial (4 mL) by liquid nitrogen and quickly added into pentane solution. The punctured gold tubes were removed from the vacuum device, cut open, and added to the vial. Then it was extracted three times by sonication in pentane and taken out. The internal standard of deuterated-nC18 was added into the solution in vial which was placed for hours after constant volume. A part of the solution (200 μL) was taken out for the quantitative analysis of light hydrocarbons. The remaining part was filtered by dichloromethane again and again. The filtrates were evaporated and weighted. The C13+ fraction then can be measured. The light hydrocarbons analysis was carried out by Finnigan Trace GC Ultra gas chromatography and Quadax column (50 m × 0.25 mm). The GC conditions were as follows: the column was kept at 50°C for 2 min, programmed to 290°C at a rate of 4°C/min, and finally kept at 290°C for 30 min. Thus, the total fraction can be accurately measured.

Organic filter papers (0.45 μm × 25 mm) were used to recover the pyrobitumen. Each of the filter papers was weighed beforehand. After being filtered, the filter papers together with pyrobitumen were dried and weighed again to determine the amount of pyrobitumen. All these operations were done in the same experimental condition (temperature 20°C; humidity 45%) and with the same electric balance (systematical error no more than 0.003 mg) to minimize experimental errors.

3. Results and Discussion

3.1. Changes in the Yields of Gas, Liquid, and Solid Pyrolysates

The total gaseous products, residual liquid hydrocarbon, and pyrobitumen generated from pyrolysis are discussed here. The volume and the quality of total gaseous products are the sum of plus H2 and CO2, and so forth. The quality of the residual liquid hydrocarbon is the sum of and C13+.

The yields of total gas (volume yield) and pyrobitumen increase continuously along with the rise of pyrolysis temperature (the two yield curves resemble each other), while the yield of residual liquid hydrocarbon decreases (Figure 1 and Table 1).

At the heating rate of 2°C/h, crude oil cracked began at about 380°C with low yields of total gas volume (20 mL/g) and pyrobitumen (9 mg/g). The two yields increased quickly to 684 mL/g and 335 mg/g, respectively, from 420°C to 500°C and then slowly increased above 500°C. It reached the maximum value (796 mL/g, 429 mg/g), respectively, by 600°C. By contrast, the yield of residual liquid hydrocarbons decreased quickly to 54 mg/g below 500°C and then went slowly to a value of zero. The result implies that the oil cracked process could be divided into two stages with possibly different mechanisms.

As shown in Figure 1, the change trends of the yield curves of all pyrolysates under two heating rates are similar. The yields under the low heating rate (2°C/h) lag behind those under the high heating rate (20°C/h). It indicates that a compensation relationship exists between time and temperature during the chemical process of oil pyrolysis [20].

A software for kinetic analysis was utilized for the data processing in order to eliminate the influence upon the oil pyrolysis experiment under different heating rates and discuss the characteristics in different evolution phases of oil cracked. The method showing thermal maturity expressed as equivalent vitrinite reflectance (%) by Sweeney was adopted [21]. The experimental data and the kinetic parameters from [21] are processed with the software for kinetic analysis. The corresponding to different temperatures can be derived. The Easy% value corresponding to different temperatures can be derived according to the formula: Easy. Then the changes in the yields of crude oil pyrolysates along with Easy (Figure 2) can be established. As shown in Figure 2, the trend of yield curves based on kinetic calculation data is approximately consistent with that derived from experimental data. The influence under different heating rates is so unconspicuous that the yield curves under different heating rates almost overlap each other. The yields of gas and pyrobitumen increase, whereas liquid hydrocarbons yield decreases quickly at about Easy  0.8%–1.6%. All the yield curves extend gently as far as Easy 4.4. The ultimate pyrolysates are methane and coke residue [8].

3.2. Yield Correlativity among Gas, Liquid, and Solid Pyrolysates

The increase or decrease of gas, liquid, and solid matter in the process of oil cracked can be shown as the correlativity among the three yields. Based on the kinetic calculation and experimental data, the correlativity among the yields of total gas volume and residual liquid hydrocarbons as well as solid pyrobitumen was regressed with high correlative coefficients in this study (Figure 3).

There is linear positive correlativity between the yields of total gas volume and pyrobitumen (Figure 3(a)). It indicates that the two yields increase continuously and simultaneously. In other words, when an amount of pyrobitumen was generated from crude oil pyrolysis, a corresponding volume of gas was generated. The rate of the two yields should be close to a constant (slope in the liner formula). The yield correlativity (quadratic function) between pyrobitumen and residual liquid hydrocarbons resemble that between total gas volume and residual liquid hydrocarbons (Figures 3(b) and 3(c)). It indicates that the residual liquid hydrocarbons decreased (ultimately with only traces left) following a quadratic function with the increase of the total gas volume and pyrobitumen. It is important to point out that the regressed correlation formula only can be applied within certain limits and the values of variable in the formula are not more than the maximum yield of all kinds of pyrolysates. Furthermore, Figure 3 shows that there is little difference between the constant coefficients in the correlation formula regressed from kinetic calculation and that from experimental data. It indicates that the correlation formula is reliable in theory.

3.3. Mass Conservation in the Process of Crude Oil Pyrolysis

The total quality was conserved in the process of oil pyrolysis, no matter how gas, liquid, and solid pyrolysates changed, respectively, in quality because a confined pyrolysis system close to the geological conditions of the plaeo-oil in reservoir cracked was adopted. At any time, the sum of the quality yields of all pyrolysates should be equal or approximately equal to 1. The volume yield of gas is converted to quality yield according to the composition of gaseous pyrolysates. The mass conservation is discussed through the changes in quality of gas, liquid, and solid pyrolysates in the process of oil pyrolysis (Figure 4).

At the beginning of oil pyrolysis, the quality of total gas and pyrobitumen increased continuously until the quality percentage of total gas went up to 65% where an inflexion appeared (Figure 4(a)). Then the quality of total gas began to decrease but the quality of pyrobitumen increased continuously. The quality percentage of the ultimate pyrolysates was 55% total gas and 45% pyrobitumen. In the later period of oil pyrolysis, the gas hydrocarbons cracked continuously to methane and pyrobitumen, and the percentage of methane in pyrolysates increased. Although the volume of total gas continued to increase, the quality reduced at the same time. The quality reduction of total gas was transferred to pyrobitumen. The quality yield curves of gas, liquid, and solid pyrolysates based on kinetic calculation data under the two heating rates coincided with each other approximately (Figure 4(b)). The curves a and (quality yield of total gas) began to descend gently (from 0.65 to ultimate value 0.55) at around over 2.0 of Easy where an inflexion appeared. The quality yield curves of pyrobitumen ( and ) and residual liquid hydrocarbons ( and ) show the same change trend with that shown in Figure 2. The value of mass balance curves ( and ) was close to 1 all the while. It indicates that the total quality is approximately constant in oil pyrolysis in confined system. Therefore, at any value of Easy the relationships and exist. The values of d and are equal or approximately equal to 1. It not only proves the mass conservation in the process of crude oil pyrolysis, but also certifies the yield correlation among gas, liquid, and solid pyrolysates.

4. Prediction Model of Oil Cracked Gas Resources

4.1. Establishment of the Prediction Model

As shown by the correlativity formulas regressed in Figure 3, the increase or decrease of gas, liquid, and solid products in the process of oil pyrolysis in confined system can be expressed by formulas with variables. If the value of quantity or yield of any phase pyrolysates is given, the values of the other two phases can be got. The correlation can be extrapolated to practical geological conditions. If the quality of solid bitumen in reservoir is given, the resources of the corresponding natural gas (oil cracked) and paleo-oil in reservoir can be determined. Therefore, a prediction model for the resource estimation of oil cracked gas is established here. The variables in the model are listed in Table 2.

The influence exerted by two different heating rates and the errors in pyrolysis can be eliminated by kinetic calculation. The results of that can be extrapolated to geological application. Therefore, the formulas regressed by kinetic calculation data in Figure 3 are adopted to express the yield correlation among gas, liquid, and solid pyrolysates in oil pyrolysis. The variables in Table 1 are substituted into them to get the formulas as follows: where .

Since the cracking of the paleo-oil in reservoir occurred in a closed geological system, mass conservation should exist in oil pyrolysis. The relevant variables in Table 1 are put into the expression of mass conservation in Figure 4(b) to get the verification formula as follows:

According to the natural gas resources (volume) derived from prediction model and its composition from GC analysis, the quality of natural gas can be calculated out and put into formula (2). If the sums are equal or approximately equal to 1, the prediction value of natural gas resources would be reliable in theory.

4.2. Factors Influencing of the Prediction Model

Because the prediction model is established based on the pyrolysis result of only one crude oil sample, some influence factors may exist as follows.(1)Property of Crude Oil. The cracking of crude oil with different properties (such as percentage of saturates, aromatics, nonhydrocarbons, and asphaltenes) may produce different yields (such as gas or pyrobitumen). For example, if a heavy crude oil with high asphaltenes cracked, the pyrobitumen yield will be high, and the quality proportion of the ultimate products (methane and coke) will be different from that of normal crude oil. Moreover, the special elements content (such as sulfur), viscosity, density and other properties of crude oil will also affect the rate of cracking and the yield of pyrolysates.(2)Measurement Errors in Experiment. Measurement errors have direct impact on the results of oil pyrolysis. For example, the factors such as the weight and water absorption of the filters will make the quantitative value higher or lower in the course of determining the amount of pyrobitumen. If the filter papers are not washed clearly by dichloromethane, the macromolecular asphaltenes retained on the filter paper will make the quantitative value of pyrobitumen higher. In addition, the volatilization of gas and light hydrocarbons will make the quality value of residual liquid hydrocarbons and gas volume lower.(3)Some secondary alterations of the natural gas or oil could produce an effect on the sensitivity of prediction model, such as TSR [24, 25] and the interaction between the minerals of reservoir rocks and crude oil.

However, the factors described above affect only the value of products yield. It cannot change the rules in pyrolysis. Perhaps, the value of constant coefficient in the correlativity formulas of prediction model can fluctuate a little because of these factors, but the linearity and quadratic function maintain.

4.3. Application Preconditions of the Prediction Model

The resources of oil cracked gas can be predicted effectively using this model with the advantages as follows. Resource evaluation can be taken according to solid bitumen in reservoir instead of source rocks; the restoring of hydrocarbon potential and the evaluation of high mature source rocks are avoided. However, the model may not be suitable for the prediction of resources of all kinds of natural gas accumulations. Some preconditions should be met with in practice. (1) The most important application precondition is that the natural gas should be generated from oil cracking. For the gas accumulations dominated by oil cracking but mixed with gas from other origins, the prediction value by this model will be lower than the factual value. (2) The generation, accumulation, and preservation of oil cracked gas should occur in a relatively closed geological system. In other words, tectonic activities did not bring large destruction on the accumulations when oil is being cracked or after gas accumulated. Meanwhile, the natural gas resources lost little in the process of accumulation. Otherwise, the prediction value by this model will be higher than the factual value of the resources of present gas accumulations.

5. A Case of Geological Application

This model was extraporlated to the oil cracked gas pools of Feixianguan Formation in NE Sichuan basin to verify the reliability and practicability of the prediction model.

5.1. Geological Setting

The NE Sichuan Basin was part of a carbonate platform during the Feixianguan (earliest Triassic). Multiple phases of tectonic deformation during the Yanshanian and Himalayan orogenies resulting in the development of the Dabashan foldbelt to the NE and the Eastern Sichuan foldbelt to the SE (Figures 5(a) and 5(b)). Carbonate platform strata have therefore been intensely deformed and the Feixianguan Formation was buried to depths of 3500 m to 6000 m [26].

Seven large natural gas accumulations have been discovered recently in the Lower Triassic Feixianguan Formation of the northeastern Sichuan Basin, located in the Luojiazhai, Dukouhe, Tieshanpo, Puguang, Jinzhuping, Huanglongchang, and Qilibei structures, respectively (Figure 5(c)). The gas pay zones are distributed under the regional caprock beds of the Lower Triassic Jialingjiang Formation, and almost all of the gas traps are associated with faults and fractures. Vertically, the highest-quality oolitic reservoirs are always within the T1f 3 to T1f 1 members. Laterally, reservoir qualities are comparable in four of the six structures, with an average 13–17% porosity. The only exception is observed in the Jinzhuping field where the average porosity is only 4.3%. The solid bitumen that occurred commonly in oolitic reservoirs is the residual product resulting from crude oil cracking [13, 22, 27, 28].

Previous studies have proved that the gas pools have experienced oil reservoir phase and formed ultimately without any large-scale destruction or dissipation, and oil cracking occurred in a relatively closed geological system [22, 23]. So the oolitic gas pools in NE Sichuan Basin are perfect area to extrapolate the resources prediction model established in our study.

5.2. Prediction of Oil Cracked Gas

As shown as Figure 6, a contour map was provided to show the content of solid bitumen in reservoir of Feixianguan Formation in NE Sichuan Basin [22]. According to the contours, weight average content of solid bitumen in reservoir is about 6.23%, with distribution area of 700 km2. Firstly, the quality of solid bitumen in gas reservoir () was estimated according to formulas where is density of solid bitumen in reservoir, general value adopted range in 1.36~1.98 g/cm3 [29], 1.65 g/cm3 was adopted here; is trap area, 700 km2 was adopted here; is efficient thickness of reservoir, average value 60 m was adopted; is volume percentage of solid bitumen in reservoir, average value 6.23% was adopted here; is efficient volume of reservoir; is volume of solid bitumen in reservoir. So the quality of solid bitumen () 3.78 × 1018 mg was derived according to formulas (3).

When the value (3.78 × 1018 mg) is inserted into formula (1), the volume of natural gas resources is equal to 6.84 × 1012 m3, oil cracked gas is 97.5 × 108 m3/km2, and the paleo-oil reserve is 97 × 108 t.

The volume () of natural gas resources was converted into the relative percentage according to where is the natural gas resources (volume), value of 6.84 × 1012 m3, is mol mass of deferent composition (such as CH4, and C2H6) in natural gas, based on an average calculation of the composition data of natural gas given in Table 3. The value of range in 5.78 × 1018~6.29 × 1018 mg was acquired and put into formula (2). The value of range in 0.957~0.985, which is approximately equal to, and a little less than, 1 and indicates that the prediction value of natural gas resources is reliable in theory.

6. Conclusions

(1)The increase or decrease of gas, liquid, and solid matter residue in the process of crude oil pyrolysis in confined high-pressure system can be shown as the correlation among the yields of three pyrolysates. There is linear positive correlativity between the yields of total gas volume and pyrobitumen. The relationship between the yields of pyrobitumen and residual liquid hydrocarbons is a quadratic negative function correlation, the same as the relationship between total gas volume and residual liquid hydrocarbons. A prediction model based on solid bitumen in reservoir is suitable for the resource estimation of oil cracked gas. The verification formula is established on the basis of mass conservation. If the amount of solid bitumen in a reservoir is given, the resource of oil cracked gas and paleo-oil accumulation can be predicted accurately to some extent.(2)The prediction model was extrapolated to the gas pools of Feixianguan formation in northeastern Sichuan basin. The prediction value of oil cracked resources is about 6.84 × 1012 m3, and generation intensity of oil cracked gas is about 97.5 × 108 m3/km2, and the paleo-oil reserves is 97 × 108 t. The verifying value for this prediction (0.957~0.985) is approximately equal to 1, indicating that the model is reliable for application in the resource estimation of oil cracked gas.(3)It is the first attempt to evaluate oil cracked gas resources through analyzing the increase or decrease of gas, liquid, and solid products in crude oil pyrolysis and the correlativity among the three yields. Some new ideas are provided for the estimation of natural gas resources and the restoration of paleo-oil accumulations in the areas with high and over mature marine source rocks in south China. It is significant for the decision-making of natural gas exploration in China.

Acknowledgments

This paper was financially supported by the National Basic Research Program of China (no. 2007CB209500), the Youth Fund of Research Institute of Petroleum Exploration and Development, PetroChina (QY012201), and the science and technology innovation Fund of PetroChina (2010D-5006-0102). The authors are grateful to Professor Sun Yongge in Guangzhou Institute of Geochemistry for providing crude oil samples and Professor Liu Jinzhong, Liao Zewen, and Senior Engineer Xu An for their help in pyrolysis experiment.