Porosity Formation and Evolution of the Deeply Buried Lower Triassic Feixianguan Formation, Puguang Gas Field, NE Sichuan Basin, China

The deeply buried shoal dolomite reservoirs of the Lower Triassic Feixianguan Formation, giant Puguang Gas Field, NE Sichuan Basin, are exceptionally porous. The influences of diageneses on pore evolution are studied. Through petrologic investigation, diagenetic phases are divided into four stages, i.e., near-surface, pre-oil window, oil window, and gas window. The Adobe Photoshop system is used to quantify the rock texture components, porosity constitutions and the influences of diageneses on reservoir porosity. Porosity evolution curves are quantitatively recovered. The Feixianguan reservoir porosities are mainly created by early meteoric dissolution, dolomitization, and organic acids dissolution. Dissolution during deep burial is insignificant. Pores are formed in near-surface and pre-oil window stages and effectively preserved till present. This result may be of great significance to the further exploration of deeply buried carbonate reservoirs not only in NE Sichuan Basin, but also around the world.


Introduction
showed a trend of carbonate porosity decreasing with burial depth and age based on Statoil's global reservoir database. However, a notable exception to that is the Triassic [1]. One of the possible reasons for this abnormity, according to the authors, is the limited data. Thirty four percent of the overall 284 Triassic carbonate reservoirs are from China, mainly the Jialingjiang Formation at a depth of <2.3 km, with porosities <5.1% [1]. However, the recently discovered Puguang Gas Field [2,3] is not included in the database. Data from this new gas-field may help evaluate the Triassic abnormity.
The giant Puguang Gas Field discovered in 2003 is the largest gas field found in marine carbonates within China with an estimated gas reserve of 350 × 10 9 m 3 (1236 BCF) [2]. The oolitic dolostone of the Lower Triassic Feixianguan Formation is one of the major reservoir rock types with the maximum porosity up to ~30%. This reservoir is exceptionally porous considering the deep burial depth (4800 -5500 m), high temperature, thermally overmature and multiple-stages of tectonic deformation [2][3][4]. How these pores formed and evolved are the major concerns of exploration as well as of this research.

Geological Setting
The Puguang Gas Field is located on the Shuangshimiao-Puguang structural belt in the NE Sichuan Basin, Southwest China (Figure 1; [2]). The Sichuan Basin has a complex sedimentary and tectonic history due to several important tectonic episodes or "movements" since the Late Proterozoic [2][3][4]. As a result of tectonic evolution and eustacy, strata within and below the middle Triassic Leikoupo Formation (T 2 l) are generally dominated by carbonate lithologies. Clastic rocks were deposited during and after the Upper Triassic Xujiahe period. The Upper Permian Changxing Formation and the Lower Triassic Feixianguan Formation are two major reservoir intervals of the Puguang Gas Field and many other gas-fields in the eastern Sichuan Basin [2].
The Feixianguan Formation is subdivided into four members, namely the first to the fourth member from base to top, of which the effective reservoirs are in the first to the third members. Sediments in the first and second members are dominated by oolites interbeded by mudstone carbonates deposited in platform-margin shoal; while the third to the fourth units are mainly micritic to fine dolomites and micritic limestones with interlayered anhydrites in the fourth member [2]. Notably, evaporates of different thicknesses (several to ~100 m thick) developed in the interior of the first to the second member in the Tongjiang-Kaixian area, which indicates an evaporatic platform was developed adjacent to the oolitic shoal (Figure 1). Dolomites are heterogeneous in thickness, which are dominantly distributed around the Puguang area. In the Well PG2, oolitic dolomite and crystalline dolomite comprise the main reservoir section of about 400 m (1312 ft.). Above the Feixianguan Formation, the Jialingjiang Formation and the Leikoupo Formation are deposited as thick and basinwide anhydrite, halite, and gypsiferous dolomite, which provide regional seals to the Puguang and other gas-fields [2].
The Feixianguan carbonates reached maximum burial depth of about 7000 m (22966 ft.) at ~120 Ma before being uplifted to the present depth of about 4800 -5700 m (15748 -18700 ft.) and the calculated present-day thermal maturity is approximately between 2.0% and 2.5% R o [4]. Although deeply buried, the Feixianguan reservoir exhibits high porosity up to 29% and high permeability 9664 md.
Natural gas produced from the Feixianguan reservoir in the Puguang field is typically dry, with approximately 80% methane, less than 0.5% C 2+ alkanes, and is rich in H 2 S (5.1% -62.2%, 19.7% in average) and CO 2 (2.6% -35.6%, 11.9% in average) [3,4]. The gas geochemistry and the presence of extensive pyrobitumens in the reservoir provide clear evidence for an early oil charge and formation of a paleo-oil reservoir during the Triassic-Jurassic late Indosinian-early Yanshanian orogenies. During the Late Jurassic-Cretaceous Yanshanian Orogeny, oil in the paleo-trap was cracked into thermal gas under higher temperatures.

Methods
Twenty eight cores and four outcrops were investigated. Six hundred and thirty samples were taken from the core and the outcrop material perpendicular to bedding surfaces at 0.2 -2.0 m intervals. Three hundred and sixty polished thin sections were prepared after the samples were impregnated with purple epoxy for easy porosity identification. The thin sections were stained with Alizarin Red-S to identify calcite and dolomite. Optical investigation was conducted conventionally under microscope. On principles of superposition and crosscutting, relative timing of diagenetic phases was determined.
Rock texture components and porosities are quantified by Adobe Photoshop Analysis [5]. This method is based on microscopic images taken from stained and impregnated thin sections. One can easily and precisely select an area of interest in those images with the selecting tools of Adobe Photoshop. In other words, grains, cements, matrix and porosities of different types can be selected separately in carbonates. Then, the pixels of the selected zone can be read in the menu of Image-Histogram-Pixels. The ratio of the pixels of the selected zone to that of the whole image equals to the percent area of the selected component. Consequently, the rock texture components and porosities of different types can be quantitatively determined. Figure 2 gives an example of the application of the Adobe Photoshop Analysis method, which shows good reproducibility of the results. In fact, this method is comparable to the commonly used methods such as estimating through naked eyes, visual comparison charts [6], point-counting [7,8], or usual image analysis [9], but more precise and easier to handle. In fact, the Adobe Photoshop Analysis method can also be used for modal analysis of granitic rocks [10] and quantitative assessment of sandstone porosity [11].
A total of 360 polished thin sections of both lime-stones and dolomites were imaged at least 6 times for each sample. Approximately 3000 images were analyzed by the Adobe Photoshop Analysis method. These data make the quantitative evaluation of porosity evolution of the reservoir rocks possible.

Results
Major diagenetic processes are identified as meteoric dissolution and cementation, early dolomitization, mechanical compaction, chemical compaction and related burial cementation, fracturing, hydrocarbon emplacement, oil cracking, pyrobitumen precipitation, gas accumulation, thermochemical sulfate reduction, and quartz crystallization during the late uplift of the Feixianguan reservoir. Porosity types are classified as fabric-selective, partial-fabric-selective, and non-fabric selective using the classification criteria of Choquette and Pray (1970) [12].
The results show that most of the pores of the Feixianguan reservoir are fabric selective intragranular (including moldic porosity), intergranular or intercrystalline, partly non-fabric selective vugs plus minor fractures, breccias and carverns. As shown in Figure 3, pores in Well PG2 are mainly intragranular and intercrystalline, partly intergranular. The reservoir quality is largely controlled by the abundance of intragranular and intercrystalline pores. While in Well MB4, pores are mainly intergranular, intercrystalline and vugs, among which intergranular and intercrystalline porosities determine the reservoir quality (Figure 3).

Diagenetic Processes
The paragenetic sequence of the Feixianguan Formation was established using the criteria of Heydari (1997) [13]. Diagenetic phases are divided into four stages: 1) nearsurface (including marine and meteoric environment); 2) pre-oil window; 3) oil window; and 4) gas window (wet gas and dry gas). Each stage is related to a distinct phase  (Figures 4 and 5).

Stage 1: Near Surface Diagenesis
The Triassic marine water had high tion and was prone to precipitate aragonite, which is called the "aragonite sea" [14,15]. Consequently, the Triassic ooids in the NE Sichuan were likely aragonitic. During and slightly after deposition, the ooids were cemented and lithified in the marine environment. Generally, the marine cements developed mainly as micrite envelopes or isopachous crusts of aragonite on the surfaces of ooids (Figure 5(a)). The marine cementation had a negative influence on the reservoir pore evolution which could be quantified as the contents of the marine cements.
In the hypersaline environment, diageneses were characterized by penecontemporaneous dolomitization, as well as the related anhydrite nodularization or cementation ( Figure 5(b)).
Due to the frequent sea-level fluctuations, the Feixianguan oolitic shoal was probably subaerially exposed after precipitation. And approximately at the end of the second period of Feixianguan, a regionally short-time exposure occurred. Unstable aragonitic ooid grains were leached or altered by neomorphism, and meteoric calcite cements precipitated (Figure 5(a)). Dissolution features were characterized by oomoldic or intragranular pores with geopetal structure (Figure 5(a)), dissolution-enlarged fractures, carvens and collapse breccia pores. However, gravitational and meniscus cements and microspar are hard to identify. d by the pre-bitumen [16] oomoldic or intragranular pores with geopetal structure, dissolution-enlarged fractures, carvens, and collapse breccias. Nevertheless, meteoric cements are hard to identify because later dolomitization altered the rock textures. It is thought that meteoric diagenesis is relatively porosity neutral [6]. However, meteoric dissolution pores constitute a major part of the present reservoir space. And the meteoric cementation not only excluded pores, but also provided hard framework of the strata and effectively resisted compaction and porosity loss.
The near surface stage generally continued to the early shallow burial and ended before the onset of pressuresolution.
When and how the first and the second members of Feixianguan oolitic dolomites formed are still controversial. Major viewpoints include: 1) dolomitization by re fluxed brines (e.g., [17]) or 2) by meteoric/mixed water (e.g., [18]) in the early diagenetic stage, and 3) burial dolomitization by formation fluids in the late diagenetic stage (e.g., [19,20]). We found several pieces of evidence for judging the relative timing of dolomitization. Important considerations are: cementation tic 1) Pyrobitumen is present in the intercrystal vuggy ores of nearly all the dolomites with different crystal 4) Oomoldic pores with geopetal structure in the oolidolomite formed during early meteoric dissolution prior to dolomite formation. Otherwise, dolomitization should have stabilized the aragonite or calcite grains to dolomite, which would not be leached fabric-selectively by meteoric water.
In conclusion, oolitic dolomite should have been formed at a shallow p size (from micritic, fine, coarse, to sucrosic), which shows that most of the dolomites formed prior to the hydrocarbon emplacement (Figure 5(g)).
2) Coarse dolomite cements developed in vugs or carvens. However, few oil inclusions are found in these cements, while oil impregnation and pyrobitumen coating are common, which reconfirms that dolomite formation, as well as dolomite recrystallization, predated oil accumulation.
3) The main features of physical or mechanical compaction and che burial depth (<500 -800 m, or 1640 -2625 ft.) postdating meteoric diagenesis and predating intensive compaction. However, dolomite may be recrystallized due to increasing burial depth to about 2000 m (6562 ft.) before petroleum charge.
After dolomitization, fractures occurred to greatly increase the reservoir permeability. Som mical compaction (also pressure solution) include dewatering and deformation or reorientation of grains, stylolites or solution seams [6]. Stylolites and seams are common in limestones while scarce in dolomites in NE Sichuan (Figure 6). Few stylolites and no intergranular microstylolites could be found in the oolitic dolomite. Moreover, individual ooid shows little deformation. Therefore, it can be deduced that dolomitization predated stylolite formation and dolomite as the rock framework resisted intensive physical and chemical compaction.
Generally, incipient stylolitization may start at a burial depth of abou e of the fractures acted as the passageways of hydrocarbon-related fluids (Figure 5(f)).
More and more evidences have shown that dolomitization may directly or indirectly influence the reservoir porosity ( [22] and references therein). It may directly cause: 1) dissolution by dolomitizing fluids which may be of little significance on pore formation; 2) mole-mole replacement of Ca 2+ by Mg 2+ which constructs porosity; 3) over dolomitization by excess fluids which destructs porosity; 4) volume-volume replacement of calcite or aragonite by dolomite which has no influence on total porosity; and 5) improvement of permeability nearly in all cases. The indirect influences are caused by the different mineral components and rock textures between t 500 m (1640 ft.), and become intense below 830 m (2723 ft.) [21]. Consequently, dolomites may have been formed before they were buried to 500 -800 m (1640 -2625 ft.). dolomites and the precursors. After dolomitization, the emaining calcite minerals are prone to be reconstructed structural configurations and burial history r selectively porosity [16,23]. However, to qu O rocks entered oil wi of the strata over geological time [3,4]. The paleo-oil ochemical sulfate reduction (TSR) [4], sh han or entered into gas wi e re d gases within the Feixianguan reservoir [4]; dissolved; while the rocks become incompatible and fissile as well. In this area, pores are mainly developed in dolomites. Consequently, dolomitization is considered as a main constructive factor for the reservoir formation. Most of the intercrystalline porosities were created by dolomitization and recrystallization and can be quantified by Adobe Photoshop Analysis. As evidence, dolomite crystal sizes are positively correlated to the intercrystalline porosities (Figure 7).
Generally, both mechanical and chemical compactions are destructive to reservoir antify the porosity loss caused by compaction is difficult, even impossible. Great attempts have been made by Heydari (2000) [16] and the citations therein. Through point counting, Heydari (2000) draw the conclusion that out of an original 40% porosity, 12.8%, 14.7%, and 12.5% were destroyed by mechanical compaction, intergranular pressure solution (IPS), and related cementation, respectively [16]. No IPS, relatively weak physical and chemical dissolution have been found in the Feixianguan oolitic dolomites; while compaction is much more intense in crystalline dolomites and especially in limestones (Figure 6). The porosity loss by compaction is qualitatively estimated based on the density of stylolite, the content of grains and many previous studies [16,23].

Stage 3: Oil Window Diagenesis
In the late Indosinian Orogeny to the early Yanshanian rogeny, the Lower Silurian source ndow [3]. Organic acids were firstly expelled and may have partially dissolved the oolitic rocks and improved the reservoir quality [12]. Oil charge followed. In Well PG2, the earliest oil migration occurred at about 198 Ma with burial depth at ~2800 m (9186 ft.) according to the field formed and was preserved for a long time until oil cracked [2,3].
An important reaction, namely, liquid-hydrocarboninvolved therm ould be noticed in this stage. Little evidence has been proposed for this reaction. However, sulfur-rich pyrobitumen and elemental sulfur in the reservoirs may indicate that liquid hydrocarbon did react with sulfate but little H 2 S gas was formed during this stage of reaction [4].

Stage 4: Gas Window Diagenesis
In the middle and late episodes of the Yans ogeny, hydrocarbon source rocks ndow for persistent subsidence and deep burial. Gases were injected into the reservoir to form a giant gas pool. At the same time, in-reservoir thermal maturation happened to the early charged oil [3]. The petroleum fluids in the paleo-oil field were altered by thermal cracking and thermochemical sulfate reduction (heavy-hydrocarbon dominated to methane dominated stages [4]) due to increasing burial depth and temperature. Abundant sulfur-rich pyrobitumens (Figure 5(g)), very dry gases with varying H 2 S and nonhydrocarbon gas contents, and cements like elemental sulfur and calcite (Figure 5(h)) were the combined results [4]. During the thermochemical sulfate reduction, the reservoir pore structures may have been modified and porosity perhaps increased [3].
During the Himalayan orogeny, the gas trap was readjusted to present configuration. Thermochemical sulfat duction might be terminated because of the decreased burial depth and temperature caused by uplift and denudation [3].
Thermochemical sulfate reduction has largely altered the fluids an Copyright © 2013 SciRes. OJG X. F. ZHANG ET AL. 308 while porosity variation during this process is controversial. Many previous studies may emphasize the importance of thermochemical sulfate reduction-related dissolution [24][25][26][27][28][29]. However, this view-point was doubted because wells with little H 2 S gas indicating weak thermochemical sulfate reduction were drilled [30]. This research provides several pieces of evidences which show that thermochemical sulfate reduction negligibly influenced the total porosity despite local minor dissolution and contemporaneously local cementation by sulfur and calcite happened.
1) Most (about 2/3) of the pores in the Feixianguan reservoirs of Well PG2 and Well MB4 are fabric-selective. Thermochemical sulfate reduction is a deep burial reaction postdating oil accumulation. Minerals are supposed to have been stabilized and should not be selectively dissolved. Consequently, pores formed by therochemical sulfate re ld m duction are relatively less important in the study area.
2) Generally, pores are partly cemented by dolomite on the inner wall. Most of these dolomite cements are euhedral without signs of later dissolution.
3) Most of the pores are rimmed by bitumen, which means that these pores are generated before oil emplacement and paleo-oil field formation.
4) The dissolution always co-exists with the cementation by thermochemical sulfate reduction byproducts. Porosity is probably produced from the dissolution of calcium sulfate to yie 2 4 SO  , the generation of CO 2 , an emical sulfate reduction. Consequ of quartz cem ( Figures 3A and 5(b)), while oolitic dolomite of the Well MB4 has abundant vugs, intergranular and intercrystaline pores (Figures 3B and 5(c)). At the near-surface stage, the Feixianguan Formation of the Well PG2 experienced intense meteoric dissolution and abundant oomoldic/ intragranular pores formed (Figure 5(b)). The meteoric dissolution would be macroscale dissolution caused by large water flux and water strongly undersaturated with respect to aragonite [15]. Porosity, as well as permeability, was then greatly intensified by dolomitization, recrystallization and organic dissolution. After oil emplacement and pyrobitumen cementation, few pores were created. In the contrary, pyrobitumen, elemental sulfur, calcite and minor quartz cements filled part of the pore space.
Even though the Feixianguan Formation of the Well MB4 was similarly exposed as did in the Well PG2, early meteoric diagenesis may have exerted different impacts l MB4. The meteoric stabilized and proceed by m he Well MB4 (Figures 3A and 5(c)). Following the meteoric diagenesis, dolomitization and crystalline pores and d fo d the precipitation of metal sulfide; while lost due to the precipitation of authigenic carbonates, elemental sulfur and/or sulfides [31]. Microscopic observation and Adobe Photoshop Analysis on the Well PG2 have shown that pores possibly created by thermochemical sulfate reduction dissolution are always less than or equals to those excluded by elemental sulfur and calcite precipitated during thermoch ently, the total porosity may stay unchanged or slightly decreased during this reaction.
In this diagenetic stage, small amounts ented in some pores (Figure 5(h)). And fracturing occurred as the latest diagenesis which has little impact on the total porosity.

Diagenetic Heterogeneity
Diagenetic features and intensity are heterogeneous in different regions of NE Sichuan Basin or even in different reservoir sections of the same well core.
The Feixianguan Formation was deposited at the similar facies in the Well PG2 and in the neighbor Well MB4 and experienced similar paragenetic sequence. However, oolitic dolomite with abundant oomoldic/intragranular and intercrystalline pores developed in the Well PG2 or fresh-marine-mixed water would be slightly undersaturated with respect to aragonite and the water flux was slow, aragonite grains may be on the reservoir sections of the Wel icroscale dissolution or replacement rather than be dissolved [15]. Consequently, oomoldic or intragranular pores are scarce, while original intergranular pores are well preserved in t recrystallization formed the inter greatly increased the permeability. Porosity slightly decreased after oil emplacement because of the pyrobitumen, sulfur, calcite and quartz cementation.
Dolomitization may be a pre-requisite to form an effective reservoir of the Feixianguan Formation in the NE Sichuan, which is proved by the Well YB2 from the opposite angle. The Feixianguan oolites in this well are not dolomitized. The present porosities are mainly less than 3% (Figure 6). Although early meteoric water dissolved some ooid grains, these oomolds were filled by burial cements. Due to the weak meteoric dissolution and the absence of dolomitization, most of the original and secondary pores were excluded by mechanical compaction, pressure solution, and burial cementation.

Porosity Evolution
To recover the porosity evolution history, an estimation of original porosity value is required. Several previous studies have estimated theoretically or experimentally the original porosity value of grainstones which is controlled by sorting and grain packing arrangement ( [16] and citations therein). An original porosity of 40% was assume r the Smackover ooid grainstones by Heydari (2000) [16]. In the Well PG2, oolites have an average grain content of about 66%, among which oolite with the best-preserved original texture and the least compaction (Well PG2, 5030.5 m or 16,504 ft.) has 60.6% of grain contents. Therefore, a 40% of original porosity seems reasonable for the Feixianguan oolites in NE Sichuan.
On the basis of the paragenetic sequence (Figure 4) and the analyses of the influences of diageneses on reservoir porosity, porosity evolution history could be (semi-)quantitatively evaluated. The recovered curves start from the present porosity (P p ) quantified by Adobe Photoshop Analysis. As the latest diagenesis is late fracturing which neatly increased porosity, the porosity before late fracturing (P p1 ) should be calculated as P p minus porosity created by late fracturing. And the porosity before quartz ce ents to get the tation (P o1 ). The P o1 will be ric dissolution and cementaes an his research may include early mesearch: mentation (P p2 ) should be computed as P p1 plus quartz contents. The P p2 plus thermochemical sulfate reduction cements (elemental sulfur and calcite) and minus porosities generated by thermochemical sulfate reduction dissolution equals to the porosity before thermochemical sulfate reduction (P p3 ). The P p3 added by bitumen contents is considered as the porosity when oil charged (P p4 ). Since the contribution of organic dissolution is hard to determine, the pre-bitumen porosity evolution could only be semi-quantitatively recovered started from the original porosity (P o ). The P o subtracts marine cem porosity after marine cemen increased to P o2 after meteo tion, therefore, P o2 equals to P o1 plus meteoric intragranular, breccia, carven and dissolutionenlarged porosities. Porosity after dolomitization and recrystallization (P o3 ) should be calculated as P o2 plus intercrystalline porosity. As stated above, the porosity loss due to compaction and increase by organic dissolution are estimated through the relative intensity of these diageneses in different reservoirs. During early diagenetic stage before intense chemical compaction, the total porosities may have been largely preserved in spite of the modification of pore systems [32].
The historical porosities are set in the time coordinate axis. Connect the dots representing the porosity valu d the porosity evolution curves are obtained. Typical reservoir sections of the Well PG2 and the Well MB4 are studied and a comparing curve of the Well YB2 is also shown in Figure 8.

Implication for Further Exploration
Zhao et al. [27], Yang et al. [33] and Ma et al. [34] have pointed out some key factors controlling the formation of the Feixianguan reservoirs including depositional facies, early meteoric dissolution, dolomitization, organic acids dissolution, burial dissolution and thermochemical sulfate reduction-related dissolution. And Xia et al. [30] provided some evidences against deep burial dissolution, especially the dissolution related to thermochemical sulfate reduction. Nevertheless, the main controlling fac-tor(s) is (are) still controversial, which greatly restricted further exploration of the Feixianguan reservoirs in other regions in the NE Sichuan Basin.
Through this research, paragenetic sequence and porosity evolution history of the Feixianguan reservoirs in the Puguang area are recovered, which clearly shows that the reservoir formation was mainly controlled by early meteoric dissolution, dolomitization, and organic acids dissolution other than deep burial dissolution. In general, the more pores were created before oil accumulation; the higher porosities are obtained at present. Porosities formed during relatively early diagenetic stages must be effectively preserved hence make great contribution to the present reservoir spaces. Consequently, the preservation mechanisms of the reservoir porosities may be of great significance to the formation of effective reservoirs in the NE Sichuan Basin. These preservation factors roughly proposed by t teoric cementation, dolomitization, thick evaporates as cap rocks leading to a closed diagenetic environments, early hydrocarbon emplacement and paleo-oil pool formation, overpressure through oil-gas cracking, etc.

Conclusions
Several conclusions can be drawn through this re 1) Diagenetic phases can be divided into near-surface (including marine and meteoric environment), pre-oil window, oil window, and gas window stages for the Feixianguan reservoirs, Puguang Gas Field, NE Sichuan Basin.
2) The main diagenetic processes affecting the Feixianguan reservoir units include: a) marine cementation; b) meteoric dissolution and cementation; c) early dolomitization and dolomite recrystallization; d) mechanical and chemical compaction; e) burial cementation; f) early fracturing; g) organic dissolution and hydrocarbon evolution; h) thermochemical sulfate reduction dissolution and cementation by pyrite, sulfur and calcite; and i) minor late fracturing. The porosity increase or loss by most of the diageneses can be quantified by the Adobe Photoshop Analysis method.
3) On the basis of the paragenetic sequence and the quantification of the influences of diageneses on reservoir porosity, porosity evolution curve could be quantitatively recovered.
4) The Feixianguan reservoirs were formed mainly by early meteoric dissolution, dolomitization, and organic acids dissolution other than deep burial dissolution. And the preservation mechanisms of the early-formed reservoir porosities may be of great significance to the further exploration of the deep buried Feixianguan reservoirs in the NE Sichuan Basin.