基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Difference in Enzyme Activity and Conformation of Glucose Oxidase Before and After Purification

EEnzyme activity of commercial glucose oxidase was enhanced after purification through a strong anionic exchange resin. In order to get a better insight into this phenomenon, surface pressure–area ( –A) isotherms and surface pressure–time ( –t) isotherms was used to study the interaction and the absorption at different pH values of the subphases between octadecylamine and glucose oxidase purified by a styrene system quaternary ammonium type strongly basic anionic exchange resin. Circular dichroism (CD), electrophoresis and enzyme activity measurements were conducted to study these phenomena. A preliminary hypothesis has been suggested to explain why the enzyme activity of purified glucose oxidase was higher than that of the commercial one. © 2002 Elsevier Science B.V. All rights reserved.

What Do We Know about Long Laminar Plasma Jets?

Silent and stable long laminar plasma jets can be generated in a rather wide range of working parameters. The laminar flow state can be maintained even if considerable parameter fluctuations exist in the laminar plasma jet or if there is an impact of laterally injected particulate matter and its carrier gas. The attractive special features of laminar plasma jets include extremely low noise level, less entrainment of ambient air, much longer and adjustable high-temperature region length, and smaller axial gradient of plasma parameters. Modeling results show that the laminar plasma jet length increases with increasing jet inlet velocity or temperature and the effect of natural convection on laminar plasma jet characteristics can be ignored, consistent with experimental observations. The large difference between laminar and turbulent plasma jet characteristics is revealed to be due to their different laws of surrounding gas entrainment. Besides the promising applications of the laminar plasma jet to remelting and cladding strengthening of the metallic surface and to thermal barrier coating preparation, it is expected that the laminar plasma jet can become a rather ideal object for the basic studies of thermal plasma science owing to the nonexistence of the complexity caused by turbulence.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

What Do We Know about Long Laminar Plasma Jets?

Silent and stable long laminar plasma jets can be generated in a rather wide range of working parameters. The laminar flow state can be maintained even if considerable parameter fluctuations exist in the laminar plasma jet or if there is an impact of laterally injected particulate matter and its carrier gas. The attractive special features of laminar plasma jets include extremely low noise level, less entrainment of ambient air, much longer and adjustable high-temperature region length, and smaller axial gradient of plasma parameters. Modeling results show that the laminar plasma jet length increases with increasing jet inlet velocity or temperature and the effect of natural convection on laminar plasma jet characteristics can be ignored, consistent with experimental observations. The large difference between laminar and turbulent plasma jet characteristics is revealed to be due to their different laws of surrounding gas entrainment. Besides the promising applications of the laminar plasma jet to remelting and cladding strengthening of the metallic surface and to thermal barrier coating preparation, it is expected that the laminar plasma jet can become a rather ideal object for the basic studies of thermal plasma science owing to the nonexistence of the complexity caused by turbulence.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Geographic variation of Chevrier's field mouse (Apodemus chevrieri) (Milne-Edwards, 1868) (Muridae : Murinae) from southwestern China based on cranial morphometric variables

Geographic variation of Chevrier's field mouse (Apodemus chevrieri) (Milne-Edwards, 1868) (Muridae: Murinae) from southwestern China based on cranial morphometric variables. Zoological Studies 47(4): 393-401. A sample of 134 specimens of Apodemus chevrieri was investigated in the present study. Individuals were divided into male and female groups, and these were respectively subjected to multivariate analysis. Results indicated that 3 geographic populations of A. chevrieri inhabit southwestern China: a Sichuan population in western Sichuan Province; a northwestern Yunnan population ranging from northwestern Yunnan Province eastward to southern Sichuan Province; and a central Yunnan population in central Yunnan Province. In addition, a coefficient of difference analysis was performed among these 3 geographic populations. The results suggested that these 3 geographical populations of A. chevrieri belonged to 2 subspecies. Furthermore, we discuss the relationships of the subspecific differentiation of A. chevneri with changes in latitude in southwestern China.

High efficiency operation of butt joint line-defect-waveguide microlaser in two-dimensional photonic crystal slab

Butt joint line-defect-waveguide microlasers are demonstrated on photonic crystal slabs with airholes in a triangular lattice. Such microlaser is designed to increase the output power from the waveguide edge directly. The output power is remarkably enhanced to 214 times higher by introducing chirped structure in the output waveguide. The lasing mode operates in the linear dispersion region of the output waveguide so that the absorption loss due to the band-edge effect is reduced. The laser resonance is illustrated theoretically using the finite difference time domain method. A practical high power efficiency of 20% is obtained in this microlaser. (C) 2008 American Institute of Physics.

Indium-Induced Effect on Polarized Electroluminescence from InGaN/GaN MQWs Light Emitting Diodes

Polarization-resolved edge-emitting electroluminescence (EL) studies of InGaN/GaN MQWs of wavelengths from near-UV (390 nm) to blue (468 nm) light-emitting diodes (LEDs) are performed. Although the TE mode is dominant in all the samples of InGaN/GaN MQW LEDs, an obvious difference of light polarization properties is found in the InGaN/GaN MQW LEDs with different wavelengths. The polarization degree decreases from 52.4% to 26.9% when light wavelength increases. Analyses of band structures of InGaN/GaN quantum wells and luminescence properties of quantum dots imply that quantum-dot-like behavior is the dominant reason for the low luminescence polarization degree of blue LEDs, and the high luminescence polarization degree of UV LEDs mainly comes from QW confinement and the strain effect. Therefore, indium induced carrier confinement (quantum-dot-like behavior) might play a major role in the polarization degree change of InGaN/GaN MQW LEDs from near violet to blue.

Integration of a photonic crystal polarization beam splitter and waveguide bend

In this work, we present the design of an integrated photonic-crystal polarization beam splitter (PC-PBS) and a low-loss photonic-crystal 60 waveguide bend. Firstly, the modal properties of the PC-PBS and the mechanism of the low-loss waveguide bend are investigated by the two-dimensional finite-difference time-domain (FDTD) method, and then the integration of the two devices is studied. It shows that, although the individual devices perform well separately, the performance of the integrated circuit is poor due to the multi-mode property of the PC-PBS. By introducing deformed airhole structures, a single-mode PC-PBS is proposed, which significantly enhance the performance of the circuit with the extinction ratios remaining above 20dB for both transverse-electric (TE) and transverse-magnetic (TM) polarizations. Both the specific result and the general idea of integration design are promising in the photonic crystal integrated circuits in the future. (C) 2009 Optical Society of America

Reduced divergence angle of photonic crystal vertical-cavity surface-emitting laser

The reduced divergence angle of the photonic crystal vertical-cavity surface-emitting laser (PC-VCSEL) was investigated in both theory and experiment. The photonic crystal waveguide possessed the weakly guiding waveguide characteristic, which accounted for the reduction of the divergence angle. The three-dimensional finite-difference time-domain method was used to simulate the designed PC-VCSEL, and a calculated divergence angle of 5.2 degrees was obtained. The measured divergence angles of our fabricated PC-VCSEL were between 5.1 degrees and 5.5 degrees over the entire drive current range, consistent with the numerical results. This is the lowest divergence angle of the fabricated PC-VCSEL ever reported.

Measurement of polar C-plane and nonpolar A-plane InN/ZnO heterojunctions band offsets by x-ray photoelectron spectroscopy

The valence band offsets of the wurtzite polar C-plane and nonpolar A-plane InN/ZnO heterojunctions are directly determined by x-ray photoelectron spectroscopy to be 1.76 +/- 0.2 eV and 2.20 +/- 0.2 eV. The heterojunctions form in the type-I straddling configuration with a conduction band offsets of 0.84 +/- 0.2 eV and 0.40 +/- 0.2 eV. The difference of valence band offsets of them mainly attributes to the spontaneous polarization effect. Our results show important face dependence for InN/ZnO heterojunctions, and the valence band offset of A-plane heterojunction is more close to the "intrinsic" valence band offset.

Band-gap bowing and p-type doping of (Zn, Mg, Be)O wide-gap semiconductor alloys: a first-principles study

Using a first-principles band-structure method and a special quasirandom structure (SQS) approach, we systematically calculate the band gap bowing parameters and p-type doping properties of (Zn, Mg, Be)O related random ternary and quaternary alloys. We show that the bowing parameters for ZnBeO and MgBeO alloys are large and dependent on composition. This is due to the size difference and chemical mismatch between Be and Zn(Mg) atoms. We also demonstrate that adding a small amount of Be into MgO reduces the band gap indicating that the bowing parameter is larger than the band-gap difference. We select an ideal N atom with lower p atomic energy level as dopant to perform p-type doping of ZnBeO and ZnMgBeO alloys. For N doped in ZnBeO alloy, we show that the acceptor transition energies become shallower as the number of the nearest neighbor Be atoms increases. This is thought to be because of the reduction of p-d repulsion. The N-O acceptor transition energies are deep in the ZnMgBeO quaternary alloy lattice-matched to GaN substrate due to the lower valence band maximum. These decrease slightly as there are more nearest neighbor Mg atoms surrounding the N dopant. The important natural valence band alignment between ZnO, MgO, BeO, ZnBeO, and ZnMgBeO quaternary alloy is also investigated.

松辽徐家围子断陷层序地层充填模式及含油气性特征

Based on the principle and method of sequence, the author describes the sequence-filling model of the rifting basin of Xujiaweizi and its gas exploration potential. The object of this paper belongs to the area around Shengping-Wangjiatun anticline. Its srtatigraphy includes Huoshiling Formation (neutral and basic volcanic rocks), Shahezi Formation (coal bedding and mud and some sandstone) and Yingcheng Formation from bottom to top. These stratigraphy units are defined by author as mesosequences respectively. The author emphasizes that the main control factors of sequence change with the types of basin and stage of basin. So the sequence is researched according to the types of basin. This viewpoint is very new, and it is consistent with the principle of sequence. Volcanic action is very frequent and acute, topography difference is obvious. Between the volcanic events, Shahezi Formation is formed, which mainly consists-of sedimentary rocks. Based on the datum from seismic section and drilling core and well-logging, the author analyzes the single unit and unit set and system tract and sedimentary fancies, then, according to the accommodation space change and marking of sequence boundary, Shahezi Formation is divided into two Third-scale sequences. The sedimentary fancies and depth distribution are described. The author also pointed out that the volcanic rocks consume the accommodation space, so volcanic rocks can influence the development of sequence. Based on the concept of accommodation space, the author put volcanic rocks into sequence frame, which normally consists of sedimentary rocks. The topography of volcanic is controlled by lithology of volcanic rocks, the pattern of volcanic eruption and the topography before volcanic eruption. The topography of volcanic can influence sedimentation and the filling pattern of sedimentary rocks. The author describes the composition and lithology fancies and depth distribution of volcanic rocks. The volcanic rocks and Volcanic building, volcanic structure is recognized on seismic section. The author paid a special attention to the relationship between sedimentation and volcanism. Finally, the author analyses the combination of source-reservoir-cover unit in sequence frame. The mudstone of Shahezi Formation has a great depth, the Kerogene in it belongs to type II and III, which tends to produce gas. The Yingcheng Formation lies between Shahezi Formation and Denglouku Formation, belonging to good reservoir. The volcanic rocks of Huoshiling Formation often formed high building, which can capture the gas produced from Shahezi Formation. The stratigraphy of rifting basin of Xujiaweizi has the great potential of gas exploration. This paper claims the following creative points: 1. The author applied the principle and method of sequence to rifting basin, greatly extending its research area and topic issues. 2. The author pointed out that basin of different type and of different stage has a different type of sequence. This is caused by the different main control factors of sequence. 3. Put volcanic rocks into the sequence frame, discussing the probability of regarding the volcanic rocks as the component of sequence, dealing with the relationship between sedimentation and volcanism and its influence to the source-reservoir-cover system. 4. The author pointed out that the filling pattern of rifting basin are determined by the filling pattern of megasequence, whose filling pattern is determined by the filling pattern of system tract and the change of accommodation space.