Spatial probabilistic modeling of slope failure using an integrated GIS Monte Carlo simulation approach

Spatial probabilistic modeling of slope failure using a combined Geographic Information System (GIS), infinite-slope stability model and Monte Carlo simulation approach is proposed and applied in the landslide-prone area of Sasebo city, southern Japan. A digital elevation model (DEM) for the study area has been created at a scale of 1/2500. Calculated results of slope angle and slope aspect derived from the DEM are discussed. Through the spatial interpolation of the identified stream network, the thickness distribution of the colluvium above Tertiary strata is determined with precision. Finally, by integrating an infinite-slope stability model and Monte Carlo simulation with GIS, and applying spatial processing, a slope failure probability distribution map is obtained for the case of both low and high water levels.     

Pn arrivals and lateral variations of Moho geometry beneath the Kaapvaal craton

Pn arrivals from mining-induced earthquakes on the edge of the Witwatersrand basin show that the P wavespeeds in the uppermost mantle are almost constant throughout most of the Kaapvaal craton. The presence of only small wavespeed variations allows the use of a simple method of estimating crustal thicknesses below the stations of the Kaapvaal broad-band network using Pn times that has been compared with results from receiver functions. One thousand three hundred thirty-seven Pn arrivals were used to derive crustal thicknesses at 46 stations on the Kaapvaal craton. The average crustal thicknesses for 19 centrally located stations on each of the northern and southern regions of the craton that yielded well-constrained thicknesses were 50.52±0.88 km and 38.07±0.85 km, respectively. In contrast, the corresponding average thicknesses determined from receiver functions were 43.58±0.57 km and 37.58±0.70 km, respectively. The systematically lower values for receiver functions in the northern part of the Kaapvaal craton that was affected by the Bushveld magmatism at 2.05 Ga, suggest that the receiver functions do not enable the petrological crust mantle boundary to be reliably resolved due to variations in composition and metamorphic grade in a mafic lower crust. The Pn times also suggest pervasive azimuthal anisotropy with maximum wavespeeds of about 8.40 km/s at azimuths of about 15° and 217° in the northern and southern regions of the craton, respectively, and minimum wavespeeds of about 8.25 km/s.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

岩溶地区桩基的有限元分析

岩溶地区长桩存在着遇洞率的问题。为了研究岩溶地区桩的受力特点,采用有限元法模拟岩溶地区桩基的工作方式,分析了岩层厚度对桩的影响以及桩的受力状态,并研究了去桩的可能性。计算结果表明:岩层厚度对桩的影响较大,岩层的位置对桩的应力不产生影响,仅对桩的位移产生影响。     

溪洛渡水电站坝址区岩体蠕变特性试验研究

介绍了溪洛渡水电站坝址区岩体柔性板法蠕变试验方法及试验结果,给出了岩体蠕变模型及其参数,为工程设计提供了参考依据。     

饱水黄土隧道变形规律研究

新松树湾隧道位于饱和的砂粘土中,土体的稳定性极差,施工难度较大。为此,在隧道施工期间,通过土工试验、隧道收敛和围岩内部位移的现场监测、隧道变形的三维有限元仿真计算等手段,对该隧道的变形规律进行了系统研究,获得了大量试验数据。研究成果对隧道的施工组织和支护参数的确定起到了重要作用,为此类隧道的施工积累了经验。     

小湾电站高边坡系统锚固与排水的优化设计方法

据小湾电站地下厂房进水口开挖高边坡的岩体结构与工程地质特征,采用基于结构面网络模拟的随机楔体稳定分析方法对边坡的稳定性进行了三维分析。进一步运用风险分析理论对边坡系统锚杆与预应力锚索的布置进行了优化设计,最后运用结构面网络模拟理论对排水方案进行了优化设计。     

土钉支护结构的优化设计

根据复形优化原理,提出了基于Bishop法的土钉支护结构的优化分析方法,编写了优化分析程序。工程实例计算表明,该方法是进行土钉优化设计的有效方法。     

宜兴抽水蓄能电站试验洞的反分析研究

结合江苏宜兴抽水蓄能电站地下厂房试验洞量测资料的分析,研究了成层各向异性岩体中洞室工程问题的反分析方法。采用分层均质模型进行反分析计算。计算区域的各层岩层均被分别视为均质各向同性体,初始地应力假设呈线性分布。有限元分析中,对倾斜成层的岩层、软弱夹层及断层节理等分别划分了单元,据以模拟成层岩体的各向异性特征。为使计算过程简化,对离洞室较远的岩层拟按地质资料对材料特性参数赋值,反分析计算的目标未知数为初始地应力和试验洞相邻围岩的弹性模量。结果表明,与工程实践符合较好。     

南京地铁联络通道冻结法施工措施分析

对南京地铁一期工程TA4标联络通道冻结法施工的成功经验进行了总结。指出地铁联络通道冻结法施工中,必须采取必要的施工技术措施,才能保证施工安全、顺利的进行。本工程的成功经验可供其他工程参考。