三峡库区巴东组泥灰岩溶蚀作用与岩体力学性质弱化研究

Abstract The karsrt erosion engineering geology became a highlight problem in recent years, in particularly, the karst erosion of marlite of Badong formation made the rock mechanics weaken in Three Gorges Reservoir area, which reduces the safety of slope. During the immigrant construction, many high slopes have been formed, whose instabilities problems pose serious threats to the safety of the people and properties. The accidents of the slope failure take place now and then. By testing, it has been found that the karst erosion pattern and dissolution rate of marlite are not weaker than that of the pure limestone. Furthermore, owning to the weathering and unloading, the karst erosion of the marlite will reach certain depth of the slope, which is named infiltrated karst erosion. The karst erosion made the rock mass quality of slope or foundation worse in a large scale. The karst erosion geological disasters, taken place or not, has become the main restrictive factors to the social stability and economic development. Thus the karst erosion process and mechanism of marlite of Badong formation are studied as the main content of this dissertation. The weakening characteristic of rock mass mechanics parameters are studied along with the rock mass structure deformation and failure processes in the course of the karst erosion. At first, the conditions and influencing factors of the karst erosion are analyzed in the investigative region, on the basis of different karst erosion phenomenon of the marlite and different failure modes of slope. Then via indoor the karst erosion tests, it is analyzed that the karst erosion will change the rock mass composition and its structure. Through test, the different karst erosion phenomena between micro and macro have been observed, and the karst erosion mechanism of the marlite has been summarized. Damage theory is introduced to explain the feature of dissolution pore and the law of crack propagation in the marlite. By microscope and the references data, it can be concluded that the karst erosion process can be divided into rock minerals damage and rock structural damage. And the percent of karst erosion volume is named damage factor, which can be used to describe the quantify karst erosion degree of marlite. Through test, the rock mechanical properties in the different period of karst erosion are studied. Based on the damage mechanics theory and the test result, the relation between the karst erosion degree of marlite and weakening degree of mechanical properties is summarized. By numerical simulations, the karst erosive rock mass mechanics is verified. The conclusion is drawn as below: to the rock mass of marlite, the karst erosion damage made mechanics parameters variation, the deformation modulus, cohesion, and inter friction angle reduce as the negative exponent with the increasing of the karst erosion volume, however, the Poisson ratio increases as the positive exponent with the karst erosion volume increasing. It should be noticed that the deduced formulations are limited to the test data and certain conditions. It is suitable to the rock mass parametric weakening process after the karst erosion of marlite in Three Gorges Reservoir area. Based on the failure types of marlite slope in the field, the karst erosion and weathering process of rock mass are analyzed. And the evolution law of deformation and failure of the marlite mass is studied. The main failure feature of the marlite slope is the karst erosive structure subsidence mode in Three Gorges Reservoir area. The karst erosive structure subsidence mode is explained as follows: the rock mass undergoes the synthetic influence, such as weathering, unloading, corrosion, and so on, many pores and cavities have been formed in the rock mass interior, the rock mass quality is worsen and the rock mass structure is changed, and then the inherent structure of rock mass is collapsed under its gravity, therefore, the failure mode of compaction and subsidence take place. Finally, two examples are used to verify the rock mass parameters in Three Gorges Reservoir area, and the relationship between the marlite slope stability and the time of karst erosion is proposed.

均质弹性条件下多洞平面问题解析算法及其应用研究

The scholars in the world have been trying to find an effective analytic algorithm of multiple hole problems usually meet in engineering designs. Though some studies on circular or elliptic holes had been achieved under specific conditions, no efforts were made to any multiple hole problems that is most significant for engineering designs. The author has made further studies on any multiple hole problems, using complex variable function method and Schwarz alternating method. After solving a series of technological difficulties, the author obtains an effective analytic algorithm, and acquires stress field and displacement field with high accuracy, which can be conducted for arbitrary many iterations according to practical accuracy requirements. In addition, th solution of stress and displacement fields, even for multiple holes of complex shapes and smaller distances. Further, the author made preliminary studies on viscoelastic displacement solution for any double holes. In terms of the obtained displacement solution of any multiple holes, this paper studies displacement back-analysis for the excavations of two tunnels, and find that the back-analysis method is accurate. Additionally, the author presents the mathematical prove of inversion uniqueness for ground stresses, elastic modulus and Poisson ratio. The author believes that the accurate analytic algorithm provided in this paper will presents an effective way to stress and displacement analysis for any multiple hole problems, optimal arrangement of multiple holes, hole shape optimization of multiple holes, etc..

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.

月表热参数模型与月壤厚度的被动微波遥感理论分析

利用被动微波遥感亮度温度数据反演月壤厚度是“嫦娥”探月工程的科学目标之一，也是人类探测月壤厚度的一种新的尝试。深入研究月表太阳辐射、月球内部热流以及月表温度的分布和变化规律，是解译遥感数据，反演月壤厚度的前提条件，也为进一步开展月球探测、开发利用月球资源乃至建立月球基地相关研究工作提供必要的参考。 本文根据月表有效太阳辐照度与太阳常数、日月距离和太阳辐射入射角的关系，建立了月表有效太阳辐照度的实时模型如下： (1) 其中， (2) (3) 通过对月表有效太阳辐照度实时模型的各个参数分析发现，影响月表有效太阳辐照度变化的主要因素是日地距离和太阳辐射入射角的变化。对模型的误差分析表明，从1950年到2050年的100年内，月表有效太阳辐照度计算结果的误差百分比小于0.28％，能更准确地反映月表有效太阳辐照度的变化情况。从2007年月表有效太阳辐照度的计算结果发现，该年内的月表有效太阳辐照度变化在1321.5～1416.6 W•m-2之间，平均为1368.0 W•m-2，一个月内的变化最小幅度为6.0 W•m-2，最大幅度为23.6 W•m-2。 在月表有效太阳辐照度的实时模型基础上，根据能量守恒和Stefan-Boltzmann定律，本文还得出了月表温度分布模型如下： (4) 其中，初始条件由下式决定， (5) 通过与月表温度实际观测结果的比较发现，当月表反射率、热发射率和热惯量分别取0.127、0.94和125 J•m-2•s-1/2•K-1时，模型的计算结果与实际观测值比较符合，能较好地预测理想条件下的月表温度。 月表热参数研究的一个重要应用就是解译对月被动微波遥感的亮度温度数据。在对月被动微波遥感探测中，辐射计获得的亮度温度反映了月球表层的热辐射特性。月球表层的热辐射与其自身的热状况紧密相关，结合文中建立的月表热参数模型，根据辐射传播理论进一步分析了对月微波遥感探测中，月球表层在不同情况下对亮度温度的贡献，确定了亮度温度随月表温度和月壤厚度的变化关系，对被动微波遥感探测月壤厚度的可能性和可能达到的精度进行了估算。 对月球表层的热辐射传播的分析发现，对月被动微波遥感探测获得的亮度温度受月球表层热辐射的控制，与月壤厚度具有指数相关性，并受到月表温度的影响。当月壤和月岩的复介电常数分别为2 + 0.005 j和9 + 1 j、相对磁导率均为1时，对应3.0GHz、7.8GHz、19.35GHz和37.0GHz四个频率的亮度温度与月壤厚度及月表温度的关系可分别近似表示为， 3.0GHz亮度温度： (6) 7.8GHz亮度温度： (7) 19.35GHz亮度温度： (8) 37.0GHz亮度温度： (9) 当月壤厚度和月表温度分别在0.5m~30m和100K~400K之间变化时，上述四个频率的亮度温度变化范围分别在212.5K～252.8K、207.4 K～266.7K、193.8 K～288.6K和174.0 K～310.9K之间。对于较低频率的被动微波遥感，亮度温度随月壤厚度的增大逐渐增大并趋于稳定；对较高频率的被动微波遥感，亮度温度随月壤厚度的增大会产生起伏波动，不利于用单波段反演月壤厚度。亮度温度梯度在频率较高时梯度较大，在很小的月壤厚度范围内很快就趋于0，不利于厚度较大时的月壤厚度反演，但对于厚度较小时的月壤厚度反演精度较高；同时，除3.0GHz外，7.8GHz、19.35GHz和37.0GHz三个频率的亮度温度梯度随月表温度的升高降幅较大，尤其是19.35GHz，适合在夜间对月壤厚度较小的地区进行更精确的探测。对于3.0GHz，其亮度温度梯度受月表温度变化的影响很小，能反映出较深层月壤厚度的信息，可以对月球进行全球全天时探测。若辐射计的分辨率为0.02K，3.0GHz频率对10m厚月壤的判别精度达到0.07m；对于20m厚月壤的精度为1.4m。当月壤厚度小于0.5m时，随着月壤厚度从0到0.5m增加，月球表层的亮度温度贡献呈先减小后增大的趋势，从而使某一亮度温度值可能对应存在两种不同的月壤厚度。因此，对于月壤厚度小于0.5m的区域，利用单波段被动微波遥感亮度温度反演月壤厚度是比较困难的。 在对月被动微波遥感探测中，可以利用月球夜晚时段的亮度温度数据判别月壤厚度是否小于0.5m。当月表温度为100K时，3.0GHz、7.8GHz、19.35GHz和37.0GHz四个频率的亮度温度判别参考值分别为212.9K、207.4K、193.5K和174.1K；月表温度为240K时，上述四个频率的亮度温度判别参考值分别为220.8K、226.8K、234.1K和237.2K。当亮度温度小于参考值时表示月壤厚度小于0.5m，反之，表示月壤厚度大于0.5m。更进一步地，可以根据月表温度的影响系数对月岩是否裸露于月表进行判断。当3.0GHz、7.8GHz、19.35GHz和37.0GHz四个频率的月表温度影响系数接近0.77、0.82、0.84和0.85时，可以认为月岩直接暴露于月表。

Parallel Simulation of Subsonic Fluid Dynamics on a Cluster of Workstations

An effective approach of simulating fluid dynamics on a cluster of non- dedicated workstations is presented. The approach uses local interaction algorithms, small communication capacity, and automatic migration of parallel processes from busy hosts to free hosts. The approach is well- suited for simulating subsonic flow problems which involve both hydrodynamics and acoustic waves; for example, the flow of air inside wind musical instruments. Typical simulations achieve $80\\%$ parallel efficiency (speedup/processors) using 20 HP-Apollo workstations. Detailed measurements of the parallel efficiency of 2D and 3D simulations are presented, and a theoretical model of efficiency is developed which fits closely the measurements. Two numerical methods of fluid dynamics are tested: explicit finite differences, and the lattice Boltzmann method.

Laser-heated high-temperature NMR - a time-resolution study

Submitted to Appl Magn Reson Sponsorship: EPSRC / EU

Simulation of multigate SOI transistors with silicon, germanium and III-V channels

In this work by employing numerical three-dimensional simulations we study the electrical performance and short channel behavior of several multi-gate transistors based on advanced SOI technology. These include FinFETs, triple-gate and gate-all-around nanowire FETs with different channel material, namely Si, Ge, and III-V compound semiconductors, all most promising candidates for future nanoscale CMOS technologies. Also, a new type of transistor called “junctionless nanowire transistor” is presented and extensive simulations are carried out to study its electrical characteristics and compare with the conventional inversion- and accumulation-mode transistors. We study the influence of device properties such as different channel material and orientation, dimensions, and doping concentration as well as quantum effects on the performance of multi-gate SOI transistors. For the modeled n-channel nanowire devices we found that at very small cross sections the nanowires with silicon channel are more immune to short channel effects. Interestingly, the mobility of the channel material is not as significant in determining the device performance in ultrashort channels as other material properties such as the dielectric constant and the effective mass. Better electrostatic control is achieved in materials with smaller dielectric constant and smaller source-to-drain tunneling currents are observed in channels with higher transport effective mass. This explains our results on Si-based devices. In addition to using the commercial TCAD software (Silvaco and Synopsys TCAD), we have developed a three-dimensional Schrödinger-Poisson solver based on the non-equilibrium Green’s functions formalism and in the framework of effective mass approximation. This allows studying the influence of quantum effects on electrical performance of ultra-scaled devices. We have implemented different mode-space methodologies in our 3D quantum-mechanical simulator and moreover introduced a new method to deal with discontinuities in the device structures which is much faster than the coupled-mode-space approach.

Dilute nitride semiconductors : band structure, scattering and high field transport

The substitution of a small fraction x of nitrogen atoms, for the group V elements in conventional III-V semiconductors such as GaAs and GaSb strongly perturbs the conduction band of the host semiconductor. In this thesis we investigate the effects of nitrogen states on the band dispersion, carrier scattering and mobility of dilute nitride alloys. In the supercell model we solve the single particle Hamiltonian for a very large supercell containing randomly placed nitrogen. This model predicts a gap in the density of states of GaNxAs1−x, where this gap is filled in the Green’s function model. Therefore we develop a self-consistent Green’s function (SCGF) approach, which provides excellent agreement with supercell calculations and reveals a gap in the DOS, in contrast with the results of previous non-self-consistent Green’s function calculations. However, including the distribution of N states destroys this gap, as seen in experiment. We then examine the high field transport of carriers by solving the steadystate Boltzmann transport equation and find that it is necessary to include the full distribution of N levels in order to account for the small, low-field mobility and the absence of a negative differential velocity regime observed experimentally with increasing x. Overall the results account well for a wide range of experimental data. We also investigate the band structure, scattering and mobility of carriers by finding the poles of the SCGF, which gives lower carrier mobility for GaNxAs1−x, compared to those already calculated, in better agreement with experiments. The calculated optical absorption spectra for InyGa1−yNxAs1−x and GaNxSb1−x using the SCGF agree well with the experimental data, confirming the validity of this approach to study the band structure of these materials.

Development of inversion-mode and junctionless Indium-Gallium-Arsenide MOSFETs

This PhD covers the development of planar inversion-mode and junctionless Al2O3/In0.53Ga0.47As metal-oxidesemiconductor field-effect transistors (MOSFETs). An implant activation anneal was developed for the formation of the source and drain (S/D) of the inversionmode MOSFET. Fabricated inversion-mode devices were used as test vehicles to investigate the impact of forming gas annealing (FGA) on device performance. Following FGA, the devices exhibited a subthreshold swing (SS) of 150mV/dec., an ION/IOFF of 104 and the transconductance, drive current and peak effective mobility increased by 29%, 25% and 15%, respectively. An alternative technique, based on the fitting of the measured full-gate capacitance vs gate voltage using a selfconsistent Poisson-Schrödinger solver, was developed to extract the trap energy profile across the full In0.53Ga0.47As bandgap and beyond. A multi-frequency inversion-charge pumping approach was proposed to (1) study the traps located at energy levels aligned with the In0.53Ga0.47As conduction band and (2) separate the trapped charge and mobile charge contributions. The analysis revealed an effective mobility (μeff) peaking at ~2850cm2/V.s for an inversion-charge density (Ninv) = 7*1011cm2 and rapidly decreasing to ~600cm2/V.s for Ninv = 1*1013 cm2, consistent with a μeff limited by surface roughness scattering. Atomic force microscopy measurements confirmed a large surface roughness of 1.95±0.28nm on the In0.53Ga0.47As channel caused by the S/D activation anneal. In order to circumvent the issue relative to S/D formation, a junctionless In0.53Ga0.47As device was developed. A digital etch was used to thin the In0.53Ga0.47As channel and investigate the impact of channel thickness (tInGaAs) on device performance. Scaling of the SS with tInGaAs was observed for tInGaAs going from 24 to 16nm, yielding a SS of 115mV/dec. for tInGaAs = 16nm. Flat-band μeff values of 2130 and 1975cm2/V.s were extracted on devices with tInGaAs of 24 and 20nm, respectively

The relationship between rural livelihoods and livestock keeping: the case of Sidama Zone, Southern Ethiopia

This study examines the relationship between rural livelihoods and livestock keeping in Sidama Zones, southern Ethiopia. The livelihood context, assets and strategies of households are the key features of rural livelihoods considered in the study; while households’ livestock ownership, dependence on livestock and livestock management are the main aspects of livestock keeping examined. The study used the sustainable livelihood approach as a framework for data collection and analysis. Describing the main features of rural livelihoods and livestock keeping, and the general pattern of relationship between them, this study mainly aims at identifying the main livelihood factors that determine livestock keeping in the study area. Descriptive statistics, pair wise correlations, mean comparisons and analysis of variance were used to describe rural livelihoods and livestock keeping as well as the relationship between them. Tobit regressions were used to examine the effect of the various livelihood factors on households’ livestock ownership and dependence; Poisson regressions are used to investigate the factors that influence the intensity of livestock management measured by the use of different technologies and inputs. The findings indicated that a number of livelihood factors - assets, livelihood strategies, livelihood shocks and institutional supports - significantly determine the different aspects of livestock keeping. These include: human assets such as age, education and family size; social assets such as membership to social groups; financial assets such as credit; natural assets such as land, and household physical assets; and livelihood strategies such as diversification into farm and nonfarm activities, and coping mechanisms. In addition the livelihood vulnerability context such as shocks and institutional support are among the main determinants of livestock keeping. The results, by and large, matched the findings of previous studies, and it is concluded that households livestock keeping depends on their livelihoods. Accordingly, it is recommended that policies aiming at livestock asset building and productivity improvement should take the livelihoods of rural households in to consideration. As such the study contribute to scholarly works in the area of rural livelihoods, in general, and livestock keeping, in particular. It also contributes to a better understanding of the problems of livestock keeping within the context of rural livelihoods in the country and to the formulation of appropriate policy for the development of the sector.

Effect of different jump distributions on the dynamics of jump processes

The paper investigates stochastic processes forced by independent and identically distributed jumps occurring according to a Poisson process. The impact of different distributions of the jump amplitudes are analyzed for processes with linear drift. Exact expressions of the probability density functions are derived when jump amplitudes are distributed as exponential, gamma, and mixture of exponential distributions for both natural and reflecting boundary conditions. The mean level-crossing properties are studied in relation to the different jump amplitudes. As an example of application of the previous theoretical derivations, the role of different rainfall-depth distributions on an existing stochastic soil water balance model is analyzed. It is shown how the shape of distribution of daily rainfall depths plays a more relevant role on the soil moisture probability distribution as the rainfall frequency decreases, as predicted by future climatic scenarios. © 2010 The American Physical Society.

Bacterial vaginosis as a risk factor for high-grade cervical lesions and cancer in HIV-seropositive women.

OBJECTIVE: To assess the effect of bacterial vaginosis (BV) on the risk of high-grade squamous intraepithelial lesions (HSIL) among HIV-seropositive women. METHODS: A hospital-based prospective cohort study of HIV-seropositive women was conducted in Johannesburg, South Africa from January 2005 to September 2009. Multivariate log-binomial and Poisson regressions were used to estimate prevalence and rate ratios, respectively. RESULTS: Among 1954 HIV-seropositive women, the baseline prevalence of HSIL was 17%. BV prevalence was high (54%) and showed no association with prevalence of HSIL (adjusted prevalence ratio, 1.12; 95% confidence intervals (CI), 0.92-1.35) nor with cervical lesion progression at follow-up visit (n=503) (adjusted rate ratio: 1.00; 95% CI, 0.65-1.53). CONCLUSION: Among HIV-seropositive women, BV was not associated with an increased risk of HSIL or cervical lesion progression.

Bacteria localization and chorion thinning among preterm premature rupture of membranes.

OBJECTIVE: Bacterial colonization of the fetal membranes and its role in pathogenesis of membrane rupture is poorly understood. Prior retrospective work revealed chorion layer thinning in preterm premature rupture of membranes (PPROM) subjects. Our objective was to prospectively examine fetal membrane chorion thinning and to correlate to bacterial presence in PPROM, preterm, and term subjects. STUDY DESIGN: Paired membrane samples (membrane rupture and membrane distant) were prospectively collected from: PPROM = 14, preterm labor (PTL = 8), preterm no labor (PTNL = 8), term labor (TL = 10), and term no labor (TNL = 8), subjects. Sections were probed with cytokeratin to identify fetal trophoblast layer of the chorion using immunohistochemistry. Fluorescence in situ hybridization was performed using broad range 16 s ribosomal RNA probe. Images were evaluated, chorion and choriodecidua were measured, and bacterial fluorescence scored. Chorion thinning and bacterial presence were compared among and between groups using Student's t-test, linear mixed effect model, and Poisson regression model (SAS Cary, NC). RESULTS: In all groups, the fetal chorion cellular layer was thinner at rupture compared to distant site (147.2 vs. 253.7 µm, p<0.0001). Further, chorion thinning was greatest among PPROM subjects compared to all other groups combined, regardless of site sampled [PPROM(114.9) vs. PTL(246.0) vs. PTNL(200.8) vs. TL(217.9) vs. TNL(246.5)]. Bacteria counts were highest among PPROM subjects compared to all other groups regardless of site sampled or histologic infection [PPROM(31) vs. PTL(9) vs. PTNL(7) vs. TL(7) vs. TNL(6)]. Among all subjects at both sites, bacterial counts were inversely correlated with chorion thinning, even excluding histologic chorioamnionitis (p<0.0001 and p = 0.05). CONCLUSIONS: Fetal chorion was uniformly thinner at rupture site compared to distant sites. In PPROM fetal chorion, we demonstrated pronounced global thinning. Although cause or consequence is uncertain, bacterial presence is greatest and inversely correlated with chorion thinning among PPROM subjects.

Lognormal and gamma mixed negative binomial regression

In regression analysis of counts, a lack of simple and efficient algorithms for posterior computation has made Bayesian approaches appear unattractive and thus underdeveloped. We propose a lognormal and gamma mixed negative binomial (NB) regression model for counts, and present efficient closed-form Bayesian inference; unlike conventional Poisson models, the proposed approach has two free parameters to include two different kinds of random effects, and allows the incorporation of prior information, such as sparsity in the regression coefficients. By placing a gamma distribution prior on the NB dispersion parameter r, and connecting a log-normal distribution prior with the logit of the NB probability parameter p, efficient Gibbs sampling and variational Bayes inference are both developed. The closed-form updates are obtained by exploiting conditional conjugacy via both a compound Poisson representation and a Polya-Gamma distribution based data augmentation approach. The proposed Bayesian inference can be implemented routinely, while being easily generalizable to more complex settings involving multivariate dependence structures. The algorithms are illustrated using real examples. Copyright 2012 by the author(s)/owner(s).

On-board Robotic Multi-pinhole SPECT System for Region-of-interest (ROI) Imaging

On-board image guidance, such as cone-beam CT (CBCT) and kV/MV 2D imaging, is essential in many radiation therapy procedures, such as intensity modulated radiotherapy (IMRT) and stereotactic body radiation therapy (SBRT). These imaging techniques provide predominantly anatomical information for treatment planning and target localization. Recently, studies have shown that treatment planning based on functional and molecular information about the tumor and surrounding tissue could potentially improve the effectiveness of radiation therapy. However, current on-board imaging systems are limited in their functional and molecular imaging capability. Single Photon Emission Computed Tomography (SPECT) is a candidate to achieve on-board functional and molecular imaging. Traditional SPECT systems typically take 20 minutes or more for a scan, which is too long for on-board imaging. A robotic multi-pinhole SPECT system was proposed in this dissertation to provide shorter imaging time by using a robotic arm to maneuver the multi-pinhole SPECT system around the patient in position for radiation therapy.

A 49-pinhole collimated SPECT detector and its shielding were designed and simulated in this work using the computer-aided design (CAD) software. The trajectories of robotic arm about the patient, treatment table and gantry in the radiation therapy room and several detector assemblies such as parallel holes, single pinhole and 49 pinholes collimated detector were investigated. The rail mounted system was designed to enable a full range of detector positions and orientations to various crucial treatment sites including head and torso, while avoiding collision with linear accelerator (LINAC), patient table and patient.

An alignment method was developed in this work to calibrate the on-board robotic SPECT to the LINAC coordinate frame and to the coordinate frames of other on-board imaging systems such as CBCT. This alignment method utilizes line sources and one pinhole projection of these line sources. The model consists of multiple alignment parameters which maps line sources in 3-dimensional (3D) space to their 2-dimensional (2D) projections on the SPECT detector. Computer-simulation studies and experimental evaluations were performed as a function of number of line sources, Radon transform accuracy, finite line-source width, intrinsic camera resolution, Poisson noise and acquisition geometry. In computer-simulation studies, when there was no error in determining angles (α) and offsets (ρ) of the measured projections, the six alignment parameters (3 translational and 3 rotational) were estimated perfectly using three line sources. When angles (α) and offsets (ρ) were provided by Radon transform, the estimation accuracy was reduced. The estimation error was associated with rounding errors of Radon transform, finite line-source width, Poisson noise, number of line sources, intrinsic camera resolution and detector acquisition geometry. The estimation accuracy was significantly improved by using 4 line sources rather than 3 and also by using thinner line-source projections (obtained by better intrinsic detector resolution). With 5 line sources, median errors were 0.2 mm for the detector translations, 0.7 mm for the detector radius of rotation, and less than 0.5° for detector rotation, tilt and twist. In experimental evaluations, average errors relative to a different, independent registration technique were about 1.8 mm for detector translations, 1.1 mm for the detector radius of rotation (ROR), 0.5° and 0.4° for detector rotation and tilt, respectively, and 1.2° for detector twist.

Simulation studies were performed to investigate the improvement of imaging sensitivity and accuracy of hot sphere localization for breast imaging of patients in prone position. A 3D XCAT phantom was simulated in the prone position with nine hot spheres of 10 mm diameter added in the left breast. A no-treatment-table case and two commercial prone breast boards, 7 and 24 cm thick, were simulated. Different pinhole focal lengths were assessed for root-mean-square-error (RMSE). The pinhole focal lengths resulting in the lowest RMSE values were 12 cm, 18 cm and 21 cm for no table, thin board, and thick board, respectively. In both no table and thin board cases, all 9 hot spheres were easily visualized above background with 4-minute scans utilizing the 49-pinhole SPECT system while seven of nine hot spheres were visible with the thick board. In comparison with parallel-hole system, our 49-pinhole system shows reduction in noise and bias under these simulation cases. These results correspond to smaller radii of rotation for no-table case and thinner prone board. Similarly, localization accuracy with the 49-pinhole system was significantly better than with the parallel-hole system for both the thin and thick prone boards. Median localization errors for the 49-pinhole system with the thin board were less than 3 mm for 5 of 9 hot spheres, and less than 6 mm for the other 4 hot spheres. Median localization errors of 49-pinhole system with the thick board were less than 4 mm for 5 of 9 hot spheres, and less than 8 mm for the other 4 hot spheres.

Besides prone breast imaging, respiratory-gated region-of-interest (ROI) imaging of lung tumor was also investigated. A simulation study was conducted on the potential of multi-pinhole, region-of-interest (ROI) SPECT to alleviate noise effects associated with respiratory-gated SPECT imaging of the thorax. Two 4D XCAT digital phantoms were constructed, with either a 10 mm or 20 mm diameter tumor added in the right lung. The maximum diaphragm motion was 2 cm (for 10 mm tumor) or 4 cm (for 20 mm tumor) in superior-inferior direction and 1.2 cm in anterior-posterior direction. Projections were simulated with a 4-minute acquisition time (40 seconds per each of 6 gates) using either the ROI SPECT system (49-pinhole) or reference single and dual conventional broad cross-section, parallel-hole collimated SPECT. The SPECT images were reconstructed using OSEM with up to 6 iterations. Images were evaluated as a function of gate by profiles, noise versus bias curves, and a numerical observer performing a forced-choice localization task. Even for the 20 mm tumor, the 49-pinhole imaging ROI was found sufficient to encompass fully usual clinical ranges of diaphragm motion. Averaged over the 6 gates, noise at iteration 6 of 49-pinhole ROI imaging (10.9 µCi/ml) was approximately comparable to noise at iteration 2 of the two dual and single parallel-hole, broad cross-section systems (12.4 µCi/ml and 13.8 µCi/ml, respectively). Corresponding biases were much lower for the 49-pinhole ROI system (3.8 µCi/ml), versus 6.2 µCi/ml and 6.5 µCi/ml for the dual and single parallel-hole systems, respectively. Median localization errors averaged over 6 gates, for the 10 mm and 20 mm tumors respectively, were 1.6 mm and 0.5 mm using the ROI imaging system and 6.6 mm and 2.3 mm using the dual parallel-hole, broad cross-section system. The results demonstrate substantially improved imaging via ROI methods. One important application may be gated imaging of patients in position for radiation therapy.

A robotic SPECT imaging system was constructed utilizing a gamma camera detector (Digirad 2020tc) and a robot (KUKA KR150-L110 robot). An imaging study was performed with a phantom (PET CT PhantomTM), which includes 5 spheres of 10, 13, 17, 22 and 28 mm in diameter. The phantom was placed on a flat-top couch. SPECT projections were acquired with a parallel-hole collimator and a single-pinhole collimator both without background in the phantom, and with background at 1/10th the sphere activity concentration. The imaging trajectories of parallel-hole and pinhole collimated detectors spanned 180 degrees and 228 degrees respectively. The pinhole detector viewed a 14.7 cm-diameter common volume which encompassed the 28 mm and 22 mm spheres. The common volume for parallel-hole was a 20.8-cm-diameter cylinder which encompassed all five spheres in the phantom. The maneuverability of the robotic system was tested by navigating the detector to trace the flat-top table while avoiding collision with the table and maintaining the closest possible proximity to the common volume. For image reconstruction, detector trajectories were described by radius-of-rotation and detector rotation angle θ. These reconstruction parameters were obtained from the robot base and tool coordinates. The robotic SPECT system was able to maneuver the parallel-hole and pinhole collimated SPECT detectors in close proximity to the phantom, minimizing impact of the flat-top couch on detector to center-of-rotation (COR) distance. In no background case, all five spheres were visible in the reconstructed parallel-hole and pinhole images. In with background case, three spheres of 17, 22 and 28 mm diameter were readily observed with the parallel-hole imaging, and the targeted spheres (22 and 28 mm diameter) were readily observed in the pinhole ROI imaging.

In conclusion, the proposed on-board robotic SPECT can be aligned to LINAC/CBCT with a single pinhole projection of the line-source phantom. Alignment parameters can be estimated using one pinhole projection of line sources. This alignment method may be important for multi-pinhole SPECT, where relative pinhole alignment may vary during rotation. For single pinhole and multi-pinhole SPECT imaging onboard radiation therapy machines, the method could provide alignment of SPECT coordinates with those of CBCT and the LINAC. In simulation studies of prone breast imaging and respiratory-gated lung imaging, the 49-pinhole detector showed better tumor contrast recovery and localization in a 4-minute scan compared to parallel-hole detector. On-board SPECT could be achieved by a robot maneuvering a SPECT detector about patients in position for radiation therapy on a flat-top couch. The robot inherent coordinate frames could be an effective means to estimate detector pose for use in SPECT image reconstruction.