The Role of Microfabric and Laminae on Pore Structure and Gas Transport Pathways of Marine Shales from Sichuan Basin, China

This study investigated the effects of microfabric and laminae on the pore structure and gas transport pathways of the Silurian Longmaxi shales from Sichuan Basin. 23 shale samples with varied lithofacies were comprehensively investigated by mineralogy, organic geochemistry, pycnometry, and low-pressure nitrogen adsorption analysis. The fabric and laminae of these samples were identified using petrographic microscope and scanning electron microscopy. Permeabilities were measured using the nonsteadystate method on both perpendicular and parallel to bedding shales. The effective pore diameter controlling gas transport was estimated from gas slippage factors obtained in permeability measurements. These values were also compared to those calculated using the Winland equation. Siliceous shales studied are faintly laminated to nonlaminated and have larger porosity and specific surface area. Argillaceous/siliceous mixed shales are well laminated, whereas argillaceous shales contain many oriented clay flakes along the lamination. Both porosity and surface area are positively correlated with TOC content. Unlike most conventional reservoirs, there is a negative correlation between porosity and permeability values of the samples studied. Permeabilities parallel to bedding, ranging from 0.4 to 76.6 μD, are in control of the oriented clay flakes and silty microlaminae. Permeability anisotropy values of the shales vary between 1.3 and 49.8. Samples rich in oriented clay flakes and microlaminated fabric have relatively larger permeability and permeability anisotropy values. The effective transport pore diameters derived from gas slippage measurements are slightly lower than those calculated from the Winland equation. However, both methods have shown that the effective transport pore diameters of argillaceous shales (averaging 552 nm) are significantly higher than siliceous shales (averaging 198 nm), which underlines the control of microfabric, rather than porosity, on gas transport pathways of the shales studied.


Introduction
Shale systems have received increased research activities since hydrocarbons have been commercially extracted from unconventional reservoirs. Shale reservoirs consist of various fine-grained sedimentary rocks that mainly include mudrocks, shales, and siltstones [1][2][3]. Due to the large variability in lithology and involvement of organic matter, reservoir characteristics of shales are very complex. To characterize the complex pore structure, innovative approaches are implemented to obtain information about the nanometer-to micrometer-sized pore systems of shales [4][5][6][7][8]. A variety of organic matter-hosted and inorganic matter-related pores with pore size ranging from micropore (<2 nm, IUPAC classification) to mesopore (2-50 nm, IUPAC classification) to macropore (>50 nm, IUPAC classification) have been reported in shales [9][10][11]. Pore structure characteristics (morphology, porosity, pore size distribution, and specific surface area) of shales are found to be associated with rock composition, total organic carbon (TOC) content, maturity, and other geological factors [12][13][14][15][16].
Permeability is one of the most important petrophysical parameters for reservoir evaluation and hydrocarbon production [17,18]. Since shale reservoirs commonly have permeabilities down to micro-to nano-Darcy range, steady and nonsteady-state flow tests are employed to determine the permeability values of shales [19]. Permeabilities are affected by many geological factors, including fabric, lithology, mineralogy, pore size distribution (PSD), and microfractures. Studies on high-maturity Horn River shales have shown that the matrix permeabilities were largely related to the connectivity between macropores and micropores [20]. Permeability values of coarser grained sandstones are commonly higher than those of finer grained mudstones by several orders of magnitude [21]. Permeabilities are further affected by diagenetic process. Mechanical compaction and mineral cementation can significantly decrease pore space and fluid transport pathways of rocks. Clay minerals are reported to cement with the framework minerals and thus decrease permeability of Eagle Ford shales [22]. It should be noted that most of the samples in these studies are lowmaturity or organic-lean, and most porosity is associated with the interparticle pores of inorganic matter. Recently, advanced scanning electron microscopy has documented diverse pores in shale systems [6,16]. However, there are relatively few detailed studies to relate pore structure to the transport properties of shales. Especially, organic matter-(OM-) hosted pores are critical for gas adsorption in many overmature marine shales [23][24][25], but the role of OM-hosted pores on fluid storage and transport needs more study in detail. Furthermore, permeability can be estimated from more easily determined petrophysical properties, such as mean pore diameter and pore throat size distribution [26,27]. However, many pore space analysis methods are conducted at atmospheric pressure. It would be much more meaningful to apply these pore structure parameters in the estimations of permeability under in situ conditions [19,28]. What is more, the effective pore diameter for gas transport can be further derived from gas slippage in gas permeability measurements [29]. These link a direct relationship between pore structure and gas transport.
Concerns about the economic potential of shale systems have led to many studies about reservoir properties of shales [12,30,31]. Significant progresses have been achieved in basin structure, sedimentology environment, organic geochemistry, and petrophysical features of shales [17,[32][33][34][35][36]. Shales are commonly composed of finegrained minerals intercalated with coarser-grained quartz and carbonate microlaminae. High-frequency laminations and strong heterogeneities probably lead to different permeability values and significant permeability anisotropy even on the same formation [37,38]. Gas transport in shales is controlled not only by geometry (pore size distribution) but also topology of pores (connectivity, preferential pathways, etc.) [39]. Although some studies have documented the heterogeneities and the pore systems associated with facies in gas shales [30,40,41], the effects of microfabric and laminae on gas flow are not yet described in detail.
We report on a suite of shale samples that represent various properties from the Jiaoshiba area, Sichuan Basin, China. This study examines the role of microfabric and laminae on pore space and transport properties of shales. With a combination of imaging observation (petrographic microscope and scanning electron microscope) and fluid invasion measurements (helium expansion porosimetry, low-pressure N 2 adsorption-desorption, and pulse-decay N 2 tests), petrophysical properties of Longmaxi shales are comprehensively characterized. Pore structures of various shales are qualitatively described. The effects of minerals and microfabric on porosity, PSD, and pore structure of shales are discussed. Furthermore, both the Winland equation and gas permeability measurements are applied to calculate the effective pore diameter. Hereby, controls of microfabric and laminae on gas transport on shales are assessed.

Geological Background
Jiaoshiba shale gas field, the first commercially developed shale reservoirs in China, is located in the Fuling District, Chongqing Municipality, in the eastern Sichuan Basin ( Figure 1). Gas production in the Jiaoshiba shale reservoir is reported to be up to 50 × 10 4 m 3 /day per well [42]. The Jiaoshiba structure is a wide and gentle anticline, constrained by two sets of northeast-trending faults [43][44][45]. As a part of Sichuan Basin, the Jiaoshiba structure was subjected to multistage tectonic movements during the formation of Yangtze platform [43,46]. In Late Ordovician-Early Silurian time, the Sichuan Basin transformed from a passive continental margin basin to a foreland basin. Because of the uplift in the southern, central, and eastern boundaries, the Sichuan Basin formed a shelf structure of being open on the northern side. Correspondingly, the sedimentary environment was undercompensated and anoxic in the study area [47][48][49]. There were two global transgressions at the end of Ordovician period and the beginning of Silurian period [50,51], which caused the deposition of Ordovician Wufeng-Silurian Longmaxi (O 3 w-S 1 l) shale system. After that, the sea level gradually declined. A mass of terrigenous detrital materials were involved in the middle and late deposition of S 1 l Formation [49].
From the lower to upper sections of O 3 w-S 1 l shales, the sedimentary facies transformed gradually from deep-water shelf to shallow-water shelf [52]. The O 3 w-S 1 l shales can be approximately divided into three members based on lithological features. The lower member of the Longmaxi shales (including the O 3 w formation), deposited in a deep-water shelf sedimentary environment, is dominated by dark gray and black siliceous shales. A lot of benthic algae, radiolarian, graptolite, and radiolarian fossils are documented in the lower section of Longmaxi shales previously [52], which provide organic matter for siliceous shales [16]. At the same time, the sedimentary environment was anoxic and euxinic when siliceous shales were deposited, which is beneficial to the preservation of organic matter in siliceous shales [47]. The middle part of Longmaxi Formation is basin-slope facies. Compared to the early deposition of Longmaxi shales, the supply of terrigenous clasts at the middle Longmaxi Formation was significantly increased [46,52]. Consequently, lithological features of the shales gradually converted from siliceous into silt-bearing shales. Silty laminae are developed well in the middle part of Longmaxi shales. The upper section was formed in shallow-water shelf environments [48]. Shales  Emeis han-L iangs han Fold Belt H u a y in g s h a n A n t ic li n e L o n g q u a n s h a n A n t i c l i n e L o n g m e n s h a n F o l d B e l t W e s t d i s t r i c t  in the upper section are dominated by light gray mudstones and siltstones [52]. The variability of lithology from the lower to the upper member of the Longmaxi shales in Jiaoshiba area offers a suitable case for assessing the relationship between shale lithofacies and reservoir characteristics. The Longmaxi Formation in the Jiaoshiba area varies in thickness from 60 to 150 m. Longmaxi shales with the vitrinite reflectance equivalent values being higher than 2.0% are highly overmatured [25,47].

Experimental Section
3.1. Sample Preparation. Marine shale samples were collected from the lower to upper section of Longmaxi Formation. Twenty-three samples with depth ranging from 2527.8 to 2615.3 m were obtained from a shale gas well drilled in the Jiaoshiba shale gas field (Table 1). Samples were selected to cover the variability of total organic carbon (TOC), inorganic minerals composition, and shale lithofacies. 21 cylindrical plugs with 25-40 mm in length and 25.4 mm (1 inch) in diameter were drilled parallel to the bedding, and 8 plugs were drilled perpendicular to bedding. Cylindrical plugs were dried in a vacuum oven (105°C) for at least 12 h before porosity and permeability measurements. Petrologic and scanning electron microscope analyses were carried out on subsamples (fragments). Mineralogical components, TOC content, and low-pressure nitrogen adsorption-desorption analyses were performed on dry particles (0.15 mm particle size).

Imaging Methods.
Petrologic observations were conducted on polished thin sections to analyze the rock fabric, texture, microorganism, and mineralogy. Polished thin sections with thickness between 0.02 and 0.06 mm were prepared and photographed using a Leica DMRX microscope. Visual pore structure was observed on a Zeiss Merlin Compact field emission-scanning electron microscope (FE-SEM). This apparatus is equipped with backscattered electron (BSE) and secondary electron (SE) detectors and achieves the highest pixel resolution of about 0.8 nm at 15 kV operating voltage. Shale samples were ion polished and coated with gold before the visual observations.
3.3. Low-Pressure Low-Temperature Nitrogen (N2) Adsorption. Low-pressure N 2 adsorption-desorption experiments at -196°C (77 K) were performed to obtain information about the specific surface area, pore size distribution (PSD), and pore volume. Powder shales were first degassed  4 Geofluids overnight at 105°C in a vacuum chamber. N 2 adsorption and desorption isotherms were obtained by collecting adsorption data on degassed samples at a pressure up to 730 mm Hg. The nitrogen adsorption data were interpreted to obtain specific surface area and PSD according to Brunauer-Emmett-Teller (BET) and Barrett-Joyner-Halenda (BJH) theories [53].

Helium Pycnometry and Gas
Permeability. Helium pycnometry was applied on dried cylindrical plugs to determine the skeletal density (ρ skeletal ). The bulk density (ρ bulk ) of these plugs was measured from the geometric dimensioning of cylindrical plugs using a caliper and the weight of samples. Then, porosity values (ϕ) under unconfined conditions were determined from the bulk density and skeletal density values according to ϕ = 1 − ρ bulk /ρ skeletal . Gas permeability coefficients were measured on dried cylindrical samples using N 2 as measuring gas at room temperature. Sample plugs were placed in a customer-designed pulse-decay permeameter (PDP-200, Core Lab) at a selected confining pressure of 2000 PSI. Small pressure pulses of about 0.5-1.0 MPa were introduced to the upstream reservoir, and the pressures of the upstream and downstream reservoirs were monitored with high-precision pressure transducers. Different average pore pressures are obtained via introducing several small-pressure pulses. The nonsteady-state nitrogen permeability at different pore pressures was determined using a linear regression performed on the pressure time data according to the fundamental flow equations [37]. The measured nitrogen permeability was corrected using the Klinkenberg method based on the permeability values at different pore pressures [29].

Mineralogy and TOC Content.
Mineralogy compositions of the studied Longmaxi samples are illustrated in Figure 2. The Jiaoshiba shale samples mainly consist of quartz (30.5-58.4%) and clay minerals (21.9-62.7%) with a minor amount of plagioclase (3.5-9.9%), dolomite (0-32.2%), calcite (0-8.4%), pyrite (0-6.6%), and potassium feldspar (0-3.0%). These agree with the results of other studies, showing that quartz and clays are the most common minerals in the Jiaoshiba shales [47]. In general, quartz content in our samples decreases from the bottom to top Longmaxi Formation 5 Geofluids whereas clay minerals present a reverse tend (Table 1). TOC content in our data set ranges from 0.3 to 5.3 wt.% with a mean value of 2.2 wt.% (Table 1). TOC content is relatively higher in the lower member of Longmaxi Formation. There is a positive correlation between quartz and TOC content ( Figure 3). The positive correlation is due to that the quartz in the lower part of Longmaxi Formation is mainly of biogenic origin [20].

Types of Shales.
Lithofacies record geological information about the rock types, texture, bedding, mineral compositions, grain size distribution, and sorting and roundness characteristics [40,41]. Initially, lithofacies were mainly applied to analyze the mineralogy and petrology of conventional reservoirs [54]. Since the shale gas development is active, lithofacies were introduced to unconventional reservoirs [24,41]. Shale lithofacies comprehensively describe reservoir characteristics and quality [41]. Based on the mineralogy composition, shales can be divided into siliceous shales (siliceous minerals content > 50%), argillaceous/siliceous mixed shales (siliceous, carbonate, and clay minerals content is less than 50%, respectively), and argillaceous shales (total clays content > 50%) ( Figure 2 and Table 1).

Siliceous Shales.
Siliceous shales (S) are mainly found in the lower member of Longmaxi Formation (Table 1). Siliceous shale samples are rich in quartz and relatively low in clay minerals. In siliceous shales, the content of siliceous minerals ranges from 50.3% to 66.2% with a mean value of 56.5%. The total clay mineral content ranges from 25.6% to 38.9% with an average of 31.9%. Illite-smectite (I/S) mixed clay minerals represent the major components of clay minerals (range from 10.5% to 19.2% with an average of 14.5%), followed by illite (range from 8.3% to 16.1% with an average of 11.7%) and chlorite (averages 4.1%). Siliceous shale samples are commonly rich in pyrite. Pyrite in siliceous shales exhibits an average of 4.8% (ranges from 3.8% to 6.6%).
Thin sections of siliceous shales show that the fabric of siliceous shales ranges from faintly laminated to nonlaminated with increasing depth (Figure 4). Abundant radiolarians are observed in siliceous shales (Figure 4(d)). Some radiolarians are replaced by pyrite. Most of the quartz particles in siliceous shales are silt-sized scale. Clays are dyed by organic matter and scattered in the minerals. Only minor lenticular clay flakes are observed along the lamination (Figure 4(c)). Siliceous shales are TOC-rich. Siliceous shales have TOC content ranging from 2.5% to 5.3% (average 3.4%).   Figure 5). The laminae mainly consist of silty quartz. The width of laminae ranges from 0.01 to 0.40 mm. Clays in argillaceous/siliceous mixed shales are commonly dyed by organic matter. TOC content of the mixed Longmaxi shales is between 1.4% and 2.3% with a mean value of 1.8%.  (Table 1). This kind of shales is characterized by the highest content of total clays and the lowest content of quarts. The total clay mineral content in CM shales studied varies from 54.6% to 62.7% (averages 58.0%), whereas the siliceous mineral content ranges between 37.3% and 42.7% (averages 39.8%). Chlorite is the highest content of clay minerals in CM shales. Chlorite exhibits a mean of 22.4% (ranges from 18.9% to 26.3%), and I/S mixed clay minerals average 18.6% (ranges from 15.8% to 21.2%). Slight of pyrite is found in CM shales (averages 1.6%).
Thin sections of CM shales show nonlaminae to intermittent laminae ( Figure 6). In argillaceous shales, many long striped and large lenticular clay flakes are oriented along the lamination. Minor organic flakes are also found in CM shales ( Figure 6(a)). Quartz is either scattered in the minerals (Figures 6(a) and 6(b)) or exhibits intermittent laminae (Figures 6(c) and 6(d)). TOC content of CM shales is really low and ranges from 0.3% to 1.7% (averages 1.0%).   In mixed shales (M), OM particles are dispersed and commonly captured in the pore spaces among rigid framework minerals (Figure 8). A handful of OM-hosted pores can be found in these residual organic matters (Figures 8(c)  and 8(d)). OM particles with large volume occur in argillaceous/siliceous mixed shales, but pores are less developed in these bulky OM particles (Figure 8(b)). Apparently, these bulky OM particles are difficult to be thermally decomposed or migrated. Dissolution pores are commonly found in argillaceous/siliceous mixed shales (Figures 8(c) and 8(d)). These dissolution pores appear to be so scattered that they are probably less connected. Dissolution-related pores/rims develop among the boundaries of feldspar grains (Figure 8(c)).
When it comes to CM shales, the amount of silty with large particles becomes less than that in mixed shales (M). Clay minerals are really common in the field of vision (Figure 9). Most clay flakes are oriented along the lamination (Figure 9(a)). There are many intraparticle pores in clay flakes, and these pores can be enhanced by bending/splitting related to compaction (Figure 9(c)). A spot of organic matter is captured in the corners of framework minerals (Figure 9(c)). However, these residual organic matters are either nonporous or less porous. Kerogen-clay aggregates are commonly observed and parallel with lamination ( Figure 9(d)). However, many pore spaces in kerogen-clay aggregates are unobservable due to strong postcompaction. Large OM particles with smooth surface can be observed occasionally, but OM pores are unobservable in these OM particles under the resolution of FE-SEM (Figure 9(b)).

4.4.
Porosity, Pore Size Distribution, and Pore Structure Parameters. Porosity of the Longmaxi shales studied ranges from 1.5% to 7.1% and averages 3.7% ( Table 2). The porosity values of shale samples are associated with shale lithofacies. Siliceous shales generally show higher porosity values than other shales (Figure 3). Porosity of siliceous shales is between 3.4% and 7.1% with an average of 4.9%. Porosity of mixed shales ranges between 1.9% and 3.3% (averages 2.4%), whereas the porosity of argillaceous shales ranges from 1.5% to 4.3% (averages 3.1%). Positive relationships exist among TOC content, siliceous mineral content, and porosity values of the Longmaxi shales studied (Figure 10).
Low-pressure N 2 adsorption-desorption isotherms for the Longmaxi Formation samples are shown in Figure 11. Samples rich in organic matter generally adsorb the high amount of nitrogen. What is more, siliceous shale samples obviously adsorb a higher amount of N 2 than other samples (Figure 11(a)). All these Longmaxi shales exhibit hysteresis loops, formed by the divergence between adsorption and (a) According to the physical adsorption theory, the shape of hysteresis loop is associated with the complexity of pore structure [53]. Large hysteresis loops qualitatively indicate that pores of the materials are mainly inkbottle-shaped or amorphous. The pore structure of siliceous shales are much more complicated than other samples, which are probably related with the high TOC content of siliceous shales [16]. The adsorption branches of nitrogen adsorption isotherms were used to calculate the PSD. All the samples exhibit broad PSDs ranging from micropores to mesopores and macropores ( Figure 11). PSD of these Longmaxi shales are dominantly bimodal. There is one maximum in the microporous ranges of PSDs. The other maximum is in the 2-3 nm on the PSDs. Comparing the peaks of different shales, it is found that the first peak in the micropore range of siliceous shales (S) is obviously higher than others, whereas the second peak of mixed (M) and argillaceous shales (CM) is slightly higher than siliceous shales (S).
Pore structure parameters measured by N 2 adsorption on the Longmaxi shales are summarized in Table 2. BET-specific surface areas of the Longmaxi shales vary from 9.7 to 26.9 m 2 /g, with a mean value 17.0 m 2 /g. The total pore volumes vary from 0.011 to 0.033 cm 3 /g. Both specific surface areas (SSA) and pore volumes of siliceous Longmaxi shales are higher than those of mixed (M) and argillaceous (CM) Longmaxi shales. The total pore volume of S, M, and CM shales averages as 0.023, 0.017, 0.015 cm 3 /g, respectively. TOC content shows obviously positive linear relationships with both the BET SSA and pore volume (Figures 12(a) and 12(b)), indicating that the pore structure of organic-rich marine shales is controlled by organic matter. The intercept of these linear fitting is approximately 9 m 2 /g and 0.01 cm 3 /g (Figures 12(a) and 12(b)), which is related to the inorganic matter. A moderate positive correlation exists between siliceous mineral content and pore structure parameters (Figures 12(c) and 12(d)). Negative relationship exists between clay mineral content and pore structure parameters of the Longmaxi shales (Figures 12(e) and 12(f)).

Permeability and Permeability
Anisotropy. Gas permeability values of overall 29 core plugs were determined. Among them, twenty-one samples were drilled parallel to bedding and eight plugs were perpendicular to bedding. Parallel to bedding, the Klinkenberg-corrected nitrogen permeability (intrinsic permeability) values measured using the pulse-decay permeameter on the confining pressure of 2000 psi ranging from 0.4 to 76.6 μD (10 -18 m 2 ) ( Table 2). The Klinkenberg-corrected permeability values measured parallel to bedding (K h ) of siliceous shales are generally lower than those of other shales. The average K h of the Longmaxi shales studied decrease in a descending order: argillaceous  (Figures 4, 5, and 6). The K h of argillaceous shales are generally more than one order of magnitude higher than that measured perpendicular to bedding (Table 2). Furthermore, samples rich in laminated fabric always show high permeability anisotropy ( Figure 4).
Gas slippage factors for nitrogen on these shales range from 0.04 to 0.88 MPa. Gas slippage factors of siliceous shales are larger than argillaceous shales. The slippage factors of siliceous minerals range from 0.14 to 0.26 MPa, while those of argillaceous shales are between 0.07 to 0.14 MPa.

The Effect of Minerals and Microfabric on Pore Structure.
Mineralogical controls on pore spaces have been extensively investigated on both marine and lacustrine shales [12,16]. Generally, there is a positive correlation between porosity and TOC content of many mature and overmature shale systems [16]. Similarly, porosity of the Longmaxi shales in our data set is positively related to TOC content. The average porosity related to inorganic matter is estimated as 1.7% according to the interception of the linear relationship ( Figure 10(a)). The positive correlation between TOC content and porosity of Longmaxi shales studied underlines the key role of organic matter on pore spaces of overmatured marine shales. With increasing TOC content, peaks in the microporous range in PSD of Longmaxi shales are higher (Figure 11(b)). TOC content also shows an obviously positive relationship with BET SSA of Longmaxi shales studied (Figure 12(a)).
Further investigations reveal that porosity of shale samples is associated with the types of shales. Porosity of siliceous shales (averages 4.9%) is commonly larger than mixed shales (averages 2.4%) and argillaceous shales (averages 3.1%) ( Table 2). High porosity values of siliceous shales (S) are ascribed to high TOC content (Figure 7). However, inorganic matter contributes a lot to the pore space of mixed shales (M) and argillaceous shales (CM). For argillaceous/siliceous mixed shales, dissolution-related porosity is commonly observed in feldspars and carbonates. Dissolutionrelated pores occur as intraparticle pores in feldspars and carbonate grains or interparticle pores at the rims of carbonates ( Figure 8). This kind of pore mainly forms by the dissolution of organic acidic during the process of hydrocarbon generation [15,16]. Moreover, silty quartz and carbonates occur in laminae as clasts ( Figure 5) and serve as rigid skeleton to support matrix. Interparticle pores exist between soft organic matter and rigid clasts (Figure 8). These kinds 10 Geofluids of interparticle pores occur as preferable transport pathways for bitumen and gas. Argillaceous shales have specific microfabrics. In argillaceous shales, clay minerals are deformed (Figure 9(c)). Intraparticle pores are situated along cleavage planes of clay flakes and enhanced by bending/splitting related to compaction. Interparticle pores also exist between soft clay flakes and rigid clasts (Figure 9(c)). The inorganic-porosity seems to be reflected in the PSD. The peaks in 2-3 nm range of mixed (M) and argillaceous shales (CM) are higher than those of siliceous shales (Figures 11(b), 11(d), and 11(f)), which is likely contributed by clay minerals [55].
Linear combination is a useful tool to quantitatively assess the petrophysical properties of rocks [12,25]. Total porosity of shales can be viewed as the sum of porosity values related to OM, framework minerals, and total clay minerals: Here, w OM , w Framework , and w Clay are the mass fractions of OM, framework minerals (quartz, carbonates, pyrite, etc.), and total clays, respectively; v OM , v Framework , and v Clay The average pore diameter was determined by low-pressure N 2 adsorption analysis. b The transport pore diameter was derived from gas slippage measurement on parallel to bedding samples. c The dominant pore-throat aperture ðr p35 Þ was estimated using the Winland equation [59]. 11 Geofluids (cm 3 /g) are the pore volume of organic matter, framework minerals, and total clays, normalized to the mass of rock, respectively.
The bulk density (ρ bulk ) and porosity (ϕ) values can be determined from helium pycnometry measurements, while the mass fractions of individual mineral components can be obtained from XRD analysis. Thus, there are three unknown parameters (specific pore volumes v OM , v Framework , and v Clay ) in the physical model, which can be simultaneously solved using a multielement nonlinear regression on a suit of shale samples. The fitted v OM , v Framework , and v Clay of the Longmaxi shales is 0.4534 cm 3 /g, 0.0002 cm 3 /g, and 0.0124 cm 3 /g, respectively. Consequently, the organic matter-hosted pores account for 79% of the total porosity while clay minerals pores account for 20% of total porosity in siliceous shales ( Figure 13). In mixed shales, organic matter-hosted pores account for 62% while the pore spaces associated with clay minerals account for 38% of the total porosity. Argillaceous shales are dominant by pores associated with clay minerals. In argillaceous shales, clay minerals account for 62% of the total porosity while the OM-hosted pores account for 37%. It should be noted that the linear combination method pro-vides fast and semiquantitative characterization of pore spaces. The effects of mineralogy on porosity and pore space of rocks are much more complicated and associated to the interaction between organic and inorganic matters. The contribution of organic matter to porosity of argillaceous shales is probably an optimistic value since it is difficult to find lots of OM-hosted pores in argillaceous shales (Figure 9).

The Effect of Microfabric and Laminae on Permeability
and Permeability Anisotropy. Shale permeability is related to many geological controls: including fabric, texture, lithology, porosity, and pore size distribution [20,38]. Geological controls on permeable property of tight shales are always rock specific. In Devonian gas shales from the Horn River and Liard Basins of Canada, high matrix permeability shales were found to contain connected transport pathways between macropores and micropores [20]. Gas shales containing balanced ratio of micro-, meso-and macropores have higher matrix permeability values. Yang and Aplin [21] analyzed the effect of lithology on permeability of deeply buried mudstones and found that coarser-grained mudstones always have higher permeability values than finer-grained The reason for that is that the clay flakes and the microlaminae of clasts oriented to the lamination/bedding noticeably influence the permeabilities (K h ) of these Longmaxi shales 13 Geofluids studied ( Figures 5 and 6). Intraparticle pores of oriented clay minerals and interparticle pores between rigid clast particles and soft minerals (clay minerals or OM) are commonly distributed along the lamination/bedding of argillaceous shales and mixed shales (Figures 8 and 9). Though the inorganic matter-related porosity is less than OM-hosted porosity, the intraparticle pores in clay minerals and interparticle pores between organic and inorganic minerals play a significant role (preferable pathways) in fluid transport processes. These can be illustrated by the residual bitumen (migrabitumen) captured in these pore spaces among grains (Figures 8(c) and 8(d) and Figure 9(c)). Conversely, siliceous shales at the bottom of Longmaxi Formation are mainly formed in euxinic sediment environment [48,49]. Siliceous shales in Longmaxi Formation are obviously less in laminae compared to mixed and argillaceous shales (Figure 4). In siliceous shales, a lot of spongy organic matter are distributed without any preferred direction in the vision of SEM (Figure 7). What is more, permeabilities (K h ) of siliceous shales with faint laminae (JYA-17, JYA-18) are also higher than nonlaminated 14 Geofluids siliceous shales (JYA-22, Figure 4). In a word, TOC content does influence the porosity of these Longmaxi shales but is not the critical factor on permeability.
The porosity-permeability relationships of shales have been assessed in previous investigations. Some investigators reported a power law relationship between porosity and permeability [21]. In this study, there is a weak negative relationship between porosity and permeability coefficients ( Figure 14). The fabric and texture highly influence the relationships between porosity and permeability. Though the porosity of siliceous shales is larger than other shales, permeability values of this kind of shales are actually lower. The intraparticle pores of oriented clay minerals and interparticle pores between organic and inorganic matters (microlaminae of clasts) develop in M and CM shales ( Figures 5 and 6), which contribute to both K h and K v .
Permeability anisotropy (the ratio of horizontal permeability K h to vertical permeability K v ) has been reported on shale systems [37,56]. Permeability anisotropy of Longmaxi shales in this study varies from 1.3 to 49.8 (confining pressure of 2000 psi), which is comparable to other shales [37]. Permeability anisotropy is related to the pores with oriented alignment, permeable layers, and preferable migration pathway of shales. In this study, permeability anisotropy of argillaceous shales is higher than that of siliceous shales ( Figure 15). The permeability anisotropy values of three siliceous shale samples are less than ten, whereas K h of argillaceous shales are generally higher than K v more than one order of Porosity associated with organic matter Porosity associated with clay Porosity associated with framework minerals Figure 13: The semiquantitative estimation of porosity related to mineralogy using the linear combination approach. magnitude. This is probably related to the microfabric of the shales. Abundant clay flakes are oriented to the lamination in argillaceous shales (Figures 6(a) and 6(c)), whereas quartz in siliceous shales is randomly orientated and clays are scattered (Figures 4(a) and 4(d)). Furthermore, samples rich in laminated fabric (sample) always show high permeability anisotropy (Figure 4(b) and Figure 5(b)).

Effective Transport
Pore Diameter of Shales. The effective transport pore diameter can be derived from gas slippage factors during permeability measurements [19,57]. By assuming cylindrical pores, the effective transport pore diameter is a function of gas slippage factors using the theory of Klinkenberg [29]: where d ðmÞ is effective pore diameter for gas transport;,c ≈ 1 is the Adzumi constant, λ ðmÞ is the gas mean free path, P m ðPaÞ is mean pore pressure, and b ðPaÞ is the gas slippage factor. Combined with the slippage factors, the effective transport pore diameter for gas transport in the shales parallel to the bedding ranges from 80 to 983 nm at the experimental confining pressure. These values are much higher than the average pore diameters derived from low-pressure N 2 adsorption analysis ( Table 2). The pore diameter obtained from low-pressure N 2 adsorption analysis is a representation of pore structure of shale matrix, while the pulse-decay permeability is directional and strongly influenced by the occurrence of preferential pathways (laminae or microfractures). The effective transport pore diameter of siliceous shales ranges between 81 and 305 nm (averaging 191 nm), while the effective transport pore diameter of argillaceous shales ranges from 218 to 826 nm (averaging 552 nm). There is a negative correlation between TOC content and effective transport pore diameters of the shales studied (Figure 16(a)), while a positive correlation exists between clay mineral content and effective transport pore diameters of shale samples (Figure 16(b)). This further implies that there are developed preferred transport paths in argillaceous shales. Through comparison with the visual pore structure characteristics, it is thought that the intraparticle pores between the oriented minerals are the preferential transport pathways.
The Winland equation, originally developed from mercury intrusion capillary pressure (MICP) measurements of conventional reservoir rocks, was recently used to identify the flow unit and estimate the dominant pore throat diameters for fluid flow in shales [58]. Though the Winland equation is an empirical method based on MICP data, it provides visual information about the pore diameter and permeability. The Winland equation describing the rock permeability (mD), porosity (decimal), and r p35 (dominant porethroat aperture, μm) can be written as [58,59] r p35 = 2:665 k 100ϕ ! 0: 45 : ð3Þ Using the permeability and porosity values, the dominant pore throat diameter (2r p35 ) for fluid transport estimated from the Winland method ranges from 71 to 1386 nm ( Table 2). These values are slightly larger than those estimated from gas slippage measurements, which is probably related to the difference of tortuosity in shales and the conventional reservoir rocks. However, both these two methods show that the effective transport pore diameters of siliceous shales are commonly lower than those of argillaceous shales ( Figure 17). The dominated pore-throat diameter of siliceous shales lies between 50 and 300 nm, while that of argillaceous shales is between 100 and 1000 nm. These values underline the controls of microfabric and texture on transport pathways of Longmaxi shales.

Conclusions
Comprehensive investigations on reservoir characteristic of Longmaxi shales from a shale gas well in the Jiaoshiba area, Sichuan Basin, were conducted. Special attentions were focused on the effects of microfabric and laminae on pore structure and gas transport pathways of shales. The following conclusions can be drawn: (1) In general, there are three different kinds of lithofacies in the Longmaxi Formation: faintly laminated to nonlaminated siliceous shales (S), laminated mixed shales (M), and argillaceous shales (CM) with oriented clay flakes. Porosity of these Longmaxi shales studied is positively correlated with TOC content and ranges from 1.5 to 7.1%. Organic matter approximately contributes 79% to the pore spaces of siliceous shale samples, whereas interparticle pores related to inorganic matter contribute 63% to the pore spaces of argillaceous shale samples (3) Permeability of these shales mainly depends on the rock microfabric, other than TOC content and porosity. There is a negative correlation between permeability values and porosities of the samples studied. Samples with well-developed laminae and oriented clay flakes show higher permeability and permeability anisotropy (4) The effective transport pore diameter was estimated using two different methods: gas slippage measurements and the Winland equation. The effective transport pore diameters of argillaceous shales (averaging 552 nm) are commonly larger than those of siliceous shales (averaging 198 nm), which is related to the preferential transport pathways of laminated fabric and intraparticle pores of oriented clay minerals

Data Availability
Data will be made available upon request.

Conflicts of Interest
The authors declare that they have no conflicts of interest.