Regional hazard prediction of rock bursts using microseismic energy attenuation tomography in deep mining

Rock burst prediction is a worldwide challenge that we have long tried to overcome. This study tentatively proposed a method to regionally predict rock burst hazards using microseismic energy attenuation. To verify the feasibility of the proposal, first, the mechanism of microseismic energy propagation and attenuation in rock medium was explored, and dominant attenuation characteristics of microseismic waves were analyzed. Second, a spatial attenuation model of microseismic energy was established, and the average energy attenuation coefficient for each wave path was defined. A 3D seismic energy attenuation inversion algorithm was put forward, and the corresponding computation matrix was developed. Third, a continuous microseismic field investigation was carried out in a deep coal mine. Seismic energy attenuation coefficient was confirmed using the calibrated focus position and energy determination. Based on data discretization processing, energy attenuation inversion and tomography, potential rock burst hazard regions were strictly zoned in mining areas. Finally, regional prediction results obtained from the microseismic energy attenuation were compared with the direct measurement results obtained from the classical drilling dust method to verify the reliability of proposed approach. It turns out that rock burst hazard regions predicted by the microseismic energy attenuation agreed well with the objective hazardous situations. Seismic energy attenuation coefficient is a significant evaluation factor that directly mirrors the inelastic performance of rock medium. Energy attenuation coefficient threshold used for determining the rock burst hazard regions was 3.0 km^?1. Reliability of the seismic energy attenuation inversion and tomography was closely related to the spatial distribution of microseisms in a localized region. The optimum spatial density of microseisms was 0.2 m^?3. Regional rock burst prediction using microseismic energy attenuation is an effective approach for revealing potential hazardous regions in deep mining conditions. This approach improves the pertinence of geological hazard prevention and provides a beforehand reference for targeted hazard management.     

Multiple-tracers-aided surface-subsurface hydrological modeling for detailed characterization of regional catchment water dynamics in Kumamoto area,southern Japan

Hydrogeology Journal - Integrated watershed modeling techniques have been applied in recent years to examine surface and subsurface interactions. Model performance is often evaluated by best fit of...     

El Niño flood deposits at Batán Grande,northern Peru

Recent construction of a 10 × 3 irrigation canal across the Batán Grande–Poma archaeological complex of northern Peru provided an opportunity for inspection of soil profiles and detailed sampling along a five km transect bisecting an important concentration of ceremonial structures built during the past 3000+ years. From laboratory analysis of 80 soil samples we conclude no long-term irrigation agriculture was practiced in this area and it was unimportant as a resource base. Stratigraphic studies indicate few sediments survive from the 1925 El Niño event, but a massive bed attributable to prehistoric slack-water flood sediments has been identified and can be bracketed between ca. A. D. 650–1000 by associated diagnostic funerary ceramics. Other flood deposits representing the same or different events are described, but we conclude intervalley correlations will be difficult until more precise dating methods are available.     

Variations of Kuroshio geostrophic transport south of Japan estimated from long-term IES observations

Two inverted echo sounders were maintained on coastal and offshore sides of the Kuroshio south of Japan from October 1993 to July 2004. Applying the gravest empirical mode method, we obtained a time series of geostrophic transport. Estimated transports generally agree well with geostrophic transports estimated from hydrography. Their agreement with the hydrographic transports is better than that of transports estimated from satellite altimetry data. The geostrophic transport is expressed as the surface transport per unit depth multiplied by the equivalent depth. The geostrophic transport varies mostly with the surface transport and fractionally with the equivalent depth. Seasonal variation of the geostrophic transport has a minimum in March and a maximum in September, with a range of about one fifth of the total transport.     

The interanual rainfall variability in Indonesia corresponding to El Niño Southern Oscillation and Indian Ocean Dipole

The Impact of the Indian Ocean Dipole(IOD) and the El Ni?o Southern Oscillation(ENSO) event for Indonesian rainfall has been investigated for the period from 1950 to 2011. Inter-annual change of IOD and ENSO indices are used to investigate their relationship with Indonesian rainfall. By using the wavelet transform method, we found a positive significant correlation between IOD and Indonesian rainfall on the time scale of nearly 2.5–4 years.Furthermore, the positive significant correlation between ENSO(sea surface temperature anomaly at Ni?o3.4 area indices) and Indonesian rainfall exists for shorter than 2 years and between 5.5 to 6.5-year time scales.     

Vertical double silicate maxima in the sea-ice reduction region of the western Arctic Ocean: Implications for an enhanced biological pump due to sea-ice reduction

The R/V Mirai conducted hydrographic surveys in the western Arctic Ocean during summer 2004 (Mirai04) over wide east-west ranges from Alaska to eastern Siberia, where sea-ice cover has been greatly reduced in recent summers. The obtained data reveal differences in silicate profiles between shelf slope areas east and west of the Chukchi Plateau, the ridge that divides the Canada Basin into the Alaskan and east Siberian sides. East of the plateau, a single silicate maximum was found in a layer of Pacific-origin winter water, as examined in many previous studies. In contrast, west of the plateau, we found vertical double silicate maxima, which are reported for the first time in this study. The shallower silicate maximum corresponded to an N** minimum, signaling denitrification at the shelf bottom. This suggests that the shallower silicate maximum was caused by the spreading of shelf water. In contrast, the deeper silicate maximum corresponded to an oxygen minimum and a maximum silicate/phosphate ratio (Si/P), suggesting that this deeper maximum resulted from the decomposition of opal-shelled organisms. We also compared a silicate profile from Mirai04 to aprofile from the Arctic Ocean Section 94 (AOS94) expedition of 1994, a heavy ice year. The results suggest that sea-ice loss has enhanced biological activities, likely resulting in the appearance of the deeper silicate maximum.     

Least squares and integral methods for the spherical harmonic analysis of theSq-field

Spherical harmonic coefficients (SHCs) for the daily magnetic variation fields (solar and lunar) and the main field of the earth are usually estimated by the method of least squares applied to a truncated spherical harmonic series. In this paper, an integral method for computing the SHCs for the solar quiet daily magnetic variation fieldSq is described and applied toSq data for May and June 1965. TheSq SHCs thus derived are then compared with the results obtained using both unweighted and weighted versions of the least squares method. The weighting used tends to orthogonalize the least squares terms. The integral and weighted least squares results agree closely for terms up to order 4 and degree 30, but both disagree considerably for the higher degree terms with the results of the unweighted least squares. Errors introduced by the numerical integration can be shown to be small, hence the disagreement between integral and unweighted least squares coefficient sets arises from improper weighting. Also, it is concluded that discrepancies between the geomagnetic northward and eastward component-derived coefficient sets arise from either time-dependent external sources that produce non-local-time, based fields or nonpotential sources and not from truncation of the spherical harmonic series as has previously been suggested.Deceased.     

Fracture strength of dry silicate rocks at high confining pressures and activity of acoustic emission

Three dry silicate rocks, gabbro, dunite and eclogite, were triaxially compressed up to a confining pressure of 3 GPa at room temperature. These rocks exhibited brittle fracture behavior up to the highest confining pressure. The change of the mechanism of fracture in the brittle region is suggested from the measurement of the compressive fracture strength and the activity of acoustic emission. The existence of the “high-pressure brittle-fracture” phase is proposed. The fracture strength increased with increase of confining pressure. The increasing rate of strength was lowered at a value of confining pressure: at about 0.8 GPa on gabbro; at about 1.0 GPa on dunite; and at about 1.5 GPa on eclogite. At lower confining pressures than the above value, the acoustic emission rate began to increase at the onset of dilatancy and increased rapidly followed by fracture as the axial stress was increased. At the higher confining pressures, however, the acoustic emission rate did not increase rapidly before final fracture, and stayed constant to the fracture. The similar behavior was shown on the granite studied previously. It is interesting that the frictional strength forms the boundary between “low- and high-pressure brittle-fracture” phases.     

Crustal structure in and around the region of the 1995 Kobe Earthquake deduced from a wide-angle and refraction seismic exploration

Abstract The 1995 Kobe (Hyogo-ken Nanbu) earthquake (M_JMA 7.2, M_w 6.9) occurred on Jan. 17, 1995, at a depth of 17 km, beneath the areas of southern part of Hyogo prefecture and Awaji Island. To investigate P-wave velocity distribution and seismological characteristics in the aftershock area of this great earthquake, a wide-angle and refraction seismic exploration was carried out by the Research Group for Explosion Seismology (RGES) . The profile including 6 shot points and 205 observations was 135 km in length, extending from Keihoku, Northern Kyoto prefecture, through Kobe, to Seidan on Awaji Island. The charge of each shot was 350–700 kg. The P-wave velocity structure model showed a complicated sedimentary layer which is shallower than 2.5 km, a 2.5 km-thick basement layer whose velocity is 5.5 km/s, overlying the crystalline upper crust, and the boundary between the upper and lower crust.
Almost all aftershock hypocenters were located in the upper crust. However, the structure model suggests that the hypocenters of the main shock and some aftershock clusters were situated deeper than the boundary between the upper and lower crust. We found that the P-velocity in the upper crust beneath the northern part of Awaji Island is 5.64 km/s which is 3% lower than that of the surrounding area. The low-velocity zone coincides with the region where the high stress moment release was observed.