CHINESE JOURNAL OF GEOPHYSICS  2015, Vol. 58 Issue (6): 599-610   PDF    
THERMAL HISTORY RECONSTRUCTION BASED ON MULTIPLE PALEO-THERMAL RECORDS OF THE YAZIHE AREA, WESTERN SICHUAN DEPRESSION, SW CHINA
ZHU Chuan-Qing1,2, QIU Nan-Sheng1,2, JIANG Qiang2, HU Sheng-Biao3, ZHANG Shuo1,2     
1 State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China;
2 College of Geosciences, China University of Petroleum(Beijing), Beijing 102249, China;
3 State Kay Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
Abstract: The thermal history of boreholes in the western Sichuan depression was reconstructed using paleotemperature indicators, such as vitrinite reflectance (Ro), apatite fission track (AFT) and the illite crystallinity index (IC index). The maximum paleo-temperatures reconstructed based on different indicators were compared. The results showed that the western Sichuan depression has experienced a continuous process of denudation and cooling since the late Cretaceous (K3), resulting in a geothermal gradient reduction from~26℃·km-1 to~22℃·km-1 and the removal of 1.3 to 1.9 km of sediment. The surface area has undergone continuous uplift and denudation from~80 Ma to 40 Ma, followed by a thermal quiet period until 2.5 Ma. Since~2.5 Ma, the geothermal gradient has risen by~5℃·km-1. In addition, there was high consistency among the paleotemperatures that were reconstructed based on different indicators. Relative to Ro and AFT, the illite crystallinity index is still controversial; more research is necessary before it can be used as a maturity and geothermal indicator for sedimentary rock.
Key words: Vitrinite reflectance     Apatite fission track     Illite crystallinity index     Thermal history reconstruction     Western Sichuan depression    
1 INTRODUCTION

The thermal history of a sedimentary basin refers to the thermal evolution of the basin in geological history; this may include the evolution of heat flow, the geothermal gradient, and the stratum temperature. Thermal history is vital to geodynamics and petroleum geology (Qiu, 2005), for the temperature is a driving force behind many geodynamic processes (Morgan, 1984; Wang, 1996). For instance, it affects the hydrocarbon accumulation process, especially the maturation of hydrocarbon source rocks and the history of hydrocarbon generation (Tissot and Welte, 1984; Tissot and Pelect, 1987; Qiu et al., 2004). The thermal history of a basin is usually reconstructed using various paleo-temperature indicators, including Ro, AFT, (U-Th)/He in apatite and zircon, reflectance of bitumen and vitrinite-like macerals, the hydrogen index, and fluid inclusions (Hu et al., 1998; Qiu et al., 2004; Qin et al., 2009a; 2009b). Ro and AFT are the most widely used indicators, and the related modeling methods are relatively well established. Recently, other indicators such as illite crystallinity (IC), Raman spectroscopy, and acoustic rock emissions have also been successfully studied and used in the quantitative reconstruction of paleo-temperatures (Xu et al., 2005; Wang et al., 2007; Qin et al., 2009a; 2009b; Li et al., 2011; Liu et al., 2013).

The marine basins in China’s mainland were formed in the Paleozoic era or even earlier and have experienced several periods of tectonic superposition and alteration. Because of this long-term, complex structuralthermal evolution, these basins have a multi-stage and complex thermal history. The continental basins lying above these marine basins have an entirely different thermal evolution, or thermal history, and thus require different methods and effective geothermal indicators to reconstruct their thermal history. In a superimposed basin composed of a continental and ocean basin (such as the representative Sichuan Basin and Tarim Basin), the Paleozoic marine strata are normally important source rock layers. These old strata contain very small quantities of vitrinite, apatite, and other conventional indicative components because they have experienced a high degree of thermal evolution. As a result, it has long been difficult to identify their maturity and reconstruct their thermal histories (Qin et al., 2009b). Thermal history reconstruction for superimposed basins requires diversified and complementary thermal history reconstruction systems (Hu et al., 2008) and, more importantly, effective geothermal indicators. The equivalence between different indicators also needs to be investigated to establish a set of geothermal indicators that applies to the reconstruction of multi-stage thermal history.

In this paper, the thermal history of the western Sichuan depression in the Sichuan Basin was reconstructed using three geothermal indicators: Ro, AFT and the IC index. The maximum paleo-temperatures reconstructed based on different indicators were compared. This research will help discover the equivalence between different geothermal indicators and establish a set of geothermal indicators for thermal history reconstruction.

2 GEOLOGICAL SETTING

Located in the western part of Sichuan Basin, the western Sichuan depression (Fig. 1) is a foreland basin occupying an area of approximately 5×104 km2. It is bounded by the Longmen Mountain front buried fault to the west and the Longquan Mountain fault to the east (Chen et al., 2010). It evolved from a Paleozoic cratonic basin. Before the Indosinian orogeny in the late Triassic (T3), the western Sichuan region was at the western boundary of the upper Yangtze Cratonic basin, collectively forming a westward sloping continental platform composed of carbonate rock (Ma et al., 2009). As the Longmen Mountain and Daba Mountain orogenic belts developed during the Indosinian orogeny (Liu et al., 2009; Li et al., 2011), seawater levels regressed, and continental sedimentation initiated in this area. An ultra-thick Upper Triassic series composed primarily of lithic sandstone formed when voluminous clasts were deposited in the western Sichuan (Chen et al., 2011). After the Indosinian orogeny, the evolution of the foreland basin began. At present, there has been considerable research on the tectono-sedimentary evolution of the foreland basin (Liu et al., 1994; Liu et al., 1995; 2011; Guo et al., 1996; Tao, 1999; Liu et al., 2010; Li et al., 2011; Tian et al., 2013). Longmen Mountain has been thrust onto the western Sichuan depression several times since the Indosinian orogeny began, resulting in regional faults and the uplift and denudation of the basin margins. The Indian Plate then collided with the Eurasian Plate during the Himalayan orogeny, inducing more violent thrusting of Longmen Mountain. These structural events created the Longmen Mountain thrust belt and the foreland basin in western Sichuan (Wang et al., 1997; Zhao et al., 2011). The strata burial history and the rate of deposition and denudation of the borehole CY92 in the Yazihe area, western Sichuan depression shows the obvious and continuous uplift since~60 Ma (Fig. 2).

Fig. 1 Geological sketch map and the stratigraphic column of the western Sichuan depression (modified from Zhao et al., 2011; Zhang et al., 2007)

Fig. 2 Strata burial history (a) and the rate of deposition and denudation of the borehole CY92 in the Yazihe area, western Sichuan depression

The Yazihe area is a semi-anticline formed in the late Himalayan orogeny, constituting the lower footwall of the Guankou fault in the front of the middle Longmenshan nappe belt. It is located in the mid-south section of the southern Sichuan depression and is a northeast faulted anticline under the control of the Guankou and Pengxian faults (Zhu et al., 2001; Liu et al., 2002). The conditions of the source-reservoir-cap rock assemblage allow the petroleum systems in this area to be divided into two parts: the upper continental petroleum system in the upper Triassic to upper Jurassic series (T3-J3) and the lower marine petroleum system in the middle Triassic (T2) and the underlying strata (Zhang et al., 2007; Zhao et al., 2011). The boreholes in the Yazihe area, such as CY92 and CY95, were drilled in the upper continental petroleum system in T3-J3 (Fig. 1). In this upper system, the T3 coal-bearing strata, formed from coal, carbonaceous shale, gray to black mudstone, and silty mudstone, are rich in hydrocarbon source rocks. The lithostratigraphic units that are composed primarily of source rocks include the Ma’antang Formation, the Xiaotangzi Formation, and the third and fifth members of the Xujiahe Formation, where the dark mudstone layers total hundreds of meters in thickness. In addition, the reservoirs in the second and fourth members of the Xujiahe Formation also contain organic-rich layers of source rocks (Yang et al., 2003; Chen et al., 2010). The sediments of alluvial fans, river deltas, and lake deltas in the sandstone strata of the Xujiahe Formation and the Jurassic System contain the main reservoirs of the continental petroleum system. The cap rocks above the Upper Triassic reservoir rocks, constituted by the third and fifth members of the Xujiahe Formation, and the mudstone and shale of the Suining Formation in J3, show large thicknesses, high breakthrough pressures, and stability. The Cretaceous gypsum layers serve as the direct cap rocks above the Jurassic reservoir rocks (Wang et al., 2002).

3 GEOTHERMAL INDICATOR DATA AND THERMAL HISTORY RECONSTRUCTION 3.1 Thermal History Reconstruction Based on Ro

Vitrinite is a type of maceral composed mainly of polycyclic aromatic hydrocarbons. With increasing thermal maturity, the benzene rings in the polycyclic aromatic hydrocarbon molecules become more concentrated, and the vitrinite structure grows denser, resulting in higher vitrinite reflectance (Chen et al., 2007). Temperature and the duration of heating of the organic matter cause the evolution of Ro values; temperature is the dominant factor. The methods used in early paleo-temperature reconstruction research, i.e., charting (Hood and Guijahr, 1975; Cooper, 1977) and TTI fitting (Lopatin, 1971; Waples, 1980; Lerche et al., 1984), have been replaced by the new method of chemical kinetic modeling of kerogen pyrolysis, as a result of more research into the thermodynamic mechanisms of Ro evolution (Armagnel et al., 1989; Burnham et al., 1989, 1987; Braun and Burnham, 1987; David and Antia, 1986; Burnham and Sweeney, 1989; Sweeney and Burnham, 1990). Currently, the EASY Ro% model constructed by Sweeney and Burnham (1990) is widely used to reconstruct thermal history.

CY92 samples for Ro analysis were mostly collected from the Triassic system approximately 2000 to 5000 m deep. Their Ro values varied from 0.8% to 2.0%, indicating that the vitrinite reflectance increased with depth. Based on the EASY Ro% model, the maximum paleo-temperatures of CY92 were reconstructed using the paleotemperature gradient method (Duddy et al., 1988; Hu et al., 1998; Qiu et al., 2004). The results showed that the paleo-temperature gradient of this borehole was~26 ℃·km−1, the heat flow was~60 mW·m−2, and the thickness of removed sediments on the surface of the unconformity between the Upper Jurassic and the Cenozoic was approximately 1900 m. As Fig. 3 shows, the paleo-temperature gradient of this borehole was larger than its present temperature gradient (~22 ℃·km−1). This indicates that the area has undergone an overall, continuing structural-thermal evolution of cooling, uplifting, and denudation since the rock strata reached the maximum temperatures.

Fig. 3 Ro data (a) and thermal history reconstruction of CY92 using the paleo-temperature gradient method (b); the rock thermal properties data refer to Xu et al. (2011)
3.2 Thermal History Reconstruction Based on AFT

Research on the thermal histories of sedimentary basins has used low-temperature thermochronological methods to analyze the fission tracks and (U-Th)/He in minerals such as apatite and zircon. The annealing temperatures of fission tracks in apatite are between 60 and 120 ℃, which is basically consistent with the hydrocarbon generation temperature range of source rocks (Green et al., 1986; 1989). Moreover, AFT has connections with various geological processes occurring at a crustal depth of 3~5 km (Li et al., 2013) and is susceptible to annealing. All of these properties make AFT analysis an important technique for modeling the uplifting history of orogenic belts, the subsidence and erosion history of sedimentary basins, and the temperature history of hydrocarbon source rocks. Considerable research on the principles and modeling methods of low-temperature thermochronology has been conducted (Gleadow, 1981; Gleadow et al., 1983; Gleadow and Fitzgerald, 1987; Fitzgerald and Gleadow, 1990; Fitzgerald et al., 1993; Armstrong et al., 1997; Armstrong, 2005; Ketcham et al., 2007). The age and length (of a closed fission track) of AFT are the frequently used data in studies on paleo-temperature based on AFT analysis.

CY92 samples for AFT analysis were collected from the strata at 2000~4500 m in depth. In CY92-1 (at 2142 m), the central AFT age was 42.3±3.0 Ma, and the average track length was 11.76±0.26 µm. In CY92-3 (at 4474 m), the central AFT age was 4.3±0.6 Ma, and no track length data were obtained (Fig. 4). To model the thermal history of CY92-1, modeling of fission-track annealing in apatite (Ketcham et al., 2007) was performed with HeFty v1.7.4 software and random AFT inversion (Corrigan, 1991). The results revealed that the basin experienced subsidence between 220 and 140 Ma and a thermal quiet period between 145 and 80 Ma. Crustal uplift and erosion occurred in this area between 80 and 40 Ma. After another thermal quiet period from 40 to 2.5 Ma, the temperatures of the strata began to rise. By the late stage of heating, the temperature had increased by approximately 10 ℃. The contribution of the burial of sediment to the temperature increase can be ignored, as the Quaternary sediment was only 23.5 m thick. Therefore, the temperature increase reflected an increase of 5 ℃·km−1 in the geothermal gradient (and heat flow). According to the modeling results, CY92-1 reached a maximum paleo-temperature of approximately 115~120 ℃ when it was buried at a maximum depth.

Fig. 4 AFT data of the samples of CY92 (a) and the radar map of AFT ages of CY92-1 (b)

Fig. 5 AFT thermal history modeling result of CY92-1 (a) and the fitting of AFT length distribution (10000 paths have been modeled using Monte Carlo method, GOF (age)=0.98; GOF (length)=0.95)

The maximum paleo-thermal temperature profile of CY92 (with a corresponding geothermal gradient of~26 ℃·km−1) shows that the maximum depth was 4040 m, compared to its present depth of 2142 m. The thickness of the removed sediments was estimated to be~1900 m, which is generally consistent with the Ro analysis results.

3.3 Thermal History Reconstruction Based on the IC Index

Illite is one of the common clay minerals in sedimentary rocks, and authigenic illite crystals are an important product of the diagenesis and metamorphic evolution of sedimentary rocks. Because illite crystallinity is closely associated with degrees of diagenesis and metamorphism and increases with temperature, it can not only be used to research diagenesis and very low-grade metamorphism but also serve as a paleo-geothermometer to reconstruct thermal history (Wang et al., 2007; Zhu and Zhu, 2006; Bignall et al., 2001; Miller and Macdonald, 2004; Aldega et al., 2003; Di, 2003; Ji and Browne, 2000). The widely used Kbler index, which is negatively correlated with illite crystallinity, uses X-ray diffraction to denote illite crystallinity in international illite crystallinity studies. Crystallinity is determined by measuring the full width at half maximum for the 10 Å(001) diffraction reflection peak of illite (Kübler, 1964; 1967; 1968; Yang, 1993; Zhu, 1995; Wang et al., 2000; You et al., 2008).

To minimize the discrepancy in test results caused by different sample preparation procedures, equipment conditions, and test methods, the Kübler index data can be calibrated based on related international standards (Wang, 1998). This improves the comparability and reliability of the data. In accordance with related international standards, this paper’s data on the illite crystallinity of the samples from CY95 (Fig. 6) have been calibrated by the Petroleum Geology & Experiment Center of the Research Institute of Petroleum Exploration & Development (You et al., 2008). According to the linear fitting, the relationship between the Kübler index and depth can be described as a decreasing linear function: Y=−9168.76X + 9988.64.

Fig. 6 KI data and the estimation of paleotemperature of CY92

Sedimentary rocks will undergo anchimetamorphism as a transition to epimetamorphism when the diagenesis ends. During this process, the rock tempe rature gradually rises, the illite crystallinity increases, and thus the Kübler index decreases. Some researchers have reported their findings regarding the relationship between crystallinity and the generation temperature of illite and the division of diagenesis to the very low-grade metamorphism process based on illite crystallinity. During the telodiagenesis stage (corresponding to the high-grade diagenetic zone), the KI for illite crystallinity varies from 1.0 to 0.42 (°∆2θ) and the temperature cap is 200 ℃ (when KI=0.42 (°∆2θ), T=~200 ℃). During burial metamorphism, the KI is between 0.25 and 0.42 (°∆2θ), and the temperature is between 200 and 300 ℃ (when KI=0.25(°∆2θ), T=~300 ℃). During epimetamorphism, which corresponds to the greenschist facies, the KI is smaller than 0.25 (°∆2θ), and the temperature is higher than 300 ℃ (Kübler, 1967; Arkai, 2002; Bi and Mo, 2004). While the application of illite crystallinity in low grade metamorphism studies is mature, the exact relationship between KI and temperature during diagenesis (KI > 0.42 (°∆2θ)) has not been determined. In light of the division of diagenesis in metamorphism processes (Arkai et al., 2002), the KI values of illite in CY95 and the maximum paleo-temperatures attained at different depths were estimated using the data in previous reports (KI=0.57 (°∆2θ), T=~141℃ (125~163 ℃) (Qin et al., 2009a); KI=0.56 (°∆2θ), T=156 ℃ (136~177 ℃) (Hu et al., 2012); KI=0.60~0.62 (°∆2θ), T=150 ℃ (Zhu and Zhu, 2006)). In CY95, at a depth of~4080 m, the KI was 0.60 (°∆2θ), and the maximum paleo-temperature was~150 ℃; at a depth of~6000 m, the KI was 0.42 (°∆2θ), and the maximum paleo-temperature was~200 ℃. The paleo-temperatures can be used to determine the maximum paleo-temperature profile of CY95 because the J3-T3 strata simultaneously reached the maximum paleo-temperatures before being uplifted and denudated.

4 RESULTS AND DISCUSSION

Figure 7 displays both the maximum paleo-temperature profile of CY92 that was reconstructed based on Ro and AFT, and the maximum paleo-temperature profile of CY95 that was reconstructed based on illite crystallinity. The maximum paleo-temperatures of CY92 reconstructed based on Ro and AFT were basically consistent, which demonstrates the strong comparability of the two widely used indicators. The maximum paleo-temperature profiles of CY92 and CY95 were nearly parallel, indicating similar geothermal gradients (~26 ℃·km−1) and paleo-heat flows (~60 mW·m−2) during their evolution. Moreover, this similarity corroborated the fact that the two adjacent boreholes were located in the same region and have similar geological settings due to similar structural and sedimentary evolutions (Liu et al., 2002).

Fig. 7 Comparison of the paleo-temperature reconstruction between CY91 and CY95

Because no Ro data of this borehole could be obtained, rock pyrolysis data1) of cores in CY95 were used. As Type Ⅲ kerogen is the primary organic matter in the source rocks in the Upper Triassic and the Jurassic of western Sichuan, the Ro data of the borehole were estimated by referring to the chart of the Tmax-Ro relationship of Type Ⅲ source rocks. When the depths varied from 2200 to 3000 m, the source rock temperatures, Tmax, were between~480 and 500 ℃, and the estimation of Ro was in the about 1.2%~1.4% range; when the depths varied from 3300 to 4000 m, Tmax varied from~480 to 500 ℃, and the estimation of Ro was in the about 1.6%~1.8% range. The estimations of Ro at different depths were all slightly smaller than the corresponding Ro values of CY92 (Fig. 3), but the temperature values were consistent with the maximum paleo-temperatures of the borehole reconstructed based on illlite crystallinity.

1)The geological report on the completed CY95 in Yazihe area, Pengxian, Sichuan Province. The No.11 Prospecting and Exploration Team of Southwest Bureau of Petroleum Geology, Ministry of Geology and Mineral Resources, 1985.

As shown in Fig. 7, the maximum paleo-temperature profiles of the two boreholes intersect with the line of depth=0 (namely, the earth’s surface) when T1=~72 ℃ and T2=~54 ℃, respectively. Given the geothermal gradient, ~26 ℃·km−1, a temperature difference of approximately 18° could indicate a difference of~690 m in the thickness of removed sediments. If this is combined with the result of AFT modeling, one can assume that the area being studied has undergone continuous cooling, uplift, and denudation since the Late Cretaceous, which resulted in a geothermal gradient reduction from about 26 ℃·km−1 to~22 ℃·km−1 and the removal of~1.3 to 1.9 km of sediments.

Estimating paleo-temperatures based on KI can yield relatively accurate results when applied to the burial metamorphism or epimetamorphism stage. However, though the KI-Ro relationship has been studied (Guthrie et al., 1986; Qin et al., 2009a; 2009b), the paleo-temperatures during diagenesis estimated with this method were unsatisfactory; the related formulas for temperature estimation using KI have been reported in related research:

(Browne and Harvey, 1992), (1)
(Ji and Browne, 2000). (2)

For example, the geothermal gradient calculated according to Formula (1) showed substantial differences from the data calculated using another geothermal indicator (Wang, 2004). In the Jiyang Depression, a Cenozoic graben basin, the present temperatures of the strata should be the maximum strata temperatures, and the paleo-temperature gradient of the basin should represent the present geothermal gradient. However, the result calculated using Formula (2), ~40 ℃·km−1 (Zhou, 2006; Jiang et al., 2008), was higher than the present geothermal gradient of~35 ℃·km−1 (Gong et al., 2003). According to the KI-Ro relationship established for the eastern Sichuan area (Qin et al., 2009a; 2009b), only a small fraction of KI data exceeded 0.42 (°∆2θ). As a result, KI=0.57(°∆2θ) deviates from the KI-Ro curve. Similar to the result of the paleo-temperature reconstruction based on Ro, the paleo-temperatures that were reconstructed based on illite crystallinity indicated that the KI data were associated with both the values and durations of maximum paleo-temperature (Qin et al., 2009b). However, differences were observed in the patterns of KI variations with temperature and time, which indicates that a particular formula for estimating temperature or an equivalent relationship cannot be applied to an area featuring a wider range of illite crystallinity. In this study, for example, when KI=0.25 was estimated according to the KI-depth relationship, the maximum paleo-temperature was not consistent with the maximum paleo-temperature profile of CY95 (Fig. 7). Research on the evolution of illite crystallinity during diagenesis is sparse compared to related research on low-grade metamorphism and the epimetamorphism process. Comparison is also difficult because research on the evolution of illite crystallinity during diagenesis usually has higher requirements for illite and test precision and often requires comparison of data from different laboratories. Further research is needed before studies can use illite crystallinity as a maturity indicator for organic matter during diagenesis.

5 CONCLUSIONS

(1) There was high consistency between the paleo-temperatures of the western Sichuan depression that were reconstructed based on Ro, AFT, and KI. Thus, these three indicators can be used as paleo-temperature indicators for the thermal history reconstruction of the area.

(2) The study area has undergone a continual process of denudation and cooling since the Late Cretaceous, resulting in a geothermal gradient reduction from~26 ℃·km−1 to~22 ℃·km−1 and the removal of~1.3 to 1.9 km of sediment. Sedimentation occurred slowly or was suspended during a thermal quiet period from 145 to 80 Ma. From~80 Ma to 40 Ma, the surface of the area experienced continuous uplift and denudation, followed by another thermal quiet period until 2.5 Ma. Since~2.5 Ma, the geothermal gradient has risen by~5 ℃·km−1 due to heating.

(3) Unlike Ro and AFT, the illite crystallinity index is still controversial due to the lack of an effective thermal evolution model, and more research is therefore necessary before it can be used as a maturity and geothermal indicator for the thermal history reconstruction of sedimentary rocks.

ACKNOWLEDGMENTS

This study was supported by the National Science Foundation of China (41102152), the PetroChina Innovation Foundation (2013D-5006-0102), the National Key Basic Research Development Plan of China (“973” Projects) (2012CB214703) and the Science Foundation of China University of Petroleum, Beijing (YJRC-2013-02).

References
[] Aldega L, Cello G, Corrado S, et al. 2003. Tectono-sedimentary evolution of the Southern Apennines (Italy):thermal constraints and modelling. Atti Ticinensi di Scienze Della Terra , 9 : 135-140.
[] Antia D D J. 1986. Kinetic method for modeling vitrinite reflectance. Geology , 14 (7) : 606-608. DOI:10.1130/0091-7613(1986)14<606:KMFMVR>2.0.CO;2
[] Arkai P, Sassi F P, Desmons J. 2002. Towards a unified nomenclature in metamorphic petrology:4. Very low-to low-grade metamorphic rocks. A proposal on behalf of the IUGS subcommission on the systematics of metamorphic rocks. Web version of 31 , 10 : 1-12.
[] Armagnel C, Bucci J, St Kendall G C, et al. 1989. Estimating the thickness of sediment removed at an unconformity using vitrinite reflectance data.//Naeser N D, McCulloh T H. Thermal History of Sedimentary Basins. New York:Springer-Verlag, Inc., 217-238.
[] Armstrong P A, Kamp P J J, Allis R G, et al. 1997. Thermal effects of intrusion below the Taranaki Basin (New Zealand):evidence from combined apatite fission track age and vitrinite reflectance data. Basin Research , 9 (2) : 151-169. DOI:10.1046/j.1365-2117.1997.00039.x
[] Armstrong P A. 2005. Thermochronometers in Sedimentary Basins. Reviews in Mineralogy and Geochemistry , 58 (1) : 499-525. DOI:10.2138/rmg.2005.58.19
[] Bi X M, Mo X X. 2004. Transition from diagenesis to low-grade metamorphism and related minerals and energy resources. Earth Science Frontiers (in Chinese) , 11 (1) : 287-294.
[] Bignall G, Tsuchiya N, Browne P R L. 2001. Use of illite crystallinity as a temperature indicator in the Orakei Korako geothermal system, New Zealand. Transactions-Geothermal Resources Council , 25 : 339-344.
[] Braun R L, Burnham A K. 1987. Analysis of chemical reaction kinetics using a distribution of activation energies and simpler models. Energy & Fuels , 1 (2) : 153-161.
[] Burnham A K, Braun R L, Gregg H R, et al. 1987. Comparison of methods for measuring kerogen pyrolysis rates and fitting kinetic parameters. Energy & Fuels , 1 (6) : 452-458.
[] Burnham A K, Oh M S, Crawford R W, et al. 1989. Pyrolysis of Argonne premium coals:activation energy distribution and related chemistry. Energy & Fuels , 3 (1) : 42-55.
[] Burnham A K, Sweeney J J. 1989. A chemical kinetic model of vitrinite maturation and reflectance. Geochimica et Cosmochimica Acta , 53 (10) : 2649-2657. DOI:10.1016/0016-7037(89)90136-1
[] Chen D X, Huang X H, Li L T, et al. 2010. Characteristics and history of hydrocarbon expulsion of the Upper Tertiary source rocks in the western Sichuan Depression. Natural Gas Industry (in Chinese) , 30 (5) : 41-45.
[] Chen Y. 2011. The formation of western Sichuan foreland basin and its significance in oil-gas exploration during late Triassic[Ph. D. thesis]. Chengdu:Chengdu University of Technology.
[] Chen Y C, Shen Z M, Luo X P. 2007. Oil and Gas Organic Geochemistry (in Chinese)[M]. Beijing: Sciences Press .
[] Cooper B S. 1977. Estimation of the maximum temperature attained in sedimentary rocks.//Developments in Petroleum Geology. London:Applied Science Publications, 127-146.
[] Corrigan J. 1991. Inversion of apatite fission track data for thermal history information. Journal of Geophysical Research , 96 (B6) : 10347-10360. DOI:10.1029/91JB00514
[] Di L P. 2003. Use of clay mineralogy in reconstructing geological processes; thermal constraints from clay minerals. AttiTicinensi di Scienze della Terra , 9 : 55-67.
[] Duddy I R, Green P F, Laslett G M. 1988. Thermal annealing of fission tracks in apatite 3. Variable temperature behaviour. Chemical Geology , 73 (1) : 25-38.
[] Fitzgerald P G, Gleadow A J W. 1990. New approaches in fission track geochronology as a tectonic tool:Examples from the transantarctic mountains. International Journal of Radiation Applications and Instrumentation. Part D. Nuclear Tracks and Radiation Measurements , 17 (3) : 351-357. DOI:10.1016/1359-0189(90)90057-5
[] Fitzgerald P G, Stump E, Redfield T F. 1993. Late Cenozoic uplift of Denali and its relation to relative plate motion and fault morphology. Science , 259 (5094) : 497-499. DOI:10.1126/science.259.5094.497
[] Gleadow A J W, Duddy I R, Lovering J F. 1983. Fission track analysis:a new tool for the evaluation of thermal histories and hydrocarbon potential. Australian Petroleum Exploration Association Journal , 23 : 93-102.
[] Gleadow A J W, Fitzgerald P G. 1987. Uplift history and structure of the Transantarctic Mountains:new evidence from fission track dating of basement apatites in the Dry Valleys area, southern Victoria Land. Earth and Planetary Science Letters , 82 (1-2) : 1-14. DOI:10.1016/0012-821X(87)90102-6
[] Gleadow A J W. 1981. Fission-track dating methods:What are the real alternatives. Nuclear Tracks , 5 (1-2) : 3-14. DOI:10.1016/0191-278X(81)90021-4
[] Gong Y L, Wang L S, Liu S W, et al. 2003. Distribution characteristics of geotemperature field in Jiyang depression, Shandong, north China. Chinese J. Geophys. (in Chinese) , 46 (5) : 652-658.
[] Green P F, Duddy I R, Gleadow A J W, et al. 1986. Thermal annealing of fission tracks in apatite:1. A qualitative description. Chemical Geology , 59 : 237-253.
[] Green P F, Duddy I R, Laslett G M, et al. 1989. Thermal annealing of fission tracks in apatite 4. Quantitative modelling techniques and extension to geological timescales. Chemical Geology , 79 (2) : 155-182.
[] Guo Z W, Deng K L, Han Y H. 1996. Formation and Evolution of the Sichuan Basin (in Chinese)[M]. Beijing: Geological Publishing House .
[] Guthrie G M, Houseknecht D W, Johns W D. 1986. Relationships among vitrinite reflectance, illite crystallinity, and organic geochemistry in Carboniferous strata, Ouachita mountains, Oklahoma and Arkansas. AAPG Bulletin , 70 (1) : 26-33.
[] Hood D, Gutjahr C C M, Heacock R L. 1975. Organic metamorphism and the generation of petroleum. AAPG Bulletin , 59 (6) : 986-996.
[] Hu D Q, Han C Y, Ma R, et al. 2012. The very low grade metamorphism in the Upper Paleozoic in Xinlingol area of Inner Mongolia, NE China:Evidence from studies of illite and vitrinite reflectance. Acta Petrologica Sinica (in Chinese) , 28 (9) : 3042-3050.
[] Hu S B, Zhang R Y, Zhou L C. 1998. Reconstruction of geothermal history in hydrocarbon basins. Petroleum Explorationist(in Chinese) , 3 (4) : 52-54.
[] Hu S B, He L J, Zhu C Q, et al. 2008. Method system of thermal reconstruction for marine basins. Oil and Gas Geology(in Chinese) , 29 (5) : 607-613.
[] Ji J F, Browne P R L. 2000. Relationship between illite crystallinity and temperature in active geothermal systems of New Zealand. Clays and Clay Minerals , 48 (1) : 139-144. DOI:10.1346/CCMN
[] Jiang H C, Xiao Y J, Zhou L. 2008. Analysis of Cenozoic subsurface temperatures of the Jiyang depression. Geology in China (in Chinese) , 35 (2) : 273-278.
[] Ketcham R A, Carter A, Donelick R A, et al. 2007. Improved modeling of fission-track annealing in apatite. Am. Mineral. , 92 (5-6) : 799-810. DOI:10.2138/am.2007.2281
[] Kübler B. 1964. Les argiles, indicateurs de métamorphisme. Rev. Inst. Franc. Petrole Ann. Combüst. , 19 : 1093-1112.
[] Kübler B. 1967. La cristallinite de I'illite et les zones tout a fait superieures du métamorphisme. Etages tectoniques. Colloque de Neuchatel, 105-121.
[] Kübler B. 1968. Évaluation qunaittative du métamorphismepar la cristallinitE de I'illité. Bulletin du Centre de Recherehes de pua-SNAP , 2 : 385-397.
[] Lerche I, Yarzab R F, St C, Kendall C G. 1984. Determination of paleoheat flux from vitrinite reflectance data. AAPG Bulletin , 68 (11) : 1704-1717.
[] Li J W, Qiu N S, Mei Q H, et al. 2011. Study on measuring the highest rock paleotemperature with thermo-acoustic emission. Chinese J. Geophys. (in Chinese) , 54 (11) : 2898-2905.
[] Li T Y, Zhou Y, Fang S, et al. 2013. A new method for apatite fission track datingthe Laser-ICPMS method. Oil & Gas Geology (in Chinese) , 34 (4) : 550-557.
[] Li Z W, Liu S G, Chen H D, et al. 2011. Structural superimposition and conjunction and its effects on hydrocarbon accumulation in the Western Sichuan Depression. Petroleum Exploration and Development (in Chinese) , 38 (5) : 538-551. DOI:10.1016/S1876-3804(11)60054-2
[] Liu D C, Li S S, Zhou Z J, et al. 2002. Oil and gas exploration target selection and evaluation of the Xujiahe formation in the Yazi River structure. Natural Gas Exploration and Development (in Chinese) , 25 (1) : 38-45.
[] Liu D H, Xiao X M, Tian H, et al. 2013. Sample maturation calculated using Raman spectroscopic parameters for solid organics:Methodology and geological applications. Chinese Science Bulletin , 58 (11) : 1285-1298. DOI:10.1007/s11434-012-5535-y
[] Liu H P, Li J M, Li X Q, et al. 2006. Evolution of cratonic basins and carbonate-evaporite sedimentary sequence hydrocarbon systems in China. Geoscience (in Chinese) , 20 (1) : 1-18.
[] Liu J H, Yang S C, Zhang S Q, et al. 2010. A study on the characters of the accommodation change and hydrocarbon in Upper Xujiahe formation in the Western Sichuan Basin. Geological Review (in Chinese) , 56 (2) : 275-282.
[] Liu S G, Yang R J, Wu X C, et al. 2009. The Late Triassic transition from marine carbonate rock to clastics in the western Sichuan Basin. Oil & Gas Geology (in Chinese) , 30 (5) : 556-565.
[] Liu S G, Deng B, Li Z W, et al. 2011. The texture of sedimentary basin-orogenic belt system and its influence on oil-gas distribution:A case study from Sichuan basin. Acta Petrologica Sinica (in Chinese) , 27 (3) : 621-635.
[] Lopatin N V. 1971. Temperature and geological time as factors of carbonification. Akademy Nauk SSSR Izvestias Serie Geologia , 3 : 95-106.
[] Ma C, Browne P R L, Harvey C C. 1992. Crystallinity of subsurface clay minerals in the Te Mihi Sector of the Wairakei geothermal system, New Zealand.//Proc. 14th New Zealand Geothermal Workshop, 267-272.
[] Ma Y S, Chen D H, Wang G L. 2009. Sequence Stratigraphy and Paleogeography in Southern China (in Chinese)[M]. Beijing: Science Press .
[] Miller S, Macdonald D I M. 2004. Metamorphic and thermal history of a fore-arc basin:the Fossil Bluff Group, Alexander Island, Antarctica. Journal of Petrology , 45 (7) : 1453-1465. DOI:10.1093/petrology/egh025
[] Morgan P. 1984. The thermal structure and thermal evolution of the continental lithosphere. Phys. Chem. Earth , 15 : 107-193. DOI:10.1016/0079-1946(84)90006-5
[] Qin J Z, Li Z M, Teng G E. 2009a. A study on paleo-geothermometer of high mature marine sequences in South China. Oil & Gas Geology (in Chinese) , 30 (5) : 608-618.
[] Qin J Z, Teng G E, Yang Q, et al. 2009b. Research on maturity indicators of high-maturity marine strata in the eastern Sichuan Basin. Acta Petrolei Sinica (in Chinese) , 30 (2) : 208-213.
[] Qiu N S, Hu S B, He L J. 2004. Theory and Application of the Sedimentary Basins Thermal Regime Research (in Chinese)[M]. Beijing: Petroleum Industry Press .
[] Qiu N S, Li H L, Jin Z J. 2005. Study of the thermal history reconstruction for Lower Paleozoic carbonate succession. Earth Science Frontiers (in Chinese) , 12 (4) : 561-567.
[] Sweeney J J, Burnham A K. 1990. Evaluation of a simple model of vitrinite reflectance based on chemical kinetics. AAPG Bulletin , 74 (10) : 1559-1570.
[] Tao X F. 1999. Evolution of nappe tectonic and foreland basin in the southern section of Longmen mountains. Journal of Chengdu University of Technology (in Chinese) , 26 (1) : 73-77.
[] Tian Y T, Kohn B P, Gleadow A J W, et al. 2013. Constructing the Longmen Shan eastern Tibetan Plateau margin:Insights from low-temperature thermochronology. Tectonics , 32 (3) : 576-592. DOI:10.1002/tect.v32.3
[] Tissot B P, Welte D H. 1984. Petroleum Formation and Occurrence. New York:Springer Verlag.
[] Tissot B P, Pelet R, Ungerer P. 1987. Thermal history of sedimentary basins, maturation indices, and kinetics of oil and gas generation. AAPG Bulletin , 71 (12) : 1445-1466.
[] Wang H J, Zhu M X, Xu Q S, et al. 2000. Relationship between slit system and the Kbler index of illite crystal unity and discussion on relevant problems. Geological Review (in Chinese) , 46 (6) : 588-593.
[] Waples D W. 1980. Time and temperature in petroleum formation:Application of Lopatin's method to petroleum exploration. AAPG Bulletin , 64 (6) : 916-926.
[] Wang H J, Tao X F, Rahn M. 2007. Some aspects of illite crystallinity and its applications in low temperature metamorphism. Earth Science Frontiers (in Chinese) , 14 (1) : 151-156.
[] Wang H J. 1998. On the error calculation of the Kbler index of illite crystallinity. Geological Review (in Chinese) , 44 (3) : 328-335.
[] Wang J. 1996. Geothermics in China[M]. Beijing: Seismological Press .
[] Wang Q. 2007. Clay minerals research and application:as geothermometers[Master's thesis] (in Chinese). Qingdao:China University of Petroleum.
[] Wang T Z. 1997. Significance of Longmenshan thrust in evolution and oil-gas exploration of western Sichuan basin. Geoscience (in Chinese) , 11 (4) : 496-500.
[] Wang Z C, Zhao W Z, Zhang L. 2002. Structural Sequence and Natural Gas Exploration of the Sichuan Basin (in Chinese)[M]. Beijing: Geological Publishing House .
[] Xu C H, Zhu G, Liu G S, et al. 2005. Application of Crystallinity of Illite to Recover Denudation Quantity:An example of cretaceous denudation quantity recovery of well Ancan 1 in Hefei basin, Anhui province. Geological Science and Technology Information (in Chinese) , 24 (1) : 41-44.
[] Xu M, Zhu C Q, Tian Y T, et al. 2011. Borehole temperature logging and characteristics of subsurface temperature in the Sichuan Basin. Chinese J. Geophys. (in Chinese) , 54 (4) : 1052-1060.
[] Yang K M, Zhu T, He L. 2003. Structural characteristics and exploration potential of the thrust block in Longmenshan. Petroleum Geology & Experiment (in Chinese) , 25 (6) : 685-700.
[] Yang X Z. 1993. Illite crystallinity and its geological significances. Acta Sedimentologica Sinaca (in Chinese) , 11 (4) : 92-98.
[] You J C, Bi X M, Hou C M. 2008. Calibration of illite crystallinity measurements by using international standards and its significance. Geosciences (in Chinese) , 22 (1) : 53-59.
[] Zhang Y G, Ma Z J, Yang K M, et al. 2007. The forecast of natural oil & gas potential in marine strata, western Sichuan basin, southwest China. Acta Geologica Sinica (in Chinese) , 81 (8) : 1041-1047.
[] Zhao J C, Liu S G, Sun W, et al. 2011. Analysis on petroleum preservation condition in the coupling area between Longmen Mountain and Sichuan Basin. Lithologic Reservoirs (in Chinese) , 23 (1) : 79-85.
[] Zhou J G. 2006. Analysis of Early-Cenozoic Prototype Basin and its structural thermal evolution characteristic study of Jiyang-Changwei depression (in Chinese)[M]. Beijing: Graduate University of China Academy of Sciences .
[] Zhu G. 1995. Grading the extreme-low metamorphic clastic sedimentary rocks by the crystallinity of the illite. Petroleum Exploration and Development (in Chinese) , 22 (1) : 33-34.
[] Zhu L, Zhu M. 2006. Thermal evolution stages for mesozoic strata in the Hefei basin:contstraints from illite crystallinity. Geology of Anhui (in Chinese) , 16 (3) : 169-172.
[] Zhu T, Liang E Y. 2001. Discussion on oil and gas exploration of Xujiahe group in Yazi River structure of the middle Longmenshan in western Sichuan. Journal of Chengdu University of Technology (in Chinese) , 28 (1) : 59-63.