Two-Point Closure Strategy in the Mapping Closure Approximation Approach

A two-point closure strategy in mapping closure approximation (MCA) approach is developed for the evolution of the probability density function (PDF) of a scalar advected by stochastic velocity fields. The MCA approach is based on multipoint statistics. We formulate a MCA modeled system using the one-point PDFs and two-point correlations. The MCA models can describe both the evolution of the PDF shape and the rate at which the PDF evolves.

改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

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

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.

Visualization of molecule interaction between antigen and antibody - One of ellipsometric imaging applications

It is to investigate molecule interactions between antigen and antibody with ellipsometric imaging technique and demonstrate some features and possibilities offered by applications of the technique. Molecule interaction is an important interest for molecule biologist and immunologist. They have used some established methods such as immufluorcence, radioimmunoassay and surface plasma resonance, etc, to study the molecule interaction. At the same time, experimentalists hope to use some updated technique with more direct visual results. Ellipsometric imaging is non-destructive and exhibits a high sensitivity to phase transitions with thin layers. It is capable of imaging local variations in the optical properties such as thickness due to the presence of different surface concentration of molecule or different deposited molecules. If a molecular mono-layer (such as antigen) with bio-activity were deposited on a surface to form a sensing surface and then incubated in a solution with other molecules (such as antibody), a variation of the layer thickness when the molecules on the sensing surface reacted with the others in the solution could be observed with ellipsometric imaging. Every point on the surface was measured at the same time with a high sensitivity to distinguish the variation between mono-layer and molecular complexes. Ellipsometric imaging is based on conventional ellipsometry with charge coupled device (CCD) as detector and images are caught with computer with image processing technique. It has advantages of high sensitivity to thickness variation (resolution in the order of angstrom), big field of view (in square centimeter), high sampling speed (a picture taken within one second), and high lateral resolution (in the order of micrometer). Here it has just shown one application in study of antigen-antibody interaction, and it is possible to observe molecule interaction process with an in-situ technique.

Iris recognition algorithm based on point covering of high-dimensional space and neural network

In this paper, we constructed a Iris recognition algorithm based on point covering of high-dimensional space and Multi-weighted neuron of point covering of high-dimensional space, and proposed a new method for iris recognition based on point covering theory of high-dimensional space. In this method, irises are trained as "cognition" one class by one class, and it doesn't influence the original recognition knowledge for samples of the new added class. The results of experiments show the rejection rate is 98.9%, the correct cognition rate and the error rate are 95.71% and 3.5% respectively. The experimental results demonstrate that the rejection rate of test samples excluded in the training samples class is very high. It proves the proposed method for iris recognition is effective.

Double proton transfer and one-electron oxidation behaviors in double H-bonded glycinamide-formamidine complex and comparison with biological base pair

The behaviors of double proton transfer (DPT) occurring in a representative glycinamide-formamidine complex have been investigated employing the B3LYP/6-311++G** level of theory. Computational results suggest that the participation of a formamidine molecule favors the proceeding of the proton transfer (PT) for glycinamide compared with that without mediator-assisted case. The DPT process proceeds with a concerted mechanism rather than a stepwise one since no zwitterionic complexes have been located during the DPT process. The barrier heights are 14.4 and 3.9 kcal/mol for the forward and reverse directions, respectively. However, both of them have been reduced by 3.1 and 2.9 kcal/mol to 11.3 and 1.0 kcal/mol with further inclusion of zero-point vibrational energy (ZPVE) corrections, where the lower reverse barrier height implies that the reverse reaction should proceed easily at any temperature of biological importance. Additionally, the one-electron oxidation process for the double H-bonded glycinamide-formamidine complex has also been investigated. The oxidated product is characterized by a distonic radical cation due to the fact that one-electron oxidation takes place on glycinamide fragment and a proton has been transferred from glycinamide to formamidine fragment spontaneously. As a result, the vertical and adiabatic ionization potentials for the neutral double H-bonded complex have been determined to be about 8.46 and 7.73 eV, respectively, where both of them have been reduced by about 0.79 and 0.87 eV relative to those of isolated glycinamide due to the formation of the intermolecular H-bond with formamidine. Finally, the differences between model system and adenine-thymine base pair have been discussed briefly.