空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

KEL-TI介质的地震波场响应与成像方法研究

To deal with the problems in multi-component converted seismic wave exploration in coal fields, the wave propagating features and imaging methods of multi-component converted waves in coal measure strata are researched in this thesis firstly. The relations between viscoelasticity and anisotropy in coal measure strata are analyzed to build KEL-TI model, and which seismic wave propagating and attenuating features are researched. The disadvantages of converted wave imaging methods based on common converted point gather are analyzed and constant velocity no NMO converted wave imaging method based on common scattering point gather is put forward, according to Huygens-Fresnel principle, which applicabilities in the elastic isotropic, elastic TI and KEL-TI situations are discussed. To different model simulation data, the common scattering point gathers’ and stacked profiles’ features are analyzed. The results show that the method can image compressional waves and converted waves with high precision. Secondly, the resolution enhancing theories and methods of converted wave are researched by Rayleigh wave suppressing, converted wave static correction and poststack inverse-Q filtering. 1) The polarization filter is designed by the instantaneous polarization information of seismic waves, and the Rayleigh wave suppressing method is researched. From the spectrum analysis before and after filtering, it can be derived that the amplitudes are kept relatively. 2) To constant velocity no NMO converted wave imaging method, the static correction method based on common equivalent offset point gather is put forward and tested to the actual converted waves. 3) The relation between equivalent quality factor of converted wave, compressional wave quality factor and the ratio of compressional to shear wave velocity is derived. The compressional wave quality inversion method by first arrivals of none-offset VSP is researched, and which is then transformed to the equivalent quality factor to perform inverse-Q filtering of actual converted waves. The result has shown that the method can recover the high frequency energy of converted waves. At last, the theories and methods researched in this thesis are practiced to the 3D3C seismic exploration in Guqiao coal mine in Huainan and achieve good results.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

生物－热催化过渡带沉积成岩演化及其成烃作用研究

The biothermocatalytic transitional zone gas is a new type of natural gas genetic theory, and also an clean, effective and high quality energy with shallow burial depth, wide distribution and few investment. Meanwhile, this puts biothermocatalytic transitional zone gas in important position to the energy resource and it is a challenging front study project. This paper introduces the concept, the present situation of study and developmental trend about biothermocatalytic transitional zone gas in detail. Then by using heat simulating of source rocks and catalysis mechanism analysis in the laboratory and studying structural evolution, sedimentation, diagenesis and the conditions of accumulation formation and so on, this paper also discusses catalytic mechanism and evolutionary model of the biothermocatalytic transitional zone gas formation, and establishes the methods of appraisal parameter and resources prediction about the biothermocatalytic transitional zone gas. At last, it shows that geochemical characteristics and differentiated mark of the biothermocatalytic transitional zone gas, and perfect natural gas genetic theory, and points out the conditions of accumulation formation, distribution characteristics and potential distribution region on the biothermocatalytic transitional zone gas m China. The paper mainly focuses on the formation mechanism and the resources potential about the biothermocatalytic transitional zone gas. Based on filed work, it is attached importance to a combination of macroscopic and microcosmic analysis, and the firsthand data are obtained to build up framework and model of the study by applying geologic theory. Based on sedimentary structure, it is expounded that structural actions have an effect on filling space and developmental cource of sediments and evolution of source rocks. Carried out sedimentary environment, sequence stratigraphy, sedimentary system and diagenesis and so on, it is concluded that diagenesis influences developmental evolution of source rocks, and basic geologic conditions of the biothermocatalytic transitional zone gas. Applying experiment simulating and catalytic simulating as well as chemical analysis, catalytic mechanism of clay minerals is discussed. Combined diagenecic dynamics with isotope fractionation dynamics, it is established that basis and method of resource appraisal about the biothermocatalytic transitional zone gas. All these results effectively assess and predict oil&gas resources about the biothermocatalytic transitional zone gas-bearing typical basin in China. I read more than 170 volumes on the biothermocatalytic transitional zone gas and complete the dissertation' summary with some 2.4 ten thousand words, draw up study contents in some detail and set up feasible experimental method and technologic course. 160 pieces of samples are obtained in oilfield such as Liaohe, Shengli, Dagang and Subei and so on, some 86 natural gas samples and more than 30 crude oil samples. Core profiles about 12 wells were observed and some 300 geologic photos were taken. Six papers were published in the center academic journal at home and abroad. Collected samples were analysised more than 1000 times, at last I complete this dissertation with more than 8 ten thousand words, and with 40 figures and 4 plates. According to these studies, it is concluded the following results and understandings. 1. The study indicates structural evolution and action of sedimentary basin influence and control the formation and accumulation the biothermocatalytic transitional zone gas. Then, the structural action can not only control accommodation space of sediments and the origin, migration and accumulation of hydrocarbon matters, but also can supply the origin of energy for hygrocarbon matters foramtion. 2. Sedimentary environments of the biothermocatalytic transitional zone gas are lake, river and swamp delta- alluvial fan sedimentary systems, having a warm, hot and humid climate. Fluctuation of lake level is from low to high., frequency, and piling rate of sedimentary center is high, which reflect a stable depression and rapidly filling sedimentary course, then resulting in source rocks with organic matter. 3. The paper perfects the natural gas genetic theory which is compound and continuous. It expounds the biothermocatalytic transitional zone gas is a special gas formation stage in continuous evolutionary sequence of organic matter, whose exogenic force is temperture and catalysis of clay minerals, at the same time, having decarbxylation, deamination and so on. 4. The methodology is established which is a combination of SEM, TEM and Engery spectrum analysis to identify microstructure of crystal morphology about clay minerals. Using differential thermal-chromatographic analysis, it can understand that hydrocarbon formation potential of different typies kerogens and catalytic method of all kinds of mineral matrix, and improve the surface acidity technology of clay minerals measured by the pyridine analytic method. 5. The experiments confirm catalysis of clay minerals to organic matter hygrocarbon formation. At low temperature (<300 ℃), there is mainly catalysis of montmorillonite, which can improve 2-3 times about produced gas of organic matters and the pyrolyzed temperature decreased 50 ℃; while at the high temperature, there is mainly catalysis of illite which can improve more than 2 times about produced gas of organic matters. 6. It is established the function relationship between organic matter (reactant) concentration and temperature, pressure, time, water and so on, that is C=f (D, t). Using Rali isotope fractionation effect to get methane isotope fractionation formula. According to the relationship between isotope fractionation of diagenesis and depth, and combined with sedimentary rate of the region, it is estimated that relict gas of the biothermocatalytic transitional zone gas in the representative basin. 7. It is revealed that hydrocarbon formation mechanism of the biothermocatalytic transitional zone gas is mainly from montmorillonite to mixed minerals during diagenesis. In interlayer, a lot of Al~(3+) substitute for Si~(4+), resulting in a imbalance between surface charge and interlayer charge of clay minerals and the occurrence of the Lewis and Bronsted acid sites, which promote to form the carbon cation. The cation can form alkene or small carbon cation. 8. It is addressed the comprehensive identification mark of the biothermo - catalytic transitional zone gas. In the temproal-spatial' distribution, its source rocks is mainly Palaeogene, secondly Cretaceous and Jurassic of Mesozoic, Triassic, having mudy rocks and coal-rich, their organic carbon being 0.2% and 0.4% respectively. The vitrinite reflection factor in source rocks Ro is 0.3-0.65%, a few up to 0.2%. The burial depth is 1000-3000m, being characterized by emerge of itself, reservoir of itself, shallow burial depth. In the transitional zone, from shallow to deep, contents of montmorillonites are progressively reduced while contents of illites increasing. Under SEM, it is observed that montmorillonites change into illite.s, firstly being mixed illite/ montmorillonite with burr-like, then itlite with silk-like. Carbon isotope of methane in the biothermocatatytic transitional zone gas , namely δ~(13)C_1-45‰- -60 ‰. 9. From the evolutionary sequence of time, distribution of the biothermocatalytic transitional zone gas is mainly oil&gas bearing basin in the Mesozoic-Neozoic Era. From the distribution region, it is mainly eastern stuctural active region and three large depressions in Bohaiwang basin. But most of them are located in evolutionary stage of the transitional zone, having the better relationship between produced, reservoir and seal layers, which is favorable about forming the biothermocatalytic transitional zone gas reservoir, and finding large gas (oil) field.

In situ study of carbon nanotube formation by C2H2 decomposition on an iron-based catalyst

The properties and formation of nanotubes have been extensively studied, but very few deal with the catalytic production mechanism of nanotubes. Two different techniques, thermogravimetric analysis and UV-Raman, have been applied to analyse the carbon deposition by catalysed decomposition of acetylene over an iron-based catalyst. The nature of the produced carbon materials depends on reaction temperature. Also, TEM allows identification of carbon nanotubes, encapsulated particles, and other nanostructures, while UV-Raman confirms its graphitic and graphite-like nature. (C) 2000 Elsevier Science Ltd. All rights reserved.

Sensitivity map of laser tweezers Raman spectroscopy for single-cell analysis of colorectal cancer

Raman spectroscopy on single, living epithelial cells captured in a laser trap is shown to have diagnostic power over colorectal cancer. This new single-cell technology comprises three major components: primary culture processing of human tissue samples to produce single-cell suspensions, Raman detection on singly trapped cells, and diagnoses of the cells by artificial neural network classifications. it is compared with DNA flow cytometry for similarities and differences. Its advantages over tissue Raman spectroscopy are also discussed. In the actual construction of a diagnostic model for colorectal cancer, real patient data were taken to generate a training set of 320 Raman spectra and, a test set of 80. By incorporating outlier corrections to a conventional binary neural classifier, our network accomplished significantly better predictions than logistic regressions, with sensitivity improved from 77.5% to 86.3% and specificity improved from 81.3% to 86.3% for the training set and moderate improvements for the test set. Most important, the network approach enables a sensitivity map analysis to quantitate the relevance of each Raman band to the normal-to-cancer transform at the cell level. Our technique has direct clinic applications for diagnosing cancers and basic science potential in the study of cell dynamics of carcinogenesis. (C) 2007 Society of Photo-Optical Instrumentation Engineers.

Sensitivity map of laser tweezers Raman spectroscopy for single-cell analysis of colorectal cancer

Raman spectroscopy on single, living epithelial cells captured in a laser trap is shown to have diagnostic power over colorectal cancer. This new single-cell technology comprises three major components: primary culture processing of human tissue samples to produce single-cell suspensions, Raman detection on singly trapped cells, and diagnoses of the cells by artificial neural network classifications. it is compared with DNA flow cytometry for similarities and differences. Its advantages over tissue Raman spectroscopy are also discussed. In the actual construction of a diagnostic model for colorectal cancer, real patient data were taken to generate a training set of 320 Raman spectra and, a test set of 80. By incorporating outlier corrections to a conventional binary neural classifier, our network accomplished significantly better predictions than logistic regressions, with sensitivity improved from 77.5% to 86.3% and specificity improved from 81.3% to 86.3% for the training set and moderate improvements for the test set. Most important, the network approach enables a sensitivity map analysis to quantitate the relevance of each Raman band to the normal-to-cancer transform at the cell level. Our technique has direct clinic applications for diagnosing cancers and basic science potential in the study of cell dynamics of carcinogenesis. (C) 2007 Society of Photo-Optical Instrumentation Engineers.

Judd-Oflet analysis of spectrum and laser performance of Ho:YAP crystal end-pumped by 1.91-mu m Tm:YLF laser

The Ho:YAP crystal is grown by the Czochralski technique. The room-temperature polarized absorption spectra of Ho:YAP crystal was measured on a c-cut sample with 1 at% holmium. According to the obtained Judd-Ofelt intensity parameters Omega(2) = 1.42 x 10(-20) cm(2), Omega(4) = 2.92 x 10(-20) cm(2), and Omega(6) = 1.71 x 10(-20) cm(2), this paper calculated the fluorescence lifetime to be 6 ms for I-5(7) -> I-5(8) transition, and the integrated emission cross section to be 2.24 x 10(-18) cm(2). It investigates the room-temperature Ho:YAP laser end-pumped by a 1.91-mu m Tm:YLF laser. The maximum output power was 4.1 W when the incident 1.91-mu m pump power was 14.4W. The slope efficiency is 40.8%, corresponding to an optical-to-optical conversion efficiency of 28.4%. The Ho:YAP output wavelength was centred at 2118 nm with full width at half maximum of about 0.8 nm.

Fast and effective analysis methods for fiber Raman amplifiers

A series of new single-step methods and their corresponding algorithms with automatic step size adjustment for model equations of fiber Raman amplifiers are proposed and compared in this paper. On the basis of the Newton-Raphson method, multiple shooting algorithms for the two-point boundary value problems involved in solving Raman amplifier propagation equations are constructed. A verified example shows that, compared with the traditional Runge-Kutta methods, the proposed methods can increase the accuracy by more than two orders of magnitude under the same conditions. The simulations for Raman amplifier propagation equations demonstrate that our methods can increase the computing speed by more than 5 times, extend the step size significantly, and improve the stability in comparison with the Dormand-Prince method. The numerical results show that the combination of the multiple shooting algorithms and the proposed methods has the capacity to rapidly and effectively solve the model equations of multipump Raman amplifiers under various conditions such as co-, counter- and bi-directionally pumped schemes, as well as dual-order pumped schemes.

Effect of rapid thermal annealing on the Raman spectrum of Si0.33Ge0.67/Si (100) alloy

The effect of thermal annealing on the Raman spectrum of Si0.33Ge0.67 alloy grown on Si (100) by molecular beam epitaxy is investigated in the temperature range of 550-800 degrees C. For annealing below 700 degrees C, interdiffusion at the interface is negligible and the residual strain plays the dominant role in the Raman shift. The strain-shift coefficients for Si-Ge and Ge-Ge phonon modes are determined to be 915 +/- 215 cm(-1) and 732 +/- 117 cm(-1), respectively. For higher temperature annealing, interdiffusion is significant and strongly affects the Raman shift and the spectral shape.

Analysis of the xuv photoabsorption spectrum of Au2+, Au3+, and Au4+

The photoabsorption processes of Au2+, Au3+, and Au4+ have been investigated experimentally and theoretically in the 70-127 eV region. Using the dual laser-produced plasma technique, the 4f and 5p photoabsorption spectrum has been recorded at 50 ns time delay and was found to be dominated by a great number of lines from 4f-5d, 6d and 5p-5d, 6s transitions, which have been identified by comparison with the aid of Hartree-Fock with configuration interaction calculations. The characteristic feature of the spectrum is that satellite lines from excited configurations containing one or two 6s electrons are more important than resonance lines, and with increasing ionization, satellite contributions from states with one 6s spectator electron gradually become more important than those with two 6s spectator electrons. Based on the assumption of a normalized Boltzmann distribution among the excited states and a steady-state collisional-radiative model, we succeeded in reproducing a spectrum which is in good agreement with experiment.

Fabrication, characterization, and application in surface-enhanced Raman spectrum of assembled type-I collagen-silver nanoparticle multilayered films

In this paper, we report a facile method for the fabrication of type-I collagen-silver nanoparticles (Ag NPs) multilayered films by utilizing type-I collagen as a medium. These samples were characterized by UV-vis spectra photometer, atomic force microscopy, scanning electron microscopy, and Fourier transform IR spectrum. Experimental results show that collagen molecules serve as effective templates to assemble Ag NPs into multilayer films. These samples exhibit high surface-enhanced Raman scattering (SERS) enhancement abilities.

Impedance spectrum measurement in liquid liquid interface electrochemistry .1. Analysis of various kinds of interfacial impedance in the four-electrode system

Applying impedance spectrum technique to liquid/liquid interfacial electrochemistry, we present the theoretical expression of the liquid/liquid interfacial impedance in the four - electrode electrolytic cell measured by alternating current impedance method. The influence of the electrolytic cell parameters and input passage parameters of the impedance apparatus on impedance spectrum curves are theoretically studied.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.