Identification of Mining Blasts at Mid- to Far-regional Distances Using Low Frequency Seismic Signals

—?This paper reports results from two recent monitoring experiments in Wyoming. Broadband seismic recordings of kiloton class delay-fired cast blasts and instantaneous calibration shots in the Black Thunder coal mine were made at four azimuths at ranges from 1° to 2°. The primary focus of this experiment was to observe and to explain low-frequency signals that can be seen at all azimuths and should routinely propagate above noise to mid-regional distances where most events will be recorded by International Monitoring System (IMS) stations.¶The recordings clearly demonstrate that large millisecond delay-fired cast blasts routinely produce seismic signals that have significant spectral modulations below 10?Hz. These modulations are independent of time, the azimuth from the source and the orientation of the sensor. Low-frequency modulations below 5?Hz are seen beyond 9°. The modulations are not due to resonance as they are not produced by the calibration shots. Linear elastic modeling of the blasts that is guided by mine-blast reports fails to reproduce the fine detail of these modulations but clearly indicates that the enhanced “spectral roughness” is due to long interrow delays and source finiteness. The mismatch between the data and the synthetics is likely due to source processes, such as nonlinear interactions between shots, that are poorly understood and to other effects, such as variations of shot time and yield from planned values, that are known to be omnipresent but cannot be described accurately. A variant of the Automated Time-Frequency Discriminant (Hedlin, 1998b), which uses low-frequency spectral modulations, effectively separates these events from the calibration shots.¶The experiment also provided evidence that kiloton class cast blasts consistently yield energetic 2–10 second surface waves. The surface waves are strongly dependent on azimuth but are seen beyond 9°. Physical modeling of these events indicates that the surface waves are due mainly to the extended source duration and to a lesser extent to the slap-down of spalled material. The directionality is largely a path effect. A discriminant that is based on the partitioning of energy between surface and body waves routinely separates these events from the calibration shots.¶The Powder River Basin has essentially no natural seismic activity. How these mining events compare to earthquake observations remains to be determined.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

Mantle Convection Modeling with Viscoelastic/Brittle Lithosphere: Numerical Methodology and Plate Tectonic Modeling

--The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.     

3-D Viscoelastic FEM Modeling of Crustal Deformation in Northeast Japan

-- As a first step toward establishing a standard earthquake cycle model in Japan, we simulate the crustal deformation during the past 100 years in northeast Japan, using a 3-D FEM based on the kinematic model. Then, we compare the computed results with the observed long-term leveling data and the recent GPS data. On the whole, although the effect of the subducting PAC is dominant, coseismic deformation of the interplate earthquakes can be clearly seen in the inland. Moreover, the postseismic deformation of the earthquakes due to the viscoelastic upper mantle seriously affects the inland movements, and continues for a few decades. Our modeling, including the effects of the interplate earthquakes and the three-dimensional viscoelastic inhomogeneity, reasonably explains the observed movement. Finally, we stress that the viscoelastic effect should be taken into consideration in the analyses, even if no earthquakes occur in the analyzed period.     

Evaluation of Subsurface Structure at Soultz Hot Dry Rock Site by the AE Reflection Method in Time-frequency Domain

—?Because the Soultz Hot Dry Rock (HDR) site, France, is to be expanded to a scientific pilot plant of greater depth, measurement of the deep area below the predeveloped artificial reservoir is gaining importance. In this paper, we present estimates of deep subsurface structure at the Soultz HDR site, obtained by a reflection method using acoustic emission (AE) signals, that is, induced seismicity, as the wave source. First, we briefly describe the AE reflection method in the time-frequency domain with wavelet transform. Then we show estimates of the subsurface structure by using 101 high-energy AE events observed in 1993. We compare and discuss the results obtained, using the AE reflection method in two wells with other independent borehole observations. Furthermore, we present the results of an investigation of the frequency dependence of reflectors identified by hodogram linearity as a possible means of further characterizing detected structures.     

Microseismicity Induced by Heavy Rainfall Around Flooded Vertical Ore Veins

—?Microseismicity (M?M?M?M??6 were also monitored at a crustal movement monitoring station located several hundred meters from the veins. It was found that the opening of the vertical ore veins primary led to significant strain and tilt, but not to seismicity, because the delay and the longer duration of the seismicity were significant. Most seismic events involve thrusting mechanisms that are consistent with the present stress state of E-W-oriented tectonic compression, but are not consistent with the opening of the deepest ore vein. Interstingly, all the events within a few months of the heavy rainfall occurred near the faults that offset the deepest ore veins, wheareas all those events located away from the deepest ore veins occurred many months after the heavy rainfall. Consequently, the delayed diffusion of water appears to have played a dominant role in reducing rock strength, which led to seismicity in the Ikuno mine.     

Modeling of Wide-Angle Seismic Attributes Using the Gaussian Beam Method

In this study, observed seismic attributes from shot gather 11 of the SAREX experiment are used to derive a preliminary velocity and attenuation model for the northern end of the profile in southern Alberta. Shot gather 11 was selected because of its prominent Pn arrivals and good signal to noise ratio. The 2-D Gaussian beam method was used to perform the modeling of the seismic attributes including travel times, peak envelope amplitudes and pulse instantaneous frequencies for selected phases. The preliminary model was obtained from the seismic attributes from shot gather 11 starting from prior tomographic results. The amplitudes and instantaneous frequencies were used to constrain the velocity and attenuation structure, with the amplitudes being more sensitive to the velocity gradients and the instantaneous frequencies more sensitive to the attenuation structure. The resulting velocity model has a velocity discontinuity between the upper and lower crust, and lower velocity gradients in the upper and lower crust compared to earlier studies. The attenuation model has Q _p ^-1 values between 0.011 and 0.004 in the upper crust, 0.0019 in the lower crust and a laterally variable Q _p ^-1 in the upper mantle. The Q _p ^-1 values are similar to those found in Archean terranes from other studies. Although the results from a single gather are non-unique, the initial model derived here provides a self-consistent starting point for a more complete seismic attribute inversion for the velocity and attenuation structure.     

Some Aspects of Resource Uncertainty and Their Economic Consequences in Assessment of the 1002 Area of the Arctic National Wildlife Refuge

Natural Resources Research - Exploration ventures in frontier areas have high risks. Before committing to them, firms prepare regional resource assessments to evaluate the potential payoffs. With...     

A method to reduce the influence of ice-rafted debris on a grain size record from northern Fram Strait, Arctic Ocean

The granulometric composition of terrigenous deep-sea sediments provides information on current speed if certain frame conditions are fulfilled. These include that current transport is the only transport process. At high latitudes this type of investigation is impaired due to the influence of ice-rafted debris (IRD) which contaminates the current-sorted grain size fractions. This study presents a new method that addresses this problem by setting the ice-rafted sand in relation to the silt of both current- and ice-transported origin. Deviations from the resulting regression function are then used to determine the behaviour of the silt mean grain size as a function of current speed largely independent from IRD bias. The study is based on sediments from the Yermak Plateau, Arctic Ocean, a region influenced by IRD brought with the south-headed Transpolar Drift and by north-directed bottom currents. The IRD correction results in displaying changes of current speed at much higher clarity; climate forcing of the currents becomes more evident. For example, the 8200 year cold event shows up as a major event in the corrected record whereas it is hardly visible in the original record.