首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 131 毫秒
1.
Sharp boundary inversion of 2D magnetotelluric data   总被引:6,自引:0,他引:6  
We consider 2D earth models consisting of laterally variable layers. Boundaries between layers are described by their depths at a set of nodes and interpolated laterally between nodes. Conductivity within each layer is described by values at a set of nodes fixed within each layer, and is interpolated laterally within each layer. Within the set of possible models of this sort, we iteratively invert magnetotelluric data for models minimizing the lateral roughness of the layer boundaries, and the lateral roughness of conductivities within layers, for a given level of data misfit. This stabilizes the inverse problem and avoids superfluous detail. This approach allows the determination of boundary positions between geological units with sharp discontinuities in properties across boundaries, while sharing the stability features of recent smooth conductivity distribution inversions.
We compare sharp boundary inversion results with smooth conductivity distribution inversion results on a numerical example, and on inversion of field data from the Columbia River flood basalts of Washington State. In the synthetic example, where true positions and resistivities are known, sharp boundary inversion results determine both layer boundary locations and layer resistivities accurately. In inversion of Columbia flood basalt data, sharp boundary inversion recovers a model with substantially less internal variation within units, and less ambiguity in both the depth to base of the basalts and depth to resistive basement.  相似文献   

2.
Abstract

The separation of sidewall boundary layers in a rotating annulus of homogeneous fluid is studied experimentally. The flow is driven by a differentially rotating lid, and a very small uniform slope of the bottom causes a weak mountain vortex pair to form in the interior, away from the sidewalls. A necessary condition for aerodynamic separation of the sidewall boundary layers is derived and compared with the experimental results. The laboratory flow separates for parameters that are just slightly more inviscid than those required by the necessary condition for the existence of adverse pressure gradients at the wall. As the bottom friction is decreased further, the flow becomes unsteady and chaotic. The most interesting aspect of this problem is that chaotic interior behavior, associated with the separated boundary layer, is observed for parameter values for which the interior topographically forced flow is, by itself, essentially linear.  相似文献   

3.
In this study, we examine the development of topography on a thin dense layer at the base of the lower mantle. The effect of the convecting mantle above is represented as a traction acting on the upper surface of the layer. Topography on the layer boundaries is predicted by a balance of dynamic flow stress and external traction. The nature of boundary topography depends on the magnitude of the driving tractions and the density variation within the layer. If we assume that the layer density is greatest beneath areas of mantle downwelling and decreases to a minimum beneath areas of mantle upwelling (the layer is thermally coupled to the convection in the overlying mantle) then its upper boundary develops a cusp-like peak beneath the upwelling mantle. The height of this peak is potentially much greater than the layer thickness. If, however, the layers are effectively coupled by viscous shear then internal density gradients of the opposite sign may be established. In this case, we observe solutions where the layer is completely swept away beneath areas of mantle downwelling leaving steep-sided ‘islands’ of dense material. This mechanism therefore provides a possible explanation for steep-sided anomalously slow regions at the base of the mantle observed by seismic methods (e.g. beneath south Africa) or for discrete ultralow velocity zones detected at the core-mantle boundary beneath locations of surface hotspots. The magnitude of the upper boundary driving tractions compared to the density gradient within the layer is the key parameter that determines the nature of flow in, and consequently boundary topography of, the layer. The deflection of the core-mantle boundary is small compared with that of the top of the dense layer, but a change in sign of the ratio of these deflections is observed as the magnitude of the driving tractions changes relative to the magnitude of the internal density gradient. We compare seismic measurements of core-mantle boundary topography and D′′ topography with the predictions of this model in an attempt to constrain model parameters, but no clear correlation seems to exist between D′′ thickness and CMB topography.  相似文献   

4.
中国构造应力场与大震复发周期关系的数值模拟   总被引:6,自引:0,他引:6       下载免费PDF全文
由古地震资料推算得到的中国境内各活动断裂带上大震复发周期值存在很大差异,对于这种差异可从构造应力产生的弹性应变能分布去探讨其原因。经过对中国及邻区边界条件和受力方式进行较详细的分析后,先以最大主压应力方向为判据,用有限元数值模拟方法反演中国及邻区的板块边界力和区内的构造应力场,然后分析计算应变能随离板块边界的距离增大所呈现的衰减规律,由此换算出要达到相同应变能(一次地震)各地所需的积累时间,并将其与用古地震方法得到的一些已知断层带大震复发周期进行对比。结果表明用应力衰减图象能解释不同地区大震复发周期的显著差异  相似文献   

5.
SUMMARY

Basing their studies on an extended Reynolds principle of similarity, the authors present in this paper a thorough theoretical investigation of the turbulent flow in pipes.

Making appropriate assumptions regarding the distribution of the turbulent apparent viscosity over the entire pipe diameter, the well-known discharge formulas, i. e. the resistance laws, of Prandtl-V. Kármán for the hydraulically smooth and the hydraulically rough region, and of Colebrook for the transition region—the accuracy of these being well established by measurements—can be clearly shown to be correct.

Thus it is also possible to determine the velocity distribution in the region near the wall for the hydraulically smooth region and the transition region, as well as for the hydraulically rough region.

The concept, that the flow with a fully developed roughness effect begins when the laminar boundary layer at the wall disappears can be seen to be in close agreement with the asssumptions made. The limits which were established for the hydraulically rough region and the transition region correspond well with the curve plotted by PRANDTL and the boundary curve of Rouse.  相似文献   

6.
This paper explores the concept of a macroscopic boundary between turbulent flows above and within rough permeable walls. The macroscopic boundary and the associated conditions for macroscopic flow variables have been thoroughly investigated for laminar, but not for turbulent flows. The literature on laminar flows follows two main conceptual models of the boundary: sharp boundary with step changes in macroscopic variables and gradual boundary with smooth changes of variables. The former approach is usually associated with the two-domain simulation models and the latter one with the single-domain models. This paper presents the derivation of the step conditions for velocity and shear stress at the macroscopic boundary between turbulent boundary layer and turbulent porous media flows. The physical meaning of the main terms in the shear stress condition is discussed in order to clarify the relationship between two-domain and single-domain simulation models.  相似文献   

7.
以实际震例的异常数据为样本,利用粗糙集理论的离散、约简和推理功能,对地震前兆异常属性进行约简,进而对地震震级进行推理分类.初步结果表明,该方法可以在一定程度上找到异常数据与震级之间的内在联系,是一种可行的地震前兆综合预测方法.  相似文献   

8.
青藏高原岩石圈多层构造应力场   总被引:3,自引:0,他引:3  
王绳祖 《地震》2002,22(3):21-26
青藏高原构造应力场可按岩石圈下层、多震层和浅层地壳区分为三层。除了震源机制解方法和井孔原地测量方法可分别用于推测多震层和浅层的应力状况外,还可根据下层塑性流动网络,采用平分网络共轭角的方法估计下层的应力方向。对比岩石圈下层与上层(多震层)的构造应力场,其结果表明: 由于板块边缘驱动力主要通过下层的网络状流动实现其远程传递,故在总体作用趋势上,上层的应力方向受控于下层;又由于高原靠近喜马拉雅驱动边界,部分驱动力直接沿上层传递,致使局部地区上、下层应力方向相差显著。  相似文献   

9.
A two-dimensional numerical model is applied to a coastal ocean wherein alongshore elevation and density gradients, normally calculated by a three-dimensional model, are instead supplied by climatologically averaged data for the California Current System between 25 and 40°N. Surface wind stress is also obtained from climatological data. Both surface and bottom boundary layers are resolved in the model calculations; a second moment turbulence closure submodel supplies vertical diffusivities. Near steady state solutions are possible when surface buoyancy flux is imposed at the surface.Model results are as follows: Southward wind stress produces a broad equatorward current with an embedded coastal jet in accordance with previous studies. Positive wind stress curl reduces the jet current and produces a poleward undercurrent which then surfaces as the curl is increased. The jet currents are reduced and poleward flow increases as bottom steepness increases; to a lesser extent, inclusion of the beta effect has a similar effect. The existence of near bottom, poleward or equatorward flow is explained rather simply in terms of the bottom stress resulting from the alongshore balance of surface wind stress and vertically integrated pressure gradient, the latter involving the alongshore surface elevation and density gradient. A further finding is that the upwelling circulation associated with wind stress is confined to the top 200 to 300 m of the ocean along the California coast.  相似文献   

10.
The determination of a local lithospheric stress pattern from low magnitude seismic activity recorded at a loca array is attempted by means of numerical simulations that include waveform inversion and annealing inversion techniques. A set of flat stratified layers are assumed for the propagating medium. The tests show that the annealing algorithm is a useful tool in this kind of study where collected seismic activity at local arrays could help us better understand local stress distributions in the lithosphere. Although real data recorded at local seismic arrays should be processed, assuming more realistic crust and upper mantle models, many characteristics of the annealing process described here could be used when studying real problems.  相似文献   

11.
The nonhydrostatic pressure effects on the generation and propagation of wind-forced internal waves are studied with a two-dimensional numerical ocean model. A one-way directed wind pulse over a stratified ocean initiates surface and internal waves in a closed basin. The studies are performed with horizontal grid sizes in the range from 1 km to 62.5 m. The experiments are performed with both a hydrostatic and a nonhydrostatic model, facilitating systematic studies of the sensitivity of the numerical model results to the grid size and to the nonhydrostatic pressure adjustments. The results show that the nonhydrostatic pressure effects are highly dependent on the grid size and grow with increased resolution. In the internal depression wave, the horizontal nonhydrostatic pressure gradients reach the same order of magnitude as the hydrostatic gradients in the high-resolution nonhydrostatic studies. In these studies, the nonhydrostatic pressure gradients approximately balance the corresponding hydrostatic pressure gradients in the internal depression wave, and the wave degenerates into a train of soliton waves. The time for the soliton form to develop agrees with the steepening timescale calculated from Korteweg-de Vries theory. In the high-resolution hydrostatic model, the internal depression wave takes the form of a single wave front. When the internal waves are generated in the boundary layers, the nonhydrostatic pressure gradients are much smaller than the hydrostatic gradients and the generation processes are not effected by the nonhydrostatic pressure with the present range of grid sizes.  相似文献   

12.
This paper presents a closed-form wave function analytic solution of two-dimensional scattering and diffraction of anti-plane SH-waves by a two-dimensional foundationless structure that corresponds to a shear wall on an elastic half-space. A wave-function expansion method is used to solve this model by first prescribing a set of wave functions with undetermined coefficients and then assembling them together based on the stress and displacement boundary conditions on the surface between the structure and half space. This results in a set of infinite equations to be solved by truncating to a finite set. The amplitudes and residuals of the displacement and stress distributions around the structure and nearby ground surface will be discussed carefully. While the solution is analytical, the computation of the numerical results involves the evaluation of complicated integrals. This analytic solution will be helpful to the understanding of propagation of seismic or other stress waves within the superstructure(s) undergoing earthquakes or other blast loads.  相似文献   

13.
The growing need for residential housing in Latin American countries has led to the construction of reinforced concrete buildings with wall thicknesses as low as 8–10 cm. Such walls have typically only a single layer of vertical rebars and are therefore particularly susceptible to out-of-plane failure. In order to investigate the response of the corresponding wall boundary elements, twelve reinforced concrete members with a single layer of vertical rebars were tested under tension–compression cycles. The objective was to gain insight into the parameters governing wall instability and out-of-plane failure, namely the thickness, reinforcement ratio, and eccentricity of the longitudinal rebars with respect to the member axis. This paper summarises the results of the test program, where the specimens' response is analysed also at the global and local levels. The results show that the crack pattern has an important influence on the out-of-plane behaviour and the conditions leading to out-of-plane failure are described. Furthermore, the differences between members with a single layer of vertical rebars and members with two layers are discussed. The influence of the parameters considered in the experimental program is addressed, showing that sections with small thickness and large reinforcement content are more prone to out-of-plane failures. Finally, predictions given by existing models are compared to the new experimental data. The entire data set is publicly available.  相似文献   

14.
Ghosh  Mithun  Das  Ratul 《Water Resources》2021,48(4):544-556
Water Resources - Present study aims to illustrate the role of form-induced shear stress on turbulent kinetic energy (TKE) budget distributions in permeable gravel bed streams. In order to achieve...  相似文献   

15.
采用多触点预压式元件的钻孔变形法地应力测量,能够测得孔壁上较多点的折算位移值,因而能够用平差术减少偶然误差,得出地应力大小和方向的最佳值。本文讨论了这种元件在绝对应力测量和相对应力测量过程中的力学作用,导出了孔壁接触点折算位移向量与元件输出向量之间的基本关系式,以及它们增量之间的关系式。本文是文献[2]工作的推广。当接触点相互之间的距离变大、接触压力的相互影响减弱时,从本文基本关系式得出的极限情形与文献[2]的一致。本文比文献[2]进一步明确了目前实际工作中所用公式的简化前提。文末指出了用多触点元件精确测定地应力尚待解决的一些问题。  相似文献   

16.
裂陷盆地基底双界面模式二维重力反演   总被引:2,自引:1,他引:1       下载免费PDF全文
裂陷盆地基底的起伏表现为非光滑的几何形态,传统的重力反演结果并不能很好地反映这种特点.此外,大多数情况下,重力观测面并不位于盆地上界面,应为单独的起伏观测面,盆地应为上界面和基底组成的双界面模式.基于此,本文研究了起伏观测面上裂陷盆地基底双界面模式二维重力反演方法.研究中假设沉积盆地的沉积层与基底的密度差随深度按双曲线规律变化.将沉积盆地的沉积层剖分成相邻的垂直柱体,其水平尺寸是已知的,顶面与沉积层上界面重合,底面深度代表基底的深度,即为要反演的参数.反演中引入全变差函数作为盆地模型的约束,使得反演结果呈现非光滑形态,符合裂陷盆地基底特征.为减小反演多解性,引入已知深度点作为约束.建立由重力数据拟合、已知深度约束及全变差函数组成的目标函数,采用非线性共轭梯度算法使目标函数最小化.模型试算结果表明该方法可反演裂陷盆地基底起伏,并通过调整正则化参数的值可反演坳陷盆地基底起伏.将该反演方法用于珠江口盆地惠州凹陷和运城-临汾裂陷盆地实际资料处理,其结果较好地反映了裂陷盆地基底起伏特征,为研究盆地构造、油气勘探等提供重要参考.  相似文献   

17.
依据最新的辽宁及其邻近地区地壳上地幔三维地震波速结构,并考虑地质构造、断裂以及地壳深部的构造环境等因素,建立研究区的三维地质模型;结合震源机制解与原地应力测量资料,分析区域现代地应力场特征,确定了研究区的边界约束和动力边界条件.在以上基础上,采用线弹性介质模型进行研究区三维构造应力场的有限元数值模拟,探讨研究区地壳应力场的基本特征.对比分析表明,本区地壳内应力场的三维有限元数值模拟结果与已发生的地震有着较好的对应关系,并指出了值得注意的有潜在发震可能的地震危险区.   相似文献   

18.
IntroductionTectonic earthquakes are the consequence of localized stress concentration. A basic way of learning the process of earthquake preparation and occurrence is to study the distribution and evolution of the tectonic stress.Since the 1980s a series of achievements in Earth stress research have been announced in the world (e.g., Zoback, 1992). In last two decades many researches on crustal tectonic stress have been carried on in China, including earthquake focal mechanism analyses, in-…  相似文献   

19.
针对川滇地区,在大量搜集、总结前人资料和结果的基础上,建立了包括上地壳、下地壳(两层)、上地幔等4层结构的三维有限元模型. 通过反复调试,确定了其边界条件,并且计算了川滇地区背景应力场、断层蠕动产生的应力场和强震触发的应力场,以及它们的动态变化. 结果表明:后续地震大多发生在前面地震引发的库仑破裂应力正值区,前一地震对后续地震有一定的触发作用,强震是在较高的应力背景下成组发生的. 这对判断、识别地震异常,提高地震预测水平有着很大的意义.   相似文献   

20.
This paper is concerned with testing the validity of the ground motions estimated by combining a boundary integral equation method to simulate dynamic rupture along finite faults with a finite difference method to compute the subsequent wave propagation. The validation exercise is conducted by comparing the calculated ground motions at about 100 hypothetical stations surrounding the pure strike-slip and pure reverse faults with those estimated by recent ground motion estimation equations derived by regression analysis of observed strong-motion data. The validity of the ground motions with respect to their amplitude, frequency content and duration is examined. It is found that the numerical simulation method adopted leads to ground motions that are mainly compatible with the magnitude and distance dependence modelled by empirical equations but that the choice of a low stress drop leads to ground motions that are smaller than generally observed. In addition, the scatter in the simulated ground motions, for which a laterally homogeneous crust and standard rock site were used, is of the same order as the scatter in observed motions therefore, close to the fault, variations in source propagation likely contribute a significant proportion of the scatter in observed motions in comparison with travel-path and site effects.  相似文献   

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

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