首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
There are several important wavenumber sampling issues associated with 2.5D seismic modelling in the frequency domain, which need careful attention if accurate results are to be obtained. At certain critical wavenumbers there exist rapid disruptions in the mainly smooth oscillatory spectra. The amplitudes of these disruptions can be very large, and this affects the accuracy of the inverse Fourier transformed frequency-space domain solution. In anisotropic elastic media there are critical wavenumbers associated with each wave mode—the quasi-P (qP) wave, and the two quasi-shear (qS1 and qS2) waves. A small wavenumber sampling interval is desirable in order to capture the highly oscillatory nature of the wavenumber spectrum, especially at increasing distance from the source. Obviously a small wavenumber sampling interval adds greatly to the computational effort because a 2D problem must be solved for every wavenumber and every frequency. The discretisation should be carried out up to some maximum wavenumber, beyond which the field becomes evanescent (exponentially decaying or diffusive). For receivers close to the source, activity persists beyond the critical wavenumber associated with the minimum shear wave velocity in the model. Fortunately, for receivers well removed from the source, the contribution from the evanescent energy is negligible and so there is no need to sample beyond this critical wavenumber. Sampling at Gauss–Legendre spacings is a satisfactory approach for acoustic media, but it is not practical in elastic media due to the difficulty of partitioning the integration around the different critical wavenumbers. We found to our surprise that in transversely isotropic media, the critical wavenumbers are independent of wave direction, but always occur at those wavenumbers corresponding to the maximum phase velocities of the three wave modes (qP, qS1 and qS2), which depend only on the elastic constants and the density. Additionally, we have observed that intermediate layers between source and receiver can filter out to a large degree, the sharp irregularities around the critical wavenumbers in the ω–k y spectra. We have found that, using the spectral element method, the singularities (poles) at the critical wavenumbers which exist with analytic solutions, do not arise. However, the troublesome spike-like behaviour still occurs and can be damped out without distorting the spectrum elsewhere, through the introduction of slight attenuation.  相似文献   

2.
We look at the large-scale dynamo properties of spatially periodic, time dependent, helical 2D flows of the form u(x, t)?=?(? y ?ψ?(x, y, t), ?? x ?ψ?(x, y, t), ?ψ (x, y, t). These flows act as kinematic fast dynamos and are able to generate a mean magnetic field uniform and constant in the xy-plane but whose direction varies periodically along z with wavenumber k. Using Mean Field Electrodynamics, the generation mechanism can be understood in terms of a k-dependent α-effect, which depends on the magnetic Reynolds number, R m . We calculate this effect for different motions and investigate how its limit as k?→?0 depends on R m and on the properties of the flows such as their spatial structure or correlation time. This work generalises earlier studies based on 2D steady flows to motions with time dependence.  相似文献   

3.
An optimization method is used to select the wavenumbers k for the inverse Fourier transform in 2.5D electrical modelling. The model tests show that with the wavenumbers k selected in this way the inverse Fourier transform performs with satisfactory accuracy.  相似文献   

4.
Abstract

We derive an equation governing the nonlinear propagation of a linearly polarized Alfvén wave in a two-dimensional, anisotropic, slightly compressible, highly magnetized, viscous plasma, where nonlinearities arise from the interaction of the Alfvén wave with fast and slow magnetoacoustic waves. The phase mixing of such a wave has been suggested as a mechanism for heating the outer solar atmosphere (Heyvaerts and Priest, 1983).

We find that cubic wave damping dominates shear linear dissipation whenever the Alfvén wave velocity amplitude δvy exceeds a few times ten metres per second. In the nonlinear regime, phase-mixed waves are marginally stable, while non-phase-mixed waves of wavenumber ka are damped over a timescale kuRe 0|δ vy/vA |?2, Re 0 being the Reynolds number corresponding to the Braginskij viscosity coefficient η0 and vA the Alfvén speed. Dissipation is most effective where β = (vs /vA) 2 ≈ 1, vs being the speed of sound.  相似文献   

5.
This study aims to realistically simulate the seismic responses of typical highway bridges in California with considerations of soil–structure interaction effects. The p‐y modeling approaches are developed and validated for embankments and pile foundations of bridges. The p‐y approach models the lateral and vertical foundation flexibility with distributed p‐y springs and associated t‐z and q‐z springs. Building upon the existing p‐y models for pile foundations, the study develops the nonlinear p‐y springs for embankments based on nonlinear 2D and 3D continuum finite element analysis under passive loading condition along both longitudinal and transverse directions. Closed‐form expressions are developed for two key parameters, the ultimate resistant force pult and the displacement y50, where 0.5pult is reached, of embankment p‐y models as functions of abutment geometry (wall width and height, embankment fill height, etc.) and soil material properties (wall‐soil friction angle, soil friction angle, and cohesion). In order to account for the kinematic and site responses, depth‐varying ground motions are derived and applied at the free‐end of p‐y springs, which reflects the amplified embankment crest motion. The modeling approach is applied to simulate the seismic responses of the Painter Street Bridge and validated through comparisons with the recorded responses during the 1992 Petrolia earthquake. It is demonstrated that the flexibility and motion amplification at end abutments are the most crucial modeling aspects. The developed p‐y models and the modeling approach can effectively predict the seismic responses of highway bridges. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

6.
This paper studies the response of pipelines to vibrations induced by the operation of a pavement breaker during the rehabilitation of concrete pavements. An efficient two-and-a-half-dimensional (2.5D) formulation is employed, where the geometry of the structure and the soil is assumed to be invariant in the longitudinal direction, allowing for a Fourier transform of the longitudinal coordinate y along the structure to the wavenumber ky. The dynamic soil–structure interaction problem is solved by means of a 2.5D coupled finite element–boundary element (FE–BE) method using a subdomain formulation. The numerical model is verified by means of results available in the literature for a buried pipeline subjected to incident P- and SV-waves with an arbitrary angle of incidence. The presented methodology is capable to incorporate any type of incident wave field induced by earthquakes, construction activities, traffic, explosions or industrial activities. The risk of damage to a high pressure steel natural gas pipeline and a concrete sewer pipe due to the operation of a pavement breaker is assessed by means of the 2.5D coupled FE–BE methodology. It is observed that the stresses in the steel pipeline due to the operation of the pavement breaker are much lower than those induced by the operating internal pressure. The steel pipeline behaves in the linear elastic range under the combined effect of the loadings, indicating that damage to steel pipelines close to the road due to the operation of a pavement breaker is unlikely. The maximum principal stress in the concrete pipe, on the other hand, remains only slightly lower than the specified tensile strength. The decision to use a pavement breaker should hence be taken with care, as its operation may induce tensile stresses in concrete sewer pipes which are of the same order of magnitude as the tensile strength of the concrete. Assessing the risk of damage by means of vibration guidelines based on the peak particle velocity (PPV) gives, for the particular cases considered, qualitatively similar results.  相似文献   

7.
A control strategy is proposed for variable damping elements (VDEs) used together with auxiliary stiffness elements (ASEs) that compose a time‐varying non‐linear Maxwell (NMW) element, considering near‐future excitation influence. The strategy first composes a state equation for the structural dynamics and the mechanical balance in the NMW elements. Next, it establishes a cost function for estimating future responses by the weighted quadratic norms of the state vector, the controlled force and the VDEs' damping coefficients. Then, the Euler equations for the optimum values are introduced, and also approximated by the first‐order terms under the autoregressive (AR) model of excitation information. Thus, at each moment tk, the strategy conducts the following steps: (1) identify the obtained seismic excitation information to an AR model, and convert it to a state equation; and (2) determine VDEs' damping coefficients under the initial conditions at tk and the final state at tk+L, using the first‐order approximation of the Euler equations. The control effects are examined by numerical experiments. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

8.
Fractal and chaotic phenomena in engineering structure are discussed in this paper, it means that the characters of fractal and chaos on dynamic system of seismic dissipated energy activity intensity E d and activity intensity of seismic dissipated energy moment I e are analyzed carefully. Based on the conceptions of the energy system of engineering structures Θ, seismic dissipated energy activity intensity E d and activity intensity of seismic dissipated energy moment I e, the chaotic phenomena of dynamic systems E d and I e are discovered by theoretic derivation, then the fractal characters of them are also discovered from theoretical inferring and numerical computation. Attractor of relative dimension d 2, Renyi entropy of the second order k 2, mean predictable time scale 1/k 2 and other parameters of the dynamic system which were constructed in light of a large number of actual measuring seismic data have been achieved in the end. These parameters are exactly what the fractal and chaotic phenomena have represented in practical dynamic system, which may be valuable for earthquake-resistant theory and analytical method in practical engineering structures.  相似文献   

9.
We develop a data set of aftershock recordings of the 1999, M = 7.4 Izmit and M = 7.2 Duzce (Turkey) earthquakes to study their source parameters. We combined seismograms from 44 stations maintained by several sources (organizations) to obtain a unified data set of events (2.1 ≤ Mw ≤ 5.5). We calculate source parameters of these small earthquakes by two methods that use different techniques to address the difficulty in obtaining source spectra for small earthquakes subject to interference from site response. One method (program NetMoment (NM), Hutchings, 2004) uses spectra of direct S waves in a simultaneous inversion of local high-frequency network data to estimate seismic moment, source corner frequency (fc), site attenuation (k) and whole-path Q. This approach takes advantage of the source commonality in all recordings for a particular earthquake by fitting a common Brune source spectrum to the data with a and individual k. The second approach (Mayeda et al., 2003) uses the coda method (CM) to obtain “nonmodel-based” source spectra and moment estimates from selected broadband recording sites. We found that both methods do well for events that allow the comparison with seismic moment estimates derived from waveform modeling. Also, source spectra obtained from the two methods are very closely matched for most of the events they have in common. We use an F test to examine the trade-off between k and fc picks identified by the direct S-wave method. About half of the events could be constrained to have less than a 50% average uncertainty in fc and k. We used these source spectra solutions to calculate energy and apparent stress and compare these to estimates from the selected “good quality” source spectra from CM. Both studies have values mutually consistent and show a similar increase in apparent stress with increasing moment. This result has added merit due to the independent approaches to calculate apparent stress. We conclude that both methods are at least partially validated by our study, and they both have usefulness for different circumstances of recording local small earthquakes. CM would work well in studies for which there is a broad magnitude range of events and NM works well for local events recorded by band-limited recorders.  相似文献   

10.
Seismic imaging is an important step for imaging the subsurface structures of the Earth. One of the attractive domains for seismic imaging is explicit frequency–space (f – x) prestack depth migration. So far, this domain focused on migrating seismic data in acoustic media, but very little work assumed visco‐acoustic media. In reality, seismic exploration data amplitudes suffer from attenuation. To tackle the problem of attenuation, new operators are required, which compensates for it. We propose the weighted L 1 ‐error minimisation technique to design visco‐acoustic f – x wavefield extrapolators. The L 1 ‐error wavenumber responses provide superior extrapolator designs as compared with the previously designed equiripple L 4 ‐norm and L∞‐norm extrapolation wavenumber responses. To verify the new compensating designs, prestack depth migration is performed on the challenging Marmousi model dataset. A reference migrated section is obtained using non‐compensating f – x extrapolators on an acoustic dataset. Then, both compensating and non‐compensating extrapolators are applied to a visco‐acoustic dataset, and both migrated sections are then compared. The final images show that the proposed weighted L 1 ‐error method enhances the resolution and results in practically stable images.  相似文献   

11.
The classical finite‐difference methods for seismic wave modelling are very accurate at low wavenumbers but suffer from inaccuracies at high wavenumbers, particularly at Nyquist wavenumber. In contrast, the optimisation finite‐difference methods reduce inaccuracies at high wavenumbers but suffer from inaccuracies at low wavenumbers, particularly at zero wavenumber when the operator length is not long and the whole range of wavenumbers is considered. Inaccuracy at zero wavenumber means that the optimisation methods only have a zeroth‐order accuracy of truncation and thus are not rigorously convergent. To guarantee the rigorous convergence of the optimisation methods, we have developed accuracy‐constrained optimisation methods. Different‐order accuracy‐constrained optimisation methods are presented. These methods not only guarantee the rigorous convergence but also reduce inaccuracies at low wavenumbers. Accuracy‐constrained optimisation methods are applied to staggered‐grid elastic wave modelling.  相似文献   

12.
A control algorithm has been developed for controlling Active Variable Stiffness (AVS) structures. This algorithm analyses information of an observed seismic excitation, estimates the future structural responses and determines how to alter the structure stiffness. An objective structure is assumed to possess N on-off elements whose states are controlled by the proposed algorithm. That is, at a given moment tk, (1) seismic excitation information is expressed by an Auto Regressive (AR) model as the identification model; (2) future excitation information is predicted using the AR model; (3) future responses due to predicted excitation are estimated; (4) based on the initial condition at tk, the responses of 2N possible structural states from tk, to tk+L are calculated; (5) the state which minimizes the input energy during tL is selected; and (6) immediately after tk, on-off elements are set up and subjected to the selected states. The effectiveness of the induced algorithm is confirmed by numerical experiments on a model of a three-storey building under sine and seismic excitations.  相似文献   

13.
The groundwater interbasin flow, Qy, from the north of Yucca Flat into Yucca Flat simulated using the Death Valley Regional Flow System (DVRFS) model greatly exceeds assessments obtained using other approaches. This study aimed to understand the reasons for the overestimation and to examine whether the Qy estimate can be reduced. The two problems were tackled from the angle of model uncertainty by considering six models revised from the DVRFS model with different recharge components and hydrogeological frameworks. The two problems were also tackled from the angle of parametric uncertainty for each model by first conducting Morris sensitivity analysis to identify important parameters and then conducting Monte Carlo simulations for the important parameters. The uncertainty analysis is general and suitable for tackling similar problems; the Morris sensitivity analysis has been utilized to date in only a limited number of regional groundwater modeling. The simulated Qy values were evaluated by using three kinds of calibration data (i.e., hydraulic head observations, discharge estimates, and constant‐head boundary flow estimates). The evaluation results indicate that, within the current DVRFS modeling framework, the Qy estimate can only be reduced to about half of the original estimate without severely deteriorating the goodness‐of‐fit to the calibration data. The evaluation results also indicate that it is necessary to develop a new hydrogeological framework to produce new flow patterns in the DVRFS model. The issues of hydrogeology and boundary flow are being addressed in a new version of the DVRFS model planned for release by the U.S. Geological Survey.  相似文献   

14.
CPT-based seismic stability assessment of a hazardous waste site   总被引：1，自引：0，他引：1  
In areas of high seismicity, seismic stability often controls hazardous and solid waste landfill closure design. The undrained shear strength (Su) of the waste mass is fundamental to seismic slope stability analyses. The value of Su for hazardous waste fill is often difficult to characterize. The physical and chemical natures of the waste fill typically preclude laboratory testing of the materials. In certain cases, Cone Penetration Test (CPT) soundings can provide a viable technique for evaluation of Su provided that the cone shear strength factor Nk can be established. If hazardous waste materials laboratory testing is not an option, Nk may be evaluated based upon results of non-intrusive in situ testing. This paper presents a case history of the seismic stability assessment of a hazardous waste site in which Nk was established from the results of non-intrusive Spectral Analysis of Surface Waves (SASW) soundings and empirical correlations to shear strength of soils. Generalization of the proposed methodology to other sites should be done with caution owing to variability among the parameters used in the analyses.  相似文献   

15.
v--vThe International Monitoring System (IMS) includes a diverse set of seismic arrays with different configurations. These configurations have apertures ranging from less than 1 to more than 25 km and minimum interelement spacings varying from 0.1 to 3.6 km. This paper presents a scheme for initial beam deployment for this variety of seismic arrays. Beamforming is equivalent to a spatiotemporal bandpass filter of which passband is defined by the minimum and maximum wavenumbers, which are functions of the geometry configuration of the array. Deployment of steered-beams for signal detection is based on the wavenumber resolution of the array, slowness and frequency distributions of seismic phases, and coherence properties of seismic signals and noises among sensors. Within the wavenumber passband, all possible slowness values are determined by the resolution for each frequency band, and those that are outside the range of seismological interest are excluded. The appropriate azimuthal distribution for each selected slowness is determined from the azimuthal resolution. Using this approach, detection beams for each array are rationally deployed in the slowness-azimuth and frequency domain.  相似文献   

16.
Nonlinear effects in seismic wave propagation are analyzed to determine the mechanical rigidity of different-order faults that thread the tectonic structures in the central part of the East European platform (Moscow syneclise and Voronezh Crystalline Massif) and the fault zones of the Balapan and Degelen mountain regions in Kazakhstan (the Degelen magmatic node in the Central Chingiz zone). The dependency of the rigidity of the fault zone on the fault’s length is obtained. The rigidity of the tectonic structures is found to experience well-expressed temporal variations with periods of 13–15 days, 27–32 days, and about one year. In the different-order fault zones, the amplitudes of both normal k n and the shear k s rigidity for semimonthly, monthly, and annual variations can span a factor of 1.3, 1.5, and 2.5, respectively.  相似文献   

17.
Following an overview of pertinent literature, this paper presents a new methodology for estimating seismic coefficients for the performance-based design of earth dams and tall embankments. The methodology is based on statistical regression of (decoupled) numerical data for 1084 potential sliding masses, originating from 110 non-linear seismic response analyses of 2D cross sections with height ranging from 20 to 120 m. At first, the methodology estimates the peak value of the seismic coefficient khmax as a function of: the peak ground acceleration at the free field, the predominant period of the seismic excitation, the non-linear fundamental period of dam vibration, the stiffness of the firm foundation soil or rock layer, as well as the geometrical characteristics and the location (upstream or downstream) of the potentially sliding mass. Then, it proceeds to the estimation of an effective value of the seismic coefficient khE, as a percentile of khmax, to be used with a requirement for pseudo-static factor of safety greater or equal to 1.0. The estimation of khE is based on allowable permanent down-slope deviatoric displacement and a conservative consideration of sliding block analysis.  相似文献   

18.
Turbulent magnetofluids appear in various geophysical and astrophysical contexts, in phenomena associated with planets, stars, galaxies and the universe itself. In many cases, large-scale magnetic fields are observed, though a better knowledge of magnetofluid turbulence is needed to more fully understand the dynamo processes that produce them. One approach is to develop the statistical mechanics of ideal (i.e. non-dissipative), incompressible, homogeneous magnetohydrodynamic (MHD) turbulence, known as “absolute equilibrium ensemble” theory, as far as possible by studying model systems with the goal of finding those aspects that survive the introduction of viscosity and resistivity. Here, we review the progress that has been made in this direction. We examine both three-dimensional (3-D) and two-dimensional (2-D) model systems based on discrete Fourier representations. The basic equations are those of incompressible MHD and may include the effects of rotation and/or a mean magnetic field B o. Statistical predictions are that Fourier coefficients of the velocity and magnetic field are zero-mean random variables. However, this is not the case, in general, for we observe non-ergodic behavior in very long time computer simulations of ideal turbulence: low wavenumber Fourier modes that have relatively large means and small standard deviations, i.e. coherent structure. In particular, ergodicity appears strongly broken when B o?=?0 and weakly broken when B o?≠?0. Broken ergodicity in MHD turbulence is explained by an eigenanalysis of modal covariance matrices. This produces a set of modal eigenvalues inversely proportional to the expected energy of their associated eigenvariables. A large disparity in eigenvalues within the same mode (identified by wavevector k ) can occur at low values of wavenumber k?=?| k |, especially when B o?=?0. This disparity breaks the ergodicity of eigenvariables with smallest eigenvalues (largest energies). This leads to coherent structure in models of ideal homogeneous MHD turbulence, which can occur at lowest values of wavenumber k for 3-D cases, and at either lowest or highest k for ideal 2-D magnetofluids. These ideal results appear relevant for unforced, decaying MHD turbulence, so that broken ergodicity effects in MHD turbulence survive dissipation. In comparison, we will also examine ideal hydrodynamic (HD) turbulence, which, in the 3-D case, will be seen to differ fundamentally from ideal MHD turbulence in that coherent structure due to broken ergodicity can only occur at maximum k in numerical simulations. However, a nonzero viscosity eliminates this ideal 3-D HD structure, so that unforced, decaying 3-D HD turbulence is expected to be ergodic. In summary, broken ergodicity in MHD turbulence leads to energetic, large-scale, quasistationary magnetic fields (coherent structures) in numerical models of bounded, turbulent magnetofluids. Thus, broken ergodicity provides a large-scale dynamo mechanism within computer models of homogeneous MHD turbulence. These results may help us to better understand the origin of global magnetic fields in astrophysical and geophysical objects.  相似文献   

19.
It is shown that the interaction of the interplanetary magnetic field (IMF), when it has southward component, with the geomagnetic field leads to the formation of an enhanced pressure layer (EPL) near the magnetopause. Currents flowing on the boundary between the EPL and the magnetosheath prevent the IMF from penetrating the magnetosphere. However, the outward boundary of the EPL is unstable. The interchange instability permanently destroys the EPL. Separate filaments of the EPL move away from the Earth. New colder plasma of the magnetosheath with a frozen magnetic field replaces the hotter EPL plasma, and the process of EPL formation and destruction repeats itself.The instability increment is calculated for various magnitudes of the azimuthal wave number, ky, and curvature radius of the magnetic field lines, Rc. The disturbances with R−1e\leqky\leq4R−1e (where Re is the Earth’s radius) and Rc\simeqRe are the most unstable.A possible result of the interchange instability of the EPL may be patchy reconnection, displayed as flux transfer events (FTEs) near the magnetopause.  相似文献   

20.
Abstract

The influence of mesoscale topography on the baroclinic instability of a two-layer model of the open ocean is considered. For westward velocities in the top layer (U), and for a sinusoidal topography independent of x or longitude (a cross-stream topography), the critical value of U (Uc ) leading to instability is the same as when there is no topography. The wavelength of the unstable perturbation corresponding to U c is shortened. For a given wavevector (k) of the perturbation the system becomes stable (as also in the absence of topography) for large values of |U|. The minimum value of the shear leading to stability is, however, significantly reduced by the topography.

For sufficiently large values of the height of the topographic features, instabilities appear which are localized within a narrow range of the shear. These instabilities are studied for a topography that depends both on x and y.

For a cross-stream topography the growth rates are somewhat smaller than those without topography and they depend only weakly on ky . For the topographies considered here which depend both on x and y, perturbations with different values of ky can again have roughly the same growth rate.

In the case of stable oscillations, variations in the eddy energy with very long periods are made possible by the coexistence of topographic modes with closely lying periods.  相似文献   

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

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