Image reconstruction by an alternating minimisation

This paper focuses on the problem of incomplete data in the applications of the circular cone-beam computed tomography. This problem is frequently encountered in medical imaging sciences and some other industrial imaging systems. For example, it is crucial when the high density region of objects can only be penetrated by X-rays in a limited angular range. As the projection data are only available in an angular range, the above mentioned incomplete data problem can be attributed to the limited angle problem, which is an ill-posed inverse problem. This paper reports a modified total variation minimisation method to reduce the data insufficiency in tomographic imaging. This proposed method is robust and efficient in the task of reconstruction by showing the convergence of the alternating minimisation method. The results demonstrate that this new reconstruction method brings reasonable performance. (C) 2010 Elsevier B.V. All rights reserved.

Comparative study of Ga2O3 : Dy3+ phosphors prepared by three methods

Trivalent dysprosium (Dy3+)-activated beta-gallium oxide (beta-Ga2O3) phosphors were synthesized by solid-state (SS), coprecipitation (CP), and citrate gel (CG) methods, respectively. The resulting beta-Ga2O3:Dy3+ phosphors were well characterized by X-ray diffraction, field-emission scanning electron microscopy (FESEM), and by photoluminescence and cathodoluminescence spectra. The phosphors prepared by different methods show different luminescence properties under the excitation of UV and low-voltage cathode rays. There exists an energy transfer from the beta-Ga2O3 host lattices to Dy3+, and the energy transfer efficiency was higher in the CG-derived phosphors than those of the SS- and CP-derived phosphors. Factors influencing the intensity of luminescence and the energy transfer efficiency from beta-Ga2O3 host to Dy3+ in beta-Ga2O3:Dy3+ phosphors were investigated in detail. It is shown that the CG method is the most respected process for the preparation of beta-Ga2O3:Dy3+ phosphors.

High strength Ni based composite reinforced by solid solution W(Al) obtained by powder metallurgy

The solid-solution-particle reinforced W(Al)-Ni composites were successfully fabricated by using mechanical alloying (MA) and hot-pressing (HP) technique when the content of Ni is between 45 wt% and 55 wt%. Besides, samples of various original component ratio of Al50W50 to Ni have been fabricated, and the corresponding microcomponents and mechanical properties such as microhardness, ultimate tensile strength and elongation were characterized and discussed. The optimum ultimate tensile strength under the experiment conditions is 1868 MPa with elongation of 10.21 % and hardness of 6.62 GPa. X-ray diffraction (XRD), FE-SEM and energy dispersive analysis of X-rays (EDS) were given to analysis the components and morphology of the composite bulk specimens.

Thermal analysis of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) irradiated under vacuum

Poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) was irradiated by Co-60 gamma-rays (doses of 50, 100 and 200kGy) under vacuum. The thermal analysis of control and irradiated PHBV, under vacuum was carried out by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). The tensile properties of control and irradiated PHBV were examined by using an Instron tensile testing machine. In the thermal degradation of control and irradiated PHBV, a one-step weight loss was observed. The derivative thermogravimetric curves of control and irradiated PHBV confirmed only one weight-loss step change. The onset degradation temperature (T-o) and the temperature of maximum weight-loss rate (T-p) of control and irradiated PHBV were in line with the heating rate (degreesC min(-1)). T-o and T-p of PHBV decreased with increasing radiation dose at the same heating rate. The DSC results showed that Co-60 gamma-radiation significantly affected the thermal properties of PHBV. With increasing radiation dose, the melting temperature (T-m) of PHBV shifted to a lower value, due to the decrease in crystal size. The tensile strength and fracture strain of the irradiated PHBV decreased, hence indicating an increased brittleness.

LDPE/carbon black conductive composites: Influence of radiation crosslinking on PTC and NTC properties

The electrical resistivity of low-density polyethylene/carbon black composites irradiated by Co-60 gamma-rays was investigated as a function of temperature. The experimental results obtained by scanning electron microscopy, solvent extraction techniques, and pressure-specific volume-temperature analysis techniques showed that the positive temperature coefficient (PTC) and negative temperature coefficient (NTC) effects of the composites were influenced by the irradiation dose, network forming (gel), and soluble fractions (Sol). The NTC effect was effectively eliminated when the radiation dose reached 400 kGy. The results showed that the elimination of the NTC effect was related to the difference in the thermal expansion of the gel and Sol regions. The thermal expansion of the sol played an important role in both increasing the PTC intensity and decreasing the NTC intensity at 400 kGy.

Effect of cross-linking of high-density polyethylene. I. On spherulitic structures

Crystallization behavior and spherulitic structure of linear high-density polyethylene (HDPE), after being irradiated in its molten state by gamma -rays, was investigated by small-angle laser scattering (SALS) and differential scanning calorimetry (DSC). Significant changes in the crystallization of HDPE during cooling in air before and after being irradiated in the melt were observed. A critical minimum average molar mass between cross-links (200 carbon-carbon bonds) for spherulite formation in such an irradiated HDPE network was obtained.

Production of active oxygen species from Taxus cuspidata induced by far-UV radiation

In order to understand the role of active oxygen species in mediating plant injuries induced by far-UV radiation, seedlings of Taxus cuspidata Sieb. et Zucc. were irradiated by far-UV rays in laboratory for 4 weeks. The production of organic free-radicals in detached needles, and the production of O-2(radical anion) and O-1(2) in isolated chloroplasts were detected weekly by electron spin resonance (ESR) to evaluate their relative importance. The results show that the cumulative effect of far-UV irradiation, is best indicated by the production of organic free radicals in the needles, O-2(radical anion) production in chloroplasts is the next. The enhancement of O-1(2) production in chloroplasts by the cumulative far-UV irradiation seems to be not so important as O-2(radical anion) in mediating injuries induced by, far-UV radiation because of its high background value.

RADIATION EFFECTS ON POLY(VINYLIDENE FLUORIDE)

Melt-crystallized poly(vinylidene fluoride)s (PVF2) with different crystallization histories were irradiated with gamma-rays within the range of irradiation doses 0-83 Mrad. The effects on the crystalline structure and mechanical properties have been measured, compared, and discussed. The degree of crystallinity of the samples was found to increase with radiation dose. The differential scanning calorimeter scans of the quenched samples indicate that there are two melting peaks, and that the area of the lower temperature peak increases while the area of higher temperature peak decreases with increasing dose. Yield stress and breaking stress for all samples are not significantly affected by irradiation but elongation at break is.

STUDIES ON RADIATION EFFECTS OF PTFE .2. RADIATION-INDUCED STRUCTURAL-CHANGES AND THE TEMPERATURE EFFECT

In this work PTFE sheets irradiated with gamma-rays at 150-degrees-C and 200-degrees-C were studied using x-ray photoelectron spectroscopy (XPS). The main structural changes in PTFE due to radiation are the formation of CF3 and CF groups. An irradiation temperature dependence of the relative content of the three kinds of groups in irradiated PTFE was observed. The CF3 groups, especially when irradiation is carried out a lower temperatures, can defluorinate in the same manner as previosly reported for CF2 groups. The CF groups, on the other hand, are observed to increase with increasing irradiation dose and irradiation temperature; the latter was explained as due to an increase in branching structures.

大云室组动量选择器

采用“预选—寄存—运算再选”的方案,设计并完成了云南站大云室组的宇宙线事例选择系统.其中动量选择器应用了数字在线计算技术.设计了小型专用计算单元,根据荷电粒子在磁场中的偏转角与动量的对应关系,对寄存下来的描迹信号进行快速运算和逻辑判断,从而实现单个荷电粒子的动量选择.

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia.

井间地震速度层析成像方法与软件系统

Crosshole Seismic tomography has been broadly studied and applied in the fields of resource exploration and engineering exploration because of its special observing manner and better resolution than normal seismic exploration. This thesis will state the theory and method of Crosshole Seismic tomography. Basing on the previous studies，the thesis studied the initial velocity model，ray-tracing method, and developed the three-dimension tomography software. All the cells that a ray passes through are of the same velocities if the paths from transmitters to receivers are straight. The cells that the each ray passes through are recorded, and rays that pass through each cell are calculated. The ray average velocity which passes through a cell is set as the cell velocity. Analogously we can make a initial node velocity model because the velocity sum is calculated on the all cells which own to a certain node, and the cell number is summed about each nodes，the ratio of the velocity sum to the all cells number is set as the node velocity. The inversion result from the initial node velocity model is better than that of the average velocity model. Ray-bending and Shortest Path for Rays (SPR) have shortcomings and limitations respectively. Using crooked rays obtained from SPR rather than straight lines as the starting point can not only avoid ray bending converging to the local minimum travel time path, but also settle the no smooth ray problem obtained by SPR. The hybrid method costs much computation time, which is roughly equal to the time that SPR expends. The Delphi development tool based on the Object Pascal language standard has an advantage of object-oriented. TDTOM (Three Dimensions Tomography) was developed by using Delphi from the DOS version. Improvement on the part of inversion was made, which bring faster convergence velocity. TDTOM can be used to do velocity tomography from the first arrival travel time of the seismic wave, and it has the good qualities of friendly user interface and convenient operation. TDTOM is used to reconstruct the velocity image for a set of crosshole data from Karamay Oil Field. The geological explanation is then given by comparing the inversion effects of different ray-tracing methods. High velocity zones mean the cover of oil reservoir, and low velocity zones correspond to the reservoir or the steam flooding layer.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

铂族元素分析方法研究及其在有关地质问题中的应用

Characterization of Platinum Group Elements (PGE) has been applied to earth, space and environmental sciences. However, all these applications are based on a basic prerequisite, i.e. their concentration or ratio in the research objects can be accurately and precisely determined. In fact, development in these related studies is a great challenge to the analytical chemistry of the PGE because their content in the geological sample (non-mineralized) is often extremely low, range from ppt (10~(-12)g/g) to ppt (10~(-9)g/g). Their distribution is highly heterogeneous, usually concentrating in single particle or phase. Therefore, the accurate determination of these elements remains a problem in analytical chemistry and it obstructs the research on geochemistry of PGE. A great effort has been made in scientific community to reliable determining of very low amounts of PGE, which has been focused on to reduce the level of background in used reagents and to solve probable heterogeneity of PGE in samples. Undoubtedly, the fire-assay method is one of the best ways for solving the heterogeneity, as a large amount of sample weight (10-50g) can be hold. This page is mainly aimed at development of the methodology on separation, concentration and determination of the ultra-trace PGE in the rock and peat samples, and then they are applied to study the trace of PGE in ophiolite suite, in Kudi, West Kunlun and Tunguska explosion in 1908. The achievements of the study are summarized as follows: 1. A PGE lab is established in the Laboratory of Lithosphere Tectonic Evolution, IGG, CAS. 2. A modified method of determination of PGE in geological samples using NiS Fire-Assay with inductively coupled plasma-mass spectrometry (ICP-MS) is set up. The technical improvements are made as following: (1) investigating the level of background in used reagents, and finding the contents of Au, Pt and Pd in carbonyl nickel powder are 30, 0.6 and 0.6ng/g, respectively and 0.35, 7.5 and 6.4ng, respectively in other flux, and the contents of Ru, Rh, Os in whole reagents used are very low (below or near the detection limits of ICP-MS); (2) measuring the recoveries of PGE using different collector (Ni+S) and finding 1.5g of carbonyl nickel is effective for recovering the PGE for 15g samples (recoveries are more than 90%), reducing the inherent blank value due to impurities reagents; (3) direct dissolving nickel button in Teflon bomb and using Te-precipitation, so reducing the loss of PGE during preconcentration process and improving the recoveries of PGE (above 60% for Os and 93.6-106.3% for other PGE, using 2g carbonyl nickel); (4) simplifying the procedure of analyzing Osmium; (5)method detection limits are 8.6, 4.8, 43, 2.4, 82pg/g for 15g sample size ofRu, Rh, Pd, Ir, Pt, respectively. 3. An analytical method is set up to determine the content of ultra-trace PGE in peat samples. The method detection limits are 0.06, 0.1, 0.001, 0.001 and 0.002ng/mL for Ru, Rh, Pd, Ir and Pt, respectively. 4. Distinct anomaly of Pd and Os are firstly found in the peat sampling near the Tunguska explosion site, using the analytical method. 5. Applying the method to the study on the origin of Tunguska explosion and making the following conclusions: (1) these excess elements were likely resulted from the Tunguska Cosmic Body (TCB) explosion of 1908. (2) The Tunguska explosive body was composed of materials (solid components) similar to C1 chondrite, and, most probably, a cometary object, which weighed more than 10~7 tons and had a radius of more than 126 m. 6. The analysis method about ultra-trace PGE in rock samples is successfully used in the study on the characteristic of PGE in Kudi ophiolite suite and the following conclusions are made: (1) The difference of the mantle normalization of PGE patterns between dunite, harzburgite and lherzolite in Kudi indicates that they are residual of multi-stage partial melt of the mantle. Their depletion of Ir at a similar degree probably indicates the existence of an upper mantle depleted Ir. (2) With the evolution of the magma produced by the partial melt of the mantle, strong differentiation has been shown between IPGE and PPGE; and the differentiation from pyroxenite to basalt would have been more and more distinct. (3) The magma forming ophiolite in Kudi probably suffered S-saturation process.