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.     

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.     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

陕西白河县月儿潭绿松石与次生磷钙铝矾研究

对陕西省白河县月儿潭风化严重的绿松石样品进行了岩石薄片观察、化学成分测试、X射线粉晶衍射和红外吸收光谱分析等研究.结果表明,绿松石风化分解由天蓝色退色成黄绿色,最后变成白色,结构也由致密细腻变得松散多孔,主要原因是发生了硫酸根对磷酸根的交代作用,最终的交代产物主要是磷钙铝矾.     

饱和粉砂不稳定性的试验研究

通过对净砂和级配良好粉砂(含10 %粉土)进行一系列三轴固结不排水试验（CU）,研究了粉土、孔隙比和围压对饱和粉砂不稳定性的影响。试验结果表明,净砂与粉砂在不排水剪切条件下均会出现应变软化现象（即不稳定性）。同一围压下脆性指数(IB)随孔隙比增加,但不稳定线的应力比随孔隙比增加而减小。引用等效粒间孔隙比(ege)后,净砂和粉砂在ege-ln p?平面上拥有基本相同的临界状态线。在临界状态理论及等效粒间孔隙比的基础上,提出在同一修正状态参数(?ge)下净砂和级配良好粉砂有相似的不稳定性。     

Re-Os同位素对峨眉山大火成岩省成因制约的探讨

峨眉山大火成岩省（ELIP）主要由玄武岩、玄武质火山碎屑岩及少量的苦橄岩（包括越南的科马提岩）、长英质岩石以及层状岩体和岩墙组成,其物质来源直接关系到其成因是否与地幔柱活动有关。Re-Os同位素体系是地核、地幔和地壳物质的最佳示踪剂。前人对ELIP内的Re-Os同位素研究表明,低Ti玄武岩的Os含量为0．006×10^-9-0．40010^-9,^187Os／^188Os初始值为0．1371～1．403,并提出其与地幔柱活动有关;而高Ti玄武岩的Os含量为0．00410^-9～0．56010^-9,^187Os／^188Os初始值为0．1271～5．19,认为起源于大陆岩石圈地幔或地幔柱上升过程中受到大量岩石圈地幔“混染”（xu JF et al．,2007）;科马提岩的0s含量为1．2410^-9～7．0010^-9,^187Os／^188Os初始值为0．1251～0．1261,苦橄岩的Os含量为0．3210^-9～2．32910^-9,^187Os／^188Os初始值为0．1233～0．1266,指示苦橄岩和科马提岩均来自亏损地幔源区（Hanski et al．,2004;陈雷等,2007）。本文利用Os含量最低、^187Os／^188Os最高的高Ti玄武岩作为地壳端员,用铁质陨石、原始上地幔（PUM）和亏损地幔（DMM）作为地核和各种地幔端员,分别做二元混合计算,结果显示绝大多数玄武岩和所有苦橄岩及科马提岩均落在地壳和DMM混合曲线附近,并且邻区特提斯洋地幔岩与DMM具有相近的Os含量和^187Os／^188Os组成,据此推测峨眉山火成岩的形成与特提斯洋的活动有关,主要受控于地壳和亏损地幔的相互作用。     

塔里木盆地寒武-奥陶系白云岩结构构造类型及其形成机理

塔里木盆地寒武-奥陶系白云岩度大、分布广,是重要的油气储集层.根据岩石学研究并结合地球化学特征将其分为4类,其特点如下:(1)泥晶-粉晶白云岩,白云石晶体小、晶形差,有序度在4类中最低(0.73);Na、K、Sr含量在4类中最高;δ13C值及Z值较高,分别为-0.94‰、122.3;Fe、Mn含量较低;阴极发光为亮桔红色;常呈纹层状,具暴露及成化标志.这些特征表明其结晶较快,形成环境蒸发强烈、盐度较高、氧化性较强,是在萨布哈环境中由浓缩海水发生白云石化所形成.(2)藻纹层白云岩,白云石晶体小,晶形较好,有序度较高(0.85);藻纹层基本保留了其原始微细结构;Na、K、Sr含量较低:δC值在4类中最高,平均为-0.50‰;含有较多黄铁矿.这些特征说明其可能形成于成岩作用的稍早期阶段且其成因与藻类活动密切相关,成岩环境盐度较低,可能是在藻类作用下,由富藻层浓集的镁与残留海水在浅埋藏环境中共同作用所形成.(3)雾心亮边白云岩,白云石普遍具雾心亮边结构,雾心含少量包裹体及灰质残余,亮边则少见或无;阴极发光下,雾心发亮桔黄色光,亮边发光暗或不发光;局部见此类白云石;呈星散状分布于亮晶颗粒灰岩或泥晶灰岩中;白云石晶体较大,晶形较好.此类白云岩的δ13C及Z值在4类中最低,K、Na、Sr、Fe含量相对较低,Mn含量最高;有序度较高,均值为0.83.这些特征表明此类白云岩中白云石的雾心可能是在早期浅埋藏环境中,由富镁孔隙水交代周围灰质所形成;亮边则形成于埋藏较深的环境,其所需镁可能来源于粘土矿物的转化.(4)重结晶白云岩:白云石晶体粗大且污浊,接触紧密;局部见白云石晶体具环带结构,局部见颗粒幻影结构;阴极发光较暗;Fe含量在4类中最高,Mn含量较高;有序度最高,平均为0.95;包裹体均一温度平均为110.2℃.这些特征表明其可能是在还原性较强的深埋藏环境中,由早期形成的白云岩经较强的重结晶作用所形成.     

海拉尔盆地贝西斜坡北部地区储层特征及影响因素分析

通过测井曲线的标准化处理及信息提取、岩心及录井岩屑观察、铸体薄片、物性分析、扫描电镜、压汞分析等技术手段,研究了海拉尔盆地贝西斜坡北部地区南屯组储层的主要岩性特征、物性特征、储集空间类型和影响因素.研究结果表明.海拉尔盆地南屯组以内陆湖相碎屑岩为主,主要包括角砾岩、砾岩、砂砾岩、粗砂岩、细砂岩、粉砂岩、泥质砂岩及泥岩等.南一段储层孔隙度平均值为6.15%,渗透率平均值为0.31 × 10-3μm2,为低孔特低渗型孔隙特征;南二段储层孔隙度平均值为12.18%,渗透率平均值为2.79× 10-3μm2,属于中孔低渗型孔隙特征.储集空间类型以粒间孔隙为主,发育一定的次生孔隙.喉道分为4种类型(Ⅰ类、Ⅱ类、Ⅲ类、Ⅳ类).扇三角洲前缘分支水道砂体、扇三角洲平原辫状河道砂体和滨浅湖砂坝微相砂体为该区有利的储层砂体类型.储层性质主要受沉积相和成岩作用影响.     

岩石破坏前后曲线分类及脆-延转换围压研究——蚀变岩常规三轴压缩试验Ⅰ

蚀变岩是工程中少见的软弱岩类,在西南某重大水电工程中,蚀变岩处于工程的重要部位,为保证工程的长期稳定性,对蚀变岩力学特性进行了深入全面的试验研究.通过对孔隙度不同的饱水蚀变岩进行系统的常规三轴压缩试验和总结分析,提出了蚀变岩三轴压缩下破坏前应力.应变曲线可分为3大类,破坏后应力-应变曲线亦可分为3大类的形态模式.并得出结论:蚀变岩的破坏类型受围压与孔隙度的共同影响,在给定的12 MPa围压下蚀变岩以脆性破坏为主,只有孔隙度大于16%且围压大于4 MPa时才有可能进入脆-延转换状态,且脆-延转换围压随孔隙度增加而降低,临界状态应力比随孔隙度增大而增加.     

Reduction of hydroelastic responses of a very-long floating structure by a floating oscillating-water-column breakwater system

The hydroelastic responses of a very-long floating structure (VLFS) placed behind a reverse T-shape freely floating breakwater with a built-in oscillating water column (OWC) chamber are analyzed in two dimensions. The Bernoulli–Euler beam equation is coupled with the equations of rigid and elastic motions of the breakwater and the VLFS. The interaction of waves between the floating rigid breakwater and the elastic VLFS is formulated in a consistent manner. It has been shown numerically that the structural deflections of the VLFS can be reduced significantly by a suitably designed reverse T-shape floating breakwater.