复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

石油二次运移机理物理模拟研究

Oil and gas migration is very important for theoretical hydrocarbon geology study and exploration practice, but related research is weak. Physical simulation is a main method to study oil migration. Systematic experiments were done to quantitatively describe the migration patterns, path characters and oil saturation by adjusting the possible dynamic factors respectively. The following conclusions were drawn. 1. Darcy velocity and pore throat diameter were calculated according to seepage cross-sectional area and glass beads arrangement. With such normalized Darcy velocity and pore throat diameter, the date from one and two dimensional experiments can be reasonably drawn in two phase diagrams. It is found that the migration pattern can be identified using only one dimensionless number L which is defined as the ration of capillary number and Bond number. 2. Oil saturated in the pores between glass beads was used as calibration and oil saturation in the path was measured by magnetic resonance imaging. The results show that oil saturation in the center of migration path can reach 100%, is higher than oil saturation in the edge of migration path. 3. Percolation backbone during secondary oil migration was identified experimentally using Hele-Shaw cell. The backbone formed mainly because of the spatial variation of the cluster conductivity caused by oil saturation heterogeneity, main resistant force change, and path shrinkage and snap-off. Percolation backbone improves hydrocarbon migration efficiency and is a favorable factor for reservoir forming. 4. In the three dimensional filling models, the thickness of the secondary migration path is mall. It is only 2.5cm even for the piston pattern. Inclination of the model is the main influencing factor of the secondary path width.

High-power laser beam shaping by inseparable two-dimensional binary-phase gratings for surface modification of stamping dies

For surface modification of stamping dies, an inseparable two-dimensional binary-phase gratings is introduced to implement the wavefront transformation of high-power laser beams. The design and fabrication of the gratings are described in detail. Two-dimensional even sampling encoding scheme is adopted to overcome the limitations of conventional Dammann grating in the design of two-dimensional output patterns. High diffractive efficiency (>70%) can be achieved through the transformation of the Gaussian laser beam into several kinds of two-dimensional arrays in focal plan. The application of the binary-phase gratings in the laser surface modification of ductile iron is investigated, and the results show that the hardness and the wear resistance of the sample surface were improved significantly by using the binary-phase gratings. (C) 2008 Elsevier Ltd. All rights reserved.

Displacive To Order-Disorder Two-Step Phase Transition Model For Para-Ferroelectric Transition

Two-step phase transition model, displacive to order-disorder, is proposed. The driving forces for these two transitions are fundamentally different. The displacive phase transition is one type of the structural phase transitions. We clearly define the structural phase transition as the symmetry broking of the unit cell and the electric dipole starts to form in the unit cell. Then the dipole-dipole interaction takes place as soon as the dipoles in unit cells are formed. We believe that the dipole-dipole interaction may cause an order-disorder phase transition following the displacive phase transition. Both structural and order-disorder phase transition can be first-order or second-order or in between. We found that the structural transition temperatures can be lower or equal or higher than the order-disorder transition temperature. The para-ferroelectric phase transition is the combination of the displacive and order-disorder phase transitions. It generates a variety of transition configurations along with confusions. In this paper, we discuss all these configurations using our displacive to order-disorder two-step phase transition model and clarified all the confusions.

Enhancement and broadening of extreme-ultraviolet supercontinuum in a relative phase controlled two-color laser field

We experimentally demonstrate the generation of an extreme-ultraviolet (XUV) supercontinuum in argon with a two-color laser field consisting of an intense 7 fs pulse at 800 nm and a relatively weak 37 fs pulse at 400 nm. By controlling the relative time delay between the two laser pulses, we observe enhanced high-order harmonic generation as well as spectral broadening of the supercontinuum. A method to produce isolated attosecond pulses with variable width and intensity is proposed. (C) 2008 Optical Society of America.

Absolute phase control of spectra effects in a two-level medium driven by two-color ultrashort laser pulses

Using a omega-3 omega combination scenario, we investigate the absolute phase control of the spectra effects for ultrashort laser pulses propagating in a two-level medium. It is found that the higher spectral components can be controlled by the absolute phases. In particular, different absolute phase combinations can lead to the buildup or split of the even harmonics. (c) 2006 Elsevier B.V. All rights reserved.

Phase-matching control of high-order harmonic generation in a two-color laser field

The phase-matching condition of high-order harmonic generation driven by intense few-cycle pulses could be controlled by adding second-harmonic pulses to change the ionization fraction of the gaseous medium. The harmonic generation efficiency could be improved by moving the phase-matching point with an all-optical control of the ionization fraction or a proper change of the confocal parameter. A specific order of harmonics could be easily controlled to reach phase matching at a fixed higher gas pressure by adding second-harmonic pulses with a suitable intensity. Such an all-optical phase-matching control was demonstrated to be dependent upon the temporal delay between the fundamental-wave and second harmonic pulses.

A two-step phase-retrieval method in Fourier-transform ghost imaging

A two-step phase-retrieval method, based on Fourier-transform ghost imaging, was demonstrated. For the complex objects, the phase-retrieval process was divided into two steps: first got the complex object's amplitude from the Fourier-transform patterns of the squared object function, then combining with the Fourier-transform patterns of the object function to get the phase. The theoretical basis of this technique is outlined, and the experimental results are presented. (C) 2008 Elsevier B.V. All rights reserved.

Phase-step calibration technique based on a two-run-times-two-frame phase-shift method

A novel phase-step calibration technique is presented on the basis of a two-run-times-two-frame phase-shift method. First the symmetry factor M is defined to describe the distribution property of the distorted phase due to phase-shifter miscalibration; then the phase-step calibration technique, in which two sets of two interferograms with a straight fringe pattern are recorded and the phase step is obtained by calculating M of the wrapped phase map, is developed. With this technique, a good mirror is required, but no uniform illumination is needed and no complex mathematical operation is involved. This technique can be carried out in situ and is applicable to any phase shifter, whether linear or nonlinear. (c) 2006 Optical Society of America.

Discussions about FFT-based two-step phase-shifting algorithm

An FFT-based two-step phase-shifting (TPS) algorithm is described in detail and implemented by use of experimental interferograms. This algorithm has been proposed to solve the TPS problem with random phase shift except pi. By comparison with the visibility-function-based TPS algorithm, it proves that the FFT-based algorithm has obvious advantages in phase extracting. Meanwhile, we present a pi-phase-shift supplement to the TPS algorithm, which combines the two interferograms and demodulates the phase map by locating the extrema of the combined fringes after removing the respective backgrounds. So combining this method and FFT-based one, one could really implement the TPS with random phase shift. Whereafter, we systematically compare the TPS with single-interferogram analysis algorithm and conventional three-step phase-shifting one. The results demonstrate that the FFT-based TPS algorithm has a satisfactory accuracy. At last, based on the polarizing interferometry, a schematic setup of two-channel TPS interferometer with random phase shift is suggested to implement the simultaneous collection of interferograms. (c) 2007 Elsevier GrnbH. All rights reserved.

Wideband two-port beam splitter of a binary fused-silica phase grating

The usual beam splitter of multilayer-coated film with a wideband spectrum is not easy to achieve. We describe the realization of a wideband transmission two-port beam splitter based on a binary fused-silica phase grating. To achieve high efficiency and equality in the diffracted 0th and -1st orders, the grating profile parameters are optimized using rigorous coupled-wave analysis at a wavelength of 1550 nm. Holographic recording and the inductively coupled plasma dry etching technique are used to fabricate the fused-silica beam splitter grating. The measured efficiency of (45% x 2) = 90% diffracted into the both orders can be obtained with the fabricated grating under Littrow mounting. The physical mechanism of such a wideband two-port beam splitter grating can be well explained by the modal method based on two-beam interference of the modes excited by the incident wave. With the high damage threshold, low coefficient of thermal expansion, and wideband high efficiency, the presented beam splitter etched in fused silica should be a useful optical element for a variety of practical applications. (C) 2008 Optical Society of America.

Phase-locking of two large-core fibre lasers with 60W of coherent output power

A phase-locking fibre laser array with up to 60 W of coherent output power based on two large-core fibre is reported. The slope efficiency of the in-phase mode is 37%. For two cases of spacings between the cores, steady high-contrast interference stripes are observed. When the whole system operates under a high pump power level, no thermal effects for the spatial filter have been observed, which means that we can increase the coherent output power further by increasing the individual fibre laser power.

Experimental demonstration of phase locking of a two-dimensional fiber laser array using a self-imaging resonator

Phase locking of a two-dimensional fiber laser array is experimentally demonstrated by using a self-imaging resonator and a spatial filter. The stable beam profiles of in-phase mode and out-of-phase mode are observed by controlling the position of spatial filter. The phase locking fiber array with in-phase mode has produced 26 W coherent output. An antisymmetric eigenmode is also observed in our experiments. The phase locking is not sensitive to power variations among the pump beams and the configuration has the ability to repair a missing element. (C) 2008 American Institute of Physics.

Classical two-site fidelity and quantum phase transitions in the Ising model and the extended Hubbard model

In this Letter, the classical two-site-ground-state fidelity (CTGF) is exploited to identify quantum phase transitions (QPTs) for the transverse field Ising model (TFIM) and the one-dimensional extended Hubbard model (EHM). Our results show that the CTGF exhibits an abrupt change around the regions of criticality and can be used to identify QPTs in spin and fermionic systems. The method is especially convenient when it is connected with the density-matrix renormalization group (DMRG) algorithm. (C) 2008 Elsevier B.V. All rights reserved.

Disentanglement of two qubits coupled to an XY spin chain: Role of quantum phase transition

We study the disentanglement evolution of two spin qubits which interact with a general XY spin-chain environment. The dynamical process of the disentanglement is numerically and analytically investigated in the vicinity of a quantum phase transition (QPT) of the spin chain in both weak and strong coupling cases. We find that the disentanglement of the two spin qubits may be greatly enhanced by the quantum critical behavior of the environmental spin chain. We give a detailed analysis to facilitate the understanding of the QPT-enhanced decaying behavior of the coherence factor. Furthermore, the scaling behavior in the disentanglement dynamics is also revealed and analyzed.