共查询到20条相似文献，搜索用时 468 毫秒
1.
Xiuying Xing Zhiqing Kou Ziyi Huang Jiin-Jen Lee 《Pure and Applied Geophysics》2013,170(6-8):1149-1168
Tsunamis waves caused by submarine earthquake or landslide might contain large wave energy, which could cause significant human loss and property damage locally as well as in distant region. The response of three harbors located at the Pacific coast (i.e. Crescent City Harbor, Los Angeles/Long Beach Port, and San Diego Harbor) to six well-known tsunamis events generated (both near-field and far-field) between 2005 and 2011 are examined and simulated using a hybrid finite element numerical model in frequency domain. The model incorporated the effects of wave refraction, wave diffraction, partial wave reflection from boundaries, entrance and bottom energy dissipation. It can be applied to harbor regions with arbitrary shapes and variable water depth. The computed resonant periods or modes of oscillation for three harbors are in good agreement with the energy spectral analysis of the time series of water surface elevations recorded at tide gauge stations inside three harbors during the six tsunamis events. The computed wave induced currents based on the present model are also in qualitative agreement with some of the reported eye-witness accounts absence of reliable current data. The simulated results show that each harbor responded differently and significantly amplified certain wave period(s) of incident wave trains according to the shape, topography, characteristic dimensions and water depth of the harbor basins. 相似文献
2.
Based on the formulation of a multiple non-isotropic scattering process, a characteristic source time is introduced to define the initial impulse width of energy density at the source. An analytical expression of the initial intensity spectral density of a seismic wave is incorporated into the integral equation of seismic wave energy density. And, a recursive formula of Green’s function is derived to obtain the higher order Green’s function, which is included to describe the stronger non-isotropic scattering process. Then, the effect of the scattering pattern on the energy density envelope is investigated by the modified scattering theory. Significant differences are found in the decay of the energy density envelopes with distances using different scattering patterns. The envelope synthesized by the forward dominated scattering pattern is larger than the results obtained by the isotropic and backward dominated scattering pattern. Different scattering patterns are also used to fit the observation data from the aftershocks of the 2008 Wenchuan earthquake. It is concluded that the envelopes synthesized by the forward scattering pattern can match the data better than the isotropic and backward dominated scattering cases, and a new interpretation of the coda wave is given. Finally, using the forward dominated scattering pattern, the envelope broadening of the observed data is reproduced. 相似文献
3.
Giant earthquakes generate rich signals that can be used to explore the characteristics of the hierarchical structure of the Earth’s interior associated with the eigenfrequencies of the Earth.We employ the spectral element method,incorporated with large-scale parallel computing technology,to investigate the characteristics of global seismic wave propagation excited by the2011 Mw9.0 Tohoku earthquake.The transversely isotropic PREM model is employed as a prototype of our numerical global Earth model.Topographic data and the effect of the oceans are taken into consideration.Wave propagation processes are simulated by solving three-dimensional elastic wave governing equations with the seismic moment tensor obtained from the Global Centroid Moment Tensor Catalog.Three-dimensional visualization of our computing results displays the nature of the global seismic wave propagation.Comparative analysis of our calculations with observations obtained from the Incorporated Research Institutions for Seismology demonstrates the reliability and feasibility of our numerical results.We compare synthetic seismograms with incorporated and unincorporated ocean models.First results show that the oceans have obvious effects on the characteristics of seismic wave propagation.The peak displacement and peak velocity of P waves become relatively small under the effect of the ocean.However,the effect of the ocean on S-waves is complex.The displacement and velocity of S waves decrease rapidly over time using an unincorporated ocean model.Therefore,the effects of the ocean should be incorporated when undertaking quantitative earthquake hazard assessments on coastal areas.In addition,we undertake comparative analysis on the characteristics of the Earth’s oscillation excited by the 2004 Sumatra-Andaman,2008 Wenchuan,and 2011Tohoku earthquakes that incorporate the effect of the Earth’s gravitational potential.A comparison of the amplitude spectra of the numerical records indicates that energy released by the three big earthquakes is different.Our comparative analysis realizes that the computing results can accurately reproduce some eigenfrequencies of the Earth,such as toroidal modes 0T2 to 0T13and spheroidal modes 0S7 to 0S31.These results demonstrate that numerical simulations can be successfully used to investigate the Earth’s oscillations.We propose that numerical simulations can be used as one of the major tools to further reveal how the Earth’s lateral heterogeneities affect the Earth’s oscillations. 相似文献
4.
5.
The 2008 Wenchuan earthquake occurred in an active earthquake zone, i.e., Longmenshan tectonic zone. Seismic waves triggered
by this earthquake can be used to explore the characteristics of the fault rupture process and the hierarchical structure
of the Earth’s interior. We employ spectral element method incorporated with large-scale parallel computing technology, to
investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms
with one-point source model and three-point source model respectively. The AK135 model is employed as a prototype of our numerical
global Earth model. The Earth’s ellipticity, Earth’s medium attenuation, and topography data are taken into consideration.
These wave propagation processes are simulated by solving three-dimensional elastic wave governing equations. Three-dimensional
visualization of our numerical results displays the profile of the seismic wave propagation. The three-point source, which
is proposed from the latest investigations through field observation and reverse estimation, can better demonstrate the spatial
and temporal characteristics of the source rupture process than the one-point source. We take comparison of synthetic seismograms
with observational data recorded at 16 observatory stations. Primary results show that the synthetic seismograms calculated
from three-point source agree well with the observations. This can further reveal that the source rupture process of Wenchuan
earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes.
Supported by National Basic Research Program of China (Grant No. 2004CB418406), National Natural Science Foundation of China
(Grant Nos. 40774049 and 40474038), and Computer Network Information Center, Chinese Academy of Sciences (Grant No. INF105-SCE-02-12) 相似文献
6.
D.C. Cuchiara E.H. Fernandes J.C. Strauch J.C. Winterwerp L.J. Calliari 《Continental Shelf Research》2009
The main objective of this study is the characterization of the wave climate in the Southern Brazilian Shelf (SBS) based on a thorough review of existing field data and on numerical modeling experiments. A quantitative knowledge about the wave climate of this area is important to understand the mechanisms driving episodic mud bank attachments to the sandy shore, and the interaction of these banks with the flow and waves. The statistical analysis of existent data on the wave climate throughout the SBS indicates that the predominant wave directions are 100° and 160° (E–SE), with wave heights varying between 1 and 1.50 m. The wave period varies between 6 and 14 s, with predominance of mean wave period of 8 s (sea conditions) and 12 s (swell conditions). The spectral wave model SWAN version 40.41 [Booij, N., Haagsma, I.J.G., Holthuijsen, L.H., Kieftenburg, A.T.M.M., Ris, R.C., van der Westhuysen, A.J., Zijlema, M., 2004. SWAN Cycle III Version 40.41 Users Manual, Delft University of Technology, Delft, The Netherlands, p. 118, http://fluidmechanics.tudelft.nl/swan/index.htm] is used to simulate the wave climate for the region. Special attention is given to Cassino Beach, describing the wave climate observed during the Cassino Experiment, carried out in 2005. The verification of the standard SWAN model was carried out based on the comparison between numerical modeling results and available data of significant wave height, peak period, mean wave direction and energy density for the period relative to February of 1998. Results showed satisfactory model predictions of significant wave height and reasonably accurate predictions of peak spectral wave period and direction. The model performance is also considered satisfactory in the representation of the wave climate of the region when the wave spectrum has only one spectral peak, but presents limitations for bimodal wave spectrum. When two spectral peaks are observed, the SWAN model agrees with the spectral level observed in the low frequency, but underestimates the spectral level in the high-frequency band. When considering the presence of mud deposits in the area, model results predict that although the presence of mud attenuates most of the wave energy on the low frequency peak, it has a smaller effect in attenuating the wave energy on the high frequency peak. 相似文献
7.
8.
M. Ioualalen H. Perfettini S. Yauri Condo C. Jimenez H. Tavera 《Pure and Applied Geophysics》2013,170(3):433-451
Following the 2007, August 15th, M w 8.0, Pisco earthquake in central Peru, Sladen et al. (J Geophys Res 115: B02405, 2010) have derived several slip models of this event. They inverted teleseismic data together with geodetic (InSAR) measurements to look for the co-seismic slip distribution on the fault plane, considering those data sets separately or jointly. But how close to the real slip distribution are those inverted slip models? To answer this crucial question, the authors generated some tsunami records based on their slip models and compared them to DART buoys, tsunami records, and available runup data. Such an approach requires a robust and accurate tsunami model (non-linear, dispersive, accurate bathymetry and topography, etc.) otherwise the differences between the data and the model may be attributed to the slip models themselves, though they arise from an incomplete tsunami simulation. The accuracy of a numerical tsunami simulation strongly depends, among others, on two important constraints: (i) A fine computational grid (and thus the bathymetry and topography data sets used) which is not always available, unfortunately, and (ii) a realistic tsunami propagation model including dispersion. Here, we extend Sladen’s work using newly available data, namely a tide gauge record at Callao (Lima harbor) and the Chilean DART buoy record, while considering a complete set of runup data along with a more realistic tsunami numerical that accounts for dispersion, and also considering a fine-resolution computational grid, which is essential. Through these accurate numerical simulations we infer that the InSAR-based model is in better agreement with the tsunami data, studying the case of the Pisco earthquake indicating that geodetic data seems essential to recover the final co-seismic slip distribution on the rupture plane. Slip models based on teleseismic data are unable to describe the observed tsunami, suggesting that a significant amount of co-seismic slip may have been aseismic. Finally, we compute the runup distribution along the central part of the Peruvian coast to better understand the wave amplification/attenuation processes of the tsunami generated by the Pisco earthquake. 相似文献
9.
Nuree Han Myung Jin Nam Hee Joon Kim Tae Jong Lee Yoonho Song Jung Hee Suh 《Journal of Applied Geophysics》2009,68(4):533-545
Magnetotelluric (MT) surveys were conducted in Pohang, Korea, for low-temperature geothermal exploration in 2002 and 2003. Pohang is located in the southeastern part of the Korean Peninsula and close to the East Sea. In the interpretation of MT data from a coastal environment, sea effects must be correctly included because seawater is a strong conductor. We first constructed a five-layered earth model with a realistic coastline and bathymetry to investigate sea effects on MT data measured in Pohang. This model clearly shows that the Pohang data are significantly influenced by sea water at frequencies blow 1 Hz at the whole measurement sites. Next, we utilized a three-dimensional inversion algorithm based on the Gauss–Newton approach to produce a reliable resistivity model. Seawater is excluded from the inversion domain to fix the resistivity, while included in the modeling domain to simulate sea effects on MT responses. Blocks for the sub-seafloor are included in unknown parameters since they are sufficiently close to the survey area to affect MT responses in Pohang. Static shifts are also considered in inversion for more accurate interpretation. The rms data misfit is smoothly reduced from 11.2 to 1.87 after 7 iterations. The resulting resistivity model shows a pattern of low–high–low resistivity with depth. The model is compatible with resistivity logs obtained from four boreholes in the survey area, and can explain major geological features in Pohang. 相似文献
10.
A discrete numerical adjoint has recently been developed for the stochastic wave model SWAN. In the present study, this adjoint code is used to construct spectral sensitivity maps for two nearshore domains. The maps display the correlations of spectral energy levels throughout the domain with the observed energy levels at a selected location or region of interest (LOI/ROI), providing a full spectrum of values at all locations in the domain. We investigate the effectiveness of sensitivity maps based on significant wave height (H s ) in determining alternate offshore instrument deployment sites when a chosen nearshore location or region is inaccessible. Wave and bathymetry datasets are employed from one shallower, small-scale domain (Duck, NC) and one deeper, larger-scale domain (San Diego, CA). The effects of seasonal changes in wave climate, errors in bathymetry, and multiple assimilation points on sensitivity map shapes and model performance are investigated. Model accuracy is evaluated by comparing spectral statistics as well as with an RMS skill score, which estimates a mean model–data error across all spectral bins. Results indicate that data assimilation from identified high-sensitivity alternate locations consistently improves model performance at nearshore LOIs, while assimilation from low-sensitivity locations results in lesser or no improvement. Use of sub-sampled or alongshore-averaged bathymetry has a domain-specific effect on model performance when assimilating from a high-sensitivity alternate location. When multiple alternate assimilation locations are used from areas of lower sensitivity, model performance may be worse than with a single, high-sensitivity assimilation point. 相似文献
11.
Soil amplification characteristics are investigated using data from the Chibaken‐Toho‐Oki earthquake and its aftershocks recorded at Chiba dense array in Japan. The frequency‐dependent amplification function of soil is calculated using uphole‐to‐downhole spectral ratio analysis, considering the horizontal components of shear wave. The identified spectral ratios consistently demonstrate the splitting of peaks in their resonance frequencies and low amplification values in comparison with a 1D model. The torsional behaviour and horizontal ground motion coupling are clarified as the reasons for these phenomena at the site. To prove the hypothesis, the torsional motion is directly evaluated using the data of the horizontal dense array in different depths at the site. The comparison between Fourier spectra of torsional motion and identified transfer functions reveals the peaks at the same frequencies. The wave equation including torsion and horizontal motion coupling is introduced and solved for the layered media by applying wave propagation theory. Using the developed model, the effects of torsional motion with horizontal motion coupling on soil transfer function are numerically examined. Splitting and low amplification at resonance frequencies are confirmed by the results of numerical analysis. Furthermore, the ground motion in two horizontal directions at the site is simulated using site geotechnical specification and optimizing the model parameters. The simulated and recorded motions demonstrate good agreement that is used to validate the hypothesis. In addition, the spectral density of torsional ground motions are compared with the calculated one and found to be well predicted by the model. Finally, the results are used to explain the overestimation of damping in back‐calculation of dynamic soil properties using vertical array data in small strain level. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
12.
This paper presents the analysis of measurements of long waves at the harbor of Marina di Carrara, Italy. Nine pressure gauges
(eight in the inner harbor and one outside the harbor) were used to record continuously the water surface elevation, for about
3 years; incoming short waves were measured using a directional buoy at about 500 m offshore the breakwaters. The analysis
is carried out splitting the energy of the long waves into two bands, i.e., very long waves (VLW; f < 0.003 Hz) and long waves (LW; 0.003 < f < 0.030 Hz); it is found that LW energy is strongly correlated with the energy of the incoming short waves while no correlation
(or very poor) is found for VLW. During energetic sea states, surface elevation spectra of VLW and LW, normalized against
the wave energy of each frequency band, appear to be autosimilar; this suggests that the spectra of the incoming long waves
have a mostly constant shape. Simple numerical computations of harbor resonance, carried out using a mild-slope equation linearized
model, indicate that the amplification of the height of the incoming LW generated by short waves (calculated using empirical
formulae from previous researches by Melito et al. 2007) inside the harbor can be reasonably predicted using such simple approach, assuming that the LW spectra are rectangular.
The shape of the amplification diagram calculated by the numerical model shows reasonable (but not perfect) agreement with
the shape of the normalized measured spectra inside the harbor. 相似文献
13.
14.
本文推导了声波方程频散函数, 分析了伪谱法的空间网格大小和采样周期对数值频散的影响, 通过数值模拟实验得到了最佳空间参数选择方法。 结果表明: 伪谱方法稳定数值模拟的最大空间采样间距选取原则是使中波长(奈奎斯特频率的一半)的采样点数为2个; 对于所有维度, 稳定性随空间采样间距的增加而增加, 但不易变化太大, 变化太大时需要适当减小震源子波的主频, 以满足空间合理采样; 空间采样间距的大小设置, 需要考虑满足采样定理和稳定性计算条件, 并且稳定性条件对空间采样间距的要求更加严格; 伪谱法数值模拟的最佳(数值频散最小)空间参数选择为中波长2个采样点, 对应主波长约6~7个采样点。 以上研究对于采用伪谱法进行声波方程数值模拟过程中, 如何合理选择模拟参数提供一些参考。 相似文献
15.
在多次各向异性散射理论的基础上,本文重新推导了方向性散射系数的球函数展开式.引入特征时间的概念,来定义震源处初始地震波脉冲宽度,并在地震波能量密度积分方程中引入任意给定频率的初始脉冲能量谱密度的解析表达.通过离散波数方法求解了修正的地震波能量密度积分方程.基于积分方程的数值解,研究了不同散射模式对S波能量密度包络曲线的影响.计算结果表明:随着震源距的增加,在S波到时之后,多次各向异性散射模式与多次各向同性散射模式合成的能量密度包络差异逐渐增大.其中通过多次前散射模式,我们可以得到不同震源距的尾波能量密度包络的同一衰减趋势,以及S波能量密度包络随着震源距的增加而出现的展宽现象.最后,利用美国内华达州Wells地震余震的台站记录验证了多次前散射模式的实用性与有效性. 相似文献
16.
O. S. Golovataya F. F. Gorbatsevich H. Kern T. Popp 《Izvestiya Physics of the Solid Earth》2006,42(11):865-876
Using core samples of the Kola ultradeep borehole (SG-3) and their surface analogues, variations in the density and elastic properties of some crystalline rocks of the Earth’s crust are estimated by modeling of in situ conditions. It is shown that the bulk density and the elastic wave velocities in the rocks have a weak depth gradient. In the SG-3 section under consideration, this gradient is negative. The resulting dependences for estimating the variations in the properties of the crystalline rocks are suitable for the depth range from the surface to 20–30 km. The initial data for the linear approximation of the characteristics can be obtained from the results of tests of surface analogue samples. It is shown that the velocity anisotropy of the metamorphic rocks can vary within wide limits. 相似文献
17.
18.
Specific features of the spatial-temporal dynamics of LF disturbances of the Earth’s magnetosphere have been analyzed by the method of numerical simulation of magnetic hydrodynamic equations taking into account plasma temperature variations. The effects related to the appearance of long-living magnetospheric jumps of density and temperature, MHD wave scattering on such inhomogeneities, Alfvén wave reflection from the near-Earth region, and magnetospheric plasma heating as a result of dissipative processes at a repeated propagation of MHD waves between magnetically conjugate regions have been considered. The problems of conformity of the discrete mathematical model with continuous equations of magnetic hydrodynamics are discussed. 相似文献
19.
本文首先从电子密度及电子温度的输运方程和考虑自作用时的电磁波波动方程出发,利用简正模展开的方法推导出泵波在反射区域激发出热自聚焦不稳定性(thermal self-focusing instabilities,TSFI)所需电场阈值以及其增长率的完整数学表达式,并估算了TSFI激发阈值及所对应的有效辐射功率(ERP)的量级.随后利用三维垂直加热的理论模型,结合国际参考电离层(IRI-2012)和中性大气模型(MSIS-E-00)给出的背景参数,数值模拟了大功率高频泵波加热电离层时泵波反射区域电子密度及电子温度因TSFI而产生的变化及发展的过程,并对比分析了不同背景参数对较热效果的影响.结果表明:当高频泵波的加热阈值达到或超过百毫伏每米的量级时,即可激发TSFI,发展出大尺度电子密度及温度不均匀体,这些不均匀体内的密度耗空约为4%~10%,而电子温度剧烈增长,到达背景温度值的1.6~2.1倍;且在相当的加热条件下,背景电子温度越低、电子密度越小,加热效果越显著;电子密度及电子温度的扰动幅度随着加热时间的推移而逐渐减小,即扰动逐渐趋于饱和,且电子温度要快于电子密度达到饱和状态.本文还对泵波反射高度处的电子密度及电子温度变化率进行采样并求得其功率谱密度,分析结果表明:TSFI发展出的大尺度不均匀体满足幂律谱结构,谱指数随着加热的进行逐渐趋于稳定,白天与夜间的幂律谱指数区别不大,但电子密度与电子温度的幂律谱有所区别. 相似文献
20.
转换波方位各向异性裂缝检测技术研究及应用 总被引：1，自引：0，他引：1
HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果. 相似文献