Diagenetic evolution of the Oligocene Huagang Formation in Xihu sag, the East China Sea Shelf Basin

The deep and ultra deep clastic reservoir is characterized by strong reservoir heterogeneity and complicated reservoir-forming characteristics for its high degree of diagenetic stage and the complexity of diagenesis. In order to better study the diagenetic evolution of deep and ultra deep reservoir in the burial process, a fine subsection scheme of 36 diagenetic micro-stage (DS) in diageneitic process was proposed based on paleotemperature (T), vitrinite reflectance (Ro%) and proportion of smectite in illite/smectite interstratified minerals (I/S-S%). Taking the Oligocene Huagang Formation in the Xihu sag of the East China Sea Basin as an example, the diagentic stage IIA1–IIA2–IIB was identified mainly by means of formation temperature data (T), homogenization temperature of fluid inclusion (Th), vitrinite analysis and clay mineral X-ray diffraction method. On this basis, diagenetic evolution using the fine subsection scheme in geological time were conducted. This fine division of diagenetic stage could bring accurate insight into porosity evolution history, hydrocarbon charging periods, diagenetic and reservoir-forming characteristics of low-permeability and tight sandstone.

With the development of petroleum exploration and exploitation, deep (3500-4500 m) to ultra-deep (> 4500 m) tight sandstone has been paid more attention because of its abundant oil and gas resources 1,2 . The deep and ultra deep clastic reservoir is characterized by strong reservoir heterogeneity and complicated reservoir-forming characteristics for its high degree of diagenetic stage and the complexity of diagenesis [3][4][5] . Considerable efforts have been made to understand the spatial and temporal distribution of diagenesis in clastic rocks, maily due to the realization that these alterations exert profound controls on reservoir quality [6][7][8][9] . It is now evident that fluid properties, diagenesis of clastic rock and hydrocarbon generating capacity are closely related with diagenetic stages, and thus research on diagenetic evolution has important significance for exploration and exploitation of low-permeability and tight sandstone 3,8,[10][11][12][13] .
The recognization that diagenetic evolution has the characteristics of stages [14][15][16][17] has already begun in the early 1970s, and Schmidt and Morad (1979) 18 first put forward the division basis and naming method of diagenetic stages. Domestic scholars have also proposed the division standard of diagenetic stages of clastic rocks on basis of previous research 11,19,20 (Table 1). In recent years, a subdividing the intermediate diagenetic stage A into four diagenetic microfacies was proposed using vitrinite reflectance 17 . Throughout those researches, the division theory of diagenetic stage from qualitative and quantitative perspectives is mainly based on the diagenetic environment parameters such as temperature, time, pressure, fluid, while the contact relationship of clastic particles, pore structure, authigenic mineral assemblage, vitrinite reflectance (Ro%), peak temperature of organic matter solution (Tmax), clay mineral conversion rate (I/S-S%), steranes horane isomerization index (SI) and quartz secondary enlargement index are integrated.
The Oligocene Huagang Formation (E 3 h) with buried depth from 2000 to 6000 m is a vital gas producing formation in the Xihu sag, East China Sea Basin [22][23][24] . In the past five decades, lots of oil and gas fields and structures have been discovered in the exploration of Xihu sag of the East China Sea Basin. For the sandstone reservoir of the E 3 h Formation, previous scholars mainly focused on the following aspects: (1) lithologic features, diagenetic environment evolution and controlling factors of the Oligocene Huagang formation tight sandstone on the central inverted structural belt 25-28 ; (2) the tectonic activities on diagenetic and reservoir-forming process

Geological setting
The Xihu sag with about 440 km long from north to south, 110 km wide from east to west and an area of 51,800 km 2 is situated in the East China Sea Shelf Basin and generally oriented in NNE direction 23,33,34 (Fig. 1). In the light of evolutionary characteristics of deposition and structures in the course of its formation and development, the Xihu sag can be divided into five structural zones as follows: Western slope belt (WSB), Western sag (WS), Central uplift belt (CUB), Eastern sag (ES) and Eastern sharp slope (ESS) 23,33,34 (Fig. 1C). The formation and evolution of the Xihusag have experienced five major tectonic movements, called Keelung movement (TG), Pinghu movement (T40), Yuquan movement (T30), Huagang movement (T20), Longjing movement (T12) and Okinawa Trough movement respectively, which formed six regional unconformities in the sag. Among them, Yuquan movement (T30) and Huagang movement (T20) divided the structural layer of Xihu sag into three parts from bottom to top: rift structural layer Synrift stage (Tg-T30), sag structural layer (T30-T12) and regional subsidence structural layer (T12-seafloor), and Yuquan movement at the end of Eocene and Longjing movement at the end of Miocene are two important periods of structural change in Xihu sag. The regional structural evolution of Xihu sag can be roughly divided into three stages: synrift stage, depression stage, regional subsidence stage 23,24 . The Cretaceous Shimentan formation, Eocene Bajiaoting formation, Baoshi formation, Pinghu Formation, Oligocene Huagang Formation, Miocene Longjing formation, Yuquan formation and Liulang formation, Pliocene Santan formation and quaternary Donghai group are developed in Xihu sag from the bottom to the top, while the Paleocene stratum is not developed in the sag 24,35 (Fig. 2). The main target reservoir in the Xihu sag is the Eocene Pinghu formation and the Oligocene Huagang formation. The Huagang formation can be subdivided into the upper and lower members (HS and HX) 22,29 . The main stratigraphic features are as follows.
(1) The lower member of the Huagang formation was deposited in the early stage of depression, but it was closed by retrogression in the late stage of the Pinghu movement and transitioned to lacustrine deposits. The thickness revealed by drilling was 80-829 m. The lithology mainly consists of a thin layer of glutenite, gravel sandstone, coarse sandstone, medium sandstone, fine sandstone, siltstone, argillaceous siltstone, mudstone, and the local entrapment of a thin coal seam, reflecting a sedimentary period with adequate sediment sources. (2) The upper member of the Huagang formation was deposited in the early stage of the depression. A lacustrine transgression and retrogression occurred again after the lacustrine regression. The lacustrine transgression was large in scale, and the thickness revealed by drilling was 105-1100 m. The lithology is coarse and mainly consists of a thin layer of glutenite, gravelly sandstone, coarse sandstone, medium sandstone, fine sandstone, siltstone, argillaceous siltstone, mudstone, shale, and oil shale, reflecting a sedimentary period with adequate sediment sources. www.nature.com/scientificreports/ tion detector. The carbon dioxide discharged from pyrolysis and generated by the heated oxidation of pyrolytic residual organic matter were detected by thermal conductivity detector. Pyrolysis peak temperature (Tmax) is the pyrolysis temperature corresponding to the highest point of S2.
The division of diagenetic stage. Based on the burial history and thermal history, with the help of Sweeney's chemical kinetic method 37 , Ro% at different historical times can be calculated. The transformation process from smectite to illite is mainly affected by temperature, time, pressure, sedimentary environment and source characteristics 21,38,39 , which can be expressed by Elliot's model formula 40 . In this study, the influence of four geological factors, temperature, pressure, fluid and time on diagenesis was considered comprehensively, while three parameters including palegeotemperature T (℃), vitrinite reflectance (Ro%) and smectite content in mixed layer (I/S-S%) were selected as the main parameters of comprehensive diagenetic evolution simulation to simulate the distribution of time domain and space domain and to divide the diagenetic stage in geological time.

Results
Formation temperature. From the shallow layer to the deep layer, the current average geothermal gradient of Xihu sag is 30.5 ℃/km, and the current formation temperature is about 167 ℃ at 5000 m buried depth. The relationship between formation temperature and depth ( Fig. 3) ( Table 2) can be fitted as T = 0.0305*H + 15, R 2 = 0.8726, in which T is formation temperature (℃), H is the buried depth (m), and the constant 15 is the average surface temperature. The present geothermal gradient in Xihu sag is higher in the Central inversion structural belt with 33.4 ℃/km and lower in the Western slope belt with 27.8 ℃/km and the Western sub sag with 30.5 ℃/km, in general lower than 3-5.5 ℃/100 m (Fig. 3A).

X-ray diffraction (XRD) of clay minerals.
It was showed that the I/S-S% of Oligocene Huagang Formation in Xihu sag obtained by X-ray diffraction ranges from 15 to 20% (Fig. 3D) ( Table 2), indicating that I/S layer has partially transformed into illite, and the source rock has reached a high maturity stage.
Vitrinite reflectance Ro%. According to the measured reflectance of vitrinite formation (Ro%) in Xihu sag ( Fig. 3C) ( Table 2), the value of Ro% has a significant increase trend with the burial depth. However, the value of Ro% is only about 1.3% when the burial depth of the stratum reaches 4500 m, and while the burial depth of the stratum is more than 5000 m, most of the value of Ro% is not more than 2%. This result showed that the rate of diagenetic evolution of Xihu sag is slow, while the degree of organic matter evolution is only in the middle diagenetic stage, and the maturity of organic matter is mostly in the low maturity maturity stage. In addition, through the analysis of the diagenetic evolution stages of different gas fields in the Western slope belt (A5), the Western sag (A10) and the Central uplift belt (A9) (Fig. 4), it is believed that the stratum is still in stage IIA at a burial depth of less than 4000 m and only start to evolve to stage IIB below 4000 m.
Homogenization temperature of inclusions. The homogenization temperature (Th) of inclusions captured in quartz grains and secondary enlarged edge ( Fig. 5M-O) can represent the paleotemperature of the reservoir, which can be well applied to the division of diagenetic stage 22 . The homogenization temperature analysis of the reservoir inclusions in the Western sag and Central uplift belt showed that most inclusions were formed www.nature.com/scientificreports/ in the stage IIA, and the maximum homogenization temperature is not more than 175 ℃ (Fig. 3B)  Maximum pyrolysis temperature Tmax. The maximum pyrolysis temperature of mudstone in Xihu sag is ranging from 402 to 479 ℃ (Fig. 3E) ( Table 2). According to the classification standard of diagenetic stage of clastic rock, it can be concluded that the reservoir at present day is mainly in stage IIA-IIB.

Discussion
Comprehensive analysis of diagenetic stage. It can be seen under the microscope that the plastic components (biotite, muscovite, feldspar) have been deformed and oriented by strong mechanical compaction (Fig. 5A,B). In addition, the contact relationships between the detrital particles (line contacts, line-to-concavity contacts, concavity-convexity contacts, and suture line contacts) (Fig. 4C) also indicate that the Huagang formation is in a relatively high diagentic stage. The carbonate cementation mainly involved the calcite, dolomite, ferro calcite, ferrodolomite ( Fig. 4D-F), and ankerite is a sign of late diagenetic stage 16 .
In consideration of formation temperature (Fig. 3A), temperature measurement of quartz mineral inclusions (Fig. 3B), vitrinite reflectance (Ro%) (Fig. 3C), clay mineral X-ray diffraction (XRD) (Fig. 3D), organic matter pyrolysis peak temperature (Tmax) (Fig. 3E), as well as the evolution characteristics of authigenic minerals, pore  www.nature.com/scientificreports/ structure, grains contact relation observed under the microscope (Fig. 5), the reservoir diagenetic stage was identified. According to the classification standard of diagenetic stage of clastic rock in China's oil and gas industry, the reservoir of Huagang Formation in Xihu sag is mainly in stage IIA-IIB. It was showed that 2300-3200 m is in the stage IIA1, while 3200-4300 m is in the stage IIA2 and 4300-5000 m is in the stage IIB (Fig. 3).   Table 3).
Analysis of thermal history. Previous researchers [43][44][45]    www.nature.com/scientificreports/ of Oligocene, ranging from 63.8 to 98.6 mW/m 2 , with an average of 83.4 mW/m 2 . It is a continuous heat flow reduction process from Oligocene to present. Based on the above understanding, the general model of terrestrial heat flow in the thermal history simulation of Xihu sag is shown in Fig. 7. It can be seen from the figure that the deposition period of Huagang Formation belongs to the early stage of the compression sag stage, with the geothermal flow of 80.5 mW/m 2 , which gradually increased later. In the rifting stage, when the Longjing formation (N1L) was deposited, the geothermal flow reached the maximum of 90.0 mW/m 2 . Then, the geothermal flow began to decline gradually. During the deposition period of Liulang formation (N1ll), the geothermal flow decreased to 71.2 mW/m 2 . In the regional subsidence stage, the fault activity is weakened, and the heat flow of the earth is further reduced, which is maintained at about 60.0 mW/m 2 (Fig. 7).
The paleotemperature in the study area is generally distributed between 85 and 180 ℃, which is manifested in the rapid thermal evolution in the CUB I. During the deposition period of Longjing formation (N1lj), the formation temperature has exceeded 120 ℃, which is in stage IIA2, and then the temperature changes slowly with temperature variation rate less than 3 ℃/Ma. Until the deposition of the upper Pleistocene Santan formation (N2st), the temperature changes rapidly with temperature variation rate more than 5 ℃/Ma, and the formation rapidly evolves to the middle diagenetic stage IIB. In the CUB II and WS, the formation temperature has exceeded 85 ℃ in the deposition period of Longjing formation (N1ll), which is in stage IIA1. Then the temperature changes slowly all the time with temperature variation rate less than 2 ℃/Ma. By the deposition of the Santan formation (N2st) of the Pliocene, the temperature changes rapidly with temperature variation rate about 2 ℃/Ma, and the formation rapidly develops to the stage IIA2. The area with the lowest degree of thermal evolution is the west slope zone, and the simulated current temperature is generally low, ranging from 85 to 110 ℃ (Fig. 8).  www.nature.com/scientificreports/ Coupling relationship between thermal history and structural evolution. The rifting process of Xihu sag lasted from paleocene to the end of oligocene, and then transformed into thermal subsidence in Miocene, and a new round of rapid subsidence started again in Pliocene Pleistocene 23,33 . The time of rapid tectonic subsidence in Xihu sag corresponds to the highest paleoheat flow and the rapid rise of paleogeotemperature, which reflects the control effect of tectonic activity on the geothermal state of the basin.
Diagenetic parameter simulation. Transformation characteristic of clay minerals. I/S mixed layer are widely distributed in mudstone and sandstone. The proportion of I/S-S% is an important means to judge the diagenetic strength and temperature of clastic rocks 21,38,39 . The transformation process from montmorillonite to illite mainly depends on time, temperature and pressure. With the increase of time, temperature and pressure, the transition sequence from montmorillonite to disordered I/S to ordered I/S and illite content increases gradually. The proportion I/S-S% decreased during burial diagenesis and geothermal alteration. The simulated value of I/S-S% in the study area ranges from 9.4 to 47.7%, with an average of 24.2%, which is in good agreement with the actual measured value in the study area. The value of I/S-S% indicates that the Huagang Formation in Xihu sag is now in stage IIA-IIB. In general, the I/S-S% value shows the following characteristics: (1) it changes slowly with a variation rate less than 2.5%/Ma in the early stage of burial and decreases rapidly with a variation rate about 15%/Ma when the depth exceeds 1500 m. The CUB I has the largest change range about 70%, followed by the CUB II with about 60% and WS with about 55%, while the Western slope has a relatively small change range with about 35%. In the later stage of burial, the change range of simulation value decreases again; (2) the older the strata in the buried history, the lower the I/S-S% value, and the I/S-S% of the CUB is lower than other areas with a value about15%, indicating relatively strong diagenesis in the CUB; (3) the uplift of the stratum does not affect the transformation of montmorillonite, but the change range is reduced. Previous studies 20,21 have shown that the transformation rate of montmorillonite to illite is closely related to the geothermal gradient. When the temperature is high, shallow burial can lead to rapid transformation, and the converse is also true. In addition, the transformation from montmorillonite to illite is also very fast in the water medium of the reservoir under the condition of rich potassium. Among various influencing factors, strong basic thermal activity and high geothermal gradient are the main reasons for low I/S-S% (Fig. 9).
Evolution characteristics of vitrinite reflectance. The simulated reflectance (Ro%) of vitrinite ranges from 0.51 to 1.45%, with an average of 0.81% (Fig. 10). Ro% indicates that the reservoir of Huagang Formation in the study area is in IIA-IIB. In general, Ro% value shows the following characteristics: (1) it changes slowly with a variation rate about 0.025%/Ma in the early stage of burial, and the simulation value increased rapidly with a variation rate about 0.075%/Ma when the depth is more than 2000 m. The CUB I has the largest change range, followed by www.nature.com/scientificreports/ CUB II and WS, and the WSB has a relatively small change range with a variation rate about 0.05%/Ma. In the later stage of burial, the change range of simulation value decreases again. (2) the older the strata in the buried history, the higher the Ro%, and the Ro% of the Central uplift belt in the study area is the largest with a value about 1.2%. The results show that the maturity of organic matter is high and diagenesis is relatively strong; (3) the uplift of strata does not affect the thermal evolution of organic matter, and the evolution of vitrinite formation is irreversible. Vitrinite reflectance is mainly controlled by temperature, and also affected by diagenetic time, tectonic activity and magmatic activity 17-37. The fine subsection of diagenetic stage (DSF). In the past decades, many scholars at home and abroad 6,10,17,18,46 have made extensive and in-depth research on the diagenesis of clastic rocks, and put forward corresponding diagenetic stage division standard of clastic reservoir. In the early 1990s, on the base of the domestic and foreign scholars 6,10,47 on diagenesis research and diagenetic stage division, Ying Fengxiang (1992) 48 combined with the specific geological characteristics of sedimentary basins in China systematically and brought forward the specification for diagenetic stage division of clastic rocks. In this division scheme, the diagenetic stage of clastic rocks was based on the thermal evolution model of organic matter, distribution characteristics of clay minerals and evolution law of sandstone authigenic minerals and divided into early diagenetic stage and late diagenetic stage. Another scheme from Wilson M. D. (1994) 10 was based on the study of clastic rock diagenesis of many scholars in the world and on account of organic matter thermal evolution, which divided the organic matter thermal evolution into immature stage, decarboxylation stage, liquid hydrocarbon generation stage and the natural gas generation stage to match the corresponding diagenetic temperature and autogenetic mineral evolution sequence for different organic matter thermal evolution stages. At the beginning of the twenty-first century, Ying revised the classification standard of the early diagenetic stage of the clastic rock and divided the diagenetic evolution of the clastic rock into three stages 19 . In this new scheme, the intermediate diagenetic stage was subdivided into two stages, A and B, and the intermediate diagenetic stage A was further divided into two sub stages, A1 and A2 (Table 1). This division standard of diagenetic stage of clastic rock has been used up to now. The accurate division of diagenetic stages brings insight into the diffidence of paleo temperature (T), vitrinite reflectance (Ro%), the highest pyrolysis peak temperature (Tmax), thermal alteration index of the proportion of smectite in illite/smectite interstratified minerals (I/S-S%) and authigenic minerals assemblage in both time domain and space domain 11,24,49 . In order to accurately describe the diagenetic evolution process of Huagang Formation, a fine division scheme of diagenetic stage of clastic reservoir (Table 4) with 36 diagenetic micro-stage (DS) in diageneitic process was proposed based on the paleotemperature, vitrinite reflectance and the proportion www.nature.com/scientificreports/ of smectite in illite/smectite interstratified minerals. DS1-DS9 is the equal of IA, while DS10-DS13, DS14-DS19, DS20-DS24, DS25-DS31 and DS32-DS36 are equivalent to the stage IB, IIA1, IIA2, IIB and III separately. This fine division scheme can be used for numerical simulation of reservoir diagenetic evolution, which can highlight the differences of reservoir evolution process and help to identify the key geological factors linked to the differential evolution of high diageneisis reservoir. The simulation results of Paleotemperature (T), vitrinite reflectance (Ro%) and I/S-S% show that the reservoir of Huagang Formation in the study area is mainly in stage IIA1, and stage IIA2, and some reservoirs in the north Central uplift belt and Western sag such as A10 well area have entered stage IIB (Fig. 11). The diagenetic stages indicated by different diagenetic parameters have the following characteristics: (1) in the shallow burial phase (H < 1000 m, T < 60 ℃), the diagenetic stages divided by paleotemperature (DS-T) and smectite in mixed layer (DS-I/S-S%) are equivalent or close and they are higher than that of diagenetic stages divided by vitrinite reflectance (DS-Ro%); (2) as the depth of stratum increased (H < 2000 m, T < 100 ℃), the difference between DS-T and DS-I/S-S% increased at the same time, which is shown as that DS-T is significantly greater than DS-I/S-S% and DS-Ro%, and DS-Ro% is minimum; (3) when the burial depth is further increased, DS-I/S-S%, DS-Ro% and DS-T tend to have the same value DS. Most of the well areas appear in stage IIA1, and some well areas such as A8 appear in stage IIB. The depth and location of DS are different in different structural areas; (4) after the poit DS, with the increase of burial depth, DS-Ro% no longer changes, or changes very little. However, DS-I/S-S% varies greatly with the increase of burial depth, and DS-I/S-S% is significantly larger than DS-T and DS-Ro%, and DS-Ro% is minimum. However, exceptional conditions occured in well A10. After DS stage, DS-I/S-S% changed little with the increase of burial depth, and DS-T was still the largest, followed by DS-Ro%, and the lowest is DS-I/S-S%. According to the measured data from the oil field (Fig. 12), there are abnormal overpressure environments in the well A4, well A11 and well A5. The initial depth values of abnormal overpressure are 3800 m, 3900 m and 4200 m respectively, and with the increase of depth, the formation pressure coefficients are measured to be up to 1.8. Previous studies 50,51 on the evolution of organic matter and the transformation of clay minerals in overpressure environment have shown that the increase of pressure can inhibit the evolution of organic matter and hydrocarbon generation, and reduce the transformation of montmorillonite to illite 52,53 . The abnormal overpressure may begin to develop at the burial depth of 3000 in well A10. The existence of the abnormal overpressure restrains the transformation of montmorillonite to illite, leading to the DS-I/S-S% value relatively smaller after DS; (5) in the diagenetic process, with the increase of temperature, the transformation of montmorillonite to illite is more sensitive than the evolution of organic matter, and DS-I/S-S% can effectively reflect the diagenetic process. In the shallow burial stage, the description of diagenetic evolution stage is relatively accurate compared with vitrinite and paleotemperature, while in the deep burial process, due to the existence of abnormal overpressure, DS-I/S-S% is relatively small.
It is generally believed that time and temperature are the main reaction processes in the evolution of diagenesis. Therefore, the rapid deposition and burial of strata can slow down diagenesis. Diagenesis is not synchronous with temperature. The current diagenetic stage is a comprehensive product of the whole geological history. When different indicators are used to judge the diagenetic stage, the prediction results are different in time and space (Fig. 11). Because the depth, temperature, duration, pressure and some chemical variables can change synchronously or asynchronously, it is difficult to determine the parameters that affect the evolution series of apparent diagenesis stage. Therefore, in order to accurately describe diagenesis and achieve the goal of successful oil and gas exploration, a variety of comprehensive prediction indexes should be needed. www.nature.com/scientificreports/

Conclusions
(1) Based on the comprehensive utilization of formation temperature data, inclusion temperature data, vitrinite reflectance test data, organic thermal solution peak temperature data and clay mineral X-ray diffraction data, combined with the analysis of rock physical properties, particle contact mode, pore structure, pore type and other characteristics, the diagenetic evolution stage of Huagang Formation reservoir was analyzed. It was showed Figure 11. Comparative analysis of diagenetic stages using different diagenetic parameter. www.nature.com/scientificreports/ that the reservoir of Huagang Formation in Xihu sag is mainly in stage IIA-IIB. Among them, 2300-3200 m is in stage IIA1, 3200-4300 m is in stage IIA2 and 4300-5000 m is in stage IIB.
(2) The burial history and thermal evolution history of the study area are restored. On this basis, paleogeotemperature, vitrinite reflectance and smectite in illite/smectite interstratified minerals were selected to simulate the evolution of Huagang clastic rocks. The simulation results showed that most of the strata in the upper part of Huagang Formation (HS) of Xihu sag are in the stage IIA1-IIA2. The CUB I is deeply buried, and some of the strata have entered the stage IIB. The lower part of Huagang Formation (HX) in Xihu sag is in the stage IIA1-IIB, in which the WSB is in the stage IIA1, while the WS and CUB III are mainly in the stage IIA2, and the CUB I and CUB II is in the stage IIB.
(3) A fine subsection scheme of 36 diagenetic micro-stage (DS) in diageneitic process was proposed based on the paleotemperature, vitrinite reflectance and the proportion of smectite in illite/smectite interstratified minerals. DS1-DS9 is the equal of IA, while DS10-DS13, DS14-DS19, DS20-DS24, DS25-DS31 and DS32-DS36 are equivalent to the stage IB, IIA1, IIA2, IIB and III separately. This fine division scheme can be used for numerical simulation of reservoir diagenetic evolution, which can highlight the differences of reservoir evolution process and help to identify the key geological factors linked to the differential evolution of high diageneisis reservoir.
(4) When judging the diagenetic stage of reservoir in geological period based on a single index, the prediction results of different parameters are quite different. This difference reflects the influence of structure and abnormal pressure on diagenetic evolution. In order to accurately describe the diagenetic evolution process of Huagang Formation, three parameters of T, Ro% and I/S-S% were used for a fine discrimination model of diagenetic stage.