AVO investigations of shallow marine sediments

Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m³, quality factor for P‐wave Q_P=15±200, quality factor for S‐wave Q_S=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m³,Q_P=790±660,Q_S=25±10.     

Characterisation and multifaceted anisotropy assessment of Corvio sandstone for geological CO2 storage studies

We present a comprehensive characterisation of the physical, mineralogical, geomechanical, geophysical, and hydrodynamic properties of Corvio sandstone. This information, together with a detailed assessment of anisotropy, is needed to establish Corvio sandstone as a useful laboratory rock‐testing standard for well‐constrained studies of thermo–hydro–mechanical–chemical coupled phenomena associated with CO₂ storage practices and for geological reservoir studies in general. More than 200 core plugs of Corvio sandstone (38.1 and 50 mm diameters, 2:1 length‐to‐diameter ratio) were used in this characterisation study, with a rock porosity of 21.7 ± 1.2%, dry density 2036 ± 32 kg m^?3, and unconfined compressive and tensile strengths of 41 ± 3.28 and 2.3 ± 0.14 MPa, respectively. Geomechanical tests show that the rock behaves elastically between ～10 and ～18 MPa under unconfined conditions with associated Young's modulus and Poisson's ratio of 11.8 ± 2.8 GPa and 0.34 ± 0.01 GPa, respectively. Permeability abruptly decreases with confining pressure up to ～10 MPa and then stabilises at ～1 mD. Ultrasonic P‐ and S‐wave velocities vary from about 2.8–3.8 km s^?1 and 1.5–2.4 km s^?1, respectively, over confining and differential pressures between 0.1 and 35 MPa, allowing derivation of associated dynamic elastic moduli. Anisotropy was investigated using oriented core plugs for electrical resistivity, elastic wave velocity and attenuation, permeability, and tracer injection tests. Corvio sandstone shows weak transverse isotropy (symmetry axis normal to bedding) of <10% for velocity and <20% for attenuation.     

A new laboratory technique for determining the compressional wave properties of marine sediments at sonic frequencies and in situ pressures

We describe a new laboratory technique for measuring the compressional wave velocity and attenuation of jacketed samples of unconsolidated marine sediments within the acoustic (sonic) frequency range 1–10 kHz and at elevated differential (confining – pore) pressures up to 2.413 MPa (350 psi). The method is particularly well suited to attenuation studies because the large sample length (up to 0.6 m long, diameter 0.069 m) is equivalent to about one wavelength, thus giving representative bulk values for heterogeneous samples. Placing a sediment sample in a water‐filled, thick‐walled, stainless steel Pulse Tube causes the spectrum of a broadband acoustic pulse to be modified into a decaying series of maxima and minima, from which the Stoneley and compressional wave, velocity and attenuation of the sample can be determined. Experiments show that PVC and copper jackets have a negligible effect on the measured values of sediment velocity and attenuation, which are accurate to better than ± 1.5% for velocity and up to ± 5% for attenuation. Pulse Tube velocity and attenuation values for sand and silty‐clay samples agree well with published data for similar sediments, adjusted for pressure, temperature, salinity and frequency using standard equations. Attenuation in sand decreases with pressure to small values below Q^?1 = 0.01 (Q greater than 100) for differential pressures over 1.5 MPa, equivalent to sub‐seafloor depths of about 150 m. By contrast, attenuation in silty clay shows little pressure dependence and intermediate Q^?1 values between 0.0206–0.0235 (Q = 49–43). The attenuation results fill a notable gap in the grain size range of published data sets. Overall, we show that the Pulse Tube method gives reliable acoustic velocity and attenuation results for typical marine sediments.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.     

Pore fluid viscosity effects on P‐ and S‐wave anisotropy in synthetic silica‐cemented sandstone with aligned fractures

Ultrasonic (500 kHz) P‐ and S‐wave velocity and attenuation anisotropy were measured in the laboratory on synthetic, octagonal‐shaped, silica‐cemented sandstone samples with aligned penny‐shaped voids as a function of pore fluid viscosity. One control (blank) sample was manufactured without fractures, another sample with a known fracture density (measured from X‐ray CT images). Velocity and attenuation were measured in four directions relative to the bedding fabric (introduced during packing of successive layers of sand grains during sample construction) and the coincident penny‐shaped voids (fractures). Both samples were measured when saturated with air, water (viscosity 1 cP) and glycerin (100 cP) to reveal poro‐visco‐elastic effects on velocity and attenuation, and their anisotropy. The blank sample was used to estimate the background anisotropy of the host rock in the fractured sample; the bedding fabric was found to show transverse isotropy with shear wave splitting (SWS) of 1.45 ± 1.18% (i.e. for S‐wave propagation along the bedding planes). In the fractured rock, maximum velocity and minimum attenuation of P‐waves was seen at 90° to the fracture normal. After correction for the background anisotropy, the fractured sample velocity anisotropy was expressed in terms of Thomsen's weak anisotropy parameters ε, γ & δ. A theory of frequency‐dependent seismic anisotropy in porous, fractured, media was able to predict the observed effect of viscosity and bulk modulus on ε and δ in water‐ and glycerin‐saturated samples, and the higher ε and δ values in air‐saturated samples. Theoretical predictions of fluid independent γ are also in agreement with the laboratory observations. We also observed the predicted polarisation cross‐over in shear‐wave splitting for wave propagation at 45° to the fracture normal as fluid viscosity and bulk modulus increases.     

Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey

A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q⁻¹_eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q⁻¹_eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q⁻¹_eff= 52 ± 15) and slowest stacking velocity. The variation of Q⁻¹_eff with azimuth φ is well fitted by 1000Q⁻¹_eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 11⁴ combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10⁻⁵, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here.     

Converted-wave imaging in anisotropic media: theory and case studies

Common‐conversion‐point binning associated with converted‐wave (C‐wave) processing complicates the task of parameter estimation, especially in anisotropic media. To overcome this problem, we derive new expressions for converted‐wave prestack time migration (PSTM) in anisotropic media and illustrate their applications using both 2D and 3D data examples. The converted‐wave kinematic response in inhomogeneous media with vertical transverse isotropy is separated into two parts: the response in horizontally layered vertical transverse isotrophy media and the response from a point‐scatterer. The former controls the stacking process and the latter controls the process of PSTM. The C‐wave traveltime in horizontally layered vertical transverse isotrophy media is determined by four parameters: the C‐wave stacking velocity V_C2, the vertical and effective velocity ratios γ₀ and γ_eff, and the C‐wave anisotropic parameter χ_eff. These four parameters are referred to as the C‐wave stacking velocity model. In contrast, the C‐wave diffraction time from a point‐scatterer is determined by five parameters: γ₀, V_P2, V_S2, η_eff and ζ_eff, where η_eff and ζ_eff are, respectively, the P‐ and S‐wave anisotropic parameters, and V_P2 and V_S2 are the corresponding stacking velocities. V_P2, V_S2, η_eff and ζ_eff are referred to as the C‐wave PSTM velocity model. There is a one‐to‐one analytical link between the stacking velocity model and the PSTM velocity model. There is also a simple analytical link between the C‐wave stacking velocities V_C2 and the migration velocity V_Cmig, which is in turn linked to V_P2 and V_S2. Based on the above, we have developed an interactive processing scheme to build the stacking and PSTM velocity models and to perform 2D and 3D C‐wave anisotropic PSTM. Real data applications show that the PSTM scheme substantially improves the quality of C‐wave imaging compared with the dip‐moveout scheme, and these improvements have been confirmed by drilling.     

Determining Q of near-surface materials from Rayleigh waves

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 (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P 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 Q_P than to Q_S. 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.     

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.     

Analytic expressions for the velocity sensitivity to the elastic moduli for the most general anisotropic media

For non‐linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density‐normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS₁, qS₂). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP‐, qS₁‐ and qS₂‐wave computations and mitigates the singularity problem for the two quasi‐shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes.     

Experimental estimation of velocities and anisotropy of a series of Swedish crystalline rocks and ores

To provide a guide for future deep (<1.5 km) seismic mineral exploration and to better understand the nature of reflections imaged by surface reflection seismic data in two mining camps and a carbonatite complex of Sweden, more than 50 rock and ore samples were collected and measured for their seismic velocities. The samples are geographically from the northern and central parts of Sweden, ranging from metallic ore deposits, meta‐volcanic and meta‐intrusive rocks to deformed and metamorphosed rocks. First, ultrasonic measurements of P‐ and S‐wave velocities at both atmospheric and elevated pressures, using 0.5 MHz P‐ and S‐wave transducers were conducted. The ultrasonic measurements suggest that most of the measured velocities show positive correlation with the density of the samples with an exception of a massive sulphide ore sample that shows significant low P‐ and S‐wave velocities. The low P‐ and S‐wave velocities are attributed to the mineral texture of the sample and partly lower pyrite content in comparison with a similar type sample obtained from Norway, which shows significantly higher P‐ and S‐wave velocities. Later, an iron ore sample from the central part of Sweden was measured using a low‐frequency (0.1–50 Hz) apparatus to provide comparison with the ultrasonic velocity measurements. The low‐frequency measurements indicate that the iron ore sample has minimal dispersion and attenuation. The iron ore sample shows the highest acoustic impedance among our samples suggesting that these deposits are favourable targets for seismic methods. This is further demonstrated by a real seismic section acquired over an iron ore mine in the central part of Sweden. Finally, a laser‐interferometer device was used to analyse elastic anisotropy of five rock samples taken from a major deformation zone in order to provide insights into the nature of reflections observed from the deformation zone. Up to 10% velocity‐anisotropy is estimated and demonstrated to be present for the samples taken from the deformation zone using the laser‐interferometery measurements. However, the origin of the reflections from the major deformation zone is attributed to a combination of anisotropy and amphibolite lenses within the deformation zone.     

Experimental verification of the fracture density and shear‐wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry

We present laboratory ultrasonic measurements of shear‐wave splitting from two synthetic silica cemented sandstones. The manufacturing process, which enabled silica cementation of quartz sand grains, was found to produce realistic sandstones of average porosity 29.7 ± 0.5% and average permeability 29.4 ± 11.3 mD. One sample was made with a regular distribution of aligned, penny‐shaped voids to simulate meso‐scale fractures in reservoir rocks, while the other was left blank. Ultrasonic shear waves were measured with a propagation direction of 90° to the coincident bedding plane and fracture normal. In the water saturated blank sample, shear‐wave splitting, the percentage velocity difference between the fast and slow shear waves, of <0.5% was measured due to the bedding planes (or layering) introduced during sample preparation. In the fractured sample, shear‐wave splitting (corrected for layering anisotropy) of 2.72 ± 0.58% for water, 2.80 ± 0.58% for air and 3.21 ± 0.58% for glycerin saturation at a net pressure of 40 MPa was measured. Analysis of X‐ray CT scan images was used to determine a fracture density of 0.0298 ± 0.077 in the fractured sample. This supports theoretical predictions that shear‐wave splitting (SWS) can be used as a good estimate for fracture density in porous rocks (i.e., SWS = 100ε_f, where ε_f is fracture density) regardless of pore fluid type, for wave propagation at 90° to the fracture normal.     

Intrinsic Qp seismic attenuation from the rise time of microearthquakes: a local scale application at Rio‐Antirrio,Western Greece

We determine the attenuation structure of a three‐dimensional medium based on first pulse‐width measurements from microearthquake data. Ninety‐five microearthquakes from a seventy stations local network were considered in this study. Measurements of the first half cycle of the wave, the so‐called rise time τ were carried out on high quality velocity seismograms and inverted to estimate the P‐waves intrinsic quality factor Q_p. The results of this investigation indicate that first pulse width data from a local microearthquake network permit retrieval with sufficient accuracy of the heterogeneous Q_p structure. The inferred attenuation variability corresponds to the known geological formations in the region.     

Physical property anisotropy of foliated fault rocks: Study from the Nobeoka Thrust,Shimanto Belt,southwest Japan

To investigate the physical property anisotropies of foliated fault rocks in subduction zones, the hanging wall phyllites and footwall cataclasites exhumed along the Nobeoka Thrust, a fossilized out‐of‐sequence‐thrust in the Shimanto Belt, Japan, was focused. Discrete physical property (electric resistivity, P‐ and S‐wave velocities, and porosity) measurements were conducted employing geologic coordinates (depth‐parallel direction, strike direction, and maximum dip direction of foliation), using the core samples obtained from the Nobeoka Thrust Drilling Project and compared the data to borehole geophysical logs. A higher sample P‐wave velocity (Vp), lower S‐wave velocity (Vs), higher Vp/Vs, and lower sample porosity and resistivity compared to the logs, are inferred to have been caused by the larger sampling scale of the logs and lower fluid saturation of the borehole. The phyllites and cataclasites exhibited substantial vertical and horizontal anisotropy of Vp (0.4–17.3 % and 2.7–13.8 %, respectively), Vs (0.5–56 % and 7.7–43 %, respectively), and resistivity (0.9–119 % and 2.0–65.9 %, respectively). The physical property anisotropies are primarily affected by the dip angles of foliation. The fault rocks that have gentler dip angles exhibit a higher Vp in the strike and maximum dip direction and a lower Vp in the depth‐parallel direction. In contrast, the fault rocks that have steeply dipping structures show a higher Vp in the strike and depth‐parallel directions with a lower velocity in the maximum dip direction. Resistivity anisotropy show a trend opposite to that of the Vp in relation to the dip angles. Our results show lower Vp anisotropy than those obtained in previous studies, which measured wave speeds perpendicular or parallel to foliation under confining pressure. This study highlights the significance of dip angles on vertical properties in geophysical surveys across foliated fault rocks.     

Dispersion and radial depth of investigation of borehole modes

Sonic techniques in geophysical prospecting involve elastic wave velocity measurements that are performed by placing acoustic transmitters and receivers in a fluid‐filled borehole. The signals recorded at the receivers are processed to obtain compressional‐ and shear‐wave velocities in the surrounding formation. These velocities are generally used in seismic surveys for the time‐to‐depth conversion and other formation parameters, such as porosity and lithology. Depending upon the type of transmitter used (e.g. monopole or dipole) and as a result of eccentering, it is possible to excite axisymmetric (n= 0) , flexural (n= 1) and quadrupole (n= 2) families of modes propagating along the borehole. We present a study of various propagating and leaky modes that includes their dispersion and attenuation characteristics caused by radiation into the surrounding formation. A knowledge of propagation characteristics of borehole modes helps in a proper selection of transmitter bandwidth for suppressing unwanted modes that create problems in the inversion for the compressional‐ and shear‐wave velocities from the dispersive arrivals. It also helps in the design of a transmitter for a preferential excitation of a given mode in order to reduce interference with drill‐collar or drilling noise for sonic measurements‐while‐drilling. Computational results for the axisymmetric family of modes in a fast formation with a shear‐wave velocity of 2032 m/s show the existence of Stoneley, pseudo‐Rayleigh and anharmonic cut‐off modes. In a slow formation with a shear‐wave velocity of 508 m/s, we find the existence of the Stoneley mode and the first leaky compressional mode which cuts in at approximately the same normalized frequency ωa/V_S= 2.5 (a is the borehole radius) as that of the fast formation. The corresponding modes among the flexural family include the lowest‐order flexural and anharmonic cut‐off modes. For both the fast and slow formations, the first anharmonic mode cuts in at a normalized frequency ωa/V_S= 1.5 approximately. Cut‐off frequencies of anharmonic modes are inversely proportional to the borehole radius in the absence of any tool. The borehole quadrupole mode can also be used for estimating formation shear slownesses. The radial depth of investigation with a quadrupole mode is marginally less than that of a flexural mode because of its higher frequency of excitation.     

Near surface velocity and QS structure of the Quaternary sediment in Bohai basin,China

Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (Q_P and Q_S) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 M_W4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average v_P and v_S of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high v_P/v_S ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that Q_S is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with Q_S modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.     

Comparison of the attenuation properties for two different areas in eastern Anatolia,Turkey

In this study, the attenuation properties of the crust and the quality factor of S wave in eastern Anatolia (Turkey) were determined by local earthquakes for two different areas, Oltu and Erzurum. Seismic wave attenuation can be changed with high pressure or structural effects. Therefore, we argued that the estimation of attenuation coefficient in seismic active zones in Eastern Anatolia is a very useful tool to determine seismic activities. It uses regional waveform data set from two stations, OLT and ERZ, for 95 events that occurred in these regions between 2001 and 2005. The attenuation has been determined using the Chobra–Alexeev model based on the epicenter distance–amplitude relations. This model allows for investigation of the effects of variations in attenuation properties for different areas. We introduced a new magnitude formula for these areas using the amplitude normalization methods for reference values M_L=4, so as to correct effects of the magnitudes. We also determined velocity of seismic waves. The average attenuation coefficient (α), average quality factor (Q_s) and P and S waves velocities were obtained with normalized amplitude values for Erzurum (ERZ) and Oltu (OLT) as 0.0135 km⁻¹, 37, 6.20 km/s and 3.38 km/s and 0.0151, 34, 6.13 and 3.48.     

Frequency-Dependent Attenuation of Coda Waves in the Crust in Southwest Anatolia (Turkey)

The attenuation of coda waves in the earth’s crust in southwest (SW) Anatolia is estimated by using the coda wave method, which is based on the decrease of coda wave amplitude in time and distance. A total of 159 earthquakes were recorded between 1997 and 2010 by 11 stations belonging to the KOERI array. The coda quality factor Q _c is determined from the properties of scattered coda waves in a heterogeneous medium. Firstly, the quality factor Q ₀ (the value of Q _c at 1 Hz.) and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves for frequencies of 1.5, 3.0, 6.0, 8.0, 12 and 20 Hz. Secondly, the attenuation coefficients (δ) are estimated. The shape of the curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The average Q _c values vary from 110 ± 15 to 1,436 ± 202 for the frequencies above. The Q ₀ and η values vary from 63 ± 7 to 95 ± 10 and from 0.87 ± 0.03 to 1.04 ± 0.09, respectively, for SW Anatolia. In this region, the average coda Q–f relation is described by Q _c = (78 ± 9)f ^0.98±0.07 and δ = 0.012 km^?1. The low Q ₀ and high η are consistent with a region characterized by high tectonic activity. The Q _c values were correlated with the tectonic pattern in SW Anatolia.     

Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.