共查询到20条相似文献，搜索用时 31 毫秒
1.
We refine the 1-D velocity model of the Central India Tectonic Zone (CITZ) using well-selected arrival times of P- and S-phases of 354 local earthquakes of magnitude (Mw) between 2.0 and 5.8, recorded by national seismic network from May 1997 to March 2016. Further, we have determined the source mechanisms of 26 selected local events using moment tensor inversion to characterize the dynamics beneath the CITZ. The best-fit simulation between observed and synthetic waveforms obtained the nodal and auxiliary planes of the each faults associated with the earthquake moment magnitude (Mw) for each event. Depth of the fault plane along the CITZ varies from 5 to 38 km. From this study, we found that the western part along the CITZ shows minimum focal depth and reaches maximum 38 kms at Jabalpur in the eastern part. This complex nature of earthquake dynamics occurrence along the CITZ. We propose that the curviplanar the CITZ dominated with sinistral curvature is subjected to compression along the longer ~E–W segments and transtension along shorter segments with ~NE–SW orientations. The occurrences of normal faulting, intrusion of mafic plutons and CLVD mechanisms for earthquakes are interpreted to be linked to the transtension zones and reverse mechanisms associated with the compressions along ~E–W segments. 相似文献
2.
In this paper, we present a seismic hazard scenario for the Garhwal region of the north-western Himalayan range, in terms of the horizontal Peak Ground Acceleration. The scenario earthquake of moment magnitude M w 8.5 has a 10% exceedance probability over the next 50 years. These estimates, the first for the region, were calculated through a stepwise process based on:
相似文献
- An estimation of the Maximum Credible Earthquake from the seismicity of the region and Global Seismic Hazard Assessment Program considerations, and
- four seismotectonic parameters abstracted from near field weak-motion data recorded at five stations installed in Chamoli District of the Garhwal region in the aftermath of the 1999 Chamoli earthquake. The latter include
- The frequency dependent power law for the shear wave quality factor, Q S
- the site amplification at each station using horizontal-to-vertical-spectral ratio and generalized inversion technique
- source parameters of various events recorded by the array and application of the resulting relations between the scalar seismic moment M 0 (dyne-cm) and moment magnitude M w and the corner frequency, ? c (Hz) and moment magnitude M w to simulate spectral acceleration due to higher magnitude events corresponding to the estimated Maximum Credible Earthquake, and
- regional and site specific local spectral attenuation relations at different geometrically central frequencies in the low, moderate and high frequency bands.
3.
Earthquake source parameters and crustal \(Q_{0}\) values for the 138 selected local events of (\(\hbox {M}_{\mathrm{w}}{:}2.5{-}4.4\)) the 2001 Bhuj earthquake sequence have been computed through inversion modelling of S-waves from three-component broadband seismometer data. SEISAN software has been used to locate the identified local earthquakes, which were recorded at least three or more stations of the Kachchh seismological network. Three component spectra of S-wave are being inverted by using the Levenberg–Marquardt non-linear inversion technique, wherein the inversion scheme is formulated based on \(\omega ^{2}\) source model. SAC Software (seismic analysis code) is being utilized for calculating three-component displacement and velocity spectra of S-wave. The displacement spectra are used for estimating corner frequency (in Hz) and long period spectral level (in nm-s). These two parameters play a key role in estimating earthquake source parameters. The crustal \({Q}_{0}\) values have been computed simultaneously for each component of three-component broadband seismograph. The estimated seismic moment (\(M_{0}\)) and source radius (r) using S-wave spectra range from 7.03E+12 to 5.36E+15 N-m and 178.56 to 565.21 m, respectively. The corner frequencies for S-wave vary from 3.025 to 7.425 Hz. We also estimated the radiated energy (\(E_{S}\)) using velocity spectra, which is varying from 2.76E+06 to 4.07E+11 Joules. The estimated apparent stress drop and static stress drop values range from 0.01 to 2.56 and 0.53 to 36.79 MPa, respectively. Our study also reveals that estimated \(Q_{0}\) values vary from 119.0 to 7229.5, with an average \(Q_{0}\) value of 701. Another important parameter, by which the earthquake rupture process can be recognized, is Zuniga parameter. It suggests that most of the Kachchh events follow the frictional overshoot model. Our estimated static stress drop values are higher than the apparent stress drop values. And the stress drop values are quite larger for intraplate earthquakes than the interplate earthquakes. 相似文献
4.
Zoya Farajpour Mehdi Zare Shahram Pezeshk Anooshiravan Ansari Esmael Farzanegan 《Arabian Journal of Geosciences》2018,11(4):80
This paper discusses a newly developed high-quality integrated dataset of shallow earthquake ground motions that occurred in Iran, from 1976 to 2013. A total of 860 three-component strong motion records are processed from 183 earthquake events, moment magnitudes 5.0?≤?M w ?≤?7.4, and rupture distances of R RUP ≤ 120 km. Strong motion data from Iran having special tectonic features and shallow earthquakes with depths less than 35 km are included. This paper presents a thorough procedure used to collect and to generate a database following the Next-Generation Attenuation-West research projects. This database can be used in the development and ranking of ground motion models and for seismological and engineering hazard and risk analyses. Unprocessed strong motion records are obtained from the Iranian Strong Motion Network (ISMN). The time series collected were thoroughly examined through several rounds of quality reviews. The newly generated database includes the peak ground acceleration, peak ground velocity, and pseudo-spectral acceleration for the 5% damped with periods ranging from 0.01 to 10 s. The database also includes ground motion information and source characterization and parameters. This study is the near-source compiled ground motion database that can be used for Iran, and it is consistent with standard worldwide databases. 相似文献
5.
Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D c) of the hypocenters, are estimated. The b-value and the D c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D c-values in both the tectonic regions. 相似文献
6.
We analyze previously published geodetic data and intensity values for the M s = 8.1 Shillong (1897), M s = 7.8 Kangra (1905), and M s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data. 相似文献
7.
The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤ b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days. 相似文献
8.
The high-pressure behavior of a vanadinite (Pb10(VO4)6Cl2, a = b = 10.3254(5), c = 7.3450(4) Å, space group P63/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V 0 = 681(1) Å3, K 0 = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a 0 = 10.3302(2) Å, K 0(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c 0 = 7.3520(3) Å, K 0(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K 0(a):K 0(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe a (0) = 35(2) GPa for the a-axis, Fe c (0) = 98(4) GPa for the c-axis and Fe V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed. 相似文献
9.
The properties of the source spectra of local shallow-focus earthquakes on Kamchatka in the range of magnitudes M w = 3.5–6.5 are studied using 460 records of S-waves obtained at the PET station. The family of average source spectra is constructed; the spectra are used to study the relationship between M w and the key quasi-dimensionless source parameters: stress drop Δσ and apparent stress σa. It is found that the parameter Δσ is almost stable, while σa grows steadily as the magnitude M w increases, indicating that the similarity is violated. It is known that at sufficiently large M w the similarity hypothesis is approximately valid: both parameters Δσ and σa do not show any noticeable magnitude dependence. It has been established that M w ≈ 5.7 is the threshold value of the magnitude when the change in regimes described occurs for the conditions on Kamchatka. 相似文献
10.
M. Yu. Piotrovich N. A. Silant’ev Yu. N. Gnedin T. M. Natsvlishvili S. D. Buliga 《Astronomy Reports》2016,60(5):486-497
The magnetic-field structure in regions of stationary, planar accretion disks around active galactic nuclei where general-relativistic effects can be neglected (from 10 to 200 gravitational radii) is considered. It is assumed that the magnetic field in the outer edges of the disk, which forms in the magnetosphere of the central black hole during the creation of the relativisitic jets, corresponds to the field of a magnetic dipole perpendicular to the plane of the disk. In this case, the azimuthal field component Bφ in the disk arises due to the presence of the radial field Bρ and the azimuthal velocity component Uφ. The value of the magnetic field at the inner radius of the disk is taken to correspond to the solution of the induction equation in a diffusion approximation. Numerical solutions of the induction equation are given for a number of cases. 相似文献
11.
On Thursday, 22 of May 2014, at 6 h 22 min 0.3.3 s (GMT?+?1) a moderate-sized earthquake struck the Mostaganem, Western Algeria, region. The main shock, recorded by many international and national seismological stations, was preceded by a foreshock, 3 hours before, on May 22, 2014 (Ml?=?4.1) at 3 h 57 min 41.4 s and followed by four well-felt aftershocks (M?>?3.0) that lasted about 1 year. The main shock did not cause loss of lives but serious panic among the population was reported. The main shock, however, caused cracks in walls and roofs, sometimes destroyed, the old non-engineered and precarious adobe dweller corresponding to I0?=?VI–VII (Msk scale). We used accelerograph records to (i) determine the epicenter location (longitude?=?0.3537 E, latitude?=?35.8598 N, (ii) perform waveforms inversion to calculate the earthquake parameters. The obtained results are, respectively, the seismic moment (M0)?=?2.71 E + 16, the Mw?=?4.9 and the focal depth?=?6 km. The obtained focal mechanism solution shows reverse faulting with small right lateral component with the following nodal plans: NP1, strike?=?193.5, dip?=?49.5, slip?=?57.6 and NP2, strike?=?57.8, dip?=?50, slip?=?122.1. On the other hand, the seismotectonic framework of the Dahra area exhibits a serie of NE-SW trending “en echelon” faulted folds that may be active as suggested by this study. 相似文献
12.
Omphacite is an important mineral component of eclogite. Single-crystal synchrotron X-ray diffraction data on natural (Ca, Na) (Mg, Fe, Al)Si2O6 omphacite have been collected at the Advanced Photon Source beamlines 13-BM-C and 13-ID-D up to 47 GPa at ambient temperature. Unit cell parameter and crystal structure refinements were carried out to constrain the isothermal equation of state and compression mechanism. The third-order Birch–Murnaghan equation of state (BM3) fit of all data gives V 0 = 423.9(3) Å3, K T0 = 116(2) GPa and K T0′ = 4.3(2). These elastic parameters are consistent with the general trend of the diopside–jadeite join. The eight-coordinated polyhedra (M2 and M21) are the most compressible and contribute to majority of the unit cell compression, while the SiO4 tetrahedra (Si1 and Si2) behave as rigid structural units and are the most incompressible. Axial compressibilities are determined by fitting linearized BM3 equation of state to pressure dependences of unit cell parameters. Throughout the investigated pressure range, the b-axis is more compressible than the c-axis. The axial compressibility of the a-axis is the largest among the three axes at 0 GPa, yet it quickly drops to the smallest at pressures above 5 GPa, which is explained by the rotation of the stiffest major compression axis toward the a-axis with the increase in pressure. 相似文献
13.
Shuangmeng Zhai Weihong Xue Daisuke Yamazaki Shuangming Shan Eiji Ito Naotaka Tomioka Akira Shimojuku Ken-ichi Funakoshi 《Physics and Chemistry of Minerals》2011,38(5):357-361
High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr3(PO4)2 has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V 0 = 498.0 ± 0.1 Å3, an isothermal bulk modulus of K T = 89.5 ± 1.7 GPa, and first pressure derivative of K T ′ = 6.57 ± 0.34. If K T ′ is fixed at 4, K T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K a = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K c = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO4 vibrations of Sr3(PO4)2 is calculated to be 0.30(2). 相似文献
14.
In this study, the seismicity rate changes that can represent an earthquake precursor were investigated along the Sagaing Fault Zone (SFZ), Central Myanmar, using the Z value technique. After statistical improvement of the existing seismicity data (the instrumental earthquake records) by removal of the foreshocks and aftershocks and man-made seismicity changes and standardization of the reported magnitude scales, 3574 earthquake events with a M w ≥ 4.2 reported during 1977–2015 were found to directly represent the seismotectonic activities of the SFZ. To find the characteristic parameters specifically suitable for the SFZ, seven known events of M w ≥ 6.0 earthquakes were recognized and used for retrospective tests. As a result, utilizing the conditions of 25 fixed earthquake events considered (N) and a 2-year time window (T w), a significantly high Z value was found to precede most of the M w ≥ 6.0 earthquakes. Therefore, to evaluate the prospective areas of upcoming earthquakes, these conditions (N = 25 and T w = 2) were applied with the most up-to-date seismicity data of 2010–2015. The results illustrate that the vicinity of Myitkyina and Naypyidaw (Z = 4.2–5.1) cities might be subject to strong or major earthquakes in the future. 相似文献
15.
N. H. Tuyen Ph. V. Phach R. B. Shakirov C. D. Trong Ph. N. Hung L. D. Anh 《Geotectonics》2018,52(3):359-381
The goal of the paper is to reveal the hierarchical block structure in the Tuan Giao area (Northwest Vietnam) and, on that base, to identify areas prone to earthquakes with M ≥ 5. Four large geoblocks of the second rank have been delineated on the basis of the joint analysis of geological, geophysical, geomorphic, and remote sensing data. The second rank geoblocks have been divided into smaller sub-blocks of the third rank. The recent geodynamics of the geoblocks have been characterized using geomorphic, seismological, gravity, and GPS data.The system of the delineated geoblocks is viewed as the Geodynamic Blocks model. The areas prone to earthquakes with M ≥ 5 have been identified using the pattern recognition algorithm CORA-3. The objects of the recognition were defined as circular areas, where boundaries of the geoblocks intersect each other. The recognition results confirm high seismic potential of the study region and provide information on potential earthquake sources for seismic hazard assessment: a number of boundary intersections have been identified as areas prone to earthquakes with M ≥ 5, where events of such magnitude have not been recorded up to date. Mordern geothermal activity and anomalies of dissolved methane (up to 10000 nL/L), helium and hydrogen in the hot mineral water springs prove high permeability of the study area and its active seismo-tectonic state. 相似文献
16.
Sankar Kumar Nath Kiran Kumar Singh Thingbaijam Abhishek Raj 《Journal of Earth System Science》2008,117(2):809-831
A comprehensive analytical as well as numerical treatment of seismological, geological, geomorphological and geotechnical concepts has been implemented through microzonation projects in the northeast Indian provinces of Sikkim Himalaya and Guwahati city, representing cases of contrasting geological backgrounds — a hilly terrain and a predominantly alluvial basin respectively. The estimated maximum earthquakes in the underlying seismic source zones, demarcated in the broad northeast Indian region, implicates scenario earthquakes of M W 8.3 and 8.7 to the respective study regions for deterministic seismic hazard assessments. The microzonation approach as undertaken in the present analyses involves multi-criteria seismic hazard evaluation through thematic integration of contributing factors. The geomorphological themes for Sikkim Himalaya include surface geology, soil cover, slope, rock outcrop and landslide integrated to achieve geological hazard distribution. Seismological themes, namely surface consistent peak ground acceleration and predominant frequency were, thereafter, overlaid on and added with the geological hazard distribution to obtain the seismic hazard microzonation map of the Sikkim Himalaya. On the other hand, the microzonation study of Guwahati city accounts for eight themes — geological and geomorphological, basement or bedrock, landuse, landslide, factor of safety for soil stability, shear wave velocity, predominant frequency, and surface consistent peak ground acceleration. The five broad qualitative hazard classifications — ‘low’, ‘moderate’, ‘high’, ‘moderate high’ and ‘very high’ could be applied in both the cases, albeit with different implications to peak ground acceleration variations. These developed hazard maps offer better representation of the local specific seismic hazard variation in the terrain. 相似文献
17.
G. Diego Gatta Nicola Rotiroti Martin Fisch Milen Kadiyski Thomas Armbruster 《Physics and Chemistry of Minerals》2008,35(9):521-533
The elastic and structural behaviour of the synthetic zeolite CsAlSi5O12 (a = 16.753(4), b = 13.797(3) and c = 5.0235(17) Å, space group Ama2, Z = 2) were investigated up to 8.5 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase-transition occurs within the P-range investigated. Fitting the volume data with a third-order Birch–Murnaghan equation-of-state gives: V 0 = 1,155(4) Å3, K T0 = 20(1) GPa and K′ = 6.5(7). The “axial moduli” were calculated with a third-order “linearized” BM-EoS, substituting the cube of the individual lattice parameter (a 3, b 3, c 3) for the volume. The refined axial-EoS parameters are: a 0 = 16.701(44) Å, K T0a = 14(2) GPa (βa = 0.024(3) GPa?1), K′ a = 6.2(8) for the a-axis; b 0 = 13.778(20) Å, K T0b = 21(3) GPa (βb = 0.016(2) GPa?1), K′ b = 10(2) for the b-axis; c 0 = 5.018(7) Å, K T0c = 33(3) GPa (βc = 0.010(1) GPa?1), K′ c = 3.2(8) for the c-axis (K T0a:K T0b:K T0c = 1:1.50:2.36). The HP-crystal structure evolution was studied on the basis of several structural refinements at different pressures: 0.0001 GPa (with crystal in DAC without any pressure medium), 1.58(3), 1.75(4), 1.94(6), 3.25(4), 4.69(5), 7.36(6), 8.45(5) and 0.0001 GPa (after decompression). The main deformation mechanisms at high-pressure are basically driven by tetrahedral tilting, the tetrahedra behaving as rigid-units. A change in the compressional mechanisms was observed at P ≤ 2 GPa. The P-induced structural rearrangement up to 8.5 GPa is completely reversible. The high thermo-elastic stability of CsAlSi5O12, the immobility of Cs at HT/HP-conditions, the preservation of crystallinity at least up to 8.5 GPa and 1,000°C in elastic regime and the extremely low leaching rate of Cs from CsAlSi5O12 allow to consider this open-framework silicate as functional material potentially usable for fixation and deposition of Cs radioisotopes. 相似文献
18.
Both abiotic and biotic factors govern distributions of estuarine vegetation, and experiments can reveal effects of these drivers under current and future conditions. In upper San Francisco Estuary (SFE), increased salinity could result from sea level rise, levee failure, or water management. We used mesocosms to test salinity effects on, as well as competition between, the native Stuckenia pectinata (sago pondweed) and invasive Egeria densa (Brazilian waterweed), species with overlapping distributions at the freshwater transition in SFE. Grown alone at a salinity of 5, E. densa decreased fivefold in biomass relative to the freshwater treatment and decomposed within 3 weeks at higher salinities. In contrast, S. pectinata biomass accumulated greatly (~4× initial) at salinities of 0 and 5, doubled at 10, and was unchanged at 15. When grown together in freshwater, S. pectinata produced 75 % less biomass than in monoculture and significantly more nodal roots (suggesting increased nutrient foraging). At a salinity of 5, a decline in E. densa performance coincided with a doubling of S. pectinata shoot density. Additional experiments on E. densa showed elevated temperature (26 and 30 °C) suppressed growth especially at higher salinities (≥5). We conclude that salinity strongly influences distributions of both species and that competition from E. densa may impose limits on S. pectinata abundance in the fresher reaches of SFE. With a salinity increase of 5, S. pectinata is likely to maintain its current distribution while spreading up-estuary at the expense of E. densa, especially if increased temperature also reduces E. densa biomass. 相似文献
19.
This paper presents the first boundary equations describing the relationship between earthquake parameters (magnitude M S and macroseismic intensity I P at the observation point on the MSK-64 scale) and clastic dikes (having maximal thickness m cd , visible height h cd , and the index of manifestation intensity of dikes in the cross section I cd ). As was expected, the maximal size of dikes grows with an increase in the earthquake magnitude and macroseismic intensity. Analysis of the dependences showed that it is better to use all three parameters for estimation of the minimal threshold M S or I P from clastic dikes, and, in the absence of data on seismogenic rupture, the maximal calculated value should be used. Some limitations in application and the advantages of the equations obtained are discussed with respect to characterizing earthquakes of the pre-instrumental period. 相似文献
20.
Channel sand acts as a stratigraphic trap for hydrocarbon accumulation in many parts of the world. Delineation of this type of reservoir is crucial as channel sand may be scarce, and inaccurate location of the drilling wells could lose a huge currency. The Hassi Messaoud (HMD) field was subjected to multiphase tectonic events, where deep-seated structures were rejuvenated leading to intensive fault complexity. The main effective tectonic events upon the studied area are the Hercynian compression and deep erosion till the Ordovician Hamra Quartzite (HQZ) oil reservoir, followed by active Triassic rifting and filling the deeply eroded areas or the graben areas by eruptive volcanic rocks at Triassic time. Hercynian erosion and volcanic rocks distribution introduce a big uncertainty to the reservoir structural model. Amplitude versus offset (AVO) method is used as a helpful tool to differentiate channel sand from surrounding formations. Several attributes (P-impedance, S-impedance, longitudinal velocity Vp, shear velocity Vs and density ρ) are estimated from pre-stack seismic inversion. They have different sensitivity to the reservoir properties. Derived attributes such as Lamé parameters, incompressibility × density (λρ) and rigidity × density (μρ) can provide key lithology and fluid indicators (Goodway et al. 1997, Goodway CSEG Rec 26(6):39-60 2001). Petrophysically relating AVO attributes both to λρ and μρ and to each other in Lambda–Mu–Rho (LMR) cross-plot space can be a good tool for AVO interpretation (Rutherford and Williams Geophysics 54:680–688 1989 and Castagna and Swan Lead Edge 16(4):337–342 1997). After proper data conditioning, simultaneous inversion of pre-stack angle gathers is performed to get acoustic wave impedance (P-impedance), elastic wave impedances (S-impedance) and density ρ, then to calculate λρ and μρ volumes. In the studied area, λρ and μρ are used as a very important key to separate reservoir sands. The λρ and μρ curves are generated at each well location. Cross plots showed a fair separation of sand in the formation, i.e. higher μρ and lower λρ can detect sand. The output λρ and μρ volumes after simultaneous inversion follow the distribution of the sand which is consistent with the wells penetrating the target reservoir. This finding on the extension of the sand reservoir in terms of λρ and μρ. 3D cross-plot zonations are used for lithology discrimination. In this study, well logs were used to constrain lithology and to control the zonation filters by reducing the limits ambiguity. Other types of advanced attributes are calculated and tested. The obtained (μρ–λρ) volume acts as a good indicator for the sand distribution. It was finally used as sand presence index in the area. Also μρ has shown a good linear relationship with porosity. To note that the porosity volume is created based on the linear relationship with μρ. A product of derived porosity and the sand presence index (μρ–λρ) provides a good tool for reservoir characterization and lead to reservoir management, future planning of the field, and setting location for new wells. 相似文献