Vertebroplasty: experimental characterization of polymethylmethacrylate bone cement spreading as a function of viscosity, bone porosity, and flow rate

STUDY DESIGN: This is an experimental study on an artificial vertebra model and human cadaveric spine. OBJECTIVE: Characterization of polymethylmethacrylate (PMMA) bone cement distribution in the vertebral body as a function of cement viscosity, bone porosity, and injection speed. Identification of relevant parameters for improved cement flow predictability and leak prevention in vertebroplasty. SUMMARY OF BACKGROUND DATA: Vertebroplasty is an efficient procedure to treat vertebral fractures and stabilize osteoporotic bone in the spine. Severe complications result from bone cement leakage into the spinal canal or the vascular system. Cement viscosity has been identified as an important parameter for leak prevention but the influence of bone structure and injection speed remain obscure. METHODS: An artificial vertebra model based on open porous aluminum foam was used to simulate bone of known porosity. Fifty-six vertebroplasties with 4 different starting viscosity levels and 2 different injection speeds were performed on artificial vertebrae of 3 different porosities. A validation on a human cadaveric spine was executed. The experiments were radiographically monitored and the shape of the cement clouds quantitatively described with the 2 indicators circularity and mean cement spreading distance. RESULTS: An increase in circularity and a decrease in mean cement spreading distance was observed with increasing viscosity, with the most striking change occurring between 50 and 100 Pas. Larger pores resulted in significantly reduced circularity and increased mean cement spreading distance whereas the effect of injection speed on the 2 indicators was not significant. CONCLUSION: Viscosity is the key factor for reducing the risk of PMMA cement leakage and it should be adapted to the degree of osteoporosis encountered in each patient. It may be advisable to opt for a higher starting viscosity but to inject the material at a faster rate.

Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.