现场大尺度摩擦实验中观测到的断层泥中的塑性变形波

本文介绍了在断层区的岩体中进行的现场大尺度摩擦实验。实验的目的是探索地壳中变形波的存在。实验结果表明,在加载时断层泥屈服时间不一致,加载一端的断层泥先屈服,然后屈服范围逐步扩大,此过程相应于塑性变形波在断层泥中的传播。实验测得试件中塑性变形波波速约为2—3 cm/min,实际相当于10—16km/year。     

A thin-layer interface model for wave propagation through filled rock joints

The present study essentially employs a thin-layer interface model for filled rock joints to analyze wave propagation across the jointed rock masses. The thin-layer interface model treats the rough-surfaced joint and the filling material as a continuum medium with a finite thickness. The filling medium is sandwiched between the adjacent rock materials. By back analysis, the relation between the normal stress and the closure of the filled joint are derived, where the effect of joint deformation process on the wave propagation through the joint is analyzed. Analytical solutions and laboratory tests are compared to evaluate the validity of the thin-layer interface model for filled rock joints with linear and nonlinear mechanical properties. The advantages and the disadvantages of the present approach are also discussed.     

地震作用下深埋隧道围岩形变的数值模拟研究

隧道围岩受地震荷载作用影响产生开裂变形,给交通、安全与经济带来极大影响,为提升隧道安全稳定性能,以某市隧道为对象,研究地震作用下深埋隧道围岩形变的数值模拟。综合考虑地震动作用与地震波扩散特性,利用FLAC软件构建三维动力模型;通过边界设置解决地震波反射问题,基于物理力学参数与地震荷载条件,通过相对变形法研究地震条件下深埋隧道围岩变形特性、能量聚集特性与安全性能。研究结果显示:深埋隧道围岩受地震荷载作用影响,形成挤让、椭圆化的变形走向;变形达到最高值时,应变能量密度高于2350 J/m^3的区域集中在隧道Ⅰ洞左侧围岩边墙和中夹岩柱上;各监测点安全系数呈非对称性,隧道右侧下角区域安全系数低于国家相关标准,计算结果与实际测量结果一致。     

Numerical Simulation of Displacement and Stress Fields Associated with the 1993 Kushiro-oki, Japan, Earthquake

—We constructed a three-dimensional finite element model to simulate coseismic and postseismic displacement and stress fields associated with the 1993 Kushiro-oki earthquake, which was a very large intermediate-depth earthquake that occurred within the subducted Pacific plate at a depth of 107 km beneath the southeastern part of Hokkaido, Japan. Taking the configuration of the subducted Pacific plate into account, we constructed a realistic model with lateral heterogeneity of viscoelastic structure. We assigned a variable slip distribution to the fault plane, which was obtained from inversion analysis of near-field seismic waveforms. The result shows that elastic deformation associated with the faulting reflects the assigned inhomogeneous slip distribution on the fault plane near the fault region, while it does not reflect the distribution on the free surface of the model. The calculated postseismic deformation does not reflect the slip distribution, but shows symmetric spatial patterns concerning the dipping direction of the fault both near the fault region and on the model surface. For the next 20 years following the earthquake, the amount of the calculated deformation is a fraction of the coseismic deformation. The calculated coeseismic deformation is large just above and below the fault plane, reaching 1 m, while the postseismic deformation is dominant near the upper and lower material boundaries between the subducted plate and the surrounding asthenosphere. The spatial distribution of maximum shear stress near the fault plane corresponds to the assigned slip distribution, amounting to 32 MPa. The directions of principal stress-change axes represent reverse fault type in the SSE region of the fault, whereas normal fault type is dominant in the NNW region with the exception of some asymmetrical spatial patterns of the principal stress-change axes on the fault due to the inhomogeneous slip distribution. Time variations both in the amount and the directions of stresses are minor, suggesting that the coseismic state of the stress would remain unchanged for two decades after the event.     

A stochastic model for tropical rainfall at a single location

Modelling data that correspond to rainfall accumulated over fixed periods of time presents the challenging problem of dealing with a random variable that has a point mass at zero which corresponds to dry periods that occur with positive probability. One way to overcome this difficulty is to assume that the data correspond to a normal variate w, that has been truncated and transformed. The dry periods correspond to the (unobserved) negative values and the wet periods correspond to some power of the positive ones. The serial structure that is present in rainfall can be modelled by imposing a serial structure to w. We use a dynamic linear model on w using a Fourier representation to allow for the seasonality of the data, which in the case of tropical rainfall is very marked. The model is fitted using a Markov chain Monte Carlo method that uses latent variables to handle both dry periods and missing values. We use the model to estimate and predict both the amount of rainfall and the probability of a dry period. The method is illustrated with data collected in the Venezuelan state of Guárico.     

Behaviour and modelling of partial‐strength beam‐to‐column composite joints for seismic applications

A refined component model is proposed to predict the inelastic monotonic response of exterior and interior beam‐to‐column joints for partial‐strength composite steel–concrete moment‐resisting frames. The joint typology is designed to exhibit ductile seismic response through plastic deformation developing simultaneously in the column web panel in shear, the bolted end‐plate connection, the column flanges in bending and the steel reinforcing bars in tension. The model can handle the large inelastic deformations consistent with high ductility moment‐resisting frames. Slip response between the concrete slab and the beams was taken into account. A fibre representation was adopted for the concrete slab to accurately capture the non‐uniform stress distribution and progressive crushing of the concrete at the interface between the concrete slab and the column flange. The model is validated against results from full‐scale subassemblages monotonic physical tests performed at the University of Pisa, Italy. A parametric study is presented to illustrate the capabilities of the model and the behaviour of the joints examined. Copyright © 2006 John Wiley & Sons, Ltd.     

Some methods of adjusting a terrestrial spatial network (Part 1)

Summary In the Part 1 and in a subsequent Part 2 to be published two methods of adjusting a spatial terrestrial network in tri-dimensional space are described. Care has been taken that the nature of the equations used, as well as of the adjustment, correspond to the same in adjustment satellite networks. The adjustment was carried out by the least-squares method according to conditioned observations. Various types of condition equations have been constructed according to the various types of adjusted quantities, and the various alternatives of the introduced errors (changes of input values) and weights. An effort was made to eliminate the ellipsoid of reference to the largest extent. The theory was applied numerically to a model of a smaller network which corresponds in position and height to usual triangulation networks with side lengths of about 30 km.Dedicated to 90th Birthday of Professor Frantiek Fiala     

Fault rock analysis of the northern part of the Chelungpu Fault and its relation to earthquake faulting of the 1999 Chi-Chi earthquake, Taiwan

Abstract   The 1999 Chi-Chi earthquake in Taiwan ( M _w = 7.6) produced a surface rupture along the north–south-striking Chelungpu thrust fault with pure dip-slip (east side up) and left lateral strike-slip displacements. Near-field strong-motion data for the northern part of the fault illustrate a distinct lack of the high-frequency seismic radiation associated with a large slip (10–15 m) and a rapid slip velocity (2–4 m/s), suggesting a smooth seismic slip associated with low dynamic frictional resistance on the fault. A drillhole was constructed at shallow depths in the possible fault zones of the northern part of the Chelungpu Fault, which may have slipped during the 1999 earthquake. One of the zones consists of a 20-cm-thick, unconsolidated fault breccia with a chaotic texture lacking both discrete slip surfaces (e.g. Riedel shears) and grain crushing. Other possible fault zones are marked by the narrow (less than a few centimeters) gouge zone in which clayey material intrudes into the damaged zone outside of the gouge zone. These characteristic fault rock textures suggest that the slip mechanisms at shallow levels during the earthquake involved either granular flow of initially unconsolidated material or slip localization under elevated pore pressure along the narrow clayey gouge zone. Because both mechanisms lead to low dynamic frictional resistance on the fault, the rapid seismic slip in the deep portions of the fault (i.e. the source region of strong-motion radiation) could have been accommodated by frictionless slip on the shallow portions of the fault. The combination of strong-motion data and fault rock analysis suggests that smooth slip associated with low dynamic friction occurred on both the deep and shallow portions of the fault, resulting in a large slip between the source region and the surface in the northern region.     

Loading Rate Dependence of Tensile Strength Anisotropy of Barre Granite

Granitic rocks usually exhibit strongly anisotropy due to pre-existing microcracks induced by long-term geological loadings. The understanding of the rock anisotropy in mechanical properties is critical to a variety of rock engineering applications. In this paper, Brazilian tests are conducted statically with a material testing machine and dynamically with a split Hopkinson pressure bar system to measure both static and dynamic tensile strength of Barre granite. To understand the anisotropy in tensile strength, samples are cored and labelled using the three principle directions of Barre granite to form six sample groups. For dynamic tests, a pulse shaping technique is used to achieve dynamic equilibrium in the samples during the dynamic test. The finite element method is then implemented to formulate equations that relate the failure load to the material tensile strength by employing an orthotropic elastic material model. For samples in the same orientation group, the tensile strength shows clear loading rate dependence. The tensile strengths also exhibit clear anisotropy under static loading while the anisotropy diminishes as the loading rate increases, which may be due to the interaction of pre-existing microcracks.     

Modelling the relationship between magnetic fabric and strain in polymineralic rocks

We present a model that is applicable to the relatively frequent case of rocks in which the magnetic fabric is dominated by uniaxial paramagnetic minerals, and in which the deformation (pure shear) corresponds to the March-Fernandez model. Borradaile et al.'s procedure for the isolation of the magnetic fabric of monocrystals of anisotropic components allows us to obtain the magnetic properties equivalent to those of all the component minerals if they should be artificially aligned. In the case of real component minerals with similar shape parameters, χ (dependent on dimension ratio), these properties will correspond to those of a theoretical equivalent mineral. Therefore, finite strain and orientation tensors may be determined from magnetic fabric measurements of such polymineralic rocks.     

含厚层软弱夹层硬质围岩的变形破坏特性

以某隧道为背景,通过地质条件分析并结合FLAC3D数值模拟手段,采用对比含有厚层软弱夹层和不含软弱夹层的2个不同模型的方法,从地应力变化、塑性区开展、变形、下覆岩体受压和位移等角度分析了含厚层软弱夹层硬质围岩的变形破坏特性。发现对于含厚层软弱夹层和不含软弱夹层的这2个不同模型,地应力的差别较小,但都为不含节理的岩体的0.5倍或更小,因此认为节理是影响含厚层软弱夹层硬质围岩变性破坏的因素之一;而2个模型在塑性区开展、变形量、下覆岩体受压和位移等方面差别较大,因此认为厚层软弱夹层对硬质围岩的变形破坏起着显著作用。同时通过创建三维模型采用分步开挖的方法,模拟了含厚层软弱夹层硬质围岩的变形破坏过程,分析认为其具有隐蔽性和发生时间不确定性的特点。研究结果可为类似地质情况的隧道建设提供参考。     

镶嵌试件表面位移场的全息显示

用激光全息干涉法研究了在单轴压力下的镶嵌试件的表面形变,显示了离面位移场的可见图象。发现在载荷达到破裂载荷的90%左右时,沿镶嵌边界出现剧烈形变带和宏观裂纹;当岩芯强度低于或等于外围介质强度时,这种剧烈形变带集中在岩芯以外的一侧;当岩芯强度高于外围介质强度时,这种剧烈形变带形成围空图象。镶嵌试件在一定程度上模拟了地壳的横向不均匀结构,观察它的表面形变过程有助于探讨震前地震活动图象的物理含义。     

RESEARCH PROGRESS ON LUMINESCENCE DATING OF ROCK SURFACES AND ITS APPLICATION

There are many examples of exposed or buried rock surfaces whose age is of interest to geologists and archaeologists. Luminescence dating is a well-established method of absolute chronology which has been successfully applied to a wide range of fine-grained sediments from hundreds of years to several hundred thousand years. Optical stimulated luminescence (OSL) has been recently proposed as a new method to date these rock surfaces (Laskariset al., 2011; Sohbatiet al., 2012a; Chapotet al., 2012; Pedersonet al., 2014; Sohbatiet al., 2015; Freieslebenet al., 2015). The basic principle is that luminescence signal of rock surface will soon decay when the rock is exposed to sunlight. When the rock surface is turned to be buried side, the OSL signal begins to accumulate again. With the variation of residual luminescence with depth, it is possible to estimate exposure and burial history of the rock. This article describes briefly the different luminescence dating methods for rock surfaces, its progress, application examples and present problems. For instance, Sohbatiet al. (2011) studied the depth dependence of the bleaching of the IRSL signal from granitic rocks, Laskaris and Liritzis (2011) proposed a mathematical function to describe the attenuation of daylight into rock surfaces, Sohbatiet al. (2012a) developed their model to include the environmental dose rate, Sohbatiet al. (2012b) overcame the problem of parameter estimation by using a known-age road-cut sample for calibration, Sohbatiet al. (2012c) further developed the OSL surface exposure dating model by including the simultaneous effect of daylight bleaching and environmental dose rate, and so on. These studies indicate that OSL dating method for rock surface can be applied widely to studies of geological and geomorphological evolution, archeology and Quaternary tectonic activity.     

Elastic Parameters of West Bohemian Granites under Hydrostatic Pressure

—The West Bohemian seismoactive region is situated near the contact of the Moldanu bian, Bohemian and Saxothuringian units in which a large volume is occupied by granitoid massifs. The spatial distribution of P-wave velocities and the rock fabric of five representative samples from these massifs were studied. The P-wave velocities were measured on spherical samples in 132 independent directions under hydrostatic pressure up to 400 MPa, using the pulse-transmission method. The pressure of 400 MPa corresponds to a depth of about 15 km in the area under study. The changes of P-wave velocity were correlated with the preferred orientations of the main rock fabric elements, i.e., rock forming minerals and microcracks. The values of the P-wave velocity from laboratory measurements on granite samples fit the velocity model used by seismologists in the West Bohemian seismoactive region.     

Experimental investigation of the strength and failure behavior of layered sandstone under uniaxial compression and Brazilian testing

As a typical inherently anisotropic rock, layered sandstones can differ from each other in several aspects, including grain size, type of material, type of cementation, and degree of compaction. An experimental study is essential to obtain and convictive evidence to characterize the mechanical behavior of such rock. In this paper, the mechanical behavior of a layered sandstone from Xuzhou, China, is investigated under uniaxial compression and Brazilian test conditions. The loading tests are conducted on 7 sets of bedding inclinations, which are defined as the angle between the bedding plane and horizontal direction. The uniaxial compression strength (UCS) and elastic modulus values show an undulatory variation when the bedding inclination increases. The overall trend of the UCS and elastic modulus values with bedding inclination is decreasing. The BTS value decreases with respect to the bedding inclination and the overall trend of it is approximating a linear variation. The 3D digital high-speed camera images reveal that the failure and fracture of a specimen are related to the surface deformation. Layered sandstone tested under uniaxial compression does not show a typical failure mode, although shear slip along the bedding plane occurs at high bedding inclinations. Strain gauge readings during the Brazilian tests indicate that the normal stress on the bedding plane transforms from compression to tension as the bedding inclination increases. The stress parallel to the bedding plane in a rock material transforms from tension to compression and agrees well with the fracture patterns; “central fractures” occur at bedding inclinations of 0°–75°, “layer activation” occurs at high bedding inclinations of 75°–90°, and a combination of the two occurs at 75°.     

Shaking table tests of 1:4 reduced‐scale models of masonry infilled reinforced concrete frame buildings

Two models of masonry infilled reinforced concrete frame buildings were tested at the shaking table. Models were built in the reduced scale 1:4 using the materials produced in accordance to modelling demands of true replica modelling technique. The first model represented a one‐storey box‐like building and the second one the two‐stories building with plan shaped in the form of a letter H. Models were shaken with the series of horizontal sine dwell motions with gradually increasing amplitude. Masonry infills of tested models were constructed of relatively strong bricks laid in weak mortar. Therefore, typical cracks developed and propagated along mortar beds without cracking of bricks or crushing of infill corners. Data collected from tests will be used in future evaluation, verification and development of computational models for prediction of in‐plane and out‐of‐plane behaviour of masonry infills. The responses of tested models can be well compared with global behaviour of real structures using the modelling rules. The similarity of local behaviour of structural elements, e.g. reinforced concrete joints, is less reliable due to limitations in modelling of steel reinforcement properties. The model responses showed that buildings designed according to Eurocodes are able to sustain relatively high dynamic excitations due to a significant level of structural overstrength. Copyright © 2001 John Wiley & Sons, Ltd.     

中朝地台东北缘3km深度地震层析成像的地质解释

利用本区丰富的地质,人工地震探测资料,对３ｋｍ深度层析成象成果进行了分析对比,结果表明：（１）３Ｋｍ深度速度值与人工地震探测所获地震波速值吻合较好;（２）３ｋｍ深度速度分布与该深度处的岩石介质分布有较好的对应关系,高速区对应地层主要为太古界－下元古界变质岩和地槽区的古生界浅变质岩,低速区对应主要为中生界火山,碎屑岩,燕山期花岗岩体和盖层性质的中元古－古代界沉积岩。（３）层析成象结果可划分深度不同的     

Comparison of seismic effects during deep tunnel excavation with different methods

The rapid release of strain energy is an important phenomenon leading to seismic events or rock failures during the excavation of deep rock.Through theoretical analysis of strain energy adjustment during blasting and mechanical excavation,and the interpretation of measured seismicity in the Jin-Ping Ⅱ Hydropower Station in China,this paper describes the characteristics of energy partition and induced seismicity corresponding to different energy release rates.The theoretical analysis indicates that part of the strain energy will be drastically released accompanied by violent crushing and fragmentation of rock under blast load,and this process will result in seismic events in addition to blasting vibration.The intensity of the seismicity induced by transient strain energy release highly depends on the unloading rate of in-situ stress.For mechanical excavation,the strain energy,which is mainly dissipated in the deformation of surrounding rock,releases smoothly,and almost no seismic events are produced in this gradual process.Field test reveals that the seismic energy transformed from the rock strain energy under high stress condition is roughly equal to that coming from explosive energy,and the two kinds of vibrations superimpose together to form the total blasting excavation-induced seismicity.In addition,the most intense seismicity is induced by the cut blasting delay; this delay contributes 50% of the total seismic energy released in a blast event.For mechanical excavation,the seismic energy of induced vibration(mainly the low intensity acoustic emission events or mechanical loading impacts),which accounts only for 1.5‰ of that caused by in-situ stress transient releasing,can be ignored in assessing the dynamic response of surrounding rock.     

循环荷载作用下粗粒料与掺砾粘土接触面单剪试验研究

土石坝内部存在较多的接触面,在静、动力荷载作用下,接触面的力学特性对坝体的稳定性具有显著的影响。基于动态单剪仪,开展了粗粒料与掺砾粘土接触面循环单剪试验。试验结果表明,初始剪应力加载方向对相对位移有较大影响：相对位移总是在沿初始剪应力加载方向一侧较大;在一个剪切循环内,正、反向加载时,若剪应力较小,发生剪缩变形,若剪应力较大,则发生剪胀变形;卸载时试样发生剪缩变形;整个剪切循环过程中,试样总体上处于剪缩状态;循环单剪过程中,接触面抗剪强度提高;接触面抗剪强度随着法向应力的增大而增大,呈线性关系。