基于小波分析的场地波速测量方法

在场地波速测量中,由于噪声等因素的影响很难准确识别P、S波的初至时刻,致使波速结果存在很大的误差。本文给出了一种基于小波变换的波速测量的新方法。该方法利用波动信号的小波变换与弹性波群速度的关系准确识别弹性波初至时刻。弹性波小波变换的峰值时刻代表着以群速度传播的弹性波的初至时刻,使P波、S波的初至时刻的确定具有明确的物理意义,波速的结果准确、可靠、稳定。此外,波动信号的小波多尺度分析还可以确定地层中传播的弹性波的频散特性。最后,该方法在场地波速测量的实测信号的应用表明该方法可准确确定P、S波速。

Research on the Seismic Wave Field of Karst Cavern Reservoirs near Deep Carbonate Weathered Crusts

Fracture and cavern hydrocarbon reservoirs in carbonates are an important pool type worldwide. The karst cavern reservoirs are easiest to identify on seismic reflection data. The prediction, exploration, and development of this type of reservoir require theoretical research on seismic wave fields reflected from complex inhomogeneous media. We compute synthetic seismic sections for fluidfilled cavern reservoirs of various heights and widths using random media models and inhomogeneous media elastic wave equations. Results indicate that even caverns significantly smaller than 1/ 4 wavelength are detectible on conventional band-width seismic sections as diffractions migrated into bead-type events. Diffraction amplitude is a function of cavern height and width. We introduce a width-amplitude factor which can be used to calculate the diffraction amplitude of a cavern with a limited width from the diffraction amplitude computed for an infinitely wide cavern.

耦合地震作用下多台阶岩质边坡动力响应分析

运用大型有限元软件ADINA对某重力坝库岸边坡,在场地波地震作用下,模拟单向地震无锚杆支护、单向地震有锚杆支护、耦合地震无锚杆支护及耦合地震有锚杆支护4种工况,从加速度、位移、应力、锚杆轴力及抗滑稳定安全系数几方面进行了动力响应分析。结果表明:锚杆在一定范围内对潜在滑移岩体有很好的嵌固作用,并能有效抑制边坡岩体对加速度的放大效应;多台阶边坡各台阶面和剪出点会产生局部拉应力,锚杆支护后可有效削减;对比单向地震,耦合地震作用下的边坡更危险,在边坡动力分析中应计入竖向地震的作用。研究结果可为同类工程抗震设计提供参考。

Direct pre-stack depth migration on rugged topography

Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.

A study of damping factors in perfectly matched layers for the numerical simulation of seismic waves

When simulating seismic wave propagation in free space, it is essential to introduce absorbing boundary conditions to eliminate reflections from artificially trtmcated boundaries. In this paper, a damping factor referred to as the Gaussian dmping factor is proposed. The Gaussian damping factor is based on the idea of perfectly matched layers （PMLs）. This work presents a detailed analysis of the theoretical foundations and advantages of the Gaussian damping factor. Additionally, numerical experiments for the simulation of seismic waves are presented based on two numerical models： a homogeneous model and a multi-layer model. The results show that the proposed factor works better. The Gaussian damping factor achieves a higher Signal-to-Noise Ratio （SNR） than previously used factors when using same number of PMLs, and requires less PMLs than other methods to achieve an identical SNR.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

Characteristics analysis on high density spatial sampling seismic data

China＇s continental deposition basins are characterized by complex geological structures and various reservoir lithologies. Therefore, high precision exploration methods are needed. High density spatial sampling is a new technology to increase the accuracy of seismic exploration. We briefly discuss point source and receiver technology, analyze the high density spatial sampling in situ method, introduce the symmetric sampling principles presented by Gijs J. O. Vermeer, and discuss high density spatial sampling technology from the point of view of wave field continuity. We emphasize the analysis of the high density spatial sampling characteristics, including the high density first break advantages for investigation of near surface structure, improving static correction precision, the use of dense receiver spacing at short offsets to increase the effective coverage at shallow depth, and the accuracy of reflection imaging. Coherent noise is not aliased and the noise analysis precision and suppression increases as a result. High density spatial sampling enhances wave field continuity and the accuracy of various mathematical transforms, which benefits wave field separation. Finally, we point out that the difficult part of high density spatial sampling technology is the data processing. More research needs to be done on the methods of analyzing and processing huge amounts of seismic data.

Radial tidal current field in a semi-enclosed rectangular basin: formation and evolution

The radial tidal current field accounts for the formation of the radial sand ridges in the South Yellow Sea. Understanding the formation and evolution of this radial tidal current field is vital to assessing the morphodynamic features in the area. A semi-enclosed rectangular basin with and without a coastal barrier was schematized from the topography of the Bohai Sea and Yellow Sea. The 2D tidal current field in this basin was simulated using the DELFT3D-FLOW model. The concept of tidal wave refraction, which highlights the effect of the sloped or stepped submarine topography on the propagation of the tidal waves, was introduced to explain the formation of the radial tidal current field. Under the effect of tidal wave refraction, co-phase lines of the counterclockwise rotating tidal wave and incident tidal wave are transformed into clockwise and counterclockwise deflections, respectively, leading to the convergence and divergence of the flow field. Regardless of whether a coastal barrier exists or not, the outer radial tidal current field might emerge over certain topography. The responses of the radial tidal current field in this basin to the environmental variations such as coastline changes and bottom erosions were discussed. Results show that local protrusion near the focal point of the radial tidal current field will have limited effects on the location of the tidal system. However, a remarkable shift of the amphidromic point toward the entrance and central axis of this basin and a movement of the focal point of the radial tidal current field toward the entrance could be caused by the significant seaward coastline advance and submarine slope erosion.

Three-dimensional transformation of magnetization direction and magnetic field component at low latitudes based on vertical relationship

The transformation of the magnetization direction and the magnetic fi eld component is one of the important methods in magnetic data processing and transformation,which can be conducted in both wavenumber and spatial domains.The transformation method in the wavenumber domain has simpler processing expression and higher processing effi ciency than in the spatial domain;however,they are unstable at low latitude.In this paper,the conclusion that the sum is 0 of two vertical magnetic fi eld components(magnetization inclinations are also perpendicular)in 2D is used for the 3D transformation of the magnetization direction and the magnetic field component.In addition,the transformation method at low latitudes based on vertical relationship(VMT)is proposed,which is an iterative algorithm that converts the transformation of the magnetization direction and the magnetic field component at the low latitude into the high latitude.This method restrains the instability of transformation of constant and variable magnetization direction and magnetic fi eld components in low latitudes.The accuracy,stability,and practicality are verifi ed from synthetic models and real data.

The Focal Mechanism and Characteristic of Stress Field in the Zipingpu Area before the M_S8.0 Wenchuan Earthquake

Based on seismic wave records of the Chengdu digital seismic network and Zipingpu reservoir digital seismic network from August 16,2004 to May 12,2008 when the Wenchuan MS8.0 earthquake occurred,the parameters of focal mechanisms of 486 earthquakes with magnitude larger than M1.6 in the Zipingpu reservoir area were calculated using amplitude ratio method,meanwhile the temporal-spatial variation characteristics of mechanisms and stress field were analyzed based on these parameters.Results show a low ratio of thrust earthquakes and an increased number of strike-slip earthquakes in the reservoir water area in the period from 2006 to 2008.While in the areas far from the reservoir waterfront,the thrust earthquakes took up a high proportion and the strike-slip ones did not increase.The direction of mean principal compressive stress field was deflected and disturbed differently in each area before the Wenchuan MS8.0 earthquake.

Discrete Wavelet Multi-scale Decomposition of the Temporal Gravity Variations in North China

On the basis of the absolute and relative gravity observations in North China,spatial dynamic variation of regional gravity fields is obtained. A multi-scale decomposition technique is used to separate anomalies at different depths,and give some explanation to gravity variation at different time space scales. Gravity variation trends in North China are improved. Based on this result and the analysis of wavelet power spectrum,the images of the depth of wavelet approximation and detail are obtained. The results obtained are of scientific significance for the deep understanding of potential seismic risk in North China from gravity variations in different time space scales.

Prices of Apartments in Relation to Noise Level in Poland

Acoustic climate of a given area ought to be a factor of considerable significance in investment processes in an urbanized area, especially in a residential real estate market, due to its extensive influence on the living standards of its inhabitants. In the following article, the authors have given an analysis of the residential market of housing units located in areas of acceptable and excessive noise levels in preselected regions of Poland. With this end in view, an entirely new source of information has been used in the research--an acoustic map which has been defined and applied to produce the outcome of the analysis. It allowed for the recognition of whether or not the noise level influences decisions made by investors existing in a local residential real estate market.

Physical random function model of ground motions for engineering purposes

A physical random function model of ground motions for engineering purposes is presented with verification of sample level. Firstly,we derive the Fourier spectral transfer form of the solution to the definition problem,which describes the one-dimensional seismic wave field. Then based on the special models of the source,path and local site,the physical random function model of ground motions is obtained whose physical parameters are random variables. The superposition method of narrow-band harmonic wave groups is improved to synthesize ground motion samples. Finally,an application of this model to simulate ground motion records in 1995 Kobe earthquake is described. The resulting accelerograms have the frequencydomain and non-stationary characteristics that are in full agreement with the realistic ground motion records.

SH-wavefield simulation for a laterally heterogeneous whole-Earth model using the pseudospectral method

We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.