首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 234 毫秒
1.
A new numerical approach to the solution of waves propagating in a fluid-saturated medium, using Biot's theory as a foundation, has important implications for oil reservoir management and earthquake prediction. A numerical scheme is developed using an exponential transformation that explicitly treats the petrophysical and fluid properties of the medium within the framework of a generalized model. The scheme accounts for wave dissipation and velocity modifications. The numerical solution is used to perform numerical experiments to study the dynamic behaviour of waves in a fluid-saturated medium at well-logging frequencies (15 kHz). The results from the numerical experiments indicate that the degree of saturation by a high-viscosity fluid (HVF) such as oil, the temperature and the porosity of a medium strongly influence the spectral power distribution, frequency content and the velocity of waves propagating through the medium. An increase in HVF saturation causes enhanced attenuation of the low-frequency components, and increases the seismic velocity. An increase in porosity, however, enriches the low-frequency components and decreases the seismic velocity. A spectral quantification procedure is suggested and used to obtain information about the petrophysical and fluid properties of the medium from the spectral characteristics of the transmitted waveform. The procedure involves segmentation of the energy or power distribution of the transmitted waveforms into specified energy bands. The energy or power in these bands is then estimated. The extracted quantification variables are found to have strong correlations with the degree of HVF saturation, and the temperature and the porosity of the medium.  相似文献   

2.
The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.  相似文献   

3.
Summary. Using a single scattering approximation, we derive equations for the scattering attenuation coefficients of P- and S -body waves. We discuss our results in the light of some recent energy renormalization approaches to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P - to S -waves and S- to P -waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average wavenumber power spectra. We approximate the power spectra with piecewise constant functions, each segment of which contributes to the net, frequency-dependent, scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behaviour associated with each spectral segment, including a 'transition' peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observed such a peak in the apparent attenuation of seismic waves. We give both the frequency and distance limits on the accuracy of the theoretical results.  相似文献   

4.
Effects of fractures on seismic-wave velocity and attenuation   总被引：1，自引：0，他引：1  
The effects of fractures on the seismic velocity and attenuation of a rock are investigated using theoretical results and experimental data. Fractures in a rock mass influence the traveltimes and amplitudes of seismic waves that have propagated through them. The displacement discontinuity model, recently employed in fracture investigations, is modified to describe the effect of fractures on seismic-wave velocity and attenuation. This new model, the modified displacement discontinuity model (MDD), is formulated in a way analogous to transmission-line analysis. The fractures are treated as transmission lines for the passage of seismic waves. The MDD takes into consideration realistic fracture parameters which include the fracture length, the fractional area of a fracture surface in contact, and the nature of the infilling material. A single fracture of varying geometric and material properties is shown to affect dramatically the transmission properties of a propagating waveform, and hence the seismic velocity and attenuation. These effects have been shown to result in a frequency-dependent velocity and attenuation. The sensitivity of the fracture parameters to seismic-wave velocity and attenuation was investigated and interesting results were obtained. Fracture parameters used in designing experimental models consisting of synthetically manufactured cracks were fed into the MDD and a well-known crack model, Hudson's model, for comparison. Velocities as a function of the incident-wave angle were obtained from both numerical models and were compared with the results from the experimental modelling. For P waves, the MDD model results show better agreement with those of the experimental model for all crack densities investigated than those from Hudson's model.  相似文献   

5.
Seismic amplitude tomography for crustal attenuation beneath China   总被引：1，自引：0，他引：1  
Amplitude tomography reconstructs seismic attenuation directly from recorded wave amplitudes. We have applied the tomography to amplitude data reported in the 'Annual Bulletin of Chinese Earthquakes' and interpreted the regionally varying crustal attenuation in terms of tectonics. The seismic amplitudes were originally recorded for determining the M L and M S magnitudes. They generally correspond to the maximum amplitudes of the horizontal components of the short-period S waves and intermediate-period Rayleigh waves. Both sets of measurements are sensitive to crustal structure. The peak amplitudes from M L amplitudes spread spherically with significant dispersion and scattering. M S amplitudes show cylindrical spreading with little dispersion. Average crustal Q values for attenuation at 1 Hz are 737 and 505 for M L and M S, respectively, with substantial regional variations. Frequency dependence in the attenuation is also indicated. Regions with the lowest attenuation (high Q values) are beneath the south China Block, Sichuan Basin, Ordos Platform, the Daxinganling and the Korea Craton. These tend to be tectonically inactive regions, which are generally dominated by intrusive and cratonic rocks in the upper crust. Regions with the highest attenuation (low Q values) are beneath Bohai Basin, Yunnan, eastern Songpan-Ganzi Terrain, margins of the Ordos platform and the Qilian Shan. These are predominantly active basins, grabens and fold belts. The continental margin also highly attenuates both S and surface waves.  相似文献   

6.
Summary. Records of 21 earthquakes with ray paths crossing south-eastern China have been examined to establish the nature of propagation of the seismic wave Sn . Sn is seen at stations within 1500km of earthquakes, but not beyond 1500km. No single localized region of attenuation can account for all the missing or weak Sn waves. The inefficient propagation of Sn through south-eastern China is surprising in light of efficient Sn propagation in other continental platforms which, like south-eastern China, appear to have been stable throughout the Cenozoic.  相似文献   

7.
《Basin Research》2018,30(5):895-925
Kilometre‐scale geobodies of diagenetic origin have been documented for the first time in a high‐resolution 3D seismic survey of the Upper Cretaceous chalks of the Danish Central Graben, North Sea Basin. Based on detailed geochemical, petrographic and petrophysical analyses, it is demonstrated that the geobodies are of an open‐system diagenetic origin caused by ascending basin fluids guided by faults and stratigraphic heterogeneities. Increased amounts of porosity‐occluding cementation, contact cement and/or high‐density/high‐velocity minerals caused an impedance contrast that can be mapped in seismic data, and represent a hitherto unrecognized, third type of heterogeneity in the chalk deposits in addition to the well‐known sedimentological and structural features. The distribution of the diagenetic geobodies is controlled by porosity/permeability contrasts of stratigraphic origin, such as hardgrounds associated with formation tops, and the feeder fault systems. One of these, the Top Campanian Unconformity at the top of the Gorm Formation, is particularly effective and created a basin‐wide barrier separating low‐porosity chalk below from high‐porosity chalk above (a Regional Porosity Marker, RPM). It is in particular in this upper high‐porosity unit (Tor and Ekofisk Formations) that the diagenetic geobodies occur, delineated by “Stratigraphy Cross‐cutting Reflectors” (SCRs) of which eight different types have been distinguished. The geobodies have been interpreted as the result of: (i) escaping pore fluids due to top seal failure, followed by local mechanical compaction of high‐porous chalks, paired with (ii) ascension of basinal diagenetic fluids along fault systems that locally triggered cementation of calcite and dolomite within the chalk, causing increased contact cements and/or reducing porosity. The migration pathway of the fluids is marked by the SCRs, which are the outlines of high‐density bodies of chalk nested in highly porous chalks. This study, thus, provides new insights into the 3D relationship between fault systems, fluid migration and diagenesis in chalks and has important applications for basin modelling and reservoir characterization.  相似文献   

8.
Ultrasonic detection technology is of great significance in the detection and evaluation of physical and mechanical properties of frozen soil, but wave propagation characteristics in frozen soil are unclear. Based on the three-phase composition of frozen saturated soil and the mixture theory, considering Bishop's effective stress formula, the wave propagation equations are establish for frozen saturated soil. In wave propagation, an entropy inequality was introduced to describe the coupling of different phases. The analytic expressions of propagation velocity and attenuation law of waves in frozen soil are obtained, and wave propagation characteristics in frozen saturated soil are discussed. Results show that four types of waves(i.e., P1, P2, P3 and S) are found in frozen saturated soil and all four wave types are dissipative waves, in which the attenuation of P3 is the maximum. The velocity of four waves increases sharply at the excitation frequency range of 10~3–10~9 Hz,but the wave velocity at high-frequency and low-frequency is almost constant. When volume ice content increases, the wave propagation velocity of P1 and S decreases dramatically, and the velocity of P2 increases gradually, but P3 velocity increases first and then decreases to zero with increasing saturation. The attenuation coefficients of P1 and S waves begins to increase gradually when the volume ice content is about 0.4, P2 increases first and then decreases with an increase of volume ice content and P3 increases with the volume ice content and decreases rapidly from extreme to zero.  相似文献   

9.
A general one-way representation of seismic data can be obtained by substituting a Green's one-way wavefield matrix into a reciprocity theorem of the convolution type for one-way wavefields. From this general one-way representation, several special cases can be derived.
By introducing a Green's one-way wavefield matrix for primaries , a generalized Bremmer series representation is obtained. Terminating this series after the first-order term yields a primary representation of seismic reflection data. According to this representation, primary seismic reflection data are proportional to a reflection operator, 'modified' by primary propagators for downgoing and upgoing waves. For seismic imaging, these propagators need to be inverted. Stable inverse primary propagators can easily be obtained from a one-way reciprocity theorem of the correlation type.
By introducing a Green's one-way wavefield matrix for generalized primaries , an alternative representation is obtained in which multiple scattering is organized quite differently (in comparison with the generalized Bremmer series representation). According to the generalized primary representation, full seismic reflection data are proportional to a reflection operator, 'modified' by generalized primary propagators for downgoing and upgoing waves. Internal multiple scattering is fully included in the generalized primary propagators {either via a series expansion or in a parametrized way). Stable inverse generalized primary propagators can be obtained from the one-way reciprocity theorem of the correlation type. These inverse propagators are the nucleus for seismic imaging techniques that take the angle-dependent dispersion effects due to fine-layering into account.  相似文献   

10.
Summary. Theoretical developments of Hudson demonstrate how to calculate the variations of velocity and attenuation of seismic waves propagating through solids containing aligned cracks. The analysis can handle a wide variety of crack configurations and crack geometries. Hudson associates the velocity variations with effective elastic constants. In this paper we associate the variation of attenuation with the imaginary parts of complex effective elastic constants. These complex elastic constants permit the simulation of wave propagation through two-phase materials by the calculation of wave propagation through homogeneous anisotropic solids.  相似文献   

11.
Numerical simulation of the propagation of P waves in fractured media   总被引：1，自引：0，他引：1  
We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q −1 of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q −1 of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) −1 in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.  相似文献   

12.
i
Displacements of Love waves generated by a two-dimensional point source in a layered medium have been studied earlier by Sezawa & Sato by the method of successive reflections at the boundaries. In this paper the same problem has been worked out by using Green's function. The paper deals with the study of attenuation of Love waves of low periods in the coastal region. Experimental observations show that Love waves of smaller periods can be obtained only in the island observing stations. A slight intervention of the continental boundary is sufficient to attenuate lower period Love waves giving a hint thereby that attenuation of lower periods takes place perhaps at the continental margin. Taking a simplified configuration for the continental boundary and using Green's function technique, the displacement of Love waves due to a point source has been obtained and it has been shown that attenuation of Love waves of smaller periods takes place in the continental margin due to the slope of the boundary.  相似文献   

13.
  相似文献   

14.
Summary. A fluid-saturated cubic packing of like elastic spheres is taken to be in equilibrium under the effect of gravity and the effects of a superimposed low-frequency elastic wave are considered. In the first place, expressions for the wave velocity, dispersion and attenuation are derived for the dry packing. This dynamic theory leads to the result that, for very low frequencies, the wave velocity is proportional to the third root of the depth and not the sixth root as is obtained by using the effective elastostatic modulus of the packing. For the fluid-saturated packing, two waves, termed respectively the 'solid wave' and the 'fluid wave', are found to propagate. The 'solid wave' has the characteristics of a wave propagating within a dry packing whose parameters differ in a specified way from those of the original packing, whereas the 'fluid wave' has those of a wave within a homogeneous fluid with similarly modified parameters.  相似文献   

15.
  相似文献   

16.
Summary Many accept thermal convection within the mantle of the Earth as the driving mechanism for continental drift. It is also of considerable interest to determine whether thermal convection is occurring within Venus, Mars, and the Moon. In this paper a systematic treatment of the stability of planetary interiors is given. The thermal stability problem for a layer of fluid heated from below is solved when the viscosity of the fluid increases exponentially with depth. For a semi-infinite fluid with exponentially increasing viscosity the critical Rayleigh number based on the surface viscosity and the scale length of the viscosity increase is found to be 30 for a fixed surface boundary condition and 23 for a free surface boundary condition. This stability analysis is also extended to include volume heat release. The thermal stability of the interiors of the Earth, Venus, Mars and the Moon is examined. Using temperature-depth profiles in the literature and a theoretical expression for the viscosity of a crystalline solid based on diffusion creep, viscosity depth profiles for planetary interiors are obtained. Because of the strong pressure effect the viscosity within the Earth and Venus increases greatly from a near surface minimum. For Mars the increase is less pronounced and within the interior of the Moon the viscosity is nearly constant. For all the cases considered the planetary interiors are found to be thermally unstable. Because of the dependence of viscosity on depth the interiors of the Moon and Mars are found to be considerably less stable than Venus and the Earth. It is concluded that thermal convection is occurring within the planetary bodies considered.  相似文献   

17.
A new and efficient method for fault seal analysis using seismic data is presented. It uses multiple seismic attributes and neural networks to enhance fluid migration pathways, including subtle features that are not detectable using single attributes only. The method may be used as a first estimate of fault seal or to calibrate results from other techniques. The results provide information about which faults and fault segments are sealing or leaking. Fluid flow along individual faults appears to be focused along zones of weakness, and fault seal research should thus be focused on finding such weak locations within fault zones, a task that is best done using three‐dimensional (3D) seismic data. Under certain conditions, it is suggested that fluids migrate along fault planes by a diapiric fluid flow mechanism. The results assist in calibrating the bulk hydraulic properties of faults and rock formations and can be used in basin modelling.  相似文献   

18.
Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= tu/Qm where t u is the travel time and Qm the quality factor in the absorption band, and τ−1 m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ m. For very small τ m , this model is shown to be very similar to previously proposed attenuation models. As τ m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.  相似文献   

19.
Summary The displacement response of an elastic half space to a plane pressure wave is examined in order to establish the conditions under which sources of this type can contribute significantly to the long-period seismic noise field. The study is restricted to pressure waves which propagate at velocities well below the seismic wave velocities characteristic of the half space. The numerical studies indicate that pressure waves with amplitudes of 100 μbar or more can contribute significantly to the long-period vertical background noise observed at the surface, provided that the detectors are located on sections of alluvial fill or poorly to moderately indurated sandstones and shales whose thicknesses are greater than about a kilometre. These same waves can also create significant tilt noise on long-period horizontal seismographs located at or near the surface, regardless of the rock type. The seismic disturbances created by pressure waves decay rapidly away from the surface. Therefore, it appears that it may be possible to eliminate the effects of atmospherically generated noise by placing the detectors at moderate depths.  相似文献   

20.
We show that seismic shear waves may be used to monitor the in situ stress state of deep inaccessible rocks in the crust. The most widespread manifestation of the stress-related behaviour of seismic waves is the shear-wave splitting (shear-wave birefringence) observed in almost all rocks, where the polarizations of the leading split shear waves are usually subparallel to the direction of the local maximum horizontal stress. It has been recognized that such shear-wave splitting is typically the result of propagation through distributions of stress-aligned fluid-filled microcracks and pores, known as extensive-dilatancy anisotropy or EDA. This paper provides a quantitative basis for the EDA hypothesis. We model the evolution of anisotropic distributions of microcracks in triaxial differential stress, where the driving mechanism is fluid migration along pressure gradients between neighbouring microcracks and pores at different orientations to the stress field. This leads to a non-linear anisotropic poroelasticity (APE) model for the stress-sensitive behaviour of fluid-saturated microcracked rocks. A companion paper shows that APE modelling matches a range of observed phenomena and is a good approximation to the equation of state of a stressed fluid-saturated rock mass.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号