共查询到20条相似文献，搜索用时 921 毫秒
1.
Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources. 相似文献
2.
Wang Jia Fujisaki‑Manome Ayumi Kessler James Cannon David Chu Philip 《Ocean Dynamics》2023,73(7):407-429
This study investigates the inertial stability properties and phase error of numerical time integration schemes in several widely-used ocean and atmospheric models. These schemes include the most widely used centered differencing (i.e., leapfrog scheme or the 3-time step scheme at n-1, n, n+1) and 2-time step (n, n+1) 1st-order Euler forward schemes, as well as 2nd-stage and 3rd- and 4th-stage Euler predictor-corrector (PC) schemes. Previous work has proved that the leapfrog scheme is neutrally stable with respect to the Coriolis force, with perfect inertial motion preservation, an amplification factor (AF) equal to unity, and a minor overestimation of the phase speed. The 1st-order Euler forward scheme, on the other hand, is known to be unconditionally inertially unstable since its AF is always greater than unity. In this study, it is shown that 3rd- and 4th-order predictor-corrector schemes 1) are inertially stable with weak damping if the Coriolis terms are equally split to n+1 (new value) and n (old value); and 2) introduce an artificial computational mode. The inevitable phase error associated with the Coriolis parameter is analyzed in depth for all numerical schemes. Some schemes (leapfrog and 2nd-stage PC schemes) overestimate the phase speed, while the others (1st-order Euler forward, 3rd- and 4th-stage PC schemes) underestimate it. To preserve phase speed as best as possible in a numerical model, alternating a scheme that overestimates the phase speed with a scheme that underestimates the phase speed is recommended. Considering all properties investigated, the leapfrog scheme is still highly recommended for a time integration scheme. As an example, a comparison between a leapfrog scheme and a 1st-order Euler forward scheme is presented to show that the leapfrog scheme reproduces much better vertical thermal stratification and circulation in the weakly-stratified Great Lakes.
相似文献3.
4.
The main aim of this work is to develop, verify and apply in simulation study an efficient hybrid approach to study seismic response of a soil-structure system taking into account all the important components as: (1) the line time-harmonic source with its specific geophysical properties; (2) the inhomogeneity and heterogeneity of the wave path from the source to the local geological region; (3) the geotechnical properties of the near-field local geological profile and finally (4) the properties of the engineering structure itself. Plane strain state is considered. The hybrid computational tool is based on the boundary element method (BEM1) for modeling the infinite far-field geological media and finite element method (FEM2) for treating the dynamic behavior of the structure and the near-field finite soil geological region. Each of the two techniques is applied in that part of the whole model where it works more efficiently. The hybrid numerical scheme is realized via the sub-structure approach, direct BEM1, conventional FEM2 and insertion of the BEM1 model of the seismically active far-field geological media as a macro-finite element (MFE3) in the FEM2 commercial program ANSYS. The accuracy and verification study of the proposed method is presented by solution of numerical test examples simulating different seismic scenarios. The obtained results show clearly that the hybrid model is able to demonstrate the sensitivity of the synthetic signals to the source properties, to the heterogeneous character of the wave path, to the relief peculiarities of the local layered geological deposit and to the specific properties of the engineering structure. 相似文献
5.
Full-wavefield inversion for distributions of acoustic velocity, density and Q on a vertical slice through a25D model is implemented for common-source gathers in a cross-hole geometry. The wavefield extrapolation used is 3D, so all geometrical spreading, scattering, reflection, and transmission effects are correctly and automatically compensated for. In order to keep the number of unknowns tractable, application was limited to 2.5D models of known geometry; the latter assurnes a prior step, such as tomography, to fix the layer geometries. With the model geometry fixed, reliable solutions are obtained using synthetic data from only two independent source locations. Solutions from data with noisy and missing traces are comparable to those from noise-free data, but with higher residuals. When the source locations are spatially widely separated, conunon-source gathers may be summed and treated as a single wavefield to yield the same model estimates as when the individual source wavefields are treated separately, at substantially reduced cost. Inversions for full 3D parameter distributions can be handled with the same software, requiring only solution for more unknowns. 相似文献
6.
L. Alvarez G. F. Panza F. Vaccari B. E. González 《Pure and Applied Geophysics》2001,158(9-10):1763-1782
— We present the results of complete P-SV and SH waves modelling, up to a maximum frequency of 1?Hz, along two profiles in Santiago de Cuba city. The seismic sources are located in the depth range from 10 to 40?km on the Oriente fault zone at distances of several tens of kilometres from the city. The calculation has been made by a hybrid method: Modal summation in the regional anelastic model (one-dimensional) where the source is buried, and finite differences in the local sedimentary anelastic models (two-dimensional). The analysis of the influence of the depth and of the distance of the source on the site effects shows that standard traditional methods, based on the deconvolution analysis of the rock outcrop motion, can lead to erroneous results. 相似文献
7.
Antoine Guitton 《Geophysical Prospecting》2006,54(2):135-152
Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries. 相似文献
8.
Stability and grid dispersion in the P-SV 4th-order in space, 2nd-order in time, displacement-stress staggered-grid finite-difference scheme is investigated in the case of a homogeneous unbounded medium. All results, however, also apply to the velocity-stress and displacement- velocity-stress finite-difference schemes.
Independent stability conditions for the P and S waves are obtained by exact separation of equations for the two types of waves.
Since the S-wave group velocity can differ from the actual velocity as much as 5% for the sampling ratio 1/5, commonly used in numerical modelling, the sampling of the minimum S wavelength by 6 grid spacings (with the velocity difference not larger than 2.5%) is recommended.
Grid dispersion is strongest for a wave propagating in a direction of a coordinate axis and weakest for a wave propagating along a plane diagonal.
Grid dispersion in the 4
th
-order scheme for the sampling ratios s = 1/5 and s = 1/6 is smaller than grid dispersion in the 2
nd
-order scheme for s = 1/10 and s = 1/12, respectively. 相似文献
9.
本文针对弹性波逆时偏移,提出稳定的激发振幅成像条件.在震源波场的正向传播过程中,计算每个网格点的能量,并保存最大能量密度的时刻和相应的波场值;在检波器波场的逆时传播过程中,在每个网格点提取最大能量密度时刻的检波器波场值,并利用保存的最大能量震源波场做归一化,获得角度依赖的反射系数成像剖面.相比于归一化互相关成像条件,该成像条件在震源波场的正向传播过程中无需存储波场快照,节省大量磁盘空间和I/O吞吐任务,提高了计算效率;相比于弹性波的激发时间成像条件,该成像条件自动校正了水平分量在震源两侧的极性反转,在多炮叠加时避免振幅损失.数值试验表明,与归一化成像条件相比,稳定激发振幅成像条件具有更小的计算量,偏移剖面的低频假象更弱,水平分量的成像能力更优,具有更高的空间分辨率. 相似文献
10.
Pedro Rau Luc Bourrel David Labat Denis Ruelland Frdric Frappart Waldo Lavado Boris Dewitte Oscar Felipe 《水文研究》2019,33(1):20-35
In a context of water scarcity in Peruvian Pacific catchments as a crucial issue for Peru, added to the paucity of data availability, we propose a methodology that provides new perspectives for freshwater availability estimation as a base reference for unimpaired conditions. Under those considerations, a regional discharge of 709 m3/s to the Pacific Ocean is estimated with a significant increasing trend of about 43 m3/s per decade over the 1970 – 2010 period. To represent the multidecadal behaviour of freshwater runoff along the region, a regional runoff analysis is proposed based on hydrological modelling at annual and monthly time step for unimpaired conditions over the whole 1970 – 2010 period. Differential Split‐Sample Tests are used to assess the hydrological modelling robustness of the GR1A and GR2M conceptual lumped models, showing a satisfactory transposability from dry to wet years inside the thresholds defined for Nash–Sutcliffe and bias criteria. This allowed relating physical catchment characteristics with calibrated and validated model parameters, thus offering a regional perspective for dryland conditions in the study area (e.g., the anticlockwise hysteresis relationship found for seasonal precipitation–runoff relationship) as well as the impacts of climate variability and catchment characteristics. 相似文献
11.
12.
The rotated Cartesian coordinate method to remove the axial singularity of cylindrical coordinates in finite‐difference schemes for elastic and viscoelastic waves 
下载免费PDF全文
下载免费PDF全文 When modelling the propagation of 3D non‐axisymmetric elastic and viscoelastic waves in cylindrical coordinates using the finite‐difference time‐domain method, a mathematical singularity occurs due to the presence of terms in the elastic and viscoelastic wave equations. For many years, this issue has been impeding the accurate numerical solution near the axis. In this work, we propose a simple but effective method for the treatment of this numerical singularity problem. By rotating the Cartesian coordinate system around the z‐axis in cylindrical coordinates, the numerical singularity problems in both 2D and 3D cylindrical coordinates can be removed. This algorithm has three advantages over the conventional treatment techniques: (i) the excitation source can be directly loaded at , (ii) the central difference scheme with second‐order accuracy is maintained, and (iii) the stability condition at the axis is consistent with the finite‐difference time‐domain in Cartesian coordinates. This method is verified by several 3D numerical examples. Results show that the rotating the Cartesian coordinate method is accurate and stable at the singularity axis. The improved finite‐difference time‐domain algorithm is also applied to sonic logging simulations in non‐axisymmetric formations and sources. 相似文献
13.
Interhemispheric transport is a key process affecting the accuracy of source quantification for species such as methane by inverse modelling, and is a source of difference among global three-dimensional chemistry transport models (CTMs). Here we use long-term observations of the atmospheric concentration of long-lived species such as CH3CCl3 and CFCl3 for testing three-dimensional chemistry transport models (CTMs); notably their ability to model the interhemispheric transport, distribution, trend, and variability of trace gases in the troposphere. The very striking contrast between the inhomogeneous source distribution and the nearly homogeneous trend, observed in the global ALE/GAGE experiments for both CH3CCl3 and CFCl3 illustrates an efficient interhemispheric transport of atmospherically long-lived chemical species. Analysis of the modelling data at two tropical stations, Barbados (13°N, 59°W) and Samoa (14°S, 124°W), show the close relationship between inter-hemispheric transport and cross-equator Hadley circulations. We found that cross-equator Hadley circulations play a key role in producing the globally homogeneous observed trends. Chemically, the most rapid interaction between CH3CCl3 and OH occurs in the northern summer troposphere; while the most rapid photolysis of CH3CCl3 and CFCl3, and the chemical reactions between CFCl3 and O(1D), take place in the southern summer stratosphere. Therefore, the cross-equator Hadley circulation plays a key role which regulates the southward flux of chemical species. The regulation by the Hadley circulations hence determines the amount of air to be processed by OH, O(1D), and ultraviolet photolysis, in both hemispheres. In summary, the dynamic regulation of the Hadley circulations, and the chemical processing (which crucially depends on the concentration of OH, O(1D), and on the intensity of solar insolation) of the air contribute to the seasonal variability and homogeneous growth rate of observed CH3CCl3 and CFCl3. 相似文献
14.
为减少直接分析三维大尺度复杂土-结构动力相互作用问题的计算量,提高计算效率,本文直接从波动方程出发,提出了较常规子结构法更简单的两步简化计算过程,即第一步简化上部复杂结构体系为集中质量杆系模型,并求基础处等效输入,第二步通过等效输入求上部结构各位置的动力反应.其中第一步计算主要采用集中质量显式有限单元法结合局部透射人工... 相似文献
15.
The present study aims to demonstrate how the dynamic behaviour of structures is affected by local inhomogeneities inside the soil. Particularly, the influence of block-shaped elastic inclusions directly beneath a vertically vibrating block foundation is considered. Comparisons between a 2D and a 3D modelling of the soil—structure system are made. It can be observed that, depending on the stiffness of the inclusion and especially on the excitation frequency, the dynamic response of the foundation may either increase or decrease. 相似文献
16.
The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image. 相似文献
17.
This study deals with the methodical aspects of k
–2(Bernard et al., 1996) kinematic strong motions modelling: (1) it is shown how to incorporate the k-dependent rise time for 2D fault geometry in the strong motion synthesis according to the representation theorem, (2) it is suggested how to produce realistic k
–2 slip models including asperity(ies), (3) modifications are introduced concerning the typeof used slip velocity function and the corner wave number in the slip distribution. High frequency effects of these generalized models are discussed.It is shown that, assuming the rise time proportional to the spatial slip wavelength at high wave numbers, the spectral decay of displacement at frequencies higher than the corner frequency is given just by the decay ofthe slip distribution spectrum, regardless of the type of slip velocity function. It is shown numerically that this model provides -squared source spectrum even in a vicinity of a 2D normal fault buried in 1D structure, which is an agreement with previous studies. 相似文献
18.
19.
G. I. Pugacheva W. N. Spjeldvik A. A. Gusev I. M. Martin N. M. Sobolevsky 《Annales Geophysicae》1998,16(8):931-939
Radial transport theory for inner radiation zone MeV ions has been extended by combining radial diffusive transport and losses due to Coulomb friction with local generation of D, T and 3He ions from nuclear reactions taking place on the inner edge of the inner radiation zone. Based on interactions between high energy trapped protons and upper atmospheric constituents we have included a nuclear reaction yield D, T and 3He flux source that was numerically derived from a nuclear reaction model code originally developed at the Institute of Nuclear Researches in Moscow, Russia. Magnetospheric transport computations have been made covering the L-shell range L=1.0–1.6. The resulting MeV energy D, T and 3He ion flux distributions show a strong influence of the local nuclear source mechanism on the inner zone energetic D, T and 3He ion content. 相似文献