Experimental Investigation of Pore Structure and Movable Fluid Traits in Tight Sandstone

Whether the variation of pore structures and movable fluid characteristics enhance, deteriorate, or have no influence on reservoir quality has long been disputed, despite their considerable implications for hydrocarbon development in tight sandstone reservoirs. To elucidate these relationships, this study systematically analyzes pore structures qualitatively and quantitatively by various kinds of direct observations, indirect methods, and imaging simulations. We found that the uncertainty of porosity measurements, caused by the complex pore-throat structure, needs to be eliminated to accurately characterize reservoir quality. Bulk water was more easily removed, while surface water tended to be retained in the pores, and the heterogeneity of pore structures was caused by the abundance of tiny pores. The rates of water saturation reduction in macropores are faster than those for tiny pores, and sandstones with poor reservoir quality show no marked descending of lower limits of movable pore radius, indicating that the movable fluid would advance exempted from the larger pores. This study suggests that the deterioration of reservoir quality is strongly affected by the reduction of larger pores and the aqueous phases tended to remain in the tiny pores in the forms of surface water.


Introduction
With depletion of conventional hydrocarbon production and increasing demand for energy, unconventional reservoirs are expected to be a vital and realistic inventory in future energy supply [1,2].As one of the major unconventional resources, tight sandstones reservoirs, to which many scholars have drawn attention, have been developed successfully in China [3][4][5][6].The definition of tight sandstone in petroliferous basins is that its porosity is no more than 10% and in suit permeability generally less than 1 mD or air permeability less than 0.1 mD, whereupon it needs specific treatments (horizontal well fracturing, correct way of completion, advanced simulation jobs, etc.) to produce commercial hydrocarbon flow [7,8].The complexity and irregularity of the pore structure are crucial microscopic factors for the deteriorated physical properties, poor reservoir quality, and complicated oil emplacement; therefore, comprehensive and accurate pore geometry evaluation has become a momentous issue in tight sandstones [9,10].Meanwhile, due to the heterogeneity of pore structure in tight sandstones, the fluid that flows in the porous media is more complex than that of conventional reservoirs because of the intricate flow path [11].Movable fluid parameters are significant in that they can provide a reliable evaluation of the characteristics of the fluid, and accurate movable fluid parameter estimation has been a point of interest because movable fluid parameters are valuable when considering the water retention and residual oil accumulation and development [12].Therefore, the assessment of pore structure (porosity, pore size distribution (PSD), connectivity, tortuosity, heterogeneity, pore shape, etc.) and movable fluid features (movable fluid percentage, movable fluid porosity, etc.) are of paramount importance in evaluating the microscopic characterization of tight sandstones [13,14].Accurate prediction of pore structure and movable fluid traits helps to assess storage capacity and evaluate percolation capacity and can provide insights into planning field exploration and development [15].
Currently, a series of optical techniques and experiments are being conducted to investigate the pore structure of tight sandstones, including radiation methods (e.g., thin section (TS), scanning electron microscope (SEM), small angle neutron scattering (SANS), computed tomography (CT), etc.), intrusive techniques (e.g., pressure-controlled mercury intrusion (PCMI), rate-controlled mercury intrusion (RCMI), etc.) and non-intrusive techniques (e.g., nuclear magnetism resonance (NMR), low-temperature N 2 adsorption (LTNA), etc.) [2,8,[10][11][12][13][15][16][17].With continuous increase in the resolution of optical microscopy and radiation methods, the field of views could be small and the costs could rise.Due to this trade-off between resolution and viewshed information, indirect techniques have been supplemented, which systematically compensate the limits of direct techniques in characterizing pore structure in tight sandstones [8,14].Although many methods can be used to characterize pore structures, no sole direct and indirect experiment can be used to determine pore structures traits due to its own limitations and strengths [18]: TS and SEM can reveal pore morphology and occurrence qualitatively, but these techniques were unable to obtain quantitative pore structure-related data [19,20].CT is an effective technique to quantitatively characterize the three-dimensional images of pore networks, but it has drawbacks in the calculation of pores larger than 3 µm and is often restricted by expense [2,21,22].PCMI is understood to be a widely-used method of estimation of pore structures because it can detect a broad spectrum of PSD on account of relatively high injection pressure according to the Washburn equation [23]; however, this technique fails to acquire an exact count of large pores because of the shielding effect of small pores [24].Unlike PCMI, RCMI injects mercury into pores at very low speed that keeps the interfacial tension and contact angle changeless and can partition the PSD into areas of pores and throats based on the periodic rising and falling of pressures, overcoming the pitfall in PCMI [2,15,[25][26][27].While it is impossible to detect pores with radius less than 0.12 µm due to the maximum injection pressure limitation (~900 psi/6.2MPa), as a result, RCMI can be envisaged as a limited technique that exaggerates the proportion of large pores [10,28].NMR has been established as an effective method to determine pore structure characterization [29,30], nevertheless, the calculation of pore radius requires calibration with other techniques (e.g., RCMI) to confirm surface relaxivity, namely, it would cause discrepancy between NMR-derived PSD and other independent pore size measurement if the process of surface relaxivity calibration is skipped and referred to an arbitrary pore size classification that is derived from other research areas [15,31,32].Therefore, in order to better evaluate pore structure characterization, a combination of each technique is necessary [33,34].
Recently, many studies have documented that NMR is a fast and nondestructive method to determine flow characteristic in tight sandstones [11][12][13][35][36][37].As for the characteristic parameters of fluid in tight sandstones, movable fluid saturation is commonly used NMR-derived parameters, which is related to the pore structures and can be regulated by T 2 cutoff values [12].T 2 cutoff values, defined here as the decay time after the external magnetic simulation disappeared, correspond to threshold partitioning of the T 2 spectrum into two subcurves, namely, the free fluid areas (corresponding to T 2 higher than the threshold) and bound fluid areas (corresponding to T 2 lower than the threshold); hence, rigorous calculation of T 2 cutoff values is vital to exactly determine the percentage of free water [12].The determination of T 2 cutoff values needs to compare the T 2 spectrum before and after centrifuging in a certain rotation speed because it can dramatically change among different specimens, even in the same research area; rather, empirically determined T 2 cutoff values ascertaining the proportion of movable fluid on the basis of predecessors is too arbitrary, and accurate determination of optimum centrifugal forces to obtain T 2 cutoff values requires an understanding of the calibration model [12,13,36,38,39].Movable fluid percentage is a significant property used to predict estimated ultimate recovery (EUR), because the irreducible fluids were not productive at all [11]; thus, it is of great importance to quantitatively determine the movable fluid parameters, in order to provide insight into evaluating the fluid storage and percolation traits and help develop proper hydrocarbon exploitation strategies.
Considering the significance of pore structure and movable fluid for hydrocarbon development in tight sandstones, a comprehensive analysis of pore structures and evaluation of movable fluid are necessary to evaluate the mechanism of hydrocarbon storage and percolation in the investigated basin [40,41].In our study, a suit of tight sandstone samples from the Upper Triassic in the Ordos basin was investigated by TS, SEM, PCMI, RCMI, NMR, and CT in order to reveal the characteristics of pore structure and movable fluid traits.We first lay out the spatial features of the pores needed for the qualitative evaluation of our work by TS and SEM.Then we propose an integrated method by a combination of PCMI, RCMI, and CT to describe PSD, and an improved PSD is elucidated.Further, we determine the optimum rotation speed and movable fluid traits and explain how pore structure and rock physical properties influence the movable fluid parameters by NMR-derived T 2 spectrum analysis before and after centrifuging at different centrifugal force, T 1 -T 2 correlation spectra and magnetic resonance imaging (MRI).Finally, we investigate the effect of pore radius lower limits on reservoir quality.This multidisciplinary research unravels the impact of pore structure and movable fluid features on reservoir quality of tight sandstones, which is critical for the future production of the Upper Triassic tight sandstone reservoirs in Ordos Basin, China and offer technological merits in similar reservoirs elsewhere.

Tight Sandstones Specimens
All specimens were taken from the Upper Triassic Yanchang tight sandstones in the Yishan slope of the Ordos Basin in NW China (Figure 1a).It is a cratonic basin with an area of around 2.5 × 10 5 km 2 , and Yishan slope is a hydrocarbon enrichment area, where dips at no more than 1 • toward the west and internal faults are rarely seen [42,43].The Yanchang formation, which is divided into ten members (Figure 1b), consists of sandstones intercalated with mudstones and shales, and Chang 2, Chang 4 + 5, Chang 6 and Chang 8 members are composed predominantly of sandstones [44,45].The Yanchang formation is dominated by fine-to moderate-grained sands deposited in a fluvial, delta and lacustrine environment; compaction, cementation, and dissolution represent the most significant diagenetic process influencing pore structure and movable fluid traits [46,47].In the present study, sandstones specimens from the Yanchang formation were collected from 14 wells, and the porosities of all specimens range from 5.53% to 14.24% with an average of 9.33%, while permeabilities vary from 0.02 mD to 2.70 mD with a mean value of 0.68 mD (Table 1), which indicates that the sandstones are typical tight sandstones [7].The tight sandstones specimens of Yanchang formation were mainly arkose according to Folk's (1980) classification [47], and quartz is predominant in specimens (av.37.38%) with a subdominant component of feldspar (av.31.86%) and low rock fragments (average contents of metamorphic, volcanic and sedimentary rock fragments equal to 7.84%, 3.82%, and 4.56%, respectively) (Table 1).The clay minerals are dominated by chlorite (av.4.81%) and illite (av.4.72%), whereas kaolinite is less dominant (av.1.50%) (Table 1).The mica and quartz overgrowth contents vary from 1.00% to 11.50% (av.4.32%) and 1.00% to 4.00% (1.90%), respectively (Table 1).

Experimental Methods
First, the rock cores selected from wells were drilled parallel to the bedding plane and cylindrical core plugs with approximately 5 cm long and 2.54 cm in diameter were obtained.The specimens were polished and cleaned with solutions of alcohol and trichloromethane before the experiments in order to remove the bitumen or drilling mud that remained in pores from the specimens, and then each specimen was put into the vacuum systems and dried at 120 °C for 24 h.Subsequently, the specimens were put in the core holder of the FYK-I testing apparatus, then the helium-based porosity was repeatedly tested five times and averaged.After the porosity tests, the pressure-transient nitrogen pulse decay permeability test was carried out using the FYK-I instrument and the specimen was tautologically measured five times for an average.

Experimental Methods
First, the rock cores selected from wells were drilled parallel to the bedding plane and cylindrical core plugs with approximately 5 cm long and 2.54 cm in diameter were obtained.The specimens were polished and cleaned with solutions of alcohol and trichloromethane before the experiments in order to remove the bitumen or drilling mud that remained in pores from the specimens, and then each specimen was put into the vacuum systems and dried at 120 • C for 24 h.Subsequently, the specimens were put in the core holder of the FYK-I testing apparatus, then the helium-based porosity was repeatedly tested five times and averaged.After the porosity tests, the pressure-transient nitrogen pulse decay permeability test was carried out using the FYK-I instrument and the specimen was tautologically measured five times for an average.
After the physical property measurement, the specimen plugs were split into several chips and plugs for direct observations and indirect tests.The TS analyses were conducted on rock slices that were impregnated by red or blue epoxy resin, using ZEISS polarized light petrographic microscope to detect the pore geometry.SEM, which was used to observe the microscopic features of the tight sandstones, was conducted by a FEI Quanta 400 FEG scanning electron microscope with working accelerating voltage set at 20 kV and the samples were covered with 10 nm thick carbon films.An OXFORD IE 350 energy dispersive spectrometer was used to collect XRD data, and each specimen was crushed and ground into powder with a grain size of 300 mesh, then dried and mixed with ethylene glycol and heated to 550 • C for 3 h before the tests.
The PCMI test was conducted using a Micromeritic Autopore IV 9420 Instrument.Injection pressure was reached at 200 MPa, which corresponded to 3.6 nm based on the Washburn equation [23], then the pressure progressively decreased to zero and the intrusion and extrusion capillary pressure curves were acquired.The ASPE 730 mercury porosimeter was adopted for the purpose of RCMI measurements.Unlike the PCMI experiment, during the injection process in RCMI, mercury was injected into the specimen at a speed of 5 × 10 −5 mL•min −1 in order to stabilize the interfacial tension and contact angle.The pressure increased as the mercury was injected into the throat, then an instantaneous drawdown was seen when the mercury entered the pores.The periodic rising and falling pressure was recorded by the software and real-time data were acquired to partition the void spaces into distributions of pores and throats [26,38].Both mercury intrusion experiments followed the SY/T 5346-2005 Chinese Oil and Gas Industry Standard, and the contact angle and interfacial tension was 140 • and 480 mN•m −1 , respectively.
Micro-CT scanning can provide a relatively high-resolution 3D reconstruction of the features of the pore network, when compared to other methods [2,[50][51][52].The 3D images of the pore characterization of tight sandstone were measured by a Phoenix Nanotom M scanner with a working voltage of 90 kV, and tomogram image spatial resolution of around 3 µm.In this method, an X-ray with conical beam penetrates the specimen and then attenuates depending on the sample density, mineral compositions and contents, and the thickness along the beam direction, and this attenuation is mainly associated with the decreasing sample density [22,53].Finally, based on the threshold value segmentation method, rock matrix with high density and void space with low density can be determined [22].
The NMR is a nondestructive technique that provides information of pore features, and insights into the fluid and porous media can be acquired by this test [11,54].The process that helped the magnetization vector recover from a non-equilibrium to equilibrium state is called relaxation, and the recovery time is called relaxation time [11].The T 1 and T 2 relaxation time, which correspond to the time that the magnetization vector is parallel to and perpendicular to the external magnetic field returned to the equilibrium state, is referred to as longitudinal and transverse relaxation time, respectively [55].This technique was performed in Niumag NMR spectrometer at a proton resonance frequency of 2 MHz and detected by a CPMG pulse sequence [55].We set the echo number for 30 and 6000 for T 1 and T 2 , respectively, and the echo spacing was varied from 10 µs to a few seconds for T 1 and the values of which were equal to 100 µs for T 2 .Each specimen was scanned 64 times to get a good signal-to-noise ratio, and the wait time between successive scans was chosen as 5 s.The specimens were saturated in brine, which was a mixture of deionized water and Calcium Chloride (CaCl 2 ) with salinity of 2.5 × 10 4 mg•L −1 to prevent the clay minerals from swelling.It was used for NMR measurement to record initial relaxation time spectra and images.After the T 1 and T 2 spectrum and MRI for the saturated water specimens were determined, the centrifugal processing with a rotation speed of 2500 r•min −1 , 2900 r•min −1 , 3500 r•min −1 , 5000 r•min −1 , 7900 r•min −1 and 9100 r•min −1 for 1 h for each speed were conducted on each specimen, and the T 2 spectra distributions were tested after different rotation speed centrifugal processing was completed.Specifically, the T 1 -T 2 correlation spectra distributions and MRI were generated after 9100 r•min −1 (corresponding to the maximum rotation speed) for centrifugal processing was implemented to determine the T 2 cutoff values, free and irreducible water distributions.

Pore Structures Observed by TS and SEM
TS and SEM can be used to determine pore types and analyses pore geometry since they provide direct characterization of the pores [17,56].Four main kinds of pore types are developed in the Upper Triassic Yanchang formation of tight sandstones, including intergranular pores, dissolution pores, intercrystalline pores, and microfractures.The intergranular pores are triangular or polygonal and remain after strong mechanical compaction and clay minerals recrystallization [20], have characteristics of relatively larger pore radius, pore-lining chlorite, and sometimes elongate euhedral quartz prisms (Figure 2a,b).The size of the intergranular pores was typically >40 µm, and the relative abundances of these pores by point counting based on TS observations were in the range from trace levels to 5.9% and averaged as 1.73%.The dissolution pores are generally a result of feldspar and unstable ferromagnesian rock fragments dissolution and show irregular geometry [2,57]; these pore systems are of three categories: intragranular dissolution pores (feldspar leaching generally developed along cleavage planes) (Figure 2c), intergranular dissolution pores (Figure 2d) and moldic pores (the pores derived from complete dissolution of grains and the shapes of the original framework were retained) (Figure 2e) [58].Point counting results of the dissolution pores reveal a wide range-from trace amounts (due to feldspars leaching) to 1.0% (resulting from the connection of moldic pores and residual intergranular pores) with an average of 0.56%.The intercrystalline pores, which refers to pores within biotite planes, clay aggregates, quartz overgrowths, and carbonate minerals, are commonly small enough (often lower than 10 µm) to be hardly observed in TS; however, the examination under SEM suggested that lamellar chlorite and fibrous illite have widespread intercrystalline pores because of its abundance (Figure 2f,g).The extensive cementation of carbonate minerals can contribute little intercrystalline pores owing to the compacted voids in the lattice (Figure 2h), and pores within sheet like biotite and elongated authigenic quartz only in minor amounts (Figure 2i,j).The microfractures, which are dominated by diagenesis and scarcely observed in cores [59], offer features of microscale in length and nanoscale in width, and can only be detected under microscopic identifications.Microscopic observations of the specimens reveal that the microfractures along the edge of grains and through the grains are the predominate microfracture type in tight sandstones (Figure 2k,l), from which acid solution dissolution and mechanical compaction was derived, respectively [60].Point counting results show that the contents of microfractures are minimal (<0.2%).

Pore Structures Characteristics Determined by PCMI
The PCMI-derived curves of mercury drainage and imbibition for the 14 specimens are presented in Figure 3 after Gane calibration [61] in order to make the curves smooth and continuous; the typical parameters for pore structures are listed in Table 2.The specimens could be categorized into four types based on intrusion and extrusion capillary curves.The first type (type α) shows low threshold pressure (P d ) (av. 0.07 MPa), low median pressure (P 50 ) (av. 0.81 MPa), high maximum mercury intrusion saturation (S Hgmax ) (av. 89.19%) and medium hysteresis between intrusion and extrusion curves, revealing that well-connected pores play a dominant role in this type of sandstone.Besides, partial percolations (plateau stage in mercury intrusion curves), corresponding to the distinct change in saturation over a small range of capillary pressure (gray boxes)), occurred in some specimens, indicating this type is better sorted (average sorting coefficient equal to 2.49) than other types (Figure 3a).For type β specimens, the average P d and P 50 rapidly increase, equal to 0.95% and 8.09%, respectively, and the average S Hgmax decreases (av.86.33%), while a relatively strong hysteresis is observed (Figure 3b).These phenomena suggest that some poorly-connected pores were developed in this type of sandstone.For the type γ specimens, the mercury intrusion curves display that the P d for all the specimens exceed 1.16 MPa, P 50 over 7.06 MPa, and S Hgmax is no more than 87.57%, that is, mercury can hardly enter the pores and the proportion of non-connected pores is relatively high, when compared to type α and β.The intrusion-extrusion cycles show more moderate mercury saturation rise and fall, corresponding to weak hysteresis, and no percolation phenomenon.However, the specimen that has strong hysteresis may result from an abundance of mica, which is deformed in the process of mercury injection and leads to the change of pore structures (Figure 3c).In terms of type δ, the average P d and P 50 are the highest, equal to 3.67 MPa and 51.47 MPa, respectively, while the average S Hgmax is no more than 83.85%.Strong hysteresis of these specimens is the evidence of poor-connected pores and minor amounts of clay minerals due to intense mechanical compaction (Figure 3d).According to the multitype pore space model proposed by Sakhaee-Pour and Bryant [40], medium hysteresis for type α suggests that intergranular pores mainly constitute the pore space, whereas the void spaces in type γ is predominantly intercrystalline pores, which corresponds to a relatively poor hysteresis; in addition, the dissolution pores may make the greatest contribution to the void spaces in type β because their diameters are midway between intergranular pores and intercrystalline pores, according to the description in Section 3.1.Figure 4 backs this and explains why more mercury was stranded in the spaces in intergranular and dissolution pore-dominated reservoirs, whereas the void spaces of intercrystalline pores can be thoroughly saturated with the wetting phase.

Pore Structures Characteristics Determined by PCMI
The PCMI-derived curves of mercury drainage and imbibition for the 14 specimens are presented in Figure 3 after Gane calibration [61] in order to make the curves smooth and continuous; the typical parameters for pore structures are listed in Table 2.The specimens could be categorized into four types based on intrusion and extrusion capillary curves.The first type (type α) shows low threshold pressure (Pd) (av.0.07 MPa), low median pressure (P50) (av.0.81 MPa), high maximum mercury intrusion saturation (SHgmax) (av.89.19%) and medium hysteresis between intrusion and corresponds to a relatively poor hysteresis; in addition, the dissolution pores may make the greatest contribution to the void spaces in type β because their diameters are midway between intergranular pores and intercrystalline pores, according to the description in Section 3.1.Figure 4 backs this and explains why more mercury was stranded in the spaces in intergranular and dissolution poredominated reservoirs, whereas the void spaces of intercrystalline pores can be thoroughly saturated with the wetting phase.  2 Pd-threshold pressure; rmax-maximum pore radius; P50-median pressure; r50-median pore radius; SHgmax-maximum mercury intrusion saturation; We-efficiency of mercury withdrawal; Sort-sorting coefficient.The PCMI-derived PSD can be calculated on the basis of the Washburn equation [23].The PSD of type α sandstones has a wide size range and possesses high peakedness over 1.0 µm with a long tail in the rests of pore radius, suggesting that specimens with high physical properties contain relatively larger pores and well-sorted PSD (Figure 5a).The PSD results of type β give two locations of the maxima, of which specimens 1# and 4# correspond to the wide range pore radius with bimodal, whereas other specimens show unimodal; the pores with radius larger than 1.0 µm are rare in this type (Figure 5b).The PSD of the type γ sandstones mainly varies from 0.01 µm to 0.4 µm, but the corresponding pore sizes of the peak pore volume values are different, from 0.008 µm to approximately 0.2 µm, indicating a heterogeneity of PSD with similar porosity and permeability as those of type γ (Figure 5c).Quantitative statistics performed on type δ sandstones indicate that the PSD mainly range from 0.003 µm to 0.15 µm, the main crests of these specimens being mainly centered at around 0.015 µm, whereas specimen 5# is accompanied by a peak larger than 0.1 µm, which might show evidence of microfractures (Figure 5d).The rapid increase of frequency in 0.003 µm represent the abundance of tiny pores (mainly the pores among clay minerals or the pores within the aluminosilicate layers); hence, porosity and permeability may be limited because of the rich clay minerals.The results of all the specimens show that the pore radius is distributed over a wide range and exhibits a multimodal PSD, and with the deterioration of physical properties, the PSD becomes more fluctuant, revealing that the pore structures are more complicated and heterogeneous.However, pores with radius over 20 µm are rare, which show a discrepancy between the PCMI and petrographic observations because many large pores are shielded by connected smaller pores.Hence, PCMI provides imperfect information of larger pores and can only uncover PSD of smaller ones, accordingly, the full range of PSD determination is required for PCMI in combination with other methods.
those of type γ (Figure 5c).Quantitative statistics performed on type δ sandstones indicate that the PSD mainly range from 0.003 μm to 0.15 μm, the main crests of these specimens being mainly centered at around 0.015 μm, whereas specimen 5# is accompanied by a peak larger than 0.1 μm, which might show evidence of microfractures (Figure 5d).The rapid increase of frequency in 0.003 μm represent the abundance of tiny pores (mainly the pores among clay minerals or the pores within the aluminosilicate layers); hence, porosity and permeability may be limited because of the rich clay minerals.The results of all the specimens show that the pore radius is distributed over a wide range and exhibits a multimodal PSD, and with the deterioration of physical properties, the PSD becomes more fluctuant, revealing that the pore structures are more complicated and heterogeneous.However, pores with radius over 20 μm are rare, which show a discrepancy between the PCMI and petrographic observations because many large pores are shielded by connected smaller pores.Hence, PCMI provides imperfect information of larger pores and can only uncover PSD of smaller ones, accordingly, the full range of PSD determination is required for PCMI in combination with other methods.

Pore Structures Characteristics Determined by RCMI
RCMI provides a reliable way to reflect detailed pores and throat structures.Since the void space can be divided into pore region and throat region, based on injection pressure fluctuation, the

Pore Structures Characteristics Determined by RCMI
RCMI provides a reliable way to reflect detailed pores and throat structures.Since the void space can be divided into pore region and throat region, based on injection pressure fluctuation, the distribution of pores and throats radius are determined [15].In order to investigate RCMI-derived PSD, six typical specimens with different physical properties were selected (Table 3).The capillary pressure curves of the specimens show that in the primary stage, the curves of total mercury intrusion saturation followed a trend of pore mercury intrusion saturation, while with the increasing pressure, the total mercury intrusion saturation was contributed to by the throat intrusion saturation and the intrusion pressure of the pore rapidly increased (Figure 6).Based on the trend of pore intrusion mercury saturation curves, the RCMI-derived PSD can be divided into a pore-dominated region corresponding to the stage that the pore mercury intrusion saturation increased exponentially, a throat-dominated region corresponding to the areas that the mercury saturation of pore remained unchanged, whereas the total mercury intrusion followed the trend of throat mercury intrusion, and a hybrid region by which both pores and throats are controlled (Figure 6).Furthermore, all specimens could be assigned to two groups based on the proportion of different regions.Group I exhibited a relatively lower threshold pressure (P d ) and medium pressure (P 50 ), while relatively higher maximum total mercury intrusion saturation (S Hgmax ); the percentage of pore-dominated and hybrid region exceeded that of throat-dominated, and shows a notable flat segment corresponding to the relatively good sorting of pores and throats (Table 3; Figure 6a,c,d).For Group II, the P d increased while the S Hgmax decreased and P 50 disappeared.Moreover, the proportion of throat-dominated region surpasses that of the pore-dominated region in this group and no hybrid region existed, showcasing characteristics of low porosity and explaining why the storage of this group was relatively poor (Table 3; Figure 6e,f).Specimen 5# has relatively low P d and S Hgmax , indicating that it may contain microfractures that lead to the reduction of entry pressure (Table 3; Figure 6b).a hybrid region by which both pores and throats are controlled (Figure 6).Furthermore, all specimens could be assigned to two groups based on the proportion of different regions.Group I exhibited a relatively lower threshold pressure (Pd) and medium pressure (P50), while relatively higher maximum total mercury intrusion saturation (SHgmax); the percentage of pore-dominated and hybrid region exceeded that of throat-dominated, and shows a notable flat segment corresponding to the relatively good sorting of pores and throats (Table 3; Figure 6a,c,d).For Group II, the Pd increased while the SHgmax decreased and P50 disappeared.Moreover, the proportion of throat-dominated region surpasses that of the pore-dominated region in this group and no hybrid region existed, showcasing characteristics of low porosity and explaining why the storage of this group was relatively poor (Table 3; Figure 6e,f).Specimen 5# has relatively low Pd and SHgmax, indicating that it may contain microfractures that lead to the reduction of entry pressure (Table 3; Figure 6b).The RCMI can test pores and throats separately; hence, the distribution of the pore radius, throat radius, and pore-throat radius ratio is recorded simultaneously [2].The pores' radius mainly ranges from 88.23 μm to 230.15 μm with main crests of approximately 105 μm, showing no obvious discrepancy between different specimens (Table 3, Figure 7a).However, the throat radius demonstrates variations among specimens with distinctive features in spectral peaks and radius spectra, mainly ranging from 0.15 μm to 2.54 μm with an average of 1.06 μm (Table 3, Figure 7b).The distribution of pore-throat radius ratio is principally in the range of 55.21-430.78 with the average value of 206.70 (Table 3, Figure 7c).Unlike PCMI, RCMI can provide more reliable results due to the quasi-statics mercury injection process; however, this method generally leads to similar pore radius distributions among various specimens.In addition, the pore radius derived from RCMI shows that tons of pores larger than 100 μm existed, while pores with radius between 10 μm and 100 μm are The RCMI can test pores and throats separately; hence, the distribution of the pore radius, throat radius, and pore-throat radius ratio is recorded simultaneously [2].The pores' radius mainly ranges from 88.23 µm to 230.15 µm with main crests of approximately 105 µm, showing no obvious discrepancy between different specimens (Table 3, Figure 7a).However, the throat radius demonstrates variations among specimens with distinctive features in spectral peaks and radius spectra, mainly ranging from 0.15 µm to 2.54 µm with an average of 1.06 µm (Table 3, Figure 7b).The distribution of pore-throat radius ratio is principally in the range of 55.21-430.78 with the average value of 206.70 (Table 3, Figure 7c).Unlike PCMI, RCMI can provide more reliable results due to the quasi-statics mercury injection process; however, this method generally leads to similar pore radius distributions among various specimens.In addition, the pore radius derived from RCMI shows that tons of pores larger than 100 µm existed, while pores with radius between 10 µm and 100 µm are seldom developed; nevertheless, it does not match with the TS and SEM observations (Figure 2).This comparison reveals that the RCMI may exaggerate the radius of the pores; therefore, necessary corrections for pore radius curves need to be made with the help of other methods.

T1 Relaxation Time before and after Centrifugal Processing
T1 measurements in a fully saturated state and centrifuged state were taken for the specimens; these spectra are shown in Figure 9.The T1 spectra of specimens show multimodal characteristics with the maximum crests being mainly distributed in the range of 51.79-719.69ms and 4.64-64.49ms in a fully saturated state and a centrifuged state, respectively.However, the major peaks show a broad massif-like peak in the fully saturated state and a narrow ridge-like peak in the centrifuged

Pore-Throat Microstructures Derived from Micro-CT
The pore-throat microstructures depicted by micro-CT allow accurate and detailed investigation of PSD; the pore-throat networks are illustrated in Figure 8.The matrix of the specimens is characterized by dark gray while the pore-throats are in color.The mutually inter-connected pores correspond to the same color, while the disconnected pores are presented by different colors or an isolated part.Due to the expense and time restriction, only two typical specimens with different permeability were selected for this measurement.Specimen 11# had relatively high porosity and permeability (ϕ: 14.24%, K: 2.70 mD); a mass of pores is mainly displayed in yellow and green, and the percentage of gray is relatively low, suggesting that this specimen had good pore-throat connectivity (Figure 8A,B).Many pores were found to be tubular or in larger spherical shapes, a few were isolated narrowed belt-like or small spherical pores, showcasing that intergranular and dissolution pores play a leading role in high permeability specimens (Figure 8A,B).For specimen 12#, which had relatively low porosity and permeability (ϕ: 12.54%, K: 0.48 mD), the pores were sparsely distributed, when compared to specimen 11#, and the same color was sporadically distributed, suggesting poor connectivity of these specimens.The pore radius distribution of these specimens show that a few narrow pores are found in 11#, whereas abundant pores with radius around 30 µm play important roles in forming major void spaces, which implies that micro-pores determined the storage and transport properties in tight sandstone (Figure 8C,F).

T1 Relaxation Time before and after Centrifugal Processing
T1 measurements in a fully saturated state and centrifuged state were taken for the specimens; these spectra are shown in Figure 9.The T1 spectra of specimens show multimodal characteristics with the maximum crests being mainly distributed in the range of 51.79-719.69ms and 4.64-64.49ms in a fully saturated state and a centrifuged state, respectively.However, the major peaks show a broad massif-like peak in the fully saturated state and a narrow ridge-like peak in the centrifuged state; the corresponding T1 relaxation time and amplitude moves from large to small (Figure 9).The spurious crests at high (>8000 ms) T1 value in Figure 9a

T 1 Relaxation Time before and after Centrifugal Processing
T 1 measurements in a fully saturated state and centrifuged state were taken for the specimens; these spectra are shown in Figure 9.The T 1 spectra of specimens show multimodal characteristics with the maximum crests being mainly distributed in the range of 51.79-719.69ms and 4.64-64.49ms in a fully saturated state and a centrifuged state, respectively.However, the major peaks show a broad massif-like peak in the fully saturated state and a narrow ridge-like peak in the centrifuged state; the corresponding T 1 relaxation time and amplitude moves from large to small (Figure 9).The spurious crests at high (>8000 ms) T 1 value in Figure 9a

T2 Relaxation Time before and after Centrifugal Processing
Unlike T1 relaxation time, tests of T2 relaxation time are preferred, as they are not timeconsuming and can provide more accurate information about the void space by liquid-state NMR [32,64].The T2 distributions in the fully saturated and centrifugated state show bimodal or multimodal features with two crests; the short T2 values represents micropores, mainly distributed in the range of 0.11-1.54ms, while the long T2 relaxation time corresponding to larger pores and cracks is mainly distributed in the range of 5.63-382.79ms (Figure 10a).The T2 spectra show similar shape before and after the centrifugal processing, while amplitude is narrowed in the irreducible water state (Figure 10).The left peak shows a relatively narrow peak distribution and larger amplitude than that of the right peaks in the fully saturated state, indicating that micropores have a

T 2 Relaxation Time before and after Centrifugal Processing
Unlike T 1 relaxation time, tests of T 2 relaxation time are preferred, as they are not time-consuming and can provide more accurate information about the void space by liquid-state NMR [32,64].The T 2 distributions in the fully saturated and centrifugated state show bimodal or multimodal features with two crests; the short T 2 values represents micropores, mainly distributed in the range of 0.11-1.54ms, while the long T 2 relaxation time corresponding to larger pores and cracks is mainly distributed in the range of 5.63-382.79ms (Figure 10a).The T 2 spectra show similar shape before and after the centrifugal processing, while amplitude is narrowed in the irreducible water state (Figure 10).The left peak shows a relatively narrow peak distribution and larger amplitude than that of the right peaks in the fully saturated state, indicating that micropores have a higher percentage, while meso-and macropores have a wider pore radius range [65].Besides, amplitude reduction with a slight drop in the left peak and distinct descending right peaks in the centrifugated state reveal that the irreducible water is mainly stranded in the micropores after centrifugal processing (Figure 10).

T2 Relaxation Time before and after Centrifugal Processing
Unlike T1 relaxation time, tests of T2 relaxation time are preferred, as they are not timeconsuming and can provide more accurate information about the void space by liquid-state NMR [32,64].The T2 distributions in the fully saturated and centrifugated state show bimodal or multimodal features with two crests; the short T2 values represents micropores, mainly distributed in the range of 0.11-1.54ms, while the long T2 relaxation time corresponding to larger pores and cracks is mainly distributed in the range of 5.63-382.79ms (Figure 10a).The T2 spectra show similar shape before and after the centrifugal processing, while amplitude is narrowed in the irreducible water state (Figure 10).The left peak shows a relatively narrow peak distribution and larger amplitude than that of the right peaks in the fully saturated state, indicating that micropores have a higher percentage, while meso-and macropores have a wider pore radius range [65].Besides, amplitude reduction with a slight drop in the left peak and distinct descending right peaks in the centrifugated state reveal that the irreducible water is mainly stranded in the micropores after centrifugal processing (Figure 10).

Full PSD Calculated by Integration of PCMI, RCMI, and NMR
As mentioned above, PCMI, RCMI, micro-CT, and NMR tests were conducted on the same specimens; however, these tests do not lead to coincident results because of the limitations and drawbacks of these experiments.Thus, no single test can be used to investigate total PSD in tight sandstones; hence, the experimental data of these experiments may be modified for the sake of precise determination.PCMI can determine the characteristics of tiny pores but acquire defective larger pore radius spectra due to the shielding effect and mercury snap-off during extrusion [24,66].RCMI could quantify the throat, but had a lower throat radius limit of 0.12 μm, while the pores' spectra detected by this method seems unreliable, because the overall pore radius of different specimens determined by RCMI are mainly in the range of 100-200 μm; observation of TS reveals few pores larger than 100 μm in the research area (Figure 2; Figure 7a).Besides, the pore radius distributions of the studied

Full PSD Calculated by Integration of PCMI, RCMI, and NMR
As mentioned above, PCMI, RCMI, micro-CT, and NMR tests were conducted on the same specimens; however, these tests do not lead to coincident results because of the limitations drawbacks of these experiments.Thus, no single test can be used to investigate total PSD in tight sandstones; hence, the experimental data of these experiments may be modified for the sake of precise determination.PCMI can determine the characteristics of tiny pores but acquire defective larger pore radius spectra due to the shielding effect and mercury snap-off during extrusion [24,66].RCMI could quantify the throat, but had a lower throat radius limit of 0.12 µm, while the pores' spectra detected by this method seems unreliable, because the overall pore radius of different specimens determined by RCMI are mainly in the range of 100-200 µm; observation of TS reveals few pores larger than 100 µm in the research area (Figure 2; Figure 7a).Besides, the pore radius distributions of the studied area are almost in accordance with that of other basins, such as the Songliao basin [15], suggesting that calibration needs to be done in order to acquire accurate pore radius spectra.Moreover, the RCMI-derived PSD shows that there is a "missing zone", namely, the pores radius range from 20 µm to 80 µm is rarely presented.These pores are commonly found in the TS and SEM micrographs (Figures 2 and 7), suggesting that this experiment may exaggerate the radius of mesopores.Micro-CT scanning may be a reasonable method to characterize PSD; however, it requires a trade-off between expense and representativeness.Moreover, the resolution of this method makes it unable to estimate the number of relatively tiny pores.The isolated pores, which have fewer contributions to storage capacity and no contributions to percolation ability, can be detected; thereupon, this test is inapplicable for the determination of PSD because commercial flow is often restricted by isolated pores [2,21,22,51].The T 2 relaxation time, derived from the NMR test, can be used to describe PSD as a nondestructive method that does not destroy in-suit pore structures; however, this method requires calibration with an independent test to calculate surface relaxivity and specific surface area and then convert the relaxation time to pore radius [15,21,22,31,32].
In order to overcome these shortcomings, an integrated method is necessary.First, fitting procedures of T 2 relaxation time and mercury intrusion need to be done.As discussed above, the PCMI-derived PSD is unreliable due to its shielding effect and induced cracks that resulted from high mercury injection pressure.There is a discrepancy between RCMI-derived pore radius distributions and optical observation results; hence, we use RCMI-derived throat radius curves to calibrate the T 2 spectra to PSD.According to previous research [16,65], the relationship between T 2 relaxation time of fluid in the porous materials and pore radius can be expressed as follows: where r t is the pore-throat radius, µm; C and n can be regarded as a constant parameter.In this way, T 2 relaxation time can be transformed to pore radius if C and n are determined.The procedures are listed below: (1) The maximum throat intrusion saturation is multiplied by specimen-measured helium porosity ϕ to obtain throat porosity, ϕ t .And then, the product of the incremental throat mercury intrusion saturation and ϕ t can be calculated as incremental throat mercury intrusion porosity.(2) The incremental amplitude of T 2 relaxation time is multiplied by helium porosity ϕ to obtain incremental porosity, ϕ t2 .(3) The cumulative porosity according to ϕ t and ϕ t2 should be calculated, and the comparison of both curves obtained.The results show that there is a good match between RCMI-and NMR-derived cumulative porosity, suggesting that these procedures and methods are reasonable and reliable for calibration (Figure 11a).Then the linear least square method is used to determine the most appropriate C and n (Figure 11b).At last, the T 2 relaxation time can be transformed to the pore radius, and the NMR-derived PSD is obtained.
storage capacity and no contributions to percolation ability, can be detected; thereupon, this test is inapplicable for the determination of PSD because commercial flow is often restricted by isolated pores [2,21,22,51].The T2 relaxation time, derived from the NMR test, can be used to describe PSD as a nondestructive method that does not destroy in-suit pore structures; however, this method requires calibration with an independent test to calculate surface relaxivity and specific surface area and then convert the relaxation time to pore radius [15,21,22,31,32].
In order to overcome these shortcomings, an integrated method is necessary.First, fitting procedures of T2 relaxation time and mercury intrusion need to be done.As discussed above, the PCMI-derived PSD is unreliable due to its shielding effect and induced cracks that resulted from high mercury injection pressure.There is a discrepancy between RCMI-derived pore radius distributions and optical observation results; hence, we use RCMI-derived throat radius curves to calibrate the T2 spectra to PSD.According to previous research [16,65], the relationship between T2 relaxation time of fluid in the porous materials and pore radius can be expressed as follows: where rt is the pore-throat radius, μm; C and n can be regarded as a constant parameter.In this way, T2 relaxation time can be transformed to pore radius if C and n are determined.The procedures are listed below: (1) The maximum throat intrusion saturation is multiplied by specimen-measured helium porosity φ to obtain throat porosity, φt.And then, the product of the incremental throat mercury intrusion saturation and φt can be calculated as incremental throat mercury intrusion porosity.(2) The incremental amplitude of T2 relaxation time is multiplied by helium porosity φ to obtain incremental porosity, φt2.(3) The cumulative porosity according to φt and φt2 should be calculated, and the comparison of both curves obtained.The results show that there is a good match between RCMI-and NMR-derived cumulative porosity, suggesting that these procedures and methods are reasonable and reliable for calibration (Figure 11a).Then the linear least square method is used to determine the most appropriate C and n (Figure 11b).At last, the T2 relaxation time can be transformed to the pore radius, and the NMR-derived PSD is obtained.The specimens were saturated with simulated formation water before it was put into the NMR apparatus; thus, some dead pores and ink-bottle pores saturated with fluid are inevitable.However, The specimens were saturated with simulated formation water before it was put into the NMR apparatus; thus, some dead pores and ink-bottle pores saturated with fluid are inevitable.However, this non-or weak-connected pores contributed less to the percolation [65].Besides, as mentioned above, the RCMI-derived throat distributions are reliable, although it cannot detect pores that are lower than 0.12 µm; hence, the integration of PCMI and NMR-derived PSD is regarded as a valid way to acquire the full range of PSD.The procedures are listed as follows: (1) the incremental porosity of improved NMR-derived PSD, which is more than 0.12 µm, is cumulated to obtain NMR cumulative porosity, ϕ N .(2) Subtract the ϕ N from helium measured porosity ϕ to get micropores porosity, ϕ M .
(3) The incremental mercury intrusion saturation of PCMI, which has radius less than 0.12 µm, is multiplied by ϕ to acquire the incremental PCMI-derived porosity.(4) The ratio of ϕ M and ϕ is multiplied by the incremental PCMI-derived porosity to obtain calibrated PCMI-derived porosity.After the improved NMR-derived porosity and PCMI-derived porosity are determined, the full range PSD can be transformed to that of Figure 12a.The improved PSD matches well with the microscopic observations (Figure 2), demonstrating that this distribution is reasonable.In the same way, the calibration equations for all specimens can be acquired and displayed in Figure 12b.0.12 μm, is multiplied by φ to acquire the incremental PCMI-derived porosity.(4) The ratio of φM and φ is multiplied by the incremental PCMI-derived porosity to obtain calibrated PCMI-derived porosity.After the improved NMR-derived porosity and PCMI-derived porosity are determined, the full range PSD can be transformed to that of Figure 12a.The improved PSD matches well with the microscopic observations (Figure 2), demonstrating that this distribution is reasonable.In the same way, the calibration equations for all specimens can be acquired and displayed in Figure 12b.

T2 Cutoff Values
T2 cutoff is a significant value in the NMR experiment because the movable fluid saturation can be determined by this parameter [13].This value is not a fixed constant, due to the variation of surface relaxivity and specific surface area of different specimens; centrifugal processing is an effective way to calculate this value.Figure 13 demonstrates how to calculate the T2 cutoff: first, the cumulative porosity of T2 spectrum in the fully saturated and irreducible water state is obtained, then a horizontal projection from the post centrifugation curve is drawn, and corresponding T2 value of the intersection of this horizontal line and cumulative curve of the fully saturated state is the T2 cutoff.We repeat this centrifugal process with various rotation speeds.In order to acquire the most accurate movable fluid saturation and determine optimum rotation speed, Coates model [65,67] was used:

T 2 Cutoff Values
T 2 cutoff is a significant value in the NMR experiment because the movable fluid saturation can be determined by this parameter [13].This value is not a fixed constant, due to the variation of surface relaxivity and specific surface area of different specimens; centrifugal processing is an effective way to calculate this value.Figure 13 demonstrates how to calculate the T 2 cutoff: first, the cumulative porosity of T 2 spectrum in the fully saturated and irreducible water state is obtained, then a horizontal projection from the post centrifugation curve is drawn, and corresponding T 2 value of the intersection of this horizontal line and cumulative curve of the fully saturated state is the T 2 cutoff.0.12 μm, is multiplied by φ to acquire the incremental PCMI-derived porosity.(4) The ratio of φM and φ is multiplied by the incremental PCMI-derived porosity to obtain calibrated PCMI-derived porosity.After the improved NMR-derived porosity and PCMI-derived porosity are determined, the full range PSD can be transformed to that of Figure 12a.The improved PSD matches well with the microscopic observations (Figure 2), demonstrating that this distribution is reasonable.In the same way, the calibration equations for all specimens can be acquired and displayed in Figure 12b.

T2 Cutoff Values
T2 cutoff is a significant value in the NMR experiment because the movable fluid saturation can be determined by this parameter [13].This value is not a fixed constant, due to the variation of surface relaxivity and specific surface area of different specimens; centrifugal processing is an effective way to calculate this value.Figure 13 demonstrates how to calculate the T2 cutoff: first, the cumulative porosity of T2 spectrum in the fully saturated and irreducible water state is obtained, then a horizontal projection from the post centrifugation curve is drawn, and corresponding T2 value of the intersection of this horizontal line and cumulative curve of the fully saturated state is the T2 cutoff.We repeat this centrifugal process with various rotation speeds.In order to acquire the most accurate movable fluid saturation and determine optimum rotation speed, Coates model [65,67] was used: IPBC-incremental porosity before centrifugation; IPAC-incremental porosity after centrifugation; CPBC-cumulative porosity before centrifugation; CPAC-cumulative porosity after centrifugation.
We repeat this centrifugal process with various rotation speeds.In order to acquire the most accurate movable fluid saturation and determine optimum rotation speed, Coates model [65,67] was used: where K C represents calculated permeability (mD), ϕ represents measured porosity (%), C c is a constant, FFI represents movable fluid saturation, and BVI represents irreducible fluid saturation.Hence, the linear regressions between ϕ 4 ( FFI BVI ) 2 under various rotation speeds and measured permeability (K a ) are performed to determine the optimum rotation speed.Figure 14 illustrate that 9600 r•min −1 can best fit Equation (2) with an R 2 value of 0.9428.The values of optimum T 2 cutoff are listed in Table 4.
where KC represents calculated permeability (mD), φ represents measured porosity (%), Cc is a constant, FFI represents movable fluid saturation, and BVI represents irreducible fluid saturation.Hence, the linear regressions between φ 4 ( ) 2   After determining the optimum rotation speed (9100 r•min −1 ), the movable fluid parameters can be calculated.In NMR-derived relaxation time studies, the movable fluid saturation (MFS) and movable fluid porosity (MFP) reflect the volumetric fraction of movable fluid occurring in the voids and the practical fluid flow ability through the sandstones, respectively; thus, there is a significant parameter that needs to be considered in evaluating reservoir quality [12].Reservoir quality is defined here as physical properties, including porosity and permeability.The results show that the relationship between the reservoir quality and movable fluid saturation is inferior to that of the movable fluid porosity; the correlation between movable fluid parameters and porosity are lower than that of permeability (Tables 1 and 4, Figure 15).These results demonstrate that the uncertainty of the porosity measurements, caused by the complex pore-throat structure, needs to be eliminated  After determining the optimum rotation speed (9100 r•min −1 ), the movable fluid parameters can be calculated.In NMR-derived relaxation time studies, the movable fluid saturation (MFS) and movable fluid porosity (MFP) reflect the volumetric fraction of movable fluid occurring in the voids and the practical fluid flow ability through the sandstones, respectively; thus, there is a significant parameter that needs to be considered in evaluating reservoir quality [12].Reservoir quality is defined here as physical properties, including porosity and permeability.The results show that the relationship between the reservoir quality and movable fluid saturation is inferior to that of the movable fluid porosity; the correlation between movable fluid parameters and porosity are lower than that of permeability (Tables 1 and 4, Figure 15).These results demonstrate that the uncertainty of the porosity measurements, caused by the complex pore-throat structure, needs to be eliminated to characterize reservoir quality accurately; in addition, the movable fluid parameters are more closely aligned with permeability.Meanwhile, with decreasing reservoir quality and movable fluid saturation, chlorite decreases from 7% to 0%, while illite increases from 0% to 7%, indicating that illite exerts a significant impact on retaining water, whereas chlorite can retard compaction (Figure 16a-c).Besides, reduction rate of reservoir quality with movable fluid saturation tends to decrease in tight sandstones with intergranular pores-dominated reservoirs to those with clay minerals, which indicates that movable fluid parameters are sensitive to a higher content of clay minerals or the proportion of tiny throats, whereas the reduction of intergranular pores play a clearly important role in the descending permeability (Figure 16a,b,d).
to characterize reservoir quality accurately; in addition, the movable fluid parameters are more closely aligned with permeability.Meanwhile, with decreasing reservoir quality and movable fluid saturation, chlorite decreases from 7% to 0%, while illite increases from 0% to 7%, indicating that illite exerts a significant impact on retaining water, whereas chlorite can retard compaction (Figure 16ac).Besides, reduction rate of reservoir quality with movable fluid saturation tends to decrease in tight sandstones with intergranular pores-dominated reservoirs to those with clay minerals, which indicates that movable fluid parameters are sensitive to a higher content of clay minerals or the proportion of tiny throats, whereas the reduction of intergranular pores play a clearly important role in the descending permeability (Figure 16a,b,d).The multidimensional correlation functions of T1 and T2, which were used to determine the characteristics of PSD, are of great interest.Appendix A presents the T1-T2 spectrum for specimens of tight sandstones; the bright red represents high amplitude, while the light red represents low amplitude.Appendix A illustrates that in either good or poor reservoir qualities specimens, the range of T1 and T2 is approximately 3 to 4 orders of magnitude because of the wide range of PSD, and the figures showed a wide distribution before centrifugation and a relatively narrow ridge-like peak in the center of the spectra after centrifugation.The bumps along the axis are likely artifacts due to noise.According to Song et al. (2002) [55], the peaks with small T1 and T2 are dominated by surface water whereas the long T1 and T2 peaks are contributed by surface and bulk water, and these are consistent with the notion that the free water (bulk water) that is in the center of the pores is easily removed after centrifugal processing [55].Besides, the zenith of the spectra before centrifugation moved from to characterize reservoir quality accurately; in addition, the movable fluid parameters are more closely aligned with permeability.Meanwhile, with decreasing reservoir quality and movable fluid saturation, chlorite decreases from 7% to 0%, while illite increases from 0% to 7%, indicating that illite exerts a significant impact on retaining water, whereas chlorite can retard compaction (Figure 16ac).Besides, reduction rate of reservoir quality with movable fluid saturation tends to decrease in tight sandstones with intergranular pores-dominated reservoirs to those with clay minerals, which indicates that movable fluid parameters are sensitive to a higher content of clay minerals or the proportion of tiny throats, whereas the reduction of intergranular pores play a clearly important role in the descending permeability (Figure 16a,b,d).

T1-T2 Correlation Spectra and Their Effects on Reservoir Quality
The multidimensional correlation functions of T1 and T2, which were used to determine the characteristics of PSD, are of great interest.Appendix A presents the T1-T2 spectrum for specimens of tight sandstones; the bright red represents high amplitude, while the light red represents low amplitude.Appendix A illustrates that in either good or poor reservoir qualities specimens, the range of T1 and T2 is approximately 3 to 4 orders of magnitude because of the wide range of PSD, and the figures showed a wide distribution before centrifugation and a relatively narrow ridge-like peak in the center of the spectra after centrifugation.The bumps along the axis are likely artifacts due to noise.According to Song et al. (2002) [55], the peaks with small T1 and T2 are dominated by surface water whereas the long T1 and T2 peaks are contributed by surface and bulk water, and these are consistent with the notion that the free water (bulk water) that is in the center of the pores is easily removed after centrifugal processing [55].Besides, the zenith of the spectra before centrifugation moved from

T 1 -T 2 Correlation Spectra and Their Effects on Reservoir Quality
The multidimensional correlation functions of T 1 and T 2 , which were used to determine the characteristics of PSD, are of great interest.Appendix A presents the T 1 -T 2 spectrum for specimens of tight sandstones; the bright red represents high amplitude, while the light red represents low amplitude.Appendix A illustrates that in either good or poor reservoir qualities specimens, the range of T 1 and T 2 is approximately 3 to 4 orders of magnitude because of the wide range of PSD, and the figures showed a wide distribution before centrifugation and a relatively narrow ridge-like peak in the center of the spectra after centrifugation.The bumps along the axis are likely artifacts due to noise.According to Song et al. (2002) [55], the peaks with small T 1 and T 2 are dominated by surface water whereas the long T 1 and T 2 peaks are contributed by surface and bulk water, and these are consistent with the notion that the free water (bulk water) that is in the center of the pores is easily removed after centrifugal processing [55].Besides, the zenith of the spectra before centrifugation moved from relatively high to low relaxation time in irreducible water conditions, suggesting that the aqueous phases have low mobility in tiny pores.Moreover, the ridges are more closely parallel to the line T 1 = T 2 , meaning that they have similar T 1 and T 2 ratio corresponding to similar properties (Appendix A).
The effect of reservoir quality is further investigated.For the specimen with greatest reservoir quality (11#) (Appendix A (A,B)), the proportion of bulk water that corresponds to larger pores is relatively high, compared to poor reservoir qualities specimens.This manifests that mill metric voids play a predominant role in optimization of reservoir quality in tight sandstone reservoirs.With permeability decreasing, the reduction rate of the water covering spectra areas of all the specimens, except specimen 5#, are decline, suggesting that the pores that allowed water to immigrate freely (bulk water) were few in the low permeability specimen, even in trace amounts in extremely tight sandstones (Appendix A (K,L)).Specimen 5# is characterized as low porosity and high permeability due the presence of microcracks (Figure 2l), which may lead to the wide range of T 2 and narrow range of T 1 before centrifugal processing (Figures 9a and 10b, Appendix A (C,D)).However, the mechanism of this phenomenon needs further investigation.Besides, for good reservoir quality specimens, peaks and shoulders are continuously distributed, whereas poor ones show separate distributions before centrifugation, revealing that the pores are relatively homogeneous, distributed in tight sandstones with high physical properties (Table 1, Appendix A (A-D)).After centrifugation, the spectra tend to show discrete distributions because the heterogeneity of the PSD caused some bulk water, which is hard to remove, especially for medium permeability specimens (Table 1, Appendix A (E-H)).For the extremely tight sandstone, however, the PSD shows a homogeneous distribution that results in the continuous distributions of T 1 -T 2 correlation spectra, because the very tiny pores generally show continuous arrangement (Table 1, Appendix A (I-L)).

NMR Imaging for Saturated and Irreducible Water Distributions
The images show distinct difference in simulated water saturation under fully saturated and irreducible water conditions; the red color represents oil and dark blue denotes deionized water, and the surrounding blue color may be caused by signal noise (Appendix B).The water saturation of the good reservoir quality specimens is always higher than that of poor reservoir quality, implying that brine can enter the majority of pore spectrum of high permeability specimens, whereas the low permeability specimens are dominated by tight section.Due to the hydrophilia of illite and I/S minerals [68] (low permeability specimens, which contain relatively abundant illite and I/S minerals), the change of water saturation before and after centrifugal processing is minor (Table 1, Appendix B).Besides, the percentage of high-saturated sections in high permeability specimens before and after centrifugal processing decline faster than that of the low permeability samples (slightly decrease) (Table 1, Appendix B).This implies that movable fluid mainly comes from a contribution of intergranular pores, while in pores.

Effects of Movable Pore Radius Lower Limit on Reservoir Quality
During NMR centrifugal processing, the specimens underwent different movable fluid loss, predominantly by centrifugal forces with minor gravity [65].In terms of the PSD being converted from the T 2 spectra and RCMI data, a series of T 2 cutoff values, which is defined by the ratio of T 2 spectra before and after centrifugal processing with various centrifugal forces, were often used to determine the lower limit of movable pore radius [11,32].In the case of high permeability specimens, a sharp decrease of macropores was observed corresponding to high T 2 relaxation time first, followed by a moderate and uniform reduction (Appendix A (A,B)).The correlation between lower limits of movable pore radius and movable fluid porosity show a distinct descending at the initial stage, followed by a moderate decrease (Figure 17).We can thus infer that the intergranular pores, which correspond to a larger pore radius, are the main result of the reduction of the lower limits of pore radius in higher permeability sandstones, and dissolution pores and intercrystalline pores are less important.For specimens with medium permeability, the amplitude of T 2 relaxation time drop uniformly (Appendix A (C,D)); however, the lower limits of movable pore radius show different changes: a distinct change in MFS and MFP over a small range of pore radius lower limits occurred in 14#, whereas the variance of 12 # is opposite (Figure 17), meaning that although these two specimens have similar reservoir quality, the percolation capacity of 14# is mainly dominated by heterogeneous pore structures with relatively low porosity, whereas the heterogeneity of the pore-throats in specimen 12# limits the improvement of permeability, even with relatively high porosity.Specimens with lower permeability resulted in no marked descending of T 2 and lower limits of movable pore radius (Appendix A (E,F)); the length of the curves for the tight specimen are generally short (except the first stage of 4# due to minor amounts of intergranular pores), indicating that pores in aqueous phases were hard to enter and were removed in poor reservoir quality sandstones (Figure 17).These results manifest that because the tight sandstones are highly heterogeneous, during the centrifugation process, rates of water saturation reduction in macropores are faster than those for tiny pores.For specimens with intergranular pores, the movable fluid advances, exempted from the larger pores.
important.For specimens with medium permeability, the amplitude of T2 relaxation time drop uniformly (Appendix C (C,D)); however, the lower limits of movable pore radius show different changes: a distinct change in MFS and MFP over a small range of pore radius lower limits occurred in 14#, whereas the variance of 12 # is opposite (Figure 17), meaning that although these two specimens have similar reservoir quality, the percolation capacity of 14# is mainly dominated by heterogeneous pore structures with relatively low porosity, whereas the heterogeneity of the porethroats in specimen 12# limits the improvement of permeability, even with relatively high porosity.Specimens with lower permeability resulted in no marked descending of T2 and lower limits of movable pore radius (Appendix C (E,F)); the length of the curves for the tight specimen are generally short (except the first stage of 4# due to minor amounts of intergranular pores), indicating that pores in aqueous phases were hard to enter and were removed in poor reservoir quality sandstones (Figure 17).These results manifest that because the tight sandstones are highly heterogeneous, during the centrifugation process, rates of water saturation reduction in macropores are faster than those for tiny pores.For specimens with intergranular pores, the movable fluid advances, exempted from the larger pores.

Conclusions
A combination of TS, SEM, X-ray diffraction, PCMI, RCMI, CT, and NMR tests were conducted on typical tight sandstone specimens from Ordos Basin to investigate pore structures and movable fluid traits.The following observations were derived: (1) Intergranular pores, dissolution pores, intercrystalline pores, and microfractures are four main kinds of pores in the research area.The PCMI-derived capillary pressure curves can be grouped into four types, and the RCMI-derived capillary pressure curves are divided into pore-dominated, throat-dominated and hybrid regions.Good pore-throat connectivity plays a vital role in enhancing physical properties based on CT images.The amplitude and incremental porosity distribution become small and narrow after centrifugal processing using NMR apparatus.

Conclusions
A combination of TS, SEM, X-ray diffraction, PCMI, RCMI, CT, and NMR tests were conducted on typical tight sandstone specimens from Ordos Basin to investigate pore structures and movable fluid traits.The following observations were derived: (1) Intergranular pores, dissolution pores, intercrystalline pores, and microfractures are four main kinds of pores in the research area.The PCMI-derived capillary pressure curves can be grouped into four types, and the RCMI-derived capillary pressure curves are divided into pore-dominated, throat-dominated and hybrid regions.Good pore-throat connectivity plays a vital role in enhancing physical properties based on CT images.The amplitude and incremental porosity distribution become small and narrow after centrifugal processing using NMR apparatus.
(2) The PSDs derived by RCMI show discrepancy with the direct observations, and the improved PSDs can be obtained by combining PCMI, RCMI and NMR and then reconstructing cumulative PSDs, which is consistent with the observation results of TS and SEM images.
(3) On the basis of Coates model, the linear regression results show that 9100 r•min −1 is the optimum rotation speed in the research area to determine appropriate T 2 cutoff values.MFP can characterize reservoir quality accurately, and a high proportion of chlorite with relatively abundant intergranular pores could improve MFS and reservoir quality.
(4) Bulk water makes great contributions to the movable fluid, whereas surface water is hard to be removed.Movable fluid mainly comes from the contribution of intergranular pores, while irreducible water is mainly stranded in intercrystalline pores, and the movable fluid could advance, exempted from the larger pores.
(5) Pore structures controls the lower limits of movable pore radius and hence affects reservoir quality.In sandstones with relatively high permeability, abundant intergranular pores can reduce the lower limits of pore radius; however, sandstones with complex pore structures have relatively poor reservoir qualities.

Figure 4 .
Figure 4. TS observations (net region with red rim represent the area that mercury occupies at this stage) and multitype model (the full line represents mercury fulfilling the network, the dashed line represents a wetting phase fulfilling the network and the line thickness represent the diameters of pore throats): (a) TS observation shows intergranular pores filled with mercury at the outset; (b) start of intrusion shows intergranular pores filled with mercury on the basis of a stick model; (c) TS observation show mercury invading the dissolution pores; (d) intrusion in this step shows that dissolution pores are occupied based on the stick model; (e) TS observation reveals that all connected pore throats are filled by mercury; (f) The end of invasion based on a stick model, mercury has occupied all the pore throats; (g) TS observation show that the pores among clay aggregates are filled with mercury; (h) The end of invasion based on a tree-like model, corresponding to the intercrystalline pores that are filled with mercury; (i) TS observation depicts the beginning of mercury withdrawal, reflecting mercury leaving the smallest pores; (j) Outset of mercury extrusion based on the tree-like model, intercrystalline pores did not leave a residual phase; (k) TS observation showing that mercury remains in weakly connected pores; (l) End of withdrawal, where the mercury remains in the tiny

Figure 4 .
Figure 4. TS observations (net region with red rim represent the area that mercury occupies at this stage) and multitype model (the full line represents mercury fulfilling the network, the dashed line represents a wetting phase fulfilling the network and the line thickness represent the diameters of pore throats): (a) TS observation shows intergranular pores filled with mercury at the outset; (b) start of intrusion shows intergranular pores filled with mercury on the basis of a stick model; (c) TS observation show mercury invading the dissolution pores; (d) intrusion in this step shows that dissolution pores are occupied based on the stick model; (e) TS observation reveals that all connected pore throats are filled by mercury; (f) The end of invasion based on a stick model, mercury has occupied all the pore throats; (g) TS observation show that the pores among clay aggregates are filled with mercury; (h) The end of invasion based on a tree-like model, corresponding to the intercrystalline pores that are filled with mercury; (i) TS observation depicts the beginning of mercury withdrawal, reflecting mercury leaving the smallest pores; (j) Outset of mercury extrusion based on the tree-like model, intercrystalline pores did not leave a residual phase; (k) TS observation showing that mercury remains in weakly connected pores; (l) End of withdrawal, where the mercury remains in the tiny pores that are less unconnected to the outside, on the basis of a stick model.IGP-intergranular pore; DP-dissolution pore; ICP-intercrystalline pore.

Figure 7 .
Figure 7. Pore-throat size distribution by RCMI of typical tight sandstones specimens: (a) pore size distribution; (b) throat size distribution; (c) pore throat radius ratio distribution.

Figure 8 .
Figure 8. (A,D) Reconstructed 3D CT images and volume of (B,E) resolved pore-throat space and (C,F) micro-CT derived pore radius distribution of tight sandstone specimen of 11# and 12#, respectively.The pore-throat clusters are shown in different colors.

Figure 7 .
Figure 7. Pore-throat size distribution by RCMI of typical tight sandstones specimens: (a) pore size distribution; (b) throat size distribution; (c) pore throat radius ratio distribution.

Figure 7 .
Figure 7. Pore-throat size distribution by RCMI of typical tight sandstones specimens: (a) pore size distribution; (b) throat size distribution; (c) pore throat radius ratio distribution.

Figure 8 .
Figure 8. (A,D) Reconstructed 3D CT images and volume of (B,E) resolved pore-throat space and (C,F) micro-CT derived pore radius distribution of tight sandstone specimen of 11# and 12#, respectively.The pore-throat clusters are shown in different colors.
are the result of optimum values of the smoothening factors arising from the low error and unstable nature of the problem [62,63].

Figure 8 .
Figure 8. (A,D) Reconstructed 3D CT images and volume of (B,E) resolved pore-throat space and (C,F) micro-CT derived pore radius distribution of tight sandstone specimen of 11# and 12#, respectively.The pore-throat clusters are shown in different colors.

Figure 11 .
Figure 11.Methods for determining the improved-NMR PSD by the integration of NMR and RCMI, taking B116 specimen in the research area as an example: (a) data selection; (b) relationship between T2 relaxation time and pore-throat radius.

Figure 11 .
Figure 11.Methods for determining the improved-NMR PSD by the integration of NMR and RCMI, taking B116 specimen in the research area as an example: (a) data selection; (b) relationship between T 2 relaxation time and pore-throat radius.

Figure 12 .
Figure 12.Full-range PSD by the integration of improved-NMR PSD and PCMI: (a) sample B116; (b) the PSDs of all specimens measured in this research.

Figure 12 .
Figure 12.Full-range PSD by the integration of improved-NMR PSD and PCMI: (a) sample B116; (b) the PSDs of all specimens measured in this research.

Figure 12 .
Figure 12.Full-range PSD by the integration of improved-NMR PSD and PCMI: (a) sample B116; (b) the PSDs of all specimens measured in this research.

Figure 16 .
Figure 16.Graphs illustrating the impacts of pore types, movable fluid saturation and clay mineral contents on reservoir quality.(a) reservoir quality; (b) clay minerals; (c) movable fluid saturations (MFS); (d) thin section images record the pore types.

Figure 16 .
Figure 16.Graphs illustrating the impacts of pore types, movable fluid saturation and clay mineral contents on reservoir quality.(a) reservoir quality; (b) clay minerals; (c) movable fluid saturations (MFS); (d) thin section images record the pore types.

Figure 16 .
Figure 16.Graphs illustrating the impacts of pore types, movable fluid saturation and clay mineral contents on reservoir quality.(a) reservoir quality; (b) clay minerals; (c) movable fluid saturations (MFS); (d) thin section images record the pore types.

Figure 17 .
Figure 17.The plot of pore radius lower limits versus (a) movable fluid saturation and (b) movable fluid porosity.

Figure 17 .
Figure 17.The plot of pore radius lower limits versus (a) movable fluid saturation and (b) movable fluid porosity.

Table 1 .
Physical properties and petrological composition of Yanchang formation tight sandstones 1 .

Table 2 .
The parameters for pore structure characterization from PCMI 2 .

Table 3 .
The parameters for pore structure characterization from RCMI.3 3 P d -threshold pressure; r max -maximum pore radius; P 50 -median pressure; r 50 -median pore radius; r t -average throat radius; r p -average pore radius; η-pore throat radius ratio; r m -mainstream throat radius; δ-sorting coefficient; S t -throat maximum mercury intrusion saturation; S p -pore maximum mercury intrusion saturation; S Hgmax -total maximum mercury intrusion saturation.

Table 4 .
The parameters for movable fluid parameters from NMR.4

3. Controls of Movable Fluid Traits on Reservoir Quality 4
.3.1.Movable Fluid Parameters and Their Effects on Reservoir Quality

Table 4 .
The parameters for movable fluid parameters from NMR.4