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.

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.

Geophysical evaluation methods for buried hill reservoirs in the Jiyang superdepression of the Bohai Bay basin, eastern China

The Jiyang superdepression is one of the richest hydrocarbon accumulations in the Bohai Bay basin, eastern China. Comprehensive seismic methods have been used in buried hill exploration in Jiyang to describe these fractured reservoirs better. Accurate seismic stratigraphic demarcation and variable-velocity mapping were applied to reveal the inner structure of the buried hills and determine the nature of the structural traps more precisely. Based on the analysis of rock properties and the characteristics of well-developed buried hill reservoirs, we have successfully linked the geology and seismic response by applying seismic forward technology. Log-constrained inversion, absorption coefficient analysis and tectonic forward-inversion with FMI loggings were applied to analyse and evaluate the buried hill reservoirs and gave satisfying results. The reservoir prediction was successful, which confirmed that the comprehensive utilization of these methods can be helpful in the exploration of buried hill reservoirs.

基于正演方法的井间地震响应特征研究

地震勘探的基本任务是测量由人工震源激发的各种地震波，通过对所得到的地震波场的数据记录进行必要的处理和解释而提取所需要的地质信息。不同的地震波场特征反应了不同的地质构造等信息，因此，对地震波场特征进行研究就显得非常重要。井间地震是在一口井中激发、另一口井中接收的一种特殊的地震勘探方法，其实际记录要比常规的地震记录复杂的多，其波场中的各种波也更丰富复杂，如何识别不同类型的波并研究它们的响应特征是一个非常必要的任务。正演模拟可以为地震数据采集、处理、解释提供理论依据，是井间地震波场特征分析的一种有效方法。 本文主要研究了井间地震的数值正演模拟问题，系统地介绍了射线追踪正演模拟和波动方程正演模拟方法，利用Visual Basic 6.0语言编写了相应的射线追踪正演模拟程序，而且分两种方法详细分析了井间地震数值模型的地震响应特征，为直观的井间地震波场识别和勘探地震数据反演提供了必要的依据。本文取得的主要成果如下： 1)在前人的VSP正演模拟的基础上，改进并总结出了一套井间地震射线追踪正演模型的制作方法。 2)编制了井间地震射线追踪正演模拟程序V1.0。 3)利用编写的井间地震正演模拟程序，作了两种射线追踪数值模拟试算；利用Tesseral声波方程正演模拟试算了含有楔形体和丘状体的理论模型；取得了较好的模拟数据。

裂缝性油气藏储层预测方法研究

In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a series researches are conducted on fracture reservoir prediction technology in general，and it especially focus on some difficult points. The technological series which integrated amplitude preserved data processing、interpretation and its comprehensive application research as a whole were developed and this new method can be applied to the other similar oilfield exploration and development. The contents and results in this paper are listed as follows: 1. An overview was given on the status and development of fracture reservoir estimation technique, compare and analyze those geophysical prediction methods. This will be very helpful to the similar reservoir researches. 2. Analyze and conclude the characters of geologies and well logging response of burial hills fracture reservoir, those conclusions are used to steer the geophysical research and get satisfying results. 3. Forward modeling anisotropy seismic response of fracture reservoir. Quantitatively describe the azimuthal amplitude variation. Amplitude ellipse at each incidence angle is used to identify the fracture orientation. 4. Numerical simulation of structure stress based on finite difference method is carried out. Quantitatively describe and analyze the direction and intensity of fracture. 5. Conventional attributes extraction of amplitude preserved seismic data、attributes with different azimuthal angle and different offset are used to determine the relationship between the results and fracture distribution. 6. With spectrum decomposition method based on wavelet transform， the author disclose the reservoir distribution in space. It is a powerful tool to display its anisotropy. 7. Integrated seismic wave impendence、elastic impendence、spectrum decomposition、attribute extraction、fracture analysis result as a whole to identify and evaluate the fracture reservoir. An optimum workflow is constructed. It is used to practical oil&gas production and good results are obtained. This can indicate the wide foreground of this technique series.

辽东湾隐蔽油气藏地球物理识别技术研究

In China and world, more than half the recent basin discovered reserves involve lithologic hydrocarbon reservoir reserves. The major target for further hydrocarbon basin exploration is the subtle reservoir. The Liaodong Bay prospect is much important in Bohai Sea, which includes Liaoxi low uplift, Liaodong uplift, Liaoxi sag and Liaozhong sag. After dozens years’ exploration in Liaodong Bay, few unexplored big-and-middle-sized favorable structural traps are remained and most of the stock structure targets are bad for fragmentary. Thus seeking for new prospect area and making a breakthrough, have become the unique way to relieve the severe exploration condition in Liaodong Bay. Technique Route Based on the petrophysical property of target area, the seismic forward inference of typical subtle trap model is expanded with analysis of logging, seismic and geologic data. According to petrophysical characteristics and forward inference and research on seismic response of actual seismic data in target area, the optimization of geophysical technique is used in subtle trap identification and the geophysical identification technique system of subtle reservoir is formed. The Key Research ① Petrophysical Model The petrophysical parameter is the basic parameter for seismic wave simulation. The seismic response difference of rocks bearing different fluids is required. With the crossplot of log data, the influence of petrophysical parameters on rock elastic properties of target area is analyzed, such as porosity, shale index, fluid property and saturation. Based on the current research on Biot-Gassmann and Kuster-Toksoz model, the petrophysical parameter calculator program which can be used for fluid substitution is established. ② S-wave evaluation based on conventional log data The shear velocity is needed during forward inference of AVO or other elastic wave field. But most of the recent conventional log data is lack of shear wave. Thus according to the research on petrophysical model, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. ③ AVO forward modeling based on well data For 6 wells in JZ31-6 block and 9 wells in LD22-1 block, the AVO forward modeling recording is made by log curve. The classification of AVO characteristics in objective interval is made by the lithologic information. ④ The 2D parameter model building and forward modeling of subtle hydrocarbon trap in target area. According to the formation interpretation of ESS03D seismic area, the 2D parameter model building and seismic wave field forward modeling are carried on the given and predicted subtle hydrocarbon trap with log curve. ⑤ The lithology and fluid identification of subtle trap in target area After study the seismic response characteristics of lithology and fluid in given target area, the optimization of geophysical technique is used for lithology identification and fluid forecast. ⑥The geophysical identification technique system of subtle reservoir The Innovative Points of this Paper ① Based on laboratory measurement and petrophysical model theory, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. Then the fluid substitution method based on B-G and K-T theory is provided. ② The method and workflow for simulating seismic wave field property of subtle hydrocarbon trap are established based on the petrophysical model building and forward modeling of wave equation. ③ The description of subtle trap structural feature is launched. According to the different reflection of frequency wave field structural attribute, the fluid property of subtle trap can be identified by wave field attenuation attribute and absorption analysis. ④ It’s the first time to identify subtle trap by geophysical technique and provide exploration drilling well location. ⑤ The technique system of subtle reservoir geophysical identification is formed to provide available workflow and research ideas for other region of interest.

大陆地壳岩石结构重建：方法与实例

The composition of the continental crust has long been a subject of interest to earth scientists as it can provide key information about the crustal growth and evolution of the continents. In this paper we make a comparative study on the lithological discrimination schemes featuring with the use of different seismic attributes, such as P-wave velocity, P- to S-wave velocity ratio, acoustic or elastic impedances, Lame impedances and high-sensitive identification factors. The results demonstrate that Lame impedances have more powerful constrains than other seismic attributes. In order to fully take the advantage of make the best of the different seismic response of crustal rock, we firstly use seismic attribute that have weak distinguish power to construct loose constrained lithological model, then use seismic attributes that have stronger distinguish power to tighten the constrains of lithological discrimination. We propose a joint scheme (chain constrain technique) by combing all available constrains to reduce the non-uniqueness in mapping rock distribution. We adopt chain constrain technique to construct lithological model beneath Tunxi-Wenzhou transect, Southeastern China, Manzhouli-Suifenhe transect, Northeastern China, and geophysical profile in Bohai Bay Basin, North China. The results can be suumarized as the follows: (1) We compare the sensitivity of different seismic factor constraints on rock types, and conclude that Lame impedances have tighter constrains than seismic velocity, Vp/Vs, density. (2) We propose chain constrains to construct lithological model from integrated geophysical data, and reduce the non-uniqueness in mapping rock distribution. (3) We reconstruct crustal lithological model beneath Tunxi-Wenzhou transect, Southeastern China. The results suggested that Jiangshan-Shaoxing fault is a crust-scale, and it is the boundary between Cathaysia and Yanthze blocks. (4) We construct crustal lithological model beneath Manzhouli-Suifenhe transect, Northeastern China. (5) We map the petrologic distribution along a geophysical profile in Bohai Bay Basin, North China, and construct a three-layered petrology model from the depth 2 km to about 10 km, consisted of basalt (the first layer), pelitic siltstone (the second layer), and silty mudstone and fine sandstone (the third layer).

AVO分析与地震资料叠前储层预测技术

In this paper we base on the anisotropic theory and Zoeppritz function of the transmission theory and the law of amplitude versus offset simplify seismic reflection coefficient of different media, analyze the characteristic of the gas or oil saturated stratum or the VTI and HTI models. Discuss the P wave reflection relationship and the meanings of the different parameters. We use measured parameters of a reservoir to simulate the characteristic of the reservoir, study the different effects of stratum saturated with gas or oil and analyze the characteristic of the seismic response of different models which change with different incident angles and different azimuths. Using the field data of logs ,analyze the rock property parameters, build the relationship of logs and parameters by Gassmann theory or empirical function. Calculate the density and the shear modulus and bulk modulus, reconstruct the log curves, calculate shear wave logs and correlate the logs affected by mud and other environmental factors. Finally perform the relationship of the seismic data log of saturated stratum and enhance the ability and reliability in reservoir prediction. Our aim is by the prestack seismic processing to get high solution and amplitude preserved seismic data. Because in incident angle gathers or azimuthal gathers, the low signal to noise ratio and low different covers affect the result of the prestack reservoir prediction. We apply prestack noise erase, cell regularization process and relatively amplitude preservation in the high solution seismic process routine to preserve the characteristic of stratum response, and erase the effects of the noise. In this paper we finished prestack invertion in the BYT survey and fractured reservoir depiction in MB survey. By the invertion and multiple attributes crossplot. we can get the stratum profiles and oil indicator profiles which can predict the distribution of the reservoir and oil. In the MB survey, we get orientation and density of fractured reservoir by the azimuthal seismic amplitude and depict the potential oil and gas reservoir. Prestak invertion works better in distinguishing oil and reservoir.

基于广义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.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Material response and failure mechanism of unidirectional metal matrix composites under impulsive shear loading

The material response and failure mechanism of unidirectional metal matrix composite under impulsive shear loading are investigated in this paper. Both experimental and analytical studies were performed. The shear strength of unidirectional C-f/A356.0 composite and A356.0 aluminum alloy at high strain rate were measured with a modified split Hopkinson torsional bar technique. The results indicated that the carbon fibers did not improve the shear strength of aluminum matrix if the fiber orientation aligned with the shear loading axis. The microscopic inspection of the fractured surface showed a multi-scale zigzag feature which implied a complicated shear failure mechanism in the composite. In addition to testing, the micromechanical stress field in the composite was analyzed by the generalized Eshelby equivalent method (GEEM). The influence of cracking in matrix on the micromechanical stress field was investigated as well. The results showed that the stress distribution in the composite is quite nonhomogeneous and very high shear stress concentrations are found in some regions in the matrix. The high shear stress concentration in the matrix induces tensile cracking at 45 degrees to the shear direction. This in turn aggravates the stress concentration at the fiber/matrix interface and finally leads to a catastrophic failure in the composite. From the correlation between the analysis and experimental results, the shear failure mechanism of unidirectional C-f/A356.0 composite can be elucidated qualitatively.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.