Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

渤海湾盆地新生界生油岩系底界面温度分布

依据渤海湾盆地2000余口测温井的地温梯度数据、地层岩性描述、分层数据以及钻井资料,计算了该盆地各生油岩系底界面的温度。统计结果表明:渤海湾盆地沙河街组大部分凹陷区地层底界面温度介于90℃至150℃,目前仍具有大量生油的温度条件,而在隆起或一些凸起地区,该地层组段温度普遍小于90℃,未能达到生油的温度指标。东营组和孔店组地层也仍具有一定的生油温度条件。研究还表明:地层温度与地层界面埋深密切相关,温度随界面埋深的增大而升高,沉积厚度大的凹陷区地层界面温度大于沉积厚度小的凸起区或斜坡地带,说明地层界面埋深是决定地温高低的主要控制因素,而地温梯度对地层界面温度的影响相对较小。     

辽东湾凝析气田超压成因及侧向传递机制初探

位于辽东湾西部低凸起北端的JZ20-2超压凝析气田,离供烃中心较远,压力系数却达到1.56~1.7,这在我国东部裂谷盆地比较少见。经过综合分析,作者认为JZ20-2凝析气田异常高压的形成与其所处的特殊地质背景紧密相关,辽中凹陷沙河街组三段超压凝析油气沿不整合面、砂体输导通道向辽东湾西部低凸起的运移、充注和能量的侧向传递是JZ20-2凝析气田异常高压形成的根本原因,而上覆东营组二段下亚段-东营组三段巨厚的强超压泥岩封盖和辽西3号边界断层侧向遮挡共同构成的优越封闭环境是超压保存的必要条件。本文分析了该气田超压特征及其远距离超压侧向传递的成因机制,不仅有助于揭示辽东湾西部低凸起超压流体运移聚集的规律和成藏作用,而且可能提供了一个超压远距离侧向传递的典型实例。     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

Implementation of reconstructed geomorphologic units in landslide susceptibility mapping: the Melen Gorge (NW Turkey)

In the international literature, although considerable amount of publications on the landslide susceptibility mapping exist, geomorphology as a conditioning factor is still used in limited number of studies. Considering this factor, the purpose of this article paper is to implement the geomorphologic parameters derived by reconstructed topography in landslide susceptibility mapping. According to the method employed in this study, terrain is generalized by the contours passed through the convex slopes of the valleys that were formed by fluvial erosion. Therefore, slope conditions before landsliding can be obtained. The reconstructed morphometric and geomorphologic units are taken into account as a conditioning parameter when assessing landslide susceptibility. Two different data, one of which is obtained from the reconstructed DEM, have been employed to produce two landslide susceptibility maps. The binary logistic regression is used to develop landslide susceptibility maps for the Melen Gorge in the Northwestern part of Turkey. Due to the high correct classification percentages and spatial effectiveness of the maps, the landslide susceptibility map comprised the reconstructed morphometric parameters exhibits a better performance than the other. Five different datasets are selected randomly to apply proper sampling strategy for training. As a consequence of the analyses, the most proper outcomes are obtained from the dataset of the reconstructed topographical parameters and geomorphologic units, and lithological variables that are implemented together. Correct classification percentage and root mean square error (RMSE) values of the validation dataset are calculated as 86.28% and 0.35, respectively. Prediction capacity of the different datasets reveal that the landslide susceptibility map obtained from the reconstructed parameters has a higher prediction capacity than the other. Moreover, the landslide susceptibility map obtained from the reconstructed parameters produces logical results.     

Analysis of the 2003 Varunawat Landslide,Uttarkashi, India using Earth Observation data

Mass movements such as landslides in mountainous terrains are natural degradation processes and one of the most important landscape-building factors. Varunawat Parbat overlooking Uttarkashi town witnessed a series of landslides on 23 September 2003 and the debris slides and rock falls continued for 2 weeks. This landslide complex was triggered due to the incessant rainfall prior to the event, and its occurrence led to the blockage of the pilgrim route to Gangotri (source of the Ganges river) and evacuation of thousands of people to safer places. Though there was no loss of lives due to timely evacuation, heavy losses to the property were reported. High-resolution stereoscopic earth observation data were acquired after the incidence to study the landslide in detail with emphasis on the cause of the landslide and mode of failure. Areas along the road and below the Varunawat foothill region are mapped for landslide risk. It was found that the foothill region of the Varunawat Parbat was highly disturbed by man-made activities and houses are dangerously located below steep slopes. The potential zones for landslides along with the existing active and old landslides are mapped. These areas are critical and their treatment with priority is required in order to minimise further landslide occurrences.     

气候变化对塔里木河来自天山的地表径流影响

塔里木河水资源主要来自天山南坡两条源流,选择西段阿克苏河和中段开都河-孔雀河作为研究区.1956-2003年研究河源山区气温呈持续升温且降水波动增加的趋势,其中1995-2003年升温强劲,升温速率高出48 a期间平均的3倍以上;降水自1986年后持续增加,20世纪90年代较80年代增幅达18%,并显示出河源山区湿岛向塔里木盆地扩展.因高山缺少气象观测,出山径流过程变化可以综合反映中高山带的气候变化.塔里木河来自天山的地表径流在1986-2003年间持续增长,以冰川融水补给为主的库玛拉克河,1994年以来年径流量增加已在前期平均值基础上提升了一个台阶;开都河以降水径流补给为主,1986-2002年出现了观测记录以来的丰水期,并使1986年后博斯腾湖水位快速上升,恢复到1958年记录的最高水位以上.两河年径流变化趋势基本相似,但也显示有西、中段的气候变化局部差异,出现丰枯水期的不一致;然而,在近16 a升温过程中,年径流增长幅度和快慢相近.     

华北克拉通北缘隆化地区S型花岗岩的独居石年龄图谱

位于华北克拉通北缘中段的隆化S型花岗岩由石榴石黑云母花岗岩、石榴石花岗岩以及片麻理化的黑云母花岗岩组成。其主体岩性石榴石黑云母花岗岩SiO_2和Al_2O_3含量分别为64.09%～69.6%以及14.6%～16.13%,K_2O/Na_2O>1.0,A/CNK>1,0,Mg~#在20.76～34.89之间变化,具有明显的Nb、Ta、P、Ti和Sr亏损以及Rb、K和Th富集。石榴石黑云母花岗岩(样品JB6031-1)采用独居石电子探针U-Th-Pb化学法进行测年,获得了2553±120Ma、2180±42Ma和1854±24Ma三个年龄峰值。一颗独居石内部成分分带上6个分析点定年结果构成2553±120Ma的峰值年龄,这一年龄与我们最新获得的2506±7Ma和2541±8Ma(继承锆石年龄)LA-ICP-MS锆石U-Pb同位素年龄相似,我们将这一独居石年龄解释为继承独居石的年龄,表明在赤城-隆化断裂以北存在太古宙陆块,并且在后期构造-热事件中发生部分熔融形成S型花岗岩。该独居石颗粒幔部成分分带上10个分析点的测年结果揭示的峰值年龄为2181±42Ma,该年龄也是出现频率最高的年龄值,我们建议2181±42Ma为S型花岗岩的结晶年龄,反映了S型花岗岩的侵位时代。独居石颗粒外部成分分带上8个分析点的测年结果构成1854±24Ma的峰值年龄,该年龄与华北克拉通中部带的变质年龄接近,我们将其解释为S型花岗岩的变质年龄,表明华北克拉通北缘的构造演化与中部带的构造演化密切相关。     

三峡库区生基包滑坡监测及成果分析

生基包滑坡监测属于三峡库区奉节县三期地质灾害监测预警项目之一,该滑坡位于长江左岸,临近人口稠密的安坪乡集镇,地理位置重要。三峡水库175 m蓄水后,其变形破坏特征有何表现？对航道安全运营有无潜在威胁？是否会对滑坡体上的重要建筑及村民生产生活构成危害？针对这些问题,首先分析了滑坡的工程地质特征及主要的影响因素;其次,确立以4种监测手段为主、人工巡查为辅的监测方案;通过对大地变形GPS、深部位移、滑坡推力等几种监测方法的运用及对其成果进行分析研究,以实例说明其在滑坡监测中的应用;再次,结合宏观人工巡查进行对照分析,以充分说明大地变形、深部位移和滑坡推力监测在实际运用中的可行性;最后,根据监测结论提出对生基包滑坡防治的建议。     

基于群桩基础沉降估算确定压缩层厚度的探讨

在现行规范中,由于桩基沉降估算模式的不同,压缩层厚度的确定方法也不完全相同。对这些方法进行归纳、总结,并通过具体的工程实例,对由不同估算模式确定的压缩层厚度及沉降估算结果进行探讨与分析,认为对于工程场地地质资料掌握较为翔实的情况下,压缩层厚度可采用变形比法确定,否则,宜根据应力比法确定。