The organopores play an important role in determining total volume of hydrocarbons in shale gas reservoir. The Lower Silurian Longmaxi Shale in southeast Chongqing was selected as a case to confirm the contribution of organopores (microscale and nanoscale pores within organic matters in shale) formed by hydrocarbon generation to total volume of hydrocarbons in shale gas reservoir. Using the material balance principle combined with chemical kinetics methods, an evaluation model of organoporosity for shale gas reservoirs was established. The results indicate that there are four important model parameters to consider when evaluating organoporosity in shale: the original organic carbon (w(TOC0)), the original hydrogen index (IH0), the transformation ratio of generated hydrocarbon (F(Ro)), and the organopore correction coefficient (C). The organoporosity of the Lower Silurian Longmaxi Shale in the Py1 well is from 0.20 to 2.76%, and the average value is 1.25%. The organoporosity variation trends and the residual organic carbon of Longmaxi Shale are consistent in section. The residual organic carbon is indicative of the relative levels of organoporosity, while the samples are in the same shale reservoirs with similar buried depths.
1. Introduction
Petrological characteristics (such as brittle mineral content), micro-nanoscale porosity, total organic carbon, and organic matter maturity of shale play important roles in the accumulation and exploration of shale gas. The nanometre pores in shale are mostly organopores (pores within organic matters in shale) created from hydrocarbon generation during formation subsidence and thermal evolution [1–3]. Many scholars studying the organoporosity of shale gas reservoirs believe that the organopores in organic matters generated during hydrocarbon generation make a significant contribution to gas reservoir space [4–7]. For example, in the Mowry Shale of the Powder River Basin, the organic matter type of the source rocks is type II. If the original organic carbon w(TOC0) is 6% and the vitrinite reflectance (Ro) goes up to 1.2%, the organoporosity is approximately 5% [4, 5]. The surface area of the organopores that approach nanoscale becomes significant and strongly affects the amount of gas stored as the free component versus the adsorbed component [8–10]. Typical methods, such as mercury injection capillary pressure (MICP) measurement and nuclear magnetic resonance (NMR) spectroscopy, can not measure organopores in shale effectively. These macroscopic methods do not involve the direct measurement of individual pores, which range in size for shale gas from a few to tens of nanometres [11, 12]. At present, the methods for measuring porosity of shale are low-temperature N2 and CO2 gas isotherm adsorption-desorption, scanning electron microscopy (SEM), 3D high-resolution imagery, and total porosity determined by bulk density coupled with skeletal density. Isotherm adsorption-desorption measures the connected pores, which range in size from 0.38 to 100 nanometres in shale, and the measured results are smaller than the actual pore values [13]. Scanning electron microscopy can image the microstructure of shale pores, especially the organopores [14–18], but the observations are limited to 2D sections and cannot be used to perform quantitative evaluation. The difference between skeletal density and bulk density measured by He pycnometry and Hg immersion, respectively, is used to calculate the total porosity of shale [19, 20]. Because of the adsorption of organic-rich shale, inconsistent crushing and sample handling procedures, the difference of relative humidity of the measurement environment, and the solvent extraction pretreatment that may potentially remove some of the solid organic matter, the measured volume has certain errors and defects [20]. Three-dimensional high-resolution imagery can be used to investigate shale microstructure and analyse the characteristics of pores and the porosity of microscale shale cubes [8, 9]. However, smaller sample diameters correspond to higher 3D-reconstruction resolutions [21]. The particle size of quantitatively evaluated shale is 1 to 5 μm in 3D-reconstruction [22]. There is strong heterogeneity in shale [23]. Therefore, smaller samples are less representative. What is more, these methods can not evaluate the organoporosity effectively.
Given that organopores play an important role in determining shale gas reservoir space, it is necessary to study organoporosity. Therefore, taking the Lower Siluriran Longmaxi Shale in the southeast Chongqing as a case study, a method was established for estimating organoporosity, which is useful in the simultaneous exploration and development of shale gas.
2. Geological Setting
The study area is in southeast Chongqing, China, and covers an area of approximately 1.98 × 104 km2. The north, south, and east of it are Hubei Province, Guizhou Province, and Hunan Province, respectively. It belongs to the Yangtze tectonic plate [24, 25] and is located in the Wuling Drape Zone and Hunan-west Hubei Drape Zone. The Xuefengshan Uplift and Sichuan Basin are in the east and northwest of southern Chongqing, respectively (Figure 1).
Stratigraphic column of Py1 well and isopach of Longmaxi Shale in southeast Chongqing.
The residual strata of Paleozoic are Cambrian, Ordovician, and Silurian, and the other layers are denuded or missing. The Lower Paleozoic marine shale mainly including the Lower Cambrian Niutitang Shale and the Lower Silurian Longmaxi Shale is widely deposited in southeast Chongqing. They have the characteristics of large thickness, deep buried depth, and well-developed fractures [26–28]. In the Longmaxi formation, the lower and middle parts belonged to deep-marine shelf sedimentary environment and the upper part to shallow-marine shelf sedimentary environment. Grey and black carbonaceous shale, shale, and grey shaly sandstone are in the lower part, grey and black mudstone is found in the middle part, grey lime mudstone, and shale and shaly sandstone are present in the upper part. Grey and black carbonaceous shale and shale in the lower part of Longmaxi formation with a range of approximately 30–100 metres called Longmaxi Shale are the target strata of shale gas exploration and development (Figure 1). In the study area, the thickness of Longmaxi shale in Py1 well is about 70 metres with a total organic carbon w(TOC) range of 0.27–4.25% and average value is 1.90%. The organic matter are types II1-I [29], with a vitrinite reflectance (Ro) range of 1.90–3.09%, and average value is 2.62%, which are calculated by bitumen reflectance [30].
3. Model Creation
Organopores are within organic matters in shale, due to hydrocarbon generation during formation subsidence and thermal evolution. These pores are important for storage and perhaps transfer of gas molecules through shale [23]. The diameter of these organopores is generally less than 1 μm. Thus, they are considered nanometre pores. The sizes of nanometre pores in organic matter mainly range from 80 to 100 nm [17, 22]. For example, in Barnett Shale, a single 10.8-μm-diameter organic matter contained more than 1,000 nanopores of various shapes and sizes [1]. The diameters of intragranular nanopores do not appear to be directly related to the size of the organic matters. Both large and small nanopores can be found in organic matters of all sizes. The total volume of the organopores is mainly related to the total organic carbon level, hydrogen index, maturity, and diagenesis (compaction effect). The paper will discuss how to evaluate the organoporosity of shale.
According to the material balance principle, under ideal conditions, the total volume of organopores is the volume of organic matter consumed by generating hydrocarbon (Formula (1)) [5]:
(1)Фorganic′=w(TOC0)·IH0·F(Ro)·ρrock/ρkerogen1000,
where Фorganic′ is the organoporosity of shale under ideal conditions (%), w(TOC0) is the weight percent of the original total organic carbon (%), IH0 is the original cracking hydrocarbon of unit quality organic carbon (mg/g), F(Ro) is the transformation ratio of oil and gas generated from organic matter (%), which is correlated with maturity, ρrock is the density of shale (g/cm3), and ρkerogen is the density of kerogen, which is approximately 1.2 g/cm3 [16, 31]. The term w(TOC0)·IH0 represents the weight of convertible organic matter. Then the term w(TOC0)·IH0·F(Ro) is the weight of organic matter consumed by generating hydrocarbon in thermal evolution. The term ρrock/ρkerogen is simply a density conversion from weight percent to volume percent. The last term 1000 is a unit conversion from mg to g.
Because of compaction and bituminization process in actual geological conditions, the organoporosity calculated under ideal conditions should be corrected. On the basis of Formula (1), a model was established for evaluating the organoporosity in shale as follows (Formula (2)):
(2)Фorganic=w(TOC0)·IH0·F(Ro)·ρrock/ρkerogen1000·C,
where Фorganic is the organoporosity of shale in actual geological conditions (%) and C is the organopore correction coefficient.
In the model, there are four important parameters: the original total organic carbon (w(TOC0)), the original hydrogen index (IH0), the transformation ratio of oil and gas generated from organic matter (F(Ro)), and the organopore correction coefficient (C). This paper will evaluate these parameters in the following part.
4. Model Parameter Evaluation4.1. Transformation Ratio of Hydrocarbon Generated from Organic Matter
To calculate the transformation ratio of hydrocarbon generated from organic matter in Longmaxi Shale of southeast Chongqing, the Xiamaling Shale in the Huabei Platform and crude oil sample generated by the Cambrian source rocks in Tarim Basin were selected for pyrolysis experiments because the shale is similar to the source rocks of the Longmaxi formation and has a relatively lower maturity. Xiamaling Shale is type II1 with a vitrinite reflectance (Ro) of approximately 0.5%, and the w(TOC) is 5.98% (Table 1). The maturity of Longmaxi Shale is relatively high with average vitrinite reflectance 2.62% in study area, and generated oil migrated to overlying strata was suffered bituminization effect. There is no proper crude oil generated from Longmaxi Shale in Sichuan Basin, leaving nature gas and asphaltite. So this paper selected the crude oil sample from the Silurian oil reservoir in the Tarim Basin (Table 2), which is generated by the Cambrian source rocks [32]. The source rocks of Cambrian and Silurian are both marine sedimentary of Paleozoic. The crude oil generated from them may have a certain degree of similarity.
Geochemical characteristics of the shale samples.
Samples
Organic matter type
Vitrinite reflectance Ro, (%)
Total organic carbon w(TOC0), (%)
Max temperature Tmax, (°C)
Free hydrocarbon S1, (mg/g)
Cracking hydrocarbon S2, (mg/g)
Hydrogen Index IH, (mg/g)
Sedimentary period
Xiamaling Shale
II1
0.5
5.98
439
0.45
19.13
320
Precambrian
Longmaxi Shale
II1-I
2.62
1.90
455
0.01
0.31
1.01
Lower Silurian
Characteristics of the crude oil sample from Cambrian source rocks in Tz62 well.
Source rock
Reservoir
Sampling depth (m)
Density g/cm3
Saturated (%)
Aromatic (%)
Nonhydrocarbon and asphaltene (%)
Cambrian
Silurian
4052.88~4073.58
0.931
53.256
33.285
13.459
Pyrolysis experiments were performed from 200°C to 600°C with different heating rates using a Rock-Eval-II type pyrolysis instrument (e.g., 30°C/h and 40°C/h). The relationship between the production of hydrocarbons and temperature (or time) was recorded. Next, the transformation ratio of oil and gas at different temperature points was obtained, and the chemical kinetics parameters were calibrated (Figure 2). Further reference studies can be consulted to understand the calibration method [33–35].
The transformation ratio of oil and gas generated from organic matter at different temperatures and heating rates. (a) The transformation ratio of oil generated from kerogen of Xiamaling Shale. (b) The transformation ratio of gas generated from kerogen of Xiamaling Shale. (c) The transformation ratio of gas generated from oil cracking.
On the basis of previous research [36–38], the burial history and thermal history of the study area were analyzed. In Py1 well area, Longmaxi Shale experienced multiple stages of tectonic activity in the southeast of Chongqing. It mainly underwent two rapid subsidence (approximately 450–400 Ma and 285–85 Ma), one slow uplift (approximately 400–285 Ma), and one rapid uplift (85 Ma to now). The maximum buried depth of the shale was approximately 4900 m, and the highest temperature was approximately 160°C (Figure 3). Based on the burial and thermal history, calibration of the chemical kinetic parameters and clarification of the relationship between the transformation ratio of hydrocarbon generated from organic matter (including the transformation ratio of oil and gas generated from kerogen and the transformation ratio of gas generated from oil cracking) and geological period (Figure 4) were combined for the model. In Py1 well area of southeast Chongqing, the main oil and gas generating periods of Longmaxi Shale are 275–190 Ma (Middle Permian to Late Triassic) and 260–80 Ma (Late Permian to Middle Cretaceous).
The burial and thermal history of the Py1 well area in southeast Chongqing.
The transformation ratio of hydrocarbon generation with geological period from Py1 well Longmaxi Shale in southeast Chongqing.
4.2. Original Hydrogen Index and Original Total Organic Carbon
Based on the restoration method of original hydrogen index and original total organic carbon [39] and using the transformation ratio of hydrocarbon generation in combination with the burial and thermal history of the study area (Figure 3), the original hydrogen index value and original total organic carbon of Longmaxi Shale from Py1 well were calculated.
The original hydrogen index value was calculated using the following formula (Formula (3)) [33]:
(3)IH0=IH+(IH0·Fo+B0-B)+IH0·Fg,
where IH0 is the original hydrogen index value (mg/g), IH is the residual hydrogen index value (mg/g), Fo is the transformation ratio of oil generation (%), Bo is the content of native asphalt in source rocks (non-thermally decomposed from kerogen, mg), B is the residue oil content calculated from chloroform “A” or the hydrocarbon index by light and heavy hydrocarbon compensation, respectively (mg), and Fg is the transformation ratio of gas generation (%).
The original total organic carbon was calculated using the following formula (Formula (4)) [33]:
(4)w(TOC0)=w(TOC)(1+ΔIH·K1000),
where w(TOC) is the weight percentage of the residual total organic carbon (%), ΔIH is the recovery content of the hydrogen index (mg/g), and K is the coefficient of organic matter conversion into organic carbon, which is approximately 0.85 [5].
The original and residual hydrogen indexes of Longmaxi Shale from Py1 well are in range of 421.31–446.03 mg/g and 0.02–4.93 mg/g, with average values 427.12 mg/g and 1.01 mg/g. The ranges of original and residual total organic carbon are 0.74–10.86% and 0.27–4.25% with average values 4.90% and 1.90%. The original organic carbon is significantly larger than residual organic carbon, but their variation trends with depth are consistent (Figure 5).
Hydrogen index and total organic carbon values of the Py1 well from Longmaxi Shale in southeast Chongqing. (a) Hydrogen index. (b) Total organic carbon.
4.3. Organopore Correction Coefficient
The organopore correction coefficient is the ratio of the volume of actual organopores to pore space, which is formed from hydrocarbon generation under ideal conditions without compaction and bituminization effects in organic matters. In the same set of shale with similar buried depth, because of the environment of temperature and pressure without difference, the process of hydrocarbon generation and expulsion and diagenesis is nearly homogeneous. So the organopore correction coefficients of them are the same. In order to obtain the organopore correction coefficient, the focused ion beam (FIB) milling and scanning electron microscopy (SEM) imaging were used on the shale samples at the same depth. The average organopore surface porosity (SPorganic) was obtained using the total area of organopores divided by the area of organic matters in several images with the same depth. Organoporosity (Фorganic) was calculated through combining with the total organic carbon, kerogen density, and bulk density of the samples using Formula (5). Then, using organoporosity (Фorganic) ratio to organoporosity (Фorganic′) under ideal conditions, the organopore correction coefficient can be obtained (Formula (1)).
Consider
(5)Фorganic=SPorganic·w(TOC)·ρrockρkerogen×100%,
where Фorganic is the organoporosity of shale in actual conditions (%) and SPorganic is the average organopore surface porosity in organic matters.
Using the focused ion beam (FIB) milling and scanning electron microscopy (SEM) imaging on the samples of Py1 well Longmaxi Shale, the surface porosity of organopores was recognized and analyzed by photo identification and forecast in two SEM images, which image two different organic matters in the same sample (Figures 6(a) and 6(c)). The depth of the sample is 2079.99 meters from Py1 well with the TOC value 1.45% and the rock density value 2.65 g/cm3. With the increasing of SEM imaging magnification, photo resolution is increased, and the smaller organopore size can be identified. The organopores smaller than the resolution of SEM were unable to be identified. The organopore surface porosity (SPorganic) is composed of the identified and unidentified organopores in SEM images. The contribution of identified organopores to SPorganic can be easily counted. For the unidentified organopores, the contribution of them to SPorganic was calculated through the following steps. (1) The distribution frequency of identified organopores in different sizes was analyzed through SEM images. (2) According to the distribution frequency of identified organopores in different sizes, the distribution frequency and size of unidentified organopores were speculated. (3) Assuming that the unidentified organopores are approximate spheres, the contribution of unidentified organopores to SPorganic can be calculated using the area formula of a circle.
The organic surface porosity and distribution frequency of the Py1 well from Longmaxi Shale in southeast Chongqing.
Combined with the SPorganic of identified and unidentified organopores, the surface porosities of different diameter organopores were obtained (Figures 6(b) and 6(d)). The column diagrams of Figures 6(b) and 6(d) are, respectively, corresponding to SEM images of Figures 6(a) and 6(c). The organopore surface porosities of the two SEM images are 32.30% and 28.07% and their average value is 30.19%. The organoporosity value is 0.97% according to Formula (5), while the organoporosity value is 3.20% using the material balance principle (Formula (1)) at the same depth of Py1 well. Finally, the organopore correction coefficient was calculated using the organoporosity value 0.97% in actual geological conditions (Фorganic) divided by the organoporosity value 3.20% under ideal conditions (Фorganic′) according to Formula (1) and Formula (2). The organopore correction coefficient is 0.303.
5. Organoporosity of Shale
Using the transformation ratio of hydrocarbon generated from organic matter, the original hydrogen index, the original total organic carbon, and the organopore correction coefficient of the Py1 well Longmaxi Shale in southeast Chongqing (Figures 4, 5, and 6), the organoporosity was calculated. The organoporosity without correction range is 0.66–9.10%, and the average value is 4.14%, whereas the organoporosity with correction range is 0.20–2.76%, and the average value is 1.25%. The organic porosities of Py1 well Longmaxi Shale before and after correction were decreased gradually from bottom to top, and both have the same variation trend with depth (Figure 7). And this variation trend and the residual organic carbon are consistent. So the residual organic carbon can indicate the relative levels of organoporosity, while the samples are in the same shale reservoirs with similar buried depths.
The organic porosities of the Py1 well from Longmaxi Shale in southeast Chongqing.
6. Conclusions
(1) Using the material balance principle and chemical kinetics, an organoporosity evaluation modelling method for shale rocks was developed. The important parameters for the method are the transformation ratio of hydrocarbon generated from organic matter, the original hydrogen index, the original total organic carbon, and the organopore correction coefficient.
(2) The organoporosity of the Lower Silurian Longmaxi Shale in the Py1 well is from 0.20 to 2.76%, and the average value is 1.25%.
(3) The organoporosity variation trends and the residual organic carbon of Longmaxi Shale are consistent in section. The residual organic carbon is indicative of the relative levels of organoporosity, while the samples are in the same shale reservoirs with similar buried depths.
Conflict of Interests
There is no conflict of interests regarding the publication of this paper.
Acknowledgments
Financial support from the National Science Foundation of China (Grant nos. 41302101 and 41330313) and the 973 Prophase Special Program of China (Grant no. 2011CB211701) is acknowledged.
LoucksR. G.ReedR. M.RuppelS. C.JarvieD. M.Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the mississippian barnett shale2009791284886110.2110/jsr.2009.0922-s2.0-72849145583LiuS. G.MaW. X.LubaJ.HuangW. M.ZengX. L.ZhangC. J.Characteristics of the shale gas reservoir rocks in the lower silurian longmaxi formation, east sichuan basin, China2011278223922522-s2.0-84863006533ChalmersG. R.BustinR. M.PowerI. M.Characterization of gas shale pore systems by porosimetry, pycnometry, surfacearea, andfield emission scanning electron microscopy/transmission electron microscopy image analyses: examples from the Barnett, Woodford, Haynesville, Marcellus, and Doig units20129661099111910.1306/101711110522-s2.0-84862742477JarvieD. M.HillR. J.RubleT. E.PollastroR. M.Unconventional shale-gas systems: the Mississippian Barnett Shale of north-central Texas as one model for thermogenic shale-gas assessment200791447549910.1306/121906060682-s2.0-34248330129ChristopherJ. M.ScottG. L.Estimation of kerogen porosity in source rocks as a function of thermal transformation: Example from the Mowry Shale in the Powder River Basin of Wyoming201296187108CurtisM. E.AmbroseR. J.SondergeldC. H.RaiC. S.Structural characterization of gas shales on the micro- and nano-scalesProceedings of the Canadian Unconventional Resources and International Petroleum ConferenceOctober 2010Calgary, CanadaSPE Paper 1376932-s2.0-79952907583HanS.ZhangJ.YangC.LinL.ZhuL.ChenY.XueB.The characteristics of nanoscale pore and its gas storage capability in the Lower Cambrian shale of southeast Chongqing2013386103810432-s2.0-84880382152AmbroseR. J.HartmanR. C.Diaz-CamposM.AkkutluI. Y.SondergeldC. H.New pore-scale considerations for shale gas in place calculationsProceedings of the Society of Petroleum Engineers Unconventional Gas ConfertuobeienceFebruary 2010Pittsburgh, Pa, USASPE Paper1671832-s2.0-77954232797CurtisM. E.SondergeldC. H.AmbroseR. J.RaiC. S.Microstructural investigation of gas shales in two and three dimensions using nanometer-scale resolution imaging201296466567710.1306/081511101882-s2.0-84860733041MosherK.HeJ. J.LiuY. Y.RuppE.WilcoxJ.Molecular simulation of methane adsorption in micro- and mesoporous carbons with applications to coal and gas shale systems2013109-110364410.1016/j.coal.2013.01.0012-s2.0-84874720743HowardJ. J.Porosimetry measurement of shale fabric and its relationship to illite/smectite diagenesis199139435536110.1346/CCMN.1991.03904032-s2.0-0026292207SondergeldC. H.AmbroseR. J.RaiC. S.MoncrieffJ.Microstructural studies of gas shaleProceedings of the Society of Petroleum Engineers Unconventional Gas ConferenceFebruary 2010Pittsburgh, Pa, USASPE Paper 131771TianH.ZhangS. C.LiuS. B.ZhangH.Determination of organic-rich shale pore features by mercury injection and gas adsorption methods20123334194272-s2.0-84863534854MoncrieffJ.Microstructure of shalesProceedings of the Annual Society of Petrophysicists and Well Log Analysts SymposiumJune 2009The Woodlands, Tex, USAWangF. P.ReedR. M.Pore networks and fluid flow in gas shales2009124253New Orleans, LouisianaSociety of Petroleum Engineers Annual Technical Conference and ExhibitionPasseyQ. R.BohacsK. M.EschW. L.KlimentidisR.SinhaS.From oil-prone source rock to gas-producing shale reservoir: geologic and petrophysical characterization of unconventional shale-gas reservoirs2010131350Beijing, ChinaChinese Petroleum Society, Society of Petroleum Engineers International Oil and Gas Conference and Exhibition in ChinaSchieberJ.Common themes in the formation and preservation of intrinsic porosity in shales and mudstones: Illustrated with examples across the PhanerozoicProceedings of the Society of Petroleum Engineers Unconventional Gas ConferenceFebruary 2010Pittsburgh, Pa, USASPE Paper 132370ZouC.ZhuR.BaiB.YangZ.WuS.SuL.DongD.LiX.First discovery of nano-pore throat in oil and gas reservoir in China and its scientific value2011276185718642-s2.0-80051702494ChalmersG. R. L.RossD. J. K.BustinR. M.Geological controls on matrix permeability of Devonian Gas Shales in the Horn River and Liard basins, northeastern British Columbia, Canada201210312013110.1016/j.coal.2012.05.0062-s2.0-84868092779UtpalenduK.DouglasK. M.ArkadiuszD.TimothyB. F.ManikaP.Total porosity measurement in gas shales by the water immersion porosimetry (WIP) method20141171115112910.1016/j.fuel.2013.09.073LongH.SwennenR.FoubertA.DierickM.JacobsP.3D quantification of mineral components and porosity distribution in Westphalian C sandstone by microfocus X-ray computed tomography20092201-211612510.1016/j.sedgeo.2009.07.0032-s2.0-69349086129MilnerM.McLinR.PetrielloJ.Imaging texture and porosity in mudstones and shales: Comparison of secondary and ion-milled backscatter SEM methodsCanadian Unconventional Resources and International Petroleum Conference 2010October 2010can305630652-s2.0-79953095784SlattR. M.O'BrienN. R.Pore types in the Barnett and Woodford gas shales: contribution to understanding gas storage and migration pathways in fine-grained rocks201195122017203010.1306/033011101452-s2.0-84856507294WangF.XuX. S.Tectonic-lithofacies palaeogeograhy of the Silurian in Sichuan-Yunnan-Guizhou-Guangxi region200352180185LiY. F.FanT. L.GaoZ. Q.ZhangJ. C.WangX. M.ZengW. T.ZhangJ. P.Sequence stratigraphy of Silurian black shale and its distribution in the southeast area of Chongqing2012233299306ZhanJ. C.JianS. L.TangX.ZhangP. X.TangY.JingT. Y.The accumulation type and characteristic of resource about shale gas in China2009291216DongD. Z.ChengK. M.WangY. M.LiX. J.WangS. J.HuangJ. L.Forming conditions and characteristics of shale gas in the Lower Paleozoic of the Upper Yangtze region, China2010313288299XieC.ZhangJ.LiY.WangX.Characteristics and gas content of the Lower Cambrian dark shale in Well Yuke 1, Southeast Chongqing201334111152-s2.0-84877041197PuB. L.JiangY. L.WangY.BaoS. J.LiuX. J.Reservoir-forming conditions and favorable exploration zones of shale gas in Lower Silurian Longmaxi Formation of Sichuan Basin20103122252302-s2.0-77952003230LiuZ. F.XiaoX. M.FuJ. M.LiuD. H.ShenJ. G.Marine vitrinite reflectance as a maturity indicator of lower Palaeozoic hydrocarbon source rocks1999286580588OkiongboK. S.AplinA. C.LarterS. R.Changes in type II Kerogen density as a function of maturity: evidence from the Kimmeridge clay formation20051962495249910.1021/ef050194+2-s2.0-29044440892XiaoZ. Y.LuY. H.SangH.PanZ. Z.LiY. F.A typical Cambrian oil reservoir: origin of oil reservoir in Well TZ62, Traim Basin2005342155160LuS. F.1996Beijing, ChinaPetroleum Industry Press(Chinese)LuS. F.LiuX. Y.FuX. T.FengY. L.WangZ. P.XueH. T.Chemical kientics study on the generation mechanism of immature to low-mature oil and its initial application2001191130135MinW.Shuang-FangL.Hai-TaoX.Kinetic simulation of hydrocarbon generation from lacustrine type I kerogen from the Songliao Basin: model comparison and geological application20112891714172610.1016/j.marpetgeo.2011.07.0042-s2.0-80051801566LiS.ZhouY.XiaoK.WoY.WangX.LiuQ.Characteristics of Silurian destroyed oil reservoir in Houtan section of Xishui area in southeastern margin of Sichuan Basin20093068498552-s2.0-72449139515JinZ. J.YuanY. S.LiuQ. Y.WoY. J.Controls of late Jurassic-early cretaceous tectonic event on source rocks and seals in marine sequences, South China2012421217911801MaW. X.LiuS. G.HuangW. M.ZhangC. J.XuG. S.YuanH. F.Characteristics of Silurian Paleo-oil reservoirs and their significance for petroleum exploration on the southeast margin of Sichuan Basin2012333432441LuS. F.FuG.WangP. Y.Beijing, ChinaPetroleum Industry Press2002 (Chinese)