Gray-scale Lattice Boltzmann – an attempt to bridge multiple length scales

Understanding and controlling the flow of fluids through porous media such as rocks, fibres, granular media and paper is of fundamental significance to a variety of industries such as oil and gas, chemical production, health and sanitary products. Numerical modelling of this physical process can be difficult not only because of the complex, three- dimensional topology of the porous medium but also because of computational limitations. For example, shale rocks which is now being intensively investigated for its oil and gas resources have porosity over a wide range of length scales from nano-metres up to millimetres. It has been shown that the micro-porosity is fundamental to the fluid movement through the rock. However, current numerical models, which work off computed tomo- graphical (CT) scans of the rock will be excessively large if they are to fully model all length scales which may span six or more orders of magnitude. Here we consider the development of a lattice Boltzmann (LB) technique which may be able to solve the fluid flow over a wide range of length scales. In the past LB techniques have proven to be ideal to model fluid flow in complex porous media since it can readily import and process digital data from CT scans. Hence the fluid flow field is quickly determined and permeabilities can be predicted. However, when the CT data contains micro- porosity, the conventional LB method is not applicable and a modified LB method needs to be developed. Here we consider a gray-scale LB method which works on voxels which are not fully void or solid but something in between, i.e. each voxel is partially resistant to fluid flow. We firstly outline the model, then validate it on test cases and then demonstrate its applicability on real porous media. We develop models not only for single phase fluid flow, but also multiphase fluid flow (i.e. a gas and a liquid) as well as a temperature model, where the temperature field is advected by the fluid flow. For all these cases the models are developed and validated and then demonstrated on realistic media. It is shown that the gray-scale LB model may be able to solve for fluid flow through multiple length scales – a difficult computational problem which is of increasing significance in many real-world applications

Transient ultrasound stimulation has lasting effects on neuronal excitability

Background Transcranial ultrasound stimulation can acutely modulate brain activity, but the lasting effects on neurons are unknown. Objective To assess the excitability profile of neurons in the hours following transient ultrasound stimulation. Methods Primary rat cortical neurons were stimulated with a 40 s, 200 kHz pulsed ultrasound stimulation or sham-stimulation. Intrinsic firing properties were investigated through whole-cell patch-clamp recording by evoking action potentials in response to somatic current injection. Recordings were taken at set timepoints following ultrasound stimulation: 0–2 h, 6–8 h, 12–14 h and 24–26 h. Transmission electron microscopy was used to assess synaptic ultrastructure at the same timepoints. Results In the 0–2 h window, neurons stimulated with ultrasound displayed an increase in the mean frequency of evoked action potentials of 32% above control cell levels (p = 0.023). After 4–6 h this increase was measured as 44% (p = 0.0043). By 12–14 h this effect was eliminated and remained absent 24–26 h post-stimulation. These changes to action potential firing occurred in conjunction with statistically significant differences between control and ultrasound-stimulated neurons in action potential half-width, depolarisation rate, and repolarisation rate, that were similarly eliminated by 24 h following stimulation. These effects occurred in the absence of alterations to intrinsic membrane properties or synaptic ultrastructure. Conclusion We report that stimulating neurons with 40 s of ultrasound enhances their excitability for up to 8 h in conjunction with modifications to action potential kinetics. This occurs in the absence of major ultrastructural change or modification of intrinsic membrane properties. These results can inform the application of transcranial ultrasound in experimental and therapeutic settings

Generation of amorphous carbon and crystallographic texture during low-temperature subseismic slip in calcite fault gouge

Identification of the nano-scale to micro-scale mechanochemical processes occurring during fault slip is of fundamental importance to understand earthquake nucleation and propagation. Here we explore the micromechanical processes occurring during fault nucleation and slip at subseismic rates (∼3 × 10−6 m s–1) in carbonate rocks. We experimentally sheared calcite-rich travertine blocks at simulated upper crustal conditions, producing a nano-grained fault gouge. Strain in the gouge is accommodated by cataclastic comminution of calcite grains and concurrent crystal-plastic deformation through twinning and dislocation glide, producing a crystallographic preferred orientation (CPO). Continued wear of fine-grained gouge particles results in the mechanical decomposition of calcite and production of amorphous carbon. We show that CPO and the production of amorphous carbon, previously attributed to frictional heating and weakening during seismic slip, can be produced at low temperature during stable slip at subseismic rates without slip weakening

Saturation effects on the joint elastic-dielectric properties of carbonates

We used a common microstructural model to investigate the cross-property relations between elastic wave velocities and dielectric permittivity in carbonate rocks. A unified model based on validated self-consistent effective medium theory was used to quantify the effects of porosity and water saturation on both elastic properties (compressional and shear wave velocities) and electromagnetic properties (dielectric permittivity). The results of the forward models are presented as a series of cross-plots covering a wide range of porosities and water saturations and for microstructures that correspond to different predominant aspect ratios. It was found that dielectric permittivity correlated approximately linearly with elastic wave velocity at each saturation stage, with slopes varying gradually from positive at low saturation conditions to negative at higher saturations. The differing sensitivities of the elastic and dielectric rock properties to changes in porosity, pore morphology and water saturation can be used to reduce uncertainty in subsurface fluid saturation estimation when co-located sonic and dielectric surveys are available. The joint approach is useful for cross-validation of rock physics models for analysing pore structure and saturation effects on elastic and dielectric responses

Evaluation of Damage Mechanisms and Skin Factor in Tight Gas Reservoirs

Tight gas reservoirs normally have production problems due to very low matrix permeability and significant damage during well drilling, completion, stimulation and production. Therefore, they may not flow gas at optimum rates without advanced production improvement techniques. The main damage mechanisms and the factors that have significant influence on total skin factor in tight gas reservoirs include mechanical damage to formation rock, plugging of natural fractures by mud solid particles invasion, relative permeability reduction around wellbore as a result of filtrate invasion, liquid leak-off into the formation during fracturing operations, water blocking, skin due to wellbore breakouts, and the damage associated with perforation. Drilling and fracturing fluids invasion mostly occurs through natural fractures and may also lead to serious permeability reduction in the rock matrix that surrounds the natural or hydraulic fractures.This study represents evaluation of different damage mechanisms in tight gas formations, and examines the factors that can have significant influence on total skin factor and well productivity. Reservoir simulation was carried out based on a typical West Australian tight gas reservoir in order to understand how well productivity is affected by each of the damage mechanisms such as natural fractures plugging, mud filtrate invasion, water blocking and perforation