Pore structure characterization of North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury intrusion

Small-angle and ultra-small-angle neutron scattering (SANS and USANS), low-pressure adsorption (N2 and CO2), and high-pressure mercury intrusion measurements were performed on a suite of North American shale reservoir samples providing the first ever comparison of all these techniques for characterizing the complex pore structure of shales. The techniques were used to gain insight into the nature of the pore structure including pore geometry, pore size distribution and accessible versus inaccessible porosity. Reservoir samples for analysis were taken from currently-active shale gas plays including the Barnett, Marcellus, Haynesville, Eagle Ford, Woodford, Muskwa, and Duvernay shales. Low-pressure adsorption revealed strong differences in BET surface area and pore volumes for the sample suite, consistent with variability in composition of the samples. The combination of CO2 and N2 adsorption data allowed pore size distributions to be created for micro–meso–macroporosity up to a limit of �1000 Å. Pore size distributions are either uni- or multi-modal. The adsorption-derived pore size distributions for some samples are inconsistent with mercury intrusion data, likely owing to a combination of grain compression during high-pressure intrusion, and the fact that mercury intrusion yields information about pore throat rather than pore body distributions. SANS/USANS scattering data indicate a fractal geometry (power-law scattering) for a wide range of pore sizes and provide evidence that nanometer-scale spatial ordering occurs in lower mesopore–micropore range for some samples, which may be associated with inter-layer spacing in clay minerals. SANS/USANS pore radius distributions were converted to pore volume distributions for direct comparison with adsorption data. For the overlap region between the two methods, the agreement is quite good. Accessible porosity in the pore size (radius) range 5 nm–10 lm was determined for a Barnett shale sample using the contrast matching method with pressurized deuterated methane fluid. The results demonstrate that accessible porosity is pore-size dependent.

In situ SANS study of the surface and pore filling adsorption of subcritical and supercritical CO2 in porous silica as a function of pore size

We applied small-angle neutron scattering (SANS) and ultra small-angle neutron scattering (USANS) to monitor evolution of the CO2 adsorption in porous silica as a function of CO2 pressure and temperature in pores of different sizes. The range of pressures (0 < P < 345 bar) and temperatures (T=18 OC, 35 OC and 60 OC) corresponded to subcritical, near critical and supercritical conditions of bulk fluid. We observed that the adsorption behavior of CO2 is fundamentally different in large and small pores with the sizes D > 100 Å and D < 30 Å, respectively. Scattering data from large pores indicate formation of a dense adsorbed film of CO2 on pore walls with the liquid-like density (ρCO2)ads≈0.8 g/cm3. The adsorbed film coexists with unadsorbed fluid in the inner pore volume. The density of unadsorbed fluid in large pores is temperature and pressure dependent: it is initially lower than (ρCO2)ads and gradually approaches it with pressure. In small pores compressed CO2 gas completely fills the pore volume. At the lowest pressures of the order of 10 bar and T=18 OC, the fluid density in smallest pores available in the matrix with D ~ 10 Å exceeds bulk fluid density by a factor of ~ 8. As pressure increases, progressively larger pores become filled with the condensed CO2. Fluid densification is only observed in pores with sizes less than ~ 25 – 30 Å. As the density of the invading fluid reaches (ρCO2)bulk~ 0.8 g/cm3, pores of all sizes become uniformly filled with CO2 and the confinement effects disappear. At higher densities the fluid in small pores appears to follow the equation of state of bulk CO2 although there is an indication that the fluid density in the inner volume of large pores may exceed the density of the adsorbed layer. The equivalent internal pressure (Pint) in the smallest pores exceeds the external pressure (Pext) by a factor of ~ 5 for both sub- and supercritical CO2. Pint gradually approaches Pext as D → 25 – 30 Å and is independent of temperature in the studied range of 18 OC ≤ T ≤ 60 OC. The obtained results demonstrate certain similarity as well as differences between adsorption of subcritical and supercritical CO2 in disordered porous silica. High pressure small angle scattering experiments open new opportunities for in situ studies of the fluid adsorption in porous media of interest to CO2 sequestration, energy storage, and heterogeneous catalysis.

Characterization of tight gas reservoir pore structure using USANS/SANS and gas adsorption analysis

Small-angle and ultra-small-angle neutron scattering (SANS and USANS) measurements were performed on samples from the Triassic Montney tight gas reservoir in Western Canada in order to determine the applicability of these techniques for characterizing the full pore size spectrum and to gain insight into the nature of the pore structure and its control on permeability. The subject tight gas reservoir consists of a finely laminated siltstone sequence; extensive cementation and moderate clay content are the primary causes of low permeability. SANS/USANS experiments run at ambient pressure and temperature conditions on lithologically-diverse sub-samples of three core plugs demonstrated that a broad pore size distribution could be interpreted from the data. Two interpretation methods were used to evaluate total porosity, pore size distribution and surface area and the results were compared to independent estimates derived from helium porosimetry (connected porosity) and low-pressure N2 and CO2 adsorption (accessible surface area and pore size distribution). The pore structure of the three samples as interpreted from SANS/USANS is fairly uniform, with small differences in the small-pore range (<2000 Å), possibly related to differences in degree of cementation, and mineralogy, in particular clay content. Total porosity interpreted from USANS/SANS is similar to (but systematically higher than) helium porosities measured on the whole core plug. Both methods were used to estimate the percentage of open porosity expressed here as a ratio of connected porosity, as established from helium adsorption, to the total porosity, as estimated from SANS/USANS techniques. Open porosity appears to control permeability (determined using pressure and pulse-decay techniques), with the highest permeability sample also having the highest percentage of open porosity. Surface area, as calculated from low-pressure N2 and CO2 adsorption, is significantly less than surface area estimates from SANS/USANS, which is due in part to limited accessibility of the gases to all pores. The similarity between N2 and CO2-accessible surface area suggests an absence of microporosity in these samples, which is in agreement with SANS analysis. A core gamma ray profile run on the same core from which the core plug samples were taken correlates to profile permeability measurements run on the slabbed core. This correlation is related to clay content, which possibly controls the percentage of open porosity. Continued study of these effects will prove useful in log-core calibration efforts for tight gas.

A USANS/SANS study of the accessibility of pores in the barnett shale to methane and water

Shale is an increasingly important source of natural gas in the United States. The gas is held in fine pores that need to be accessed by horizontal drilling and hydrofracturing techniques. Understanding the nature of the pores may provide clues to making gas extraction more efficient. We have investigated two Mississippian Barnett Shale samples, combining small-angle neutron scattering (SANS) and ultrasmall-angle neutron scattering (USANS) to determine the pore size distribution of the shale over the size range 10 nm to 10 μm. By adding deuterated methane (CD4) and, separately, deuterated water (D2O) to the shale, we have identified the fraction of pores that are accessible to these compounds over this size range. The total pore size distribution is essentially identical for the two samples. At pore sizes >250 nm, >85% of the pores in both samples are accessible to both CD4 and D2O. However, differences in accessibility to CD4 are observed in the smaller pore sizes (∼25 nm). In one sample, CD4 penetrated the smallest pores as effectively as it did the larger ones. In the other sample, less than 70% of the smallest pores (<25 nm) were accessible to CD4, but they were still largely penetrable by water, suggesting that small-scale heterogeneities in methane accessibility occur in the shale samples even though the total porosity does not differ. An additional study investigating the dependence of scattered intensity with pressure of CD4 allows for an accurate estimation of the pressure at which the scattered intensity is at a minimum. This study provides information about the composition of the material immediately surrounding the pores. Most of the accessible (open) pores in the 25 nm size range can be associated with either mineral matter or high reflectance organic material. However, a complementary scanning electron microscopy investigation shows that most of the pores in these shale samples are contained in the organic components. The neutron scattering results indicate that the pores are not equally proportioned in the different constituents within the shale. There is some indication from the SANS results that the composition of the pore-containing material varies with pore size; the pore size distribution associated with mineral matter is different from that associated with organic phases.

Pore size distribution and accessible pore size distribution in bituminous coals

The porosity and pore size distribution of coals determine many of their properties, from gas release to their behavior on carbonization, and yet most methods of determining pore size distribution can only examine a restricted size range. Even then, only accessible pores can be investigated with these methods. Small-angle neutron scattering (SANS) and ultra small-angle neutron scattering (USANS) are increasingly used to characterize the size distribution of all of the pores non-destructively. Here we have used USANS/SANS to examine 24 well-characterized bituminous and subbituminous coals: three from the eastern US, two from Poland, one from New Zealand and the rest from the Sydney and Bowen Basins in Eastern Australia, and determined the relationships of the scattering intensity corresponding to different pore sizes with other coal properties. The range of pore radii examinable with these techniques is 2.5nm to 7μm. We confirm that there is a wide range of pore sizes in coal. The pore size distribution was found to be strongly affected by both rank and type (expressed as either hydrogen or vitrinite content) in the size range 250nm to 7μm and 5 to 10nm, but weakly in intermediate regions. The results suggest that different mechanisms control coal porosity on different scales. Contrast-matching USANS and SANS were also used to determine the size distribution of the fraction of the pores in these coals that are inaccessible to deuterated methane, CD4, at ambient temperature. In some coals most of the small (~10nm) pores were found to be inaccessible to CD4 on the time scale of the measurement (~30min–16h). This inaccessibility suggests that in these coals a considerable fraction of inherent methane may be trapped for extended periods of time, thus reducing the effectiveness of methane release from (or sorption by) these coals. Although the number of small pores was less in higher rank coals, the fraction of total pores that was inaccessible was not rank dependent. In the Australian coals, at the 10nm to 50nm size scales the pores in inertinites appeared to be completely accessible to CD4, whereas the pores in the vitrinite were about 75% inaccessible. Unlike the results for total porosity that showed no regional effects on relationships between porosity and coal properties, clear regional differences in the relationships between fraction of closed porosity and coal properties were found. The 10 to 50nm-sized pores of inertinites of the US and Polish coals examined appeared less accessible to methane than those of the inertinites of Australian coals. This difference in pore accessibility in inertinites may explain why empirical relationships between fluidity and coking properties developed using Carboniferous coals do not apply to Australian coals.

Accessibility of pores in coal to methane and carbon dioxide

Fluid–solid interactions in natural and engineered porous solids underlie a variety of technological processes, including geological storage of anthropogenic greenhouse gases, enhanced coal bed methane recovery, membrane separation, and heterogeneous catalysis. The size, distribution and interconnectivity of pores, the chemical and physical properties of the solid and fluid phases collectively dictate how fluid molecules migrate into and through the micro- and meso-porous media, adsorb and ultimately react with the solid surfaces. Due to the high penetration power and relatively short wavelength of neutrons, smallangle neutron scattering (SANS) as well as ultra small-angle scattering (USANS) techniques are ideally suited for assessing the phase behavior of confined fluids under pressure as well as for evaluating the total porosity in engineered and natural porous systems including coal. Here we demonstrate that SANS and USANS can be also used for determining the fraction of the pore volume that is actually accessible to fluids as a function of pore sizes and study the fraction of inaccessible pores as a function of pore size in three coals from the Illinois Basin (USA) and Bowen Basin (Australia). Experiments were performed at CO2 and methane pressures up to 780 bar, including pressures corresponding to zero average contrast condition (ZAC), which is the pressure where no scattering from the accessible pores occurs. Scattering curves at the ZAC were compared with the scattering from same coals under vacuum and analysed using a newly developed approach that shows that the volume fraction of accessible pores in these coals varies between �90% in the macropore region to �30% in the mesopore region and the variation is distinctive for each of the examined coals. The developed methodology may be also applied for assessing the volume of accessible pores in other natural underground formations of interest for CO2 sequestration, such as saline aquifers as well as for estimating closed porosity in engineered porous solids of technological importance.

Monitoring phase behavior of sub- and supercritical CO2 confined in porous fractal silica with 85% porosity

Phase behavior of CO2 confined in porous fractal silica with volume fraction of SiO2 φs = 0.15 was investigated using small-angle neutron scattering (SANS) and ultrasmall-angle neutron scattering (USANS) techniques. The range of fluid densities (0<(FCO2)bulk<0.977 g/cm3) and temperatures (T=22 °C, 35 and 60 °C) corresponded to gaseous, liquid, near critical and supercritical conditions of the bulk fluid. The results revealed formation of a dense adsorbed phase in small pores with sizes D<40 A° at all temperatures. At low pressure (P <55 bar, (FCO2)bulk <0.2 g/cm3) the average fluid density in pores may exceed the density of bulk fluid by a factor up to 6.5 at T=22 °C. This “enrichment factor” gradually decreases with temperature, however significant fluid densification in small pores still exists at temperature T=60°C, i.e., far above the liquid-gas critical temperature of bulk CO2 (TC=31.1 °C). Larger pores are only partially filled with liquid-like adsorbed layer which coexists with unadsorbed fluid in the pore core. With increasing pressure, all pores become uniformly filled with the fluid, showing no measurable enrichment or depletion of the porous matrix with CO2.

Stability of the bituminous coal microstructure upon exposure to high pressures of helium

Small-angle neutron scattering (SANS) and ultra-small-angle neutron scattering (USANS) measurements of the structure of two Australian bituminous coals (particle size of 1-0.5 mm) before, during, and after exposure to 155 bar of helium were made to identify any effects of pressure alone on the pore size distribution of coal and any irreversible effects upon exposure to high pressures of helium in the pore size range from 3 nm to 10 μm. No irreversible effects upon exposure were identified for any pore size. No effects of pressure on pore size distribution were observed, except for a small effect at a pore size of about 2 μm for one coal. This study provides a convenient baseline for SANS and USANS investigations on sorption of gases at elevated pressures on coals, by distinguishing between the effect of pressure alone on coal pore size distribution and against the effect of the gas to be investigated.

Transition from disc to rod-like shape of 16-3-16 dimeric micelles in aqueous solutions

Small-angle neutron scattering (SANS) measurements from bis-cationic C16H33N+(CH3)(2)-(CH2)(3)-N+ (CH3)(2)C16H33 2Br(-) dimeric surfactant, referred to as 16-3-16, at different concentrations and temperatures, are reported. It is seen that micelles are disc-like for concentrations C = 2.5 and 10 mM at temperature T = 30 degrees C. At low concentration C = 0.5 mM micelles are rod-like. Similarly, there is a disc to rod-like transition of micelles on increasing the temperature. For C = 2.5 mM, micelles are rod-like at T = 45 and 70 degrees C.

Structure and Rheology of Cationic Molecular Hydrogels of Quinuclidine Grafted Bile Salts. Influence of the Ionic Strength and Counter-Ion type

Quinuclidine grafted cationic bile salts are forming salted hydrogels. An extensive investigation of the effect of the electrolyte and counterions on the gelation has been envisaged. The special interest of the quinuclidine grafted bile salt is due to its broader experimental range of gelation to study the effect of electrolyte. Rheological features of the hydrogels are typical of enthalpic networks exhibiting a scaling law of the elastic shear modulus with the concentration (scaling exponent 2.2) modeling cellular solids in which the bending modulus is the dominant parameter. The addition of monovalent salt (NaCl) favors the formation of gels in a first range (0.00117 g cm-3 (0.02 M) < TNaCl < 0.04675 g cm-3 (0.8 M)). At larger salt concentrations, the gels become more heterogeneous with nodal zones in the micron scale. Small-angle neutron scattering experiments have been used to characterize the rigid fibers ( ≈ 68 Å) and the nodal zones. Stress sweep and creeprecovery measurements are used to relate the lack of linear viscoelastic domain to a mechanism of disentanglement of the fibers from their associations into fagots. The electrostatic interactions can be screened by addition of salt to induce a progressive evolution toward flocculation. SEM, UV absorbance, and SAXS study of the Bragg peak at large Q-values complete the investigation.

Structure and properties of two component hydrogels comprising lithocholic acid and organic amines

We demonstrate the aptitude of supramolecular hydrogel formation using simple bile acid such as lithocholic acid in aqueous solution in the presence of various dimeric or oligomeric amines. By variation of the choice of the amines in such mixtures the gelation properties could be modulated. However, the replacement of lithocholic acid (LCA) by cholic acid or deoxycholic acid resulted in no hydrogel formation. FT-IR studies confirm that the carboxylate and ammonium residues of the two components are involved in the salt (ion-pair) formation. This promotes further assembly of the components reinforced by a continuous hydrogen bonded network leading to gelation. Electron microscopy shows the morphology of the internal organization of gels of two component systems which also depends significantly on the amine part. Variation of the amine component from the simple 1,2-ethanediamine (EDA) to oligomeric amines in such gels of lithocholic acid changes the morphology of the assembly from long one-dimensional nanotubes to three-dimensional complex structures. Single crystal X-ray diffraction analysis with one of the amine-LCA complexes suggested the motif of fiber formation where the amines interact with the carboxylate and hydroxyl moieties through electrostatic forces and hydrogen bonding. From small angle neutron scattering study, it becomes clear that the weak gel from LCA-EDA shows scattering oscillation due to the presence of non-interacting nanotubules while for gels of LCA with oligomeric amines the individual fibers come together to form complex three-dimensional organizations of higher length scale. The rheological properties of this class of two component system provide clear evidence that the flow behavior can be modulated varying the acid-amine ratio.

When using small samples to evaluate hydrocarbon reservoirs, proceed with caution

Small-angle neutron scattering (SANS) and ultra-small angle neutron scattering (USANS) with contrast matching techniques (Melnichenko and others, 2012) were used to investigate size distribution and gas accessibility in pores in an approximately 10.6 cm long Mississippian Barnett Shale butt core from the Fort Worth Basin, Texas, USA. SANS and USANS measurements record scattering from all pores, both open and closed, in the size range 10nm - ~10 μ. The techniques can also be used to determine the material that contains pores and the number of pores as a function of size. By injecting deuterated methane gas (CD4) at contrast matching pressure it is possible to distinguish which pores are accessible, or open, to fluids and which ones are not.

Structure of polyamidoamide dendrimers up to limiting generations: A mesoscale description

The polyamidoamide (PAMAM) class of dendrimers was one of the first dendrimers synthesized by Tomalia and co-workers at Dow. Since its discovery the PAMAMs have stimulated many discussions on the structure and dynamics of such hyperbranched polymers. Many questions remain open because the huge conformation disorder combined with very similar local symmetries have made it difficult to characterize experimentally at the atomistic level the structure and dynamics of PAMAM dendrimers. The higher generation dendrimers have also been difficult to characterize computationally because of the large size (294852 atoms for generation 11) and the huge number of conformations. To help provide a practical means of atomistic computational studies, we have developed an atomistically informed coarse-grained description for the PAMAM dendrimer. We find that a two-bead per monomer representation retains the accuracy of atomistic simulations for predicting size and conformational complexity, while reducing the degrees of freedom by tenfold. This mesoscale description has allowed us to study the structural properties of PAMAM dendrimer up to generation 11 for time scale of up to several nanoseconds. The gross properties such as the radius of gyration compare very well with those from full atomistic simulation and with available small angle x-ray experiment and small angle neutron scattering data. The radial monomer density shows very similar behavior with those obtained from the fully atomistic simulation. Our approach to deriving the coarse-grain model is general and straightforward to apply to other classes of dendrimers.

Phenyl-Ring-Bearing Cationic Surfactants: Effect of Ring Location on the Micellar Structure

A series of isomeric cationic surfactants (S1-S5) bearing a long alkyl chain that carries a 1,4-phenylene unit and a trimethyl ammonium headgroup was synthesized; the location of the phenyl ring within the alkyl tail was varied in an effort to understand its influence on the amphiphilic properties of the surfactants. The cmc's of the surfactants were estimated using ionic conductivity measurements and isothermal calorimetric titrations (ITC); the values obtained by the two methods were found to be in excellent agreement. The ITC measurements provided additional insight into the various thermodynamic parameters associated with the micellization process. Although all five surfactants have exactly the same molecular formula, their micellar properties were seen to vary dramatically depending on the location of the phenyl ring; the cmc was seen to decrease by almost an order of magnitude when the phenyl ring was moved from the tail end (cmc of S1 is 23 mM) to the headgroup region (cmc of S5 is 3 mM). In all cases, the enthalpy of micellization was negative but the entropy of micellization was positive, suggesting that in all of these systems the formation of micelles is both enthalpically and entropically favored. As expected, the decrease in cmc values upon moving the phenyl ring from the tail end to he headgroup region is accompanied by an increase in the thermodynamic driving force (Delta G) for micellization. To understand further the differences in the micellar structure of these surfactants, small-angle neutron scattering (SANS) measurements were carried out; these measurements reveal that the aggregation number of the micelles increases as the cmc decreases. This increase in the aggregation number is also accompanied by an increase in the asphericity of the micellar aggregate and a decrease in the fractional charge. Geometric packing arguments are presented to account for these changes in aggregation behavior as a function of phenyl ring location.

Colloidal dispersions in a liquid crystalline medium

We have prepared stable colloidal suspensions in a lyotropic liquid crystal exhibiting an isotropic-nematic-lamellar phase sequence. Small angle neutron scattering (SANS) and dynamic light scattering (DLS) studies show the existence of attractive interparticle interactions in the nematic phase, which lead to a gas-liquid transition of the particles. The resulting liquid phase is weakly anisotropic. Further, the nematic-lamellar transition of the liquid crystal is found to be accompanied by a liquid-solid transition of the particles.