各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

高斯束射线追踪

Using the approximate high-frequency asymptotic methods to solve the scalar wave equation, we can get the eikonal equation and transport equation. Solving the eikonal equation by the method of characteristics provides a mathematical derivation of ray tracing equations. So, the ray tracing system is folly based on the approximate high-frequency asymptotic methods. If the eikonal is complex, more strictly, the eikonal is real value at the rays and complex outside rays, we can derive the Gaussian beam. This article mainly concentrates on the theory of Gaussian beam. To classical ray tracing theory, the Gaussina beam method (GBM) has many advantages. First, rays are no longer required to stop at the exact position of the receivers; thus time-consuming two-point ray tracing can be avoided. Second, the GBM yields stable results in regions of the wavefield where the standard ray theory fails (e.g., caustics, shadows zones and critical distance). Third, unlike seismograms computed by conventional ray tracing techniques, the GBM synthetic data are less influenced by minor details in the model representation. Here, I realize kinematical and dynamical system, and based on this, realize the GBM. Also, I give some mathematical examples. From these examples, we can find the importance and feasibility of the ray tracing system. Besides, I've studied about the reflection coefficient of inhomogeneous S-electromagnetic wave at the interface of conductive media. Basing on the difference of directions of phase shift constant and attenuation constant when the electromagnetic wave propagates in conductive medium, and using the boundary conditions of electromagnetic wave at the interface of conductive media, we derive the reflection coefficient of inhomogeneous S-electromagnetic wave, and draw the curves of it. The curves show that the quasi total reflection will occur when the electromagnetic wave incident from the medium with greater conductivity to the medium with smaller conductivity. There are two peak, values at the points of the critical angles of phase shift constant and attenuation constant, and the reflection coefficient is smaller than 1. This conclusion is different from that of total reflection light obviously.

Analysis of electromagnetic propagation in birefringent thin film

In general, the propagating behavior of extraordinary wave in anisotropic materials is different from that in isotropic materials. With the tangential continuity of Maxwell's equations, the electromagnetic propagating behaviors have been investigated at the incident and exit interfaces of the uniaxial anisotropic thin film. The emphasis was placed on two interesting optical phenomena such as homolateral refraction behavior and wide-angle Brewster's phenomenon, which occurred at the interfaces of uniaxial anisotropic thin film.

Three-Dimensional Propagation of the Transmitted Shock Wave in a Square Cross-Sectional Chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.

Distinct Element Analysis on Propagation Characteristics of P-Wave in Rock Pillar with Finite length

以节理岩体等效刚度的概念为基础,讨论了离散元刚性块体模型中节理刚度的选取问题。采用面-面接触模型模拟了纵波在一维岩体中的传播,给出了纵波波形；研究了阻尼比、软弱夹层以及节理间是否可拉对波传播规律的影响。

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle Hows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained.

Anisotropic-source-induced changes of the degree of polarization of stochastic electromagnetic beams on propagation

It is shown that stochastic electromagnetic beams may have different degrees of polarization on propagation, even though they have the same coherence properties in the source plane. This fact is due to a possible difference in the anisotropy of the field in the source plane. The result is illustrated by some examples.

Effects of coherence on anisotropic electromagnetic Gaussian-Schell model beams on propagation

An analytical formula for the cross-spectral density matrix of the electric field of anisotropic electromagnetic Gaussian-Schell model beams propagating in free space is derived by using a tensor method. The effects of coherence on those beams are studied. It is shown that two anisotropic stochastic electromagnetic beams that propagate from the source plane z = 0 into the half-space z > 0 may have different beam shapes (i.e., spectral density) and states of polarization in the half-space, even though they have the same beam shape and states of polarization in the source plane. This fact is due to a difference in the coherence properties of the field in the source plane. (C) 2007 Optical Society of America.

Dressed multi-wave mixing in a V-type four-level atomic system

The dressed four- and six-wave mixings in a V-type four-level system are considered. Under two different dressed conditions, two- and three-photon resonant Autler-Townes splittings, accompanied by enhancement and suppression of wave mixing signal, are obtained analytically. Meanwhile, an electromagnetic induced transparency of multi-wave mixing is presented, which shows multiple peaks and asymmetric effects caused by one-photon, two-photon and three-photon resonances, separately. The slow light propagation multiple region of multi-wave mixing signal is also obtained.

Numerical simulation of an adaptive optics system with laser propagation in the atmosphere

A comprehensive model of laser propagation in the atmosphere with a complete adaptive optics (AO) system for phase compensation is presented, and a corresponding computer program is compiled. A direct wave-front gradient control method is used to reconstruct the wave-front phase. With the long-exposure Strehl ratio as the evaluation parameter, a numerical simulation of an AO system in a stationary state with the atmospheric propagation of a laser beam was conducted. It was found that for certain conditions the phase screen that describes turbulence in the atmosphere might not be isotropic. Numerical experiments show that the computational results in imaging of lenses by means of the fast Fourier transform (FFT) method agree well with those computed by means of an integration method. However, the computer time required for the FFT method is 1 order of magnitude less than that of the integration method. Phase tailoring of the calculated phase is presented as a means to solve the problem that variance of the calculated residual phase does not correspond to the correction effectiveness of an AO system. It is found for the first time to our knowledge that for a constant delay time of an AO system, when the lateral wind speed exceeds a threshold, the compensation effectiveness of an AO system is better than that of complete phase conjugation. This finding indicates that the better compensation capability of an AO system does not mean better correction effectiveness. (C) 2000 Optical Society of America.

Shock wave reflection from a wedge in a dusty gas

The present paper contains a detailed study of shock wave reflection from a wedge placed in various suspensions. In past works, the incident shock propagated initially in pure gas and the suspension started only at the leading edge of the deflecting wedge. However, in the present case the entire flow field is filled with a gas-dust suspension and the initial shock wave has steady-state structure relative to the shock front. In former studies the transmitted shock wave starts its propagation into the suspension and is reflected from the wedge at the same time. It is therefore obvious that the two unrelated processes of (2D) reflection and (1D) "transitional" relaxation occur simultaneously. In the present case the suspension behind the incident shock wave has reached steady state (i.e., it is a traveling wave) before the shock reaches the wedge leading edge. The reflection process from the deflecting wedge is studied for different dust mass loadings and different dust-particle diameter. It is shown that when the dust loading is low and the dust particle diameter is small the wave reflection pattern is similar to that observed in a similar pure gas case. In addition, an equilibrium state is reached, behind the evolved waves, very quickly. On the other hand, when the dust loading is relatively high and/or the dust particle diameter is relatively large, the observed reflection wave pattern is very different from that seen in a similar pure gas case. In such cases it takes much longer time to reach an equilibrium state behind the reflecting waves. It is also shown that the dust presence significantly affects the (gas) pressure on the wedge surface. The higher the dust loading is, the higher the pressure on the wedge surface. Suspensions composed of solid particle of different size, but having the same dust mass loading, will approach the same equilibrium pressure. However, it will take longer time to reach an equilibrium state for suspensions having large diameter particles. (C) 2004 Elsevier Ltd. All rights reserved.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

Experimental investigation on low-temperature fatigue crack propagation in offshore structural steel A131 under random sea ice loading

An apparatus of low-temperature controlling for fatigue experiments and its crack measuring system were developed and used for offshore structural steel A131 under conditions of both low temperature and random sea ice. The experimental procedures and data processing were described, and a universal random data processing software for FCP under spectrum loading was written. Many specific features of random ice-induced FCP which differed with constant amplitude FCP behaviours were proposed and temperature effect on ice-induced FCP was pointed out with an easily neglected aspect in designing for platforms in sea ice emphasized. In the end, differences of FCP behaviours between sea ice and ocean wave were presented.