弹性波数值模拟的混合边界与频散抑制

针对刚性边界条件、衰减边界条件与透射边界条件对弹性波场的影响及计算中产生的数值频散现象,将衰减边界与透射边界双重作用于人工边界,引入了混合边界条件,并在计算中采用了微量阻尼的方法,结果表明,混合边界条件能更好地吸收人工边界的反射;数值频散的阻尼算法能大大减弱由网格剖分、计算误差造成的频散现象。该方法可有效提高波场模拟精度。

An intersection method for locating earthquakes in complex velocity models

The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.

The Precise Finite Difference Method for Seismic Modeling

D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.

Three-dimensional inversion of borehole-surface electrical data based on quasi-analytical approximation

3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation （QA） and re-weighted regularized conjugate gradient method （RRCG） algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.

Improvements in seismic event locations in a deep western U.S. coal mine using tomographic velocity models and an evolutionary search algorithm

Methods of improving seismic event locations were investigated as part of a research study aimed at reducing ground control safety hazards. Seismic event waveforms collected with a 23-station three-dimensional sensor array during longwall coal mining provide the data set used in the analyses. A spatially variable seismic velocity model is constructed using seismic event sources in a passive tomographic method. The resulting three-dimensional velocity model is used to relocate seismic event positions. An evolutionary optimization algorithm is implemented and used in both the velocity model development and in seeking improved event location solutions. Results obtained using the different velocity models are compared. The combination of the tomographic velocity model development and evolutionary search algorithm provides improvement to the event locations.

Deformation Patterns and Shortening Rates in the East Part of the Kalpin Thrust System in the Southwest Tianshan Mountains During the Late Quaternary

Deformation patterns, shortening amounts and rates in the late Quaternary across the Kalpin thrust system have received tittle attention in the past. This paper attempts to discuss them, mainly in the eastern part of the thrust system by doing field investigation along the faults and folds, measuring geomorphic deformation, excavating trenches in several important sites where young alluvial fans were obviously displaced and dating young deposits of alluvial terraces. There are two types of deformation in the surface and near surface for the Kalpin thrust system in the late Quaternary. They are movement of thrust faults on lower angles and bending of young folds. Both kinds of deformation are shown by shortening and uplifting of young geomorphic surfaces. The surface ages of 3 stages are calculated by dating 20 examples using the TL method in the study area and comparing the results of our predecessors on the deposition and incision times of alluvial terraces in the Tianshan mountain which are 100ka B. P., 33 - 18ka B.P. and 6.6 - 8.2ka B.P. respectively for the large-scale deformed alluvial surfaces： T3, T2 and T1 in the Kalpin region. Then, 19 sets of shortening amounts and rates are obtained in 13 sites along 4 rows of anticlines in front of the Kalpin thrust system and Piqiang fold. The shortening amounts and rates show that there are two sections where deformation is stronger than others. The two sections consist of two arcs that are towards the south. The shortening rates near the top of arcs are 1.32mm/a in the west and 1.39mm/a in the east across the thrust system, respectively. In addition, deformation is stronger in the front rows than the rear ones for bifurcate folds.

Influence of Inhomogeneity on Critical Behavior of Earthquake Model on Random Graph

We consider the earthquake model on a random graph. A detailed analysis of the probability distribution of the size of the avalanches will be given. The model with different inhomogeneities is studied in order to compare the critical behavior of different systems. The results indicate that with the increase of the inhomogeneities, the avalanche exponents reduce, i.e., the different numbers of defects cause different critical behaviors of the system. This is virtually ascribed to the dynamical perturbation.

Influence of Structure Plane Size on Seismic Response of Soil-Structure Interaction

The influence of the change of structure plane size on seismic response was studied for a soil-structure interaction system.Based on the finite element method,a soil-structure interaction calculation model was established to analyze the seismic response by changing the structure plane size and choosing different earthquake waves for different soil fields.The results show that when the natural periods of vibration for different structure plane sizes are close,under the same earthquake wave,the total displacement on the top layer of the structure and the foundation rotation displacement decrease with the increase of structure plane size,and the proportion of superstructure elastic selfdeformation displacement to the total displacement increases with the increase of structure plane size.While for different types of sites and seismic waves,under the horizontal and vertical seismic waves,the seismic responses of different plane sizes have a similar change rule.

Modeling seismic wave propagation in a coal-bearing porous medium by a staggered-grid finite difference method

A staggered-grid finite difference method is used to model seismic wave records in a coal bearing, porous medium. The variables analyzed include the order of the difference calculations, the use of a perfect match layer to provide absorbing boundary conditions, the source location, the stability conditions, and dispersion in the medium. The results show that the location of the first derivative of the dynamic variable with respect to space is coincident with the location of the first derivative of the kinematic varable with respect to time. Outgoing waves are effectively absorbed and reflection at the boundary is very weak when more than 20 perfect match layer cells are used. Blot theory considers the liquid phase to be homogeneous so the ratio of liquid to solid exposure of the seismic source depends upon the medium porosity. Numerical dispersion and generation of false frequencies is reduced by increasing the accuracy of the difference calculations and by reducing the grid size and time step. Temporal second order accuracy, a tenth order spatial accuracy, and a wavelength over more than ten grid points gave acceptable numerical results. Larger grid step sizes in the lateral direction and smaller grid sizes in the vertical direction allow control of dispersion when the medium is a low speed body. This provides a useful way to simulate seismic waves in a porous coal bearing medium.