人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.

火成岩储层裂缝识别、预测理论与方法研究－以济阳坳陷火成岩油藏为例

Since 1970s, igneous reservoirs such as Shang741, Bin674 and Luol51 have been found in Jiyang depression, which are enrichment and heavy-producing. Showing good prospect of exploration and development, igneous reservoirs have been the main part of increasing reserves and production in Shengli oilfield. As fracture igneous reservoir being an extraordinary complex concealed reservoir and showing heavy heterogeneity in spatial distribution, the study of recognition, prediction, formation mechanism and the law of distribution of fracture is essential to develop the reservoir. Guided by multiple discipline theory such as sedimentology, geophysics, mineralogy, petroleum geology, structural geology and reservoir engineering, a set of theories and methods of recognition and prediction of fractured igneous rock reservoir are formed in this paper. Rock data, three-dimensional seismic data, log data, borehole log data, testing data and production data are combined in these methods by the means of computer. Based on the research of igneous rock petrography and reservoir formation mechanism, emphasized on the assessment and forecast of igneous rock reservoir, aimed at establishing a nonhomogeneity quantification model of fractured igneous rock reservoir, the creativity on the fracture recognition, prediction and formation mechanism are achieved. The research result is applied to Jiyang depression, suggestion of exploration and development for fractured igneous rock reservoir is supplied and some great achievement and favourable economic effect are achieved. The main achievements are gained as follows: 1. The main facies models of igneous rock reservoir in JiYang depression are summarized. Based on data and techniques of seism, well log and logging,started from the research of single well rock facies, proceeded by seismic and log facies research, from point to line and line to face, the regional igneous facies models are established. And hypabyssal intrusion allgovite facies model, explosion volcaniclastic rock facies model and overfall basaltic rocks facies model are the main facies models of igneous rock reservoir in JiYang depression. 2. Four nonhomogenous reservoir models of igneous reservoirs are established, which is the base of fracture prediction and recognition. According to characteristics of igneous petrology and spatial types of reservoir, igneous reservoirs of Jiyang depression are divided into four categories: fractured irruptive rock reservoir, fracture-pore thermocontact metamorphic rock and irruptive rock compound reservoir, pore volcanic debris cone reservoir and fracture-pore overfall basaltic rock reservoir. The spatial distribution of each model's reservoir has its features. And reservoirs can be divided into primary ones and secondary ones, whose mechanism of formation and laws of distribution are studied in this paper. 3. Eight geologic factors which dominate igneous reservoirs are presented. The eight geologic factors which dominates igneous reservoirs are igneous facies, epigenetic tectonics deformation, fracture motion, intensity of intrusive effect and adjoining-rock characters, thermo-contact metamorphic rock facies, specific volcano-tectonic position, magmatic cyclicity and epigenetic diagenetic evolution. The interaction of the eight factors forms the four types nonhomogenous reservoir models of igneous reservoirs in Jiyang depression. And igneous facies and fracture motion are the most important and primary factors. 4. Identification patterns of seismic, well log and logging facies of igneous rocks are established. Igneous rocks of Jiyang depression show typical reflecting features on seismic profile. Tabular reflection seismic facies, arc reflection seismic facies and hummocky or mushroom reflection seismic facies are the three main facies. Logging response features of basic basalt and diabase are shown as typical "three low and two high", which means low natural gamma value, low interval transit-time, low neutron porosity, high resistivity and high density. Volcaniclastic rocks show "two high and three low"-high neutron porosity, high interval transit-time, low density, low-resistance and low natural gamma value. Thermo-contact metamorphic rocks surrounding to diabase show "four high and two low" on log data, which is high natural gamma value, high self-potential anomaly, high neutron porosity, high interval transit-time and low density and low-resistance. Based on seismic, well log and logging data, spatial shape of Shang 741 igneous rock is described. 5. The methods of fracture prediction and recognition for fractured igneous reservoir are summarized. Adopting FMI image log and nuclear magnetic resonance log to quantitative analysis of fractured igneous reservoir and according to formation mechanism and shape of fracture, various fractures are recognized, such as high-angle fracture, low-angle fracture, vertical fracture, reticulated fracture, induced fracture, infilling fracture and corrosion vug. Shang 741 intrusive rock reservoir can be divided into pore-vug compound type, pore fracture type, micro-pore and micro-fracture type. Physical properties parameters of the reservoir are computed and single-well fracture model and reservoir parameters model are established. 6. Various comprehensive methods of fracture prediction and recognition for fractured igneous reservoir are put forward. Adopting three-element (igneous facies, fracture motion and rock bending) geologic comprehensive reservoir evaluation technique and deep-shallow unconventional laterolog constrained inversion technique, lateral prediction of fractured reservoir such as Shang 741 is taken and nonhomogeneity quantification models of reservoirs are established.

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.

Interferometric Measurement of Density in Nonstationary Shock Wave Reflection Flow and Comparison with CFD

We present density measurements from the application of interferometry and Fourier transform fringe analysis to the problem of nonstationary shock wave reflection over a semicircular cylinder and compare our experimental measurements to theoretical results from a CFD simulation of the same problem. The experimental results demonstrate our ability to resolve detailed structure in this complex shock wave reflection problem, allowing visualization of multiple shocks in the vicinity of the triple point, plus visualization of the shear layer and an associated vortical structure. Comparison between CFD and experiment show significant discrepancies with experiment producing a double Mach Reflection when CFD predicts a transitional Mach reflection.

Experimental Study of the Earthquake Recurrence Period and the Trend of Post-Seismic Development

In order to study the earthquake recurrence and the characteristics of earthquake series, rupture tests of rock samples and plexiglass samples were made. On rock samples, a number of acoustic emission (AE) and strain measuring points were deployed; the load was one side direct shear. The variation characteristics of AE and strain at different detecting points around the extra large fracture were observed and studied. On plexiglass samples, a series of inclined cracks were prefabricated by a small-scale compressive testing machine. The samples were then loaded on a shockproof platen, when the samples were loaded, the stress intensity factor (SIF) was determined by the laser interferometric technique and shadow optical method of caustics. The fracture conditions such as material toughness around the extra large fracture were also studied. From those experimental results and the theory of fracture mechanics, the earthquake recurrence period and the trend of post-seismic development were studied.

Biosensor with Total Internal reflection Imaging Ellipsometry

In order to monitor multiple protein reaction processes simultaneously, a biosensor based on imaging ellipsometry operated in the total internal reflection mode is proposed. It could be realised as an automatic analysis for protein interaction processes with real-time label-free method. Its principle and methodology as well as a demonstration for its applications are presented.

The Mach reflection of a detonation based on soot track measurements

This paper presents a series of soot tracks formed by gaseous detonation waves diffracting around wedges with different wedge angles. These cellular structure patterns describe the Mach-reflection processes of a detonation and reveal some unique characteristics. They can be used to analyze the relationship between the trajectory angle of the triple point, wedge angle, and initial pressure in Mach reflection. Compared to the Mach-reflected one-dimensional shock wave in nonreactive air, all these unique characteristics for a Mach-reflected detonation should be attributed to the transverse-wave structure of the detonation front; meanwhile, the precursor shock wave and transverse wave influence the Mach-reflected detonation, respectively. The experimental results support the recently published numerical simulation of this complex phenomenon.

Wave Dynamic Processes in Cellular Detonation Reflection from Wedges

When the cell width of the incident detonation wave (IDW) is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H (2)/O (2) diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes.

Hydrothermal Wave in a Shallow Liquid Layer

The oscillatory thermocapillary convection and hydrothermal wave in a shallow liquid layer, where a temperature difference is applied between two parallel sidewalls, have been numerically investigated in a two-dimensional model. The oscillatory thermocapillary convection and hydrothermal wave appear if the Marangoni number is larger than a critical value. The critical phase speed and critical wave number of the hydrothermal wave agree with the ones given analytically by Smith and Davis in the microgravity environment, and it travels in the direction opposed to the surface flow. Another wave traveled downstream in addition to the hydrothermal wave traveled upstream was observed in the case of earth gravity condition.

Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

Shock wave diffraction and reflection around a dusty square cavity

The diffraction and reflection of planar shock wave around a dusty square cavity is investigated numerically, which is embedded in the net bottom surface of a two-dimensional channel, and the induced gas-particle two-phase now. The wave patterns at different times are obtained for three different values of the particle diameter. The computational results show that the existence of particles affects appreciably the shock wave diffraction and cavity flow.

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.

The Hydrothermal Wave Of Large-Prandtl-Number Fluid In A Shallow Cavity

The hydrothermal wave was investigated numerically for large-Prandtl-number fluid (Pr = 105.6) in a shallow cavity with different heated sidewalls. The traveling wave appears and propagates in the direction opposite to the surface flow (upstream) in the case of zero gravity when the applied temperature difference grows and over the critical value. The phase relationships of the disturbed velocity, temperature and pressure demonstrate that the traveling wave is driven by the disturbed temperature, which is named hydrothermal wave. The hydrothermal wave is so weak that the oscillatory flow field and temperature distribution can hardly be observed in the liquid layer. The exciting mechanism of hydrothermal wave is analyzed and discussed in the present paper.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.