首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Two peculiar characteristics of the VAN method of short-term earthquake prediction are: (1) the existence of sensitive and insensitive sites to seismic electric signal (SES) and (2) the selective sensitivity of a sensitive station to SES from specific focal area(s). The process of their discovery, which so far has not been published adequately, is examined for the Ioannina (IOA) station. The sensitive site was discovered simply through repeatedly moving the temporary observation network. The selectivity map of the IOA station was made by first identifying the earthquake source areas with SESs which were considered to have been detected at IOA, and second by examining the direction of SESs from different such source areas. For the latter purpose, information from long dipoles appears to have been helpful. The selectivity map is empirically constructed and is updated regularly. We have independently tried to construct a selectivity map of IOA by examining VAN's raw SES records and by using the same criteria as VAN. The selectivity map we obtained was not exactly, but roughly similar to that of VAN.  相似文献   

2.
The effects of using different mesh sizes on interpretation of seasonal changes in the fauna of a sandy station off the north-east coast of England was investigated. Retention of smaller animals depended both upon absolute size, and the tendency of some species to fragment during routine processing of samples in the field. Seasonal effects may lead to difficulties in establishing trends in abundance using a larger mesh size, and the implications for certain types of investigation of this and other sources of sampling error are discussed.  相似文献   

3.
地震电信号选择性数值模拟及可能影响因素   总被引:18,自引:10,他引:8       下载免费PDF全文
本文探讨了将商业有限元软件COMSOL Multiphysics应用于地震电信号选择性三维数值模拟方面可能面临的一些问题.在检验算法可靠性的同时,通过选取不同类型的边界条件、不同大小的计算区域、不同的介质电性参数以及地表介质电性非均匀性和地下特定导电通道等不同电性结构模型等,考察了影响地震电信号选择性数值模拟结果的一些可能因素.结果表明,在开展地震电信号选择性数值模拟研究时需对模型参数的选取进行必要的分析,而地表介质电性非均匀性以及地下特定导电通道等可以较好地解释地震电信号选择性现象,其贡献的大小与模型具体的电性结构和对比度等参数有关.  相似文献   

4.
The Newmark average acceleration method is non-dissipative and unconditionally stable, but its accuracy in transient analysis of wave propagation problems depends not only on the spatial discretization but also on the temporal discretization. It has been found that the effects of spatial and temporal discretization when considered separately as commonly done, are far from adequate for most transient analysis. A better criterion manipulating the interdependent relationship between mesh size and time-step magnitude is imperative to achieve sufficiently accurate results of analysis. In this paper, the accuracy of the Newmark method is investigated by considering the two basic sources of errors, namely, numerical amplitude dissipation and velocity dispersion. The effects of both spatial and temporal discretizations are considered. A new technique to describe the characteristics of various frequency spectra is established. A criterion for mesh design and the selection of time-step magnitude is also proposed based on the combined effects of the amplitude dissipation and velocity dispersion. The efficiency and effectiveness of the proposed criterion are demonstrated using two one-dimensional wave propagation problems. A two-dimensional application shows that this criterion is equally applicable to multidimensional problems.  相似文献   

5.
Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

6.
Groundwater flow in a 3-D domain with fracture planes is numerically investigated using the finite element method. A flexible mesh generation method for discretization is proposed in this paper. The method, based on Delaunay triangulation, divides the whole domain into subdomains separated by fracture planes. It then triangulates each subdomain independently into tetrahedra which are further subdivided into hexahedra. By putting together all of the meshes of the subdomain, a finite element mesh of the whole domain is obtained. The appropriateness of the mesh generation method is topologically proved. Several applications of the proposed method are given, and the numerical solutions are in good agreement with those obtained with a structured grid. It is concluded that the proposed mesh generation method can replace the structured grid.  相似文献   

7.
The turbulent flow in the local scour hole around a single non-submerged spur dyke is investigated with both experimental and numerical methods. The experiments are conducted under clear-water scour regime with an impermeable spur dyke. The scour geometry and flow velocities are measured in details with a high-resolution laser displacement meter, electro-magnetic velocimetries and PIV (Particle image velocimetry). A 3D non-linear k-ε model is developed to simulate the complex local flow field around the scour area. The numerical model is formulated using FVM (Finite volume method) on a collocated unstructured mesh, capable of resolving complex geometries and boundaries. It is found that the simulation results are reasonably consistent with those of the experimental measurements. Based on the study results, the nature of the flow structure around a spur dyke with local scour hole is analyzed.  相似文献   

8.
The Varotsos-Alexopoulos-Nomicos(VAN) method of short-term earthquake prediction was introduced in the 1980s. The VAN method enables estimation of the epicenter, magnitude and occurrence time of an impending earthquake by observing transient changes of the electric field of the Earth termed seismic electric signals(SES). Here, we present a few examples of SES observed in various earthquake prone areas worldwide.  相似文献   

9.
对了解地下构造和描绘地下油气分布来说,用地面地震勘探技术可以分辩15~30 m的地层.低分辨率和缺乏详细的速度信息会降低地面地震勘探的效果,在研究更薄细的油气储层和研究已研究过的老区时,我们要采用一些更先进的勘探技术,对储层进行更精确的预测和描绘.本文分析了 VSP、测井和地面地震勘探各自的优势、不足和局限,采用VSP方法,我们只能得到井周围一定空问内的丰富速度信息,而详细的速度模型可以极大地提高地面地震勘探的效果.据此分析,组合VSP和地面地震方法优势的井-地观测系统(3D VSP),可有效地增加信息量,并可通过引入VSP数据处理时所用的基于模型的处理技术,来使井-地观测系统的优势进一步扩大,在此基础上,如果与精细测井相结合,可以更好的寻找和开发剩余油气,这将是我们今后重要的发展方向.  相似文献   

10.
A mesh grading approach based on investigated lump method has been presented for simulating wave propagation in high velocity-contrast media. Unstructured fine grids are used to discretize the low wave-velocity medium in order to ensure the accuracy of numerical computation, and unstructured coarse grids are used for the high wave-velocity medium in order to substantially reduce the computational cost. On the interface, one coarse grid can match the fine grids of arbitrary odd number. The key feature of the proposed method is the constructions of investigated lumps on the interfaces of media. The transition zone, which is commonly used in the discontinuous grid scheme based on the staggered-grid finite-difference method, will not be used any more. Moreover, the computational instability that the discontinuous grid schemes frequently encountered does not arise in the proposed method. The comparisons with the analytical solutions and the application in studying the effects of sedimentary basin demonstrated that the mesh grading approach is a valid, accurate, convenient and flexible algorithm in simulating wave propagations in high velocity-contrast media with irregular interfaces.  相似文献   

11.
Comprehensive studies on the dynamic behaviours of a high rock slope of the Three Gorges shiplock and a blocky arch structure are performed using the distinct element method (DEM). Some additional functions have been incorporated into the current DEM code, including a scheme of non-uniform seismic input on the boundaries and an improved triangular mesh generator for deformable elements. Parameters from the field and laboratory tests of the project are used in the analysis of the shiplock slope. Comparisons between the DEM results and the field measurements under the excavation unloading process are made and the good agreement obtained demonstrates the effectiveness of the DEM for predicting the deformation process of the slope. The above examples of the high rock slope and the blocky arch demonstrate that the effects of different earthquake input, input mechanism and parameters on the dynamic response behaviours and collapse pattern of the systems are significant.  相似文献   

12.
随着计算机技术的飞速发展,大地电磁测深数据的各种反演方法得到空前发展。对于反演过程中的网格剖分问题,前人研究有限。本文就此进行研究,确定影响网格剖分方式的因素,找到一种合适的网格剖分方法。为此,设计一个复杂模型,对正演模拟得到的数据加入随机噪声作为反演的输入数据,选择不同的网格剖分方式,使用DASOCC与NLCG反演方法,以相同参数分别进行反演,对所得结果进行对比分析。结果显示,在设定的模型及参数条件下,针对DASOCC反演方法,纵向及横向网格均不加密,反演效果比较理想,即横向采用一个测点对应一个网格的剖分方式,纵向采用表层加密,100 m以下后一网格为前一网格厚度1.1倍的剖分方式;针对NLCG反演方法,纵向与DASOCC反演方法相同,采用不加密网格的剖分方式得到的结果较好,而横向则建议采用一个测点对应一个网格,并在2个测点之间插入一个网格的剖分方式。  相似文献   

13.
The finite element method is widely employed in numerical analysis of seismic liquefaction at present. However, due to the mesh distortion in strong dynamic action, the results are generally inaccurate and can even stop the calculation. Therefore, an adaptive mesh refinement (AMR) scheme is introduced to improve the numerical results based on a finite element and finite difference coupled dynamic method (FEM–FDM). The changing law about the errors due to the linear and the bilinear recovered fields of coarse meshes, AMR meshes and fine meshes is discussed. The seismic response on soil is compared between different AMR meshes. The AMR strategy can better simulate seismic liquefaction and offer more accurate results than the normal finite element simulation.  相似文献   

14.
The Laplace Transform Galerkin (LTG) method and the Arnoldi modal reduction method (AMRM) have been implemented in finite element schemes designed to solve mass transport problems in porous media by Sudicky [Sudicky, E.A., Water Resour. Res., 25(8) (1989) 1833-46] and Woodbury et al [Woodbury, A.D., Dunbar, W.S., & Nour-Omid, B., Water Resour. Res., 26(10) (1990) 2579-90]. In this work, a comparative analysis of the two methods is performed with attention focused on efficiency and accuracy. The analysis is performed over one- and two-dimensional domains composed of homogeneous and heterogeneous material properties. The results obtained using homogeneous material properties indicate that for a given mesh design the LTG method maintains a higher degree of accuracy than does the AMRM. However, in terms of efficiency, the Arnoldi attains a pre-defined level of accuracy faster than does the LTG method. It is also shown that for problems involving homogeneous material properties the solution obtained using the LTG method on a coarse mesh is comparable in terms of solution time and accuracy to that obtained using the AMRM on a fine mesh. Comparisons similar to those performed using homogeneous material properties are also performed for the case where the hydraulic conductivity field is heterogeneous. For this case, the level of accuracy achieved by the AMRM and the LTG method are similar. However, as with the analysis involving homogeneous material properties, the AMRM is found to be more efficient than the LTG method. It is also shown that for heterogeneous material properties, use of the LTG method under high grid Peclet conditions can be potentially problematic.  相似文献   

15.
三角谱元法及其在地震正演模拟中的应用   总被引:2,自引:2,他引:0       下载免费PDF全文
李琳  刘韬  胡天跃 《地球物理学报》2014,57(4):1224-1234
谱元法(SEM)是基于有限元(FEM)的一种算法,在地震正演模拟中应用广泛,但是大部分研究都是基于四边形网格下的谱元法.本文给出了2阶谱元法在三角网格中(TSEM)的基本原理,包括Lagrange形函数的构建,数值积分公式的选取.在此基础上,分析了2阶TSEM方法的数值频散特性以及稳定性条件,并引入三角网格下3阶有限元方法进行分析对比,数值算例的结果证明2阶TSEM相比于3阶FEM具有更高的计算精度,以及更宽松的稳定条件.最后,本文将TSEM方法应用于中国西部地区的两个含溶洞介质的地质模型中,数值模拟结果表明TSEM方法能够有效地模拟复杂结构的介质,有助于对地震波场传播特性的认识.  相似文献   

16.
Damping solvent extraction is a finite element method for the analysis of unbounded (visco-)elastic media which was suggested by Wolf and Song in 1994. It was originally recommended that the method should be employed with a variable domain size depending on excitation frequency. Furthermore, other researchers who have utilized this method in the context of constant domain size, have often imposed strict conditions on the mesh size for the whole domain, which reduces the effectiveness of the approach. Considering the effect of artificial damping on mesh density selection, the present study introduces damping solvent as a method in which one can relax typical mesh density requirements to a large extent by utilizing a large value of artificial damping. Therefore, it makes the finite element mesh to benefit from a large domain size which improves the results for low frequency range. Moreover, good results are simultaneously obtained for high frequency range due to employing a large value of artificial damping. To illustrate the point, a rigid strip foundation with a cross section of rectangle embedded in half plane is considered. According to some comparison between the results obtained from several finite element meshes, the best one which takes full advantages of a large value of artificial damping for dynamic stiffness coefficients of strip foundation is introduced. These comparisons are carried out on domain size, mesh density and artificial damping.  相似文献   

17.
The finite difference method (FDM) is an important numerical approach for simulating the propagation of seismic waves, and some FDMs can be used to study the impact of the Earth’s curvature and topography over large distances. To efficiently model the effects of the Earth’s irregular topography on the propagation of seismic waves, here we optimize a previously proposed grid mesh method and develop a novel two-dimensional boundary-conforming FDM based on a curvilinear polar coordinate system. This method efficiently simulates the propagation of seismic waves in an arc-shaped model with large variations in surface topography. Our method was benchmarked against other reported methods using several global-scale models. The consistency of the results confirms the validity of our proposed optimization strategy. Furthermore, our findings indicate that the proposed optimization strategy improves computational efficiency.  相似文献   

18.
Implementation of ADI-FDTD subgrids in ground penetrating radar FDTD models   总被引:1,自引:0,他引:1  
Realistic numerical modeling of ground penetrating radar (GPR) using the finite-difference time-domain (FDTD) method could greatly benefit from the implementation of subgrids – supporting finer spatial resolution – into the conventional FDTD mesh. This is particularly important, when parts of the computational domain need to be modeled in detail or when there are features or regions in the overall computational mesh with values of high dielectric constant supporting propagation of waves at very short wavelengths. A scheme that simplifies the process of implementing these subgrids into the traditional FDTD method is presented. This scheme is based on the combination of the standard FDTD method and the unconditionally stable alternating-direction implicit (ADI) FDTD technique. Because ADI-FDTD is unconditionally stable its time-step can be set to any value that facilitates the accurate calculation of the fields. By doing so, the two grids can efficiently communicate information across their boundary without requiring to use a costly time-interpolation scheme. This paper discusses the performance of ADI-FDTD subgrids when implemented into the traditional FDTD method, using different communication schemes for the information exchange at the boundary of the two grids. The developed algorithm, can handle cases where the subgrid crosses dielectrically inhomogeneous media. In addition, results from the comparison between the proposed scheme and a commonly employed purely FDTD subgridding technique are presented.  相似文献   

19.
— Accurate simulation of seismic ground motion for three-dimensionally complex topography and structures is one of the most important goals of strong motion seismology. The finite-element method (FEM) is well suited for this kind of simulation, since traction-free conditions are already included in the formulation, and the Courant condition is less strict than for the finite-difference method (FDM). However, the FEM usually requires both large memory and computation time. These limitations can be overcome by using a mesh consisting of voxels (rectangular prisms) with isotropy built into the explicit formulation of the dynamic matrix equation. Since operators in the voxel FEM are the combinations of ordinary FDM operators and additional terms, the method keeps accuracy of the same order as FDM and the terms relax the Courant condition. The voxel FEM requires a similar amount of memory and only takes 1.2~1.4 times longer computation time. The voxel mesh can be generated considerably faster than the popular tetrahedral mesh. Both ground motions and static displacements due to a point or line source can be calculated using the voxel FEM approach. Comparisons with the reflectivity method and theoretical solutions demonstrate the successful implementation of the method, which is then applied to more complex problems.  相似文献   

20.
《国际泥沙研究》2020,35(5):540-549
Land use is an important factor influencing soil erosion and sediment yield (SESY). Regressing soil erosion intensity (SEI) and sediment yield (SY) to land use characteristics can provide necessary information for controlling soil loss. However, current simple regression methods emphasize cross sectional parameters, with less emphasis on temporal variability of relevant land use parameters so that the derived effects of land use change on SESY can be biased. Here, a panel data method was applied to quantify the impact of land use change on SESY in 1954, 1975, and 2015, based on the WaTEM/SEDEM model and seven landscape metrics for 25 reservoir catchments in northeastern China. The results indicate that SEI and area-specific SY (SSY) continuously decreased from 1954 to 2015, which were significantly correlated with landscape metrics such as area-edge metrics of mean patch area (AREA_MN), shape index of the mean related circumscribing circle (CIRCLE_MN), aggregation index of effective mesh size (MESH), patch cohesion index (COHESION), and diversity metrics such as Shannon's diversity index (SHDI), patch richness density (PRD), and modified Simpson's evenness index (MSIEI). The results suggested that catchment SESY can be reduced through decreasing mean patch area, patch mesh size, and physical connectivity of patches, enriching landscape types, and elongating land use patches. These findings are helpful to effectively implement soil conservation measures in northeastern China and similar regions worldwide. The current study also implies that the panel data approach will have beneficial potential applications in earth-science research fields.  相似文献   

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

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