首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 296 毫秒
1.
Saturation of porous rocks with a mixture of two fluids (known as partial saturation) has a substantial effect on the seismic waves propagating through these rocks. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves, due to wave-induced fluid flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. As partial fluid saturation can occur on different length scales, attenuation due to wave-induced fluid flow is ubiquitous. In particular, mesoscopic fluid flow due to heterogeneities occurring on a scale greater than porescale, but less than wavelength scale, is responsible for significant attenuation in the frequency range from 10 to 1000 Hz.Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a periodic/regular way. In 1D this corresponds to periodically alternating layering, in 3D as periodically distributed inclusions of a given shape (usually spheres). All these models yield very similar estimates of attenuation and dispersion.Experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Therefore, theoretical models are required which would simulate the effect of more general and realistic fluid distributions.We have developed two theoretical models which simulate the effect of random distributions of mesoscopic fluid heterogeneities. The first model assumes that one fluid forms a random ensemble of spherical inclusions in a porous medium saturated by the other fluid. The attenuation and dispersion predicted by this model are very similar to those predicted for 3D periodic distribution. Attenuation (inverse quality factor) is proportional to ω at low frequencies for both distributions. This is in contrast to the 1D case, where random and periodically alternating layering shows different attenuation behaviour at low frequencies. The second model, which assumes a 3D continuous distribution of fluid heterogeneities, also predicts the same low-frequency asymptote of attenuation. However, the shapes of the frequency dependencies of attenuation are different. As the 3D continuous random approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental results. Further research is required in order to uncover how to relate the random functions to experimentally significant parameters.  相似文献   

2.
The increasing popularity of simplified nonlinear methods in seismic design has recently led to many proposals for procedures aimed at extending pushover analysis to plan asymmetric structures. In terms of practical applications, one particularly promising approach is based on combining pushover analysis of a 3D structural model with the results of linear (modal) dynamic analysis. The effectiveness of such procedure, however, is contingent on one fundamental requirement: the elastic prediction of the envelope of lateral displacements must be conservative with respect to the actual inelastic one. This paper aims at verifying the above assumption through an extensive parametric analysis conducted with simplified single‐storey models. The main structural parameters influencing torsional response in the elastic and inelastic range of behaviour are varied, while devoting special attention to the system stiffness eccentricity and radius. The analysis clarifies the main features of inelastic torsional response of different types of building structures; in this manner, it is found that the above‐mentioned method is generally suitable for structures characterized by moderate to large torsional stiffness, whereas it cannot be recommended for extremely torsionally stiff structures, as their inelastic torsional response almost always exceeds the elastic one. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

3.
The high level of seismic activity in eastern Turkey is thought to be mainly associated with the continuing collision of the Arabian and Eurasian tectonic plates. The determination of a detailed three-dimensional (3D) structure is crucial for a better understanding of this on-going collision or subduction process; therefore, a body wave tomographic inversion technique was performed on the region. The tomographic inversion used high quality arrival times from earthquakes occurring in the region from 1999 to 2001 recorded by a temporary 29 station broadband IRIS-PASSCAL array operated by research groups from the Universities of Bo?azi?i (Turkey) and Cornell (USA). The data was inverted and consisted of 3,114 P- and 2,298 S-wave arrival times from 252 local events with magnitudes (M D) ranging from 2.5 to 4.8. The stability and resolution of the results were qualitatively assessed by two synthetic tests: a spike test and checkerboard resolution test and it was found that the models were well resolved for most parts of the imaged domain. The tomographic inversion results reveal significant lateral heterogeneities in the study area to a depth of ~20?km. The P- and S-wave velocity models are consistent with each other and provide evidence for marked heterogeneities in the upper crustal structure beneath eastern Turkey. One of the most important features in the acquired tomographic images is the high velocity anomalies, which are generally parallel to the main tectonic units in the region, existing at shallow depths. This may relate to the existence of ophiolitic units at shallow depths. The other feature is that low velocities are widely dispersed through the 3D structure beneath the region at deeper crustal depths. This feature can be an indicator of the mantle upwelling or support the hypothesis that the Anatolian Plateau is underlain by a partially molten uppermost mantle.  相似文献   

4.
We present a new, fast and versatile method, the lateral parameter correlation method, of invoking lateral smoothness in model sections of one-dimensional (1D) models. Modern, continuous electrical and electromagnetic methods are capable of recording very large data sets and except for a few cases, standard inversion methodology still relies on 1D models. In environments where the lateral rate of change of resistivity is small, 1D inversion can be justified but model sections of concatenated 1D models do not necessarily display the expected lateral smoothness.
The lateral parameter correlation method has three steps. First, all sounding data are inverted individually. Next, a laterally smooth version of each model parameter, one at a time, is found by solving a simple constrained inversion problem. Identity is postulated between the uncorrelated and correlated parameters and the equations are solved including a model covariance matrix. As a last step, all sounding data are inverted again to produce models that better fit the data, now subject to constraints by including the correlated parameter values as a priori values. Because the method separates the inversion from the correlation it is much faster than methods where the inversion and correlation are solved simultaneously, typically with a factor of 200–500.
Theoretical examples show that the method produces laterally smooth model sections where the main influence comes from the well-determined parameters in such a way that problems with equivalence and poor resolution are alleviated. A field example is presented, demonstrating the improved resolution obtained with the lateral parameter correlation method. The method is very flexible and is capable of coupling models from inversion of different data types and information from boreholes.  相似文献   

5.
Seismic studies of the lowermost mantle suggest that the core-mantle boundary (CMB) region is strongly laterally heterogeneous over both local and global scales. These heterogeneities are likely to be associated with significant lateral viscosity variations that may influence the shape of the long-wavelength non-hydrostatic geoid. In the present paper we investigate the effect of these lateral viscosity variations on the solution of the inverse problem known as the inferences of viscosity from the geoid. We find that the presence of lateral viscosity variations in the CMB region can significantly improve the percentage fit of the predicted data with observations (from 42 to 70% in case of free-air gravity) while the basic characterisics of the mantle viscosity model, namely the viscosity increase with depth and the rate of layering, remain more or less the same as in the case of the best-fitting radially symmetric viscosity models. Assuming that viscosity is laterally dependent in the CMB region, and radially dependent elsewhere, we determine the largescale features of the viscosity structure in the lowermost mantle. The viscosity pattern found for the CMB region shows a high density of hotspots above the regions of higher-than-average viscosity. This result suggests an important role for petrological heterogeneities in the lowermost mantle, potentially associated with a post-perovskite phase transition. Another potential interpretation is that the lateral viscosity variations derived for the CMB region correspond in reality to lateral variations in the mechanical conditions at the CMB boundary or to large-scale undulations of a chemically distinct layer at the lowermost mantle.  相似文献   

6.
Summary This paper studies the propagation of Love waves in the following two-layered models: (i) a heterogeneous layer overlying a heterogeneous halfspace, (ii) a heterogeneous layer overlying homogeneous halfspace, (iii) a homogeneous layer overlying a homogeneous halfspace. The heterogeneities of the layer and halfspace are specified by different exponential functions of the depth. The dispersion curves have been drawn taking account of the actual variation of elastic parameters with depth in the Earth. The effects of the heterogeneities of the layer (crust) and halfspace (mantle) on dispersion curves have been discussed.  相似文献   

7.
Several two-dimensional structures are modelled for vertical electrical soundings in arrays parallel and perpendicular to the strike of the structure. The models are a horst and a graben within a three-layer medium, a cliff over two layers, and heterogeneities around the electrodes. Apparent resistivity curves are shown for different model parameters and different distances to the two-dimensional structures. Some of the features on the shape of these curves are inflections that may be misinterpreted as fictitious layers, of slopes greater than 45°; some features are simply anomalous peaks. One-dimensional interpretation of a two-dimensional graben model has been performed, in order to show typical errors when 2D structures are interpreted as one-dimensional. A real case corresponding to a landfill near Barcelona and showing strong lateral and topographic effects is presented. This represents combined effects of the above theoretical 2D models.  相似文献   

8.
M.A. Soofi  P. Wu   《Journal of Geodynamics》2008,46(1-2):38-47
The region of Alaska and adjacent northwest Canada is tectonically active and is subjected to multiple tectonic processes including plate subduction and terrane accretion. These tectonic processes and the forces originating thereof are responsible for high seismicity in the region and deformation of the crust. In the present-day tectonic setting, the Yakutat terrane is obliquely colliding with Alaska along the Aleutian Trench. Also, flat subduction due to under thrusting of a thickened crust, probably of oceanic affinity, is contributing to the tectonic evolution of this region in a basal traction collision style. This study uses the 2D, planform, thin-viscous-sheet model to investigate the effect of the Yakutat terrane colliding with Alaska and adjacent northwest Canada. Along with the obliquity and velocity of convergence, the lateral strength heterogeneities in the crust are considered in this investigation. The results of the numerical model are constrained with the observed topography and stress orientation in Alaska. It is shown that the Alaska–Yakutat collision is producing asymmetric deformation of the crust with respect to the normal to the collision boundary and that lateral strength heterogeneities contribute significantly to the deformation of the crust. Also, the influence of this collision can be observed up to a distance of 700 km inland from the collision boundary.  相似文献   

9.
The authors have presented the results of inelastic earthquake response analysis of a class of asymmetric building models and concluded that regularly asymmetric buildings excited well into the inelastic range may not be conservatively designed via modal analysis of their elastic models. The purpose of this short communication is to refer to the results in contrast to the authors' and to make some comments on their paper.  相似文献   

10.
The paper deals with the applicability of the extended N2 method to base‐isolated plan‐asymmetric building structures. The results obtained by the usual pushover analysis of a 3D structural model are further combined with the aid of linear dynamic (spectral) analysis to account for the dynamic effects caused by structural asymmetry. In the paper, the method has been applied to the seismic analysis of a mass‐eccentric four‐storeyed RC frame building isolated with lead rubber bearings. Three different positions of the center of isolation system (CI) with respect to the center of mass (CM) and the center of stiffness of the superstructure (CS) were considered. The response was analyzed for three different eccentricities, three different torsional to lateral frequency ratios of the superstructure, and two ground motion intensities. The stiffness of the isolators was selected for three different protection levels, which resulted in elastic as well as moderately to excessively damaged superstructure performance levels. The results are presented in terms of the top, base and relative displacements, as well as the stiff/flexible side amplification factors. A more detailed insight into the nonlinear behavior of the superstructure is given in a form of ductility factors for the flexible and stiff side frames. The results of the extended N2 method for selected lateral load distributions are compared with the average results of nonlinear dynamic analyses. It was concluded that the extended N2 method could, with certain limitations, provide a reasonable prediction of the torsional influences in minor to moderately asymmetric base‐isolated structures. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

11.
Understanding the processes that occur in the transition from the Pacific Ocean to Eurasia is key to constructing the tectonic models of the Earth’s shells and the convection models of the upper mantle. The electromagnetic methods permit estimating the temperature and fluid content (and/or carbon (graphite) content) in the Earth’s interior. These estimates are independent of the traditionally used estimates based on seismic methods because the dependence of electrical conductivity on the physical properties of the rock is based on different principles than the behavior of the elastic waves. The region is characterized by a complicated geological structure with intense three-dimensional (3D) surface heterogeneities, which significantly aggravate the retrieval of the information about the deep horizons in the structure of the Earth’s mantle from the observed electromagnetic (EM) fields. The detailed analysis of the nature of the deep electrical conductivity and structural features of the transition from the Pacific to Eurasia included numerical modeling of the typical two- and three-dimensional models has been carried out. Based on this analysis, the approaches that increase the reliability of the interpretation of the results of the EM studies are suggested.  相似文献   

12.
Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.  相似文献   

13.
本文以宁夏区域地震台网为例,分析了波形互相关法在判定重复地震中可能存在的问题并讨论了相应的处理方法,通过构建三维非均匀体模型并利用谱元法数值模拟地震波的传播,统计了不同台站观测到的地震对波形互相关系数的分布,进一步研究了互相关系数与非均匀体性质及震源机制解之间的关系。结果表明:针对宁夏区域地震台网,利用波形互相关法判定重复地震比地震定位方法更有效;互相关系数在不同台站的取值与震源附近三维非均匀体强度和直达波与尾波的振幅比有关,对于相同的震源间距,较弱震源、较弱非均匀体或者较强振幅的直达波均会导致波形互相关系数变高,因此应选取更高的互相关系数阀值来判定重复地震。宁夏区域地震台网平均台间距为30—50 km,通过选取直达波较弱的台站或只截取尾波窗口计算互相关系数并设定较高的阀值,利用波形互相关法可有效地判定ML1.0—3.0重复地震,进而为重复地震的监测与研究提供依据。   相似文献   

14.
Strong ground motion observed at an instrumented hill site is first analysed through the standard (SSR) and the horizontal-to-vertical (HVSR) spectral ratio techniques. A reasonable agreement is found between these approaches. The observations are then compared with 3D numerical simulations, performed with a highly efficient numerical code based on a spectral method, that allowed for reasonable computer times also on a PC. The observed amplification is significantly higher than that computed with a 3D homogeneous model of the mountain, suggesting that local response is governed by large-scale and small-scale soil heterogeneities rather than by topographic site effects. The introduction of a local near-surface inclusion of nonhomogeneous soil material under one of the recording stations has not significantly improved the numerical results. The observed data are also compared with the results of simplified simulations, either using 2D homogeneous models or coupling the 3D response with a 1D local soil profile. The results of such simplified approaches are discussed and their usefulness is emphasised.  相似文献   

15.
The purpose of this work is to present three methods of analysis for elastic waves propagating in two dimensional, elastic nonhomogeneous media. The first step, common to all methods, is a transformation of the governing equations of motion so that derivatives with respect to the material parameters no longer appear in the differential operator. This procedure, however, restricts analysis to a very specific class of nonhomogeneous media, namely those for which Poisson's ratio is equal to 0.25 and the elastic parameters are quadratic functions of position. Subsequently, fundamental solutions are evaluated by: (i) conformal mapping in conjunction with wave decomposition, which in principle allows for both vertical and lateral heterogeneities; (ii) wave decomposition into pseudo-dilatational and pseudo-rotational components, which results in an Euler-type equation for the transformed solution if medium heterogeneity is a function of one coordinate only; and (iii) Fourier transformation followed by a first order differential equation system solution, where the final step involving inverse transformation from the wavenumber domain is accomplished numerically. Finally, in the companion paper numerical examples serve to illustrate the above methodologies and to delineate their range of applicability.  相似文献   

16.
利用大地测量和历史相对海平面变化数据,结合地震剪切波层析模型,联合确定了新的末次冰期冰川均衡调整（GIA）模型,其中地幔黏滞度不仅沿径向而且沿横向变化.研究思路是,先尝试性地选择比例系数β,利用与地震剪切波速异常的线性关系,计算地幔黏滞度横向扰动,并与横向均匀的黏滞度参考模型叠加给出3D地幔黏滞度模型;再利用耦合拉普拉斯方程的有限元法算法进行GIA预测;然后,重复该过程,直到预测与观测之间的吻合满意为止.主要结论有：（1）给出了横向非均匀的地幔黏滞度模型（RF3L20(β=0.4)）,发现了黏滞度显著的横向非均匀性和其对GIA预测的显著影响,指出横向非均匀不完全是由热效应引起的,可能还与化学组分等其他因素有关,该模型可用于地幔动力学研究.（2）给出了全球现今多种GIA预测速率,可为板块运动、陆地水储量、海水质量变化和冰川冰雪质量非平衡监测提供重要的改正.  相似文献   

17.
The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case.  相似文献   

18.
This investigation deals with the torsional balance of the earthquake response and design of elastic asymmetric structures with frictional dampers. Plan asymmetry leads to an uneven lateral deformation demand among structural members and to unbalanced designs with larger capacities in some resisting planes. Frictional dampers are capable of controlling lateral‐torsional coupling by placing the so‐called empirical center of balance (ECB) of the structure at equal distance from all edges of the building. This rule is developed for single‐story systems with linear and inelastic behavior. However, recently obtained theoretical and experimental results demonstrate that this rule carries over to multistory structures. Results show that the peak displacement demand at the building edges and that of resisting planes equidistant from the geometric center may be similar if the damper is optimally placed. It is also shown that torsional amplification of the edge displacements of arbitrary asymmetric structures relative to the displacement of the symmetric counterparts are approximately bound by a factor of 2. Furthermore, frictional dampers are equally effective in controlling lateral‐torsional coupling of torsionally flexible as well as stiff structures. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

19.
A 3D relocation technique permits precise locations of induced earthquakes. Geostatistical processing using the data of 87 boreholes provides the basis of a precise 3D structure, with a dome geometry. Conventional laboratory mechanical tests performed on deep rock samples (1000 m to 5000 m) define the rock properties at depths similar to those of the seismic events (1<M L<4.2) that range from 1 to 7 km.In the studied period, most (85%), of the events were located above the gas reservoir, with very few located in the reservoir itself. Because the production parameters (50 MPa depletion of the gas pressure reservoir) are homogeneous throughout the gas field, the lateral inhomogeneity of the seismic rupture locations are a consequence of variations in the rheological response of the dome to the deformation induced by gas production.Here a ratio of two is found between the elastic modulus of the seismic rock matrix and the elastic modulus of the aseismic rock matrix. The contrast in strength is at least as great, if not greater. Repeated measured surface deformations involve the whole structure. Spatial and temporal deformations indicate that aseismic deformation is quantitatively the main process of this structural deformation. The heterogeneous stress pattern inferred fromP-axes of induced earthquakes disagrees with the tectonic regional stress field. The radial distribution ofP-axes towards the gas reservoir probably reflects the production induced deformation. The inferred deformation of the dome occurs in response to weak induced stresses.  相似文献   

20.
《Journal of Geodynamics》2009,47(3-5):104-117
Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.  相似文献   

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

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