浅层地震技术在隧道勘察中的应用

论述了在深圳市东部沿海高速公路隧道勘察中,应用浅层地震勘探方法查明地下地质构造、风化带、分层等,并根据资料解释成果、钻探资料和工区地质资料,对隧道开挖过程中可能遇到的断层及塌陷等情况进行了讨论。     

Constraints from melt inclusions and their host olivines on the petrogenesis of Oligocene-Early Miocene Xindian basalts, Chifeng area, North China Craton

The geochemical characteristics of melt inclusions and their host olivines provide important information on the processes that create magmas and the nature of their mantle and crustal source regions. We report chemical compositions of melt inclusions, their host olivines and bulk rocks of Xindian basalts in Chifeng area, North China Craton. Compositions of both bulk rocks and melt inclusions are tholeiitic. Based on petrographic observations and compositional variation of melt inclusions, the crystallizing sequence of Xindian basalts is as follows: olivine (at MgO > ~5.5 wt%), plagioclase (beginning at MgO = ~5.5 wt%), clinopyroxene and ilmenite (at MgO < 5.0 wt%). High Ni contents and Fe/Mn ratios, and low Ca and Mn contents in olivine phenocrysts, combining with low CaO contents of relatively high MgO melt inclusions (MgO > 6 wt%), indicate that Xindian basalts are possibly derived from a pyroxenite source rather than a peridotite source. In the CS-MS-A diagram, all the high MgO melt inclusions (MgO > 6.0 wt%) project in the field between garnet + clinopyroxene + liquid and garnet + clinopyroxene + orthopyroxene + liquid near 3.0 GPa, further suggesting that residual minerals are mainly garnet and clinopyroxene, with possible presence of orthopyroxene, but without olivine. Modeling calculations using MELTS show that the water content of Xindian basalts is 0.3–0.7 wt% at MgO = 8.13 wt%. Using 20–25 % of partial melting estimated by moderately incompatible element ratios, the water content in the source of Xindian basalts is inferred to be ≥450 ppm, much higher than 6–85 ppm in dry lithospheric mantle. The melting depth is inferred to be ~3.0 GPa, much deeper than that of tholeiitic lavas (<2.0 GPa), assuming a peridotite source with a normal mantle potential temperature. Such melting depth is virtually equal to the thickness of lithosphere beneath Chifeng area (~100 km), suggesting that Xindian basalts are derived from the asthenospheric mantle, if the lithospheric lid effect model is assumed.     

Analysis of mass concentration of atmospheric particulate matter in a sandstorm course and its affecting factors in the Taklimakan Desert

During the course of a major sandstorm from April 17 to April 23, 2008 in the Taklimakan Desert, data pertaining to the mass concentrations of different-sized atmospheric particulate matter were observed continuously with Grimm 1.108, Thermo RP 1400a, TSP, and CAWS-600 instruments. The results showed that: (1) during the entire sandstorm process there were some differences between the daily mean particle concentration peaks and the hourly mean particle concentration peaks because the actual sandstorm lasted for only about 4 hr, whereas more particles were accumulated in the floating dust days before and after the actual sandstorm; (2) the intensity of the sandstorm was enhanced with the increase of wind speed, and this was related to the peak mass concentrations of atmospheric particulate matter; the wind speed directly affected the concentration of atmospheric particulate matter: the higher the wind speed, the higher the mass concentration (>0.23 μm was 39,496.5 μg/m³, and >20.0 μm was 5,390.7 μg/m³); (3) the concentration changes of PM₁₀ and TSP were also related to the course and intensity of the sandstorm; and (4) the mass concentration of atmospheric particulate matter had the following sequence during the dust weather: clear day < floating dust < floating and blowing dust < sandstorm. Temperature, relative humidity, and barometric pressure are important factors affecting the strength of storms, which could also indirectly affect the concentration change of atmospheric particulate matter.     

Comparison of spectral and spatial methods for a Moho recovery from gravity and vertical gravity-gradient data

In global studies investigating the Earth’s lithospheric structure, the spectral expressions for the gravimetric forward and inverse modeling of the global gravitational and crustal structure models are preferably used, because of their numerical efficiency. In regional studies, the applied numerical schemes typically utilize the expressions in spatial form. Since the gravity-gradient observations have a more localized support than the gravity measurements, the gravity-gradient data (such as products from the Gravity field and steady-state Ocean Circulation Explorer - GOCE - gravity-gradiometry satellite mission) could preferably be used in regional studies, because of reducing significantly the spatial data-coverage required for a regional inversion or interpretation. In this study, we investigate this aspect in context of a regional Moho recovery. In particular, we compare the numerical performance of solving the Vening Meinesz-Moritz’s (VMM) inverse problem of isostasy in spectral and spatial domains from the gravity and (vertical) gravity-gradient data. We demonstrate that the VMM spectral solutions from the gravity and gravity-gradient data are (almost) the same, while the VMM spatial solutions differ from the corresponding spectral solutions, especially when using the gravity-gradient data. The validation of the VMM solutions, however, reveals that the VMM spatial solution from the gravity-gradient data has a slightly better agreement with seismic models. A more detailed numerical analysis shows that the VMM spatial solution formulated for the gravity gradient is very sensitive to horizontal spatial variations of the vertical gravity gradient, especially in vicinity of the computation point. Consequently, this solution provides better results in regions with a relatively well-known crustal structure, while suppressing errors caused by crustal model uncertainties from distant zones. Based on these findings we argue that the gravity-gradient data are more suitable than the gravity data for a regional Moho recovery.     

Assessment of shallow landslides from Hurricane Mitch in central America using a physically based model

Shallow landslides induced by heavy rainfall events represent one of the most disastrous hazards in mountainous regions because of their high frequency and rapid mobility. Recent advancements in the availability and accessibility of remote sensing data, including topography, land cover and precipitation products, allow landslide hazard assessment to be considered at larger spatial scales. A theoretical framework for a landslide forecasting system was prototyped in this study using several remotely sensed and surface parameters. The applied physical model SLope-Infiltration-Distributed Equilibrium (SLIDE) takes into account some simplified hypotheses on water infiltration and defines a direct relation between factor of safety and the rainfall depth on an infinite slope. This prototype model is applied to a case study in Honduras during Hurricane Mitch in 1998. Two study areas were selected where a high density of shallow landslides occurred, covering approximately 1,200 km². The results were quantitatively evaluated using landslide inventory data compiled by the United States Geological Survey (USGS) following Hurricane Mitch’s landfall. The agreement between the SLIDE modeling results and landslide observations demonstrates good predictive skill and suggests that this framework could serve as a potential tool for the future early landslide warning systems. Results show that within the two study areas, the values of rates of successful estimation of slope failure locations reached as high as 78 and 75%, while the error indices were 35 and 49%. Despite positive model performance, the SLIDE model is limited by several assumptions including using general parameter calibration rather than in situ tests and neglecting geologic information. Advantages and limitations of this physically based model are discussed with respect to future applications of landslide assessment and prediction over large scales.     

Rapid estimation of maximum and minimum void ratios of granular soils

Maximum and minimum void ratios (e_max and e_min) of granular soils are commonly used as indicators of many engineering properties. However, few methods, apart from laboratory tests, are available to provide a rapid estimation of both e_max and e_min. In this study, we present a theoretical model to map the densest and the loosest packing configurations of granular soils onto the void space. A corresponding numerical procedure that can predict both e_max and e_min of granular soils with arbitrary grain size distributions is proposed. The capacity of the proposed method is evaluated by predicting the maximum and minimum void ratios of medium to fine mixed graded sands with different contents of fines. The influence of the grain size distribution, characterized quantitatively by uniformity parameter and the fractal dimension, on e_max and e_min is discussed using the proposed method. Moreover, application of this method in understanding the controlling mechanism for the void ratio change during grain crushing is presented.

     

Temporal vP/vS variation characteristics in different zones of China’s Capital Circle area before and after Wen’an earthquake

On July 4, 2006, an earthquake of MS5.1 took place in Wen’an, Hebei Province, just at the south center of China’s Capital Circle area digital seismograph network. It is the strongest event recorded ever since the network went into operation in 2002. We processed the vast amounts of phase data yielded by the 107 digital seismic stations between 2002~2007 using Wadati method. In order to improve the precision and stability of shear and compressional wave velocities (vP/vS) calculation, we impose a number of restrictions on the computation environment and condition, e.g., the earthquakes are densely concentrated, selected stations are limited in range, the number of stations in- volved in the computation is larger than 5 and linear fitting features high precision and small error. Under these restrictions, the study shows that vP/vS in and around Wen’an and Tangshan underwent a normal-low-normal proc- ess one year before Wen’an earthquake, vP/vS became obviously low and the low ratio lasted for about one year, meanwhile, little variation of vP/vS was seen in Xingtai, northwest of Beijing, southwest of Beijing, Beijing-Tianjin and Beijing; after the quake, the vP/vS returned normal in Wen’an and Tangshan. Error and stability analysis of the calculated result for vP/vS shows it is convincible that anomaly appeared in and around Wen’an and Tangshan be- fore Wen’an earthquake.     

1989 Datong-Yanggao Earthquake with Ms=6.1

In this paper,the basic conditions and characteristics of the Datong-Yanggao earthquake are described in detail.The type of the earthquake,the aftershock activity,the earthquake trend,the characteristics of damage are analyzed and discussed.This earthquake belongs to a sequence of earthquake swarm type and is of the concentratedly occurring characteristics of several main shocks in a time sequence,and the aftershock activity was normally attenuated.Because of the unreasonable design of buildings and repeated vibration by several destructive main shocks within a short period,the damage in the high intensity region was aggravated.The analysis and the prediction procedure of the 1991 Datong-Yanggao earthquake(Ms=5.8)are also briefly given.     

Meteoric diagenesis in Pleistocene reef limestones of Xisha Islands, China

This paper focuses on a borehole, Xichen-1 well, drilled on the Chenhang Island, Xisha Islands in the South China Sea. Mineralogical, petrographic, stable isotopic and minor-element data from the Holocene to Pleistocene interval (0–179 m ) in the Xichen-1 well are discussed in detail. The 400-m-long core is divisible into four mineralogical facies: a high-Mg calcitic aragonite facies (0–16.91 m, Holocene), an aragonitic low- Mg calcite facies (16.91–30.60 m, Late Pleistocene), a low- Mg calcitic facies (30.60–179 m, Middle-Early Pleistocene) and a low- Mg calcitic and dolomitic facies (179–400 m, Early Pleistocene–Late Miocene). The Holocene section has much higher whole-rock δ¹⁸O and δ¹³C values and Mg and Sr content than the non-dolomitized Pleistocene limestones (16.91–179 m). The 16.91–165 m interval is characterized by a relatively invariant oxygen isotopic composition and very heterogeneous carbon isotopic composition. Between 165 and 179 m, there is a positively correlated increase of whole-rock δ¹⁸O and δ¹⁸C with depth, and Mg content also shows a gradual increase with depth. Petrographic data demonstrate that the Pleistocene reef sequence has been extensively affected by meteoric waters. We conclude that the Late Pleistocene section (16.91–30.60 m) and the Middle-Early Pleistocene section (30.60–165 m) have suffered incomplete and complete meteoric diagenesis, respectively, and that the Early Pleistocene interval (165–179 m) was diagenetically altered in a meteoric–marine mixing environment.     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.