由多个震源机制解分析民乐Ms6.1级地震前后应力场演化

2003年10月25日在甘肃省民乐县和山丹县境内发生了Ms6．1与Ms5．8级地震，对震区200km范围、震前2年内发生的地震，及对震后部分余震矩张量反演得到它们的震源机制解，并对其进行统计分析得到该区的背景应力场及应力轴的整体特征，根据发震区地质构造背景和对平均应力场变化分析，寻求民乐地震的发震机制。     

Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling

Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (M_W = 5.8), and the large August 30th, 1986 (M_W = 7.1) and March 4th, 1977 (M_W = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ M_W ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

利用地基土应力分布原理确定砂垫层厚度

软弱土层场地若采用天然地基不能满足建筑物所需强度与沉降要求,采用换土撼砂法就可以提高地基的承载力。运用应力分布原理合理确定撼砂层厚度不仅可保证工程稳固性,并可节约投资和缩短工期。     

A hybrid discrete-continuum approach to model hydro-mechanical behaviour of jointed rocks

A hybrid discrete-continuum approach has been presented in this paper to simulate water flow in the near and far fields of deformable fractured rocks. In the near field, the discrete model is used; while in the far field, the equivalent continuum model is employed. The discrete element method (with the static relaxation algorithm) is used in the near field and the boundary integral equation method in the far field. Along the interface of these two domains, both mechanical and hydraulic compatibility conditions are satisfied. Fully coupled hydro-mechanical analysis can be conducted in the combined near and far fields. Application to a dam foundation problem has demonstrated the capability of the developed approach.     

东北地区5 000 aB.P.～4 700 aB.P.左右的降温事件及对考古文化的影响

距今5000年左右世界各地普遍存在一次幅度较大的短期降温气候事件，通过对金川剖面泥炭地层14C测年，孢粉分析，认为东北地区此次降温事件结束的14C年龄在4700aB.P.左右，降温事件出现之前气候温暖湿润，降温事件过后温度持续下降，金川剖面的记录与东北地区其它类型沉积物所记录的降温事件出现的时间基本吻合，降温事件对燕辽地区古文化演变产生了深刻的影响。     

Earthquake strong motion simulated by stress drop models

Shallow strike slip earthquakes on vertical faults are modelled as two-dimensional antiplane strain ruptures in a uniformly prestressed homogeneous halfspace. Behind the rupture front, which is specified, the stress drops to a lower value. The elastodynamic boundary value problem is solved with a finite difference method. Several cases are studied, which include symmetric and one-directional rupture propagation, surface faulting, multiple events, variable rupture velocity, sticking and rebreaking of the fault plane. The time function of displacement, velocity and acceleration are interpreted as signals generated by events in the focus, namely starting, stopping and breaking through the surface of the rupture. The model explains peak velocity and peak acceleration in the near field of M5.5–6 earthquakes; which are typically about 0.2 m/s and 5 m/s² at 10 km epicentral distance, if the rupture velocity is close to the shear wave velocity. Sticking of the fault does not alter the accelerograms significantly, but it increases the seismic moment in simple events and decreases it in multiple events.Contribution No. 226, Geophysical Institute, University of Karlsruhe.     

The coalescence of water drops II. The coalescence problem and coalescence efficiency

By the use of the model of approaching drops (Arbel and Levin, 1977) the coalescence efficiencies of drops are computed. It is found that for interactions of drops at their terminal velocities the coalescence depends both on the size of the large drop and on the size ratio of the interacting drops in agreement with the experimental results of Whelpdale and List (1971) and Levin and Machnes (1977).The results were found to be sensitive to the assumption of the drops deformation and to the critical separation distance. This distance is defined as the distance at which the drops begin to merge. The variations of the coalescence efficiency with these parameters is discussed.Appendix: List of symbols D distance between the deformed surfaces of the drops - D _o initial value ofD - D _s stop distance, the distance at which the impact velocity vanishes - D _c critical coalescence distance - E collection efficiency - E ₁ collision efficiency - E ₂ coalescence efficiency - E _2R coalescence efficiency for collisions with stationary targets - F _c centrifugal force - p ratio of the radii of the interacting drops - r _o initial distance between drops' centers - R _L radius of larger drop - R _s radius of smaller drop - R _D radius of deformation - v approach velocity of two deformed surfaces - v _o initial value ofv - V _i impact velocity (given negative sign when drops approach each other) - V _c critical impact velocity - W _i velocity of the smaller drop at infinity for it to reachD _o with velocityv _o - x _i impact distance, the distance between the trajectories of the two drops - x _c critical impact distance for coalescence -  average critical impact distance for coalescence - X _c critical impact distance for collisions - coefficient of deformation given in equation 1 - _i impact angle defined byWhelpdale andList (1971) given also inArbel andLevin (1977) - coefficient of deformation given in equation 2 - viscosity of air - _i impact angle used inArbel andLevin (1977) and here - _c critical angle for coalescence - average critical angle for coalescence On sabbatical leave (1976–77) from the Department of Geophysics and Planetary Sciences, Tel Aviv University, Ramat Aviv, Israel.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Stress in the lithosphere: Inferences from steady state flow of rocks

Mechanical data and flow processes from steady state deformation experiments may be used to infer the state of stress in the lithosphere and asthenosphere. Extrapolations of flow equations to a representative geologic strain rate of 10^–14/sec. for halite, marble, quartzite, dolomite, dunite and enstatolite are now warranted because the steady state flow processes in the experiments are identical to those in rocks and because the geotherms are reasonably well established. More direct estimates are obtained from free dislocation densities, subgrain sizes and recrystallized grain sizes all of which are functions only of stress. Using the last of these techniques, we have estimated stress profiles as a function of depth from xenoliths in basalts and kimberlites, whose depths of equilibration were determined by pyroxene techniques, from four different areas of subcontinental and suboceanic upper mantle. The results are similar and indicate stress differences of about 200 to 300 bars at 40 to 50 km, decaying to a few tens of bars at depths betow 100 km. These stresses are reasonable and are in accord with extrapolations of the mechanical data provided that allowance is made for a general increase in strain rate and decrease in viscosity with depth.     

Shear and tension hydraulic fractures in low permeability rocks

Laboratory hydrofracture experiments were performed on triaxially stressed specimens of oil shale and low-permeability granite. The results show that either shear or tension fractures could develop depending on the level of differentials stress, even in specimens containing preexisting fractures. With 1 kb of confining pressure and differential stress greater than 2kb, hydraulic fluid diffusion into the specimens reduced the effective confining pressure until failure occurred by shear fracture. Below 2kb of differential stress, tension fractures occurred. These results suggest that hydraulic fracturing in regions of significant tectonic stress may produce shear rather than tension fractures. In this casein situ stress determinations based on presumed tension fractures would lead to erroneous results.