粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

Thermocapillary Motion Of Droplets At Large Marangoni Numbers

In this paper, the thermocapillary motion problem of drops is investigated using the axisymmetric model. The front-tracking method is employed to capture the drop interface. We find that the migration velocity of the drop is greatly influenced by the temperature field in the drop when Ma is fairly large (>100), which leads to an increase-decrease migration velocity at the beginning of our simulations. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

神光Ⅱ装置中空间滤波器系统的稳定性设计与分析

考虑到空间滤波器系统在激光装置中的稳定性要求，分析了影响稳定性的各种因素。对现行装置中的SF4空间滤波器系统进行振动模态测试，初步确定系统的工作状况，并结合振动理论明确系统稳定性的设计参数。利用有限元软件ANSYS对空间滤波器系统及相关的设计进行仿真验证计算。结果表明：SF4系统的一阶固有频率达到83Hz，与实验结果相符合，满足系统稳定性的要求；在系统中引入支撑不但能有效增加系统的刚度，同时能转移系统的振型；系统基础支撑平台的选择要求平台独立，并且选择较低的高度以便于控制上层光学元件的变形；在系统中的引入

数值模拟探针诱导表面等离子体共振耦合纳米光刻

采用有损耗介质和色散介质的二维时域有限差分方法,数值模拟了以光波长514.5 nm的p偏振基模高斯光束为入射光源,激发Kretschmann型表面等离子体共振,并通过探针的局域场增强效应实现纳米光刻的新方法——探针诱导表面等离子体共振耦合纳米光刻.分别就探针与记录层的间距以及探针针尖大小,模拟分析了不同情况下探针的局域场增强效应和记录层表面的相对电场强度振幅分布.结果表明,探针工作在接触模式时,探针的局域场增强效应最明显,记录层表面的相对电场强度振幅的对比度最大;当探针针尖距记录层5 nm时,针尖下方记录层表面的相对电场强度振幅大于光刻临界值的分布宽度与针尖尺寸相近.

Impact of symmetrized and Burt-Foreman Hamiltonians on spurious solutions and energy levels of InAs/GaAs quantum dots

We present a systematic investigation of calculating quantum dots (QDs) energy levels using finite element method in the frame of eight-band k . p method. Numerical results including piezoelectricity, electron and hole levels, as yell as wave functions are achieved. In the calculation of energy levels, we do observe spurious solutions (SSs) no matter Burt-Foreman or symmetrized Hamiltonians are used. Different theories are used to analyse the SSs, we find that the ellipticity theory can give a better explanation for the origin of SSs and symmetrized Hamiltonian is easier to lead to SSs. The energy levels simulated with the two Hamiltonians are compared to each other after eliminating SSs, different Hamiltonians cause a larger difference on electron energy levels than that on hole energy levels and this difference decreases with the increase of QD size.

Effects of multistage coupling and disk flexibility on mistuned bladed disk dynamics

The effects ofdisk flexibility and multistage coupling on the dynamics of bladed disks with and without blade mistuning are investigated. Both free and forced responses are examined using finite element representations of example single and two-stage rotor models. The reported work demonstrates the importance of proper treatment of interstage (stage-to-stage) boundaries in order to yield adequate capture of disk-blade modal interaction in eigenfrequency veering regions. The modified disk-blade modal interactions resulting from interstage-coupling-induced changes in disk flexibility are found to have a significant impact on (a) tuned responses due to excitations passing through eigenfrequency veering regions, and (b) a design's sensitivity to blade mistuning. Hence, the findings in this paper suggest that multistage analyses may be required when excitations are expected to fall in or near eigenfrequency veering regions or when the sensitivity to blade mistuning is to be accounted for Conversely, the observed sensitivity to disk flexibility also indicates that the severity of unfavorable structural interblade coupling may be reduced significantly by redesigning the disk(s) and stage-to-stage connectivity. The relatively drastic effects of such modifications illustrated in this work indicate that the design modifications required to alleviate veering-related response problems may be less comprehensive than what might have been expected.

Trapping of low-refractive-index particles with azimuthally polarized beam

Azimuthally polarized beams, focused by a high-numerical-aperture (NA) objective lens, form a hollow intensity distribution near the focus, which is appropriate for trapping low-refractive-index particles, in contrast to common linearly polarized or radially polarized beams. In this paper, the field distribution of the azimuthally polarized beam focused by a high-NA objective is described by the vectorial diffraction integral, and then the radiation forces on spherical particles with different parameters such as radius and refractive index are calculated by the T-matrix method. Numerical results show that the azimuthally polarized beam not only can steadily trap low-refractive-index particles at the focus center but also can trap multiple high-refractive-index particles around the focus center by virtue of the hollow-ring configuration. The range of the sizes of low-refractive-index particles that can be trapped steadily are presented, corresponding to different parameters such as the NA of the objective and the relative refractive index, based on which the NA of the objective can be selected to trap the appropriate size of particles. (C) 2009 Optical Society of America

Exact description of a cylindrically symmetrical complex-argument Laguerre-Gauss beam

Based on the perturbative series representation of a complex-source-point spherical wave an expression for cylindrically symmetrical complex-argument Laguerre-Gauss beams of radial order n is derived. This description acquires the accuracy up to any order of diffraction angle, and its first three corrected terms are in accordance with those given by Seshadri [Opt. Lett. 27, 1872 (2002)] based on the virtual source method. Numerical results show that on the beam axis the number of orders of nonvanishing nonparaxial corrections is equal to n. Meanwhile a higher radial mode number n leads to a smaller convergent domain of radius. (C) 2008 Optical Society of America.

Design and optimization of highly nonlinear low-dispersion crystal fiber with high birefringence for four-wave mixing

We have proposed a novel type of photonic crystal fiber (PCF) with low dispersion and high nonlinearity for four-wave mixing. This type of fiber is composed of a solid silica core and a cladding with a squeezed hexagonal lattice elliptical airhole along the fiber length. Its dispersion and nonlinearity coefficient are investigated simultaneously by using the full vectorial finite element method. Numerical results show that the proposed highly nonlinear low-dispersion fiber has a total dispersion as low as +/- 2.5 ps nm(-1) km(-1) over an ultrabroad wavelength range from 1.43 to 1.8 mu m, and the corresponding nonlinearity coefficient and birefringence are about 150 W-1 km(-1) and 2.5 x 10(-3) at 1.55 mu m, respectively. The proposed PCF with low ultraflattened dispersion, high nonlinearity, and high birefringence can have important application in four-wave mixing. (C) 2010 Optical Society of America

Progress in luminescent films prepared by sol-gel process

In this paper, state of the art for luminescent films prepared by sol-gel process has been reviewed. The basic process and characteristics for the synthesis of luminescent films via sol-gel method, Characterization methods for the films and the current status for the development and application of luminescent films are discussed in the context, An elucidation has been made on the luminescent films classified by composition, including inorganic luminescent films and organic/inorganic hybrid luminescent films. The sol-gel derived luminescent films have found applications in the display devices for photoluminescence, electroluminescence, cathodoluminescence and field emission etc. Finally, the future development tendency for the luminescent films are forecasted.