Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.

Line focusing for soft x-ray laser-plasma lasing

A computational study of line-focus generation was done using a self-written ray-tracing code and compared to experimental data. Two line-focusing geometries were compared, i.e., either exploiting the sagittal astigmatism of a tilted spherical mirror or using the spherical aberration of an off-axis- illuminated spherical mirror. Line focusing by means of astigmatism or spherical aberration showed identical results as expected for the equivalence of the two frames of reference. The variation of the incidence angle on the target affects the line-focus length, which affects the amplification length such that as long as the irradiance is above the amplification threshold, it is advantageous to have a longer line focus. The amplification threshold is physically dependent on operating parameters and plasma-column conditions and in the present study addresses four possible cases.

Assessment of the mineral density and mineral content of the equine third metacarpal and first phalanx bone by dual energy x-ray absorptiometry

In the first part of this methodological study eleven metacarpi of 9 skeletally normal horses were examined from 4 directions by dual energy x-ray absorptiometry (DXA). The differences between the dorsopalmar-palmarodorsal and lateromedial-mediolateral (opposite sites) bone mineral density (BMD) values were found to be nonsignificant. In the second part of the study the precision of the Norland XR-26 densitometer was tested by measuring 34 metacarpal bones and 34 proximal phalanges, each of them three times, from a single direction. The difference between the individual measurements of the first phalanges and of the metacarpal bones originating from the right or the left side of the same horse were not significant, nor did the age or breed have a significant effect on BMD or bone mineral content (BMC). However, both BMD and BMC are greater in the metacarpal bones than in the proximal phalanges and are higher in geldings than in mares or to stallions, while the BMD or BMC values of mares and stallions did not differ from each other significantly. These data point to the necessity of further BMD studies in a higher number of patients.

[Exposure of the Swiss population to radiation from diagnostic x-ray studies]

A new study is presently being conducted on the exposure of the Swiss population to radiation by diagnostic measures. This study is performed by the Department of Medical Radiation Physics of the University of Berne in collaboration with the Federal Health Bureau and the Swiss Institute for Health and Hospital Matters. In earlier studies the genetically significant exposure of the population and subsequently the median exposure of the red bone marrow had been investigated, whereas now the risk exposure to radiation of as far as possible practically all the risk-relevant organs will be studied. Prior to the initiation of the study, all results of earlier investigations during 1957, 1971 and 1978 were collected and analysed. It was found that the published results are hardly comparable, since the first study was based on individual X-ray examinations and the two subsequent studies on the localised X-ray examinations. To ensure that all data are now comparable, the results of the three studies were appropriately recalculated. Although certain assumptions had to be made that cannot be fully verified any more in view of the time that has elapsed, the collected results will provide a fairly reliable overview of the present-day state of knowledge in this particular field.

Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography

The alveolated structure of the pulmonary acinus plays a vital role in gas exchange function. Three-dimensional (3D) analysis of the parenchymal region is fundamental to understanding this structure-function relationship, but only a limited number of attempts have been conducted in the past because of technical limitations. In this study, we developed a new image processing methodology based on finite element (FE) analysis for accurate 3D structural reconstruction of the gas exchange regions of the lung. Stereologically well characterized rat lung samples (Pediatr Res 53: 72-80, 2003) were imaged using high-resolution synchrotron radiation-based X-ray tomographic microscopy. A stack of 1,024 images (each slice: 1024 x 1024 pixels) with resolution of 1.4 mum(3) per voxel were generated. For the development of FE algorithm, regions of interest (ROI), containing approximately 7.5 million voxels, were further extracted as a working subunit. 3D FEs were created overlaying the voxel map using a grid-based hexahedral algorithm. A proper threshold value for appropriate segmentation was iteratively determined to match the calculated volume density of tissue to the stereologically determined value (Pediatr Res 53: 72-80, 2003). The resulting 3D FEs are ready to be used for 3D structural analysis as well as for subsequent FE computational analyses like fluid dynamics and skeletonization.

Could full-body digital X-ray (LODOX-Statscan) screening in trauma challenge conventional radiography?

ATLS Guidelines recommend single plain radiography of the chest and pelvis as part of the primary survey. Such isolated radiographs, usually obtained by bedside machines, can result in limited, low-quality studies that can adversely affect management. A new digital, low-radiation imaging device, the "Lodox Statscan" (LS), provides full-body anterior and lateral views based on enhanced linear slot-scanning technology in just over 5 minutes. We have the first LS in Europe at our facility. The aim of this study was to compare LS with computed tomographic (CT) scanning, as the gold standard, to determine the sensitivity of LS investigation in detecting injuries to the chest, thoracolumbar spine, and pelvis from our own experience, and to compare our findings with those of conventional radiography in the literature.

Differences between dual X-ray absorptiometry using pencil beam and fan beam modes and their determinants in vivo and in vitro

The aim of this study was to determine the influence of individual factors on differences in bone mineral density (BMD) using dual X-ray absorptiometry pencil beam (PB) and fan beam (FB) modes in vivo and in vitro. PB.BMD and FB.BMD of 63 normal Caucasian females ages 21-80 yr were measured at the lumbar spine and hip. Residuals of the FB/PB regression were used to assess the impact of height, weight, adiposity index (AI) (= weight/height(3/2)), back tissue thickness, and PB.BMD, respectively, on FB/PB difference. The Hologic Anthropomorphic Spine Phantom (ASP) was measured using the PB and FB modes at two different levels to assess the impact of scanning mode and focus distance. The European Spine Phantom (ESP) prototype, a geometrically well-defined phantom with known vertebral densities, was measured using PB and FB modes and analyzed manually to determine the impact of bone density on FB/PB difference and automatically to determine the impact of edge detection on FB/PB difference. Population BMD results were perfectly correlated, but significantly overestimated by 1.5% at the lumbar spine and underestimated by 0.7% at the neck, 1.8% at the trochanter, and 2.0% at the total hip, respectively, when using the FB compared with PB mode. At the lumbar spine, the FB/PB residual correlated negatively with height (r = 0.34, p < 0.01) and PB.BMD (r = 0.48, p <: 0. 0001) and positively with AI (r = 0.26, p < 0.05). At the hip, residual of trochanter correlated positively with weight (r = 0.36, p < 0.01) and AI (r = 0.36, p < 0.01). The FB mode significantly increased ASP BMD by 0.7% compared with PB. Using the FB mode, increasing focus distance significantly (p < 0.001) decreased area and bone mineral content, but not BMD. By contrast, increasing focus distance significantly decreased PB.BMD by 0.7%. With the ESP, the PB mode supplied accurate projected are of the bone (AREA) results but significant underestimation of specified BMD in the manual analysis. The FB mode significantly underestimated PB. AREA by 2.9% but fitted specified BMD quite well. FB/PB overestimation was larger for the low-density (+8.7%) than for the high-density vertebra (+4. 9%). The automated analysis resulted in more than 14% underestimation of PB. AREA (low-density vertebra) and an almost 13% overestimation of PB.BMD (high-density vertebra) using FB. In conclusion, FB and PB measurements are highly correlated at the lumbar spine and hip with small but significant BMD differences related to height, adiposity, and BMD. In clinical practice, it can be erroneous to switch from one method to another, especially in women with low bone density.

Mössbauer, X-ray and Magnetic Studies of Black Sand from the Italian Mediterranean Sea

The study of natural magnetic sands is instrumental to investigate the geological aspects of their formation and of the origin of their territory. In particular, Mössbauer spectroscopy provides unique information on their iron content and on the oxidation state of iron in their mineral composition. The Italian coast on the Mediterranean Sea near Rome is known for the presence of highly magnetic black sands of volcanic origin. A study of the room temperature Mössbauer spec- trum, powder X-ray diffraction, energy dispersive X-ray spectroscopy, and magnetic measurements of a sample of black magnetic sand collected on the seashore of the town of Ladispoli is performed. This study reveals magnetite as main constituent with iron in both tetrahedral and octahedral sites. Minor constituents are the iron minerals hematite and ilmenite, the iron containing minerals diopsite, gossular, and allanite, as well as ubiquitous sanidine, quartz, and calcite.

3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

Synchrotron x ray induced axonal transections in the brain of rats assessed by high-field diffusion tensor imaging tractography

Since approximately two thirds of epileptic patients are non-eligible for surgery, local axonal fiber transections might be of particular interest for them. Micrometer to millimeter wide synchrotron-generated X-ray beamlets produced by spatial fractionation of the main beam could generate such fiber disruptions non-invasively. The aim of this work was to optimize irradiation parameters for the induction of fiber transections in the rat brain white matter by exposure to such beamlets. For this purpose, we irradiated cortex and external capsule of normal rats in the antero-posterior direction with a 4 mm×4 mm array of 25 to 1000 µm wide beamlets and entrance doses of 150 Gy to 500 Gy. Axonal fiber responses were assessed with diffusion tensor imaging and fiber tractography; myelin fibers were examined histopathologically. Our study suggests that high radiation doses (500 Gy) are required to interrupt axons and myelin sheaths. However, a radiation dose of 500 Gy delivered by wide minibeams (1000 µm) induced macroscopic brain damage, depicted by a massive loss of matter in fiber tractography maps. With the same radiation dose, the damage induced by thinner microbeams (50 to 100 µm) was limited to their paths. No macroscopic necrosis was observed in the irradiated target while overt transections of myelin were detected histopathologically. Diffusivity values were found to be significantly reduced. A radiation dose ≤ 500 Gy associated with a beamlet size of < 50 µm did not cause visible transections, neither on diffusion maps nor on sections stained for myelin. We conclude that a peak dose of 500 Gy combined with a microbeam width of 100 µm optimally induced axonal transections in the white matter of the brain.

In-situ X-ray micro-diffraction studies of heterogeneous interfaces between cementitious materials and geological formations

In the present study the challenge of analyzing complex micro X-ray diffraction (microXRD) patterns from cement–clay interfaces has been addressed. In order to extract the maximum information concerning both the spatial distribution and the crystal structure type associated with each of the many diffracting grains in heterogeneous, polycrystalline samples, an approach has been developed in which microXRD was applied to thin sections which were rotated in the X-ray beam. The data analysis, performed on microXRD patterns collected from a filled vein of a cement–clay interface from the natural analogue in Maqarin (Jordan), and a sample from a two-year-old altered interface between cement and argillaceous rock, demonstrate the potential of this method.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

Hohlraum Target for Overcoming Refractive Losses in Plasma X-ray Lasers

Refractive losses in laser-produced plasmas used as gain media are caused by electron density gradients, and limit the energy transport range. The pump pulse is thus deflected from the high-gain region and the short wavelength laser signal also steers away, causing loss of collimation. A Hohlraum used as a target makes the plasma homogeneous and can mitigate refractive losses by means of wave-guiding. A computational study combining a hydrodynamics code and an atomic physics code is presented, which includes a ray-tracing modeling based on the eikonal theory of the trajectory equation. This study presents gain calculations based on population inversion produced by free-electron collisions exciting bound electrons into metastable levels in the 3d94d1(J = 0) → 3d94p1(J = 1) transition of Ni-like Sn. Further, the Hohlraum suggests a dramatic enhancement of the conversion efficiency of collisionally excited x-ray lasing for Ni-like Sn.

Automatic X-ray landmark detection and shape segmentation via data-driven joint estimation of image displacements

In this paper, we propose a new method for fully-automatic landmark detection and shape segmentation in X-ray images. To detect landmarks, we estimate the displacements from some randomly sampled image patches to the (unknown) landmark positions, and then we integrate these predictions via a voting scheme. Our key contribution is a new algorithm for estimating these displacements. Different from other methods where each image patch independently predicts its displacement, we jointly estimate the displacements from all patches together in a data driven way, by considering not only the training data but also geometric constraints on the test image. The displacements estimation is formulated as a convex optimization problem that can be solved efficiently. Finally, we use the sparse shape composition model as the a priori information to regularize the landmark positions and thus generate the segmented shape contour. We validate our method on X-ray image datasets of three different anatomical structures: complete femur, proximal femur and pelvis. Experiments show that our method is accurate and robust in landmark detection, and, combined with the shape model, gives a better or comparable performance in shape segmentation compared to state-of-the art methods. Finally, a preliminary study using CT data shows the extensibility of our method to 3D data.

Evaluation and Comparison of Anatomical Landmark Detection Methods for Cephalometric X-Ray Images: A Grand Challenge

Cephalometric analysis is an essential clinical and research tool in orthodontics for the orthodontic analysis and treatment planning. This paper presents the evaluation of the methods submitted to the Automatic Cephalometric X-Ray Landmark Detection Challenge, held at the IEEE International Symposium on Biomedical Imaging 2014 with an on-site competition. The challenge was set to explore and compare automatic landmark detection methods in application to cephalometric X-ray images. Methods were evaluated on a common database including cephalograms of 300 patients aged six to 60 years, collected from the Dental Department, Tri-Service General Hospital, Taiwan, and manually marked anatomical landmarks as the ground truth data, generated by two experienced medical doctors. Quantitative evaluation was performed to compare the results of a representative selection of current methods submitted to the challenge. Experimental results show that three methods are able to achieve detection rates greater than 80% using the 4 mm precision range, but only one method achieves a detection rate greater than 70% using the 2 mm precision range, which is the acceptable precision range in clinical practice. The study provides insights into the performance of different landmark detection approaches under real-world conditions and highlights achievements and limitations of current image analysis techniques.