基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

Structures on sliding base subject to horizontal and vertical motions

In this paper particular investigation is directed towards the combined effects of horizontal and vertical motions of real earthquakes to structures resting on sliding base. A simplified method is presented to treat the nonlinear effects of time dependent frictional force of the sliding base as a function of the vertical reaction produced by the foundation. As an example, the El Centro 1940 earthquake record is used on a structural model to show the structural responses due to a sliding base with different frictional and stiffness characteristics. The study shows that vertical ground motion does affect both the superstructure response and the base sliding displacement. Nevertheless, the sliding base isolator is shown to be effective for the reduction of seismic response of a superstructure.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

Hill Instability Analysis Of Tlp Tether Subjected To Combined Platform Surge And Heave Motions

This paper presents the Hill instability analysis of Tension Leg Platform (TLP) tether it, deep sea. The 2-D nonlinear beam model which is Undergoing Coupled axial and transverse vibrations, is applied. The governing equations are reduced to nonlinear Hill equation by use of the Galerkin's method and the modes superposition principle. The Hill instability charted Lip to large parameters is obtained. An important parameter M is defined and can he expressed as the functions of tether length, the platform surge and heave motion amplitudes. Some example studies are performed for various environmental conditions. The results demonstrate that the nonlinear coupling between the axial and transverse vibrations has a significant effect on the response of structure.. It needs to be considered for the accurate dynamic analysis of long TLP tether subjected to the combined platform surge and heave motions.

Soil microbial community analysis using denaturing gradient gel electrophoresis

The most biological diversity on this planet is probably harbored in soils. Understanding the diversity and function of the microbiological component of soil poses great challenges that are being overcome by the application of molecular biological approaches. This review covers one of many approaches being used: separation of polymerase chain reaction (PCR) amplicons using denaturing gradient gel electrophoresis (DGGE). Extraction of nucleic acids directly from soils allows the examination of a community without the limitation posed by cultivation. Polymerase chain reaction provides a means to increase the numbers of a target for its detection on gels. Using the rRNA genes as a target for PCR provides phylogenetic information on populations comprising communities. Fingerprints produced by this method have allowed spatial and temporal comparisons of soil communities within and between locations or among treatments. Numerous samples can be compared because of the rapid high throughput nature of this method. Scientists now have the means to begin addressing complex ecological questions about the spatial, temporal, and nutritional interactions faced by microbes in the soil environment.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

大别山岩石圈有效弹性厚度初步研究

As powerful tools to study the lithosphere dynamics, the effective elastic thickness (Te) as well as the envelope of yielding stress of lithosphere have been attracted great attention of geoscientists in the past thirty years. The oceanic lithosphere, contrary to the continental lithosphere, has more fruits for its simple structures and evolution process. In continent, the lithosphere commonly is complex and variable in the rheological, thermal structures, and has a complicated history. Therefore, the application of the effective elastic thickness in continent is still a subject to learn in a long time. Te, with the definition of the thickness of an elastic plate in theory flexured by the equal benging of the real stress in the lithosphere plate (Turcotte, 1982), marks the depth of transition between elastic and fluid behaviors of rocks subjected to stress exceeding 100 MPa over the geological timescales (McNutt, 1990). There are three methods often adapted: admittance or isostatic response function, coherence and forwarding. In principle, the models of Te consist of thermal-rheological, non-linear Maxwell, non-linear work hardening and rheological layered models. There is a tentative knowledge of Te that it is affected by the following factors: crustal thickness, crust-mantle decoupling, plate bending, boundary conditions of plate (end forces and bending moments), stress state, sedimentary layer, faulting effect, variation in the mountain belts' strike, foreland basin, inheritance of tectonic evolution, convection of mantle, seismic depth and lithosphere strength. In this thesis, the author introduces the geological sketch of the Dabie collisional orogenic belt and the Hefei Basin. The Dabie Mts. is famous for the ultra-high pressure metamorphism. The crustal materials subducted down to the depth of at least 100 km and exhumed. So that the front subjects arise such as the deeply subduction of continent, and the post-collisional crust-mantle interaction. In a geological journey at June of 1999, the author found the rarely variolitic basaltic andesite in the Dabie Mts. It occurs in Susong Group, near Zhifenghe Countryside, Susong County, Anhui Province. It is just to the south of the boundary between the high-grade Susong melange and the ultra-high grade South Dabie melange. It has a noticeable knobby or pitted appearance in the surface. The size of the varioles is about 1-4 mm. In hand-specimen and under microscope, there are distinct contacts between the varioles and the matrice. The mineralogy of the varioles is primarily radiate plagioclase, with little pyroxene, hornblende and quartz. The pyroxene, hornblende and quartz are in the interstices between plagioclase. The matrix is consisted of glass, and micro-crystals of chlorite, epidote and zoisite. It is clearly subjected and extensive alteration. The andesite has an uncommon chemical composition. The SiO_2 content is about 56.8%, TiO_2 = 0.9%, MgO = 6.4%, (Fe_2O_3)_(Total) = 6.7% ～ 7.6%, 100 Mg/(Mg+Fe) = 64.1 ～ 66.2. Mg# is significantly high. The andesite has higher abundances of large-lithophile trace elements (e.g. K, Ba, Sr, LREE), e.g. La/Nd = 5.56-6.07, low abundances of high-strength-field elements (HFSE, e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subducted-related magmas (Pearcce, 1982; Sun and McDonaugh, 1989). In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a characteristic of the continental margins (Pearce, 1982). There has a strongly enrichment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)_N is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)_N = 28.63-36.74, (La/Y)_N = 70.33 - 82.84. The elements Y and Yb depleted greatly, Y < 20 ppm, Y_N = 2.74-2.84, Yb_N = 2.18 - 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is -18.7 ～ -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the basaltic andesite was resulted from the post-collisional crust-mantle interaction between the subducted Yangze carton and the mantle of Sino-Korea carton. To obtain the Te of the lithosphere beneath the Dabie Mts. and the Hefei Basin, the author applies the coherence method in this thesis. The author makes two topography-gravity profiles (profiles 7~(th) and 9~(th)) across the Dabie Mts. and the Hefei Basin, and calculates the auto-coherence, across coherence, power spectrum, across power spectrum of the topography and gravity of the two profiles. From the relationships between the coherence and the wave-number of profiles. From the relationships between the coherence and the wave-number of profiles 7~(th) and 9~(th), it is obtained that the characteristic wavelengths respectively are 157 km and 126 km. Consequently the values of effective elastic thickness are 6.5 km and 4.8 km, respectively. However, the Te values merely are the minimum value of the lithosphere because the coherencemethod in a relative small region will generate a systemic underestimation. Why there is a so low Te value? In order to check the strength of the lithosphere beneath the Dabie Mts., the authore tries to outline the yielding-stress envelope of the lithosphere. It is suggested that the elastic layers in the crust and upper mantle are 18 km and 35 km, respectively. Since there exist a low viscosity layer about 3-5 km thickness, so it is reasonable that the decoupling between the crust and mantle occurred. So the effective thickness of the lithosphere can be estimated from the two elastic layers. Te is about 34 km. This is the maximum strength of the lithosphere. We can make an approximately estimation about the strength of the lithosphere beneath the Dabie Mts.: Te is about 20-30 km. The author believes that the following factors should be responsible for the low Te value: (1) the Dabie Mts. has elevated strongly since K_3-J_1. The north part of the Dabie Mts. elevates faster than the south part today; (2) there occur large active striking faults in this area. And in the east, the huge Tan-Lu striking fault anyway tends to decrease the lithosphere strength; (3) the lithosphere beneath the Dabie Mts. is heter-homogeneous in spatio-temporal; (4) the study area just locates in the adjacent region between the eastern China where the lithosphere thickness is significantly reduced and the normal western China. These factors will decrease the lithosphere strength.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Effect of desert soil algae on the stabilization of fine sands

Four filamentous cyanobacteria, Microcoleus vaginatus, Phormidium tenue, Scytonema javanicum (Kutz.) and Nostoc sp., and a single-celled green alga, Desmococcus olivaceus, all isolated from Shapotou (Ningxia Hui Autonomous Region of China), were batch cultured and inoculated onto unconsolidated sand in greenhouse and field experiments. Their ability to reduce wind erosion in sands was quantified by using a wind tunnel laboratory. The major factors related to cohesion of algal crusts, such as biomass, species, species combinations, bioactivity, niche, growth phase of algae, moisture, thickness of the crusts, dust accretion (including dust content and manner of dust added) and other cryptogams (lichens, fungi and mosses) were studied. The best of the five species were M. vaginatus and P. tenue, while the best mix was a blend of 80% M. vaginatus and 5% each of P. tenue, S. javanicum, Nostoc sp. and D. olivaceus. The threshold friction velocity was significantly increased by the presence of all of the cyanobacterial species, while the threshold impact velocity was notably increased only by the filamentous species. Thick crusts were less easily eroded than thin crusts, while biomass was more effective than thickness. Dust was incorporated best into Microcoleus crust when added in small amounts over time, and appeared to increase growth of the cyanobacterium as well as strengthen the cohesion of the crust. Microbial crust cohesion was mainly attributed to algal aggregation, while lichens, fungi and mosses affected more the soil structure and physico-chemical properties.