首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 125 毫秒
1.
Seismic monitoring of reservoir and overburden performance during subsurface CO2 storage plays a key role in ensuring efficiency and safety. Proper interpretation of monitoring data requires knowledge about the rock physical phenomena occurring in the subsurface formations. This work focuses on rock stiffness and elastic velocity changes of a shale overburden formation caused by both reservoir inflation induced stress changes and leakage of CO2 into the overburden. In laboratory experiments, Pierre shale I core plugs were loaded along the stress path representative for the in situ stress changes experienced by caprock during reservoir inflation. Tests were carried out in a triaxial compaction cell combining three measurement techniques and permitting for determination of (i) ultrasonic velocities, (ii) quasistatic rock deformations, and (iii) dynamic elastic stiffness at seismic frequencies within a single test, which allowed to quantify effects of seismic dispersion. In addition, fluid substitution effects connected with possible CO2 leakage into the caprock formation were modelled by the modified anisotropic Gassmann model. Results of this work indicate that (i) stress sensitivity of Pierre shale I is frequency dependent; (ii) reservoir inflation leads to the increase of the overburden Young's modulus and Poisson's ratio; (iii) in situ stress changes mostly affect the P‐wave velocities; (iv) small leakage of the CO2 into the overburden may lead to the velocity changes, which are comparable with one associated with geomechanical influence; (v) non‐elastic effects increase stress sensitivity of an acoustic waves; (iv) and both geomechanical and fluid substitution effects would create significant time shifts, which should be detectable by time‐lapse seismic.  相似文献   

2.
Acoustic emissions (AE), compressional (P), shear (S) wave velocities, and volumetric strain of Etna basalt and Aue granite were measured simultaneously during triaxial compression tests. Deformation-induced AE activity and velocity changes were monitored using twelve P-wave sensors and eight orthogonally polarized S-wave piezoelectric sensors; volumetric strain was measured using two pairs of orthogonal strain gages glued directly to the rock surface. P-wave velocity in basalt is about 3 km/s at atmospheric pressure, but increases by > 50% when the hydrostatic pressure is increased to 120 MPa. In granite samples initial P-wave velocity is 5 km/s and increases with pressure by < 20%. The pressure-induced changes of elastic wave speed indicate dominantly compliant low-aspect ratio pores in both materials, in addition Etna basalt also contains high-aspect ratio voids. In triaxial loading, stress-induced anisotropy of P-wave velocities was significantly higher for basalt than for granite, with vertical velocity components being faster than horizontal velocities. However, with increasing axial load, horizontal velocities show a small increase for basalt but a significant decrease for granite. Using first motion polarity we determined AE source types generated during triaxial loading of the samples. With increasing differential stress AE activity in granite and basalt increased with a significant contribution of tensile events. Close to failure the relative contribution of tensile events and horizontal wave velocities decreased significantly. A concomitant increase of double-couple events indicating shear, suggests shear cracks linking previously formed tensile cracks.  相似文献   

3.
Seismologically determined properties of the 400 km discontinuity may be compared to experimentally determined properties of the associated phase transformation in order to place constraints upon upper mantle bulk composition. Disagreement among previous studies is commonly ascribed to differences in elastic equations of state (especially to assumptions about pressure and temperature derivatives) between studies. However, much of the disparity between studies is actually due to the selection of different seismic data functionals (P-wave velocity,S-wave velocity, etc.) for comparison to minnral clasticity calculations, rather than to the differences in elasticity data sets and equations of state. Within any given study, bulk sound velocity comparisons generally yield more olivine-rich compositional estimates than doP-wave velocity comparisons, which in turn indicate more olivine thanS-wave velocities. Indeed, such variation in compositional estimates within a given study (arising from choice of data functional) exceeds the variation between studies (arising from elastic equation of state approx mations). it can be argued that bulk sound velocities are better constrained seismologically than densities and, being independent of assumptions about shear moduli, should provide more reliable compositional estimates thanP-orS-wave velocities.Using recently measured bulk and shear moduli equations of state, mutually consistent estimates of upper mantle olivine content can be obtained fromP-wave,S-wave, and bulk sound velocity contrasts at 400 km only if ln /T of has a value of about–2×10–4K–1, yielding approximately 52% olivine by volume. A value of ln /T smaller in magnitude would require reassessment of several underlying assumptions.  相似文献   

4.
We develop a data set of aftershock recordings of the 1999, M = 7.4 Izmit and M = 7.2 Duzce (Turkey) earthquakes to study their source parameters. We combined seismograms from 44 stations maintained by several sources (organizations) to obtain a unified data set of events (2.1 ≤ Mw ≤ 5.5). We calculate source parameters of these small earthquakes by two methods that use different techniques to address the difficulty in obtaining source spectra for small earthquakes subject to interference from site response. One method (program NetMoment (NM), Hutchings, 2004) uses spectra of direct S waves in a simultaneous inversion of local high-frequency network data to estimate seismic moment, source corner frequency (fc), site attenuation (k) and whole-path Q. This approach takes advantage of the source commonality in all recordings for a particular earthquake by fitting a common Brune source spectrum to the data with a and individual k. The second approach (Mayeda et al., 2003) uses the coda method (CM) to obtain “nonmodel-based” source spectra and moment estimates from selected broadband recording sites. We found that both methods do well for events that allow the comparison with seismic moment estimates derived from waveform modeling. Also, source spectra obtained from the two methods are very closely matched for most of the events they have in common. We use an F test to examine the trade-off between k and fc picks identified by the direct S-wave method. About half of the events could be constrained to have less than a 50% average uncertainty in fc and k. We used these source spectra solutions to calculate energy and apparent stress and compare these to estimates from the selected “good quality” source spectra from CM. Both studies have values mutually consistent and show a similar increase in apparent stress with increasing moment. This result has added merit due to the independent approaches to calculate apparent stress. We conclude that both methods are at least partially validated by our study, and they both have usefulness for different circumstances of recording local small earthquakes. CM would work well in studies for which there is a broad magnitude range of events and NM works well for local events recorded by band-limited recorders.  相似文献   

5.
The source parameters of the M W = 7.6 Olyutorskii earthquake were estimated using the moments of the slip rate function with degrees 1 and 2. The moments were estimated from broadband P-wave records at 52 stations of the worldwide network. The first step was to find a function S(t) for each station; this function is an apparent source time function, i.e., the P-wave slip as radiated by the source toward a station under consideration. The method of empirical Green’s functions was used to estimate S(t). The next step was to calculate the moments of S(t) of degrees 1 and 2 over time and to set up relevant equations to be solved by least squares for the unknown source moments. The horizontal linear source was used as a nonparametric model for calculating the source moments. Haskell’s parametric model was used for further interpretation of the source moments. The resulting estimates are as follows: the source centroid was 13–25 km southwest of the epicenter, the source was 105–120 km long, the source strike was 222°–228°, the rupture velocity was 2.7–3.0 km/s, and the total radiation duration was 24–27 s. These estimates indicate a bilateral rupture dominated by a southwestward sense of rupture propagation. The source characteristics are consistent with the aftershock area geometry and with the focal mechanism, as well as with surface breakage as observed by geologists in the field.  相似文献   

6.
We analyzed the local earthquakes waveform recorded on a broadband seismic network in the northwestern Himalayan Region to compute lapse time and frequency dependence of coda Q (Q c). The observed Q c values increase with increasing lapse time at all frequency bands. The increase in Q c values with lapse time is attributed to an increase in Q c with depth. This implies that attenuation decreases with increasing depth. The approximate radius of medium contributing to coda generation varies from 55 to 130 km. By comparing the Q c values with those from other regions of the world, we find that they are similar to those obtained from tectonically active regions. The estimated Q c values show a frequency-dependent relationship, Q c = Q 0 f n , where Q 0 is Q c at 1 Hz and n represents degree of frequency dependence. They represent the level of heterogeneity and tectonic activity in an area. Our results show that northwest Himalayas are highly heterogeneous and tectonically very active. Q 0 increases from 113 ± 7 to 243 ± 10 and n decreases from 1.01 ± 0.05 to 0.85 ± 0.03 when lapse time increases from 30 to 70 s. As larger time window sees the effect of deeper part of the Earth, it is concluded that Q 0 increases and n decreases with increasing depth; i.e., heterogeneity decreases with depth in the study area.  相似文献   

7.
The western part of Anatolia is one of the most seismically and tectonically active continental regions in the world, and much of it has been undergoing NS-directed extensional deformation since the Early Miocene. In this study, we determine 3-D tomographic images of the crust under the southwestern part of the North Anatolian Fault Zone by inverting a large number of arrival time data of P and S waves. From the obtained P- and S-wave velocity models, we estimated the Poisson’s ratio structures for a more reliable interpretation of the obtained anomalies. Our tomographic results confirmed the major tectonic features detected by previous studies and revealed new structural heterogeneities related to the active seismotectonics of the studied area. High P-wave velocity anomalies are recognized near the surface, while at deeper crustal layers, low P-wave velocities are widely distributed. The crustal S-wave velocity and Poisson’s ratio exhibit more structural heterogeneities compared to the P-wave velocity structure. Microearthquake activity is intense along highly heterogeneous zones in the southwestern part, which is characterized by low to high P-wave velocity, low S-wave velocity, and high Poisson’s ratio anomalies. Large earthquakes are also concentrated in zones dominated by low velocities and low to high Poisson’s ratios. Results of the checkerboard and synthetic tests indicate that the imaged anomalies are reliable features down to a depth of 25 km. Moreover, they are consistent with many geological and geophysical results obtained by other researchers along the southwestern part of the North Anatolian Fault Zone. An erratum to this article can be found at  相似文献   

8.
A layeredP- andS-wave velocity model is obtained for the Friuli seismic area using the arrival time data ofP- andS-waves from local earthquakes. A damped least-squares method is applied in the inversion.The data used are 994P-wave arrival times for 177 events which have epicenters in the region covered by the Friuli seismic network operated by Osservatorio Geofisico sperimentale (OGS) di Trieste, which are jointly inverted for the earthquake hypocenters andP-wave velocity model. TheS-wave velocity model is estimated on the basis of 978S-wave arrival times and the hypocenters obtained from theP-wave arrival time inversion. We also applied an approach thatP- andS-wave arrival time data are jointly used in the inversion (Roecker, 1982). The results show thatS-wave velocity structures obtained from the two methods are quite consistent, butP-wave velocity structures have obvious differences. This is apparent becauseP-waves are more sensitive to the hypocentral location thanS-waves, and the reading errors ofS-wave arrival times, which are much larger than those ofP-waves, bring large location errors in the joint inversion ofP- andS-wave arrival time. The synthetic data tests indicated that when the reading errors ofS-wave arrivals are larger than four times that ofP-wave arrivals, the method proposed in this paper seems more valid thanP- andS-wave data joint inversion. Most of the relocated events occurred in the depth range between 7 and 11 km, just above the biggest jump in velocity. This jump might be related to the detachment line hypothesized byCarulli et al. (1982). From the invertedP- andS-wave velocities, we obtain an average value 1.82 forV p /V s in the first 16 km depth.  相似文献   

9.
Numerical algorithms are suggested for reconstructing the potential field in the region Ω = {0 ≤ R 0 ≤ ρ ≤ R, 0 ≤ θ ≤ π, 0 ≤ φ ≤ 2π} from its values on the surface Γ1 = {R,0 ≤ α1 ≤ θ ≤ α2 ≤ π, 0 ≤ β1 ≤ φ ≤ β2 ≤ π}.  相似文献   

10.
We analyze the strong motion accelerograms of the moderate (M w = 6.1), March 31, 2006, Darb-e-Astane earthquake of western Iran and also those of one of its prominently recorded, large (M w = 5.1) foreshock and (M w = 4.9) aftershock. (1) Using derived SH-wave spectral data, we first objectively estimate the parameters W o\mathit{\Omega} _{\rm o} (long period spectral level), f c (corner frequency) and Q(f) (frequency dependent, average shear wave quality factor), appropriate for the best-fit Brune ω  − 2 spectrum of each of these three events. We then perform a non-linear least square analysis of the SH-wave spectral data to provide approximate near-field estimates of the strike, dip, and rake of the causative faults and also the seismic moment, moment magnitude, source size, and average stress drop of these three events. (2) In the next step, we use these approximate values and an empirical Green’s function approach, in an iterative manner, to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration/velocity/displacement and duration of ground shaking and thereby provide improved, more reliable estimates of the causative fault parameters of the main event and its asperities. Our near-field estimates for both the main moderate event and the two smaller events are in good conformity with the corresponding far-field estimates reported by other studies.  相似文献   

11.
We analysed the complications in laboratory velocity anisotropy measurement on shales. There exist significant uncertainties in the laboratory determination of c13 and Thomsen parameter δ. These uncertainties are primarily related to the velocity measurement in the oblique direction. For reliable estimation of c13 and δ, it is important that genuine phase velocity or group velocity be measured with minimum uncertainty. The uncertainties can be greatly reduced if redundant oblique velocities are measured. For industrial applications, it is impractical to make multiple oblique velocity measurements on multiple core plugs. We demonstrated that it is applicable to make multiple genuine oblique group velocity measurements on a single horizontal core plug. The measurement results show that shales can be classified as a typical transversely isotropic medium. There is a coupling relation between c44 and c13 in determining the directional dependence of the seismic velocities. The quasi‐P‐wave or quasi‐S‐wave velocities can be approximated by three elastic parameters.  相似文献   

12.
P- and S-wave travel times from local volcanic earthquakes recorded in the North group of volcanoes area during the 2005–2009 period were treated by the “reverse wave” method to calculate the V P velocity field and the TAU parameter, which is an analogue of the P- to S-wave velocity ratio. We constructed 3D velocity distributions along the line traversing the volcanic group along the direction from Ploskii Tolbachik Volcano in the southwest toward Shiveluch Volcano in the northeast. Dynamic changes in the velocity field were identified, both over time and depth. We examine the relationships of these dynamic changes to the evolution of volcanic activity during the period indicated.  相似文献   

13.
The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulav CMB=r c /p (v CMB=velocity at the CMB;r c =core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for Pdiff and Sdiff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe.  相似文献   

14.
Fracturing and hydrothermal alteration in normal fault zones   总被引：9，自引：0，他引：9  
Large normal fault zones are characterized by intense fracturing and hydrothermal alteration. Displacement is localized in a slip zone of cataclasite, breccia and phyllonite surrounding corrugated and striated fault surfaces. Slip zone rock grades into fractured, but less comminuted and hydrothermally altered rock in the transition zone, which in turn grades abruptly into the wall rock. Fracturing and fluid flow is episodic, because permeability generated during earthquakes is destroyed by hydrothermal processes during the time between earthquakes.Fracture networks are described by a fracture fabric tensor (F). The permeability tensor (k) is used to estimate fluid transport properties if the trace of F is sufficiently large. Variations in elastic moduli and seismic velocities between fault zone and wall rock are estimated as a function of fracture density (). Fracturing decreases elastic moduli in the transition zone by 50–100% relative to the country rock, and similar or even greater changes presumably occur in the slip zone.P-andS-wave velocity decrease, andV p /V s increases in the fault zone relative to the wall rock. Fracture permeability is highly variable, ranging between 10–13 m2 and 10–19 m2 at depths near 10 km. Changes in permeability arise from variations in effective stress and fracture sealing and healing.Hydrothermal alteration of quartzo-feldspathic rock atT>300°C creates mica, chlorite, epidote and alters the quartz content. Alteration changes elastic moduli, but the changes are much less than those caused by fracturing.P-andS-wave velocities also decrease in the hydrothermally altered fault rock relative to the country rock, and there is a slight decrease inV p /V s , which partially offsets the increase inV p /V s caused by fracturing.Fracturing and hydrothermal alteration affect fault mechanics. Low modulus rock surrounding fault surfaces increases the probability of exceeding the critical slip distance required for the onset of unstable slip during rupture initiation. Boundaries between low modulus fault rock and higher modulus wall rock also act as rupture guides and enhance rupture acceleration to dynamic velocity. Hydrothermal alteration at temperatures in excess of 300°C weakens the deeper parts of the fault zone by producingphyllitic mineral assemblages. Sealing of fracture in time periods between large earthquakes generates pods of abnormally pressured fluid which may play a fundamental role in the initiation of large earthquakes.  相似文献   

15.
A method for point estimation of coseismic volume deformation from amplitudes of earthquake-induced jumps in the water level and from the tidal sensitivity of the water level is presented. Application of this method is illustrated by six Kamchatka earthquakes of 1997–2004 with M w = 6.2–7.8 that occurred 128–316 km from the YuZ-5 well. Estimates of coseismic deformation provided by the level gage observations and by the model agree in amplitude and sign. The model reproduces an extended dislocation source in a homogeneous elastic isotropic half-space. Source rupture parameters correspond to focal mechanisms from the CMT international catalog. Consistence in experimental and theoretical estimates of the volume coseismic deformation indicates the high strain-resistance properties of YuZ-5 borehole and the possibility of quantitative estimation of seismotectonic deformation from high-precision water level measurements with time resolution not lower than 10 minutes.  相似文献   

16.
Pore-pressure depletion causes changes in the triaxial stress state. Pore-pressure depletion in a flat reservoir, for example, can be reasonably approximated as uniaxial compaction, in which the horizontal effective stress change is smaller than the vertical effective stress. Furthermore, the stress sensitivity of velocities can be angle-dependent. Therefore, time-lapse changes in reservoir elastic anisotropy are expected as a consequence of production, which can complicate the interpretation of the 4D seismic response. The anisotropic 4D seismic response caused by pore-pressure depletion was investigated using existing core velocity measurements. To make a direct comparison between the anisotropic 4D seismic response and the isotropic response based only on vertical velocities, pseudoisotropic elastic properties were utilized, and the two responses were compared in terms of a dynamic rock physics template. A comparison of the dynamic rock physics templates indicates that time-lapse changes in reservoir elastic anisotropy have a noticeable impact on the interpretation of 4D seismic data. Changes in anisotropy as a result of pore-pressure depletion cause a time-lapse amplitude variation with offset response as if there is a reduction in VP/VS (i.e., pseudoisotropic VP/VS decreases), although the vertical VP/VS increases. The impact of time-lapse changes in anisotropy on the amplitude variation with offset gradient was also investigated, and the time-lapse anisotropy was found to enhance changes in the amplitude variation with offset gradient for a given case.  相似文献   

17.
High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (VS), compressional (P)-wave velocity (VP), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor QP and S-wave quality factor QS in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting QS from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from QP cannot be ignored when Vs/VP reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert QP from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to QP than to QS. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.  相似文献   

18.
Erosion of cohesive soils in fluvial environments is dependent on physical, geochemical and biological properties, which govern inter‐particle attraction forces and control detachment rates from stream beds and banks. Most erosion rate models are based on the excess shear stress equation where the soil erodibility coefficient (kd) is multiplied by the difference between the boundary hydraulic shear stress (τb) and the soil critical shear stress (τc). Both kd and τc are a function of soil properties and must be obtained through in situ field or laboratory testing. Many studies have generated predictive relationships for kd and τc derived from various soil properties. These studies typically were conducted in watersheds within a single physiographic region with a common surficial geology and/or investigated a limited number of soil properties, particularly geochemical properties. With widely reported differences in relationships between τc and soil properties, this study investigated differences in predictive relationships for τc among different physiographic provinces in Tennessee, USA. Erodibility parameters were determined in the field using a mini‐jet test device. Among these provinces, statistically four unique clusters were identified from a dataset of 128 observations and these data clusters were used to develop predictive models for τc to identify dominant properties governing erosion. In these clusters, 16 significant physical and geochemical soil properties were identified for τc prediction. Among these soil properties, water content and passing #200 sieve (percentage soil less than 75 μm) were the dominant controlling parameters to predict τc in addition to clay percentage (< 2 μm), bulk density, and soil pore water chemistry. This study suggests that unique relationships exist for physiographic provinces that are likely due to soil physical‐geochemical processes associated with surficial geology that determine minerology of the cohesive soil. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

19.
The paper presents a model of the kinetics of electronically excited O2(c1Σu−,v), O2(A′3Δu,v), O2(A3Σu+,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O2(A3Σu+,v) and O2(A′3Δu,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.  相似文献   

20.
Velocity as well as attenuation factorQ –1 ofP-wave in a dry granitic rock sample under uniaxial compressions were measured in the range of frequency between 100 kHz and 710 kHz by using the pulse transmission technique. Above the stress of 0.5 f , where f is the fracture stress, theP-wave velocity decreases with increasing axial stress, whereasQ –1 increases. Particularly, the change ofQ –1 is greater for high frequency than for low frequency. At a given stress level, the higher the frequency, the higher theP-wave velocity and the largerQ –1. This result means that the velocity decrease with increasing stress is smaller for higher frequency. Because of this frequency-dependence of velocity decrease, theP-wave in the rock under dilatant state shows dispersion. The body wave dispersion is more remarkable at higher stress, and is not found in a homogeneous material with no cracks. Thus the disperison is attributed to the generation of cracks. When the frequency-dependence ofQ –1 is approximated asf n in the present frequency range, the exponentn takes a value from 0.63 to 0.77.  相似文献   

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

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