Temperature-Dependent Oxygen Effect on NMR D-[Formula: see text] Relaxation-Diffusion Correlation of n-Alkanes.

Nuclear magnetic resonance (NMR) diffusion-relaxation correlation experiments (D-[Formula: see text]) are widely used for the petrophysical characterisation of rocks saturated with petroleum fluids both in situ and for laboratory analyses. The encoding for both diffusion and relaxation offers increased fluid typing contrast by discriminating fluids based on their self-diffusion coefficients, while relaxation times provide information about the interaction of solid and fluid phases and associated confinement geometry (if NMR responses of pure fluids at particular temperature and pressure are known). Petrophysical interpretation of D-[Formula: see text] correlation maps is typically assisted by the "standard alkane line"-a relaxation-diffusion correlation valid for pure normal alkanes and their mixtures in the absence of restrictions to diffusing molecules and effects of internal gradients. This correlation assumes fluids are free from paramagnetic impurities. In situations where fluid samples cannot be maintained at air-free state the diffusion-relaxation response of fluids shift towards shorter relaxation times due to oxygen paramagnetic relaxation enhancement. Interpretation of such a response using the "standard alkane line" would be erroneous and is further complicated by the temperature-dependence of oxygen solubility for each component of the alkane mixture. We propose a diffusion-relaxation correlation suitable for interpretation of low-field NMR D-[Formula: see text] responses of normal alkanes and their mixtures saturating rocks over a broad temperature range, in equilibrium with atmospheric air. We review and where necessary revise existing viscosity-relaxation correlations. Findings are applied to diffusion-relaxation dependencies taking into account the temperature dependence of oxygen solubility and solvent vapour pressure. The effect is demonstrated on a partially saturated carbonate rock.

Lattice Boltzmann framework for accurate NMR simulation in porous media.

Nuclear magnetic resonance (NMR) responses of fluids saturating porous media arise from complex relaxation-diffusion dynamics of polarized spins. These constitute a sensitive probe of the microstructure and are described by the Bloch-Torrey equations. An NMR simulation framework based on an augmented lattice Boltzmann method aimed at the fine-scale resolution of nuclear polarization density is presented. The approach encapsulates the time evolution of the full magnetization vector and naturally incorporates the mechanisms of diffusional transport. Spin dephasing mechanisms are fully resolved at tomogram voxel scale to account for magnetic field inhomogeneity. The approach is validated against analytical solutions of spin-echo decays for simple pore geometries. An application to a nano-computed-tomography image of chalk with inhomogeneous internal fields yields T_{2} spectral measures in good agreement with experiment and illustrates the spatial pore-scale dynamics of net magnetization. Findings establish the feasibility of the framework for pure diffusion and present an approach vector to modeling the evolution of magnetization under flow conditions.

A comparison between the characteristics of a biochar-NPK granule and a commercial NPK granule for application in the soil.

Continual application of nitrogen (N), phosphorous (P) and potassium (K) fertilizer may not return a profit to farmers due to the costs of application and the loss of NPK from soil in various ways. Thus, a combination of NPK granule with a porous biochar (termed here as BNPK) appears to offer multiple benefits resulting from the excellent properties of biochar. Given the lack of information on the properties of NPK and BNPK fertilizers, it is necessary to investigate the characteristics of both to achieve a good understanding of why BNPK granule is superior to NPK granule. Therefore, this study aims to investigate the characteristics of a maize straw biochar mixed with NPK granule, before and after application to soil, and compare them to those for a commercial NPK granule. The BNPK granule, with a greater surface area and porosity, showed a higher capacity to store and donate electrons than the NPK granule. Relatively lower concentrations of Ca, P, K, Si and Mg were dissolved from the BNPK, indicating the ability of the BNPK granule to maintain these mineral elements and reduce dissolution rate. To study the nutrient storage mechanism of the BNPK granule in the soil, short- and long-term leaching experiments were conducted. During the experiments, organo-mineral clusters, comprising C, P, K, Si, Al and Fe, were formed on the surface and inside the biochar pores. However, BNPK was not effective in reducing N leaching, in the absence of plants, in a red chromosol soil.

Fast Fourier transform and support-shift techniques for pore-scale microstructure classification using additive morphological measures.

The Minkowski functionals, as the full set of additive morphological measures in three dimensions (3D) consisting of volume, surface area, mean curvature, and total curvature, can be calculated directly by evaluating the local contributions of vertices of a discrete structure. They are sensitive measures of microstructure, and for microstructures generated by a Boolean process, relate to their physical properties. In this work we introduce fast numerical techniques based on the additivity of the Minkowski functionals to derive fields of regional Minkowski measures over large regional support for large 3D data sets as generated, e.g., from x-ray tomography techniques. We demonstrate the application of these 3D feature fields to microstructure classification for a set of heterogeneous microstructures using a multivariate Gaussian mixture model and a thin-bedded sandstone. It is shown that for the case of a spatially heterogeneous Boolean process the internal boundaries of the generating process are recovered with high accuracy, while for the thin-bedded sandstone, compact partitions with clear layering are extracted.

Application of low-field, 1H/13C high-field solution and solid state NMR for characterisation of oil fractions responsible for wettability change in sandstones.

Asphaltene adsorption on solid surfaces is a standing problem in petroleum industry. It has an adverse effect on reservoir production and development by changing rock wettability, plugging pore throats, and affects oil transport through pipelines. Asphaltene chemistry constitutes important part of the ageing process as part of petrophysical studies and core analysis. The mechanisms and contribution of various oil components to adsorption processes is not fully understood. To investigate the kinetics of the ageing process and address the relative contribution of different oil components, we prepared three sets of sandstone core plugs aged in different oil mixtures over various time intervals. Cores were then re-saturated with decane to evaluate their wetting state using low-field NMR relaxometry by monitoring a change of surface relaxivity. Adsorbed deposits were then extracted from cores for solution-state NMR analysis. Their 1H and 1H-13C correlation spectra obtained using heteronuclear single quantum coherence (HSQC) technique were matched to spectra of four SARA (saturates, aromatics, resins and asphaltenes) components of oil mixtures to deduce components of deposits and inter-component interactions. We notice that wettability reversal of rock is inversely proportional to initial asphaltene concentration. Analysis of deposits reveals an increase in their aliphatic content over ageing time, which is accompanied by a change of the morphology of the pore space due to cluster aggregates forming a network. Results suggest that the ageing process in respect to the wetting state of rock samples consists of three distinctive stages: (i) an early-time period, when the fraction of most polar asphaltenes creates a discontinuous layer corresponding to mixed-wet state; (ii) an intermediate-time interval, at which the full grain coverage may be achieved (at favourable chemical environment) corresponding to strong oil-wetting; (iii) a late-time stage, where intense macro-aggregates accumulation occurs, changing the pore space integrity. It is likely asphaltene-aliphatic interactions leading to growth of sub-micron size macro-aggregates.

About the connectivity of dual-scale media based on micro-structure based regional analysis of NMR flow propagators.

The characterisation of heterogeneous porous media at multiple length scales typically requires the classification of structure at some scale to allow the calculation of effective transport properties at a scale relevant for macroscopic description. While such a classification may be derived from various imaging methods, a shortcoming is often the simultaneous characterisation of the connectivity between regions representing different micro-structure. In this work we combine NMR based flow propagators with the simulations performed on corresponding reconstructed structure, and relate the NMR measurements to their simulated global and local representations to study fluid transport locally and the exchange between micro- and macro-porous regions. This is achieved by carrying out detailed lattice Boltzmann simulations and random walk method to track the displacements of tracers in each kind of region. Using Euclidean distance maps (EDT) we analyse the fluid invasion to regions of different scale and relate it to the connectivity of the system. We demonstrate that numerical simulation has great flexibility in providing additional sensitivity to the inference of region-region connectivity.

Theoretical investigation of heterogeneous wettability in porous media using NMR.

It is highly important to understand the heterogeneous wettability properties of porous media for enhanced oil recovery (EOR). However, wettability measurements are still challenging in directly investigating the wettability of porous media. In this paper, we propose a multidimensional nuclear magnetic resonance (NMR) method and the concept of apparent contact angles to characterize the heterogeneous wettability of porous media. The apparent contact angle, which is determined by both the wetting surface coverage and the local wettability (wetting contact angles of each homogeneous wetting regions or wetting patches), is first introduced as an indicator of the heterogeneous wettability of porous media using the NMR method. For homogeneously wetting patches, the relaxation time ratio T1/T2 is employed to probe the local wettabiity of wetting patches. The T2 - D is introduced to obtain the wetting surface coverage using the effective relaxivity. Numerical simulations are conducted to validate this method.

Assessment of bone ingrowth into porous biomaterials using MICRO-CT.

The three-dimensional (3D) structure and architecture of biomaterial scaffolds play a critical role in bone formation as they affect the functionality of the tissue-engineered constructs. Assessment techniques for scaffold design and their efficacy in bone ingrowth studies require an ability to accurately quantify the 3D structure of the scaffold and an ability to visualize the bone regenerative processes within the scaffold structure. In this paper, a 3D micro-CT imaging and analysis study of bone ingrowth into tissue-engineered scaffold materials is described. Seven specimens are studied in this paper; a set of three specimens with a cellular structure, varying pore size and implant material, and a set of four scaffolds with two different scaffold designs investigated at early (4 weeks) and late (12 weeks) explantation times. The difficulty in accurately phase separating the multiple phases within a scaffold undergoing bone regeneration is first highlighted. A sophisticated three-phase segmentation approach is implemented to develop high-quality phase separation with minimal artifacts. A number of structural characteristics and bone ingrowth characteristics of the scaffolds are quantitatively measured on the phase separated images. Porosity, pore size distributions, pore constriction sizes, and pore topology are measured on the original pore phase of the scaffold volumes. The distribution of bone ingrowth into the scaffold pore volume is also measured. For early explanted specimens we observe that bone ingrowth occurs primarily at the periphery of the scaffold with a constant decrease in bone mineralization into the scaffold volume. Pore size distributions defined by both the local pore geometry and by the largest accessible pore show distinctly different behavior. The accessible pore size is strongly correlated to bone ingrowth. In the specimens studied a strong enhancement of bone ingrowth is observed for pore diameters>100 microm. Little difference in bone ingrowth is measured with different scaffold design. This result illustrates the benefits of microtomography for analyzing the 3D structure of scaffolds and the resultant bone ingrowth.

Recent Fourier and Laplace perspectives for multidimensional NMR in porous media.

Multidimensional NMR techniques used in the measurement of molecular displacements, whether by diffusion or advection, and in the measurement of nuclear spin relaxation times are categorised. Fourier-Fourier, Fourier-Laplace and Laplace-Laplace methods are identified, and recent developments discussed in terms of the separation, correlation and exchange perspective of multidimensional NMR spectroscopy.

Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis - part 2: Calculation of the evolution of percolation and transport properties.

Percolation of reactive fluids in carbonate rocks affects the rock microstructure and hence changes the rock macroscopic properties. In Part 1 paper, we examined the voxel-wise evolution of microstructure of the rock in terms of mineral dissolution/detachment, mineral deposition, and unchanged regions. In the present work, we investigate the relationships between changes in two characteristic transport properties, i.e. permeability and electrical conductivity and two critical parameters of the pore phase, i.e. the fraction of the pore space connecting the inlet and outlet faces of the core sample and the critical pore-throat diameter. We calculate the aforementioned properties on the images of the sample, wherein a homogeneous modification of pore structure occurred in order to ensure the representativeness of the calculated transport properties at the core scale. From images, the evolution of pore connectivity and the potential role of micropores on the connectivity are quantified. It is found that the changing permeability and electrical conductivity distributions along the core length are generally in good agreement with the longitudinal evolution of macro-connected macroporosity and the critical pore-throat diameter. We incorporate microporosity into critical length and permeability calculations and show how microporosity locally plays a role in permeability. It is shown that the Katz-Thompson model reasonably predicts the post-alteration permeability in terms of pre-alteration simulated parameters. This suggests that the evolution of permeability and electrical conductivity of the studied complex carbonate core are controlled by the changes in the macro-connected macroporosity as well as the smallest pore-throats between the connected macropores.

High-fidelity replication of thermoplastic microneedles with open microfluidic channels.

Development of microneedles for unskilled and painless collection of blood or drug delivery addresses the quality of healthcare through early intervention at point-of-care. Microneedles with submicron to millimeter features have been fabricated from materials such as metals, silicon, and polymers by subtractive machining or etching. However, to date, large-scale manufacture of hollow microneedles has been limited by the cost and complexity of microfabrication techniques. This paper reports a novel manufacturing method that may overcome the complexity of hollow microneedle fabrication. Prototype microneedles with open microfluidic channels are fabricated by laser stereolithography. Thermoplastic replicas are manufactured from these templates by soft-embossing with high fidelity at submicron resolution. The manufacturing advantages are (a) direct printing from computer-aided design (CAD) drawing without the constraints imposed by subtractive machining or etching processes, (b) high-fidelity replication of prototype geometries with multiple reuses of elastomeric molds, (c) shorter manufacturing time compared to three-dimensional stereolithography, and (d) integration of microneedles with open-channel microfluidics. Future work will address development of open-channel microfluidics for drug delivery, fluid sampling and analysis.

Structure and properties of clinical coralline implants measured via 3D imaging and analysis.

The development and design of advanced porous materials for biomedical applications requires a thorough understanding of how material structure impacts on mechanical and transport properties. This paper illustrates a 3D imaging and analysis study of two clinically proven coral bone graft samples (Porites and Goniopora). Images are obtained from X-ray micro-computed tomography (micro-CT) at a resolution of 16.8 microm. A visual comparison of the two images shows very different structure; Porites has a homogeneous structure and consistent pore size while Goniopora has a bimodal pore size and a strongly disordered structure. A number of 3D structural characteristics are measured directly on the images including pore volume-to-surface-area, pore and solid size distributions, chord length measurements and tortuosity. Computational results made directly on the digitized tomographic images are presented for the permeability, diffusivity and elastic modulus of the coral samples. The results allow one to quantify differences between the two samples. 3D digital analysis can provide a more thorough assessment of biomaterial structure including the pore wall thickness, local flow, mechanical properties and diffusion pathways. We discuss the implications of these results to the development of optimal scaffold design for tissue ingrowth.

Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis- part 1: Assessment of pore-scale mineral dissolution and deposition.

The application of X-ray micro-computed tomography (µ-CT) for quantitatively characterizing reactive-flow induced pore structure evolution including local particle detachment, displacement and deposition in carbonate rocks is investigated. In the studies conducted in this field of research, the experimental procedure has involved alternating steps of imaging and ex-situ core sample alteration. Practically, it is impossible to return the sample, with micron precision, to the same position and orientation. Furthermore, successive images of a sample in pre- and post-alteration states are usually taken at different conditions such as different scales, resolutions and signal-to-noise ratios. These conditions accompanying with subresolution features in the images make voxel-by-voxel comparisons of successive images problematic. In this paper, we first address the respective challenges in voxel-wise interpretation of successive images of carbonate rocks subject to reactive flow. Reactive coreflood in two carbonate cores with different rock types are considered. For the first rock, we used the experimental and imaging results published by Qajar et al. (2013) which showed a quasi-uniform dissolution regime. A similar reactive core flood was conducted in the second rock which resulted in wormhole-like dissolution regime. We particularly examine the major image processing operations such as transformation of images to the same grey-scale, noise filtering and segmentation thresholding and propose quantitative methods to evaluate the effectiveness of these operations in voxel-wise analysis of successive images of a sample. In the second part, we generalize the methodology based on the three-phase segmentation of normalized images, microporosity assignment and 2D histogram of image intensities to estimate grey-scale changes of individual image voxels for a general case where the greyscale images are segmented into arbitrary number of phases. The results show that local (voxel-based) porosity changes can be decomposed into local mineral dissolution and deposition. Moreover, it is found that the microporosity evolutions are differently distributed in the samples after the reactive coreflood experiments. In the last part of the paper, for the case of quasi-uniform dissolution, we combine the tomographic images with numerical calculations of permeability along the core to characterize the relationship between changes in permeability and the fractions of the mineral dissolved and deposited. A consistency is found between the calculated longitudinal permeability changes and the quantified distribution of mineral dissolved and deposited along the sample.

Cross-property correlations and permeability estimation in sandstone.

We computationally investigate cross-property correlations linking fluid permeability to conductive properties in sedimentary rock for a number of pore size parameters based on three-dimensional digitized rock images. In particular, we focus on correlations based on the pore volume-to-surface-area ratio (V(p)/S), a critical channel diameter (c) associated with mercury porosimetry measurements, length scales associated with the nuclear magnetic resonance relaxation time T2, as well as the mean survival time (tau). Differences between the length scales are discussed. All these correlations yield good agreement with our simulations, but permeability estimates based on the critical diameter (c) are found to be most reliable.

The heterogeneity in femoral neck structure and strength.

Most measures of femoral neck strength derived using dual-energy X-ray absorptiometry or computed tomography (CT) assume the femoral neck is a cylinder with a single cortical thickness. We hypothesized that these simplifications introduce errors in estimating strength and that detailed analyses will identify new parameters that more accurately predict femoral neck strength. High-resolution CT data were used to evaluate 457 cross-sectional slices along the femoral neck of 12 postmortem specimens. Cortical morphology was measured in each cross-section. The distribution of cortical thicknesses was evaluated to determine whether the mean or median better estimated central tendency. Finite-element models were used to calculate the stresses in each cross-section resulting from the peak hip joint forces created during a sideways fall. The relationship between cortical morphology and peak bone stress along the femoral neck was analyzed using multivariate regression analysis. In all cross-sections, cortical thicknesses were non-normally distributed and skewed toward smaller thicknesses (p < 0.0001). The central tendency of cortical thickness was best estimated by the median, not the mean. Stress increased as the median cortical thickness decreased along the femoral neck. The median, not mean, cortical thickness combined with anterior-posterior diameter best predicted peak bone stress generated during a sideways fall (R(2) = 0.66, p < 0.001). Heterogeneity in the structure of the femoral neck determines the diversity of its strength. The median cortical thickness best predicted peak femoral neck stress and is likely to be a relevant predictor of femoral neck fragility.

Morphology and linear-elastic moduli of random network solids.

The effective linear-elastic moduli of disordered network solids are analyzed by voxel-based finite element calculations. We analyze network solids given by Poisson-Voronoi processes and by the structure of collagen fiber networks imaged by confocal microscopy. The solid volume fraction ϕ is varied by adjusting the fiber radius, while keeping the structural mesh or pore size of the underlying network fixed. For intermediate ϕ, the bulk and shear modulus are approximated by empirical power-laws K(phi)proptophin and G(phi)proptophim with n≈1.4 and m≈1.7. The exponents for the collagen and the Poisson-Voronoi network solids are similar, and are close to the values n=1.22 and m=2.11 found in a previous voxel-based finite element study of Poisson-Voronoi systems with different boundary conditions. However, the exponents of these empirical power-laws are at odds with the analytic values of n=1 and m=2, valid for low-density cellular structures in the limit of thin beams. We propose a functional form for K(ϕ) that models the cross-over from a power-law at low densities to a porous solid at high densities; a fit of the data to this functional form yields the asymptotic exponent n≈1.00, as expected. Further, both the intensity of the Poisson-Voronoi process and the collagen concentration in the samples, both of which alter the typical pore or mesh size, affect the effective moduli only by the resulting change of the solid volume fraction. These findings suggest that a network solid with the structure of the collagen networks can be modeled in quantitative agreement by a Poisson-Voronoi process.

Minimal surface scaffold designs for tissue engineering.

Triply-periodic minimal surfaces are shown to be a more versatile source of biomorphic scaffold designs than currently reported in the tissue engineering literature. A scaffold architecture with sheetlike morphology based on minimal surfaces is discussed, with significant structural and mechanical advantages over conventional designs. These sheet solids are porous solids obtained by inflation of cubic minimal surfaces to sheets of finite thickness, as opposed to the conventional network solids where the minimal surface forms the solid/void interface. Using a finite-element approach, the mechanical stiffness of sheet solids is shown to exceed that of conventional network solids for a wide range of volume fractions and material parameters. We further discuss structure-property relationships for mechanical properties useful for custom-designed fabrication by rapid prototyping. Transport properties of the scaffolds are analyzed using Lattice-Boltzmann computations of the fluid permeability. The large number of different minimal surfaces, each of which can be realized as sheet or network solids and at different volume fractions, provides design flexibility essential for the optimization of competing design targets.

The correlation of pore morphology, interconnectivity and physical properties of 3D ceramic scaffolds with bone ingrowth.

In the design of tissue engineering scaffolds, design parameters including pore size, shape and interconnectivity, mechanical properties and transport properties should be optimized to maximize successful inducement of bone ingrowth. In this paper we describe a 3D micro-CT and pore partitioning study to derive pore scale parameters including pore radius distribution, accessible radius, throat radius, and connectivity over the pore space of the tissue engineered constructs. These pore scale descriptors are correlated to bone ingrowth into the scaffolds. Quantitative and visual comparisons show a strong correlation between the local accessible pore radius and bone ingrowth; for well connected samples a cutoff accessible pore radius of approximately 100 microM is observed for ingrowth. The elastic properties of different types of scaffolds are simulated and can be described by standard cellular solids theory: (E/E(0))=(rho/rho(s))(n). Hydraulic conductance and diffusive properties are calculated; results are consistent with the concept of a threshold conductance for bone ingrowth. Simple simulations of local flow velocity and local shear stress show no correlation to in vivo bone ingrowth patterns. These results demonstrate a potential for 3D imaging and analysis to define relevant pore scale morphological and physical properties within scaffolds and to provide evidence for correlations between pore scale descriptors, physical properties and bone ingrowth.