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

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.

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.

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.

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.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Cauchy singular integral equation method for transient antiplane dynamic problems

In this paper, by use of the boundary integral equation method and the techniques of Green basic solution and singularity analysis, the dynamic problem of antiplane is investigated. The problem is reduced to solving a Cauchy singular integral equation in Laplace transform space. This equation is strictly proved to be equivalent to the dual integral equations obtained by Sih [Mechanics of Fracture, Vol. 4. Noordhoff, Leyden (1977)]. On this basis, the dynamic influence between two parallel cracks is also investigated. By use of the high precision numerical method for the singular integral equation and Laplace numerical inversion, the dynamic stress intensity factors of several typical problems are calculated in this paper. The related numerical results are compared to be consistent with those of Sih. It shows that the method of this paper is successful and can be used to solve more complicated problems. Copyright (C) 1996 Elsevier Science Ltd

Singular behaviour near the tip of a sharp V-notch in a power law hardening material

This paper presents an asymptotic analysis of the near-tip stress and strain fields of a sharp V-notch in a power law hardening material. First, the asymptotic solutions of the HRR type are obtained for the plane stress problem under symmetric loading. It is found that the angular distribution function of the radial stress sigma(r) presents rapid variation with the polar angle if the notch angle beta is smaller than a critical notch angle; otherwise, there is no such phenomena. Secondly, the asymptotic solutions are developed for antisymmetric loading in the cases of plane strain and plane stress. The accurate calculation results and the detailed comparisons are given as well. All results show that the singular exponent s is changeable for various combinations of loading condition and plane problem.

NK-lysin of channel catfish: Gene triplication, sequence variation, and expression analysis

Antimicrobial peptides (AMPs) are important components of the host innate immune response against microbial invasion. In addition to the previously known four classes of antimicrobial peptides, a fifth class of antimicrobial peptides has been recently identified to include NK-lysins that have a globular three-dimensional structure and are larger with 74-78 amino acid residues. NK-lysin has been shown to harbor antimicrobial activities against a wide spectrum of microorganisms including bacteria, fungi, protozoa, and parasites. To date, NK-lysin genes have been reported from only a limited number of organisms. We previously identified a NK-lysin cDNA in channel catfish. Here we report the identification of two noveltypes of NK-lysin transcripts in channel catfish. Altogether, three distinct NK-lysin transcripts exist in channel catfish. In this work, their encoding genes were identified, sequenced, and characterized. We provide strong evidence that the catfish NK-lysin gene is tripled in the same genomic neighborhood. All three catfish NK-lysin genes are present in the same genomic region and are tightly linked on the same chromosome, as the same BAC clones harbor all three copies of the NK-lysin genes. All three NK-lysin genes are expressed, but exhibit distinct expression profiles in various tissues. In spite of the existence of a single copy of NK-lysin gene in the human genome, and only a single hit from the pufferfish,genome, there are two tripled clusters of NK-lysin genes on chromosome 17 of zebrafish in addition to one more copy on its chromosome 5. The similarity in the genomic arrangement of the tripled NK-lysin genes in channel catfish and zebrafish suggest similar evolution of NK-lysin genes. (c) 2005 Elsevier Ltd. All rights reserved.

Lineshape Analysis of the Beat Signal Between Optical Carrier and Delayed Sidebands

This paper presents the lineshape analysis of the beat signal between the optical carrier and the shifted and delayed side-bands produced by sinusoidal amplitude modulation. It is shown that the beat signal has a typical lineshape with a very narrow delta-peak superposed on a quasi-Lorentzian profile. Theoretical explanation for the appearance of this peak has been given based on optical spectral structure constructed by a large number of optical wave trains. It is predicted that the delta-peak is originated from the beat between the wave trains in the carrier and those in the delayed sidebands when their average coherence length is longer than the delay line. Experiments carried out using different delay lines clearly show that the delta-peak is always located at the modulation frequency and decreases with the increasing delay line. Our analysis explicitly indicates that the linewidth is related to the observation time. It is also suggested that the disappearance of the delta-peak can be used as the criterion of coherence elimination.

Numerical analysis of gain saturation, noise figure, and carrier distribution for quantum-dot semiconductor-optical amplifiers

The gain saturation behaviors and noise figure are numerically analyzed for quantum-dot semiconductor optical amplifiers (QD-SOAs). The carrier and photon distributions in the longitudinal direction as well as the photon energy dependent facet reflectivity are accounted in the rate equations, which are solved with output amplified spontaneous emission spectrum as iterative variables. The longitudinal distributions of the occupation probabilities and spectral-hole burning are presented for electrons in the excited and ground states of quantum dots. The saturation output power 19.7 dBm and device gain 20.6 dB are obtained for a QD-SOA with the cavity length of 6 rum at the bias current of 500 mA. The influences of them electron intradot relaxation time and the QD capture time on the gain spectrum are simulated with the relaxation time of 1, 30, and 60 ps and capture time of 1, 5, and 10 ps. The noise figure as low as 3.5 dB is expected due to the strong polarization sensitive spontaneous emission. The characteristics of gain saturation and noise figure versus input signal power for QD-SOAs are similar to that of semiconductor. linear optical amplifiers with gain clamping by vertical laser fields.

Analysis of the additional modulation of semiconductor lasers using optical heterodyne techniques

In this letter, the power spectrum of a cooled distributed feedback laser module is measured using the self-heterodyne technique. Periodical oscillation peaks have been observed in the measurement. Further investigation shows that the additional modulation signal is coupled from the thermal electric cooler (TEC) controller to the laser driver, and then applied to the laser diode. The additional modulation can be eliminated by properly isolating the laser driving source from the TEC controller.

Thermal stress analysis for GaInAsP multiple quantum well wafer chemically bonded to Si (100)

An n-InP-based InGaAsP multiple-quantum-well wafer was bonded with p-Si by chemical surface activated bonding at 70 degrees C, and then annealed at 450 degrees C. Different thermal expansion coefficients between InP and Si will induce thermal stresses in the bonded wafer. Planar and cross-sectional distributions of thermal stress in the bonded InP-Si pairs were analyzed by a two-dimensional finite element method. In addition, the normal, peeling, and shear stresses were calculated by an analytic method. Furthermore, x-ray double crystalline diffraction was applied to measure the thermal strain and the strain caused by the mismatching of the crystalline orientation between InP (100) and Si (100). The wavelength redshift of the photoluminescence (PL) spectrum due to thermal strain was investigated via the calculation of the band structure, which is in agreement with the measured PL spectra.