Total Organic Carbon Enrichment and Its Impact on Pore Characteristics : A Case Study from the Niutitang Formation Shales in Northern Guizhou

This study analyzes samples from the Lower Cambrian Niutitang Formation in northern Guizhou Province to enable a better understanding of total organic carbon (TOC) enrichment and its impact on the pore characteristics of over-mature marine shale. Organic geochemical analysis, X-ray diffraction, scanning electron microscopy, helium porosity, and low-temperature nitrogen adsorption experiments were conducted on shale samples. Their original TOC (TOCo) content and organic porosity were estimated by theoretical calculation, and fractal dimension D was computed with the fractal Frenkel–Halsey–Hill model. The results were then used to consider which factors control TOC enrichment and pore characteristics. The samples are shown to be dominated by type-I kerogen with a TOC content of 0.29–9.36% and an equivalent vitrinite reflectance value of 1.72–2.72%. The TOCo content varies between 0.64% and 18.17%, and the overall recovery coefficient for the Niutitang Formation was 2.16. Total porosity of the samples ranged between 0.36% and 6.93%. TOC content directly controls porosity when TOC content lies in the range 1.0% to 6.0%. For samples with TOC < 1.0% and TOC > 6.0%, inorganic pores are the main contributors to porosity. Additionally, pore structure parameters show no obvious trends with TOC, quartz, and clay mineral content. The fractal dimension D1 is between 2.619 and 2.716, and D2 is between 2.680 and 2.854, illustrating significant pore surface roughness and structural heterogeneity. No single constituent had a dominant effect on the fractal characteristics.


Introduction
Organic matter (OM) in shale refers to material abundant in organic carbon generated or retained in the shale through deep burial and increasing thermal maturity, including kerogen, bitumen, solid bitumen, residual OM, pyrobitumen, and char [1,2].OM receives considerable attention in shale gas exploration and exploitation because of its essential role in hydrocarbon generation and in organic pore formation and development [3,4].
Total organic carbon (TOC) content is a key parameter for characterizing a shale reservoir and is an indicator of shale gas potential.However, previous studies adopted various values for the lower limit for TOC in a viable resource.For example, the principal gas-productive section of the Ohio Shale had at least 2.0% TOC, which enabled the shale gas system to become the first source of commercial gas production in North America [5].The minimum TOC of an exploration target shale pay should be 2.5-3.0%[6].The source rock in the core producing area of the Barnett Shale system had a present-day The study area is in northern Guizhou, and well YX1 is located in northwestern Fenggang County (Figure 1).Tectonically, the northern Guizhou area evolved along with the Yangtze Platform, which experienced multiple periods of tectonic movement (from the Xuefeng movement to the Himalayan movement), resulting in complex geology [31].These tectonic movements resulted in northeast (NE)and north-northeast (NNE)-orientated folds and faults.The study area mainly features trough-like folds striking NE and NNE.The formation of S-shaped or Z-shaped single folds is attributed to tectonic deformation dominated by compression and facilitated by strike-slip.Their presence indicates that the geological conditions under which the area developed are complicated [31,32].Differently trending faults intersect each other.The dip angles of faults are relatively steep, including some well-preserved upright fault planes [33].The Lower Cambrian succession in the Upper Yangtze Platform was deposited during a major marine transgression from southeast to northwest.It is mainly composed of carbonaceous shale, and calcareous shale, carbonate, and siltstone.The black shale thickens southeastward; the present thick black shale in Sichuan, western Hunan and Hubei, Chongqing, and Guizhou province can be interpreted as being deposited in an open marine platform to marine shelf environment, while the deposits at the western margin of the Upper Yangtze Platform, which are dominated by terrestrial clastics and carbonates, are indicative of a littoral environment [28,35].
Occurrences of the Niutitang Formation in the study area can be subdivided into two members: (1) the lower member (LM), mainly consisting of carbonaceous shale, siliceous shale, and black shale with relatively high TOC content.It lies at depths of 1996-2031 m and is 35 m thick; (2) the upper member (UM), mainly consisting of organic-poor argillaceous mudstone and silty mudstone.It lies at depths of 1921-1996 m and is 75 m thick (Figure 2).The Lower Cambrian succession in the Upper Yangtze Platform was deposited during a major marine transgression from southeast to northwest.It is mainly composed of carbonaceous shale, and calcareous shale, carbonate, and siltstone.The black shale thickens southeastward; the present thick black shale in Sichuan, western Hunan and Hubei, Chongqing, and Guizhou province can be interpreted as being deposited in an open marine platform to marine shelf environment, while the deposits at the western margin of the Upper Yangtze Platform, which are dominated by terrestrial clastics and carbonates, are indicative of a littoral environment [28,35].
Occurrences of the Niutitang Formation in the study area can be subdivided into two members: (1) the lower member (LM), mainly consisting of carbonaceous shale, siliceous shale, and black shale with relatively high TOC content.It lies at depths of 1996-2031 m and is 35 m thick; (2) the upper member (UM), mainly consisting of organic-poor argillaceous mudstone and silty mudstone.It lies at depths of 1921-1996 m and is 75 m thick (Figure 2).

Samples and Experiments
Sixty-two samples were collected from the Lower Cambrian Niutitang Formation in YX1 well (Figure 1) with the depth from 1921 m to 2031 m (Figure 2), including 31 samples from the lower member (LM) and 31 samples from the upper member (UM).The core samples of the LM are dominated by black siliceous shale and carbonaceous shale, while the UM is mainly gray to dark-gray argillaceous shale and silty shale [34].Measurements including organic geochemistry, optical microscopy, X-ray diffraction, scanning electron microscopy (SEM), nitrogen adsorption, and helium pycnometry were performed on samples.The methods used in this study were described thoroughly in Reference [36]; thus, they are not fully presented in this paper.Table 1 lists experiments and measurements conducted on each sample.

Samples and Experiments
Sixty-two samples were collected from the Lower Cambrian Niutitang Formation in YX1 well (Figure 1) with the depth from 1921 m to 2031 m (Figure 2), including 31 samples from the lower member (LM) and 31 samples from the upper member (UM).The core samples of the LM are dominated by black siliceous shale and carbonaceous shale, while the UM is mainly gray to dark-gray argillaceous shale and silty shale [34].Measurements including organic geochemistry, optical microscopy, X-ray diffraction, scanning electron microscopy (SEM), nitrogen adsorption, and helium pycnometry were performed on samples.The methods used in this study were described thoroughly in Reference [36]; thus, they are not fully presented in this paper.Table 1 lists experiments and measurements conducted on each sample.

Calculation for Original TOC and Organic Porosity
Claypool's Equation is considered to be useful [3,7,11], but the original production index (PI o ) used for calculation was assumed by the researchers; thus, the transfer ratio (TR) of kerogen was unreliable.Both hydrogen index (HI) and original hydrogen index (HI o ) were used to compute the TR, which was more reliable [13].However, this calculation method ignores the mass change of source rock triggered by hydrocarbon generation and expulsion.Thus, rock mass ratio was taken into consideration [12], but the HI o determination was not applicable for over-mature shale.As a result, the first step of their method was revised, and the method is presented in detail below.
Firstly, the original hydrogen index (HI o ) was calculated by Equation (1).
where S, L, V, and I are volume percentages of sapropelinite, liptinite, vitrinite, and inertinite, respectively, on a mineral matter-free basis [11].
The TR refers to the degree of kerogen transformation due to hydrocarbon generation.Hydrogen index (HI) and original hydrogen index (HI o ) can be used to determine TR using Equation (2) [37].
The expulsion efficiency, f, refers to the ratio of the hydrocarbon expulsion amount to the total hydrocarbon generation amount, which can be expressed by Equation (3).
The definition of mass conversion factor w is the ratio of present-day rock mass to original rock mass [12]; the factor can be obtained from Equation (4).
Then, TOC o is determined using Equation (5).
where the scaling parameter α refers to the proportion of convertible carbon to total carbon, and α = HI o /1200.Finally, organic porosity Φ org is obtained from Equation (6).
where t stands for the carbon equivalent mass of kerogen, t = 1.2; and ρ b and ρ k are the bulk density of source rock and kerogen density, respectively.The type-I kerogen density is 1.45 g/cm 3 [38].

Fractal Theory
A fractal dimension (D) can characterize pore geometry by evaluating complexity and irregularity of both pore surface and structure [39][40][41].In this study, fractal dimensions can be determined using the Frenkel-Halsey-Hill (FHH) method, together with the N 2 adsorption isotherm data.The FHH model can be expressed by the following equation: where V is the volume of adsorbed gas molecules at the equilibrium pressure P, V 0 is the volume of monolayer coverage, A is the power-law exponent depending on D and the mechanism of adsorption, and P 0 is the saturation pressure of the gas.A is derived from the slope of linear regression equation in the plot of lnV versus ln(ln(P 0 /P)).The fractal dimension D derives from the two following expressions: D = A + 3 and D = 3A + 3.

Organic Geochemistry and Petrography
Figure 2 shows that the TOC contents of the Niutitang shale samples have a wide range of 0.29% to 9.36%; the upper and lower members average 0.69% and 5.93%, respectively.Only one sample from the upper member has a TOC value of over 2.0%.The bitumen reflectance (BR o ) value ranges from 1.47% to 2.56%, corresponding to vitrinite reflectance between 1.72% and 2.72% with a mean of 2.19% as calculated with the conversion equation [42].The majority (60-84%) of the organic matter is sapropelinite, while liptinite accounts for 15-37%.The amounts of vitrinite and inertinite are quite small, 0-3% and 0-2%, respectively.Table 2 presents the calculated type index (TI) and shows that only two samples have a TI lower than 80, indicating that the Niutitang shale is dominated by type-I kerogen.This finding is consistent with previous studies [28,34].Table 3 illustrates some of the main characteristics of the relative mineral composition obtained with X-ray diffraction (XRD) analysis.The samples mainly contain quartz, feldspar, carbonate, pyrite, and clay (illite, mixed-layer illite/smectite, and chlorite).The quartz content ranges from 20.3% to 54.5% with an average of 37.8%.The clay mineral content ranges from 8.7% to 47.1% and comprises illite (50-100%), mixed-layer illite/smectite (0-29%), and chlorite (0-31%).In eight of the 11 samples from the lower member of the Niutitang Formation, the only clay mineral is illite, revealing that the lower member entered the post-mature stage, with vitrinite reflectance >2.0% [43].The other samples contain illite, chlorite, and mixed-layer illite/smectite, in which 10% comprises smectite layers, indicating a highly mature stage [43,44].The relationships between thermal maturity and clay mineral type and content ensure that the equivalent vitrinite reflectance (EqVR o ) values calculated by the fitting equation [42] are relatively reasonable.Furthermore, the pyrite content lies in the range of 1.3% to 13.8%; the mean pyrite contents of the upper and the lower member are 3.3% and 10.2%, respectively.Generally, pyrite forms in an anoxic environment [15]; thus, it can be taken as evidence that the lower member of Niutitang Formation had a stagnant, anoxic depositional environment [34].

Estimation of Original TOC and Organic Porosity
Table 2 presents the results of original TOC reconstruction and organic porosity estimation.The HI o ranges from 640.75 mg/gTOC to 689.5 mg/gTOC with a mean of 660.75 mg/gTOC.The TR values are over 95%, indicating that the organic matter present is approaching the end of hydrocarbon generation [3].The low S2 values suggest that the source rock has almost no hydrocarbon generation potential (Table 4).Additionally, the high expulsion efficiency reveals that there is a minute quantity of remaining hydrocarbon.The TOC o values for the upper member are between 0.64% and 1.91% with a mean of 1.15%, and those for the lower member range from 10.66% to 18.17% and average 14.02%, significantly higher than the upper member.S2 was between 4.14 mg/g (rock) and 125.28 mg/g (rock).The parameters in Table 5 demonstrate that a considerable amount of hydrocarbon was generated and expelled during thermal maturation.Furthermore, the recovery coefficient (the ratio of TOC o to TOC pd ) is between 2.02 and 2.28 with an average of 2.16.TOC o has a good positive correlation with TOC pd , indicating that all of the samples experienced similar changes in TOC o .The organic porosity derived from TOC o is shown in Table 5.Values for the upper member range from 0.85% to 2.67% and average 1.55%, whereas those for the lower member are between 12.84% and 22.59% with a mean of 16.94%.The two members had different capacities for generating organic pores due to the marked discrepancy in TOC o between them.However, the calculated organic porosity is merely a theoretical value.

Total Porosity
The helium porosity of the Niutitang shale in the study area ranges from 0.36% to 6.93% with a mean of 2.61%.That of the upper member lies in the range of 0.36% to 2.49% with a mean of 1.01%, while that of the lower member varies from 2.26% to 6.93% with a mean of 4.36%.Figure 3a shows how porosity changes with TOC.When the TOC is lower than 1.0%, the points are scattered; when the TOC is between 1.0% and 6.0%, porosity is positively correlated with TOC; when the TOC content exceeds 6.0%, the positive trend between these two parameters becomes negative, which is consistent with the findings of some previous studies [20,24,25,45].However, it is noteworthy that sample 47 had a TOC content of 6.86% and a porosity of 6.93%, which differs from the negative trend illustrated by Figure 3a.A possible explanation for this phenomenon is that sample 47 has a relatively high brittle mineral content, as evidenced by the well-developed fractures seen in its SEM image (Figure 3b); fractures are less developed in the other samples.

Total Porosity
The helium porosity of the Niutitang shale in the study area ranges from 0.36% to 6.93% with a mean of 2.61%.That of the upper member lies in the range of 0.36% to 2.49% with a mean of 1.01%, while that of the lower member varies from 2.26% to 6.93% with a mean of 4.36%.Figure 3a shows how porosity changes with TOC.When the TOC is lower than 1.0%, the points are scattered; when the TOC is between 1.0% and 6.0%, porosity is positively correlated with TOC; when the TOC content exceeds 6.0%, the positive trend between these two parameters becomes negative, which is consistent with the findings of some previous studies [20,24,25,45].However, it is noteworthy that sample 47 had a TOC content of 6.86% and a porosity of 6.93%, which differs from the negative trend illustrated by Figure 3a.A possible explanation for this phenomenon is that sample 47 has a relatively high brittle mineral content, as evidenced by the well-developed fractures seen in its SEM image (Figure 3b); fractures are less developed in the other samples.In Figure 3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

N2 Adsorption-Desorption Isotherm
The nitrogen adsorption-desorption isotherms from this study are presented in Figure 4.The isotherms of all the Niutitang Formation shale samples can be classed as type IV in the International Union of Pure and Applied Chemistry (IUPAC) classification system [46,47].These isotherms have no plateau within the high relative pressure interval and are characterized by pronounced hysteresis loops due to capillary condensation within mesopores [46].The adsorption volume is low and rises slowly in the low relative pressure range, but the isotherms became steep in the high relative pressure interval (p/p0 > 0.9).Although the degree of development of hysteresis loops varies, samples 7, 15, and 33 can be considered type H3 and the rest of samples type H4 (Figure 4) according to the IUPAC classification In Figure 3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

N 2 Adsorption-Desorption Isotherm
The nitrogen adsorption-desorption isotherms from this study are presented in Figure 4.The isotherms of all the Niutitang Formation shale samples can be classed as type IV in the International Union of Pure and Applied Chemistry (IUPAC) classification system [46,47].These isotherms have no plateau within the high relative pressure interval and are characterized by pronounced hysteresis loops due to capillary condensation within mesopores [46].The adsorption volume is low and rises slowly in the low relative pressure range, but the isotherms became steep in the high relative pressure interval (p/p 0 > 0.9).

Total Porosity
The helium porosity of the Niutitang shale in the study area ranges from 0.36% to 6.93% with a mean of 2.61%.That of the upper member lies in the range of 0.36% to 2.49% with a mean of 1.01%, while that of the lower member varies from 2.26% to 6.93% with a mean of 4.36%.Figure 3a shows how porosity changes with TOC.When the TOC is lower than 1.0%, the points are scattered; when the TOC is between 1.0% and 6.0%, porosity is positively correlated with TOC; when the TOC content exceeds 6.0%, the positive trend between these two parameters becomes negative, which is consistent with the findings of some previous studies [20,24,25,45].However, it is noteworthy that sample 47 had a TOC content of 6.86% and a porosity of 6.93%, which differs from the negative trend illustrated by Figure 3a.A possible explanation for this phenomenon is that sample 47 has a relatively high brittle mineral content, as evidenced by the well-developed fractures seen in its SEM image (Figure 3b); fractures are less developed in the other samples.In Figure 3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

N2 Adsorption-Desorption Isotherm
The nitrogen adsorption-desorption isotherms from this study are presented in Figure 4.The isotherms of all the Niutitang Formation shale samples can be classed as type IV in the International Union of Pure and Applied Chemistry (IUPAC) classification system [46,47].These isotherms have no plateau within the high relative pressure interval and are characterized by pronounced hysteresis loops due to capillary condensation within mesopores [46].The adsorption volume is low and rises slowly in the low relative pressure range, but the isotherms became steep in the high relative pressure interval (p/p0 > 0.9).Although the degree of development of hysteresis loops varies, samples 7, 15, and 33 can be considered type H3 and the rest of samples type H4 (Figure 4) according to the IUPAC classification Although the degree of development of hysteresis loops varies, samples 7, 15, and 33 can be considered type H3 and the rest of samples type H4 (Figure 4) according to the IUPAC classification system, suggesting the existence of slit-or plate-like pores [46].The hysteresis loops also indicate that ink-bottle pores are rare in the samples.Cylindrical pores, slit-like pores, and wedge-like pores are open pores that are conducive to gas flow, while ink-bottle pores can adsorb gas [47].Consequently, the sorption capacity of these samples may be low.

Pore Structure
Values for pore structure parameters derived from the nitrogen adsorption isotherms are shown in Table 6.The Brunauer-Emmett-Teller (BET) specific surface area ranges from 5.08 m 2 /g to 25.31 m 2 /g with a mean of 16.21 m 2 /g.The micropore surface area computed by the t-plot method ranges from 0.96 m 2 /g to 7.78 m 2 /g with an average of 4.85 m 2 /g.The calculated Barrett-Joyner-Halenda (BJH) total pore volume varies between 0.013 cm 3 /g and 0.029 cm 3 /g and averages 0.019 cm 3 /g, while the micropore volume lies in the range 0.0004-0.0037cm 3 /g with a mean of 0.0024 cm 3 /g.The average pore diameter varies between 6.18 nm and 12.62 nm with an average of 7.45 nm, indicating that the Niutitang Formation shale is dominated by mesopores.Figure 5 shows all of the samples to have similar pore size distributions except for sample 7, which shows no significant peak and, thus, has a relatively uniform pore distribution.As a result, this sample exhibits the smallest BET surface area and pore volume and the largest average pore size.Other samples tend to be unimodal with a major peak in the pore size of 2-3 nm, revealing that they are dominated by mesopores.Additionally, a small but non-negligible amount of macropores are present.
Energies 2019, 12, 1480 11 of 23 system, suggesting the existence of slit-or plate-like pores [46].The hysteresis loops also indicate that ink-bottle pores are rare in the samples.Cylindrical pores, slit-like pores, and wedge-like pores are open pores that are conducive to gas flow, while ink-bottle pores can adsorb gas [47].Consequently, the sorption capacity of these samples may be low.

Pore Structure
Values for pore structure parameters derived from the nitrogen adsorption isotherms are shown in Table 6.The Brunauer-Emmett-Teller (BET) specific surface area ranges from 5.08 m 2 /g to 25.31 m 2 /g with a mean of 16.21 m 2 /g.The micropore surface area computed by the t-plot method ranges from 0.96 m 2 /g to 7.78 m 2 /g with an average of 4.85 m 2 /g.The calculated Barrett-Joyner-Halenda (BJH) total pore volume varies between 0.013 cm 3 /g and 0.029 cm 3 /g and averages 0.019 cm 3 /g, while the micropore volume lies in the range 0.0004-0.0037cm 3 /g with a mean of 0.0024 cm 3 /g.The average pore diameter varies between 6.18 nm and 12.62 nm with an average of 7.45 nm, indicating that the Niutitang Formation shale is dominated by mesopores.Figure 5 shows all of the samples to have similar pore size distributions except for sample 7, which shows no significant peak and, thus, has a relatively uniform pore distribution.As a result, this sample exhibits the smallest BET surface area and pore volume and the largest average pore size.Other samples tend to be unimodal with a major peak in the pore size of 2-3 nm, revealing that they are dominated by mesopores.Additionally, a small but non-negligible amount of macropores are present.

Fractal Characteristics
The FHH plots for the shale samples are presented in Figure 6.The presence of two distinct linear segments in these plots at P/P0 ranges of 0-0.5 and 0.5-1 suggests that the fractal characteristics vary within those two relative pressure ranges; thus, the fractal dimensions D1 and D2 were acquired from the two linear intervals, respectively.Table 7 gives the slopes of the linear regression equations A1 and A2 and the fractal dimension values D1 and D2.By using the equation D = 3A + 3, the D1 values of samples 7 and 52 were calculated to be less than 2, which is

Fractal Characteristics
The FHH plots for the shale samples are presented in Figure 6.The presence of two distinct linear segments in these plots at P/P 0 ranges of 0-0.5 and 0.5-1 suggests that the fractal characteristics vary within those two relative pressure ranges; thus, the fractal dimensions D1 and D2 were acquired from the two linear intervals, respectively.Table 7 gives the slopes of the linear regression equations A 1 and A 2 and the fractal dimension values D1 and D2.By using the equation D = 3A + 3, the D1 values of samples 7 and 52 were calculated to be less than 2, which is unreasonable because the fractal dimension lies in the range 2 to 3 with increasing irregularity and complexity [39].As a result, fractal dimensions calculated with the equation D = A + 3 are used for further analysis and discussion.
Energies 2019, 12, 1480 12 of 23 unreasonable because the fractal dimension lies in the range 2 to 3 with increasing irregularity and complexity [39].As a result, fractal dimensions calculated with the equation D = A + 3 are used for further analysis and discussion.The D1 values range from 2.619 to 2.716 with an average value of 2.684, and the D2 values range from 2.680 to 2.854 with a mean of 2.813.The results illustrate that both the pore surfaces and structures of the shales are highly irregular and heterogeneous.Moreover, a good positive correlation between D1 and D2 can be observed in Figure 7, revealing that D1 and D2 were both suitable for characterizing pore structures in the Niutitang Formation shale.This finding is consistent with a previous study focusing on the Lower Silurian Longmaxi Formation shale [27].The higher D2 values suggest that, with rising pore diameter, the pore surface of the shale becomes rougher and the pore structure more complicated.The D1 values range from 2.619 to 2.716 with an average value of 2.684, and the D2 values range from 2.680 to 2.854 with a mean of 2.813.The results illustrate that both the pore surfaces and structures of the shales are highly irregular and heterogeneous.Moreover, a good positive correlation between D1 and D2 can be observed in Figure 7, revealing that D1 and D2 were both suitable for characterizing pore structures in the Niutitang Formation shale.This finding is consistent with a previous study focusing on the Lower Silurian Longmaxi Formation shale [27].The higher D2 values suggest that, with rising pore diameter, the pore surface of the shale becomes rougher and the pore structure more complicated.unreasonable because the fractal dimension lies in the range 2 to 3 with increasing irregularity and complexity [39].As a result, fractal dimensions calculated with the equation D = A + 3 are used for further analysis and discussion.The D1 values range from 2.619 to 2.716 with an average value of 2.684, and the D2 values range from 2.680 to 2.854 with a mean of 2.813.The results illustrate that both the pore surfaces and structures of the shales are highly irregular and heterogeneous.Moreover, a good positive correlation between D1 and D2 can be observed in Figure 7, revealing that D1 and D2 were both suitable for characterizing pore structures in the Niutitang Formation shale.This finding is consistent with a previous study focusing on the Lower Silurian Longmaxi Formation shale [27].The higher D2 values suggest that, with rising pore diameter, the pore surface of the shale becomes rougher and the pore structure more complicated.

TOC Recovery
Hydrocarbon generation is a result of organic matter decomposition with increasing thermal maturity [6]; thus, original TOC decreases gradually due to the loss of the convertible part of the organic matter.As more than three-quarters of the samples entered the over-mature stage, the HI values of the LM samples dropped to below 1 mg/g (Table 4), illustrating that a large amount of hydrocarbon was generated and expelled.
In addition, 13 samples have a recovery coefficient between 2.02 and 2.28, showing little variation.However, the present-day TOC of the LM is markedly higher than that of the UM.Four LM samples with high present-day TOC had much higher TOC o , while nine UM samples with low present-day TOC (<1.0%) had relatively low TOC o .This indicates that TOC enrichment was dependent on the sedimentary environment.Previous studies characterized paleo-environments on the basis of productivity proxies and redox condition proxies [34,48,49].For example, upwelling and hydrothermal events triggered a remarkable increase in paleo-productivity during LM deposition and showed that the redox conditions evolved from anoxic with euxinic intervals during LM deposition to oxic conditions with anoxic intervals upward [34].In this study, mineral compositions are used to characterize the sedimentary environment.Firstly, quartz content is positively correlated with TOC in the LM samples but has a weak negative correlation with TOC in the UM samples (Figure 8).Quartz of biogenic origin has a positive correlation with TOC, whereas the correlation between terrigenous detrital origin quartz and TOC is negative [50,51].Thus, the quartz source of the LM was dominantly authigenic, while that in the UM had both terrigenous inputs and an authigenic source.Secondly, the average pyrite contents of the UM and the LM shale are 3.4% and 10.3%, respectively.These results also indicate that the LM was deposited in a reducing environment with relatively stagnant water mass and that this transformed into an oxidizing environment upward.Hydrocarbon generation is a result of organic matter decomposition with increasing thermal maturity [6]; thus, original TOC decreases gradually due to the loss of the convertible part of the organic matter.As more than three-quarters of the samples entered the over-mature stage, the HI values of the LM samples dropped to below 1 mg/g (Table 4), illustrating that a large amount of hydrocarbon was generated and expelled.
In addition, 13 samples have a recovery coefficient between 2.02 and 2.28, showing little variation.However, the present-day TOC of the LM is markedly higher than that of the UM.Four LM samples with high present-day TOC had much higher TOCo, while nine UM samples with low present-day TOC (<1.0%) had relatively low TOCo.This indicates that TOC enrichment was dependent on the sedimentary environment.Previous studies characterized paleo-environments on the basis of productivity proxies and redox condition proxies [34,48,49].For example, upwelling and hydrothermal events triggered a remarkable increase in paleo-productivity during LM deposition and showed that the redox conditions evolved from anoxic with euxinic intervals during LM deposition to oxic conditions with anoxic intervals upward [34].In this study, mineral compositions are used to characterize the sedimentary environment.Firstly, quartz content is positively correlated with TOC in the LM samples but has a weak negative correlation with TOC in the UM samples (Figure 8).Quartz of biogenic origin has a positive correlation with TOC, whereas the correlation between terrigenous detrital origin quartz and TOC is negative [50,51].Thus, the quartz source of the LM was dominantly authigenic, while that in the UM had both terrigenous inputs and an authigenic source.Secondly, the average pyrite contents of the UM and the LM shale are 3.4% and 10.3%, respectively.These results also indicate that the LM was deposited in a reducing environment with relatively stagnant water mass and that this transformed into an oxidizing environment upward.

Porosity and Controlling Factors
Organic matter decomposes during thermal maturation, leading to the formation of hydrocarbons and intraparticle organic pores [4], and migrated organic matter develops secondary pores as it cracks at high maturity [52].Nanoporosity evolutionary sequence includes the formation stage (0.60% < VRo < 2.0%), the development stage (2.0%< VRo < 3.5%), and the conversion or destruction stage (VRo > 3.5%) [53].However, no general trend between organic porosity and maturity was observed when VRo > 0.9% in the Woodford shale [14].In this study, organoporosity was calculated with the method presented in Reference [12].There is no obvious relationship between the calculated organoporosity and EqVRo (Figure 9).Because of the marked differences in the TOC of the two members, samples from the LM had higher organoporosity (over 10%) than those from the UM.Thus, both total porosity and organic porosity had no direct relationship with thermal maturity, and other factors like TOC content should be given more consideration when characterizing shale porosity.

Porosity and Controlling Factors
Organic matter decomposes during thermal maturation, leading to the formation of hydrocarbons and intraparticle organic pores [4], and migrated organic matter develops secondary pores as it cracks at high maturity [52].Nanoporosity evolutionary sequence includes the formation stage (0.60% < VR o < 2.0%), the development stage (2.0%< VR o < 3.5%), and the conversion or destruction stage (VR o > 3.5%) [53].However, no general trend between organic porosity and maturity was observed when VR o > 0.9% in the Woodford shale [14].In this study, organoporosity was calculated with the method presented in Reference [12].There is no obvious relationship between the calculated organoporosity and EqVR o (Figure 9).Because of the marked differences in the TOC of the two members, samples from the LM had higher organoporosity (over 10%) than those from the UM.Thus, both total porosity and organic porosity had no direct relationship with thermal maturity, and other factors like TOC content should be given more consideration when characterizing shale porosity.Figure 3a reveals the relationship between total porosity and TOC.Total porosity is positively correlated with TOC when 1.0% < TOC < 6.0%, after which the trend changes.A similar pattern was reported in other studies; for example, there was a change in a positive correlation at TOC = 5.5% in the Marcellus shale from North America [24], TOC = 5.0% in the Qiongzhusi Formation from southern Sichuan Basin [45], and TOC = 12% in the Permian shale from eastern China [25].
The contribution of per-gram organic carbon content to total porosity can be derived from the slope of the linear regression equation in a plot of total porosity versus TOC.For samples with 1.0% < TOC < 6.0%, the TOC contribution is 0.93% (Figure 3a), which is greater than in Lower Silurian shale in southern China (0.72%) [23] and higher than in Lower Cambrian shale in the southern Sichuan Basin (0.40%) [45].However, these values may be invalid because they are based on limited samples; further study is needed to determine the contribution made by TOC content.For samples with TOC < 1.0% and TOC > 6.0%, organic pores make less of a contribution to porosity.Figure 10 shows that the recovery coefficient and calculated organoporosity have no obvious relationship with porosity, indicating that inorganic pores are the main contributors to porosity.When TOC < 1.0%, the volume of organic pores is too small to influence the porosity.The porosity decreasing when TOC > 6.0% may be attributable to mechanical compaction [24,25].A previous study showed that Young's modulus values of kerogen, clay, and quartz were 10 GPa, 15-45 GPa, and 50-70 GPa, respectively [54], and an association between significant lowering of the kerogen modulus and organic pore development was illustrated.The immature kerogen modulus was 15-20 GPa, whereas kerogen in the gas window had a smaller modulus of 7-12 GPa [54].This is one reason why fewer organic pores can be observed in the SEM images (Figure 11).Another possible explanation is that the organic matter seen in the SEM images is depositional organic matter, which mainly consists of kerogen and seldom yields pores [52].Lastly, some micropores and fine mesopores cannot be clearly observed; thus, the insufficient magnification of the SEM images may be a factor [7].  Figure 3a reveals the relationship between total porosity and TOC.Total porosity is positively correlated with TOC when 1.0% < TOC < 6.0%, after which the trend changes.A similar pattern was reported in other studies; for example, there was a change in a positive correlation at TOC = 5.5% in the Marcellus shale from North America [24], TOC = 5.0% in the Qiongzhusi Formation from southern Sichuan Basin [45], and TOC = 12% in the Permian shale from eastern China [25].
The contribution of per-gram organic carbon content to total porosity can be derived from the slope of the linear regression equation in a plot of total porosity versus TOC.For samples with 1.0% < TOC < 6.0%, the TOC contribution is 0.93% (Figure 3a), which is greater than in Lower Silurian shale in southern China (0.72%) [23] and higher than in Lower Cambrian shale in the southern Sichuan Basin (0.40%) [45].However, these values may be invalid because they are based on limited samples; further study is needed to determine the contribution made by TOC content.For samples with TOC < 1.0% and TOC > 6.0%, organic pores make less of a contribution to porosity.Figure 10 shows that the recovery coefficient and calculated organoporosity have no obvious relationship with porosity, indicating that inorganic pores are the main contributors to porosity.When TOC < 1.0%, the volume of organic pores is too small to influence the porosity.The porosity decreasing when TOC > 6.0% may be attributable to mechanical compaction [24,25].A previous study showed that Young's modulus values of kerogen, clay, and quartz were 10 GPa, 15-45 GPa, and 50-70 GPa, respectively [54], and an association between significant lowering of the kerogen modulus and organic pore development was illustrated.The immature kerogen modulus was 15-20 GPa, whereas kerogen in the gas window had a smaller modulus of 7-12 GPa [54].This is one reason why fewer organic pores can be observed in the SEM images (Figure 11).Another possible explanation is that the organic matter seen in the SEM images is depositional organic matter, which mainly consists of kerogen and seldom yields pores [52].Lastly, some micropores and fine mesopores cannot be clearly observed; thus, the insufficient magnification of the SEM images may be a factor [7].  Figure 3a reveals the relationship between total porosity and TOC.Total porosity is positively correlated with TOC when 1.0% < TOC < 6.0%, after which the trend changes.A similar pattern was reported in other studies; for example, there was a change in a positive correlation at TOC = 5.5% in the Marcellus shale from North America [24], TOC = 5.0% in the Qiongzhusi Formation from southern Sichuan Basin [45], and TOC = 12% in the Permian shale from eastern China [25].
The contribution of per-gram organic carbon content to total porosity can be derived from the slope of the linear regression equation in a plot of total porosity versus TOC.For samples with 1.0% < TOC < 6.0%, the TOC contribution is 0.93% (Figure 3a), which is greater than in Lower Silurian shale in southern China (0.72%) [23] and higher than in Lower Cambrian shale in the southern Sichuan Basin (0.40%) [45].However, these values may be invalid because they are based on limited samples; further study is needed to determine the contribution made by TOC content.For samples with TOC < 1.0% and TOC > 6.0%, organic pores make less of a contribution to porosity.Figure 10 shows that the recovery coefficient and calculated organoporosity have no obvious relationship with porosity, indicating that inorganic pores are the main contributors to porosity.When TOC < 1.0%, the volume of organic pores is too small to influence the porosity.The porosity decreasing when TOC > 6.0% may be attributable to mechanical compaction [24,25].A previous study showed that Young's modulus values of kerogen, clay, and quartz were 10 GPa, 15-45 GPa, and 50-70 GPa, respectively [54], and an association between significant lowering of the kerogen modulus and organic pore development was illustrated.The immature kerogen modulus was 15-20 GPa, whereas kerogen in the gas window had a smaller modulus of 7-12 GPa [54].This is one reason why fewer organic pores can be observed in the SEM images (Figure 11).Another possible explanation is that the organic matter seen in the SEM images is depositional organic matter, which mainly consists of kerogen and seldom yields pores [52].Lastly, some micropores and fine mesopores cannot be clearly observed; thus, the insufficient magnification of the SEM images may be a factor [7].The relationships between quartz, clay minerals, and total porosity are presented in Figure 12.Total porosity has a slight negative correlation with quartz content in the lower member but no obvious relationship with it in the Upper member.Similar findings were obtained from analysis of Devonian-Mississippian shale in western Canada, showing that total porosity had a negative correlation with biogenic quartz due to the loss of pores with a diameter of less than 10 nm [50].Furthermore, quartz of detrital origin with angular grain boundaries is more likely to form interparticle pores, which are mainly macropores and mesopores, leading to a positive correlation between detrital quartz and total porosity [23,31].
The total porosity has a good negative correlation with clay content.Previous studies proposed that clay minerals increase total porosity [51,55].Additionally, the lattice of randomly oriented clay mineral platelets forms a "house of cards" structure [14], contributing to the preservation of linear and triangular pores.However, for shale layers, higher clay content may weaken their capacity to resist compaction, which will cause a pronounced decrease in porosity.Clay-associated pores form with increasing clay content, but more pores, including organic pores and mineral pores, vanish during the mechanical compaction resulting from high clay content.As a result, the contribution of clay-associated pores to total porosity may be negligible [25].Hypotheses regarding this phenomenon are preliminary and need further confirmation.The relationships between quartz, clay minerals, and total porosity are presented in Figure 12.Total porosity has a slight negative correlation with quartz content in the lower member but no obvious relationship with it in the Upper member.Similar findings were obtained from analysis of Devonian-Mississippian shale in western Canada, showing that total porosity had a negative correlation with biogenic quartz due to the loss of pores with a diameter of less than 10 nm [50].Furthermore, quartz of detrital origin with angular grain boundaries is more likely to form interparticle pores, which are mainly macropores and mesopores, leading to a positive correlation between detrital quartz and total porosity [23,31].
The total porosity has a good negative correlation with clay content.Previous studies proposed that clay minerals increase total porosity [51,55].Additionally, the lattice of randomly oriented clay mineral platelets forms a "house of cards" structure [14], contributing to the preservation of linear and triangular pores.However, for shale layers, higher clay content may weaken their capacity to resist compaction, which will cause a pronounced decrease in porosity.Clay-associated pores form with increasing clay content, but more pores, including organic pores and mineral pores, vanish during the mechanical compaction resulting from high clay content.As a result, the contribution of clay-associated pores to total porosity may be negligible [25].Hypotheses regarding this phenomenon are preliminary and need further confirmation.The relationships between quartz, clay minerals, and total porosity are presented in Figure 12.Total porosity has a slight negative correlation with quartz content in the lower member but no obvious relationship with it in the Upper member.Similar findings were obtained from analysis of Devonian-Mississippian shale in western Canada, showing that total porosity had a negative correlation with biogenic quartz due to the loss of pores with a diameter of less than 10 nm [50].Furthermore, quartz of detrital origin with angular grain boundaries is more likely to form interparticle pores, which are mainly macropores and mesopores, leading to a positive correlation between detrital quartz and total porosity [23,31].
The total porosity has a good negative correlation with clay content.Previous studies proposed that clay minerals increase total porosity [51,55].Additionally, the lattice of randomly oriented clay mineral platelets forms a "house of cards" structure [14], contributing to the preservation of linear and triangular pores.However, for shale layers, higher clay content may weaken their capacity to resist compaction, which will cause a pronounced decrease in porosity.Clay-associated pores form with increasing clay content, but more pores, including organic pores and mineral pores, vanish during the mechanical compaction resulting from high clay content.As a result, the contribution of clay-associated pores to total porosity may be negligible [25].Hypotheses regarding this phenomenon are preliminary and need further confirmation.Figure 13 illustrates the relationships among the pore structure parameters of the shale samples.Specific surface area is significantly correlated with total pore volume (R 2 = 0.81) and micropore volume (R 2 = 0.96), which agrees with previous studies on highly mature marine shale [27,40,51,56].Additionally, the average pore size exhibits a negative correlation with total pore volume and specific Energies 2019, 12, 1480 16 of 23 surface area (Figure 13a,b).Samples with smaller average pore diameters tend to possess more micropores and mesopores, leading to a higher pore volume and specific surface area.Micropores may be a more significant contributor to the specific surface area than are mesopores.
The relationships between shale constituents and pore structure parameters are shown in Figure 13.The slightly positive relationships between TOC content and total surface area, micropore surface area, total pore volume, and micropore volume seen in Figure 13c,d suggest that TOC is not the major contributor to these structural parameters in the samples.However, previous studies drew different conclusions, finding that specific surface area and pore volume have a significantly positive correlation with TOC content [27,40,57].An examination of the SEM images of the samples used in those studies reveals that this correlation may be attributable to the presence of a large amount of well-preserved organic pores.Organic pores are limited in the samples considered in the current study, as seen in the SEM images, yielding a different outcome.
The moderate negative trends between quartz content and total surface area and total pore volume shown in Figure 13e,f reveal that quartz content does not exert a significant influence on pore structure in this study.Biogenetic quartz has smooth rims, which tend to form less interparticle porosity.In addition, quartz tends to provide less pore volume [58].
Specific surface area and total pore volume have a slightly negative relationship with clay content (Figure 13g,h), which is in agreement with some previous studies [25,26,58].On the one hand, the specific surface area values of clay minerals vary considerably.For example, montmorillonite and mixed-layer illite/smectite have relatively high specific surface areas of 76.4 m 2 /g and 30.8 m 2 /g, respectively, but chlorite and illite have smaller surface areas of 11.7 m 2 /g and 7.1 m 2 /g [59].As presented in Section 4.1, the clay content of these samples is dominated by illite (average 78%), which makes less of a contribution to the surface area.On the other hand, parts of clay-associated pores that are not protected from brittle minerals will be deformed by mechanical compaction.Thus, the specific surface area and total pore volume show little correlation with clay mineral content.To conclude, no simple relationship is seen between pore structure and shale composition in this study.
volume and specific surface area (Figure 13a,b).Samples with smaller average pore diameters tend to possess more micropores and mesopores, leading to a higher pore volume and specific surface area.Micropores may be a more significant contributor to the specific surface area than are mesopores.
The relationships between shale constituents and pore structure parameters are shown in Figure 13.The slightly positive relationships between TOC content and total surface area, micropore surface area, total pore volume, and micropore volume seen in Figures 13c,d suggest that TOC is not the major contributor to these structural parameters in the samples.However, previous studies drew different conclusions, finding that specific surface area and pore volume have a significantly positive correlation with TOC content [27,40,57].An examination of the SEM images of the samples used in those studies reveals that this correlation may be attributable to the presence of a large amount of well-preserved organic pores.Organic pores are limited in the samples considered in the current study, as seen in the SEM images, yielding a different outcome.
The moderate negative trends between quartz content and total surface area and total pore volume shown in Figure 13e,f reveal that quartz content does not exert a significant influence on pore structure in this study.Biogenetic quartz has smooth rims, which tend to form less interparticle porosity.In addition, quartz tends to provide less pore volume [58].
Specific surface area and total pore volume have a slightly negative relationship with clay content (Figure 13g,h), which is in agreement with some previous studies [25,26,58].On the one hand, the specific surface area values of clay minerals vary considerably.For example, montmorillonite and mixed-layer illite/smectite have relatively high specific surface areas of 76.4 m 2 /g and 30.8 m 2 /g, respectively, but chlorite and illite have smaller surface areas of 11.7 m 2 /g and 7.1 m 2 /g [59].As presented in Section 4.1, the clay content of these samples is dominated by illite (average 78%), which makes less of a contribution to the surface area.On the other hand, parts of clay-associated pores that are not protected from brittle minerals will be deformed by mechanical compaction.Thus, the specific surface area and total pore volume show little correlation with clay mineral content.To conclude, no simple relationship is seen between pore structure and shale composition in this study.

Characteristics and Controlling Factors of Fractal Dimension
Figure 14 shows the relationship between the fractal dimensions D1 and D2 and the constituents of the shale.Neither D1 nor D2 shows an obvious correlation with quartz content.This is because the quartz particles have smooth surfaces, resulting in a relatively smooth and regular overall surface and microstructure (Figure 14a).D1 and D2 exhibit no clear correlation with clay content (Figure 14b).The layered and flocculent structures of clay minerals increase the specific surface area of pores [58]; however, in this study, the samples have high illite content due to high thermal maturity, and the contribution is not significant.The lack of a relationship between mineral composition and fractal dimension indicates that mineral composition is not the major factor influencing fractal dimension.

Characteristics and Controlling Factors of Fractal Dimension
Figure 14 shows the relationship between the fractal dimensions D1 and D2 and the constituents of the shale.Neither D1 nor D2 shows an obvious correlation with quartz content.This is because the quartz particles have smooth surfaces, resulting in a relatively smooth and regular overall surface and microstructure (Figure 14a).D1 and D2 exhibit no clear correlation with clay content (Figure 14b).The layered and flocculent structures of clay minerals increase the specific surface area of pores [58]; however, in this study, the samples have high illite content due to high thermal maturity, and the contribution is not significant.The lack of a relationship between mineral composition and fractal dimension indicates that mineral composition is not the major factor influencing fractal dimension.
Figure 14c shows a slight increase in D1 with rising TOC content and a moderately positive correlation between D2 and TOC content.Organic-rich shale tends to develop more nanopores due to kerogen decomposition during thermal maturation [3,52].These nanopores form a complicated pore network, resulting in higher fractal dimensions [16,26].However, limited organic pore development occurred in these samples due to their kerogen type; thus, the positive correlation between fractal dimension and TOC content is not significant in this study.In sum, the sedimentary environment and thermal maturation control the shale constituents, and no single constituent dominates the fractal characteristics.Figure 14c shows a slight increase in D1 with rising TOC content and a moderately positive correlation between D2 and TOC content.Organic-rich shale tends to develop more nanopores due to kerogen decomposition during thermal maturation [3,52].These nanopores form a complicated pore network, resulting in higher fractal dimensions [16,26].However, limited organic pore development occurred in these samples due to their kerogen type; thus, the positive correlation between fractal dimension and TOC content is not significant in this study.In sum, the sedimentary environment and thermal maturation control the shale constituents, and no single constituent dominates the fractal characteristics.The relationship between the fractal dimensions and the pore structure parameters is also shown in Figure 14.D1 and D2 have positive linear correlations with surface area, and D1 has a better correlation with both total surface area (Figure 14d) and micropore surface area than does D2 (Figure 14e).This is because D1 represents the surface fractal dimension, which is dependent on the pore surface area [39].Figure 14f exhibits a good linear correlation between fractal dimension and average pore size.As previous studies illustrated, a smaller pore size will not only increase the complexity and heterogeneity of the pore structure, but will also generate a rougher surface [16,27].
Moreover, the fractal dimension shows no obvious correlation with total pore volume (Figure 14g) but a better correlation with micropore volume (Figure 14h).On the one hand, the specific surface area is well correlated with pore volume; thus, samples with a higher fractal dimension have a larger pore volume.On the other hand, micropores are the major contributors to pore structure complexity and surface roughness.Thus, the correlation of the fractal dimension with micropores will be better than with other pores.Furthermore, the observation that both D1 and D2 have a pronounced correlation with the pore structure parameters reveals that fractal dimensions can be applied to characterize pore structure.

Conclusions
In this paper, the TOC distribution, pore structure, and evolutionary characteristics of Niutitang Formation shales from the northern Guizhou Province in southwest China were investigated using a series of experiments, mass balance calculations, and FHH theory.The following conclusions were drawn: (1) The organic matter was dominated by type-I kerogen with a total organic carbon (TOC) content of 0.29-9.36%,averaging 0.69% and 5.93% in the upper and the lower members, respectively.The equivalent vitrinite reflectance (EqVR o ) values ranged from 1.72% to 2.72% with a mean of 2.19%.The original TOC content values of the upper member were between 0.64% to 1.91% with a mean of 1.15%, while those of the lower member ranged from 10.66% to 18.17% and averaged 14.02%.The overall TOC recovery coefficient for the Niutitang formation was 2.16, indicating good hydrocarbon generation potential and that it is a promising exploration prospect.
The sedimentary environment was shown to be the dominant factor in both present-day TOC and the original TOC distribution.(2) The total porosity of the Niutitang shale ranged from 0.36% to 6.93% with a mean of 2.61%.
This ranged from 0.36% to 2.49% with a mean of 1.01% in samples from the upper member and 2.26% to 6.93% with a mean of 4.36% in samples from the lower member.When TOC was below 1.0%, porosity had no correlation with it; when TOC was between 1.0% and 6.0%, porosity was positively correlated with it; when TOC < 1.0% and TOC > 6.0%, porosity had no obvious correlation with it.This indicates that, when within a certain range, TOC was a significant controlling factor in porosity.The calculated organic porosity was higher than the measured total porosity, and the relationship between porosity and mineral composition was found to be indirect.(3) The nitrogen adsorption isotherms of all of the Niutitang Formation shale samples were type IV, with hysteresis loops identified as type H3 and type H4, suggesting the existence of slit-or plate-like pores.The BET specific surface area ranged from 5.08 m 2 /g to 25.31 m 2 /g, total pore volume from 0.013 cm 3 /g to 0.029 cm 3 /g, and average pore diameter from 6.18 nm to 12.62 nm.Pore structure parameters showed no obvious correlation with the composition of the shale.(4) The D1 values ranged from 2.619 to 2.716, with an average value of 2.684, and the D2 values ranged from 2.680 to 2.854, with a mean of 2.813.The irregularity of both the pore surface and structure indicated good gas sorption capacity.Sedimentary environment and thermal maturation were shown to control the constituents of the shale, and no single constituent dominated its fractal characteristics.

Figure 1 .
Figure 1.(a) Location of Guizhou Province; (b) location of Northern Guizhou area; (c) location of the study area and well location (Modified from Reference [34]).

Figure 1 .
Figure 1.(a) Location of Guizhou Province; (b) location of Northern Guizhou area; (c) location of the study area and well location (Modified from Reference [34]).

Figure 2 .
Figure 2. Lithological profile and total organic carbon (TOC) content of well YX1 section.

Figure 2 .
Figure 2. Lithological profile and total organic carbon (TOC) content of well YX1 section.

Figure 3 .
Figure 3. Relationship between total porosity and TOC content (a) and SEM image of sample 47 (b).In Figure3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

Figure 3 .
Figure 3. Relationship between total porosity and TOC content (a) and SEM image of sample 47 (b).In Figure3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

Figure 3 .
Figure 3. Relationship between total porosity and TOC content (a) and SEM image of sample 47 (b).In Figure3a, red dots represent samples with TOC values over 6.0%, green dots represent samples with TOC values between 1.0% and 6.0%, and blue dots represent samples with TOC values below 1.0%.

Figure 6 .
Figure 6.Plots of lnV vs. ln(ln(P 0 /P)) obtained from the N 2 adsorption isotherms of two typical shale samples: (a) sample 7 and (b) sample 15.In Figure 6a,b, red dots indicate the P/P 0 range of 0.5-1, and blue dots indicate the P/P 0 range of 0-0.5.

Figure 8 .
Figure 8. Relationship between TOC and quartz content.

Figure 8 .
Figure 8. Relationship between TOC and quartz content.

5. 3 .Figure 11 .
Figure13illustrates the relationships among the pore structure parameters of the shale samples.Specific surface area is significantly correlated with total pore volume (R 2 = 0.81) and micropore volume (R 2 = 0.96), which agrees with previous studies on highly mature marine shale[27,40,51,56].Additionally, the average pore size exhibits a negative correlation with total pore

5. 3 .Figure 12 .
Figure13illustrates the relationships among the pore structure parameters of the shale samples.Specific surface area is significantly correlated with total pore volume (R 2 = 0.81) and micropore volume (R 2 = 0.96), which agrees with previous studies on highly mature marine shale[27,40,51,56].Additionally, the average pore size exhibits a negative correlation with total pore

Figure 13 .
Figure 13.Relationship among structure parameters (a,b) and effects of shale composition on structure parameters (c-h); in (b) and (e-h), specific surface area (SSA) refers to total SSA, pore volume (PV) refers to total PV.

Figure 13 .
Figure 13.Relationship among structure parameters (a,b) and effects of shale composition on structure parameters (c-h); in (b) and (e-h), specific surface area (SSA) refers to total SSA, pore volume (PV) refers to total PV.

Figure 14 .
Figure 14.Relationship between fractal dimensions and shale composition (a-c) and relationship between fractal dimensions and structure parameters (d-h).

Figure 14 .
Figure 14.Relationship between fractal dimensions and shale composition (a-c) and relationship between fractal dimensions and structure parameters (d-h).

Table 1 .
Experiments and measurements conducted on samples.

Table 1 .
Experiments and measurements conducted on samples.

Table 3 .
Mineralogical composition of the shale samples.

Table 4 .
Rock-Eval data and HI of samples.
Note: T max : maturity parameter based on the temperature at which the maximum amount of pyrolyzate (S2) is generated from the kerogen in a rock sample; S1: free hydrocarbons present in the rock (mg HC/g of rock); S2: remaining generation potential (mg HC/g of rock); S3: oxidizable carbon, (mg CO 2 /g rock); HI: hydrogen index.

Table 5 .
TOC o calculation, total porosity, and organic porosity estimation.

Table 7 .
Fractal dimensions of pores obtained from N2 adsorption of the shale samples.

Table 7 .
Fractal dimensions of pores obtained from N 2 adsorption of the shale samples.

Table 7 .
Fractal dimensions of pores obtained from N2 adsorption of the shale samples.