Kinking of an interface crack between 2 dissimilar anisotropic elastic solids

A detailed analysis of kinking of an interface crack between two dissimilar anisotropic elastic solids is presented in this paper. The branched crack is considered as a distributed dislocation. A set of the singular integral equations for the distribution function of the dislocation density is developed. Explicit formulas of the stress intensity factors and the energy release rates for the branched crack are given for orthotropic bimaterials and misoriented orthotropic bicrystals. The role of the stress parallel to the interface, sigma0 is taken into account in these formulas. The interface crack can advance either by continued extension along the interface or by kinking out of the interface into one of the adjoining materials. This competition depends on the ratio of the energy release rates for interface cracking and for kinking out of the interface and the ratio of interface toughness to substrate toughness. Throughout the paper, the influences of the inplane stress sigma0 on the stress intensity factors and the energy release rates for the branched crack, which can significantly alter the conditions for interface cracking, are emphasized.

Pulse shaping properties of volume holographic gratings in anisotropic media

Based on a modified coupled wave theory, the pulse shaping properties of volume holographic gratings (VHGs) in anisotropic media VHGs are studied systematically. Taking photorefractive LiNbO3 crystals as an example, the combined effect that the grating parameters, the dispersion and optical anisotropy of the crystal, the pulse width, and the polarization state of the input ultrashort pulsed beam (UPB) have on the pulse shaping properties are considered when the input UPB with arbitrary polarization state propagates through the VHG. Under the combined effect, the diffraction bandwidth, pulse profiles of the diffracted and transmitted pulsed beams, and the total diffraction efficiency are shown. The studies indicate that the properties of the shaping of the o and e components of the input UPB in the crystal are greatly different; this difference can be used for pulse shaping applications. (c) 2006 Optical Society of America.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.

Design of anisotropic reflector with birefringent thin films

A novel design for dielectric anisotropic mirrors with birefringent thin films for normal incidence is presented. This mirror consists of a stack of quarter-wave biaxial layers. The biaxial anisotropic layers can be fabricated by oblique deposition. The reflectance is different for two linear polarizations of light incidence on the mirrors. As a numerical example, the design is carried out on glass with TiO2 and ZrO2. These thin films could be applied to anisotropic reflective devices for lasers.

Analysis of characteristic matrix for a uniaxially birefringent thin film

As the principal section consistent with the principal plane, electromagnetic propagation in a uniaxially birefingent thin film can be described with a concise 2 x 2 characteristic matrix, in which the refractive indices of the forward and backward propagating extraordinary rays are different and dependent. In this letter, based on Huygen's construction, the refractive indices and effective optical admittances of the forward and backward propagating extraordinary rays are discussed further, and the 2 x 2 characteristic matrix is simplified. Furthermore, the input optical admittance, reflectance and transmittance of assembly is presented, just as an isotropic thin film. The result can be extended to the general case of multilayer uniaxially birefringent thin films with their optic axes in the incident plane. Copyright (c) EPLA, 2007.

Propagation properties of beams generated by Gaussian mirror resonator in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of beams generated by Gaussian mirror resonator (GMR) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident beams generated by GMR propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of the incident beams generated by GMR do not remain during propagating in crystals. (c) 2006 Elsevier GmbH. All rights reserved.

Propagation and polarization properties of hollow Gaussian beams in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of hollow Gaussian beams (HGBs) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident HGBs propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of incident HGBs do not remain during propagating in crystals. (c) 2007 Published by Elsevier B.V.

地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

旋转轴对称弱各向异性介质模型的射线追踪方法

Ray tracing is a rapid and effective method for wave field calculation. Not only in the field of seismic-wave theory, but also in the field of seismic inversion and migration imaging,the seismic ray tracing method has become one of the most important methods. In anisotropic media, group velocity and phase velocity have different propagation directions. The seismic wave propagates along the direction of group velocity , it does not depend on the direction of phase velocity. Ray angle is a complex function with respect to phase angle, it is difficult to measure and calculate. But most rocks are weak anisotropic, so the expression of phase velocity can be simplified greatly. Based on the approximate expression of phase velocity this thesis for rotating axisymmetric weak anisotropic media deduces an expression of the partial derivative of phase velocity and an expression of group velocity with the method of linear approximation. This paper uses the fourth order Runge-Kutta method together with the two-dimensional interpolation and linear interpolation to obtain the parameters of the physical locations. At last the paths of seismic wave in rotating axisymmetric weak anisotropic media are computed. According to the analysis of the computational results, it indicates that the method developed in this paper has strong adaptability, high computational efficiency and high accuracy for rotating axisymmetric weak anisotropic media.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity