首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 93 毫秒
1.
This paper reviews the theoretical principles of subspace system identification as applied to the problem of estimating black‐box state‐space models of support‐excited structures (e.g., structures exposed to earthquakes). The work distinguishes itself from past studies by providing readers with a powerful geometric interpretation of subspace operations that relates directly to theoretical structural dynamics. To validate the performance of subspace system identification, a series of experiments are conducted on a multistory steel frame structure exposed to moderate seismic ground motions; structural response data is used off‐line to estimate black‐box state‐space models. Ground motions and structural response measurements are used by the subspace system identification method to derive a complete input–output state‐space model of the steel frame system. The modal parameters of the structure are extracted from the estimated input–output state‐space model. With the use of only structural response data, output‐only state‐space models of the system are also estimated by subspace system identification. The paper concludes with a comparison study of the modal parameters extracted from the input–output and output‐only state‐space models in order to quantify the uncertainties present in modal parameters extracted from output‐only models. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

2.
大气边界层湍流的动力非平稳性的验证   总被引：7，自引：1，他引：6       下载免费PDF全文
首次用验证时间序列中是否存在动力非平稳性的一种简单图示方法——space time index法来分析大气边界层湍流的动力平稳性特征.本文以取自淮河流域和威斯康星森林下垫面条件下的三维高精度风速和温度、湿度湍流脉动资料对大气边界层湍流的平稳性特征进行了分析.结果表明space time index方法能有效地检验大气边界层湍流信号中是否存在动力平稳性.另外,均匀下垫面条件（水稻田）及复杂下垫面条件（森林）下的大气边界层湍流信号中几乎都存在动力非平稳性,大气湍流动力学非平稳性可能是边界层湍流信号相当普遍具有的一种特性.大气边界层湍流中的间歇性和相干结构使得其非平稳性图形的特征不同于一般时间序列非平稳性图形的“V”型特征;森林下垫面条件下的湍流信号比相对均匀下垫面（水稻田）下的湍流信号更有组织性,相干结构更强.  相似文献   

3.
The aim of this study is the refinement of the dynamics of a recent (1994?C1999) minor, slow-inflation episode of the Santorini (Thera) volcano, famous for the Minoan (??3600 B.C.) eruption and the identification of the parameters of the magmatic source responsible for the inflation. Based on the Mogi source equations, on geodetic observations of base-line changes, on a topological, grid-search approach and on the reasonable assumption that the magma source remained practically stable in map view during the inflation period, we have been able to refine the location and depth (approximately 2.7 km) of the magma center. A tendency for increase of the magma pressure with time, roughly corresponding to a sphere with radius between 30 and 60 m, and a short deflation interval were also documented. The overall modeling was based on a topological method of inversion in two steps and for a selected 4-D grid. At a first step the system of Mogi-source equations was approximated by the intersection of the 4-D subspaces (defined by sets of grid points) each satisfying one observation equation on the basis of a grid-search procedure. At a second step, the best estimate of the Mogi source solution and its full variance-covariance matrix were defined using a common stochastic approach. The overall approach leads to a solution of a system of equations focusing on a 4-D space bounding significant minima in the misfits between model and observed values, and not on solutions focusing on single points, usually trapped in local minima. This study is important to understand a new phase of volcanic unrest since January 2011, while the proposed methodology, inspired from traditional navigation methods may be useful for other inversion problems leading to redundant systems of highly non-linear equations with n unknowns (i.e. topological solutions in the n-D space).  相似文献   

4.
Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.  相似文献   

5.
During our recent work with 3-D dynamic ray-tracing and velocity inversion problems, a new 3-D model generation system has been developed using a so-called “solid modeling” technique. The term “solid modeling” refers to the fact that the logical system governing the internal geometrical properties of the model describes the model as a combination of “solids” or “volumes” in 3-D space. In each of these volumes the physical parameters (such as seismic velocity, density) vary continuously. Discontinuous changes occur only across the model interfaces separating the volumes. The model is constructed by firstly forming a number of “simple volumes” from the given interfaces and then combining these simple volumes into more complex volumes which represent the physical volumes of the model. It is easy to make changes to the model, by adding volumes or subtracting volumes and perform more composite operations, all by simple use of Boolean expressions. Every time a model has been specified (or changed), the internal logic automatically carries out a check of physical consistency of the 3-D model space (no overlapping volumes, no holes). By including various types of coordinate transformations, different kinds of complex structures can be handled, such as salt domes and vertical and near vertical faulting.  相似文献   

6.
A complex and highly dynamical ocean region, the Agulhas Current System plays an important role in the transfer of energy, nutrients and organic material from the Indian to the Atlantic Ocean. Its dynamics are not only important locally, but affect the global ocean-atmosphere system. In working towards improved ocean reanalysis and forecasting capabilities, it is important that numerical models simulate mesoscale variability accurately—especially given the scarcity of coherent observational platforms in the region. Data assimilation makes use of scarce observations, a dynamical model and their respective error statistics to estimate a new, improved model state that minimises the distance to the observations whilst preserving dynamical consistency. Qualitatively, it is unclear whether this minimisation directly translates to an improved representation of mesoscale dynamics. In this study, the impact of assimilating along-track sea-level anomaly (SLA) data into a regional Hybrid Coordinate Ocean Model (HYCOM) is investigated with regard to the simulation of mesoscale eddy characteristics. We use an eddy-tracking algorithm and compare the derived eddy characteristics of an assimilated (ASSIM) and an unassimilated (FREE) simulation experiment in HYCOM with gridded satellite altimetry-derived SLA data. Using an eddy tracking algorithm, we are able to quantitatively evaluate whether assimilation updates the model state estimate such that simulated mesoscale eddy characteristics are improved. Additionally, the analysis revealed limitations in the dynamical model and the data assimilation scheme, as well as artefacts introduced from the eddy tracking scheme. With some exceptions, ASSIM yields improvements over FREE in eddy density distribution and dynamics. Notably, it was found that FREE significantly underestimates the number of eddies south of Madagascar compared to gridded altimetry, with only slight improvements introduced through assimilation, highlighting the models’ limitation in sustaining mesoscale activity in this region. Interestingly, it was found that the threshold for the maximum eddy propagation velocity in the eddy detection scheme is often exceeded when data assimilation relocates an eddy, causing the algorithm to interpret the discontinuity as eddy genesis, which directly influences the eddy count, lifetime and propagation velocity, and indirectly influences other metrics such as non-linearity. Finally, the analysis allowed us to separate eddy kinetic energy into contributions from detected mesoscale eddies and meandering currents, revealing that the assimilation of SLA has a greater impact on mesoscale eddies than on meandering currents.  相似文献   

7.
历史地震资料与地震的中长期预测   总被引：1，自引：0，他引：1  
故障诊断模型的基本内容是根据动态系统的外部特征来判断系统内部是否发生故障及确定故障发生的部位、时间和大小。由于故障诊断技术在监测及诊断故障的思路上和地震预报有很多相似之处，因此，把故障诊断技术应用于地震预报是可行的。由于该模型与其他数学模型一样，需要有较多的学习过程，所以历史地震资料在该模型中有重要作用。而其中的鲁棒性故障诊断模型在抑制各子模型的个性，凸现其共性方面有其特有的性质，所以把它作为一种综合模型，能抑制各子模型的个性，突出在地震预测方面的共性，从而提高预测的精确性。本文根据一个实例，说明了这种综合性模型的可行性。在文章的最后，由信息量的分析，说明了模型的鲁棒性特征。  相似文献   

8.
本文采用欧拉反褶积、场源参数成像(SPI)、场源边界提取(SED)、莫霍面反演、地壳三维可视化等多源方法,对青藏高原东北缘地区的布格重力场进行反演与分析,深入研究该地区的深部结构与变形特征,探讨区域深部孕震环境及动力学机制.研究表明,青藏高原东北缘的布格重力场整体呈负异常值,具有明显的分区性,表现出鄂尔多斯盆地异常值相对偏高、阿拉善块体次之、青藏高原块体极低的特点,其中海源断裂系形成了一条宽缓的弧形重力梯度条带,梯度值达1.2 mGal·km^-1.欧拉结果显示,鄂尔多斯盆地相比于青藏高原块体而言,场源点具有较强的均一性,场源强度值高(密度值高)且深度稳定在25~32 km范围内,而高原块体的中下地壳尺度广泛分布着低密度异常体.SPI图可知,海源弧形断裂系位于“浅源异常”弧形区,反映其地壳较为活跃,易发生中强地震.SED图揭示青藏高原地壳向东北扩展,经过几大断裂系的调节后运动矢量向东或东南转化,SED与GPS、SKS运动特征大致相同,说明地表-地壳-地幔的运动特征有着较强的一致性.青藏高原东北缘地区壳幔变形是连贯的,加之莫霍面由北向南、由东向西是逐渐加深的,因此属于垂向连贯变形机制,不符合下地壳管道流动力学模式.区域形成了似三联点构造格局,其中海源弧形断裂系的深部地壳结构复杂,高低密度异常体复杂交汇,是青藏高原、阿拉善、鄂尔多斯三大块体相互作用的重要枢纽,其运动学特征总体为中段走滑尾端逆冲,而断裂系正处于大型的弧形莫霍面斜坡带之上,具备强震的深部孕震环境,因此大尺度的运动调节与深部孕震条件共同促使了该地区中强震的多发.  相似文献   

9.
Abstract

Severe unidirectional Fourier truncation of the equations for 2-D incompressible flow leads to a system of three coupled PDEs in one space dimension with the same quadratic invariants as the original set (i.e. energy and enstrophy). Numerically generated equilibria for inviscid, truncated versions of the reduced system are well approximated by Kraichnan's energy-enstrophy equipartition spectra. Viscous calculations for decaying turbulence at moderate resolution (1024 degrees of freedom) also appear to be consistent with a direct, k ?3, enstrophy cascading inertial range when the dissipation is small. Dissipation range intermittency in the form of spatially intermittent enstrophy dissipation with occasional strong bursts producing linear phase locking is also observed. In contrast to full 2-D simulations, no tendency towards the emergence of isolated, coherent vorticity structures is observed. The model consequently mimics some, but not all, of the properties of the full 2-D set.  相似文献   

10.
We present results on evolving geometrical and material properties of large strike-slip fault zones and associated deformation fields, using 3-D numerical simulations in a rheologically-layered model with a seismogenic upper crust governed by a continuum brittle damage framework over a viscoelastic substrate. The damage healing parameters we employ are constrained using results of test models and geophysical observations of healing along active faults. The model simulations exhibit several results that are likely to have general applicability. The fault zones form initially as complex segmented structures and evolve overall with continuing deformation toward contiguous, simpler structures. Along relatively-straight mature segments, the models produce flower structures with depth consisting of a broad damage zone in the top few kilometers of the crust and highly localized damage at depth. The flower structures form during an early evolutionary stage of the fault system (before a total offset of about 0.05 to 0.1 km has accumulated), and persist as continued deformation localizes further along narrow slip zones. The tectonic strain at seismogenic depths is concentrated along the highly damaged cores of the main fault zones, although at shallow depths a small portion of the strain is accommodated over a broader region. This broader domain corresponds to shallow damage (or compliant) zones which have been identified in several seismic and geodetic studies of active faults. The models produce releasing stepovers between fault zone segments that are locations of ongoing interseismic deformation. Material within the fault stepovers remains damaged during the entire earthquake cycle (with significantly reduced rigidity and shear-wave velocity) to depths of 10 to 15 km. These persistent damage zones should be detectable by geophysical imaging studies and could have important implications for earthquake dynamics and seismic hazard.  相似文献   

11.
利用曲线网格有限差分方法,研究了三维倾斜断层的破裂传播过程.基于断层面生成贴体曲线网格,并通过坐标变换将含曲线网格的物理空间转换到含均匀直角网格的计算空间,实现了有限差分方法对复杂界面的处理.通过模拟地震断层的自相似破裂和自发破裂,并与已有发表的结果对比,发现拟合程度较高,验证了本方法的有效性和精确性.重点研究了不同倾角的倾斜断层破裂,最后展望了今后用本方法对非均匀介质中和任意起伏地表下的任意非平面断层破裂动力学的进一步研究.  相似文献   

12.
—Large, extended fault systems such as those in California demonstrate complex space-time seismicity patterns, which include repetitive events, precursory activity and quiescence, and aftershock sequences. Although the characteristics of these patterns can be qualitatively described, a systematic quantitative analysis remains elusive. Our research suggests that a new pattern dynamics methodology can be used to define a unique, finite set of seismicity patterns for a given fault system. In addition, while a long-sought goal of earthquake research has been the reliable forecasting of these events, very little progress has been made in developing a successful, consistent methodology. In this report, we document the discovery of systematic space-time variations in seismicity from southern California using a new technique. Here we present examples of this analysis technique on data obtained prior to events in seismically active areas that show coherent regions associated with the future occurrence of major earthquakes in the same areas. These results strongly support the hypothesis that seismic activity is highly correlated across many space and time scales within large volumes of the earth's crust.  相似文献   

13.
We present a model for earthquake failure at intermediate scales (space: 100 m–100 km, time: 100 m/v shear- 1000’s of years). The model consists of a segmented strike–slip fault embedded in a 3-D elastic solid as in the framework of Ben-Zion and Rice (1993). The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of the regions around the fault, static/kinetic friction laws with possible gradual healing, and stress transfer based on the solution of Chinnery (1963) for static dislocations in an elastic half-space. As a new ingredient, we approximate the dynamic rupture on a continuous time scale using a finite stress propagation velocity (quasi–dynamic model) instead of instantaneous stress transfer (quasi–static model). We compare the quasi–dynamic model with the quasi–static version and its mean field approximation, and discuss the conditions for the occurrence of frequency-size statistics of the Gutenberg–Richter type, the characteristic earthquake type, and the possibility of a spontaneous mode switching from one distribution to the other. We find that the ability of the system to undergo a spontaneous mode switching depends on the range of stress transfer interaction, the cell size, and the level of strength heterogeneities. We also introduce time-dependent log (t) healing and show that the results can be interpreted in the phase diagram framework. To have a flexible computational environment, we have implemented the model in a modular C++ class library.  相似文献   

14.
The capacity spectrum method (CSM), capable of predicting the demands of forces and deformations of the inelastic system, has been applied in the ATC and FEMA guidelines. The deformation of an inelastic system is solved iteratively by using the equivalent linearization for CSM, which actually forms a nonlinear map or discrete dynamical system. However, the iterative procedure of CSM did not converge for some inelastic systems, and the complicated dynamical phenomena for the solutions such as the periodic oscillation, period-doubling bifurcation and chaos may occur, which were shown in the bifurcation plots of iterative map of the simplified CSM in ATC40 and FEMA440. This paper presents a novel method to analyze and control the non-convergence of the iterative procedure of CSM from the perspective of chaotic dynamics. The Lyapunov exponent of the dynamical system is employed to identify the evolutional state and stability of solutions. Finally, the stability transformation method as a simple, versatile and effective chaos feedback control approach is applied to control the convergent failure of CSM in ATC40 and FEMA440. The numerical results illustrate that the stability transformation method can capture the desired fixed points of the dynamical system and obtain the stable convergent solutions of CSM.  相似文献   

15.
— Using 3-D numerical modeling of seismic wave propagation we investigate the possibility of generating fault zone (FZ) trapped wave energy from sources well outside a fault. The FZ is represented by a O(200 m) wide vertical low velocity layer in a half space. We find that FZ trapped waves can be excited from sources well outside the fault if (1) the low-velocity structure extemds only to shallow depth and the source is located at greater depth or (2) the structure of the low-velocity zone is such that only the shallow part of the FZ traps energy. FZ trapped waves are not excited from sources well outside a FZ continuous with depth. The results support, in conjunction with recent observational evidence, a model for natural faults with shallow trapping structures rather than ones that span the entire seismogenic zone. This may have implications for fault mechanics as well as for aspects of shaking hazard near faults.  相似文献   

16.
  相似文献   

17.
  相似文献   

18.
While rainfall intermittency is a dynamical phenomenon, little progress has been made in the literature on the link between rainfall intermittency and atmospheric dynamics. We present the basic dynamical models of intermittency that are phenomenologically most similar to rainfall: Pomeau–Manneville Type-III and On–Off. We then illustrate each type with both a 1-D iterative map and a corresponding stochastic process stressing the appearance of these dynamics in high-dimensional (stochastic) systems as opposed to low-dimensional chaotic systems. We show that the pdf of rainfall intensities, the pdf of “laminar phases” (periods of zero rainfall intensity), and the spectrum of the rainfall series all have power-law behavior that is broadly consistent with intermittency in the classic types. Using a seasonal analysis, we find that summer convective rainfall at daily and sub-daily scales seems consistent with features of Type-III intermittency. The correspondence with Type-III intermittency and a preliminary entropic analysis further suggest that rainfall may be an example of sporadic randomness, blending deterministic and stochastic components.  相似文献   

19.
This paper describes a methodology, based on dynamical systems theory, to model and predict streamflow at the daily scale. The model is constructed by developing a multidimensional phase-space map from observed streamflow signals. Predictions are made by examining trajectories on the reconstructed phase space. Prediction accuracy is used as a diagnostic tool to characterize the nature, which ranges from low-order deterministic to stochastic, of streamflow signals. To demonstrate the utility of this diagnostic tool, the proposed method is first applied to a time series with known characteristics. The paper shows that the proposed phase-space model can be used to make a tentative distinction between a noisy signal and a deterministic chaotic signal.The proposed phase-space model is then applied to daily streamflow records for 28 selected stations from the Continental United States covering basin areas between 31 and 35 079 km2. Based on the analyses of these 28 streamflow time series and 13 artificially generated signals with known characteristics, no direct relationship between the nature of underling stream flow characteristics and basin area has been found. In addition, there does not appear to be any physical threshold (in terms of basin area, average flow rate and yield) that controls the change in streamflow dynamics at the daily scale. These results suggest that the daily streamflow signals span a wide dynamical range between deterministic chaos and periodic signal contaminated with additive noise.  相似文献   

20.
In this article, we first reviewed the method of boundary integral equation (BIEM) for modelling rupture dynamics of a planar fault embedded in a 3-D elastic half space developed recently (ZHANG and CHEN, 2005a,b). By incorporating the half-space Green's function, we successfully extended the BIEM, which is a powerful tool to study earthquake rupture dynamics on complicated fault systems but limited to full-space model to date, to half-space model. In order to effectively compute the singular integrals in the kernels of the fundamental boundary integral equation, we proposed a regularization procedure consisting of the generalized Apsel-Luco correction and the Karami-Derakhshan algorithm to remove all the singularities, and developed an adaptive integration scheme to efficiently deal with those nonsingular while slowly convergent integrals. The new BIEM provides a powerful tool for investigating the physics of earthquake dynamics. We then applied the new BIEM to investigate the influences of geometrical and physical parameters, such as the dip angle (δ) and depth (h) of the fault, radius of the nucleation region (Rasp), slip-weakening distance (Dc), and stress inside (Ti) and outside (Te) the nucleation region, on the dynamic rupture processes on the fault embedded in a 3-D half space, and found that (1) overall pattern of the rupture depends on whether the fault runs up to the free surface or not, especially for strike-slip, (2) although final slip distribution is influenced by the dip angle of the fault, the dip angle plays a less important role in the major feature of the rupture progress, (3) different value of h, δ, Rasp, Te, Ti and Dc may influence the balance of energy and thus the acceleration time of the rupture, but the final rupture speed is not controlled by these parameters.  相似文献   

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

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