首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 611 毫秒
1.
Estimates of core density deficit (cdd) of the Earth's outer core recently reported by Anderson and Isaak [Another look at the core density deficit of Earth's outer core, Phys. Earth Planet Int. 131 (2002) 19-27] are questionable in view of the serious errors in the pressure-volume and bulk modulus data due to an inadequacy in the calibration process used by Mao et al. [Static compression of iron to 300 GPa and Fe0.8Ni0.2 alloy to 200 GPa: implications for the core, J. Geophys. Res. 94 (1990) 21737-21742]. The data used by Anderson and Isaak deviate significantly from the corresponding values derived from seismology. In the present study we have used the input data on density, isothermal bulk modulus and its pressure derivative from Stacey and Davis [High pressure equations of state with application to lower mantle and core, Phys. Earth Planet Int. 142 (2004) 137-184] which are consistent with the seismological data. Volumes of hexagonal close-packed iron have been calculated at different temperatures under isobaric conditions at P = 330 GPa, the inner core boundary (ICB) pressure using the relationship between thermal pressure and volume expansion based on the lattice potential theory originally due to Born and Huang [Dynamical Theory of Crystal Lattices, Oxford University Press, Oxford, 1954, p. 50]. The formulation for thermal pressure used by Anderson and Isaak has been modified by taking into account the variations of thermal expansivity α and isothermal bulk modulus KT with temperature. Values of cdd are then estimated corresponding to different temperatures ranging from 4000 to 8000 K. The results for cdd at different temperatures obtained in the present study are significantly higher than those estimated by Anderson and Isaak suggesting that the cdd for the Earth's outer core is nearly 10%. The effects of nickel when an Fe-Ni alloy replaces Fe are estimated and found to be insignificant.  相似文献   

2.
The results of solving the inverse problem of forced nutations and free oscillations of the Earth by decomposing the Q-factor and small depth variations in density in a system of orthogonal functions are considered. These functions are determined by orthogonalization of the functional derivatives of the observed parameters with respect to the depth distributions of the sought parameters (assuming there are no distributions of the velocities of body seismic waves V p and V S with depth and unchanged total mass M and inertia moments I of the Earth). The examples are presented to illustrate the numerical solution of the inverse problem on finding the density distributions in the mantle and core of the Earth using orthogonalization of the integral constraints for the probable depth distributions of density describing the conditions of unchanged M and I, as well as the constraints posed by the data on the periods of the free low-order oscillations of the Earth.  相似文献   

3.
The models of the mechanical Q factor of the inner solid core of the Earth are reconstructed from the amplitudes and phases of forced nutation and the periods and damping constants of the high-order overtones of free radial modes. The admissible range of the Q-factor in the solid core is estimated and the stability of the obtained distributions is analyzed. The real accuracy of the obtained model distributions is estimated on the basis of the previous conclusions concerning the uncertainty in the solution of the inverse problem of reconstructing the internal structure of the Earth in the low-frequency range.  相似文献   

4.
The question of ambiguity in the solution of the inverse problem for determining the Brünt-Väisäla frequency in the Earth’s mantle from the entire set of the up-to-date data on seismicity, free oscillations, and forced nutations of the Earth, as well as the data on the Earth’s total mass and total moment of inertia, is considered. Based on the results of a series of numerical experiments, the band of admissible distributions of the Brünt-Väisäla frequency and mantle density with depth is calculated. This estimate is used for investigating the convective and gravitational stability of the different regions of the mantle against relatively small adiabatic and nonadiabatic perturbations. The generalization of the known Rayleigh criterion of convective stability of homogeneous and a nonself-gravitating incompressible viscous fluid for the case of a compressible self-gravitating fluid is given. A system of the ordinary eight-order differential equations with complex coefficients and homogeneous boundary conditions, whose eigenvalues determine the transition from the stable state to instability, is obtained. Examples of the numerical determination of these eignevalues are presented. For interpreting the data about the band of the admissible distributions of the Brünt-Väisäla frequency with depth, the notion of the effective bulk modulus of the medium at different depths is introduced. This quantity governs the depth changes in temperature in a convecting mantle and allows us to make a conclusion about the role of heat conduction and the radial heterogeneity of the mantle composition without imposing any constraints on the convection mechanism. It is shown that within the present-day observation errors in the frequencies of the Earth’s free oscillations, the simplest reasonable model is that in which the ratio of the effective bulk modulus to its adiabatic value in the lower and middle mantle is 1.043 ± 0.05. The closeness of this value to unity indicates that convection in the lower and middle mantle is fairly close to adiabatic. At the same time, when the analysis only relies on seismic data and on the information about the periods of the free oscillations of the Earth, there is a significant uncertainty in the models of the effective bulk modulus distribution in the upper mantle and crust. This uncertainty precludes us from making purely empirically derived conclusions that reliably and unambiguously describe the role of the effects of heat conduction and radially heterogeneous material composition in the convection in the upper mantle.  相似文献   

5.
We analyze the present-day data on the periods of free oscillations and amplitudes of the forced nutations of the Earth for evaluating the admissible range of the mass and moment of inertia for the liquid core. The initial model for this study is taken in the form of the model distribution of density and mechanical Q parameters of the mantle suggested in (Molodenskii, 2010; 2011a; 2011b). This model was constructed by the steepest descent method in the space of 64 parameters, which determine the distribution of density and parameters of mechanical Q in the mantle, liquid outer core, and solid inner core of the Earth. We assumed the Q parameter of the mantle and inner solid core to be constant and sought for the density variations for the simplest two-parameter model of the piecewise-linear functions with the jumps on the boundary between the liquid core and the mantle and at the olivine-spinel phase transition at a depth of 670 km in the mantle. After this, the computations were repeated for the other distributions of Q (which were also assumed to be unchanged) that correspond to their limiting admissible values. Using this approach, we managed to find the most probable values of the mass and moment of inertia of the liquid core and determine the admissible range of their values. According to our estimates, the ratios of the mass and moments of inertia of the liquid core to the mass and moment of inertia of the whole Earth fall in the intervals 0.317996 ± 0.00065 and 0.110319 ± 0.00022, respectively. These values are lower than the corresponding values for the PREM model (0.322757 and 0.112297) by (1.48 ± 0.30)% and (1.76 ± 0.35)%, respectively. The interpretation of these results requires the revision and thorough analysis of the data on the admissible temperature range of the liquid core and (or) its chemical composition.  相似文献   

6.
—We investigate a new nonlinear inversion method for low frequencies to determine the bulk and shear modulus as well as the material density and the location of subsurface inhomogeneities. The solution is a direct exact nonlinear inversion of single scattered waves containing near- and far-field terms for incident P and scattered P and S waves, allowing for inversion of parameters in the vicinity and at distance from the sources and receivers. Because the approach is based on single scattering theory, the range of application includes single strong scattering anomalies of various sizes like magma chambers, gas- or fluid-filled cavities, or buried near-surface obstacles. The replacement of the material properties by a new set of parameters, referred to as scattering factors, allows the inverse problem to be solved analytically. The nonlinear nature of the scattering problem is investigated and implications for the inversion process are discussed. The deviations in the elastic parameters as a function of the scattering factors show a strong asymmetry about zero, and therefore linearized approximations will perform differently, depending on the sign of the perturbation. Based on the low frequency (Rayleigh) approximation, we introduce and evaluate a pair of approximations (Mie) derived by numerical and analytical integration of the Rayleigh approximation. Both approximations are based on the underlying principle of subdividing the inhomogeneities into a number of small noninteracting parts and subsequent integration over the total volume, thus increasing the Rayleigh limit and producing better resolution of the parameter estimates during the inversion. The two Mie approximations, when evaluated as a function of scattering angle and distance, produce similar results in the mid- and far-field of the inhomogeneity and reveal better resolution than the Rayleigh approximation. For three anomalies of ± 50% in bulk modulus, shear modulus, and density, the relative error between the exact solution and the two Mie approximations remains below 10%, 20%, and 30%, respectively, for values of k p R < 3.0, where R is the radius of the heterogeneity. However, smaller errors for individual cases are found for values up to k p R≈ 4.5. The performance of the inversion based on the analytically and the numerically integrated Mie approximation is tested for single parameter perturbations, revealing reliable and stable inversion results for the bulk and the shear modulus, reasonable results for the density, and crosstalk between the shear modulus and the density. The results show well-defined locations of the anomalies and slight deviations in the estimates of their magnitudes, which can be explained by amplitude and phase deviations between the analytical solution used for forward modeling and the approximations used for the inversion. The analytical Mie approximation provides a fast means to estimate elastic parameters compared to the more time consuming numerically integrated approximation, while the latter can be applied to more arbitrarily shaped inhomogeneities.  相似文献   

7.
The most abundant mineral on Earth has a perovskite crystal structure and a chemistry that is dominated by MgSiO3 with the next most abundant cations probably being aluminum and ferric iron. The dearth of experimental elasticity data for this chemically complex mineral limits our ability to calculate model seismic velocities for the lower mantle. We have calculated the single crystal elastic moduli (cij) for (Mg, Fe3 +)(Si, Al)O3 perovskite using density functional theory in order to investigate the effect of chemical variations and spin state transitions of the Fe3+ ions. Considering the favored coupled substitution of Mg2+-Si4 + by Fe3+-Al3+, we find that the effect of ferric iron on seismic properties is comparable with the same amount of ferrous iron. Ferric iron lowers the elastic moduli relative to the Al charge-coupled substitution. Substitution of Fe3+ for Al3+, giving rise to an Fe/Mg ratio of 6%, causes 1.8% lower longitudinal velocity and 2.5% lower shear velocity at ambient pressure and 1.1% lower longitudinal velocity and 1.8% lower shear velocity at 142 GPa. The spin state of the iron for this composition has a relatively small effect (< 0.5% variation) on both bulk modulus and shear modulus.  相似文献   

8.
Increasing amount of crystals tends to reduce the mobility of magmas and modifies its elastic characteristics (e.g. [Caricchi, L. et al., 2007. Non-Newtonian rheology of crystal-bearing magmas and implications for magma ascent dynamics. Earth and Planetary Science Letters, 264: 402–419.; Bagdassarov, N., Dingwell, D.B. and Webb, S.L., 1994. Viscoelasticity of crystal- and bubble-bearing rhyolite melts. Physics of the Earth and Planetary Interior, 83: 83–99.]). To quantify the effect of crystals on the elastic properties of magmas the propagation speed of shear and compressional waves have been measured at pressure and temperatures relevant for natural magmatic reservoirs. The measurements have been performed in aggregates at variable particle fractions (? = 0–0.7). The measurements were carried out at 200 MPa confining pressure and temperatures between 300 K and 1273 K (i.e. across the glass transition temperature (Tg) from glass to melt). The specimens were mixtures of a haplogranitic melt containing 5.25 wt.% H2O and variable amounts of sub-spherical alumina particles. Additional experiments were carried out on a sample containing both, crystals and air bubbles. The temperature derivatives of the shear (dVs/dT) and compressional wave (dVp/dT) velocities for pure glass and samples with a crystal fraction of 0.5 are different below and above the glass transition temperature. For a crystal fraction 0.7, only dVp/dT changed above the Tg. In the presence of gas bubbles, Vp and Vs decrease constantly with increasing temperature. The bubble-bearing material yields a lower bulk modulus relative to its shear modulus. The propagation velocities of compressional and shear waves increase non-linearly with increasing crystal fraction with a prominent raise in the range 0.5 < ? < 0.7. The speed variations are only marginally related to the density increase due to the presence of crystals, but are dominantly related to the achievement of a continuous crystal framework. The experimental data set presented here can be utilized to estimate the relative proportions of crystals and melt present in a magmatic reservoir, which, in turn, is one of the fundamental parameters determining the mobility of magma and, consequently, exerting a prime control on the likelihood of an eruption from a sub-surficial magma reservoir.  相似文献   

9.
A Gibson half-space model (a non-layered Earth model) has the shear modulus varying linearly with depth in an inhomogeneous elastic half-space. In a half-space of sedimentary granular soil under a geostatic state of initial stress, the density and the Poisson’s ratio do not vary considerably with depth. In such an Earth body, the dynamic shear modulus is the parameter that mainly affects the dispersion of propagating waves. We have estimated shear-wave velocities in the compressible Gibson half-space by inverting Rayleigh-wave phase velocities. An analytical dispersion law of Rayleigh-type waves in a compressible Gibson half-space is given in an algebraic form, which makes our inversion process extremely simple and fast. The convergence of the weighted damping solution is guaranteed through selection of the damping factor using the Levenberg-Marquardt method. Calculation efficiency is achieved by reconstructing a weighted damping solution using singular value decomposition techniques. The main advantage of this algorithm is that only three parameters define the compressible Gibson half-space model. Theoretically, to determine the model by the inversion, only three Rayleigh-wave phase velocities at different frequencies are required. This is useful in practice where Rayleigh-wave energy is only developed in a limited frequency range or at certain frequencies as data acquired at manmade structures such as dams and levees. Two real examples are presented and verified by borehole S-wave velocity measurements. The results of these real examples are also compared with the results of the layered-Earth model.  相似文献   

10.
In the first part of the paper, we obtained the refined estimates for the periods and Q-factors of the fundamental modes and overtones of spherical and toroidal oscillations with periods longer than 3 min from the data on the free oscillations of the Earth, which were excited by the earthquakes with magnitude 9 that occurred in Sumatra, Japan, and the Sea of Okhotsk. In (Molodenskii et al., 2013), we analyzed the limits of the admissible density distributions in the mantle and liquid core of the Earth, using the data on the amplitudes and phases of the forced nutations, as well as the periods and attenuation factors of the fundamental modes of the free spheroidal and toroidal oscillations of the Earth. These studies were conducted with the fixed values of the total mass and total moment of inertia of the Earth and the fixed distributions of the body seismic waves in the mantle and in the core. The solution was obtained by orthogonalizing the kernels of the integral equations that link the residuals of the observed frequencies and attenuation factors of the free oscillations, as well as the amplitudes and phases of the forced nutations, with the sought densities and Q-factors of the mantle and liquid core. Below, we present the solution of the same problem with allowance for the results obtained in the first part of this paper, namely, the new data on the periods and attenuation factors of the fundamental modes of free oscillations of the Earth and on the periods of the first four overtones of the free spheroidal and toroidal oscillations. Despite the involvement of the new data on the overtones, which have not been considered in our calculations, the weighted root mean square deviations of the theoretical predictions from the observed periods and attenuation factors of the free oscillations, as well as the amplitudes and phases of the forced nutations, have significantly decreased. This is due to (1) the noticeable reduction of the real errors in estimating the parameters of the free oscillations described in the first part of the paper and (2) the inclusion of the quantities determining the depth- and frequency dependences of the Q-factor in the mantle in the set of the independently varied parameters.  相似文献   

11.
Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. A shear wave propagating in an azimuthally anisotropic medium splits into two components with different polarizations if the source polarization is not aligned with the principal axes of the medium. For vertical propagation of shear waves in a horizontally layered medium containing vertical fractures, the shear‐wave splitting depends on the shear compliance of the fractures, but is independent of their normal compliance. If the fractures are not perfectly vertical, the shear‐wave splitting also depends on the normal compliance of the fractures. The normal compliance of a fluid‐filled fracture decreases with increasing fluid bulk modulus. For dipping fractures, this results in a decrease in shear‐wave splitting and an increase in shear‐wave velocity with increasing fluid bulk modulus. The sensitivity of the shear‐wave splitting to fluid bulk modulus depends on the interconnectivity of the fracture network, the permeability of the background medium and on whether the fracture is fully or partially saturated.  相似文献   

12.
Ambiguity in the inverse problem of retrieval of the mechanical parameters of the Earth’s shell and core from the set of data on the velocities V p and V S , of longitudinal and transverse seismic body waves, the frequencies f i and quality factors Q i , of free oscillations, and the amplitudes and phases of forced nutation is considered. The numerical experiments show that the inverse problem of simultaneous retrieval of the density profile ρ in the mantle-liquid core system and the mechanical quality factor Q μ of the mantle (if the total mass M and the total mean moment of inertia I of the Earth, and V p and V S are constant at all depths) has most unstable solutions. An example of depth distributions of ρ and Q μ which are alternative to the well-known PREM model is given. In these distributions, the values of M and I and the velocities V p and V S at all depths for the period of oscillations T = 1 s exactly coincide with their counterparts yielded by PREM model (T = 1 s); however, the maximum deviations of the ρ and Q μ profiles from those in the PREM model are about 3% and 40%, respectively; the mass and the moment of inertia of the liquid core are smaller than those for the PREM model by 0.75% and 0.63%, respectively. In this model, the root mean square (rms) deviations of all the measured values of f i and Q i from their values predicted by theory are half to third the corresponding values in the PREM model; the values of Δ for natural frequencies of the fundamental tone and overtones of radial oscillations, the fundamental tones of torsional oscillations, and the fundamental tones of spheroidal oscillations, which are measured with the highest relative accuracy, are smaller by a factor of 30, 6.6, and 2 than those in the PREM model, respectively. Such a large ambiguity in the solution of the inverse problem indicates that the current models of the depth distribution of density have relatively low accuracy, and the models of the depth distribution of the mechanical Q in the mantle are extremely unreliable. It is shown that the ambiguity in the models of depth distribution of density considerably decreases after the new data on the amplitudes and phases of the forced nutation of the Earth are taken into account. Using the same data, one may also refine by several times the recent estimates of the creep function for the lower mantle within a wide interval of periods ranging from a second to a day.  相似文献   

13.
Existing data supporting or disputing the validity of the Hashin-Shtrikman bounds on the elastic properties of multiphase aggregates often do not consider porosity, elastic anisotropy, or experimental errors. In this experiment, two-phase aggregates of KCl + (NH4Br, TlBr, CsCl, NaCl, Cu, and LiF) at every 20% volume fraction were vacuum hot-pressed and the compressional and shear velocities were measured with a computer-controlled ultrasonic interferometer to ±0.2%. The ratio of the shear moduli, μ, (phase 2/KCl) varied from about 1 to 5, producing a range of separations between the theoretical two-phase Hashin-Shtrikman bounds for the composites. Samples were generally 99% or better of the theoretical density, with less than 1% velocity anisotropy. Porosity corrections were applied assuming spherical pores, based on the observed velocity-pressure behaviour. Velocities agreed with the HS bounds calculated from the end-member single-crystal stiffnesses when anisotropy was taken into account.The velocity data were also used to estimate the bulk modulus, K, and shear modulus of the second phase by means of the matrix method — taking the K and μ of KCl as known and calculating the moduli of the other phase assuming that the measured velocities were the two-phase Hashin-Shtrikman bounds or the Voigt-Reuss-Hill average. A narrow range of moduli estimates results only if the μ's of both phases are fairly closely matched. For μ's mismatched by a factor of 5, the theoretical uncertainty in the estimates can be 10 times larger than the experimental uncertainty. Estimates using the VRH average can lie outside the HS-based results.  相似文献   

14.
A low strain shear modulus plays a fundamental role in earthquake geotechnical engineering to estimate the ground response parameters for seismic microzonation. A large number of site response studies are being carried out using the standard penetration test (SPT) data, considering the existing correlation between SPT N values and shear modulus. The purpose of this paper is to review the available empirical correlations between shear modulus and SPT N values and to generate a new correlation by combining the new data obtained by the author and the old available data. The review shows that only few authors have used measured density and shear wave velocity to estimate shear modulus, which were related to the SPT N values. Others have assumed a constant density for all the shear wave velocities to estimate the shear modulus. Many authors used the SPT N values of less than 1 and more than 100 to generate the correlation by extrapolation or assumption, but practically these N values have limited applications, as measuring of the SPT N values of less than 1 is not possible and more than 100 is not carried out. Most of the existing correlations were developed based on the studies carried out in Japan, where N values are measured with a hammer energy of 78%, which may not be directly applicable for other regions because of the variation in SPT hammer energy. A new correlation has been generated using the measured values in Japan and in India by eliminating the assumed and extrapolated data. This correlation has higher regression coefficient and lower standard error. Finally modification factors are suggested for other regions, where the hammer energy is different from 78%.  相似文献   

15.
A critical porosity model establishes the empirical relationship between a grain matrix and a dry rock by the concept of critical porosity. The model is simple and practical and widely used. But the critical porosity in the model is a fixed value that cannot relate to pore structure. The aim of this paper is to establish the theoretical relationship between critical porosity and pore structure by combining Kuster–Toksöz theory with the critical porosity model. Different from the traditional critical porosity model, critical porosity is not an empirical value but varied with pore shape and the ratio of bulk modulus versus shear modulus of the grain matrix. The substitution of the theoretical relationship into Kuster–Toksöz theory will generate the formulae for the bulk and shear moduli of multiple-porosity dry rocks, which is named the multiple-porosity variable critical porosity model. The new model has been used to predict elastic moduli for sandstone and carbonate rock. We compare the modelling results for P- and S-wave velocities and elastic moduli with the experimental data. The comparison shows that the new model can be used to describe the elastic properties for the rocks with multiple pore types.  相似文献   

16.
耿煜  王君恒 《地球物理学报》2015,58(10):3530-3539
地球形成初期,构成地球的物质在组成上是大致均一的.目前地球的地核-地幔-地壳圈层结构,是由分异作用形成的.分异过程释放的能量称为分异能.Sorokhtin和Chilingarian等人从行星吸积的定义出发,导出了基于地球内部密度分布的势能计算公式,计算出的分异能大小为1.698×1031J.本文采用计算球体势能的思路,导出分异能计算的解析公式和数值计算公式,通过求取原始地球模型与均匀分层模型、PREM模型的势能差计算分异能.两种方法的计算结果分别为1.535×1031J和1.698×1031J.前者与Sorokhtin等的结果相近,后者与之相同.本文初步分析了方法间的异同以及造成结果偏差的主要原因.  相似文献   

17.
深内部地球结构对内核平动振荡本征周期的影响   总被引：2，自引：2，他引：0       下载免费PDF全文
地球固态内核的平动振荡是地球的基本简正模之一,又称Slichter模,其本征周期大约为几个小时,与地球内部结构密切相关.为了研究影响内核平动振荡的本征周期与内部结构的依赖关系,本文利用球对称、非自转、弹性和各向同性地球模型（SNREI）,通过自由振荡运动方程的数值积分,以地球模型PREM为基础,理论上系统研究了地球内部介质（包括密度、地震波速等）分布异常对Slichter模本征周期的影响.数值结果表明,Slichter模周期随着内外核边界（ICB）密度差的增加以类似于双曲线的特征显著减小,当ICB密度差从597 kg·m-3减小到200 kg·m-3时,周期增大66.44%,当ICB密度差从597 kg·m-3增大到1000 kg·m-3时,周期减小21.48%;Slichter模周期随着核幔边界（CMB）密度差的增大而缓慢增大;相对于PREM,地球模型1066A在ICB和CMB的密度差分别相差45.321%和1.132%,内部地震波速度和密度梯度也存在差异,但是,当密度差减小到1066A模型提供的数值时,得到的Slichter模周期与基于1066A获得的结果（4.599 h）非常接近,差异分别只有3.762%和0.037%;表明Slichter模本征周期与地球内部介质的精细结构关系不大,而对ICB的密度差非常敏感.内、外核P波波速分布异常对Slichter模周期的影响基本相当,当内核和外核P波波速均增加5%时,Slichter周期分别减小1.02%和1.69%,P波波速分别减小5%时,Slichter模周期分别增加1.27%和1.847%,内核S波波速分布异常比P波波速分布异常对Slichter模周期的影响小1个量级;与地核相比,地幔中的地震波速异常对Slichter模本征周期的影响小1~2个量级;表明地核中地震波速异常对Slichter模周期的影响很小,目前有关Slichter模周期理论计算的差异主要来自于所采用的地球模型中内核边界的密度差的差异,本文结果可以为Slichter模的研究、探测及其对地球深内部结构的约束提供理论依据.  相似文献   

18.
The aim of the present study is twofold. Firstly, the paper investigates the undrained cyclic and post-cyclic behaviour of two silica sands by means of multi-stage cyclic triaxial tests. Secondly, based on the post-cyclic response observed in the element test, the authors formulate a simplified stress–strain relationship that can be conveniently used for the construction of p–y curves for liquefiable soils. The multi-stage loading condition consists of an initial cyclic loading applied to cause liquefaction, followed by undrained monotonic loading that aimed to investigate the post-cyclic response of the liquefied sample. It was found that due to the tendency of the liquefied soil to dilate upon undrained shearing, the post-liquefaction strain–stress response was characterised by a distinct strain–hardening behaviour. The latter is idealized by means of a bi-linear stress–strain model, which can be conveniently formulated in terms of three parameters, i.e.: (i) take-off shear strain, γto, i.e. shear strain required to mobilize 1 kPa of shear strength; (b) initial secant shear modulus, G1, defined as 1/γto; (c) post-liquefied shear modulus at large strain, G2 (γ⪢γto). Based on the experimental results, it is concluded that these parameters are strongly influenced by the initial relative density of the sample, whereby γto decreases with increasing relative density. Differently both shear moduli (G1 and G2) increases with increasing relative density. Lastly, the construction of new p–y curves for liquefiable soils based on the idealized bi-linear model is described.  相似文献   

19.
基于叠前地震纵横波模量直接反演的流体检测方法   总被引：15，自引：9，他引：6       下载免费PDF全文
流体因子是储层流体识别的重要方法,而叠前地震反演是获得流体因子的有效途径之一.本文从流体因子的构建出发,基于多孔弹性介质岩石物理模型,建立了流体因子与纵横波模量之间的直接关系,避免了流体因子计算所需的密度参数无法准确求取的问题,通过推导基于纵横波模量的Zeoppritz近似方程及弹性阻抗方程,探讨了基于弹性阻抗的纵横波模量直接反演方法,模型与实际应用表明,基于弹性阻抗的纵横波模量直接反演方法合理、可靠,减少了常规方法间接计算纵横波模量带来的累积误差,基于纵横波模量的流体因子计算方法有较好的实际应用效果.  相似文献   

20.
Short-period seismograms are synthesized for PKP phases in anelastic Earth models. The synthetics were constructed using a synthetic technique valid at grazing incidence, a source-time function appropriate for deep-focus earthquakes, and an instrument response for either a short-period WWSSN or SRO seismograph. The agreement between predicted and observed amplitudes and spectral ratios requires neither a low-Qα zone at 0.2–2 Hz nor a low or negative P-velocity gradient at the bottom of the outer core. Thin low-Qα zones beneath the inner core boundary fit spectral ratio data that sample the upper 200 km of the inner core but fail to fit data that sample the lower inner core. Only a model having Qα?1?[0.003, 0.004] at 0.2–2 Hz, nearly constant with depth in the inner core, satisfies all of the spectral ratio and amplitude data. The assumption of a bulk viscosity of 10-103 Pa s for the liquid phase of a partially molten inner core combined with the observation of low shear attenuation in the inner core at frequencies less than 0.005 Hz limit the physical parameters associated with two possible attenuation mechanisms: (1) fluid flow and viscous relaxation due to ellipsoidally shaped inclusions of melt, and (2) the solid-liquid phase transformation induced by the stress change during the passage of a seismic wave. Both mechanisms require an order of 0.1% partial melt to reproduce the observed Qα?1. In the outer core, the time constant of the mechanism of phase transformation is predicted to be 104–106 s. Confirmation of small shear attenuation in the inner core in the frequency band of seismic body waves would favor the mechanism of phase transformation.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号