基于广义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 urinary 8-hydroxydeoxyguanosine by capillary electrophoresis and solid-phase extraction

Urinary 8-hydroxydeoxyguanosine (80HdG) has been considered as an excellent marker of individuals at high risk of developing cancer. Until now, urinary 80HdG has largely been measured by high-performance liquid chromatography with electrochemical detection. A new method for the analysis of urinary 80HdG by high-performance capillary electrophoresis has been developed and optimized in our laboratory. A single step solid-phase extraction procedure was optimized and used for extracting 80HdG from human urine. Separations were performed in an uncoated silica capillary (50 cm x 50 tm i.d.) using a P/ACE MDQ system with UV detection. The separation of 80HdG from interfering urinary matrix components is optimized with regard to pH, applied voltage, pressure injection time and concentration of SDS in running buffer. The detection limit of this method is 0.4 mug/ml, the linear range is 0.8-500 mug/ml, the correlation coefficients levels is better than 0.999. The developed method is simple, fast and good reproducibility, furthermore, it requires a very small injection volumes and low costs of analysis, which makes it possible to provide a new noninvasive assay for an indirect measurement of oxidative DNA damage.

Fast formation of NaA zeolite membrane in the microwave field

NaA zeolite membrane was successfully synthesized on the porous alpha-Al2O3 support by microwave heating. The synthesis of NaA zeolite membrane in the microwave field only needs 15 min and the synthesis time is 10 times shorter than that by conventional heating. SEM characterization indicates that the zeolite crystals in the NaA zeolite membrane synthesized by microwave heating are uniform in size; the membrane thickness is about 4 mu m and is thinner than that of the NaA zeolite membrane synthesized by conventional heating. Gas permeation studies indicate that the permeances of the NaA zeolite membrane synthesized by microwave heating are 3-4 times higher than those of the NaA zeolite membrane synthesized by conventional heating, while their permselectivities are comparable.

Fast analysis of capillary electrochromatography with non-porous ODS as the stationary phase

High-speed capillary electrochromatography was developed on both short and long packed columns with 2 mu m non-porous ODS as the stationary phase. Factors that affect the analysis time of samples, such as voltage, electrolyte concentration, pH and organic modifier concentration in the mobile phase, were studied systematically. Fast analysis of aromatic compounds within 13 seconds was realized with column efficiency of 573,000 plates/m and a R.S.D.% of the retention times of all components in 8 consecutive injections below 1.0%. which demonstrated the high efficiency and high reproducibility of such a technique. In addition, DNPH derived aldehydes and ketones in both standards and environmental samples were separated with high speed.

Soil column chromatography for correlation between capacity factors and soil organic partition coefficients for eight pesticides

A soil column chromatographic method was developed to measure the capacity factors (k') of pesticides, in which soil acted as a stationary phase and methanol-water mixture as an eluent. The k' values of eight pesticides, including three insecticides (methiocarb, azinphos-methyl, fenthion), four fungicides (triadimenol, fuberidazole, tebuconazole, pencycuron), and one herbicide (atrazine), were found to be well fitted to a retention equation, ln k'=ln k(w)'-S-phi. Due to similar interactions of solutes with soil and solvent in both sorption determination and retention experiment, log k' has a good linear correlation with log K-oc for the eight pesticides from different classes, in contrast with poor correlation between log k' from C-18 column and log K-oc. So the method provides a tool for rapid estimation of K-oc from experimental k'. (C) 1999 Elsevier Science Ltd. All rights reserved.