环渤海花岗岩体地下油库建设的工程地质力学研究

As a kind of strategic resource，petroleum play an very important role in current social stability, economic development and state safety. Since 1993 China has turned from a net oil exporter into a net oil importer, the figure of imported oil increased from then on. In 2004 China's total energy consumption exceeded Japan’s, and ranked in the second place, just inferior to America. Today China is the world’s third-largest importing nation, accounting for 6% of world imports and 8% of world consumption. Comparing with other strategic petroleum reserve schemes, underground oil storage possess many advantages, such as security, economy, less pollution, save land, suited for strategic reserve and so on, so it is the most ideal form for strategic petroleum reserve. In the background of China Strategic Petroleum Reserve Program started just now, this paper choose Circum-Bo sea region as a study area, and do some system study on the underground oil storage caverns constructed in inter-large granite rock masses in Circum-Bo sea region. On the foundation of a great amount of information come from both home and abroad, firstly this paper analysed the principle, economy, cavern shape, profile dimension, and gain some cognizances and logos, as follows: ①Hard rock mass such as granite is the major rock, in which underground oil storage are constructed; ②Unlined underground oil storage caverns had been wide spread used as a sort of oil storage form abroad, there already exist a suit of skilled experience and technologies to prevent oil product from leaking; ③Compared with surface tanks, underground oil storage cavern possess predominance in economy clearly. In general, it will be more economical when the storage capacity exceed 50000m3. The quality of rock mass is the most important factor for underground storage cost, however such as hydrogeology, storage capacity, the number of storage galleries, the length, storage product, mechanical equipments, geographic location also influent the cost. In designed depth of the underground storage, the rock mass of Jinzhou mainly belong to class Ⅱ, but parts with dykes, clayization alteration, and dense joints are Ⅲ, Ⅳ; ④Now, there are few underground oil storages span more than 25m in both abroad and home. The examples of some ancient underground works and modern underground excavation with wide span surely give us many precious elicitations to construct more great unlined storage caverns, when the rock mass quality is good, cavern shape and construction method also are proper, it is quite possible to construct underground oil storage cavern with span more than 30m . The main axis orientation of Jinzhou underground oil storage cavern is NW direction, the cavern's elevation locate between -53msl and -76msl. The storage's total volume is about 3×106m3, composed of 8 parallel galleries with 950m length, the pillars between them are 45m, and every two of galleries form one unit, which can deposit 75×104m3 for each unit. The product will be stored are Saudi light and Saudi medium crude oil, the main cavern's section is 411.5m2, with 23m height and 19m width. According to the principle and technique of engineering geomechanics, this study supply a sort of system scientific thinking and method for sitting location of underground oil storage in granite region: ① On the foundation of the earth crust stability sub-zone appraise of Circum-Bo sea region, farther research concerning granite distribution, genesis, geological period and fault structure are conducted in stable areas, generally, this paper select Liaoxi, east shore of Liaotung peninsula and Jiaotung peninsula as target areas for underground oil storage regions, where Mesozoic granite is magnitude; ②After roundly comparison in facts of geologic structure, engineering geology, hydrogeology, topography, transportation and so on of three granite distributed areas, at last, selecting Jingzhou granite zone in Liaoxi out as an ideal construction area; ③ Detailed investigation is conducted in the southeast of Baimashi in Jingzhou development district, the final field. Ultrasonic Borehole Television, as a major way to collect original information of borehole rock mass were used, which is very effective to appraise the quality of deep rock mass; ④ According to the field data of tectonic stress, rock mass quality, the spatial distribution of fracture water, some optimum designs in cross section, axial direction and cavern span have been designed for the underground oil storage cavern layout in Jinzhou. To understand the characteristics of swelling alteration rock in Jinzhou granite mass, collected abundant swelling alteration rock engineering examples in granite, which study them in detail, concluded the swelling alteration rock distribute nearly everywhere in China, intruded medium-basic dykes alteration, along discontinuities and mineral hydrothermal alteration with genesis of granite are three main forms clayization alteration rock in granite rock mass. In Jinzhou field, from macro to micro studied the swelling rock which induced by mid-basic dyke intrusion, with weak swelling. In conclusion, this paper conclude the distribution rule and features of expansion alteration rock in filed, and advise some technical suggestions for excavation at swelling alteration rock part. The main features of this paper: ①In the process of site selection, investigation and design, a suit of technique and method of engineering geomechanics metasynthesis were formed, which is significative to guide the large scale underground oil storage cavern sitting location, investigation and design in granite rock mass; ②The detailed discussion on the engineering geology problems in granite mass, such as weathering crust, faults, dykes and clayization alteration rock, are useful for other projects in aspects of site selection, engineering geology evaluation and stability estimation; ③The summary and integration of the genesis, type, countermeasure relate to swelling alteration rock, also is likely to be used for other underground oil storage caverns constructed in swelling alteration granite. In conclusion, this study is meaningful for guiding the large scale underground oil storage for site selection, investigation and design in granite rock mass.

基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

大别山硬玉石英岩带中的水

The Dabie Mountains is a collisional orogenic belt between the North China and Yantze Continental plates. It is the eastern elongation of the Tongbai and Qingling orogen, and is truncated at its east end by the Tan-Lu fault. Jadeite-quartzite belt occurs in the eastern margin of UHPMB from the Dabie Mountains. Geochemical features indicate that the protoliths of the jadeite-quartzite and associated eclogite to be supracrustal rocks. The occurrence of micro-inclusions of coesite in jadeite and garnet confirmed that the continental crust can be subducted to great depth (8 0-100km) and then exhumed rapidly with its UHP mineral signature fairly preserved. Therefore, study of UHP jadeite-quartzite provides important information on subduction of continental crustal rocks and their exhumation histories, as well as the dynamics of plate tectonic processes at convergent margins. The purpose of this paper is to investigate the presence of hydrous component in the jadeite-quartzite belt, significant natural variations in the hydrous component content of UHP minerals and to discuss the role of water in petrology, geochemistry and micro-tectonic. On the basis of our previous studies, some new geological evidences have been found in the jadeite-quartzite belt by researches on petrography, mineralogy, micro-tectonic, hydrous component content of UHP minerals and combined with the study on rheology of materials using microprob, ER, TEM. By research and analysis of these phenomenona, the results obtained are as follows: 1. The existence of fluid during ultra-high pressure metamorphic process. Jadeites, omphacite, garnet, rutile, coesite and quartz from the jadeite-quartzite belt have been investigated by Fourier transform infrared spectrometer and TEM. Results show that all of these minerals contain trace amount of water which occur as hydroxyl and free-water in these minerals. The two-type hydrous components in UHP minerals are indicated stable in the mantle-depth. The results demonstrated that these ultra-high pressure metamorphic minerals, which were derived from continental crust protoliths, they could bring water into the mantle depth during the ultra-high pressure metamorphism. The clusters of water molecules within garnet are very important evidence of the existence of fluid during ultra-high pressure metamorphic process. It indicated that the metamorphic system was not "dry"during the ultra-high pressure stage. 2.The distribution of hydrous component in UHP minerals of jadeite-quartzite. The systematic distribution of hydrous components in UHP minerals are a strong indication that water in these minerals, are controlled by some factors and that the observed variations are not of a random nature. The distribution and concentration of hydrous component is not only correlated with composition of minerals, but also a function of geological environment. Therefore, the hydrous component in the minerals can not only take important part in the UHP metamorphic fluid during subduction of continental crustal rocks, but also their hydroxyl transported water molecules with decreasing pressure during their exhumation. And these water molecules can not only promote the deformation of jadeite through hydrolytic weakening, but also may be the part of the retrograde metamorphic fluid. 3.The role of water in the deformed UHP minerals. The jadeite, omphacite, garnet are strong elongated deformation in the jadeite-quartzite from the Dabie Mountains. They are (1) they are developed strong plastic deformation; (2) developed dislocation loop, dislocation wall; (3) the existence of clusters of water molecular in the garnet; and (4) the evolution of micero-tectonic from clusters of water molecular-dislocation loop in omphacite. That indicated that the water weakening controlled the mechanism of deformed minerals. Because the data presented here are not only the existence of clusters of water molecular in the garnet, but also developed strong elongation, high density of dislocation and high aspect ratios, adding microprobe data demonstrate the studied garnet crystals no compositional zoning. Therefore, this indicates that the diffusion process of the grain boundary mobility did not take place in these garnets. On the basis of above features, we consider that it can only be explained by plastic deformation of the garnets. The clusters of water molecules present in garnet was directly associated with mechanical weakening and inducing in plastic deformation of garnet by glissile dislocations. Investigate of LPO, strain analysis, TEM indicated that these clinopyroxenes developed strong elongation, high aspect ratios, and developed dislocation loop, dislocation wall and free dislocations. These indicated that the deformation mechanism of the clinopyroxenes plastically from the Dabie Mountains is dominant dislocation creep under the condition of the UHP metamorphic conditions. There are some bubbles with dislocation loops attached to them in the omphacite crystal. The bubbles attached to the dislocation loops sometimes form a string of bubble beads and some loops are often connected to one another via a common bubble. The water present in omphacite was directly associated with hydrolitic weakening and inducing in plastic deformation of omphacite by dislocations. The role of water in brittle deformation. Using microscopy, deformation has been identified as plastic deformation and brittle deformation in UHP minerals from the Dabie Mountains. The study of micro-tectonic on these minerals shows that the brittle deformation within UHP minerals was related to local stresses. The brittle deformation is interpreted as being caused by an interaction of high fluid pressure, volume changes. The hydroxyl within UHP minerals transported water molecules with decreasing pressure due to their exhumation. However, under eclogite facies conditions, the litho-static pressure is extreme, but a high fluid pressure will reduce the effective stress and make brittle deformation possible. The role of water in prograde metamorphism. Geochemical research on jadeite-quartzite and associated eclogite show that the protoliths of these rocks are supracrustal rocks. With increasing of temperature and pressure, the chlorite, biotite, muscovite was dehydrous reaction and released hydrous component during the subduction of continental lithosphere. The supracrustal rocks were transformed UHP rocks and formed UHP facies assemblage promotely by water introduction, and was retained in UHP minerals as hydrous component. The water within UHP minerals may be one of the retrograde metamorphic fluids. Petrological research on UHP rocks of jadeite-quartzite belt shows that there was existence of local fluids during early retrograde metamorphism. That are: (1) coronal textures and symplectite around relict UHP minerls crystals formed from UHP minerls by hydration reactions; (2) coronal textures of albite around ruitle; and (3) micro-fractures in jadeite or garnet were filled symplectite of Amp + PI + Mt. That indicated that the reactions of early retrograde metamorphism dependent on fluid introduction. These fluids not only promoted retrograde reaction of UHP minerals, but also were facilitate to diffuse intergranular and promote growth in minerals. Therefore, the hydrous component in the UHP minerals can not only take important part in the UHP metamorphic fluid during subduction of continental crustal rocks, but also their hydroxyl transport water molecules with decreasing pressure and may take part in the retrograde metamorphic fluid during their exhumation. 7. The role of water in geochemistry of UHP jadeite-quartzite. Geochemical research show that there are major, trace and rare earth element geochemical variations in the jadeite-quartzite from the Changpu district of Dabie Mountains, during retrograde metamorphic processes from the jadeite-quartzite--gneiss. The elements such as SiO_2、FeO、Ba、Zr、Ga、La、Ce、PTN Nd% Sm and Eu increase gradually from the jadeite-quartzite to retrograded jadeite-quartzite and to gneiss, whilst TiO_2. Na_2CK Fe2O_3、Rb、Y、Nb、Gd、Tb、Dy、Ho、Er、Tm、Yb decrease gradually. And its fO_2 keep nearly unchanged during early retrograde metamorphism, but decreased obviously during later retrograde metamorphism. These indicate that such changes are not only controlled by element transformation between mineralogical phases, but also closely relative to fluid-rock interaction in the decompression retrograde metamorphic processes.

A Binocular, Foveated Active Vision System

This report documents the design and implementation of a binocular, foveated active vision system as part of the Cog project at the MIT Artificial Intelligence Laboratory. The active vision system features a three degree of freedom mechanical platform that supports four color cameras, a motion control system, and a parallel network of digital signal processors for image processing. To demonstrate the capabilities of the system, we present results from four sample visual-motor tasks.

Quantitative Inference in a Mechanical Design Compiler

This paper presents the ideas underlying a program that takes as input a schematic of a mechanical or hydraulic power transmission system, plus specifications and a utility function, and returns catalog numbers from predefined catalogs for the optimal selection of components implementing the design. It thus provides the designer with a high level "language" in which to compose new designs, then performs some of the detailed design process for him. The program is based on a formalization of quantitative inferences about hierarchically organized sets of artifacts and operating conditions, which allows design compilation without the exhaustive enumeration of alternatives.

Complex System Classification

The use of terms such as “Engineering Systems”, “System of systems” and others have been coming into greater use over the past decade to denote systems of importance but with implied higher complexity than for the term systems alone. This paper searches for a useful taxonomy or classification scheme for complex Systems. There are two aspects to this problem: 1) distinguishing between Engineering Systems (the term we use) and other Systems, and 2) differentiating among Engineering Systems. Engineering Systems are found to be differentiated from other complex systems by being human-designed and having both significant human complexity as well as significant technical complexity. As far as differentiating among various engineering systems, it is suggested that functional type is the most useful attribute for classification differentiation. Information, energy, value and mass acted upon by various processes are the foundation concepts underlying the technical types.

Control of Vibration in Mechanical Systems Using Shaped Reference Inputs

Dynamic systems which undergo rapid motion can excite natural frequencies that lead to residual vibration at the end of motion. This work presents a method to shape force profiles that reduce excitation energy at the natural frequencies in order to reduce residual vibration for fast moves. Such profiles are developed using a ramped sinusoid function and its harmonics, choosing coefficients to reduce spectral energy at the natural frequencies of the system. To improve robustness with respect to parameter uncertainty, spectral energy is reduced for a range of frequencies surrounding the nominal natural frequency. An additional set of versine profiles are also constructed to permit motion at constant speed for velocity-limited systems. These shaped force profiles are incorporated into a simple closed-loop system with position and velocity feedback. The force input is doubly integrated to generate a shaped position reference for the controller to follow. This control scheme is evaluated on the MIT Cartesian Robot. The shaped inputs generate motions with minimum residual vibration when actuator saturation is avoided. Feedback control compensates for the effect of friction Using only a knowledge of the natural frequencies of the system to shape the force inputs, vibration can also be attenuated in modes which vibrate in directions other than the motion direction. When moving several axes, the use of shaped inputs allows minimum residual vibration even when the natural frequencies are dynamically changing by a limited amount.

Quantum behavior in mesoscopic systems

In this thesis I present the work done during my PhD. The Thesis is divided into two parts; in the first one I present the study of mesoscopic quantum systems whereas in the second one I address the problem of the definition of Markov regime for quantum system dynamics. The first work presented is the study of vortex patterns in (quasi) two dimensional rotating Bose Einstein condensates (BECs). I consider the case of an anisotropy trapping potential and I shall show that the ground state of the system hosts vortex patterns that are unstable. In a second work I designed an experimental scheme to transfer entanglement from two entangled photons to two BECs. This work is meant to propose a feasible experimental set up to bring entanglement from microscopic to macroscopic systems for both the study of fundamental questions (quantum to classical transition) and technological applications. In the last work of the first part another experimental scheme is presented in order to detect coherences of a mechanical oscillator which is assumed to have been previously cooled down to the quantum regime. In this regime in fact the system can rapidly undergo decoherence so that new techniques have to be employed in order to detect and manipulate their states. In the scheme I propose a micro-mechanical oscillator is coupled to a BEC and the detection is performed by monitoring the BEC with a negligible back-action on the cantilever. In the second part of the thesis I give a definition of Markov regime for open quantum dynamics. The importance of such definition comes from both the mathematical description of the system dynamics and from the understanding of the role played by the environment in the evolution of an open system. In the Markov regime the mathematical description can be simplified and the role of the environment is a passive one.

Investigation of micro-devices for neurobiological applications

The aim of this project is to integrate neuronal cell culture with commercial or in-house built micro-electrode arrays and MEMS devices. The resulting device is intended to support neuronal cell culture on its surface, expose specific portions of a neuronal population to different environments using microfluidic gradients and stimulate/record neuronal electrical activity using micro-electrode arrays. Additionally, through integration of chemical surface patterning, such device can be used to build neuronal cell networks of specific size, conformation and composition. The design of this device takes inspiration from the nervous system because its development and regeneration are heavily influenced by surface chemistry and fluidic gradients. Hence, this device is intended to be a step forward in neuroscience research because it utilizes similar concepts to those found in nature. The large part of this research revolved around solving technical issues associated with integration of biology, surface chemistry, electrophysiology and microfluidics. Commercially available microelectrode arrays (MEAs) are mechanically and chemically brittle making them unsuitable for certain surface modification and micro-fluidic integration techniques described in the literature. In order to successfully integrate all the aspects into one device, some techniques were heavily modified to ensure that their effects on MEA were minimal. In terms of experimental work, this thesis consists of 3 parts. The first part dealt with characterization and optimization of surface patterning and micro-fluidic perfusion. Through extensive image analysis, the optimal conditions required for micro-contact printing and micro-fluidic perfusion were determined. The second part used a number of optimized techniques and successfully applied these to culturing patterned neural cells on a range of substrates including: Pyrex, cyclo-olefin and SiN coated Pyrex. The second part also described culturing neurons on MEAs and recording electrophysiological activity. The third part of the thesis described integration of MEAs with patterned neuronal culture and microfluidic devices. Although integration of all methodologies proved difficult, a large amount of data relating to biocompatibility, neuronal patterning, electrophysiology and integration was collected. Original solutions were successfully applied to solve a number of issues relating to consistency of micro printing and microfluidic integration leading to successful integration of techniques and device components.

A low-cost, portable, and quantitative spectral imaging system for application to biological tissues.

The ability of diffuse reflectance spectroscopy to extract quantitative biological composition of tissues has been used to discern tissue types in both pre-clinical and clinical cancer studies. Typically, diffuse reflectance spectroscopy systems are designed for single-point measurements. Clinically, an imaging system would provide valuable spatial information on tissue composition. While it is feasible to build a multiplexed fiber-optic probe based spectral imaging system, these systems suffer from drawbacks with respect to cost and size. To address these we developed a compact and low cost system using a broadband light source with an 8-slot filter wheel for illumination and silicon photodiodes for detection. The spectral imaging system was tested on a set of tissue mimicking liquid phantoms which yielded an optical property extraction accuracy of 6.40 +/- 7.78% for the absorption coefficient (micro(a)) and 11.37 +/- 19.62% for the wavelength-averaged reduced scattering coefficient (micro(s)').

Micro-Anatomical Quantitative Imaging Towards Enabling Automated Diagnosis of Thick Tissues at the Point of Care

Histopathology is the clinical standard for tissue diagnosis. However, histopathology has several limitations including that it requires tissue processing, which can take 30 minutes or more, and requires a highly trained pathologist to diagnose the tissue. Additionally, the diagnosis is qualitative, and the lack of quantitation leads to possible observer-specific diagnosis. Taken together, it is difficult to diagnose tissue at the point of care using histopathology.

Several clinical situations could benefit from more rapid and automated histological processing, which could reduce the time and the number of steps required between obtaining a fresh tissue specimen and rendering a diagnosis. For example, there is need for rapid detection of residual cancer on the surface of tumor resection specimens during excisional surgeries, which is known as intraoperative tumor margin assessment. Additionally, rapid assessment of biopsy specimens at the point-of-care could enable clinicians to confirm that a suspicious lesion is successfully sampled, thus preventing an unnecessary repeat biopsy procedure. Rapid and low cost histological processing could also be potentially useful in settings lacking the human resources and equipment necessary to perform standard histologic assessment. Lastly, automated interpretation of tissue samples could potentially reduce inter-observer error, particularly in the diagnosis of borderline lesions.

To address these needs, high quality microscopic images of the tissue must be obtained in rapid timeframes, in order for a pathologic assessment to be useful for guiding the intervention. Optical microscopy is a powerful technique to obtain high-resolution images of tissue morphology in real-time at the point of care, without the need for tissue processing. In particular, a number of groups have combined fluorescence microscopy with vital fluorescent stains to visualize micro-anatomical features of thick (i.e. unsectioned or unprocessed) tissue. However, robust methods for segmentation and quantitative analysis of heterogeneous images are essential to enable automated diagnosis. Thus, the goal of this work was to obtain high resolution imaging of tissue morphology through employing fluorescence microscopy and vital fluorescent stains and to develop a quantitative strategy to segment and quantify tissue features in heterogeneous images, such as nuclei and the surrounding stroma, which will enable automated diagnosis of thick tissues.

To achieve these goals, three specific aims were proposed. The first aim was to develop an image processing method that can differentiate nuclei from background tissue heterogeneity and enable automated diagnosis of thick tissue at the point of care. A computational technique called sparse component analysis (SCA) was adapted to isolate features of interest, such as nuclei, from the background. SCA has been used previously in the image processing community for image compression, enhancement, and restoration, but has never been applied to separate distinct tissue types in a heterogeneous image. In combination with a high resolution fluorescence microendoscope (HRME) and a contrast agent acriflavine, the utility of this technique was demonstrated through imaging preclinical sarcoma tumor margins. Acriflavine localizes to the nuclei of cells where it reversibly associates with RNA and DNA. Additionally, acriflavine shows some affinity for collagen and muscle. SCA was adapted to isolate acriflavine positive features or APFs (which correspond to RNA and DNA) from background tissue heterogeneity. The circle transform (CT) was applied to the SCA output to quantify the size and density of overlapping APFs. The sensitivity of the SCA+CT approach to variations in APF size, density and background heterogeneity was demonstrated through simulations. Specifically, SCA+CT achieved the lowest errors for higher contrast ratios and larger APF sizes. When applied to tissue images of excised sarcoma margins, SCA+CT correctly isolated APFs and showed consistently increased density in tumor and tumor + muscle images compared to images containing muscle. Next, variables were quantified from images of resected primary sarcomas and used to optimize a multivariate model. The sensitivity and specificity for differentiating positive from negative ex vivo resected tumor margins was 82% and 75%. The utility of this approach was further tested by imaging the in vivo tumor cavities from 34 mice after resection of a sarcoma with local recurrence as a bench mark. When applied prospectively to images from the tumor cavity, the sensitivity and specificity for differentiating local recurrence was 78% and 82%. The results indicate that SCA+CT can accurately delineate APFs in heterogeneous tissue, which is essential to enable automated and rapid surveillance of tissue pathology.

Two primary challenges were identified in the work in aim 1. First, while SCA can be used to isolate features, such as APFs, from heterogeneous images, its performance is limited by the contrast between APFs and the background. Second, while it is feasible to create mosaics by scanning a sarcoma tumor bed in a mouse, which is on the order of 3-7 mm in any one dimension, it is not feasible to evaluate an entire human surgical margin. Thus, improvements to the microscopic imaging system were made to (1) improve image contrast through rejecting out-of-focus background fluorescence and to (2) increase the field of view (FOV) while maintaining the sub-cellular resolution needed for delineation of nuclei. To address these challenges, a technique called structured illumination microscopy (SIM) was employed in which the entire FOV is illuminated with a defined spatial pattern rather than scanning a focal spot, such as in confocal microscopy.

Thus, the second aim was to improve image contrast and increase the FOV through employing wide-field, non-contact structured illumination microscopy and optimize the segmentation algorithm for new imaging modality. Both image contrast and FOV were increased through the development of a wide-field fluorescence SIM system. Clear improvement in image contrast was seen in structured illumination images compared to uniform illumination images. Additionally, the FOV is over 13X larger than the fluorescence microendoscope used in aim 1. Initial segmentation results of SIM images revealed that SCA is unable to segment large numbers of APFs in the tumor images. Because the FOV of the SIM system is over 13X larger than the FOV of the fluorescence microendoscope, dense collections of APFs commonly seen in tumor images could no longer be sparsely represented, and the fundamental sparsity assumption associated with SCA was no longer met. Thus, an algorithm called maximally stable extremal regions (MSER) was investigated as an alternative approach for APF segmentation in SIM images. MSER was able to accurately segment large numbers of APFs in SIM images of tumor tissue. In addition to optimizing MSER for SIM image segmentation, an optimal frequency of the illumination pattern used in SIM was carefully selected because the image signal to noise ratio (SNR) is dependent on the grid frequency. A grid frequency of 31.7 mm-1 led to the highest SNR and lowest percent error associated with MSER segmentation.

Once MSER was optimized for SIM image segmentation and the optimal grid frequency was selected, a quantitative model was developed to diagnose mouse sarcoma tumor margins that were imaged ex vivo with SIM. Tumor margins were stained with acridine orange (AO) in aim 2 because AO was found to stain the sarcoma tissue more brightly than acriflavine. Both acriflavine and AO are intravital dyes, which have been shown to stain nuclei, skeletal muscle, and collagenous stroma. A tissue-type classification model was developed to differentiate localized regions (75x75 µm) of tumor from skeletal muscle and adipose tissue based on the MSER segmentation output. Specifically, a logistic regression model was used to classify each localized region. The logistic regression model yielded an output in terms of probability (0-100%) that tumor was located within each 75x75 µm region. The model performance was tested using a receiver operator characteristic (ROC) curve analysis that revealed 77% sensitivity and 81% specificity. For margin classification, the whole margin image was divided into localized regions and this tissue-type classification model was applied. In a subset of 6 margins (3 negative, 3 positive), it was shown that with a tumor probability threshold of 50%, 8% of all regions from negative margins exceeded this threshold, while over 17% of all regions exceeded the threshold in the positive margins. Thus, 8% of regions in negative margins were considered false positives. These false positive regions are likely due to the high density of APFs present in normal tissues, which clearly demonstrates a challenge in implementing this automatic algorithm based on AO staining alone.

Thus, the third aim was to improve the specificity of the diagnostic model through leveraging other sources of contrast. Modifications were made to the SIM system to enable fluorescence imaging at a variety of wavelengths. Specifically, the SIM system was modified to enabling imaging of red fluorescent protein (RFP) expressing sarcomas, which were used to delineate the location of tumor cells within each image. Initial analysis of AO stained panels confirmed that there was room for improvement in tumor detection, particularly in regards to false positive regions that were negative for RFP. One approach for improving the specificity of the diagnostic model was to investigate using a fluorophore that was more specific to staining tumor. Specifically, tetracycline was selected because it appeared to specifically stain freshly excised tumor tissue in a matter of minutes, and was non-toxic and stable in solution. Results indicated that tetracycline staining has promise for increasing the specificity of tumor detection in SIM images of a preclinical sarcoma model and further investigation is warranted.

In conclusion, this work presents the development of a combination of tools that is capable of automated segmentation and quantification of micro-anatomical images of thick tissue. When compared to the fluorescence microendoscope, wide-field multispectral fluorescence SIM imaging provided improved image contrast, a larger FOV with comparable resolution, and the ability to image a variety of fluorophores. MSER was an appropriate and rapid approach to segment dense collections of APFs from wide-field SIM images. Variables that reflect the morphology of the tissue, such as the density, size, and shape of nuclei and nucleoli, can be used to automatically diagnose SIM images. The clinical utility of SIM imaging and MSER segmentation to detect microscopic residual disease has been demonstrated by imaging excised preclinical sarcoma margins. Ultimately, this work demonstrates that fluorescence imaging of tissue micro-anatomy combined with a specialized algorithm for delineation and quantification of features is a means for rapid, non-destructive and automated detection of microscopic disease, which could improve cancer management in a variety of clinical scenarios.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.

Monte Carlo Analysis and Physics Characterization of a Novel Nanoparticle Detector for Medical and Micro-dosimetry Applications

The outcomes for both (i) radiation therapy and (ii) preclinical small animal radio- biology studies are dependent on the delivery of a known quantity of radiation to a specific and intentional location. Adverse effects can result from these procedures if the dose to the target is too high or low, and can also result from an incorrect spatial distribution in which nearby normal healthy tissue can be undesirably damaged by poor radiation delivery techniques. Thus, in mice and humans alike, the spatial dose distributions from radiation sources should be well characterized in terms of the absolute dose quantity, and with pin-point accuracy. When dealing with the steep spatial dose gradients consequential to either (i) high dose rate (HDR) brachytherapy or (ii) within the small organs and tissue inhomogeneities of mice, obtaining accurate and highly precise dose results can be very challenging, considering commercially available radiation detection tools, such as ion chambers, are often too large for in-vivo use.

In this dissertation two tools are developed and applied for both clinical and preclinical radiation measurement. The first tool is a novel radiation detector for acquiring physical measurements, fabricated from an inorganic nano-crystalline scintillator that has been fixed on an optical fiber terminus. This dosimeter allows for the measurement of point doses to sub-millimeter resolution, and has the ability to be placed in-vivo in humans and small animals. Real-time data is displayed to the user to provide instant quality assurance and dose-rate information. The second tool utilizes an open source Monte Carlo particle transport code, and was applied for small animal dosimetry studies to calculate organ doses and recommend new techniques of dose prescription in mice, as well as to characterize dose to the murine bone marrow compartment with micron-scale resolution.

Hardware design changes were implemented to reduce the overall fiber diameter to <0.9 mm for the nano-crystalline scintillator based fiber optic detector (NanoFOD) system. Lower limits of device sensitivity were found to be approximately 0.05 cGy/s. Herein, this detector was demonstrated to perform quality assurance of clinical 192Ir HDR brachytherapy procedures, providing comparable dose measurements as thermo-luminescent dosimeters and accuracy within 20% of the treatment planning software (TPS) for 27 treatments conducted, with an inter-quartile range ratio to the TPS dose value of (1.02-0.94=0.08). After removing contaminant signals (Cerenkov and diode background), calibration of the detector enabled accurate dose measurements for vaginal applicator brachytherapy procedures. For 192Ir use, energy response changed by a factor of 2.25 over the SDD values of 3 to 9 cm; however a cap made of 0.2 mm thickness silver reduced energy dependence to a factor of 1.25 over the same SDD range, but had the consequence of reducing overall sensitivity by 33%.

For preclinical measurements, dose accuracy of the NanoFOD was within 1.3% of MOSFET measured dose values in a cylindrical mouse phantom at 225 kV for x-ray irradiation at angles of 0, 90, 180, and 270˝. The NanoFOD exhibited small changes in angular sensitivity, with a coefficient of variation (COV) of 3.6% at 120 kV and 1% at 225 kV. When the NanoFOD was placed alongside a MOSFET in the liver of a sacrificed mouse and treatment was delivered at 225 kV with 0.3 mm Cu filter, the dose difference was only 1.09% with use of the 4x4 cm collimator, and -0.03% with no collimation. Additionally, the NanoFOD utilized a scintillator of 11 µm thickness to measure small x-ray fields for microbeam radiation therapy (MRT) applications, and achieved 2.7% dose accuracy of the microbeam peak in comparison to radiochromic film. Modest differences between the full-width at half maximum measured lateral dimension of the MRT system were observed between the NanoFOD (420 µm) and radiochromic film (320 µm), but these differences have been explained mostly as an artifact due to the geometry used and volumetric effects in the scintillator material. Characterization of the energy dependence for the yttrium-oxide based scintillator material was performed in the range of 40-320 kV (2 mm Al filtration), and the maximum device sensitivity was achieved at 100 kV. Tissue maximum ratio data measurements were carried out on a small animal x-ray irradiator system at 320 kV and demonstrated an average difference of 0.9% as compared to a MOSFET dosimeter in the range of 2.5 to 33 cm depth in tissue equivalent plastic blocks. Irradiation of the NanoFOD fiber and scintillator material on a 137Cs gamma irradiator to 1600 Gy did not produce any measurable change in light output, suggesting that the NanoFOD system may be re-used without the need for replacement or recalibration over its lifetime.

For small animal irradiator systems, researchers can deliver a given dose to a target organ by controlling exposure time. Currently, researchers calculate this exposure time by dividing the total dose that they wish to deliver by a single provided dose rate value. This method is independent of the target organ. Studies conducted here used Monte Carlo particle transport codes to justify a new method of dose prescription in mice, that considers organ specific doses. Monte Carlo simulations were performed in the Geant4 Application for Tomographic Emission (GATE) toolkit using a MOBY mouse whole-body phantom. The non-homogeneous phantom was comprised of 256x256x800 voxels of size 0.145x0.145x0.145 mm3. Differences of up to 20-30% in dose to soft-tissue target organs was demonstrated, and methods for alleviating these errors were suggested during whole body radiation of mice by utilizing organ specific and x-ray tube filter specific dose rates for all irradiations.

Monte Carlo analysis was used on 1 µm resolution CT images of a mouse femur and a mouse vertebra to calculate the dose gradients within the bone marrow (BM) compartment of mice based on different radiation beam qualities relevant to x-ray and isotope type irradiators. Results and findings indicated that soft x-ray beams (160 kV at 0.62 mm Cu HVL and 320 kV at 1 mm Cu HVL) lead to substantially higher dose to BM within close proximity to mineral bone (within about 60 µm) as compared to hard x-ray beams (320 kV at 4 mm Cu HVL) and isotope based gamma irradiators (137Cs). The average dose increases to the BM in the vertebra for these four aforementioned radiation beam qualities were found to be 31%, 17%, 8%, and 1%, respectively. Both in-vitro and in-vivo experimental studies confirmed these simulation results, demonstrating that the 320 kV, 1 mm Cu HVL beam caused statistically significant increased killing to the BM cells at 6 Gy dose levels in comparison to both the 320 kV, 4 mm Cu HVL and the 662 keV, 137Cs beams.

An integrated approach to flow, thermal and mechanical modeling of electronics devices

The future success of many electronics companies will depend to a large extent on their ability to initiate techniques that bring schedules, performance, tests, support, production, life-cycle-costs, reliability prediction and quality control into the earliest stages of the product creation process. Earlier papers have discussed the benefits of an integrated analysis environment for system-level thermal, stress and EMC prediction. This paper focuses on developments made to the stress analysis module and presents results obtained for an SMT resistor. Lifetime predictions are made using the Coffin-Manson equation. Comparison with the creep strain energy based models of Darveaux (1997) shows the shear strain based method to underestimate the solder joint life. Conclusions are also made about the capabilities of both approaches to predict the qualitative and quantitative impact of design changes.