共查询到20条相似文献，搜索用时 35 毫秒
1.
Ernian Pan 《Physics of the Earth and Planetary Interiors》1990,60(1-4):254-264
A general method is developed for the study of transient thermoelastic deformation in a transversely isotropic and layered half-space by surface loads and internal sources. A Laplace transform is first applied to the field quantities; Cartesian and cylindrical systems of vector functions are then introduced for reducing the basic equations to three sets of simultaneous linear differential equations. General solutions are obtained from these sets, and propagator matrices from the solutions by a partitioned matrix method.
Source functions for a variety of sources are derived in the Cartesian and cylindrical systems, and the Laplace transformed expressions of the field variables at the surface presented explicitly in the two systems in terms of a layer matrix. The effect of gravity is included by multiplying simply an effect matrix resulting from the modification of continuity conditions at the surface and the layer interfaces.
It should be noted that the present analytical method has great advantages over either the classical thin plate approach or the finite element method, and that the present result can be reduced directly to the solutions of the corresponding isotropic case. 相似文献
2.
Tarkeshwar Lal 《Pure and Applied Geophysics》1970,82(1):259-269
Summary The theory of resistivity sounding on a three layer earth comprising an inhomogeneous interstratum has been presented. The electrical conductivity in the intermediate layer, embedded between two layers of uniform conductivity, is assumed to follow either generalised power law or exponential variation. A point source of current is placed on the surface. The electric potential functions on the ground surface have been obtained. The expressions for apparent resistivity over the above models for the Wenner and Schlumberger electrode configurations have been derived. 相似文献
3.
Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, have not been simultaneously satisfied by an isotropic structure. In this paper available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. We assume that the surface wave data represents an azimuthal average of actual velocities. Thus, we can treat the mantle as transversely isotropic. The resulting models for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SH > SV in the vicinity of the low-velocity zone. Neat the East Pacific Rise, however, SV > SH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of ～50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures. A new model for average Earth is obtained which includes a thin LID. This model extends the fit of a PREM, type model to shorter period surface waves. 相似文献
4.
An earth model with a transition layer (anisotropic inhomogeneous) is considered. The inhomogeneity in σv (vertical conductivity) of the transition layer is represented by a power law variation. Expressions for potential distribution in the upper layer, transition layer and bottom layer are obtained by solving appropriate differential equation for each layer. By utilizing the boundary conditions, expressions of apparent resistivity for Wenner and Schlumberger configurations are derived. Numerical analysis is performed for linear and quadratic variation of σv. The results are presented in the form of theoretical apparent resistivity curves for both configurations. Negative apparent resistivities are the interesting feature of this analysis. 相似文献
5.
Akinori Hashima Yukitoshi Fukahata Chihiro Hashimoto Mitsuhiro Matsu’ura 《Pure and Applied Geophysics》2014,171(8):1669-1693
We derived explicit expressions in the time domain for 3-D quasi-static strain and stress fields, due to a point moment tensor source in an elastic surface layer overlying viscoelastic half-space under gravity. The expressions of strain in the elastic surface layer were directly obtained from the expressions of displacement in our previous paper. The conversion of strain into stress is easy, because the stress–strain relation of elastic material is linear. In the viscoelastic substratum, the expressions of strain were obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. The strain–stress conversion is not straightforward, as the stress–strain relation of viscoelastic material is usually given in a differential form. To convert strain into stress, we used an integral form of the stress–strain relation instead of the usual differential form. The expressions give the responses of elastic half-space at \( t = 0 \) , and the responses of an elastic plate floating on non-viscous liquid at \( t = \infty \) . The moment tensor is rationally decomposed into the three independent force systems, corresponding to isotropic expansion, shear faulting and crack opening, and so the expressions include the strain and stress fields for these force systems as special cases. As the first numerical example, we computed the temporal changes in strain and stress fields after the sudden opening of an infinitely long vertical crack cutting the elastic surface layer. Here, we observe that the stress changes caused by the sudden crack opening gradually decay with time and vanish at \( t = \infty \) everywhere. After the completion of stress relaxation, a characteristic pattern of shear strain remains in the viscoelastic substratum. Since the strain and stress fields at \( t = \infty \) can be read as the strain- and stress-rate fields caused by steady crack opening, respectively, this numerical example demonstrates the realization of a steady stress state supported by steady viscous flow in the asthenosphere, associated with steady seafloor spreading at mid-ocean ridges. For the second numerical example, we computed the temporal changes in strain and stress fields after the 2011 Tohoku-oki mega-thrust earthquake, which occurred at the North American-Pacific plate interface. In this numerical example, the stress changes caused by coseismic fault slip vanish at \( t = \infty \) in the viscoelastic substratum, but remain in the elastic surface layer. The coseismic stress changes (and also strain changes) in the elastic surface layer diffuse away from the source region with time, due to gradual stress relaxation in the viscoelastic substratum. 相似文献
6.
The electrical potential due to a point source of current placed on the ground surface is studied for a multi-layered earth consisting of homogeneous overburden of constant conductivity over a stack of transition layers where conductivity varies with depth according to power and exponential laws in even and odd layers, respectively. The general recursion relations are derived and are used to obtain expressions for the apparent resistivities for Schlumberger and Wenner electrode arrays. Their asymptotic behavior has been studied. The solutions for some particular cases are given: (i) odd layers with conductivities exponentially varying with depth while all even layers (and the first) have constant conductivities; (ii) even layer conductivities varying as a power law while odd number layers are of constant conductivity; and (iii) any two successive layers as transition layers and all others having constant conductivities. Further it is shown that Patella's theory is a particular case of the present study. It is concluded, therefore, that the present treatment is more general as all earlier models consisting of trasition layers can be derived from this study. 相似文献
7.
Dynamic analysis of slab track on multi-layered transversely isotropic saturated soils subjected to train loads 总被引：1，自引：0，他引：1
The dynamic responses of a slab track on transversely isotropic saturated soils subjected to moving train loads are investigated by a semi-analytical approach. The track model is described as an upper Euler beam to simulate the rails and a lower Euler beam to model the slab. Rail pads between the rails and slab are represented by a continuous layer of springs and dashpots. A series of point loads are formulated to describe the moving train loads. The governing equations of track-ground systems are solved using the double Fourier transform, and the dynamic responses in the time domain are obtained by the inverse Fourier transform. The results show that a train load with high velocity will generate a larger response in transversely isotropic saturated soil than the lower velocity load, and special attention should be paid on the pore pressure in the vicinity of the ground surface. The anisotropic parameters of a surface soil layer will have greater influence on the displacement and excess pore water pressure than those of the subsoil layer. The traditional design method taking ground soil as homogeneous isotropic soil is unsafe for the case of RE 1 and RG 1, so a transversely isotropic foundation model is of great significance to the design for high train velocities. 相似文献
8.
S. K. Ghosh 《Pure and Applied Geophysics》1973,105(1):781-801
Summary In this paper approximate expressions for the surface displacements due to a buried spherical source inside an isotropic elastic layer with one surface stress-free and the other (i) rigid, (ii) stress-free, have been obtained. Graphical representations for the surface displacements at a large distance have also been found. 相似文献
9.
A. Sidorov A. Bakulin B. Kashtan S. Ziatdinov D. Alexandrov 《Geophysical Prospecting》2009,57(5):863-882
Many tasks in geophysics and acoustics require estimation of mode velocities in cylindrically layered media. For example, acoustic logging or monitoring in open and cased boreholes need to account for radial inhomogeneity caused by layers inside the borehole (sand screen, gravel pack, casing) as well as layers outside (cement, altered and unaltered formation layers). For these purposes it is convenient to study a general model of cylindrically layered media with inner fluid layer and free surface on the outside. Unbounded surrounding media can be described as a limiting case of this general model when thickness of the outer layer is infinite. At low frequencies such composite media support two symmetric modes called Stoneley (tube) and plate (extensional) wave. Simple expressions are obtained for these two mode velocities valid at zero frequency. They are written in a general form using elements of a propagator matrix describing axisymmetric waves in the entire layered composite. This allows one to apply the same formalism and compute velocities for n -layered composites as well as anisotropic pipes. It is demonstrated that the model of periodical cylindrical layers is equivalent to a homogeneous radially transversely isotropic media when the number of periods increases to infinity, whereas their thickness goes to zero. Numerical examples confirm good validity of obtained expressions and suggest that even small number of periods may already be well described by equivalent homogeneous anisotropic media. 相似文献
10.
Summary Theoretical expressions are derived for the displacement, velocity and stress in the time domain induced by an axially symmetric shearing stress applied at the inner surface of a spherical cavity in a hornogeneous, isotropic, elastic medium of infinite extent. Theoretical seismograms are computed for a step source and for three sources with exponential decay in time. A satisfactory time-dependence of the source can be obtained by combining the step source with one or more exponentially decaying sources. 相似文献
11.
For curves over a horizontally stratified earth where any of the layers has a conductivity variation proportional to (1 +β·z)N, where β and N are arbitrary constants and z is the depth to the layer, expressions for apparent resistivity for Wenner and Schlumberger sounding are derived. No assumption has been made about the continuity of conductivity at the interface. It is shown that most of the previous investigations in this connection can be regarded as particular cases of the present study. 相似文献
12.
一、引言 地震面波的频散性质、地震波辐射的方向性等特性已经广泛地用于地壳和上地幔结构以及震源机制的研究中,并且取得了许多有用的成果.研究地震波如何从震源辐射出来、如何在实际介质中传播和衰减的这一问题,对于利用地震波确定地壳和上地幔结构以及震源的参数,是很有必要的.关于这一问题的研究,已经作过许多工作。已往的工作中,为了分析方便,往往采用简单的地壳-上地幔模型或简单的震源模型,或两者都相 相似文献
13.
Perturbation methods are common tools for describing wave propagation in weakly anisotropic media. The anisotropic medium is replaced by an average isotropic medium where wave propagation can be treated analytically and the correction for the effect of anisotropy is computed by perturbation techniques. This works well for anisotropies of up to 10%. Some materials (e.g. shales), however, can exhibit a much stronger anisotropy. In this case a background is required which still can be treated analytically but is applicable to stronger P-wave anisotropy. We present an averaging technique to compute a best-fitting ellipsoidal medium to an arbitrary anisotropic medium. Ellipsoidal media are sufficiently simple for analytical expressions to be available for many applications and allow consideration of strong P-wave anisotropy. The averaging of the arbitrary anisotropic medium can be carried out globally (i.e. for the whole sphere) or sectorially (e.g. for seismic waves propagating predominantly in the vertical direction). We derive linear relationships for the coefficients of the ellipsoid which depend on the elastic coefficients of the anisotropic medium. We also provide specifications for best-fitting elliptical and best-fitting isotropic media. Numerical examples for different rocks demonstrate the improved approximation of the anisotropic model obtained using the formulae derived, compared with the conventionally used average isotropic medium. 相似文献
14.
Dispersion of Rayleigh-type surface wave is studied in a homogeneous transversely isotropic elastic layer overlying a nondissipative liquid-saturated porous solid half-space and lying under a uniform layer of homogeneous liquid. The frequency equation in the form of ninth-order determinant is obtained.Special cases have been deduced by reducing the depth of the layers to zero and by changing the transverse isotropic layer to an isotropic layer. Dispersion curves for the phase velocity have been plotted for a particular model. 相似文献
15.
Robert Bloor 《Geophysical Prospecting》1996,44(2):197-213
The near-surface layer is modelled as a constant-velocity layer with varying thickness. The base of the layer is described by a B-spline curve. The optimum model is calculated by minimizing, with respect to the model parameters, the difference between traveltimes predicted by the model and those observed in the data. Once a model has been produced, corrections that are dependent on the raypath geometry through the near-surface layer can be calculated. The effect of the near-surface layer is normally considered to be consistent at each shot or geophone station for all traveltimes arriving at that location (the surface-consistent approximation). This assumption linearizes the problem, allowing timeshifts to be calculated and the traveltimes corrected to a chosen datum, representing static corrections. The single correction at each point is an averaged correction, based on an assumption that is particularly inaccurate in the presence of lateral variations of velocity or thickness of the surface layer, in the presence of large surface layer velocities or in the presence of a thick surface layer. The method presented considers the non-linear relationship between data and model explicitly, hence the correction that is dependent on the raypath. Linearization removes this dependence and reduces the problem to a surface-consistent approximation. The method is applied to synthetic data calculated from a model with surface layer variations. Comparisons are made between the corrected data resulting from the method described here and the conventional surface-consistent approach. From these results it becomes apparent that the near-surface layer inversion method presented here can reproduce accurate models and correct for near-surface layer effects in cases where conventional methods encounter difficulties. Additionally the method can be readily extended to 3D. 相似文献
16.
台址下表层介质电阻率变化对视电阻率观测值影响特征的初步研究 总被引：2，自引：0，他引：2
气候的交替变化将导致大地表层含水量发生变化，从而使该层介质的电阻率发生变化，进而使得观测结果(即视电阻率)发生变化，显然这对地震前兆观测而言是一种干扰。目前我国绝大多数前兆地电阻率台采用单极距观测方式，因此观测到的视电阻率变化应是表层变化(主要是干扰)和观测范围内深部变化(地下介质应力与应变状态发生变化)的综合反映。作先依据台址下的水平层状模型进行了模拟计算研究，发现大地表层含水量的变化对单极距四极装置视电阻率观测值的影响有两个特征：(1)无论台址下电阻率结构是下伏低阻还是下伏高阻，表层介质电阻率升高对观测值的影响都很小；(2)表层电阻率降低对观测值的影响很明显，且对下伏高阻结构的影响大于对下伏低阻结构的影响。然后用视电阻率实际观测资料说明了将其特征用于预报的可能性。 相似文献
17.
18.
THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA 
下载免费PDF全文
下载免费PDF全文 In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research. 相似文献
19.
可控源音频大地电磁测深法(CSAMT)基本上是基于电性各向同性理论的假设,而实际地层普遍表现为电性各向异性现象,因此研究各向异性地层条件下CSAMT方法的电磁响应具有重要的理论意义和应用价值.本文从有源Maxwell方程出发,推导出了二维各向异性地层条件下无限长线源频率域的电磁场响应,得到了一组平行于线源方向的电场分量E_x和磁场分量H_x的偏微分方程,并采用有限元法进行求解.模型的计算结果表明:不论是对均匀半空间还是一维层状介质中的各向异性二维体模型的计算,均能从结果中明显地看出各向同性背景场和各向异性异常体的存在;在计算中通过改变各向异性系数或各向异性旋转角,其结果在视电阻率和相位曲线上均存在显著区别.计算结果表明了在各向异性地层中,仍然采用各向同性的假设进行CSAMT资料的处理与解释,可能会带来比较严重的误差,这种情况下必须采用基于各向异性的地质模型对CSAMT资料进行处理和解释.研究结果对于提高CSAMT方法的勘探效果和应用水平具有重要的指导意义. 相似文献
20.
Min Liang Denis Marcotte 《Stochastic Environmental Research and Risk Assessment (SERRA)》2016,30(3):973-987
This article describes the use of non-stationary covariance functions with compact support to estimate and simulate a random function. Based on the kernel convolution theory, the functions are derived by convolving hyperspheres in \(\mathbb{R}^n\) followed by a Radon transform. The order of the Radon transform controls the differentiability of the covariance functions. By varying spatially the hyperspheres radius one defines non-stationary isotropic versions of the spherical, the cubic and the penta-spherical models. Closed-form expressions for the non-stationary covariances are derived for the isotropic spherical, cubic, and penta-spherical models. Simulation of the different non-stationary models is easily obtained by weighted average of independent standard Gaussian variates in both the isotropic and the anisotropic case. The non-stationary spherical covariance model is applied to estimate the overburden thickness over an area composed of two different geological domains. The results are compared to the estimation with a single stationary model and the estimation with two stationary models, one for each geological domain. It is shown that the non-stationary model enables a reduction of the mean square error and a more realistic transition between the two geological domains. 相似文献