地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

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.

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic–plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a “hold” period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and “hold-at-the-maximum-indenter-displacement” for determining the instantaneous modulus using spherical indenters.

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic-plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a "hold" period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and "hold-at-the-maximum-indenter-displacement" for determining the instantaneous modulus using spherical indenters.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Fast Discrimination of Commerical Corn Varieties Based on Near Infrared Spectra

The existing methods for the discrimination of varieties of commodity corn seed are unable to process batch data and speed up identification, and very time consuming and costly. The present paper developed a new approach to the fast discrimination of varieties of commodity corn by means of near infrared spectral data. Firstly, the experiment obtained spectral data of 37 varieties of commodity corn seed with the Fourier transform near infrared spectrometer in the wavenurnber range from 4 000 to 12 000 cm (1). Secondly, the original data were pretreated using statistics method of normalization in order to eliminate noise and improve the efficiency of models. Thirdly, a new way based on sample standard deviation was used to select the characteristic spectral regions, and it can search very different wavenumbers among all wavenumbers and reduce the amount of data in part. Fourthly, principal component analysis (PCA) was used to compress spectral data into several variables, and the cumulate reliabilities of the first ten components were more than 99.98%. Finally, according to the first ten components, recognition models were established based on BPR. For every 25 samples in each variety, 15 samples were randomly selected as the training set. The remaining 10 samples of the same variety were used as the first testing set, and all the 900 samples of the other varieties were used as the second testing set. Calculation results showed that the average correctness recognition rate of the 37 varieties of corn seed was 94.3%. Testing results indicate that the discrimination method had higher precision than the discrimination of various kinds of commodity corn seed. In short, it is feasible to discriminate various varieties of commodity corn seed based on near infrared spectroscopy and BPR.

Splitting the Fast and Slow Motions in Molecular Dynamics Simulations Based on the Change of Cold Potential Well Bottom

The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Fast Haar Transform Based Feature Extraction for Face Representation and Recognition

Subspace learning is the process of finding a proper feature subspace and then projecting high-dimensional data onto the learned low-dimensional subspace. The projection operation requires many floating-point multiplications and additions, which makes the projection process computationally expensive. To tackle this problem, this paper proposes two simple-but-effective fast subspace learning and image projection methods, fast Haar transform (FHT) based principal component analysis and FHT based spectral regression discriminant analysis. The advantages of these two methods result from employing both the FHT for subspace learning and the integral vector for feature extraction. Experimental results on three face databases demonstrated their effectiveness and efficiency.

Ultrasonic relaxation kinetics on fast deuteron transfer reaction

Propylamine has been selected to investigate the isotope effect of a fast deuteron transfer reaction by ultrasonic relaxation method. Ultrasonic absorption coefficients of propylamine in heavy water (D2O) at 25 degrees C in the concentration range from 0.0107 to 0.6300 mol dm(-3) have been measured by pulse and resonance methods over the frequency range from 0.8 to 220 MHz. A Debye-type single relaxation absorption has been observed in the solution. From the dependence of the ultrasonic relaxation parameters on the concentration and solution pH, the source of the observed relaxation has been attributed to a perturbation of the chemical equilibrium associated with the deuteron transfer reaction. The rate and equilibrium constants have been determined by the measurement of the deuteroxyl ion concentration dependence of the relaxation frequency. Also the standard volume change of the reaction has been determined from the concentration dependence of the maximum absorption per wavelength and the adiabatic compressibility has been calculated from the density and the sound velocity in the solution. These results have then been compared with those obtained for propylamine in light water (H2O). The forward rate constant is greater and the reverse rate constant is smaller in DO than in H2O. The standard volume change for deuteron transfer is greater than that for proton transfer reaction, and the adiabatic compressibility shows a similar trend. These data support an argument that there exists a stronger hydrogen bond in D2O than in H2O. The difference of the stability in the intermediate states, R-ND3+... OD- and R-NH3+... OH-, has also been considered from the results of the isotope effects.

DETERMINATION OF PARAMETERS OF FAST CHEMICAL REACTION OF ALIZARIN RED S BY SPECTROELECTROCHEMICAL METHOD WITH LOPTLC

The rate constant of very fast chemical reaction generally can be measured by electrochemical methods, but can not by the thin layer electrochemical methods because of the influence of diffusion effect. Long optical path length thin layer cell (LOPTLC) with large ratio of electrode area to solution volume can be used to monitor the fist chemical reaction in situ with high sensitivity and accuracy. It enable the adsorption spectra to be measured without the influence of diffusion effect. In the present paper, a fast chemical reaction of Alizarin Red S (ARS) with its oxidative state has been studied. The reaction equilibrium constant (K) under different potentials can be determined by single step potential-absorption spectra in LOPTLC. An equilibrium constant of 7.94 x 10(5) l.mol(-1) for the chemical reaction has been obtained from the plot of lgK vs. (E - E-1(0)'). Rate constant (k) under different potentials can be measured by single step potential-chronoabsorptiometry. A rate constant of 426.6 l.mol(-1).s(-1) for the chemical reaction has been obtained from the plot of lgK vs. (E - E-1(0)') with (E - E-1(0)') = 0.